Allosteric sites: remote control in regulation of protein activity

GETdb: A comprehensive database for genetic and evolutionary features of drug targets

The development of an innovative drug is complex and time-consuming, and the drug target identification is one of the critical steps in drug discovery process. Effective and accurate identification of drug targets can accelerate the drug development process. According to previous research, evolutionary and genetic information of genes has been found to facilitate the identification of approved drug targets. In addition, allosteric proteins have great potential as targets due to their structural diversity. However, this information that could facilitate target identification has not been collated in existing drug target databases. Here, we construct a comprehensive drug target database named Genetic and Evolutionary features of drug Targets database (GETdb, http://zhanglab.hzau.edu.cn/GETdb/page/index.jsp). This database not only integrates and standardizes data from dozens of commonly used drug and target databases, but also innovatively includes the genetic and evolutionary information of targets. Moreover, this database features an effective allosteric protein prediction model. GETdb contains approximately 4000 targets and over 29,000 drugs, and is a user-friendly database for searching, browsing and downloading data to facilitate the development of novel targets.

Megalin-related mechanism of hemolysis-induced acute kidney injury and the therapeutic strategy

Hemolysis-induced acute kidney injury (AKI) is attributed to heme-mediated proximal tubule epithelial cell (PTEC) injury and tubular cast formation due to intratubular protein condensation. Megalin is a multiligand endocytic receptor for proteins, peptides, and drugs in PTECs and mediates the uptake of free hemoglobin and the heme-scavenging protein α1-microglobulin. However, understanding of how megalin is involved in the development of hemolysis-induced AKI remains elusive. Here, we investigated the megalin-related pathogenesis of hemolysis-induced AKI and a therapeutic strategy using cilastatin, a megalin blocker. A phenylhydrazine-induced hemolysis model developed in kidney-specific mosaic megalin knockout (MegKO) mice confirmed megalin-dependent PTEC injury revealed by the co-expression of kidney injury molecule-1 (KIM-1). In the hemolysis model in kidney-specific conditional MegKO mice, the uptake of hemoglobin and α1-microglobulin as well as KIM-1 expression in PTECs was suppressed, but tubular cast formation was augmented, likely due to the nonselective inhibition of protein reabsorption in PTECs. Quartz crystal microbalance analysis revealed that cilastatin suppressed the binding of megalin with hemoglobin and α1-microglobulin. Cilastatin also inhibited the specific uptake of fluorescent hemoglobin by megalin-expressing rat yolk sac tumor-derived L2 cells. In a mouse model of hemolysis-induced AKI, repeated cilastatin administration suppressed PTEC injury by inhibiting the uptake of hemoglobin and α1-microglobulin and also prevented cast formation. Hemopexin, another heme-scavenging protein, was also found to be a novel ligand of megalin, and its binding to megalin and uptake by PTECs in the hemolysis model were suppressed by cilastatin. Mass spectrometry-based semiquantitative analysis of urinary proteins in cilastatin-treated C57BL/6J mice indicated that cilastatin suppressed the reabsorption of a limited number of megalin ligands in PTECs, including α1-microglobulin and hemopexin. Collectively, cilastatin-mediated selective megalin blockade is an effective therapeutic strategy to prevent both heme-mediated PTEC injury and cast formation in hemolysis-induced AKI. © 2024 The Pathological Society of Great Britain and Ireland.

Ethaverine and Papaverine Target Cyclin-Dependent Kinase 5 and Inhibit Lung Cancer Cell Proliferation and Migration

CDK5 kinase plays a central role in the regulation of neuronal functions, and its hyperactivation has been associated with neurodegenerative pathologies and more recently with several human cancers, in particular lung cancer. However, ATP-competitive inhibitors targeting CDK5 are poorly selective and suffer limitations, calling for new classes of inhibitors. In a screen for allosteric modulators of CDK5, we identified ethaverine and closely related derivative papaverine and showed that they inhibit cell proliferation and migration of non small cell lung cancer cell lines. Moreover the efficacy of these compounds is significantly enhanced when combined with the ATP-competitive inhibitor roscovitine, suggesting an additive dual mechanism of inhibition targeting CDK5. These compounds do not affect CDK5 stability, but thermodenaturation studies performed with A549 cell extracts infer that they interact with CDK5 in cellulo. Furthermore, the inhibitory potentials of ethaverine and papaverine are reduced in A549 cells treated with siRNA directed against CDK5. Taken together, our results provide unexpected and novel evidence that ethaverine and papaverine constitute promising leads that can be repurposed for targeting CDK5 in lung cancer.

Back in time to the Gly-rich prototype of the phosphate binding elementary function

Binding of nucleotides and their derivatives is one of the most ancient elementary functions dating back to the Origin of Life. We review here the works considering one of the key elements in binding of (di)nucleotide-containing ligands - phosphate binding. We start from a brief discussion of major participants, conditions, and events in prebiotic evolution that resulted in the Origin of Life. Tracing back to the basic functions, including metal and phosphate binding, and, potentially, formation of primitive protein-protein interactions, we focus here on the phosphate binding. Critically assessing works on the structural, functional, and evolutionary aspects of phosphate binding, we perform a simple computational experiment reconstructing its most ancient and generic sequence prototype. The profiles of the phosphate binding signatures have been derived in form of position-specific scoring matrices (PSSMs), their peculiarities depending on the type of the ligands have been analyzed, and evolutionary connections between them have been delineated. Then, the apparent prototype that gave rise to all relevant phosphate-binding signatures had also been reconstructed. We show that two major signatures of the phosphate binding that discriminate between the binding of dinucleotide- and nucleotide-containing ligands are GxGxxG and GxxGxG, respectively. It appears that the signature archetypal for dinucleotide-containing ligands is more generic, and it can frequently bind phosphate groups in nucleotide-containing ligands as well. The reconstructed prototype's key signature GxGGxG underlies the role of glycine residues in providing flexibility and interactions necessary for binding the phosphate groups. The prototype also contains other ancient amino acids, valine, and alanine, showing versatility towards evolutionary design and functional diversification.

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.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Allosteric inhibition of myosin by phenamacril: a synergistic mechanism revealed by computational and experimental approaches

BACKGROUND

Myosin plays a crucial role in cellular processes, while its dysfunction can lead to organismal malfunction. Phenamacril (PHA), a highly species-specific and non-competitive inhibitor of myosin I (FgMyoI) from Fusarium graminearum, has been identified as an effective fungicide for controlling plant diseases caused by partial Fusarium pathogens, such as wheat scab and rice bakanae. However, the molecular basis of its action is still unclear.

RESULTS

This study used multiple computational approaches first to elucidate the allosteric inhibition mechanism of FgMyoI by PHA at the atomistic level. The results indicated the increase of adenosine triphosphate (ATP) binding affinity upon PHA binding, which might impede the release of hydrolysis products. Furthermore, simulations revealed a broadened outer cleft and a significantly more flexible interface for actin binding, accompanied by a decrease in signaling transduction from the catalytic center to the actin-binding interface. These various effects might work together to disrupt the actomyosin cycle and hinder the ability of motor to generate force. Our experimental results further confirmed that PHA reduces the enzymatic activity of myosin and its binding with actin.

CONCLUSION

Therefore, our findings demonstrated that PHA might suppress the function of myosin through a synergistic mechanism, providing new insights into myosin allostery and offering new avenues for drug/fungicide discovery targeting myosin. © 2023 Society of Chemical Industry.

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.

Quest for selective MMP9 inhibitors: a computational approach

Matrix Metalloproteinases-9 (MMP-9) is one of the important targets that play a vital role in various diseases such as cancer, Alzheimer's, arthritis, etc. Traditionally, MMP-9 inhibitors have been unable to achieve selectivity to get around this target; thereby, novel mechanisms such as inhibition of activated MMP-9 zymogen (pro-MMP-9) have been discovered. The JNJ0966 was one of the few compounds that attained the requisite selectivity by inhibiting the activation of MMP-9 zymogen (pro-MMP-9). Since JNJ0966, no other small molecules have been identified. Herein, extensive in silico studies were called upon to bolster the prospect of exploring potential candidates. The key objective of this research is to identify the potential hits from the ChEMBL database via molecular docking and dynamics approach. Protein with PDB ID: 5UE4, having a unique inhibitor in an allosteric binding pocket of MMP-9, was chosen for the study. Structure-based virtual screening and MMGBSA binding affinity calculations were performed, and five potential hits were finalized. Detailed analysis of the best-scoring molecules was performed with ADMET analysis and molecular dynamics (MD) simulation. All five hits outperformed JNJ0966 in the docking assessment, ADMET analysis, and molecular dynamics simulation. Accordingly, our research findings imply that these hits can be investigated for in vitro and in vivo studies against proMMP9 and might be explored as potential anticancer drugs. The outcome of our research might contribute in expediting the exploration of drugs that inhibits proMMP-9.Communicated by Ramaswamy H. Sarma.

Structural Analysis, Multi-Conformation Virtual Screening and Molecular Simulation to Identify Potential Inhibitors Targeting pS273R Proteases of African Swine Fever Virus

The African Swine Fever virus (ASFV) causes an infectious viral disease in pigs of all ages. The development of antiviral drugs primarily aimed at inhibition of proteases required for the proteolysis of viral polyproteins. In this study, the conformation of the pS273R protease in physiological states were investigated, virtually screened the multi-protein conformation of pS273R target proteins, combined various molecular docking scoring functions, and identified five potential drugs from the Food and Drug Administration drug library that may inhibit pS273R. Subsequent validation of the dynamic interactions of pS273R with the five putative inhibitors was achieved using molecular dynamics simulations and binding free energy calculations using the molecular mechanics/Poison-Boltzmann (Generalized Born) (MM/PB(GB)SA) surface area. These findings demonstrate that the arm domain and Thr159-Lys167 loop region of pS273R are significantly more flexible compared to the core structural domain, and the Thr159-Lys167 loop region can serve as a "gatekeeper" in the substrate channel. Leucovorin, Carboprost, Protirelin, Flavin Mononucleotide, and Lovastatin Acid all have Gibbs binding free energies with pS273R that were less than -20 Kcal/mol according to the MM/PBSA analyses. In contrast to pS273R in the free energy landscape, the inhibitor and drug complexes of pS273R showed distinct structural group distributions. These five drugs may be used as potential inhibitors of pS273R and may serve as future drug candidates for treating ASFV.

In Vitro, In Silico and Network Pharmacology Mechanistic Approach to Investigate the α-Glucosidase Inhibitors Identified by Q-ToF-LCMS from Phaleria macrocarpa Fruit Subcritical CO2 Extract

The fruit of Phaleria macrocarpa have been traditionally used as an antidiabetic remedy in Malaysia and neighbouring countries. Despite its potential for diabetes treatment, no scientific study has ever been conducted to predict the inhibitor interaction of the protein α-glucosidase identified in an extract prepared with a non-conventional extraction technique. Hence, the major aim of this research was to evaluate the in vitro antioxidant, the α-glucosidase inhibitors, and the molecular dynamic simulations of the α-glucosidase inhibitors identified by Quadrupole Time-of-Flight Liquid Chromatography Mass Spectrometry (Q-ToF-LCMS) analysis. Initially, dry fruit were processed using non-conventional and conventional extraction methods to obtain subcritical carbon dioxide extracts (SCE-1 and SCE-2) and heating under reflux extract (HRE), respectively. Subsequently, all extracts were evaluated for their in vitro antioxidative and α-glucosidase inhibitory potentials. Subsequently, the most bioactive extract (SCE-2) was subjected to Q-ToF-LCMS analysis to confirm the presence of α-glucosidase inhibitors, which were then analysed through molecular dynamic simulations and network pharmacology approaches to confirm their possible mechanism of action. The highest inhibitory effects of the 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical and α-glucosidase on SCE-2 was found as 75.36 ± 0.82% and 81.79 ± 0.82%, respectively, compared to the SCE-1 and HRE samples. The Q-ToF-LCMS analysis tentatively identified 14 potent α-glucosidase inhibitors. Finally, five identified compounds, viz., lupenone, swertianolin, m-coumaric acid, pantothenic acid, and 8-C-glucopyranosyleriodictylol displayed significant stability, compactness, stronger protein-ligand interaction up to 100 ns further confirming their potential as α-glucosidase inhibitors. Consequently, it was concluded that the SCE-2 possesses a strong α-glucosidase inhibitory effect due to the presence of these compounds. The findings of this study might prove useful to develop these compounds as alternative safe α-glucosidase inhibitors to manage diabetes more effectively.

Variable Regions of p53 Isoforms Allosterically Hard Code DNA Interaction

Allosteric regulation of protein activity pervades biology as the "second secret of life." We have been examining the allosteric regulation and mutant reactivation of the tumor suppressor protein p53. We have found that generalizing the definition of allosteric effector to include entire proteins and expanding the meaning of binding site to include the interface of a transcription factor with its DNA to be useful in understanding the modulation of protein activity. Here, we cast the variable regions of p53 isoforms as allosteric regulators of p53 interactions with its consensus DNA. We implemented molecular dynamics simulations and our lab's new techniques of molecular dynamics (MD) sectors and MD-Markov state models to investigate the effects of nine naturally occurring splice variant isoforms of p53. We find that all of the isoforms differ from wild type in their dynamic properties and how they interact with the DNA. We consider the implications of these findings on allostery and cancer treatment.

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.

Mutagenic Activation of Glutathione Peroxidase-4: Approaches toward Rational Design of Allosteric Drugs

Glutathione peroxidase 4 (GPX4) reduces lipid hydroperoxides in lipid membranes, effectively inhibiting iron-dependent cell death or ferroptosis. The upregulation of the enzyme by the mutations at residues D21 and D23 has been suggested to be associated with higher protein activity, which confers more protection against neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases. Therefore, it has become an attractive target for treating and preventing neurodegenerative diseases. However, identifying means of mimicking the beneficial effects of these mutations distant from the active site constitutes a formidable challenge in moving toward therapeutics. In this study, we explore using molecular dynamics simulations to computationally map the conformational and energetic landscape of the wild-type GPX4 protein and three mutant variants to identify the allosteric networks of the enzyme. We present the conformational dynamic profile providing the desired signature behavior of the enzyme. We also discuss the implications of these findings for drug design efforts.

AlphaFold, Artificial Intelligence (AI), and Allostery

AlphaFold has burst into our lives. A powerful algorithm that underscores the strength of biological sequence data and artificial intelligence (AI). AlphaFold has appended projects and research directions. The database it has been creating promises an untold number of applications with vast potential impacts that are still difficult to surmise. AI approaches can revolutionize personalized treatments and usher in better-informed clinical trials. They promise to make giant leaps toward reshaping and revamping drug discovery strategies, selecting and prioritizing combinations of drug targets. Here, we briefly overview AI in structural biology, including in molecular dynamics simulations and prediction of microbiota–human protein–protein interactions. We highlight the advancements accomplished by the deep-learning-powered AlphaFold in protein structure prediction and their powerful impact on the life sciences. At the same time, AlphaFold does not resolve the decades-long protein folding challenge, nor does it identify the folding pathways. The models that AlphaFold provides do not capture conformational mechanisms like frustration and allostery, which are rooted in ensembles, and controlled by their dynamic distributions. Allostery and signaling are properties of populations. AlphaFold also does not generate ensembles of intrinsically disordered proteins and regions, instead describing them by their low structural probabilities. Since AlphaFold generates single ranked structures, rather than conformational ensembles, it cannot elucidate the mechanisms of allosteric activating driver hotspot mutations nor of allosteric drug resistance. However, by capturing key features, deep learning techniques can use the single predicted conformation as the basis for generating a diverse ensemble.

Learning about allosteric drugs and ways to design them

While the accelerating quest for precision medicine requires new individually targeting and selective drugs, and the ability to work with so-called undruggable targets, the realm of allosteric drugs meeting this need remains largely uncharted. Generalizing the observations on two major drug targets with widely observed inherent allostery, GPCRs and kinases, we describe and discuss basic allosteric modes of action that are universally applicable in all types of structures and functions. Using examples of Class A GPCRs and CMGC protein kinases, we show how Allosteric Signalling and Probing Fingerprints can be used to identify potential allosteric sites and reveal effector-leads that may serve as a starting point for the development of allosteric drugs targeting these regulatory sites. A set of distinct characteristics of allosteric ligands was established, which highlights the versatility of their design and make them advantageous before their orthosteric counterparts in personalized medicine. We argue that rational design of allosteric drugs should begin with the search for latent sites or design of non-natural binding sites followed by fragment-based design of allosteric ligands and by the mutual adjustment of the site-ligand pair in order to achieve required effects. On the basis of the perturbative nature and reversibility of allosteric communication, we propose a generic protocol for computational design of allosteric effectors, enabling also the allosteric tuning of biologics, in obtaining allosteric control over protein functions.

Cooperative Protein Allosteric Transition Mediated by a Fluctuating Transmission Network

Allosteric communication between distant protein sites represents a key mechanism of biomolecular regulation and signal transduction. Compared to other processes such as protein folding, however, the dynamical evolution of allosteric transitions is still not well understood. As an example of allosteric coupling between distant protein regions, we consider the global open-closed motion of the two domains of T4 lysozyme, which is triggered by local motion in the hinge region. Combining extensive molecular dynamics simulations with a correlation analysis of interresidue contacts, we identify a network of interresidue distances that move in a concerted manner. The cooperative process originates from a cogwheel-like motion of the hydrophobic core in the hinge region, which constitutes an evolutionary conserved and flexible transmission network. Through rigid contacts and the protein backbone, the small local changes of the hydrophobic core are passed on to the distant terminal domains and lead to the emergence of a rare global conformational transition. As in an Ising-type model, the cooperativity of the allosteric transition can be explained via the interaction of local fluctuations.

Insights into Phosphorylation-Induced Protein Allostery and Conformational Dynamics of Glycogen Phosphorylase via Integrative Structural Mass Spectrometry and In Silico Modeling

Allosteric regulation plays a fundamental role in innumerable biological processes. Understanding its dynamic mechanism and impact at the molecular level is of great importance in disease diagnosis and drug discovery. Glycogen phosphorylase (GP) is a phosphoprotein responding to allosteric regulation and has significant biological importance to glycogen metabolism. Although the atomic structures of GP have been previously solved, the conformational dynamics of GP related to allostery regulation remain largely elusive due to its macromolecular size (∼196 kDa). Here, we integrated native top-down mass spectrometry (nTD-MS), hydrogen-deuterium exchange MS (HDX-MS), protection factor (PF) analysis, molecular dynamics (MD) simulations, and allostery signaling analysis to examine the structural basis and dynamics for the allosteric regulation of GP by phosphorylation. nTD-MS reveals differences in structural stability as well as oligomeric state between the unphosphorylated (GPb) and phosphorylated (GPa) forms. HDX-MS, PF analysis, and MD simulations further pinpoint the structural differences between GPb and GPa involving the binding interfaces (the N-terminal and tower-tower helices), catalytic site, and PLP-binding region. More importantly, it also allowed us to complete the missing link of the long-range communication process from the N-terminal tail to the catalytic site caused by phosphorylation. This integrative MS and in silico-based platform is highly complementary to biophysical approaches and yields valuable insights into protein structures and dynamic regulation.

Allosteric regulation of the soluble epoxide hydrolase by nitro fatty acids using a combined experimental and computational approach

The human soluble epoxide hydrolase (hsEH) is a key regulator of epoxy fatty acid (EpFA) metabolism. Inhibition of sEH can maintain endogenous levels of beneficial EpFAs and reduce the levels of their corresponding diol products, thus ameliorating a variety of pathological conditions including cardiovascular, central nervous system and metabolic diseases. The quest for orthosteric drugs that bind directly to the catalytic crevice of hsEH has been prolonged and sustained over the past decades, but the disappointing outcome of clinical trials to date warrants alternative pharmacological approaches. Previously, we have shown that hsEH can be allosterically inhibited by the endogenous electrophilic lipid 15-deoxy-Δ^12,14-Prostaglandin-J₂, via covalent adduction to two cysteines, C423 and C522. In this study, we explore the properties and behaviour of three electrophilic lipids belonging to the class of the nitro fatty acids, namely 9- and 10-nitrooleate and 10-nitrolinoleate. Biochemical and biophysical investigations revealed that, in addition to C423 and C522, nitro fatty acids can covalently bind to additional nucleophilic residues in hsEH C-terminal domain (CTD), two of which predicted in this study to be latent allosteric sites. Systematic mapping of the protein mutational space and evaluation of possible propagation pathways delineated selected residues, both in the allosteric patches and in other regions of the enzyme, envisaged to play a role on allosteric signalling. The responses elicited by the ligands on the covalent adduction sites supports future fragment-based design studies of new allosteric effectors for hsEH with increased efficacy and selectivity.

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.

Mechanistic Insights Into Co-Administration of Allosteric and Orthosteric Drugs to Overcome Drug-Resistance in T315I BCR-ABL1

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm, driven by the BCR-ABL1 fusion oncoprotein. The discovery of orthosteric BCR-ABL1 tyrosine kinase inhibitors (TKIs) targeting its active ATP-binding pocket, such as first-generation Imatinib and second-generation Nilotinib (NIL), has profoundly revolutionized the therapeutic landscape of CML. However, currently targeted therapeutics still face considerable challenges with the inevitable emergence of drug-resistant mutations within BCR-ABL1. One of the most common resistant mutations in BCR-ABL1 is the T315I gatekeeper mutation, which confers resistance to most current TKIs in use. To resolve such conundrum, co-administration of orthosteric TKIs and allosteric drugs offers a novel paradigm to tackle drug resistance. Remarkably, previous studies have confirmed that the dual targeting BCR-ABL1 utilizing orthosteric TKI NIL and allosteric inhibitor ABL001 resulted in eradication of the CML xenograft tumors, exhibiting promising therapeutic potential. Previous studies have demonstrated the cooperated mechanism of two drugs. However, the conformational landscapes of synergistic effects remain unclear, hampering future efforts in optimizations and improvements. Hence, extensive large-scale molecular dynamics (MD) simulations of wide type (WT), WT-NIL, T315I, T315I-NIL, T315I-ABL001 and T315I-ABL001-NIL systems were carried out in an attempt to address such question. Simulation data revealed that the dynamic landscape of NIL-bound BCR-ABL1 was significantly reshaped upon ABL001 binding, as it shifted from an active conformation towards an inactive conformation. The community network of allosteric signaling was analyzed to elucidate the atomistic overview of allosteric regulation within BCR-ABL1. Moreover, binding free energy analysis unveiled that the affinity of NIL to BCR-ABL1 increased by the induction of ABL001, which led to its favorable binding and the release of drug resistance. The findings uncovered the in-depth structural mechanisms underpinning dual-targeting towards T315I BCR-ABL1 to overcome its drug resistance and will offer guidance for the rational design of next generations of BCR-ABL1 modulators and future combinatory therapeutic regimens.

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

Recent developments in the SARS-CoV-2 pandemic point to its inevitable transformation into an endemic disease, urging both refinement of diagnostics for emerging variants of concern (VOCs) and design of variant-specific drugs in addition to vaccine adjustments. Exploring the structure and dynamics of the SARS-CoV-2 Spike protein, we argue that the high-mutability characteristic of RNA viruses coupled with the remarkable flexibility and dynamics of viral proteins result in a substantial involvement of allosteric mechanisms. While allosteric effects of mutations should be considered in predictions and diagnostics of new VOCs, allosteric drugs advantageously avoid escape mutations via non-competitive inhibition originating from alternative distal locations. The exhaustive allosteric signaling and probing maps presented herein provide a comprehensive picture of allostery in the spike protein, making it possible to locate potential mutations that could work as new VOC “drivers” and to determine binding patches that may be targeted by newly developed allosteric drugs.

Harnessing Reversed Allosteric Communication: A Novel Strategy for Allosteric Drug Discovery

Allostery is a fundamental and extensive mechanism of intramolecular signal transmission. Allosteric drugs possess several unique pharmacological advantages over traditional orthosteric drugs, including greater selectivity, better physicochemical properties, and lower off-target toxicity. However, owing to the complexity of allosteric regulation, experimental approaches for the development of allosteric modulators are traditionally serendipitous. Recently, the reversed allosteric communication theory has been proposed, providing a feasible tool for the unbiased detection of allosteric sites. Herein, we review the latest research on the reversed allosteric communication effect using the examples of sirtuin 6, epidermal growth factor receptor, 3-phosphoinositide-dependent protein kinase 1, and Related to A and C kinases (RAC) serine/threonine protein kinase B and recapitulate the methodologies of reversed allosteric communication strategy. The novel reversed allosteric communication strategy greatly expands the horizon of allosteric site identification and allosteric mechanism exploration and is expected to accelerate an end-to-end framework for drug discovery.

Hidden allosteric sites and De-Novo drug design

INTRODUCTION

Hidden allosteric sites are not visible in apo-crystal structures, but they may be visible in holo-structures when a certain ligand binds and maintains the ligand intended conformation. Several computational and experimental techniques have been used to investigate these hidden sites but identifying them remains a challenge.

AREAS COVERED

This review provides a summary of the many theoretical approaches for predicting hidden allosteric sites in disease-related proteins. Furthermore, promising cases have been thoroughly examined to reveal the hidden allosteric site and its modulator.

EXPERT OPINION

In the recent past, with the development in scientific techniques and bioinformatics tools, the number of drug targets for complex human diseases has significantly increased but unfortunately most of these targets are undruggable due to several reasons. Alternative strategies such as finding cryptic (hidden) allosteric sites are an attractive approach for exploitation of the discovery of new targets. These hidden sites are difficult to recognize compared to allosteric sites, mainly due to a lack of visibility in the crystal structure. In our opinion, after many years of development, MD simulations are finally becoming successful for obtaining a detailed molecular description of drug-target interaction.

Allosteric pockets and dynamic residue network hubs of falcipain 2 in mutations including those linked to artemisinin resistance

Continually emerging resistant strains of malarial parasites to current drugs present challenges. Understanding the underlying resistance mechanisms, especially those linked to allostery is, thus, highly crucial for drug design. This forms the main concern of the paper through a case study of falcipain 2 (FP-2) and its mutations, some of which are linked to artemisinin (ART) drug resistance. Here, we applied a variety of in silico approaches and tools that we developed recently, together with existing computational tools. This included novel essential dynamics and dynamic residue network (DRN) analysis algorithms. We identified six pockets demonstrating dynamic differences in the presence of some mutations. We observed striking allosteric effects in two mutant proteins. In the presence of M245I, a cryptic pocket was detected via a unique mechanism in which Pocket 2 fused with Pocket 6. In the presence of the A353T mutation, which is located at Pocket 2, the pocket became the most rigid among all protein systems analyzed. Pocket 6 was also highly stable in all cases, except in the presence of M245I mutation. The effect of ART linked mutations was more subtle, and the changes were at residue level. Importantly, we identified an allosteric communication path formed by four unique averaged BC hubs going from the mutated residue to the catalytic site and passing through the interface of three identified pockets. Collectively, we established and demonstrated that we have robust tools and a pipeline that can be applicable to the analysis of mutations.

Synthesis and kinetic characterization of hyperbolic inhibitors of human cathepsins K and S based on a succinimide scaffold

Cathepsins K and S are closely related papain-like cysteine peptidases and potential therapeutic targets for metabolic and inflammatory diseases such as osteoporosis and arthritis. Here we describe the reduction of a previously characterized succinimide (2,5-dioxopyrrolidine)-containing hyperbolic inhibitor of cathepsin K (methyl (RS)-N-[1-(4-methoxyphenyl)-2,5-dioxopyrrolidin-3-yl]glycinate), to obtain a better and more selective compound (compound 4a - methyl (2,5-dioxopyrrolidin-3-yl)glycinate), which acted as a hyperbolic mixed inhibitor/activator similar to already known allosteric effectors of cathepsin K. We then investigated the potential of the succinimide scaffold as inhibitors of cathepsins K and/or S and synthesized a library of such compounds by 1,4-addition of α-amino acid esters and related compounds to N-substituted maleimides. From the generated library, we identified the first small molecule hyperbolic inhibitors of cathepsin S (methyl ((R)-2,5-dioxopyrrolidin-3-yl)-l-threoninate (compound R-4c) and 3-{[(1S,2R,3'S)-2-hydroxycyclohexyl]amino}pyrrolidine-2,5-dione (compound (1S,2R,3'S-10)). The former acted via a similar mechanism to compound 4a, while the latter was a hyperbolic specific inhibitor of cathepsin S. Given the versatility of the scaffold, the identified compounds will be used as the basis for the development of high-affinity hyperbolic inhibitors of the individual peptidases and to explore the potential of hyperbolic inhibitors for the inhibition of cysteine cathepsins in in vitro models.

SCOUR: a stepwise machine learning framework for predicting metabolite-dependent regulatory interactions

BACKGROUND

The topology of metabolic networks is both well-studied and remarkably well-conserved across many species. The regulation of these networks, however, is much more poorly characterized, though it is known to be divergent across organisms-two characteristics that make it difficult to model metabolic networks accurately. While many computational methods have been built to unravel transcriptional regulation, there have been few approaches developed for systems-scale analysis and study of metabolic regulation. Here, we present a stepwise machine learning framework that applies established algorithms to identify regulatory interactions in metabolic systems based on metabolic data: stepwise classification of unknown regulation, or SCOUR.

RESULTS

We evaluated our framework on both noiseless and noisy data, using several models of varying sizes and topologies to show that our approach is generalizable. We found that, when testing on data under the most realistic conditions (low sampling frequency and high noise), SCOUR could identify reaction fluxes controlled only by the concentration of a single metabolite (its primary substrate) with high accuracy. The positive predictive value (PPV) for identifying reactions controlled by the concentration of two metabolites ranged from 32 to 88% for noiseless data, 9.2 to 49% for either low sampling frequency/low noise or high sampling frequency/high noise data, and 6.6-27% for low sampling frequency/high noise data, with results typically sufficiently high for lab validation to be a practical endeavor. While the PPVs for reactions controlled by three metabolites were lower, they were still in most cases significantly better than random classification.

CONCLUSIONS

SCOUR uses a novel approach to synthetically generate the training data needed to identify regulators of reaction fluxes in a given metabolic system, enabling metabolomics and fluxomics data to be leveraged for regulatory structure inference. By identifying and triaging the most likely candidate regulatory interactions, SCOUR can drastically reduce the amount of time needed to identify and experimentally validate metabolic regulatory interactions. As high-throughput experimental methods for testing these interactions are further developed, SCOUR will provide critical impact in the development of predictive metabolic models in new organisms and pathways.

ProteinLens: a web-based application for the analysis of allosteric signalling on atomistic graphs of biomolecules

The investigation of allosteric effects in biomolecular structures is of great current interest in diverse areas, from fundamental biological enquiry to drug discovery. Here we present ProteinLens, a user-friendly and interactive web application for the investigation of allosteric signalling based on atomistic graph-theoretical methods. Starting from the PDB file of a biomolecule (or a biomolecular complex) ProteinLens obtains an atomistic, energy-weighted graph description of the structure of the biomolecule, and subsequently provides a systematic analysis of allosteric signalling and communication across the structure using two computationally efficient methods: Markov Transients and bond-to-bond propensities. ProteinLens scores and ranks every bond and residue according to the speed and magnitude of the propagation of fluctuations emanating from any site of choice (e.g. the active site). The results are presented through statistical quantile scores visualised with interactive plots and adjustable 3D structure viewers, which can also be downloaded. ProteinLens thus allows the investigation of signalling in biomolecular structures of interest to aid the detection of allosteric sites and pathways. ProteinLens is implemented in Python/SQL and freely available to use at: www.proteinlens.io.

Repurposing of FDA-approved drugs against active site and potential allosteric drug-binding sites of COVID-19 main protease

The novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) still has serious negative effects on health, social life, and economics. Recently, vaccines from various companies have been urgently approved to control SARS-CoV-2 infections. However, any specific antiviral drug has not been confirmed so far for regular treatment. An important target is the main protease (M^pro ), which plays a major role in replication of the virus. In this study, Gaussian and residue network models are employed to reveal two distinct potential allosteric sites on M^pro that can be evaluated as drug targets besides the active site. Then, Food and Drug Administration (FDA)-approved drugs are docked to three distinct sites with flexible docking using AutoDock Vina to identify potential drug candidates. Fourteen best molecule hits for the active site of M^pro are determined. Six of these also exhibit high docking scores for the potential allosteric regions. Full-atom molecular dynamics simulations with MM-GBSA method indicate that compounds docked to active and potential allosteric sites form stable interactions with high binding free energy (∆G_bind ) values. ∆G_bind values reach -52.06 kcal/mol for the active site, -51.08 kcal/mol for the potential allosteric site 1, and - 42.93 kcal/mol for the potential allosteric site 2. Energy decomposition calculations per residue elucidate key binding residues stabilizing the ligands that can further serve to design pharmacophores. This systematic and efficient computational analysis successfully determines ivermectine, diosmin, and selinexor currently subjected to clinical trials, and further proposes bromocriptine, elbasvir as M^pro inhibitor candidates to be evaluated against SARS-CoV-2 infections.

Mechanism of allosteric activation of SIRT6 revealed by the action of rationally designed activators

The recent discovery of activator compounds binding to an allosteric site on the NAD⁺-dependent protein lysine deacetylase, sirtuin 6 (SIRT6) has attracted interest and presents a pharmaceutical target for aging-related and cancer diseases. However, the mechanism underlying allosteric activation of SIRT6 by the activator MDL-801 remains largely elusive because no major conformational changes are observed upon activator binding. By combining molecular dynamics simulations with biochemical and kinetic analyses of wild-type SIRT6 and its variant M136A, we show that conformational rotation of 2-methyl-4-fluoro-5-bromo substituent on the right phenyl ring (R-ring) of MDL-801, which uncovers previously unseen hydrophobic interactions, contributes to increased activating deacetylation activity of SIRT6. This hypothesis is further supported by the two newly synthesized MDL-801 derivatives through the removal of the 5-Br atom on the R-ring (MDL-801-D1) or the restraint of the rotation of the R-ring (MDL-801-D2). We further propose that the 5-Br atom serves as an allosteric driver that controls the ligand allosteric efficacy. Our study highlights the effect of allosteric enzyme catalytic activity by activator binding and provides a rational approach for enhancing deacetylation activity.

Molecular Dynamics-Based Allosteric Prediction Method to Design Key Residues in Threonine Dehydrogenase for Amino-Acid Production

Allosteric proteins are considered as one of the most critical targets to design cell factories via synthetic biology approaches. Here, we proposed a molecular dynamics-based allosteric prediction method (MBAP) to screen indirect-binding sites and potential mutations for protein re-engineering. Using this MBAP method, we have predicted new sites to relieve the allosteric regulation of threonine dehydrogenase (TD) by isoleucine. An obtained mutation P441L has been verified with the ability to significantly reduce the allosteric regulation of TD in vitro assays and with the fermentation application in vivo for amino-acid production. These findings have proved the MBAP method as an effective and efficient predicting tool to find new positions of the allosteric enzymes, thus opening a new path to constructing cell factories in synthetic biology.

Exploring the Allosteric Territory of Protein Function

While the pervasiveness of allostery in proteins is commonly accepted, we further show the generic nature of allosteric mechanisms by analyzing here transmembrane ion-channel viroporin 3a and RNA-dependent RNA polymerase (RdRp) from SARS-CoV-2 along with metabolic enzymes isocitrate dehydrogenase 1 (IDH1) and fumarate hydratase (FH) implicated in cancers. Using the previously developed structure-based statistical mechanical model of allostery (SBSMMA), we share our experience in analyzing the allosteric signaling, predicting latent allosteric sites, inducing and tuning targeted allosteric response, and exploring the allosteric effects of mutations. This, yet incomplete list of phenomenology, forms a complex and unique allosteric territory of protein function, which should be thoroughly explored. We propose a generic computational framework, which not only allows one to obtain a comprehensive allosteric control over proteins but also provides an opportunity to approach the fragment-based design of allosteric effectors and drug candidates. The advantages of allosteric drugs over traditional orthosteric compounds, complemented by the emerging role of the allosteric effects of mutations in the expansion of the cancer mutational landscape and in the increased mutability of viral proteins, leave no choice besides further extensive studies of allosteric mechanisms and their biomedical implications.

Computational Ways to Enhance Protein Inhibitor Design

Two new computational approaches are described to aid in the design of new peptide-based drugs by evaluating ensembles of protein structures from their dynamics and through the assessing of structures using empirical contact potential. These approaches build on the concept that conformational variability can aid in the binding process and, for disordered proteins, can even facilitate the binding of more diverse ligands. This latter consideration indicates that such a design process should be less restrictive so that multiple inhibitors might be effective. The example chosen here focuses on proteins/peptides that bind to hemagglutinin (HA) to block the large-scale conformational change for activation. Variability in the conformations is considered from sets of experimental structures, or as an alternative, from their simple computed dynamics; the set of designe peptides/small proteins from the David Baker lab designed to bind to hemagglutinin, is the large set considered and is assessed with the new empirical contact potentials.

Gas Chromatography-Mass Spectrometry Coupled with Multivariate Statistical Analysis to Identify the Alpha Glucosidase Inhibitors from Flesh of Salacca zalacca Fruits and Their Molecular Docking Studies

Fruit of salak (Salaaca zalacca) is traditionally used and commercialized as an antidiabetic agent. However, scientific evidence to prove this folk claim is quite lacking. Therefore, this research was aimed to evaluate the α-glucosidase inhibition activity of S. zalacca fruit and identify the bioactive compounds. The fruits were extracted by different ratios of ethanol and water (0, 20, 40, 60, 80, 100%, v/v) to get E0 (100% water), E20 (20% ethanol), E40 (40% ethanol), E60 (60% ethanol), E80 (80% ethanol), and E100 (100% ethanol) extracts. The extracts obtained were subjected to the α-glucosidase inhibitory assay. Gas chromatography-mass spectrometry- (GC-MS-) based metabolomics approach was used in profiling the bioactive metabolites present in the extracts. Orthogonal partial least square (OPLS) was used to correlate GC-MS data and α-glucosidase assay results to identify the possible chemical markers. All active compounds identified were subjected to molecular docking. The extracts from the S. zalacca fruit showed potent inhibition activity against α-glucosidase. The IC50 values from the α-glucosidase inhibitory assay ranged between 16 and 275 µg/ml. Overall, E60 displayed significantly higher α-glucosidase inhibition activity, while E0 showed the lowest α-glucosidase inhibition activity. Major compounds detected in S. zalacca fruits were sugars, fatty acids, and sterols, including myo-inositol, palmitic acid, stearic acid, and β-sitosterol. Moreover, the results obtained from molecular docking indicated that palmitic acid and β-sitosterol were close to the active side of the enzyme. Some of the residues that interacted include HID295, ASN259, LEU313, LYS125, PHE159, VAL216, PHE178, TYR72, TYR158, HIE315, ARG315, and PHE303. The bioassay result strongly suggests that E60 extract from S. zalacca fruits has potential α-glucosidase inhibitory activity. The hydrophobic compounds, including palmitic acid and β-sitosterol, were found to induce the α-glucosidase inhibition activity.

Discovery of allosteric binding sites by crystallographic fragment screening

Understanding allosteric regulation of proteins is fundamental to our study of protein structure and function. Moreover, allosteric binding pockets have become a major target of drug discovery efforts in recent years. However, even though the function of almost every protein can be influenced by allostery, it remains a challenge to discover, rationalise and validate putative allosteric binding pockets. This review examines how the discovery and analysis of putative allosteric binding sites have been influenced by the availability of centralised facilities for crystallographic fragment screening, along with newly developed computational methods for modelling low occupancy features. We discuss the experimental parameters required for success, and how new methods could influence the field in the future. Finally, we reflect on the general problem of how to translate these findings into actual ligand development programs.

Discovery of cryptic allosteric sites using reversed allosteric communication by a combined computational and experimental strategy

Allostery, which is one of the most direct and efficient methods to fine-tune protein functions, has gained increasing recognition in drug discovery. However, there are several challenges associated with the identification of allosteric sites, which is the fundamental cornerstone of drug design. Previous studies on allosteric site predictions have focused on communication signals propagating from the allosteric sites to the orthosteric sites. However, recent biochemical studies have revealed that allosteric coupling is bidirectional and that orthosteric perturbations can modulate allosteric sites through reversed allosteric communication. Here, we proposed a new framework for the prediction of allosteric sites based on reversed allosteric communication using a combination of computational and experimental strategies (molecular dynamics simulations, Markov state models, and site-directed mutagenesis). The desirable performance of our approach was demonstrated by predicting the known allosteric site of the small molecule MDL-801 in nicotinamide dinucleotide (NAD⁺)-dependent protein lysine deacetylase sirtuin 6 (Sirt6). A potential novel cryptic allosteric site located around the L116, R119, and S120 residues within the dynamic ensemble of Sirt6 was identified. The allosteric effect of the predicted site was further quantified and validated using both computational and experimental approaches. This study proposed a state-of-the-art computational pipeline for detecting allosteric sites based on reversed allosteric communication. This method enabled the identification of a previously uncharacterized potential cryptic allosteric site on Sirt6, which provides a starting point for allosteric drug design that can aid the identification of candidate pockets in other therapeutic targets.

Using reversed allosteric communication, we performed MD simulations, MSMs, and mutagenesis experiments, to discover allosteric sites. It reproduced the known allosteric site for MDL-801 on Sirt6 and uncovered a novel cryptic allosteric Pocket X.

AlloSigMA 2: paving the way to designing allosteric effectors and to exploring allosteric effects of mutations

The AlloSigMA 2 server provides an interactive platform for exploring the allosteric signaling caused by ligand binding and/or mutations, for analyzing the allosteric effects of mutations and for detecting potential cancer drivers and pathogenic nsSNPs. It can also be used for searching latent allosteric sites and for computationally designing allosteric effectors for these sites with required agonist/antagonist activity. The server is based on the implementation of the Structure-Based Statistical Mechanical Model of Allostery (SBSMMA), which allows one to evaluate the allosteric free energy as a result of the perturbation at per-residue resolution. The Allosteric Signaling Map (ASM) providing a comprehensive residue-by-residue allosteric control over the protein activity can be obtained for any structure of interest. The Allosteric Probing Map (APM), in turn, allows one to perform the fragment-based-like computational design experiment aimed at finding leads for potential allosteric effectors. The server can be instrumental in elucidating of allosteric mechanisms and actions of allosteric mutations, and in the efforts on design of new elements of allosteric control. The server is freely available at: http://allosigma.bii.a-star.edu.sg.

Disorder driven allosteric control of protein activity

Studies of protein allostery increasingly reveal an involvement of the back and forth order-disorder transitions in this mechanism of protein activity regulation. Here, we investigate the allosteric mechanisms mediated by structural disorder using the structure-based statistical mechanical model of allostery (SBSMMA) that we have previously developed. We show that SBSMMA accounts for the energetics and causality of allosteric communication underlying dimerization of the BirA biotin repressor, activation of the sortase A enzyme, and inhibition of the Rac1 GTPase. Using the SBSMMA, we also show that introducing structural order or disorder in various regions of esterases can originate tunable allosteric modulation of the catalytic triad. On the basis of obtained results, we propose that operating with the order-disorder continuum allows one to establish an allosteric control scale for achieving desired modulation of the protein activity.

•

Back and forth order-disorder transitions can induce allosteric signaling.

•

Allosteric signaling originated by order/disorder follow universal rules.

•

Allosteric control scale facilitates engineering of the protein activity regulation.

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Molecular dynamics (MD) simulation is a rigorous theoretical tool that when used efficiently could provide reliable answers to questions pertaining to the structure-function relationship of proteins. Data collated from protein dynamics can be translated into useful statistics that can be exploited to sieve thermodynamics and kinetics crucial for the elucidation of mechanisms responsible for the modulation of biological processes such as protein-ligand binding and protein-protein association. Continuous modernization of simulation tools enables accurate prediction and characterization of the aforementioned mechanisms and these qualities are highly beneficial for the expedition of drug development when effectively applied to structure-based drug design (SBDD). In this review, current all-atom MD simulation methods, with focus on enhanced sampling techniques, utilized to examine protein structure, dynamics, and functions are discussed. This review will pivot around computer calculations of protein-ligand and protein-protein systems with applications to SBDD. In addition, we will also be highlighting limitations faced by current simulation tools as well as the improvements that have been made to ameliorate their efficiency.

Zinc-mediated conformational preselection mechanism in the allosteric control of DNA binding to the zinc transcriptional regulator (ZitR)

The zinc transcriptional regulator (ZitR) functions as a metalloregulator that fine tunes transcriptional regulation through zinc-dependent DNA binding. However, the molecular mechanism of zinc-driven allosteric control of the DNA binding to ZitR remains elusive. Here, we performed enhanced sampling accelerated molecular dynamics simulations to figure out the mechanism, revealing the role of protein dynamics in the zinc-induced allosteric control of DNA binding to ZitR. The results suggest that zinc-free ZitR samples distinct conformational states, only a handful of which are compatible with DNA binding. Remarkably, zinc binding reduces the conformational plasticity of the DNA-binding domain of ZitR, promoting the population shift in the ZitR conformational ensemble towards the DNA binding-competent conformation. Further co-binding of DNA to the zinc–ZitR complex stabilizes this competent conformation. These findings suggest that ZitR–DNA interactions are allosterically regulated in a zinc-mediated conformational preselection manner, highlighting the importance of conformational dynamics in the regulation of transcription factor family.

Modelling and Recognition of Protein Contact Networks by Multiple Kernel Learning and Dissimilarity Representations

Multiple kernel learning is a paradigm which employs a properly constructed chain of kernel functions able to simultaneously analyse different data or different representations of the same data. In this paper, we propose an hybrid classification system based on a linear combination of multiple kernels defined over multiple dissimilarity spaces. The core of the training procedure is the joint optimisation of kernel weights and representatives selection in the dissimilarity spaces. This equips the system with a two-fold knowledge discovery phase: by analysing the weights, it is possible to check which representations are more suitable for solving the classification problem, whereas the pivotal patterns selected as representatives can give further insights on the modelled system, possibly with the help of field-experts. The proposed classification system is tested on real proteomic data in order to predict proteins' functional role starting from their folded structure: specifically, a set of eight representations are drawn from the graph-based protein folded description. The proposed multiple kernel-based system has also been benchmarked against a clustering-based classification system also able to exploit multiple dissimilarities simultaneously. Computational results show remarkable classification capabilities and the knowledge discovery analysis is in line with current biological knowledge, suggesting the reliability of the proposed system.

The role of cooperativity in a p53-miR34 dynamical mathematical model

The objective of this study is to evaluate the role of cooperativity, captured by the Hill coefficient, in a minimal mathematical model describing the interactions between p53 and miR-34a. The model equations are analyzed for negative, none and normal cooperativity using a specific version of bifurcation theory and they are solved numerically. Special attention is paid to the sign of so-called first Lyapunov value. Interpretations of the results are given, both according to dynamic theory and in biological terms. In terms of cell signaling, we propose the hypothesis that when the outgoing signal of a system spends a physiologically significant amount of time outside of its equilibrium state, then the value of that signal can be sampled at any point along the trajectory towards that equilibrium and indeed, at multiple points. Coupled with non-linear behavior, such as that caused by cooperativity, this feature can account for a complex and varied response, which p53 is known for. From dynamical point of view, we found that when cooperativity is negative, the system has only one stable equilibrium point. In the absence of cooperativity, there is a single unstable equilibrium point with a critical boundary of stability. In the case with normal cooperativity, the system can have one, two, or three steady states with both, bi-stability and bi-instability occurring.

Unveiling functional motions based on point mutations in biased signaling systems: A normal mode study on nerve growth factor bound to TrkA

Many receptors elicit signal transduction by activating multiple intracellular pathways. This transduction can be triggered by a non-specific ligand, which simultaneously activates all the signaling pathways of the receptors. However, the binding of one biased ligand preferentially trigger one pathway over another, in a process called biased signaling. The identification the functional motions related to each of these distinct pathways has a direct impact on the development of new effective and specific drugs. We show here how to detect specific functional motions by considering the case of the NGF/TrkA-Ig2 complex. NGF-mediated TrkA receptor activation is dependent on specific structural motions that trigger the neuronal growth, development, and survival of neurons in nervous system. The R221W mutation in the ngf gene impairs nociceptive signaling. We discuss how the large-scale structural effects of this mutation lead to the suppression of collective motions necessary to induce TrkA activation of nociceptive signaling. Our results suggest that subtle changes in the NGF interaction network due to the point mutation are sufficient to inhibit the motions of TrkA receptors putatively linked to nociception. The methodological approach presented in this article, based jointly on the normal mode analysis and the experimentally observed functional alterations due to point mutations provides an essential tool to reveal the structural changes and motions linked to the disease, which in turn could be necessary for a drug design study.