Browse
Submit Data
Databases
API
Help

Dataset Information

13 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Electronic continuum model for molecular dynamics simulations.

ABSTRACT: A simple model for accounting for electronic polarization in molecular dynamics (MD) simulations is discussed. In this model, called molecular dynamics electronic continuum (MDEC), the electronic polarization is treated explicitly in terms of the electronic continuum (EC) approximation, while the nuclear dynamics is described with a fixed-charge force field. In such a force-field all atomic charges are scaled to reflect the screening effect by the electronic continuum. The MDEC model is rather similar but not equivalent to the standard nonpolarizable force-fields; the differences are discussed. Of our particular interest is the calculation of the electrostatic part of solvation energy using standard nonpolarizable MD simulations. In a low-dielectric environment, such as protein, the standard MD approach produces qualitatively wrong results. The difficulty is in mistreatment of the electronic polarizability. We show how the results can be much improved using the MDEC approach. We also show how the dielectric constant of the medium obtained in a MD simulation with nonpolarizable force-field is related to the static (total) dielectric constant, which includes both the nuclear and electronic relaxation effects. Using the MDEC model, we discuss recent calculations of dielectric constants of alcohols and alkanes, and show that the MDEC results are comparable with those obtained with the polarizable Drude oscillator model. The applicability of the method to calculations of dielectric properties of proteins is discussed.

SUBMITTER: Leontyev IV

PROVIDER: S-EPMC3910273 | biostudies-other | 2009 Feb

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

GEM*: A Molecular Electronic Density-Based Force Field for Molecular Dynamics Simulations.

Project description:GEM*, a force field that combines Coulomb and Exchange terms calculated with Hermite Gaussians with the polarization, bonded, and modified van der Waals terms from AMOEBA is presented. GEM* is tested on an initial water model fitted at the same level as AMOEBA. The integrals required for the evaluation of the intermolecular Coulomb interactions are efficiently evaluated by means of reciprocal space methods. The GEM* water model is tested by comparing energies and forces for a series of water oligomers and MD simulations. Timings for GEM* compared to AMOEBA are presented and discussed.

| S-EPMC5207213 | biostudies-literature

Squeezing protein shells: how continuum elastic models, molecular dynamics simulations, and experiments coalesce at the nanoscale.

Project description:The current rapid growth in the use of nanosized particles is fueled in part by our increased understanding of their physical properties and ability to manipulate them, which is essential for achieving optimal functionality. Here we report detailed quantitative measurements of the mechanical response of nanosized protein shells (viral capsids) to large-scale physical deformations and compare them with theoretical descriptions from continuum elastic modeling and molecular dynamics (MD). Specifically, we used nanoindentation by atomic force microscopy to investigate the complex elastic behavior of Hepatitis B virus capsids. These capsids are hollow, approximately 30 nm in diameter, and conform to icosahedral (5-3-2) symmetry. First we show that their indentation behavior, which is symmetry-axis-dependent, cannot be reproduced by a simple model based on Föppl-von Kármán thin-shell elasticity with the fivefold vertices acting as prestressed disclinations. However, we can properly describe the measured nonlinear elastic and orientation-dependent force response with a three-dimensional, topographically detailed, finite-element model. Next, we show that coarse-grained MD simulations also yield good agreement with our nanoindentation measurements, even without any fitting of force-field parameters in the MD model. This study demonstrates that the material properties of viral nanoparticles can be correctly described by both modeling approaches. At the same time, we show that even for large deformations, it suffices to approximate the mechanical behavior of nanosized viral shells with a continuum approach, and ignore specific molecular interactions. This experimental validation of continuum elastic theory provides an example of a situation in which rules of macroscopic physics can apply to nanoscale molecular assemblies.

| S-EPMC2920642 | biostudies-literature

Molecular dynamics simulations of indolicidin association with model lipid bilayers.

Project description:Identifying the mechanisms responsible for the interaction of peptides with cell membranes is critical to the design of new antimicrobial peptides and membrane transporters. We report here the results of a computational simulation of the interaction of the 13-residue peptide indolicidin with single-phase lipid bilayers of dioleoylphosphatidylcholine, distearoylphosphatidylcholine, dioleoylphosphatidylglycerol, and distearoylphosphatidylglycerol. Ensemble analysis of the membrane-bound peptide revealed that, in contrast to the extended, linear backbone structure reported for indolicidin in sodium dodecyl sulphate detergent micelles, the peptide adopts a boat-shaped conformation in both phosphatidylglycerol and phosphatidylcholine lipid bilayers, similar to that reported for dodecylphosphocholine micelles. In agreement with fluorescence and NMR experiments, simulations confirmed that the peptide localizes in the membrane interface, with the distance between phosphate headgroups of each leaflet being reduced in the presence of indolicidin. These data, along with a concomitant decrease in lipid order parameters for the upper-tail region, suggest that indolicidin binding results in membrane thinning, consistent with recent in situ atomic force microscopy studies.

| S-EPMC1877785 | biostudies-literature

Entropy Rules: Molecular Dynamics Simulations of Model Oligomers for Thermoresponsive Polymers.

Project description:We attempted to attain atomic-scale insights into the mechanism of the heat-induced phase transition of two thermoresponsive polymers containing amide groups, poly(N-isopropylacrylamide) (PNIPAM) and poly(2-isopropyl-2-oxazoline) (PIPOZ), and we succeeded in reproducing the existence of lower critical solution temperature (LCST). The simulation data are in accord with experimental findings. We found out that the entropy has an important contribution to the thermodynamics of the phase separation transition. Moreover, after decomposing further the entropy change to contributions from the solutes and from the solvent, it appeared out that the entropy of the solvent has the decisive share for the lowering of the free energy of the system when increasing the temperature above the LCST. Our conclusion is that the thermoresponsive behavior is driven by the entropy of the solvent. The water molecules structured around the functional groups of the polymer that are exposed to contact with the solvent in the extended conformation lower the enthalpy of the system, but at certain temperature the extended conformation of the polymer collapses as a result of dominating entropy gain from "released" water molecules. We stress also on the importance of using more than one reference molecule in the simulation box at the setup of the simulation.

| S-EPMC7597358 | biostudies-literature

Polarizable molecular dynamics in a polarizable continuum solvent.

Project description:We present, for the first time, scalable polarizable molecular dynamics (MD) simulations within a polarizable continuum solvent with molecular shape cavities and exact solution of the mutual polarization. The key ingredients are a very efficient algorithm for solving the equations associated with the polarizable continuum, in particular, the domain decomposition Conductor-like Screening Model (ddCOSMO), which involves a rigorous coupling of the continuum with the polarizable force field achieved through a robust variational formulation and an effective strategy to solve the coupled equations. The coupling of ddCOSMO with nonvariational force fields, including AMOEBA, is also addressed. The MD simulations are feasible, for real-life systems, on standard cluster nodes; a scalable parallel implementation allows for further acceleration in the context of a newly developed module in Tinker, named Tinker-HP. NVE simulations are stable, and long-term energy conservation can be achieved. This paper is focused on the methodological developments, the analysis of the algorithm, and the stability of the simulations; a proof-of-concept application is also presented to attest to the possibilities of this newly developed technique.

| S-EPMC4621208 | biostudies-other

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

Conformational model of the Holliday junction transition deduced from molecular dynamics simulations.

Project description:Homologous recombination plays a key role in the restart of stalled replication forks and in the generation of genetic diversity. During this process, two homologous DNA molecules undergo strand exchange to form a four-way DNA (Holliday) junction. In the presence of metal ions, the Holliday junction folds into the stacked-X structure that has two alternative conformers. Experiments have revealed the spontaneous transitions between these conformers, but their detailed pathways are not known. Here, we report a series of molecular dynamics simulations of the Holliday junction at physiological and elevated (400 K) temperatures. The simulations reveal new tetrahedral intermediates and suggest a schematic framework for conformer transitions. The tetrahedral intermediates bear resemblance to the junction conformation in complex with a junction-resolving enzyme, T7 endonuclease I, and indeed, one intermediate forms a stable complex with the enzyme as demonstrated in one simulation. We also describe free energy minima for various states of the Holliday junction system, which arise during conformer transitions. The results show that magnesium ions stabilize the stacked-X form and destabilize the open and tetrahedral intermediates. Overall, our study provides a detailed dynamic model of the Holliday junction undergoing a conformer transition.

