Molecular dynamics simulation using enhanced sampling methods is one of the powerful computational tools used to explore protein conformations and free energy landscapes. Enhanced sampling methods often employ either an increase in temperature or a flattening of the potential energy surface to rapidly sample phase space, and a corresponding reweighting algorithm is used to recover the Boltzmann statistics. However, potential energies of complex biomolecules usually involve large fluctuations on a magnitude of hundreds of kcal/mol despite minimal structural changes during simulation. This leads to noisy reweighting statistics and complicates the obtainment of accurate final results. To overcome this common issue in enhanced conformational sampling, we propose a scaled molecular dynamics method, which modifies the biomolecular potential energy surface and employs a reweighting scheme based on configurational populations. Statistical mechanical theory is applied to derive the reweighting formula, and the canonical ensemble of simulated structures is recovered accordingly. Test simulations on alanine dipeptide and the fast folding polypeptide Chignolin exhibit sufficiently enhanced conformational sampling and accurate recovery of free energy surfaces and thermodynamic properties. The results are comparable to long conventional molecular dynamics simulations and exhibit better recovery of canonical statistics over methods which employ a potential energy term in reweighting.

An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038-2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pK(a) values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pK(a) values in closer agreement to those obtained experimentally for the more slowly converging residues.

We use thermodynamic integration (TI) and explicit solvent molecular dynamics (MD) simulation to estimate the absolute free energy of host-guest binding. In the unbound state, water molecules visit all of the internally accessible volume of the host, which is fully hydrated on all sides. Upon binding of an apolar guest, the toroidal host cavity is fully dehydrated; thus, during the intermediate λ stages along the integration, the hydration of the host fluctuates between hydrated and dehydrated states. Estimating free energies by TI can be especially challenging when there is a considerable difference in hydration between the two states of interest. We investigate these aspects using the popular TIP3P and TIP4P water models. TI free energy estimates through MD largely depend on water-related interactions, and water dynamics significantly affect the convergence of binding free energy calculations. Our results indicate that wetting/dewetting transitions play a major role in slowing the convergence of free energy estimation. We employ two alternative approaches-one analytical and the other empirically based on actual MD sampling-to correct for the standard state free energy. This correction is sizable (up to 4 kcal/mol), and the two approaches provide corrections that differ by about 1 kcal/mol. For the system considered here, the TIP4P water model combined with an analytical correction for the standard state free energy provides higher overall accuracy. This observation might be transferable to other systems in which water-related contributions dominate the binding process.

The synthetic host cucurbit[7]uril (CB[7]) binds aromatic guests or metal complexes with ultrahigh affinity compared with that typically displayed in protein-ligand binding. Due to its small size, CB[7] serves as an ideal receptor-ligand system for developing computational methods for molecular recognition. Here, we apply the recently developed variational implicit-solvent model (VISM), numerically evaluated by the level-set method, to study hydration effects in the high-affinity binding of the B2 bicyclo[2.2.2]octane derivative to CB[7]. For the unbound host, we find that the host cavity favors the hydrated state over the dry state due to electrostatic effects. For the guest binding, we find reasonable agreement to experimental binding affinities. Dissection of the individual VISM free-energy contributions shows that the major driving forces are water-mediated hydrophobic interactions and the intrinsic (vacuum) host-guest van der Waals interactions. These findings are in line with recent experiments and molecular dynamics simulations with explicit solvent. It is expected that the level-set VISM, with further refinement on the electrostatic descriptions, can efficiently predict molecular binding and recognition in a wide range of future applications.