Discovery of hidden allosteric sites as novel targets for allosteric drug design

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.

Protein Allostery Study in Cells Using NMR Spectroscopy

Protein allostery is commonly observed in vitro. But how protein allostery behaves in cells is unknown. In this work, a protein monomer-dimer equilibrium system was built with the allosteric effect on the binding characterized using NMR spectroscopy through mutations away from the dimer interface. A chemical shift linear fitting method was developed that enabled us to accurately determine the dissociation constant. A total of 28 allosteric mutations were prepared and grouped to negative allosteric, nonallosteric, and positive allosteric modulators. ∼ 50% of mutations displayed the allosteric-state changes when moving from a buffered solution into cells. For example, there were no positive allosteric modulators in the buffered solution but eight in cells. The change in protein allostery is correlated with the interactions between the protein and the cellular environment. These interactions presumably drive the surrounding macromolecules in cells to transiently bind to the monomer and dimer mutational sites and change the free energies of the two species differently which generate new allosteric effects. These surrounding macromolecules create a new protein allostery pathway that is only present in cells.

Review on allosteric modulators of dopamine receptors so far

Background Contemporary research is predominantly directed towards allosteric modulators, a class of compounds designed to interact with specific sites distinct from the orthosteric site on G protein-coupled receptors. These allosteric modulators play a pivotal role in influencing diverse pharmacological effects, such as agonism/inverse agonism, efficacy modulation, and affinity modulation. One particularly intriguing aspect is the demonstrated capacity of allosteric modulation to enhance drug selectivity for therapeutic purposes, potentially leading to a reduction in serious side effects associated with traditional approaches. Allosteric ligands, a majority of which fall into the categories of negative allosteric modulators or positive allosteric modulators, exhibit the unique ability to either diminish or enhance the effects of endogenous ligands. Negative allosteric modulators weaken the response, while positive allosteric modulators intensify it. Additionally, silent allosteric modulators represent a distinct class that neither activates nor blocks the effects of endogenous ligands, adding complexity to the spectrum of allosteric modulation. In the broader context of central nervous system disorders, allosteric modulation takes center stage, particularly in the realm of dopamine receptors specifically, D1, D2, and D3 receptors. These receptors hold immense therapeutic potential for a range of conditions spanning neurodegenerative disorders to neurobehavioral and psychiatric disorders. The intricate modulation of dopamine receptors through allosteric mechanisms offers a nuanced and versatile approach to drug development. As research endeavors continue to unfold, the exploration of allosteric modulation stands as a promising frontier, holding the potential to reshape the landscape of drug discovery and therapeutic interventions in the field of neurology and psychiatry.

Identifying novel therapeutic inhibitors to target FMS-like tyrosine kinase-3 (FLT3) against acute myeloid leukemia: a molecular docking, molecular dynamics, and DFT study

Around 30% of acute myeloid leukemia (AML) patients have triggering mutations in Feline McDonough Sarcoma (FMS)-like tyrosine kinase 3 (FLT3), which has been suggested as a possible therapeutic candidate for AML therapy. Many tyrosine kinase inhibitors are available and have a wide variety of applications in the treatment of cancer by inhibiting subsequent steps of cell proliferation. Therefore, our study aims to identify effective antileukemic agents against FLT3 gene. Initially, well-known antileukemic drug candidates have been chosen to generate a structure-based pharmacophore model to assist the virtual screening of 217,77,093 compounds from the Zinc database. The final hits compounds were retrieved and evaluated by docking against the target protein, where the top four compounds have been selected for the analysis of ADMET. Based on the density functional theory (DFT), the geometry optimization, frontier molecular orbital (FMO), HOMO-LUMO, and global reactivity descriptor values have been evaluated that confirming a satisfactory profile and reactivity order for the selected candidates. In comparison to control compounds, the docking results revealed that the four compounds had substantial binding energies (-11.1 to -11.5 kcal/mol) with FLT3. The physicochemical and ADMET (adsorption, distribution, metabolism, excretion, toxicity) prediction results corresponded to the bioactive and safe candidates. Molecular dynamics (MD) confirmed the better binding affinity and stability compared to gilteritinib as a potential FLT3 inhibitor. In this study, a computational approach has been performed that found a better docking and dynamics score against target proteins, indicating potent and safe antileukemic agents, furthermore in-vivo and in-vitro investigations are recommended.Communicated by Ramaswamy H. Sarma.

Thrombin - A Molecular Dynamics Perspective

Thrombin is a crucial enzyme involved in blood coagulation, essential for maintaining circulatory system integrity and preventing excessive bleeding. However, thrombin is also implicated in pathological conditions such as thrombosis and cancer. Despite the application of various experimental techniques, including X-ray crystallography, NMR spectroscopy, and HDXMS, none of these methods can precisely detect thrombin's dynamics and conformational ensembles at high spatial and temporal resolution. Fortunately, molecular dynamics (MD) simulation, a computational technique that allows the investigation of molecular functions and dynamics in atomic detail, can be used to explore thrombin behavior. This review summarizes recent MD simulation studies on thrombin and its interactions with other biomolecules. Specifically, the 17 studies discussed here provide insights into thrombin's switch between 'slow' and 'fast' forms, active and inactive forms, the role of Na+ binding, the effects of light chain mutation, and thrombin's interactions with other biomolecules. The findings of these studies have significant implications for developing new therapies for thrombosis and cancer. By understanding thrombin's complex behavior, researchers can design more effective drugs and treatments that target thrombin.

