Allosteric perspective on the mutability and druggability of the SARS-CoV-2 Spike protein

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Mechanisms of allostery at the viral surface through the eyes of molecular simulation

The outermost surface layer of any virus is formed by either a capsid shell or envelope. Such layers have traditionally been thought of as immovable structures, but it is becoming apparent that they cannot be viewed exclusively as static architectures protecting the viral genome. A limited number of proteins on the virion surface must perform a multitude of functions in order to orchestrate the viral life cycle, and allostery can regulate their structures at multiple levels of organization, spanning individual molecules, protomers, large oligomeric assemblies, or entire viral surfaces. Here, we review recent contributions from the molecular simulation field to viral surface allostery, with a particular focus on the trimeric spike glycoprotein emerging from the coronavirus surface, and the icosahedral flaviviral envelope complex. As emerging viral pathogens continue to pose a global threat, an improved understanding of viral dynamics and allosteric regulation will prove crucial in developing novel therapeutic strategies.

Allosteric drugs: New principles and design approaches

Focusing on an important biomedical implication of allostery - design of allosteric drugs, we describe characteristics of allosteric sites, effectors, and their modes of actions distinguishing them from the orthosteric counterparts and calling for new principles and protocols in the quests for allosteric drugs. We show the importance of considering both binding affinity and allosteric signaling in establishing the structure-activity relationships (SARs) toward design of allosteric effectors, arguing that pairs of allosteric sites and their effector ligands - the site-effector pairs - should be generated and adjusted simultaneously in the framework of what we call directed design protocol. Key ideas and approaches for designing allosteric effectors including reverse perturbation, targeted and agnostic analysis are also discussed here. Several promising computational approaches are highlighted, along with the need for and potential advantages of utilizing generative models to facilitate discovery/design of new allosteric drugs.

Allosteric Signal within the Receptor-Binding Domain of the SARS-CoV-2 Spike Protein Mediated by a Class 3 Monoclonal Antibody Revealed through Molecular Dynamics Simulations and Protein Residue Networks

This study investigated the allosteric action within the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein caused by class 3 monoclonal antibody (mAb) binding. As the emergence of SARS-CoV-2 variants has raised concerns about the effectiveness of treatments by antibodies, targeting the highly conserved class 3 epitopes has become an alternative strategy of antibody design. Simulations of explicitly solvated RBD of the BA.2.75 omicron subvariants were carried out both in the presence and in the absence of bebtelovimab, as a model example of class 3 monoclonal antibodies against the RBD of the SARS-CoV-2 spike protein. The comparative analysis showed that bebtelovimab's binding on two α helices at the epitope region disrupted the nearby interaction network, which triggered a denser interaction network formation on the opposite side of the receptor-binding motif (RBM) region and resulted in a "close" conformation that could prevent the ACE2 binding. A better understanding of this allosteric action could lead to the development of alternative mAbs for further variants of concern. In terms of computational techniques, the communicability matrix could serve as a tool to visualize the effects of allostery, as the pairs of amino acids or secondary structures with high communicability could pinpoint the possible sites to transfer the allosteric signal. Additionally, the communicability gain/loss matrix could help elucidate the consequences of allosteric actions, which could be employed along with other allostery quantification techniques in some previous studies.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 Spike (S) protein. With this approach, we first identified Candidate Adaptive Polymorphisms (CAPs) of the SARS-CoV-2 Spike protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

ASD2023: towards the integrating landscapes of allosteric knowledgebase

