Two pathways mediate interdomain allosteric regulation in pin1

Dynamic Allostery Modulates Catalytic Activity by Modifying the Hydrogen Bonding Network in the Catalytic Site of Human Pin1

Allosteric communication among domains in modular proteins consisting of flexibly linked domains with complimentary roles remains poorly understood. To understand how complementary domains communicate, we have studied human Pin1, a representative modular protein with two domains mutually tethered by a flexible linker: a WW domain for substrate recognition and a peptidyl-prolyl isomerase (PPIase) domain. Previous studies of Pin1 showed that physical contact between the domains causes dynamic allostery by reducing conformation dynamics in the catalytic domain, which compensates for the entropy costs of substrate binding to the catalytic site and thus increases catalytic activity. In this study, the S138A mutant PPIase domain, a mutation that mimics the structural impact of the interdomain contact, was demonstrated to display dynamic allostery by rigidification of the α2-α3 loop that harbors the key catalytic residue C113. The reduced dynamics of the α2-α3 loop stabilizes the C113-H59 hydrogen bond in the hydrogen-bonding network of the catalytic site. The stabilized hydrogen bond between C113 and H59 retards initiation of isomerization, which explains the reduced isomerization rate by ~20% caused by the S138A mutation. These results provide new insight into the interdomain allosteric communication of Pin1.

Fine-tuning the extent and dynamics of binding cleft opening as a potential general regulatory mechanism in parvulin-type peptidyl prolyl isomerases

Parvulins or rotamases form a distinct group within peptidyl prolyl cis-trans isomerases. Their exact mode of action as well as the role of conserved residues in the family are still not unambiguously resolved. Using backbone S² order parameters and NOEs as restraints, we have generated dynamic structural ensembles of three distinct parvulins, SaPrsA, TbPin1 and CsPinA. The resulting ensembles are in good agreement with the experimental data but reveal important differences between the three enzymes. The largest difference can be attributed to the extent of the opening of the substrate binding cleft, along which motional mode the three molecules occupy distinct regions. Comparison with a wide range of other available parvulin structures highlights structural divergence along the bottom of the binding cleft acting as a hinge during the opening-closing motion. In the prototype WW-domain containing parvulin, Pin1, this region is also important in forming contacts with the WW domain known to modulate enzymatic activity of the catalytic domain. We hypothesize that modulation of the extent and dynamics of the identified ‘breathing motion’ might be one of the factors responsible for functional differences in the distinct parvulin subfamilies.

Long-Range Signaling in MutS and MSH Homologs via Switching of Dynamic Communication Pathways

Allostery is conformation regulation by propagating a signal from one site to another distal site. This study focuses on the long-range communication in DNA mismatch repair proteins MutS and its homologs where intramolecular signaling has to travel over 70 Å to couple lesion detection to ATPase activity and eventual downstream repair. Using dynamic network analysis based on extensive molecular dynamics simulations, multiple preserved communication pathways were identified that would allow such long-range signaling. The pathways appear to depend on the nucleotides bound to the ATPase domain as well as the type of DNA substrate consistent with previously proposed functional cycles of mismatch recognition and repair initiation by MutS and homologs. A mechanism is proposed where pathways are switched without major conformational rearrangements allowing for efficient long-range signaling and allostery.

Unidirectional allostery in the regulatory subunit RIα facilitates efficient deactivation of protein kinase A

The holoenzyme complex of protein kinase A is in an inactive state; activation involves ordered cAMP binding to two tandem domains of the regulatory subunit and release of the catalytic subunit. Deactivation has been less studied, during which the two cAMPs unbind from the regulatory subunit to allow association of the catalytic subunit to reform the holoenzyme complex. Unbinding of the cAMPs appears ordered as indicated by a large difference in unbinding rates from the two sites, but the cause has remained elusive given the structural similarity of the two tandem domains. Even more intriguingly, NMR data show that allosteric communication between the two domains is unidirectional. Here, we present a mechanism for the unidirectionality, developed from extensive molecular dynamics simulations of the tandem domains in different cAMP-bound forms. Disparate responses to cAMP releases from the two sites (A and B) in conformational flexibility and chemical shift perturbation confirmed unidirectional allosteric communication. Community analysis revealed that the A-site cAMP, by forming across-domain interactions, bridges an essential pathway for interdomain communication. The pathway is impaired when this cAMP is removed but remains intact when only the B-site cAMP is removed. Specifically, removal of the A-site cAMP leads to the separation of the two domains, creating room for binding the catalytic subunit. Moreover, the A-site cAMP, by maintaining interdomain coupling, retards the unbinding of the B-site cAMP and stalls an unproductive pathway of cAMP release. Our work expands the perspective on allostery and implicates functional importance for the directionality of allostery.

