Controlling Allosteric Networks in Proteins

Mapping allosteric linkage to channel gating by extracellular domains in the human epithelial sodium channel

The epithelial sodium channel (ENaC) mediates sodium absorption in lung, kidney, and colon epithelia. Channels in the ENaC/degenerin family possess an extracellular region that senses physicochemical changes in the extracellular milieu and allosterically regulates the channel opening. Proteolytic cleavage activates the ENaC opening, by the removal of specific segments in the finger domains of the α- and γ ENaC-subunits. Cleavage causes perturbations in the extracellular region that propagate to the channel gate. However, it is not known how the channel structure mediates the propagation of activation signals through the extracellular sensing domains. Here, to identify the structure-function determinants that mediate allosteric ENaC activation, we performed MD simulations, thiol modification of residues substituted by cysteine, and voltage-clamp electrophysiology recordings. Our simulations of an ENaC heterotetramer, α₁βα₂γ, in the proteolytically cleaved and uncleaved states revealed structural pathways in the α-subunit that are responsible for ENaC proteolytic activation. To validate these findings, we performed site-directed mutagenesis to introduce cysteine substitutions in the extracellular domains of the α-, β-, and γ ENaC-subunits. Insertion of a cysteine at the α-subunit Glu⁵⁵⁷ site, predicted to stabilize a closed state of ENaC, inhibited ENaC basal activity and retarded the kinetics of proteolytic activation by 2-fold. Our results suggest that the lower palm domain of αENaC is essential for ENaC activation. In conclusion, our integrated computational and experimental approach suggests key structure-function determinants for ENaC proteolytic activation and points toward a mechanistic model for the allosteric communication in the extracellular domains of the ENaC/degenerin family channels.

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.

Effects of Distal Mutations on the Structure, Dynamics and Catalysis of Human Monoacylglycerol Lipase

An understanding of how conformational dynamics modulates function and catalysis of human monoacylglycerol lipase (hMGL), an important pharmaceutical target, can facilitate the development of novel ligands with potential therapeutic value. Here, we report the discovery and characterization of an allosteric, regulatory hMGL site comprised of residues Trp-289 and Leu-232 that reside over 18 Å away from the catalytic triad. These residues were identified as critical mediators of long-range communication and as important contributors to the integrity of the hMGL structure. Nonconservative replacements of Trp-289 or Leu-232 triggered concerted motions of structurally distinct regions with a significant conformational shift toward inactive states and dramatic loss in catalytic efficiency of the enzyme. Using a multimethod approach, we show that the dynamically relevant Trp-289 and Leu-232 residues serve as communication hubs within an allosteric protein network that controls signal propagation to the active site, and thus, regulates active-inactive interconversion of hMGL. Our findings provide new insights into the mechanism of allosteric regulation of lipase activity, in general, and may provide alternative drug design possibilities.

In silico identification of rescue sites by double force scanning

Motivation

A deleterious amino acid change in a protein can be compensated by a second-site rescue mutation. These compensatory mechanisms can be mimicked by drugs. In particular, the location of rescue mutations can be used to identify protein regions that can be targeted by small molecules to reactivate a damaged mutant.

Results

We present the first general computational method to detect rescue sites. By mimicking the effect of mutations through the application of forces, the double force scanning (DFS) method identifies the second-site residues that make the protein structure most resilient to the effect of pathogenic mutations. We tested DFS predictions against two datasets containing experimentally validated and putative evolutionary-related rescue sites. A remarkably good agreement was found between predictions and experimental data. Indeed, almost half of the rescue sites in p53 was correctly predicted by DFS, with 65% of remaining sites in contact with DFS predictions. Similar results were found for other proteins in the evolutionary dataset.

Availability and implementation

The DFS code is available under GPL at https://fornililab.github.io/dfs/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Allosteric Modulation of Intact γ-Secretase Structural Dynamics

