Efficient and accurate calculation of proline cis/trans isomerization free energies from Hamiltonian replica exchange molecular dynamics simulations

Proline cis/trans isomerization plays an important role in many biological processes but occurs on time scales not accessible to brute-force molecular dynamics (MD) simulations. We have designed a new Hamiltonian replica exchange scheme, ω-bias potential replica exchange molecular dynamics (ωBP-REMD), to efficiently and accurately calculate proline cis/trans isomerization free energies. ωBP-REMD is applied to various proline-containing tripeptides and a biologically important proline residue in the N2-domain of the gene-3-protein of phage fd in the wildtype and mutant variants of the protein. Excellent cis/trans transition rates are obtained. Reweighting of the sampled probability distribution along the peptide bond dihedral angle allows construction of the corresponding free-energy profile and calculation of the cis/trans isomerization free energy with high statistical precision. Very good agreement with experimental data is obtained. ωBP-REMD outperforms standard umbrella sampling in terms of convergence and agreement with experiment and strongly reduces perturbation of the local structure near the proline residue.

Evolutionarily Conserved Proline Residues Impede the Misfolding of the Mouse Prion Protein by Destabilizing an Aggregation-competent Partially Unfolded Form

The misfolding of the prion protein has been linked to several neurodegenerative diseases. Despite extensive studies, the mechanism of the misfolding process remains poorly understood. The present study structurally delineates the role of the conserved proline residues present in the structured C-terminal domain of the mouse prion protein (moPrP) in the misfolding process. It is shown that mutation of these Pro residues to Ala leads to destabilization of the native (N) state, and also to rapid misfolding. Using hydrogen-deuterium exchange (HDX) studies coupled with mass spectrometry (MS), it has been shown that the N state of moPrP is in rapid equilibrium with a partially unfolded form (PUF2*) at pH 4. It has been shown that the Pro to Ala mutations make PUF2* energetically more accessible from the N state by stabilizing it relative to the unfolded (U) state. The apparent rate constant of misfolding is found to be linearly proportional to the extent to which PUF2* is populated in equilibrium with the N state, strongly indicating that misfolding commences from PUF2*. It has also been shown that the Pro residues restrict the boundary of the structural core of the misfolded oligomers. Overall, this study highlights how the conserved proline residues control misfolding of the prion protein by modulating the stability of the partially unfolded form from which misfolding commences.

Protein folding and unfolding: proline cis-trans isomerization at the c subunits of F1 FO -ATPase might open a high conductance ion channel

The c subunits, which constitute the c-ring apparatus of the F1 FO -ATPase, could be the main components of the mitochondrial permeability transition pore (mPTP). The well-known modulator of the mPTP formation and opening is the cyclophilin D (CyPD), a peptidyl-prolyl cis-trans isomerase. On the loop, which connects the two hairpin α-helix of c subunit, is present the unique proline residue (Pro40 ) that could be a biological target of CyPD. Indeed, the proline cis-trans isomerization might provide the switch that interconverts the open/closed states of the pore by pulling out the c-ring lipid plug.

Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

Photopharmacology of Ion Channels through the Light of the Computational Microscope

The optical control and investigation of neuronal activity can be achieved and carried out with photoswitchable ligands. Such compounds are designed in a modular fashion, combining a known ligand of the target protein and a photochromic group, as well as an additional electrophilic group for tethered ligands. Such a design strategy can be optimized by including structural data. In addition to experimental structures, computational methods (such as homology modeling, molecular docking, molecular dynamics and enhanced sampling techniques) can provide structural insights to guide photoswitch design and to understand the observed light-regulated effects. This review discusses the application of such structure-based computational methods to photoswitchable ligands targeting voltage- and ligand-gated ion channels. Structural mapping may help identify residues near the ligand binding pocket amenable for mutagenesis and covalent attachment. Modeling of the target protein in a complex with the photoswitchable ligand can shed light on the different activities of the two photoswitch isomers and the effect of site-directed mutations on photoswitch binding, as well as ion channel subtype selectivity. The examples presented here show how the integration of computational modeling with experimental data can greatly facilitate photoswitchable ligand design and optimization. Recent advances in structural biology, both experimental and computational, are expected to further strengthen this rational photopharmacology approach.

The folding space of proteinβ2-microglobulin is modulated by a single disulfide bridge

Protein beta-2-microglobulin (β2m) is classically considered the causative agent of dialysis related amyloidosis, a conformational disorder that affects patients undergoing long-term hemodialysis. The wild type (WT) form, the ΔN6 structural variant, and the D76N mutant have been extensively used as model systems ofβ2m aggregation. In all of them, the native structure is stabilized by a disulfide bridge between the sulphur atoms of the cysteine residues 25 (at B strand) and 80 (at F strand), which has been considered fundamental inβ2m fibrillogenesis. Here, we use extensive discrete molecular dynamics simulations of a full atomistic structure-based model to explore the role of this disulfide bridge as a modulator of the folding space ofβ2m. In particular, by considering different models for the disulfide bridge, we explore the thermodynamics of the folding transition, and the formation of intermediate states that may have the potential to trigger the aggregation cascade. Our results show that the dissulfide bridge affects folding transition and folding thermodynamics of the considered model systems, although to different extents. In particular, when the interaction between the sulphur atoms is stabilized relative to the other intramolecular interactions, or even locked (i.e. permanently established), the WT form populates an intermediate state featuring a well preserved core and two unstructured termini, which was previously detected only for the D76N mutant. The formation of this intermediate state may have important implications in our understanding ofβ2m fibrillogenesis.