Multiscale Simulations on the Catalytic Plasticity of CYP76AH1.
Ontology highlight
ABSTRACT:
SUBMITTER: Qiu Y
PROVIDER: S-EPMC8207200 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In this article we review the key modeling tools available for simulating biomolecular systems. We consider recent developments and representative applications of mixed quantum mechanics/molecular mechanics (QM/MM), elastic network models (ENMs), coarse-grained molecular dynamics, and grid-based tools for calculating interactions between essentially rigid protein assemblies. We consider how the different length scales can be coupled, both in a sequential fashion (e.g. a coarse-grained or grid model using parameterization from MD simulations), and via concurrent approaches, where the calculations are performed together and together control the progression of the simulation. We suggest how the concurrent coupling approach familiar in the context of QM/MM calculations can be generalized, and describe how this has been done in the CHARMM macromolecular simulation package.
| S-EPMC6407689 | biostudies-literature
Project description:Multiscale simulations model phenomena across natural scales using monolithic or component-based code, running on local or distributed resources. In this work, we investigate the performance of distributed multiscale computing of component-based models, guided by six multiscale applications with different characteristics and from several disciplines. Three modes of distributed multiscale computing are identified: supplementing local dependencies with large-scale resources, load distribution over multiple resources, and load balancing of small- and large-scale resources. We find that the first mode has the apparent benefit of increasing simulation speed, and the second mode can increase simulation speed if local resources are limited. Depending on resource reservation and model coupling topology, the third mode may result in a reduction of resource consumption.
| S-EPMC4084531 | biostudies-literature
Project description:Smoldyn is a software package for stochastic modelling of spatial biochemical networks and intracellular systems. It was originally developed with an accurate off-lattice particle-based model at its core. This has recently been enhanced with the addition of a computationally efficient on-lattice model, which can be run stand-alone or coupled together for multiscale simulations using both models in regions where they are most required, increasing the applicability of Smoldyn to larger molecule numbers and spatial domains. Simulations can switch between models with only small additions to their configuration file, enabling users with existing Smoldyn configuration files to run the new on-lattice model with any reaction, species or surface descriptions they might already have.Source code and binaries freely available for download at www.smoldyn.org, implemented in C/C++ and supported on Linux, Mac OSX and MS Windows.
| S-EPMC4495299 | biostudies-other
Project description:Protein functions require specific structures frequently coupled with conformational changes. The scale of the structural dynamics of proteins spans from the atomic to the molecular level. Theoretically, all-atom molecular dynamics (MD) simulation is a powerful tool to investigate protein dynamics because the MD simulation is capable of capturing conformational changes obeying the intrinsically structural features. However, to study long-timescale dynamics, efficient sampling techniques and coarse-grained (CG) approaches coupled with all-atom MD simulations, termed multiscale MD simulations, are required to overcome the timescale limitation in all-atom MD simulations. Here, we review two examples of rotary motor proteins examined using free energy landscape (FEL) analysis and CG-MD simulations. In the FEL analysis, FEL is calculated as a function of reaction coordinates, and the long-timescale dynamics corresponding to conformational changes is described as transitions on the FEL surface. Another approach is the utilization of the CG model, in which the CG parameters are tuned using the fluctuation matching methodology with all-atom MD simulations. The long-timespan dynamics is then elucidated straightforwardly by using CG-MD simulations.
| S-EPMC5899732 | biostudies-literature
Project description:When photoactive molecules interact strongly with confined light modes in optical cavities, new hybrid light-matter states form. They are known as polaritons and correspond to coherent superpositions of excitations of the molecules and of the cavity photon. The polariton energies and thus potential energy surfaces are changed with respect to the bare molecules, such that polariton formation is considered a promising paradigm for controlling photochemical reactions. To effectively manipulate photochemistry with confined light, the molecules need to remain in the polaritonic state long enough for the reaction on the modified potential energy surface to take place. To understand what determines this lifetime, we have performed atomistic molecular dynamics simulations of room-temperature ensembles of rhodamine chromophores strongly coupled to a single confined light mode with a 15 fs lifetime. We investigated three popular experimental scenarios and followed the relaxation after optically pumping (i) the lower polariton, (ii) the upper polariton, or (iii) uncoupled molecular states. The results of the simulations suggest that the lifetimes of the optically accessible lower and upper polaritons are limited by (i) ultrafast photoemission due to the low cavity lifetime and (ii) reversible population transfer into the "dark" state manifold. Dark states are superpositions of molecular excitations but with much smaller contributions from the cavity photon, decreasing their emission rates and hence increasing their lifetimes. We find that population transfer between polaritonic modes and dark states is determined by the overlap between the polaritonic and molecular absorption spectra. Importantly, excitation can also be transferred "upward" from the lower polariton into the dark-state reservoir due to the broad absorption spectra of the chromophores, contrary to the common conception of these processes as a "one-way" relaxation from the dark states down to the lower polariton. Our results thus suggest that polaritonic chemistry relying on modified dynamics taking place within the lower polariton manifold requires cavities with sufficiently long lifetimes and, at the same time, strong light-matter coupling strengths to prevent the back-transfer of excitation into the dark states.
| S-EPMC6914212 | biostudies-literature
Project description:A model of myocardial electromechanics is suggested. It combines modified and simplified versions of previously published models of cardiac electrophysiology, excitation-contraction coupling, and mechanics. The mechano-calcium and mechano-electrical feedbacks, including the strain-dependence of the propagation velocity of the action potential, are also accounted for. The model reproduces changes in the twitch amplitude and Ca2+-transients upon changes in muscle strain including the slow response. The model also reproduces the Bowditch effect and changes in the twitch amplitude and duration upon changes in the interstimulus interval, including accelerated relaxation at high stimulation frequency. Special efforts were taken to reduce the stiffness of the differential equations of the model. As a result, the equations can be integrated numerically with a relatively high time step making the model suitable for multiscale simulation of the human heart and allowing one to study the impact of myocardial mechanics on arrhythmias.
| S-EPMC8297763 | biostudies-literature
Project description:Within VECMAtk platform we perform Uncertainty Quantification (UQ) for multiscale fusion plasmas simulations. The goal of VECMAtk is to enable modular and automated tools for a wide range of applications to archive robust and actionable results. Our aim in the current paper is to incorporate suitable features to build UQ workflow over the existing fusion codes and to tackle simulations on high performance parallel computers.
| S-EPMC7304739 | biostudies-literature
Project description:Interactions of membrane proteins with lipid molecules are central to their stability and function. We have used multiscale molecular dynamics simulations to determine the extent to which interactions with lipids are conserved across the aquaporin (Aqp) family of membrane proteins. Simulation-based assessment of the lipid interactions made by Aqps when embedded within a simple phospholipid bilayer agrees well with the protein-lipid contacts determined by electron diffraction from 2D crystals. Extending this simulation-based analysis to all Aqps of known structure reveals a degree of conservation of such interactions across the Aqp structural proteome. Despite similarities in the binding orientations and interactions of the lipids, there do not appear to be distinct, high-specificity lipid binding sites on the surface of Aqps. Rather Aqps exhibit a more broadly conserved protein/lipid interface, suggestive of interchange between annular and bulk lipids, instead of a fixed annular "shell" of lipids.
| S-EPMC3746155 | biostudies-literature
Project description:Integrins are large cell-surface adhesion receptors that can be activated to a high affinity state by the formation of an intracellular complex between the integrin β-subunit tail, the membrane, and talin. The F2 and F3 subdomains of the talin head play a key role in formation of this complex. Here, activation of the integrin αIIb/β3 dimer by the talin head domain was probed using multiscale molecular dynamics simulations. A number of novel insights emerge from these studies, including (i) the importance of the integrin αIIb subunit F992 and F993 residues in stabilizing the "off" state of the αIIb/β3 dimer, (ii) a crucial role for negatively charged groups in the F2-F3/membrane interaction, (iii) binding of the talin F2-F3 domain to negatively charged lipid headgroups in the membrane induces a reorientation of the β transmembrane (TM) domain, (iv) an increase in the tilt angle of the β TM domain relative to the bilayer normal helps to destabilize the α/β TM interaction and promote a scissor-like movement of the integrin TM helices. These results, combined with various published experimental observations, suggest a model for the mechanism of inside-out activation of integrins by talin.
| S-EPMC3141943 | biostudies-literature
Project description:Multiscale modeling and simulations in neuroscience is gaining scientific attention due to its growing importance and unexplored capabilities. For instance, it can help to acquire better understanding of biological phenomena that have important features at multiple scales of time and space. This includes synaptic plasticity, memory formation and modulation, homeostasis. There are several ways to organize multiscale simulations depending on the scientific problem and the system to be modeled. One of the possibilities is to simulate different components of a multiscale system simultaneously and exchange data when required. The latter may become a challenging task for several reasons. First, the components of a multiscale system usually span different spatial and temporal scales, such that rigorous analysis of possible coupling solutions is required. Then, the components can be defined by different mathematical formalisms. For certain classes of problems a number of coupling mechanisms have been proposed and successfully used. However, a strict mathematical theory is missing in many cases. Recent work in the field has not so far investigated artifacts that may arise during coupled integration of different approximation methods. Moreover, in neuroscience, the coupling of widely used numerical fixed step size solvers may lead to unexpected inefficiency. In this paper we address the question of possible numerical artifacts that can arise during the integration of a coupled system. We develop an efficient strategy to couple the components comprising a multiscale test problem in neuroscience. We introduce an efficient coupling method based on the second-order backward differentiation formula (BDF2) numerical approximation. The method uses an adaptive step size integration with an error estimation proposed by Skelboe (2000). The method shows a significant advantage over conventional fixed step size solvers used in neuroscience for similar problems. We explore different coupling strategies that define the organization of computations between system components. We study the importance of an appropriate approximation of exchanged variables during the simulation. The analysis shows a substantial impact of these aspects on the solution accuracy in the application to our multiscale neuroscientific test problem. We believe that the ideas presented in the paper may essentially contribute to the development of a robust and efficient framework for multiscale brain modeling and simulations in neuroscience.
| S-EPMC5018489 | biostudies-literature