As a protease complex involved in the cleavage of amyloid precursor proteins that lead to the formation of amyloid β fibrils implicated in Alzheimer's disease, γ-secretase is an important target for developing therapeutics against Alzheimer's disease. γ-secretase is composed of four subunits: nicastrin (NCT) in the extracellular (EC) domain, presenilin-1 (PS1), anterior pharynx defective 1, and presenilin enhancer 2 in the transmembrane (TM) domain. NCT and PS1 play important roles in binding amyloid β precursor proteins and modulating PS1 catalytic activity. Yet, the molecular mechanisms of coupling between substrate/modulator binding and catalytic activity remain to be elucidated. Recent determination of intact human γ-secretase cryo-electron microscopy structure has opened the way for a detailed investigation of the structural dynamics of this complex. Our analysis, based on a membrane-coupled anisotropic network model, reveals two types of NCT motions, bending and twisting, with respect to PS1. These underlie the fluctuations between the "open" and "closed" states of the lid-like NCT with respect to a hydrophilic loop 1 (HL1) on PS1, thus allowing or blocking access of the substrate peptide (EC portion) to HL1 and to the neighboring helix TM2. In addition to this alternating access mechanism, fluctuations in the volume of the PS1 central cavity facilitate the exposure of the catalytic site for substrate cleavage. Druggability simulations show that γ-secretase presents several hot spots for either orthosteric or allosteric inhibition of catalytic activity, consistent with experimental data. In particular, a hinge region at the interface between the EC and TM domains, near the interlobe groove of NCT, emerges as an allo-targeting site that would impact the coupling between HL1/TM2 and the catalytic pocket, opening, to our knowledge, new avenues for structure-based design of novel allosteric modulators of γ-secretase protease activity.

Single Binding Pockets Versus Allosteric Binding

An orthosteric site is commonly viewed as the primary, functionally binding pocket on a receptor. Signal molecules, endogenous agonists, and substrates are recognized by and bind to the orthosteric site of a specific target, resulting in a biological effect. A malfunctioning active site on a crucial receptor has been confirmed as the culprit that causes many metabolic disturbances, neurologic disorders, and genetic diseases. A competitive inhibitor that has a stronger binding affinity can outcompete an orthosteric ligand. An allosteric site, which is nonoverlapping and topographically distinct from the active pocket, can emerge as a potential regulatory site on the protein surface. An allosteric modulator interacts with a specific binding site, affecting the atoms of nearby residues, thus eliciting a series of conformational changes in the residues at the active site through propagation pathways. Allosteric regulation can potentiate or inhibit function instead of blocking it, and this is a promising strategy for drug design. In this chapter, we describe the tools and protocols for allosteric site analysis and allosteric ligand design.

Ensemble- and Rigidity Theory-Based Perturbation Approach To Analyze Dynamic Allostery

Allostery describes the functional coupling between sites in biomolecules. Recently, the role of changes in protein dynamics for allosteric communication has been highlighted. A quantitative and predictive description of allostery is fundamental for understanding biological processes. Here, we integrate an ensemble-based perturbation approach with the analysis of biomolecular rigidity and flexibility to construct a model of dynamic allostery. Our model, by definition, excludes the possibility of conformational changes, evaluates static, not dynamic, properties of molecular systems, and describes allosteric effects due to ligand binding in terms of a novel free-energy measure. We validated our model on three distinct biomolecular systems: eglin c, protein tyrosine phosphatase 1B, and the lymphocyte function-associated antigen 1 domain. In all cases, it successfully identified key residues for signal transmission in very good agreement with the experiment. It correctly and quantitatively discriminated between positively or negatively cooperative effects for one of the systems. Our model should be a promising tool for the rational discovery of novel allosteric drugs.

Decomposing Dynamical Couplings in Mutated scFv Antibody Fragments into Stabilizing and Destabilizing Effects

