The mechanism of cis-trans isomerization of prolyl peptides by cyclophilin

Analytic Gradients for the Electrostatic Embedding QM/MM Model in Periodic Boundary Conditions Using Particle-Mesh Ewald Sums and Electrostatic Potential Fitted Charge Operators

Long-range electrostatic effects are fundamental for describing chemical reactivity in the condensed phase. Here, we present the methodology of an efficient quantum mechanical/molecular mechanical (QM/MM) model in periodic boundary conditions (PBC) compatible with QM/MM boundaries at chemical bonds. The method combines electrostatic potential fitted charge operators and electrostatic potentials derived from the smooth particle-mesh Ewald (PME) sum approach. The total energy and its analytic first derivatives with respect to QM, MM, and lattice vectors allow QM/MM molecular dynamics (MD) in the most common thermodynamic ensembles. We demonstrate the robustness of the method by performing a QM/MM MD equilibration of methanol in water. We simulate the cis/trans isomerization free-energy profiles in water of proline amino acid and a proline-containing oligopeptide, showing a correct description of the reaction barrier. Our PBC-compatible QM/MM model can efficiently be used to study the chemical reactivity in the condensed phase and enzymatic catalysis.

Average Electron Density: A Quantitative Tool for Evaluating Non-Classical Bioisosteres of Amides

Bioisosterism is strategically used in drug design to enhance the pharmacokinetic and pharmacodynamic properties of therapeutic molecules. The average electron density (AED) tool has been used in several studies to quantify similarities among nonclassical bioisosteres of carboxylic acid. In this study, the AED tool is used to quantify the similarities among nonclassical bioisosteres of an amide group. In particular, amide-to-1,2,3-triazole bioisosterism is considered. To evaluate the AED differences exhibited by isomers of nonclassical bioisosteres, both isomers of amide (cis and trans) and 1,2,3-triazole (1,4 and 1,5 disubstituted moieity) were considered. The amide and 1,2,3-triazole bioisosteric moieties were capped with various R groups (R= methyl, hydrogen, and chloro) to account for changes in their environment. Amide-to-triazole bioisosteric substitutions were then explored in a more realistic environment, that is, within the FDA-approved anticancer imatinib drug. The AED tool effectively identified similarities between substantially different moieties, 1,2,3-triazole and amide, showing AED differences of no more than 4%. The AED tool was also proven to be useful in evaluating the contribution of various factors affecting triazole-amide bioisosterism including isomerism and changes in their environment. The AED values of each bioisostere were transferable within a maximum difference of 2.6%, irrespective of the change in environment. The 1,4- and 1,5-disubstituted isomers of 1,2,3-triazole have AED values that differ by less than unity, 0.52%. Similarly, the AED values of the cis- and trans-amide isomers differ by only 1.31%. Overall, the AED quantitative tool not only replicated experimental observations regarding similarities in bioisosteres, but also explained and quantified each contributing factor. This demonstrates the extended utility of the AED tool from nonclassical carboxylic acid bioisosteres to amide equivalents.On the contrary, electrostatic potential maps, usually used in the literature to qualitatively evaluate bioisosterism, were not similar for the 1,2,3-triazole and amide bioisosteres, under different environments. Overall, the AED tool proves to be powerful in quantitatively evaluating and predicting bioisosterism across diverse moieties considering structural and environmental variations, making it valuable in drug design.

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.

Dynamic cross-linking of an alginate-acrylamide tough hydrogel system: time-resolved in situ mapping of gel self-assembly

Hydrogels are a popular class of biomaterial that are used in a number of commercial applications (e.g.; contact lenses, drug delivery, and prophylactics). Alginate-based tough hydrogel systems, interpenetrated with acrylamide, reportedly form both ionic and covalent cross-links, giving rise to their remarkable mechanical properties. In this work, we explore the nature, onset and extent of such hybrid bonding interactions between the complementary networks in a model double-network alginate-acrylamide system, using a host of characterisation techniques (e.g.; FTIR, Raman, UV-vis, and fluorescence spectroscopies), in a time-resolved manner. Further, due to the similarity of bonding effects across many such complementary, interpenetrating hydrogel networks, the broad bonding interactions and mechanisms observed during gelation in this model system, are thought to be commonly replicated across alginate-based and broader double-network hydrogels, where both physical and chemical bonding effects are present. Analytical techniques followed real-time bond formation, environmental changes and re-organisational processes that occurred. Experiments broadly identified two phases of reaction; phase I where covalent interaction and physical entanglements predominate, and; phase II where ionic cross-linking effects are dominant. Contrary to past reports, ionic cross-linking occurred more favourably via mannuronate blocks of the alginate chain, initially. Evolution of such bonding interactions was also correlated with the developing tensile and compressive properties. These structure-property findings provide mechanistic insights and future synthetic intervention routes to manipulate the chemo-physico-mechanical properties of dynamically-forming tough hydrogel structures according to need (i.e.; durability, biocompatibility, adhesion, etc.), allowing expansion to a broader range of more physically and/or environmentally demanding biomaterials applications.

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.

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.

Unfolding of CPR3 Gets Initiated at the Active Site and Proceeds via Two Intermediates

