Analysis of phi and chi 1 torsion angles for hen lysozyme in solution from 1H NMR spin-spin coupling constants

Three-bond 3JHN alpha coupling constants have been determined for 106 residues and 3J alpha beta coupling constants have been measured for 57 residues of the 129-residue protein hen egg white lysozyme. These NMR data have been compared with torsion angles defined in the tetragonal and the triclinic crystal forms of the protein. For most residues the measured 3JHN alpha values were consistent with the phi torsion angles found in both crystal forms; the RMS difference between the coupling constants calculated by using the tetragonal crystal structure phi angles and the experimental 3JHN alpha values is 0.88 Hz. Thus there appears to be no significant averaging of the phi torsion angle either in the interior or at the surface of the protein. For 41 of the residues where 3J alpha beta coupling constants have been determined, the values are consistent with a single staggered conformation about the chi 1 torsion angle and there is complete agreement between the NMR data in solution and the torsion angles defined in the crystalline state. In contrast, for the other 16 residues where 3J alpha beta coupling constant values have been measured, the data indicate extensive motional averaging about the chi 1 torsion angle. These residues occur largely on the surface of the protein and examination of the crystal structures shows that many of these residues adopt a different conformation in the triclinic and tetragonal crystal forms and have high crystallographic temperature factors. It appears, however, that in solution conformational flexibility of the side chains of surface residues is significantly more pronounced than in individual crystal structures.     

Graphical representation of the salient conformational features of protein residues.

A composite plot for depicting in two dimensions the conformation and the secondary structural features of protein residues has been developed. Instead of presenting the exact values of the main- and side-chain torsion angles (φ, psi and chi(1)), it indicates the region in the three-dimensional conformational space to which a residue belongs. Other structural aspects, like the presence of a cis peptide bond and disulfide linkages, are also displayed. The plot may be used to recognize patterns in the backbone and side-chain conformation along a polypeptide chain and to compare protein structures derived from X-ray crystallography, NMR spectroscopy or molecular modelling studies and also to highlight the effect of mutation on structure.     

CAP: conformation angles package-displaying the conformation angles of side chains in proteins

SUMMARY: A graphics package has been developed to display all side chain conformation angles of the user selected residue in a given protein structure. The proposed package is incorporated with all the protein structures (solved using X-ray diffraction and NMR spectroscopy) available in the Protein Data Bank. The package displays the multiple conformations adopted by a single amino acid residue whose structure is solved and refined at very high resolution. In addition, it shows the percentage distribution of the side chain conformation angles in different rotameric states. AVAILABILITY: http://144.16.71.146/cap/     

CADB: Conformation Angles DataBase of proteins

Conformation Angles DataBase (CADB) provides an online resource to access data on conformation angles (both main-chain and side-chain) of protein structures in two data sets corresponding to 25% and 90% sequence identity between any two proteins, available in the Protein Data Bank. In addition, the database contains the necessary crystallographic parameters. The package has several flexible options and display facilities to visualize the main-chain and side-chain conformation angles for a particular amino acid residue. The package can also be used to study the interrelationship between the main-chain and side-chain conformation angles. A web based JAVA graphics interface has been deployed to display the user interested information on the client machine. The database is being updated at regular intervals and can be accessed over the World Wide Web interface at the following URL: http://144.16.71.148/cadb/.     

Conformational properties of biologically active melanostatin analogs

On the basis of the two-dimensional 1H NMR studies the conformation of melanocyte inhibiting hormone (HCl.Pro-Leu-Gly-NH2, MIF) and its five analogues with p-substituted phenylalanine has been determined. Structure-antidepressant activity relationship (examined by the Porsolt test) of MIF and its analogues has been estimated by means of the multivariate statistical analysis, the vicinal spin coupling constants, which determine phi and chi dihedral angles of the second amino acid and phi and dihedral angle of glycine, being selected as structural parameters. It is shown that a biologically active conformation (10-membered beta-turn II) is realized for the considered peptides.     

Statistical and energetic analysis of side-chain conformations in oligopeptides

