How Does Pin1 Catalyze the Cis–Trans Prolyl Peptide Bond Isomerization? A QM/MM and Mean Reaction Force Study

The transition state region in nonsynchronous concerted reactions

The critical and vanishing points of the reaction force F(ξ) = -dV(ξ)/dξ yield five important coordinates (ξ_R, ξ_R^* , ξ_TS, ξ_P^* , ξ_P) along the intrinsic reaction coordinate (IRC) for a given concerted reaction or reaction step. These points partition the IRC into three well-defined regions, reactants (ξ_R→ξ_R^* ), transition state (ξ_R^* →ξ_P^* ), and products (ξ_P^* →ξ_P), with traditional roles of mostly structural changes associated with the reactants and products regions and mostly electronic activity associated with the transition state (TS) region. Following the evolution of chemical bonding along the IRC using formal descriptors of synchronicity, reaction electron flux, Wiberg bond orders, and their derivatives (or, more precisely, the intensity of the electron activity) unambiguously indicates that for nonsynchronous reactions, electron activity transcends the TS region and takes place well into the reactants and products regions. Under these circumstances, an extension of the TS region toward the reactants and products regions may occur.

Molecular Insights into the Intrinsic Dynamics and Their Roles During Catalysis in Pin1 Peptidyl-prolyl Isomerase

Proteins are intrinsically dynamic and change conformations over a wide range of time scales. While the conformational dynamics have been realized to be important for protein functions, e.g., in activity-stability trade-offs, how they play a role during enzyme catalysis has been of debate over decades. By studying Pin1 peptidyl-prolyl isomerase using extensive molecular dynamics simulations, here we discuss how the slow intrinsic dynamics of Pin1 observed in the NMR relaxation dispersion experiment occur and couple to isomerization reactions in molecular detail. In particular, we analyze the angular correlation functions of the backbone N-H bonds and find that slow conformational transitions occur at about the 3₁₀ helix in the apo state. These events at the helical region further affect the residues at about the ligand binding site. Unfolding of this helix leads to a tight hydrogen bond between the helical region and the ligand binding loop, thus forming a stable coiled structure. The helical and coiled structures are found to be characteristic of the Pin1-ligand complex with the ligand in the trans and cis states, respectively. These results indicate that the changes in the slow dynamics of Pin1 by the isomerization reaction occur via the shift in populations of the helical and coiled states, where the balance is dependent on the ligand isomerization states.

Dissecting the Dynamics during Enzyme Catalysis: A Case Study of Pin1 Peptidyl-Prolyl Isomerase

Free energy surfaces have played a central role in studying protein conformational changes and enzymatic reactions over decades. Yet, free energy barriers and kinetics are highly dependent on the coordinates chosen to define the surface, and furthermore, the dynamics during the reactions are often overlooked. Our recent study on the Pin1-catalyzed isomerization reaction has indicated that the isomerization transition events remarkably deviate from the free energy path, highlighting the need to understand the reaction dynamics in more detail. To this end, here we investigate the reaction coordinates that describe the transition states of the free energy and transition pathways by minimizing the cross-entropy function. We show that the isomerization transition events can be expressed by the concerted changes in the improper torsion angle ζ and nearby backbone torsional angles of the ligand, whereas the transition state of the free energy surface involves changes in a broad range of coordinates including multiple protein-ligand interactions. The current result supports the previous finding that the isomerization transitions occur quickly from the conformational excited states, which is in sharp contrast to the slow and collective changes suggested from the free energy path. Our results further indicate that the coordinates derived from the transition trajectories are not sufficient for finding the transition states on the free energy surfaces due to the lack of information from conformational excited states.

Conformational Excitation and Nonequilibrium Transition Facilitate Enzymatic Reactions: Application to Pin1 Peptidyl-Prolyl Isomerase