Conformational fluctuations within scFv antibodies are characterized by a novel perturbation-response decomposition of molecular dynamics trajectories. Both perturbation and response profiles are stratified into stabilizing and destabilizing conditions. The linker between the VH and VL domains exhibits the dominant dynamical response by being coupled to nearly the entire protein, responding to both stabilizing and destabilizing perturbations. Perturbations within complementarity-determining regions (CDR) induce rich behavior in dynamic response. Among many effects, stabilizing any CDR loop in the VH domain triggers a destabilizing response in all CDR loops in the VL domain and vice versa. Destabilizing residues within the VL domain are likely to stabilize all CDR loops in the VH domain, and, when these residues are not buried, the CDR loops in the VL domain are also likely to be stabilized. These effects, described by shifts in normal mode characteristics, initiate a propensity for dynamic allostery with possible functional implications in bispecific antibodies.

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BACKGROUND

Accumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.

METHODS/RESULTS

We employ two computationally efficient approaches - graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.

CONCLUSIONS

Contact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.

GENERAL SIGNIFICANCE

This study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.

Time-resolved observation of protein allosteric communication

Allostery represents a fundamental mechanism of biological regulation that is mediated via long-range communication between distant protein sites. Although little is known about the underlying dynamical process, recent time-resolved infrared spectroscopy experiments on a photoswitchable PDZ domain (PDZ2S) have indicated that the allosteric transition occurs on multiple timescales. Here, using extensive nonequilibrium molecular dynamics simulations, a time-dependent picture of the allosteric communication in PDZ2S is developed. The simulations reveal that allostery amounts to the propagation of structural and dynamical changes that are genuinely nonlinear and can occur in a nonlocal fashion. A dynamic network model is constructed that illustrates the hierarchy and exceeding structural heterogeneity of the process. In compelling agreement with experiment, three physically distinct phases of the time evolution are identified, describing elastic response ([Formula: see text] ns), inelastic reorganization ([Formula: see text] ns), and structural relaxation ([Formula: see text]s). Issues such as the similarity to downhill folding as well as the interpretation of allosteric pathways are discussed.

Leveraging Reciprocity to Identify and Characterize Unknown Allosteric Sites in Protein Tyrosine Phosphatases

Drug-like molecules targeting allosteric sites in proteins are of great therapeutic interest; however, identification of potential sites is not trivial. A straightforward approach to identify hidden allosteric sites is demonstrated in protein tyrosine phosphatases (PTP) by creation of single alanine mutations in the catalytic acid loop of PTP1B and VHR. This approach relies on the reciprocal interactions between an allosteric site and its coupled orthosteric site. The resulting NMR chemical shift perturbations (CSPs) of each mutant reveal clusters of distal residues affected by acid loop mutation. In PTP1B and VHR, two new allosteric clusters were identified in each enzyme. Mutations in these allosteric clusters altered phosphatase activity with changes in k_cat/K_M ranging from 30% to nearly 100-fold. This work outlines a simple method for identification of new allosteric sites in PTP, and given the basis of this method in thermodynamics, it is expected to be generally useful in other systems.

Identifying Protein Allosteric Transitions for Drug Discovery with 1D NMR

Allosteric drugs present many advantages over orthosteric drugs and are therefore an attractive approach in drug discovery, despite being highly challenging. First, the binding of ligands in protein allosteric pockets do not ensure an allosteric effect, and second, allosteric ligands can possess diverse modes of pharmacology even within a compound family. Herein we report a new method to: 1) detect allosteric communication between protein binding sites, and 2) compare the effect of allosteric ligands on the allosteric transitions of the protein target. The method, illustrated with glycogen phosphorylase, consists of comparing 1D saturation transfer difference (STD) NMR spectra of a molecular spy (here fragments) in the absence and presence of allosteric ligands. The modification of the STD NMR spectrum of the fragment indicates whether the protein dynamics/conformations have been changed in the presence of the allosteric modulator, thereby highlighting allosteric coupling between the binding pocket of the reference compound (in this case the fragment) and the allosteric pocket.

Designing allostery-inspired response in mechanical networks

