Exploring gene knockout strategies to identify potential drug targets using genome-scale metabolic models

Semi-Automatic Detection of Errors in Genome-Scale Metabolic Models

Background

Genome-Scale Metabolic Models (GSMMs) are used for numerous tasks requiring computational estimates of metabolic fluxes, from predicting novel drug targets to engineering microbes to produce valuable compounds. A key limiting step in most applications of GSMMs is ensuring their representation of the target organism's metabolism is complete and accurate. Identifying and visualizing errors in GSMMs is complicated by the fact that they contain thousands of densely interconnected reactions. Furthermore, many errors in GSMMs only become apparent when considering pathways of connected reactions collectively, as opposed to examining reactions individually.

Results

We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a collection of algorithms for detecting errors in GSMMs. The relative frequencies of errors we detect in manually curated GSMMs appear to reflect the different approaches used to curate them. Changing the method used to automatically create a GSMM from a particular organism's genome can have a larger impact on the kinds of errors in the resulting GSMM than using the same method with a different organism's genome. Our algorithms are particularly capable of identifying errors that are only apparent at the pathway level, including loops, and nontrivial cases of dead ends.

Conclusions

MACAW is capable of identifying inaccuracies of varying severity in a wide range of GSMMs. Correcting these errors can measurably improve the predictive capacity of a GSMM. The relative prevalence of each type of error we identify in a large collection of GSMMs could help shape future efforts for further automation of error correction and GSMM creation.

Incidence of chemotherapy-related cardiac dysfunction in cancer patients

BACKGROUND

Cancer patients are increasingly affected by chemotherapy-related cardiac dysfunction. The reported incidence of this condition vary significantly across different studies.

HYPOTHESIS

A better comprehensive understanding of chemotherapy-related cardiac dysfunction incidence in cancer patients is imperative. Therefore, we performed a meta-analysis to establish the overall incidence of chemotherapy-related cardiac dysfunction in cancer patients.

METHODS

We searched articles in PubMed and EMBASE from database inception to May 1, 2023. Studies that reported the incidence of chemotherapy-related cardiac dysfunction in cancer patients were included.

RESULTS

A total of 53 studies involving 35 651 individuals were finally included in the meta-analysis. The overall pooled incidence of chemotherapy-related cardiac dysfunction in cancer patients was 63.21 per 1000 person-years (95% CI: 57.28-69.14). The chemotherapy-related cardiac dysfunction incidence increased steeply within half a year of cancer chemotherapy. Also, the trend of chemotherapy-related cardiac dysfunction incidence appeared to have plateaued after a longer duration of follow-up. In addition, chemotherapy-related cardiac dysfunction incidence rates are significantly higher among patients with age ≥50 years versus patients with age <50 years (99.96 vs. 34.48 per 1000 person-years). The incidence rate of cardiac dysfunction was higher among breast cancer patients (72.97 per 1000 person-years), leukemia patients (65.21 per 1000 person-years), and lymphoma patients (55.43 per 1000 person-years).

CONCLUSION

Our meta-analysis unveiled a definitive overall incidence rate of chemotherapy-related cardiac dysfunction in cancer patients. In addition, it was found that the risk of developing this condition escalates within the initial 6 months postchemotherapy, subsequently tapering off to become statistically insignificant after a duration of 6 years.

A reconstructed genome-scale metabolic model of Helicobacter pylori for predicting putative drug targets in clarithromycin and rifampicin resistance conditions

BACKGROUND

Helicobacter pylori is considered a true human pathogen for which rising drug resistance constitutes a drastic concern globally. The present study aimed to reconstruct a genome-scale metabolic model (GSMM) to decipher the metabolic capability of H. pylori strains in response to clarithromycin and rifampicin along with identification of novel drug targets.

MATERIALS AND METHODS

The iIT341 model of H. pylori was updated based on genome annotation data, and biochemical knowledge from literature and databases. Context-specific models were generated by integrating the transcriptomic data of clarithromycin and rifampicin resistance into the model. Flux balance analysis was employed for identifying essential genes in each strain, which were further prioritized upon being nonhomologs to humans, virulence factor analysis, druggability, and broad-spectrum analysis. Additionally, metabolic differences between sensitive and resistant strains were also investigated based on flux variability analysis and pathway enrichment analysis of transcriptomic data.

