Concurrent Coupling of Electronic-Density-Functional, Molecular Dynamics, and Coarse-Grained Particles Schemes for Multiscale Simulation of Nanostructured Materials

2005 ◽  
Vol 502 ◽  
pp. 33-38 ◽  
Author(s):  
Shuji Ogata ◽  
Takahiro Igarashi
Keyword(s):  
Molecular Dynamics ◽  
Nanostructured Materials ◽  
Density Functional ◽  
Dft Method ◽  
Md Simulations ◽  
Multiscale Simulation ◽  
Coarse Grained ◽  
Dynamic Simulations ◽  
Electronic Density ◽  
Wide Range

Feature sizes of useful electronic devices are becoming smaller and reaching nanometer ranges. There is increasing demand to perform dynamic simulations of such nano-devices with realistic sizes. To date, various kinds of simulation methods have been used for materials and devices including the density-functional theory (DFT) and the molecular dynamics (MD) for atomistic mechanics and the finite element method for continuum mechanics. We review recent progresses in our multiscale, hybrid simulation schemes that combine those methods. The coarse-grained particles (CG) method originally proposed by Rudd and Broughton [Phys. Rev. B58 (1998), p. R5893] has features suitable to such hybridization. We improve the CG method so that it is applicable to realistic nanostructured materials with large deformations. A novel hybridization scheme that couples the DFT method with the MD method is presented, which is applicable to virtually any selection of the DFT region in a wide range of materials. Hybrid DFT-MD simulations of the H2O reaction with nanostructured Si and alumina systems under stresses are performed, to demonstrate significant effects of stress on the chemical reaction.

scholarly journals Solvent-induced depletion interactions in multiparticle collision dynamic simulations

2019 ◽  
Vol 30 (10) ◽  
pp. 1941008 ◽  
Author(s):  
Martin Wagner ◽  
Marisol Ripoll
Keyword(s):  
Molecular Dynamics ◽  
Md Simulations ◽  
Coarse Grained ◽  
Simulation Methods ◽  
Mesoscale Simulation ◽  
Dynamic Simulations ◽  
Collision Dynamic ◽  
Versatile Tool ◽  
Depletion Interactions

Molecular-dynamics-coupled multiparticle collision dynamic (MPC-MD) simulations have emerged to be an efficient and versatile tool in the description of mesoscale colloidal dynamics. However, the compressibility of the coarse-grained fluid leads to this method being prone to spurious depletion interactions that may dominate the colloidal dynamics. In this paper, we review the existing methodology to deal with these interactions, establish and report depletion measurements, and present a method to avoid artificial depletion in mesoscale simulation methods.

scholarly journals Binding of Ca2+-Independent C2 Domains to Lipid Membranes: a Multi-Scale Molecular Dynamics Study

2020 ◽  
Author(s):  
Andreas Haahr Larsen ◽  
Mark S.P. Sansom
Keyword(s):  
Molecular Dynamics ◽  
Lipid Bilayers ◽  
Md Simulations ◽  
Multiscale Simulation ◽  
Coarse Grained ◽  
Type I ◽  
Type Ii ◽  
Binding Modes ◽  
C2 Domains ◽  
Anionic Lipids

AbstractC2 domains facilitate protein-lipid interaction in cellular recognition and signalling processes. They possess a β-sandwich structure, with either type I or type II topology. C2 domains can interact with anionic lipid bilayers in either a Ca2+-dependent or a Ca2+-independent manner. The mechanism of recognition of anionic lipids by Ca2+-independent C2 domains is incompletely understood. We have used molecular dynamics (MD) simulations to explore the membrane interactions of six Ca2+– independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse grained MD simulations these C2 domains bound to lipid bilayers, forming transient interactions with zwitterionic (phosphatidylcholine, PC) bilayers compared to long lived interactions with anionic bilayers also containing either phosphatidylserine (PS) or PS and phosphatidylinositol bisphosphate (PIP2). Type I C2 domains bound non-canonically via the front, back or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops (as is typically the case for Ca2+-dependent C2 domains). C2 domains interacted strongly (up to 120 kJ/mol) with membranes containing PIP2 causing the bound anionic lipids to clustered around the protein. The C2 domains bound less strongly to anionic membranes without PIP2 (<50 kJ/mol), and most weakly to neutral membranes (<33 kJ/mol). Productive binding modes were identified and further analysed in atomistic simulations. For PTEN and SHIP2, CG simulations were also performed of the intact enzymes (i.e. phosphatase domain plus C2 domain) with PIP2-contating bilayers and the roles of the two domains in membrane localization were compared. From a methodological perspective, these studies establish a multiscale simulation protocol for studying membrane binding/recognition proteins, capable of revealing binding modes alongside details of lipid binding affinity and specificity.

