Selective Enhanced Sampling in Dihedral Energy Facilitating to Overcome the Dihedral Energy Increase in Protein Folding and Accelerating the Searching for Protein Native Structure

Qiang Shao,Lijiang Yang,Weiliang Zhu
DOI: https://doi.org/10.1039/C9CP00615J
IF: 3.3
2019-01-01
Physical Chemistry Chemical Physics
Abstract:The dihedral energy function is the most primarily involved in the molecular mechanics (MM) force field parameter optimization. A selective enhanced sampling on dihedral energy could effectively reflect the influence of dihedral energy settings on protein secondary structure representation, which in turn testifies the availability of the force field in folding simulation. Here, a Dihedral based Selective Integrated-Tempering-Sampling Molecular Dynamics (D-SITSMD) simulation method is shown to provide a selective enhanced sampling on dihedral energy without introducing large energetic noise. Its capabilities of searching for protein natively folded structure and providing underlying folding pathway are evaluated through the folding test of three peptides (chignolin, TC5b, and HP35) with multiple AMBER force fields (FF14SBonlysc, FF99SBildn, or FF03) and the comparison to presenting experimental data and REMD simulations. Both abovementioned capabilities are improved, displaying the potential of D-SITSMD in the studies of in silico protein folding and structure refinement. Additionally, it is commonly observed among the test simulation systems that their folding processes are thermodynamically favorable for non-bonded vdW and electrostatic energies but unfavorable for dihedral energy, giving the folding barrier height correlated with the dihedral energy increase from the unfolded to folded state whereas the unfolding free energy barrier correlated with the combined increase of vdW and electrostatic energies in the unfolding process. It is speculated that the influence of a force field on the folding barrier of a protein is fulfilled mainly through regulating the contribution of dihedral energy to determine the secondary structure formation in the global folding process.
What problem does this paper attempt to address?