RESULTS

The reconstructed GSMM was named as HpM485 model. Pathway enrichment and flux variability analyses demonstrated reduced activity in the ribosomal pathway in both clarithromycin- and rifampicin-resistant strains. Also, a significant decrease was detected in the activity of metabolic pathways of clarithromycin-resistant strain. Moreover, 23 and 16 essential genes were exclusively detected in clarithromycin- and rifampicin-resistant strains, respectively. Based on prioritization analysis, cyclopropane fatty acid synthase and phosphoenolpyruvate synthase were identified as putative drug targets in clarithromycin- and rifampicin-resistant strains, respectively.

CONCLUSIONS

We present a robust and reliable metabolic model of H. pylori. This model can predict novel drug targets to combat drug resistance and explore the metabolic capability of H. pylori in various conditions.

PEI functionalized cell membrane for tumor targeted and glutathione responsive gene delivery

Polyethylenimine (PEI) is a broadly exploited cationic polymer due to its remarkable gene-loading capacity. However, the high cytotoxicity caused by its high surface charge density has been reported in many cell lines, limiting its application significantly. In this study, two different molecular weights of PEI (PEI10k and PEI25k) were crosslinked with red blood cell membranes (RBCm) via disulfide bonds to form PEI derivatives (RMPs) with lower charge density. Furthermore, the targeting molecule folic acid (FA) molecules were further grafted onto the polymers to obtain FA-modified PEI-RBCm copolymers (FA-RMP25k) with tumor cell targeting and glutathione response. In vitro experiments showed that the FA-RMP25k/DNA complex had satisfactory uptake efficiency in both HeLa and 293T cells, and did not cause significant cytotoxicity. Furthermore, the uptake and transfection efficiency of the FA-RMP25k/DNA complex was significantly higher than that of the PEI25k/DNA complex, indicating that FA grafting can increase transfection efficiency by 15 %. These results suggest that FA-RMP25k may be a promising non-viral gene vector with potential applications in gene therapy.

VirulentPred 2.0: An improved method for prediction of virulent proteins in bacterial pathogens

Virulence proteins in pathogens are essential for causing disease in a host. They enable the pathogen to invade, survive and multiply within the host, thus enhancing its potential to cause disease while also causing evasion of host defense mechanisms. Identifying these factors, especially potential vaccine candidates or drug targets, is critical for vaccine or drug development research. In this context, we present an improved version of VirulentPred 1.0 for rapidly identifying virulent proteins. The VirulentPred 2.0 is based on training machine learning models with experimentally validated virulent protein sequences. VirulentPred 2.0 achieved 84.71% accuracy with the validation dataset and 85.18% on an independent test dataset. The models are trained and evaluated with the latest sequence datasets of virulent proteins, which are three times greater in number than the proteins used in the earlier version of VirulentPred. Moreover, a significant improvement of 11% in the prediction accuracy over the earlier version is achieved with the best position-specific scoring matrix (PSSM)-based model for the latest test dataset. VirulentPred 2.0 is available as a user-friendly web interface at https://bioinfo.icgeb.res.in/virulent2/ and a standalone application suitable for bulk predictions. With higher efficiency and availability as a standalone tool, VirulentPred 2.0 holds immense potential for high throughput yet efficient identification of virulent proteins in bacterial pathogens.

Chemical genetic approaches for the discovery of bacterial cell wall inhibitors

Antimicrobial resistance (AMR) in bacterial pathogens is a worldwide health issue. The innovation gap in discovering new antibiotics has remained a significant hurdle in combating the AMR problem. Currently, antibiotics target various vital components of the bacterial cell envelope, nucleic acid and protein biosynthesis machinery and metabolic pathways essential for bacterial survival. The critical role of the bacterial cell envelope in cell morphogenesis and integrity makes it an attractive drug target. While a significant number of in-clinic antibiotics target peptidoglycan biosynthesis, several components of the bacterial cell envelope have been overlooked. This review focuses on various antibacterial targets in the bacterial cell wall and the strategies employed to find their novel inhibitors. This review will further elaborate on combining forward and reverse chemical genetic approaches to discover antibacterials that target the bacterial cell envelope.

