The energetic and allosteric landscape for KRAS inhibition

Thousands of proteins have been validated genetically as therapeutic targets for human diseases1. However, very few have been successfully targeted, and many are considered 'undruggable'. This is particularly true for proteins that function via protein-protein interactions-direct inhibition of binding interfaces is difficult and requires the identification of allosteric sites. However, most proteins have no known allosteric sites, and a comprehensive allosteric map does not exist for any protein. Here we address this shortcoming by charting multiple global atlases of inhibitory allosteric communication in KRAS. We quantified the effects of more than 26,000 mutations on the folding of KRAS and its binding to six interaction partners. Genetic interactions in double mutants enabled us to perform biophysical measurements at scale, inferring more than 22,000 causal free energy changes. These energy landscapes quantify how mutations tune the binding specificity of a signalling protein and map the inhibitory allosteric sites for an important therapeutic target. Allosteric propagation is particularly effective across the central β-sheet of KRAS, and multiple surface pockets are genetically validated as allosterically active, including a distal pocket in the C-terminal lobe of the protein. Allosteric mutations typically inhibit binding to all tested effectors, but they can also change the binding specificity, revealing the regulatory, evolutionary and therapeutic potential to tune pathway activation. Using the approach described here, it should be possible to rapidly and comprehensively identify allosteric target sites in many proteins.

Genetic switching by the Lac repressor is based on two-state Monod-Wyman-Changeux allostery

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod-Wyman-Changeux model.

Disease-causing cystathionine β-synthase linker mutations impair allosteric regulation

Cystathionine β-synthase (CBS) catalyzes the committing step in the transsulfuration pathway, which is important for clearing homocysteine and furnishing cysteine. The transsulfuration pathway also generates H2S, a signaling molecule. CBS is a modular protein with a heme and pyridoxal phosphate-binding catalytic core, which is separated by a linker region from the C-terminal regulatory domain that binds S-adenosylmethionine (AdoMet), an allosteric activator. Recent cryo-EM structures reveal that CBS exists in a fibrillar form and undergoes a dramatic architectural rearrangement between the basal and AdoMet-bound states. CBS is the single most common locus of mutations associated with homocystinuria, and, in this study, we have characterized three clinical variants (K384E/N and M391I), which reside in the linker region. The native fibrillar form is destabilized in the variants, and differences in their limited proteolytic fingerprints also reveal conformational alterations. The crystal structure of the truncated K384N variant, lacking the regulatory domain, reveals that the overall fold of the catalytic core is unperturbed. M391I CBS exhibits a modest (1.4-fold) decrease while the K384E/N variants exhibit a significant (∼8-fold) decrease in basal activity, which is either unresponsive to or inhibited by AdoMet. Pre-steady state kinetic analyses reveal that the K384E/N substitutions exhibit pleiotropic effects and that the differences between them are expressed in the second half reaction, that is, homocysteine binding and reaction with the aminoacrylate intermediate. Together, these studies point to an important role for the linker in stabilizing the higher-order oligomeric structure of CBS and enabling AdoMet-dependent regulation.

A critical inter-subunit interaction for the transmission of the allosteric signal in the Agrobacterium tumefaciens ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase is a key regulatory enzyme involved in starch and glycogen synthesis in plants and bacteria, respectively. It has been hypothesized that inter-subunit communications are important for the allosteric effect in this enzyme. However, no specific interactions have been identified as part of the regulatory signal. The enzyme from Agrobacterium tumefaciens is a homotetramer allosterically regulated by fructose 6-phosphate and pyruvate. Three pairs of distinct subunit-subunit interfaces are present. Here we focus on an interface that features two symmetrical interactions between Arg11 and Asp141 from one subunit with residues Asp141 and Arg11 of the neighbor subunit, respectively. Previously, scanning mutagenesis showed that a mutation at the Arg11 position disrupted the activation of the enzyme. Considering the distance of these residues from the allosteric and catalytic sites, we hypothesized that the interaction between Arg11 and Asp141 is critical for allosteric signaling rather than effector binding. To prove our hypothesis, we mutated those two sites (D141A, D141E, D141N, D141R, R11D, and R11K) and performed kinetic and binding analysis. Mutations that altered the charge affected the regulation the most. To prove that the interaction per se (rather than the presence of specific residues) is critical, we partially rescued the R11D protein by introducing a second mutation (R11D/D141R). This could not restore the activator effect on kcat , but it did rescue the effect on substrate affinity. Our results indicate the critical functional role of Arg11 and Asp141 to relay the allosteric signal in this subunit interface.

A complete allosteric map of a GTPase switch in its native cellular network

Allosteric regulation is central to protein function in cellular networks. A fundamental open question is whether cellular regulation of allosteric proteins occurs only at a few defined positions or at many sites distributed throughout the structure. Here, we probe the regulation of GTPases-protein switches that control signaling through regulated conformational cycling-at residue-level resolution by deep mutagenesis in the native biological network. For the GTPase Gsp1/Ran, we find that 28% of the 4,315 assayed mutations show pronounced gain-of-function responses. Twenty of the sixty positions enriched for gain-of-function mutations are outside the canonical GTPase active site switch regions. Kinetic analysis shows that these distal sites are allosterically coupled to the active site. We conclude that the GTPase switch mechanism is broadly sensitive to cellular allosteric regulation. Our systematic discovery of new regulatory sites provides a functional map to interrogate and target GTPases controlling many essential biological processes.

Coevolution-based prediction of key allosteric residues for protein function regulation

Allostery is fundamental to many biological processes. Due to the distant regulation nature, how allosteric mutations, modifications, and effector binding impact protein function is difficult to forecast. In protein engineering, remote mutations cannot be rationally designed without large-scale experimental screening. Allosteric drugs have raised much attention due to their high specificity and possibility of overcoming existing drug-resistant mutations. However, optimization of allosteric compounds remains challenging. Here, we developed a novel computational method KeyAlloSite to predict allosteric site and to identify key allosteric residues (allo-residues) based on the evolutionary coupling model. We found that protein allosteric sites are strongly coupled to orthosteric site compared to non-functional sites. We further inferred key allo-residues by pairwise comparing the difference of evolutionary coupling scores of each residue in the allosteric pocket with the functional site. Our predicted key allo-residues are in accordance with previous experimental studies for typical allosteric proteins like BCR-ABL1, Tar, and PDZ3, as well as key cancer mutations. We also showed that KeyAlloSite can be used to predict key allosteric residues distant from the catalytic site that are important for enzyme catalysis. Our study demonstrates that weak coevolutionary couplings contain important information of protein allosteric regulation function. KeyAlloSite can be applied in studying the evolution of protein allosteric regulation, designing and optimizing allosteric drugs, and performing functional protein design and enzyme engineering.

Allostery: Allosteric Cancer Drivers and Innovative Allosteric Drugs

Here, we discuss the principles of allosteric activating mutations, propagation downstream of the signals that they prompt, and allosteric drugs, with examples from the Ras signaling network. We focus on Abl kinase where mutations shift the landscape toward the active, imatinib binding-incompetent conformation, likely resulting in the high affinity ATP outcompeting drug binding. Recent pharmacological innovation extends to allosteric inhibitor (GNF-5)-linked PROTAC, targeting Bcr-Abl1 myristoylation site, and broadly, allosteric heterobifunctional degraders that destroy targets, rather than inhibiting them. Designed chemical linkers in bifunctional degraders can connect the allosteric ligand that binds the target protein and the E3 ubiquitin ligase warhead anchor. The physical properties and favored conformational state of the engineered linker can precisely coordinate the distance and orientation between the target and the recruited E3. Allosteric PROTACs, noncompetitive molecular glues, and bitopic ligands, with covalent links of allosteric ligands and orthosteric warheads, increase the effective local concentration of productively oriented and placed ligands. Through covalent chemical or peptide linkers, allosteric drugs can collaborate with competitive drugs, degrader anchors, or other molecules of choice, driving innovative drug discovery.

Known allosteric proteins have central roles in genetic disease

Allostery is a form of protein regulation, where ligands that bind sites located apart from the active site can modify the activity of the protein. The molecular mechanisms of allostery have been extensively studied, because allosteric sites are less conserved than active sites, and drugs targeting them are more specific than drugs binding the active sites. Here we quantify the importance of allostery in genetic disease. We show that 1) known allosteric proteins are central in disease networks, contribute to genetic disease and comorbidities much more than non-allosteric proteins, and there is an association between being allosteric and involvement in disease; 2) they are enriched in many major disease types like hematopoietic diseases, cardiovascular diseases, cancers, diabetes, or diseases of the central nervous system; 3) variants from cancer genome-wide association studies are enriched near allosteric proteins, indicating their importance to polygenic traits; and 4) the importance of allosteric proteins in disease is due, at least partly, to their central positions in protein-protein interaction networks, and less due to their dynamical properties.

Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9

CRISPR-Cas9 (clustered regularly interspaced short palindromic repeat and associated Cas9 protein) is a molecular tool with transformative genome editing capabilities. At the molecular level, an intricate allosteric signaling is critical for DNA cleavage, but its role in the specificity enhancement of the Cas9 endonuclease is poorly understood. Here, multi-microsecond molecular dynamics is combined with solution NMR and graph theory-derived models to probe the allosteric role of key specificity-enhancing mutations. We show that mutations responsible for increasing the specificity of Cas9 alter the allosteric structure of the catalytic HNH domain, impacting the signal transmission from the DNA recognition region to the catalytic sites for cleavage. Specifically, the K855A mutation strongly disrupts the allosteric connectivity of the HNH domain, exerting the highest perturbation on the signaling transfer, while K810A and K848A result in more moderate effects on the allosteric communication. This differential perturbation of the allosteric signal correlates to the order of specificity enhancement (K855A > K848A ~ K810A) observed in biochemical studies, with the mutation achieving the highest specificity most strongly perturbing the signaling transfer. These findings suggest that alterations of the allosteric communication from DNA recognition to cleavage are critical to increasing the specificity of Cas9 and that allosteric hotspots can be targeted through mutational studies for improving the system's function.

The novel GCK variant p.Val455Leu associated with hyperinsulinism is susceptible to allosteric activation and is conducive to weight gain and the development of diabetes

AIMS/HYPOTHESIS

The mammalian enzyme glucokinase (GK), expressed predominantly in liver and pancreas, plays an essential role in carbohydrate metabolism. Monogenic GK disorders emphasise the role of GK in determining the blood glucose set point.

METHODS

A family with congenital hyperinsulinism (CHI) was examined for GCK gene variants by Sanger sequencing. A combined approach, involving kinetic analysis (also using GK activators and inhibitors), intracellular translocation assays, insulin secretion measurements and structural modelling, was used to investigate the novel variant compared with known variants.

RESULTS

We report on the novel gain-of-function GCK variant p.Val455Leu (V455L), inherited as an autosomal dominant trait in a German family with CHI and concomitant obesity (fasting blood glucose 2.1 mmol/l, BMI 45.0 kg/m², HOMA-IR 1.5 in an adult female family member); one male family member developed type 2 diabetes until age 35 years (with fasting glucose 2.8-3.7 mmol/l, BMI 38.9 kg/m², HOMA-IR 4.6). Kinetic characterisation of the V455L variant revealed a significant increase in glucose affinity (glucose concentration at which reaction rate is half its maximum rate [S_0.5]: mutant 2.4 ± 0.3 mmol/l vs wild-type 7.6 ± 1.0 mmol/l), accompanied by a distinct additive susceptibility to both the endogenous activator fructose 2,6-bisphosphatase and the synthetic allosteric activator RO-28-1675. The effect of RO-28-1675 was more pronounced when compared with the previously known GK variants V455M and V455E. Binding to the inhibitor glucokinase regulatory protein was unimpaired for V455L and V455E but was reduced for V455M, whereas mannoheptulose inhibited all GK variants and the wild-type enzyme. Structural analyses suggested a role for residue 455 in rearrangements between the inactive and active conformations of GK and also in allosteric activation. Comparison with V455M and V455E and an overview of activating GK variants provided a context for the novel sequence aberration in terms of altered GK enzyme characteristics caused by single amino acid changes.

CONCLUSION/INTERPRETATION

We provide new knowledge on the structure-function relationship of GK, with special emphasis on enzyme activation, potentially yielding fresh strategic insights into breaking the vicious circle of fluctuating blood glucose levels and the attendant risk of long-lasting metabolic changes in both CHI and type 2 diabetes.

Insights into Conformational Dynamics and Allostery in DNMT1-H3Ub/USP7 Interactions

DNA methyltransferases (DNMTs) including DNMT1 are a conserved family of cytosine methylases that play crucial roles in epigenetic regulation. The versatile functions of DNMT1 rely on allosteric networks between its different interacting partners, emerging as novel therapeutic targets. In this work, based on the modeling structures of DNMT1-ubiquitylated H3 (H3Ub)/ubiquitin specific peptidase 7 (USP7) complexes, we have used a combination of elastic network models, molecular dynamics simulations, structural residue perturbation, network modeling, and pocket pathway analysis to examine their molecular mechanisms of allosteric regulation. The comparative intrinsic and conformational dynamics analysis of three DNMT1 systems has highlighted the pivotal role of the RFTS domain as the dynamics hub in both intra- and inter-molecular interactions. The site perturbation and network modeling approaches have revealed the different and more complex allosteric interaction landscape in both DNMT1 complexes, involving the events caused by mutational hotspots and post-translation modification sites through protein-protein interactions (PPIs). Furthermore, communication pathway analysis and pocket detection have provided new mechanistic insights into molecular mechanisms underlying quaternary structures of DNMT1 complexes, suggesting potential targeting pockets for PPI-based allosteric drug design.

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.

Anti-apoptotic BCL-2 regulation by changes in dynamics of its long unstructured loop

BCL-2, a key protein in inhibiting apoptosis, has a 65-residue-long highly flexible loop domain (FLD) located on the opposite side of its ligand-binding groove. In vivo phosphorylation of the FLD enhances the affinity of BCL-2 for pro-apoptotic ligands, and consequently anti-apoptotic activity. However, it remains unknown as to how the faraway, unstructured FLD modulates the affinity. Here we investigate the protein-ligand interactions by fluorescence techniques and monitor protein dynamics by DEER and NMR spectroscopy tools. We show that phosphomimetic mutations on the FLD lead to a reduction in structural flexibility, hence promoting ligand access to the groove. The bound pro-apoptotic ligands can be displaced by the BCL-2-selective inhibitor ABT-199 efficiently, and thus released to trigger apoptosis. We show that changes in structural flexibility on an unstructured loop can activate an allosteric protein that is otherwise structurally inactive.

Molecular Simulations and Network Modeling Reveal an Allosteric Signaling in the SARS-CoV-2 Spike Proteins

The development of computational strategies for the quantitative characterization of the functional mechanisms of SARS-CoV-2 spike proteins is of paramount importance in efforts to accelerate the discovery of novel therapeutic agents and vaccines combating the COVID-19 pandemic. Structural and biophysical studies have recently characterized the conformational landscapes of the SARS-CoV-2 spike glycoproteins in the prefusion form, revealing a spectrum of stable and more dynamic states. By employing molecular simulations and network modeling approaches, this study systematically examined functional dynamics and identified the regulatory centers of allosteric interactions for distinct functional states of the wild-type and mutant variants of the SARS-CoV-2 prefusion spike trimer. This study presents evidence that the SARS-CoV-2 spike protein can function as an allosteric regulatory engine that fluctuates between dynamically distinct functional states. Perturbation-based modeling of the interaction networks revealed a key role of the cross-talk between the effector hotspots in the receptor binding domain and the fusion peptide proximal region of the SARS-CoV-2 spike protein. The results have shown that the allosteric hotspots of the interaction networks in the SARS-CoV-2 spike protein can control the dynamic switching between functional conformational states that are associated with virus entry to the host receptor. This study offers a useful and novel perspective on the underlying mechanisms of the SARS-CoV-2 spike protein through the lens of allosteric signaling as a regulatory apparatus of virus transmission that could open up opportunities for targeted allosteric drug discovery against SARS-CoV-2 proteins and contribute to the rapid response to the current and potential future pandemic scenarios.

Allosteric regulation of U1 snRNP by splicing regulatory proteins controls spliceosomal assembly