Cyclophilin catalyzes the ubiquitous process "peptidyl-prolyl cis-trans isomerization," which plays a key role in protein folding, regulation, and function. Here, we present a detailed characterization of the unfolding of yeast mitochondrial cyclophilin (CPR3) induced by urea. It is seen that CPR3 unfolding is reversible and proceeds via two intermediates, I1 and I2. The I1 state has native-like secondary structure and shows strong anilino-8-naphthalenesulphonate binding due to increased exposure of the solvent-accessible cluster of non-polar groups. Thus, it has some features of a molten globule. The I2 state is more unfolded, but it retains some residual secondary structure, and shows weak anilino-8-naphthalenesulphonate binding. Chemical shift perturbation analysis by ¹H-¹⁵N heteronuclear single quantum coherence spectra reveals disruption of the tertiary contacts among the regions close to the active site in the first step of unfolding, i.e., the N-I1 transition. Both of the intermediates, I1 and I2, showed a propensity to self-associate under stirring conditions, but their kinetic profiles are different; the native protein did not show any such tendency under the same conditions. All these observations could have significant implications for the function of the protein.

Chemotactic Activity of Cyclophilin A in the Skin Mucus of Yellow Catfish (Pelteobagrus fulvidraco) and Its Active Site for Chemotaxis

Fish skin mucus is a dynamic barrier for invading pathogens with a variety of anti-microbial enzymes, including cyclophilin A (CypA), a multi-functional protein with peptidyl-prolyl cis/trans isomerase (PPIase) activity. Beside various other immunological functions, CypA induces leucocytes migration in vitro in teleost. In the current study, we have discovered several novel immune-relevant proteins in yellow catfish skin mucus by mass spectrometry (MS). The CypA present among them was further detected by Western blot. Moreover, the CypA present in the skin mucus displayed strong chemotactic activity for yellow catfish leucocytes. Interestingly, asparagine (like arginine in mammals) at position 69 was the critical site in yellow catfish CypA involved in leucocyte attraction. These novel efforts do not only highlight the enzymatic texture of skin mucus, but signify CypA to be targeted for anti-inflammatory therapeutics.

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.

The TOR signaling pathway regulates vegetative development and virulence in Fusarium graminearum

The target of rapamycin (TOR) signaling pathway plays critical roles in controlling cell growth in a variety of eukaryotes. However, the contribution of this pathway in regulating virulence of plant pathogenic fungi is unknown. We identified and characterized nine genes encoding components of the TOR pathway in Fusarium graminearum. Biological, genetic and biochemical functions of each component were investigated. The FgFkbp12-rapamycin complex binds to the FgTor kinase. The type 2A phosphatases FgPp2A, FgSit4 and FgPpg1 were found to interact with FgTap42, a downstream component of FgTor. Among these, we determined that FgPp2A is likely to be essential for F. graminearum survival, and FgSit4 and FgPpg1 play important roles in cell wall integrity by positively regulating the phosphorylation of FgMgv1, a key MAP kinase in the cell wall integrity pathway. In addition, the FgPpg1 interacting protein, FgTip41, is involved in regulating mycelial growth and virulence. Notably, FgTip41 does not interact with FgTap42 but with FgPpg1, suggesting the existence of FgTap42:FgPpg1:FgTip41 heterotrimer in F. graminearum, a complex not observed in the yeast model. Collectively, we defined a genetic regulatory framework that elucidates how the TOR pathway regulates virulence and vegetative development in F. graminearum.

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.

Characterization of three Arabidopsis thaliana immunophilin genes involved in the plant defense response against Pseudomonas syringae

Plant immunophilins are a broadly conserved family of proteins, which carry out a variety of cellular functions. In this study, we investigated three immunophilin genes involved in the Arabidopsis thaliana response to Pseudomonas syringae infection: a cytoplasmic localized AtCYP19, a cytoplasmic and nuclear localized AtCYP57, and one nucleus directed FKBP known as AtFKBP65. Arabidopsis knock-out mutations in these immunophilins result in an increased susceptibility to P. syringae, whereas overexpression of these genes alters the transcription profile of pathogen-related defense genes and led to enhanced resistance. Histochemical analysis revealed local gene expression of AtCYP19, AtCYP57, and AtFKBP65 in response to pathogen infection. AtCYP19 was shown to be involved in reactive oxygen species production, and both AtCYP57 and AtFKBP65 provided callose accumulation in plant cell wall. Identification of the involvement of these genes in biotic stress response brings a new set of data that will advance plant immune system research and can be widely used for further investigation in this area.

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.

Structural insights into substrate binding by PvFKBP35, a peptidylprolyl cis-trans isomerase from the human malarial parasite Plasmodium vivax

