Investigation of slow molecular dynamics using R-CODEX.
Ontology highlight
ABSTRACT:
SUBMITTER: Li W
PROVIDER: S-EPMC3572235 | biostudies-literature | 2012 Sep
REPOSITORIES: biostudies-literature
Publications
Investigation of slow molecular dynamics using R-CODEX.
Journal of magnetic resonance (San Diego, Calif. : 1997) 20120607
A solid state NMR experiment is introduced for probing motions on the millisecond time scale, based on dephasing and refocusing (1)H-(13)C or (1)H-(15)N dipolar couplings. The method is related to the previously described Centerband-Only Detection of Exchange or CODEX experiment. The use of an R-type dipolar recoupling sequence takes advantage of the strong (1)H-(13)C or (1)H-(15)N dipolar coupling, while suppressing the effect of (1)H-(1)H homonuclear coupling. This approach paves the way to de ...[more]
Similar Datasets
Project description:Individual DNA bases are known to be able to flip out of the helical stack, providing enzymes with access to the genetic information otherwise hidden inside the helix. Consequently, base flipping is a necessary first step to many more complex biological processes such as DNA transcription or replication. Much remains unknown about this elementary step, despite a wealth of experimental and theoretical studies. From the theoretical point of view, the involved timescale of milliseconds or longer requires the use of enhanced sampling techniques. In contrast to previous theoretical studies employing umbrella sampling along a predefined flipping coordinate, this study attempts to induce flipping without prior knowledge of the pathway, using information from a molecular dynamics simulation of a B-DNA fragment and the conformational flooding method. The relevance to base flipping of the principal components of the simulation is assayed, and a combination of modes optimally related to the flipping of the base through either helical groove is derived for each of the two bases of the central guanine-cytosine basepair. By applying an artificial flooding potential along these collective coordinates, the flipping mechanism is accelerated to within the scope of molecular dynamics simulations. The associated free energy surface is found to feature local minima corresponding to partially flipped states, particularly relevant to flipping in isolated DNA; further transitions from these minima to the fully flipped conformation are accelerated by additional flooding potentials. The associated free energy profiles feature similar barrier heights for both bases and pathways; the flipped state beyond is a broad and rugged attraction basin, only a few kcal/mol higher in energy than the closed conformation. This result diverges from previous works but echoes some aspects of recent experimental findings, justifying the need for novel approaches to this difficult problem: this contribution represents a first step in this direction. Important structural factors involved in flipping, both local (sugar-phosphate backbone dihedral angles) and global (helical axis bend), are also identified.
| S-EPMC1913169 | biostudies-literature
Project description:The synchronization of native state motions as they transition between microstates influences catalysis kinetics, mediates allosteric interactions, and reduces the conformational entropy of proteins. However, it has proven difficult to describe native microstates because they are usually minimally frustrated and may interconvert on the micro- to millisecond time scale. Direct observation of concerted equilibrium fluctuations would therefore be an important tool for describing protein native states. Here we propose a strategy that relates NMR cross-correlated relaxation (CCR) rates between dipolar interactions to residual dipolar couplings (RDCs) of individual consecutive H(N)-N and H(?)-C(?) bonds, which act as a proxy for the peptide planes and the side chains, respectively. Using Xplor-NIH ensemble structure calculations restrained with the RDC and CCR data, we observe collective motions on time scales slower than nanoseconds in the backbone for GB3. To directly access the correlations from CCR, we develop a structure-free data analysis. The resulting dynamic correlation map is consistent with the ensemble-restrained simulations and reveals a complex network. In general, we find that the bond motions are on average slightly correlated and that the local environment dominates many observations. Despite this, some patterns are typical over entire secondary structure elements. In the ?-sheet, nearly all bonds are weakly correlated, and there is an approximately binary alternation in correlation intensity corresponding to the solvent exposure/shielding alternation of the side chains. For ?-helices, there is also a weak correlation in the H(N)-N bonds. The degree of correlation involving H(?)-C(?) bonds is directly affected by side-chain fluctuations, whereas loops show complex and nonuniform behavior.
| S-EPMC5055379 | biostudies-literature
Project description:In this study, we investigate the atomistic details of Keap1-Nrf2 inhibitors by in-depth modeling techniques, including molecular dynamics (MD) simulations, and the path-based free energy method of umbrella sampling (US). The protein-protein interaction (PPI) of Keap1-Nrf2 is implicated in several neurodegenerative diseases like cancer, diabetes, and cardiomyopathy. A better understanding of the five sub-pocket binding sites for Nrf2 (ETGE and DLG motifs) inside the Kelch domain would expedite the inhibitor design process. We selected four protein-ligand complexes with distinct co-crystal ligands and binding occupancies inside the Nrf2 binding site. We performed 100 ns of MD simulation for each complex and analyzed the trajectories. From the results, it is evident that one ligand (1VV) has flipped inside the binding pocket, whereas the remaining three were stable. We found that Coulombic (Arg483, Arg415, Ser363, Ser508, and Ser602) and Lennard-Jones (Tyr525, Tyr334, and Tyr572) interactions played a significant role in complex stability. The obtained binding free energy values from US simulations were consistent with the potencies of simulated ligands. US simulation highlight the importance of basic and aromatic residues in the binding pocket. A detailed description of the dissociation process brings valuable insight into the interaction of the four selected protein-ligand complexes, which could help in the future to design more potent PPI inhibitors.
| S-EPMC6891428 | biostudies-literature
Project description:The molecular flexibility of an inhibitor in ligand-binding process has been investigated by the mass-weighted molecular-dynamics simulation, a computational method adopted from the standard molecular-dynamics simulation and one by which the conformational space of a biomolecular system over potential energy barriers can be sampled effectively. The bimolecular complex of the aspartyl proteinase from Rhizopus chinensis, rhizopuspepsin, and an octapeptide inhibitor was previously studied in a mass-weighted molecular-dynamics simulation; the study has been extended for investigating the molecular flexibility in ligand binding. A series of mass-weighted molecular-dynamics simulations was carried out in which libration of the inhibitor dihedral angles was parametrically controlled, and threshold values of dihedral angle libration amplitudes were observed from monitoring the sampling of the enzyme binding pocket by the inhibitor in the simulations. The computational results are consistent with the general notion of molecular-flexibility requirement for ligand binding; the freedom of dihedral rotations of side-chain groups was found to be particularly important for ligand binding. Thus the critical degree of molecular flexibility which would contribute to effective enzyme inhibition can be obtained precisely from the modified molecular-dynamics simulations; the procedure described herein represents a first step toward providing quantitative measures of such a molecular-flexibility index for inhibitor molecules that have been otherwise targeted for optimal protein-ligand interactions.
| S-EPMC1132086 | biostudies-other
Project description:Current studies of cell signaling dynamics that use live cell fluorescent biosensors routinely yield thousands of single-cell, heterogeneous, multi-dimensional trajectories. Typically, the extraction of relevant information from time series data relies on predefined, human-interpretable features. Without a priori knowledge of the system, the predefined features may fail to cover the entire spectrum of dynamics. Here we present CODEX, a data-driven approach based on convolutional neural networks (CNNs) that identifies patterns in time series. It does not require a priori information about the biological system and the insights into the data are built through explanations of the CNNs' predictions. CODEX provides several views of the data: visualization of all the single-cell trajectories in a low-dimensional space, identification of prototypic trajectories, and extraction of distinctive motifs. We demonstrate how CODEX can provide new insights into ERK and Akt signaling in response to various growth factors, and we recapitulate findings in p53 and TGFβ-SMAD2 signaling.
| S-EPMC8034356 | biostudies-literature
Project description:Vinculin activation plays a critical role in focal adhesion initiation and formation. In its native state, vinculin is in an autoinhibitory conformation in which domain 1 prevents interaction of the vinculin tail domain with actin by steric hindrance. Once activated, vinculin is able to interact with both actin and talin. Several hypotheses have been put forth addressing the mechanisms of vinculin activation. One set of studies suggests that vinculin interaction with talin is sufficient to cause activation, whereas another set of studies suggests that a simultaneous interaction with several binding partners is necessary to achieve vinculin activation. Using molecular-dynamics (MD) simulations, we investigate the mechanisms of vinculin activation and suggest both a trajectory of conformational changes leading to vinculin activation, and key structural features that are likely involved in stabilizing the autoinhibited conformation. Assuming that the simultaneous interaction of vinculin with both actin and talin causes a stretching force on vinculin, and that vinculin activation results from a removal of steric hindrance blocking the actin-binding sites, we simulate with MD the stretching and activation of vinculin. The MD simulations are further confirmed by normal-mode analysis and simulation after residue modification. Taken together, the results of these simulations suggest that bending of the vinculin-binding-site region in vinculin away from the vinculin tail is the likely trajectory of vinculin activation.
| S-EPMC2920635 | biostudies-literature
Project description:The protein kinase C? (PKC?) enzyme is a member of a broad family of serine/threonine kinases, which are involved in varied cellular signaling pathways. The initial step of PKC? activation involves the C2 subunit docking with the cell membrane, which is followed by interactions of the C1 domains with diacylglycerol (DAG) in the membrane. Notably, the molecular mechanisms of these interactions remain poorly understood, especially what effects, if any, DAG may have on the initial C2 docking. To further understand this process, we have performed a series of conventional molecular dynamics simulations to systematically investigate the interaction between PKC?-C2 domains and lipid bilayers with different compositions to examine the effects of POPS, PIP2, and 1-palmitoyl-2-oleoyl-sn-glycerol (POG) on domain docking. Our results show that the PKC?-C2 domain does not interact with the bilayer surface in the absence of POPS and PIP2. In contrast, the inclusion of POPS and PIP2 to the bilayer resulted in strong domain docking in both perpendicular and parallel orientations, whereas the further inclusion of POG resulted in only parallel domain docking. In addition, lysine residues in the C2 domain formed hydrogen bonds with PIP2 molecule bilayers containing POG. These effects were further explored with umbrella sampling calculations to estimate the free energy of domain docking to the lipid bilayer in the presence of one or two PIP2 molecules. The results show that the binding of one or two PIP2 molecules is thermodynamically favorable, although stronger in bilayers lacking POG. However, in POG-containing bilayers, the binding mode of the C2 domain appears to be more flexible, which may have implications for activation of full-length PKC?. Together, our results shed new insights into the process of C2 bilayer binding and suggest new mechanisms for the roles of different phospholipids in the activation process of PKC?.
| S-EPMC5582975 | biostudies-literature
Project description:Lactate dehydrogenase A (LDHA) is an important enzyme in fermentative glycolysis, generating most energy for cancer cells that rely on anaerobic respiration even under normal oxygen concentrations. This renders LDHA a promising molecular target for the treatment of various cancers. Several efforts have been made recently to develop LDHA inhibitors with nanomolar inhibition and cellular activity, some of which have been studied in complex with the enzyme by X-ray crystallography. In this work, we present a molecular dynamics (MD) study of the binding interactions of selected ligands with human LDHA. Conventional MD simulations demonstrate different binding dynamics of inhibitors with similar binding affinities, whereas steered MD simulations yield discrimination of selected LDHA inhibitors with qualitative correlation between the in silico unbinding difficulty and the experimental binding strength. Further, our results have been used to clarify ambiguities in the binding modes of two well-known LDHA inhibitors.
| S-EPMC3895040 | biostudies-literature
Project description:Recent experiments with undersaturated aqueous glycine solutions have repeatedly exhibited the presence of giant liquid-like clusters or nanodroplets around 100 nm in diameter. These nanodroplets re-appear even after careful efforts for their removal and purification of the glycine solution. The composition of these clusters is not clear, although it has been suggested that they are mainly composed of glycine, a small and very soluble amino acid. To gain insights into this phenomenon, we study the aggregation of glycine in aqueous solutions at concentrations below the experimental solubility limit using large-scale molecular dynamics simulations under ambient conditions. Three protonation states of glycine (zwitterion = GLZ, anion = GLA, and cation = GLC) are simulated using molecular force fields based on the 1.14*CM1A partial charge scheme, which incorporates the OPLS all-atom force field and TIP3P water. When initiated from dispersed states, we find that giant clusters do not form in our simulations unless salt impurities are present. Moreover, if simulations are initiated from giant cluster states, we find that they tend to dissolve in the absence of salt impurities. Therefore, the simulation results provide little support for the possibility that the giant clusters seen in experiments are composed purely of glycine (and water). Considering that strenuous efforts are made in experiments to remove impurities such as salt, we propose that the giant clusters observed might instead result from the aggregation of reaction products of aqueous glycine, such as diketopiperazine or other oligoglycines which may be difficult to separate from glycine using conventional methods, or their co-aggregation with glycine.
| S-EPMC9251761 | biostudies-literature
Project description:A molecular-level understanding of the function of a protein requires knowledge of both its structural and dynamic properties. NMR spectroscopy allows the measurement of generalized order parameters that provide an atomistic description of picosecond and nanosecond fluctuations in protein structure. Molecular dynamics (MD) simulation provides a complementary approach to the study of protein dynamics on similar time scales. Comparisons between NMR spectroscopy and MD simulations can be used to interpret experimental results and to improve the quality of simulation-related force fields and integration methods. However, apparent systematic discrepancies between order parameters extracted from simulations and experiments are common, particularly for elements of noncanonical secondary structure. In this paper, results from a 1.2 micros explicit solvent MD simulation of the protein ubiquitin are compared with previously determined backbone order parameters derived from NMR relaxation experiments [Tjandra, N.; Feller, S. E.; Pastor, R. W.; Bax, A. J. Am. Chem. Soc. 1995, 117, 12562-12566]. The simulation reveals fluctuations in three loop regions that occur on time scales comparable to or longer than that of the overall rotational diffusion of ubiquitin and whose effects would not be apparent in experimentally derived order parameters. A coupled analysis of internal and overall motion yields simulated order parameters substantially closer to the experimentally determined values than is the case for a conventional analysis of internal motion alone. Improved agreement between simulation and experiment also is encouraging from the viewpoint of assessing the accuracy of long MD simulations.
| S-EPMC2805408 | biostudies-literature