Conformational flexibility of protein is essential for enzyme catalysis. Yet, how protein's conformational rearrangements and dynamics contribute to catalysis remains highly controversial. To unravel protein's role in catalysis, it is inevitable to understand the static and dynamic mechanisms simultaneously. To this end, here the Pin1-catalyzed isomerization reaction is studied from the two perspectives. The static view indicates that the hydrogen bonds involving Pin1 rearrange in a tightly coupled manner with isomerization. In sharp contrast, the isomerization dynamics are found to be very rapid; protein's slow conformational rearrangements thus cannot occur simultaneously with isomerization, and the reaction proceeds in a nonequilibrium manner. The distinctive protein conformations necessary to stabilize the transition state are prepared a priori, i.e., as conformational excited states. The present result suggests that enzymatic reaction is not a simple thermal activation from equilibrium directly to the transition state, thus adding a novel perspective to Pauling's view.

An irreversible inhibitor of peptidyl-prolyl cis/trans isomerase Pin1 and evaluation of cytotoxicity

Pin1 (protein interacting with never in mitosis A-1) is a member of the peptidyl prolyl isomerase (PPIase) family, and catalyzes cis-trans isomerization of pThr/Ser-Pro amide bonds. Because Pin1 is overexpressed in various cancer cell lines and promotes cell growth, it is considered a target for anticancer agents. Here, we designed and synthesized a covalently binding Pin1 inhibitor (S)-2 to target Pin1's active site. This compound inhibited Pin1 in protease-coupled assay, and formed a covalent bond with Cys113 of Pin1, as determined by ESI-MS. The acetoxymethyl ester of (S)-2, i.e., 6, suppressed cyclin D1 expression in human prostate cancer PC-3 cells, and exhibited cytotoxicity. Pin1-knockdown experiments indicated that a target for the cytotoxicity of 6 is Pin1.

Structure and function of the human parvulins Pin1 and Par14/17

Parvulins belong to the family of peptidyl-prolyl cis/trans isomerases (PPIases) assisting in protein folding and in regulating the function of a broad variety of proteins in all branches of life. The human representatives Pin1 and Par14/17 are directly involved in processes influencing cellular maintenance and cell fate decisions such as cell-cycle progression, metabolic pathways and ribosome biogenesis. This review on human parvulins summarizes the current knowledge of these enzymes and intends to oppose the well-studied Pin1 to its less well-examined homolog human Par14/17 with respect to structure, catalytic and cellular function.

Amide rotation trajectories probed by symmetry

Amide rotation of peptidyl-prolyl fragments is an important factor in backbone structure organization of proteins. Computational studies have indicated that this rotation preferentially proceeds through a defined transition-state structure (syn/exo). Here, we complement the computational findings by determining the amide bond rotation barriers for derivatives of the two symmetric proline analogues, meso and racemic pyrrolidine-2,5-dicarboxylic acids. The rotations around these residues represent syn/exo-syn/exo and anti/endo-syn/exo hybrid transition states for the meso and racemic diastereomer, respectively. The rotation barriers are lower for the former rotation by about 9 kJ mol^-1 (aqueous medium), suggesting a strong preference for the syn/exo (clockwise) rotation over the anti/endo (anticlockwise) rotation. The results show that both hybrid rotation processes are enthalpically driven but respond differently to solvent polarity changes due to the different transition state dipole-dipole interactions.

Initiation of prolyl cis-trans isomerisation in the CDR-H3 loop of an antibody in response to antigen binding

Proline cis-trans isomerisation is a regulatory mechanism used in a range of biological processes, and is related to various diseases such as Alzheimers disease and cancer. However, the details of the exact molecular mechanism by which it occurs are not known. Using X-ray crystallography, proline isomerisation has been shown to occur following formation of an antigen-antibody complex between the target epiregulin (EPR) and the antibody 9E5, at proline (Pro103), located in the third complementarity-determining region (CDR) of the heavy chain of 9E5. To obtain an accurate description of the pathway involved in cis-trans isomerisation in this system, we performed ten independent long molecular dynamics (MD) simulations starting at a stable transient bound structure obtained from many short binding MD simulations. As a result, we were able to describe the process by which cis-trans isomerisation is initiated, and suggest a catalysis mechanism for cis-trans isomerization in this antigen-antibody system. We found that Asp102, which is immediately adjacent to Pro103, rotates while changing its interacting partner residues in the light chain of 9E5, and at the same time EPR polar residues help to stabilise the intermediate states in the isomerisation process by interacting strongly with Asp102.