scholarly journals Ligand-induced structural changes analysis of ribose-binding protein as studied by molecular dynamics simulations

2021 ◽  
pp. 1-12
Author(s):  
Haiyan Li ◽  
Zanxia Cao ◽  
Guodong Hu ◽  
Liling Zhao ◽  
Chunling Wang ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Structural Changes ◽  
Binding Protein ◽  
Network Models ◽  
Md Simulations ◽  
Protein Domain ◽  
Opening Angle ◽  
Conformation Changes ◽  
Wide Range ◽  
Ribose Binding Protein

BACKGROUND: The ribose-binding protein (RBP) from Escherichia coli is one of the representative structures of periplasmic binding proteins. Binding of ribose at the cleft between two domains causes a conformational change corresponding to a closure of two domains around the ligand. The RBP has been crystallized in the open and closed conformations. OBJECTIVE: With the complex trajectory as a control, our goal was to study the conformation changes induced by the detachment of the ligand, and the results have been revealed from two computational tools, MD simulations and elastic network models. METHODS: Molecular dynamics (MD) simulations were performed to study the conformation changes of RBP starting from the open-apo, closed-holo and closed-apo conformations. RESULTS: The evolution of the domain opening angle θ clearly indicates large structural changes. The simulations indicate that the closed states in the absence of ribose are inclined to transition to the open states and that ribose-free RBP exists in a wide range of conformations. The first three dominant principal motions derived from the closed-apo trajectories, consisting of rotating, bending and twisting motions, account for the major rearrangement of the domains from the closed to the open conformation. CONCLUSIONS: The motions showed a strong one-to-one correspondence with the slowest modes from our previous study of RBP with the anisotropic network model (ANM). The results obtained for RBP contribute to the generalization of robustness for protein domain motion studies using either the ANM or PCA for trajectories obtained from MD.

scholarly journals Molecular Dynamics Scoring of Protein–Peptide Models Derived from Coarse-Grained Docking

Molecules ◽  
2021 ◽  
Vol 26 (11) ◽  
pp. 3293
Author(s):  
Mateusz Zalewski ◽  
Sebastian Kmiecik ◽  
Michał Koliński
Keyword(s):  
Molecular Dynamics ◽  
Geometry Optimization ◽  
Computational Prediction ◽  
Md Simulations ◽  
Coarse Grained ◽  
Ligand Interaction ◽  
Vast Number ◽  
Docking Simulations ◽  
Peptide Models ◽  
Peptide Complexes

One of the major challenges in the computational prediction of protein–peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein–peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor–ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein–peptide complexes show that the proposed MD-based scoring approach can be used to identify protein–peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.

scholarly journals Theoretical Insights into the Electron Capture Behavior of H2SO4···N2O Complex: A DFT and Molecular Dynamics Study

Molecules ◽  
2018 ◽  
Vol 23 (9) ◽  
pp. 2349 ◽  
Author(s):  
Wei-Hua Wang ◽  
Wen-Ling Feng ◽  
Wen-Liang Wang ◽  
Ping Li
Keyword(s):  
Molecular Dynamics ◽  
Electron Capture ◽  
Atmospheric Chemistry ◽  
Density Functional ◽  
Molecular Complexes ◽  
Dft Method ◽  
Dynamics Simulation ◽  
Oh Radical ◽  
Before And After ◽  
Reduced Density

Both sulfuric acid (H2SO4) and nitrous oxide (N2O) play a central role in the atmospheric chemistry in regulating the global environment and climate changes. In this study, the interaction behavior between H2SO4 and N2O before and after electron capture has been explored using the density functional theory (DFT) method as well as molecular dynamics simulation. The intermolecular interactions have been characterized by atoms in molecules (AIM), natural bond orbital (NBO), and reduced density gradient (RDG) analyses, respectively. It was found that H2SO4 and N2O can form two transient molecular complexes via intermolecular H-bonds within a certain timescale. However, two molecular complexes can be transformed into OH radical, N2, and HSO4− species upon electron capture, providing an alternative formation source of OH radical in the atmosphere. Expectedly, the present findings not only can provide new insights into the transformation behavior of H2SO4 and N2O, but also can enable us to better understand the potential role of the free electron in driving the proceeding of the relevant reactions in the atmosphere.

