A tool set to map allosteric networks through the NMR chemical shift covariance analysis

NMR tools to detect protein allostery

Allostery is a fundamental mechanism of cellular homeostasis by intra-protein communication between distinct functional sites. It is an internal process of proteins to steer interactions not only with each other but also with other biomolecules such as ligands, lipids, and nucleic acids. In addition, allosteric regulation is particularly important in enzymatic activities. A major challenge in structural and molecular biology today is unraveling allosteric sites in proteins, to elucidate the detailed mechanism of allostery and the development of allosteric drugs. Here we summarize the recently developed tools and approaches which enable the elucidation of regulatory hotspots and correlated motion in biomolecules, focusing primarily on solution-state nuclear magnetic resonance spectroscopy (NMR). These tools open an avenue towards a rational understanding of the mechanism of allostery and provide essential information for the design of allosteric drugs.

Allostery, engineering and inhibition of tryptophan synthase

The final two steps of tryptophan biosynthesis are catalyzed by the enzyme tryptophan synthase (TS), composed of alpha (αTS) and beta (βTS) subunits. Recently, experimental and computational methods have mapped "allosteric networks" that connect the αTS and βTS active sites. In αTS, allosteric networks change across the catalytic cycle, which might help drive the conformational changes associated with its function. Directed evolution studies to increase catalytic function and expand the substrate profile of stand-alone βTS have also revealed the importance of αTS in modulating the conformational changes in βTS. These studies also serve as a foundation for the development of TS inhibitors, which can find utility against Mycobacterium tuberculosis and other bacterial pathogens.

Toward a molecular mechanism for the interaction of ATP with alpha-synuclein

The ability of Adenosine Triphosphate (ATP) to modulate protein solubility establishes a critical link between ATP homeostasis and proteinopathies, such as Parkinson's (PD). The most significant risk factor for PD is aging, and ATP levels decline dramatically with age. However, the mechanism by which ATP interacts with alpha-synuclein (αS), whose aggregation is characteristic of PD, is currently not fully understood, as is ATP's effect on αS aggregation. Here, we use nuclear magnetic resonance spectroscopy as well as fluorescence, dynamic light scattering and microscopy to show that ATP affects multiple species in the αS self-association cascade. The triphosphate moiety of ATP disrupts long-range electrostatic intramolecular contacts in αS monomers to enhance initial aggregation, while also inhibiting the formation of late-stage β-sheet fibrils by disrupting monomer-fibril interactions. These effects are modulated by magnesium ions and early onset PD-related αS mutations, suggesting that loss of the ATP hydrotropic function on αS fibrillization may play a role in PD etiology.

Analysis of coordinated NMR chemical shifts to map allosteric regulatory networks in proteins

The exquisite sensitivity of the NMR chemical shift to local environment makes it an ideal probe to assess atomic level perturbations in proteins of all sizes and structural compositions. Recent advances in solution and solid-state NMR spectroscopy of biomolecules have leveraged the chemical shift to report on short- and long-range couplings between individual amino acids to establish "networks" of residues that form the basis of allosteric pathways that transmit chemical signals through the protein matrix to induce functional responses. The simple premise that thermodynamically and functionally coupled regions of a protein (i.e. active and allosteric sites) should be reciprocally sensitive to structural or dynamic perturbations has enabled NMR spectroscopy, the premier method for molecular resolution of protein structural fluctuations, to occupy a place at the forefront of investigations into protein allostery. Here, we detail several key methods of NMR chemical shift analysis to extract mechanistic information about long-range chemical signaling in a protein, focusing on practical methodological aspects and the circumstances under which a given approach would be relevant. We also detail some of the experimental considerations that should be made when applying these methods to specific protein systems.

MptpA Kinetics Enhanced by Allosteric Control of an Active Conformation

Understanding allostery in the Mycobacterium tuberculosis low molecular weight protein tyrosine phosphatase (MptpA) is a subject of great interest since MptpA is one of two protein tyrosine phosphatases (PTPs) from the pathogenic organism Mycobacterium tuberculosis expressed during host cell infection. Here, we combine computational modeling with solution NMR spectroscopy and we find that Q75 is an allosteric site. Removal of the polar side chain of Q75 by mutation to leucine results in a cascade of events that reposition the acid loop over the active site and relocates the catalytic aspartic acid (D126) at an optimal position for proton donation to the leaving aryl group of the substrate and for subsequent hydrolysis of the thiophosphoryl intermediate. The computational analysis is consistent with kinetic data, and NMR spectroscopy, showing that the Q75L mutant exhibits enhanced reaction kinetics with similar substrate binding affinity. We anticipate that our findings will motivate further studies on the possibility that MptpA remains passivated during the chronic state of infection and increases its activity as part of the pathogenic life cycle of M. tuberculosis possibly via allosteric means.

Identification of Core Allosteric Sites through Temperature- and Nucleus-Invariant Chemical Shift Covariance