A Structure-Based Allosteric Modulator Design Paradigm

Importance: Allosteric drugs bound to topologically distal allosteric sites hold a substantial promise in modulating therapeutic targets deemed undruggable at their orthosteric sites. Traditionally, allosteric modulator discovery has predominantly relied on serendipitous high-throughput screening. Nevertheless, the landscape has undergone a transformative shift due to recent advancements in our understanding of allosteric modulation mechanisms, coupled with a significant increase in the accessibility of allosteric structural data. These factors have extensively promoted the development of various computational methodologies, especially for machine-learning approaches, to guide the rational design of structure-based allosteric modulators. Highlights: We here presented a comprehensive structure-based allosteric modulator design paradigm encompassing 3 critical stages: drug target acquisition, allosteric binding site, and modulator discovery. The recent advances in computational methods in each stage are encapsulated. Furthermore, we delve into analyzing the successes and obstacles encountered in the rational design of allosteric modulators. Conclusion: The structure-based allosteric modulator design paradigm holds immense potential for the rational design of allosteric modulators. We hope that this review would heighten awareness of the use of structure-based computational methodologies in advancing the field of allosteric drug discovery.

Allosteric modulation of SHP2: Quest from known to unknown

Src homology-2 domain-containing protein tyrosine phosphatase-2 (SHP2) is a key regulatory factor in the cell cycle and its activating mutations play an important role in the development of various cancers, making it an important target for antitumor drugs. Due to the highly conserved amino acid sequence and positively charged nature of the active site of SHP2, it is difficult to discover inhibitors with high affinity for the catalytic site of SHP2 and sufficient cell permeability, making it considered an "undruggable" target. However, the discovery of allosteric regulation mechanisms provides new opportunities for transforming undruggable targets into druggable ones. Given the limitations of orthosteric inhibitors, SHP2 allosteric inhibitors have become a more selective and safer research direction. In this review, we elucidate the oncogenic mechanism of SHP2 and summarize the discovery methods of SHP2 allosteric inhibitors, providing new strategies for the design and improvement of SHP2 allosteric inhibitors.

Combining structural and coevolution information to unveil allosteric sites

Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.

Multisite Labeling of Proteins Using the Ligand-Directed Reactivity of Triggerable Michael Acceptors

Targeted modification of endogenous proteins without genetic manipulation of protein expression machinery has a range of applications from chemical biology to drug discovery. Despite being demonstrated to be effective in various applications, target-specific protein labeling using ligand-directed strategies is limited by stringent amino acid selectivity. Here, we present highly reactive ligand-directed triggerable Michael acceptors (LD-TMAcs) that feature rapid protein labeling. Unlike previous approaches, the unique reactivity of LD-TMAcs enables multiple modifications on a single target protein, effectively mapping the ligand binding site. This capability is attributed to the tunable reactivity of TMAcs that enable the labeling of several amino acid functionalities via a binding-induced increase in local concentration while remaining fully dormant in the absence of protein binding. We demonstrate the target selectivity of these molecules in cell lysates using carbonic anhydrase as the model protein. Furthermore, we demonstrate the utility of this method by selectively labeling membrane-bound carbonic anhydrase XII in live cells. We envision that the unique features of LD-TMAcs will find use in target identification, investigation of binding/allosteric sites, and studying membrane proteins.

Zafirlukast Is a Promising Scaffold for Selectively Inhibiting TNFR1 Signaling

Tumor necrosis factor (TNF) plays an important role in the pathogenesis of inflammatory and autoimmune diseases such as rheumatoid arthritis and Crohn's disease. The biological effects of TNF are mediated by binding to TNF receptors, TNF receptor 1 (TNFR1), or TNF receptor 2 (TNFR2), and this coupling makes TNFR1-specific inhibition by small-molecule therapies essential to avoid deleterious side effects. Recently, we engineered a time-resolved fluorescence resonance energy transfer biosensor for high-throughput screening of small molecules that modulate TNFR1 conformational states and identified zafirlukast as a compound that inhibits receptor activation, albeit at low potency. Here, we synthesized 16 analogues of zafirlukast and tested their potency and specificity for TNFR1 signaling. Using cell-based functional assays, we identified three analogues with significantly improved efficacy and potency, each of which induces a conformational change in the receptor (as measured by fluorescence resonance energy transfer (FRET) in cells). The best analogue decreased NF-κB activation by 2.2-fold, IκBα efficiency by 3.3-fold, and relative potency by two orders of magnitude. Importantly, we showed that the analogues do not block TNF binding to TNFR1 and that binding to the receptor's extracellular domain is strongly cooperative. Despite these improvements, the best candidate's maximum inhibition of NF-κB is only 63%, leaving room for further improvements to the zafirlukast scaffold to achieve full inhibition and prove its potential as a therapeutic lead. Interestingly, while we find that the analogues also bind to TNFR2 in vitro, they do not inhibit TNFR2 function in cells or cause any conformational changes upon binding. Thus, these lead compounds should also be used as reagents to study conformational-dependent activation of TNF receptors.

