Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Molecular Dynamics Simulations of Molecules in Uniform Flow.

ABSTRACT: Flow at the molecular level induces shear-induced unfolding of single proteins and can drive their assembly, the mechanisms of which are not completely understood. To be able to analyze the role of flow on molecules, we present uniform-flow molecular dynamics simulations at atomic level. The pull module of the GRoningen MAchine for Chemical Simulations package was extended to be able to force-group atoms within a defined layer of the simulation box. Application of this external enforcement to explicit water molecules, together with the coupling to a thermostat, led to a uniform terminal velocity of the solvent water molecules. We monitored the density of the whole system to establish the conditions under which the simulated flow is well-behaved. A maximal velocity of 1.3 m/s can be generated if a pull slice of 8 nm is used, and high velocities would require larger pull slices to still maintain a stable density. As expected, the target velocity increases linearly with the total external force applied. Finally, we suggest an appropriate setup to stretch a protein by uniform flow, in which protein extensions depend on the flow conditions. Our implementation provides an efficient computational tool to investigate the effect of the flow at the molecular level.

SUBMITTER: Herrera-Rodriguez AM

PROVIDER: S-EPMC6506623 | biostudies-literature | 2019 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Exploring Structures and Dynamics of Protamine Molecules through Molecular Dynamics Simulations.

Project description:Protamines are arginine-rich proteins that condense DNA in sperm. Despite their importance in reproduction, information on protamine structure is scarce. We, therefore, used molecular dynamics to examine the structures of salmon, bull P1, and human P1 protamines. The sizes and shapes of each protamine varied widely, indicating that they were disordered with structures covering a broad conformational landscape, from hairpin loop structures to extended coils. Despite their general disorder, the protamines did form secondary structures, including helices and hairpin loops. In eutherians, hairpins may promote disulfide bonding that facilitates protamine-DNA condensation, but the specifics of this bonding is not well established. We examined inter-residue distances in the simulations to predict residue pairs likely to form intramolecular bonds, leading to the identification of bonding pairs consistent with previous results in bull and human. These results support a model for eutherian protamine structures where a highly charged center is surrounded by disulfide-bond-stabilized loops.

| S-EPMC9685783 | biostudies-literature

Molecular Dynamics Simulations of Kir2.2 Interactions with an Ensemble of Cholesterol Molecules.

Project description:Cholesterol is a major regulator of multiple types of ion channels, but the specific mechanisms and the dynamics of its interactions with the channels are not well understood. Kir2 channels were shown to be sensitive to cholesterol through direct interactions with "cholesterol-sensitive" regions on the channel protein. In this work, we used Martini coarse-grained simulations to analyze the long (μs) timescale dynamics of cholesterol with Kir2.2 channels embedded into a model membrane containing POPC phospholipid with 30 mol% cholesterol. This approach allows us to simulate the dynamic, unbiased migration of cholesterol molecules from the lipid membrane environment to the protein surface of Kir2.2 and explore the favorability of cholesterol interactions at both surface sites and recessed pockets of the channel. We found that the cholesterol environment surrounding Kir channels forms a complex milieu of different short- and long-term interactions, with multiple cholesterol molecules concurrently interacting with the channel. Furthermore, utilizing principles from network theory, we identified four discrete cholesterol-binding sites within the previously identified cholesterol-sensitive region that exist depending on the conformational state of the channel-open or closed. We also discovered that a twofold decrease in the cholesterol level of the membrane, which we found earlier to increase Kir2 activity, results in a site-specific decrease of cholesterol occupancy at these sites in both the open and closed states: cholesterol molecules at the deepest of these discrete sites shows no change in occupancy at different cholesterol levels, whereas the remaining sites showed a marked decrease in occupancy.

| S-EPMC6170799 | biostudies-literature

On achieving experimental accuracy from molecular dynamics simulations of flexible molecules: aqueous glycerol.

Project description:The rotational isomeric states (RIS) of glycerol at infinite dilution have been characterized in the aqueous phase via a 1 micros conventional molecular dynamics (MD) simulation, a 40 ns enhanced sampling replica exchange molecular dynamics (REMD) simulation, and a reevaluation of the experimental NMR data. The MD and REMD simulations employed the GLYCAM06/AMBER force field with explicit treatment of solvation. The shorter time scale of the REMD sampling method gave rise to RIS and theoretical scalar 3J(HH) coupling constants that were comparable to those from the much longer traditional MD simulation. The 3J(HH) coupling constants computed from the MD methods were in excellent agreement with those observed experimentally. Despite the agreement between the computed and the experimental J-values, there were variations between the rotamer populations computed directly from the MD data and those derived from the experimental NMR data. The experimentally derived populations were determined utilizing limiting J-values from an analysis of NMR data from substituted ethane molecules and may not be completely appropriate for application in more complex molecules, such as glycerol. Here, new limiting J-values have been derived via a combined MD and quantum mechanical approach and were used to decompose the experimental 3J(HH) coupling constants into population distributions for the glycerol RIS.

| S-EPMC4201037 | biostudies-literature

Multi-scale molecular dynamics simulations of enhanced energy transfer in organic molecules under strong coupling.

Project description:Exciton transport can be enhanced in the strong coupling regime where excitons hybridize with confined light modes to form polaritons. Because polaritons have group velocity, their propagation should be ballistic and long-ranged. However, experiments indicate that organic polaritons propagate in a diffusive manner and more slowly than their group velocity. Here, we resolve this controversy by means of molecular dynamics simulations of Rhodamine molecules in a Fabry-Pérot cavity. Our results suggest that polariton propagation is limited by the cavity lifetime and appears diffusive due to reversible population transfers between polaritonic states that propagate ballistically at their group velocity, and dark states that are stationary. Furthermore, because long-lived dark states transiently trap the excitation, propagation is observed on timescales beyond the intrinsic polariton lifetime. These insights not only help to better understand and interpret experimental observations, but also pave the way towards rational design of molecule-cavity systems for coherent exciton transport.

