Gamma estimator of Jarzynski equality for recovering binding energies from noisy dynamic data sets

A fundamental problem in thermodynamics is the recovery of macroscopic equilibrated interaction energies from experimentally measured single-molecular interactions. The Jarzynski equality forms a theoretical basis in recovering the free energy difference between two states from exponentially averaged work performed to switch the states. In practice, the exponentially averaged work value is estimated as the mean of finite samples. Numerical simulations have shown that samples having thousands of measurements are not large enough for the mean to converge when the fluctuation of external work is above 4 kBT, which is easily observable in biomolecular interactions. We report the first example of a statistical gamma work distribution applied to single molecule pulling experiments. The Gibbs free energy of surface adsorption can be accurately evaluated even for a small sample size. The values obtained are comparable to those derived from multi-parametric surface plasmon resonance measurements and molecular dynamics simulations

Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations

Molecular dynamics (MD) is one of the most common tools in computational chemistry. Recently, our group has employed accelerated molecular dynamics (aMD) to improve the conformational sampling over conventional molecular dynamics techniques. In the original aMD implementation, sampling is greatly improved by raising energy wells below a predefined energy level. Recently, our group presented an alternative aMD implementation where simulations are accelerated by lowering energy barriers of the potential energy surface. When coupled with thermodynamic integration simulations, this implementation showed very promising results. However, when applied to large systems, such as proteins, the simulation tends to be biased to high energy regions of the potential landscape. The reason for this behavior lies in the boost equation used since the highest energy barriers are dramatically more affected than the lower ones. To address this issue, in this work, we present a new boost equation that prevents oversampling of unfavorable high energy conformational states. The new boost potential provides not only better recovery of statistics throughout the simulation but also enhanced sampling of statistically relevant regions in explicit solvent MD simulations

Gating of a pH-Sensitive K2P Potassium Channel by an Electrostatic Effect of Basic Sensor Residues on the Selectivity Filter

K+ channels share common selectivity characteristics but exhibit a wide diversity in how they are gated open. Leak K2P K+ channels TASK-2, TALK-1 and TALK-2 are gated open by extracellular alkalinization. The mechanism for this alkalinization-dependent gating has been proposed to be the neutralization of the side chain of a single arginine (lysine in TALK-2) residue near the pore of TASK-2, which occurs with the unusual pKa of 8.0. We now corroborate this hypothesis by transplanting the TASK-2 extracellular pH (pHo) sensor in the background of a pHo-insensitive TASK-3 channel, which leads to the restitution of pHo-gating. Using a concatenated channel approach, we also demonstrate that for TASK-2 to open, pHo sensors must be neutralized in each of the two subunits forming these dimeric channels with no apparent cross-talk between the sensors. These results are consistent with adaptive biasing force analysis of K+ permeation using a model selectivity filter in wild-type and mutated channels. The underlying free-energy profiles confirm that either a doubly or a singly charged pHo sensor is sufficient to abolish ion flow. Atomic detail of the associated mechanism reveals that, rather than a collapse of the pore, as proposed for other K2P channels gated at the selectivity filter, an increased height of the energetic barriers for ion translocation accounts for channel blockade at acid pHo. Our data, therefore, strongly suggest that a cycle of protonation/deprotonation of pHo-sensing arginine 224 side chain gates the TASK-2 channel by electrostatically tuning the conformational stability of its selectivity filter

Parameterization of a coarse-grained model of cholesterol with point-dipole electrostatics

© 2018, Springer Nature Switzerland AG. We present a new coarse-grained (CG) model of cholesterol (CHOL) for the electrostatic-based ELBA force field. A distinguishing feature of our CHOL model is that the electrostatics is modeled by an explicit point dipole which interacts through an ideal vacuum permittivity. The CHOL model parameters were optimized in a systematic fashion, reproducing the electrostatic and nonpolar partitioning free energies of CHOL in lipid/water mixtures predicted by full-detailed atomistic molecular dynamics simulations. The CHOL model has been validated by comparison to structural, dynamic and thermodynamic properties with experimental and atomistic simulation reference data. The simulation of binary DPPC/cholesterol mixtures covering the relevant biological content of CHOL in mammalian membranes is shown to correctly predict the main lipid behavior as observed experimentally