Modeling Conformational Ensembles of Slow Functional Motions in Pin1-WW

Protein-protein interactions are often mediated by flexible loops that experience conformational dynamics on the microsecond to millisecond time scales. NMR relaxation studies can map these dynamics. However, defining the network of inter-converting conformers that underlie the relaxation data remains generally challenging. Here, we combine NMR relaxation experiments with simulation to visualize networks of inter-converting conformers. We demonstrate our approach with the apo Pin1-WW domain, for which NMR has revealed conformational dynamics of a flexible loop in the millisecond range. We sample and cluster the free energy landscape using Markov State Models (MSM) with major and minor exchange states with high correlation with the NMR relaxation data and low NOE violations. These MSM are hierarchical ensembles of slowly interconverting, metastable macrostates and rapidly interconverting microstates. We found a low population state that consists primarily of holo-like conformations and is a “hub” visited by most pathways between macrostates. These results suggest that conformational equilibria between holo-like and alternative conformers pre-exist in the intrinsic dynamics of apo Pin1-WW. Analysis using MutInf, a mutual information method for quantifying correlated motions, reveals that WW dynamics not only play a role in substrate recognition, but also may help couple the substrate binding site on the WW domain to the one on the catalytic domain. Our work represents an important step towards building networks of inter-converting conformational states and is generally applicable

Molecular Dynamics Simulations Towards The Understanding of the Cis-Trans Isomerization of Proline As A Conformational Switch For The Regulation of Biological Processes

Pin1 is an enzyme central to cell signaling pathways because it catalyzes the cis–trans isomerization of the peptide ω-bond in phosphorylated serine/threonine-proline motifs in many proteins. This regulatory function makes Pin1 a drug target in the treatment of various diseases. The effects of phosphorylation on Pin1 substrates and the basis for Pin1 recognition are not well understood. The conformational consequences of phosphorylation on Pin1 substrate analogues and the mechanism of recognition by the catalytic domain of Pin1 were determined using molecular dynamics simulations. Phosphorylation perturbs the backbone conformational space of Pin1 substrate analogues. It is also shown that Pin1 recognizes specific conformations of its substrate by conformational selection. Dynamical correlated motions in the free Pin1 enzyme are present in the enzyme of the enzyme–substrate complex when the substrate is in the transition state configuration. This suggests that these motions play a significant role during catalysis. These results provide a detailed mechanistic understanding of Pin1 substrate recognition that can be exploited for drug design purposes and further our understanding of the subtleties of post-translational phosphorylation and cis–trans isomerization. Results from accelerated molecular dynamics simulations indicate that catalysis occurs along a restricted path of the backbone configuration of the substrate, selecting specific subpopulations of the conformational space of the substrate in the active site of Pin1. The simulations show that the enzyme–substrate interactions are coupled to the state of the prolyl peptide bond during catalysis. The transition-state configuration of the substrate binds better than the cis and trans states to the catalytic domain of Pin1. This suggests that Pin1 catalyzes its substrate by noncovalently stabilizing the transition state. These results suggest an atomistic detail understanding of the catalytic mechanism of Pin1 that is necessary for the design of novel inhibitors and the treatment of several diseases. Additionally, a set of constant force biased molecular dynamics simulations are presented to explore the kinetic properties of a Pin1 substrate and its unphosphorylated analogue. The simulations indicate that the phosphorylated Pin1 substrate isomerizes slower than the unphosphorylated analogue. This is due to the lower diffusion constant for the phosphorylated Pin1 substrate

Investigating the Roles of the C-terminal Tail and the Metal-dependent C2 Domain in PKCa Regulation