Allosteric regulation is essential to control biological function. In addition, allosteric sites offer a promising venue for selective drug targeting. However, accurate mapping of allosteric sites remains challenging since allostery relies on often subtle, yet functionally relevant, structural and dynamical changes. A viable approach proposed to overcome such challenge is the chemical shift covariance analysis (CHESCA). Although CHESCA offers an exhaustive map of allosteric networks, it is critical to define the core allosteric sites to be prioritized in subsequent functional studies or the design of allosteric drugs. Here we propose two new CHESCA-based methodologies, called temperature CHESCA (T-CHESCA) and CLASS-CHESCA, aimed at narrowing down allosteric maps to the core allosteric residues. Both T- and CLASS-CHESCAs rely on the invariance of core inter-residue correlations to changes in the chemical shifts of the active and inactive conformations interconverting in fast exchange. In the T-CHESCA the chemical shifts of such states are modulated through temperature changes, while in the CLASS-CHESCA through variations in the spin-active nuclei involved in pairwise correlations. The T- and CLASS-CHESCAs as well as complete-linkage CHESCA were applied to the cAMP-binding domain of the exchange protein directly activated by cAMP (EPAC), which serves as a prototypical allosteric switch. Residues consistently identified by the three CHESCA methods were found in previously identified EPAC allosteric core sites. Hence, the T-, CLASS- and CL-CHESCA provide a toolset to establish allosteric site hierarchy and triage allosteric sites to be further analyzed by mutations and functional assays. Furthermore, the core allosteric networks selectively revealed through T- and CLASS-CHESCA are expected to facilitate the mechanistic understanding of disease-related mutations and the design of selective allosteric modulators.

Allosteric pluripotency: challenges and opportunities

Allosteric pluripotency arises when the functional response of an allosteric receptor to an allosteric stimulus depends on additional allosteric modulators. Here, we discuss allosteric pluripotency as observed in the prototypical Protein Kinase A (PKA) as well as in other signaling systems, from typical multidomain signaling proteins to bacterial enzymes. We identify key drivers of pluripotent allostery and illustrate how hypothesizing allosteric pluripotency may solve apparent discrepancies currently present in the literature regarding the dual nature of known allosteric modulators. We also outline the implications of allosteric pluripotency for cellular signaling and allosteric drug design, and analyze the challenges and opportunities opened by the pluripotent nature of allostery.

Allosteric Regulation of Cyclic Nucleotide Dependent Protein Kinases

Kinases include a wide variety of valuable drug targets, but full therapeutic exploitation requires a high degree of selectivity. A promising avenue to engineer the desired kinase selectivity relies on allosteric sites. Here we provide a focused minireview of recent progress in allosteric modulation of cyclic nucleotide-dependent kinases, including protein kinases A and G. We show how apparently diverse emerging concepts such as allosteric pluripotency, allosteric non-additive binding and uncompetitive allosteric inhibition are all manifestations of complex conformational ensembles. Such ensembles include not only the typical apo-inactive and effector-bound-active states, but also mixed intermediates that feature attributes of the former states within a single molecule. We also discuss how allosteric responses are amplified by aggregation processes, thus establishing a novel interface between the signaling and amyloid fields. Finally, we critically evaluate the challenges and opportunities for clinical translation opened by these emerging allosteric concepts.

Is Disrupted Nucleotide-Substrate Cooperativity a Common Trait for Cushing's Syndrome Driving Mutations of Protein Kinase A?

Somatic mutations in the PRKACA gene encoding the catalytic α subunit of protein kinase A (PKA-C) are responsible for cortisol-producing adrenocortical adenomas. These benign neoplasms contribute to the development of Cushing's syndrome. The majority of these mutations occur at the interface between the two lobes of PKA-C and interfere with the enzyme's ability to recognize substrates and regulatory (R) subunits, leading to aberrant phosphorylation patterns and activation. Rarely, patients with similar phenotypes carry an allosteric mutation, E31V, located at the C-terminal end of the αA-helix and adjacent to the αC-helix, but structurally distinct from the PKA-C/R subunit interface mutations. Using a combination of solution NMR, thermodynamics, kinetic assays, and molecular dynamics simulations, we show that the E31V allosteric mutation disrupts central communication nodes between the N- and C- lobes of the enzyme as well as nucleotide-substrate binding cooperativity, a hallmark for kinases' substrate fidelity and regulation. For both orthosteric (L205R and W196R) and allosteric (E31V) Cushing's syndrome mutants, the loss of binding cooperativity is proportional to the density of the intramolecular allosteric network. This structure-activity relationship suggests a possible common mechanism for Cushing's syndrome driving mutations in which decreased nucleotide/substrate binding cooperativity is linked to loss in substrate fidelity and dysfunctional regulation.

Biochemical and structural characterization of analogs of MRE11 breast cancer-associated mutant F237C

The MRE11-RAD50-NBS1 (MRN) protein complex plays a vital role in DNA double strand break sensing, signaling, and repair. Mutation in any component of this complex may lead to disease as disrupting DNA double strand break repair has the potential to cause translocations and loss of genomic information. Here, we have investigated an MRE11 mutation, F237C, identified in a breast cancer tumor. We found that the analogous mutant of Pyrococcus furiosus Mre11 diminishes both the exonuclease and endonuclease activities of Mre11 in vitro. Solution state NMR experiments show that this mutant causes structural changes in the DNA-bound Mre11 for both exo- and endonuclease substrates and causes the protein to become generally more rigid. Moreover, by comparing the NMR data for this cancer-associated mutant with two previously described Mre11 separation-of-nuclease function mutants, a potential allosteric network was detected within Mre11 that connects the active site to regions responsible for recognizing the DNA ends and for dimerization. Together, our data further highlight the dynamics required for Mre11 nuclease function and illuminate the presence of allostery within the enzyme.