Recent advances in designing metamaterials have demonstrated that global mechanical properties of disordered spring networks can be tuned by selectively modifying only a small subset of bonds. Here, using a computationally efficient approach, we extend this idea to tune more general properties of networks. With nearly complete success, we are able to produce a strain between any two target nodes in a network in response to an applied source strain on any other pair of nodes by removing only ∼1% of the bonds. We are also able to control multiple pairs of target nodes, each with a different individual response, from a single source, and to tune multiple independent source/target responses simultaneously into a network. We have fabricated physical networks in macroscopic 2D and 3D systems that exhibit these responses. This work is inspired by the long-range coupled conformational changes that constitute allosteric function in proteins. The fact that allostery is a common means for regulation in biological molecules suggests that it is a relatively easy property to develop through evolution. In analogy, our results show that long-range coupled mechanical responses are similarly easy to achieve in disordered networks.

Allosteric Transmission along a Loosely Structured Backbone Allows a Cardiac Troponin C Mutant to Function with Only One Ca2+ Ion

Hypertrophic cardiomyopathy (HCM) is one of the most common cardiomyopathies and a major cause of sudden death in young athletes. The Ca²⁺ sensor of the sarcomere, cardiac troponin C (cTnC), plays an important role in regulating muscle contraction. Although several cardiomyopathy-causing mutations have been identified in cTnC, the limited information about their structural defects has been mapped to the HCM phenotype. Here, we used high-resolution electron-spray ionization mass spectrometry (ESI-MS), Carr-Purcell-Meiboom-Gill relaxation dispersion (CPMG-RD), and affinity measurements of cTnC for the thin filament in reconstituted papillary muscles to provide evidence of an allosteric mechanism in mutant cTnC that may play a role to the HCM phenotype. We showed that the D145E mutation leads to altered dynamics on a μs-ms time scale and deactivates both of the divalent cation-binding sites of the cTnC C-domain. CPMG-RD captured a low populated protein-folding conformation triggered by the Glu-145 replacement of Asp. Paradoxically, although D145E C-domain was unable to bind Ca²⁺, these changes along its backbone allowed it to attach more firmly to thin filaments than the wild-type isoform, providing evidence for an allosteric response of the Ca²⁺-binding site II in the N-domain. Our findings explain how the effects of an HCM mutation in the C-domain reflect up into the N-domain to cause an increase of Ca²⁺ affinity in site II, thus opening up new insights into the HCM phenotype.

DNA-binding protects p53 from interactions with cofactors involved in transcription-independent functions

Binding-induced conformational changes of a protein at regions distant from the binding site may play crucial roles in protein function and regulation. The p53 tumour suppressor is an example of such an allosterically regulated protein. Little is known, however, about how DNA binding can affect distal sites for transcription factors. Furthermore, the molecular details of how a local perturbation is transmitted through a protein structure are generally elusive and occur on timescales hard to explore by simulations. Thus, we employed state-of-the-art enhanced sampling atomistic simulations to unveil DNA-induced effects on p53 structure and dynamics that modulate the recruitment of cofactors and the impact of phosphorylation at Ser215. We show that DNA interaction promotes a conformational change in a region 3 nm away from the DNA binding site. Specifically, binding to DNA increases the population of an occluded minor state at this distal site by more than 4-fold, whereas phosphorylation traps the protein in its major state. In the minor conformation, the interface of p53 that binds biological partners related to p53 transcription-independent functions is not accessible. Significantly, our study reveals a mechanism of DNA-mediated protection of p53 from interactions with partners involved in the p53 transcription-independent signalling. This also suggests that conformational dynamics is tightly related to p53 signalling.

Computational approaches to investigating allostery

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

Structural basis for KCNE3 modulation of potassium recycling in epithelia

