A study of the preferred environment of amino acid residues in globular proteins

In the native folded conformation of a globular protein, amino acid residues distant along the polypeptide chain come together to form the compact structure. This spatial structure is such that most of the polar residues are on the surface and have contact with the solvent medium and the nonpolar residues buried in the interior which have contact with similar nonpolar side chains. This cooperativity and mutual interaction among the randomly aligned amino acid residues suggest that each type of residue may prefer to have a specific environment. To gain more insight into this aspect of residue-residue cooperativity, a detailed analysis of the preferred environment associated with each of the 20 different amino acid residues in a number of protein crystals has been carried out. The variation of nonpolar nature computed for different sizes of spheres shows that the spatial region between radii of 6 and 8 Å is more favored for hydrophobic interactions and indicates that the influence of each residue over the surrounding medium extends predominantly up to a distance of 8 Å. The analysis of the surrounding amino acid residues associated with each type of residue shows that there is a definite tendency for each type of residue to have association with specific residues. The variation in environment is found even within the polar group as well as in the nonpolar group of residues. The surrounding residues associated with isoleucine, leucine, and valine are purely nonpolar. Proline, a nonpolar residue, is often surrounded by polar residues. The surrounding nonpolar nature of the tryptophan and tyrosine residues implies that even a single polar atom in a nonpolar side chain is sufficient to reduce their hydrophobic environment. There exists a high degree of mutual residue-residue cooperativity between the pairs glutamic acid-lysine, methionine-arginine, asparagine-tryptophan, and glutamine-proline, and the mutual residue-residue noncooperativity is high for the pairs methionine-aspartic acid, cysteine-glutamic acid, histidine-glutamine, and leucine-asparagine. The formation of secondary and tertiary structures is discussed in terms of the preferred environment and mutual cooperativity among various types of amino acid residues.     

残基相互作用网络特征与木聚糖酶耐热性的关系

残基相互作用网络是体现蛋白质中残基与残基之间协同和制约关系的重要形式。残基相互作用网络的拓扑性质以及社团结构与蛋白质的功能和性质有密切的关系。本文在构建一系列耐热木聚糖酶和常温木聚糖酶的残基相互作用网络后,通过计算网络的度、聚类系数、连接强度、特征路径长度、接近中心性、介数中心性等拓扑参数来确定网络拓扑结构与木聚糖酶耐热性的关系。识别残基相互作用网络的hub点,分析hub点的亲疏水性、带电性以及各种氨基酸在hub点中所占的比例。进一步使用GA-Net算法对网络进行社团划分,并计算社团的规模、直径和密度。网络的高平均度、高连接强度、以及更短的最短路径等表明耐热木聚糖酶残基相互作用网络的结构更加紧密;耐热木聚糖酶网络中的hub节点比常温木聚糖酶网络hub节点具有更多的疏水性残基,hub点中Phe、Ile、Val的占比更高。社团检测后发现,耐热木聚糖酶网络拥有更大的社团规模、较小的社团直径和较大的社团密度。社团规模越大表明耐热木聚糖酶的氨基酸残基更倾向于形成大的社团,而较小的社团直径和较大的社团密度则表明社团内部氨基酸残基的相互作用比常温木聚糖酶更强。     

残基相互作用网络特征与木聚糖酶耐热性的关系

残基相互作用网络是体现蛋白质中残基与残基之间协同和制约关系的重要形式。残基相互作用网络的拓扑性质以及社团结构与蛋白质的功能和性质有密切的关系。本文在构建一系列耐热木聚糖酶和常温木聚糖酶的残基相互作用网络后,通过计算网络的度、聚类系数、连接强度、特征路径长度、接近中心性、介数中心性等拓扑参数来确定网络拓扑结构与木聚糖酶耐热性的关系。识别残基相互作用网络的hub点,分析hub点的亲疏水性、带电性以及各种氨基酸在hub点中所占的比例。进一步使用GA-Net算法对网络进行社团划分,并计算社团的规模、直径和密度。网络的高平均度、高连接强度、以及更短的最短路径等表明耐热木聚糖酶残基相互作用网络的结构更加紧密;耐热木聚糖酶网络中的hub节点比常温木聚糖酶网络hub节点具有更多的疏水性残基,hub点中Phe、Ile、Val的占比更高。社团检测后发现,耐热木聚糖酶网络拥有更大的社团规模、较小的社团直径和较大的社团密度。社团规模越大表明耐热木聚糖酶的氨基酸残基更倾向于形成大的社团,而较小的社团直径和较大的社团密度则表明社团内部氨基酸残基的相互作用比常温木聚糖酶更强。     

Role of β-lactamase residues in a common interface for binding the structurally unrelated inhibitory proteins BLIP and BLIP-II

The β-lactamase inhibitory proteins (BLIPs) are a model system for examining molecular recognition in protein-protein interactions. BLIP and BLIP-II are structurally unrelated proteins that bind and inhibit TEM-1 β-lactamase. Both BLIPs share a common binding interface on TEM-1 and make contacts with many of the same TEM-1 surface residues. BLIP-II, however, binds TEM-1 over 150-fold tighter than BLIP despite the fact that it has fewer contact residues and a smaller binding interface. The role of eleven TEM-1 amino acid residues that contact both BLIP and BLIP-II was examined by alanine mutagenesis and determination of the association (k_on) and dissociation (k_off) rate constants for binding each partner. The substitutions had little impact on association rates and resulted in a wide range of dissociation rates as previously observed for substitutions on the BLIP side of the interface. The substitutions also had less effect on binding affinity for BLIP than BLIP-II. This is consistent with the high affinity and small binding interface of the TEM-1-BLIP-II complex, which predicts per residue contributions should be higher for TEM-1 binding to BLIP-II versus BLIP. Two TEM-1 residues (E104 and M129) were found to be hotspots for binding BLIP while five (L102, Y105, P107, K111, and M129) are hotspots for binding BLIP-II with only M129 as a common hotspot for both. Thus, although the same TEM-1 surface binds to both BLIP and BLIP-II, the distribution of binding energy on the surface is different for the two target proteins, that is, different binding strategies are employed.     

Serum albumin: Search for new sites of interaction with organophosphorus compounds by the example of soman

Albumin is known to be able to interact with organophosphorus compounds (OPCs), but neither amino acid residues of albumin that are responsible for this interaction nor the nature of the forming bonds have been finally established. Catalytic and pseudocatalytic functions of albumin are under consideration. Possible sites of interaction of albumin with soman have been elucidated by the methods of molecular modeling. Structures of the soman-albumin complexes have been determined by molecular docking. Stability of the obtained complexes has been evaluated by the method of molecular dynamics. The chemical bond between soman and the tyrosine-411 residue has been found to form only after deprotonation of the latter. The tyrosine-150 residue of albumin binds soman more effectively than tyrosine-411, and the tyrosine-150-deprotonation does not determine the efficacy of the binding (sorption) of soman, but affects the stability of the formed bound. It was proposed that the albumin residues of tyrosine-150 and serine-193 could serve as sites of the catalytic interaction with soman. We hypothesized that the deprotonation of an amino acid residue in one albumin site influenced initiation of the ligand binding in the other albumin site (allosteric albumin regulation).     

Functional mapping of the anti-idiotypic antibody anti-TS1 scFv using site-directed mutagenesis and kinetic analysis

Recombinant antibodies may be engineered to obtain improved functional properties. Functional mapping of the residues in the binding surfaces is of importance for predicting alterations needed to yield the desired properties. In this investigation, 17 single mutation mutant single-chain variable fragments (scFvs) of the anti-idiotypic antibody anti-TS1 were generated in order to functionally map amino acid residues important for the interaction with its idiotype TS1. Residues in anti-TS1 determined to be very important for the interaction were identified, Y32L, K50L, K33H, and Y52H, and they were distributed adjacent to a centrally located hydrophobic area, and contributed extensively to the interaction energy (≥2.5 kcal/mol) in the interaction. Quantitative ELISA assays, BIAcore technologies and three-dimensional surface analysis by modeling were employed to visualize the consequences of the mutations. The expression levels varied between 2 - 1,800 nM as determined by ELISA. All the 17 scFvs displayed higher dissociation rates (60 - 1,300 times) and all but two of them also faster association rates (1.3 - 56 times). The decrease in affinity was determined to be 1.6 - 12,200 times. Two of the mutants displayed almost identical affinity with the wild type anti-TS1, but with a change in both association and dissociation rates. The present investigation demonstrates that it is possible to generate a large panorama of anti-idiotypic antibodies, and single out a few that might be of potential use for future clearing and pre-targeting purposes of idiotypic-anti-idiotypic interactions.     