Defective internal allosteric network imparts dysfunctional ATP/substrate-binding cooperativity in oncogenic chimera of protein kinase A

An aberrant fusion of the DNAJB1 and PRKACA genes generates a chimeric protein kinase (PKA-CDNAJB1) in which the J-domain of the heat shock protein 40 is fused to the catalytic α subunit of cAMP-dependent protein kinase A (PKA-C). Deceivingly, this chimeric construct appears to be fully functional, as it phosphorylates canonical substrates, forms holoenzymes, responds to cAMP activation, and recognizes the endogenous inhibitor PKI. Nonetheless, PKA-CDNAJB1 has been recognized as the primary driver of fibrolamellar hepatocellular carcinoma and is implicated in other neoplasms for which the molecular mechanisms remain elusive. Here we determined the chimera's allosteric response to nucleotide and pseudo-substrate binding. We found that the fusion of the dynamic J-domain to PKA-C disrupts the internal allosteric network, causing dramatic attenuation of the nucleotide/PKI binding cooperativity. Our findings suggest that the reduced allosteric cooperativity exhibited by PKA-CDNAJB1 alters specific recognitions and interactions between substrates and regulatory partners contributing to dysregulation.

Distinct conformational dynamics and allosteric networks in alpha tryptophan synthase during active catalysis

Experimental observations of enzymes under active turnover conditions have brought new insight into the role of protein motions and allosteric networks in catalysis. Many of these studies characterize enzymes under dynamic chemical equilibrium conditions, in which the enzyme is actively catalyzing both the forward and reverse reactions during data acquisition. We have previously analyzed conformational dynamics and allosteric networks of the alpha subunit of tryptophan synthase under such conditions using NMR. We have proposed that this working state represents a four to one ratio of the enzyme bound with the indole-3-glycerol phosphate substrate (E:IGP) to the enzyme bound with the products indole and glyceraldehyde-3-phosphate (E:indole:G3P). Here, we analyze the inactive D60N variant to deconvolute the contributions of the substrate- and products-bound states to the working state. While the D60N substitution itself induces small structural and dynamic changes, the D60N E:IGP and E:indole:G3P states cannot entirely account for the conformational dynamics and allosteric networks present in the working state. The act of chemical bond breakage and/or formation, or possibly the generation of an intermediate, may alter the structure and dynamics present in the working state. As the enzyme transitions from the substrate-bound to the products-bound state, millisecond conformational exchange processes are quenched and new allosteric connections are made between the alpha active site and the surface which interfaces with the beta subunit. The structural ordering of the enzyme and these new allosteric connections may be important in coordinating the channeling of the indole product into the beta subunit.

The dynamic nature of the Mre11-Rad50 DNA break repair complex

The Mre11-Rad50-Nbs1/Xrs2 protein complex plays a pivotal role in the detection and repair of DNA double strand breaks. Through traditional and emerging structural biology techniques, various functional structural states of this complex have been visualized; however, relatively little is known about the transitions between these states. Indeed, it is these structural transitions that are important for Mre11-Rad50-mediated DNA unwinding at a break and the activation of downstream repair signaling events. Here, we present a brief overview of the current understanding of the structure of the core Mre11-Rad50 complex. We then highlight our recent studies emphasizing the contributions of solution state NMR spectroscopy and other biophysical techniques in providing insight into the structures and dynamics associated with Mre11-Rad50 functions.

Driving Protein Conformational Cycles in Physiology and Disease: "Frustrated" Amino Acid Interaction Networks Define Dynamic Energy Landscapes: Amino Acid Interaction Networks Change Progressively Along Alpha Tryptophan Synthase's Catalytic Cycle

A general framework by which dynamic interactions within a protein will promote the necessary series of structural changes, or "conformational cycle," required for function is proposed. It is suggested that the free-energy landscape of a protein is biased toward this conformational cycle. Fluctuations into higher energy, although thermally accessible, conformations drive the conformational cycle forward. The amino acid interaction network is defined as those intraprotein interactions that contribute most to the free-energy landscape. Some network connections are consistent in every structural state, while others periodically change their interaction strength according to the conformational cycle. It is reviewed here that structural transitions change these periodic network connections, which then predisposes the protein toward the next set of network changes, and hence the next structural change. These concepts are illustrated by recent work on tryptophan synthase. Disruption of these dynamic connections may lead to aberrant protein function and disease states.

Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning

Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.

Mechanism of allosteric inhibition in the Plasmodium falciparum cGMP-dependent protein kinase

Most malaria deaths are caused by the protozoan parasite Plasmodium falciparum Its life cycle is regulated by a cGMP-dependent protein kinase (PfPKG), whose inhibition is a promising antimalaria strategy. Allosteric kinase inhibitors, such as cGMP analogs, offer enhanced selectivity relative to competitive kinase inhibitors. However, the mechanisms underlying allosteric PfPKG inhibition are incompletely understood. Here, we show that 8-NBD-cGMP is an effective PfPKG antagonist. Using comparative NMR analyses of a key regulatory domain, PfD, in its apo, cGMP-bound, and cGMP analog-bound states, we elucidated its inhibition mechanism of action. Using NMR chemical shift analyses, molecular dynamics simulations, and site-directed mutagenesis, we show that 8-NBD-cGMP inhibits PfPKG not simply by reverting a two-state active versus inactive equilibrium, but by sampling also a distinct inactive "mixed" intermediate. Surface plasmon resonance indicates that the ability to stabilize a mixed intermediate provides a means to effectively inhibit PfPKG, without losing affinity for the cGMP analog. Our proposed model may facilitate the rational design of PfPKG-selective inhibitors for improved management of malaria.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Cushing's syndrome driver mutation disrupts protein kinase A allosteric network, altering both regulation and substrate specificity

