Van der Waals interaction and the stability of helical polypeptide chains

The characterization of Thermotoga maritima Arginine Binding Protein variants demonstrates that minimal local strains have an important impact on protein stability

The Ramachandran plot is a versatile and valuable tool that provides fundamental information for protein structure determination, prediction, and validation. The structural/thermodynamic effects produced by forcing a residue to adopt a conformation predicted to be forbidden were here explored using Thermotoga maritima Arginine Binding Protein (TmArgBP) as model. Specifically, we mutated TmArgBP Gly52 that assumes a conformation believed to be strictly disallowed for non-Gly residues. Surprisingly, the crystallographic characterization of Gly52Ala TmArgBP indicates that the structural context forces the residue to adopt a non-canonical conformation never observed in any of the high-medium resolution PDB structures. Interestingly, the inspection of this high resolution structure demonstrates that only minor alterations occur. Nevertheless, experiments indicate that Gly52 replacements in TmArgBP produce destabilizations comparable to those observed upon protein truncation or dissection in domains. Notably, we show that force-fields commonly used in computational biology do not reproduce this non-canonical state. Using TmArgBP as model system we here demonstrate that the structural context may force residues to adopt conformations believed to be strictly forbidden and that barely detectable alterations produce major destabilizations. Present findings highlight the role of subtle strains in governing protein stability. A full understanding of these phenomena is essential for an exhaustive comprehension of the factors regulating protein structures.

Peptide folding driven by Van der Waals interactions

Contrary to the widespread view that hydrogen bonding and its entropy effect play a dominant role in protein folding, folding into helical and hairpin-like structures is observed in molecular dynamics (MD) simulations without hydrogen bonding in the peptide-solvent system. In the widely used point charge model, hydrogen bonding is calculated as part of the interaction between atomic partial charges. It is removed from these simulations by setting atomic charges of the peptide and water to zero. Because of the structural difference between the peptide and water, van der Waals (VDW) interactions favor peptide intramolecular interactions and are a major contributing factor to the structural compactness. These compact structures are amino acid sequence dependent and closely resemble standard secondary structures, as a consequence of VDW interactions and covalent bonding constraints. Hydrogen bonding is a short range interaction and it locks the approximate structure into the specific secondary structure when it is included in the simulation. In contrast to standard molecular simulations where the total energy is dominated by charge-charge interactions, these simulation results will give us a new view of the folding mechanism.

Approximate treatment of the conformational characteristics of a cyclic nonapeptide, cyclolinopeptide a

Possible conformations for a cyclic nonapeptide that are consistent with conformation-dependent information obtained from an NMR investigation of the peptide in solution are presented. These several conformations are deduced from the myriad of possible conformations by eliminating from consideration all cyclic species having one or more residues in a conformation that does not correspond to the vicinal coupling constants observed by NMR between the amide and alpha-protons. A Karplus-like relation connecting the dihedral angle [unk]' and the vicinal coupling J(Nalpha) between N-H and C(alpha)-H(alpha) is used to test this correspondence. A further reduction in the number of cyclic conformations under consideration is made possible by rejecting the conformations that have a high intramolecular conformational energy. The intramolecular conformational energy of the cyclic nonapeptide is estimated by summing the independent residue energies. These have been calculated by others with approximate potential functions to account for the intrinsic torsional potentials and the nonbonded steric (6-12 potential) and electrostatic (monopole-monopole) interactions solely dependent upon one or both of the residue rotations, [unk] and Psi, about the N-C(alpha) and C(alpha)-C bonds, respectively.

Efficient minimization of angle-dependent potentials for polypeptides in internal coordinates

Angular potentials play an important role in the refinement of protein structures through angle-dependent restraints (e.g., those determined by cross-correlated relaxations, residual dipolar couplings, and hydrogen bonds). Analytic derivatives of such angular potentials with respect to the dihedral angles of proteins would be useful for optimizing such restraints and other types of angular potentials (i.e., such as we are now introducing into protein structure prediction) but have not been described. In this article, analytic derivatives are calculated for four types of angular potentials and integrated with the efficient recursive derivative calculation methods of Gō and coworkers. The formulas are implemented in publicly available software and illustrated by refining a low-resolution protein structure with idealized vector-angle, dipolar-coupling, and hydrogen-bond restraints. The method is now being used routinely to optimize hydrogen-bonding potentials in ROSETTA.

Flexible-geometry conformational energy maps for the amino acid residue preceding a proline

