ABSTRACT
Molecular dynamics (MD) simulations have become increasingly popular in studying the motions and functions of biomolecules. The accuracy of the simulation, however, is highly determined by the molecular mechanics (MM) force field (FF), a set of functions with adjustable parameters to compute the potential energies from atomic positions. However, the overall quality of the FF, such as our previously published ff99SB and ff14SB, can be limited by assumptions that were made years ago. In the updated model presented here (ff19SB), we have significantly improved the backbone profiles for all 20 amino acids. We fit coupled φ/ψ parameters using 2D φ/ψ conformational scans for multiple amino acids, using as reference data the entire 2D quantum mechanics (QM) energy surface. We address the polarization inconsistency during dihedral parameter fitting by using both QM and MM in aqueous solution. Finally, we examine possible dependency of the backbone fitting on side chain rotamer. To extensively validate ff19SB parameters, and to compare to results using other Amber models, we have performed a total of â¼5 ms MD simulations in explicit solvent. Our results show that after amino-acid-specific training against QM data with solvent polarization, ff19SB not only reproduces the differences in amino-acid-specific Protein Data Bank (PDB) Ramachandran maps better but also shows significantly improved capability to differentiate amino-acid-dependent properties such as helical propensities. We also conclude that an inherent underestimation of helicity is present in ff14SB, which is (inexactly) compensated for by an increase in helical content driven by the TIP3P bias toward overly compact structures. In summary, ff19SB, when combined with a more accurate water model such as OPC, should have better predictive power for modeling sequence-specific behavior, protein mutations, and also rational protein design. Of the explicit water models tested here, we recommend use of OPC with ff19SB.
Subject(s)
Amino Acids/chemistry , Peptides/chemistry , Proteins/chemistry , Water/chemistry , Molecular Dynamics Simulation , Protein Conformation , Protein Stability , Quantum Theory , ThermodynamicsABSTRACT
Carbon-donated hydrogen bonds (CDHBs) are weak forms of hydrogen bonding (0.5-1.0â kcal mol(-1) ) that are difficult to detect, and thus their roles in the structure and functionality of chemical systems often go unrecognized. Utilizing a computational approach, the existence of a structurally significant CDHB in the medically relevant protein Streptococcus pneumoniae hyaluronate lyase (SpnHL) is affirmed. The structure of a tetrapeptide fragment model containing the CDHB was optimized with second-order perturbation theory. From this, a CDHB with bond distance and angle consistent with previously discovered CDHBs and comparable to neighboring traditional HBs in the fragment model was found. The CDHB competes with another donor T253 OH, whereby the two alternate in strength between protein conformations, imbuing αHelixâ 3 appreciable flexibility. The CDHB seems to exist in spite of torsional and steric strain on the donor methyl group. It is postulated that the CDHB could aid in either counteracting the macrodipole of αHelixâ 3 or protecting the A249 CO from destabilizing interactions with the adjacent solvent. Employing the energy gradients from the optimization, the torque generated by the fragment model was computed, which accurately predicts the direction of rotation of αHelixâ 3 observed from experiment. A strongly correlated motion between αHelixâ 3 and αHelicesâ 2, 4, and 5 was noted, which the interactions of the fragment model help drive by generating a torque much larger than necessary to rotate just αHelixâ 3. Considering these results, we conclude that CDHBs should be considered as possible beneficial components of chemical and biological phenomena.