The immunosuppressive drug FK506 binding proteins (FKBPs), an immunophilin family with the immunosuppressive drug FK506 binding property, exhibit peptidylprolyl cis-trans isomerase (PPIase) activity. While the cyclophilin-catalyzed peptidylprolyl isomerization of X-Pro peptide bonds has been extensively studied, the mechanism of the FKBP-mediated peptidylprolyl isomerization remains uncharacterized. Thus, to investigate the binding of FKBP with its substrate and the underlying catalytic mechanism of the FKBP-mediated proline isomerization, here we employed the FK506 binding domain (FKBD) of the human malarial parasite Plasmodium vivax FK506 binding protein 35 (PvFKBP35) and examined the details of the molecular interaction between the isomerase and a peptide substrate. The crystallographic structures of apo PvFKBD35 and its complex with the tetrapeptide substrate succinyl-Ala-Leu-Pro-Phe-p-nitroanilide (sALPFp) determined at 1.4 Å and 1.65 Å resolutions, respectively, showed that the substrate binds to PvFKBD35 in a cis conformation. Nuclear magnetic resonance (NMR) studies demonstrated the chemical shift perturbations of D55, H67, V73, and I74 residues upon the substrate binding. In addition, the X-ray crystal structure, along with the mutational studies, shows that Y100 is a key residue for the catalytic activity. Taken together, our results provide insights into the catalytic mechanism of PvFKBP35-mediated cis-trans isomerization of substrate and ultimately might aid designing substrate mimetic inhibitors targeting the malarial parasite FKBPs.

Cyclohexyl ketone inhibitors of Pin1 dock in a trans-diaxial cyclohexane conformation

Cyclohexyl ketone substrate analogue inhibitors (Ac–pSer-Ψ[C = OCH]-Pip–tryptamine) of Pin1, the cell cycle regulatory peptidyl-prolyl isomerase (PPIase), were designed and synthesized as potential electrophilic acceptors for the Pin1 active site Cys113 nucleophile to test a proposed nucleophilic addition-isomerization mechanism. Because they were weak inhibitors, models of all three stereoisomers were docked into the active site of Pin1. Each isomer consistently minimized to a trans-diaxial cyclohexane conformation. From this, we hypothesize that Pin1 stretches substrates into a trans-pyrrolidine conformation to lower the barrier to isomerization. Our reduced amide inhibitor of Pin1 adopted a similar trans-pyrrolidine conformation in the crystal structure. The molecular model of 1, which mimics the l-Ser-l-Pro stereochemistry, in the Pin1 active site showed a distance of 4.4 Å, and an angle of 31° between Cys113-S and the ketone carbon. The computational models suggest that the mechanism of Pin1 PPIase is not likely to proceed through nucleophilic addition.

Rapamycin exerts antifungal activity in vitro and in vivo against Mucor circinelloides via FKBP12-dependent inhibition of Tor

The zygomycete Mucor circinelloides is an opportunistic fungal pathogen that commonly infects patients with malignancies, diabetes mellitus, and solid organ transplants. Despite the widespread use of antifungal therapy in the management of zygomycosis, the incidence of infections continues to rise among immunocompromised individuals. In this study, we established that the target and mechanism of antifungal action of the immunosuppressant rapamycin in M. circinelloides are mediated via conserved complexes with FKBP12 and a Tor homolog. We found that spontaneous mutations that disrupted conserved residues in FKBP12 conferred rapamycin and FK506 resistance. Disruption of the FKBP12-encoding gene, fkbA, also conferred rapamycin and FK506 resistance. Expression of M. circinelloides FKBP12 (McFKBP12) complemented a Saccharomyces cerevisiae mutant strain lacking FKBP12 to restore rapamycin sensitivity. Expression of the McTor FKBP12-rapamycin binding (FRB) domain conferred rapamycin resistance in S. cerevisiae, and McFKBP12 interacted in a rapamycin-dependent fashion with the McTor FRB domain in a yeast two-hybrid assay, validating McFKBP12 and McTor as conserved targets of rapamycin. We showed that in vitro, rapamycin exhibited potent growth inhibitory activity against M. circinelloides. In a Galleria mellonella model of systemic mucormycosis, rapamycin improved survival by 50%, suggesting that rapamycin and nonimmunosuppressive analogs have the potential to be developed as novel antifungal therapies for treatment of patients with mucormycosis.

A reduced-amide inhibitor of Pin1 binds in a conformation resembling a twisted-amide transition state

The mechanism of the cell cycle regulatory peptidyl prolyl isomerase (PPIase), Pin1, was investigated using reduced-amide inhibitors designed to mimic the twisted-amide transition state. Inhibitors, R-pSer-Ψ[CH(2)N]-Pro-2-(indol-3-yl)ethylamine, 1 [R = fluorenylmethoxycarbonyl (Fmoc)] and 2 (R = Ac), of Pin1 were synthesized and bioassayed. Inhibitor 1 had an IC(50) value of 6.3 μM, which is 4.5-fold better for Pin1 than our comparable ground-state analogue, a cis-amide alkene isostere-containing inhibitor. The change of Fmoc to Ac in 2 improved aqueous solubility for structural determination and resulted in an IC(50) value of 12 μM. The X-ray structure of the complex of 2 bound to Pin1 was determined to 1.76 Å resolution. The structure revealed that the reduced amide adopted a conformation similar to the proposed twisted-amide transition state of Pin1, with a trans-pyrrolidine conformation of the prolyl ring. A similar conformation of substrate would be destabilized relative to the planar amide conformation. Three additional reduced amides, with Thr replacing Ser and l- or d-pipecolate (Pip) replacing Pro, were slightly weaker inhibitors of Pin1.

