On using oscillating time-dependent restraints in MD simulation

The use of time-dependent restraints in molecular simulation in order to generate a conformational ensemble for molecules that is in accordance with measured ensemble averages for particular observable quantities is investigated. Using a model system consisting of liquid butane and the cyclic peptide antamanide the reproduction of particular average 3 J-coupling constant values in a molecular dynamics simulation is analysed. It is shown that the multiple-valuedness and the sizeable gradients of the Karplus curve relating 3 J-coupling constants measured in NMR experiments to the corresponding torsional-angle values cause severe problems when trying to restrain a 3 J-coupling constant to a value close to the extrema of the Karplus curve. The introduction of a factor oscillating with time into the restraining penalty function alleviates this problem and enhances the restrained conformational samplin

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand-solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand-DNA bindin

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding

Validation of the 53A6 GROMOS force field

The quality of biomolecular dynamics simulations relies critically on the force field that is used to describe the interactions between particles in the system. Force fields, which are generally parameterized using experimental data on small molecules, can only prove themselves in realistic simulations of relevant biomolecular systems. In this work, we begin the validation of the new 53A6 GROMOS parameter set by examining three test cases. Simulations of the well-studied 129 residue protein hen egg-white lysozyme, of the DNA dodecamer d(CGCGAATTCGCG)2, and a proteinogenic β3-dodecapeptide were performed and analysed. It was found that the new parameter set performs as well as the previous parameter sets in terms of protein (45A3) and DNA (45A4) stability and that it is better at describing the folding-unfolding balance of the peptide. The latter is a property that is directly associated with the free enthalpy of hydration, to which the 53A6 parameter set was parameterize

A history of Missouri's counties, county seats, and courthouse squares (1983)

Series information taken from "1888-1984, Missouri Agricultural Experiment Station and Cooperative Extension Service publications." May not represent this version.Includes bibliographical references (pages 139-142) and index

On the use of one-step perturbation to investigate the dependence of NOE-derived atom-atom distance bound violations of peptides upon a variation of force-field parameters

The method of one-step perturbation can be used to predict from a single molecular dynamics simulation the values of observable quantities as functions of variations in the parameters of the Hamiltonian or biomolecular force field used in the simulation. The method is used to predict violations of nuclear overhauser effect (NOE) distance bounds measured in nuclear magnetic resonance (NMR) experiments by atom-atom distances of the NOE atom pairs when varying force-field parameters. Predictions of NOE distance bound violations between different versions of the GROMOS force field for a hexa-β-peptide in solution show that the technique works for rather large force-field parameter changes as well as for very different NOE bound violation patterns. The effect of changing individual force-field parameters on the NOE distance bound violations of the β-peptide and an α-peptide was investigated too. One-step perturbation, which in this case is equivalent to reweighting configurations, constitutes an efficient technique to predict many values of different quantities from a single conformational ensemble for a particular system, which makes it a powerful force-field development technique that easily reduces the number of required separate simulations by an order of magnitude

Cellular levels and molecular dynamics simulations of estragole DNA adducts point at inefficient repair resulting from limited distortion of the double-stranded DNA helix

Estragole, naturally occurring in a variety of herbs and spices, can form DNA adducts after bioactivation. Estragole DNA adduct formation and repair was studied in in vitro liver cell models, and a molecular dynamics simulation was used to investigate the conformation dependent (in)efficiency of N2-(trans-isoestragol-3′-yl)-2′-deoxyguanosine (E-3′-N2-dG) DNA adduct repair. HepG2, HepaRG cells, primary rat hepatocytes and CHO cells (including CHO wild-type and three NER-deficient mutants) were exposed to 50 μM estragole or 1′-hydroxyestragole and DNA adduct formation was quantified by LC–MS immediately following exposure and after a period of repair. Results obtained from CHO cell lines indicated that NER plays a role in repair of E-3′-N2-dG adducts, however, with limited efficiency since in the CHO wt cells 80% DNA adducts remained upon 24 h repair. Inefficiency of DNA repair was also found in HepaRG cells and primary rat hepatocytes. Changes in DNA structure resulting from E-3′-N2-dG adduct formation were investigated by molecular dynamics simulations. Results from molecular dynamics simulations revealed that conformational changes in double-stranded DNA by E-3′-N2-dG adduct formation are small, providing a possible explanation for the restrained repair, which may require larger distortions in the DNA structure. NER-mediated enzymatic repair of E-3′-N2-dG DNA adducts upon exposure to estragole will be limited, providing opportunities for accumulation of damage upon repeated daily exposure. The inability of this enzymatic repair is likely due to a limited distortion of the DNA double-stranded helix resulting in inefficient activation of nucleotide excision repair.</p