Protein kinase C (PKC) isoenzymes sit in the crossroad of numerous signaling pathways involved in cellular functions such as proliferation, differentiation, migration and survival. The dysregulation of PKCs have been shown to associate with human diseases including cancers, cardiovascular diseases and neurodegenerative diseases. However, the knowledge of PKC regulation is still limited. The objective of this dissertation is to determine the roles of the PKCa C-terminal tail and the C2 regulatory domain in its maturation, activation and down-regulation. The C-terminal V5 domain of PKC contains the least conserved sequence among the isoenzymes. In this study, nuclear magnetic resonance (NMR) and circular dichroism are used to show that the isolated V5 domain is intrinsically disordered in solution. A detailed characterization is provided for the V5 domain’s secondary structural preference, dynamic properties and propensity to interact with a hydrophobic environment. NMR techniques are used to demonstrate that the PKCa C2 domain interacts with the V5 domain. A structural model for the C2–V5 complex is determined. In addition, NMR-detected binding studies reveal that V5 and calcium interact with C2 cooperatively. Mutations that disrupt the C2–V5 interface altered both the conformation of full-length PKCa and the kinetics of membrane translocation. These results indicate that C2-V5 interaction plays an essential role in PKC regulation, through its contribution to both autoinhibition and activation. Furthermore, the V5 domain is shown to directly interact with the peptidyl-prolyl isomerase Pin1. The V5–Pin1 interaction is observed to be highly specific, non-catalytic, and is enhanced by the avidity from bivalent binding. These data provide insights into a novel mechanism for the Pin1-mediated down-regulation process of PKCs. The metal-dependent membrane interactions of the C2 domain are studied with cadmium as structural surrogate and the results compared to other divalent metals. Conformational dynamics change induced by calcium binding is detected for regions connecting to other PKC domains. These results reveal specific roles of calcium ion during membrane interaction and conformational rearrangement of PKC. Together, these data contribute to our understanding of PKC regulation with concerted intramolecular contacts and complex intermolecular interactions, which can help to develop isoenzyme-specific agents to modulate PKC activity

Insights from Coarse-Grained Gō Models for Protein Folding and Dynamics

Exploring the landscape of large scale conformational changes such as protein folding at atomistic detail poses a considerable computational challenge. Coarse-grained representations of the peptide chain have therefore been developed and over the last decade have proved extremely valuable. These include topology-based Gō models, which constitute a smooth and funnel-like approximation to the folding landscape. We review the many variations of the Gō model that have been employed to yield insight into folding mechanisms. Their success has been interpreted as a consequence of the dominant role of the native topology in folding. The role of local contact density in determining protein dynamics is also discussed and is used to explain the ability of Gō-like models to capture sequence effects in folding and elucidate conformational transitions

Computational Perspective on Intricacies of Interactions, Enzyme Dynamics and Solvent Effects in the Catalytic Action of Cyclophilin A

Cyclophilin A (CypA) is the well-studied member of a group of ubiquitous and evolutionarily conserved families of enzymes called peptidyl–prolyl isomerases (PPIases). These enzymes catalyze the cis-trans isomerization of peptidyl-prolyl bond in many proteins. The distinctive functional path triggered by each isomeric state of peptidyl-prolyl bond renders PPIase-catalyzed isomerization a molecular switching mechanism to be used on physiological demand. PPIase activity has been implicated in protein folding, signal transduction, and ion channel gating as well as pathological condition such as cancer, Alzheimer’s, and microbial infections. The more than five order of magnitude speed-up in the rate of peptidyl–prolyl cis–trans isomerization by CypA has been the target of intense research. Normal and accelerated molecular dynamic simulations were carried out to understand the catalytic mechanism of CypA in atomistic details. The results reaffirm transition state stabilization as the main factor in the astonishing enhancement in isomerization rate by enzyme. The ensuing intramolecular polarization, as a result of the loss of pseudo double bond character of the peptide bond at the transition state, was shown to contribute only about −1.0 kcal/mol to stabilizing the transition state. This relatively small contribution demonstrates that routinely used fixed charge classical force fields can reasonably describe these types of biological systems. The computational studies also revealed that the undemanding exchange of the free substrate between β- and α-helical regions is lost in the active site of the enzyme, where it is mainly in the β-region. The resultant relative change in conformational entropy favorably contributes to the free energy of stabilizing the transition state by CypA. The isomerization kinetics is strongly coupled to the enzyme motions while the chemical step and enzyme–substrate dynamics are in turn buckled to solvent fluctuations. The chemical step in the active site of the enzyme is therefore not separated from the fluctuations in the solvent. Of special interest is the nature of catalysis in a more realistic crowded environment, for example, the cell. Enzyme motions in such complicated medium are subjected to different viscosities and hydrodynamic properties, which could have implications for allosteric regulation and function

Hydrogen Exchange Mass Spectrometry for Studying Protein-Ligand Interactions