A mechanistic and structural analysis of the inhibition of the 90-kDa heat shock protein by the benzoquinone and hydroquinone ansamycins

The benzoquinone ansamycins inhibit the ATPase activity of the 90-kDa heat shock protein (Hsp90), disrupting the function of numerous client proteins involved in oncogenesis. In this study, we examine the role of NAD(P)H:quinone oxidoreductase 1 (NQO1) in the metabolism of trans- and cis-amide isomers of the benzoquinone ansamycins and their mechanism of Hsp90 inhibition. Inhibition of purified human Hsp90 by a series of benzoquinone ansamycins was examined in the presence and absence of NQO1, and their relative rate of NQO1-mediated reduction was determined. Computational-based molecular docking simulations indicated that the trans- but not the cis-amide isomers of the benzoquinone ansamycins could be accommodated by the NQO1 active site, and the ranking order of binding energies correlated with the relative reduction rate using purified human NQO1. The trans-cis isomerization of the benzoquinone ansamycins in Hsp90 inhibition has been disputed in recent reports. Previous computational studies have used the closed or cocrystallized Hsp90 structures in an attempt to explore this isomerization step; however, we have successfully docked both the trans- and cis-amide isomers of the benzoquinone ansamycins into the open Hsp90 structure. The results of these studies indicate that both trans- and cis-amide isomers of the hydroquinone ansamycins exhibited increased binding affinity for Hsp90 relative to their parent quinones. Our data support a mechanism in which trans- rather than cis-amide forms of benzoquinone ansamycins are metabolized by NQO1 to hydroquinone ansamycins and that Hsp90-mediated trans-cis isomerization via tautomerization plays an important role in subsequent Hsp90 inhibition.

Cyclophilin A (CyPA) induces chemotaxis independent of its peptidylprolyl cis-trans isomerase activity: direct binding between CyPA and the ectodomain of CD147

Cyclophilin A (CyPA) is a ubiquitously distributed peptidylprolyl cis-trans isomerase (PPIase) that possesses diverse biological functions. Extracellular CyPA is a potent chemokine, which can directly induce leukocyte chemotaxis and contribute to the pathogenesis of inflammation-mediated diseases. Although it has been identified that the chemotaxis activity of CyPA is mediated through its cell surface signaling receptor CD147, the role of CyPA PPIase activity in this process is disputable, and the underlying molecular mechanism is still poorly understood. In this study, we present the first evidence that CyPA induces leukocyte chemotaxis through a direct binding with the ectodomain of CD147 (CD147(ECT)), independent of its PPIase activity. Although NMR study indicates that the CD147(ECT) binding site on CyPA overlaps with the PPIase active site, the PPIase inactive mutant CyPA(R55A) exhibits similar CD147(ECT) binding ability and chemotaxis activity to those of CyPA(WT). Furthermore, we have identified three key residues of CyPA involved in CD147(ECT) binding and found that mutations H70A, T107A, and R69A result in similar levels of reduction in CD147(ECT) binding ability and chemotaxis activity for CyPA, without affecting the PPIase activity. Our findings indicate that there exists a novel mechanism for CyPA to regulate cellular signaling processes, shedding new light on its applications in drug development and providing a new targeting site for drug design.

Molecular dynamics simulation of β₂-microglobulin in denaturing and stabilizing conditions

β₂-Microglobulin has been a model system for the study of fibril formation for 20 years. The experimental study of β₂-microglobulin structure, dynamics, and thermodynamics in solution, at atomic detail, along the pathway leading to fibril formation is difficult because the onset of disorder and aggregation prevents signal resolution in Nuclear Magnetic Resonance experiments. Moreover, it is difficult to characterize conformers in exchange equilibrium. To gain insight (at atomic level) on processes for which experimental information is available at molecular or supramolecular level, molecular dynamics simulations have been widely used in the last decade. Here, we use molecular dynamics to address three key aspects of β₂-microglobulin, which are known to be relevant to amyloid formation: (1) 60 ns molecular dynamics simulations of β₂-microglobulin in trifluoroethanol and in conditions mimicking low pH are used to study the behavior of the protein in environmental conditions that are able to trigger amyloid formation; (2) adaptive biasing force molecular dynamics simulation is used to force cis-trans isomerization at Proline 32 and to calculate the relative free energy in the folded and unfolded state. The native-like trans-conformer (known as intermediate 2 and determining the slow phase of refolding), is simulated for 10 ns, detailing the possible link between cis-trans isomerization and conformational disorder; (3) molecular dynamics simulation of highly concentrated doxycycline (a molecule able to suppress fibril formation) in the presence of β₂-microglobulin provides details of the binding modes of the drug and a rationale for its effect.

Amide cis-trans isomerization in aqueous solutions of methyl N-formyl-D-glucosaminides and methyl N-acetyl-D-glucosaminides: chemical equilibria and exchange kinetics