Recent advancement in targeted therapy and role of emerging technologies to treat cancer

Cancer is a complex disease that causes abnormal cell growth and spread. DNA mutations, chemical or environmental exposure, viral infections, chronic inflammation, hormone abnormalities, etc., are underlying factors that can cause cancer. Drug resistance and toxicity complicate cancer treatment. Additionally, the variability of cancer makes it difficult to establish universal treatment guidelines. Next-generation sequencing has made genetic testing inexpensive. This uncovers genetic mutations that can be treated with specialty drugs. AI (artificial intelligence), machine learning, biopsy, next-generation sequencing, and digital pathology provide personalized cancer treatment. This allows for patient-specific biological targets and cancer treatment. Monoclonal antibodies, CAR-T, and cancer vaccines are promising cancer treatments. Recent trial data incorporating these therapies have shown superiority in clinical outcomes and drug tolerability over conventional chemotherapies. Combinations of these therapies with new technology can change cancer treatment and help many. This review discusses the development and challenges of targeted therapies like monoclonal antibodies (mAbs), bispecific antibodies (BsAbs), bispecific T cell engagers (BiTEs), dual variable domain (DVD) antibodies, CAR-T therapy, cancer vaccines, oncolytic viruses, lipid nanoparticle-based mRNA cancer vaccines, and their clinical outcomes in various cancers. We will also study how artificial intelligence and machine learning help find new cancer treatment targets.

SL-scan identifies synthetic lethal interactions in cancer using metabolic networks

Exploiting synthetic lethality is a promising strategy for developing targeted cancer therapies. However, identifying clinically significant synthetic lethal (SL) interactions among a large number of gene combinations is a challenging computational task. In this study, we developed the SL-scan pipeline based on metabolic network modeling to discover SL interaction. The SL-scan pipeline identifies the association between simulated Flux Balance Analysis knockout scores and mutation data across cancer cell lines and predicts putative SL interactions. We assessed the concordance of the SL pairs predicted by SL-scan with those of obtained from analysis of the CRISPR, shRNA, and PRISM datasets. Our results demonstrate that the SL-scan pipeline outperformed existing SL prediction approaches based on metabolic networks in identifying SL pairs in various cancers. This study emphasizes the importance of integrating multiple data sources, particularly mutation data, when identifying SL pairs for targeted cancer therapies. The findings of this study may lead to the development of novel targeted cancer therapies.

Latent tuberculosis and computational biology: A less-talked affair

Tuberculosis (TB) is a pervasive and devastating air-borne disease caused by the organisms belonging to the Mycobacterium tuberculosis (Mtb) complex. Currently, it is the global leader in infectious disease-related death in adults. The proclivity of TB to enter the latent state has become a significant impediment to the global effort to eradicate TB. Despite decades of research, latent tuberculosis (LTB) mechanisms remain poorly understood, making it difficult to develop efficient treatment methods. In this review, we seek to shed light on the current understanding of the mechanism of LTB, with an accentuation on the insights gained through computational biology. We have outlined various well-established computational biology components, such as omics, network-based techniques, mathematical modelling, artificial intelligence, and molecular docking, to disclose the crucial facets of LTB. Additionally, we highlighted important tools and software that may be used to conduct a variety of systems biology assessments. Finally, we conclude the article by addressing the possible future directions in this field, which might help a better understanding of LTB progression.

Genome-scale models as a vehicle for knowledge transfer from microbial to mammalian cell systems