Alternative splicing is responsible for much of the transcriptomic and proteomic diversity observed in eukaryotes and involves combinatorial regulation by many cis-acting elements and trans-acting factors. SR and hnRNP splicing regulatory proteins often have opposing effects on splicing efficiency depending on where they bind the pre-mRNA relative to the splice site. Position-dependent splicing repression occurs at spliceosomal E-complex, suggesting that U1 snRNP binds but cannot facilitate higher order spliceosomal assembly. To test the hypothesis that the structure of U1 snRNA changes during activation or repression, we developed a method to structure-probe native U1 snRNP in enriched conformations that mimic activated or repressed spliceosomal E-complexes. While the core of U1 snRNA is highly structured, the 5' end of U1 snRNA shows different SHAPE reactivities and psoralen crosslinking efficiencies depending on where splicing regulatory elements are located relative to the 5' splice site. A motif within the 5' splice site binding region of U1 snRNA is more reactive toward SHAPE electrophiles when repressors are bound, suggesting U1 snRNA is bound, but less base-paired. These observations demonstrate that splicing regulators modulate splice site selection allosterically.

Cooperative Binding of the Cationic Porphyrin Tris-T4 Enhances Catalytic Activity of 20S Proteasome Unveiling a Complex Distribution of Functional States

The present study provides new evidence that cationic porphyrins may be considered as tunable platforms to interfere with the structural “key code” present on the 20S proteasome α-rings and, by consequence, with its catalytic activity. Here, we describe the functional and conformational effects on the 20S proteasome induced by the cooperative binding of the tri-cationic 5-(phenyl)-10,15,20-(tri N-methyl-4-pyridyl) porphyrin (Tris-T4). Our integrated kinetic, NMR, and in silico analysis allowed us to disclose a complex effect on the 20S catalytic activity depending on substrate/porphyrin concentration. The analysis of the kinetic data shows that Tris-T4 shifts the relative populations of the multiple interconverting 20S proteasome conformations leading to an increase in substrate hydrolysis by an allosteric pathway. Based on our Tris-T4/h20S interaction model, Tris-T4 is able to affect gating dynamics and substrate hydrolysis by binding to an array of negatively charged and hydrophobic residues present on the protein surface involved in the 20S molecular activation by the regulatory proteins (RPs). Accordingly, despite the fact that Tris-T4 also binds to the α3ΔN mutant, allosteric modulation is not observed since the molecular mechanism connecting gate dynamics with substrate hydrolysis is impaired. We envisage that the dynamic view of the 20S conformational equilibria, activated through cooperative Tris-T4 binding, may work as a simplified model for a better understanding of the intricate network of 20S conformational/functional states that may be mobilized by exogenous ligands, paving the way for the development of a new generation of proteasome allosteric modulators.

Allosteric regulation of lysosomal enzyme recognition by the cation-independent mannose 6-phosphate receptor

The cation-independent mannose 6-phosphate receptor (CI-MPR, IGF2 receptor or CD222), is a multifunctional glycoprotein required for normal development. Through the receptor's ability to bind unrelated extracellular and intracellular ligands, it participates in numerous functions including protein trafficking, lysosomal biogenesis, and regulation of cell growth. Clinically, endogenous CI-MPR delivers infused recombinant enzymes to lysosomes in the treatment of lysosomal storage diseases. Although four of the 15 domains comprising CI-MPR's extracellular region bind phosphorylated glycans on lysosomal enzymes, knowledge of how CI-MPR interacts with ~60 different lysosomal enzymes is limited. Here, we show by electron microscopy and hydroxyl radical protein footprinting that the N-terminal region of CI-MPR undergoes dynamic conformational changes as a consequence of ligand binding and different pH conditions. These data, coupled with X-ray crystallography, surface plasmon resonance and molecular modeling, allow us to propose a model explaining how high-affinity carbohydrate binding is achieved through allosteric domain cooperativity.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Binding and allosteric transmission of histone H3 Lys-4 trimethylation to the recombinase RAG-1 are separable functions of the RAG-2 plant homeodomain finger

V(D)J recombination is initiated by the recombination-activating gene protein (RAG) recombinase, consisting of RAG-1 and RAG-2 subunits. The susceptibility of gene segments to cleavage by RAG is associated with gene transcription and with epigenetic marks characteristic of active chromatin, including histone H3 trimethylated at lysine 4 (H3K4me3). Binding of H3K4me3 by a plant homeodomain (PHD) in RAG-2 induces conformational changes in RAG-1, allosterically stimulating substrate binding and catalysis. To better understand the path of allostery from the RAG-2 PHD finger to RAG-1, here we employed phylogenetic substitution. We observed that a chimeric RAG-2 protein in which the mouse PHD finger is replaced by the corresponding domain from the shark Chiloscyllium punctatum binds H3K4me3 but fails to transmit an allosteric signal, indicating that binding of H3K4me3 by RAG-2 is insufficient to support recombination. By substituting residues in the C. punctatum PHD with the corresponding residues in the mouse PHD and testing for rescue of allostery, we demonstrate that H3K4me3 binding and transmission of an allosteric signal to RAG-1 are separable functions of the RAG-2 PHD finger.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

Allosteric Cooperativity in Proton Energy Conversion in A1-Type Cytochrome c Oxidase

Cytochrome c oxidase (CcO), the Cu_A, heme a, heme a₃, Cu_B enzyme of respiratory chain, converts the free energy released by aerobic cytochrome c oxidation into a membrane electrochemical proton gradient (ΔμH⁺). ΔμH⁺ derives from the membrane anisotropic arrangement of dioxygen reduction to two water molecules and transmembrane proton pumping from a negative (N) space to a positive (P) space separated by the membrane. Spectroscopic, potentiometric, and X-ray crystallographic analyses characterize allosteric cooperativity of dioxygen binding and reduction with protonmotive conformational states of CcO. These studies show that allosteric cooperativity stabilizes the favorable conformational state for conversion of redox energy into a transmembrane ΔμH⁺.

Two PKA RIα holoenzyme states define ATP as an isoform-specific orthosteric inhibitor that competes with the allosteric activator, cAMP

Protein kinase A (PKA) holoenzyme, comprised of a cAMP-binding regulatory (R)-subunit dimer and 2 catalytic (C)-subunits, is the master switch for cAMP-mediated signaling. Of the 4 R-subunits (RIα, RIβ, RIIα, RIIβ), RIα is most essential for regulating PKA activity in cells. Our 2 RIα2C2 holoenzyme states, which show different conformations with and without ATP, reveal how ATP/Mg2+ functions as a negative orthosteric modulator. Biochemical studies demonstrate how the removal of ATP primes the holoenzyme for cAMP-mediated activation. The opposing competition between ATP/cAMP is unique to RIα. In RIIβ, ATP serves as a substrate and facilitates cAMP-activation. The isoform-specific RI-holoenzyme dimer interface mediated by N3A-N3A' motifs defines multidomain cross-talk and an allosteric network that creates competing roles for ATP and cAMP. Comparisons to the RIIβ holoenzyme demonstrate isoform-specific holoenzyme interfaces and highlights distinct allosteric mechanisms for activation in addition to the structural diversity of the isoforms.

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.

Perspective: The promises of a holistic view of proteins-impact on antibody engineering and drug discovery

The reductionist approach is prevalent in biomedical science. However, increasing evidence now shows that biological systems cannot be simply considered as the sum of its parts. With experimental, technological, and computational advances, we can now do more than view parts in isolation, thus we propose that an increasing holistic view (where a protein is investigated as much as a whole as possible) is now timely. To further advocate this, we review and discuss several studies and applications involving allostery, where distant protein regions can cross-talk to influence functionality. Therefore, we believe that an increasing big picture approach holds great promise, particularly in the areas of antibody engineering and drug discovery in rational drug design.

Directed Evolution Mimics Allosteric Activation by Stepwise Tuning of the Conformational Ensemble