Amide cis-trans isomerization (CTI) in methyl 2-deoxy-2-acylamido-d-glucopyranosides was investigated by (1)H and (13)C NMR spectroscopy. Singly (13)C-labeled methyl 2-deoxy-2-formamido-d-glucopyranoside (MeGlcNFm) anomers provided standard (1)H and (13)C chemical shifts and (1)H-(1)H and (13)C-(13)C spin-coupling constants for cis and trans amides that are detected readily in aqueous solution. Equipped with this information, doubly (13)C-labeled methyl 2-deoxy-2-acetamido-d-glucopyranoside (MeGlcNAc) anomers were investigated, leading to the detection and quantification of cis and trans amides in this biologically important aminosugar. In comparison to MeGlcNFm anomers, the percentage of cis amide in aqueous solutions of MeGlcNAc anomers is small ( approximately 23% for MeGlcNFm versus approximately 1.8% for MeGlcNAc at 42 degrees C) but nevertheless observable with assistance from (13)C-labeling. Temperature studies gave thermodynamic parameters DeltaG degrees , DeltaH degrees , and DeltaS degrees for cis-trans interconversion in MeGlcNFm and MeGlcNAc anomers. Cis/trans equilibria depended on anomeric configuration, with solutions of alpha-anomers containing less cis amide than those of beta-anomers. Confirmation of the presence of cis amide in MeGlcNAc solutions derived from quantitative (13)C saturation transfer measurements of CTI rate constants as a function of solution temperature, yielding activation parameters E(act), DeltaG degrees (), DeltaH degrees (), and DeltaS degrees () for saccharide CTI. Rate constants for the conversion of trans to cis amide in MeGlcNFm and MeGlcNAc anomers ranged from 0.02 to 3.59 s(-1) over 31-85 degrees C, compared to 0.24-80 s(-1) for the conversion of cis to trans amide over the same temperature range. Energies of activation ranged from 16-19 and 19-20 kcal/mol for the cis --> trans and trans --> cis processes, respectively. Complementary DFT calculations on MeGlcNFm and MeGlcNAc model structures were conducted to evaluate the effects of an acyl side chain and anomeric structure, as well as C2-N2 bond rotation, on CTI energetics. These studies show that aqueous solutions of GlcNAc-containing structures contain measurable amounts of both cis and trans amides, which may influence their biological properties.

Identification and comparative analysis of the peptidyl-prolyl cis/trans isomerase repertoires of H. sapiens, D. melanogaster, C. elegans, S. cerevisiae and Sz. pombe

The peptidyl-prolyl cis/trans isomerase (PPIase) class of proteins comprises three member families that are found throughout nature and are present in all the major compartments of the cell. Their numbers appear to be linked to the number of genes in their respective genomes, although we have found the human repertoire to be smaller than expected due to a reduced cyclophilin repertoire. We show here that whilst the members of the cyclophilin family (which are predominantly found in the nucleus and cytoplasm) and the parvulin family (which are predominantly nuclear) are largely conserved between different repertoires, the FKBPs (which are predominantly found in the cytoplasm and endoplasmic reticulum) are not. It therefore appears that the cyclophilins and parvulins have evolved to perform conserved functions, while the FKBPs have evolved to fill ever-changing niches within the constantly evolving organisms. Many orthologous subgroups within the different PPIase families appear to have evolved from a distinct common ancestor, whereas others, such as the mitochondrial cyclophilins, appear to have evolved independently of one another. We have also identified a novel parvulin within Drosophila melanogaster that is unique to the fruit fly, indicating a recent evolutionary emergence. Interestingly, the fission yeast repertoire, which contains no unique cyclophilins and parvulins, shares no PPIases solely with the budding yeast but it does share a majority with the higher eukaryotes in this study, unlike the budding yeast. It therefore appears that, in comparison with Schizosaccharomyces pombe, Saccharomyces cerevisiae is a poor representation of the higher eukaryotes for the study of PPIases.

A nonessential role for Arg 55 in cyclophilin18 for catalysis of proline isomerization during protein folding

The protein folding process is often in vitro rate-limited by slow cis-trans proline isomerization steps. Importantly, the rate of this process in vivo is accelerated by prolyl isomerases (PPIases). The archetypal PPIase is the human cyclophilin 18 (Cyp18 or CypA), and Arg 55 has been demonstrated to play a crucial role when studying short peptide substrates in the catalytic action of Cyp18 by stabilizing the transition state of isomerization. However, in this study we show that a R55A mutant of Cyp18 is as efficient as the wild type to accelerate the refolding reaction of human carbonic anhydrase II (HCA II). Thus, it is evident that the active-site located Arg 55 is not required for catalysis of the rate-limiting prolyl cis-trans isomerization steps during the folding of a protein substrate as HCA II. Nevertheless, catalysis of cis-trans proline isomerization in HCA II occurs in the active-site of Cyp18, since binding of the inhibitor cyclosporin A abolishes rate acceleration of the refolding reaction. Obviously, the catalytic mechanisms of Cyp18 can differ when acting upon a simple model peptide, four residues long, with easily accessible Pro residues compared with a large protein molecule undergoing folding with partly or completely buried Pro residues. In the latter case, the isomerization kinetics are significantly slower and simpler mechanistic factors such as desolvation and/or strain might operate during folding-assisted catalysis, since binding to the hydrophobic active site is still a prerequisite for catalysis.

Structural and dynamic implications of an effector-induced backbone amide cis-trans isomerization in cytochrome P450cam