Genetic alterations in the PRKACA gene coding for the catalytic α subunit of the cAMP-dependent protein kinase A (PKA-C) are linked to cortisol-secreting adrenocortical adenomas, resulting in Cushing's syndrome. Among those, a single mutation (L205R) has been found in up to 67% of patients. Because the x-ray structures of the wild-type and mutant kinases are essentially identical, the mechanism explaining aberrant function of this mutant remains under active debate. Using NMR spectroscopy, thermodynamics, kinetic assays, and molecular dynamics simulations, we found that this single mutation causes global changes in the enzyme, disrupting the intramolecular allosteric network and eliciting losses in nucleotide/pseudo-substrate binding cooperativity. Remarkably, by rewiring its internal allosteric network, PKA-CL205R is able to bind and phosphorylate non-canonical substrates, explaining its changes in substrate specificity. Both the lack of regulation and change in substrate specificity reveal the complex role of this mutated kinase in the formation of cortisol-secreting adrenocortical adenomas.

Inter-Active Site Communication Mediated by the Dimer Interface β-Sheet in the Half-the-Sites Enzyme, Thymidylate Synthase

Thymidylate synthase (TS) is a dimeric enzyme conserved in all life forms that exhibits the allosteric feature of half-the-sites activity. Neither the reason for nor the mechanism of this phenomenon is understood. We used a combined nuclear magnetic resonance (NMR) and molecular dynamics approach to study a stable intermediate preceding hydride transfer, which is the rate-limiting and half-the-sites step. In NMR titrations with ligands leading to this intermediate, we measured chemical shifts of the apoenzyme (lig₀), the saturated holoenzyme (lig₂), and the typically elusive singly bound (lig₁) states. Approximately 40 amides showed quartet patterns providing direct NMR evidence of coupling between the active site and probes >30 Å away in the distal subunit. Quartet peak patterns have symmetrical character, indicating reciprocity in communicating the first and second binding events to the distal protomer. Quartets include key catalytic residues and map to the dimer interface β-sheet, which also represents the shortest path between the two active sites. Simulations corroborate the coupling observed in solution in that there is excellent overlap between quartet residues and main-chain atoms having intersubunit cross-correlated motions. Simulations identify five hot spot residues, three of which lie at the kink in the unique β-bulge abutting the active sites on either end of the sheet. Interstrand cross-correlated motions become more organized and pronounced as the enzyme progresses from lig₀ to lig₁ and ultimately lig₂. Coupling in the apparently symmetrical complex has implications for half-the-sites reactivity and potentially resolves the paradox of inequivalent TS active sites despite the vast majority of X-ray structures appearing to be symmetrical.

Mechanisms of Specific versus Nonspecific Interactions of Aggregation-Prone Inhibitors and Attenuators

A common source of false positives in drug discovery is ligand self-association into large colloidal assemblies that nonspecifically inhibit target proteins. However, the mechanisms of aggregation-based inhibition (ABI) and ABI-attenuation by additives, such as Triton X-100 (TX) and human serum albumin (HSA), are not fully understood. Here, we investigate the molecular basis of ABI and ABI-attenuation through the lens of NMR and coupled thermodynamic cycles. We unexpectedly discover a new class of aggregating ligands that exhibit negligible interactions with proteins but act as competitive sinks for the free inhibitor, resulting in bell-shaped dose-response curves. TX attenuates ABI by converting inhibitory, protein-binding aggregates into nonbinding coaggregates, whereas HSA minimizes nonspecific ligand interactions by functioning as a reservoir for free inhibitor and preventing self-association. Hence, both TX and HSA are useful tools to minimize false positives arising from nonspecific binding but at the cost of potentially introducing false negatives due to suppression of specific interactions.

Globally correlated conformational entropy underlies positive and negative cooperativity in a kinase's enzymatic cycle

Enzymes accelerate the rate of chemical transformations by reducing the activation barriers of uncatalyzed reactions. For signaling enzymes, substrate recognition, binding, and product release are often rate-determining steps in which enthalpy-entropy compensation plays a crucial role. While the nature of enthalpic interactions can be inferred from structural data, the molecular origin and role of entropy in enzyme catalysis remains poorly understood. Using thermocalorimetry, NMR, and MD simulations, we studied the conformational landscape of the catalytic subunit of cAMP-dependent protein kinase A, a ubiquitous phosphoryl transferase involved in a myriad of cellular processes. Along the enzymatic cycle, the kinase exhibits positive and negative cooperativity for substrate and nucleotide binding and product release. We found that globally coordinated changes of conformational entropy activated by ligand binding, together with synchronous and asynchronous breathing motions of the enzyme, underlie allosteric cooperativity along the kinase's cycle.

Identifying coupled clusters of allostery participants through chemical shift perturbations

