GROMOS Force Field (GROMO + force_field)

Distribution by Scientific Domains


Selected Abstracts


A new GROMOS force field for hexopyranose-based carbohydrates

JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 13 2005
Roberto D. Lins
Abstract A new parameter set (referred to as 45A4) is developed for the explicit-solvent simulation of hexopyranose-based carbohydrates. This set is compatible with the most recent version of the GROMOS force field for proteins, nucleic acids, and lipids, and the SPC water model. The parametrization procedure relies on: (1) reassigning the atomic partial charges based on a fit to the quantum-mechanical electrostatic potential around a trisaccharide; (2) refining the torsional potential parameters associated with the rotations of the hydroxymethyl, hydroxyl, and anomeric alkoxy groups by fitting to corresponding quantum-mechanical profiles for hexopyranosides; (3) adapting the torsional potential parameters determining the ring conformation so as to stabilize the (experimentally predominant) 4C1 chair conformation. The other (van der Waals and nontorsional covalent) parameters and the rules for third and excluded neighbors are taken directly from the most recent version of the GROMOS force field (except for one additional exclusion). The new set is general enough to define parameters for any (unbranched) hexopyranose-based mono-, di-, oligo- or polysaccharide. In the present article, this force field is validated for a limited set of monosaccharides (,- and ,-D-glucose, ,- and ,-D-galactose) and disaccharides (trehalose, maltose, and cellobiose) in solution, by comparing the results of simulations to available experimental data. More extensive validation will be the scope of a forthcoming article. 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1400,1412, 2005 [source]


An improved nucleic acid parameter set for the GROMOS force field

JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 7 2005
Thereza A. Soares
Abstract Over the past decades, the GROMOS force field for biomolecular simulation has primarily been developed for performing molecular dynamics (MD) simulations of polypeptides and, to a lesser extent, sugars. When applied to DNA, the 43A1 and 45A3 parameter sets of the years 1996 and 2001 produced rather flexible double-helical structures, in which the Watson,Crick hydrogen-bonding content was more limited than expected. To improve on the currently available parameter sets, the nucleotide backbone torsional-angle parameters and the charge distribution of the nucleotide bases are reconsidered based on quantum-chemical data. The new 45A4 parameter set resulting from this refinement appears to perform well in terms of reproducing solution NMR data and canonical hydrogen bonding. The deviation between simulated and experimental observables is now of the same order of magnitude as the uncertainty in the experimental values themselves. 2005 Wiley Periodicals, Inc. J Comput Chem 26: 725,737, 2005 [source]


Temperature and urea induced denaturation of the TRP-cage mini protein TC5b: A simulation study consistent with experimental observations

PROTEIN SCIENCE, Issue 10 2009
Z. Gattin
Abstract The effects of temperature and urea denaturation (6M urea) on the dominant structures of the 20-residue Trp-cage mini-protein TC5b are investigated by molecular dynamics simulations of the protein at different temperatures in aqueous and in 6M urea solution using explicit solvent degrees of freedom and the GROMOS force-field parameter set 45A3. In aqueous solution at 278 K, TC5b is stable throughout the 20 ns of MD simulation and the trajectory structures largely agree with the NMR-NOE atom,atom distance data available. Raising the temperature to 360 K and to 400 K, the protein denatures within 22 ns and 3 ns, showing that the denaturation temperature is well below 360 K using the GROMOS force field. This is 40,90 K lower than the denaturation temperatures observed in simulations using other much used protein force fields. As the experimental denaturation temperature is about 315 K, the GROMOS force field appears not to overstabilize TC5b, as other force fields and the use of continuum solvation models seem to do. This feature may directly stem from the GROMOS force-field parameter calibration protocol, which primarily involves reproduction of condensed phase thermodynamic quantities such as energies, densities, and solvation free energies of small compounds representative for protein fragments. By adding 6M urea to the solution, the onset of denaturation is observed in the simulation, but is too slow to observe a particular side-chain side-chain contact (Trp6-Ile4) that was experimentally observed to be characteristic for the denatured state. Interestingly, using temperature denaturation, the process is accelerated and the experimental data are reproduced. [source]