Experimental evidence has been provided for a functionally relevant cis-trans isomerization of the Ile88-Pro89 peptide bond in cytochrome P450(cam) (CYP101). The isomerization is proposed to be a key element of the structural reorganization leading to the catalytically competent form of CYP101 upon binding of the effector protein putidaredoxin (Pdx). A detailed comparison of the results of molecular dynamics simulations on the cis and trans conformations of substrate- and carbonmonoxy-bound ferrous CYP101 with sequence-specific Pdx-induced structural perturbations identified by nuclear magnetic resonance is presented, providing insight into the structural and dynamic consequences of the isomerization. The mechanical coupling between the Pdx binding site on the proximal face of CYP101 and the site of isomerization is described.

Mechanistic insight into the role of transition-state stabilization in cyclophilin A

Peptidyl prolyl cis-trans isomerases (PPIases) are ubiquitous enzymes in biology that catalyze the cis-trans isomerization of the proline imide peptide bond in many cell signaling pathways. The local change of the isomeric state of the prolyl peptide bond acts as a switching mechanism in altering the conformation of proteins. A complete understanding of the mechanism of PPIases is still lacking, and current experimental techniques have not been able to provide a detailed atomistic picture. Here we have carried out several accelerated molecular dynamics simulations with explicit solvent, and we have provided a detailed description of cis-trans isomerization of the free and cyclophilin A-catalyzed process. We show that the catalytic mechanism of cyclophilin is due mainly to the stabilization and preferential binding of the transition state that is achieved by a favorable hydrogen bond interaction with a backbone NH group. We also show that the substrate in the transition state interacts more favorably with the enzyme than the cis isomer, which in turn interacts more favorably than the trans isomer. The stability of the enzyme-substrate complex is directly correlated with the interaction the substrate makes with a highly conserved arginine residue. Finally, we show that catalysis is achieved through the rotation of the carbonyl oxygen on the N-terminal of the prolyl peptide bond in a predominately unidirectional fashion.

Mechanism of action of cyclophilin a explored by metadynamics simulations

Trans/cis prolyl isomerisation is involved in several biological processes, including the development of numerous diseases. In the HIV-1 capsid protein (CA), such a process takes place in the uncoating and recruitment of the virion and is catalyzed by cyclophilin A (CypA). Here, we use metadynamics simulations to investigate the isomerization of CA's model substrate HAGPIA in water and in its target protein CypA. Our results allow us to propose a novel mechanistic hypothesis, which is finally consistent with all of the available molecular biology data.

Peptidyl prolyl isomerases are ubiquitous enzymes whose actions are crucial in several biological processes, such as, for instance, in cellular signalling and in the onset of several diseases, e.g., HIV infection. Therefore, these isomerases are promising targets for the design of new drugs. For this purpose, we need to understand their molecular mechanism of action. One of the most characterized peptidyl prolyl isomerases is cyclophilin A. Previous studies characterized the roles of several protein regions in isomerase function. However, there are still experimentally identified important portions of the protein whose specific actions in the mechanism are still not known. Here, we address this problem by an extensive computational study of cyclophilin A and a substrate peptide that is part of the HIV-1 capside protein. We present a novel four-step mechanism of the whole enzymatic process, which is consistent with all of the available experimental data. Moreover, these steps can be used as targets for the development of drugs, e.g., for HIV-1 infection.

A functional proline switch in cytochrome P450cam

The two-protein complex between putidaredoxin (Pdx) and cytochrome P450(cam) (CYP101) is the catalytically competent species for camphor hydroxylation by CYP101. We detected a conformational change in CYP101 upon binding of Pdx that reorients bound camphor appropriately for hydroxylation. Experimental evidence shows that binding of Pdx converts a single X-proline amide bond in CYP101 from trans or distorted trans to cis. Mutation of proline 89 to isoleucine yields a mixture of both bound camphor orientations, that seen in Pdx-free and that seen in Pdx-bound CYP101. A mutation in CYP101 that destabilizes the cis conformer of the Ile 88-Pro 89 amide bond results in weaker binding of Pdx. This work provides direct experimental evidence for involvement of X-proline isomerization in enzyme function.

Functional drift of sequence attributes in the FK506-binding proteins (FKBPs)

Diverse members of the FK506-binding proteins (FKBPs) group and their complexes with different macrocyclic ligands of fungal origins such as FK506, rapamycin, ascomycin, and their immunosuppressive and nonimmunosuppressive derivatives display a variety of cellular and biological activities. The functional relatedness of the FKBPs was estimated from the following attributes of their aligned sequences: 1 degrees conservation of the consensus sequence; 2 degrees sequence similarity; 3 degrees pI; 4 degrees hydrophobicity; 5 degrees amino acid hydrophobicity and bulkiness profiles. Analyses of the multiple sequence alignments and intramolecular interaction networks calculated from a series of structures of the FKBPs revealed some variations in the interaction clusters formed by the AA residues that are crucial for sustaining peptidylprolyl cis/trans isomerases (PPIases) activity and binding capacity of the FKBPs. Fine diversification of the sequences of the multiple paralogues and orthologues of the FKBPs encoded in different genomes alter the intramolecular interaction patterns of their structures and allowed them to gain some selectivity in binding to diverse targets (functional drift).