Allosteric couplings underlie many cellular signaling processes and provide an exciting avenue for development of new diagnostics and therapeutics. A general method for identifying important residues in allosteric mechanisms would be very useful, but remains elusive due to the complexity of long-range phenomena. Here, we introduce an NMR method to identify residues involved in allosteric coupling between two ligand-binding sites in a protein, which we call chemical shift detection of allostery participants (CAP). Networks of functional groups responding to each ligand are defined through correlated NMR perturbations. In this process, we also identify allostery participants, groups that respond to both binding events and likely play a role in the coupling between the binding sites. Such residues exhibit multiple functional states with distinct NMR chemical shifts, depending on binding status at both binding sites. Such a strategy was applied to the prototypical ion channel KcsA. We had previously shown that the potassium affinity at the extracellular selectivity filter is strongly dependent on proton binding at the intracellular pH sensor. Here, we analyzed proton and potassium binding networks and identified groups that depend on both proton and potassium binding (allostery participants). These groups are viewed as candidates for transmitting information between functional units. The vital role of one such identified amino acid was validated through site-specific mutagenesis, electrophysiology functional studies, and NMR-detected thermodynamic analysis of allosteric coupling. This strategy for identifying allostery participants is likely to have applications for many other systems.

Engineering Allostery into Proteins

Our ability to engineer protein structure and function has grown dramatically over recent years. Perhaps the next level in protein design is to develop proteins whose function can be regulated in response to various stimuli, including ligand binding, pH changes, and light. Endeavors toward these goals have tested and expanded on our understanding of protein function and allosteric regulation. In this chapter, we provide examples from different methods for developing new allosterically regulated proteins. These methods range from whole insertion of regulatory domains into new host proteins, to covalent attachment of photoswitches to generate light-responsive proteins, and to targeted changes to specific amino acid residues, especially to residues identified to be important for relaying allosteric information across the protein framework. Many of the examples we discuss have already found practical use in medical and biotechnology applications.

Mechanism of Selective Enzyme Inhibition through Uncompetitive Regulation of an Allosteric Agonist

Classical uncompetitive inhibitors are potent pharmacological modulators of enzyme function. Since they selectively target enzyme-substrate complexes (E:S), their inhibitory potency is amplified by increasing substrate concentrations. Recently, an unconventional uncompetitive inhibitor, called CE3F4R, was discovered for the exchange protein activated by cAMP isoform 1 (EPAC1). Unlike conventional uncompetitive inhibitors, CE3F4R is uncompetitive with respect to an allosteric effector, cAMP, as opposed to the substrate (i.e., CE3F4R targets the E:cAMP rather than the E:S complex). However, the mechanism of CE3F4R as an uncompetitive inhibitor is currently unknown. Here, we elucidate the mechanism of CE3F4R's action using NMR spectroscopy. Due to limited solubility and line broadening, which pose major challenges for traditional structural determination approaches, we resorted to a combination of protein- and ligand-based NMR experiments to comparatively analyze EPAC mutations, inhibitor analogs, and cyclic nucleotide derivatives that trap EPAC at different stages of activation. We discovered that CE3F4R binds within the EPAC cAMP-binding domain (CBD) at a subdomain interface distinct from the cAMP binding site, acting as a wedge that stabilizes a cAMP-bound mixed-intermediate. The mixed-intermediate includes attributes of both the apo/inactive and cAMP-bound/active states. In particular, the intermediate targeted by CE3F4R traps a CBD's hinge helix in its inactive conformation, locking EPAC into a closed domain topology that restricts substrate access to the catalytic domain. The proposed mechanism of action also explains the isoform selectivity of CE3F4R in terms of a single EPAC1 versus EPAC2 amino acid difference that destabilizes the active conformation of the hinge helix.

Dynamic Allostery in PLCγ1 and Its Modulation by a Cancer Mutation Revealed by MD Simulation and NMR

Phosphatidylinositol phospholipase Cγ (PLCγ) is an intracellular membrane-associated second-messenger signaling protein activated by tyrosine kinases such as fibroblast growth factor receptor 1. PLCγ contains the regulatory γ-specific array (γSA) comprising a tandem Src homology 2 (SH2) pair, an SH3 domain, and a split pleckstrin homology domain. Binding of an activated growth factor receptor to γSA leads to Tyr783 phosphorylation and consequent PLCγ activation. Several disease-relevant mutations in γSA have been identified; all lead to elevated phospholipase activity. In this work, we describe an allosteric mechanism that connects the Tyr783 phosphorylation site to the nSH2-cSH2 junction and involves dynamic interactions between the cSH2-SH3 linker and cSH2. Molecular dynamics simulations of the tandem SH2 protein suggest that Tyr783 phosphorylation is communicated to the nSH2-cSH2 junction by modulating cSH2 binding to sections of the cSH2-SH3 linker. NMR chemical shift perturbation analyses for designed tandem SH2 constructs reveal combined fast and slow dynamic processes that can be attributed to allosteric communication involving these regions of the protein, establishing an example in which complex N-site exchange can be directly inferred from 1H,15N-HSQC spectra. Furthermore, in tandem SH2 and γSA constructs, molecular dynamics and NMR results show that the Arg687Trp mutant in PLCγ1 (equivalent to the cancer mutation Arg665Trp in PLCγ2) perturbs the dynamic allosteric pathway. This combined experimental and computational study reveals a rare example of multistate kinetics involved in a dynamic allosteric process that is modulated in the context of a disease-relevant mutation. The allosteric influences and the weakened binding of the cSH2-SH3 linker to cSH2 should be taken into account in any more holistic investigation of PLCγ regulation.