Glycomimetics for the inhibition and modulation of lectins

Carbohydrates are essential mediators of many processes in health and disease. They regulate self-/non-self- discrimination, are key elements of cellular communication, cancer, infection and inflammation, and determine protein folding, function and life-times. Moreover, they are integral to the cellular envelope for microorganisms and participate in biofilm formation. These diverse functions of carbohydrates are mediated by carbohydrate-binding proteins, lectins, and the more the knowledge about the biology of these proteins is advancing, the more interfering with carbohydrate recognition becomes a viable option for the development of novel therapeutics. In this respect, small molecules mimicking this recognition process become more and more available either as tools for fostering our basic understanding of glycobiology or as therapeutics. In this review, we outline the general design principles of glycomimetic inhibitors (Section 2). This section is then followed by highlighting three approaches to interfere with lectin function, i.e. with carbohydrate-derived glycomimetics (Section 3.1), novel glycomimetic scaffolds (Section 3.2) and allosteric modulators (Section 3.3). We summarize recent advances in design and application of glycomimetics for various classes of lectins of mammalian, viral and bacterial origin. Besides highlighting design principles in general, we showcase defined cases in which glycomimetics have been advanced to clinical trials or marketed. Additionally, emerging applications of glycomimetics for targeted protein degradation and targeted delivery purposes are reviewed in Section 4.

Structural insights into angiotensin receptor signaling modulation by balanced and biased agonists

The peptide hormone angiotensin II regulates blood pressure mainly through the type 1 angiotensin II receptor AT1 R and its downstream signaling proteins Gq and β-arrestin. AT1 R blockers, clinically used as antihypertensive drugs, inhibit both signaling pathways, whereas AT1 R β-arrestin-biased agonists have shown great potential for the treatment of acute heart failure. Here, we present a cryo-electron microscopy (cryo-EM) structure of the human AT1 R in complex with a balanced agonist, Sar1 -AngII, and Gq protein at 2.9 Å resolution. This structure, together with extensive functional assays and computational modeling, reveals the molecular mechanisms for AT1 R signaling modulation and suggests that a major hydrogen bond network (MHN) inside the receptor serves as a key regulator of AT1 R signal transduction from the ligand-binding pocket to both Gq and β-arrestin pathways. Specifically, we found that the MHN mutations N1113.35 A and N2947.45 A induce biased signaling to Gq and β-arrestin, respectively. These insights should facilitate AT1 R structure-based drug discovery for the treatment of cardiovascular diseases.

NMR applications to GPCR recognition by peptide ligands

Peptides form the largest group of ligands that modulate the activity of more than 120 different GPCRs. Among which linear disordered peptide ligands usually undergo significant conformational changes upon binding that is essential for receptor recognition and activation. Conformational selection and induced fit are the extreme mechanisms of coupled folding and binding that can be distinguished by analysis of binding pathways by methods that include NMR. However, the large size of GPCRs in membrane-mimetic environments limits NMR applications. In this review, we highlight advances in the field that can be adopted to address coupled folding and binding of peptide ligands to their cognate receptors.

How protein topology controls allosteric regulations

Allostery is an important regulatory mechanism of protein functions. Among allosteric proteins, certain protein structure types are more observed. However, how allosteric regulation depends on protein topology remains elusive. In this study, we extracted protein topology graphs at the fold level and found that known allosteric proteins mainly contain multiple domains or subunits and allosteric sites reside more often between two or more domains of the same fold type. Only a small fraction of fold-fold combinations are observed in allosteric proteins, and homo-fold-fold combinations dominate. These analyses imply that the locations of allosteric sites including cryptic ones depend on protein topology. We further developed TopoAlloSite, a novel method that uses the kernel support vector machine to predict the location of allosteric sites on the overall protein topology based on the subgraph-matching kernel. TopoAlloSite successfully predicted known cryptic allosteric sites in several allosteric proteins like phosphopantothenoylcysteine synthetase, spermidine synthase, and sirtuin 6, demonstrating its power in identifying cryptic allosteric sites without performing long molecular dynamics simulations or large-scale experimental screening. Our study demonstrates that protein topology largely determines how its function can be allosterically regulated, which can be used to find new druggable targets and locate potential binding sites for rational allosteric drug design.

Mapping cryptic binding sites of drug targets to overcome drug resistance

The emergence of drug resistance is a primary obstacle for successful chemotherapy. Drugs that target cryptic binding sites (CBSs) represent a novel strategy for overcoming drug resistance. In this short communication, we explain and discuss how the discovery of CBSs and their inhibitors can overcome drug resistance.

Comprehensive Survey of Consensus Docking for High-Throughput Virtual Screening

The rapid advances of 3D techniques for the structural determination of proteins and the development of numerous computational methods and strategies have led to identifying highly active compounds in computer drug design. Molecular docking is a method widely used in high-throughput virtual screening campaigns to filter potential ligands targeted to proteins. A great variety of docking programs are currently available, which differ in the algorithms and approaches used to predict the binding mode and the affinity of the ligand. All programs heavily rely on scoring functions to accurately predict ligand binding affinity, and despite differences in performance, none of these docking programs is preferable to the others. To overcome this problem, consensus scoring methods improve the outcome of virtual screening by averaging the rank or score of individual molecules obtained from different docking programs. The successful application of consensus docking in high-throughput virtual screening highlights the need to optimize the predictive power of molecular docking methods.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Evolving Experimental Techniques for Structure-Based Drug Design