Allosteric communication pathways routed by Ca2+/Mg2+ exchange in GCAP1 selectively switch target regulation modes

GCAP1 is a neuronal calcium sensor protein that regulates the phototransduction cascade in vertebrates by switching between activator and inhibitor of the target guanylate cyclase (GC) in a Ca²⁺-dependent manner. We carried out exhaustive molecular dynamics simulations of GCAP1 and determined the intramolecular communication pathways involved in the specific GC activator/inhibitor switch. The switch was found to depend on the Mg²⁺/Ca²⁺ loading states of the three EF hands and on the way the information is transferred from each EF hand to specific residues at the GCAP1/GC interface. Post-translational myristoylation is fundamental to mediate long range allosteric interactions including the EF2-EF4 coupling and the communication between EF4 and the GC binding interface. Some hubs in the identified protein network are the target of retinal dystrophy mutations, suggesting that the lack of complete inhibition of GC observed in many cases is likely due to the perturbation of intra/intermolecular communication routes.

Computational approaches to investigating allostery

Allosteric regulation plays a key role in many biological processes, such as signal transduction, transcriptional regulation, and many more. It is rooted in fundamental thermodynamic and dynamic properties of macromolecular systems that are still poorly understood and are moreover modulated by the cellular context. Here we review the computational approaches used in the investigation of allosteric processes in protein systems. We outline how the models of allostery have evolved from their initial formulation in the sixties to the current views, which more fully account for the roles of the thermodynamic and dynamic properties of the system. We then describe the major classes of computational approaches employed to elucidate the mechanisms of allostery, the insights they have provided, as well as their limitations. We complement this analysis by highlighting the role of computational approaches in promising practical applications, such as the engineering of regulatory modules and identifying allosteric binding sites.

Allosteric activation of SENP1 by SUMO1 β-grasp domain involves a dock-and-coalesce mechanism

Small ubiquitin-related modifiers (SUMOs) are conjugated to proteins to regulate a variety of cellular processes. SENPs are cysteine proteases with a catalytic center located within a channel between two subdomains that catalyzes SUMO C-terminal cleavage for processing of SUMO precursors and de-SUMOylation of target proteins. The β-grasp domain of SUMOs binds to an exosite cleft, and allosterically activates SENPs via an unknown mechanism. Our molecular dynamics simulations showed that binding of the β-grasp domain induces significant conformational and dynamic changes in SENP1, including widening of the exosite cleft and quenching of nanosecond dynamics in all but a distal region. A dock-and-coalesce mechanism emerges for SENP-catalyzed SUMO cleavage: the wedging of the β-grasp domain enables the docking of the proximal portion of the C-terminus and the strengthened cross-channel motional coupling initiates inter-subdomain correlated motions to allow for the distal portion to coalesce around the catalytic center.

DOI:http://dx.doi.org/10.7554/eLife.18249.001

The Dynamic Basis for Signal Propagation in Human Pin1-WW

Allostery is the structural manifestation of information transduction in biomolecules. Its hallmark is conformational change induced by perturbations at a distal site. An increasing body of evidence demonstrates the presence of allostery in very flexible and even disordered proteins, encouraging a thermodynamic description of this phenomenon. Still, resolving such processes at atomic resolution is difficult. Here we establish a protocol to determine atomistic thermodynamic models of such systems using high-resolution solution state nuclear magnetic resonance data and extensive molecular simulations. Using this methodology, we study information transduction in the WW domain of a key cell-cycle regulator Pin1. Pin1 binds promiscuously to phospho-Ser/Thr-Pro motifs, however, disparate structural and dynamic responses have been reported upon binding different ligands. Our model consists of two topologically distinct states whose relative population may be specifically skewed by an incoming ligand. This model provides a canonical basis for the understanding of multi-functionality in Pin1.

