Skip to main content
  • Methodology article
  • Open access
  • Published:

Stereochemical errors and their implications for molecular dynamics simulations

Abstract

Background

Biological molecules are often asymmetric with respect to stereochemistry, and correct stereochemistry is essential to their function. Molecular dynamics simulations of biomolecules have increasingly become an integral part of biophysical research. However, stereochemical errors in biomolecular structures can have a dramatic impact on the results of simulations.

Results

Here we illustrate the effects that chirality and peptide bond configuration flips may have on the secondary structure of proteins throughout a simulation. We also analyze the most common sources of stereochemical errors in biomolecular structures and present software tools to identify, correct, and prevent stereochemical errors in molecular dynamics simulations of biomolecules.

Conclusions

Use of the tools presented here should become a standard step in the preparation of biomolecular simulations and in the generation of predicted structural models for proteins and nucleic acids.

Background

Biomolecules often feature asymmetries in stereochemistry. Many biologically active molecules are chiral, i.e., they exist in two forms, called enantiomers, which are non-superimposable mirror images of each other. Of particular relevance to biological compounds is the carbon atom as a chiral center: a carbon atom is chiral if it carries four nonequivalent substituents. Thus, all amino acids save glycine have at least one chiral center at C α (see Figure 1A). Threonine and isoleucine have an additional chiral center at C β . Interestingly, only one of the two enantiomers is widely used in nature: according to the D-/L-naming convention, most naturally occurring amino acids are found in the L-configuration. Note, however, that D-amino acids do occur in biology, e.g., in cell walls of bacteria [1–3]. Nucleic acids also have chiral centers. For example, in DNA the atoms C1', C3', and C4' of the sugar moiety are chiral, while in RNA the presence of an additional OH group renders also C2' of the ribose chiral (see Figure 1B). Although the origin of the homochirality is not understood, the asymmetry due to preferential use of one enantiomer in biological systems has wide consequences. In particular, recognition processes of chiral molecules are impacted as can be demonstrated, e.g., by the different smell of D- and L-carvone or by the inhibition of proteases by D-amino acids [4]. On the level of protein structure, the occurrence of D-amino acids in an L-amino-acid environment is known to disrupt secondary structure [5, 6]. The crucial role of the correct enantiochemistry can be also seen in the fact that organisms tightly control the use of D-amino acids. Some of the aminoacyl-tRNA synthetases, which implement the genetic code by loading a tRNA with the corresponding amino acid have proofreading capabilities to exclude D-amino acids [7]. Additionally, organisms have evolved specific enzymes, deacylases, which are able to recognize tRNAs acylated with a D-amino acid and cleave the bond between the amino acid and the tRNA, thus preventing a D-acylated tRNA from entering into protein synthesis [8].

Figure 1
figure 1

Selected stereochemical parameters of amino acids and nucleosides. (A) Chirality configuration at C α of an amino acid; (B) chiral centers at the ribose moiety of guanosine monophosphate; (C) cis and trans isomers of a peptide group. Carbon atoms are shown in green, nitrogen in blue, oxygen in red, and hydrogen in white. Chiral centers are surrounded by a transparent purple sphere.

Another type of asymmetry is encountered in the conformation of the peptide bond connecting the carboxy end of one amino acid to the amino end of the next one in a peptide or protein. Due to the partial double-bond character of the C n -Nn+1bond, the atoms Cα, n, C n , On, Cα, n+1, Nn+1and its hydrogen are in a plane (see, however, Ref [9]) and the rotation around the C n -Nn+1bond is restricted by a barrier of about 20 kcal/mol [10, 11]. Depending on the value of the dihedral angle ω described by Cα, n, C n , Nn+1and Cα, n+1, one can distinguish cis (ω ≈ 0°) and trans (ω ≈ 180°) isomers [12] (see Figure 1C). For sterical reasons, the trans isomer is energetically more stable and, thus, is the prevalent form in proteins. Additionally, the rather high rotational barrier makes the interconversion of the two isomers a very slow process at room temperature. Nevertheless, cis peptide bonds can be found in nature [13, 14]. The vast majority of cis bonds are observed before a proline residue, Xaa-Pro, with Xaa being any amino acid. The formation of these cis peptides is catalyzed by special enzymes, prolyl-cis/trans isomerases [15, 16]. Prolyl-cis/trans isomerization is an important molecular switch [17] and the occurrence of enzymes specialized for this particular isomerization underpins its biological significance. Non-prolyl cis peptide bonds can also be found in proteins, but much less frequently than Xaa-Pro [13, 14] and only one protein, DnaK, is known to promote peptide isomerization of non-prolyl peptide bonds [18]. In particular, DnaK was found to accelerate the isomerization of Ala-Xaa bonds, Xaa being Ala, Gly, Glu, Ile, Leu, Lys, Met, or Ser. Metal ions can also play a role both in stabilizing the cis isomer and in promoting isomerization [19–22].