Previously calculated conformational energy maps suggest that the alpha-helical conformation for the residue preceding a proline is disfavored relative to the extended conformation by more than 7 kcal/mol. In known protein structures this conformation is observed, however, to occur for about 9% of all prolines. In addition, introduction or removal of prolines at theoretically unfavorable positions in proteins and peptides can have modest effects on stability and structure. To investigate the discrepancy between calculation and experiment, we have determined how the conformation of the proline affects the calculated energy. We have also explored the effect of bond length and bond angle relaxation on the conformational energy map. The conformational energy of the preceding residue is found to be unaffected by the conformation of the proline, but the effect of allowing covalent bond relaxation is dramatic. If bond lengths and angles, and dihedral angles within the pyrrolidine ring, are allowed to relax, a calculated energy difference between the alpha and beta conformations of 1.1 kcal/mol is obtained, in reasonable agreement with experiment. The detailed shape of the calculated energy surface is also in excellent agreement with the observed conformational distributions in known protein structures.

Theoretical prediction of the gel electrophoretic retardation changes due to point mutations in a tract of SV40 DNA

The changes of gel electrophoretic retardation due to single base substitutions in a 173 bp fragment of Sv40 DNA were predicted by using a theoretical model based on conformational energy calculations. As described in previous papers, this model allows successful prediction of the gel electrophoretic retardation of synthetic as well as natural DNAs reported in literature. The experimental retardations related to 195 point-mutated DNAs were reproduced with a standard deviation of 0.05 comparable with the experimental one of 0.04. This result, which represents a very critical test for the proposed model, indicates that DNA superstructures can be satisfactorily predicted on the simple physical basis of the integration of the nearest-neighbour perturbations in the dinucleotide steps. Thus, cooperative effects appear, in the majority of cases investigated, to play a second order role.

Protein conformational prediction

The prediction of the secondary and tertiary structure of globular and membrane proteins is reviewed. Prospects are encouraging for future developments, but present algorithms require cautious interpretation.

A theoretical model of DNA curvature

Distortions from the uniform idealized B-DNA structure are investigated in terms of differential interactions between adjacent nucleotide pairs on the basis of conformational energy calculations. A theoretical model of DNA curvature is proposed based on the evaluation of the curvature vector defined in the complex plane and the corresponding variance. The model appears to contain the basic physical features for translating the deterministic fluctuations of DNA sequences in superstructure elements. It allows the quantitative reproduction of all the available gel electrophoresis experiments on both periodical polynucleotides and tracts of DNAs as well as the theoretical prediction of the sequence dependent DNA writhing in good agreement with the experimental data. The general pattern of agreement between the theoretical and experimental data and the biological significance of the results obtained allow an extensive application of the model for the screening of DNA regions which are possible candidates for protein recognition.

Helix geometry in proteins

In this report we describe a general survey of all helices found in 57 of the known protein crystal structures, together with a detailed analysis of 48 alpha-helices found in 16 of the structures that are determined to high resolution. The survey of all helices reveals a total of 291 alpha-helices, 71 3(10)-helices and no examples of pi-helices. The conformations of the observed helices are significantly different from the "ideal" linear structures. The mean phi, psi angles for the alpha- and 3(10)-helices found in proteins are, respectively, (-62 degrees, -41 degrees) and (-71 degrees, -18 degrees). A computer program, HBEND, is used to characterize and to quantify the different types of helix distortion. alpha-Helices are classified as regular or irregular, linear, curved or kinked. Of the 48 alpha-helices analysed, only 15% are considered to be linear; 17% are kinked, and 58% are curved. The curvature of helices is caused by differences in the peptide hydrogen bonding on opposite faces of the helix, reflecting carbonyl-solvent/side-chain interactions for the exposed residues, and packing constraints for residues involved in the hydrophobic core. Kinked helices arise either as a result of included proline residues, or because of conflicting requirements for the optimal packing of the helix side-chains. In alpha-helices where there are kinks caused by proline residues, we show that the angle of kink is relatively constant (approximately 26 degrees), and that there is minimal disruption of the helix hydrogen bonding. The proline residues responsible for the kinks are highly conserved, suggesting that these distortions may be structurally/functionally important.

Studies on hydrogen bonds. Part V--Hydrogen bonding in energy minimization studies of peptides

