Replica Exchange Molecular Dynamics Study Research Articles

Can Free Energy Perturbation Simulations Coupled with Replica-Exchange Molecular Dynamics Study Ligands with Distributed Binding Sites?

Free energy perturbation coupled with replica exchange with solute tempering (FEP/REST) offers a rigorous approach to compute relative free energy changes for ligands. To determine the applicability of FEP/REST for the ligands with distributed binding poses, we considered two alchemical transformations involving three putative inhibitors I0, I1, and I2 of the Venezuelan equine encephalitis virus nuclear localization signal sequence binding to the importin-α (impα) transporter protein. I0 → I1 and I0 → I2 transformations, respectively, increase or decrease the polarity of the parent molecule. Our objective was three-fold─(i) to verify FEP/REST technical performance and convergence, (ii) to estimate changes in binding free energy ΔΔG, and (iii) to determine the utility of FEP/REST simulations for conformational binding analysis. Our results are as follows. First, our FEP/REST implementation properly follows FEP/REST formalism and produces converged ΔΔG estimates. Due to ligand inherent unbinding, the better FEP/REST strategy lies in performing multiple independent trajectories rather than extending their length. Second, I0 → I1 and I0 → I2 transformations result in overall minor changes in inhibitor binding free energy, slightly strengthening the affinity of I1 and weakening that of I2. Electrostatic interactions dominate binding interactions, determining the enthalpic changes. The two transformations cause opposite entropic changes, which ultimately govern binding affinities. Importantly, we confirm the validity of FEP/REST free energy estimates by comparing them with our previous REST simulations, directly probing binding of three ligands to impα. Third, we established that FEP/REST simulations can sample binding ensembles of ligands. Thus, FEP/REST can be applied (i) to study the energetics of the ligand binding without defined poses and showing minor differences in affinities |ΔΔG| ≲ 0.5 kcal/mol and (ii) to collect ligand binding conformational ensembles.

The influences of E22Q mutant on solvated 3Aβ11-40 peptide: A REMD study

The residue E22 plays a critical role in the aggregation process of Amyloid beta (Aβ) peptides. The effect of E22Q mutant on the shapes of the solvated Aβ11-40 trimer is clarified using a replica exchange molecular dynamics (REMD) simulation employing ∼20.6 μs of MD simulations with 48 disparate replicas. The increase of intramolecular polar contacts and salt bridge between the residue D23 to residues (24–29) was observed. The residual secondary structure of the mutated trimer is shifted in a similar way to the picture observed in previous investigations of F19W mutant. The free energy surface (FES) of the mutated E22Q system has a fewer number of minima in comparison with the wild-type trimer. The optimized shapes of the mutated E22Q form a significant increase in beta structure (47%) and serious decrease in coil content (46%) compared with the wild-type (of 36 and 56%, respectively). The binding affinity of constituting chains to the rest is of −43.7 ± 6.5 kcal/mol, implying that the representative structure of E22Q is more stable than the wild-type one. Furthermore, the E22Q mutant increases the size of stable structures due to larger collision cross section (CCS) and solvent accessible area (SASA). The observed results may enhance the Aβ inhibition throughout the contribution to the knowledge of the Aβ oligomerization/aggregation.

Aqueous ionic liquids influence the disulfide bond isoform equilibrium in conotoxin AuIB: a consequence of the Hofmeister effect?

The appearance of several disulfide bond isoforms in multiple cysteine containing venom peptides poses a significant challenge in their synthesis and purification under laboratory conditions. Recent experiments suggest that careful tuning of solvent and temperature conditions can propel the disulfide bond isoform equilibrium in favor of the most potent, native form. Certain aqueous ionic liquids (ILs) have proven significantly useful as solvents for this purpose, while exceptions have also been noted. To elucidate the molecular level origin behind such a preference, we report a detailed explicit solvent replica exchange molecular dynamics study of a conotoxin, AuIB, in pure water and four different aqueous IL solutions (~45-60% v/v). The ILs studied here are comprised of cations like 1-ethyl-3-methyl-imidazolium (Im21+) or 1-butyl-3-methyl-imidazolium (Im41+) coupled with either acetate (OAc-) or chloride (Cl-) as the counter anion. Our simulations unfold interesting features of the conformational spaces sampled by the peptide and its solvation in pure water and aqueous IL solutions. Detailed investigation into populations of the globular disulfide bond isoform of AuIB in aqueous IL solutions reveal distinct trends which might be related to the Hofmeister effect of the cation and anion of the IL and of specific interactions of the aqueous IL solutions with the peptide. In accordance with experimental observations, the aqueous [Im21][OAc] solution is found to promote the highest globular isoform population in AuIB.

Replica exchange molecular dynamics study of the truncated amyloid beta (11-40) trimer in solution.