Correct stereochemistry of a structural model is important for its interpretation and critical if the model is to be subject to a molecular dynamics simulation. Force fields typically employed in biomolecular simulations do not contain terms to enforce stereochemistry and support either enantiomer or peptide isomer. Errors in the input structure usually persist throughout the simulation and, as will be shown below, can propagate and lead to severe artifacts. Even in cases where stereochemical errors do not lead to such large-scale problems, such errors must be avoided since they represent deviations of the simulated system from the biological reality that is to be modeled. There is a steady trend in the field of biomolecular simulations toward the study of large biomolecular assemblies and the usage of models based on structure prediction. Our recent experience [23–25] shows that stereochemical errors often arise in the preparation of large systems for simulation, particularly when some components must be modeled prior to simulation.

A variety of servers and programs are available for structure validation, a vital stage in the preparation of files for deposition in a coordinate database. One example is the SAVES server [26], which provides an interface to tools such as PROCHECK [27], WHAT_CHECK [28], and other programs to detect irregularities in geometry and structure such as chirality, bond angles, close contacts, or rotamer states of amino acid side chains. Other examples with similar functionality include the MolProbity [29] server as well as the PDB validation service [30].

All of the aforementioned tools for structure validation are primarily designed to validate experimentally obtained models and not to ensure proper stereochemistry in simulations. Additionally, although irregularities are reported, none of the tools we are aware of allow the user to inspect and correct stereochemical errors easily and immediately. We have thus written software tools to help researchers easily detect, correct, and avoid stereochemical errors in simulations.

In the following, we start by illustrating how errors in chirality and peptide bond configuration can affect secondary structure of proteins in molecular dynamics simulations. We then present software tools to identify, visually inspect, and interactively correct stereochemical errors in structural models of proteins and nucleic acids. Next, we discuss the most common sources of stereochemical errors in simulations, alongside a systematic analysis of the entire Protein Data Bank [31]. Finally, we provide a recommended workflow to avoid stereochemical errors in biomolecular simulations.

Results and Discussion

Consequences of stereochemical errors in biomolecular simulations

In order to understand more directly the effects that errors in chirality or peptide bond configuration have on secondary structure, consider three simulations involving a 15-amino-acid-long α-helix AAQAAAAQAAAAQAA solvated in water. In the first simulation, the helix is stereochemically correct. For the second simulation, a chirality error is introduced at C α of Gln8 (Figure 2A). In the third case, the peptide bond between Gln8 and Ala9 is present as a cis isomer (Figure 2B).

Figure 2
figure 2

Stereochemical manipulations used for the illustrative simulations. Stereochemical manipulations used for the illustrative simulations. (A) Different configurations at the chiral C α at Gln8. (B) Different isomers at the Gln8-Ala9 peptide bond.

Figure 3 shows the conformations of the peptides after 32ns for each of the three scenarios, as well as the secondary structure content throughout each trajectory. By the end of the simulation the stereochemically correct helix is intact; in fact, the secondary structure is maintained during the entire simulation (Figure 3A). In contrast, a flip in chirality introduces a kink of almost 90° into the helix (Figure 3B). Interestingly, the two pieces separated by the kink remain helical. This conformational change can be understood considering that a D-amino acid in an environment with the backbone dihedrals φ/ψ of an L-amino acid helix features unfavourable steric interaction between the side chain and the peptide oxygen atom. The relaxation leading to the very prominent kink brings the side chain of Gln8 out of the eclipsed conformation relative to the peptide oxygen atom. In case of the cis peptide bond, the initial disturbance in the hydrogen bond network stabilizing the helix leads to a rearrangement of the network and to a complete loss of helicity downstream of Gln8 (Figure 3C). The configuration of a peptide bond is central to the types of secondary structure the peptide chain can assume: only the trans isomer accepts and donates hydrogen bonds in opposite directions allowing for formation of α-helices and β-sheets.

