Analysis of a Novel pncA Mutation for Susceptibility to Pyrazinamide Therapy

Chemical genetic interactions elucidate pathways controlling tuberculosis antibiotic efficacy during infection

Successful tuberculosis therapy requires treatment with an unwieldy multidrug combination for several months. Thus, there is a growing need to identify novel genetic vulnerabilities that can be leveraged to develop new, more effective antitubercular drugs. Consequently, recent efforts to optimize tuberculosis (TB) therapy have exploited Mycobacterium tuberculosis (Mtb) chemical genetics to identify pathways influencing antibiotic efficacy, novel mechanisms of antibiotic action, and new targets for TB drug discovery. However, the influence of the complex host environment on these interactions remains largely unknown, leaving the therapeutic potential of the identified targets unclear. In this study, we leveraged a library of conditional mutants targeting 467 essential Mtb genes to characterize the chemical-genetic interactions (CGIs) with TB drugs directly in the mouse infection model. We found that these in vivo CGIs differ significantly from those identified in vitro. Both drug-specific and drug-agnostic effects were identified, and many were preserved during treatment with a multidrug combination, suggesting numerous strategies for enhancing therapy. This work also elucidated the complex effects of pyrazinamide (PZA), a drug that relies on aspects of the infection environment for efficacy. Specifically, our work supports the importance of coenzyme A synthesis- inhibition during infection, as well as the antagonistic effect of iron limitation on PZA activity. In addition, we found that inhibition of thiamine and purine synthesis increases PZA efficacy, suggesting additional therapeutically exploitable metabolic dependencies. Our findings present a map of the unique in vivo CGIs, characterizing the mechanism of PZA activity in vivo and identifying potential targets for TB drug development.

Exploring the effects of missense mutations on protein thermodynamics through structure-based approaches: findings from the CAGI6 challenges

Missense mutations are known contributors to diverse genetic disorders, due to their subtle, single amino acid changes imparted on the resultant protein. Because of this, understanding the impact of these mutations on protein stability and function is crucial for unravelling disease mechanisms and developing targeted therapies. The Critical Assessment of Genome Interpretation (CAGI) provides a valuable platform for benchmarking state-of-the-art computational methods in predicting the impact of disease-related mutations on protein thermodynamics. Here we report the performance of our comprehensive platform of structure-based computational approaches to evaluate mutations impacting protein structure and function on 3 challenges from CAGI6: Calmodulin, MAPK1 and MAPK3. Our stability predictors have achieved correlations of up to 0.74 and AUCs of 1 when predicting changes in ΔΔG for MAPK1 and MAPK3, respectively, and AUC of up to 0.75 in the Calmodulin challenge. Overall, our study highlights the importance of structure-based approaches in understanding the effects of missense mutations on protein thermodynamics. The results obtained from the CAGI6 challenges contribute to the ongoing efforts to enhance our understanding of disease mechanisms and facilitate the development of personalised medicine approaches.

Chemical genetic interactions elucidate pathways controlling tuberculosis antibiotic efficacy during infection

Successful tuberculosis therapy requires treatment with an unwieldy multidrug combination for several months. Thus, there is a growing need to identify novel genetic vulnerabilities that can be leveraged to develop new, more effective antitubercular drugs. Consequently, recent efforts to optimize TB therapy have exploited Mtb chemical genetics to identify pathways influencing antibiotic efficacy, novel mechanisms of antibiotic action, and new targets for TB drug discovery. However, the influence of the complex host environment on these interactions remains largely unknown, leaving the therapeutic potential of the identified targets unclear. In this study, we leveraged a library of conditional mutants targeting 467 essential Mtb genes to characterize the chemical-genetic interactions (CGIs) with TB drugs directly in the mouse infection model. We found that these in vivo CGIs differ significantly from those identified in vitro . Both drug-specific and drug-agnostic effects were identified, and many were preserved during treatment with a multidrug combination, suggesting numerous strategies for enhancing therapy. This work also elucidated the complex effects of pyrazinamide (PZA), a drug that relies on aspects of the infection environment for efficacy. Specifically, our work supports the importance of coenzyme A synthesis inhibition during infection, as well as the antagonistic effect of iron limitation on PZA activity. In addition, we found that inhibition of thiamine and purine synthesis increases PZA efficacy, suggesting novel therapeutically exploitable metabolic dependencies. Our findings present a map of the unique in vivo CGIs, characterizing the mechanism of PZA activity in vivo and identifying novel targets for TB drug development.

Significance

The inevitable rise of multi-drug-resistant tuberculosis underscores the urgent need for new TB drugs and novel drug targets while prioritizing synergistic drug combinations. Chemical-genetic interaction (CGI) studies have delineated bacterial pathways influencing antibiotic efficacy and uncovered druggable pathways that synergize with TB drugs. However, most studies are conducted in vitro , limiting our understanding of how the host environment influences drug-mutant interactions. Using an inducible mutant library targeting essential Mtb genes to characterize CGIs during infection, this study reveals that CGIs are both drug-specific and drug-agnostic and differ significantly from those observed in vitro . Synergistic CGIs comprised distinct metabolic pathways mediating antibiotic efficacy, revealing novel drug mechanisms of action, and defining potential drug targets that would synergize with frontline antitubercular drugs.

Engineering G protein-coupled receptors for stabilization

G protein-coupled receptors (GPCRs) are one of the most important families of targets for drug discovery. One of the limiting steps in the study of GPCRs has been their stability, with significant and time-consuming protein engineering often used to stabilize GPCRs for structural characterization and drug screening. Unfortunately, computational methods developed using globular soluble proteins have translated poorly to the rational engineering of GPCRs. To fill this gap, we propose GPCR-tm, a novel and personalized structurally driven web-based machine learning tool to study the impacts of mutations on GPCR stability. We show that GPCR-tm performs as well as or better than alternative methods, and that it can accurately rank the stability changes of a wide range of mutations occurring in various types of class A GPCRs. GPCR-tm achieved Pearson's correlation coefficients of 0.74 and 0.46 on 10-fold cross-validation and blind test sets, respectively. We observed that the (structural) graph-based signatures were the most important set of features for predicting destabilizing mutations, which points out that these signatures properly describe the changes in the environment where the mutations occur. More specifically, GPCR-tm was able to accurately rank mutations based on their effect on protein stability, guiding their rational stabilization. GPCR-tm is available through a user-friendly web server at https://biosig.lab.uq.edu.au/gpcr_tm/.

Characterizing and predicting ccRCC-causing missense mutations in Von Hippel-Lindau disease

BACKGROUND

Mutations within the Von Hippel-Lindau (VHL) tumor suppressor gene are known to cause VHL disease, which is characterized by the formation of cysts and tumors in multiple organs of the body, particularly clear cell renal cell carcinoma (ccRCC). A major challenge in clinical practice is determining tumor risk from a given mutation in the VHL gene. Previous efforts have been hindered by limited available clinical data and technological constraints.

METHODS

To overcome this, we initially manually curated the largest set of clinically validated VHL mutations to date, enabling a robust assessment of existing predictive tools on an independent test set. Additionally, we comprehensively characterized the effects of mutations within VHL using in silico biophysical tools describing changes in protein stability, dynamics and affinity to binding partners to provide insights into the structure-phenotype relationship. These descriptive properties were used as molecular features for the construction of a machine learning model, designed to predict the risk of ccRCC development as a result of a VHL missense mutation.

RESULTS

Analysis of our model showed an accuracy of 0.81 in the identification of ccRCC-causing missense mutations, and a Matthew's Correlation Coefficient of 0.44 on a non-redundant blind test, a significant improvement in comparison to the previous available approaches.

CONCLUSION

This work highlights the power of using protein 3D structure to fully explore the range of molecular and functional consequences of genomic variants. We believe this optimized model will better enable its clinical implementation and assist guiding patient risk stratification and management.

AI-Driven Enhancements in Drug Screening and Optimization

The greatest challenge in drug discovery remains the high rate of attrition across the different phases of the process, which cost the industry billions of dollars every year. While all phases remain crucial to ensure pharmaceutical-level safety, quality, and efficacy of the end product, streamlining these efforts toward compounds with success potential is pivotal for a more efficient and cost-effective process. The use of artificial intelligence (AI) within the pharmaceutical industry aims at just this, and has applications in preclinical screening for biological activity, optimization of pharmacokinetic properties for improved drug formulation, early toxicity prediction which reduces attrition, and pre-emptively screening for genetic changes in the biological target to improve therapeutic longevity. Here, we present a series of in silico tools that address these applications in small molecule development and describe how they can be embedded within the current pharmaceutical development pipeline.

Identifying Innate Resistance Hotspots for SARS-CoV-2 Antivirals Using In Silico Protein Techniques

The development and approval of antivirals against SARS-CoV-2 has further equipped clinicians with treatment strategies against the COVID-19 pandemic, reducing deaths post-infection. Extensive clinical use of antivirals, however, can impart additional selective pressure, leading to the emergence of antiviral resistance. While we have previously characterized possible effects of circulating SARS-CoV-2 missense mutations on proteome function and stability, their direct effects on the novel antivirals remains unexplored. To address this, we have computationally calculated the consequences of mutations in the antiviral targets: RNA-dependent RNA polymerase and main protease, on target stability and interactions with their antiviral, nucleic acids, and other proteins. By analyzing circulating variants prior to antiviral approval, this work highlighted the inherent resistance potential of different genome regions. Namely, within the main protease binding site, missense mutations imparted a lower fitness cost, while the opposite was noted for the RNA-dependent RNA polymerase binding site. This suggests that resistance to nirmatrelvir/ritonavir combination treatment is more likely to occur and proliferate than that to molnupiravir. These insights are crucial both clinically in drug stewardship, and preclinically in the identification of less mutable targets for novel therapeutic design.

Using Graph-Based Signatures to Guide Rational Antibody Engineering

Antibodies are essential experimental and diagnostic tools and as biotherapeutics have significantly advanced our ability to treat a range of diseases. With recent innovations in computational tools to guide protein engineering, we can now rationally design better antibodies with improved efficacy, stability, and pharmacokinetics. Here, we describe the use of the mCSM web-based in silico suite, which uses graph-based signatures to rapidly identify the structural and functional consequences of mutations, to guide rational antibody engineering to improve stability, affinity, and specificity.

Structure activity relationship of pyrazinoic acid analogs as potential antimycobacterial agents

Tuberculosis (TB) remains a leading cause of infectious disease-related mortality and morbidity. Pyrazinamide (PZA) is a critical component of the first-line TB treatment regimen because of its sterilizing activity against non-replicating Mycobacterium tuberculosis (Mtb), but its mechanism of action has remained enigmatic. PZA is a prodrug converted by pyrazinamidase encoded by pncA within Mtb to the active moiety, pyrazinoic acid (POA) and PZA resistance is caused by loss-of-function mutations to pyrazinamidase. We have recently shown that POA induces targeted protein degradation of the enzyme PanD, a crucial component of the coenzyme A biosynthetic pathway essential in Mtb. Based on the newly identified mechanism of action of POA, along with the crystal structure of PanD bound to POA, we designed several POA analogs using structure for interpretation to improve potency and overcome PZA resistance. We prepared and tested ring and carboxylic acid bioisosteres as well as 3, 5, 6 substitutions on the ring to study the structure activity relationships of the POA scaffold. All the analogs were evaluated for their whole cell antimycobacterial activity, and a few representative molecules were evaluated for their binding affinity, towards PanD, through isothermal titration calorimetry. We report that analogs with ring and carboxylic acid bioisosteres did not significantly enhance the antimicrobial activity, whereas the alkylamino-group substitutions at the 3 and 5 position of POA were found to be up to 5 to 10-fold more potent than POA. Further development and mechanistic analysis of these analogs may lead to a next generation POA analog for treating TB.

Drug Degradation Caused by mce3R Mutations Confers Contezolid (MRX-I) Resistance in Mycobacterium tuberculosis

Contezolid (MRX-I), a safer antibiotic of the oxazolidinone class, is a promising new antibiotic with potent activity against Mycobacterium tuberculosis (MTB) both in vitro and in vivo. To identify resistance mechanisms of contezolid in MTB, we isolated several in vitro spontaneous contezolid-resistant MTB mutants, which exhibited 16-fold increases in the MIC of contezolid compared with the parent strain but were still unexpectedly susceptible to linezolid. Whole-genome sequencing revealed that most of the contezolid-resistant mutants bore mutations in the mce3R gene, which encodes a transcriptional repressor. The mutations in mce3R led to markedly increased expression of a monooxygenase encoding gene Rv1936. We then characterized Rv1936 as a putative flavin-dependent monooxygenase that catalyzes the degradation of contezolid into its inactive 2,3-dihydropyridin-4-one (DHPO) ring-opened metabolites, thereby conferring drug resistance. While contezolid is an attractive drug candidate with potent antimycobacterial activity and low toxicity, the occurrence of mutations in Mce3R should be considered when designing combination therapy using contezolid for treating tuberculosis.

Systematic evaluation of computational tools to predict the effects of mutations on protein stability in the absence of experimental structures

Changes in protein sequence can have dramatic effects on how proteins fold, their stability and dynamics. Over the last 20 years, pioneering methods have been developed to try to estimate the effects of missense mutations on protein stability, leveraging growing availability of protein 3D structures. These, however, have been developed and validated using experimentally derived structures and biophysical measurements. A large proportion of protein structures remain to be experimentally elucidated and, while many studies have based their conclusions on predictions made using homology models, there has been no systematic evaluation of the reliability of these tools in the absence of experimental structural data. We have, therefore, systematically investigated the performance and robustness of ten widely used structural methods when presented with homology models built using templates at a range of sequence identity levels (from 15% to 95%) and contrasted performance with sequence-based tools, as a baseline. We found there is indeed performance deterioration on homology models built using templates with sequence identity below 40%, where sequence-based tools might become preferable. This was most marked for mutations in solvent exposed residues and stabilizing mutations. As structure prediction tools improve, the reliability of these predictors is expected to follow, however we strongly suggest that these factors should be taken into consideration when interpreting results from structure-based predictors of mutation effects on protein stability.

Estimating tuberculosis drug resistance amplification rates in high-burden settings

Background

Antimicrobial resistance develops following the accrual of mutations in the bacterial genome, and may variably impact organism fitness and hence, transmission risk. Classical representation of tuberculosis (TB) dynamics using a single or two strain (DS/MDR-TB) model typically does not capture elements of this important aspect of TB epidemiology. To understand and estimate the likelihood of resistance spreading in high drug-resistant TB incidence settings, we used epidemiological data to develop a mathematical model of Mycobacterium tuberculosis (Mtb) transmission.

Methods

A four-strain (drug-susceptible (DS), isoniazid mono-resistant (INH-R), rifampicin mono-resistant (RIF-R) and multidrug-resistant (MDR)) compartmental deterministic Mtb transmission model was developed to explore the progression from DS- to MDR-TB in The Philippines and Viet Nam. The models were calibrated using data from national tuberculosis prevalence (NTP) surveys and drug resistance surveys (DRS). An adaptive Metropolis algorithm was used to estimate the risks of drug resistance amplification among unsuccessfully treated individuals.

Results

The estimated proportion of INH-R amplification among failing treatments was 0.84 (95% CI 0.79–0.89) for The Philippines and 0.77 (95% CI 0.71–0.84) for Viet Nam. The proportion of RIF-R amplification among failing treatments was 0.05 (95% CI 0.04–0.07) for The Philippines and 0.011 (95% CI 0.010–0.012) for Viet Nam.

Conclusion

The risk of resistance amplification due to treatment failure for INH was dramatically higher than RIF. We observed RIF-R strains were more likely to be transmitted than acquired through amplification, while both mechanisms of acquisition were important contributors in the case of INH-R. These findings highlight the complexity of drug resistance dynamics in high-incidence settings, and emphasize the importance of prioritizing testing algorithms which allow for early detection of INH-R.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07067-1.

HGDiscovery: An online tool providing functional and phenotypic information on novel variants of homogentisate 1,2- dioxigenase

Alkaptonuria (AKU), a rare genetic disorder, is characterized by the accumulation of homogentisic acid (HGA) in the body. Affected individuals lack functional levels of an enzyme required to breakdown HGA. Mutations in the homogentisate 1,2-dioxygenase (HGD) gene cause AKU and they are responsible for deficient levels of functional HGD, which, in turn, leads to excess levels of HGA. Although HGA is rapidly cleared from the body by the kidneys, in the long term it starts accumulating in various tissues, especially cartilage. Over time (rarely before adulthood), it eventually changes the color of affected tissue to slate blue or black. Here we report a comprehensive mutation analysis of 111 pathogenic and 190 non-pathogenic HGD missense mutations using protein structural information. Using our comprehensive suite of graph-based signature methods, mCSM complemented with sequence-based tools, we studied the functional and molecular consequences of each mutation on protein stability, interaction and evolutionary conservation. The scores generated from the structure and sequence-based tools were used to train a supervised machine learning algorithm with 89% accuracy. The empirical classifier was used to generate the variant phenotype for novel HGD missense mutations. All this information is deployed as a user friendly freely available web server called HGDiscovery (https://biosig.lab.uq.edu.au/hgdiscovery/).

•

Functional and phenotypic consequences of HGD non-synonymous variations.

•

Biophysical, structural and evolutionary analysis of novel and known clinical variants.

•

Pathogenic mutations affected protein stability and conformational flexibility.

•

Pathogenic mutations associated with deleterious scores for sequence-based features.

•

HGDiscovery (http://biosig.unimelb.edu.au/hgdiscovery/) – webserver.

Structure-guided machine learning prediction of drug resistance mutations in Abelson 1 kinase

Kinases play crucial roles in cellular signalling and biological processes with their dysregulation associated with diseases, including cancers. Kinase inhibitors, most notably those targeting ABeLson 1 (ABL1) kinase in chronic myeloid leukemia, have had a significant impact on cancer survival, yet emergence of resistance mutations can reduce their effectiveness, leading to therapeutic failure. Limited effort, however, has been devoted to developing tools to accurately identify ABL1 resistance mutations, as well as providing insights into their molecular mechanisms. Here we investigated the structural basis of ABL1 mutations modulating binding affinity of eight FDA-approved drugs. We found mutations impair affinity of type I and type II inhibitors differently and used this insight to developed a novel web-based diagnostic tool, SUSPECT-ABL, to pre-emptively predict resistance profiles and binding free-energy changes (ΔΔG) of all possible ABL1 mutations against inhibitors with different binding modes. Resistance mutations in ABL1 were successfully identified, achieving a Matthew's Correlation Coefficient of up to 0.73 and the resulting change in ligand binding affinity with a Pearson's correlation of up to 0.77, with performances consistent across non-redundant blind tests. Through an in silico saturation mutagenesis, our tool has identified possibly emerging resistance mutations, which offers opportunities for in vivo experimental validation. We believe SUSPECT-ABL will be an important tool not just for improving precision medicine efforts, but for facilitating the development of next-generation inhibitors that are less prone to resistance. We have made our tool freely available at http://biosig.unimelb.edu.au/suspect_abl/.

Structural and Genomic Insights Into Pyrazinamide Resistance in Mycobacterium tuberculosis Underlie Differences Between Ancient and Modern Lineages

Resistance to drugs used to treat tuberculosis disease (TB) continues to remain a public health burden, with missense point mutations in the underlying Mycobacterium tuberculosis bacteria described for nearly all anti-TB drugs. The post-genomics era along with advances in computational and structural biology provide opportunities to understand the interrelationships between the genetic basis and the structural consequences of M. tuberculosis mutations linked to drug resistance. Pyrazinamide (PZA) is a crucial first line antibiotic currently used in TB treatment regimens. The mutational promiscuity exhibited by the pncA gene (target for PZA) necessitates computational approaches to investigate the genetic and structural basis for PZA resistance development. We analysed 424 missense point mutations linked to PZA resistance derived from ∼35K M. tuberculosis clinical isolates sourced globally, which comprised the four main M. tuberculosis lineages (Lineage 1-4). Mutations were annotated to reflect their association with PZA resistance. Genomic measures (minor allele frequency and odds ratio), structural features (surface area, residue depth and hydrophobicity) and biophysical effects (change in stability and ligand affinity) of point mutations on pncA protein stability and ligand affinity were assessed. Missense point mutations within pncA were distributed throughout the gene, with the majority (>80%) of mutations with a destabilising effect on protomer stability and on ligand affinity. Active site residues involved in PZA binding were associated with multiple point mutations highlighting mutational diversity due to selection pressures at these functionally important sites. There were weak associations between genomic measures and biophysical effect of mutations. However, mutations associated with PZA resistance showed statistically significant differences between structural features (surface area and residue depth), but not hydrophobicity score for mutational sites. Most interestingly M. tuberculosis lineage 1 (ancient lineage) exhibited a distinct protein stability profile for mutations associated with PZA resistance, compared to modern lineages.

mmCSM-PPI: predicting the effects of multiple point mutations on protein-protein interactions

Protein-protein interactions play a crucial role in all cellular functions and biological processes and mutations leading to their disruption are enriched in many diseases. While a number of computational methods to assess the effects of variants on protein-protein binding affinity have been proposed, they are in general limited to the analysis of single point mutations and have been shown to perform poorly on independent test sets. Here, we present mmCSM-PPI, a scalable and effective machine learning model for accurately assessing changes in protein-protein binding affinity caused by single and multiple missense mutations. We expanded our well-established graph-based signatures in order to capture physicochemical and geometrical properties of multiple wild-type residue environments and integrated them with substitution scores and dynamics terms from normal mode analysis. mmCSM-PPI was able to achieve a Pearson's correlation of up to 0.75 (RMSE = 1.64 kcal/mol) under 10-fold cross-validation and 0.70 (RMSE = 2.06 kcal/mol) on a non-redundant blind test, outperforming existing methods. Our method is freely available as a user-friendly and easy-to-use web server and API at http://biosig.unimelb.edu.au/mmcsm_ppi.

Distinguishing between PTEN clinical phenotypes through mutation analysis

Phosphate and tensin homolog on chromosome ten (PTEN) germline mutations are associated with an overarching condition known as PTEN hamartoma tumor syndrome. Clinical phenotypes associated with this syndrome range from macrocephaly and autism spectrum disorder to Cowden syndrome, which manifests as multiple noncancerous tumor-like growths (hamartomas), and an increased predisposition to certain cancers. It is unclear, however, the basis by which mutations might lead to these very diverse phenotypic outcomes. Here we show that, by considering the molecular consequences of mutations in PTEN on protein structure and function, we can accurately distinguish PTEN mutations exhibiting different phenotypes. Changes in phosphatase activity, protein stability, and intramolecular interactions appeared to be major drivers of clinical phenotype, with cancer-associated variants leading to the most drastic changes, while ASD and non-pathogenic variants associated with more mild and neutral changes, respectively. Importantly, we show via saturation mutagenesis that more than half of variants of unknown significance could be associated with disease phenotypes, while over half of Cowden syndrome mutations likely lead to cancer. These insights can assist in exploring potentially important clinical outcomes delineated by PTEN variation.

Genome analyses of 174 strains of Mycobacterium tuberculosis provide insight into the evolution of drug resistance and reveal potential drug targets

Mycobacterium tuberculosis is a known human pathogen that causes the airborne infectious disease tuberculosis (TB). Every year TB infects millions of people worldwide. The emergence of multi-drug resistant (MDR), extensively drug resistant (XDR) and totally drug resistant (TDR) M. tuberculosis strains against the first- and second-line anti-TB drugs has created an urgent need for the development and implementation of new drug strategies. In this study, the complete genomes of 174 strains of M. tuberculosis are analysed to understand the evolution of molecular drug target (MDT) genes. Phylogenomic placements of M. tuberculosis strains depicted close association and temporal clustering. Selection pressure analysis by deducing the ratio of non-synonymous to synonymous substitution rates (dN/dS) in 51 MDT genes of the 174 M. tuberculosis strains led to categorizing these genes into diversifying (D, dN/dS>0.70), moderately diversifying (MD, dN/dS=0.35-0.70) and stabilized (S, dN/dS<0.35) genes. The genes rpsL, gidB, pncA and ahpC were identified as diversifying, and Rv0488, kasA, ndh, ethR, ethA, embR and ddn were identified as stabilized genes. Furthermore, sequence similarity networks were drawn that supported these divisions. In the multiple sequence alignments of diversifying and stabilized proteins, previously reported resistance mutations were checked to predict sensitive and resistant strains of M. tuberculosis. Finally, to delineate the potential of stabilized or least diversified genes/proteins as anti-TB drug targets, protein-protein interactions of MDT proteins with human proteins were analysed. We predict that kasA (dN/dS=0.29), a stabilized gene that encodes the most host-interacting protein, KasA, should serve as a potential drug target for the treatment of TB.

ThermoMutDB: a thermodynamic database for missense mutations

Proteins are intricate, dynamic structures, and small changes in their amino acid sequences can lead to large effects on their folding, stability and dynamics. To facilitate the further development and evaluation of methods to predict these changes, we have developed ThermoMutDB, a manually curated database containing >14,669 experimental data of thermodynamic parameters for wild type and mutant proteins. This represents an increase of 83% in unique mutations over previous databases and includes thermodynamic information on 204 new proteins. During manual curation we have also corrected annotation errors in previously curated entries. Associated with each entry, we have included information on the unfolding Gibbs free energy and melting temperature change, and have associated entries with available experimental structural information. ThermoMutDB supports users to contribute to new data points and programmatic access to the database via a RESTful API. ThermoMutDB is freely available at: http://biosig.unimelb.edu.au/thermomutdb.

HARP: a database of structural impacts of systematic missense mutations in drug targets of Mycobacterium leprae

Computational Saturation Mutagenesis is an in-silico approach that employs systematic mutagenesis of each amino acid residue in the protein to all other amino acid types, and predicts changes in thermodynamic stability and affinity to the other subunits/protein counterparts, ligands and nucleic acid molecules. The data thus generated are useful in understanding the functional consequences of mutations in antimicrobial resistance phenotypes. In this study, we applied computational saturation mutagenesis to three important drug-targets in Mycobacterium leprae (M. leprae) for the drugs dapsone, rifampin and ofloxacin namely Dihydropteroate Synthase (DHPS), RNA Polymerase (RNAP) and DNA Gyrase (GYR), respectively. M. leprae causes leprosy and is an obligate intracellular bacillus with limited protein structural information associating mutations with phenotypic resistance outcomes in leprosy. Experimentally solved structures of DHPS, RNAP and GYR of M. leprae are not available in the Protein Data Bank, therefore, we modelled the structures of these proteins using template-based comparative modelling and introduced systematic mutations in each model generating 80,902 mutations and mutant structures for all the three proteins. Impacts of mutations on stability and protein-subunit, protein-ligand and protein-nucleic acid affinities were computed using various in-house developed and other published protein stability and affinity prediction software. A consensus impact was estimated for each mutation using qualitative scoring metrics for physicochemical properties and by a categorical grouping of stability and affinity predictions. We developed a web database named HARP (a database of Hansen's Disease Antimicrobial Resistance Profiles), which is accessible at the URL - https://harp-leprosy.org and provides the details to each of these predictions.

Combining structure and genomics to understand antimicrobial resistance

Antimicrobials against bacterial, viral and parasitic pathogens have transformed human and animal health. Nevertheless, their widespread use (and misuse) has led to the emergence of antimicrobial resistance (AMR) which poses a potentially catastrophic threat to public health and animal husbandry. There are several routes, both intrinsic and acquired, by which AMR can develop. One major route is through non-synonymous single nucleotide polymorphisms (nsSNPs) in coding regions. Large scale genomic studies using high-throughput sequencing data have provided powerful new ways to rapidly detect and respond to such genetic mutations linked to AMR. However, these studies are limited in their mechanistic insight. Computational tools can rapidly and inexpensively evaluate the effect of mutations on protein function and evolution. Subsequent insights can then inform experimental studies, and direct existing or new computational methods. Here we review a range of sequence and structure-based computational tools, focussing on tools successfully used to investigate mutational effect on drug targets in clinically important pathogens, particularly Mycobacterium tuberculosis. Combining genomic results with the biophysical effects of mutations can help reveal the molecular basis and consequences of resistance development. Furthermore, we summarise how the application of such a mechanistic understanding of drug resistance can be applied to limit the impact of AMR.

Prediction of rifampicin resistance beyond the RRDR using structure-based machine learning approaches

Rifampicin resistance is a major therapeutic challenge, particularly in tuberculosis, leprosy, P. aeruginosa and S. aureus infections, where it develops via missense mutations in gene rpoB. Previously we have highlighted that these mutations reduce protein affinities within the RNA polymerase complex, subsequently reducing nucleic acid affinity. Here, we have used these insights to develop a computational rifampicin resistance predictor capable of identifying resistant mutations even outside the well-defined rifampicin resistance determining region (RRDR), using clinical M. tuberculosis sequencing information. Our tool successfully identified up to 90.9% of M. tuberculosis rpoB variants correctly, with sensitivity of 92.2%, specificity of 83.6% and MCC of 0.69, outperforming the current gold-standard GeneXpert-MTB/RIF. We show our model can be translated to other clinically relevant organisms: M. leprae, P. aeruginosa and S. aureus, despite weak sequence identity. Our method was implemented as an interactive tool, SUSPECT-RIF (StrUctural Susceptibility PrEdiCTion for RIFampicin), freely available at https://biosig.unimelb.edu.au/suspect_rif/ .

mCSM-membrane: predicting the effects of mutations on transmembrane proteins

Significant efforts have been invested into understanding and predicting the molecular consequences of mutations in protein coding regions, however nearly all approaches have been developed using globular, soluble proteins. These methods have been shown to poorly translate to studying the effects of mutations in membrane proteins. To fill this gap, here we report, mCSM-membrane, a user-friendly web server that can be used to analyse the impacts of mutations on membrane protein stability and the likelihood of them being disease associated. mCSM-membrane derives from our well-established mutation modelling approach that uses graph-based signatures to model protein geometry and physicochemical properties for supervised learning. Our stability predictor achieved correlations of up to 0.72 and 0.67 (on cross validation and blind tests, respectively), while our pathogenicity predictor achieved a Matthew's Correlation Coefficient (MCC) of up to 0.77 and 0.73, outperforming previously described methods in both predicting changes in stability and in identifying pathogenic variants. mCSM-membrane will be an invaluable and dedicated resource for investigating the effects of single-point mutations on membrane proteins through a freely available, user friendly web server at http://biosig.unimelb.edu.au/mcsm_membrane.

mCSM-AB2: guiding rational antibody design using graph-based signatures

MOTIVATION

A lack of accurate computational tools to guide rational mutagenesis has made affinity maturation a recurrent challenge in antibody (Ab) development. We previously showed that graph-based signatures can be used to predict the effects of mutations on Ab binding affinity.

RESULTS

Here we present an updated and refined version of this approach, mCSM-AB2, capable of accurately modelling the effects of mutations on Ab-antigen binding affinity, through the inclusion of evolutionary and energetic terms. Using a new and expanded database of over 1800 mutations with experimental binding measurements and structural information, mCSM-AB2 achieved a Pearson's correlation of 0.73 and 0.77 across training and blind tests, respectively, outperforming available methods currently used for rational Ab engineering.

AVAILABILITY AND IMPLEMENTATION

mCSM-AB2 is available as a user-friendly and freely accessible web server providing rapid analysis of both individual mutations or the entire binding interface to guide rational antibody affinity maturation at http://biosig.unimelb.edu.au/mcsm_ab2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DynaMut2: Assessing changes in stability and flexibility upon single and multiple point missense mutations

Predicting the effect of missense variations on protein stability and dynamics is important for understanding their role in diseases, and the link between protein structure and function. Approaches to estimate these changes have been proposed, but most only consider single‐point missense variants and a static state of the protein, with those that incorporate dynamics are computationally expensive. Here we present DynaMut2, a web server that combines Normal Mode Analysis (NMA) methods to capture protein motion and our graph‐based signatures to represent the wildtype environment to investigate the effects of single and multiple point mutations on protein stability and dynamics. DynaMut2 was able to accurately predict the effects of missense mutations on protein stability, achieving Pearson's correlation of up to 0.72 (RMSE: 1.02 kcal/mol) on a single point and 0.64 (RMSE: 1.80 kcal/mol) on multiple‐point missense mutations across 10‐fold cross‐validation and independent blind tests. For single‐point mutations, DynaMut2 achieved comparable performance with other methods when predicting variations in Gibbs Free Energy (ΔΔG) and in melting temperature (ΔT_m). We anticipate our tool to be a valuable suite for the study of protein flexibility analysis and the study of the role of variants in disease. DynaMut2 is freely available as a web server and API at http://biosig.unimelb.edu.au/dynamut2.

Identifying Genotype-Phenotype Correlations via Integrative Mutation Analysis

Mutations in protein-coding regions can lead to large biological changes and are associated with genetic conditions, including cancers and Mendelian diseases, as well as drug resistance. Although whole genome and exome sequencing help to elucidate potential genotype-phenotype correlations, there is a large gap between the identification of new variants and deciphering their molecular consequences. A comprehensive understanding of these mechanistic consequences is crucial to better understand and treat diseases in a more personalized and effective way. This is particularly relevant considering estimates that over 80% of mutations associated with a disease are incorrectly assumed to be causative. A thorough analysis of potential effects of mutations is required to correctly identify the molecular mechanisms of disease and enable the distinction between disease-causing and non-disease-causing variation within a gene. Here we present an overview of our integrative mutation analysis platform, which focuses on refining the current genotype-phenotype correlation methods by using the wealth of protein structural information.

mmCSM-AB: guiding rational antibody engineering through multiple point mutations

While antibodies are becoming an increasingly important therapeutic class, especially in personalized medicine, their development and optimization has been largely through experimental exploration. While there have been many efforts to develop computational tools to guide rational antibody engineering, most approaches are of limited accuracy when applied to antibody design, and have largely been limited to analysing a single point mutation at a time. To overcome this gap, we have curated a dataset of 242 experimentally determined changes in binding affinity upon multiple point mutations in antibody-target complexes (89 increasing and 153 decreasing binding affinity). Here, we have shown that by using our graph-based signatures and atomic interaction information, we can accurately analyse the consequence of multi-point mutations on antigen binding affinity. Our approach outperformed other available tools across cross-validation and two independent blind tests, achieving Pearson's correlations of up to 0.95. We have implemented our new approach, mmCSM-AB, as a web-server that can help guide the process of affinity maturation in antibody design. mmCSM-AB is freely available at http://biosig.unimelb.edu.au/mmcsm_ab/.

Artificial Intelligence and Machine learning based prediction of resistant and susceptible mutations in Mycobacterium tuberculosis

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis (M.tb), causes highest number of deaths globally for any bacterial disease necessitating novel diagnosis and treatment strategies. High-throughput sequencing methods generate a large amount of data which could be exploited in determining multi-drug resistant (MDR-TB) associated mutations. The present work is a computational framework that uses artificial intelligence (AI) based machine learning (ML) approaches for predicting resistance in the genes rpoB, inhA, katG, pncA, gyrA and gyrB for the drugs rifampicin, isoniazid, pyrazinamide and fluoroquinolones. The single nucleotide variations were represented by several sequence and structural features that indicate the influence of mutations on the target protein coded by each gene. We used ML algorithms - naïve bayes, k nearest neighbor, support vector machine, and artificial neural network, to build the prediction models. The classification models had an average accuracy of 85% across all examined genes and were evaluated on an external unseen dataset to demonstrate their application. Further, molecular docking and molecular dynamics simulations were performed for wild type and predicted resistance causing mutant protein and anti-TB drug complexes to study their impact on the conformation of proteins to confirm the observed phenotype.

Structure guided prediction of Pyrazinamide resistance mutations in pncA

Pyrazinamide plays an important role in tuberculosis treatment; however, its use is complicated by side-effects and challenges with reliable drug susceptibility testing. Resistance to pyrazinamide is largely driven by mutations in pyrazinamidase (pncA), responsible for drug activation, but genetic heterogeneity has hindered development of a molecular diagnostic test. We proposed to use information on how variants were likely to affect the 3D structure of pncA to identify variants likely to lead to pyrazinamide resistance. We curated 610 pncA mutations with high confidence experimental and clinical information on pyrazinamide susceptibility. The molecular consequences of each mutation on protein stability, conformation, and interactions were computationally assessed using our comprehensive suite of graph-based signature methods, mCSM. The molecular consequences of the variants were used to train a classifier with an accuracy of 80%. Our model was tested against internationally curated clinical datasets, achieving up to 85% accuracy. Screening of 600 Victorian clinical isolates identified a set of previously unreported variants, which our model had a 71% agreement with drug susceptibility testing. Here, we have shown the 3D structure of pncA can be used to accurately identify pyrazinamide resistance mutations. SUSPECT-PZA is freely available at: http://biosig.unimelb.edu.au/suspect_pza/.

A Comprehensive Computational Platform to Guide Drug Development Using Graph-Based Signature Methods

High-throughput computational techniques have become invaluable tools to help increase the overall success, process efficiency, and associated costs of drug development. By designing ligands tailored to specific protein structures in a disease of interest, an understanding of molecular interactions and ways to optimize them can be achieved prior to chemical synthesis. This understanding can help direct crucial chemical and biological experiments by maximizing available resources on higher quality leads. Moreover, predicting molecular binding affinity within specific biological contexts, as well as ligand pharmacokinetics and toxicities, can aid in filtering out redundant leads early on within the process. We describe a set of computational tools which can aid in drug discovery at different stages, from hit identification (EasyVS) to lead optimization and candidate selection (CSM-lig, mCSM-lig, Arpeggio, pkCSM). Incorporating these tools along the drug development process can help ensure that candidate leads are chemically and biologically feasible to become successful and tractable drugs.

Computational saturation mutagenesis to predict structural consequences of systematic mutations in the beta subunit of RNA polymerase in Mycobacterium leprae

Rifampin resistance in leprosy may remain undetected due to the lack of rapid and effective diagnostic tools. A quick and reliable method is essential to determine the impacts of emerging detrimental mutations in the drug targets. The functional consequences of missense mutations in the β-subunit of RNA polymerase (RNAP) in Mycobacterium leprae (M. leprae) contribute to phenotypic resistance to rifampin in leprosy. Here, we report in-silico saturation mutagenesis of all residues in the β-subunit of RNAP to all other 19 amino acid types (generating 21,394 mutations for 1126 residues) and predict their impacts on overall thermodynamic stability, on interactions at subunit interfaces, and on β-subunit-RNA and rifampin affinities (only for the rifampin binding site) using state-of-the-art structure, sequence and normal mode analysis-based methods. Mutations in the conserved residues that line the active-site cleft show largely destabilizing effects, resulting in increased relative solvent accessibility and a concomitant decrease in residue-depth (the extent to which a residue is buried in the protein structure space) of the mutant residues. The mutations at residue positions S437, G459, H451, P489, K884 and H1035 are identified as extremely detrimental as they induce highly destabilizing effects on the overall protein stability, and nucleic acid and rifampin affinities. Destabilizing effects were predicted for all the clinically/experimentally identified rifampin-resistant mutations in M. leprae indicating that this model can be used as a surveillance tool to monitor emerging detrimental mutations that destabilise RNAP-rifampin interactions and confer rifampin resistance in leprosy.

Author summary

The emergence of primary and secondary drug resistance to rifampin in leprosy is a growing concern and poses a threat to the leprosy control and elimination measures globally. In the absence of an effective in-vitro system to detect and monitor phenotypic resistance to rifampin in leprosy, diagnosis mainly relies on the presence of mutations in drug resistance determining regions of the rpoB gene that encodes the β-subunit of RNAP in M. leprae. Few labs in the world perform mouse food pad propagation of M. leprae in the presence of drugs (rifampin) to determine growth patterns and confirm resistance, however the duration of these methods lasts from 8 to 12 months making them impractical for diagnosis. Understanding molecular mechanisms of drug resistance is vital to associating mutations to clinically detected drug resistance in leprosy. Here we propose an in-silico saturation mutagenesis approach to comprehensively elucidate the structural implications of any mutations that exist or that can arise in the β-subunit of RNAP in M. leprae. Most of the predicted mutations may not occur in M. leprae due to fitness costs but the information thus generated by this approach help decipher the impacts of mutations across the structure and conversely enable identification of stable regions in the protein that are least impacted by mutations (mutation coolspots) which can be a potential choice for small molecule binding and structure guided drug discovery.

Molecular investigation against the resistant mechanism of PncA mutated pyrazinamide resistance and insight into the role of pH environment for pyrazinamide activation

Pyrazinamide (PZA), a crucial component of anti-TB therapy, is a prodrug. PZA interacts with PncA protein to be converted into its functional form i.e. pyrazinoic acid (POA). It has unique feature to kill dormant tubercle bacilli of acidic environment. Although significance of pH environment in PZA activation has been investigated in several of previous studies, insight into the significant atomistic variations in the interaction pattern of PZA with PncA, at different pH environments, are still required to be explored. On the other hand, continuously emerging PncA mutants, associated with PZA resistance, have also become a serious threat for global TB control program. Therefore, the current study was designed to understand the role of pH environment in the PZA activation and to explore the PZA resistance mechanism in various PncA mutants. The study included various in silico experiments like molecular docking, MD simulation, binding free energy estimation, PCA and FEL. In our study, we have found pH-3 and pH-5 environment as a highly significant environment for PZA activation. It was found that protonation or deprotonation of PZA activation site (PAS) residues, majorly K48, D56, K96 and E107, resulted in rearrangement of the PAS according to the pH conditions. It has also been observed that positioning of PZA binding near to Fe²⁺ and residues of catalytic triad (i.e. D8, K96 and C138) also play a very crucial role in the activation of PZA. The overall insight from the current study may help to develop new therapeutics against PncA mutated PZA resistance.Communicated by Ramaswamy H. Sarma.

Empirical ways to identify novel Bedaquiline resistance mutations in AtpE

Clinical resistance against Bedaquiline, the first new anti-tuberculosis compound with a novel mechanism of action in over 40 years, has already been detected in Mycobacterium tuberculosis. As a new drug, however, there is currently insufficient clinical data to facilitate reliable and timely identification of genomic determinants of resistance. Here we investigate the structural basis for M. tuberculosis associated bedaquiline resistance in the drug target, AtpE. Together with the 9 previously identified resistance-associated variants in AtpE, 54 non-resistance-associated mutations were identified through comparisons of bedaquiline susceptibility across 23 different mycobacterial species. Computational analysis of the structural and functional consequences of these variants revealed that resistance associated variants were mainly localized at the drug binding site, disrupting key interactions with bedaquiline leading to reduced binding affinity. This was used to train a supervised predictive algorithm, which accurately identified likely resistance mutations (93.3% accuracy). Application of this model to circulating variants present in the Asia-Pacific region suggests that current circulating variants are likely to be susceptible to bedaquiline. We have made this model freely available through a user-friendly web interface called SUSPECT-BDQ, StrUctural Susceptibility PrEdiCTion for bedaquiline (http://biosig.unimelb.edu.au/suspect_bdq/). This tool could be useful for the rapid characterization of novel clinical variants, to help guide the effective use of bedaquiline, and to minimize the spread of clinical resistance.

Synthesis and Structure-Activity relationship of 1-(5-isoquinolinesulfonyl)piperazine analogues as inhibitors of Mycobacterium tuberculosis IMPDH

Tuberculosis (TB) is a major infectious disease associated increasingly with drug resistance. Thus, new anti-tubercular agents with novel mechanisms of action are urgently required for the treatment of drug-resistant TB. In prior work, we identified compound 1 (cyclohexyl(4-(isoquinolin-5-ylsulfonyl)piperazin-1-yl)methanone) and showed that its anti-tubercular activity is attributable to inhibition of inosine-5′-monophosphate dehydrogenase (IMPDH) in Mycobacterium tuberculosis. In the present study, we explored the structure–activity relationship around compound 1 by synthesizing and evaluating the inhibitory activity of analogues against M. tuberculosis IMPDH in biochemical and whole-cell assays. X-ray crystallography was performed to elucidate the mode of binding of selected analogues to IMPDH. We establish the importance of the cyclohexyl, piperazine and isoquinoline rings for activity, and report the identification of an analogue with IMPDH-selective activity against a mutant of M. tuberculosis that is highly resistant to compound 1. We also show that the nitrogen in urea analogues is required for anti-tubercular activity and identify benzylurea derivatives as promising inhibitors that warrant further investigation.

•

Forty-eight analogues of 1-(5-isoquinolinesulfonyl)piperazine were synthesized.

•

Biochemical, whole-cell, and X-ray studies were performed to elucidate the IMPDH inhibition.

•

Piperazine and isoquinoline rings were essential for target-selective whole-cell activity.

•

Compound 47 showed improved IC₅₀ against the MtbIMPDH and maintained on-target whole-cell activity.

•

Compound 21 showed activity against IMPDH in both wild type M. tuberculosis and a resistant mutant of compound 1.

Exploring Protein Supersecondary Structure Through Changes in Protein Folding, Stability, and Flexibility

The ability to predict how mutations affect protein structure, folding, and flexibility can elucidate the molecular mechanisms leading to disruption of supersecondary structures, the emergence of phenotypes, as well guiding rational protein engineering. The advent of fast and accurate computational tools has enabled us to comprehensively explore the landscape of mutation effects on protein structures, prioritizing mutations for rational experimental validation.Here we describe the use of two complementary web-based in silico methods, DUET and DynaMut, developed to infer the effects of mutations on folding, stability, and flexibility and how they can be used to explore and interpret these effects on protein supersecondary structures.

Kinact: a computational approach for predicting activating missense mutations in protein kinases

Protein phosphorylation is tightly regulated due to its vital role in many cellular processes. While gain of function mutations leading to constitutive activation of protein kinases are known to be driver events of many cancers, the identification of these mutations has proven challenging. Here we present Kinact, a novel machine learning approach for predicting kinase activating missense mutations using information from sequence and structure. By adapting our graph-based signatures, Kinact represents both structural and sequence information, which are used as evidence to train predictive models. We show the combination of structural and sequence features significantly improved the overall accuracy compared to considering either primary or tertiary structure alone, highlighting their complementarity. Kinact achieved a precision of 87% and 94% and Area Under ROC Curve of 0.89 and 0.92 on 10-fold cross-validation, and on blind tests, respectively, outperforming well established tools (P < 0.01). We further show that Kinact performs equally well on homology models built using templates with sequence identity as low as 33%. Kinact is freely available as a user-friendly web server at http://biosig.unimelb.edu.au/kinact/.