Search

While the construction of a dependable force field for performing classical molecular dynamics (MD) simulation is crucial for elucidating the structure and function of biomolecular systems, the attempts to do this for glycans are relatively sparse compared to those for proteins and nucleic acids. Currently, the use of GLYCAM06 force field is the most popular, but there have been a number of concerns about its accuracy in the systematic description of structural changes. In the present work, we focus on the improvement of the GLYCAM06 force field for β-d-glucose, a simple and the most abundant monosaccharide molecule, with the aid of machine learning techniques implemented with the TensorFlow library. Following the pre-sampling over a wide range of configuration space generated by MD simulation, the atomic charge and dihedral angle parameters in the GLYCAM06 force field were re-optimized to accurately reproduce the relative energies of β-d-glucose obtained by the density functional theory (DFT) calculations according to the structural changes. The validation for the newly proposed force-field parameters was then carried out by verifying that the relative energy errors compared to the DFT value were significantly reduced and that some inconsistencies with experimental (e.g., NMR) results observed in the GLYCAM06 force field were resolved relevantly.

A new derivation of the GLYCAM06 force field, which removes its previous specificity for carbohydrates, and its dependency on the AMBER force field and parameters, is presented. All pertinent force field terms have been explicitly specified and so no default or generic parameters are employed. The new GLYCAM is no longer limited to any particular class of biomolecules, but is extendible to all molecular classes in the spirit of a small-molecule force field. The torsion terms in the present work were all derived from quantum mechanical data from a collection of minimal molecular fragments and related small molecules. For carbohydrates, there is now a single parameter set applicable to both α- and β-anomers and to all monosaccharide ring sizes and conformations. We demonstrate that deriving dihedral parameters by fitting to QM data for internal rotational energy curves for representative small molecules generally leads to correct rotamer populations in molecular dynamics simulations, and that this approach removes the need for phase corrections in the dihedral terms. However, we note that there are cases where this approach is inadequate. Reported here are the basic components of the new force field as well as an illustration of its extension to carbohydrates. In addition to reproducing the gas-phase properties of an array of small test molecules, condensed-phase simulations employing GLYCAM06 are shown to reproduce rotamer populations for key small molecules and representative biopolymer building blocks in explicit water, as well as crystalline lattice properties, such as unit cell dimensions, and vibrational frequencies. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2008 [ABSTRACT FROM AUTHOR]

While the construction of a dependable force field for performing classical molecular dynamics (MD) simulation is crucial for elucidating the structure and function of biomolecular systems, the attempts to do this for glycans are relatively sparse compared to those for proteins and nucleic acids. Currently, the use of GLYCAM06 force field is the most popular, but there have been a number of concerns about its accuracy in the systematic description of structural changes. In the present work, we focus on the improvement of the GLYCAM06 force field for β -d-glucose, a simple and the most abundant monosaccharide molecule, with the aid of machine learning techniques implemented with the TensorFlow library. Following the pre-sampling over a wide range of configuration space generated by MD simulation, the atomic charge and dihedral angle parameters in the GLYCAM06 force field were re-optimized to accurately reproduce the relative energies of β -d-glucose obtained by the density functional theory (DFT) calculations according to the structural changes. The validation for the newly proposed force-field parameters was then carried out by verifying that the relative energy errors compared to the DFT value were significantly reduced and that some inconsistencies with experimental (e.g., NMR) results observed in the GLYCAM06 force field were resolved relevantly. [ABSTRACT FROM AUTHOR]