Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often “cryptic” without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein–ligand systems.     

Studies on fragment-based design of allosteric inhibitors of human factor XIa

Human factor XIa (hFXIa) has emerged as an attractive target for development of new anticoagulants that promise higher level of safety. Different strategies have been adopted so far for the design of anti-hFXIa molecules including competitive and non-competitive inhibition. Of these, allosteric dysfunction of hFXIa’s active site is especially promising because of the possibility of controlled reduction in activity that may offer a route to safer anticoagulants. In this work, we assess fragment-based design approach to realize a group of novel allosteric hFXIa inhibitors. Starting with our earlier discovery that sulfated quinazolinone (QAO) bind in the heparin-binding site of hFXIa, we developed a group of two dozen dimeric sulfated QAOs with intervening linkers that displayed a progressive variation in inhibition potency. In direct opposition to the traditional wisdom, increasing linker flexibility led to higher potency, which could be explained by computational studies. Sulfated QAO 19S was identified as the most potent and selective inhibitor of hFXIa. Enzyme inhibition studies revealed that 19S utilizes a non-competitive mechanism of action, which was supported by fluorescence studies showing a classic sigmoidal binding profile. Studies with selected mutants of hFXIa indicated that sulfated QAOs bind in heparin-binding site of the catalytic domain of hFXIa. Overall, the approach of fragment-based design offers considerable promise for designing heparin-binding site-directed allosteric inhibitors of hFXIa.     

Neutralizing Antibodies Against Allosteric Proteins: Insights From a Bacterial Adhesin

Allosteric proteins transition between ‘inactive’ and ‘active’ states. In general, such proteins assume distinct conformational states at the level of secondary, tertiary and/or quaternary structure. Different conformers of an allosteric protein can be antigenically dissimilar and induce antibodies with a highly distinctive specificities and neutralizing functional effects. Here we summarize studies on various functional types of monoclonal antibodies obtained against different allosteric conformers of the mannose-specific bacterial adhesin FimH – the most common cell attachment protein of Escherichia coli and other enterobacterial pathogens. Included are types of antibodies that activate the FimH function via interaction with ligand-induced binding sites or by wedging between domains as well as antibodies that inhibit FimH through orthosteric, parasteric, or novel dynasteric mechanisms. Understanding the molecular mechanism of antibody action against allosteric proteins provides insights on how to design antibodies with a desired functional effect, including those with neutralizing activity against bacterial and viral cell attachment proteins.     

Design and engineering of allosteric communications in proteins

Allostery in proteins plays an important role in regulating protein activities and influencing many biological processes such as gene expression, enzyme catalysis, and cell signaling. The process of allostery takes place when a signal detected at a site on a protein is transmitted via a mechanical pathway to a functional site and, thus, influences its activity. The pathway of allosteric communication consists of amino acids that form a network with covalent and non-covalent bonds. By mutating residues in this allosteric network, protein engineers have successfully established novel allosteric pathways to achieve desired properties in the target protein. In this review, we highlight the most recent and state-of-the-art techniques for allosteric communication engineering. We also discuss the challenges that need to be overcome and future directions for engineering protein allostery.     

Triggering Closure of a Sialic Acid TRAP Transporter Substrate Binding Protein through Binding of Natural or Artificial Substrates

The pathogens Vibrio cholerae and Haemophilus influenzae use tripartite ATP-independent periplasmic transporters (TRAPs) to scavenge sialic acid from host tissues. They use it as a nutrient or to evade the innate immune system by sialylating surface lipopolysaccharides. An essential component of TRAP transporters is a periplasmic substrate binding protein (SBP). Without substrate, the SBP has been proposed to rest in an open-state, which is not recognised by the transporter. Substrate binding induces a conformational change of the SBP and it is thought that this closed state is recognised by the transporter, triggering substrate translocation. Here we use real time single molecule FRET experiments and crystallography to investigate the open- to closed-state transition of VcSiaP, the SBP of the sialic acid TRAP transporter from V. cholerae. We show that the conformational switching of VcSiaP is strictly substrate induced, confirming an important aspect of the proposed transport mechanism. Two new crystal structures of VcSiaP provide insights into the closing mechanism. While the first structure contains the natural ligand, sialic acid, the second structure contains an artificial peptide in the sialic acid binding site. Together, the two structures suggest that the ligand itself stabilises the closed state and that SBP closure is triggered by physically bridging the gap between the two lobes of the SBP. Finally, we demonstrate that the affinity for the artificial peptide substrate can be substantially increased by varying its amino acid sequence and by this, serve as a starting point for the development of peptide-based inhibitors of TRAP transporters.     

