Browse
Submit Data
Databases
API
Help

Dataset Information

15 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Low-mass molecular dynamics simulation for configurational sampling enhancement: More evidence and theoretical explanation.

ABSTRACT: It has been reported recently that classical, isothermal-isobaric molecular dynamics (NTP MD) simulations at a time step of 1.00 fs of the standard-mass time (Δt=1.00 fssmt) and a temperature of ≤340 K using uniformly reduced atomic masses by tenfold offers better configurational sampling than standard-mass NTP MD simulations at the same time step. However, it has long been reported that atomic masses can also be increased to improve configurational sampling because higher atomic masses permit the use of a longer time step. It is worth investigating whether standard-mass NTP MD simulations at Δt=2.00 or 3.16 fssmt can offer better or comparable configurational sampling than low-mass NTP MD simulations at Δt=1.00 fssmt. This article reports folding simulations of two β-hairpins showing that the configurational sampling efficiency of NTP MD simulations using atomic masses uniformly reduced by tenfold at Δt=1.00 fssmt is statistically equivalent to and better than those using standard masses at Δt=3.16 and 2.00 fssmt, respectively. The results confirm that, relative to those using standard masses at routine Δt=2.00 fssmt, the low-mass NTP MD simulations at Δt=1.00 fssmt are a simple and generic technique to enhance configurational sampling at temperatures of ≤340 K.

SUBMITTER: Pang YP

PROVIDER: S-EPMC5668912 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Theoretical Explanation for the Rarity of Antibody-Dependent Enhancement of Infection (ADE) in COVID-19.

Project description:Global vaccination against the SARS-CoV-2 virus has proved to be highly effective. However, the possibility of antibody-dependent enhancement of infection (ADE) upon vaccination remains underinvestigated. Here, we aimed to theoretically determine conditions for the occurrence of ADE in COVID-19. We developed a series of mathematical models of antibody response: model Ab-a model of antibody formation; model Cv-a model of infection spread in the body; and a complete model, which combines the two others. The models describe experimental data on SARS-CoV and SARS-CoV-2 infections in humans and cell cultures, including viral load dynamics, seroconversion times and antibody concentration kinetics. The modelling revealed that a significant proportion of macrophages can become infected only if they bind antibodies with high probability. Thus, a high probability of macrophage infection and a sufficient amount of pre-existing antibodies are necessary for the development of ADE in SARS-CoV-2 infection. However, from the point of view of the dynamics of pneumocyte infection, the two cases where the body has a high concentration of preexisting antibodies and a high probability of macrophage infection and where there is a low concentration of antibodies in the body and no macrophage infection are indistinguishable. This conclusion could explain the lack of confirmed ADE cases for COVID-19.

| S-EPMC9569501 | biostudies-literature

On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations.

Project description:We have implemented the accelerated molecular dynamics approach (Hamelberg, D.; Mongan, J.; McCammon, J. A. J. Chem. Phys. 2004, 120 (24), 11919) in the framework of ab initio MD (AIMD). Using three simple examples, we demonstrate that accelerated AIMD (A-AIMD) can be used to accelerate solvent relaxation in AIMD simulations and facilitate the detection of reaction coordinates: (i) We show, for one cyclohexane molecule in the gas phase, that the method can be used to accelerate the rate of the chair-to-chair interconversion by a factor of ?1 × 10(5), while allowing for the reconstruction of the correct canonical distribution of low-energy states; (ii) We then show, for a water box of 64 H(2)O molecules, that A-AIMD can also be used in the condensed phase to accelerate the sampling of water conformations, without affecting the structural properties of the solvent; and (iii) The method is then used to compute the potential of mean force (PMF) for the dissociation of Na-Cl in water, accelerating the convergence by a factor of ?3-4 compared to conventional AIMD simulations.(2) These results suggest that A-AIMD is a useful addition to existing methods for enhanced conformational and phase-space sampling in solution. While the method does not make the use of collective variables superfluous, it also does not require the user to define a set of collective variables that can capture all the low-energy minima on the potential energy surface. This property may prove very useful when dealing with highly complex multidimensional systems that require a quantum mechanical treatment.

| S-EPMC3074571 | biostudies-literature

Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation.

Project description:Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.

| S-EPMC6789927 | biostudies-literature

