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

Binding mechanism of oligopeptides on solid surface: assessing the significance of single-molecule approach

This paper addresses the complementarity and potential disparities between single-molecule and ensemble-average approaches to probe the binding mechanism of oligopeptides on inorganic solids. Specifically, we explore the peptide/gold interface owing to its significance in various topics and its suitability to perform experiments both in model and real conditions. Experimental results show that the studied peptide adopts a lying configuration upon adsorption on the gold surface and interacts through its peptidic links and deprotonated thiolate extremities, in agreement with theoretical predictions. Single-molecule force spectroscopy (SMFS) measurements revealed the existence of a wide panel of adhesion forces, resulting from the interaction between individual peptide moieties and the abundant surface sites. We therefore propose methodological developments for sorting the events of interest to understand the peptide adsorption mechanism. Thermodynamic and kinetic aspects of the peptide adsorption are probed using both static and dynamic force spectroscopy measurements. Specifically, we show the possibility of providing a reasonable estimate of the peptide free energy of adsorption ΔadsG° by exploring the fluctuations of the adhesion work, based on the Jarzynski equality, and by using a parametric Gamma estimator. The proposed approach offers a relevant method for studying the different factors influencing the peptide adsorption and evaluating their impact on ΔadsG° as an alternative to exploring adhesion forces that may lead to misinterpretations. This is illustrated by the comparison of the adsorption of two peptides with specific amino acids substitution. Our method provides insights into the overall mechanism by which peptides interact with the surface and allows an integration of the single-molecule versus ensemble-average points of view.

Conformational dynamics and putative substrate extrusion pathways of the N-glycosylated outer membrane factor CmeC from Campylobacter jejuni

The outer membrane factor CmeC of the efflux machinery CmeABC plays an important role in conferring antibiotic and bile resistance to Campylobacter jejuni. Curiously, the protein is N-glycosylated, with the glycans playing a key role in the effective function of this system. In this work we have employed atomistic equilibrium molecular dynamics simulations of CmeC in a representative model of the C. jejuni outer membrane to characterise the dynamics of the protein and its associated glycans. We show that the glycans are more conformationally labile than had previously been thought. The extracellular loops of CmeC visit the open and closed states freely suggesting the absence of a gating mechanism on this side, while the narrow periplasmic entrance remains tightly closed, regulated via coordination to solvated cations. We identify several cation binding sites on the interior surface of the protein. Additionally, we used steered molecular dynamics simulations to elucidate translocation pathways for a bile acid and a macrolide antibiotic. These, and additional equilibrium simulations suggest that the anionic bile acid utilises multivalent cations to climb the ladder of acidic residues that line the interior surface of the protein.

Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules

Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.

Molecular Simulations of Surfactant Adsorption on Iron Oxide from Hydrocarbon Solvents

The performance of organic friction modifiers (OFMs) depends on their ability to adsorb onto surfaces and form protective monolayers. Understanding the relationship between OFM concentration in the base oil and the resulting surface coverage is important for improving lubricant formulations. Here, we use molecular dynamics (MD) simulations to study the adsorption of three OFMs─stearic acid (SA), glycerol monoostearate (GMS), and glycerol monooleate (GMO)─onto a hematite surface from two hydrocarbon solvents─n-hexadecane and poly(α-olefin) (PAO). We calculate the potential of mean force of the adsorption process using the adaptive biasing force algorithm, and the adsorption strength increases in the order SA < GMS < GMO. We estimate the minimum area occupied by OFM molecules on the surface using annealing MD simulations and obtained a similar hard-disk area for GMS and GMO but a lower value for SA. Using the MD results, we determine the adsorption isotherms using the molecular thermodynamic theory (MTT), which agree well with one previous experimental data set for SA on hematite. For two other experimental data sets for SA, lateral interactions between surfactant molecules need to be accounted for within the MTT framework. SA forms monolayers with lower surface coverage than GMO and GMS at low concentrations but also has the highest plateau coverage. We validate the adsorption energies from the MD simulations using high-frequency reciprocating rig friction experiments with different concentrations of the OFMs in PAO. For OFMs with saturated tailgroups (SA and GMS), we obtain good agreement between the simulations and the experiments. The results deviate for OFMs containing Z-unsaturated tailgroups (GMO) due to the additional steric hindrance, which is not accounted for in the current simulation framework. This study demonstrates that MD simulations, alongside MTT, are an accurate and efficient tool to predict adsorption isotherms at solid-liquid interfaces.