The effect of the environment on the methyl transfer reaction mechanism between trimethylsulfonium and phenolate

The methyl transfer reaction mechanism in different molecular environments were studied by electronic structure methods and QM/MM molecular dynamics simulations.

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.

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.

Cysteine-mediated dynamic hydrogen-bonding network in the active site of Pin1

Enzymes catalyze a plethora of chemical reactions that are tightly regulated and intricately coupled in biology. Catalysis of phosphorylation-dependent cis-trans isomerization of peptidyl-prolyl bonds, which act as conformational switches in regulating many post-phosphorylation processes, is considered to be one of the most critical. Pin1 is a cis-trans isomerase of peptidyl-prolyl(ω-) bonds of phosphorylated-Ser/Thr-Pro motifs and has been implicated in many diseases. Structural and experimental studies are still unable to resolve the mechanistic role and protonation states of two adjacent histidines (His59 and His157) and a cysteine (Cys113) in the active site of Pin1. Here, we show that the protonation state of Cys113 mediates a dynamic hydrogen-bonding network in the active site of Pin1, involving the two adjacent histidines and several other residues that are highly conserved and necessary for catalysis. We have used detailed free energy calculations and molecular dynamics simulations, complementing previous experiments, to resolve the ambiguities in the orientations of the histidines and protonation states of these key active site residues, details that are critical for fully understanding the mechanism of Pin1 and necessary for developing potent inhibitors. Importantly, Cys113 is shown to alternate between the unprotonated and neutral states, unprotonated in free Pin1 and neutral in substrate-bound Pin1. Our results are consistent with experiments and provide an explanation for the chemical reactivity of free Pin1 that is suggested to be necessary for the regulation of the enzyme.

Molecular dynamics of the proline switch and its role in Crk signaling

The Crk adaptor proteins play a central role as a molecular timer for the formation of protein complexes including various growth and differentiation factors. The loss of regulation of Crk results in many kinds of cancers. A self-regulatory mechanism for Crk was recently proposed, which involves domain–domain rearrangement. It is initiated by a cis–trans isomerization of a specific proline residue (Pro238 in chicken Crk II) and can be accelerated by Cyclophilin A. To understand how the proline switch controls the autoinhibition at the molecular level, we performed large-scale molecular dynamics and metadynamics simulations in the context of short peptides and multidomain constructs of chicken Crk II. We found that the equilibrium and kinetic properties of the macrostates are regulated not only by the local environments of specified prolines but also by the global organization of multiple domains. We observe the two macrostates (cis closed/autoinhibited and trans open/uninhibited) consistent with NMR experiments and predict barriers. We also propose an intermediate state, the trans closed state, which interestingly was reported to be a prevalent state in human Crk II. The existence of this macrostate suggests that the rate of switching off the autoinhibition by Cyp A may be limited by the relaxation rate of this intermediate state.

The Ess1 prolyl isomerase: traffic cop of the RNA polymerase II transcription cycle

Ess1 is a prolyl isomerase that regulates the structure and function of eukaryotic RNA polymerase II. Ess1 works by catalyzing the cis/trans conversion of pSer5-Pro6 bonds, and to a lesser extent pSer2-Pro3 bonds, within the carboxy-terminal domain (CTD) of Rpb1, the largest subunit of RNA pol II. Ess1 is conserved in organisms ranging from yeast to humans. In budding yeast, Ess1 is essential for growth and is required for efficient transcription initiation and termination, RNA processing, and suppression of cryptic transcription. In mammals, Ess1 (called Pin1) functions in a variety of pathways, including transcription, but it is not essential. Recent work has shown that Ess1 coordinates the binding and release of CTD-binding proteins that function as co-factors in the RNA pol II complex. In this way, Ess1 plays an integral role in writing (and reading) the so-called CTD code to promote production of mature RNA pol II transcripts including non-coding RNAs and mRNAs.