Network visualization of conformational sampling during molecular dynamics simulation.

Project description:Effective data reduction methods are necessary for uncovering the inherent conformational relationships present in large molecular dynamics (MD) trajectories. Clustering algorithms provide a means to interpret the conformational sampling of molecules during simulation by grouping trajectory snapshots into a few subgroups, or clusters, but the relationships between the individual clusters may not be readily understood. Here we show that network analysis can be used to visualize the dominant conformational states explored during simulation as well as the connectivity between them, providing a more coherent description of conformational space than traditional clustering techniques alone. We compare the results of network visualization against 11 clustering algorithms and principal component conformer plots. Several MD simulations of proteins undergoing different conformational changes demonstrate the effectiveness of networks in reaching functional conclusions.

| S-EPMC3893665 | biostudies-literature

Fine grained sampling of residue characteristics using molecular dynamics simulation.

Project description:In a fine-grained computational analysis of protein structure, we investigated the relationships between a residue's backbone conformations and its side-chain packing as well as conformations. To produce continuous distributions in high resolution, we ran molecular dynamics simulations over a set of protein folds (dynameome). In effect, the dynameome dataset samples not only the states well represented in the PDB but also the known states that are not well represented in the structural database. In our analysis, we characterized the mutual influence among the backbone phi,psi angles with the first side-chain torsion angles (chi(1)) and the volumes occupied by the side-chains. The dependencies of these relationships on side-chain environment and amino acids are further explored. We found that residue volumes exhibit dependency on backbone 2 degrees structure conformation: side-chains pack more densely in extended beta-sheet than in alpha-helical structures. As expected, residue volumes on the protein surface were larger than those in the interior. The first side-chain torsion angles are found to be dependent on the backbone conformations in agreement with previous studies, but the dynameome dataset provides higher resolution of rotamer preferences based on the backbone conformation. All three gauche(-), gauche(+), and trans rotamers show different patterns of phi,psi dependency, and variations in chi(1) value are skewed from their canonical values to relieve the steric strains. By demonstrating the utility of dynameomic modeling on the native state ensemble, this study reveals details of the interplay among backbone conformations, residue volumes and side-chain conformations.

| S-EPMC2916028 | biostudies-literature

Blood oxygenation level-dependent functional MRI signal turbulence caused by ultrahigh spatial resolution: numerical simulation and theoretical explanation.

Project description:High-spatial-resolution functional MRI (fMRI) can enhance image contrast and improve spatial specificity for brain activity mapping. As the voxel size is reduced, an irregular magnetic fieldmap will emerge as a result of less local averaging, and will lead to abnormal fMRI signal evolution with respect to the image acquisition TE. In this article, we report this signal turbulence phenomenon observed in simulations of ultrahigh-spatial-resolution blood oxygenation level-dependent (BOLD) fMRI (voxel size of less than 50 × 50 × 50 µm³). We present a four-level coarse-to-fine multiresolution BOLD fMRI signal simulation. Based on the statistical histogram of an intravoxel fieldmap, we reformulate the intravoxel dephasing summation (a form of Riemann sum) into a new formula that is a discrete Fourier transformation of the intravoxel fieldmap histogram (a form of Lebesgue sum). We interpret the BOLD signal formation by relating its magnitude (phase) to the even (odd) symmetry of the fieldmap histogram. Based on multiresolution BOLD signal simulation, we find that the signal turbulence mainly emerges at the vessel boundary, and that there are only a few voxels (less than 10%) in an ultrahigh-resolution image that reveal turbulence in the form of sparse point noise. Our simulation also shows that, for typical human brain imaging of the cerebral cortex with millimeter resolution, TE < 30 ms and B₀ = 3 T, we are unlikely to observe BOLD signal turbulence. Overall, the main causes of voxel signal turbulence include a high spatial resolution, high field, long TE and large vessel.

| S-EPMC4435714 | biostudies-other

SARS-CoV-2 Omicron Variant Binds to Human Cells More Strongly than the Wild Type: Evidence from Molecular Dynamics Simulation.

