Order parameters for macromolecules: application to multiscale simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Singharoy A
PROVIDER: S-EPMC3045216 | biostudies-other | 2011 Jan
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Coarse-grained features of macromolecular assemblies are understood via a set of order parameters (OPs) constructed in terms of their all-atom configuration. OPs are shown to be slowly changing in time and capture the large-scale spatial features of macromolecular assemblies. The relationship of these variables to the classic notion of OPs based on symmetry breaking phase transitions is discussed. OPs based on space warping transformations are analyzed in detail as they naturally provide a connection between overall structure of an assembly and all-atom configuration. These OPs serve as the basis of a multiscale analysis that yields Langevin equations for OP dynamics. In this context, the characteristics of OPs and PCA modes are compared. The OPs enable efficient all-atom multiscale simulations of the dynamics of macromolecular assemblies in response to changes in microenvironmental conditions, as demonstrated on the structural transitions of cowpea chlorotic mottle virus capsid (CCMV) and RNA of the satellite tobacco mosaic virus (STMV).
| S-EPMC4937887 | biostudies-literature
Project description:While densely packed DNA arrays are known to exhibit hexagonal and orthorhombic local packings, the detailed mechanism governing the associated phase transition remains rather elusive. Furthermore, at high densities the atomistic resolution is paramount to properly account for fine details, encompassing the DNA molecular order, the contingent ordering of counterions and the induced molecular ordering of the bathing solvent, bringing together electrostatic, steric, thermal and direct hydrogen-bonding interactions, resulting in the observed osmotic equation of state. We perform a multiscale simulation of dense DNA arrays by enclosing a set of 16 atomistically resolved DNA molecules within a semi-permeable membrane, allowing the passage of water and salt ions, and thus mimicking the behavior of DNA arrays subjected to external osmotic stress in a bathing solution of monovalent salt and multivalent counterions. By varying the DNA density, local packing symmetry, and counterion type, we obtain osmotic equation of state together with the hexagonal-orthorhombic phase transition, and full structural characterization of the DNA subphase in terms of its positional and angular orientational fluctuations, counterion distributions, and the solvent local dielectric response profile with its order parameters that allow us to identify the hydration force as the primary interaction mechanism at high DNA densities.
| S-EPMC5500594 | biostudies-literature
Project description:Soft materials (e.g., enveloped viruses, liposomes, membranes and supercooled liquids) simultaneously deform or display collective behaviors, while undergoing atomic scale vibrations and collisions. While the multiple space-time character of such systems often makes traditional molecular dynamics simulation impractical, a multiscale approach has been presented that allows for long-time simulation with atomic detail based on the co-evolution of slowly-varying order parameters (OPs) with the quasi-equilibrium probability density of atomic configurations. However, this approach breaks down when the structural change is extreme, or when nearest-neighbor connectivity of atoms is not maintained. In the current study, a self-consistent approach is presented wherein OPs and a reference structure co-evolve slowly to yield long-time simulation for dynamical soft-matter phenomena such as structural transitions and self-assembly. The development begins with the Liouville equation for N classical atoms and an ansatz on the form of the associated N-atom probability density. Multiscale techniques are used to derive Langevin equations for the coupled OP-configurational dynamics. The net result is a set of equations for the coupled stochastic dynamics of the OPs and centers of mass of the subsystems that constitute a soft material body. The theory is based on an all-atom methodology and an interatomic force field, and therefore enables calibration-free simulations of soft matter, such as macromolecular assemblies.
| S-EPMC3650908 | biostudies-literature
Project description:This paper focuses on the geometric modeling and computational algorithm development of biomolecular structures from two data sources: Protein Data Bank (PDB) and Electron Microscopy Data Bank (EMDB) in the Eulerian (or Cartesian) representation. Molecular surface (MS) contains non-smooth geometric singularities, such as cusps, tips and self-intersecting facets, which often lead to computational instabilities in molecular simulations, and violate the physical principle of surface free energy minimization. Variational multiscale surface definitions are proposed based on geometric flows and solvation analysis of biomolecular systems. Our approach leads to geometric and potential driven Laplace-Beltrami flows for biomolecular surface evolution and formation. The resulting surfaces are free of geometric singularities and minimize the total free energy of the biomolecular system. High order partial differential equation (PDE)-based nonlinear filters are employed for EMDB data processing. We show the efficacy of this approach in feature-preserving noise reduction. After the construction of protein multiresolution surfaces, we explore the analysis and characterization of surface morphology by using a variety of curvature definitions. Apart from the classical Gaussian curvature and mean curvature, maximum curvature, minimum curvature, shape index, and curvedness are also applied to macromolecular surface analysis for the first time. Our curvature analysis is uniquely coupled to the analysis of electrostatic surface potential, which is a by-product of our variational multiscale solvation models. As an expository investigation, we particularly emphasize the numerical algorithms and computational protocols for practical applications of the above multiscale geometric models. Such information may otherwise be scattered over the vast literature on this topic. Based on the curvature and electrostatic analysis from our multiresolution surfaces, we introduce a new concept, the polarized curvature, for the prediction of protein binding sites.
| S-EPMC3855405 | biostudies-literature
Project description:Geometric modeling of biomolecules plays an essential role in the conceptualization of biolmolecular structure, function, dynamics, and transport. Qualitatively, geometric modeling offers a basis for molecular visualization, which is crucial for the understanding of molecular structure and interactions. Quantitatively, geometric modeling bridges the gap between molecular information, such as that from X-ray, NMR, and cryo-electron microscopy, and theoretical/mathematical models, such as molecular dynamics, the Poisson-Boltzmann equation, and the Nernst-Planck equation. In this work, we present a family of variational multiscale geometric models for macromolecular systems. Our models are able to combine multiresolution geometric modeling with multiscale electrostatic modeling in a unified variational framework. We discuss a suite of techniques for molecular surface generation, molecular surface meshing, molecular volumetric meshing, and the estimation of Hadwiger's functionals. Emphasis is given to the multiresolution representations of biomolecules and the associated multiscale electrostatic analyses as well as multiresolution curvature characterizations. The resulting fine resolution representations of a biomolecular system enable the detailed analysis of solvent-solute interaction, and ion channel dynamics, whereas our coarse resolution representations highlight the compatibility of protein-ligand bindings and possibility of protein-protein interactions.
| S-EPMC3760017 | biostudies-literature
Project description:Proteins and other macromolecules have coupled dynamics over multiple time scales (from femtosecond to millisecond and beyond) that make resolving molecular dynamics challenging. We present an approach based on periodically decomposing the dynamics of a macromolecule into slow and fast modes based on a scalable coarse-grained normal mode analysis. A Langevin equation is used to propagate the slowest degrees of freedom while minimizing the nearly instantaneous degrees of freedom. We present numerical results showing that time steps of up to 1000 fs can be used, with real speedups of up to 200 times over plain molecular dynamics. We present results of successfully folding the Fip35 mutant of WW domain.
| S-EPMC4308582 | biostudies-literature
Project description:Multiscale techniques integrating detailed atomistic information on materials and reactions to predict the performance of heterogeneous catalytic full-scale reactors have been suggested but lack seamless implementation. The largest challenges in the multiscale modeling of reactors can be grouped into two main categories: catalytic complexity and the difference between time and length scales of chemical and transport phenomena. Here we introduce the Automated MUltiscale Simulation Environment AMUSE, a workflow that starts from Density Functional Theory (DFT) data, automates the analysis of the reaction networks through graph theory, prepares it for microkinetic modeling, and subsequently integrates the results into a standard open-source Computational Fluid Dynamics (CFD) code. We demonstrate the capabilities of AMUSE by applying it to the unimolecular iso-propanol dehydrogenation reaction and then, increasing the complexity, to the pre-commercial Pd/In2O3 catalyst employed for the CO2 hydrogenation to methanol. The results show that AMUSE allows the computational investigation of heterogeneous catalytic reactions in a comprehensive way, providing essential information for catalyst design from the atomistic to the reactor scale level.
| S-EPMC10694852 | biostudies-literature
Project description:Hopping-type charge transport in an amorphous thin layer composed of organic molecules is simulated by the combined use of molecular dynamics, quantum chemical, and Monte Carlo calculations. By explicitly considering the molecular structure and the disordered intermolecular packing, we reasonably reproduce the experimental hole and electron mobilities and their applied electric field dependence (Poole-Frenkel behaviour) without using any adjustable parameters. We find that the distribution of the density-of-states originating from the amorphous nature has a significant impact on both the mobilities and Poole-Frenkel behaviour. Detailed analysis is also provided to reveal the molecular-level origin of the charge transport, including the origin of Poole-Frenkel behaviour.
| S-EPMC5175281 | biostudies-literature
Project description:Optogenetics has become a key tool for understanding the function of neural circuits and controlling their behavior. An array of directly light driven opsins have been genetically isolated from several families of organisms, with a wide range of temporal and spectral properties. In order to characterize, understand and apply these opsins, we present an integrated suite of open-source, multi-scale computational tools called PyRhO. The purpose of developing PyRhO is three-fold: (i) to characterize new (and existing) opsins by automatically fitting a minimal set of experimental data to three-, four-, or six-state kinetic models, (ii) to simulate these models at the channel, neuron and network levels, and (iii) provide functional insights through model selection and virtual experiments in silico. The module is written in Python with an additional IPython/Jupyter notebook based GUI, allowing models to be fit, simulations to be run and results to be shared through simply interacting with a webpage. The seamless integration of model fitting algorithms with simulation environments (including NEURON and Brian2) for these virtual opsins will enable neuroscientists to gain a comprehensive understanding of their behavior and rapidly identify the most suitable variant for application in a particular biological system. This process may thereby guide not only experimental design and opsin choice but also alterations of the opsin genetic code in a neuro-engineering feed-back loop. In this way, we expect PyRhO will help to significantly advance optogenetics as a tool for transforming biological sciences.
| S-EPMC4834562 | biostudies-other
Project description:The method presented here refines molecular conformations directly against projections of single particles measured by electron microscopy. By optimizing the orientation of the projection at the same time as the conformation, the method is well-suited to two-dimensional class averages from cryoelectron microscopy. Such direct use of two-dimensional images circumvents the need for a three-dimensional density map, which may be difficult to reconstruct from projections due to structural heterogeneity or preferred orientations of the sample on the grid. Our refinement protocol exploits Natural Move Monte Carlo to model a macromolecule as a small number of segments connected by flexible loops, on multiple scales. After tests on artificial data from lysozyme, we applied the method to the Methonococcus maripaludis chaperonin. We successfully refined its conformation from a closed-state initial model to an open-state final model using just one class-averaged projection. We also used Natural Moves to iteratively refine against heterogeneous projection images of Methonococcus maripaludis chaperonin in a mix of open and closed states. Our results suggest a general method for electron microscopy refinement specially suited to macromolecules with significant conformational flexibility. The algorithm is available in the program Methodologies for Optimization and Sampling In Computational Studies.
| S-EPMC3382478 | biostudies-literature