An energy term, representing the N-H...O type of hydrogen bond, which is a function of the hydrogen bond length (R) and angle (theta) has been introduced in an energy minimization program, taking into consideration its interpolation with the non-bonded energy for borderline values of R and theta. The details of the mathematical formulation of the derivatives of the hydrogen bond function as applicable to the energy minimization have been given. The minimization technique has been applied to hydrogen bonded two and three linked peptide units (gamma-turns and beta-turns), and having Gly, Ala and Pro side chains. Some of the conformational highlights of the resulting minimum energy conformations are a) the occurrence of the expected 4----1 hydrogen bond in all of the burn-turn tripeptide sequences and b) the presence of an additional 3----1 hydrogen bond in some of the type I and II tripeptides with the hydrogen bonding scheme in such type I beta-turns occurring in a bifurcated form. These and other conformational features have been discussed in the light of experimental evidence and theoretical predictions of other workers.

Side chain characteristic main chain conformations of amino acid residues

Main chain conformation characteristics of the respective side chains of the 20 naturally occurring amino acid residues were obtained by the analysis of (phi, psi)-data from the crystal structures of 38 different globular proteins. The following observations are of interest: (i) For amino acids other than Glu, Thr and aliphatic amino acids, at least one main chain conformation is stabilised solely by side chain atoms. Such conformations are listed. (ii) In globular proteins, the main chain conformations which are significantly destabilised by side chain atom interactions are observed. The stabilising force for those conformations seems to come from main chain atom interactions. (iii) A set of conformations for which the main chain and side chain interactions are almost equal but in opposite directions, is listed and these conformations will be taken by residues mainly due to the effect of surroundings. The results can be used to study the folding of polypeptide chains and they also provide an insight into the role of the respective side chains on main chin conformations of amino acid residues.

Conformational studies on somatostatin. II. The C-terminal hexapeptide fragment

The results of a conformational study on the C terminal hexapeptide of Somatostatin are presented. Semi-empirical energy calculations and high resolution NMR methods have been used to obtain information on the conformational properties of SRIF9-14 in [2H6]dimethylsulfoxide and 2H2O. It is concluded from the energy calculations that the peptide has an averaged conformation in which semi extended and folded structures are important. Only some of the folded conformations can explain the chemical shift differences between the amino acid residues Thr10 and Thr12 as a ring current shift by the Phe11 aromatic ring on Thr10. The nonequivalence is more pronounced in dimethyl-sulfoxide (0.23--0.15 ppm) where it decreases with increasing temperature towards the temperature independent value in 2H2O (0.03 ppm). This suggests that the folded conformations are somewhat predominant in dimethylsulfoxide solutions. In 2H2O the semi extended and folded structures are statistically equally important and the peptide is more flexible. A comparison with a study on the smaller fragments SRIF10-12 and SRIF10-13 which have similar conformational properties, demonstrates the usefulness of the fragment approach in conformational studies of peptides.

The conformational properties of some fragments of the peptide hormone somatostatin

The conformational properties of some somatostatin fragments have been studied using high resolution NMR and semi-empirical calculations. The fragments are Thr10-Phe11, Phe11-Thr12, Thr10-Phe11-Thr12 and Thr10-Phe11-Thr12-Ser13. The results of high resolution 1H and 13C NMR using dimethylsulfoxide and 2H2O as solvents, combined with a new method for determining dihedral angles phi and psi from 13C and 1H spin lattice relaxation times are presented. A marked inequivalence of the two Thr10,12 residues is attributed to a shielding of the the Thr10 side chain by the Phe11 aromatic ring. The calculations show the existence of extended and folded low energy conformations in the tri and tetrapeptide. Only the folded structures give the observed shielding of Thr10. A temperature study of the tri and the tetrapeptide indicates that the folded structures are energetically the most favorable conformations at room temperature in dimethylsulfoxide. Increasing temperature reduces the nonequivalence of the Thr residues towards the differences that are observed in 2H2O. A detailed comparison of 3JalphaNH coupling constants and relaxation time measurements with the calculated conformations gives in general good agreement between both approaches. It is concluded that in these linear peptides, although several quite different low energy conformations exist, some of them are predominant. The continuity of both NMR parameters and calculated low energy conformations, when going from the smaller to the larger peptides, demonstrates the existence of structural properties far from the "random conformation".

A program for generation and selection of possible conformations of cyclic molecules

An algorithm has been developped for the computer generation of the possible conformations of non-planar cyclic molecules. This algorithm has been used successfully to select the known conformations of the ribose ring, of substituted cyclohexane and of disulfide bridged peptides. Its use may be extended to open molecules when the end-to-end distance is known or imposed.

Prediction of protein structure from amino acid sequence

In principle, it is possible to predict theoretically the three-dimensional structure of a protein from its amino acid sequence. Recently substantial progress towards this goal has been made by the use of simple models to represent protein conformation and interatomic interactions, together with the knowledge gained form analyses of know protein structures.