Construction of enzyme-substrate complexes between hen egg-white lysozyme and N-acetyl-D-glucosamine hexamer by systematic conformational search and molecular dynamics simulation

We constructed the complexes between HEWL and (GlcNAc)6 oligomer in order to investigate the amino acid residues related to substrate binding in the productive and nonproductive complexes, and the relationship between the distortion of the GlcNAc residue D and the formation of the productive complexes. We obtained 49 HEWL-(GlcNAc)6 complexes by a systematic conformational search and classified the each one to the three binding modes; left side, center, or right side. Furthermore we performed the molecular dynamics simulation against 20 HEWL-(GlcNAc)6 complexes (8: chair model, 12 : half-chair model) selected from the 49 complexes to investigate the interaction between HEWL and (GlcNAc)6. As results, we confirmed that it is necessary for GlcNAc residue D to be half-chaired form to bind toward the right side to form productive complexes. We found newly that eight amino acid residues interact with the (GlcNAc)6 oligomer, as follows, Arg73, Gly102, Asn103 for GlcNAc residue A; Asn103 for GlcNAc residues B and C; Leu56, Ala107, Val109 for GlcNAc residue D; Ala110 for GlcNAc residue E; and Lys33 for GlcNAc residue F. We also indicated that GlcNAc residue F does not interact with Thr47 and rarely interacts with Phe34 and Asn37.     

Analysis of beta-tubulin sequences reveals highly conserved, coordinated amino acid substitutions. Evidence that these 'hot spots' are directly involved in the conformational change required for dynamic instability

Vertebrate beta-tubulins have been classified into six classes on the basis of their C-terminal sequences [(1987) J. Cell Biol. 105, 1707-1720]. In particular, the sequences starting at residue 430 differ between isotypes of the same animal but are conserved between species. We extend this analysis and show that there are three 'hot spots', at residues 35, 55-57 and 124 which exhibit intra-species heterogeneity but inter-species conservation. There is a remarkable correlation between the identity of these residues and the C-terminal sequences, and suggests that the vertebrate beta-tubulins fall into three broad types. This correlation extends to those non-vertebrate organisms which have the Type 1 C-terminal sequence. We propose that these three 'hot spots' and the C-terminal peptide interact in the tertiary structure. We have also noted that the C-terminal peptide almost always contains a single phenylalanine or tyrosine residue, and that there is a strong correlation between this residue and the amino acids at positions 217/218, in both the vertebrate and non-vertebrate sequences. We propose that the C-terminal aromatic amino acid interacts with residues 217/218 in the tertiary structure. Analysis of conditions which stabilise microtubules and/or lower the steady state critical concentration strongly suggests that these two sets of coordinated amino acid substitutions are directly involved in effecting the conformational change associated with GTP hydrolysis which results in dynamic instability. We propose that there is an interaction between the highly acidic sequence between residue 430 and the aromatic amino acid (termed peptide A) and conserved basic amino acids located close to the 'hot spots'. We suggest that this interaction is altered in response to the assembly-dependent GTP hydrolysis, with the consequential increase in the subunit dissociation rate constant.     

Protein modeling combined with spectroscopic techniques: an attractive quick alternative to obtain structural information