The single-span membrane protein KCNE3 modulates a variety of voltage-gated ion channels in diverse biological contexts. In epithelial cells, KCNE3 regulates the function of the KCNQ1 potassium ion (K(+)) channel to enable K(+) recycling coupled to transepithelial chloride ion (Cl(-)) secretion, a physiologically critical cellular transport process in various organs and whose malfunction causes diseases, such as cystic fibrosis (CF), cholera, and pulmonary edema. Structural, computational, biochemical, and electrophysiological studies lead to an atomically explicit integrative structural model of the KCNE3-KCNQ1 complex that explains how KCNE3 induces the constitutive activation of KCNQ1 channel activity, a crucial component in K(+) recycling. Central to this mechanism are direct interactions of KCNE3 residues at both ends of its transmembrane domain with residues on the intra- and extracellular ends of the KCNQ1 voltage-sensing domain S4 helix. These interactions appear to stabilize the activated "up" state configuration of S4, a prerequisite for full opening of the KCNQ1 channel gate. In addition, the integrative structural model was used to guide electrophysiological studies that illuminate the molecular basis for how estrogen exacerbates CF lung disease in female patients, a phenomenon known as the "CF gender gap."

A Panel of TrpB Biocatalysts Derived from Tryptophan Synthase through the Transfer of Mutations that Mimic Allosteric Activation

Naturally occurring enzyme homologues often display highly divergent activity with non-natural substrates. Exploiting this diversity with enzymes engineered for new or altered function, however, is laborious because the engineering must be replicated for each homologue. A small set of mutations of the tryptophan synthase β-subunit (TrpB) from Pyrococcus furiosus, which mimics the activation afforded by binding of the α-subunit, was demonstrated to have a similar activating effect in different TrpB homologues with as little as 57 % sequence identity. Kinetic and spectroscopic analyses indicate that the mutations function through the same mechanism: mimicry of α-subunit binding. From these enzymes, we identified a new TrpB catalyst that displays a remarkably broad activity profile in the synthesis of 5-substituted tryptophans. This demonstrates that allosteric activation can be recapitulated throughout a protein family to explore natural sequence diversity for desirable biocatalytic transformations.

Allosteric Inhibition of Factor XIIIa. Non-Saccharide Glycosaminoglycan Mimetics, but Not Glycosaminoglycans, Exhibit Promising Inhibition Profile

Factor XIIIa (FXIIIa) is a transglutaminase that catalyzes the last step in the coagulation process. Orthostery is the only approach that has been exploited to design FXIIIa inhibitors. Yet, allosteric inhibition of FXIIIa is a paradigm that may offer a key advantage of controlled inhibition over orthosteric inhibition. Such an approach is likely to lead to novel FXIIIa inhibitors that do not carry bleeding risks. We reasoned that targeting a collection of basic amino acid residues distant from FXIIIa’s active site by using sulfated glycosaminoglycans (GAGs) or non-saccharide GAG mimetics (NSGMs) would lead to the discovery of the first allosteric FXIIIa inhibitors. We tested a library of 22 variably sulfated GAGs and NSGMs against human FXIIIa to discover promising hits. Interestingly, although some GAGs bound to FXIIIa better than NSGMs, no GAG displayed any inhibition. An undecasulfated quercetin analog was found to inhibit FXIIIa with reasonable potency (efficacy of 98%). Michaelis-Menten kinetic studies revealed an allosteric mechanism of inhibition. Fluorescence studies confirmed close correspondence between binding affinity and inhibition potency, as expected for an allosteric process. The inhibitor was reversible and at least 9-fold- and 26-fold selective over two GAG-binding proteins factor Xa (efficacy of 71%) and thrombin, respectively, and at least 27-fold selective over a cysteine protease papain. The inhibitor also inhibited the FXIIIa-mediated polymerization of fibrin in vitro. Overall, our work presents the proof-of-principle that FXIIIa can be allosterically modulated by sulfated non-saccharide agents much smaller than GAGs, which should enable the design of selective and safe anticoagulants.

The dynamics of the β-propeller domain in Kelch protein KLHL40 changes upon nemaline myopathy-associated mutation

The nemaline myopathy-associated E528K mutation in the KLHL40 alters the communication between the Kelch propeller blades.