Structure-based drug design (SBDD) is a prominent method in rational drug development and has traditionally benefitted from the atomic models of protein targets obtained using X-ray crystallography at cryogenic temperatures. In this perspective, we highlight recent advances in the development of structural techniques that are capable of probing dynamic information about protein targets. First, we discuss advances in the field of X-ray crystallography including serial room-temperature crystallography as a method for obtaining high-resolution conformational dynamics of protein-inhibitor complexes. Next, we look at cryogenic electron microscopy (cryoEM), another high-resolution technique that has recently been used to study proteins and protein complexes that are too difficult to crystallize. Finally, we present small-angle X-ray scattering (SAXS) as a potential high-throughput screening tool to identify inhibitors that target protein complexes and protein oligomerization.

Structural-Guided Identification of Small Molecule Inhibitor of UHRF1 Methyltransferase Activity

Ubiquitin-like containing plant homeodomain Ring Finger 1 (UHRF1) protein is recognized as a cell-cycle-regulated multidomain protein. UHRF1 importantly manifests the maintenance of DNA methylation mediated by the interaction between its SRA (SET and RING associated) domain and DNA methyltransferase-1 (DNMT1)-like epigenetic modulators. However, overexpression of UHRF1 epigenetically responds to the aberrant global methylation and promotes tumorigenesis. To date, no potential molecular inhibitor has been studied against the SRA domain. Therefore, this study focused on identifying the active natural drug-like candidates against the SRA domain. A comprehensive set of in silico approaches including molecular docking, molecular dynamics (MD) simulation, and toxicity analysis was performed to identify potential candidates. A dataset of 709 natural compounds was screened through molecular docking where chicoric acid and nystose have been found showing higher binding affinities to the SRA domain. The MD simulations also showed the protein ligand interaction stability of and in silico toxicity analysis has also showed chicoric acid as a safe and nontoxic drug. In addition, chicoric acid possessed a longer interaction time and higher LD50 of 5000 mg/kg. Moreover, the global methylation level (%5 mC) has been assessed after chicoric acid treatment was in the colorectal cancer cell line (HCT116) at different doses. The result showed that 7.5 µM chicoric acid treatment reduced methylation levels significantly. Thus, the study found chicoric acid can become a possible epidrug-like inhibitor against the SRA domain of UHRF1 protein.

Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study

Small ubiquitin-related modifier (SUMO)-specific protease 1 (SENP1) is a cysteine protease that catalyzes the cleavage of the C-terminus of SUMO1 for the processing of SUMO precursors and deSUMOylation of target proteins. SENP1 is considered to be a promising target for the treatment of hepatocellular carcinoma (HCC) and prostate cancer. SENP1 Gln597 is located at the unstructured loop connecting the helices α4 to α5. The Q597A mutation of SENP1 allosterically disrupts the hydrolytic reaction of SUMO1 through an unknown mechanism. Here, extensive multiple replicates of microsecond molecular dynamics (MD) simulations, coupled with principal component analysis, dynamic cross-correlation analysis, community network analysis, and binding free energy calculations, were performed to elucidate the detailed mechanism. Our MD simulations showed that the Q597A mutation induced marked dynamic conformational changes in SENP1, especially in the unstructured loop connecting the helices α4 to α5 which the mutation site occupies. Moreover, the Q597A mutation caused conformational changes to catalytic Cys603 and His533 at the active site, which might impair the catalytic activity of SENP1 in processing SUMO1. Moreover, binding free energy calculations revealed that the Q597A mutation had a minor effect on the binding affinity of SUMO1 to SENP1. Together, these results may broaden our understanding of the allosteric modulation of the SENP1−SUMO1 complex.

Unmasking allosteric binding sites: Novel targets for GPCR drug discovery

INTRODUCTION

Unexpected non-apparent and hidden allosteric binding sites are non-classical and non-apparent allosteric centers in 3-D X-ray protein structures until orthosteric or allosteric ligands bind to them. The orthosteric center of one protomer that modulates binding centers of the other protomers within an oligomer is also an unexpected allosteric site. Furthermore, another partner protein can also produce these effects, acting as an unexpected allosteric modulator.

AREAS COVERED

This review summarizes both classical and non-classical allosterism. The authors focus on G protein-coupled receptor (GPCR) oligomers as a paradigm of allosteric molecules. Moreover, they show several examples of unexpected allosteric sites such as hidden allosteric sites in a protomer that appear after the interaction with other molecules and the allosterism exerted between orthosteric sites within GPCR oligomer, emphasizing on the allosteric modulations that can occur between binding sites.

EXPERT OPINION

The study of these new non-classical allosteric sites will expand the diversity of allosteric control on the function of orthosteric sites within proteins, whether GPCRs or other receptors, enzymes or transporters. Moreover, the design of new drugs targeting these hidden allosteric sites or already known orthosteric sites acting as allosteric sites in protein homo- or hetero-oligomers will increase the therapeutic potential of allosterism.