Figure 3
figure 3

Impact of stereochemical errors on the structure of an α -helix. Impact of stereochemical errors on the structure of an α-helix. The figure shows the starting conformations, the conformations after 32 ns of equilibrium simulation, as well as the secondary structure content throughout the simulation for (A) a stereochemically correct helix, (B) a helix with a D-Gln8, and (C) a helix with a cis peptide Gln8-Ala9. The axis on the left shows residue numbers and the red dashed line indicates the position of the introduced stereochemical errors. The colors in the bottom panels represent secondary structure content: helix (pink), turn (cyan) and coil (white).

As demonstrated above, the impact of errors in chirality or peptide bond isomerization on secondary structure can be dramatic. Note, however, that the chosen example represents the worst case scenario in terms of severity of the structural disturbance - in a real protein, tertiary interactions may provide stabilization of the native structure and, thus, dampen the effect or extend the time scale on which the structural disturbance becomes apparent.

Tools to identify, inspect, and correct stereochemical errors

Having demonstrated the impact of stereochemical errors on structure, the question of how to ensure stereochemical correctness in simulation naturally arises. In general, biological systems can contain amino acids or sugars of different chirality, as well as both peptide isomers. Thus, an automatic procedure to "correct" the structure is not appropriate unless one is absolutely certain that only one enantiomer and isomer occurs. Therefore, we designed a semi-automatic four-step protocol to correct errors and to ensure stereochemical integrity of a simulated system. For both chirality and peptide bond conformation, the protocol was implemented into easy-to-use plugins for the molecular visualization and analysis program VMD [32], referred to as Chirality and Cispeptide plugins, respectively. The plugins make use of the molecular dynamics simulation package NAMD [33] in the correction step. Both software packages are open source and freely available. The current implementation provides both a graphical and a command-line user interface. Use of each plugin follows a similar 4-step process, namely: (1) identify stereochemical anomalies; (2) visually inspect each anomaly and decide if it should be corrected; (3) move selected atoms as to change the stereochemical configuration; and (4) locally optimize the structure. In the following, we describe the four steps using the Cispeptide plugin as an example. Full documentation is available on the VMD website [34]. An additional step-by-step practical guide can be found in a tutorial describing both plugins [35].

In step (1), unusual stereochemical configurations (cis peptide bonds in case of the Cispeptide plugin) are identified. A peptide bond is considered to be an irregularity if the value of the dihedral angle ω' formed by O n , C n , Nn+1and Cα, n+1is larger than 85°. The slightly different, but equivalent, definition of the angle with respect to ω ensures unique atom names used for the angle measurement (exploited for computational efficiency) and the threshold is chosen so as to avoid false positives while maximizing the sensitivity when all the structures of the PDB were tested (see "Sources of stereochemical error below"). For chiral centers, the improper angle made up of non-hydrogen atoms with the chiral atom in the center is used to define non-standard configurations: a negative value of the considered improper indicates an unusual chirality. All identified irregularities are displayed in the corresponding panel (see "Identified cis peptide bonds" in Figure 4A). In step (2), the irregularities are inspected and a decision is made whether to keep or to modify each stereochemical configuration. Clicking on each stereochemical anomaly generates a corresponding molecular visualization (Figure 4B). As mentioned before, there are tools [36] that may help decide if the cis configuration is correct or if it should be switched to trans. In step (3), for each cis peptide bond to be manipulated, the user selects an atom (hydrogen or oxygen) to be moved (see Figure 4A). C α atoms are not allowed to move in order to minimize the structural impact. The movements are reflections of the positions of hydrogen or oxygen relative to the peptide nitrogen or carbon, respectively (see Figure 4C). Which of the two atoms should be moved highly depends on the environment, and is thus not automated. By visual inspection, the user can quite confidently determine which atom (hydrogen or oxygen) should be moved in order to optimize the hydrogen bonding network. Finally, in step (4) a local structure optimization is carried out using NAMD through the interactive molecular dynamics [37, 38] interface in VMD (Figure 4D). To enforce the target isomer, local restraints on the peptide dihedral are employed. The local nature of the manipulation and relaxation, which only involves residues within a given cutoff, prevents widespread influence of the procedure on the system. Finally, once the system is stereochemically correct, restraints can be constructed and used during further relaxation steps, prior to production simulations, to prevent changes in stereochemistry.

