Search

We have previously proposed a method for refining force-field parameters of protein systems, which consists of minimising the summation of the square of the force acting on each atom in the proteins with the structures from the protein data bank (PDB). The results showed that the modified force-field parameters for all-atom model gave structures more consistent with the experimental implications than the original force fields. In this work, we applied this method and a new method to the OPLS-UA force field. In the new method, we perform a minimisation of the average of the root-mean-square deviation of various protein structures from the native structure. We selected some torsion-energy parameters for this optimisation, and 100 molecules from the PDB were used. The results imply that the new force-field parameters gave structures of two peptides more consistent with the experimental implications for the secondary structure-forming tendencies than the original OPLS-UA force field. [ABSTRACT FROM AUTHOR]

We performed folding simulations of three proteins using four force fields, AMBER parm96, AMBER parm99, CHARMM 27 and OPLS-AA/L, in order to examine the features of these force fields. We studied three proteins, protein A (all α-helix), cold-shock protein (all β-strand) and protein G (α/β-structures), for the folding simulations. For the simulation, we used the simulated annealing molecular dynamics method, which was performed 50 times for each protein using the four force fields. The results showed that the secondary-structure-forming tendencies are largely different among the four force fields. AMBER parm96 favours β-bridge structures and extended β-strand structures, and AMBER parm99 favours α-helix structures and 310-helix structures. CHARMM 27 slightly favours α-helix structures, and there are also π-helix and β-bridge structures. OPLS-AA/L favours α-helix structures and 310-helix structures. [ABSTRACT FROM AUTHOR]

We examined a new backbone torsion-energy term proposed by us in the force field for protein systems. This torsion-energy term is represented by a double Fourier series in two variables, namely the backbone dihedral angles φ and ψ. It gives a natural representation of the torsion energy in the Ramachandran space in the sense that any two-dimensional energy surface periodic in both φ and ψ can be expanded by the double Fourier series. We can then easily control secondary-structure-forming tendencies by modifying the torsion-energy surface. For instance, we can increase or decrease the α-helix-forming-tendencies by lowering or raising the torsion-energy surface in the α-helix region and likewise increase or decrease the β-sheet-forming tendencies by lowering or raising the surface in the β-sheet region in the Ramachandran space. We applied this torsion-energy modification method to six force fields, AMBER parm94, AMBER parm96, AMBER parm99, CHARMM27, OPLS-AA and OPLS-AA/L, and demonstrated that our modifications of the torsion-energy terms resulted in the expected changes of secondary-structure-forming tendencies by performing folding simulations of α-helical and β-hairpin peptides. [ABSTRACT FROM AUTHOR]

We propose a new method of optimisation of backbone torsion-energy parameters in the force field for molecular simulations of protein systems. This method is based on the idea of balancing the secondary-structure-forming tendencies, namely, those of α-helix and β-sheet structures. We perform a minimisation of the backbone dihedral angle-based root-mean-square deviation of the helix and β structure regions in many protein structures. As an example, we optimised the backbone torsion-energy parameters of AMBER parm96 force field using 100 protein molecules from the Protein Data Bank. We then performed folding simulations of α-helical and β-hairpin peptides, using the optimised force field. The results imply that the new force-field parameters give structures more consistent with the experimental implications than the original AMBER parm96 force field. [ABSTRACT FROM AUTHOR]