Understanding the P-Loop Conformation in the Determination of Inhibitor Selectivity Toward the Hepatocellular Carcinoma-Associated Dark Kinase STK17B

As a member of the death-associated protein kinase family of serine/threonine kinases, the STK17B has been associated with diverse diseases such as hepatocellular carcinoma. However, the conformational dynamics of the phosphate-binding loop (P-loop) in the determination of inhibitor selectivity profile to the STK17B are less understood. Here, a multi-microsecond length molecular dynamics (MD) simulation of STK17B in the three different states (ligand-free, ADP-bound, and ligand-bound states) was carried out to uncover the conformational plasticity of the P-loop. Together with the analyses of principal component analysis, cross-correlation and generalized correlation motions, secondary structural analysis, and community network analysis, the conformational dynamics of the P-loop in the different states were revealed, in which the P-loop flipped into the ADP-binding site upon the inhibitor binding and interacted with the inhibitor and the C-lobe, strengthened the communication between the N- and C-lobes. These resulting interactions contributed to inhibitor selectivity profile to the STK17B. Our results may advance our understanding of kinase inhibitor selectivity and offer possible implications for the design of highly selective inhibitors for other protein kinases.

Computational Identification of Druggable Bioactive Compounds from Catharanthus roseus and Avicennia marina against Colorectal Cancer by Targeting Thymidylate Synthase

Colorectal cancer (CRC) is the second most common cause of death worldwide, affecting approximately 1.9 million individuals in 2020. Therapeutics of the disease are not yet available and discovering a novel anticancer drug candidate against the disease is an urgent need. Thymidylate synthase (TS) is an important enzyme and prime precursor for DNA biosynthesis that catalyzes the methylation of deoxyuridine monophosphate (dUMP) to deoxythymidine monophosphate (dTMP) that has emerged as a novel drug target against the disease. Elevated expression of TS in proliferating cells promotes oncogenesis as well as CRC. Therefore, this study aimed to identify potential natural anticancer agents that can inhibit the activity of the TS protein, subsequently blocking the progression of colorectal cancer. Initially, molecular docking was implied on 63 natural compounds identified from Catharanthus roseus and Avicennia marina to evaluate their binding affinity to the desired protein. Subsequently, molecular dynamics (MD) simulation, ADME (Absorption, Distribution, Metabolism, and Excretion), toxicity, and quantum chemical-based DFT (density-functional theory) approaches were applied to evaluate the efficacy of the selected compounds. Molecular docking analysis initially identified four compounds (PubChem CID: 5281349, CID: 102004710, CID: 11969465, CID: 198912) that have better binding affinity to the target protein. The ADME and toxicity properties indicated good pharmacokinetics (PK) and toxicity ability of the selected compounds. Additionally, the quantum chemical calculation of the selected molecules found low chemical reactivity indicating the bioactivity of the drug candidate. The global descriptor and HOMO-LUMO energy gap values indicated a satisfactory and remarkable profile of the selected molecules. Furthermore, MD simulations of the compounds identified better binding stability of the compounds to the desired protein. To sum up, the phytoconstituents from two plants showed better anticancer activity against TS protein that can be further developed as an anti-CRC drug.

Investigation of Unprecedented Sites and Proposition of New Ligands for Programmed Cell Death Protein I through Molecular Dynamics with Probes and Virtual Screening

Cancer immunotherapy has attracted increasing attention over the last few years. Programmed cell death protein 1 (PD-1) promotes self-tolerance and inhibits immune responses by modulating the T-cell function. The interaction between PD-1 and programmed cell death ligand-1 (PD-L1) leads to immune exhaustion, protecting cancer cells from destruction. Here, we computationally designed a novel ligand named 1508 that binds to an unprecedented PD-1 cavity identified by MixMD and defined by amino acid residues Lys78 to Val97. We showed through a set of MD simulations totaling 12.5 μs that ligand 1508 establishes frequent cation-π and hydrogen bonding interactions with amino acid residues Lys78 and Arg86, respectively, and stabilizes the PD-1 C'D loop in a conformation that does not favor PD-1-PD-L1 complex formation. This study highlights the power of MixMD in exposing new cavities prone to protein-protein complex inhibition and establishes the basis for the design of new molecules that target the PD-1 C'D cavity as an alternative for exploring the modulation of the PD-1-PD-L1 complex in cancer therapy.

Probing GPCR Dimerization Using Peptides

G protein-coupled receptors (GPCRs) are the largest class of membrane proteins and the most common and extensively studied pharmacological target. Numerous studies over the last decade have confirmed that GPCRs do not only exist and function in their monomeric form but in fact, have the ability to form dimers or higher order oligomers with other GPCRs, as well as other classes of receptors. GPCR oligomers have become increasingly attractive to investigate as they have the ability to modulate the pharmacological responses of the receptors which in turn, could have important functional roles in diseases, such as cancer and several neurological & neuropsychiatric disorders. Despite the growing evidence in the field of GPCR oligomerisation, the lack of structural information, as well as targeting the 'undruggable' protein-protein interactions (PPIs) involved in these complexes, has presented difficulties. Outside the field of GPCRs, targeting PPIs has been widely studied, with a variety of techniques being investigated; from small-molecule inhibitors to disrupting peptides. In this review, we will demonstrate several physiologically relevant GPCR dimers and discuss an array of strategies and techniques that can be employed when targeting these complexes, as well as provide ideas for future development.

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.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Strategies for targeting undruggable targets

INTRODUCTION

Undruggable targets refer to clinically meaningful therapeutic targets that are 'difficult to drug' or 'yet to be drugged' via traditional approaches. Featuring characteristics of lacking defined ligand-binding pockets, non-catalytic protein-protein interaction functional modes and less-investigated 3D structures, these undruggable targets have been targeted with novel therapeutic entities developed with the progress of unconventional drug discovery approaches, such as targeted degradation molecules and display technologies.

AREA COVERED

This review first presents the concept of 'undruggable' exemplified by RAS and other targets. Next, detailed strategies are illustrated in two aspects: innovation of therapeutic entities and development of unconventional drug discovery technologies. Finally, case studies covering typical undruggable targets (Bcl-2, p53, and RAS) are depicted to further demonstrate the feasibility of the strategies and entities above.

EXPERT OPINION

Targeting the undruggable expands the scope of therapeutically reachable targets. Consequently, it represents the drug discovery frontier. Biomedical studies are capable of dissecting disease mechanisms, thus broadening the list of undruggable targets. Encouraged by the recent approval of the KRAS inhibitor Sotorasib, we believe that merging multiple discovery approaches and exploiting various novel therapeutic entities would pave the way for dealing with more 'undruggable' targets in the future.

Mechanistic Insights into the Allosteric Inhibition of Androgen Receptors by Binding Function 3 Antagonists from an Integrated Molecular Modeling Study

An androgen receptor (AR) is an intensively studied treatment target for castration-resistant prostate cancer that is irresponsive to conventional antiandrogen therapeutics. Binding function 3 (BF3) inhibitors with alternative modes of action have emerged as a promising approach to overcoming antiandrogen resistance. However, how these BF3 inhibitors modulate AR function remains elusive, hindering the development of BF3-targeting agents. Here, we performed an integrated computational study to interrogate the binding mechanism of several known BF3 inhibitors with ARs. Our results show that the inhibitory effect of the BF3 antagonists arises from their allosteric modulation of the activation function (AF2) site, which alters the dynamic coupling between the BF3 and AF2 sites as well as the AF2-coactivator (SRC2-3) interaction. Moreover, the per-residue binding energy analyses reveal the "anchor" role of the linker connecting the phenyl ring and benzimidazole/indole in these BF3 inhibitors. Furthermore, the allosteric driver-interacting residues are found to include both "positive", e.g., Phe673 and Asn833, and "negative" ones, e.g., Phe826, and the differential interactions with these residues provide an explanation why stronger binding does not necessarily result in higher inhibitory activities. Finally, our allosteric communication pathway analyses delineate how the allosteric signals triggered by BF3 binding are propagated to the AF2 pocket through multiple short- and/or long-ranged transmission pathways. Collectively, our combined computational study provides a comprehensive structural mechanism underlying how the selected set of BF3 inhibitors modulate AR function, which will help guide future development of BF3 antagonists.

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.

Inhibition of Nonfunctional Ras

Intuitively, functional states should be targeted; not nonfunctional ones. So why could drugging the inactive K-Ras4BG12Cwork-but drugging the inactive kinase will likely not? The reason is the distinct oncogenic mechanisms. Kinase driver mutations work by stabilizing the active state and/or destabilizing the inactive state. Either way, oncogenic kinases are mostly in the active state. Ras driver mutations work by quelling its deactivation mechanisms, GTP hydrolysis, and nucleotide exchange. Covalent inhibitors that bind to the inactive GDP-bound K-Ras4BG12C conformation can thus work. By contrast, in kinases, allosteric inhibitors work by altering the active-site conformation to favor orthosteric drugs. From the translational standpoint this distinction is vital: it expedites effective pharmaceutical development and extends the drug classification based on the mechanism of action. Collectively, here we postulate that drug action relates to blocking the mechanism of activation, not to whether the protein is in the active or inactive state.

Mechanistic insights into the effect of phosphorylation on Ras conformational dynamics and its interactions with cell signaling proteins

Ras undergoes interconversion between the active GTP-bound state and the inactive GDP-bound state. This GTPase cycle, which controls the activities of Ras, is accelerated by Ras GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors (SOS). Oncogenic Ras mutations could affect the GTPase cycle and impair Ras functions. Additionally, Src-induced K-Ras Y32/64 dual phosphorylation has been reported to disrupt GTPase cycle and hinder Ras downstream signaling. However, the underlying mechanisms remain unclear. To address this, we performed molecular dynamics simulations (~30 μs in total) on unphosphorylated and phosphorylated K-Ras4B in GTP- and GDP-bound states, and on their complexes with GTPase cycle regulators (GAP and SOS) and the effector protein Raf. We found that K-Ras4B dual phosphorylation mainly alters the conformation at the nucleotide binding site and creates disorder at the catalytic site, resulting in the enlargement of GDP binding pocket and the retard of Ras-GTP intrinsic hydrolysis. We observed phosphorylation-induced shift in the distribution of Ras-GTP inactive-active sub-states and recognized potential druggable pockets in the phosphorylated Ras-GTP. Moreover, decreased catalytic competence or signal delivery abilities due to reduced binding affinities and/or distorted catalytic conformations of GAP, SOS and Raf were observed. In addition, the allosteric pathway from Ras/Raf interface to the distal Raf L4 loop was compromised by Ras phosphorylation. These results reveal the mechanisms by which phosphorylation influences the intrinsic or GAP/SOS catalyzed transformations between GTP- and GDP-bound states of Ras and its signal transduction to Raf. Our findings project Ras phosphorylation as a target for cancer drug discovery.