Molecular dynamics and simulation analysis against superoxide dismutase (SOD) target of Micrococcus luteus with secondary metabolites from Bacillus licheniformis recognized by genome mining approach

Micrococcus luteus, also known as M. luteus, is a bacterium that inhabits mucous membranes, human skin, and various environmental sources. It is commonly linked to infections, especially among individuals who have compromised immune systems. M. luteus is capable of synthesizing the enzyme superoxide dismutase (SOD) as a component of its protective response to reactive oxygen species (ROS). This enzyme serves as a promising target for drug development in various diseases. The current study utilized a subtractive genomics approach to identify potential therapeutic targets from M. luteus. Additionally, genome mining was employed to identify and characterize the biosynthetic gene clusters (BGCs) responsible for the production of secondary metabolites in Bacillus licheniformis (B. licheniformis), a bacterium known for its production of therapeutically relevant secondary metabolites. Subtractive genomics resulted in identification of important extracellular protein SOD as a drug target that plays a crucial role in shielding cells from damage caused by ROS. Genome mining resulted in identification of five potential ligands (secondary metabolites) from B. licheniformis such as, Bacillibactin (BAC), Paenibactin (PAE), Fengycin (FEN), Surfactin (SUR) and Lichenysin (LIC). Molecular docking was used to predict and analyze the binding interactions between these five ligands and target protein SOD. The resulting protein–ligand complexes were further analyzed for their motions and interactions of atoms and molecules over 250 ns using molecular dynamics (MD) simulation analysis. The analysis of MD simulations suggests, Bacillibactin as the probable candidate to arrest the activities of SOD. All the five compounds reported in this study were found to act by directly/indirectly interacting with ROS molecules, such as superoxide radicals (O₂–) and hydrogen peroxide (H₂O₂), and transforming them into less reactive species. This antioxidant activity contributes to its protective effects against oxidative stress-induced damage in cells making them likely candidate for various applications, including in the development of antioxidant-based therapies, nutraceuticals, and functional foods.     

Impact on S. aureus and E. coli Membranes of Treatment with Chlorhexidine and Alcohol Solutions: Insights from Molecular Simulations and Nuclear Magnetic Resonance

Membranes form the first line of defence of bacteria against potentially harmful molecules in the surrounding environment. Understanding the protective properties of these membranes represents an important step towards development of targeted anti-bacterial agents such as sanitizers. Use of propanol, isopropanol and chlorhexidine can significantly decrease the threat imposed by bacteria in the face of growing anti-bacterial resistance via mechanisms that include membrane disruption. Here we have employed molecular dynamics simulations and nuclear magnetic resonance to explore the impact of chlorhexidine and alcohol on the S. aureus cell membrane, as well as the E. coli inner and outer membranes. We identify how sanitizer components partition into these bacterial membranes, and show that chlorhexidine is instrumental in this process.     

Identification of potent pyrazole based APELIN receptor (APJ) agonists

The apelinergic system comprises the apelin receptor and its cognate apelin and elabela peptide ligands of various lengths. This system has become an increasingly attractive target for pulmonary and cardiometabolic diseases. Small molecule regulators of this receptor with good drug-like properties are needed. Recently, we discovered a novel pyrazole based small molecule agonist 8 of the apelin receptor (EC₅₀ = 21.5 µM, K_i = 5.2 µM) through focused screening which was further optimized to initial lead 9 (EC₅₀ = 0.800 µM, K_i = 1.3 µM). In our efforts to synthesize more potent agonists and to explore the structural features important for apelin receptor agonism, we carried out structural modifications at N1 of the pyrazole core as well as the amino acid side-chain of 9. Systematic modifications at these two positions provided potent small molecule agonists exhibiting EC₅₀ values of <100 nM. Recruitment of β-arrestin as a measure of desensitization potential of select compounds was also investigated. Functional selectivity was a feature of several compounds with a bias towards calcium mobilization over β-arrestin recruitment. These compounds may be suitable as tools for in vivo studies of apelin receptor function.     

