Dataset Information

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

ABSTRACT: 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.

SUBMITTER: Yang M

PROVIDER: S-EPMC5809005 | biostudies-literature | 2018 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

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

Yang Mingjun M Aytenfisu Asaminew H AH MacKerell Alexander D AD

Carbohydrate research 20180131

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 mon ...[more]

PMID: 29422120

Similar Datasets

Project description:The CHARMM Drude-2013 polarizable force field (FF) was developed to include the explicit treatment of induced electronic polarizability, resulting in a more accurate description of the electrostatic interactions in molecular dynamics (MD) simulations. While the Drude-2013 protein FF has shown success in improving the folding properties of α-helical peptides and to reproduce experimental observables in simulations up to 1 μs, some limitations were noted regarding the stability of β-sheet structures in simulations longer than 100 ns as well as larger deviations from crystal structures in simulations of a number of proteins compared to the additive CHARMM36 protein FF. The origin of the instability has been identified and appears to be primarily due to overestimated atomic polarizabilities and induced dipole-dipole interactions on the Cβ, Cγ, and Cδ side chain atoms. To resolve this and other issues, a number of aspects of the model were revisited, resulting in Drude-2019 protein FF. Backbone parameters were optimized targeting the conformational properties of the (Ala)5 peptide in solution along with gas phase properties of the alanine dipeptide. Dipeptides that contain N-acetylated and N'-methylamidated termini, excluding Gly, Pro, and Ala, were used as models to optimize the atomic polarizabilities and Thole screening factors on selected Cβ, Cγ, and Cδ carbons by targeting quantum mechanical (QM) dipole moments and molecular polarizabilities. In addition, to obtain better conformational properties, side chain χ1 and χ2 dihedral parameters were optimized targeting QM data for the respective side chain dipeptide conformations as well as Protein Data Bank survey data based on the χ1, χ2 sampling from Hamiltonian replica-exchange MD simulations of (Ala)4-X-(Ala)4 in solution, where X is the amino acid of interest. Further improvements include optimizing nonbonded interactions between charged residues to reproduce QM interaction energies of the charged-protein model compounds and experimental osmotic pressures. Validation of the optimized Drude protein FF includes MD simulations of a collection of peptides and proteins including β-sheet structures, as well as transmembrane ion channels. Results showed that the updated Drude-2019 protein FF yields smaller overall root-mean-square differences of proteins as compared to the additive CHARMM36m and Drude-2013 FFs as well as similar or improved agreement with experimental NMR properties, allowing for long time scale simulation studies of proteins and more complex biomolecular systems in conjunction with the remainder of the Drude polarizable FF.

Project description:Molecular dynamics (MD) simulations based on atomic models play an important role in the drug-discovery process to screen molecules, estimate binding free energies, and optimize lead compounds in chemical space. Accurate computations of thermodynamic and kinetic properties using MD simulations are highly dependent on the accuracy of the underlying atomic force field. In this context, going beyond the nonpolarizable fixed-charge model by accounting explicitly for induced polarization is highly desirable. The CHARMM polarizable force field based on classical Drude oscillators, in which an auxiliary charged particle is attached via a harmonic spring to its parent nucleus, offers both a computationally convenient and rigorous framework to model explicitly induced electronic polarization in MD simulations. For any molecule of interest, electrostatic partial charges, atomic polarizabilities, and Thole shielding factors, as well as bonded parameters can either be determined from ab initio calculations or ascribed from the knowledge-based library of the CHARMM Generalized force field. While this approach is fairly reliable in general, it is well understood that the overall accuracy of the models with respect to thermodynamic properties such as bulk density, enthalpies, and solvation free energies is particularly sensitive to the nonbonded Lennard-Jones (LJ) parameters. In the present study, we systematically refined the set of LJ parameters for the atom types available in the Drude force field to best match the experimental thermodynamic properties for 416 small druglike organic molecules. To further test the transferability of the optimized parameters, the hydration free energy of 372 molecules was computed. The calculations resulted in a small average error of 0.46 kcal/mol and a Pearson R of 0.9, representing a significant improvement over the additive GAFF force field in our previous study, where an average error of ∼2 kcal/mol was obtained. Such an improvement is consistent with the ability of the polarizable Drude model to more accurately model interactions in different environments. The effort provides a roadmap for the global optimization of force field parameters using experimental data. It is hoped that the present effort will further the application of the Drude polarizable force field in molecular simulations including drug design and discovery.

Dataset Information

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

Publications

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets