Conformational Selection in the Recognition of Phosphorylated Substrates by the Catalytic Domain of Human Pin1

Bispidine as a promising scaffold for designing molecular machines

The development of artificial molecular machines is a challenging endeavor. Herein, we have synthesized a series of bispidine diamides D1-D6 that exhibit rotation reminiscent of a motor motion. Dynamic NMR, X-ray diffraction, quantum mechanical calculations, and molecular dynamics simulations provided insights into their rotational dynamics. All the diamides D1-D6 exhibited mutually independent rotation around the two bispidine arms. However, the rate of rotation and the presence or absence of directionality in amide bond rotation were found to depend on the solvent, temperature, and nature of substitution on the amide carbonyl. These engineered systems may aid in the development of biologically relevant synthetic molecular motors. Studies on homochiral and heterochiral bispidine-peptides revealed that the direction of rotation can be controlled by chirality and the nature of the amino acid.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Proteomics Analysis of SsNsd1-Mediated Compound Appressoria Formation in Sclerotinia sclerotiorum

Sclerotinia sclerotiorum (Lib.) de Bary is a devastating necrotrophic fungal pathogen attacking a broad range of agricultural crops. In this study, although the transcript accumulation of SsNsd1, a GATA-type IVb transcription factor, was much lower during the vegetative hyphae stage, its mutants completely abolished the development of compound appressoria. To further elucidate how SsNsd1 influenced the appressorium formation, we conducted proteomics-based analysis of the wild-type and ΔSsNsd1 mutant by two-dimensional electrophoresis (2-DE). A total number of 43 differentially expressed proteins (≥3-fold change) were observed. Of them, 77% were downregulated, whereas 14% were upregulated. Four protein spots fully disappeared in the mutants. Further, we evaluated these protein sequences by mass spectrometry analysis of the peptide mass and obtained functionally annotated 40 proteins, among which only 17 proteins (38%) were identified to have known functions including energy production, metabolism, protein fate, stress response, cellular organization, and cell growth and division. However, the remaining 23 proteins (56%) were characterized as hypothetical proteins among which four proteins (17%) were predicted to contain the signal peptides. In conclusion, the differentially expressed proteins identified in this study shed light on the ΔSsNsd1 mutant-mediated appressorium deficiency and can be used in future investigations to better understand the signaling mechanisms of SsNsd1 in S. sclerotiorum.

Exploring the molecular basis of RNA recognition by the dimeric RNA-binding protein via molecular simulation methods

RNA-binding protein with multiple splicing (RBPMS) is critical for axon guidance, smooth muscle plasticity, and regulation of cancer cell proliferation and migration. Recently, different states of the RNA-recognition motif (RRM) of RBPMS, one in its free form and another in complex with CAC-containing RNA, were determined by X-ray crystallography. In this article, the free RRM domain, its wild type complex and 2 mutant complex systems are studied by molecular dynamics (MD) simulations. Through comparison of free RRM domain and complex systems, it's found that the RNA binding facilitates stabilizing the RNA-binding interface of RRM domain, especially the C-terminal loop. Although both R38Q and T103A/K104A mutations reduce the binding affinity of RRM domain and RNA, the underlining mechanisms are different. Principal component analysis (PCA) and Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) methods were used to explore the dynamical and recognition mechanisms of RRM domain and RNA. R38Q mutation is positioned on the homodimerization interface and mainly induces the large fluctuations of RRM domains. This mutation does not directly act on the RNA-binding interface, but some interfacial hydrogen bonds are weakened. In contrast, T103A/K104A mutations are located on the RNA-binding interface of RRM domain. These mutations obviously break most of high occupancy hydrogen bonds in the RNA-binding interface. Meanwhile, the key interfacial residues lose their favorable energy contributions upon RNA binding. The ranking of calculated binding energies in 3 complex systems is well consistent with that of experimental binding affinities. These results will be helpful in understanding the RNA recognition mechanisms of RRM domain.

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 C113D mutation in human Pin1 causes allosteric structural changes in the phosphate binding pocket of the PPIase domain through the tug of war in the dual-histidine motif