Expression profiles and functional characterization of common carp (Cyprinus carpio) T2Rs

Bitter taste perception is mediated by a family of G protein-coupled receptors (T2Rs) in vertebrates. Common carp (Cyprinus carpio), which has experienced an additional round of whole genome duplication during the course of evolution, has a small number of T2R genes similar to zebrafish, a closely related cyprinid fish species, and their expression pattern at the cellular level or their cognate ligands have not been elucidated yet. Here, we showed through in situ hybridization experiments, that three common carp T2R (ccT2R) genes encoding ccT2R200-1, ccT2R202-1, and ccT2R202-2, were specifically expressed in the subsets of taste receptor cells in the lips and gill rakers. ccT2R200-1 was co-expressed with genes encoding downstream signal transduction molecules, such as PLC-β2 and Gαia. Heterologous expression system revealed that each ccT2R showed narrowly, intermediately, or broadly tuned ligand specificity, as in the case of zebrafish T2Rs. However, ccT2Rs showed different ligand profiles from their orthologous zebrafish T2Rs previously reported. Finally, we identified three ccT2Rs, namely ccT2R200-1, ccT2R200-2, and ccT2R203-1, to be activated by natural bitter compounds, andrographolide and/or picrotoxinin, which elicited no response to zebrafish T2Rs, in a dose-dependent manner. These results suggest that some ccT2Rs may have evolved to function in the oral cavity as taste receptors for natural bitter compounds found in the habitats in a species-specific manner.     

Generation of highly potent DYRK1A-dependent inducers of human β-Cell replication via Multi-Dimensional compound optimization

Small molecule stimulation of β-cell regeneration has emerged as a promising therapeutic strategy for diabetes. Although chemical inhibition of dual specificity tyrosine-phosphorylation-regulated kinase 1A (DYRK1A) is sufficient to enhance β-cell replication, current lead compounds have inadequate cellular potency for in vivo application. Herein, we report the clinical stage anti-cancer kinase inhibitor OTS167 as a structurally novel, remarkably potent DYRK1A inhibitor and inducer of human β-cell replication. Unfortunately, OTS167’s target promiscuity and cytotoxicity curtails utility. To tailor kinase selectivity towards DYRK1A and reduce cytotoxicity we designed a library of fifty-one OTS167 derivatives based upon a modeled structure of the DYRK1A-OTS167 complex. Indeed, derivative characterization yielded several leads with exceptional DYRK1A inhibition and human β-cell replication promoting potencies but substantially reduced cytotoxicity. These compounds are the most potent human β-cell replication-promoting compounds yet described and exemplify the potential to purposefully leverage off-target activities of advanced stage compounds for a desired application.     

Integrating structure-based approaches in generative molecular design

Generative molecular design for drug discovery and development has seen a recent resurgence promising to improve the efficiency of the design-make-test-analyse cycle; by computationally exploring much larger chemical spaces than traditional virtual screening techniques. However, most generative models thus far have only utilized small-molecule information to train and condition de novo molecule generators. Here, we instead focus on recent approaches that incorporate protein structure into de novo molecule optimization in an attempt to maximize the predicted on-target binding affinity of generated molecules. We summarize these structure integration principles into either distribution learning or goal-directed optimization and for each case whether the approach is protein structure-explicit or implicit with respect to the generative model. We discuss recent approaches in the context of this categorization and provide our perspective on the future direction of the field.     

Inhibition of Porphyromonas gingivalis peptidyl arginine deiminase,a virulence factor,by antioxidant-rich Cratoxylum cochinchinense: In vitro and in silico evaluation