Allosteric regulation, induced by perturbations at an allosteric site topographically distinct from the orthosteric site, is one of the most direct and efficient ways to fine-tune macromolecular function. The Allosteric Database (ASD; accessible online at http://mdl.shsmu.edu.cn/ASD) has been systematically developed since 2009 to provide comprehensive information on allosteric regulation. In recent years, allostery has seen sustained growth and wide-ranging applications in life sciences, from basic research to new therapeutics development, while also elucidating emerging obstacles across allosteric research stages. To overcome these challenges and maintain high-quality data center services, novel features were curated in the ASD2023 update: (i) 66 589 potential allosteric sites, covering > 80% of the human proteome and constituting the human allosteric pocketome; (ii) 748 allosteric protein-protein interaction (PPI) modulators with clear mechanisms, aiding protein machine studies and PPI-targeted drug discovery; (iii) 'Allosteric Hit-to-Lead,' a pioneering dataset providing panoramic views from 87 well-defined allosteric hits to 6565 leads and (iv) 456 dualsteric modulators for exploring the simultaneous regulation of allosteric and orthosteric sites. Meanwhile, ASD2023 maintains a significant growth of foundational allosteric data. Based on these efforts, the allosteric knowledgebase is progressively evolving towards an integrated landscape, facilitating advancements in allosteric target identification, mechanistic exploration and drug discovery.

Discovery of a Potential Allosteric Site in the SARS-CoV-2 Spike Protein and Targeting Allosteric Inhibitor to Stabilize the RBD Down State using a Computational Approach

BACKGROUND

The novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to a worldwide public health crisis. At present, the development of effective drugs and/or related therapeutics is still the most urgent and important task for combating the virus. The viral entry and associated infectivity mainly rely on its envelope spike protein to recognize and bind to the host cell receptor angiotensin-converting enzyme 2 (ACE2) through a conformational switch of the spike receptor binding domain (RBD) from inactive to active state. Thus, it is of great significance to design an allosteric inhibitor targeting spike to lock it in the inactive and ACE2-inaccessible state.

OBJECTIVE

This study aims to discover the potential broad-spectrum allosteric inhibitors capable of binding and stabilizing the diverse spike variants, including the wild type, Delta, and Omicron, in the inactive RBD down state.

METHODS

In this work, we first detected a potential allosteric pocket within the SARS-CoV-2 spike protein. Then, we performed large-scale structure-based virtual screening by targeting the putative allosteric pocket to identify allosteric inhibitors that could stabilize the spike inactive state. Molecular dynamics simulations were further carried out to evaluate the effects of compound binding on the stability of spike RBD.

RESULTS

Finally, we identified three potential allosteric inhibitors, CPD3, CPD5, and CPD6, against diverse SARS-CoV-2 variants, including Wild-type, Delta, and Omicron variants. Our simulation results showed that the three compounds could stably bind the predicted allosteric site and effectively stabilize the spike in the inactive state.

CONCLUSION

The three compounds provide novel chemical structures for rational drug design targeting spike protein, which is expected to greatly assist in the development of new drugs against SARS-CoV-2.

Progressive computational approaches to facilitate decryption of allosteric mechanism and drug discovery

Allostery is a ubiquitous biological phenomenon where perturbation at topologically distal areas of a protein serves as a trigger to fine-tune the orthosteric site and thus regulate protein function. The investigation of allosteric regulation greatly enhances our understanding of human diseases and broadens avenue for drug discovery. For decades, owing to the difficulty in allostery characterization through serendipitous experimental screening, researchers have developed several innovative computational approaches, which proves to accelerate the elucidation of allostery. Herein, we review the state-of-the-art advance of computational methodologies for allostery study, with particular emphasis on promising trends emerging over the past two years. We expect this review will outline the latest landscape of allostery study and inspire researchers to further facilitate this field.

Inconspicuous Yet Indispensable: The Coronavirus Spike Transmembrane Domain

Membrane-spanning portions of proteins' polypeptide chains are commonly known as their transmembrane domains (TMDs). The structural organisation and dynamic behaviour of TMDs from proteins of various families, be that receptors, ion channels, enzymes etc., have been under scrutiny on the part of the scientific community for the last few decades. The reason for such attention is that, apart from their obvious role as an "anchor" in ensuring the correct orientation of the protein's extra-membrane domains (in most cases functionally important), TMDs often actively and directly contribute to the operation of "the protein machine". They are capable of transmitting signals across the membrane, interacting with adjacent TMDs and membrane-proximal domains, as well as with various ligands, etc. Structural data on TMD arrangement are still fragmentary at best due to their complex molecular organisation as, most commonly, dynamic oligomers, as well as due to the challenges related to experimental studies thereof. Inter alia, this is especially true for viral fusion proteins, which have been the focus of numerous studies for quite some time, but have provoked unprecedented interest in view of the SARS-CoV-2 pandemic. However, despite numerous structure-centred studies of the spike (S) protein effectuating target cell entry in coronaviruses, structural data on the TMD as part of the entire spike protein are still incomplete, whereas this segment is known to be crucial to the spike's fusogenic activity. Therefore, in attempting to bring together currently available data on the structure and dynamics of spike proteins' TMDs, the present review aims to tackle a highly pertinent task and contribute to a better understanding of the molecular mechanisms underlying virus-mediated fusion, also offering a rationale for the design of novel efficacious methods for the treatment of infectious diseases caused by SARS-CoV-2 and related viruses.

Synthesis and in silico study of chenodeoxycholic acid and its analogues as an alternative inhibitor of spike glycoprotein of SARS-CoV-2

COVID-19, caused by SARS-CoV-2, is a viral infection that has generated one of the most significant health problems in the world. Spike glycoprotein is a crucial enzyme in viral replication and transcription mediation. There are reports in the literature on using bile acid in the fight against this virus through in vitro tests. This work presents the synthesis of nine chenodeoxycholic acid derivatives (1-9), which were prepared by oxidation, acetylation, formylation, and esterification reactions, and the analogs 6-9 have not yet been reported in the literature and the possibility of conducting an in silico study of bile acid derivatives as a therapeutic alternative to combat the virus using glycoprotein as a macromolecular target. As a result, five compounds (1, 6-9) possessed favorable competitive interactions with the lowest energies compared to the native ligand (BLA), and the highlighted compound 9 got the best scores. At the same time, analog 1 presented the best ADME filter result. Molecular dynamics also simulated these compounds to verify their stability within the active protein site to seek new therapeutic propositions to fight against the pandemic. Physical and spectroscopic data have fully characterized all the compounds.Communicated by Ramaswamy H. Sarma.

An innovative strategy to investigate microbial protein modifications in a reliable fast and sensitive way: A therapy oriented proof of concept based on UV-C irradiation of SARS-CoV-2 spike protein

The characterization of modifications of microbial proteins is of primary importance to dissect pathogen lifecycle mechanisms and could be useful in identifying therapeutic targets. Attempts to solve this issue yielded only partial and non-exhaustive results. We developed a multidisciplinary approach by coupling in vitro infection assay, mass spectrometry (MS), protein 3D modelling, and surface plasma resonance (SPR). As a proof of concept, the effect of low UV-C (273 nm) irradiation on SARS-CoV-2 spike (S) protein was investigated. Following UV-C exposure, MS analysis identified, among other modifications, the disruption of a disulphide bond within the conserved S2 subunit of S protein. Computational analyses revealed that this bond breakage associates with an allosteric effect resulting in the generation of a closed conformation with a reduced ability to bind the ACE2 receptor. The UV-C-induced reduced affinity of S protein for ACE2 was further confirmed by SPR analyses and in vitro infection assays. This comprehensive approach pinpoints the S2 domain of S protein as a potential therapeutic target to prevent SARS-CoV-2 infection. Notably, this workflow could be used to screen a wide variety of microbial protein domains, resulting in a precise molecular fingerprint and providing new insights to adequately address future epidemics.

SARS-CoV-2 variant evasion of monoclonal antibodies based on in vitro studies

Monoclonal antibodies (mAbs) offer a treatment option for individuals with severe COVID-19 and are especially important in high-risk individuals where vaccination is not an option. Given the importance of understanding the evolution of resistance to mAbs by SARS-CoV-2, we reviewed the available in vitro neutralization data for mAbs against live variants and viral constructs containing spike mutations of interest. Unfortunately, evasion of mAb-induced protection is being reported with new SARS-CoV-2 variants. The magnitude of neutralization reduction varied greatly among mAb-variant pairs. For example, sotrovimab retained its neutralization capacity against Omicron BA.1 but showed reduced efficacy against BA.2, BA.4 and BA.5, and BA.2.12.1. At present, only bebtelovimab has been reported to retain its efficacy against all SARS-CoV-2 variants considered here. Resistance to mAb neutralization was dominated by the action of epitope single amino acid substitutions in the spike protein. Although not all observed epitope mutations result in increased mAb evasion, amino acid substitutions at non-epitope positions and combinations of mutations also contribute to evasion of neutralization. This Review highlights the implications for the rational design of viral genomic surveillance and factors to consider for the development of novel mAb therapies.

Urgent Need for Next-Generation COVID-19 Vaccines

This Viewpoint argues that the development of a distinctly improved generation of SARS-CoV-2 vaccines is paramount to offering a greater breadth and depth of protection for a longer duration against COVID-19 disease.

Green and efficient one-pot three-component synthesis of novel drug-like furo[2,3–d]pyrimidines as potential active site inhibitors and putative allosteric hotspots modulators of both SARS-CoV-2 MPro and PLPro

In this paper, an environmentally benign, convenient, and efficient one-pot three-component reaction has been developed for the regioselective synthesis of novel 5-aroyl(or heteroaroyl)-6-(alkylamino)-1,3-dimethylfuro[2,3-d]pyrimidine-2,4(1H,3H)-diones (4a‒n) through the sequential condensation of aryl(or heteroaryl)glyoxal monohydrates (1a‒g), 1,3-dimethylbarbituric acid (2), and alkyl(viz. cyclohexyl or tert-butyl)isocyanides (3a or 3b) catalyzed by ultra-low loading ZrOCl₂•8H₂O (just 2 mol%) in water at 50 ˚C. After synthesis and characterization of the mentioned furo[2,3-d]pyrimidines (4a‒n), their multi-targeting inhibitory properties were investigated against the active site and putative allosteric hotspots of both SARS-CoV-2 main protease (M^Pro) and papain-like protease (PL^Pro) based on molecular docking studies and compare the attained results with various medicinal compounds which approximately in three past years were used, introduced, and or repurposed to fight against COVID-19. Furthermore, drug-likeness properties of the mentioned small heterocyclic frameworks (4a‒n) have been explored using in silico ADMET analyses. Interestingly, the molecular docking studies and ADMET-related data revealed that the novel series of furo[2,3-d]pyrimidines (4a‒n), especially 5-(3,4-methylendioxybenzoyl)-6-(cyclohexylamino)-1,3-dimethylfuro[2,3-d]pyrimidine-2,4(1H,3H)-dione (4g) as hit one is potential COVID-19 drug candidate, can subject to further in vitro and in vivo studies. It is worthwhile to note that the protein-ligand-type molecular docking studies on the human body temperature-dependent M^Pro protein that surprisingly contains zinc^II (Zn^II) ion between His41/Cys145 catalytic dyad in the active site, which undoubtedly can make new plans for designing novel SARS-CoV-2 M^Pro inhibitors, is performed for the first time in this paper, to the best of our knowledge.

AlloMAPS 2: allosteric fingerprints of the AlphaFold and Pfam-trRosetta predicted structures for engineering and design

AlloMAPS 2 is an update of the Allosteric Mutation Analysis and Polymorphism of Signalling database, which contains data on allosteric communication obtained for predicted structures in the AlphaFold database (AFDB) and trRosetta-predicted Pfam domains. The data update contains Allosteric Signalling Maps (ASMs) and Allosteric Probing Maps (APMs) quantifying allosteric effects of mutations and of small probe binding, respectively. To ensure quality of the ASMs and APMs, we performed careful and accurate selection of protein sets containing high-quality predicted structures in both databases for each organism/structure, and the data is available for browsing and download. The data for remaining structures are available for download and should be used at user's discretion and responsibility. We believe these massive data can facilitate both diagnostics and drug design within the precision medicine paradigm. Specifically, it can be instrumental in the analysis of allosteric effects of pathological and rescue mutations, providing starting points for fragment-based design of allosteric effectors. The exhaustive character of allosteric signalling and probing fingerprints will be also useful in future developments of corresponding machine learning applications. The database is freely available at: http://allomaps.bii.a-star.edu.sg.

Allosteric Hotspots in the Main Protease of SARS-CoV-2

Inhibiting the main protease of SARS-CoV-2 is of great interest in tackling the COVID-19 pandemic caused by the virus. Most efforts have been centred on inhibiting the binding site of the enzyme. However, considering allosteric sites, distant from the active or orthosteric site, broadens the search space for drug candidates and confers the advantages of allosteric drug targeting. Here, we report the allosteric communication pathways in the main protease dimer by using two novel fully atomistic graph-theoretical methods: Bond-to-bond propensity, which has been previously successful in identifying allosteric sites in extensive benchmark data sets without a priori knowledge, and Markov transient analysis, which has previously aided in finding novel drug targets in catalytic protein families. Using statistical bootstrapping, we score the highest ranking sites against random sites at similar distances, and we identify four statistically significant putative allosteric sites as good candidates for alternative drug targeting.

Evolutionary progression of collective mutations in Omicron sub-lineages towards efficient RBD-hACE2: Allosteric communications between and within viral and human proteins

The interaction between the Spike (S) protein of SARS-CoV-2 and the human angiotensin converting enzyme 2 (hACE2) is essential for infection, and is a target for neutralizing antibodies. Consequently, selection of mutations in the S protein is expected to be driven by the impact on the interaction with hACE2 and antibody escape. Here, for the first time, we systematically characterized the collective effects of mutations in each of the Omicron sub-lineages (BA.1, BA.2, BA.3 and BA.4) on both the viral S protein receptor binding domain (RBD) and the hACE2 protein using post molecular dynamics studies and dynamic residue network (DRN) analysis. Our analysis suggested that Omicron sub-lineage mutations result in altered physicochemical properties that change conformational flexibility compared to the reference structure, and may contribute to antibody escape. We also observed changes in the hACE2 substrate binding groove in some sub-lineages. Notably, we identified unique allosteric communication paths in the reference protein complex formed by the DRN metrics betweenness centrality and eigencentrality hubs, originating from the RBD core traversing the receptor binding motif of the S protein and the N-terminal domain of the hACE2 to the active site. We showed allosteric changes in residue network paths in both the RBD and hACE2 proteins due to Omicron sub-lineage mutations. Taken together, these data suggest progressive evolution of the Omicron S protein RBD in sub-lineages towards a more efficient interaction with the hACE2 receptor which may account for the increased transmissibility of Omicron variants.

A Uniquely Stable Trimeric Model of SARS-CoV-2 Spike Transmembrane Domain

Understanding fusion mechanisms employed by SARS-CoV-2 spike protein entails realistic transmembrane domain (TMD) models, while no reliable approaches towards predicting the 3D structure of transmembrane (TM) trimers exist. Here, we propose a comprehensive computational framework to model the spike TMD only based on its primary structure. We performed amino acid sequence pattern matching and compared the molecular hydrophobicity potential (MHP) distribution on the helix surface against TM homotrimers with known 3D structures and selected an appropriate template for homology modeling. We then iteratively built a model of spike TMD, adjusting “dynamic MHP portraits” and residue variability motifs. The stability of this model, with and without palmitoyl modifications downstream of the TMD, and several alternative configurations (including a recent NMR structure), was tested in all-atom molecular dynamics simulations in a POPC bilayer mimicking the viral envelope. Our model demonstrated unique stability under the conditions applied and conforms to known basic principles of TM helix packing. The original computational framework looks promising and could potentially be employed in the construction of 3D models of TM trimers for a wide range of membrane proteins.

Uncovering cryptic pockets in the SARS-CoV-2 spike glycoprotein

The COVID-19 pandemic has prompted a rapid response in vaccine and drug development. Herein, we modeled a complete membrane-embedded SARS-CoV-2 spike glycoprotein and used molecular dynamics simulations with benzene probes designed to enhance discovery of cryptic pockets. This approach recapitulated lipid and host metabolite binding sites previously characterized by cryo-electron microscopy, revealing likely ligand entry routes, and uncovered a novel cryptic pocket with promising druggable properties located underneath the 617-628 loop. A full representation of glycan moieties was essential to accurately describe pocket dynamics. A multi-conformational behavior of the 617-628 loop in simulations was validated using hydrogen-deuterium exchange mass spectrometry experiments, supportive of opening and closing dynamics. The pocket is the site of multiple mutations associated with increased transmissibility found in SARS-CoV-2 variants of concern including Omicron. Collectively, this work highlights the utility of the benzene mapping approach in uncovering potential druggable sites on the surface of SARS-CoV-2 targets.

Optimising Elastic Network Models for Protein Dynamics and Allostery: Spatial and Modal Cut-offs and Backbone Stiffness

The family of coarse-grained models for protein dynamics known as Elastic Network Models (ENMs) require careful choice of parameters to represent well experimental measurements or fully-atomistic simulations. The most basic ENM that represents each protein residue by a node at the position of its C-alpha atom, all connected by springs of equal stiffness, up to a cut-off in distance. Even at this level a choice is required of the optimum cut-off distance and the upper limit of elastic normal modes taken in any sum for physical properties, such as dynamic correlation or allosteric effects on binding. Additionally, backbone-enhanced ENM (BENM) may improve the model by allocating a higher stiffness to springs that connect along the protein backbone. This work reports on the effect of varying these three parameters (distance and mode cutoffs, backbone stiffness) on the dynamical structure of three proteins, Catabolite Activator Protein (CAP), Glutathione S-transferase (GST), and the SARS-CoV-2 Main Protease (M pro ). Our main results are: (1) balancing B-factor and dispersion-relation predictions, a near-universal optimal value of 8.5 Å is advisable for ENMs; (2) inhomogeneity in elasticity brings the first mode containing spatial structure not well-resolved by the ENM typically within the first 20; (3) the BENM only affects modes in the upper third of the distribution, and, additionally to the ENM, is only able to model the dispersion curve better in this vicinity; (4) BENM does not typically affect fluctuation-allostery, which also requires careful treatment of the effector binding to the host protein to capture.

Conservation and diversity in allosteric fingerprints of proteins for evolutionary-inspired engineering and design

Hand-in-hand work of physics and evolution delivered protein universe with diversity of forms, sizes, and functions. Pervasiveness and advantageous traits of allostery made it an important component of the protein function regulation, calling for thorough investigation of its structural determinants and evolution. Learning directly from nature, we explored here allosteric communication in several major folds and repeat proteins, including α/β and β-barrels, β-propellers, Ig-like fold, ankyrin and α/β leucine-rich repeat proteins, which provide structural platforms for many different enzymatic and signalling functions. We obtained a picture of conserved allosteric communication characteristic in different fold types, modifications of the structure-driven signalling patterns via sequence-determined divergence to specific functions, as well as emergence and potential diversification of allosteric regulation in multi-domain proteins and oligomeric assemblies. Our observations will be instrumental in facilitating the engineering and de novo design of proteins with allosterically regulated functions, including development of therapeutic biologics. In particular, results described here may guide the identification of the optimal structural platforms (e.g. fold type, size, and oligomerization states) and the types of diversifications/perturbations, such as mutations, effector binding, and order-disorder transition. The tunable allosteric linkage across distant regions can be used as a pivotal component in the design/engineering of modular biological systems beyond the traditional scaffolding function.