Role of Dimers in the cAMP-Dependent Activation of Hyperpolarization-Activated Cyclic-Nucleotide-Modulated (HCN) Ion Channels

Hyperpolarization-activated cyclic-nucleotide-modulated (HCN) ion channels control rhythmicity in neurons and cardiomyocytes. Cyclic AMP (cAMP) modulates HCN activity through the cAMP-induced formation of a tetrameric gating ring spanning the intracellular region (IR) of HCN. Although evidence from confocal patch-clamp fluorometry indicates that the cAMP-dependent gating of HCN occurs through a dimer of dimers, the structural and dynamical basis of cAMP allostery in HCN dimers has so far remained elusive. Thus, here we examine how dimers influence IR structural dynamics, and the role that such structural dynamics play in HCN allostery. To this end, we performed molecular dynamics (MD) simulations of HCN4 IR dimers in their fully apo, fully holo, and partially cAMP-bound states, resulting in a total simulated time of 1.2 μs. Comparative analyses of these MD trajectories, as well as previous monomer and tetramer simulations utilized as benchmarks for comparison, reveal that dimers markedly sensitize the HCN IR to cAMP-modulated allostery. Our results indicate that dimerization fine-tunes the IR dynamics to enhance, relative to both monomers and tetramers, the allosteric intra- and interprotomer coupling between the cAMP-binding domain and tetramerization domain components of the IR. The resulting allosteric model provides a viable rationalization of electrophysiological data on the role of IR dimers in HCN activation.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

A dynamic allosteric pathway underlies Rad50 ABC ATPase function in DNA repair

The Mre11-Rad50 protein complex is an initial responder to sites of DNA double strand breaks. Many studies have shown that ATP binding to Rad50 causes global changes to the Mre11-Rad50 structure, which are important for DNA repair functions. Here we used methyl-based NMR spectroscopy on a series of mutants to describe a dynamic allosteric pathway within Rad50. Mutations result in changes in the side chain methyl group chemical environment that are correlated with altered nanosecond timescale dynamics. We also observe striking relationships between the magnitude of chemical shift perturbations and Rad50 and Mre11 activities. Together, these data suggest an equilibrium between a ground state and an "active" dimerization competent state of Rad50 that has locally altered structure and dynamics and is poised for ATP-induced dimerization and eventual ATP hydrolysis. Thus, this sparsely populated intermediate is critical for Mre11-Rad50-directed DNA double strand break repair.

A solution NMR toolset to probe the molecular mechanisms of amyloid inhibitors

A chemical exchange-based solution NMR toolset to probe the molecular mechanisms of amyloid inhibitors.

Implementation of the NMR CHEmical Shift Covariance Analysis (CHESCA): A Chemical Biologist's Approach to Allostery

Mapping allosteric sites is emerging as one of the central challenges in physiology, pathology, and pharmacology. Nuclear Magnetic Resonance (NMR) spectroscopy is ideally suited to map allosteric sites, given its ability to sense at atomic resolution the dynamics underlying allostery. Here, we focus specifically on the NMR CHEmical Shift Covariance Analysis (CHESCA), in which allosteric systems are interrogated through a targeted library of perturbations (e.g., mutations and/or analogs of the allosteric effector ligand). The atomic resolution readout for the response to such perturbation library is provided by NMR chemical shifts. These are then subject to statistical correlation and covariance analyses resulting in clusters of allosterically coupled residues that exhibit concerted responses to the common set of perturbations. This chapter provides a description of how each step in the CHESCA is implemented, starting from the selection of the perturbation library and ending with an overview of different clustering options.

The PROSECCO server for chemical shift predictions in ordered and disordered proteins

The chemical shifts measured in solution-state and solid-state nuclear magnetic resonance (NMR) are powerful probes of the structure and dynamics of protein molecules. The exploitation of chemical shifts requires methods to correlate these data with the protein structures and sequences. We present here an approach to calculate accurate chemical shifts in both ordered and disordered proteins using exclusively the information contained in their sequences. Our sequence-based approach, protein sequences and chemical shift correlations (PROSECCO), achieves the accuracy of the most advanced structure-based methods in the characterization of chemical shifts of folded proteins and improves the state of the art in the study of disordered proteins. Our analyses revealed fundamental insights on the structural information carried by NMR chemical shifts of structured and unstructured protein states.

A Universal Pattern in the Percolation and Dissipation of Protein Structural Perturbations

Understanding the extent to which information is transmitted through the intramolecular interaction network of proteins upon a perturbation, that is, an allosteric effect, has long remained an unsolved problem. Through an analysis of high-resolution NMR data from the literature on 28 different proteins and 49 structural perturbations, we show that the extent of induced structural changes through mutations and molecular events including protein-protein, protein-peptide, protein-ligand binding, and post-translational modifications exhibit a near-universal exponential functional form. The extent of percolation into the protein structures can be up to 20-25 Å despite no apparent change in the 3D structures. These observations are also consistent with theoretical expectations, elementary graph theoretic analysis of protein structures, detailed molecular dynamics simulations, and experimental double-mutant cycles. Our analysis highlights that most molecular events would contribute to allosteric effects independent of protein structure, topology, or identity and provides a simple avenue to test and potentially model their effects.

