Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Analysis of helix-helix interactions of bacteriorhodopsin by replica-exchange simulations.

ABSTRACT: We performed long-time replica-exchange Monte Carlo simulations of bacteriorhodopsin transmembrane helices, which made it possible that wide conformational space was sampled. Using only the helix-helix interactions and starting from random initial configurations, we obtained the nativelike helix arrangement successfully and predicted a part of the configurations (three helices out of seven) precisely. By the principal component analysis we classified low-energy structures into some clusters of similar structures, and we showed that the above nativelike three-helix configuration is reproduced properly in most clusters and that not only the van der Waals interactions but also the electrostatic interactions contributed to the stabilization of the native structures.

SUBMITTER: Kokubo H

PROVIDER: S-EPMC2716567 | biostudies-literature | 2009 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Analysis of helix-helix interactions of bacteriorhodopsin by replica-exchange simulations.

Kokubo Hironori H Okamoto Yuko Y

Biophysical journal 20090201 3

We performed long-time replica-exchange Monte Carlo simulations of bacteriorhodopsin transmembrane helices, which made it possible that wide conformational space was sampled. Using only the helix-helix interactions and starting from random initial configurations, we obtained the nativelike helix arrangement successfully and predicted a part of the configurations (three helices out of seven) precisely. By the principal component analysis we classified low-energy structures into some clusters of s ...[more]

PMID: 18835905

Similar Datasets

Multidimensional umbrella sampling and replica-exchange molecular dynamics simulations for structure prediction of transmembrane helix dimers.

Project description:Structural information of a transmembrane (TM) helix dimer is useful in understanding molecular mechanisms of important biological phenomena such as signal transduction across the cell membrane. Here, we describe an umbrella sampling (US) scheme for predicting the structure of a TM helix dimer in implicit membrane using the interhelical crossing angle and the TM-TM relative rotation angles as the reaction coordinates. This scheme conducts an efficient conformational search on TM-TM contact interfaces, and its robustness is tested by predicting the structures of glycophorin A (GpA) and receptor tyrosine kinase EphA1 (EphA1) TM dimers. The nuclear magnetic resonance (NMR) structures of both proteins correspond to the global free-energy minimum states in their free-energy landscapes. In addition, using the landscape of GpA as a reference, we also examine the protocols of temperature replica-exchange molecular dynamics (REMD) simulations for structure prediction of TM helix dimers in implicit membrane. A wide temperature range in REMD simulations, for example, 250-1000 K, is required to efficiently obtain a free-energy landscape consistent with the US simulations. The interhelical crossing angle and the TM-TM relative rotation angles can be used as reaction coordinates in multidimensional US and be good measures for conformational sampling of REMD simulations.

| S-EPMC4083740 | biostudies-literature

Side chain interactions can impede amyloid fibril growth: replica exchange simulations of Abeta peptide mutant.

Project description:Using replica exchange molecular dynamics, we study the effect of Asp23Tyr mutation on Abeta(10-40) fibril growth. The effect of this mutation is revealed through the computation of free energy landscapes, the distributions of peptide-fibril interactions, and by comparison with the wild-type Abeta(10-40) peptide. Asp23Tyr mutation has a relatively minor influence on the docking of Abeta peptides to the fibril. However, it has a strong impact on the locking stage due to profound stabilization of the parallel in-registry beta-sheets formed by the peptides on the fibril edge. The enhanced stability of parallel beta-sheets results from the deletion of side chain interactions formed by Asp23, which are incompatible with the fibril-like conformers. Consequently, Asp23Tyr mutation is expected to promote fibril growth. We argue that strong off-registry side chain interactions may slow down fibril assembly as it occurs for the wild-type Abeta peptide. The analysis of experimental data offers support to our in silico conclusions.

| S-EPMC2765228 | biostudies-literature

Large-scale asynchronous and distributed multidimensional replica exchange molecular simulations and efficiency analysis.

