Browse
Submit Data
Databases
API
Help

Dataset Information

10 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Evidence for percolation diffusion of cations and reordering in disordered pyrochlore from accelerated molecular dynamics.

ABSTRACT: Diffusion in complex oxides is critical to ionic transport, radiation damage evolution, sintering, and aging. In complex oxides such as pyrochlores, anionic diffusion is dramatically affected by cation disorder. However, little is known about how disorder influences cation transport. Here, we report results from classical and accelerated molecular dynamics simulations of vacancy-mediated cation diffusion in Gd₂Ti₂O₇ pyrochlore, on the microsecond timescale. We find that diffusion is slow at low levels of disorder, while higher disorder allows for fast diffusion, which is then accompanied by antisite annihilation and reordering, and thus a slowing of cation transport. Cation diffusivity is therefore not constant, but decreases as the material reorders. We also show that fast cation diffusion is triggered by the formation of a percolation network of antisites. This is in contrast with observations from other complex oxides and disordered media models, suggesting a fundamentally different relation between disorder and mass transport.Diffusion plays an important role in sintering, damage tolerance and transport. Here authors perform classical and accelerated molecular dynamics simulations of vacancy-mediated cation diffusion in Gd₂Ti₂O₇ pyrochlore and report non-monotonic evolution of cation diffusivity.

SUBMITTER: Perriot R

PROVIDER: S-EPMC5606988 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics.

Project description:Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ?1.5 ?s simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans-gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2-3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems.

| S-EPMC3191728 | biostudies-literature

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics.

Project description:An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038-2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pK(a) values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pK(a) values in closer agreement to those obtained experimentally for the more slowly converging residues.

| S-EPMC2817915 | biostudies-literature

Coulomb spin liquid in anion-disordered pyrochlore Tb2Hf2O7.

Project description:The charge ordered structure of ions and vacancies characterizing rare-earth pyrochlore oxides serves as a model for the study of geometrically frustrated magnetism. The organization of magnetic ions into networks of corner-sharing tetrahedra gives rise to highly correlated magnetic phases with strong fluctuations, including spin liquids and spin ices. It is an open question how these ground states governed by local rules are affected by disorder. Here we demonstrate in the pyrochlore Tb2Hf2O7, that the vicinity of the disordering transition towards a defective fluorite structure translates into a tunable density of anion Frenkel disorder while cations remain ordered. Quenched random crystal fields and disordered exchange interactions can therefore be introduced into otherwise perfect pyrochlore lattices of magnetic ions. We show that disorder can play a crucial role in preventing long-range magnetic order at low temperatures, and instead induces a strongly fluctuating Coulomb spin liquid with defect-induced frozen magnetic degrees of freedom.Experimental studies of frustrated spin systems such as pyrochlore magnetic oxides test our understanding of quantum many-body physics. Here the authors show experimentally that Tb2Hf2O7 may be a model material for investigating how structural disorder can stabilize a quantum spin liquid phase.

| S-EPMC5638941 | biostudies-other

Tumor proliferation and diffusion on percolation clusters.

Project description:We study in silico the influence of host tissue inhomogeneity on tumor cell proliferation and diffusion by simulating the mobility of a tumor on percolation clusters with different homogeneities of surrounding tissues. The proliferation and diffusion of a tumor in an inhomogeneous tissue could be characterized in the framework of the percolation theory, which displays similar thresholds (0.54, 0.44, and 0.37, respectively) for tumor proliferation and diffusion in three kinds of lattices with 4, 6, and 8 connecting near neighbors. Our study reveals the existence of a critical transition concerning the survival and diffusion of tumor cells with leaping metastatic diffusion movement in the host tissues. Tumor cells usually flow in the direction of greater pressure variation during their diffusing and infiltrating to a further location in the host tissue. Some specific sites suitable for tumor invasion were observed on the percolation cluster and around these specific sites a tumor can develop into scattered tumors linked by some advantage tunnels that facilitate tumor invasion. We also investigate the manner that tissue inhomogeneity surrounding a tumor may influence the velocity of tumor diffusion and invasion. Our simulation suggested that invasion of a tumor is controlled by the homogeneity of the tumor microenvironment, which is basically consistent with the experimental report by Riching et al. as well as our clinical observation of medical imaging. Both simulation and clinical observation proved that tumor diffusion and invasion into the surrounding host tissue is positively correlated with the homogeneity of the tissue.

| S-EPMC5059599 | biostudies-literature

Gaussian Accelerated Molecular Dynamics in NAMD.

