Browse
Submit Data
Databases
API
Help

Dataset Information

28 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

On the Use of Molecular Dynamics Simulations for Probing Allostery through DNA.

ABSTRACT: A recent study described an allosteric effect in which the binding of a protein to DNA is influenced by another protein bound nearby. The effect shows a periodicity of ?10 basepairs and decays with increasing protein-protein distance. As a mechanistic explanation, the authors reported a similar periodic, decaying pattern of the correlation coefficient between major groove widths inferred from a shorter molecular dynamics simulation. Here we show that in a state-of-the-art, microsecond-long simulation of the same DNA sequence, the periodicity of the correlation coefficient is not observed. To study the problem further, we extend an earlier mechanical model of DNA allostery based on constrained minimization of effective quadratic deformation energy of the DNA. We demonstrate that, if the constraints mimicking the bound proteins are properly applied, the periodicity in the binding energy is indeed recovered.

SUBMITTER: Drsata T

PROVIDER: S-EPMC4776040 | biostudies-literature | 2016 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

On the Use of Molecular Dynamics Simulations for Probing Allostery through DNA.

Dršata Tomáš T Zgarbová Marie M Jurečka Petr P Šponer Jiří J Lankaš Filip F

Biophysical journal 20160127 4

A recent study described an allosteric effect in which the binding of a protein to DNA is influenced by another protein bound nearby. The effect shows a periodicity of ∼10 basepairs and decays with increasing protein-protein distance. As a mechanistic explanation, the authors reported a similar periodic, decaying pattern of the correlation coefficient between major groove widths inferred from a shorter molecular dynamics simulation. Here we show that in a state-of-the-art, microsecond-long simul ...[more]

PMID: 26827073

Similar Datasets

Probing allostery through DNA.

Project description:Allostery is well documented for proteins but less recognized for DNA-protein interactions. Here, we report that specific binding of a protein on DNA is substantially stabilized or destabilized by another protein bound nearby. The ternary complex's free energy oscillates as a function of the separation between the two proteins with a periodicity of ~10 base pairs, the helical pitch of B-form DNA, and a decay length of ~15 base pairs. The binding affinity of a protein near a DNA hairpin is similarly dependent on their separation, which-together with molecular dynamics simulations-suggests that deformation of the double-helical structure is the origin of DNA allostery. The physiological relevance of this phenomenon is illustrated by its effect on gene expression in live bacteria and on a transcription factor's affinity near nucleosomes.

| S-EPMC3586787 | biostudies-literature

Allostery of actin filaments: molecular dynamics simulations and coarse-grained analysis.

Project description:The structural and mechanical properties of monomeric actin (G-actin), the trimer nucleus, and actin filaments (F-actins) are determined as a function of the conformation of the DNase I-binding loop (DB loop) by using all-atom molecular dynamics simulations and coarse-grained (CG) analysis. Recent x-ray structures of ADP-bound G-actin (G-ADP) by Otterbein et al. [Otterbein, L. R., Graceffa, P. & Dominguez, R. (2001) Science 293, 708-711] and ATP-bound G-actin (G-ATP) by Graceffa and Dominguez [Graceffa, P. & Dominguez, R. (2003) J. Biol. Chem. 278, 34172-34180] indicate that the DB loop of actin does not have a well defined secondary structure in the ATP state but folds into an alpha-helix in the ADP state. MD simulations and CG analysis indicate that such a helical DB loop significantly weakens the intermonomer interactions of actin assemblies and thus leads to a wider, shorter, and more disordered filament. The computed persistence lengths of F-actin composed of G-ATP (16 microm) and of G-ADP (8.5 microm) agree well with the experimental values for the two states. Therefore, the loop-to-helix transition of the DB loop may be one of the factors that lead to the changes in structural and mechanical properties of F-actin after ATP hydrolysis. This result may provide a direct connection between the conformational changes of an actin monomer and the structural and mechanical properties of the cytoskeleton. The information provided by MD simulations also helps to understand the possible origin of the special features of actin dynamics.

| S-EPMC1201585 | biostudies-literature

Molecular dynamics simulations of aptamer-binding reveal generalized allostery in thrombin.

