Entropic and Surprisingly Small Intramolecular Polarization Effects in the Mechanism of Cyclophilin A

Dissecting the Allosteric Fine-Tuning of Enzyme Catalysis

Fully understanding the mechanism of allosteric regulation in biomolecules requires separating and examining all of the involved factors. In enzyme catalysis, allosteric effector binding shifts the structure and dynamics of the active site, leading to modified energetic (e.g., energy barrier) and dynamical (e.g., diffusion coefficient) factors underlying the catalyzed reaction rate. Such modifications can be subtle and dependent on the type of allosteric effector, representing a fine-tuning of protein function. The microscopic description of allosteric regulation at the level of function-dictating factors has prospective applications in fundamental and pharmaceutical sciences, which is, however, largely missing so far. Here, we characterize the allosteric fine-tuning of enzyme catalysis, using human Pin1 as an example, by performing more than half-millisecond all-atom molecular dynamics simulations. Changes of reaction kinetics and the dictating factors, including the free energy surface along the reaction coordinate and the diffusion coefficient of the reaction dynamics, under various enzyme and allosteric effector binding conditions are examined. Our results suggest equal importance of the energetic and dynamical factors, both of which can be modulated allosterically, and the combined effect determines the final allosteric output. We also reveal the potential dynamic basis for allosteric modulation using an advanced statistical technique to detect function-related conformational dynamics. Methods developed in this work can be applied to other allosteric systems.

Conserved Conformational Dynamics Reveal a Key Dynamic Residue in the Gatekeeper Loop of Human Cyclophilins

Cyclophilins are ubiquitous human enzymes that catalyze peptidyl-prolyl cis-trans isomerization in protein substrates. Of the 17 unique isoforms, five closely related isoforms (CypA-E) are found in various environments and participate in diverse cellular processes, yet all have similar structures and the same core catalytic function. The question is what key residues are behind the conserved function of these enzymes. Here, conformational dynamics are compared across these isoforms to detect conserved dynamics essential for the catalytic activity of cyclophilins. A set of key dynamic residues, defined by the most dynamically conserved positions, are identified in the gatekeeper 2 region. The highly conserved glycine (Gly80) in this region is predicted to underlie the local flexibility, which is further tested by molecular dynamics simulations performed on mutants (G80A) of CypE and CypA. The mutation leads to decreased flexibility of CypE and CypA during substrate binding but increased flexibility during catalysis. Dynamical changes occur in the mutated region and a distal loop downstream of the mutation site in sequence. Examinations of the mutational effect on catalysis show that both mutated CypE and CypA exhibit shifted binding free energies of the substrate under distinct isomer conformations. The results suggest a loss of function in the mutated CypE and CypA. These catalytic changes by the mutation are likely independent of the substrate sequence, at least in CypA. Our work presents a method to identify function-related key residues in proteins.

Post-translational Modifications of Cyclophilin D Fine-Tune Its Conformational Dynamics and Activity: Implications for Its Mitochondrial Function

Mitochondria are the powerhouse of a cell, whose disruption due to mitochondrial pore opening can cause cell death, leading to necrosis and many other diseases. The peptidyl-prolyl cis-trans isomerase cyclophilin D (CypD) is a key player in the regulation of the mitochondrial pore. The activity of CypD can be modulated by the post-translational modification (PTM). However, the detailed mechanism of this functional modulation is not well understood. Here, we investigate the catalytic mechanism of unmodified and modified CypD by calculating the reaction free energy profiles and characterizing the function-related conformational dynamics using molecular dynamics simulations and associated analyses. Our results show that unmodified and modified CypD considerably lower the isomerization free energy barrier compared to a free peptide substrate, supporting the catalytic activity of CypD in the simulation systems. The unmodified CypD reduces the free energy difference between the cis and trans states of the peptide substrate, suggesting a stronger binding affinity of CypD toward cis, consistent with experiments. In contrast, phosphorylated CypD further stabilizes trans, leading to a lower catalytic rate in the trans-to-cis direction. The differential catalytic activities of the unmodified and phosphorylated CypD are due to a significant shift of the conformational ensemble upon phosphorylation under different functional states. Interestingly, the local flexibility is both reduced and enhanced at distinct regions by phosphorylation, which is explained by a "seesaw" model of flexibility modulation. The allosteric pathway between the phosphorylation site and a distal site displaying substantial conformational changes upon phosphorylation is also identified, which is influenced by the presence of the substrate or the substrate conformation. Similar conclusions are obtained for the acetylation of CypD using the same peptide substrate and the influence of substrate sequence is also examined. Our work may serve as the basis for the understanding of other PTMs and PTM-initiated allosteric regulations in CypD.

p53 Is Potentially Regulated by Cyclophilin D in the Triple-Proline Loop of the DNA Binding Domain