Beside of the protein crystallography or NMR, another attractive option in protein structure analysis has recently appeared: computer modeling of the protein structure based on homology and similarity with proteins of already known structures. We have used the combination of computer modeling with spectroscopic techniques, such as steady-state or time-resolved fluorescence spectroscopy, and with molecular biology techniques. This method could provide useful structural information in the cases where crystal or NMR structure is not available. Molecular modeling of the ATP site within the H4-H5-loop revealed eight amino acids residues, namely besides the previously reported amino acids Asp443, Lys480, Lys501, Gly502 and Arg544, also Glu446, Phe475 and Gln482, which form the complete ATP recognition site. Moreover, we have proved that a hydrogen bond between Arg423 and Glu472 supports the connection of two opposite halves of the ATP-binding pocket. Similarly, the conserved residue Pro489 is important for the proper interaction of the third and fourth beta-strands, which both contain residues that take part in the ATP-binding. Alternatively, molecular dynamics simulation combined with dynamic fluorescence spectroscopy revealed that 14-3-3 zeta C-terminal stretch is directly involved in the interaction of 14-3-3 protein with the ligand. Phosphorylation at Thr232 induces a conformational change of the C-terminus, which is presumably responsible for observed inhibition of binding abilities. Phosphorylation at Thr232 induces more extended conformation of 14-3-3zeta C-terminal stretch and changes its interaction with the rest of the 14-3-3 molecule. This could explain negative regulatory effect of phosphorylation at Thr232 on 14-3-3 binding properties.     

Structure-wise discrimination of cytosine,thymine, and uracil by proteins in terms of their nonbonded interactions

The molecular recognition and discrimination of very similar ligand moieties by proteins are important subjects in protein–ligand interaction studies. Specificity in the recognition of molecules is determined by the arrangement of protein and ligand atoms in space. The three pyrimidine bases, viz. cytosine, thymine, and uracil, are structurally similar, but the proteins that bind to them are able to discriminate them and form interactions. Since nonbonded interactions are responsible for molecular recognition processes in biological systems, our work attempts to understand some of the underlying principles of such recognition of pyrimidine molecular structures by proteins. The preferences of the amino acid residues to contact the pyrimidine bases in terms of nonbonded interactions; amino acid residue–ligand atom preferences; main chain and side chain atom contributions of amino acid residues; and solvent-accessible surface area of ligand atoms when forming complexes are analyzed. Our analysis shows that the amino acid residues, tyrosine and phenyl alanine, are highly involved in the pyrimidine interactions. Arginine prefers contacts with the cytosine base. The similarities and differences that exist between the interactions of the amino acid residues with each of the three pyrimidine base atoms in our analysis provide insights that can be exploited in designing specific inhibitors competitive to the ligands.     

The DNA-bound orientation of Cu(II)-Xaa-Gly-His metallopeptides

The DNA-bound orientations of Cu(II) x Xaa-Gly-L-His metallopeptides (where Xaa is Gly, L-Lys or L-Arg) were investigated by DNA fiber EPR spectroscopy and molecular modeling. Observed and calculated EPR spectra indicated that the g// axes of 1:1 Cu(II) complexes of the tripeptides tilted about 50 degrees from the DNA fiber axis. These results suggest that the complexes are stereospecifically oriented in the DNA minor groove. Although the side chain of the N-terminal amino acid residue did not affect the orientation of the DNA-bound complexes, it contributed to their stability in the presence of DNA; the Cu(II) complex of Gly-Gly-L-His was found to dissociate to hydrated Cu(II) ion more extensively than the respective L-Lys-Gly-L-His and L-Arg-Gly-L-His complexes. The ionic interaction between the positively charged lysine or arginine residues and the negatively charged DNA phosphodiester backbone may result in the reduced dissociation of these complexes when bound to the DNA minor groove.     

Conformational fluctuations versus constraints in amino acid side chains: the evolution of information content from free amino acids to proteins

Like all other complex biological systems, proteins exhibit properties not found in free amino acids (i.e., emergent properties). Here, we explore top-down constraints experienced by the residue side chains in proteins compared to amino acids in increasingly complex molecular environments: free amino acids, end-capped amino acids, and the central residue in an alpha-helical nonapeptide. The crystalline structure of the contractile protein profilin Ib and the enzyme trypsin were chosen as objects of study, and submitted to 10 ns molecular dynamics (MD) simulations. The results revealed increased conformational constraints on the side chains when going from the simpler to the more complex compounds. A Shannon entropy (SE) analysis of the conformational behavior of the side chains showed in most cases a progressive and marked decrease in the SE of the chi1 and chi2 dihedral angles. This is equivalent to stating that conformational constraints on the side chain of residues increase their information content and, hence, recognition specificity compared to free amino acids. In other words, the vastly increased information content of a protein relative to its free monomers is embedded not only in the tertiary structure of the backbone, but also in the conformational behavior of the side chains. The postulated implication is that both backbone and side chains, by virtue of being conformationally constrained, contribute to the protein's recognition specificity toward other macromolecules and ligands.     

