Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Krylov subspace methods for computing hydrodynamic interactions in brownian dynamics simulations.

ABSTRACT: Hydrodynamic interactions play an important role in the dynamics of macromolecules. The most common way to take into account hydrodynamic effects in molecular simulations is in the context of a brownian dynamics simulation. However, the calculation of correlated brownian noise vectors in these simulations is computationally very demanding and alternative methods are desirable. This paper studies methods based on Krylov subspaces for computing brownian noise vectors. These methods are related to Chebyshev polynomial approximations, but do not require eigenvalue estimates. We show that only low accuracy is required in the brownian noise vectors to accurately compute values of dynamic and static properties of polymer and monodisperse suspension models. With this level of accuracy, the computational time of Krylov subspace methods scales very nearly as O(N(2)) for the number of particles N up to 10 000, which was the limit tested. The performance of the Krylov subspace methods, especially the "block" version, is slightly better than that of the Chebyshev method, even without taking into account the additional cost of eigenvalue estimates required by the latter. Furthermore, at N = 10,000, the Krylov subspace method is 13 times faster than the exact Cholesky method. Thus, Krylov subspace methods are recommended for performing large-scale brownian dynamics simulations with hydrodynamic interactions.

SUBMITTER: Ando T

PROVIDER: S-EPMC3427343 | biostudies-other | 2012 Aug

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

Predicting protein interactions by Brownian dynamics simulations.

Project description:We present a newly adapted Brownian-Dynamics (BD)-based protein docking method for predicting native protein complexes. The approach includes global BD conformational sampling, compact complex selection, and local energy minimization. In order to reduce the computational costs for energy evaluations, a shell-based grid force field was developed to represent the receptor protein and solvation effects. The performance of this BD protein docking approach has been evaluated on a test set of 24 crystal protein complexes. Reproduction of experimental structures in the test set indicates the adequate conformational sampling and accurate scoring of this BD protein docking approach. Furthermore, we have developed an approach to account for the flexibility of proteins, which has been successfully applied to reproduce the experimental complex structure from the structure of two unbounded proteins. These results indicate that this adapted BD protein docking approach can be useful for the prediction of protein-protein interactions.

| S-EPMC3303761 | biostudies-literature

A molecule-centered method for accelerating the calculation of hydrodynamic interactions in Brownian dynamics simulations containing many flexible biomolecules.

Project description:Inclusion of hydrodynamic interactions (HIs) is essential in simulations of biological macromolecules that treat the solvent implicitly if the macromolecules are to exhibit correct translational and rotational diffusion. The present work describes the development and testing of a simple approach aimed at allowing more rapid computation of HIs in coarse-grained Brownian dynamics simulations of systems that contain large numbers of flexible macromolecules. The method combines a complete treatment of intramolecular HIs with an approximate treatment of the intermolecular HIs which assumes that the molecules are effectively spherical; all of the HIs are calculated at the Rotne-Prager-Yamakawa level of theory. When combined with Fixman's Chebyshev polynomial method for calculating correlated random displacements, the proposed method provides an approach that is simple to program but sufficiently fast that it makes it computationally viable to include HIs in large-scale simulations. Test calculations performed on very coarse-grained models of the pyruvate dehydrogenase (PDH) E2 complex and on oligomers of ParM (ranging in size from 1 to 20 monomers) indicate that the method reproduces the translational diffusion behavior seen in more complete HI simulations surprisingly well; the method performs less well at capturing rotational diffusion but its discrepancies diminish with increasing size of the simulated assembly. Simulations of residue-level models of two tetrameric protein models demonstrate that the method also works well when more structurally detailed models are used in the simulations. Finally, test simulations of systems containing up to 1024 coarse-grained PDH molecules indicate that the proposed method rapidly becomes more efficient than the conventional BD approach in which correlated random displacements are obtained via a Cholesky decomposition of the complete diffusion tensor.

| S-EPMC3731167 | biostudies-literature

Absolute protein-protein association rate constants from flexible, coarse-grained Brownian dynamics simulations: the role of intermolecular hydrodynamic interactions in barnase-barstar association.

