Proper balance of solvent-solute and solute-solute interactions in the treatment of the diffusion of glucose using the Drude polarizable force field.

Motivated by underestimation of the diffusion constant of glucose in aqueous solution at high glucose concentrations we performed additional optimization of the Drude polarizable hexopyranose monosaccharide force field. This indicated aggregation of the glucose at higher concentrations, which is a concern for studies of complex glycan systems such as the HIV Envelope where high effective concentrations of sugars are present. High-level quantum mechanical calculations were undertaken on water monohydrate-glucose interactions, on water cluster-glucose interactions and on glucose-glucose dimers in stacked (parallel) and perpendicular orientations. Optimization of the nonbond and dihedral parameters targeting these data yielded a revised model that showed improved agreement with experimental aqueous diffusion data. However, limitations in the diffusion constants were still present. These were due to the SWM4-NDP inherently overestimating the diffusion constant of water, a problem that was validated by calculation of the aqueous diffusion constants using the SWM6-NDP water model. In addition, results show the water diffusion constant to be significantly overestimated at high glucose concentrations though the glucose diffusion is in satisfactory agreement with experiment. These results indicate the subtle balance of water-sugar, water-water and sugar-sugar interactions that needs to be properly modeled to account for the full range of aqueous behavior of sugars in aqueous solution.

[1]  Alexander D. MacKerell,et al.  Polarizable Force Field for DNA Based on the Classical Drude Oscillator: I. Refinement Using Quantum Mechanical Base Stacking and Conformational Energetics. , 2017, Journal of chemical theory and computation.

[2]  E. Carballo,et al.  Thermodynamic Properties of Aqueous Carbohydrate Solutions , 1997 .

[3]  Alexander D. MacKerell,et al.  A comparative Kirkwood-Buff study of aqueous methanol solutions modeled by the CHARMM additive and Drude polarizable force fields. , 2013, The journal of physical chemistry. B.

[4]  Wonpil Im,et al.  Effects of N-Glycan Composition on Structure and Dynamics of IgG1 Fc and Their Implications for Antibody Engineering , 2017, Scientific Reports.

[5]  Alexander D. MacKerell,et al.  Differential Impact of the Monovalent Ions Li+, Na+, K+, and Rb+ on DNA Conformational Properties , 2014, The journal of physical chemistry letters.

[6]  Martin Frank,et al.  Bioinformatics and molecular modeling in glycobiology , 2010, Cellular and Molecular Life Sciences.

[7]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[8]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[9]  Alexander D. MacKerell,et al.  All‐atom polarizable force field for DNA based on the classical drude oscillator model , 2014, J. Comput. Chem..

[10]  Alexander D. MacKerell,et al.  Induced Polarization Influences the Fundamental Forces in DNA Base Flipping , 2014, The journal of physical chemistry letters.

[11]  Peter Pulay,et al.  Systematic AB Initio Gradient Calculation of Molecular Geometries, Force Constants, and Dipole Moment Derivatives , 1979 .

[12]  Benoît Roux,et al.  A polarizable force field of dipalmitoylphosphatidylcholine based on the classical Drude model for molecular dynamics simulations of lipids. , 2013, The journal of physical chemistry. B.

[13]  Robin Taylor,et al.  Life-science applications of the Cambridge Structural Database. , 2002, Acta crystallographica. Section D, Biological crystallography.

[14]  Alexander D. MacKerell,et al.  Additive empirical force field for hexopyranose monosaccharides , 2008, J. Comput. Chem..

[15]  Alexander D. MacKerell,et al.  Drude polarizable force field for aliphatic ketones and aldehydes, and their associated acyclic carbohydrates , 2017, Journal of Computer-Aided Molecular Design.

[16]  Philippe H. Hünenberger,et al.  A reoptimized GROMOS force field for hexopyranose‐based carbohydrates accounting for the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers , 2011, J. Comput. Chem..

[17]  Alexander D. MacKerell,et al.  Induction of peptide bond dipoles drives cooperative helix formation in the (AAQAA)3 peptide. , 2014, Biophysical journal.

[18]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[19]  Alexander D. MacKerell,et al.  Automated conformational energy fitting for force-field development , 2008, Journal of molecular modeling.

[20]  Alexander D. MacKerell,et al.  Polarizable empirical force field for the primary and secondary alcohol series based on the classical Drude model. , 2007, Journal of chemical theory and computation.

[21]  Alexander D. MacKerell,et al.  A polarizable model of water for molecular dynamics simulations of biomolecules , 2006 .

[22]  Alexander D. MacKerell,et al.  Induced Dipole–Dipole Interactions Influence the Unfolding Pathways of Wild-Type and Mutant Amyloid β-Peptides , 2015, Journal of Physical Chemistry B.

[23]  Gerhard Hummer,et al.  System-Size Dependence of Diffusion Coefficients and Viscosities from Molecular Dynamics Simulations with Periodic Boundary Conditions , 2004 .

[24]  Malcolm Dole,et al.  Diffusion in Supersaturated Solutions. II. Glucose Solutions , 1953 .

[25]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[26]  Alexander D. MacKerell,et al.  Polarizable empirical force field for nitrogen‐containing heteroaromatic compounds based on the classical Drude oscillator , 2009, J. Comput. Chem..

[27]  Udo Kaatze,et al.  Molecular Dynamics of Carbohydrate Aqueous Solutions. Dielectric Relaxation as a Function of Glucose and Fructose Concentration , 2001 .

[28]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[29]  Alexander D. MacKerell,et al.  Polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator. , 2013, Biopolymers.

[30]  Alexander D. MacKerell,et al.  An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications , 2016, Chemical reviews.

[31]  Victor H Rusu,et al.  GROMOS 53A6GLYC, an Improved GROMOS Force Field for Hexopyranose-Based Carbohydrates. , 2012, Journal of chemical theory and computation.

[32]  Alexander D. MacKerell,et al.  Recent Advances in Polarizable Force Fields for Macromolecules: Microsecond Simulations of Proteins Using the Classical Drude Oscillator Model , 2014, The journal of physical chemistry letters.

[33]  Xiao Zhu,et al.  Polarizable empirical force field for sulfur‐containing compounds based on the classical Drude oscillator model , 2010, J. Comput. Chem..

[34]  Alexander D. MacKerell,et al.  Polarizable empirical force field for aromatic compounds based on the classical drude oscillator. , 2007, The journal of physical chemistry. B.

[35]  Alexander D. MacKerell,et al.  CHARMM Drude Polarizable Force Field for Aldopentofuranoses and Methyl-aldopentofuranosides. , 2015, The journal of physical chemistry. B.

[36]  Alexander D. MacKerell,et al.  Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field. , 2010, Journal of chemical theory and computation.

[37]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[38]  Alexander D. MacKerell,et al.  Robustness in the fitting of molecular mechanics parameters , 2015, J. Comput. Chem..

[39]  Alexander D. MacKerell,et al.  Drude Polarizable Force Field for Molecular Dynamics Simulations of Saturated and Unsaturated Zwitterionic Lipids. , 2017, Journal of chemical theory and computation.

[40]  Lai-Xi Wang,et al.  Conformational Heterogeneity of the HIV Envelope Glycan Shield , 2017, Scientific Reports.

[41]  D. Cremer,et al.  General definition of ring puckering coordinates , 1975 .

[42]  C. Talón,et al.  Dynamics of Water Molecules in Glucose Solutions , 2004 .

[43]  Alexander D. MacKerell,et al.  Polarizable Force Field for DNA Based on the Classical Drude Oscillator: II. Microsecond Molecular Dynamics Simulations of Duplex DNA. , 2017, Journal of chemical theory and computation.

[44]  Alexander D. MacKerell,et al.  Polarizable Empirical Force Field for Hexopyranose Monosaccharides Based on the Classical Drude Oscillator , 2014, The journal of physical chemistry. B.

[45]  Alexander D. MacKerell,et al.  Six-site polarizable model of water based on the classical Drude oscillator. , 2013, The Journal of chemical physics.

[46]  Alexander D. MacKerell,et al.  Force Field for Peptides and Proteins based on the Classical Drude Oscillator. , 2013, Journal of chemical theory and computation.