Allosteric Network Research Articles

Using Optical Tweezers to Dissect Allosteric Communication Networks in Protein Kinases.

Mutations in protein kinases are often associated with the development of cancer, and application of mutant-specific inhibitors as therapeutic measures have shown a remarkable improvement in prolonging patient survival. However, it has also been observed that tumors bearing certain mutation types are more resistant to current approved drugs. Importantly, many resistant mutations are located in regions outside substrate or inhibitor binding sites, indicating allosteric effects. Understanding how mutations trigger effects over a distant site of the protein requires a deeper investigation ofthe molecular origin of allosteric regulation networks in kinases. In this chapter, we show the application ofsingle-molecule optical tweezers to selectively manipulate specific regions of proteins to trace allosteric signals, thereby allowing theelucidation ofallosteric communication networks. We illustrate this approach using as model system the regulatory subunit of protein kinase A. This single-molecule optical tweezers approach, however, can be readily applicable to studyother kinases, and can be further expanded to screen potential allosteric drugs for future therapeutics.

Equilibria between conformational states of the Ras oncogene protein revealed by high pressure crystallography.

In this work, we experimentally investigate the allosteric transitions between conformational states on the Ras oncogene protein using high pressure crystallography. Ras protein is a small GTPase involved in central regulatory processes occurring in multiple conformational states. Ras acts as a molecular switch between active GTP-bound, and inactive GDP-bound states, controlling essential signal transduction pathways. An allosteric network of interactions between the effector binding regions and the membrane interacting regions is involved in Ras cycling. The conformational states which coexist simultaneously in solution possess higher Gibbs free energy than the ground state. Equilibria between these states can be shifted by applying pressure favouring conformations with lower partial molar volume, and has been previously analyzed by high-pressure NMR spectroscopy. High-pressure macromolecular crystallography (HPMX) is a powerful tool perfectly complementary to high-pressure NMR, allowing characterization at the molecular level with a high resolution the different allosteric states involved in the Ras cycling. We observe a transition above 300 MPa in the crystal leading to more stable conformers. Thus, we compare the crystallographic structures of Ras(wt)·Mg2+·GppNHp and Ras(D33K)·Mg2+·GppNHp at various high hydrostatic pressures. This gives insight into per-residue descriptions of the structural plasticity involved in allosteric equilibria between conformers. We have mapped out at atomic resolution the different segments of Ras protein which remain in the ground-state conformation or undergo structural changes, adopting excited-energy conformations corresponding to transient intermediate states. Such in crystallo phase transitions induced by pressure open the possibility to finely explore the structural determinants related to switching between Ras allosteric sub-states without any mutations nor exogenous partners.

Client binding shifts the populations of dynamic Hsp90 conformations through an allosteric network.

Hsp90 is a molecular chaperone that interacts with a specific set of client proteins and assists their folding. The underlying molecular mechanisms, involving dynamic transitions between open and closed conformations, are still enigmatic. Combining nuclear magnetic resonance, small-angle x-ray scattering, and biochemical experiments, we have identified a key intermediate state of Hsp90 induced by adenosine triphosphate (ATP) binding, in which rotation of the Hsp90 N-terminal domain (NTD) yields a domain arrangement poised for closing. This ATP-stabilized NTD rotation is allosterically communicated across the full Hsp90 dimer, affecting distant client sites. By analyzing the interactions of four distinct clients, i.e., steroid hormone receptors (glucocorticoid receptor and mineralocorticoid receptor), p53, and Tau, we show that client-specific interactions with Hsp90 select and enhance the NTD-rotated state and promote closing of the full-length Hsp90 dimer. The p23 co-chaperone shifts the population of Hsp90 toward the closed state, thereby enhancing client interaction and processing.

Mutation-Induced Long-Range Allosteric Interactions in the Spike Protein Determine the Infectivity of SARS-CoV-2 Emerging Variants.