Fighting against bacterial lipopolysaccharide-caused infections through molecular dynamics simulations: a review

Lipopolysaccharide (LPS) is the primary component of the outer leaflet of Gram-negative bacterial outer membranes. LPS elicits an overwhelming immune response during infection, which can lead to life-threatening sepsis or septic shock for which no suitable treatment is available so far. As a result of the worldwide expanding multidrug-resistant bacteria, the occurrence and frequency of sepsis are expected to increase; thus, there is an urge to develop novel strategies for treating bacterial infections. In this regard, gaining an in-depth understanding about the ability of LPS to both stimulate the host immune system and interact with several molecules is crucial for fighting against LPS-caused infections and allowing for the rational design of novel antisepsis drugs, vaccines and LPS sequestration and detection methods. Molecular dynamics (MD) simulations, which are understood as being a computational microscope, have proven to be of significant value to understand LPS-related phenomena, driving and optimizing experimental research studies. In this work, a comprehensive review on the methods that can be combined with MD simulations, recently applied in LPS research, is provided. We focus especially on both enhanced sampling methods, which enable the exploration of more complex systems and access to larger time scales, and free energy calculation approaches. Thereby, apart from outlining several strategies for surmounting LPS-caused infections, this work reports the current state-of-the-art of the methods applied with MD simulations for moving a step forward in the development of such strategies.Financial support from the Spanish Ministry of Science, Innovation and Universities under the project RTI2018- 093310-B-I00 is gratefully acknowledged. C.G.F. and A.B. are also thankful for the FPU (FPU18/03525) and FPI (BES-2016-077206) postgraduate research grants, respectively

Redox thermodynamics of B-class dye-decolorizing peroxidases

With&gt;5000 annotated genes dye-decolorizing peroxidases (DyPs) represent a heme b peroxidase family of broad functional diversity. Bacterial B-class DyPs are poor peroxidases of unknown physiological function. Hydrogen peroxide efficiently mediates the rapid formation of Compound I in B-class DyPs, which, however, is stable and shows modest reactivity towards organic and inorganic electron donors. To understand these characteristics, we have investigated the redox thermodynamics of the one-electron reduction of the ferric high-spin form of wild-type B-class DyP from the pathogenic bacterium Klebsiella pneumoniae (KpDyP) and the variants D143A, R232A and D143A/R232A. These distal amino acids are fully conserved in all DyPs and play important roles in Compound I formation and maintenance of the heme cavity architecture and substrate access route(s). The E°′ values of the respective redox couples Fe(III)/Fe(II) varied from −350 mV (wild-type KpDyP) to −299 mV (D143A/R232A) at pH 7.0. Variable-temperature spectroelectrochemical experiments revealed that the reduction reaction of B-class DyPs is enthalpically unfavored but entropically favored with significant differences in enthalpic and entropic contributions to E°′ between the four proteins. Molecular dynamics simulations demonstrated the impact of solvent reorganization on the entropy change during reduction reaction and revealed the dynamics and restriction of substrate access channels. Obtained data are discussed with respect to the poor peroxidase activities of B-class DyPs and compared with heme peroxidases from other (super)families as well as with chlorite dismutases, which do not react with hydrogen peroxide but share a similar fold and heme cavity architecture