Porphyromonas gingivalis, the cause of periodontitis, is also linked to many systemic disorders due to its citrullination capability from a unique peptidyl arginine deiminase (PPAD). Protein citrullination is able to trigger an autoimmune response, increasing the severity of rheumatoid arthritis. The main objective of this study is to evaluate the inhibitory activity of Cratoxylym cochinchinense leaves extract towards the PPAD in vitro and in silico. Methanolic extract of Cratoxylum cochinchinense (CCM) was tested for total phenolic and flavonoid contents along with antioxidative assays. Inhibition of PPAD activities was conducted thereafter using recombinant PPAD in cell lysate. Phytocompounds postulated present in the CCM such as mangiferin, vismiaquinone A, δ-tocotrienol and α-tocotrienol and canophyllol were used as ligands in a simulated docking study against PPAD. Results obtained indicated high antioxidant potential in CCM while recording abundant phenolic (129.0 ± 2.5495 mg GA/g crude extract) and flavonoid (159.0 ± 2.1529 mg QE/g crude extract) contents. A dose-dependent inhibition of PPAD was observed when CCM was evaluated at various concentrations. CCM at 1 mg/mL exhibited citrulline concentration of 24.37 ± 3.25 mM which was 5 times lower than the negative control (114.23 ± 3.31 mM). Molecular docking simulation revealed that mangiferin and vismiaquinone A engaged in H-bonding and pi-pi interactions with important active site residues (Asp130, Arg152, Arg154 and Trp127) of PPAD and could be the potential phytochemicals that accounted for the inhibitory activities observed in the methanolic leaves extract. As such, CCM could be further explored for its therapeutic properties not only for periodontitis, but also for other systemic diseases like rheumatoid arthritis.     

Subsets of Slow Dynamic Modes Reveal Global Information Sources as Allosteric Sites

Allostery is a key biological control mechanism, and dynamic information flow provides a perspective to describe allosteric interactions in causal relationships. Here, as a novel implementation of the Gaussian Network Model (GNM) based Transfer Entropy (TE) calculations, we show that the dissection of dynamic information into subsets of slow dynamic modes discloses different layers of multi-directional allosteric pathways inherent in a given protein structure. In these subsets of slow modes, the degree of collectivity (Col) in the information transfer of residues with their TE values (TECol score) identifies distinct residues as powerful effectors, global information sources; showing themselves with a high dynamic capacity to collectively disseminate information to others. As exemplified on aspartate transcarbamoylase (ATCase), Na⁺/K⁺-adenosine triphosphatase (Na⁺/K⁺-ATPase), and human transient receptor potential melastatin 2 (TRPM2) along with a dataset of 20 proteins, these specific residues are associated with known active and allosteric sites. These information source residues, which collectively control others and lead allosteric communication pathways, hint at plausible binding sites for structure-based rational drug design.     

Utilizing in silico and in vitro methods to identify possible binding sites of a novel ligand against Pseudomonas aeruginosa phospholipase toxin ExoU

Multi-drug resistant infections caused by the opportunistic pathogen, Pseudomonas aeruginosa (P. aeruginosa), are a continuing problem that contribute to morbidity and mortality in immunocompromised hosts such as cystic fibrosis (CF), wound and burn patients. The bacterial toxin ExoU is one of four potent toxins that P. aeruginosa secretes into the epithelial cells of hosts. In this study, NMR Saturation Transfer Difference (STD) and in silico Schrödinger Computational Modeling were used to identify a possible binding site of a novel ligand methoctramine targeting ExoU. Future project goals will be to design a structure activity relationship (SAR) study of methoctramine and ExoU and lead to a new drug solving ExoU toxicity P. aeruginosa exerts in the clinical environment.     

Structure-Function Studies of Two Yeast Homing Endonucleases that Evolved to Cleave Identical Targets with Dissimilar Rates and Specificities