Project description:Theory and computation have long been used to rationalize the experimental association rate constants of protein-protein complexes, and Brownian dynamics (BD) simulations, in particular, have been successful in reproducing the relative rate constants of wild-type and mutant protein pairs. Missing from previous BD studies of association kinetics, however, has been the description of hydrodynamic interactions (HIs) between, and within, the diffusing proteins. Here we address this issue by rigorously including HIs in BD simulations of the barnase-barstar association reaction. We first show that even very simplified representations of the proteins--involving approximately one pseudoatom for every three residues in the protein--can provide excellent reproduction of the absolute association rate constants of wild-type and mutant protein pairs. We then show that simulations that include intermolecular HIs also produce excellent estimates of association rate constants, but, for a given reaction criterion, yield values that are decreased by ∼35-80% relative to those obtained in the absence of intermolecular HIs. The neglect of intermolecular HIs in previous BD simulation studies, therefore, is likely to have contributed to the somewhat overestimated absolute rate constants previously obtained. Consequently, intermolecular HIs could be an important component to include in accurate modeling of the kinetics of macromolecular association events.

| S-EPMC2965997 | biostudies-literature

Brownian dynamics simulations of interactions between aldolase and G- or F-actin.

Project description:Compartmentation of proteins in cells is important to proper cell function. Interactions of F-actin and glycolytic enzymes is one mechanism by which glycolytic enzymes can compartment. Brownian dynamics (BD) simulations of the binding of the muscle form of the glycolytic enzyme fructose-1,6-bisphosphate aldolase (aldolase) to F- or G-actin provide first-encounter snapshots of these interactions. Using x-ray structures of aldolase, G-actin, and three-dimensional models of F-actin, the electrostatic potential about each protein was predicted by solving the linearized Poisson-Boltzmann equation for use in BD simulations. The BD simulations provided solution complexes of aldolase with F- or G-actin. All complexes demonstrate the close contacts between oppositely charged regions of the protein surfaces. Positively charged surface regions of aldolase (residues Lys 13, 27, 288, 293, and 341 and Arg 257) are attracted to the negatively charged amino terminus (Asp 1 and Glu 2 and 4) and other patches (Asp 24, 25, and 363 and Glu 361, 364, 99, and 100) of actin subunits. According to BD results, the most important factor for aldolase binding to actin is the quaternary structure of aldolase and actin. Two pairs of adjacent aldolase subunits greatly add to the positive electrostatic potential of each other creating a region of attraction for the negatively charged subdomain 1 of the actin subunit that is exposed to solvent in the quaternary F-actin structure.

| S-EPMC1302496 | biostudies-other

Brownian reservoir computing realized using geometrically confined skyrmion dynamics.

Project description:Reservoir computing (RC) has been considered as one of the key computational principles beyond von-Neumann computing. Magnetic skyrmions, topological particle-like spin textures in magnetic films are particularly promising for implementing RC, since they respond strongly nonlinearly to external stimuli and feature inherent multiscale dynamics. However, despite several theoretical proposals that exist for skyrmion reservoir computing, experimental realizations have been elusive until now. Here, we propose and experimentally demonstrate a conceptually new approach to skyrmion RC that leverages the thermally activated diffusive motion of skyrmions. By confining the electrically gated and thermal skyrmion motion, we find that already a single skyrmion in a confined geometry suffices to realize nonlinearly separable functions, which we demonstrate for the XOR gate along with all other Boolean logic gate operations. Besides this universality, the reservoir computing concept ensures low training costs and ultra-low power operation with current densities orders of magnitude smaller than those used in existing spintronic reservoir computing demonstrations. Our proposed concept is robust against device imperfections and can be readily extended by linking multiple confined geometries and/or by including more skyrmions in the reservoir, suggesting high potential for scalable and low-energy reservoir computing.

| S-EPMC9666653 | biostudies-literature

Brownian dynamics simulations of ion transport through the VDAC.

Project description:It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant was used in the GCMC/BD simulations, the current-voltage characteristics and ion number profiles inside the pore showed excellent agreement with those calculated from all-atom molecular-dynamics (MD) simulations, thereby validating the GCMC/BD approach. Of the 20 NMR models of hVDAC1 currently available, the third one (NMR03) best reproduces both experimental single-channel conductance and ion selectivity (i.e., the reversal potential). In addition, detailed analyses of the ion trajectories, one-dimensional multi-ion potential of mean force, and protein charge distribution reveal that electrostatic interactions play an important role in the channel structure and ion transport relationship. Finally, the GCMC/BD simulations of various mutants based on NMR03 show good agreement with experimental ion selectivity. The difference in ion selectivity between the wild-type and the mutants is the result of altered potential of mean force profiles that are dominated by the electrostatic interactions.

| S-EPMC3030170 | biostudies-literature

Computing 3D Chromatin Configurations from Contact Probability Maps by Inverse Brownian Dynamics.

Project description:The three-dimensional (3D) organization of chromatin, on the length scale of a few genes, is crucial in determining the functional state-accessibility and amount of gene expression-of the chromatin. Recent advances in chromosome conformation capture experiments provide partial information on the chromatin organization in a cell population, namely the contact count between any segment pairs, but not on the interaction strength that leads to these contact counts. However, given the contact matrix, determining the complete 3D organization of the whole chromatin polymer is an inverse problem. In this work, a novel inverse Brownian dynamics method based on a coarse-grained bead-spring chain model has been proposed to compute the optimal interaction strengths between different segments of chromatin such that the experimentally measured contact count probability constraints are satisfied. Applying this method to the α-globin gene locus in two different cell types, we predict the 3D organizations corresponding to active and repressed states of chromatin at the locus. We show that the average distance between any two segments of the region has a broad distribution and cannot be computed as a simple inverse relation based on the contact probability alone. The results presented for multiple normalization methods suggest that all measurable quantities may crucially depend on the nature of normalization. We argue that by experimentally measuring predicted quantities, one may infer the appropriate form of normalization.

| S-EPMC7203009 | biostudies-literature

Motion of a nano-spheroid in a cylindrical vessel flow: Brownian and hydrodynamic interactions.

Project description:We study the motion of a buoyant or a nearly neutrally buoyant nano-sized spheroid in a fluid filled tube without or with an imposed pressure gradient (weak Poiseuille flow). The fluctuating hydrodynamics approach and the deterministic method are both employed. We ensure that the fluctuation-dissipation relation and the principle of thermal equipartition of energy are both satisfied. The major focus is on the effect of the confining boundary. Results for the velocity and the angular velocity autocorrelations (VACF and AVACF), the diffusivities and the drag and the lift forces as functions of the shape, the aspect ratio, the inclination angle and the proximity to the wall are presented. For the parameters considered, the boundary modifies the VACF and AVACF such that three distinct regimes are discernible - an initial exponential decay followed by an algebraic decay culminating in a second exponential decay. The first is due to the thermal noise, the algebraic regime is due both to the thermal noise and the hydrodynamic correlations, while the second exponential decay shows the effect of momentum reflection from the confining wall. Our predictions display excellent comparison with published results for the algebraic regime (the only regime for which earlier results exist). We also discuss the role of the off-diagonal elements of the mobility and the diffusivity tensors that enable the quantifications of the degree of lift and margination of the nanocarrier. Our study covers a range of parameters that are of wide applicability in nanotechnology, microrheology and in targeted drug delivery.

| S-EPMC5669124 | biostudies-literature

Coarse-Grained Brownian Dynamics Simulations of the 10-23 DNAzyme.

Project description:Deoxyribozymes (DNAzymes) are single-stranded DNA that catalyze nucleic acid biochemistry. Although a number of DNAzymes have been discovered by in vitro selection, the relationship between their tertiary structure and function remains unknown. We focus here on the well-studied 10-23 DNAzyme, which cleaves mRNA with a catalytic efficiency approaching that of RNase A. Using coarse-grained Brownian dynamics simulations, we find that the DNAzyme bends its substrate away from the cleavage point, exposing the reactive site and buckling the DNAzyme catalytic core. This hypothesized transition state provides microscopic insights into experimental observations concerning the size of the DNAzyme/substrate complex, the impact of the recognition arm length, and the sensitivity of the enzymatic activity to point mutations of the catalytic core. Upon cleaving the pertinent backbone bond in the substrate, we find that the catalytic core of the DNAzyme unwinds and the overall complex rapidly extends, in agreement with experiments on the related 8-17 DNAzyme. The results presented here provide a starting point for interpreting experimental data on DNAzyme kinetics, as well as developing more detailed simulation models. The results also demonstrate the limitations of using a simple physical model to understand the role of point mutations.

| S-EPMC2776252 | biostudies-literature

Brownian-dynamics simulations of metal-ion binding to four-way junctions.

Project description:Four-way junctions (4Hs) are important intermediates in DNA rearrangements such as genetic recombination. Under the influence of multivalent cations these molecules undergo a conformational change, from an extended planar form to a quasi-continuous stacked X-structure. Recently, a number of X-ray structures and a nuclear magnetic resonance (NMR) structure of 4Hs have been reported and in three of these the position of multivalent cations is revealed. These structures belong to two main families, characterized by the angle between the two co-axial stacked helices, which is either around +40 to +55 degrees or around -70 to -80 degrees. To investigate the role of metal-ion binding on the conformation of folded 4Hs we performed Brownian-dynamics simulations on the set of available structures. The simulations confirm the proposed metal-ion binding sites in the NMR structure and in one of the X-ray structures. Furthermore, the calculations suggest positions for metal-ion binding in the other X-ray structures. The results show a striking dependence of the ion density on the helical environment (B-helix or A-helix) and the structural family.

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