Allosteric enzymes contain a wealth of catalytic diversity that remains distinctly underutilized for biocatalysis. Tryptophan synthase is a model allosteric system and a valuable enzyme for the synthesis of noncanonical amino acids (ncAA). Previously, we evolved the β-subunit from Pyrococcus furiosus, PfTrpB, for ncAA synthase activity in the absence of its native partner protein PfTrpA. However, the precise mechanism by which mutation activated TrpB to afford a stand-alone catalyst remained enigmatic. Here, we show that directed evolution caused a gradual change in the rate-limiting step of the catalytic cycle. Concomitantly, the steady-state distribution of the intermediates shifts to favor covalently bound Trp adducts, which have increased thermodynamic stability. The biochemical properties of these evolved, stand-alone TrpBs converge on those induced in the native system by allosteric activation. High-resolution crystal structures of the wild-type enzyme, an intermediate in the lineage, and the final variant, encompassing five distinct chemical states, show that activating mutations have only minor structural effects on their immediate environment. Instead, mutation stabilizes the large-scale motion of a subdomain to favor an otherwise transiently populated closed conformational state. This increase in stability enabled the first structural description of Trp covalently bound in a catalytically active TrpB, confirming key features of catalysis. These data combine to show that sophisticated models of allostery are not a prerequisite to recapitulating its complex effects via directed evolution, opening the way to engineering stand-alone versions of diverse allosteric enzymes.

Biased signaling of the proton-sensing receptor OGR1 by benzodiazepines

GPCRs have diverse signaling capabilities, based on their ability to assume various conformations. Moreover, it is now appreciated that certain ligands can promote distinct receptor conformations and thereby bias signaling toward a specific pathway to differentially affect cell function. The recently deorphanized G protein-coupled receptor OGR1 [ovarian cancer G protein-coupled receptor 1 ( GPR68)] exhibits diverse signaling events when stimulated by reductions in extracellular pH. We recently demonstrated airway smooth muscle cells transduce multiple signaling events, reflecting a diverse capacity to couple to multiple G proteins. Moreover, we recently discovered that the benzodiazepine lorazepam, more commonly recognized as an agonist of the γ-aminobutyric acid A (GABAA) receptor, can function as an allosteric modulator of OGR1 and, similarly, can promote multiple signaling events. In this study, we demonstrated that different benzodiazepines exhibit a range of biases for OGR1, with sulazepam selectively activating the canonical Gs of the G protein signaling pathway, in heterologous expression systems, as well as in several primary cell types. These findings highlight the potential power of biased ligand pharmacology for manipulating receptor signaling qualitatively, to preferentially activate pathways that are therapeutically beneficial.-Pera, T., Deshpande, D. A., Ippolito, M., Wang, B., Gavrila, A., Michael, J. V., Nayak, A. P., Tompkins, E., Farrell, E., Kroeze, W. K., Roth, B. L., Panettieri, R. A. Jr Benovic, J. L., An, S. S., Dulin, N. O., Penn, R. B. Biased signaling of the proton-sensing receptor OGR1 by benzodiazepines.

Allosteric control of a bacterial stress response system by an anti-σ factor

Bacterial signal transduction systems commonly use receiver (REC) domains, which regulate adaptive responses to the environment as a function of their phosphorylation state. REC domains control cell physiology through diverse mechanisms, many of which remain understudied. We have defined structural features that underlie activation of the multi-domain REC protein, PhyR, which functions as an anti-anti-σ factor and regulates transcription of genes required for stress adaptation and host-microbe interactions in Alphaproteobacteria. Though REC phosphorylation is necessary for PhyR function in vivo, we did not detect expected changes in inter-domain interactions upon phosphorylation by solution X-ray scattering. We sought to understand this result by defining additional molecular requirements for PhyR activation. We uncovered specific interactions between unphosphorylated PhyR and an intrinsically disordered region (IDR) of the anti-σ factor, NepR, by solution NMR spectroscopy. Our data support a model whereby nascent NepR(IDR)-PhyR interactions and REC phosphorylation coordinately impart the free energy to shift PhyR to an open, active conformation that binds and inhibits NepR. This mechanism ensures PhyR is activated only when NepR and an activating phosphoryl signal are present. Our study provides new structural understanding of the molecular regulatory logic underlying a conserved environmental response system.

Using competition assays to quantitatively model cooperative binding by transcription factors and other ligands

BACKGROUND

The affinities of DNA binding proteins for target sites can be used to model the regulation of gene expression. These proteins can bind to DNA cooperatively, strongly impacting their affinity and specificity. However, current methods for measuring cooperativity do not provide the means to accurately predict binding behavior over a wide range of concentrations.

METHODS

We use standard computational and mathematical methods, and develop novel methods as described in Results.

RESULTS

We explore some complexities of cooperative binding, and develop an improved method for relating in vitro measurements to in vivo function, based on ternary complex formation. We derive expressions for the equilibria among the various complexes, and explore the limitations of binding experiments that model the system using a single parameter. We describe how to use single-ligand binding and ternary complex formation in tandem to determine parameters that have thermodynamic relevance. We develop an improved method for finding both single-ligand dissociation constants and concentrations simultaneously. We show how the cooperativity factor can be found when only one of the single-ligand dissociation constants can be measured.

CONCLUSIONS

The methods that we develop constitute an optimized approach to accurately model cooperative binding.

GENERAL SIGNIFICANCE

The expressions and methods we develop for modeling and analyzing DNA binding and cooperativity are applicable to most cases where multiple ligands bind to distinct sites on a common substrate. The parameters determined using these methods can be fed into models of higher-order cooperativity to increase their predictive power.

Site-to-site interdomain communication may mediate different loss-of-function mechanisms in a cancer-associated NQO1 polymorphism

Disease associated genetic variations often cause intracellular enzyme inactivation, dysregulation and instability. However, allosteric communication of mutational effects to distant functional sites leading to loss-of-function remains poorly understood. We characterize here interdomain site-to-site communication by which a common cancer-associated single nucleotide polymorphism (c.C609T/p.P187S) reduces the activity and stability in vivo of NAD(P)H:quinone oxidoreductase 1 (NQO1). NQO1 is a FAD-dependent, two-domain multifunctional stress protein acting as a Phase II enzyme, activating cancer pro-drugs and stabilizing p53 and p73α oncosuppressors. We show that p.P187S causes structural and dynamic changes communicated to functional sites far from the mutated site, affecting the FAD binding site located at the N-terminal domain (NTD) and accelerating proteasomal degradation through dynamic effects on the C-terminal domain (CTD). Structural protein:protein interaction studies reveal that the cancer-associated polymorphism does not abolish the interaction with p73α, indicating that oncosuppressor destabilization largely mirrors the low intracellular stability of p.P187S. In conclusion, we show how a single disease associated amino acid change may allosterically perturb several functional sites in an oligomeric and multidomain protein. These results have important implications for the understanding of loss-of-function genetic diseases and the identification of novel structural hot spots as targets for pharmacological intervention.

Thermal activation of 'allosteric-like' large-scale motions in a eukaryotic Lactate Dehydrogenase

Conformational changes occurring during the enzymatic turnover are essential for the regulation of protein functionality. Individuating the protein regions involved in these changes and the associated mechanical modes is still a challenge at both experimental and theoretical levels. We present here a detailed investigation of the thermal activation of the functional modes and conformational changes in a eukaryotic Lactate Dehydrogenase enzyme (LDH). Neutron Spin Echo spectroscopy and Molecular Dynamics simulations were used to uncover the characteristic length- and timescales of the LDH nanoscale motions in the apo state. The modes involving the catalytic loop and the mobile region around the binding site are activated at room temperature, and match the allosteric reorganisation of bacterial LDHs. In a temperature window of about 15 degrees, these modes render the protein flexible enough and capable of reorganising the active site toward reactive configurations. On the other hand an excess of thermal excitation leads to the distortion of the protein matrix with a possible anti-catalytic effect. Thus, the temperature activates eukaryotic LDHs via the same conformational changes observed in the allosteric bacterial LDHs. Our investigation provides an extended molecular picture of eukaryotic LDH's conformational landscape that enriches the static view based on crystallographic studies alone.

Computational Analysis of Residue Interaction Networks and Coevolutionary Relationships in the Hsp70 Chaperones: A Community-Hopping Model of Allosteric Regulation and Communication

Allosteric interactions in the Hsp70 proteins are linked with their regulatory mechanisms and cellular functions. Despite significant progress in structural and functional characterization of the Hsp70 proteins fundamental questions concerning modularity of the allosteric interaction networks and hierarchy of signaling pathways in the Hsp70 chaperones remained largely unexplored and poorly understood. In this work, we proposed an integrated computational strategy that combined atomistic and coarse-grained simulations with coevolutionary analysis and network modeling of the residue interactions. A novel aspect of this work is the incorporation of dynamic residue correlations and coevolutionary residue dependencies in the construction of allosteric interaction networks and signaling pathways. We found that functional sites involved in allosteric regulation of Hsp70 may be characterized by structural stability, proximity to global hinge centers and local structural environment that is enriched by highly coevolving flexible residues. These specific characteristics may be necessary for regulation of allosteric structural transitions and could distinguish regulatory sites from nonfunctional conserved residues. The observed confluence of dynamics correlations and coevolutionary residue couplings with global networking features may determine modular organization of allosteric interactions and dictate localization of key mediating sites. Community analysis of the residue interaction networks revealed that concerted rearrangements of local interacting modules at the inter-domain interface may be responsible for global structural changes and a population shift in the DnaK chaperone. The inter-domain communities in the Hsp70 structures harbor the majority of regulatory residues involved in allosteric signaling, suggesting that these sites could be integral to the network organization and coordination of structural changes. Using a network-based formalism of allostery, we introduced a community-hopping model of allosteric communication. Atomistic reconstruction of signaling pathways in the DnaK structures captured a direction-specific mechanism and molecular details of signal transmission that are fully consistent with the mutagenesis experiments. The results of our study reconciled structural and functional experiments from a network-centric perspective by showing that global properties of the residue interaction networks and coevolutionary signatures may be linked with specificity and diversity of allosteric regulation mechanisms.

The diversity of allosteric mechanisms in the Hsp70 proteins could range from modulation of the inter-domain interactions and conformational dynamics to fine-tuning of the Hsp70 interactions with co-chaperones. The goal of this study is to present a systematic computational analysis of the dynamic and evolutionary factors underlying allosteric structural transformations of the Hsp70 proteins. We investigated the relationship between functional dynamics, residue coevolution, and network organization of residue interactions in the Hsp70 proteins. The results of this study revealed that conformational dynamics of the Hsp70 proteins may be linked with coevolutionary propensities and mutual information dependencies of the protein residues. Modularity and connectivity of allosteric interactions in the Hsp70 chaperones are coordinated by stable functional sites that feature unique coevolutionary signatures and high network centrality. The emergence of the inter-domain communities that are coordinated by functional centers and include highly coevolving residues could facilitate structural transitions through cooperative reorganization of the local interacting modules. We determined that the differences in the modularity of the residue interactions and organization of coevolutionary networks in DnaK may be associated with variations in their allosteric mechanisms. The network signatures of the DnaK structures are characteristic of a population-shift allostery that allows for coordinated structural rearrangements of local communities. A dislocation of mediating centers and insufficient coevolutionary coupling between functional regions may render a reduced cooperativity and promote a limited entropy-driven allostery in the Sse1 chaperone that occurs without structural changes. The results of this study showed that a network-centric framework and a community-hopping model of allosteric communication pathways may provide novel insights into molecular and evolutionary principles of allosteric regulation in the Hsp70 proteins.

RET Functions as a Dual-Specificity Kinase that Requires Allosteric Inputs from Juxtamembrane Elements

Receptor tyrosine kinases exhibit a variety of activation mechanisms despite highly homologous catalytic domains. Such diversity arises through coupling of extracellular ligand-binding portions with highly variable intracellular sequences flanking the tyrosine kinase domain and specific patterns of autophosphorylation sites. Here, we show that the juxtamembrane (JM) segment enhances RET catalytic domain activity through Y687. This phospho-site is also required by the JM region to rescue an otherwise catalytically deficient RET activation-loop mutant lacking tyrosines. Structure-function analyses identified interactions between the JM hinge, αC helix, and an unconventional activation-loop serine phosphorylation site that engages the HRD motif and promotes phospho-tyrosine conformational accessibility and regulatory spine assembly. We demonstrate that this phospho-S909 arises from an intrinsic RET dual-specificity kinase activity and show that an equivalent serine is required for RET signaling in Drosophila. Our findings reveal dual-specificity and allosteric components for the mechanism of RET activation and signaling with direct implications for drug discovery.

M1 muscarinic allosteric modulators slow prion neurodegeneration and restore memory loss

The current frontline symptomatic treatment for Alzheimer's disease (AD) is whole-body upregulation of cholinergic transmission via inhibition of acetylcholinesterase. This approach leads to profound dose-related adverse effects. An alternative strategy is to selectively target muscarinic acetylcholine receptors, particularly the M1 muscarinic acetylcholine receptor (M1 mAChR), which was previously shown to have procognitive activity. However, developing M1 mAChR-selective orthosteric ligands has proven challenging. Here, we have shown that mouse prion disease shows many of the hallmarks of human AD, including progressive terminal neurodegeneration and memory deficits due to a disruption of hippocampal cholinergic innervation. The fact that we also show that muscarinic signaling is maintained in both AD and mouse prion disease points to the latter as an excellent model for testing the efficacy of muscarinic pharmacological entities. The memory deficits we observed in mouse prion disease were completely restored by treatment with benzyl quinolone carboxylic acid (BQCA) and benzoquinazoline-12 (BQZ-12), two highly selective positive allosteric modulators (PAMs) of M1 mAChRs. Furthermore, prolonged exposure to BQCA markedly extended the lifespan of diseased mice. Thus, enhancing hippocampal muscarinic signaling using M1 mAChR PAMs restored memory loss and slowed the progression of mouse prion disease, indicating that this ligand type may have clinical benefit in diseases showing defective cholinergic transmission, such as AD.

WD40 domain of Apc1 is critical for the coactivator-induced allosteric transition that stimulates APC/C catalytic activity

The anaphase-promoting complex/cyclosome (APC/C) is a large multimeric cullin-RING E3 ubiquitin ligase that orchestrates cell-cycle progression by targeting cell-cycle regulatory proteins for destruction via the ubiquitin proteasome system. The APC/C assembly comprises two scaffolding subcomplexes: the platform and the TPR lobe that together coordinate the juxtaposition of the catalytic and substrate-recognition modules. The platform comprises APC/C subunits Apc1, Apc4, Apc5, and Apc15. Although the role of Apc1 as an APC/C scaffolding subunit has been characterized, its specific functions in contributing toward APC/C catalytic activity are not fully understood. Here, we report the crystal structure of the N-terminal domain of human Apc1 (Apc1N) determined at 2.2-Å resolution and provide an atomic-resolution description of the architecture of its WD40 (WD40 repeat) domain (Apc1(WD40)). To understand how Apc1(WD40) contributes to APC/C activity, a mutant form of the APC/C with Apc1(WD40) deleted was generated and evaluated biochemically and structurally. We found that the deletion of Apc1(WD40) abolished the UbcH10-dependent ubiquitination of APC/C substrates without impairing the Ube2S-dependent ubiquitin chain elongation activity. A cryo-EM structure of an APC/C-Cdh1 complex with Apc1(WD40) deleted showed that the mutant APC/C is locked into an inactive conformation in which the UbcH10-binding site of the catalytic module is inaccessible. Additionally, an EM density for Apc15 is not visible. Our data show that Apc1(WD40) is required to mediate the coactivator-induced conformational change of the APC/C that is responsible for stimulating APC/C catalytic activity by promoting UbcH10 binding. In contrast, Ube2S activity toward APC/C substrates is not dependent on the initiation-competent conformation of the APC/C.

Metabolic Actions of the Type 1 Cholecystokinin Receptor: Its Potential as a Therapeutic Target

Cholecystokinin (CCK) regulates appetite and reduces food intake by activating the type 1 CCK receptor (CCK1R). Attempts to develop CCK1R agonists for obesity have yielded active agents that have not reached clinical practice. Here we discuss why, along with new strategies to target CCK1R more effectively. We examine signaling events and the possibility of developing agents that exhibit ligand-directed bias, to dissociate satiety activity from undesirable side effects. Potential allosteric sites of modulation are also discussed, along with desired properties of a positive allosteric modulator (PAM) without intrinsic agonist action as another strategy to treat obesity. These new types of CCK1R-active drugs could be useful as standalone agents or as part of a rational drug combination for management of obesity.