With the plethora of omics data becoming available for mammalian cell and, increasingly, human cell systems, Genome-scale metabolic models (GEMs) have emerged as a useful tool for their organisation and analysis. The systems biology community has developed an array of tools for the solution, interrogation and customisation of GEMs as well as algorithms that enable the design of cells with desired phenotypes based on the multi-omics information contained in these models. However, these tools have largely found application in microbial cells systems, which benefit from smaller model size and ease of experimentation. Herein, we discuss the major outstanding challenges in the use of GEMs as a vehicle for accurately analysing data for mammalian cell systems and transferring methodologies that would enable their use to design strains and processes. We provide insights on the opportunities and limitations of applying GEMs to human cell systems for advancing our understanding of health and disease. We further propose their integration with data-driven tools and their enrichment with cellular functions beyond metabolism, which would, in theory, more accurately describe how resources are allocated intracellularly.

Fuzzy optimization for identifying anti-cancer targets with few side effects in constraint-based models of head and neck cancer

Computer-aided methods can be used to screen potential candidate targets and to reduce the time and cost of drug development. In most of these methods, synthetic lethality is used as a therapeutic criterion to identify drug targets. However, these methods do not consider the side effects during the identification stage. This study developed a fuzzy multi-objective optimization for identifying anti-cancer targets that not only evaluated cancer cell mortality, but also minimized side effects due to treatment. We identified potential anti-cancer enzymes and antimetabolites for the treatment of head and neck cancer (HNC). The identified one- and two-target enzymes were primarily involved in six major pathways, namely, purine and pyrimidine metabolism and the pentose phosphate pathway. Most of the identified targets can be regulated by approved drugs; thus, these drugs are potential candidates for drug repurposing as a treatment for HNC. Furthermore, we identified antimetabolites involved in pathways similar to those identified using a gene-centric approach. Moreover, HMGCR knockdown could not block the growth of HNC cells. However, the two-target combinations of (UMPS, HMGCR) and (CAD, HMGCR) could achieve cell mortality and improve metabolic deviation grades over 22% without reducing the cell viability grade.

Elucidating the metabolic characteristics of pancreatic β-cells from patients with type 2 diabetes (T2D) using a genome-scale metabolic modeling

Diabetes is a global health problem caused primarily by the inability of pancreatic β-cells to secrete adequate insulin. Despite extensive research, the identity of factors contributing to the dysregulated metabolism-secretion coupling in the β-cells remains elusive. The present study attempts to capture some of these factors responsible for the impaired β-cell metabolism-secretion coupling that contributes to diabetes pathogenesis. The metabolic-flux profiles of pancreatic β-cells were predicted using genome-scale metabolic modeling for ten diabetic patients and ten control subjects. Analysis of these flux states shows reduction in the mitochondrial fatty acid oxidation and mitochondrial oxidative phosphorylation pathways, that leads to decreased insulin secretion in diabetes. We also observed elevated reactive oxygen species (ROS) generation through peroxisomal fatty acid β-oxidation. In addition, cellular antioxidant defense systems were found to be attenuated in diabetes. Our analysis also uncovered the possible changes in the plasma metabolites in diabetes due to the β-cells failure. These efforts subsequently led to the identification of seven metabolites associated with cardiovascular disease (CVD) pathogenesis, thus establishing its link as a secondary complication of diabetes.

Recombinant COL6 α2 as a Self-Organization Factor That Triggers Orderly Nerve Regeneration Without Guidance Cues

Collagen VI (COL6) in the microenvironment was recently identified as an extracellular signal that bears the function of promoting orderly axon bundle formation. However, the large molecular weight of COL6 (≈2,000 kDa) limits its production and clinical application. It remains unclear whether the smaller subunit α chains of COL6 can exert axon bundling and ordering effects independently. Herein, based on a dorsal root ganglion (DRG) ex vivo model, the contributions of three main COL6 α chains on orderly nerve bundle formation were analyzed, and COL6 α2 showed the largest contribution weight. A recombinant COL6 α2 chain was produced and demonstrated to promote the formation of orderly axon bundles through the NCAM1-mediated pathway. The addition of COL6 α2 in conventional hydrogel triggered orderly nerve regeneration in a rat sciatic nerve defect model. Immunogenicity assessment showed weaker immunogenicity of COL6 α2 compared to that of the COL6 complex. These findings suggest that recombinant COL6 α2 is a promising material for orderly nerve regeneration.