Amyloid beta (Aβ) oligomers are neurotoxic compounds that destroy the brain of Alzheimer's disease patients. Recent studies indicated that the trimer is one of the most cytotoxic forms of low molecular weight Aβ oligomers. As there was limited information about the structure of the Aβ trimer, either by experiment or by computation, we determined in this work the structure of the 3Aβ11-40 oligomer for the first time using the temperature replica exchange molecular dynamics simulations in the presence of an explicit solvent. More than 20.0 μs of MD simulations were performed. The probability of the β-content and random coil structure of the solvated trimer amounts to 42 ± 6 and 49 ± 7% which is in good agreement with experiments. Intermolecular interactions in central hydrophobic cores play a key role in stabilizing the oligomer. Intermolecular polar contacts between D23 and residues 24-29 replace the salt bridge D23-K28 to secure the loop region. The hydrophilic region of the N-terminus is maintained by the intermolecular polar crossing contacts H13A-Q15B and H13B-Q15C. The difference in the free energy of binding between the constituting monomers and the others amounts to -36 ± 8 kcal mol-1. The collision cross section of the representative structures of the trimer was computed to be 1330 ± 47 Å2, which is in good agreement with previous experiments.

Replica exchange molecular dynamics study of the amyloid beta (11–40) trimer penetrating a membrane

The transmembrane Aβ11–40 trimer is investigated for the first time using REMD and FEP.

Replica Exchange Molecular Dynamics Study of Dimerization in Prion Protein: Multiple Modes of Interaction and Stabilization.

The pathological forms of prions are known to be a result of misfolding, oligomerization, and aggregation of the cellular prion. While the mechanism of misfolding and aggregation in prions has been widely studied using both experimental and computational tools, the structural and energetic characterization of the dimer form have not garnered as much attention. On one hand dimerization can be the first step toward a nucleation-like pathway to aggregation, whereas on the other hand it may also increase the conformational stability preventing self-aggregation. In this work, we have used extensive all-atom replica exchange molecular dynamics simulations of both monomer and dimer forms of a mouse prion protein to understand the structural, dynamic, and thermodynamic stability of dimeric prion as compared to the monomeric form. We show that prion proteins can dimerize spontaneously being stabilized by hydrophobic interactions as well as intermolecular hydrogen bonding and salt bridge formation. We have computed the conformational free energy landscapes for both monomer and dimer forms to compare the thermodynamic stability and misfolding pathways. We observe large conformational heterogeneity among the various modes of interactions between the monomers and the strong intermolecular interactions may lead to as high as 20% β-content. The hydrophobic regions in helix-2, surrounding coil regions, terminal regions along with the natively present β-sheet region appear to actively participate in prion-prion intermolecular interactions. Dimerization seems to considerably suppress the inherent dynamic instability observed in monomeric prions, particularly because the regions of structural frustration constitute the dimer interface. Further, we demonstrate an interesting reversible coupling between the Q160-G131 interaction (which leads to inhibition of β-sheet extension) and the G131-V161 H-bond formation.

The intrinsic helical propensities of the helical fragments in prion protein under neutral and low pH conditions: a replica exchange molecular dynamics study

Replica exchange molecular dynamics simulations in neutral and acidic aqueous solutions were employed to study the intrinsic helical propensities of three helices in both Syrian hamster (syPrP) and human (huPrP) prion proteins. The helical propensities of syPrP HA and huPrP HA are very high under both pH conditions, which implies that HA is barely involved in the helix-to-β transition. The SyPrP HB chain has a strong tendency to adopt an extended conformation, which is possibly involved in the mechanism of infectious prion diseases in Syrian hamster. HuPrP HC has more of a preference for the extended conformation than huPrP HA and huPrP HB do, which leads to the conjecture that it is more likely to be the source of β-rich structure for human prion protein. We also noticed that the presence of salt bridges is not correlated with helical propensity, indicating that salt bridges do not stabilize helices.

Investigation of the Polymeric Properties of α-Synuclein and Comparison with NMR Experiments: A Replica Exchange Molecular Dynamics Study.

Intrinsically disordered proteins (IDPs) have been shown to be involved in a number of cellular functions, in addition to their predominance in diseased states. α-synuclein may be described as one such IDP implicated in the pathology of Parkinson's disease. Understanding the conformational characteristics of the monomeric state of α-synuclein is necessary for understanding the role of the monomer conformation in aggregation. Polymer theories have been applied to investigate the statistical properties of homopolymeric IDPs. Here we use Replica Exchange Molecular Dynamics (REMD) simulations using temperature as a proxy for solvent quality to examine how well these theories developed for homopolymeric chains describe heteropolymeric α-synuclein. Our results indicate that α-synuclein behaves like a homopolymer at the extremes of solvent quality, while in the intermediate solvent regime, the uneven distribution of charged residues along the sequence strongly influences the conformations adopted by the chain. We refine the ensemble extracted from the REMD simulations of α-synuclein, which shows the best qualitative agreement with experiment, by fitting to the experimental NMR Residual Dipolar Couplings (RDCs) and Paramagnetic Relaxation Enhancements (PREs). Our results demonstrate that the detailed shape of the RDC patterns are sensitive to the angular correlations that are local in sequence while longer range anti-correlations which arise from packing constraints affect the RDC magnitudes.

Low molecular weight oligomers of amyloid peptides display β-barrel conformations: A replica exchange molecular dynamics study in explicit solvent

