On the Use of Molecular Dynamics Simulations for Probing Allostery through DNA.
Ontology highlight
ABSTRACT:
SUBMITTER: Drsata T
PROVIDER: S-EPMC4776040 | biostudies-literature | 2016 Feb
REPOSITORIES: biostudies-literature
Similar Datasets
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
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
Project description:Many biological processes are regulated by allosteric mechanisms that communicate with distant sites in the protein responsible for functionality. The binding of a small molecule at an allosteric site typically induces conformational changes that propagate through the protein along allosteric pathways regulating enzymatic activity. Elucidating those communication pathways from allosteric sites to orthosteric sites is, therefore, essential to gain insights into biochemical processes. Targeting the allosteric pathways by mutagenesis can allow the engineering of proteins with desired functions. Furthermore, binding small molecule modulators along the allosteric pathways is a viable approach to target reactions using allosteric inhibitors/activators with temporal and spatial selectivity. Methods based on network theory can elucidate protein communication networks through the analysis of pairwise correlations observed in molecular dynamics (MD) simulations using molecular descriptors that serve as proxies for allosteric information. Typically, single atomic descriptors such as α-carbon displacements are used as proxies for allosteric information. Therefore, allosteric networks are based on correlations revealed by that descriptor. Here, we introduce a Python software package that provides a comprehensive toolkit for studying allostery from MD simulations of biochemical systems. MDiGest offers the ability to describe protein dynamics by combining different approaches, such as correlations of atomic displacements or dihedral angles, as well as a novel approach based on the correlation of Kabsch-Sander electrostatic couplings. MDiGest allows for comparisons of networks and community structures that capture physical information relevant to allostery. Multiple complementary tools for studying essential dynamics include principal component analysis, root mean square fluctuation, as well as secondary structure-based analyses.
| S-EPMC10769569 | biostudies-literature
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
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
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
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
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
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
Project description:The anti-Pneumocystis carinii activity of 13 synthetic pentamidine analogs was analyzed. The experimental differences in melting points of DNA dodecamer 5'-(CGCGAATTCGCG)2-3' complexes (ΔTm), and in the biological activity measured using ATP bioluminescence assay (IC50) together with the theoretical free energy of DNA-ligand binding estimated by the proposed computational protocol, showed that the experimental activity of the tested pentamidines appeared to be due to the binding to the DNA minor groove with extended AT sequences. The effect of heteroatoms in the aliphatic linker, and the sulfonamide or methoxy substituents on the compound inducing changes in the interactions with the DNA minor groove was examined and was correlated with biological activity. In computational analysis, the explicit solvent approximation with the discrete water molecules was taken into account, and the role of water molecules in the DNA-ligand complexes was defined.
| S-EPMC6272165 | biostudies-literature