Investigation, using density function theory, of coverage of the kaolinite (001) surface during hydrogen adsorption

Clay Minerals ◽  
2018 ◽  
Vol 53 (3) ◽  
pp. 393-402 ◽  
Author(s):  
Jian Zhao ◽  
Wei Gao ◽  
Zhi-Gang Tao ◽  
Hong-Yun Guo ◽  
Man-Chao He
Keyword(s):  
Density Functional ◽  
Hydrogen Adsorption ◽  
Function Theory ◽  
Lattice Relaxation ◽  
Electronic Density ◽  
Functional Theory ◽  
Adsorption Sites ◽  
Coverage Dependence ◽  
Hydrogen Molecules ◽  
Wide Range

ABSTRACTKaolinite can be used for many applications, including the underground storage of gases. Density functional theory was employed to investigate the adsorption of hydrogen molecules on the kaolinite (001) surface. The coverage dependence of the adsorption sites and energetics was studied systematically for a wide range of coverage, Θ (from 1/16 to 1 monolayer). The three-fold hollow site is the most stable, followed by the bridge, top-z and top sites. The adsorption energy of H2 decreased with increasing coverage, thus indicating the lower stability of surface adsorption due to the repulsion of neighbouring H2 molecules. The coverage has obvious effects on hydrogen adsorption. Other properties of the H2/kaolinite (001) system, including the lattice relaxation and changes of electronic density of states, were also studied and are discussed in detail.

scholarly journals A hierarchical Bayesian framework for force field selection in molecular dynamics simulations

2016 ◽  
Vol 374 (2060) ◽  
pp. 20150032 ◽  
Author(s):  
S. Wu ◽  
P. Angelikopoulos ◽  
C. Papadimitriou ◽  
R. Moser ◽  
P. Koumoutsakos
Keyword(s):  
Experimental Data ◽  
Molecular Dynamics ◽  
Force Field ◽  
Computational Cost ◽  
Md Simulations ◽  
Bayesian Framework ◽  
Underlying Structure ◽  
Hierarchical Bayesian ◽  
Wide Range ◽  
Selection Of

We present a hierarchical Bayesian framework for the selection of force fields in molecular dynamics (MD) simulations. The framework associates the variability of the optimal parameters of the MD potentials under different environmental conditions with the corresponding variability in experimental data. The high computational cost associated with the hierarchical Bayesian framework is reduced by orders of magnitude through a parallelized Transitional Markov Chain Monte Carlo method combined with the Laplace Asymptotic Approximation. The suitability of the hierarchical approach is demonstrated by performing MD simulations with prescribed parameters to obtain data for transport coefficients under different conditions, which are then used to infer and evaluate the parameters of the MD model. We demonstrate the selection of MD models based on experimental data and verify that the hierarchical model can accurately quantify the uncertainty across experiments; improve the posterior probability density function estimation of the parameters, thus, improve predictions on future experiments; identify the most plausible force field to describe the underlying structure of a given dataset. The framework and associated software are applicable to a wide range of nanoscale simulations associated with experimental data with a hierarchical structure.

scholarly journals Antimicrobial action of the cationic peptide, chrysophsin-3: a coarse-grained molecular dynamics study

Soft Matter ◽  
2018 ◽  
Vol 14 (15) ◽  
pp. 2796-2807 ◽  
Author(s):  
Andrea Catte ◽  
Mark R. Wilson ◽  
Martin Walker ◽  
Vasily S. Oganesyan
Keyword(s):  
Molecular Dynamics ◽  
Large Scale ◽  
Md Simulations ◽  
Antimicrobial Action ◽  
Coarse Grained ◽  
Cationic Peptide ◽  
Coarse Grained Molecular Dynamics

Antimicrobial action of a cationic peptide is modelled by large scale MD simulations.

Sign in / Sign up

Export Citation Format

Share Document