Optogenetic Analysis of Allosteric Control in Protein Tyrosine Phosphatases

Allosteric regulation enables dynamic adjustments to protein function that permit tight control over cellular biochemistry. Discrepancies in the allosteric systems of related proteins can thus reveal important differences in their susceptibilities to influential stimuli (e.g., allosteric ligands, mutations, or post-translational modifications). This study uses an optogenetic actuator as a tool to compare the allosteric systems of two structurally related regulatory proteins: protein tyrosine phosphatase 1B (PTP1B) and T-cell protein tyrosine phosphatase (TCPTP). It begins with an interesting observation: The fusion of a protein light switch to the allosterically influential α7 helix of PTP1B permits optical modulation of its catalytic activity, but a similar fusion to TCPTP does not. A subsequent analysis of different PTP chimeras shows that replacing regions of TCPTP with homologous regions from PTP1B can enhance photocontrol; as TCPTP becomes more "PTP1B-like", its photosensitivity increases. Interestingly, the structural changes required for photocontrol also enhance the sensitivity of TCPTP to other allosteric inputs, notably, an allosteric inhibitor and a newly reported activating mutation. Our findings indicate that the allosteric functionality of the α7 helix of PTP1B is not conserved across the PTP family and highlight residues necessary to transfer this functionality to other PTPs. More broadly, our results suggest that simple gene fusion events can strengthen allosteric communication within individual protein domains and describe an intriguing application for optogenetic actuators as structural probes-a sort of physically disruptive "ratchet"-for studying protein allostery.

Molecular Dynamics Simulation Techniques as Tools in Drug Discovery and Pharmacology: A Focus on Allosteric Drugs

Allosteric drugs are ligands that when bound to an allosteric site modify the conformational state of the pharmacological target, leading then to a modification of functional response upon binding of the endogenous ligand. Pharmacological targets are defined as biological entities, to which a ligand/drug binds and leads to a functional effect. Pharmacological targets can be proteins or nucleic acids. Computational approaches such as molecular dynamics (MD) sped up discovery and identification of allosteric binding sites and allosteric ligands. Classical all-atom and hybrid classical/quantum MD simulations can be generalized as simulation techniques aimed at analysis of atoms and molecular motion. Main limitations of MD simulations are related to high computational costs, that in turn limit the conformational sampling of biological systems. Indeed, other techniques have been developed to overcome limitations of MD, such as enhanced sampling MD simulations. In this chapter, classical MD and enhanced sampling MD simulations will be described, along with their application to drug discovery, with a focus on allosteric drugs.

Protein folding stability and binding interactions through the lens of evolution: a dynamical perspective

While the function of a protein depends heavily on its ability to fold into a correct 3D structure, billions of years of evolution have tailored proteins from highly stable objects to flexible molecules as they adapted to environmental changes. Nature maintains the fine balance of protein folding and stability while still evolving towards new function through generations of fine-tuning necessary interactions with other proteins and small molecules. Here we focus on recent computational and experimental studies that shed light onto how evolution molds protein folding and the functional landscape from a conformational dynamics' perspective. Particularly, we explore the importance of dynamic allostery throughout protein evolution and discuss how the protein anisotropic network can give rise to allosteric and epistatic interactions.

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.

Spatiotemporal identification of druggable binding sites using deep learning

Identification of novel protein binding sites expands druggable genome and opens new opportunities for drug discovery. Generally, presence or absence of a binding site depends on the three-dimensional conformation of a protein, making binding site identification resemble the object detection problem in computer vision. Here we introduce a computational approach for the large-scale detection of protein binding sites, that considers protein conformations as 3D-images, binding sites as objects on these images to detect, and conformational ensembles of proteins as 3D-videos to analyze. BiteNet is suitable for spatiotemporal detection of hard-to-spot allosteric binding sites, as we showed for conformation-specific binding site of the epidermal growth factor receptor, oligomer-specific binding site of the ion channel, and binding site in G protein-coupled receptor. BiteNet outperforms state-of-the-art methods both in terms of accuracy and speed, taking about 1.5 minutes to analyze 1000 conformations of a protein with ~2000 atoms.

Kozlovskii and Popov present BiteNet, a new computational method utilizing deep learning principles for rapid detection of binding sites. BiteNet considers proteins as 3D images, enabling rapid detection of allosteric sites from either static protein structures or its dynamic ensembles.

Backbone resonance assignment and dynamics of 110 kDa hexameric inorganic pyrophosphatase from Mycobacterium tuberculosis

