24 research outputs found
Secondary-Structure Design of Proteins by a Backbone Torsion Energy
We propose a new backbone-torsion-energy term in the force field for protein
systems. This torsion-energy term is represented by a double Fourier series in
two variables, the backbone dihedral angles phi and psi. 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 phi and psi 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/decrease the alpha-helix-forming-tendencies by
lowering/raising the torsion-energy surface in the alpha-helix region and
likewise increase/decrease the beta-sheet-forming tendencies by
lowering/raising the surface in the beta-sheet region in the Ramachandran
space. We applied our approach to AMBER parm94 and AMBER parm96 force fields
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 alpha-helical and beta-hairpin peptides.Comment: 13 pages, (Revtex4), 5 figure
Energetics of Base Pairs in B-DNA in Solution: An Appraisal of Potential Functions and Dielectric Treatments
Generalization of the Gaussian electrostatic model: Extension to arbitrary angular momentum, distributed multipoles, and speedup with reciprocal space methods
GRIFFIN: A Versatile Methodology for Optimization of Protein−Lipid Interfaces for Membrane Protein Simulations
Free Energy Barriers for the N-Terminal Asparagine to Succinimide Conversion: Quantum Molecular Dynamics Simulations for the Fully Solvated Model
The onset of Taylor-Görtler vortices in the time-dependent Couette flow induced by an impulsively imposed shear stress
Recommended from our members
Molecular dynamics simulations of proteins: can the explicit water model be varied?
In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties