Characterization of the water defect at the HIV-1 gp41 membrane spanning domain in bilayers with and without cholesterol using molecular simulations

The membrane spanning domain (MSD) of human immunodeficiency virus 1 (HIV-1) envelope glycoprotein gp41 is important for fusion and infection. We used molecular dynamics (MD) simulations (3.4 μs total) to relate membrane and peptide properties that lead to water solvation of the α-helical gp41 MSD's midspan arginine in pure dipalmitoylphosphatidylcholine (DPPC) and in 50/50 DPPC/cholesterol membranes. We find that the midspan arginine is solvated by water that penetrates the inner leaflet, leading to a so-called water defect. The water defect is surprisingly robust across initial conditions and membrane compositions, but the presence of cholesterol modulates its behavior in several key ways. In the cholesterol-containing membranes, fluctuations in membrane thickness and water penetration depth are localized near the midspan arginine, and the MSD helices display a tightly regulated tilt angle. In the cholesterol-free membranes, thickness fluctuations are not as strongly correlated to the peptide position and tilt angles vary significantly depending on protein position relative to boundaries between domains of differing thickness. Cholesterol in an HIV-1 viral membrane is required for infection. Therefore, this work suggests that the colocalized water defect and membrane thickness fluctuations in cholesterol-containing viral membranes play an important role in fusion by bringing the membrane closer to a stability limit that must be crossed for fusion to occur.     

Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides

Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl α-d-glucopyranoside and methyl α-d-galactopyranoside, employing the glycam06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the ω (O-5–C-5–C-6–O-6) and χ (C-3–C-4–O-4–HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4–HO-4 and O-6–HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2 kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics.     

Modelling of β-d-glucopyranose ring distortion in different force fields: a metadynamics study

Modelling of carbohydrate conformations is a challenging task for force field developers. Three carbohydrate force fields, namely GLYCAM06, GROMOS 45a4 and OPLS were evaluated. Free energies of different ring conformations of β-d-glucopyranose were calculated using metadynamics in vacuum as well as in explicitly modelled water. All three force fields model the ⁴C₁ conformation as the most stable by at least 6 kJ/mol, as compared to other conformations. Interconversion from the ⁴C₁ to any other conformation is associated with a barrier of no lower than 26 kJ/mol. The free energy surface calculated in the GLYCAM06 force field is in remarkably good agreement with the recent Car-Parrinello metadynamics study. The effect of a water environment is relatively low and analogous in all tested force fields. Namely, the presence of water stabilizes the upper-left (^3,OB) versus bottom-right (B_3,O) area of Stoddard’s plot, relative to the situation in vacuum. Comparison of free and potential surfaces is also provided for vacuum calculations.     

Continuous metadynamics in essential coordinates as a tool for free energy modelling of conformational changes

Modelling of conformational changes in biopolymers is one of the greatest challenges of molecular biophysics. Metadynamics is a recently introduced free energy modelling technique that enhances sampling of configurational (e.g. conformational) space within a molecular dynamics simulation. This enhancement is achieved by the addition of a history-dependent bias potential, which drives the system from previously visited regions. Discontinuous metadynamics in the space of essential dynamics eigenvectors (collective motions) has been proposed and tested in conformational change modelling. Here, we present an implementation of two continuous formulations of metadynamics in the essential subspace. The method was performed in a modified version of the molecular dynamics package GROMACS. These implementations were tested on conformational changes in cyclohexane, alanine dipeptide (terminally blocked alanine, Ace-Ala-Nme) and SH3 domain. The results illustrate that metadynamics in the space of essential coordinates can accurately model free energy surfaces associated with conformational changes. Figure The conformational free energy surface of cyclohexane in the space of the two most intensive collective motions.

---

     

A metadynamic approach to understand the recognition mechanism of the histone H3 tail with the ATRXADD domain

The binding affinity between the histone 3 (H3) tail and the ADD domain of ATRX (ATRX_ADD) increases with the subsequent addition of methyl groups on lysine 9 on H3. To improve our understanding of how the difference in methylation state affects binding between H3 and the ATRX_ADD, we adopted a metadynamic approach to explore the recognition mechanism between the two proteins and identify the key intermolecular interactions that mediate this protein-peptide interaction (PPI). The non-methylated H3 peptide is recognized only by the PHD finger of ATRX_ADD while mono-, di-, and trimethylated H3 is recognized by both the PHD and GATA-like zinc finger of the domain. Furthermore, water molecules play an important role in orienting the lysine 9 anchor towards the GATA-like zinc finger, which results in stabilizing the lysine 9 binding pocket on ATRX_ADD. We compared our computational results against experimentally determined NMR and X-ray structures by demonstrating the RMSD, order parameter S² and hydration number of the complex. The metadynamics data provide new insight into roles of water-bridges and the mechanisms through which K9 hydration stabilizes the H3K9me3:ATRX_ADD PPI, providing context for the high affinity demonstrated between this protein and peptide.     

Modeling catalytic reaction mechanisms in glycoside hydrolases

Modeling catalysis in carbohydrate-active enzymes is a daunting challenge because of the high flexibility and diversity of both enzymes and carbohydrates. Glycoside hydrolases (GHs) are an illustrative example, where conformational changes and subtle interactions have been shown to be critical for catalysis. GHs have pivotal roles in industry (e.g. biofuel or detergent production) and biomedicine (e.g. targets for cancer and diabetes), and thus, a huge effort is devoted to unveil their molecular mechanisms. Besides experimental techniques, computational methods have served to provide an in-depth understanding of GH mechanisms, capturing complex reaction coordinates and the conformational itineraries that substrates follow during the whole catalytic pathway, providing a framework that ultimately may assist the engineering of these enzymes and the design of new inhibitors.     

Molecular modeling study on the dynamical structural features of human smoothened receptor and binding mechanism of antagonist LY2940680 by metadynamics simulation and free energy calculation

Background

The smoothened (SMO) receptor, one of the Class F G protein coupled receptors (GPCRs), is an essential component of the canonical hedgehog signaling pathway which plays a key role in the regulation of embryonic development in animals. The function of the SMO receptor can be modulated by small-molecule agonists and antagonists, some of which are potential antitumour agents. Understanding the binding mode of an antagonist in the SMO receptor is crucial for the rational design of new antitumour agents.

Methods

Molecular dynamics (MD) simulation and dynamical network analysis are used to study the dynamical structural features of SMO receptor. Metadynamics simulation and free energy calculation are employed to explore the binding mechanism between the antagonist and SMO receptor.

Results

The MD simulation results and dynamical network analysis show that the conserved KTXXXW motif in helix VIII has strong interaction with helix I. The α-helical extension of transmembrane 6 (TM6) is detected as part of the ligand-binding pocket and dissociation pathway of the antagonist. The metadynamics simulation results illustrate the binding mechanism of the antagonist in the pocket of SMO receptor, and free energy calculation shows the antagonist needs to overcome about 38 kcal/mol of energy barrier to leave the binding pocket of SMO receptor.

Conclusions

The unusually long TM6 plays an important role on the binding behavior of the antagonist in the pocket of SMO receptor.

General significance

The results can not only profile the binding mechanism between the antagonist and Class F GPCRs, but also supply the useful information for the rational design of a more potential small molecule antagonist bound to SMO receptor.     

Unbinding free energy of acetylcholinesterase bound oxime drugs along the gorge pathway from metadynamics‐umbrella sampling investigation

Because of the pivotal role that the nerve enzyme, acetylcholinesterase plays in terminating nerve impulses at cholinergic synapses. Its active site, located deep inside a 20 Å gorge, is a vulnerable target of the lethal organophosphorus compounds. Potent reactivators of the intoxicated enzyme are nucleophiles, such as bispyridinium oxime that binds to the peripheral anionic site and the active site of the enzyme through suitable cation–π interactions. Atomic scale molecular dynamics and free energy calculations in explicit water are used to study unbinding pathways of two oxime drugs (Ortho‐7 and Obidoxime) from the gorge of the enzyme. The role of enzyme‐drug cation–π interactions are explored with the metadynamics simulation. The metadynamics discovered potential of mean force (PMF) of the unbinding events is refined by the umbrella sampling (US) corrections. The bidimensional free energy landscape of the metadynamics runs are further subjected to finite temperature string analysis to obtain the transition tube connecting the minima and bottlenecks of the unbinding pathway. The PMF is also obtained from US simulations using the biasing potential constructed from the transition tube and are found to be consistent with the metadynamics‐US corrected results. Although experimental structural data clearly shows analogous coordination of the two drugs inside the gorge in the bound state, the PMF of the drug trafficking along the gorge pathway point, within an equilibrium free energy context, to a multistep process that differs from one another. Routes, milestones and subtlety toward the unbinding pathway of the two oximes at finite temperature are identified. Proteins 2014; 82:1799–1818. © 2014 Wiley Periodicals, Inc.     

Uracil/H+ Symport by FurE Refines Aspects of the Rocking-bundle Mechanism of APC-type Transporters

Transporters mediate the uptake of solutes, metabolites and drugs across the cell membrane. The eukaryotic FurE nucleobase/H⁺ symporter of Aspergillus nidulans has been used as a model protein to address structure–function relationships in the APC transporter superfamily, members of which are characterized by the LeuT-fold and seem to operate by the so-called ‘rocking-bundle’ mechanism. In this study, we reveal the binding mode, translocation and release pathway of uracil/H⁺ by FurE using path collective variable, funnel metadynamics and rational mutational analysis. Our study reveals a stepwise, induced-fit, mechanism of ordered sequential transport of proton and uracil, which in turn suggests that FurE, functions as a multi-step gated pore, rather than employing ‘rocking’ of compact domains, as often proposed for APC transporters. Finally, our work supports that specific residues of the cytoplasmic N-tail are involved in substrate translocation, in line with their essentiality for FurE function.