Functional dynamics in cyclic nucleotide signaling and amyloid inhibition

It is now established that understanding the molecular basis of biological function requires atomic resolution maps of both structure and dynamics. Here, we review several illustrative examples of functional dynamics selected from our work on cyclic nucleotide signaling and amyloid inhibition. Although fundamentally diverse, a central aspect common to both fields is that function can only be rationalized by considering dynamic equilibria between distinct states of the accessible free energy landscape. The dynamic exchange between ground and excited states of signaling proteins is essential to explain auto-inhibition and allosteric activation. The dynamic exchange between non-toxic monomeric species and toxic oligomers of amyloidogenic proteins provides a foundation to understand amyloid inhibition. NMR ideally probes both types of dynamic exchange at atomic resolution. Specifically, we will show how NMR was utilized to reveal the dynamical basis of cyclic nucleotide affinity, selectivity, agonism and antagonism in multiple eukaryotic cAMP and cGMP receptors. We will also illustrate how NMR revealed the mechanism of action of plasma proteins that act as extracellular chaperones and inhibit the self-association of the prototypical amyloidogenic Aβ peptide. The examples outlined in this review illustrate the widespread implications of functional dynamics and the power of NMR as an indispensable tool in molecular pharmacology and pathology.

Calcium-Mediated Control of S100 Proteins: Allosteric Communication via an Agitator/Signal Blocking Mechanism

Allosteric proteins possess dynamically coupled residues for the propagation of input signals to distant target binding sites. The input signals usually correspond to "effector is present" or "effector is not present". Many aspects of allosteric regulation remain incompletely understood. This work focused on S100A11, a dimeric EF-hand protein with two hydrophobic target binding sites. An annexin peptide (Ax) served as the target. Target binding is allosterically controlled by Ca²⁺ over a distance of ∼26 Å. Ca²⁺ promotes formation of a [Ca₄ S100 Ax₂] complex, where the Ax peptides are accommodated between helices III/IV and III'/IV'. Without Ca²⁺ these binding sites are closed, precluding interactions with Ax. The allosteric mechanism was probed by microsecond MD simulations in explicit water, complemented by hydrogen exchange mass spectrometry (HDX/MS). Consistent with experimental data, MD runs in the absence of Ca²⁺ and Ax culminated in target binding site closure. In simulations on [Ca₄ S100] the target binding sites remained open. These results capture the essence of allosteric control, revealing how Ca²⁺ prevents binding site closure. Both HDX/MS and MD data showed that the metalation sites become more dynamic after Ca²⁺ loss. However, these enhanced dynamics do not represent the primary trigger of the allosteric cascade. Instead, a labile salt bridge acts as an incessantly active "agitator" that destabilizes the packing of adjacent residues, causing a domino chain of events that culminates in target binding site closure. This agitator represents the starting point of the allosteric signal propagation pathway. Ca²⁺ binding rigidifies elements along this pathway, thereby blocking signal transmission. This blocking mechanism does not conform to the commonly held view that allosteric communication pathways generally originate at the sites where effectors interact with the protein.

Applications of NMR and computational methodologies to study protein dynamics

Overwhelming evidence now illustrates the defining role of atomic-scale protein flexibility in biological events such as allostery, cell signaling, and enzyme catalysis. Over the years, spin relaxation nuclear magnetic resonance (NMR) has provided significant insights on the structural motions occurring on multiple time frames over the course of a protein life span. The present review article aims to illustrate to the broader community how this technique continues to shape many areas of protein science and engineering, in addition to being an indispensable tool for studying atomic-scale motions and functional characterization. Continuing developments in underlying NMR technology alongside software and hardware developments for complementary computational approaches now enable methodologies to routinely provide spatial directionality and structural representations traditionally harder to achieve solely using NMR spectroscopy. In addition to its well-established role in structural elucidation, we present recent examples that illustrate the combined power of selective isotope labeling, relaxation dispersion experiments, chemical shift analyses, and computational approaches for the characterization of conformational sub-states in proteins and enzymes.

Free energy landscape remodeling of the cardiac pacemaker channel explains the molecular basis of familial sinus bradycardia

The hyperpolarization-activated and cyclic nucleotide-modulated ion channel (HCN) drives the pacemaker activity in the heart, and its malfunction can result in heart disorders. One such disorder, familial sinus bradycardia, is caused by the S672R mutation in HCN, whose electrophysiological phenotypes include a negative shift in the channel activation voltage and an accelerated HCN deactivation. The outcomes of these changes are abnormally low resting heart rates. However, the molecular mechanism underlying these electrophysiological changes is currently not fully understood. Crystallographic investigations indicate that the S672R mutation causes limited changes in the structure of the HCN intracellular gating tetramer, but its effects on protein dynamics are unknown. Here, we utilize comparative S672R versus WT NMR analyses to show that the S672R mutation results in extensive perturbations of the dynamics in both apo- and holo-forms of the HCN4 isoform, reflecting how S672R remodels the free energy landscape for the modulation of HCN4 by cAMP, i.e. the primary cyclic nucleotide modulator of HCN channels. We show that the S672R mutation results in a constitutive shift of the dynamic auto-inhibitory equilibrium toward inactive states of HCN4 and broadens the free-energy well of the apo-form, enhancing the millisecond to microsecond dynamics of the holo-form at sites critical for gating cAMP binding. These S672R-induced variations in dynamics provide a molecular basis for the electrophysiological phenotypes of this mutation and demonstrate that the pathogenic effects of the S672R mutation can be rationalized primarily in terms of modulations of protein dynamics.