The LAGLIDADG family of homing endonucleases (LHEs) bind to and cleave their DNA recognition sequences with high specificity. Much of our understanding for how these proteins evolve their specificities has come from studying LHE homologues. To gain insight into the molecular basis of LHE specificity, we characterized I-WcaI, the homologue of the Saccharomyces cerevisiae I-SceI LHE found in Wickerhamomyces canadensis. Although I-WcaI and I-SceI cleave the same recognition sequence, expression of I-WcaI, but not I-SceI, is toxic in bacteria. Toxicity suppressing mutations frequently occur at I-WcaI residues critical for activity and I-WcaI cleaves many more non-cognate sequences in the Escherichia coli genome than I-SceI, suggesting I-WcaI endonuclease activity is the basis of toxicity. In vitro, I-WcaI is a more active and a less specific endonuclease than I-SceI, again accounting for the observed toxicity in vivo. We determined the X-ray crystal structure of I-WcaI bound to its cognate target site and found that I-WcaI and I-SceI use residues at different positions to make similar base-specific contacts. Furthermore, in some regions of the DNA interface where I-WcaI specificity is lower, the protein makes fewer DNA contacts than I-SceI. Taken together, these findings demonstrate the plastic nature of LHE site recognition and suggest that I-WcaI and I-SceI are situated at different points in their evolutionary pathways towards acquiring target site specificity.     

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.     

In vitro Selection of High Affinity DNA and RNA Aptamers that Detect Hepatitis C Virus Core Protein of Genotypes 1 to 4 and Inhibit Virus Production in Cell Culture

Hepatitis C virus (HCV) core is a highly conserved and multifunctional protein that forms the viral capsid, making it an attractive target for HCV detection and inhibition. Aptamers are in vitro selected, single-stranded nucleic acids (RNA or ssDNA) with growing applicability in viral diagnostics and therapy. We have carried out DNA and RNA in vitro selection against six different variants of HCV core protein: two versions of the full-length protein of genotype 1, and the hydrophilic domain of genotypes 1 to 4. The aptamer populations obtained were analyzed by means of Ultra-Deep Sequencing (UDS), the most abundant sequences were identified and a number of highly represented sequence motifs were unveiled. Affinity (measured as the dissociation constant, Kd) of the most abundant DNA and RNA aptamers were quantified using Enzyme-Linked OligoNucleotide Assay (ELONA)-based methods. Some aptamers with nanomolar or subnanomolar Kd values (as low as 0.4 nM) were the common outcome of DNA and RNA selections against different HCV core variants. They were tested in sandwich and competitive biosensor assays, reaching a limit of detection for HCV core of 2 pM. Additionally, the two most prevalent and high affinity aptamers were assayed in Huh-7.5 reporter cell lines infected with HCV, where they decreased both the viral progeny titer and the extracellular viral RNA level, while increasing the amount of intracellular viral RNA. Our results suggest that these aptamers inhibit HCV capsid assembly and virion formation, thus making them good candidate molecules for the design of novel therapeutic approaches for hepatitis C.     

Anti-aging and antioxidant of four traditional malaysian plants using simplex centroid mixture design approach

Aging is a naturally biological process with adverse effects. The continuous accumulation of reactive oxygen species (ROS) trigger cellular and tissue damage by activating several aging enzymes. The antioxidant properties of traditional medicinal plants used by Jakun aborigine’s community are a promising approach to alleviate aging process and prevent Alzheimer. The aim of the current investigation was to optimize a novel anti-aging formulation from traditional plants (Cnestis palala stem, Urceola micrantha stem, Marantodes pumilum stem and Microporus xanthopus fruiting bodies) using simplex centroid mixture design (SCMD). After selecting the optimal formulations based on desirability function of antioxidant activity (DPPḢ, ABTS^̇+ and FRAP), they were further examined against the activity of aging-related-enzymes (collagenase, tyrosinase, acetyl- and butyrylcholinesterase). The single extracts of C. palala, U. micrantha and the binary mixture of C. palala and U. micrantha were the optimal formulations with high antioxidant activities. Single extract of U. micrantha showed the highest inhibition towards matrix metalloproteinase-1 (49.44 ± 4.11 %), while C. palala water extract showed highest inhibitions towards tyrosinase (14.06 ± 0.31%), acetylcholinesterase (32.92 ± 2.13%) and butyrylcholinesterase (34.89 ± 2.84%) enzymes. The single extracts of C. palala and U. micrantha displayed better activity as compared to the binary mixture formulation. In conclusion, these findings could be a baseline for further exploration of novel anti-aging agents from natural resources.