| S-EPMC545448 | biostudies-literature

Comparative molecular dynamics simulations of the antimicrobial peptide CM15 in model lipid bilayers.

Project description:We report altogether 3-?s molecular dynamics (MD) simulations of the antimicrobial peptide CM15 to systematically investigate its interaction with two model lipid bilayers, pure POPC and mixed POPG:POPC (1:2). Starting with either an ?-helical or a random-coil conformation, CM15 is found to insert into both bilayers. Peptide-lipid interaction is stronger with the anionic POPG:POPC than the zwitterionic POPC, which is largely attributed to the electrostatic attraction between CM15 and the negatively charged POPG. Simulations initiated with CM15 as a random coil allowed us to study peptide folding at the lipid-water interface. Interestingly, CM15 folding appears to be faster in POPC than POPG:POPC, which may be explained by a lower activation energy barrier of structural rearrangement in the former system. Our data also suggest that compared with the random-coil conformation, CM15 in a pre-folded ?-helix has significantly reduced interactions with the lipids, indicating that peptide initial structures may bias the simulation results considerably on the 100-ns timescale. The implications of this result should be considered when preparing and interpreting future AMP simulations.

| S-EPMC3329965 | biostudies-literature

A kinetic model of trp-cage folding from multiple biased molecular dynamics simulations.

Project description:Trp-cage is a designed 20-residue polypeptide that, in spite of its size, shares several features with larger globular proteins.Although the system has been intensively investigated experimentally and theoretically, its folding mechanism is not yet fully understood. Indeed, some experiments suggest a two-state behavior, while others point to the presence of intermediates. In this work we show that the results of a bias-exchange metadynamics simulation can be used for constructing a detailed thermodynamic and kinetic model of the system. The model, although constructed from a biased simulation, has a quality similar to those extracted from the analysis of long unbiased molecular dynamics trajectories. This is demonstrated by a careful benchmark of the approach on a smaller system, the solvated Ace-Ala3-Nme peptide. For theTrp-cage folding, the model predicts that the relaxation time of 3100 ns observed experimentally is due to the presence of a compact molten globule-like conformation. This state has an occupancy of only 3% at 300 K, but acts as a kinetic trap.Instead, non-compact structures relax to the folded state on the sub-microsecond timescale. The model also predicts the presence of a state at Calpha-RMSD of 4.4 A from the NMR structure in which the Trp strongly interacts with Pro12. This state can explain the abnormal temperature dependence of the Pro12-delta3 and Gly11-alpha3 chemical shifts. The structures of the two most stable misfolded intermediates are in agreement with NMR experiments on the unfolded protein. Our work shows that, using biased molecular dynamics trajectories, it is possible to construct a model describing in detail the Trp-cage folding kinetics and thermodynamics in agreement with experimental data.

| S-EPMC2711228 | biostudies-literature

Molecular dynamics simulations of a helicase.

Project description:Helicases are ubiquitous enzymes involved in nucleic acid metabolism. The PcrA DNA helicase is an essential bacterial protein involved in rolling circle plasmid replication and DNA repair. Recent crystal structures of PcrA bound to DNA indicate that a flexible loop mediates a functionally important rigid-body-domain rotation. In this study, we report stochastic boundary molecular dynamics simulations focused on this region for wild-type and mutants designed to increase the rigidity of the region. Residues in the region that were helix-disfavoring, such as glycine, threonine, and others, were mutated to alanine. The simulated dynamics, analyzed with a variety of measures of structure and mobility, indicate that a few point mutations will substantially increase helix formation in this region. Subnanosecond stochastic boundary molecular dynamics simulations at several temperatures offer a rapid protocol for assessing large numbers of mutants and provides a novel strategy for the design of experiments to test the role of this flexible loop region in the function of PcrA.

| S-EPMC3034217 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

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