Kinetic isotope effects support the twisted amide mechanism of Pin1 peptidyl-prolyl isomerase

The Pin1 peptidyl-prolyl isomerase catalyzes isomerization of pSer/pThr-Pro motifs in regulating the cell cycle. Peptide substrates, Ac-Phe-Phe-phosphoSer-Pro-Arg-p-nitroaniline, were synthesized in unlabeled form, and with deuterium-labeled Ser-d3 and Pro-d7 amino acids. Kinetic data were collected as a function of Pin1 concentration to measure kinetic isotope effects (KIEs) on catalytic efficiency (kcat/Km). The normal secondary (2°) KIE value measured for the Ser-d3 substrate (kH/kD = 1.6 ± 0.2) indicates that the serine carbonyl does not rehybridize from sp(2) to sp(3) in the rate-determining step, ruling out a nucleophilic addition mechanism. The normal 2° KIE can be explained by hyperconjugation between Ser α-C-H/D and C═O and release of steric strain upon rotation of the amide bond from cis to syn-exo. The inverse 2° KIE value (kH/kD = 0.86 ± 0.08) measured for the Pro-d7 substrate indicates rehybridization of the prolyl nitrogen from sp(2) to sp(3) during the rate-limiting step of isomerization. No solvent kinetic isotope was measured by NMR exchange spectroscopy (kH2O/kD2O = 0.92 ± 0.12), indicating little or no involvement of exchangeable protons in the mechanism. These results support the formation of a simple twisted amide transition state as the mechanism for peptidyl prolyl isomerization catalyzed by Pin1. A model of the reaction mechanism is presented using crystal structures of Pin1 with ground state analogues and an inhibitor that resembles a twisted amide transition state.

Pin1 inhibitors: Pitfalls, progress and cellular pharmacology

Compelling data supports the hypothesis that Pin1 inhibitors will be useful for the therapy of cancer: Pin1 deficient mice resist the induction of breast cancers normally evoked by expression of MMTV-driven Ras or Erb2 alleles. While Pin1 poses challenges for drug discovery, several groups have identified potent antagonists by structure based drug design, significant progress has been made designing peptidic inhibitors and a number of natural products have been found that blockade Pin1, notably epigallocatchechin gallate (EGCG), a major flavonoid in green tea. Here we critically discuss the modes of action and likely specificity of these compounds, concluding that a suitable chemical biology tool for probing the function of Pin1 has yet to be found. We conclude by outlining some open questions regarding the target validation of Pin1 and the prospects for identification of improved inhibitors in the future.

Cis-trans isomerization of omega dihedrals in proteins

Peptide bonds in protein structures are mainly found in trans conformation with a torsion angle ω close to 180°. Only a very low proportion is observed in cis conformation with ω angle around 0°. Cis-trans isomerization leads to local conformation changes which play an important role in many biological processes. In this paper, we reviewed the recent discoveries and research achievements in this field. First, we presented some interesting cases of biological processes in which cis-trans isomerization is directly implicated. It is involved in protein folding and various aspect of protein function like dimerization interfaces, autoinhibition control, channel gating, membrane binding. Then we reviewed conservation studies of cis peptide bonds which emphasized evolution constraints in term of sequence and local conformation. Finally we made an overview of the numerous molecular dynamics studies and prediction methodologies already developed to take into account this structural feature in the research area of protein modeling. Many cis peptide bonds have not been recognized as such due to the limited resolution of the data and to the refinement protocol used. Cis-trans proline isomerization reactions represents a vast and promising research area that still needs to be further explored for a better understanding of isomerization mechanism and improvement of cis peptide bond predictions.