The distributions of side-chain conformations in 258 crystal structures of oligopeptides have been analyzed. The sample contains 321 residues having side chains that extend beyond the C beta atom. Statistically observed preferences of side-chain dihedral angles are summarized and correlated with stereochemical and energetic constraints. The distributions are compared with observed distributions in proteins of known X-ray structures and with computed minimum-energy conformations of amino acid derivatives. The distributions are similar in all three sets of data, and they appear to be governed primarily by intraresidue interactions. In side chains with no beta-branching, the most important interactions that determine chi 1 are those between the C gamma H2 group and atoms of the neighboring peptide groups. As a result, the g- conformation (chi 1 congruent to -60 degrees) occurs most frequently for rotation around the C alpha-C beta bond in oligopeptides, followed by the t conformation (chi 1 congruent to 180 degrees), while the g+ conformation (chi 1 congruent to 60 degrees) is least favored. In residues with beta-branching, steric repulsions between the C gamma H2 or C gamma H3 groups and backbone atoms govern the distribution of chi 1. The extended (t) conformation is highly favored for rotation around the C beta-C gamma and C gamma-C delta bonds in unbranched side chains, because the t conformer has a lower energy than the g+ and g- conformers in hydrocarbon chains. This study of the observed side-chain conformations has led to a refinement of one of the energy parameters used in empirical conformational energy computations.     

Analysis of the relationship between side-chain conformation and secondary structure in globular proteins

The relationship between the preferred side-chain dihedral angles and the secondary structure of a residue was examined. The structures of 61 proteins solved to a resolution of 2.0 A (1 A = 0.1 nm) or better were analysed using a relational database to store the information. The strongest feature observed was that the chi 1 distribution for most side-chains in an alpha-helix showed an absence of the g- conformation and a shift towards the t conformation when compared to the non-alpha/beta structures. The exceptions to this tendency were for short polar side-chains that form hydrogen bonds with the main-chain which prefer g+. Shifts in the chi 1 preferences for residues in the beta-sheet were observed. Other side-chain dihedral angles (chi 2, chi 3, chi 4) were found to be influenced by the main-chain. This paper presents more accurate distributions for the side-chain dihedral angles which were obtained from the increased number of proteins determined to high resolution. The means and standard deviations for chi 1 and chi 2 angles are presented for all residues according to the secondary structure of the main-chain. The means and standard deviations are given for the most popular conformations for side-chains in which chi 3 and chi 4 rotations affect the position of C atoms.     

Ramachandran plot on the web (2.0)

The Ramachandran plot displays the main chain conformation angles (Phi and Psi) of the polypeptide chain of a protein molecule. The paper reports the updated version of the Ramachandran plot web server and has several improved options for displaying the conformation angles in various regions. In addition, options are provided to display the conformation angles in various secondary structural elements and regions within the user specified Phi and Psi values in the plot. The updated version is accessible at the following URL: http://dicsoft1.physics.iisc.ernet.in/rp/.     

An algorithm for determining the conformation of polypeptide segments in proteins by systematic search

The feasibility of determining the conformation of segments of a polypeptide chain up to six residues in length in globular proteins by means of a systematic search through the possible conformations has been investigated. Trial conformations are generated by using representative sets of phi, psi, and chi angles that have been derived from an examination of the distributions of these angles in refined protein structures. A set of filters based on simple rules that protein structures obey is used to reduce the number of conformations to a manageable total. The most important filters are the maintenance of chain integrity and the avoidance of too-short van der Waals contacts with the rest of the protein and with other portions of the segment under construction. The procedure is intended to be used with approximate models so that allowance is made throughout for errors in the rest of the structure. All possible main chains are first constructed and then all possible side-chain conformations are built onto each of these. The electrostatic energy, including a solvent screening term, and the exposed hydrophobic area are evaluated for each accepted conformation. The method has been tested on two segments of chain in the trypsin like enzyme from Streptomyces griseus. It is found that there is a wide spread of energies among the accepted conformations, and the lowest energy ones have satisfactorily small root mean square deviations from the X-ray structure.     

Serine and threonine residues bend alpha-helices in the chi(1) = g(-) conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.     

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.     

Geometry of interaction of metal ions with sulfur-containing ligands in protein structures

An analysis of the geometry of binding of metal ions by cysteine and methionine residues in protein structures has been made by using the Protein Data Bank. Metal ions have a distinct model of binding to each of these residues, and this is independent of the nature of the metal center or the type of protein. Metal ions tend to approach the sulfur of Met roughly 38 degrees from the perpendicular to the plane through atoms C gamma-S delta-C epsilon. For Cys, the approach direction is such that the M...S gamma-C beta-C alpha torsional angle is about +/- 90 or 180 degrees. The side-chain conformation of the cysteine residue is affected by the presence of the metal ion; there is a shift from the g+ conformation toward g- and mainly t conformations. When two Cys residues at positions i-3 and i bind to the same metal center, there appears to be some restriction on the geometry of metal binding by the residue i; for such a residue chi 1 and M...S gamma-C beta-C alpha angles are likely to be around 60 degrees and 270 degrees, respectively. Met and Cys residues coordinating to a metal ion are usually from coil or turn regions of the protein structure.     