Figure 4
figure 4

The Cispeptide plugin. (A) The main window of the plugin showing the selection options in the top form, identified cis peptide bonds in the center form, and options to move atoms and to relax the structure in the bottom form. (B) Molecular visualization highlighting the selected bond generated by the plugin. (C) The configuration generated upon moving the oxygen atom along the C-O bond. (D) The control window for interactive molecular dynamics used in the correction of unusual configurations.

Sources of stereochemical errors

Most molecular dynamics simulations of single proteins are relatively safe from stereochemical errors, given that molecular structures are validated upon deposition into the PDB. The irregularities detected in the validation step are usually either corrected or reported in the header of the PDB file. Using a structure from the PDB for simulations, together with a careful examination of the file header, usually means that the stereochemical integrity of the simulation is secure (barring special cases discussed below), because during an equilibrium simulation the force field will preserve stereochemistry. However, since force fields used for molecular dynamics simulations support both types of chirality as well as both peptide bond isomers, an existing error will persist throughout the simulation.

To quantify stereochemical anomalies within the Protein Data Bank, all the structures contained in the PDB were analyzed using the Chirality and Cispeptide VMD plugins introduced above. As shown in Table 1 around 4,000 unusual configurations with respect to chirality could be found in protein and nucleic acid structures, as well as more than 100,000 cis peptide bonds. Analysis of the complete dataset shows that one chirality error is expected for every 10,000 residues, whereas for every few hundred residues a cis peptide bond can be expected. The likelihood of stereochemical anomalies varies with the structure determination method, in particular when chirality is concerned. For example, for structures derived from electron microscopy, one chirality error can be expected for every 500 residues (see Table 1). The higher error rate seen in electron microscopy-derived structures are likely due to the much lower resolution on average, combined with use of a wide range of modeling approaches, including manual manipulation of structures.

Table 1 Chirality errors and cis peptide bonds in the Protein Data Bank

In a subset of structures, chirality errors (153 structures) and cis peptide bonds (62 structures) were reported in the PDB file, but were not detected by the plugins. To understand the discrepancies, about half of the structures in each subset was visually inspected. In the case of cis peptides, the differences were due to missing atoms in the deposited model, errors being present in only one alternative conformation or model (see Methods), pathologically distorted structures, and, most frequently, incorrect entries in the PDB header. Similarly, the reasons for the discrepancies in chirality were that errors occurred in cofactors not checked by the plugin, the anomalies were present in only one alternative conformation or model (see Methods), or the annotations in the PDB header were incorrect.

The plugins also identified many structures with stereochemical anomalies not reported in the PDB header. According to the current PDB format standard, each cis peptide bond should be reported in a separate CISPEP record. It is then possible to compare each identified cis peptide bond with the PDB header. There were in total 2,518 cis peptide bonds identified by the Cispeptide plugin but not reported in the PDB files. Visual inspection of 20 such files did not reveal any incorrect identification by the plugin. Chirality issues should be reported in CAVEAT records of the PDB header. Unfortunately, since CAVEAT records are free format, rarely are chirality errors individually reported. Furthermore, many PDB structures have chiral errors inconsistently documented in REMARK 500 records instead. Thus, a similar comparison on a per-error basis cannot be performed for the Chirality plugin; instead, we can only report that 348 PDB files contain chirality anomalies according to the Chirality plugin, but lack the corresponding annotation in their PDB headers. Upon visual inspection of 20 files, the structures fell into two categories: either there should be a D-amino acid (non-ribosomal peptides) or the PDB header is simply missing the required annotation.