Analysis of catalytic residues in enzyme active sites

We present an analysis of the residues directly involved in catalysis in 178 enzyme active sites. Specific criteria were derived to define a catalytic residue, and used to create a catalytic residue dataset, which was then analysed in terms of properties including secondary structure, solvent accessibility, flexibility, conservation, quaternary structure and function. The results indicate the dominance of a small set of amino acid residues in catalysis and give a picture of a general active site environment. It is hoped that this information will provide a better understanding of the molecular mechanisms involved in catalysis and a heuristic basis for predicting catalytic residues in enzymes of unknown function.     

Functional mapping of the anti-idiotypic antibody anti-TS1 scFv using site-directed mutagenesis and kinetic analysis

Recombinant antibodies may be engineered to obtain improved functional properties. Functional mapping of the residues in the binding surfaces is of importance for predicting alterations needed to yield the desired properties. In this investigation, 17 single mutation mutant single-chain variable fragments (scFvs) of the anti-idiotypic antibody anti-TS1 were generated in order to functionally map amino acid residues important for the interaction with its idiotype TS1. Residues in anti-TS1 determined to be very important for the interaction were identified, Y32_L, K50_L, K33_H and Y52_H, and they were distributed adjacent to a centrally located hydrophobic area and contributed extensively to the interaction energy (≥2.5 kcal/mol) in the interaction. Quantitative ELISA assays, BIAcore technologies and three-dimensional surface analysis by modeling were employed to visualize the consequences of the mutations. The expression levels varied between 2–1,800 nM as determined by ELISA. All the 17 scFvs displayed higher dissociation rates (60–1,300 times) and all but two of them also displayed faster association rates (1.3–56 times). The decrease in affinity was determined to be 1.6–12,200 times. Two of the mutants displayed almost identical affinity with the wild-type anti-TS1, but with a change in both association and dissociation rates. The present investigation demonstrates that it is possible to generate a large panorama of anti-idiotypic antibodies and single out a few that might be of potential use for future clearing and pre-targeting purposes of idiotypic-anti-idiotypic interactions.     

Proteochemometrics analysis of substrate interactions with dengue virus NS3 proteases

The prime side specificity of dengue protease substrates was investigated by use of proteochemometrics, a technology for drug target interaction analysis. A set of 48 internally quenched peptides were designed using statistical molecular design (SMD) and assayed with proteases of four subtypes of dengue virus (DEN-1-4) for Michaelis (K(m)) and cleavage rate constants (k(cat)). The data were subjected to proteochemometrics modeling, concomitantly modeling all peptides on all the four dengue proteases, which yielded highly predictive models for both activities. Detailed analysis of the models then showed that considerably differing physico-chemical properties of amino acids contribute independently to the K(m) and k(cat) activities. For k(cat), only P1' and P2' prime side residues were important, while for K(m) all four prime side residues, P1'-P4', were important. The models could be used to identify amino acids for each P' substrate position that are favorable for, respectively, high substrate affinity and cleavage rate.     

To be folded or to be unfolded?

The lack of ordered structure in “natively unfolded” proteins raises a general question: Are there intrinsic properties of amino acid residues that are responsible for the absence of fixed structure at physiological conditions? In this article, we demonstrate that the competence of a protein to be folded or to be unfolded may be determined by the property of amino acid residues to form a sufficient number of contacts in a globular state. The expected average number of contacts per residue calculated from the amino acid sequence alone (using the average number of contacts for 20 amino acid residues in globular proteins) can be used as one of the simple indicators of natively unfolded proteins. The prediction accuracy for the sets of 80 folded and 90 natively unfolded proteins reaches 89% if the expected average number of contacts is used as a parameter and 83% in the case of hydrophobicity. An optimal set of artificial parameters for 20 amino acid residues obtained by Monte Carlo algorithm to maximally separate the sets of 90 natively unfolded and 80 folded proteins demonstrates the upper limit for prediction accuracy, which is 95%.     

Proline-specific endopeptidases from microbial sources: isolation of an enzyme from a Xanthomonas sp.

