Browse
Submit Data
Databases
API
Help

Dataset Information

22 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca²⁺-Calmodulin.

ABSTRACT: Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM's affinity for the calcium ions (Ca²⁺), or vice versa, may provide insight into how Ca²⁺-CaM selects its target binding proteins. However, modeling of Ca²⁺-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca²⁺ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca²⁺. This feature that underlies the reciprocal relation between Ca²⁺ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca²⁺ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a "community model" of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca²⁺ binding in future studies.

SUBMITTER: Nde J

PROVIDER: S-EPMC8421859 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Coarse grained molecular dynamics simulations of transmembrane protein-lipid systems.

Project description:Many biological cellular processes occur at the micro- or millisecond time scale. With traditional all-atom molecular modeling techniques it is difficult to investigate the dynamics of long time scales or large systems, such as protein aggregation or activation. Coarse graining (CG) can be used to reduce the number of degrees of freedom in such a system, and reduce the computational complexity. In this paper the first version of a coarse grained model for transmembrane proteins is presented. This model differs from other coarse grained protein models due to the introduction of a novel angle potential as well as a hydrogen bonding potential. These new potentials are used to stabilize the backbone. The model has been validated by investigating the adaptation of the hydrophobic mismatch induced by the insertion of WALP-peptides into a lipid membrane, showing that the first step in the adaptation is an increase in the membrane thickness, followed by a tilting of the peptide.

| S-EPMC2904924 | 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

Exploring TRPC3 Interaction with Cholesterol through Coarse-Grained Molecular Dynamics Simulations.

Project description:Transient receptor potential canonical 3 (TRPC3) channel belongs to the superfamily of transient receptor potential (TRP) channels which mediate Ca2+ influx into the cell. These channels constitute essential elements of cellular signalling and have been implicated in a wide range of diseases. TRPC3 is primarily gated by lipids and its surface expression has been shown to be dependent on cholesterol, yet a comprehensive exploration of its interaction with this lipid has thus far not emerged. Here, through 80 µs of coarse-grained molecular dynamics simulations, we show that cholesterol interacts with multiple elements of the transmembrane machinery of TRPC3. Through our approach, we identify an annular binding site for cholesterol on the pre-S1 helix and a non-annular site at the interface between the voltage-sensor-like domain and pore domains. Here, cholesterol interacts with exposed polar residues and possibly acts to stabilise the domain interface.

| S-EPMC9313397 | biostudies-literature

Capturing RNA Folding Free Energy with Coarse-Grained Molecular Dynamics Simulations.

Project description:We introduce a coarse-grained RNA model for molecular dynamics simulations, RACER (RnA CoarsE-gRained). RACER achieves accurate native structure prediction for a number of RNAs (average RMSD of 2.93 Å) and the sequence-specific variation of free energy is in excellent agreement with experimentally measured stabilities (R2 = 0.93). Using RACER, we identified hydrogen-bonding (or base pairing), base stacking, and electrostatic interactions as essential driving forces for RNA folding. Also, we found that separating pairing vs. stacking interactions allowed RACER to distinguish folded vs. unfolded states. In RACER, base pairing and stacking interactions each provide an approximate stability of 3-4 kcal/mol for an A-form helix. RACER was developed based on PDB structural statistics and experimental thermodynamic data. In contrast with previous work, RACER implements a novel effective vdW potential energy function, which led us to re-parameterize hydrogen bond and electrostatic potential energy functions. Further, RACER is validated and optimized using a simulated annealing protocol to generate potential energy vs. RMSD landscapes. Finally, RACER is tested using extensive equilibrium pulling simulations (0.86 ms total) on eleven RNA sequences (hairpins and duplexes).

| S-EPMC5385882 | biostudies-literature

Coarse-Grained Structural Modeling of Molecular Motors Using Multibody Dynamics.

Project description:Experimental and computational approaches are needed to uncover the mechanisms by which molecular motors convert chemical energy into mechanical work. In this article, we describe methods and software to generate structurally realistic models of molecular motor conformations compatible with experimental data from different sources. Coarse-grained models of molecular structures are constructed by combining groups of atoms into a system of rigid bodies connected by joints. Contacts between rigid bodies enforce excluded volume constraints, and spring potentials model system elasticity. This simplified representation allows the conformations of complex molecular motors to be simulated interactively, providing a tool for hypothesis building and quantitative comparisons between models and experiments. In an example calculation, we have used the software to construct atomically detailed models of the myosin V molecular motor bound to its actin track. The software is available at www.simtk.org.

| S-EPMC2860290 | biostudies-literature

Insertion of the Ca²?-independent phospholipase A? into a phospholipid bilayer via coarse-grained and atomistic molecular dynamics simulations.

