Transfer learning predicts species-specific drug interactions in emerging pathogens.

Machine learning (ML) algorithms are necessary to efficiently identify potent drug combinations within a large candidate space to combat drug resistance. However, existing ML approaches cannot be applied to emerging and under-studied pathogens with limited training data. To address this, we developed a transfer learning and crowdsourcing framework (TACTIC) to train ML models on data from multiple bacteria. TACTIC was built using 2,965 drug interactions from 12 bacterial strains and outperformed traditional ML models in predicting drug interaction outcomes for species that lack training data. Top TACTIC model features revealed genetic and metabolic factors that influence cross- species and species-specific drug interaction outcomes. Upon analyzing ∼600,000 predicted drug interactions across 9 metabolic environments and 18 bacterial strains, we identified a small set of drug interactions that are selectively synergistic against Gram- negative (e.g., A. baumannii ) and non-tuberculous mycobacteria (NTM) pathogens. We experimentally validated synergistic drug combinations containing clarithromycin, ampicillin, and mecillinam against M. abscessus , an emerging pathogen with growing levels of antibiotic resistance. Lastly, we leveraged TACTIC to propose selectively synergistic drug combinations to treat bacterial eye infections (endophthalmitis).

Integrative analysis of multimodal patient data identifies personalized predictors of tuberculosis treatment prognosis.

Tuberculosis (TB) afflicted 10.6 million people in 2021, and its global burden is increasing due to multidrug-resistant TB (MDR-TB) and extensively resistant TB (XDR-TB). Here, we analyze multi-domain information from 5,060 TB patients spanning 10 countries with high burden of MDR-TB from the NIAID TB Portals database to determine predictors of TB treatment outcome. Our analysis revealed significant associations between radiological, microbiological, therapeutic, and demographic data modalities. Our machine learning model, built with 203 features across modalities outperforms models built using each modality alone in predicting treatment outcomes, with an accuracy of 83% and area under the curve of 0.84. Notably, our analysis revealed that the drug regimens Bedaquiline-Clofazimine-Cycloserine-Levofloxacin-Linezolid and Bedaquiline-Clofazimine-Linezolid-Moxifloxacin were associated with treatment success and failure, respectively, for MDR non-XDR-TB. Drug combinations predicted to be synergistic by the INDIGO algorithm performed better than antagonistic combinations. Our prioritized set of features predictive of treatment outcomes can ultimately guide the personalized clinical management of TB.

A flux-based machine learning model to simulate the impact of pathogen metabolic heterogeneity on drug interactions.

Drug combinations are a promising strategy to counter antibiotic resistance. However, current experimental and computational approaches do not account for the entire complexity involved in combination therapy design, such as the effect of pathogen metabolic heterogeneity, changes in the growth environment, drug treatment order, and time interval. To address these limitations, we present a comprehensive approach that uses genome-scale metabolic modeling and machine learning to guide combination therapy design. Our mechanistic approach (a) accommodates diverse data types, (b) accounts for time- and order-specific interactions, and (c) accurately predicts drug interactions in various growth conditions and their robustness to pathogen metabolic heterogeneity. Our approach achieved high accuracy (area under the receiver operating curve (AUROC) = 0.83 for synergy, AUROC = 0.98 for antagonism) in predicting drug interactions for Escherichia coli cultured in 57 metabolic conditions based on experimental validation. The entropy in bacterial metabolic response was predictive of combination therapy outcomes across time scales and growth conditions. Simulation of metabolic heterogeneity using population FBA identified two subpopulations of E. coli cells defined by the levels of three proteins (eno, fadB, and fabD) in glycolysis and lipid metabolism that influence cell tolerance to a broad range of antibiotic combinations. Analysis of the vast landscape of condition-specific drug interactions revealed a set of 24 robustly synergistic drug combinations with potential for clinical use.

Integrative analysis of clinical health records, imaging and pathogen genomics identifies personalized predictors of disease prognosis in tuberculosis.