Project description:Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that provides efficient free energy calculations of biomolecules. Like the previous accelerated molecular dynamics (aMD), GaMD allows for "unconstrained" enhanced sampling without the need to set predefined collective variables and so is useful for studying complex biomolecular conformational changes such as protein folding and ligand binding. Furthermore, because the boost potential is constructed using a harmonic function that follows Gaussian distribution in GaMD, cumulant expansion to the second order can be applied to recover the original free energy profiles of proteins and other large biomolecules, which solves a long-standing energetic reweighting problem of the previous aMD method. Taken together, GaMD offers major advantages for both unconstrained enhanced sampling and free energy calculations of large biomolecules. Here, we have implemented GaMD in the NAMD package on top of the existing aMD feature and validated it on three model systems: alanine dipeptide, the chignolin fast-folding protein, and the M3 muscarinic G protein-coupled receptor (GPCR). For alanine dipeptide, while conventional molecular dynamics (cMD) simulations performed for 30 ns are poorly converged, GaMD simulations of the same length yield free energy profiles that agree quantitatively with those of 1000 ns cMD simulation. Further GaMD simulations have captured folding of the chignolin and binding of the acetylcholine (ACh) endogenous agonist to the M3 muscarinic receptor. The reweighted free energy profiles are used to characterize the protein folding and ligand binding pathways quantitatively. GaMD implemented in the scalable NAMD is widely applicable to enhanced sampling and free energy calculations of large biomolecules.

| S-EPMC5743237 | biostudies-literature

A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations.

Project description:Four different molecular dynamics (MD) simulations have been performed for ordered DNA decamers, d(5'-ATGCAGTCAG)*d(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine3+ (Spd3+) and putrescine2+ (Put2+), the synthetic polyamine diaminopropane2+ (DAP2+) and Na+. The simulation set-up corresponds to an infinite array of parallel DNA mimicking the state in oriented DNA fibers or crystals. This work describes general properties of polyamine and Na+ binding to DNA. Simulated diffusion coefficients show satisfactory agreement with experimental NMR diffusion data of comparable systems. The interaction of the polyamines with DNA is dynamic in character and the cations mostly form short-lived contacts with the electronegative binding sites of DNA. Polyamines, Na+ and water interact most frequently with the charged phosphate atoms with preference for association from the minor groove side with O1P over O2P. There is a strong anti-correlation in the cation binding to the electronegative groups of DNA, i.e. the presence of a cation near one of the DNA sites repels other cations from binding to this and to the other sites separated by <7.5 A from each other. In contrast to the other polyamines, DAP2+ is able to form 'bridges' connecting neighboring phosphate groups along the DNA strand. A small fraction of DAP2+ and Put2+ can be found in the major grooves, while Spd3+ is absent there. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural (Spd3+, Put2+ and spermine4+) and synthetic (DAP2+) polyamines.

| S-EPMC219480 | biostudies-literature

Implementation of Accelerated Molecular Dynamics in NAMD.

Project description:Accelerated molecular dynamics (aMD) is an enhanced-sampling method that improves the conformational space sampling by reducing energy barriers separating different states of a system. Here we present the implementation of aMD in the parallel simulation program NAMD. We show that aMD simulations performed with NAMD have only a small overhead compared with classical MD simulations. Through example applications to the alanine dipeptide, we discuss the choice of acceleration parameters, the interpretation of aMD results, as well as the advantages and limitations of the aMD method.

| S-EPMC3115733 | biostudies-literature

Accelerated molecular dynamics simulations of protein folding.

Project description:Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native NMR or X-ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second-order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein-folding studies.

| S-EPMC4487363 | biostudies-literature

Gaussian Accelerated Molecular Dynamics: Theory, Implementation, and Applications.

Project description:A novel Gaussian Accelerated Molecular Dynamics (GaMD) method has been developed for simultaneous unconstrained enhanced sampling and free energy calculation of biomolecules. Without the need to set predefined reaction coordinates, GaMD enables unconstrained enhanced sampling of the biomolecules. Furthermore, by constructing a boost potential that follows a Gaussian distribution, accurate reweighting of GaMD simulations is achieved via cumulant expansion to the second order. The free energy profiles obtained from GaMD simulations allow us to identify distinct low energy states of the biomolecules and characterize biomolecular structural dynamics quantitatively. In this chapter, we present the theory of GaMD, its implementation in the widely used molecular dynamics software packages (AMBER and NAMD), and applications to the alanine dipeptide biomolecular model system, protein folding, biomolecular large-scale conformational transitions and biomolecular recognition.

| S-EPMC5927394 | biostudies-literature

Gaussian accelerated molecular dynamics (GaMD): principles and applications.

Project description:Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

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