Fluctuations of water near extended hydrophobic and hydrophilic surfaces

We use molecular dynamics simulations of the SPC-E model of liquid water to derive probability distributions for water density fluctuations in probe volumes of different shapes and sizes, both in the bulk as well as near hydrophobic and hydrophilic surfaces. To obtain our results, we introduce a biased sampling of coarse-grained densities, which in turn biases the actual solvent density. The technique is easily combined with molecular dynamics integration algorithms. Our principal result is that the probability for density fluctuations of water near a hydrophobic surface, with or without surface-water attractions, is akin to density fluctuations at the water-vapor interface. Specifically, the probability of density depletion near the surface is significantly larger than that in bulk. In contrast, we find that the statistics of water density fluctuations near a model hydrophilic surface are similar to that in the bulk

Reaction profiling of a set of acrylamide-based human tissue transglutaminase inhibitors

The major function of the enzyme human tissue transglutaminase (TG2) is the crosslinking of proteins via a transamidation between the γ-carboxamide of a glutamine and the ε-amino group of a lysine. Overexpression of TG2 can lead to undesirable outcomes and has been linked to conditions such as fibrosis, celiac disease and neurodegenerative diseases. Accordingly, TG2 is a tempting drug target. The most effective TG2 inhibitors to date are small-molecule peptidomimetics featuring electrophilic warheads that irreversibly modify the active site catalytic cysteine (CYS277). In an effort to facilitate the design of such TG2 inhibitors, we undertook a quantum mechanical reaction profiling of the Michael reaction between a set of six acrylamide-based known TG2 inhibitors and the TG2 CYS277. The inhibitors were docked into the active site and the coordinates were refined by MD simulations prior to modelling the covalent modification of the CYS277 thiolate. The results of QM/MM MD umbrella sampling applied to reaction coordinates driving the Michael reaction are presented for two approximations of the Michael reaction: a concerted reaction (simultaneous thiolate attack onto the acrylamide warhead and pronation from the adjacent HIS335) and a two-stage reaction (consecutive thiolate attack and protonation). The two-stage approximation of the Michael reaction gave the better results for the evaluation of acrylamide-based potential TG2 inhibitors in silico. Good correlations were observed between the experimental TG2 IC50 data and the calculated activation energies over the range 0.0061 – 6.3 µM (three orders of magnitude) and we propose that this approach may be used to evaluate acrylamide-based potential TG2 inhibitors

Analyzing and Biasing Simulations with PLUMED

This chapter discusses how the PLUMED plugin for molecular dynamics can be used to analyze and bias molecular dynamics trajectories. The chapter begins by introducing the notion of a collective variable and by then explaining how the free energy can be computed as a function of one or more collective variables. A number of practical issues mostly around periodic boundary conditions that arise when these types of calculations are performed using PLUMED are then discussed. Later parts of the chapter discuss how PLUMED can be used to perform enhanced sampling simulations that introduce simulation biases or multiple replicas of the system and Monte Carlo exchanges between these replicas. This section is then followed by a discussion on how free-energy surfaces and associated error bars can be extracted from such simulations by using weighted histogram and block averaging techniques

Characterizing Structural Transitions Using Localized Free Energy Landscape Analysis

Structural changes in molecules are frequently observed during biological processes like replication, transcription and translation. These structural changes can usually be traced to specific distortions in the backbones of the macromolecules involved. Quantitative energetic characterization of such distortions can greatly advance the atomic-level understanding of the dynamic character of these biological processes.Molecular dynamics simulations combined with a variation of the Weighted Histogram Analysis Method for potential of mean force determination are applied to characterize localized structural changes for the test case of cytosine (underlined) base flipping in a GTCAGCGCATGG DNA duplex. Free energy landscapes for backbone torsion and sugar pucker degrees of freedom in the DNA are used to understand their behavior in response to the base flipping perturbation. By simplifying the base flipping structural change into a two-state model, a free energy difference of upto 14 kcal/mol can be attributed to the flipped state relative to the stacked Watson-Crick base paired state. This two-state classification allows precise evaluation of the effect of base flipping on local backbone degrees of freedom.The calculated free energy landscapes of individual backbone and sugar degrees of freedom expectedly show the greatest change in the vicinity of the flipping base itself, but specific delocalized effects can be discerned upto four nucleotide positions away in both 5' and 3' directions. Free energy landscape analysis thus provides a quantitative method to pinpoint the determinants of structural change on the atomic scale and also delineate the extent of propagation of the perturbation along the molecule. In addition to nucleic acids, this methodology is anticipated to be useful for studying conformational changes in all macromolecules, including carbohydrates, lipids, and proteins

Direct observation of topoisomerase IA gate dynamics

Type IA topoisomerases cleave single-stranded DNA and relieve negative supercoils in discrete steps corresponding to the passage of the intact DNA strand through the cleaved strand. Although type IA topoisomerases are assumed to accomplish this strand passage via a protein-mediated DNA gate, opening of this gate has never been observed. We developed a single-molecule assay to directly measure gate opening of the Escherichia coli type IA topoisomerases I and III. We found that after cleavage of single-stranded DNA, the protein gate opens by as much as 6.6 nm and can close against forces in excess of 16 pN. Key differences in the cleavage, ligation, and gate dynamics of these two enzymes provide insights into their different cellular functions. The single-molecule results are broadly consistent with conformational changes obtained from molecular dynamics simulations. These results allowed us to develop a mechanistic model of interactions between type IA topoisomerases and single-stranded DNA

Minimum Free Energy Path of Ligand-Induced Transition in Adenylate Kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme