Two Dimensional Window Exchange Umbrella Sampling for Transmembrane Helix Assembly.
Ontology highlight
ABSTRACT:
SUBMITTER: Park S
PROVIDER: S-EPMC3593060 | biostudies-literature | 2013
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A method of window exchange umbrella sampling molecular dynamics simulation is employed for transmembrane helix assembly. An analytical expression for the average acceptance probability between neighboring windows is derived and combined with the first passage time optimization method to predetermine a parameter set in an optimal range. With the parameter set, the method provides a substantially more efficient sampling of helix-helix interfaces together with the potential of mean force along the helix-helix distance of a transmembrane helix-dimer model, compared to the umbrella sampling method.
| S-EPMC3657132 | biostudies-literature
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
Project description:Intrinsically disordered regions of proteins can gain structure by binding to a partner. This process, of coupled folding and binding (CFaB), is a fundamental part of many important biological processes. Structure-based models have proven themselves capable of revealing fundamental aspects of how CFaB occurs, however, typical methods to enhance the sampling of these transitions, such as replica exchange, do not adequately sample the transition state region of this extremely rare process. Here, we use a variant of Umbrella Sampling to enforce sampling of the transition states of CFaB of HdeA monomers at neutral pH, an extremely rare process that occurs over timescales ranging from seconds to hours. Using high resolution sampling in the transition state region, we cluster states along the principal transition path to obtain a detailed description of coupled binding and folding for the HdeA dimer, revealing new insight into the ensemble of states that are accessible to client recognition. We then demonstrate that exchanges between umbrella sampling windows, as done in previous work, significantly improve relaxation in variables orthogonal to the restraints used. Altogether, these results show that Window-Exchange Umbrella Sampling is a promising approach for systems that exhibit flexible binding, and can reveal transition state ensembles of these systems in high detail. © 2015 Wiley Periodicals, Inc.
| S-EPMC4744578 | biostudies-literature
Project description:Molecular dynamics (MD) simulations are increasingly used to study various biological processes such as protein folding, conformational changes, and ligand binding. These processes generally involve slow dynamics that occur on the millisecond or longer timescale, which are difficult to simulate by conventional atomistic MD. Recently, we applied a two-dimensional (2D) replica-exchange MD (REMD) method, which combines the generalized replica exchange with solute tempering (gREST) with the replica-exchange umbrella sampling (REUS) in kinase-inhibitor binding simulations, and successfully observed multiple ligand binding/unbinding events. To efficiently apply the gREST/REUS method to other kinase-inhibitor systems, we establish modified, practical protocols with non-trivial simulation parameter tuning. The current gREST/REUS simulation protocols are tested for three kinase-inhibitor systems: c-Src kinase with PP1, c-Src kinase with Dasatinib, and c-Abl kinase with Imatinib. We optimized the definition of kinase-ligand distance as a collective variable (CV), the solute temperatures in gREST, and replica distributions and umbrella forces in the REUS simulations. Also, the initial structures of each replica in the 2D replica space were prepared carefully by pulling each ligand from and toward the protein binding sites for keeping stable kinase conformations. These optimizations were carried out individually in multiple short MD simulations. The current gREST/REUS simulation protocol ensures good random walks in 2D replica spaces, which are required for enhanced sampling of inhibitor dynamics around a target kinase.
| S-EPMC9099257 | biostudies-literature
Project description:Umbrella sampling efficiently yields equilibrium averages that depend on exploring rare states of a model by biasing simulations to windows of coordinate values and then combining the resulting data with physical weighting. Here, we introduce a mathematical framework that casts the step of combining the data as an eigenproblem. The advantage to this approach is that it facilitates error analysis. We discuss how the error scales with the number of windows. Then, we derive a central limit theorem for averages that are obtained from umbrella sampling. The central limit theorem suggests an estimator of the error contributions from individual windows, and we develop a simple and computationally inexpensive procedure for implementing it. We demonstrate this estimator for simulations of the alanine dipeptide and show that it emphasizes low free energy pathways between stable states in comparison to existing approaches for assessing error contributions. Our work suggests the possibility of using the estimator and, more generally, the eigenvector method for umbrella sampling to guide adaptation of the simulation parameters to accelerate convergence.
| S-EPMC5010559 | biostudies-other
Project description:In this study, we report binding free energy calculations of various drugs-of-abuse to Cucurbit-[8]-uril as part of the SAMPL8 blind challenge. Force-field parameters were obtained from force-matching with different quantum mechanical levels of theory. The Replica Exchange Umbrella Sampling (REUS) approach was used with a cylindrical restraint to enhance the sampling of host-guest binding. Binding free energy was calculated by pulling the guest molecule from one side of the symmetric and cylindrical host, then into and through the host, and out the other side (bidirectional) as compared to pulling only to the bound pose inside the cylindrical host (unidirectional). The initial results with force-matched MP2 parameter set led to RMSE of 4.68 [Formula: see text] from experimental values. However, the follow-up study with CHARMM generalized force field parameters and force-matched PM6-D3H4 parameters resulted in RMSEs from experiment of [Formula: see text] and [Formula: see text], respectively, which demonstrates the potential of REUS for accurate binding free energy calculation given a more suitable description of energetics. Moreover, we compared the free energies for the so called bidirectional and unidirectional free energy approach and found that the binding free energies were highly similar. However, one issue in the bidirectional approach is the asymmetry of profile on the two sides of the host. This is mainly due to the insufficient sampling for these larger systems and can be avoided by longer sampling simulations. Overall REUS shows great promise for binding free energy calculations.
| S-EPMC8131287 | biostudies-literature
Project description:The channel Orai1 requires Ca2+ store depletion in the endoplasmic reticulum and an interaction with the Ca2+ sensor STIM1 to mediate Ca2+ signaling. Alterations in Orai1-mediated Ca2+ influx have been linked to several pathological conditions including immunodeficiency, tubular myopathy, and cancer. We screened large-scale cancer genomics data sets for dysfunctional Orai1 mutants. Five of the identified Orai1 mutations resulted in constitutively active gating and transcriptional activation. Our analysis showed that certain Orai1 mutations were clustered in the transmembrane 2 helix surrounding the pore, which is a trigger site for Orai1 channel gating. Analysis of the constitutively open Orai1 mutant channels revealed two fundamental gates that enabled Ca2+ influx: Arginine side chains were displaced so they no longer blocked the pore, and a chain of water molecules formed in the hydrophobic pore region. Together, these results enabled us to identify a cluster of Orai1 mutations that trigger Ca2+ permeation associated with gene transcription and provide a gating mechanism for Orai1.
| S-EPMC6433236 | biostudies-literature
Project description:Methods that predict membrane helices have become increasingly useful in the context of analyzing entire proteomes, as well as in everyday sequence analysis. Here, we analyzed 27 advanced and simple methods in detail. To resolve contradictions in previous works and to reevaluate transmembrane helix prediction algorithms, we introduced an analysis that distinguished between performance on redundancy-reduced high- and low-resolution data sets, established thresholds for significant differences in performance, and implemented both per-segment and per-residue analysis of membrane helix predictions. Although some of the advanced methods performed better than others, we showed in a thorough bootstrapping experiment based on various measures of accuracy that no method performed consistently best. In contrast, most simple hydrophobicity scale-based methods were significantly less accurate than any advanced method as they overpredicted membrane helices and confused membrane helices with hydrophobic regions outside of membranes. In contrast, the advanced methods usually distinguished correctly between membrane-helical and other proteins. Nonetheless, few methods reliably distinguished between signal peptides and membrane helices. We could not verify a significant difference in performance between eukaryotic and prokaryotic proteins. Surprisingly, we found that proteins with more than five helices were predicted at a significantly lower accuracy than proteins with five or fewer. The important implication is that structurally unsolved multispanning membrane proteins, which are often important drug targets, will remain problematic for transmembrane helix prediction algorithms. Overall, by establishing a standardized methodology for transmembrane helix prediction evaluation, we have resolved differences among previous works and presented novel trends that may impact the analysis of entire proteomes.
| S-EPMC2373751 | biostudies-literature
Project description:Oligosaccharides and polysaccharides exert numerous functional roles in biology through their structural diversity and conformational properties. To investigate their conformational properties using computational methods, Hamiltonian replica exchange (H-REX) combined with two-dimensional grid-based correction maps as biasing potentials (bpCMAP) significantly improves the sampling efficiency about glycosidic linkages. In the current study, we extend the application of H-REX with bpCMAP to complex saccharides and establish systematic procedures for bpCMAP construction, determination of replica distribution, and data analysis. Our main findings are that (1) the bpCMAP for each type of glycosidic linkage can be constructed from the corresponding disaccharide using gas-phase umbrella sampling simulations, (2) the replica distribution can be conveniently determined following the exact definition of the average acceptance ratio based on the assigned distribution of biasing potentials, and (3) the extracted free energy surface (or potential of mean force (PMF)) can be improved using the weighted histogram analysis method (WHAM) allowing for the inclusion of data from the excited state replicas in the calculated probability distribution. The method is applied to a branched N-glycan found on the HIV gp120 protein, and a linear N-glycan. Considering the general importance of N-glycans and the wide appreciation of the sampling problem, the present method represents an efficient procedure for the conformational sampling of complex oligo- and polysaccharides under explicit solvent conditions. More generally, the use of WHAM is anticipated to be of general utility for the calculation of PMFs from H-REX simulations in a wide range of macromolecular systems.
| S-EPMC4334390 | biostudies-literature
Project description:We have previously studied the unfolding equilibrium of bacterioopsin in a single phase solvent, using Förster mechanism fluorescence resonance energy transfer (FRET) as a probe, from tryptophan donors to a dansyl acceptor. We observed an apparent unfolding transition in bacterioopsin perturbed by increasing ethanol concentrations [Nannepaga et al. (2004) Biochemistry 43, 50-59]. We have further investigated this transition and find that the unfolding is pH-dependent. We have now measured the apparent pK of acid-induced unfolding of bacterioopsin in 90% ethanol. When the acceptor is on helix B (Lys 41), the apparent pK for unfolding is 4.75; on the EF connecting loop (Cys 163), 5.15; and on helix G (Cys 222), 5.75. Five-helix proteolytic fragments are less stable. The apparent unfolding pKs are 5.46 for residues 72-248 (Cys 163) and 7.36 for residues 1-166 (Lys 41). When interpreted in terms of a simple equilibrium model for unfolding, the apparent pKs give relative free energies of unfolding in the range of -0.54 to -3.5 kcal/mol. The results suggest that the C-terminal helix of bacterioopsin is less stably folded than the N-terminal helices. We analyzed the pairwise helix-helix interaction surfaces of bacteriorhodopsin and three other seven-transmembrane-helix proteins on the basis of crystal structures. The results show that the interaction surfaces are smoother and the helix axis separations are closer in the amino-terminal two-thirds of the proteins compared with the carboxyl-terminal one-third. However, the F helix is important in stabilizing the folded structure, as shown by the instability of the 1-166 fragment. Considering the high-resolution crystal structure of bacteriorhodopsin, there are no obvious helix-helix interactions involving protein side chains which would be destabilized by protonation at the estimated pH of the unfolding transitions. However, a number of helix-bridging water molecules could become protonated, thereby weakening the helix-helix interactions.
| S-EPMC2519877 | biostudies-literature