Hydrogen deuterium exchange (HDX) coupled with mass spectrometry is widely used for probing protein structure and dynamics. Protein-ligand interactions usually induce a reduction in the measured HDX rates an effect that may be ascribed to stabilization of the protein structure. This work aims to improve the general understanding of the changes in HDX patterns associated with ligand binding. We initially applied HDX for studying differences between oxy-hemoglobin (Oxy-Hb) and aquomet-hemoglobin (Chapter 2). The results show that the α and β subunits respond differently to the oxy to aquomet transition with the heme binding pocket being destabilized in both cases. The results suggest that enhanced structural dynamics in the heme binding pocket may have adverse effects on heme-protein interactions. Chapter 3 focuses on the different scenarios that can be encountered in an HDX experiment upon ligand binding. Myoglobin and hemoglobin were used as model systems, focusing on the oxy and deoxy states of both proteins. Our results demonstrate that ligand binding can be stabilizing or destabilizing, leading to decreased or increased HDX rates respectively. In Chapters 4 HDX was used to probe the changes in structural dynamics of caseinolytic protease P (ClpP), an antibiotic drug target, after binding ADEP antibiotics. The mechanism of ADEP binding and the N-terminal structure of ClpP is not well understood with conflicting x-ray structures reported in literature. Our findings demonstrate that the N-terminus of ClpP remains quite unstructured after ADEP binding, while belt region undergoes tightening. Pin 1, a peptidyl prolyl isomerase, binding to a cyclic peptide inhibitor was studied in Chapter 5. Characterization of Pin1-CRYPEVEIC interactions by other techniques has been iv difficult. This study demonstrates that binding of the inhibitor triggers an overall stabilization of Pin 1. We identify a loop that interacts with basic sites of the ligand and that becomes destabilized upon ligand binding. This destabilization is ascribed to steric clashes between the peptide inhibitor and the protei

Application and Optimization of Contact-Guided Replica Exchange Molecular Dynamics