The emergence of a variety of highly transmissible SARS-CoV-2 variants, the causative agent of COVID-19, with multiple spike mutations poses serious challenges in overcoming the ongoing deadly pandemic. It is, therefore, essential to understand how these variants gain enhanced ability to evade immune responses with a higher rate of spreading infection. To address this question, here we have individually assessed the effects of SARS-CoV-2 variant-specific spike (S) protein receptor-binding domain (RBD) mutations E484K, K417N, L452Q, L452R, N501Y, and T478K that characterize and differentiate several emerging variants. Despite the hundreds of apparently neutral mutations observed in the domains other than the RBD, we have shown that each RBD mutation site is differentially engaged in an interdomain allosteric network involving mutation sites from a distant domain, affecting interactions with the human receptor angiotensin-converting enzyme-2 (ACE2). This allosteric network couples the residues of the N-terminal domain (NTD) and the RBD, which are modulated by the RBD-specific mutations and are capable of propagating mutation-induced perturbations between these domains through a combination of structural changes and effector-dependent modulations of dynamics. One key feature of this network is the inclusion of compensatory mutations segregated into three characteristically different clusters, where each cluster residue site is allosterically coupled with specific RBD mutation sites. Notably, each RBD mutation acted like a positive allosteric modulator; nevertheless, K417N was shown to have the largest effects among all of the mutations on the allostery and thereby holds the highest binding affinity with ACE2. This result will be useful for designing the targeted control measure and therapeutic efforts aiming at allosteric modulators.

Molecular reconstruction of recurrent evolutionary switching in olfactory receptor specificity.

Olfactory receptor repertoires exhibit remarkable functional diversity, but how these proteins have evolved is poorly understood. Through analysis of extant and ancestrally reconstructed drosophilid olfactory receptors from the Ionotropic receptor (Ir) family, we investigated evolution of two organic acid-sensing receptors, Ir75a and Ir75b. Despite their low amino acid identity, we identify a common 'hotspot' in their ligand-binding pocket that has a major effect on changing the specificity of both Irs, as well as at least two distinct functional transitions in Ir75a during evolution. Moreover, we show that odor specificity is refined by changes in additional, receptor-specific sites, including those outside the ligand-binding pocket. Our work reveals how a core, common determinant of ligand-tuning acts within epistatic and allosteric networks of substitutions to lead to functional evolution of olfactory receptors.

Perturbing dimer interactions and allosteric communication modulates the immunosuppressive activity of human galectin-7

The design of allosteric modulators to control protein function is a key objective in drug discovery programs. Altering functionally essential allosteric residue networks provides unique protein family subtype specificity, minimizes unwanted off-target effects, and helps avert resistance acquisition typically plaguing drugs that target orthosteric sites. In this work, we used protein engineering and dimer interface mutations to positively and negatively modulate the immunosuppressive activity of the proapoptotic human galectin-7 (GAL-7). Using the PoPMuSiC and BeAtMuSiC algorithms, mutational sites and residue identity were computationally probed and predicted to either alter or stabilize the GAL-7 dimer interface. By designing a covalent disulfide bridge between protomers to control homodimer strength and stability, we demonstrate the importance of dimer interface perturbations on the allosteric network bridging the two opposite glycan-binding sites on GAL-7, resulting in control of induced apoptosis in Jurkat T cells. Molecular investigation of G16X GAL-7 variants using X-ray crystallography, biophysical, and computational characterization illuminates residues involved in dimer stability and allosteric communication, along with discrete long-range dynamic behaviors involving loops 1, 3, and 5. We show that perturbing the protein–protein interface between GAL-7 protomers can modulate its biological function, even when the overall structure and ligand-binding affinity remains unaltered. This study highlights new avenues for the design of galectin-specific modulators influencing both glycan-dependent and glycan-independent interactions.

Allosteric Control of Structural Mimicry and Mutational Escape in the SARS-CoV-2 Spike Protein Complexes with the ACE2 Decoys and Miniprotein Inhibitors: A Network-Based Approach for Mutational Profiling of Binding and Signaling.