Project description:Thrombin is an attractive target for antithrombotic therapy due to its central role in thrombosis and hemostasis as well as its role in inducing tumor growth, metastasis, and tumor invasion. The thrombin-binding DNA aptamer (TBA), is under investigation for anticoagulant drugs. Although aptamer binding experiments have been revealed various effects on thrombin's enzymatic activities, the detailed picture of the thrombin's allostery from TBA binding is still unclear. To investigate thrombin's response to the aptamer-binding at the molecular level, we compare the mechanical properties and free energy landscapes of the free and aptamer-bound thrombin using microsecond-scale all-atom GPU-based molecular dynamics simulations. Our calculations on residue fluctuations and coupling illustrate the allosteric effects of aptamer-binding at the atomic level, highlighting the exosite II, 60s, ? and the sodium loops, and the alpha helix region in the light chains involved in the allosteric changes. This level of details clarifies the mechanisms of previous experimentally demonstrated phenomena, and provides a prediction of the reduced autolysis rate after aptamer-binding. The shifts in thrombin's ensemble of conformations and free energy surfaces after aptamer-binding demonstrate that the presence of bound-aptamer restricts the conformational freedom of thrombin suggesting that conformational selection, i.e. generalized allostery, is the dominant mechanism of thrombin-aptamer binding. The profound perturbation on thrombin's mechanical and thermodynamic properties due to the aptamer-binding, which was revealed comprehensively as a generalized allostery in this work, may be exploited in further drug discovery and development.

| S-EPMC6876308 | biostudies-literature

Probing antibody internal dynamics with fluorescence anisotropy and molecular dynamics simulations.

Project description:The solution dynamics of antibodies are critical to antibody function. We explore the internal solution dynamics of antibody molecules through the combination of time-resolved fluorescence anisotropy experiments on IgG1 with more than two microseconds of all-atom molecular dynamics (MD) simulations in explicit water, an order of magnitude more than in previous simulations. We analyze the correlated motions with a mutual information entropy quantity, and examine state transition rates in a Markov-state model, to give coarse-grained descriptors of the motions. Our MD simulations show that while there are many strongly correlated motions, antibodies are highly flexible, with F(ab) and F(c) domains constantly forming and breaking contacts, both polar and non-polar. We find that salt bridges break and reform, and not always with the same partners. While the MD simulations in explicit water give the right time scales for the motions, the simulated motions are about 3-fold faster than the experiments. Overall, the picture that emerges is that antibodies do not simply fluctuate around a single state of atomic contacts. Rather, in these large molecules, different atoms come in contact during different motions.

| S-EPMC3893241 | biostudies-literature

Probing the disordered domain of the nuclear pore complex through coarse-grained molecular dynamics simulations.

Project description:The distribution of disordered proteins (FG-nups) that line the transport channel of the nuclear pore complex (NPC) is investigated by means of coarse-grained molecular dynamics simulations. A one-bead-per-amino-acid model is presented that accounts for the hydrophobic/hydrophilic and electrostatic interactions between different amino acids, polarity of the solvent, and screening of free ions. The results indicate that the interaction of the FG-nups forms a high-density, doughnut-like distribution inside the NPC, which is rich in FG-repeats. We show that the obtained distribution is encoded in the amino-acid sequence of the FG-nups and is driven by both electrostatic and hydrophobic interactions. To explore the relation between structure and function, we have systematically removed different combinations of FG-nups from the pore to simulate inviable and viable NPCs that were previously studied experimentally. The obtained density distributions show that the maximum density of the FG-nups inside the pore does not exceed 185 mg/mL in the inviable NPCs, whereas for the wild-type and viable NPCs, this value increases to 300 mg/mL. Interestingly, this maximum density is not correlated to the total mass of the FG-nups, but depends sensitively on the specific combination of essential Nups located in the central plane of the NPC.

| S-EPMC4167297 | biostudies-literature

In situ structure and dynamics of DNA origami determined through molecular dynamics simulations.

Project description:The DNA origami method permits folding of long single-stranded DNA into complex 3D structures with subnanometer precision. Transmission electron microscopy, atomic force microscopy, and recently cryo-EM tomography have been used to characterize the properties of such DNA origami objects, however their microscopic structures and dynamics have remained unknown. Here, we report the results of all-atom molecular dynamics simulations that characterized the structural and mechanical properties of DNA origami objects in unprecedented microscopic detail. When simulated in an aqueous environment, the structures of DNA origami objects depart from their idealized targets as a result of steric, electrostatic, and solvent-mediated forces. Whereas the global structural features of such relaxed conformations conform to the target designs, local deformations are abundant and vary in magnitude along the structures. In contrast to their free-solution conformation, the Holliday junctions in the DNA origami structures adopt a left-handed antiparallel conformation. We find the DNA origami structures undergo considerable temporal fluctuations on both local and global scales. Analysis of such structural fluctuations reveals the local mechanical properties of the DNA origami objects. The lattice type of the structures considerably affects global mechanical properties such as bending rigidity. Our study demonstrates the potential of all-atom molecular dynamics simulations to play a considerable role in future development of the DNA origami field by providing accurate, quantitative assessment of local and global structural and mechanical properties of DNA origami objects.

| S-EPMC3864285 | biostudies-literature

Probing energetics of Abeta fibril elongation by molecular dynamics simulations.