The multifunctional protein p53 is the central molecular sensor of cellular stresses. The canonical function of p53 is to transcriptionally activate target genes in response to, for example, DNA damage that may trigger apoptosis. Recently, p53 was also found to play a role in the regulation of necrosis, another type of cell death featured by the mitochondrial permeability transition (mPT). In this process, p53 directly interacts with the mPT regulator cyclophilin D, the detailed mechanism of which however remains poorly understood. Here, we report a comprehensive computational investigation of the p53-cyclophilin D interaction using molecular dynamics simulations and associated analyses. We have identified the specific cyclophilin D binding site on p53 that is located at proline 151 in the DNA binding domain. As a peptidyl-prolyl isomerase, cyclophilin D binds p53 and catalyzes the cis-trans isomerization of the peptide bond preceding proline 151. We have also characterized the effect of such an isomerization and found that the p53 domain in the cis state is overall more rigid than the trans state except for the local region around proline 151. Dynamical changes upon isomerization occur in both local and distal regions, indicating an allosteric effect elicited by the isomerization. We present potential allosteric communication pathways between proline 151 and distal sites, including the DNA binding surface. Our work provides, for the first time, a model for how cyclophilin D binds p53 and regulates its activity by switching the configuration of a specific site.

Computational perspective and evaluation of plausible catalytic mechanisms of peptidyl-prolyl cis-trans isomerases

BACKGROUND

Peptidyl prolyl cis-trans isomerization of the protein backbone is involved in the regulation of many biological processes. Cis-trans isomerization is notoriously slow and is catalyzed by a family of cis-trans peptidyl prolyl isomerases (PPIases) that have been implicated in many diseases. A general consensus on how these enzymes speed up prolyl isomerization has not been reached after decades of both experimental and computational studies.

SCOPE OF REVIEW

Computational studies carried out to understand the catalytic mechanism of the prototypical FK506 binding protein 12, Cyclophilin A and peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) are reviewed. A summary and an evaluation of the implications of the proposed mechanisms from computational studies are presented.

MAJOR CONCLUSIONS

The analysis of computational studies and evaluation of the proposed mechanisms provide a general consensus and a better understanding of PPIase catalysis. The speedup of the rate of peptidyl-prolyl isomerization by PPIases can be best described by a catalytic mechanism in which the substrate in transition state configuration is stabilized. The enzymes preferentially bind the transition state configuration of the substrate relative to the cis conformation, which in most cases is bound better than the trans conformation of the substrate. Stabilization of the transition state configuration of the substrate leads to a lower free energy barrier and a faster rate of isomerization when compared to the uncatalyzed isomerization reaction.

GENERAL SIGNIFICANCE

Fully understanding the catalytic mechanism of PPIases has broad implications for drug design, elucidation of the molecular basis of many diseases, protein engineering, and enzyme catalysis in general. This article is part of a Special Issue entitled Proline-directed Foldases: Cell Signaling Catalysts and Drug Targets.

Towards fast, rigorous and efficient conformational sampling of biomolecules: Advances in accelerated molecular dynamics

BACKGROUND

Accelerated molecular dynamics (aMD) has been proven to be a powerful biasing method for enhanced sampling of biomolecular conformations on general-purpose computational platforms. Biologically important long timescale events that are beyond the reach of standard molecular dynamics can be accessed without losing the detailed atomistic description of the system in aMD. Over other biasing methods, aMD offers the advantages of tuning the level of acceleration to access the desired timescale without any advance knowledge of the reaction coordinate.

SCOPE OF REVIEW

Recent advances in the implementation of aMD and its applications to small peptides and biological macromolecules are reviewed here along with a brief account of all the aMD variants introduced in the last decade.

MAJOR CONCLUSIONS

In comparison to the original implementation of aMD, the recent variant in which all the rotatable dihedral angles are accelerated (RaMD) exhibits faster convergence rates and significant improvement in statistical accuracy of retrieved thermodynamic properties. RaMD in conjunction with accelerating diffusive degrees of freedom, i.e. dual boosting, has been rigorously tested for the most difficult conformational sampling problem, protein folding. It has been shown that RaMD with dual boosting is capable of efficiently sampling multiple folding and unfolding events in small fast folding proteins.

GENERAL SIGNIFICANCE

RaMD with the dual boost approach opens exciting possibilities for sampling multiple timescales in biomolecules. While equilibrium properties can be recovered satisfactorily from aMD-based methods, directly obtaining dynamics and kinetic rates for larger systems presents a future challenge. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Mechanistic insights into Pin1 peptidyl-prolyl cis-trans isomerization from umbrella sampling simulations

The peptidyl-proyl isomerase Pin1 plays a key role in the regulation of phospho(p)-Ser/Thr-Pro proteins, acting as a molecular timer of the cell cycle. After recognition of these motifs, Pin1 catalyzes the rapid cis-trans isomerization of proline amide bonds of substrates, contributing to maintain the equilibrium between the two conformations. Although a great interest has arisen on this enzyme, its catalytic mechanism has long been debated. Here, the cis-trans isomerization of a model peptide system was investigated by means of umbrella sampling simulations in the Pin1-bound and unbound states. We obtained free energy barriers consistent with experimental data, and identified several enzymatic features directly linked to the acceleration of the prolyl bond isomerization. In particular, an enhanced autocatalysis, the stabilization of perturbed ground state conformations, and the substrate binding in a procatalytic conformation were found as main contributions to explain the lowering of the isomerization free energy barrier.