Project description:Group VI Ca²?-independent phospholipase A? (iPLA?) is a water-soluble enzyme that is active when associated with phospholipid membranes. Despite its clear pharmaceutical relevance, no X-ray or NMR structural information is currently available for the iPLA? or its membrane complex. In this paper, we combine homology modeling with coarse-grained (CG) and all-atom (AA) molecular dynamics (MD) simulations to build structural models of iPLA? in association with a phospholipid bilayer. CG-MD simulations of the membrane insertion process were employed to provide a starting point for an atomistic description. Six AA-MD simulations were then conducted for 60 ns, starting from different initial CG structures, to refine the membrane complex. The resulting structures are shown to be consistent with each other and with deuterium exchange mass spectrometry (DXMS) experiments, suggesting that our approach is suitable for the modeling of iPLA? at the membrane surface. The models show that an anchoring region (residues 710-724) forms an amphipathic helix that is stabilized by the membrane. In future studies, the proposed iPLA? models should provide a structural basis for understanding the mechanisms of lipid extraction and drug-inhibition. In addition, the dual-resolution approach discussed here should provide the means for the future exploration of the impact of lipid diversity and sequence mutations on the activity of iPLA? and related enzymes.

| S-EPMC3723492 | biostudies-literature

Peptide nanopores and lipid bilayers: interactions by coarse-grained molecular-dynamics simulations.

Project description:A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (alpha-helical and beta-barrel), and the seven different bilayer systems range in thickness from approximately 28 to approximately 43 A. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.

| S-EPMC2711413 | biostudies-literature

Coarse-grained molecular dynamics simulations reveal lipid access pathways in P-glycoprotein.

Project description:P-glycoprotein (P-gp) exports a broad range of dissimilar compounds, including drugs, lipids, and lipid-like molecules. Because of its substrate promiscuity, P-gp is a key player in the development of cancer multidrug resistance. Although P-gp is one of the most studied ABC transporters, the mechanism by which its substrates access the cavity remains unclear. In this study, we perform coarse-grained molecular dynamics simulations to explore possible lipid access pathways in the inward-facing conformation of P-gp embedded in bilayers of different lipid compositions. In the inward-facing orientation, only lipids from the lower leaflet access the cavity of the transporter. We identify positively charged residues at the portals of P-gp that favor lipid entrance to the cavity, as well as lipid-binding sites at the portals and within the cavity, which is in good agreement with previous experimental studies. This work includes several examples of lipid pathways for phosphatidylcholine and phosphatidylethanolamine lipids that help elucidate the molecular mechanism of lipid binding in P-gp.

| S-EPMC5839720 | biostudies-literature

COARSE-GRAINED MODELING OF PROTEIN UNFOLDING DYNAMICS.

Project description:We present a new dynamic elastic network model (DENM) that describes the unfolding process of a force-loaded protein. The protein interaction network and its potentials are constructed based on information of its native-state structure obtained from the Protein Data Bank, with network nodes positioned at the C? coordinates of the protein backbone. Specifically, to mimic the unfolding process, i.e., to simulate the process of overcoming the local energy barrier on the free energy landscape with force loading, the noncovalent protein network bonds (i.e., hydrogen bonds, salt bridges, hydrophobic contacts, etc.) are broken one-by-one with a certain probability, while the strong covalent bonds along the backbone (i.e., peptide bonds, disulfide bonds, etc.) are kept intact. The jumping event from local energy minima (bonds breaking rate) are chosen according to Kramer's theory and the Bell model. Moreover, we exploit the self-similar structure of proteins at different scales to design an effective coarse-graining procedure for DENM with optimal parameter selection. The robustness of DENM is validated by coarse-grained molecular dynamics (MD) simulation against atomistic MD simulation of force-extension processes of the Fibrinogen and Titin Immunoglobulin proteins. We observe that the native structure of the proteins determines the total unfolding dynamics (including large deviations) and not just the fluctuations around the native state.

| S-EPMC4230303 | biostudies-literature

Disorder Mediated Oligomerization of DISC1 Proteins Revealed by Coarse-Grained Molecular Dynamics Simulations.

Project description:Disrupted-in-schizophrenia-1 (DISC1) is a scaffold protein of significant importance for neuro-development and a prominent candidate protein in the etiology of mental disorders. In this work, we investigate the role of conformational heterogeneity and local structural disorder in the oligomerization pathway of the full-length DISC1 and of two truncation variants. Through extensive coarse-grained molecular dynamics simulations with a predictive energy landscape-based model, we shed light on the interplay of local and global disorder which lead to different oligomerization pathways. We found that both global conformational heterogeneity and local structural disorder play an important role in shaping distinct oligomerization trends of DISC1. This study also sheds light on the differences in oligomerization pathways of the full-length protein compared to the truncated variants produced by a chromosomal translocation associated with schizophrenia. We report that oligomerization of full-length DISC1 sequence works in a nonadditive manner with respect to truncated fragments that do not mirror the conformational landscape or binding affinities of the full-length unit.

| S-EPMC7255695 | biostudies-literature