Project description:The emergence of the variant of concern Omicron (B.1.1.529) of the severe acute respiratory syndrome coronavirus 2 has aggravated the Covid-19 pandemic due to its very contagious ability. The high infection rate may be due to the high binding affinity of Omicron to human cells, but both experimental and computational studies have yielded conflicting results on this issue. Some studies have shown that the Omicron variant binds to human angiotensin-converting enzyme 2 (hACE2) more strongly than the wild type (WT), but other studies have reported comparable binding affinities. To shed light on this open problem, in this work, we calculated the binding free energy of the receptor binding domain (RBD) of the WT and Omicron spike protein to hACE2 using all-atom molecular dynamics simulation and the molecular mechanics Poisson-Boltzmann surface area method. We showed that Omicron binds to human cells more strongly than the WT due to increased RBD charge, which enhances electrostatic interaction with negatively charged hACE2. N440K, T478K, E484A, Q493R, and Q498R mutations in the RBD have been found to play a critical role in the stability of the RBD-hACE2 complex. The effect of homogeneous and heterogeneous models of glycans coating the viral RBD and the peptidyl domain of hACE2 was examined. Although the total binding free energy is not sensitive to the glycan model, the distribution of per-residue interaction energies depends on it. In addition, glycans have a little effect on the binding affinity of the WT RBD to hACE2.

| S-EPMC9235043 | biostudies-literature

Enhancing Biomolecular Sampling with Reinforcement Learning: A Tree Search Molecular Dynamics Simulation Method.

Project description:This paper proposes a novel molecular simulation method, called tree search molecular dynamics (TS-MD), to accelerate the sampling of conformational transition pathways, which require considerable computation. In TS-MD, a tree search algorithm, called upper confidence bounds for trees, which is a type of reinforcement learning algorithm, is applied to sample the transition pathway. By learning from the results of the previous simulations, TS-MD efficiently searches conformational space and avoids being trapped in local stable structures. TS-MD exhibits better performance than parallel cascade selection molecular dynamics, which is one of the state-of-the-art methods, for the folding of miniproteins, Chignolin and Trp-cage, in explicit water.

| S-EPMC6714528 | biostudies-literature

Configurational and Constitutional Dynamics of Enamine Molecular Switches.

Project description:Dual configurational and constitutional dynamics in systems based on enamine molecular switches has been systematically studied. pH-responsive moieties, such as 2-pyridyl and 2-quinolinyl units, were required on the "stator" part, also providing enamine stability through intramolecular hydrogen-bonding (IMHB) effects. Upon protonation or deprotonation, forward and backward switching could be rapidly achieved. Extension of the stator π-system in the 2-quinolinyl derivative provided a higher E-isomeric equilibrium ratio under neutral conditions, pointing to a means to achieve quantitative forward/backward isomerization processes. The "rotor" part of the enamine switches exhibited constitutional exchange ability with primary amines. Interestingly, considerably higher exchange rates were observed with amines containing ester groups, indicating potential stabilization of the transition state through IMHB. Acids, particularly BiIII , were found to efficiently catalyze the constitutional dynamic processes. In contrast, the enamine and the formed dynamic enamine system showed excellent stability under basic conditions. This coupled configurational and constitutional dynamics expands the scope of dynamic C-C and C-N bonds and potentiates further studies and applications in the fields of molecular machinery and systems chemistry.

| S-EPMC7756271 | biostudies-literature

A theoretical explanation of grain size distributions in explosive rock fragmentation.

Project description:We have measured grain size distributions of the results of laboratory decompression explosions of volcanic rock. The resulting distributions can be approximately represented by gamma distributions of weight per cent as a function of [Formula: see text], where d is the grain size in millimetres measured by sieving, with a superimposed long tail associated with the production of fines. We provide a description of the observations based on sequential fragmentation theory, which we develop for the particular case of 'self-similar' fragmentation kernels, and we show that the corresponding evolution equation for the distribution can be explicitly solved, yielding the long-time lognormal distribution associated with Kolmogorov's fragmentation theory. Particular features of the experimental data, notably time evolution, advection, truncation and fines production, are described and predicted within the constraints of a generalized, 'reductive' fragmentation model, and it is shown that the gamma distribution of coarse particles is a natural consequence of an assumed uniform fragmentation kernel. We further show that an explicit model for fines production during fracturing can lead to a second gamma distribution, and that the sum of the two provides a good fit to the observed data.

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