Proteine sind komplexe Makromoleküle, die in lebenden Organismen eine große Vielfalt an wichtigen Aufgaben erfüllen. Proteine können beispielsweise Gene regulieren, Struktur stabilisieren, Zellsignale übertragen, Substanzen transportieren und vieles mehr. Typischerweise sind umfassende Kenntnisse von Struktur und Dynamik eines Proteins erforderlich um dessen physiologische Funktion und Interaktionsmechanismen vollständig zu verstehen. Gewonnene Erkenntnisse sind für Biowissenschaften unerlässlich und können auf viele Bereiche angewendet werden, wie z.B. für Arzneimitteldesign oder zur Krankheitsbehandlung. Trotz des unfassbaren Fortschritts experimenteller Techniken bleibt die Bestimmung einer Proteinstruktur immer noch eine herausfordernde Aufgabe. Außerdem können Experimente nur Teilinformationen liefern und Messdaten können mehrdeutig und schwer zu interpretieren sein. Aus diesem Grund werden häufig Computersimulationen durchgeführt um weitere Erkenntnisse zu liefern und die Lücke zwischen Theorie und Experiment zu schließen. Heute sind viele in-silico Methoden in der Lage genaue Protein Strukturmodelle zu erzeugen, sei es mit einem de novo Ansatz oder durch Verbesserung eines anfänglichen Modells unter Berücksichtigung experimenteller Daten. In dieser Dissertation erforsche ich die Möglichkeiten von Replica Exchange Molekulardynamik (REX MD) als ein physikbasierter Ansatz zur Erzeugung von physikalisch sinnvollen Proteinstrukturen. Dabei lege ich den Fokus darauf möglichst nativähnliche Strukturen zu erhalten und untersuche die Stärken und Schwächen der angewendeten Methode. Ich erweitere die Standardanwendung, indem ich ein kontaktbasiertes Bias-Potential integriere um die Leistung und das Endergebnis von REX zu verbessern. Die Einbeziehung nativer Kontaktpaare, die sowohl aus theoretischen als auch aus experimentellen Quellen abgeleitet werden können, treibt die Simulation in Richtung gewünschter Konformationen und reduziert dementsprechend den notwendigen Rechenaufwand. Während meiner Arbeit führte ich mehrere Studien durch mit dem Ziel, die Anreicherung von nativ-ähnlichen Strukturen zu maximieren, wodurch der End-to-End Prozess von geleitetem REX MD optimiert wird. Jede Studie zielt darauf ab wichtige Aspekte der verwendeten Methode zu untersuchen und zu verbessern: 1) Ich studiere die Auswirkungen verschiedener Auswahlen von Bias-Kontakten, insbesondere die Reichweitenabhängigkeit und den negativen Einfluss von fehlerhaften Kontakten. Dadurch kann ich ermitteln, welche Art von Bias zu einer signifikanten Anreicherung von nativ-ähnlichen Konformationen führen im Vergleich zu regulärem REX. 2) Ich führe eine Parameteroptimierung am verwendeten Bias-Potential durch. Der Vergleich von Ergebnissen aus REX-Simulationen unter Verwendung unterschiedlicher sigmoidförmiger Potentiale weist mir sinnvolle Parameter Bereiche auf, wodurch ich ein ideales Bias-Potenzial für den allgemeinen Anwendungsfall ableiten kann. 3) Ich stelle eine de novo Faltungsmethode vor, die möglichst schnell viele einzigartige Startstrukturen für REX generieren kann. Dabei untersuche ich ausführlich die Leistung dieser Methode und vergleiche zwei verschiedene Ansätze zur Auswahl der Startstruktur. Das Ergebnis von REX wird stark verbessert, falls Strukturen bereits zu Beginn eine große Bandbreite des Konformationsraumes abdecken und gleichzeitig eine geringe Distanz zum angestrebten Zustand aufweisen. 4) Ich untersuche vier komplexe Algorithmusketten, die in der Lage sind repräsentative Strukturen aus großen biomolekularen Ensembles zu extrahieren, welche durch REX erzeugt wurden. Dabei studiere ich ihre Robustheit und Zuverlässigkeit, vergleiche sie miteinander und bewerte ihre erbrachte Leistung numerisch. 5) Basierend auf meiner Erfahrung mit geleitetem REX MD habe ich ein Python-Paket entwickelt um REX-Projekte zu automatisieren und zu vereinfachen. Es ermöglicht einem Benutzer das Entwerfen, Ausführen, Analysieren und Visualisieren eines REX-Projektes in einer interaktiven und benutzerfreundlichen Umgebung

Computational Methods in Biophysics and Medicinal Chemistry: Applications and Challenges

In this thesis I described the theory and application of several computational methods in solving medicinal chemistry and biophysical tasks. I pointed out to the valuable information which could be achieved by means of computer simulations and to the possibility to predict the outcome of traditional experiments. Nowadays, computer represents an invaluable tool for chemists. In particular, the main topics of my research consisted in the development of an automated docking protocol for the voltage-gated hERG potassium channel blockers, and the investigation of the catalytic mechanism of the human peptidyl-prolyl cis-trans isomerase Pin1

A Role for Both Conformational Selection and Induced Fit in Ligand Binding by the LAO Protein

Molecular recognition is determined by the structure and dynamics of both a protein and its ligand, but it is difficult to directly assess the role of each of these players. In this study, we use Markov State Models (MSMs) built from atomistic simulations to elucidate the mechanism by which the Lysine-, Arginine-, Ornithine-binding (LAO) protein binds to its ligand. We show that our model can predict the bound state, binding free energy, and association rate with reasonable accuracy and then use the model to dissect the binding mechanism. In the past, this binding event has often been assumed to occur via an induced fit mechanism because the protein's binding site is completely closed in the bound state, making it impossible for the ligand to enter the binding site after the protein has adopted the closed conformation. More complex mechanisms have also been hypothesized, but these have remained controversial. Here, we are able to directly observe roles for both the conformational selection and induced fit mechanisms in LAO binding. First, the LAO protein tends to form a partially closed encounter complex via conformational selection (that is, the apo protein can sample this state), though the induced fit mechanism can also play a role here. Then, interactions with the ligand can induce a transition to the bound state. Based on these results, we propose that MSMs built from atomistic simulations may be a powerful way of dissecting ligand-binding mechanisms and may eventually facilitate a deeper understanding of allostery as well as the prediction of new protein-ligand interactions, an important step in drug discovery