Pin1 peptidyl-prolyl isomerase (PPIase) catalyzes specifically the pSer/pThr-Pro motif. The cis-trans isomerization mechanism has been studied by various approaches, including X-ray crystallography, site-directed mutagenesis, and the kinetic isotope effect on isomerization. However, a complete picture of the reaction mechanism remains elusive. On the basis of the X-ray structure of Pin1, residue C113 was proposed to play a nucleophile attacker to catalyze the isomerization. The controversial result that the C113D Pin1 mutant retains the activity, albeit at a reduced level, challenges the importance of C113 as a catalyst. To facilitate our understanding of the Pin1 isomerization process, we compared the structures and dynamics of the wild type with those of the C113D mutant Pin1 PPIase domains (residues 51-163). We found the C113D mutation disturbed the hydrogen bonds between the conserved histidine residues, H59 and H157 ("dual-histidine motif"); H59 imidazole forms a stable hydrogen bond to H157 in the wild type, whereas it has a strong hydrogen bond to D113 with weakened bonding to H157 in the C113D mutant. The C113D mutation unbalanced the hydrogen bonding tug of war for H59 between C113/D113 and H157 and destabilized the catalytic site structure, which eventually resulted in an altered conformation of the basic triad (K63, R68, and R69) that binds to the phosphate group in a substrate. The change in the basic triad structure could explain the severely weakened substrate binding ability of the C113D mutant. Overall, this work demonstrated that C113 plays a role in keeping the catalytic site in an active fold, which has never before been described.

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 modelling and simulations in cancer research

The complexity of cancer and the vast amount of experimental data available have made computer-aided approaches necessary. Biomolecular modelling techniques are becoming increasingly easier to use, whereas hardware and software are becoming better and cheaper. Cross-talk between theoretical and experimental scientists dealing with cancer-research from a molecular approach, however, is still uncommon. This is in contrast to other fields, such as amyloid-related diseases, where molecular modelling studies are widely acknowledged. The aim of this review paper is therefore to expose some of the more common approaches in molecular modelling to cancer scientists in simple terms, illustrating success stories while also revealing the limitations of computational studies at the molecular level.

Solution structural analysis of the single-domain parvulin TbPin1

Background

Pin1-type parvulins are phosphorylation-dependent peptidyl-prolyl cis-trans isomerases. Their functions have been widely reported to be involved in a variety of cellular responses or processes, such as cell division, transcription, and apoptosis, as well as in human diseases including Alzheimer's disease and cancers. TbPin1 was identified as a novel class of Pin1-type parvulins from Trypanosoma brucei, containing a unique PPIase domain, which can catalyze the isomerization of phosphorylated Ser/Thr-Pro peptide bond.

Methodology/Principal Findings

We determined the solution structure of TbPin1 and performed ¹⁵N relaxation measurements to analyze its backbone dynamics using multi-dimensional heteronuclear NMR spectroscopy. The average RMSD values of the 20 lowest energy structures are 0.50±0.05 Å for backbone heavy atoms and 0.85±0.08 Å for all heavy atoms. TbPin1 adopts the typical catalytic tertiary structure of Pin1-type parvulins, which comprises a globular fold with a four-stranded anti-parallel β-sheet core surrounded by three α-helices and one 3₁₀-helix. The global structure of TbPin1 is relatively rigid except the active site. The 2D EXSY spectra illustrate that TbPin1 possesses a phosphorylation-dependent PPIase activity. The binding sites of TbPin1 for a phosphorylated peptide substrate {SSYFSG[p]TPLEDDSD} were determined by the chemical shift perturbation approach. Residues Ser15, Arg18, Asn19, Val21, Ser22, Val32, Gly66, Ser67, Met83, Asp105 and Gly107 are involved in substantial contact with the substrate.

Conclusions/Significance

The solution structure of TbPin1 and the binding sites of the phosphorylated peptide substrate on TbPin1 were determined. The work is helpful for further understanding the molecular basis of the substrate specificity for Pin1-type parvulin family and enzyme catalysis.