Project description:We describe methods to perform replica exchange molecular dynamics (REMD) simulations asynchronously (ASyncRE). The methods are designed to facilitate large scale REMD simulations on grid computing networks consisting of heterogeneous and distributed computing environments as well as on homogeneous high-performance clusters. We have implemented these methods on NSF (National Science Foundation) XSEDE (Extreme Science and Engineering Discovery Environment) clusters and BOINC (Berkeley Open Infrastructure for Network Computing) distributed computing networks at Temple University and Brooklyn College at CUNY (the City University of New York). They are also being implemented on the IBM World Community Grid. To illustrate the methods, we have performed extensive (more than 60 ms in aggregate) simulations for the beta-cyclodextrin-heptanoate host-guest system in the context of one- and two-dimensional ASyncRE, and we used the results to estimate absolute binding free energies using the binding energy distribution analysis method. We propose ways to improve the efficiency of REMD simulations: these include increasing the number of exchanges attempted after a specified molecular dynamics (MD) period up to the fast exchange limit and/or adjusting the MD period to allow sufficient internal relaxation within each thermodynamic state. Although ASyncRE simulations generally require long MD periods (>picoseconds) per replica exchange cycle to minimize the overhead imposed by heterogeneous computing networks, we found that it is possible to reach an efficiency similar to conventional synchronous REMD, by optimizing the combination of the MD period and the number of exchanges attempted per cycle.

| S-EPMC4512903 | biostudies-literature

Molecular dynamics simulations using temperature-enhanced essential dynamics replica exchange.

Project description:Today's standard molecular dynamics simulations of moderately sized biomolecular systems at full atomic resolution are typically limited to the nanosecond timescale and therefore suffer from limited conformational sampling. Efficient ensemble-preserving algorithms like replica exchange (REX) may alleviate this problem somewhat but are still computationally prohibitive due to the large number of degrees of freedom involved. Aiming at increased sampling efficiency, we present a novel simulation method combining the ideas of essential dynamics and REX. Unlike standard REX, in each replica only a selection of essential collective modes of a subsystem of interest (essential subspace) is coupled to a higher temperature, with the remainder of the system staying at a reference temperature, T(0). This selective excitation along with the replica framework permits efficient approximate ensemble-preserving conformational sampling and allows much larger temperature differences between replicas, thereby considerably enhancing sampling efficiency. Ensemble properties and sampling performance of the method are discussed using dialanine and guanylin test systems, with multi-microsecond molecular dynamics simulations of these test systems serving as references.

| S-EPMC1877756 | biostudies-literature

Replica exchange simulations of the thermodynamics of Abeta fibril growth.

Project description:Replica exchange molecular dynamics and an all-atom implicit solvent model are used to probe the thermodynamics of deposition of Alzheimer's Abeta monomers on preformed amyloid fibrils. Consistent with the experiments, two deposition stages have been identified. The docking stage occurs over a wide temperature range, starting with the formation of the first peptide-fibril interactions at 500 K. Docking is completed when a peptide fully adsorbs on the fibril edge at the temperature of 380 K. The docking transition appears to be continuous, and occurs without free energy barriers or intermediates. During docking, incoming Abeta monomer adopts a disordered structure on the fibril edge. The locking stage occurs at the temperature of approximately 360 K and is characterized by the rugged free energy landscape. Locking takes place when incoming Abeta peptide forms a parallel beta-sheet structure on the fibril edge. Because the beta-sheets formed by locked Abeta peptides are typically off-registry, the structure of the locked phase differs from the structure of the fibril interior. The study also reports that binding affinities of two distinct fibril edges with respect to incoming Abeta peptides are different. The peptides bound to the concave edge have significantly lower free energy compared to those bound on the convex edge. Comparison with the available experimental data is discussed.

| S-EPMC2716483 | biostudies-other

Coupling of replica exchange simulations to a non-Boltzmann structure reservoir.

Project description:Computing converged ensemble properties remains challenging for large biomolecules. Replica exchange molecular dynamics (REMD) can significantly increase the efficiency of conformational sampling by using high temperatures to escape kinetic traps. Several groups, including ours, introduced the idea of coupling replica exchange to a pre-converged, Boltzmann-populated reservoir, usually at a temperature higher than that of the highest temperature replica. This procedure reduces computational cost because the long simulation times needed for extensive sampling are only carried out for a single temperature. However, a weakness of the approach is that the Boltzmann-weighted reservoir can still be difficult to generate. We now present the idea of employing a non-Boltzmann reservoir, whose structures can be generated through more efficient conformational sampling methods. We demonstrate that the approach is rigorous and derive a correct statistical mechanical exchange criterion between the reservoir and the replicas that drives Boltzmann-weighted probabilities for the replicas. We test this approach on the trpzip2 peptide and demonstrate that the resulting thermal stability profile is essentially indistinguishable from that obtained using very long (>100 ns) standard REMD simulations. The convergence of this reservoir-aided REMD is significantly faster than for regular REMD. Furthermore, we demonstrate that modification of the exchange criterion is essential; REMD simulations using a standard exchange function with the non-Boltzmann reservoir produced incorrect results.

| S-EPMC4819981 | biostudies-literature

