## Abstract Computer simulations offer a valuable way to study membrane systems, from simple lipid bilayers to large transmembrane protein complexes and lipidβnucleic acid complexes for drug delivery. Their accuracy depends on the quality of the force field parameters used to describe the component
Parametrization of aliphatic CHn united atoms of GROMOS96 force field
β Scribed by Daura, Xavier; Mark, Alan E.; Van Gunsteren, Wilfred F.
- Publisher
- John Wiley and Sons
- Year
- 1998
- Tongue
- English
- Weight
- 203 KB
- Volume
- 19
- Category
- Article
- ISSN
- 0192-8651
No coin nor oath required. For personal study only.
β¦ Synopsis
The derivation of the van der Waals parameters for the aliphatic CH united atoms of the GROMOS96 force field is presented. The parameters n have been adjusted to reproduce the experimental enthalpies of vaporization and vapor pressures or densities of a set of nine alkanes in the liquid state at 298 Ε½ . K or at the boiling point in the case of methane , using a cutoff radius for the van der Waals interactions of 1.6 nm. Force fields to be used in molecular simulations are bound to the conditions chosen for their parametrization, for example, the temperature, the densities of the systems included in the calibration set, or the cutoff radius used for the nonbonded interactions. Van der Waals parameters for the CH united atoms of earlier GROMOS force fields were n developed using a cutoff radius of 0.8 nm for the van der Waals interactions. Because the van der Waals interaction energy between aliphatic groups separated by distances between 0.8 and 1.4 nm is not negligible at liquid densities, the use of these parameters in combination with longer cutoffs leads to an overestimation of the attractive van der Waals interaction energy. The relevance of this excess attraction depends on the size of the groups that are interacting, as well as on their local densities. Free energies of hydration have been calculated for five alkanes.
π SIMILAR VOLUMES
Biomolecular force fields for use in molecular dynamics (MD) simulations of proteins, DNA, or membranes are generally parametrized against ab initio quantum-chemical and experimental data for small molecules. The application of a force field in a simulation of a biomolecular system, such as a protei
We refined the united atom field for the simulations of phospholipid membranes. To validate this potential we performed 1000-ps Ε½ . constant pressure simulation of a dipalmitoylphosphatidicholine DPPC bilayer Λ2 Ε½ . at T s 50Π C. The average area per head group 61.6 " 0.6 A obtained in our Λ2 Ε½ . si
A molecular mechanics force field was developed for systems bearing the N -C -0 unit on the basis of 6-31G\* and 4-21G %b initio" calculations with full optimization of the geometry and experimental heats of formation. The parameters used, which implicitly included the anomeric effect, provided good