We developed a computational framework for comprehensive and rapid mutational scanning of binding energetics and residue interaction networks in the SARS-CoV-2 spike protein complexes. Using this approach, we integrated atomistic simulations and conformational landscaping of the SARS-CoV-2 spike protein complexes with ensemble-based mutational screening and network modeling to characterize mechanisms of structure-functional mimicry and resilience toward mutational escape by the ACE2 protein decoy and de novo designed miniprotein inhibitors. A detailed analysis of structural plasticity of the SARS-CoV-2 spike proteins obtained from atomistic simulations of conformational landscapes and sequence-based profiling of the disorder propensities revealed the intrinsically flexible regions that harbor key functional sites targeted by circulating variants. The conservation of collective dynamics in the SARS-CoV-2 spike protein complexes showed that mutational escape positions are important for modulation of functional motions and that mutational changes in these sites can alter allosteric interaction networks. Through mutational profiling of binding and allosteric propensities in the SARS-CoV-2 spike protein complexes, we identified the key binding and regulatory hotspots that collectively determine functional response and resilience of miniproteins to mutational variants. The results suggest that binding affinities and allosteric signatures of the SARS-CoV-2 complexes can be determined by dynamic crosstalk between structurally stable regulatory centers and conformationally adaptable allosteric hotspots that collectively control the resilience toward mutational escape. This may underlie a mechanism in which moderate perturbations in the mutational escape positions can induce global allosteric changes and alter functional protein response by modulating signaling in the residue interaction networks.

Epistasis shapes the fitness landscape of an allosteric specificity switch

Epistasis is a major determinant in the emergence of novel protein function. In allosteric proteins, direct interactions between inducer-binding mutations propagate through the allosteric network, manifesting as epistasis at the level of biological function. Elucidating this relationship between local interactions and their global effects is essential to understanding evolution of allosteric proteins. We integrate computational design, structural and biophysical analysis to characterize the emergence of novel inducer specificity in an allosteric transcription factor. Adaptive landscapes of different inducers of the designed mutant show that a few strong epistatic interactions constrain the number of viable sequence pathways, revealing ridges in the fitness landscape leading to new specificity. The structure of the designed mutant shows that a striking change in inducer orientation still retains allosteric function. Comparing biophysical and functional properties suggests a nonlinear relationship between inducer binding affinity and allostery. Our results highlight the functional and evolutionary complexity of allosteric proteins.

Probing Universal Protein Dynamics Using Hydrogen-Deuterium Exchange Mass Spectrometry-Derived Residue-Level Gibbs Free Energy.

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is a powerful technique to monitor protein intrinsic dynamics. The technique provides high-resolution information on how protein intrinsic dynamics are altered in response to biological signals, such as ligand binding, oligomerization, or allosteric networks. However, identification, interpretation, and visualization of such events from HDX-MS data sets is challenging as these data sets consist of many individual data points collected across peptides, time points, and experimental conditions. Here, we present PyHDX, an open-source Python package and webserver, that allows the user to batch extract the universal quantity Gibbs free energy at residue levels over multiple protein conditions and homologues. The output is directly visualized on a linear map or 3D structures or is exported as .csv files or PyMOL scripts.

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.

Letting go: Deep computational modeling insights into pH-dependent calcium affinity

Calcium and other cofactors can feature as key additions to a molecular interface, to the extent that the cofactor is completely buried in the bound state. How can such an interaction be regulated then? The answer: By facilitating a switch through an allosteric network. Although a number of unbinding mechanisms are being characterized, an extensive computational study by Joswig et al. reveals a detailed model for the pattern recognition receptor langerin.

Allosteric Mechanism of Human Mitochondrial Phenylalanyl-tRNA Synthetase: An Atomistic MD Simulation and a Mutual Information-Based Network Study

Aminoacyl-tRNA synthetases (aaRSs), a family of ubiquitous and essential enzymes, can bind target tRNAs and catalyze the aminoacylation reaction in genetic code translation. In this work, we explore the dynamic properties and allosteric communication of human mitochondrial phenylalanyl-tRNA synthetase (hmPheRS) in free and bound states to understand the mechanisms of its tRNAPhe recognition and allostery using molecular dynamics simulations combined with the torsional mutual information-based network model. Our results reveal that hmPheRS's residue mobility and inter-residue motional coupling are significantly enhanced by tRNAPhe binding, and there occurs a strong allosteric communication which is critical for the aminoacylation reaction, suggesting the vital role of tRNAPhe binding in the enzyme's function. The identified signaling pathways mainly make the connections between the anticodon binding domain (ABD) and catalytic domain (CAD), as well as within the CAD composed of many functional fragments and active sites, revealing the co-regulation role of them to act coordinately and achieve hmPheRS's aminoacylation function. Besides, several key residues along the communication pathways are identified to be involved in mediating the coordinated coupling between anticodon recognition at the ABD and activation process at the CAD, showing their pivotal role in the allosteric network, which are well consistent with the experimental observation. This study sheds light on the allosteric communication mechanism in hmPheRS and can provide important information for the structure-based drug design targeting aaRSs.

Structural insights into the mechanism of human methyltransferase hPRMT4

Human PRMT4, also known as CARM1, is a type I arginine methyltransferase protein that catalyse the formation of asymmetrical dimethyl arginine product. Structural studies done to date on PRMT4 have shown that the N-terminal region, Rossmann fold and dimerization arm play an important role in PRMT4 activity. Elucidating the functions of these regions in catalysis remains to be explored. Studies have shown the existence of communication pathways in PRMT4, which need further elucidation. The molecular dynamics (MD) simulations performed in this study show differences in different monomeric and dimeric forms of hPRMT4, revealing the role of the N-terminal region, Rossmann fold and dimerization arm. The study shows the conformational changes that occur during dimerization and SAM binding. Our cross-correlation analysis showed a correlation between these regions. Further, we performed PSN and network analysis to establish the existence of communication networks and an allosteric pathway. This study shows the use of MD simulations and network analysis to explore the aspects of PRMT4 dimerization, SAM binding and demonstrates the existence of an allosteric network. These findings shed novel insights into the conformational changes associated with hPRMT4, the mechanism of its dimerization, SAM binding and clues for better inhibitor designs. Communicated by Ramaswamy H. Sarma

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

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

Local rules for fabricating allosteric networks

Mechanical properties of disordered networks can be significantly tailored by modifying a small fraction of their bonds. This procedure has been used to design and build mechanical metamaterials with a variety of responses. A long-range 'allosteric' response, where a localized input strain at one site gives rise to a localized output strain at a distant site, has been of particular interest. This work presents a novel approach to incorporating allosteric responses in experimental systems by pruning disordered networks $\textit{in-situ}$. Previous work has relied on computer simulations to design and predict the response of such systems using a cost function where the response of the entire network to each bond removal is used at each step to determine which bond to prune. It is not feasible to follow such a design protocol in experiments where one has access only to local response at each site. This paper presents design algorithms that allow determination of what bonds to prune based purely on the local stresses in the network without employing a cost function; using only local information, allosteric networks are designed in simulations and then built out of real materials. The results show that some pruning strategies work better than others when translated into an experimental system. A method is presented to measure local stresses experimentally in disordered networks. This approach is then used to implement pruning methods to design desired responses $\textit{in-situ}$. Results from these experiments confirm that the pruning methods are robust and work in a real laboratory material.

Free energy calculations of ALS-causing SOD1 mutants reveal common perturbations to stability and dynamics along the maturation pathway.

With over 150 heritable mutations identified as disease-causative, superoxide dismutase 1 (SOD1) has been a main target of amyotrophic lateral sclerosis (ALS) research and therapeutic efforts. However, recent evidence has suggested that neither loss of function nor protein aggregation is responsible for promoting neurotoxicity. Furthermore, there is no clear pattern to the nature or the location of these mutations that could suggest a molecular mechanism behind SOD1-linked ALS. Here, we utilize reliable and accurate computational techniques to predict the perturbations of 10 such mutations to the free energy changes of SOD1 as it matures from apo monomer to metallated dimer. We find that the free energy perturbations caused by these mutations strongly depend on maturational progress, indicating the need for state-specific therapeutic targeting. We also find that many mutations exhibit similar patterns of perturbation to native and non-native maturation, indicating strong thermodynamic coupling between the dynamics at various sites of maturation within SOD1. These results suggest the presence of an allosteric network in SOD1 which is vulnerable to disruption by these mutations. Analysis of these perturbations may contribute to uncovering a unifying molecular mechanism which explains SOD1-linked ALS and help to guide future therapeutic efforts.

Activation mechanism of Drosophila cryptochrome through an allosteric switch.