Project description:Using replica exchange molecular dynamics simulations and an all-atom implicit solvent model, we probed the energetics of Abeta(10-40) fibril growth. The analysis of the interactions between incoming Abeta peptides and the fibril led us to two conclusions. First, considerable variations in fibril binding propensities are observed along the Abeta sequence. The peptides in the fibril and those binding to its edge interact primarily through their N-termini. Therefore, the mutations affecting the Abeta positions 10-23 are expected to have the largest impact on fibril elongation compared with those occurring in the C-terminus and turn. Second, we performed weak perturbations of the binding free energy landscape by scanning partial deletions of side-chain interactions at various Abeta sequence positions. The results imply that strong side-chain interactions--in particular, hydrophobic contacts--impede fibril growth by favoring disordered docking of incoming peptides. Therefore, fibril elongation may be promoted by moderate reduction of Abeta hydrophobicity. The comparison with available experimental data is presented.

| S-EPMC2711452 | biostudies-literature

Probing mutation-induced conformational transformation of the GTP/M-RAS complex through Gaussian accelerated molecular dynamics simulations.

Project description:Mutations highly affect the structural flexibility of two switch domains in M-RAS considered an important target of anticancer drug design. Gaussian accelerated molecular dynamics (GaMD) simulations were applied to probe the effect of mutations P40D, D41E, and P40D/D41E/L51R on the conformational transition of the switch domains from the GTP-bound M-RAS. The analyses of free energy landscapes (FELs) not only reveal that three mutations induce less energetic states than the wild-type (WT) M-RAS but also verify that the switch domains are extremely disordered. Principal component analysis (PCA) and dynamics analysis suggest that three mutations greatly affect collective motions and structural flexibility of the switch domains that mostly overlap with binding regions of M-RAS to its effectors, which in turn disturbs the activity of M-RAS. The analyses of the interaction network between GTP and M-RAS show that the high instability in hydrogen bonding interactions (HBIs) of GTP with residue 41 and Y42 in the switch domain I drives the disordered states of the switch domains. This work is expected to provide a molecular mechanism for deeply understanding the function of M-RAS and future drug design towards the treatment of cancers.

| S-EPMC10114971 | biostudies-literature

Probing inhibition mechanisms of adenosine deaminase by using molecular dynamics simulations.

Project description:Adenosine deaminase (ADA) catalyzes the deamination of adenosine, which is important in purine metabolism. ADA is ubiquitous to almost all human tissues, and ADA abnormalities have been reported in various diseases, including rheumatoid arthritis. ADA can be divided into two conformations based on the inhibitor that it binds to: open and closed forms. Here, we chose three ligands, namely, FR117016 (FR0), FR221647 (FR2) (open form), and HDPR (PRH, closed form), to investigate the inhibition mechanism of ADA and its effect on ADA through molecular dynamics simulations. In open forms, Egap and electrostatic potential (ESP) indicated that electron transfer might occur more easily in FR0 than in FR2. Binding free energy and hydrogen bond occupation revealed that the ADA-FR0 complex had a more stable structure than ADA-FR2. The probability of residues Pro159 to Lys171 of ADA-FR0 and ADA-FR2 to form a helix moderately increased compared with that in nonligated ADA. In comparison with FR0 and FR2 PRH could maintain ADA in a closed form to inhibit the function of ADA. The ?7 helix (residues Thr57 to Ala73) of ADA in the closed form was mostly unfastened because of the effect of PRH. The number of H bonds and the relative superiority of the binding free energy indicated that the binding strength of PRH to ADA was significantly lower than that of an open inhibitor, thereby supporting the comparison of the inhibitory activities of the three ligands. Alanine scanning results showed that His17, Gly184, Asp295, and Asp296 exerted the greatest effects on protein energy, suggesting that they played crucial roles in binding to inhibitors. This study served as a theoretical basis for the development of new ADA inhibitors.

| S-EPMC6239307 | biostudies-literature

Probing the druggability of membrane-bound Rab5 by molecular dynamics simulations.

Project description:Rab5 is a small GTPase and a key regulator in early endosomal trafficking. Rab5 and its effectors are involved in a large number of infectious diseases and certain types of cancer. We performed µs atomistic molecular dynamics simulations of inactive and active full-length Rab5 anchored to a complex model bilayer with composition of the early endosome membrane. Direct interactions between the Rab5?G domain and the bilayer were observed. We found two dominant nucleotide-dependent orientations characterised by a different accessibility of the switch regions. The "buried switch" orientation was mainly associated with inactive Rab5 accompanied with a rather extended structure of the hypervariable C-terminal region. Active Rab5 preferred an orientation in which the switch regions are accessible to effector proteins. These structural differences may provide an opportunity to selectively target one Rab5 state and lead to new approaches in the development of Rab5-specific therapies.

| S-EPMC6010109 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data