All-Atom Molecular Dynamics of Virus Capsids as Drug Targets

Virus capsids are protein shells that package the viral genome. Although their morphology and biological functions can vary markedly, capsids often play critical roles in regulating viral infection pathways. A detailed knowledge of virus capsids, including their dynamic structure, interactions with cellular factors, and the specific roles that they play in the replication cycle, is imperative for the development of antiviral therapeutics. The following Perspective introduces an emerging area of computational biology that focuses on the dynamics of virus capsids and capsid-protein assemblies, with particular emphasis on the effects of small-molecule drug binding on capsid structure, stability, and allosteric pathways. When performed at chemical detail, molecular dynamics simulations can reveal subtle changes in virus capsids induced by drug molecules a fraction of their size. Here, the current challenges of performing all-atom capsid-drug simulations are discussed, along with an outlook on the applicability of virus capsid simulations to reveal novel drug targets.

Rapid Characterization of Allosteric Networks with Ensemble Normal Mode Analysis

Allosteric regulation is a primary means of controlling protein function. By definition, allostery involves the propagation of structural dynamic changes between distal protein sites that yields a functional change. Gaining improved knowledge of these fundamental mechanisms is important for understanding many biomolecular processes and for guiding protein engineering and drug design efforts. In this work we compare and contrast a range of normal mode analysis (NMA) approaches together with network analysis for the prediction of structural dynamics and allosteric sites. Application to heterotrimeric G proteins, hemoglobin, and caspase 7 indicates that atomistic elastic network models provide improved predictions of experimental allosteric mutation sites. Results for G proteins also display an improved consistency with those derived from more computationally demanding MD simulations. Application of this approach across available experimental structures for a given protein family in a unified manner, that we refer to as ensemble NMA, yields the best overall predictive performance. We propose that this atomistic ensemble NMA approach represents an efficient and powerful tool for guiding the exploration of coupled motions and allosteric mechanisms in cases where multiple structures are available and where MD may prove prohibitively expensive.

Emerging Computational Methods for the Rational Discovery of Allosteric Drugs

Allosteric drug development holds promise for delivering medicines that are more selective and less toxic than those that target orthosteric sites. To date, the discovery of allosteric binding sites and lead compounds has been mostly serendipitous, achieved through high-throughput screening. Over the past decade, structural data has become more readily available for larger protein systems and more membrane protein classes (e.g., GPCRs and ion channels), which are common allosteric drug targets. In parallel, improved simulation methods now provide better atomistic understanding of the protein dynamics and cooperative motions that are critical to allosteric mechanisms. As a result of these advances, the field of predictive allosteric drug development is now on the cusp of a new era of rational structure-based computational methods. Here, we review algorithms that predict allosteric sites based on sequence data and molecular dynamics simulations, describe tools that assess the druggability of these pockets, and discuss how Markov state models and topology analyses provide insight into the relationship between protein dynamics and allosteric drug binding. In each section, we first provide an overview of the various method classes before describing relevant algorithms and software packages.

Protein Allostery and Conformational Dynamics

The functions of many proteins are regulated through allostery, whereby effector binding at a distal site changes the functional activity (e.g., substrate binding affinity or catalytic efficiency) at the active site. Most allosteric studies have focused on thermodynamic properties, in particular, substrate binding affinity. Changes in substrate binding affinity by allosteric effectors have generally been thought to be mediated by conformational transitions of the proteins or, alternatively, by changes in the broadness of the free energy basin of the protein conformational state without shifting the basin minimum position. When effector binding changes the free energy landscape of a protein in conformational space, the change affects not only thermodynamic properties but also dynamic properties, including the amplitudes of motions on different time scales and rates of conformational transitions. Here we assess the roles of conformational dynamics in allosteric regulation. Two cases are highlighted where NMR spectroscopy and molecular dynamics simulation have been used as complementary approaches to identify residues possibly involved in allosteric communication. Perspectives on contentious issues, for example, the relationship between picosecond-nanosecond local and microsecond-millisecond conformational exchange dynamics, are presented.

Dynamic Coupling and Allosteric Networks in the α Subunit of Heterotrimeric G Proteins

G protein α subunits cycle between active and inactive conformations to regulate a multitude of intracellular signaling cascades. Important structural transitions occurring during this cycle have been characterized from extensive crystallographic studies. However, the link between observed conformations and the allosteric regulation of binding events at distal sites critical for signaling through G proteins remain unclear. Here we describe molecular dynamics simulations, bioinformatics analysis, and experimental mutagenesis that identifies residues involved in mediating the allosteric coupling of receptor, nucleotide, and helical domain interfaces of Gαi. Most notably, we predict and characterize novel allosteric decoupling mutants, which display enhanced helical domain opening, increased rates of nucleotide exchange, and constitutive activity in the absence of receptor activation. Collectively, our results provide a framework for explaining how binding events and mutations can alter internal dynamic couplings critical for G protein function.

Pre-Anchoring of Pin1 to Unphosphorylated c-Myc in a Fuzzy Complex Regulates c-Myc Activity

Hierarchic phosphorylation and concomitant Pin1-mediated proline isomerization of the oncoprotein c-Myc controls its cellular stability and activity. However, the molecular basis for Pin1 recognition and catalysis of c-Myc and other multisite, disordered substrates in cell regulation and disease is unclear. By nuclear magnetic resonance, surface plasmon resonance, and molecular modeling, we show that Pin1 subdomains jointly pre-anchor unphosphorylated c-Myc1-88 in the Pin1 interdomain cleft in a disordered, or "fuzzy", complex at the herein named Myc Box 0 (MB0) conserved region N-terminal to the highly conserved Myc Box I (MBI). Ser62 phosphorylation in MBI intensifies previously transient MBI-Pin1 interactions in c-Myc1-88 binding, and increasingly engages Pin1PPIase and its catalytic region with maintained MB0 interactions. In cellular assays, MB0 mutated c-Myc shows decreased Pin1 interaction, increased protein half-life, but lowered rates of Myc-driven transcription and cell proliferation. We propose that dynamic Pin1 recognition of MB0 contributes to the regulation of c-Myc activity in cells.

Negative Regulation of Peptidyl-Prolyl Isomerase Activity by Interdomain Contact in Human Pin1

Pin1 is a modular peptidyl-prolyl isomerase specific for phosphorylated Ser/Thr-Pro (pS/T-P) motifs, typically within intrinsically disordered regions of signaling proteins. Pin1 consists of two flexibly linked domains: an N-terminal WW domain for substrate binding and a larger C-terminal peptidyl-prolyl isomerase (PPIase) domain. Previous studies showed that binding of phosphopeptide substrates to Pin1 could alter Pin1 interdomain contact, strengthening or weakening it depending on the substrate sequence. Thus, substrate-induced changes in interdomain contact may act as a trigger within the Pin1 mechanism. Here, we investigate this possibility via nuclear magnetic resonance studies of several Pin1 mutants. Our findings provide new mechanistic insights for those substrates that reduce interdomain contact. Specifically, the reduced interdomain contact can allosterically enhance PPIase activity relative to that when the contact is sustained. These findings suggest Pin1 interdomain contact can negatively regulate its activity.

Cyclophilin40 isomerase activity is regulated by a temperature-dependent allosteric interaction with Hsp90

Cyclophilin 40 (Cyp40) comprises an N-terminal cyclophilin domain with peptidyl-prolyl isomerase (PPIase) activity and a C-terminal tetratricopeptide repeat (TPR) domain that binds to the C-terminal-EEVD sequence common to both heat shock protein 70 (Hsp70) and Hsp90. We show in the present study that binding of peptides containing the MEEVD motif reduces the PPIase activity by ∼30%. CD and fluorescence assays show that the TPR domain is less stable than the cyclophilin domain and is stabilized by peptide binding. Isothermal titration calorimetry (ITC) shows that the affinity for the-MEEVD peptide is temperature sensitive in the physiological temperature range. Results from these biophysical studies fit with the MD simulations of the apo and holo (peptide-bound) structures which show a significant reduction in root mean square (RMS) fluctuation in both TPR and cyclophilin domains when-MEEVD is bound. The MD simulations of the apo-protein also highlight strong anti-correlated motions between residues around the PPIase-active site and a band of residues running across four of the seven helices in the TPR domain. Peptide binding leads to a distortion in the shape of the active site and a significant reduction in these strongly anti-correlated motions, providing an explanation for the allosteric effect of ligand binding and loss of PPIase activity. Together the experimental and MD results suggest that on heat shock, dissociation of Cyp40 from complexes mediated by the TPR domain leads to an increased pool of free Cyp40 capable of acting as an isomerase/chaperone in conditions of cellular stress.