Tuberculosis (TB) afflicts over 10 million people every year and its global burden is projected to increase dramatically due to multidrug-resistant TB (MDR-TB). The Covid-19 pandemic has resulted in reduced access to TB diagnosis and treatment, reversing decades of progress in disease management globally. It is thus crucial to analyze real-world multi-domain information from patient health records to determine personalized predictors of TB treatment outcome and drug resistance. We conduct a retrospective analysis on electronic health records of 5060 TB patients spanning 10 countries with high burden of MDR-TB including Ukraine, Moldova, Belarus and India available on the NIAID-TB portals database. We analyze over 200 features across multiple host and pathogen modalities representing patient social demographics, disease presentations as seen in cChest X rays and CT scans, and genomic records with drug susceptibility features of the pathogen strain from each patient. Our machine learning model, built with diverse data modalities outperforms models built using each modality alone in predicting treatment outcomes, with an accuracy of 81% and AUC of 0.768. We determine robust predictors across countries that are associated with unsuccessful treatmentclinical outcomes, and validate our predictions on new patient data from TB Portals. Our analysis of drug regimens and drug interactions suggests that synergistic drug combinations and those containing the drugs Bedaquiline, Levofloxacin, Clofazimine and Amoxicillin see more success in treating MDR and XDR TB. Features identified via chest imaging such as percentage of abnormal volume, size of lung cavitation and bronchial obstruction are associated significantly with pathogen genomic attributes of drug resistance. Increased disease severity was also observed in patients with lower BMI and with comorbidities. Our integrated multi-modal analysis thus revealed significant associations between radiological, microbiological, therapeutic, and demographic data modalities, providing a deeper understanding of personalized responses to aid in the clinical management of TB.

Acetyl-CoA metabolism drives epigenome change and contributes to carcinogenesis risk in fatty liver disease.

BACKGROUND: The incidence of non-alcoholic fatty liver disease (NAFLD)-associated hepatocellular carcinoma (HCC) is increasing worldwide, but the steps in precancerous hepatocytes which lead to HCC driver mutations are not well understood. Here we provide evidence that metabolically driven histone hyperacetylation in steatotic hepatocytes can increase DNA damage to initiate carcinogenesis. METHODS: Global epigenetic state was assessed in liver samples from high-fat diet or high-fructose diet rodent models, as well as in cultured immortalized human hepatocytes (IHH cells). The mechanisms linking steatosis, histone acetylation and DNA damage were investigated by computational metabolic modelling as well as through manipulation of IHH cells with metabolic and epigenetic inhibitors. Chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) and transcriptome (RNA-seq) analyses were performed on IHH cells. Mutation locations and patterns were compared between the IHH cell model and genome sequence data from preneoplastic fatty liver samples from patients with alcohol-related liver disease and NAFLD. RESULTS: Genome-wide histone acetylation was increased in steatotic livers of rodents fed high-fructose or high-fat diet. In vitro, steatosis relaxed chromatin and increased DNA damage marker γH2AX, which was reversed by inhibiting acetyl-CoA production. Steatosis-associated acetylation and γH2AX were enriched at gene clusters in telomere-proximal regions which contained HCC tumour suppressors in hepatocytes and human fatty livers. Regions of metabolically driven epigenetic change also had increased levels of DNA mutation in non-cancerous tissue from NAFLD and alcohol-related liver disease patients. Finally, genome-scale network modelling indicated that redox balance could be a key contributor to this mechanism. CONCLUSIONS: Abnormal histone hyperacetylation facilitates DNA damage in steatotic hepatocytes and is a potential initiating event in hepatocellular carcinogenesis.

Lipid exposure activates gene expression changes associated with estrogen receptor negative breast cancer.

Improved understanding of local breast biology that favors the development of estrogen receptor negative (ER-) breast cancer (BC) would foster better prevention strategies. We have previously shown that overexpression of specific lipid metabolism genes is associated with the development of ER- BC. We now report results of exposure of MCF-10A and MCF-12A cells, and mammary organoids to representative medium- and long-chain polyunsaturated fatty acids. This exposure caused a dynamic and profound change in gene expression, accompanied by changes in chromatin packing density, chromatin accessibility, and histone posttranslational modifications (PTMs). We identified 38 metabolic reactions that showed significantly increased activity, including reactions related to one-carbon metabolism. Among these reactions are those that produce S-adenosyl-L-methionine for histone PTMs. Utilizing both an in-vitro model and samples from women at high risk for ER- BC, we show that lipid exposure engenders gene expression, signaling pathway activation, and histone marks associated with the development of ER- BC.

Machine learning to design antimicrobial combination therapies: Promises and pitfalls.

Combination therapies can overcome antimicrobial resistance (AMR) and repurpose existing drugs. However, the large combinatorial space to explore presents a daunting challenge. In response, machine learning (ML) algorithms are being applied to identify novel synergistic drug interactions from millions of potential combinations. Here, we compare ML-based approaches for combination therapy design based on the type of input information used, specifically: drug properties, microbial response and infection microenvironment. We also provide a compilation of publicly available drug interaction datasets relevant to AMR. Finally, we discuss limitations of current ML-based methods and propose new strategies for designing efficacious combination therapies. These include consideration of in vivo conditions, design of sequential combinations, enhancement of model interpretability and application of deep learning algorithms.

Next-Generation Genome-Scale Metabolic Modeling through Integration of Regulatory Mechanisms.

Genome-scale metabolic models (GEMs) are powerful tools for understanding metabolism from a systems-level perspective. However, GEMs in their most basic form fail to account for cellular regulation. A diverse set of mechanisms regulate cellular metabolism, enabling organisms to respond to a wide range of conditions. This limitation of GEMs has prompted the development of new methods to integrate regulatory mechanisms, thereby enhancing the predictive capabilities and broadening the scope of GEMs. Here, we cover integrative models encompassing six types of regulatory mechanisms: transcriptional regulatory networks (TRNs), post-translational modifications (PTMs), epigenetics, protein-protein interactions and protein stability (PPIs/PS), allostery, and signaling networks. We discuss 22 integrative GEM modeling methods and how these have been used to simulate metabolic regulation during normal and pathological conditions. While these advances have been remarkable, there remains a need for comprehensive and widespread integration of regulatory constraints into GEMs. We conclude by discussing challenges in constructing GEMs with regulation and highlight areas that need to be addressed for the successful modeling of metabolic regulation. Next-generation integrative GEMs that incorporate multiple regulatory mechanisms and their crosstalk will be invaluable for discovering cell-type and disease-specific metabolic control mechanisms.

The deacylase SIRT5 supports melanoma viability by influencing chromatin dynamics.

Cutaneous melanoma remains the most lethal skin cancer, and ranks third among all malignancies in terms of years of life lost. Despite the advent of immune checkpoint and targeted therapies, only roughly half of patients with advanced melanoma achieve a durable remission. Sirtuin 5 (SIRT5) is a member of the sirtuin family of protein deacylases that regulates metabolism and other biological processes. Germline Sirt5 deficiency is associated with mild phenotypes in mice. Here we showed that SIRT5 was required for proliferation and survival across all cutaneous melanoma genotypes tested, as well as uveal melanoma, a genetically distinct melanoma subtype that arises in the eye and is incurable once metastatic. Likewise, SIRT5 was required for efficient tumor formation by melanoma xenografts and in an autochthonous mouse Braf Pten-driven melanoma model. Via metabolite and transcriptomic analyses, we found that SIRT5 was required to maintain histone acetylation and methylation levels in melanoma cells, thereby promoting proper gene expression. SIRT5-dependent genes notably included MITF, a key lineage-specific survival oncogene in melanoma, and the c-MYC proto-oncogene. SIRT5 may represent a druggable genotype-independent addiction in melanoma.

Dissecting Murine Muscle Stem Cell Aging through Regeneration Using Integrative Genomic Analysis.

During aging, there is a progressive loss of volume and function in skeletal muscle that impacts mobility and quality of life. The repair of skeletal muscle is regulated by tissue-resident stem cells called satellite cells (or muscle stem cells [MuSCs]), but in aging, MuSCs decrease in numbers and regenerative capacity. The transcriptional networks and epigenetic changes that confer diminished regenerative function in MuSCs as a result of natural aging are only partially understood. Herein, we use an integrative genomics approach to profile MuSCs from young and aged animals before and after injury. Integration of these datasets reveals aging impacts multiple regulatory changes through significant differences in gene expression, metabolic flux, chromatin accessibility, and patterns of transcription factor (TF) binding activities. Collectively, these datasets facilitate a deeper understanding of the regulation tissue-resident stem cells use during aging and healing.