Theoretical and experimental investigation of the energetics of cis-trans proline isomerization in peptide models

The energetics of cis-trans proline isomerization in small peptide models have been investigated using the hybrid density functional theory method B3LYP with a 6-31+G* basis set. The molecules studied are models for the phospho-Ser/Thr-Pro substrate for Pin-1, a peptidyl-prolyl isomerase (PPIase) involved in cell division. Pin-1 requires phosphorylation of a Ser or Thr residue adjacent to a Pro residue in the substrate and catalyzes cis-trans isomerization about the proline amide bond. The dihedral angle that would correspond to the reaction coordinate for isomerization of the omega peptide bond was investigated for several small models. Relaxed potential energy scans for this dihedral angle in N-methylacetamide, 1, N,N-dimethylacetamide, 2, acetylpyrrolidine, 3 and acetylproline, 4, were carried out in 20 degrees steps using the B3LYP/6-31+G* level of theory. In addition, similar scans were carried out for 1-4 protonated on the acetylamide carbonyl oxygen. Optimized structures for 1-4 protonated on the amide nitrogen were also obtained at B3LYP/6-31+G*. Relative proton affinities were determined for each site at various angles along the reaction coordinate for isomerization. The relative proton affinities were anchored to experimental gas phase proton affinities, which were taken from the literature for 1 and 2, or determined in an electrospray ionization-quadrupole ion trap instrument using the extended kinetic method for 3 and 4. Proton affinities of 925 +/- 10 and 911 +/- 12 kJ/mol were determined for 3 and 4, respectively. These studies suggest that the nitrogen atom in these amides becomes the most basic site in the molecule at a dihedral angle of ca. 130 degrees . In addition, the nitrogen atoms in 2-4 are predicted to attain basicities in the range 920-950 kJ/mol, making them basic enough to be the preferred site for hydrogen bonding in the Pin-1 active site, in support of the proposed mechanism for PPIases.

Catalytic mechanism of cyclophilin as observed in molecular dynamics simulations: pathway prediction and reconciliation of X-ray crystallographic and NMR solution data

Cyclophilins are proteins that catalyze X-proline cis-trans interconversion, where X represents any amino acid. Its mechanism of action has been investigated over the past years but still generates discussion, especially because until recently structures of the ligand in the cis and trans conformations for the same system were lacking. X-ray crystallographic structures for the complex cyclophilin A and HIV-1 capsid mutants with ligands in the cis and trans conformations suggest a mechanism where the N-terminal portion of the ligand rotates during the cis-trans isomerization. However, a few years before, a C-terminal rotating ligand was proposed to explain NMR solution data. In the present study we use molecular dynamics (MD) simulations to generate a trans structure starting from the cis structure. From simulations starting from the cis and trans structures obtained through the rotational pathways, the seeming contradiction between the two sets of experimental data could be resolved. The simulated N-terminal rotated trans structure shows good agreement with the equivalent crystal structure and, moreover, is consistent with the NMR data. These results illustrate the use of MD simulation at atomic resolution to model structural transitions and to interpret experimental data.

A molecular dynamics study of Cyclophilin A free and in complex with the Ala-Pro dipeptide

Six different molecular dynamics simulations of Cyclophilin A, three with the protein free in water and three with the Ala-Pro dipeptide bound to the protein, have been performed, and analysed with respect to structure and hydration of the active site. The water structure in the binding pocket of the free Cyclophilin A was found to mimic the experimentally obtained binding cis conformation of the dipeptide. Cyclophilin A is a peptidyl-prolyl cis-trans isomerase (PPIase), but the mechanism of the cis/trans isomerization is not exactly clear. This study was performed to understand better the binding between dipeptide and Cyclophilin A, but also two previously proposed isomerization mechanisms are discussed.

The cyclophilin repertoire of the fission yeast Schizosaccharomyces pombe

The cyclophilin repertoire of the fission yeast Schizosaccharomyces pombe is comprised of nine members that are distributed over all three of its chromosomes and range from small single-domain to large multi-domain proteins. Each cyclophilin possesses only a single prolyl-isomerase domain, and these vary in their degree of consensus, including at positions that are likely to affect their drug-binding ability and catalytic activity. The additional identified motifs are involved in putative protein or RNA interactions, while a novel domain that is specific to SpCyp7 and its orthologues may have functions that include an interaction with hnRNPs. The Sz. pombe cyclophilins are found throughout the cell but appear to be absent from the mitochondria, which is unique among the characterized eukaryotic repertoires. SpCyp5, SpCyp6 and SpCyp8 have exhibited significant upregulation of their expression during the meiotic cycle and SpCyp5 has exhibited significant upregulation of its expression during heat stress. All nine have identified members in the repertoires of H. sapiens, D. melanogaster and A. thaliana. However, only three identified members in the cyclophilin repertoire of S. cerevisiae with SpCyp7 identifying a fourth protein that is not a member of the recognized repertoire due to its possession of a degenerate prolyl-isomerase domain. The cyclophilin repertoire of Sz. pombe therefore represents a better model group for the study of cyclophilin function in the higher eukaryotes.

Cis/trans isomerization in HIV-1 capsid protein catalyzed by cyclophilin A: insights from computational and theoretical studies