Two Deafness-Causing Actin Mutations (DFNA20/26) Have Allosteric Effects on the Actin Structure

Point mutations in γ-cytoplasmic actin have been shown to result in autosomal-dominant, nonsyndromic, early-onset deafness. Two mutations at the same site, K118M and K118N, provide a unique opportunity to compare the effects of two dissimilar amino acid substitutions that produce a similar phenotype in humans. K118 resides in a helix that runs from K113 to T126, and mutations that alter the position, dynamics, and/or biochemistry of this helix can result in a wide range of pathologies. Using a combination of computational and experimental studies, both employing yeast actin, we find that these mutations at K118 result in changes in the structure and dynamics of the DNase-I loop, alterations in the structure of the H73 loop as well as the side-chain orientations of W79 and W86, changes in nucleotide exchange rates, and significant shifts in the twist of the actin monomer. Interestingly, in the case of K118N, the twist of the monomer is nearly identical to that of the F-actin protomer, and in vitro polymerization assays show that this mutation results in faster polymerization. Taken together, these results indicate that mutations at this site give rise to a series of small changes that can be tolerated in vivo but result in misregulation of actin assembly and dynamics.

Allosteric Regulation of Fibronectin/α5β1 Interaction by Fibronectin-Binding MSCRAMMs

Adherence of microbes to host tissues is a hallmark of infectious disease and is often mediated by a class of adhesins termed MSCRAMMs (Microbial Surface Components Recognizing Adhesive Matrix Molecules). Numerous pathogens express MSCRAMMs that specifically bind the heterodimeric human glycoprotein fibronectin (Fn). In addition to roles in adhesion, Fn-binding MSCRAMMs exploit physiological Fn functions. For example, several pathogens can invade host cells by a mechanism whereby MSCRAMM-bound Fn bridges interaction with α₅β₁ integrin. Here, we investigate two Fn-binding MSCRAMMs, FnBPA (Staphylococcus aureus) and BBK32 (Borrelia burgdorferi) to probe structure-activity relationships of MSCRAMM-induced Fn/α₅β₁integrin activation. Circular dichroism, fluorescence resonance energy transfer, and dynamic light scattering techniques uncover a conformational rearrangement of Fn involving domains distant from the MSCRAMM binding site. Surface plasmon resonance experiments demonstrate a significant enhancement of Fn/α₅β₁ integrin affinity in the presence of FnBPA or BBK32. Detailed kinetic analysis of these interactions reveal that this change in affinity can be attributed solely to an increase in the initial Fn/α₅β₁ on-rate and that this rate-enhancement is dependent on high-affinity Fn-binding by MSCRAMMs. These data implicate MSCRAMM-induced perturbation of specific intramolecular contacts within the Fn heterodimer resulting in activation by exposing previously cryptic α₅β₁ interaction motifs. By correlating structural changes in Fn to a direct measurement of increased Fn/α₅β₁ affinity, this work significantly advances our understanding of the structural basis for the modulation of integrin function by Fn-binding MSCRAMMs.

Structural Mechanisms of Peptide Recognition and Allosteric Modulation of Gene Regulation by the RRNPP Family of Quorum-Sensing Regulators

The members of RRNPP family of bacterial regulators sense population density-specific secreted oligopeptides and modulate the expression of genes involved in cellular processes, such as sporulation, competence, virulence, biofilm formation, conjugative plasmid transfer and antibiotic resistance. Signaling by RRNPP regulators include several steps: generation and secretion of the signaling oligopeptides, re-internalization of the signaling molecules into the cytoplasm, signal sensing by the cytosolic RRNPP regulators, signal-specific allosteric structural changes in the regulators, and interaction of the regulators with their respective regulatory target and gene regulation. The recently determined structures of the RRNPP regulators provide insight into the mechanistic aspects for several steps in this signaling circuit. In this review, we discuss the structural principles underlying peptide specificity, regulatory target recognition, and ligand-induced allostery in RRNPP regulators and its impact on gene regulation. Despite the conserved tertiary structure of these regulators, structural analyses revealed unexpected diversity in the mechanism of activation and molecular strategies that couple the peptide-induced allostery to gene regulation. Although these structural studies provide a sophisticated understanding of gene regulation by RRNPP regulators, much needs to be learned regarding the target DNA binding by yet-to-be characterized RNPP regulators and the several aspects of signaling by Rgg regulators.

mGlu5 positive allosteric modulation normalizes synaptic plasticity defects and motor phenotypes in a mouse model of Rett syndrome

Rett syndrome (RS) is a neurodevelopmental disorder that shares many symptomatic and pathological commonalities with idiopathic autism. Alterations in protein synthesis-dependent synaptic plasticity (PSDSP) are a hallmark of a number of syndromic forms of autism; in the present work, we explore the consequences of disruption and rescue of PSDSP in a mouse model of RS. We report that expression of a key regulator of synaptic protein synthesis, the metabotropic glutamate receptor 5 (mGlu5) protein, is significantly reduced in both the brains of RS model mice and in the motor cortex of human RS autopsy samples. Furthermore, we demonstrate that reduced mGlu5 expression correlates with attenuated DHPG-induced long-term depression in the hippocampus of RS model mice, and that administration of a novel mGlu5 positive allosteric modulator (PAM), termed VU0462807, can rescue synaptic plasticity defects. Additionally, treatment of Mecp2-deficient mice with VU0462807 improves motor performance (open-field behavior and gait dynamics), corrects repetitive clasping behavior, as well as normalizes cued fear-conditioning defects. Importantly, due to the rationale drug discovery approach used in its development, our novel mGlu5 PAM improves RS phenotypes and synaptic plasticity defects without evoking the overt adverse effects commonly associated with potentiation of mGlu5 signaling (i.e. seizures), or affecting cardiorespiratory defects in RS model mice. These findings provide strong support for the continued development of mGlu5 PAMs as potential therapeutic agents for use in RS, and, more broadly, for utility in idiopathic autism.

Human Secreted Ly-6/uPAR Related Protein-1 (SLURP-1) Is a Selective Allosteric Antagonist of α7 Nicotinic Acetylcholine Receptor

SLURP-1 is a secreted toxin-like Ly-6/uPAR protein found in epithelium, sensory neurons and immune cells. Point mutations in the slurp-1 gene cause the autosomal inflammation skin disease Mal de Meleda. SLURP-1 is considered an autocrine/paracrine hormone that regulates growth and differentiation of keratinocytes and controls inflammation and malignant cell transformation. The majority of previous studies of SLURP-1 have been made using fusion constructs containing, in addition to the native protein, extra polypeptide sequences. Here we describe the activity and pharmacological profile of a recombinant analogue of human SLURP-1 (rSLURP-1) differing from the native protein only by one additional N-terminal Met residue. rSLURP-1 significantly inhibited proliferation (up to ~ 40%, EC₅₀ ~ 4 nM) of human oral keratinocytes (Het-1A cells). Application of mecamylamine and atropine,—non-selective inhibitors of nicotinic acetylcholine receptors (nAChRs) and muscarinic acetylcholine receptors, respectively, and anti-α7-nAChRs antibodies revealed α7 type nAChRs as an rSLURP-1 target in keratinocytes. Using affinity purification from human cortical extracts, we confirmed that rSLURP-1 binds selectively to the α7-nAChRs. Exposure of Xenopus oocytes expressing α7-nAChRs to rSLURP-1 caused a significant non-competitive inhibition of the response to acetylcholine (up to ~ 70%, IC₅₀ ~ 1 μM). It was shown that rSLURP-1 binds to α7-nAChRs overexpressed in GH₄C_l cells, but does not compete with ¹²⁵I-α-bungarotoxin for binding to the receptor. These findings imply an allosteric antagonist-like mode of SLURP-1 interaction with α7-nAChRs outside the classical ligand-binding site. Contrary to rSLURP-1, other inhibitors of α7-nAChRs (mecamylamine, α-bungarotoxin and Lynx1) did not suppress the proliferation of keratinocytes. Moreover, the co-application of α-bungarotoxin with rSLURP-1 did not influence antiproliferative activity of the latter. This supports the hypothesis that the antiproliferative activity of SLURP-1 is related to ‘metabotropic’ signaling pathway through α7-nAChR, that activates intracellular signaling cascades without opening the receptor channel.