The Mechanism of ATP-Dependent Allosteric Protection of Akt Kinase Phosphorylation

Kinases use ATP to phosphorylate substrates; recent findings underscore the additional regulatory roles of ATP. Here, we propose a mechanism for allosteric regulation of Akt1 kinase phosphorylation by ATP. Our 4.7-μs molecular dynamics simulations of Akt1 and its mutants in the ATP/ADP bound/unbound states revealed that ATP occupancy of the ATP-binding site stabilizes the closed conformation, allosterically protecting pT308 by restraining phosphatase access and key interconnected residues on the ATP→pT308 allosteric pathway. Following ATP→ADP hydrolysis, pT308 is exposed and readily dephosphorylated. Site-directed mutagenesis validated these predictions and indicated that the mutations do not impair PDK1 and PP2A phosphatase recruitment. We further probed the function of residues around pT308 at the atomic level, and predicted and experimentally confirmed that Akt1(H194R/R273H) double mutant rescues pathology-related Akt1(R273H). Analysis of classical Akt homologs suggests that this mechanism can provide a general model of allosteric kinase regulation by ATP; as such, it offers a potential avenue for allosteric drug discovery.

Investigating Dynamic Interdomain Allostery in Pin1

Signaling proteins often sequester complementary functional sites in separate domains. How do the different domains communicate with one another? An attractive system to address this question is the mitotic regulator, human Pin1 (Lu et al. 1996). Pin-1 consists of two tethered domains: a WW domain for substrate binding, and a catalytic domain for peptidyl-prolyl isomerase (PPIase) activity. Pin1 accelerates the cis-trans isomerization of phospho-Ser/Thr-Pro (pS/T-P) motifs within proteins regulating the cell cycle and neuronal development. The early x-ray (Ranganathan et al. 1997; Verdecia et al. 2000) and solution NMR studies (Bayer et al. 2003; Jacobs et al. 2003) of Pin1 indicated inter- and intradomain motion. We became interested in exploring how such motions might affect interdomain communication, using NMR. Our accumulated results indicate substrate binding to Pin1 WW domain changes the intra/inter domain mobility, thereby altering substrate activity in the distal PPIase domain catalytic site. Thus, Pin1 shows evidence of dynamic allostery, in the sense of Cooper and Dryden (Cooper and Dryden 1984). We highlight our results supporting this conclusion, and summarize them via a simple speculative model of conformational selection.

Mechanism-Based Mathematical Model for Gating of Ionotropic Glutamate Receptors

We present a mathematical model for ionotropic glutamate receptors (iGluR's) that is built on mechanistic understanding and yields a number of thermodynamic and kinetic properties of channel gating. iGluR's are ligand-gated ion channels responsible for the vast majority of fast excitatory neurotransmission in the central nervous system. The effects of agonist-induced closure of the ligand-binding domain (LBD) are transmitted to the transmembrane channel (TMC) via interdomain linkers. Our model demonstrates that, relative to full agonists, partial agonists may reduce either the degree of LBD closure or the curvature of the LBD free energy basin, leading to less stabilization of the channel open state and hence lower channel open probability. A rigorous relation is derived between the channel closed-to-open free energy difference and the tension within the linker. Finally, by treating LBD closure and TMC opening as diffusive motions, we obtain gating trajectories that resemble stochastic current traces from single-channel recordings and calculate the rate constants for transitions between the channel open and closed states. Our model can be implemented by molecular dynamics simulations to realistically depict iGluR gating and may guide functional experiments in gaining deeper insight into this essential family of channel proteins.