Impact of interfacial high-density water layer on accurate estimation of adsorption free energy by Jarzynski's equality

Entropy-Enthalpy Compensation in Peptide Adsorption on Solid Surfaces: Dependence on Surface Hydration

Although protein adsorption at the solid-water interface is of immense importance, understanding the crucial role of the water phase in mediating protein-surface interactions is lacking, particularly due to the lack of fundamental thermodynamic data. Herein, we have performed complicated free energy calculations and successfully extracted the entropy and enthalpy changes of molecular adsorption on solids. Using the gold and graphene as the surface models with distinct affinities to the water phase, we successfully unravel the sharply opposite manners of entropy-enthalpy compensation in driving water and tripeptide adsorptions on two surfaces. Though the thermodynamic features of water adsorption on surface are enthalpically dominated based on the positions of free energy barriers and minima, the favorable entropy term significantly decreases the free energy barrier and further stabilizes the adsorbate at the adsorption site on the graphene surface. For the peptide, the shape of the adsorption free energy profile is jointly determined by the enthalpy and entropy changes, which, however, alternatively act the driving force to promote the peptide adsorption on the Au surface and graphene surface. The distinct structural and dynamic properties of solid-liquid interfaces account for the special role of the interfacial water phase in regulating the competitive relationship between the entropy and enthalpy variations.

Computational Study for the Unbinding Routes of β-N-Acetyl-d-Hexosaminidase Inhibitor: Insight from Steered Molecular Dynamics Simulations

β-N-Acetyl-d-hexosaminidase from Ostrinia furnacalis (OfHex1) is a new target for the design of insecticides. Although some of its inhibitors have been found, there is still no commercial drug available at present. The residence time of the ligand may be important for its pharmacodynamic effect. However, the unbinding routes of ligands from OfHex1 still remain largely unexplored. In the present study, we first simulated the six dissociation routes of N,N,N-trimethyl-d-glucosamine-chitotriomycin (TMG-chitotriomycin, a highly selective inhibitor of OfHex1) from the active pocket of OfHex1 by steered molecular dynamics simulations. By comparing the potential of mean forces (PMFs) of six routes, Route 1 was considered as the most possible route with the lowest energy barrier. Furthermore, the structures of six different states for Route 1 were snapshotted, and the key amino acid residues affecting the dissociated time were analyzed in the unbinding pathway. Moreover, we also analyzed the "open⁻close" mechanism of Glu368 and Trp448 and found that their conformational changes directly affected the dissociation of TMG-chitotriomycin. Our findings would be helpful to understanding and identifying novel inhibitors against OfHex1 from virtual screening or lead-optimization.

Elastic Barrier Dynamical Freezing in Free Energy Calculations: A Way To Speed Up Nonequilibrium Molecular Dynamics Simulations by Orders of Magnitude

An important issue concerning computer simulations addressed to free energy estimates via nonequilibrium work theorems, such as the Jarzynski equality [Phys. Rev. Lett. 1997, 78, 2690], is the computational effort required to achieve results with acceptable accuracy. In this respect, the dynamical freezing approach [Phys. Rev. E 2009, 80, 041124] has been shown to improve the efficiency of this kind of simulations, by blocking the dynamics of particles located outside an established mobility region. In this report, we show that dynamical freezing produces a systematic spurious decrease of the particle density inside the mobility region. As a consequence, the requirements to apply nonequilibrium work theorems are only approximately met. Starting from these considerations, we have developed a simulation scheme, called "elastic barrier dynamical freezing", according to which a stiff potential-energy barrier is enforced at the boundaries of the mobility region, preventing the particles from leaving this region of space during the nonequilibrium trajectories. The method, tested on the calculation of the distance-dependent free energy of a dimer immersed into a Lennard-Jones fluid, provides an accuracy comparable to the conventional steered molecular dynamics, with a computational speedup exceeding a few orders of magnitude.