| S-EPMC10587084 | biostudies-literature

Molecular dynamics simulations in photosynthesis.

Project description:Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular dynamics (MD) simulations provide a powerful toolkit to investigate dynamical processes in (bio)molecular ensembles from the (sub)picosecond to the (sub)millisecond regime and from the Å to hundreds of nm length scale. Therefore, MD is well suited to address a variety of questions arising in the field of photosynthesis research. In this review, we provide an introduction to the basic concepts of MD simulations, at atomistic and coarse-grained level of resolution. Furthermore, we discuss applications of MD simulations to model photosynthetic systems of different sizes and complexity and their connection to experimental observables. Finally, we provide a brief glance on which methods provide opportunities to capture phenomena beyond the applicability of classical MD.

| S-EPMC7203591 | biostudies-literature

Flexible Fitting of Small Molecules into Electron Microscopy Maps Using Molecular Dynamics Simulations with Neural Network Potentials.

Project description:Despite significant advances in resolution, the potential for cryo-electron microscopy (EM) to be used in determining the structures of protein-drug complexes remains unrealized. Determination of accurate structures and coordination of bound ligands necessitates simultaneous fitting of the models into the density envelopes, exhaustive sampling of the ligand geometries, and, most importantly, concomitant rearrangements in the side chains to optimize the binding energy changes. In this article, we present a flexible-fitting pipeline where molecular dynamics flexible fitting (MDFF) is used to refine structures of protein-ligand complexes from 3 to 5 Å electron density data. Enhanced sampling is employed to explore the binding pocket rearrangements. To provide a model that can accurately describe the conformational dynamics of the chemically diverse set of small-molecule drugs inside MDFF, we use QM/MM and neural-network potential (NNP)/MM models of protein-ligand complexes, where the ligand is represented using the QM or NNP model, and the protein is represented using established molecular mechanical force fields (e.g., CHARMM). This pipeline offers structures commensurate to or better than recently submitted high-resolution cryo-EM or X-ray models, even when given medium to low-resolution data as input. The use of the NNPs makes the algorithm more robust to the choice of search models, offering a radius of convergence of 6.5 Å for ligand structure determination. The quality of the predicted structures was also judged by density functional theory calculations of ligand strain energy. This strain potential energy is found to systematically decrease with better fitting to density and improved ligand coordination, indicating correct binding interactions. A computationally inexpensive protocol for computing strain energy is reported as part of the model analysis protocol that monitors both the ligand fit as well as model quality.

| S-EPMC7311632 | biostudies-literature

Molecular dynamics simulations and drug discovery.

Project description:This review discusses the many roles atomistic computer simulations of macromolecular (for example, protein) receptors and their associated small-molecule ligands can play in drug discovery, including the identification of cryptic or allosteric binding sites, the enhancement of traditional virtual-screening methodologies, and the direct prediction of small-molecule binding energies. The limitations of current simulation methodologies, including the high computational costs and approximations of molecular forces required, are also discussed. With constant improvements in both computer power and algorithm design, the future of computer-aided drug design is promising; molecular dynamics simulations are likely to play an increasingly important role.

| S-EPMC3203851 | biostudies-literature

Molecular dynamics simulations: advances and applications.

Project description:Molecular dynamics simulations have evolved into a mature technique that can be used effectively to understand macromolecular structure-to-function relationships. Present simulation times are close to biologically relevant ones. Information gathered about the dynamic properties of macromolecules is rich enough to shift the usual paradigm of structural bioinformatics from studying single structures to analyze conformational ensembles. Here, we describe the foundations of molecular dynamics and the improvements made in the direction of getting such ensemble. Specific application of the technique to three main issues (allosteric regulation, docking, and structure refinement) is discussed.

| S-EPMC4655909 | biostudies-literature

Lipid Configurations from Molecular Dynamics Simulations.

Project description:The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na+ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the 31P chemical shift anisotropy and 14N and CD3 NMR quadrupole splittings.

| S-EPMC5937052 | biostudies-literature

Molecular dynamics simulations of biotin carboxylase.

Project description:Biotin carboxylase catalyzes the ATP-dependent carboxylation of biotin and is one component of the multienzyme complex acetyl-CoA carboxylase that catalyzes the first committed step in fatty acid synthesis in all organisms. Biotin carboxylase from Escherichia coli, whose crystal structures with and without ATP bound have been determined, has served as a model system for this component of the acetyl-CoA carboxylase complex. The two crystal structures revealed a large conformational change of one domain relative to the other domains when ATP is bound. Unfortunately, the crystal structure with ATP bound was obtained with an inactive site-directed mutant of the enzyme. As a consequence the structure with ATP bound lacked key structural information such as for the Mg2+ ions and contained altered conformations of key active-site residues. Therefore, nanosecond molecular dynamics studies of the wild-type biotin carboxylase were undertaken to supplant and amend the results of the crystal structures. Specifically, the protein-metal interactions of the two catalytically critical Mg2+ ions bound in the active site are presented along with a reevaluation of the conformations of active-site residues bound to ATP. In addition, the regions of the polypeptide chain that serve as hinges for the large conformational change were identified. The results of the hinge analysis complemented a covariance analysis that identified the individual structural elements of biotin carboxylase that change their conformation in response to ATP binding.

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