NMR reveals a dynamic allosteric pathway in thrombin

Although serine proteases are found ubiquitously in both eukaryotes and prokaryotes, and they comprise the largest of all of the peptidase families, their dynamic motions remain obscure. The backbone dynamics of the coagulation serine protease, apo-thrombin (S195M-thrombin), were compared to the substrate-bound form (PPACK-thrombin). R₁, R₂, ¹⁵N-{¹H}NOEs, and relaxation dispersion NMR experiments were measured to capture motions across the ps to ms timescale. The ps-ns motions were not significantly altered upon substrate binding. The relaxation dispersion data revealed that apo-thrombin is highly dynamic, with μs-ms motions throughout the molecule. The region around the N-terminus of the heavy chain, the Na⁺-binding loop, and the 170 s loop, all of which are implicated in allosteric coupling between effector binding sites and the active site, were dynamic primarily in the apo-form. Most of the loops surrounding the active site become more ordered upon PPACK-binding, but residues in the N-terminal part of the heavy chain, the γ-loop, and anion-binding exosite 1, the main allosteric binding site, retain μs-ms motions. These residues form a dynamic allosteric pathway connecting the active site to the main allosteric site that remains in the substrate-bound form.

Elucidation of a four-site allosteric network in fibroblast growth factor receptor tyrosine kinases

Receptor tyrosine kinase (RTK) signaling is tightly regulated by protein allostery within the intracellular tyrosine kinase domains. Yet the molecular determinants of allosteric connectivity in tyrosine kinase domain are incompletely understood. By means of structural (X-ray and NMR) and functional characterization of pathogenic gain-of-function mutations affecting the FGF receptor (FGFR) tyrosine kinase domain, we elucidated a long-distance allosteric network composed of four interconnected sites termed the 'molecular brake', 'DFG latch', 'A-loop plug', and 'αC tether'. The first three sites repress the kinase from adopting an active conformation, whereas the αC tether promotes the active conformation. The skewed design of this four-site allosteric network imposes tight autoinhibition and accounts for the incomplete mimicry of the activated conformation by pathogenic mutations targeting a single site. Based on the structural similarity shared among RTKs, we propose that this allosteric model for FGFR kinases is applicable to other RTKs.

NMR Characterization of Information Flow and Allosteric Communities in the MAP Kinase p38γ

The intramolecular network structure of a protein provides valuable insights into allosteric sites and communication pathways. However, a straightforward method to comprehensively map and characterize these pathways is not currently available. Here we present an approach to characterize intramolecular network structure using NMR chemical shift perturbations. We apply the method to the mitogen activated protein kinase (MAPK) p38γ. p38γ contains allosteric sites that are conserved among eukaryotic kinases as well as unique to the MAPK family. How these regulatory sites communicate with catalytic residues is not well understood. Using our method, we observe and characterize for the first time information flux between regulatory sites through a conserved kinase infrastructure. This network is accessed, reinforced, and broken in various states of p38γ, reflecting the functional state of the protein. We demonstrate that the approach detects critical junctions in the network corresponding to biologically significant allosteric sites and pathways.

Biophysical and computational methods to analyze amino acid interaction networks in proteins

Globular proteins are held together by interacting networks of amino acid residues. A number of different structural and computational methods have been developed to interrogate these amino acid networks. In this review, we describe some of these methods, including analyses of X-ray crystallographic data and structures, computer simulations, NMR data, and covariation among protein sequences, and indicate the critical insights that such methods provide into protein function. This information can be leveraged towards the design of new allosteric drugs, and the engineering of new protein function and protein regulation strategies.

Solution NMR Spectroscopy for the Study of Enzyme Allostery

Allostery is a ubiquitous biological regulatory process in which distant binding sites within a protein or enzyme are functionally and thermodynamically coupled. Allosteric interactions play essential roles in many enzymological mechanisms, often facilitating formation of enzyme-substrate complexes and/or product release. Thus, elucidating the forces that drive allostery is critical to understanding the complex transformations of biomolecules. Currently, a number of models exist to describe allosteric behavior, taking into account energetics as well as conformational rearrangements and fluctuations. In the following Review, we discuss the use of solution NMR techniques designed to probe allosteric mechanisms in enzymes. NMR spectroscopy is unequaled in its ability to detect structural and dynamical changes in biomolecules, and the case studies presented herein demonstrate the range of insights to be gained from this valuable method. We also provide a detailed technical discussion of several specialized NMR experiments that are ideally suited for the study of enzymatic allostery.

Functionalized N,N-Diphenylamines as Potent and Selective EPAC2 Inhibitors

N,N-Diphenylamines were discovered as potent and selective EPAC2 inhibitors. A study was conducted to determine the structure-activity relationships in a series of inhibitors of which several compounds displayed submicromolar potencies. Selectivity over the related EPAC1 protein was also demonstrated. Computational modeling reveals an allosteric site that is distinct from the cAMP binding domain shared by both EPAC isoforms, providing a theory with regards to subtype selectivity.