Method of determining protein conformations by the two-dimensional nuclear Overhauser enhancement spectroscopy data

A method is suggested to determine the most probable values of the angles phi, psi of the protein backbone by the data on the availability and absence of d connectivities in the two-dimensional nuclear Overhauser enhancement spectra. In view of this, the dependences of the proton-proton distances in dipeptide units of L-amino acid residues on the dihedral angles phi, psi, chi1 are considered and the conformational states of amino acid residues of the proteins with the known spatial structure are analysed statistically. The potentialities of the method are assessed with the aid of model spectral nuclear magnetic resonance (NMR) parameters obtained from the X-ray data for the bovine pancreatic trypsin inhibitor and avian pancreatic polypeptide. It is shown that the developed procedure of structural interpretation of the NMR data allows one to correctly reproduce the local conformation of the protein backbone. The obtained backbone conformation may serve as a starting point to build and refine molecular three-dimensional structure.     

Minimum energy conformations of proline-containing helices.

Proline occurs frequently in transmembrane alpha-helices of transport and receptor proteins even though statistical surveys demonstrate the overwhelming preference of this residue for a non-alpha-helical, hydrophilic environment. As a result, membrane-buried proline has been proposed to be functionally important, with function arising from structural discontinuity or destabilization of the helix. Destabilization may occur by Pro-mediated conformational transitions between discrete states, and may be manifested in membrane protein systems through reversible processes such as channel opening and closing or signal transduction. In this study, computer modeling of a model transmembrane alpha-helix, (Ala)8-Leu-Pro-Phe-(Ala)8, in a medium of low polarity (dielectric = 2), is used to examine the occurrence and energetic accessibility of Pro-mediated conformational interconversions. Leu psi and chi 1, Pro psi, and Phe phi and chi 1 torsion angles were assigned random values so that a data base of 200 conformations for each of the cis and trans states was generated. The conformations were minimized and low-energy structures organized into families. This analysis demonstrated that the most populated lowest energy family is the Trans-I conformation, corresponding to proline in a kinked alpha-helix. Two additional trans structures, Trans-II and Trans-III, as well as a cis conformation, Cis-I, are also energetically competitive. Interconversions between the trans states could thus be mediated by changes at a single torsion angle, accompanied by minor local hydrogen-bonding rearrangements. This work substantiates that membrane-buried proline can provide the basis for conformational transitions between discrete alpha-helix-based structures in a nonpolar environment.     

Side-chain structure and dynamics at the lipid-protein interface: Val1 of the gramicidin A channel.

High resolution dynamics and structural information has been resolved from 2H solid-state NMR spectra of the Val-1 side-chain of the gramicidin channel in a lipid bilayer. Both powder pattern lineshapes and spectra from uniformly aligned samples of gramicidin in lipid bilayers have been analyzed to achieve a fully consistant interpretation of the data. Torsional motions about the C alpha C beta axis (chi 1) are shown to be three-state jumps in which the occupancy of the states is given by the ratio, 75:15:10 for the chi 1 angles of 184 degrees:304 degrees:64 degrees. The dominant conformer is also the most common conformation observed for valines in well defined protein structures. The distribution of conformational substates that represents the chi 1 dynamics appears to be largely independent of the lipid phase transition and the hydration of the sample. However, there is evidence that the residence time between jumps is dependent on the lipid phase transition. Although this time is shown to be approximately 1 microseconds below the phase transition temperature, it is in the fast exchange limit above the transition temperature.     

Conformational Angles DataBase (CADB-3.0)

Transitions in amino-acid conformation angles tend to accompany various structural modifications in protein structures. Thus, to benefit the modeling of protein structures, the Conformation Angles DataBase (CADB-3.0) has been updated to visualize the conformational angles in varied regions (fully, generously, additionally and disallowed regions). In addition, options are provided to display the angles in the secondary structural elements (alpha-helix, beta-sheet and 3(10)-helix) of the Ramachandran plot. The database is being updated periodically and can be accessed over the World Wide Web at the following URL: http://cluster.physics.iisc.ernet.in/cadb/.     

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.     