Understanding folding and design: replica-exchange simulations of "Trp-cage" miniproteins.

Project description:Replica-exchange molecular dynamics simulations in implicit solvent have been carried out to study the folding thermodynamics of a designed 20-residue peptide, or "miniprotein." The simulations in this study used the amber (parm94) force field along with the generalized Born/solvent-accessible surface area implicit solvent model, and they spanned a range of temperatures from 273 to 630 K. Starting from a completely extended initial conformation, simulations of one peptide sequence sample conformations that are <1.0 A Calpha rms positional deviation from structures in the corresponding NMR ensemble. These folded states are thermodynamically stable with a simulated melting temperature of approximately 400 K, and they satisfy the majority of experimentally observed NMR restraints. Simulations of a related mutant peptide show a degenerate ensemble of states at low temperature, in agreement with experimental results.

| S-EPMC164630 | biostudies-literature

Replica exchange simulations of transient encounter complexes in protein-protein association.

Project description:Recent paramagnetic relaxation enhancement (PRE) studies on several weakly interacting protein complexes have unequivocally demonstrated the existence of transient encounter complexes. Here, we present a computational method to study protein-protein binding by creating equilibrium ensembles that include both specific and nonspecific protein complexes. In a joint analysis of simulation and experiment we explore the physical nature and underlying physicochemical characteristics of encounter complexes involving three protein-protein interactions of the bacterial phosphotransferase system. Replica exchange Monte Carlo simulations using a coarse-grained energy function recover the structures of the specific complexes and produce binding affinities in good agreement with experiment. Together with the specific complex, a relatively small number of distinct nonspecific complexes largely accounts for the measured PRE data. The combined relative population of the latter is less than approximately 10%. The binding interfaces of the specific and nonspecific complexes differ primarily in size but exhibit similar amino acid compositions. We find that the overall funnel-shaped energy landscape of complex formation is dominated by the specific complex, a small number of structured nonspecific complexes, and a diffuse cloud of loosely bound complexes connecting the specific and nonspecific binding sites with each other and the unbound state. Nonspecific complexes may not only accelerate the binding kinetics by enhancing the rate of success of random diffusional encounters but also play a role in protein function as alternative binding modes.

| S-EPMC2529066 | biostudies-literature

Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations.

Project description:2-Hydroxypropyl-β-cyclodextrin (HPβCD) has unique properties to enhance the stability and the solubility of low water-soluble compounds by inclusion complexation. An understanding of the structural properties of HPβCD and its derivatives, based on the number of 2-hydroxypropyl (HP) substituents at the α-d-glucopyranose subunits is rather important. In this work, replica exchange molecular dynamics simulations were performed to investigate the conformational changes of single- and double-sided HP-substitution, called 6-HPβCDs and 2,6-HPβCDs, respectively. The results show that the glucose subunits in both 6-HPβCDs and 2,6-HPβCDs have a lower chance of flipping than in βCD. Also, HP groups occasionally block the hydrophobic cavity of HPβCDs, thus hindering drug inclusion. We found that HPβCDs with a high number of HP-substitutions are more likely to be blocked, while HPβCDs with double-sided HP-substitutions have an even higher probability of being blocked. Overall, 6-HPβCDs with three and four HP-substitutions are highlighted as the most suitable structures for guest encapsulation, based on our conformational analyses, such as structural distortion, the radius of gyration, circularity, and cavity self-closure of the HPβCDs.

| S-EPMC6401915 | biostudies-literature

Exploring Protocols to Build Reservoirs to Accelerate Temperature Replica Exchange MD Simulations.

Project description:Temperature replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method for accelerating biomolecular simulations. During the past 2 decades, several variants of REMD have been developed to further improve the rate of conformational sampling of REMD. One such variant, reservoir REMD (RREMD), was shown to improve the rate of conformational sampling by around 5-20×. Despite the significant increase in the sampling speed, RREMD methods have not been widely used because of the difficulties in building the reservoir and also because of the code not being available on the graphics processing units (GPUs). In this work, we ported the Amber RREMD code onto GPUs making it 20× faster than the central processing unit code. Then, we explored protocols for building Boltzmann-weighted reservoirs as well as non-Boltzmann reservoirs and tested how each choice affects the accuracy of the resulting RREMD simulations. We show that, using the recommended protocols outlined here, RREMD simulations can accurately reproduce Boltzmann-weighted ensembles obtained by much more expensive conventional temperature-based REMD simulations, with at least 15× faster convergence rates even for larger proteins (>50 amino acids) compared to conventional REMD.

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