We used microsecond time scale molecular dynamics simulations to compute, at high precision, binding enthalpies for cucurbit[7]uril (CB7) with eight guests in aqueous solution. The results correlate well with experimental data from previously published isothermal titration calorimetry studies, and decomposition of the computed binding enthalpies by interaction type provides plausible mechanisti...

Calcium silicate hydrate, or C–S–H, is the chief hydration product of Portland cement. The structure of the C–S–H phase within cement has been proposed and developed via molecular simulations. In such simulations, empirical interatomic potentials for water molecules within C–S–H are adopted to govern the position and relative motion of this key constituent. Initial simulations and force fields ...

Contactins are modular extracellular cell matrix proteins that are present in the brain, and they are responsible for the proper development and functioning of neurons. They contain six immunoglobulin-like IgC2 domains and four fibronectin type III repeats. The interactions of contactin with other proteins are poorly understood. The mechanical properties of all IgC2 domains of human contactin 4...

We have estimated free energies for the binding of nine cyclic carboxylate guest molecules to the octa-acid host in the SAMPL4 blind-test challenge with four different approaches. First, we used standard free-energy perturbation calculations of relative binding affinities, performed at the molecular-mechanics (MM) level with TIP3P waters, the GAFF force field, and two different sets of charges ...

In this article, we systematically apply a novel implicit-solvent model, the variational implicit-solvent model (VISM) together with the Coulomb-Field Approximation (CFA), to calculate the hydration free energy of a large set of small organic molecules. Because these molecules have been studied in detail by molecular dynamics simulations and other implicit-solvent models, they provide a good be...

Indole-3-glycerol phosphate synthase catalyzes the terminal ring closure step in tryptophan biosynthesis. In this paper, we compare the results from molecular dynamics (MD) simulations of enzyme-bound substrate at 298, 333, 363, and 385 K and the enzyme-bound intermediate at 385 K, solvated in TIP3P water box with a CHARMM force field. Results from MD simulations agree with experimental studies...

In this work, we have evaluated how well the general assisted model building with energy refinement (AMBER) force field performs in studying the dynamic properties of liquids. Diffusion coefficients (D) have been predicted for 17 solvents, five organic compounds in aqueous solutions, four proteins in aqueous solutions, and nine organic compounds in nonaqueous solutions. An efficient sampling st...

Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against 1H-1H NOESY distances and 3J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are repo...

Does aqueous solvent discriminate among peptide conformers? To address this question, we computed the solvation free energy of a blocked, 12-residue polyalanyl-peptide in explicit water and analyzed its solvent structure. The peptide was modeled in each of 4 conformers: alpha-helix, antiparallel beta-strand, parallel beta-strand, and polyproline II helix (P(II)). Monte Carlo simulations in the ...

We have performed long molecular dynamics simulations of water using four popular water models, namely simple point charge ~SPC!, extended simple point charge ~SPC/E!, and the three point ~TIP3P! and four point ~TIP4P! transferable intermolecular potentials. System sizes of 216 and 820 molecules were used to study the dependence of properties on the system size. All systems were simulated at 30...