A network of protein vibrations has recently been identified in the enzyme cyclophilin A (CypA) that is associated with its peptidyl-prolyl cis/trans isomerization activity of small peptide substrates. It has been suggested that this network may have a role in promoting the catalytic step during the isomerization reaction. This work presents the results from the characterization of this network during the isomerization of the Gly89-Pro90 peptide bond in the N-terminal domain of the capsid protein (CA(N)) from human immunodeficiency virus type 1 (HIV-1), which is a naturally occurring, biologically relevant protein substrate for CypA. A variety of computational and theoretical studies are utilized to investigate the protein dynamics of the CypA-CA(N) complex, at multiple time scales, during the isomerization step. The results provide insights into the detailed mechanism of isomerization and confirm the presence of previously reported network of protein vibrations coupled to the reaction. Conserved CypA residues at the complex interface and at positions distal to the interface form parts of this network. There is HIV-1 related medical interest in CypA; incorporation of CypA, complexed with the capsid protein, into the virion is required for the infectious activity of HIV-1. Interaction energy and dynamical cross-correlation calculations are used for a detailed investigation of the protein-protein interactions in the CypA-CA(N) complex. The results show that CA(N) residues His87-Ala-Gly-Pro-Ile-Ala92 form the majority of the interactions with CypA residues. New protein-protein interactions distal to the active site (CypA Arg148-CA(N) Gln95 and CypA Arg148-CA(N) Asn121) are also identified.

The immunophilin-interacting protein AtFIP37 from Arabidopsis is essential for plant development and is involved in trichome endoreduplication

The FKBP12 (FK506-binding protein 12 kD) immunophilin interacts with several protein partners in mammals and is a physiological regulator of the cell cycle. In Arabidopsis, only one specific partner of AtFKBP12, namely AtFIP37 (FKBP12 interacting protein 37 kD), has been identified but its function in plant development is not known. We present here the functional analysis of AtFIP37 in Arabidopsis. Knockout mutants of AtFIP37 show an embryo-lethal phenotype that is caused by a strong delay in endosperm development and embryo arrest. AtFIP37 promoter::beta-glucuronidase reporter gene constructs show that the gene is expressed during embryogenesis and throughout plant development, in undifferentiating cells such as meristem or embryonic cells as well as highly differentiating cells such as trichomes. A translational fusion with the enhanced yellow fluorescent protein indicates that AtFIP37 is a nuclear protein localized in multiple subnuclear foci that show a speckled distribution pattern. Overexpression of AtFIP37 in transgenic lines induces the formation of large trichome cells with up to six branches. These large trichomes have a DNA content up to 256C, implying that these cells have undergone extra rounds of endoreduplication. Altogether, these data show that AtFIP37 is critical for life in Arabidopsis and implies a role for AtFIP37 in the regulation of the cell cycle as shown for FKBP12 and TOR (target of rapamycin) in mammals.

Transition state stabilization by general acid catalysis, water expulsion, and enzyme reorganization in Medicago savita chalcone isomerase

In aqueous solution, Medicago savita chalcone isomerase (CHI) enhances the reaction rate for the unimolecular rearrangement of chalcone (CHN) into flavanone by seven orders of magnitude. Conformations of CHN and their relative free energies in water and CHI were investigated by the thermodynamic perturbation method. In water, CHN adopts two conformations (I and II) with conformation I being higher in energy than conformation II by 3 kcal/mol. Only I can give rise to a near attack conformer (NAC) where the nucleophile O2' and the electrophile C9 are placed in proximity. In CHI, I binds less tightly than II by approximately 2 kcal/mol, resulting in the free energy for NAC formation being approximately 2 kcal/mol higher in the enzyme than in water. This unfavorable feature in the ground state of the CHI reaction requires the predominant catalytic advantage to be taken in the step of NAC --> transition state (TS). From the molecular dynamics simulations of apo-CHI, CHI complexed with CHN (CHI.CHN) and CHI.TS, we found: (i) Lys-97-general-acid catalysis of the O2'(-) nucleophilic addition; (ii) expulsion of three water molecules in the process of TS formation; (iii) release of enzyme structural distortion on TS formation. In the conclusion, CHI's remarkable efficiency of stabilizing the TS and its relatively poor ability in organizing the ground state is compared with chorismate mutase whose catalytic prowess, when compared with water, originates predominantly from the enhanced NAC population at the active site.

Structural insights into the catalytic mechanism of cyclophilin A

Cyclophilins constitute a ubiquitous protein family whose functions include protein folding, transport and signaling. They possess both sequence-specific binding and proline cis-trans isomerase activities, as exemplified by the interaction between cyclophilin A (CypA) and the HIV-1 CA protein. Here, we report crystal structures of CypA in complex with HIV-1 CA protein variants that bind preferentially with the substrate proline residue in either the cis or the trans conformation. Cis- and trans-Pro substrates are accommodated within the enzyme active site by rearrangement of their N-terminal residues and with minimal distortions in the path of the main chain. CypA Arg55 guanidinium group probably facilitates catalysis by anchoring the substrate proline oxygen and stabilizing sp3 hybridization of the proline nitrogen in the transition state.