An extensive screening among microorganisms for the presence of post-proline-specific endopeptidase activity was performed. This activity was found among ordinary bacteria from soil samples but not among fungi and actinomycetes. This result is in contrast to the previous notion that this activity is confined to the genus Flavobacterium. A proline endopeptidase was isolated from a Xanthomonas sp. and characterized with respect to physicochemical and enzymatic properties. The enzyme is composed of a single peptide chain with a molecular weight of 75,000. The isoelectric point is 6.2. It is inhibited by diisopropylfluorophosphate and may therefore be classified as a serine endopeptidase. The activity profile is bell shaped with an optimum at pH 7.5. By using synthetic peptide substrates and intramolecular fluorescence quenching it was possible to study the influence of substrate structure on the rate of hydrolysis. The enzyme specifically hydrolyzed Pro-X peptide bonds. With Glu at position X, low rates of hydrolysis were observed; otherwise the enzyme exhibited little preference for particular amino acid residues at position X. A similar substrate preference was observed with respect to the amino acid residue preceding the prolyl residue in the substrate. The enzyme required a minimum of two amino acid residues toward the N terminus from the scissile bond, but further elongation of the peptide chain by up to six amino acid residues caused only a threefold increase in the rate of hydrolysis. Attempts to cleave at the prolyl residues in oxidized RNase failed, indicating that the enzyme does not hydrolyze long peptides, a peculiar property it shares with other proline-specific endopeptidases.     

Structural requirements of N-glycosylation of proteins. Studies with proline peptides as conformational probes

Conformational aspects of N-glycosylation have been investigated with a series of proline-containing peptides as molecular probes. The results demonstrate that, depending on the position of the imino acid in the peptide chain, dramatic alterations of glycosylation rates are produced, pointing to a critical contribution of the amino acids framing the 'marker sequence' triplet Asn-Xaa-Thr(Ser) on the formation of a potential sugar-attachment site. No glycosyl transfer at all was detectable to those peptides containing a proline residue either in position Xaa or in the next position beyond the threonine of the Asn-sequon on the C-terminal side, whereas the hexapeptide Pro-Asn-Gly-Thr-Ala-Val was glycosylated at a high rate. (Emboldened residues denote the 'marker sequence' that is identical in all the peptides; italicized residues distinguish the positions of proline in the various peptides.) Studies with space-filling models reveal that the lack of glycosyl-acceptor capabilities of Ala(Pro)-Asn-Gly-Thr-Pro-Val might be directly related to their inability to adopt and/or stabilize a turn or loop conformation which permits the catalytically essential interaction between the hydroxy amino acid and the asparagine residue within the 'marker sequence' [Bause & Legler (1981) Biochem. J. 195, 639-644]. This conclusion is supported by circular-dichroism spectroscopic data, which suggest structure-forming potentials in this type of non-acceptor peptides dominating over those that favour the induction of an appropriate sugar-attachment site in the acceptor peptides. The lack of acceptor properties of Tyr-Asn-Pro-Thr-Ser-Val indicates that even small modifications in the 'recognition' pattern are not tolerated by the N-glycosyltransferases.     

Mapping Interactions between Myosin Relay and Converter Domains That Power Muscle Function

Intramolecular communication within myosin is essential for its function as motor, but the specific amino acid residue interactions required are unexplored within muscle cells. Using Drosophila melanogaster skeletal muscle myosin, we performed a novel in vivo molecular suppression analysis to define the importance of three relay loop amino acid residues (Ile⁵⁰⁸, Asn⁵⁰⁹, and Asp⁵¹¹) in communicating with converter domain residue Arg⁷⁵⁹. We found that the N509K relay mutation suppressed defects in myosin ATPase, in vitro motility, myofibril stability, and muscle function associated with the R759E converter mutation. Through molecular modeling, we define a mechanism for this interaction and suggest why the I508K and D511K relay mutations fail to suppress R759E. Interestingly, I508K disabled motor function and myofibril assembly, suggesting that productive relay-converter interaction is essential for both processes. We conclude that the putative relay-converter interaction mediated by myosin residues 509 and 759 is critical for the biochemical and biophysical function of skeletal muscle myosin and the normal ultrastructural and mechanical properties of muscle.