Browse
Submit Data
Databases
API
Help

Dataset Information

0 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

Electronic Dynamics of a Molecular System Coupled to a Plasmonic Nanoparticle Combining the Polarizable Continuum Model and Many-Body Perturbation Theory.

Project description:The efficiency of plasmonic metallic nanoparticles in harvesting and concentrating light energy in their proximity triggers a wealth of important and intriguing phenomena. For example, spectroscopies are able to reach single-molecule and intramolecule sensitivities, and important chemical reactions can be effectively photocatalyzed. For the real-time description of the coupled dynamics of a molecule's electronic system and of a plasmonic nanoparticle, a methodology has been recently proposed (J. Phys. Chem. C. 120, 2016, 28774-28781) which combines the classical description of the nanoparticle as a polarizable continuum medium with a quantum-mechanical description of the molecule treated at the time-dependent configuration interaction (TDCI) level. In this work, we extend this methodology by describing the molecule using many-body perturbation theory: the molecule's excitation energies, transition dipoles, and potentials computed at the GW/Bethe-Salpeter equation (BSE) level. This allows us to overcome current limitations of TDCI in terms of achievable accuracy without compromising on the accessible molecular sizes. We illustrate the developed scheme by characterizing the coupled nanoparticle/molecule dynamics of two prototype molecules, LiCN and p-nitroaniline.

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

Nonaqueous Ion Pairing Exemplifies the Case for Including Electronic Polarization in Molecular Dynamics Simulations.

Project description:The inclusion of electronic polarization is of crucial importance in molecular simulations of systems containing charged moieties. When neglected, as often done in force field simulations, charge-charge interactions in solution may become severely overestimated, leading to unrealistically strong bindings of ions to biomolecules. The electronic continuum correction introduces electronic polarization in a mean-field way via scaling of charges by the reciprocal of the square root of the high-frequency dielectric constant of the solvent environment. Here, we use ab initio molecular dynamics simulations to quantify the effect of electronic polarization on pairs of like-charged ions in a model nonaqueous environment where electronic polarization is the only dielectric response. Our findings confirm the conceptual validity of this approach, underlining its applicability to complex aqueous biomolecular systems. Simultaneously, the results presented here justify the potential employment of weaker charge scaling factors in force field development.

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

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