Apart from errors already present in experimentally determined structures, our experience shows that there are three main sources of stereochemical errors in simulations. The first source is found in modeling steps, particularly homology modeling: any regions of a structure that were modeled de novo, especially at the junctions of the known and modeled part of a protein, are prone to peptide isomerization errors. The second common source is found in the setup protocol for a simulation. This includes the preparation of the system and its initial relaxation. In particular, sterical clashes between atoms can lead to errors in chirality and isomerization state of peptide bonds. Again, such behavior is particularly likely at the interface between known and modeled portions of a structure, which may contain severe distortions prior to equilibration. Although the barriers for isomerization or a flip in chirality are large enough to prevent these events in an equilibrium simulation at physiological temperatures (e.g., 21 kcal/mol in CHARMM22 [39]), forces arising during the initial structure optimization, necessary to relax possible clashes, may be large enough to introduce errors into the structure. This source of errors becomes increasingly important as the field of computational biophysics moves towards multi-component assemblies, the structures of which are often modeled based on high-resolution models of their constituents and low-resolution data of the whole complex. Finally, errors can also be introduced during structure optimizations if additional forces are applied on the system. Of particular interest at this point is the molecular dynamics flexible fitting (MDFF) method [40], which flexibly fits atomic-resolution structures into low-resolution density maps. In some (rare) cases, it was observed that forces arising from the MDFF method during initial structure optimization were large enough to cause stereochemical errors. System setup protocols where the system is simulated at very high temperatures (e.g., to obtain heat-denatured structures) may also allow incorrect isomerization events.

The Chirality and Cispeptide plugins can be used to generate harmonic restraints designed to preserve the current isomerization state of each chiral center and peptide bond. The restraints can be used in simulations with NAMD, effectively preventing stereochemical errors from arising during simulation. These restraints should be removed prior to production equilibrium simulations, since they are not required and would represent an unnecessary modification of the force field employed.

Recommended workflow

In order to avoid stereochemical errors, the following simple workflow is recommended for standard MD simulations:

  1. 1.

    Build system for MD simulation (model missing components, assemble structure, embed system in a water box, and add counterions).

  2. 2.

    Energy-minimize structure.

  3. 3.

    Check stereochemistry with the Chirality and Cispeptide plugins, correcting errors if applicable. Repeat until no further errors are detected. Make sure that the detected irregularities are indeed errors and not naturally occurring.

  4. 4.

    Proceed to production simulation.

For simulations in which large forces are expected (e.g., flexible fitting with the MDFF method [40] or temperature-induced denaturation), it is recommended that, in addition to the workflow above, harmonic restraints generated by the Chirality and Cispeptide plugins are applied throughout the simulation.

Conclusions

The simulations presented here illustrate the drastic effects that stereochemical errors can have in biomolecular simulations. Experimentally determined structures may contain stereochemical errors, and various modeling approaches can further increase the number of such errors. As the community moves toward simulation of large, multi-component complexes and uses to an increasing extent models based on structure prediction, the issue of stereochemical correctness becomes even more relevant. We thus developed tools to identify, inspect, and correct stereochemical errors in protein and nucleic acid structures. In particular, chirality and the isomerization state of a peptide bond are examined. The main advantage of the offered tools is the possibility to immediately inspect and correct the detected errors. The tools are implemented as plugins to the molecular visualization and analysis program VMD. The recommended workflow presented above effectively avoids artifacts in simulations due to stereochemical errors. We hope that checks for stereochemical correctness become a standard step of any biomolecular simulation or generation of predicted structural models for proteins and nucleic acids.

Availability and Requirements

Methods

Molecular dynamics simulations

Molecular dynamics simulations were conducted using NAMD 2.7 [33]. The system consisted of the 15-amino-acid-long α-helix AAQAAAAQAAAAQAA solvated in TIP3P water. The N- and C-terminus were acetylated and amidated, respectively. The system was set up in VMD [32]. In particular, the helix was constructed using the molefacture plugin, after which the full system was built with the solvate and psfgen plugins.

All simulations were performed in the NpT ensemble (T = 310 K, p = 1 atm). The stereochemically correct system was equilibrated using a 2-step protocol. First, water and side chains were equilibrated for 400 ps while the backbone of the peptide was restrained, after which all restraints were removed and the system was equilibrated for an additional 3.6 ns. The resulting structure was used to manually introduce a chirality flip at Gln8 and to isomerize the peptide bond between Gln8 and Ala9 into the cis form (see Figure 2). Starting from the three obtained systems, the simulations were continued for 32 ns in each case. The equations of motion were integrated using periodic boundary conditions and a 2-fs time step, with bonded interactions calculated every 2 fs. Nonbonded, short-range interactions were calculated every 4 fs using a distance cut-off of 10 Ã… with a switching function applied at 9 Ã…. Long-range electrostatic interactions were updated every 6 fs using the Particle Mesh Ewald (PME) method and the PME grid density was never less than 1/Ã…3. The simulations presented here were performed with the CHARMM27 force field [39, 41] with the CMAP correction [42].