Positive and Negative Allosteric Modulation of an α1β3γ2 γ-Aminobutyric Acid Type A (GABAA) Receptor by Binding to a Site in the Transmembrane Domain at the γ+-β- Interface

In the process of developing safer general anesthetics, isomers of anesthetic ethers and barbiturates have been discovered that act as convulsants and inhibitors of γ-aminobutyric acid type A receptors (GABAARs) rather than potentiators. It is unknown whether these convulsants act as negative allosteric modulators by binding to the intersubunit anesthetic-binding sites in the GABAAR transmembrane domain (Chiara, D. C., Jayakar, S. S., Zhou, X., Zhang, X., Savechenkov, P. Y., Bruzik, K. S., Miller, K. W., and Cohen, J. B. (2013) J. Biol. Chem. 288, 19343-19357) or to known convulsant sites in the ion channel or extracellular domains. Here, we show that S-1-methyl-5-propyl-5-(m-trifluoromethyl-diazirynylphenyl) barbituric acid (S-mTFD-MPPB), a photoreactive analog of the convulsant barbiturate S-MPPB, inhibits α1β3γ2 but potentiates α1β3 GABAAR responses. In the α1β3γ2 GABAAR, S-mTFD-MPPB binds in the transmembrane domain with high affinity to the γ(+)-β(-) subunit interface site with negative energetic coupling to GABA binding in the extracellular domain at the β(+)-α(-) subunit interfaces. GABA inhibits S-[(3)H]mTFD-MPPB photolabeling of γ2Ser-280 (γM2-15') in this site. In contrast, within the same site GABA enhances photolabeling of β3Met-227 in βM1 by an anesthetic barbiturate, R-[(3)H]methyl-5-allyl-5-(m-trifluoromethyl-diazirynylphenyl)barbituric acid (mTFD-MPAB), which differs from S-mTFD-MPPB in structure only by chirality and two hydrogens (propyl versus allyl). S-mTFD-MPPB and R-mTFD-MPAB are predicted to bind in different orientations at the γ(+)-β(-) site, based upon the distance in GABAAR homology models between γ2Ser-280 and β3Met-227. These results provide an explanation for S-mTFD-MPPB inhibition of α1β3γ2 GABAAR function and provide a first demonstration that an intersubunit-binding site in the GABAAR transmembrane domain binds negative and positive allosteric modulators.

Comparative sequence analysis suggests a conserved gating mechanism for TRP channels

The transient receptor potential (TRP) channel superfamily plays a central role in transducing diverse sensory stimuli in eukaryotes. Although dissimilar in sequence and domain organization, all known TRP channels act as polymodal cellular sensors and form tetrameric assemblies similar to those of their distant relatives, the voltage-gated potassium (Kv) channels. Here, we investigated the related questions of whether the allosteric mechanism underlying polymodal gating is common to all TRP channels, and how this mechanism differs from that underpinning Kv channel voltage sensitivity. To provide insight into these questions, we performed comparative sequence analysis on large, comprehensive ensembles of TRP and Kv channel sequences, contextualizing the patterns of conservation and correlation observed in the TRP channel sequences in light of the well-studied Kv channels. We report sequence features that are specific to TRP channels and, based on insight from recent TRPV1 structures, we suggest a model of TRP channel gating that differs substantially from the one mediating voltage sensitivity in Kv channels. The common mechanism underlying polymodal gating involves the displacement of a defect in the H-bond network of S6 that changes the orientation of the pore-lining residues at the hydrophobic gate.

Subtype-dependent N-methyl-D-aspartate receptor amino-terminal domain conformations and modulation by spermine

The N-methyl-d-aspartate (NMDA) subtype of the ionotropic glutamate receptors is the primary mediator of calcium-permeable excitatory neurotransmission in the central nervous system. Subunit composition and binding of allosteric modulators to the amino-terminal domain determine the open probability of the channel. By using luminescence resonance energy transfer with functional receptors expressed in CHO cells, we show that the cleft of the amino-terminal domain of the GluN2B subunit, which has a lower channel open probability, is on average more closed than the GluN2A subunit, which has a higher open probability. Furthermore, the GluN1 amino-terminal domain adopts a more open conformation when coassembled with GluN2A than with GluN2B. Binding of spermine, an allosteric potentiator, opens the amino-terminal domain cleft of both the GluN2B subunit and the adjacent GluN1 subunit. These studies provide direct structural evidence that the inherent conformations of the amino-terminal domains vary based on the subunit and match the reported open probabilities for the receptor.

Evolution of sparsity and modularity in a model of protein allostery

The sequence of a protein is not only constrained by its physical and biochemical properties under current selection, but also by features of its past evolutionary history. Understanding the extent and the form that these evolutionary constraints may take is important to interpret the information in protein sequences. To study this problem, we introduce a simple but physical model of protein evolution where selection targets allostery, the functional coupling of distal sites on protein surfaces. This model shows how the geometrical organization of couplings between amino acids within a protein structure can depend crucially on its evolutionary history. In particular, two scenarios are found to generate a spatial concentration of functional constraints: high mutation rates and fluctuating selective pressures. This second scenario offers a plausible explanation for the high tolerance of natural proteins to mutations and for the spatial organization of their least tolerant amino acids, as revealed by sequence analysis and mutagenesis experiments. It also implies a faculty to adapt to new selective pressures that is consistent with observations. The model illustrates how several independent functional modules may emerge within the same protein structure, depending on the nature of past environmental fluctuations. Our model thus relates the evolutionary history of proteins to the geometry of their functional constraints, with implications for decoding and engineering protein sequences.

Allosteric activation of Bordetella pertussis adenylyl cyclase by calmodulin: molecular dynamics and mutagenesis studies

Adenylyl cyclase (AC) toxin is an essential toxin that allows Bordetella pertussis to invade eukaryotic cells, where it is activated after binding to calmodulin (CaM). Based on the crystal structure of the AC catalytic domain in complex with the C-terminal half of CaM (C-CaM), our previous molecular dynamics simulations (Selwa, E., Laine, E., and Malliavin, T. (2012) Differential role of calmodulin and calcium ions in the stabilization of the catalytic domain of adenyl cyclase CyaA from Bordetella pertussis. Proteins 80, 1028–1040) suggested that three residues (i.e. Arg(338), Asn(347), and Asp(360)) might be important for stabilizing the AC/CaM interaction. These residues belong to a loop-helix-loop motif at the C-terminal end of AC, which is located at the interface between CaM and the AC catalytic loop. In the present study, we conducted the in silico and in vitro characterization of three AC variants, where one (Asn(347); ACm1A), two (Arg(338) and Asp(360); ACm2A), or three residues (Arg(338), Asn(347), and Asp(360); ACm3A) were substituted with Ala. Biochemical studies showed that the affinities of ACm1A and ACm2A for CaM were not affected significantly, whereas that of ACm3A was reduced dramatically. To understand the effects of these modifications, molecular dynamics simulations were performed based on the modified proteins. The molecular dynamics trajectories recorded for the ACm3AC-CaM complex showed that the calcium-binding loops of C-CaM exhibited large fluctuations, which could be related to the weakened interaction between ACm3A and its activator. Overall, our results suggest that the loop-helix-loop motif at the C-terminal end of AC is crucial during CaM binding for stabilizing the AC catalytic loop in an active configuration.

Hotspot mutations in KIT receptor differentially modulate its allosterically coupled conformational dynamics: impact on activation and drug sensitivity

Receptor tyrosine kinase KIT controls many signal transduction pathways and represents a typical allosterically regulated protein. The mutation-induced deregulation of KIT activity impairs cellular physiological functions and causes serious human diseases. The impact of hotspots mutations (D816H/Y/N/V and V560G/D) localized in crucial regulatory segments, the juxtamembrane region (JMR) and the activation (A-) loop, on KIT internal dynamics was systematically studied by molecular dynamics simulations. The mutational outcomes predicted in silico were correlated with in vitro and in vivo activation rates and drug sensitivities of KIT mutants. The allosteric regulation of KIT in the native and mutated forms is described in terms of communication between the two remote segments, JMR and A-loop. A strong correlation between the communication profile and the structural and dynamical features of KIT in the native and mutated forms was established. Our results provide new insight on the determinants of receptor KIT constitutive activation by mutations and resistance of KIT mutants to inhibitors. Depiction of an intra-molecular component of the communication network constitutes a first step towards an integrated description of vast communication pathways established by KIT in physiopathological contexts.

Receptor tyrosine kinase KIT plays a crucial role in the regulation of cell signaling. This allosterically controlled activity may be affected by gain-of-function mutations that promote the development of several cancers. Identification of the molecular basis of KIT constitutive activation and allosteric regulation has inspired computational study of KIT hotspot mutations. In the present contribution, we investigated the mutation-induced effects on KIT conformational dynamics and intra-protein communication conditionally on the mutation location and the nature of the substituting amino acid. Our data elucidate that all studied mutations stabilize an inactive non-autoinhibited state of KIT over the inactive auto-inhibited state prevalent for the native protein. This shift in the protein conformational landscape promotes KIT constitutive activation. Our in silico analysis established correlations between the structural and dynamical effects induced by oncogenic mutations and the mutants auto-activation rates and drug sensitivities measured in vitro and in vivo. Particularly, the A-loop mutations stabilize the drug-resistant forms, while the JMR mutations may facilitate inhibitors binding to the active site. Cross-correlations established between local and long-range structural and dynamical effects demonstrate the allosteric character of the gain-of-function mutations mode of action.

Allosteric inhibition of human immunodeficiency virus integrase: late block during viral replication and abnormal multimerization involving specific protein domains

HIV-1 replication in the presence of antiviral agents results in evolution of drug-resistant variants, motivating the search for additional drug classes. Here we report studies of GSK1264, which was identified as a compound that disrupts the interaction between HIV-1 integrase (IN) and the cellular factor lens epithelium-derived growth factor (LEDGF)/p75. GSK1264 displayed potent antiviral activity and was found to bind at the site occupied by LEDGF/p75 on IN by x-ray crystallography. Assays of HIV replication in the presence of GSK1264 showed only modest inhibition of the early infection steps and little effect on integration targeting, which is guided by the LEDGF/p75-IN interaction. In contrast, inhibition of late replication steps was more potent. Particle production was normal, but particles showed reduced infectivity. GSK1264 promoted aggregation of IN and preformed LEDGF/p75-IN complexes, suggesting a mechanism of inhibition. LEDGF/p75 was not displaced from IN during aggregation, indicating trapping of LEDGF/p75 in aggregates. Aggregation assays with truncated IN variants revealed that a construct with catalytic and C-terminal domains of IN only formed an open polymer associated with efficient drug-induced aggregation. These data suggest that the allosteric inhibitors of IN are promising antiviral agents and provide new information on their mechanism of action.

Structure of fully liganded Hb ζ2β2s trapped in a tense conformation

A variant Hb ζ2β2(s) that is formed from sickle hemoglobin (Hb S; α2β2(s)) by exchanging adult α-globin with embryonic ζ-globin subunits shows promise as a therapeutic agent for sickle-cell disease (SCD). Hb ζ2β2(s) inhibits the polymerization of deoxygenated Hb S in vitro and reverses characteristic features of SCD in vivo in mouse models of the disorder. When compared with either Hb S or with normal human adult Hb A (α2β2), Hb ζ2β2(s) exhibits atypical properties that include a high oxygen affinity, reduced cooperativity, a weak Bohr effect and blunted 2,3-diphosphoglycerate allostery. Here, the 1.95 Å resolution crystal structure of human Hb ζ2β2(s) that was expressed in complex transgenic knockout mice and purified from their erythrocytes is presented. When fully liganded with carbon monoxide, Hb ζ2β2(s) displays a central water cavity, a ζ1-β(s)2 (or ζ2-β(s)1) interface, intersubunit salt-bridge/hydrogen-bond interactions, C-terminal βHis146 salt-bridge interactions, and a β-cleft, that are highly unusual for a relaxed hemoglobin structure and are more typical of a tense conformation. These quaternary tense-like features contrast with the tertiary relaxed-like conformations of the ζ1β(s)1 dimer and the CD and FG corners, as well as the overall structures of the heme cavities. This crystallographic study provides insights into the altered oxygen-transport properties of Hb ζ2β2(s) and, moreover, decouples tertiary- and quaternary-structural events that are critical to Hb ligand binding and allosteric function.

Endogenous positive allosteric modulation of GABA(A) receptors by diazepam binding inhibitor

Benzodiazepines (BZs) allosterically modulate γ-aminobutyric acid type-A receptors (GABAARs) to increase inhibitory synaptic strength. Diazepam binding inhibitor (DBI) protein is a BZ site ligand expressed endogenously in the brain, but functional evidence for BZ-mimicking positive modulatory actions has been elusive. We demonstrate an endogenous potentiation of GABAergic synaptic transmission and responses to GABA uncaging in the thalamic reticular nucleus (nRT) that is absent in both nm1054 mice, in which the Dbi gene is deleted, and mice in which BZ binding to α3 subunit-containing GABAARs is disrupted. Viral transduction of DBI into nRT is sufficient to rescue the endogenous potentiation of GABAergic transmission in nm1054 mice. Both mutations enhance thalamocortical spike-and-wave discharges characteristic of absence epilepsy. Together, these results indicate that DBI mediates endogenous nucleus-specific BZ-mimicking ("endozepine") roles to modulate nRT function and suppress thalamocortical oscillations. Enhanced DBI signaling might serve as a therapy for epilepsy and other neurological disorders.

Broadly neutralizing antibody PGT121 allosterically modulates CD4 binding via recognition of the HIV-1 gp120 V3 base and multiple surrounding glycans

New broad and potent neutralizing HIV-1 antibodies have recently been described that are largely dependent on the gp120 N332 glycan for Env recognition. Members of the PGT121 family of antibodies, isolated from an African donor, neutralize ∼70% of circulating isolates with a median IC₅₀ less than 0.05 µg ml⁻¹. Here, we show that three family members, PGT121, PGT122 and PGT123, have very similar crystal structures. A long 24-residue HCDR3 divides the antibody binding site into two functional surfaces, consisting of an open face, formed by the heavy chain CDRs, and an elongated face, formed by LCDR1, LCDR3 and the tip of the HCDR3. Alanine scanning mutagenesis of the antibody paratope reveals a crucial role in neutralization for residues on the elongated face, whereas the open face, which accommodates a complex biantennary glycan in the PGT121 structure, appears to play a more secondary role. Negative-stain EM reconstructions of an engineered recombinant Env gp140 trimer (SOSIP.664) reveal that PGT122 interacts with the gp120 outer domain at a more vertical angle with respect to the top surface of the spike than the previously characterized antibody PGT128, which is also dependent on the N332 glycan. We then used ITC and FACS to demonstrate that the PGT121 antibodies inhibit CD4 binding to gp120 despite the epitope being distal from the CD4 binding site. Together, these structural, functional and biophysical results suggest that the PGT121 antibodies may interfere with Env receptor engagement by an allosteric mechanism in which key structural elements, such as the V3 base, the N332 oligomannose glycan and surrounding glycans, including a putative V1/V2 complex biantennary glycan, are conformationally constrained.

An estimated 33 million adults and children currently live with the human immunodeficiency virus type 1 (HIV-1), which represents a global prevalence of 0.8%. In the absence of a cure, the development of a protective vaccine is the long sought-after goal in containment of the pandemic. HIV-1 Env is the sole viral surface glycoprotein and mediates viral engagement and entry into host cells, which constitutes the first step of the virus life cycle. Recently, a plethora of exciting new antibodies have been discovered that interact with HIV-1 Env and inhibit infection of target cells (i.e. neutralize the virus). Here, we structurally characterize the interaction of a recombinant HIV-1 Env with one class of such antibodies, namely antibodies of the PGT121 family. These studies have uncovered a novel mode of HIV-1 Env recognition. By interacting with key structural elements of HIV-1 Env near the apex at its membrane-distal end, these antibodies can interfere with binding to CD4, the receptor on T cells that is required for HIV-1 infection. These observations further delineate a glycan-dependent site of vulnerability on HIV-1 Env that can be used in vaccine design efforts.