The effect of S-substitution at the O6-guanine site on the structure and dynamics of a DNA oligomer containing a G:T mismatch

The effect of S-substitution on the O6 guanine site of a 13-mer DNA duplex containing a G:T mismatch is studied using molecular dynamics. The structure, dynamic evolution and hydration of the S-substituted duplex are compared with those of a normal duplex, a duplex with Ssubstitution on guanine, but no mismatch and a duplex with just a G:T mismatch. The S-substituted mismatch leads to cell death rather than repair. One suggestion is that the G:T mismatch recognition protein recognises the S-substituted mismatch (GS:T) as G:T. This leads to a cycle of futile repair ending in DNA breakage and cell death. We find that some structural features of the helix are similar for the duplex with the G:T mismatch and that with the S-substituted mismatch, but differ from the normal duplex, notably the helical twist. These differences arise from the change in the hydrogen-bonding pattern of the base pair. However a marked feature of the S-substituted G:T mismatch duplex is a very large opening. This showed considerable variability. It is suggested that this enlarged opening would lend support to an alternative model of cell death in which the mismatch protein attaches to thioguanine and activates downstream damage-response pathways. Attack on the sulphur by reactive oxygen species, also leading to cell death, would also be aided by the large, variable opening

An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints

In this work, we introduce an algorithm to compute the derivatives of physical observables along the constrained subspace when flexible constraints are imposed on the system (i.e., constraints in which the hard coordinates are fixed to configuration-dependent values). The presented scheme is exact, it does not contain any tunable parameter, and it only requires the calculation and inversion of a sub-block of the Hessian matrix of second derivatives of the function through which the constraints are defined. We also present a practical application to the case in which the sought observables are the Euclidean coordinates of complex molecular systems, and the function whose minimization defines the constraints is the potential energy. Finally, and in order to validate the method, which, as far as we are aware, is the first of its kind in the literature, we compare it to the natural and straightforward finite-differences approach in three molecules of biological relevance: methanol, N-methyl-acetamide and a tri-glycine peptideComment: 13 pages, 8 figures, published versio

A-Site Residues Move Independently from P-Site Residues in all-Atom Molecular Dynamics Simulations of the 70S Bacterial Ribosome

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation

Binding mode analyses and pharmacophore model development for stilbene derivatives as a novel and competitive class of α-glucosidase inhibitors

Stilbene urea derivatives as a novel and competitive class of non-glycosidic α-glucosidase inhibitors are effective for the treatment of type II diabetes and obesity. The main purposes of our molecular modeling study are to explore the most suitable binding poses of stilbene derivatives with analyzing the binding affinity differences and finally to develop a pharmacophore model which would represents critical features responsible for α-glucosidase inhibitory activity. Three-dimensional structure of S. cerevisiae α-glucosidase was built by homology modeling method and the structure was used for the molecular docking study to find out the initial binding mode of compound 12, which is the most highly active one. The initial structure was subjected to molecular dynamics (MD) simulations for protein structure adjustment at compound 12-bound state. Based on the adjusted conformation, the more reasonable binding modes of the stilbene urea derivatives were obtained from molecular docking and MD simulations. The binding mode of the derivatives was validated by correlation analysis between experimental Ki value and interaction energy. Our results revealed that the binding modes of the potent inhibitors were engaged with important hydrogen bond, hydrophobic, and π-interactions. With the validated compound 12-bound structure obtained from combining approach of docking and MD simulation, a proper four featured pharmacophore model was generated. It was also validated by comparison of fit values with the Ki values. Thus, these results will be helpful for understanding the relationship between binding mode and bioactivity and for designing better inhibitors from stilbene derivatives

Disulfide Bridges Remain Intact while Native Insulin Converts into Amyloid Fibrils

Amyloid fibrils are β-sheet-rich protein aggregates commonly found in the organs and tissues of patients with various amyloid-associated diseases. Understanding the structural organization of amyloid fibrils can be beneficial for the search of drugs to successfully treat diseases associated with protein misfolding. The structure of insulin fibrils was characterized by deep ultraviolet resonance Raman (DUVRR) and Nuclear Magnetic Resonance (NMR) spectroscopy combined with hydrogen-deuterium exchange. The compositions of the fibril core and unordered parts were determined at single amino acid residue resolution. All three disulfide bonds of native insulin remained intact during the aggregation process, withstanding scrambling. Three out of four tyrosine residues were packed into the fibril core, and another aromatic amino acid, phenylalanine, was located in the unordered parts of insulin fibrils. In addition, using all-atom MD simulations, the disulfide bonds were confirmed to remain intact in the insulin dimer, which mimics the fibrillar form of insulin

Transmembrane but not soluble helices fold inside the ribosome tunnel

Integral membrane proteins are assembled into the ER membrane via a continuous ribosome-translocon channel. The hydrophobicity and thickness of the core of the membrane bilayer leads to the expectation that transmembrane (TM) segments minimize the cost of harbouring polar polypeptide backbones by adopting a regular pattern of hydrogen bonds to form α-helices before integration. Co-translational folding of nascent chains into an α-helical conformation in the ribosomal tunnel has been demonstrated previously, but the features governing this folding are not well understood. In particular, little is known about what features influence the propensity to acquire α-helical structure in the ribosome. Using in vitro translation of truncated nascent chains trapped within the ribosome tunnel and molecular dynamics simulations, we show that folding in the ribosome is attained for TM helices but not for soluble helices, presumably facilitating SRP (signal recognition particle) recognition and/or a favourable conformation for membrane integration upon translocon entry

Crystal structure of nucleotide-free dynamin

Dynamin is a mechanochemical GTPase that oligomerizes around the neck of clathrin-coated pits and catalyses vesicle scission in a GTP-hydrolysis-dependent manner. The molecular details of oligomerization and the mechanism of the mechanochemical coupling are currently unknown. Here we present the crystal structure of human dynamin 1 in the nucleotide-free state with a four-domain architecture comprising the GTPase domain, the bundle signalling element, the stalk and the pleckstrin homology domain. Dynamin 1 oligomerized in the crystals via the stalks, which assemble in a criss-cross fashion. The stalks further interact via conserved surfaces with the pleckstrin homology domain and the bundle signalling element of the neighbouring dynamin molecule. This intricate domain interaction rationalizes a number of disease-related mutations in dynamin 2 and suggests a structural model for the mechanochemical coupling that reconciles previous models of dynamin function

Computer-Aided Lead Optimization: Improved Small-Molecule Inhibitor of the Zinc Endopeptidase of Botulinum Neurotoxin Serotype A

Optimization of a serotype-selective, small-molecule inhibitor of botulinum neurotoxin serotype A (BoNTA) endopeptidase is a formidable challenge because the enzyme-substrate interface is unusually large and the endopeptidase itself is a large, zinc-binding protein with a complex fold that is difficult to simulate computationally. We conducted multiple molecular dynamics simulations of the endopeptidase in complex with a previously described inhibitor (Kiapp of 7±2.4 µM) using the cationic dummy atom approach. Based on our computational results, we hypothesized that introducing a hydroxyl group to the inhibitor could improve its potency. Synthesis and testing of the hydroxyl-containing analog as a BoNTA endopeptidase inhibitor showed a twofold improvement in inhibitory potency (Kiapp of 3.8±0.8 µM) with a relatively small increase in molecular weight (16 Da). The results offer an improved template for further optimization of BoNTA endopeptidase inhibitors and demonstrate the effectiveness of the cationic dummy atom approach in the design and optimization of zinc protease inhibitors