The self-assembly of proteins and peptides into amyloid fibrils is connected to over 40 pathological conditions including neurodegenerative diseases and systemic amyloidosis. Diffusible, low molecular weight protein and peptide oligomers that form in the early steps of aggregation appear to be the harmful cytotoxic species in the molecular etiology of these diseases. So far, the structural characterization of these oligomers has remained elusive owing to their transient and dynamic features. We here address, by means of full atomistic replica exchange molecular dynamics simulations, the energy landscape of heptamers of the amyloidogenic peptide NHVTLSQ from the beta-2 microglobulin protein. The simulations totaling 5 micros show that low molecular weight oligomers in explicit solvent consist of beta-barrels in equilibrium with amorphous states and fibril-like assemblies. The results, also accounting for the influence of the pH on the conformational properties, provide a strong evidence of the formation of transient beta-barrel assemblies in the early aggregation steps of amyloid-forming systems. Our findings are discussed in terms of oligomers cytotoxicity.

How Is cis-trans Isomerization Controlled in Dronpa Mutants? A Replica Exchange Molecular Dynamics Study.

The reversibly photoactivatable green fluorescent protein analog Dronpa holds great promise as a marker for various new cellular imaging applications. Using a replica exchange method which combines both Hamiltonian and temperature exchanges, the ground-state dynamics of Dronpa and two mutants with increased switching kinetics, Val157Gly and Met159Thr, were compared. The dominant chromophore state was found to be the cis isomer in all three proteins. The simulation data suggest that both mutations strongly increase the chromophore flexibility and cis-trans isomerization rate. We identify three key amino acids, Val157, Met159, and Phe173, which are able to impede the bottom hula-twist transition path, depending on their position and rotameric state. We believe our insights will help to understand the switching process and provide useful information for the design of new variants with improved fluorescence properties.

Simulation of the Pressure and Temperature Folding/Unfolding Equilibrium of a Small RNA Hairpin

We report molecular dynamics simulations of the equilibrium folding/unfolding thermodynamics of the RNA tetraloop in explicit solvent. A replica exchange molecular dynamics study of the r(CGUUGCCG) oligomer that forms a hairpin is performed for 226 ns per replica, using 52 replicas. We are able to show the unbiased folding of all replicas starting from extended conformations. The equilibrium pressure-temperature free energy of folding, DeltaG(P,T), is calculated from the averaged energy, pressure, and specific volume change upon folding of the oligomer as a function of T at constant volume. We find that this oligomer is destabilized by increasing hydrostatic pressure, similar to the behavior of globular proteins.

Free energy landscape and folding mechanism of a β‐hairpin in explicit water: A replica exchange molecular dynamics study

The free energy landscape and the folding mechanism of the C-terminal beta-hairpin of protein G is studied by extensive replica exchange molecular dynamics simulations (40 replicas and 340 ns total simulation time), using the GROMOS96 force field and the SPC explicit water solvent. The study reveals that the system preferentially adopts a beta-hairpin structure at biologically important temperatures, and that the helix content is low at all temperatures studied. Representing the free energy landscape as a function of several types of reaction coordinates, four local minima corresponding to the folded, partially folded, molten globule, and unfolded states are identified. The findings suggest that the folding of the beta-hairpin occurs as the sequence: collapse of hydrophobic core --> formation of H-bond --> formation of the turn. Identifying the folded and molten globule states as the main conformations, the free energy landscape of the beta-hairpin is consistent with a two-state behavior with a broad transition state. The temperature dependence of the folding-unfolding transition is investigated in some detail. The enthalpy and entropy jumps at the folding transition temperature are found to be about three times lower than the experimental estimates, indicating that the folding-unfolding transition in silico is less cooperative than its in vitro counterpart.

Structure and energy landscape of a photoswitchable peptide: A replica exchange molecular dynamics study

A replica exchange molecular dynamics (REMD) simulation of a bicyclic azobenzene peptide in explicit dimethyl sulfoxide solution is presented in order to characterize the conformational structures and energy landscape of a photoswitchable peptide. It is shown that an enhanced-sampling technique such as the REMD method is essential to obtain a converged conformational sampling of the peptide at room temperature. This is because conventional MD simulations of less than approximately 100-ns length are either trapped in local minima (at 295 K) or-if run at high temperature-do not resemble the room-temperature REMD results. Calculating various nuclear Overhauser effects (NOEs) and (3)J-couplings, a good overall agreement between the REMD simulations and the NMR experiments of Renner et al. (Biopolymers 2000;54:501-514) is found. In particular, the REMD study confirms the general picture drawn by Renner et al. that the trans-isomer of the azobenzene peptide exhibits a well-defined structure, while the cis-isomer is a conformational heterogeneous system; that is, the trans-isomer occurs in 2 well-defined conformers, while the cis-isomer represents an energetically frustrated system that leads to an ensemble of conformational structures. Employing a principal component analysis of the REMD data, the free energy landscape of the systems is studied at various temperatures. The implications for the folding and unfolding pathways of the system are discussed.