Analysis of the PDB

Each structure available in the PDB as of mid-September 2010 was analyzed using the Chirality and Cispeptide VMD plugins. In case of multiple models (typically in structures solved by NMR), only the first model was considered in the analysis. When multiple alternative conformations of certain residues were present, as indicated by the altLoc field in the PDB file, only the first conformation was considered. For comparison and validation purposes, each PDB file header was parsed for reported cis peptide bonds (CISPEP records) and unusual chirality configurations (CAVEAT and REMARK 500 records).

References

  1. Höltje JV: Growth of the Stress-Bearing and Shape-Maintaining Murein Sacculus of Escherichia coli. Microbiol Mol Biol Rev 1998, 62: 181–203.

    PubMed Central  PubMed  Google Scholar 

  2. Vollmer W, Blanot D, de Pedro MA: Peptidoglycan structure and architecture. FEMS Microbiol Rev 2008, 32: 149–167. 10.1111/j.1574-6976.2007.00094.x

    Article  CAS  PubMed  Google Scholar 

  3. Lam H, Oh DC, Cava F, Takacs CN, Clardy J, de Pedr MA, Waldor MK: D-Amino Acids Govern Stationary Phase Cell Wall Remodeling in Bacteria. Science 2009.

    Google Scholar 

  4. Annedi SC, Biabani F, Poduch E, Mannargudi BM, Majumder K, Wei L, Khayat R, Tong L, Kotra LP: Engineering D-amino acid containing novel protease inhibitors using catalytic site architecture. Bioorg Med Chem 2006, 14: 214–236. 10.1016/j.bmc.2005.08.019

    Article  CAS  PubMed  Google Scholar 

  5. Krause E, Beyermann M, Dathe M, Rothemund S, Bienert M: Location of an amphipathic alpha-helix in peptides using reversed-phase HPLC retention behavior of D-amino acid analogs. Anal Chem 1995, 67: 252–258. 10.1021/ac00098a003

    Article  CAS  PubMed  Google Scholar 

  6. Lee DL, Powers JP, Pflegerl K, Vasil ML, Hancock RE, Hodges RS: Effects of single D-amino acid substitutions on disruption of beta-sheet structure and hydrophobicity in cyclic 14-residue antimicrobial peptide analogs related to gramicidin S. J Pept Res 2004, 63: 69–84. 10.1046/j.1399-3011.2003.00106.x

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Hussain T, Kruparani S, Pal B, Dock-Bregeon AC, Dwivedi S, Shekar M, Sureshbabu K, Sankaranarayanan R: Post-transfer editing mechanism of a D-aminoacyl-tRNA deacylase-like domain in threonyl-tRNA synthetase from archaea. EMBO J 2006, 25: 4152–4162. 10.1038/sj.emboj.7601278

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Soutourina J, Plateau P, Blanquet S: Metabolism of D-Aminoacyl-tRNAs in Escherichia coli and Saccharomyces cerevisiae Cells. J Biol Chem 2000, 275: 32535–32542.

    Article  CAS  PubMed  Google Scholar 

  9. MacArthur MW, Thornton JM: Deviations from Planarity of the Peptide Bond in Peptides and Proteins. J Mol Biol 1996, 264: 1180–1195. 10.1006/jmbi.1996.0705

    Article  CAS  PubMed  Google Scholar 

  10. Scherer G, Kramer ML, Schutkowski M, Reimer U: Barriers to rotation in secondary amide peptide bonds. J Am Chem Soc 1998, 120: 5568–5574. 10.1021/ja980181t

    Article  CAS  Google Scholar 

  11. Schiene-Fischer C, Fischer G: Direct Measurement Indicates a Slow Cis/Trans Isomerization at the Secondary Amide Peptide Bond of Glycylglycine. J Am Chem Soc 2001, 123: 6227–6231. 10.1021/ja0042480

    Article  CAS  PubMed  Google Scholar 

  12. on Biochemical Nomenclature IIC: Abbreviations and symbols for the description of the conformation of polypeptide chains. J Mol Biol 1970, 52: 1–17.

    Article  Google Scholar 

  13. Pal D, Chakrabarti P: Cis peptide bondsnext term in proteins: residues involved, their conformations, interactions and locations1. J Mol Biol 1999, 294: 271–288. 10.1006/jmbi.1999.3217

    Article  CAS  PubMed  Google Scholar 

  14. Jabs A, Weiss MS, Hilgenfeld R: Non-proline cis peptide bonds in proteins. J Mol Biol 1999, 286: 291–304. 10.1006/jmbi.1998.2459

    Article  CAS  PubMed  Google Scholar 

  15. Fischer G, Bang H, Mech C: Detection of enzyme catalysis for cis-trans isomerization of peptide binding using proline-containing peptides as substrates. Biochim Biophys Acta 1984, 43: 1101–1111.

    CAS  Google Scholar 

  16. Schmid FX, Mayr LM, Mücke M, Schönbrunner E: Prolyl isomerases: role in protein folding. Adv Prot Chem 1993, 44: 25–66.

    Article  CAS  Google Scholar 

  17. Eckert B, Martin A, Balbach J, Schmid FX: Prolyl isomerization as a molecular timer in phage infection. Nat Struct Mol Biol 2005, 12: 619–623. 10.1038/nsmb946

    Article  CAS  PubMed  Google Scholar 

  18. Schiene-Fischer C, Habazettl J, Schmid FX, Fischer G: The Hsp70 chaperone DnaK is a secondary amide peptide bond cis-trans isomerase. Nat Struct Biol 2002, 9: 419–424. 10.1038/nsb804

    Article  CAS  PubMed  Google Scholar 

  19. Hohenester E, Maurer P, Hohenadl C, Timpl R, Jansonius JN, Engel J: Structure of a novel extracellular Ca2+-binding module in BM-40. Nat Struct Biol 1996, 3: 67–73. 10.1038/nsb0196-67

    Article  CAS  PubMed  Google Scholar 

  20. Bouckaert J, Dewallef Y, Poortmans F, Wyns L, Loris R: The structural features of concanavalin A governing non-proline peptide isomerization. J Biol Chem 2000, 275: 19778–19787. 10.1074/jbc.M001251200

    Article  CAS  PubMed  Google Scholar 

  21. Wilson JJ, Matsushita O, Okab A, Sakon J: A bacterial collagen-binding domain with novel calcium-binding motif controls domain orientation. EMBO J 2003, 22: 1743–1752. 10.1093/emboj/cdg172

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Spiriti‡ J, van der Vaart A: Mechanism of the Calcium-Induced trans-cis Isomerization of a Non-Prolyl Peptide Bond in Clostridium histolyticum Collagenase. Biochemistry 2010, 49: 5314–5320. 10.1021/bi1007047

    Article  Google Scholar 

  23. Villa E, Sengupta J, Trabuco LG, LeBarron J, Baxter WT, Shaikh TR, Grassucci RA, Nissen P, Ehrenberg M, Schulten K, Frank J: Ribosome-induced Changes in Elongation Factor Tu Conformation Control GTP Hydrolysis. Proc Natl Acad Sci USA 2009, 106: 1063–1068. 10.1073/pnas.0811370106

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Gumbart J, Trabuco LG, Schreiner E, Villa E, Schulten K: Regulation of the protein-conducting channel by a bound ribosome. Structure 2009, 17: 1453–1464. 10.1016/j.str.2009.09.010

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Trabuco LG, Schreiner E, Eargle J, Cornish P, Ha T, Luthey-Schulten Z, Schulten K: The role of L1 stalk-tRNA interaction in the ribosome elongation cycle. J Mol Biol 2010, 402: 741–760. 10.1016/j.jmb.2010.07.056

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. SAVES[Http://nihserver.mbi.ucla.edu/SAVES/]

  27. Laskowski RA, MacArthur MW, Moss DS, Thornton JM: PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Cryst 1993, 26: 283–291. 10.1107/S0021889892009944

    Article  CAS  Google Scholar 

  28. Hooft RWW, Vriend G, Sander C, Abola EE: Errors in protein structures. Nature 1996, 381: 272–272.

    Article  CAS  PubMed  Google Scholar 

  29. Chen VB, Arendall WB III, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, Murray LW, Richardson JS, Richardson DC: MolProbity: all-atom structure validation for macromolecular crystallography. Acta Cryst D 2010, 66: 12–21. 10.1107/S0907444909042073

    Article  CAS  Google Scholar 

  30. PDB Validation Server[Http://deposit.pdb.org/validate/]

  31. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res 2000, 28: 235–242. 10.1093/nar/28.1.235

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Humphrey W, Dalke A, Schulten K: VMD - Visual Molecular Dynamics. J Mol Graphics 1996, 14: 33–38. 10.1016/0263-7855(96)00018-5

    Article  CAS  Google Scholar 

  33. Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel RD, Kale L, Schulten K: Scalable Molecular Dynamics with NAMD. J Comp Chem 2005, 26: 1781–1802. 10.1002/jcc.20289

    Article  CAS  Google Scholar 

  34. http://www.ks.uiuc.edu/Research/vmd/plugins/cispeptide; http://www.ks.uiuc.edu/Research/vmd/plugins/chirality

  35. [http://www.ks.uiuc.edu/Training/Tutorials/#structurecheck].

  36. Exarchos KP, Exarchos TP, Papaloukas C, Troganis AN, Fotiadis DI: PBOND: Web Server for the Prediction of Proline and Non-Proline cis/trans Isomerization. Gen Prot Bioinf 2009, 7: 138–142.

    Article  CAS  Google Scholar 

  37. Stone J, Gullingsrud J, Grayson P, Schulten K: A System for Interactive Molecular Dynamics Simulation. In 2001 ACM Symposium on Interactive 3D Graphics. Edited by: Hughes JF, Séquin CH. New York: ACM SIGGRAPH; 2001:191–194.

    Chapter  Google Scholar 

  38. Grayson P, Tajkhorshid E, Schulten K: Mechanisms of Selectivity in Channels and Enzymes Studied with Interactive Molecular Dynamics. Biophys J 2003, 85: 36–48. 10.1016/S0006-3495(03)74452-X

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. MacKerell AD Jr, Bashford D, Bellott M, Dunbrack RL Jr, Evanseck JD, Field MJ, Fischer S, Gao J, Guo H, Ha S, Joseph D, Kuchnir L, Kuczera K, Lau FTK, Mattos C, Michnick S, Ngo T, Nguyen DT, Prodhom B, Reiher IWE, Roux B, Schlenkrich M, Smith J, Stote R, Straub J, Watanabe M, Wiorkiewicz-Kuczera J, Yin D, Karplus M: All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B 1998, 102: 3586–3616. 10.1021/jp973084f

    Article  CAS  PubMed  Google Scholar 

  40. Trabuco LG, Villa E, Mitra K, Frank J, Schulten K: Flexible Fitting of Atomic Structures into Electron Microscopy Maps Using Molecular Dynamics. Structure 2008, 16: 673–683. 10.1016/j.str.2008.03.005

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Foloppe N, MacKerell AD Jr: All-atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data. J Comp Chem 2000, 21: 86–104. 10.1002/(SICI)1096-987X(20000130)21:2<86::AID-JCC2>3.0.CO;2-G

    Article  CAS  Google Scholar 

  42. MacKerell AD Jr, Feig M, Brooks CL III: Extending the Treatment of Backbone Energetics in Protein Force Fields: Limitations of Gas-Phase Quantum Mechanics in Reproducing Protein Conformational Distributions in Molecular Dynamics Simulations. J Comp Chem 2004, 25: 1400–1415. 10.1002/jcc.20065

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Institutes of Health (P41-RR005969) and the National Science Foundation (PHY0822613). E.S. and L.G.T. were supported by fellowships from the Humboldt Foundation and the European Molecular Biology Organization, respectively.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Klaus Schulten.

Additional information

Authors' contributions

ES and LGT conceived the ideas behind the paper. LGT wrote and tested the software. ES assisted in writing and testing of the software. PLF assisted in testing the software. ES performed and analyzed molecular dynamics simulations. LGT analyzed structures of the Protein Data Bank. All authors were involved in writing of the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Schreiner, E., Trabuco, L.G., Freddolino, P.L. et al. Stereochemical errors and their implications for molecular dynamics simulations. BMC Bioinformatics 12, 190 (2011). https://doi.org/10.1186/1471-2105-12-190

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-12-190

Keywords