Family I soluble inorganic pyrophosphatases (PPases; EC 3.6.1.1) are enzymes essential for all organisms. They hydrolyze inorganic pyrophosphate, thus providing the driving force for numerous biosynthetic reactions. Soluble PPases retain enzymatic activity only in multimeric forms. PPases from various organisms are extensively studied by X-ray crystallography but until now there was no information on their structure and dynamics in solution. Hexameric 110 kDa (6 × 18.3 kDa) PPase from Mycobacterium tuberculosis (Mt-PPase) is a promising target for the rational design of potential anti-tuberculosis agents. In order to use NMR techniques in functional studies of Mt-PPase and rational design of the inhibitors for this enzyme, it is necessary to have information on the backbone ¹H, ¹³C and ¹⁵N resonance assignments. Samples of Mt-PPase enriched with 99% of ¹³C and ¹⁵N isotopes, and 95% of ²H were obtained using recombinant protein expression in an isotopically-labeled medium and effective heat-shock protocol for the deuterium-to-hydrogen exchange of the amide groups. Backbone resonance assignment was achieved for more than 95% of the residues. It was found that the secondary structure of Mt-PPase in solution corresponds well to the crystal structure of this protein. Protein backbone dynamics were studied using ¹⁵N NMR relaxation experiments. Determined resonance assignments and dynamic properties provide the basis for the subsequent structure-based design of novel inhibitors of Mt-PPase-potential anti-tuberculosis drugs.

Recent progress on cheminformatics approaches to epigenetic drug discovery

The ability of epigenetic markers to affect genome function has enabled transformative changes in drug discovery, especially in cancer and other emerging therapeutic areas. Concordant with the introduction of the term 'epi-informatics', the size of the epigenetically relevant chemical space has grown substantially and so did the number of applications of cheminformatic methods to epigenetics. Recent progress in epi-informatics has improved our understanding of the structure-epigenetic activity relationships and boosted the development of models predicting novel epigenetic agents. Herein, we review the advances in computational approaches to drug discovery of small molecules with epigenetic modulation profiles, summarize the current chemogenomics data available for epigenetic targets, and provide a perspective on the greater utility of biomedical knowledge mining as a means to advance the epigenetic drug discovery.

Bcr-Abl Allosteric Inhibitors: Where We Are and Where We Are Going to

The fusion oncoprotein Bcr-Abl is an aberrant tyrosine kinase responsible for chronic myeloid leukemia and acute lymphoblastic leukemia. The auto-inhibition regulatory module observed in the progenitor kinase c-Abl is lost in the aberrant Bcr-Abl, because of the lack of the N-myristoylated cap able to bind the myristoyl binding pocket also conserved in the Bcr-Abl kinase domain. A way to overcome the occurrence of resistance phenomena frequently observed for Bcr-Abl orthosteric drugs is the rational design of allosteric ligands approaching the so-called myristoyl binding pocket. The discovery of these allosteric inhibitors although very difficult and extremely challenging, represents a valuable option to minimize drug resistance, mostly due to the occurrence of mutations more frequently affecting orthosteric pockets, and to enhance target selectivity with lower off-target effects. In this perspective, we will elucidate at a molecular level the structural bases behind the Bcr-Abl allosteric control and will show how artificial intelligence can be effective to drive the automated de novo design towards off-patent regions of the chemical space.

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.

Driving Protein Conformational Cycles in Physiology and Disease: "Frustrated" Amino Acid Interaction Networks Define Dynamic Energy Landscapes: Amino Acid Interaction Networks Change Progressively Along Alpha Tryptophan Synthase's Catalytic Cycle

A general framework by which dynamic interactions within a protein will promote the necessary series of structural changes, or "conformational cycle," required for function is proposed. It is suggested that the free-energy landscape of a protein is biased toward this conformational cycle. Fluctuations into higher energy, although thermally accessible, conformations drive the conformational cycle forward. The amino acid interaction network is defined as those intraprotein interactions that contribute most to the free-energy landscape. Some network connections are consistent in every structural state, while others periodically change their interaction strength according to the conformational cycle. It is reviewed here that structural transitions change these periodic network connections, which then predisposes the protein toward the next set of network changes, and hence the next structural change. These concepts are illustrated by recent work on tryptophan synthase. Disruption of these dynamic connections may lead to aberrant protein function and disease states.

Fractional-Order Susceptible-Infected Model: Definition and Applications to the Study of COVID-19 Main Protease

We propose a model for the transmission of perturbations across the amino acids of a protein represented as an interaction network. The dynamics consists of a Susceptible-Infected (SI) model based on the Caputo fractional-order derivative. We find an upper bound to the analytical solution of this model which represents the worse-case scenario on the propagation of perturbations across a protein residue network. This upper bound is expressed in terms of Mittag-Leffler functions of the adjacency matrix of the network of inter-amino acids interactions. We then apply this model to the analysis of the propagation of perturbations produced by inhibitors of the main protease of SARS CoV-2. We find that the perturbations produced by strong inhibitors of the protease are propagated far away from the binding site, confirming the long-range nature of intra-protein communication. On the contrary, the weakest inhibitors only transmit their perturbations across a close environment around the binding site. These findings may help to the design of drug candidates against this new coronavirus.

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.