Extending the accuracy limits of prediction for side-chain conformations

Current techniques for the prediction of side-chain conformations on a fixed backbone have an accuracy limit of about 1.0-1.5 A rmsd for core residues. We have carried out a detailed and systematic analysis of the factors that influence the prediction of side-chain conformation and, on this basis, have succeeded in extending the limits of side-chain prediction for core residues to about 0.7 A rmsd from native, and 94 % and 89 % of chi(1) and chi(1+2 ) dihedral angles correctly predicted to within 20 degrees of native, respectively. These results are obtained using a force-field that accounts for only van der Waals interactions and torsional potentials. Prediction accuracy is strongly dependent on the rotamer library used. That is, a complete and detailed rotamer library is essential. The greatest accuracy was obtained with an extensive rotamer library, containing over 7560 members, in which bond lengths and bond angles were taken from the database rather than simply assuming idealized values. Perhaps the most surprising finding is that the combinatorial problem normally associated with the prediction of the side-chain conformation does not appear to be important. This conclusion is based on the fact that the prediction of the conformation of a single side-chain with all others fixed in their native conformations is only slightly more accurate than the simultaneous prediction of all side-chain dihedral angles.     

Excluded volume in protein side-chain packing

The excluded volume occupied by protein side-chains and the requirement of high packing density in the protein interior should severely limit the number of side-chain conformations compatible with a given native backbone. To examine the relationship between side-chain geometry and side-chain packing, we use an all-atom Monte Carlo simulation to sample the large space of side-chain conformations. We study three models of excluded volume and use umbrella sampling to effectively explore the entire space. We find that while excluded volume constraints reduce the size of conformational space by many orders of magnitude, the number of allowed conformations is still large. An average repacked conformation has 20 % of its chi angles in a non-native state, a marked reduction from the expected 67 % in the absence of excluded volume. Interestingly, well-packed conformations with up to 50 % non-native chi angles exist. The repacked conformations have native packing density as measured by a standard Voronoi procedure. Entropy is distributed non-uniformly over positions, and we partially explain the observed distribution using rotamer probabilities derived from the Protein Data Bank database. In several cases, native rotamers that occur infrequently in the database are seen with high probability in our simulation, indicating that sequence-specific excluded volume interactions can stabilize rotamers that are rare for a given backbone. In spite of our finding that 65 % of the native rotamers and 85 % of chi(1) angles can be predicted correctly on the basis of excluded volume only, 95 % of positions can accommodate more than one rotamer in simulation. We estimate that, in order to quench the side-chain entropy observed in the presence of excluded volume interactions, other interactions (hydrophobic, polar, electrostatic) must provide an additional stabilization of at least 0.6 kT per residue in order to single out the native state.     

1H-NMR studies on nucleotide binding to the sarcoplasmic reticulum Ca2+ ATPase. Determination of the conformations of bound nucleotides by the measurement of proton-proton transferred nuclear Overhauser enhancements

The glycosidic bond torsion angles and the conformations of the ribose of Mg2+ATP, Mg2+ADP and Mg2+AdoPP[NH]P (magnesium adenosine 5'-[beta, gamma-imido]triphosphate) bound to Ca2+ATPase, both native and modified with fluorescein isothiocyanate (FITC), in intact sarcoplasmic reticulum have been determined by the measurement of proton-proton transferred nuclear Overhauser enhancements by 1H-NMR spectroscopy. This method shows clearly the existence of a low-affinity ATP binding site after modification of the high-affinity site with FITC. For all three nucleotides bound to both the high-affinity (catalytic) site and the low-affinity site, we find that the conformation about the glycosidic bond is anti, the conformation of the ribose 3'-endo of the N type and the conformation about the ribose C4'-C5' bond either gauche-trans or trans-gauche. The values for the glycosidic bond torsion angles chi (O4'-C1'-N9-C4) for Mg2+ATP, Mg2+ADP and Mg2+AdoPP[NH]P bound to the low-affinity site of FITC-modified Ca2+ATPase are approximately equal to 270 degrees, approximately equal to 260 degrees and approximately equal to 240 degrees respectively. In the case of the nucleotides bound to the high-affinity (catalytic) site of native Ca2+ATPase, chi lies in the range 240-280 degrees.