Cryptochromes are signaling proteins activated by photoexcitation of the flavin adenine dinucleotide (FAD) cofactor. Although extensive research has been performed, the mechanism for this allosteric process is still unknown. We constructed three computational models, corresponding to different redox states of the FAD cofactor in Drosophila cryptochrome (dCRY). Analyses of the dynamics trajectories reveal that the activation process occurs in the semiquinone state FAD-●, resulting from excited-state electron transfer. The Arg381-Asp410 salt bridge acts as an allosteric switch, regulated by the change in the redox state of FAD. In turn, Asp410 forms new hydrogen bonds, connecting allosteric networks of the amino-terminal and carboxyl-terminal domains initially separated in the resting state. The expansion to a global dynamic network leads to enhanced protein fluctuations, an increase in the radius of gyration, and the expulsion of the carboxyl-terminal tail. These structural features are in accord with mutations and spectroscopic experiments.

A putative structural mechanism underlying the antithetic effect of homologous RND1 and RhoD GTPases in mammalian plexin regulation.

Plexins are semaphorin receptors that play essential roles in mammalian neuronal axon guidance and in many other important mammalian biological processes. Plexin signaling depends on a semaphorin-induced dimerization mechanism and is modulated by small GTPases of the Rho family, of which RND1 serves as a plexin activator yet its close homolog RhoD an inhibitor. Using molecular dynamics (MD) simulations, we showed that RND1 reinforces the plexin dimerization interface, whereas RhoD destabilizes it due to their differential interaction with the cell membrane. Upon binding plexin at the Rho-GTPase-binding domain (RBD), RND1 and RhoD interact differently with the inner leaflet of the cell membrane and exert opposite effects on the dimerization interface via an allosteric network involving the RBD, RBD linkers, and a buttress segment adjacent to the dimerization interface. The differential membrane interaction is attributed to the fact that, unlike RND1, RhoD features a short C-terminal tail and a positively charged membrane interface.

Substitution of a Surface-Exposed Residue Involved in an Allosteric Network Enhances Tryptophan Synthase Function in Cells.

Networks of noncovalent amino acid interactions propagate allosteric signals throughout proteins. Tryptophan synthase (TS) is an allosterically controlled bienzyme in which the indole product of the alpha subunit (αTS) is transferred through a 25 Å hydrophobic tunnel to the active site of the beta subunit (βTS). Previous nuclear magnetic resonance and molecular dynamics simulations identified allosteric networks in αTS important for its function. We show here that substitution of a distant, surface-exposed network residue in αTS enhances tryptophan production, not by activating αTS function, but through dynamically controlling the opening of the indole channel and stimulating βTS activity. While stimulation is modest, the substitution also enhances cell growth in a tryptophan-auxotrophic strain of Escherichia coli compared to complementation with wild-type αTS, emphasizing the biological importance of the network. Surface-exposed networks provide new opportunities in allosteric drug design and protein engineering, and hint at potential information conduits through which the functions of a metabolon or even larger proteome might be coordinated and regulated.

Computational analysis of protein stability and allosteric interaction networks in distinct conformational forms of the SARS-CoV-2 spike D614G mutant: reconciling functional mechanisms through allosteric model of spike regulation

In this study, we used an integrative computational approach to examine molecular mechanisms underlying functional effects of the D614G mutation by exploring atomistic modeling of the SARS-CoV-2 spike proteins as allosteric regulatory machines. We combined coarse-grained simulations, protein stability and dynamic fluctuation communication analysis with network-based community analysis to examine structures of the native and mutant SARS-CoV-2 spike proteins in different functional states. Through distance fluctuations communication analysis, we probed stability and allosteric communication propensities of protein residues in the native and mutant SARS-CoV-2 spike proteins, providing evidence that the D614G mutation can enhance long-range signaling of the allosteric spike engine. By combining functional dynamics analysis and ensemble-based alanine scanning of the SARS-CoV-2 spike proteins we found that the D614G mutation can improve stability of the spike protein in both closed and open forms, but shifting thermodynamic preferences towards the open mutant form. Our results revealed that the D614G mutation can promote the increased number of stable communities and allosteric hub centers in the open form by reorganizing and enhancing the stability of the S1-S2 inter-domain interactions and restricting mobility of the S1 regions. This study provides atomistic-based view of allosteric communications in the SARS-CoV-2 spike proteins, suggesting that the D614G mutation can exert its primary effect through allosterically induced changes on stability and communications in the residue interaction networks. Communicated by Ramaswamy H. Sarma