CONTRABASS: exploiting flux constraints in genome-scale models for the detection of vulnerabilities

Motivation: Despite the fact that antimicrobial resistance is an increasing health concern, the pace of production of new drugs is slow due to the high cost and uncertain success of the process. The development of high-throughput technologies has allowed the integration of biological data into detailed genome-scale models of multiple organisms. Such models can be exploited by means of computational methods to identify system vulnerabilities such as chokepoint reactions and essential reactions. These vulnerabilities are appealing drug targets that can lead to novel drug developments. However, the current approach to compute these vulnerabilities is only based on topological data and ignores the dynamic information of the model. This can lead to misidentified drug targets. Results: This work computes flux constraints that are consistent with a certain growth rate of the modelled organism, and integrates the computed flux constraints into the model to improve the detection of vulnerabilities. By exploiting these flux constraints, we are able to obtain a directionality of the reactions of metabolism consistent with a given growth rate of the model, and consequently, a more realistic detection of vulnerabilities can be performed. Several sets of reactions that are system vulnerabilities are defined and the relationships among them are studied. The approach for the detection of these vulnerabilities has been implemented in the Python tool CONTRABASS. Such tool, for which an online web server has also been implemented, computes flux constraints and generates a report with the detected vulnerabilities

Multispecies reconstructions uncover widespread conservation, and lineage-specific elaborations in eukaryotic mRNA metabolism.

The degree of conservation and evolution of cytoplasmic mRNA metabolism pathways across the eukaryotes remains incompletely resolved. In this study, we describe a comprehensive genome and transcriptome-wide analysis of proteins involved in mRNA maturation, translation, and mRNA decay across representative organisms from the six eukaryotic super-groups. We demonstrate that eukaryotes share common pathways for mRNA metabolism that were almost certainly present in the last eukaryotic common ancestor, and show for the first time a correlation between intron density and a selective absence of some Exon Junction Complex (EJC) components in eukaryotes. In addition, we identify pathways that have diversified in individual lineages, with a specific focus on the unique gene gains and losses in members of the Excavata and SAR groups that contribute to their unique gene expression pathways compared to other organisms

Integrated human/SARS-CoV-2 metabolic models present novel treatment strategies against COVID-19.

The coronavirus disease 2019 (COVID-19) pandemic caused by the new coronavirus (SARS-CoV-2) is currently responsible for more than 3 million deaths in 219 countries across the world and with more than 140 million cases. The absence of FDA-approved drugs against SARS-CoV-2 has highlighted an urgent need to design new drugs. We developed an integrated model of the human cell and SARS-CoV-2 to provide insight into the virus' pathogenic mechanism and support current therapeutic strategies. We show the biochemical reactions required for the growth and general maintenance of the human cell, first, in its healthy state. We then demonstrate how the entry of SARS-CoV-2 into the human cell causes biochemical and structural changes, leading to a change of cell functions or cell death. A new computational method that predicts 20 unique reactions as drug targets from our models and provides a platform for future studies on viral entry inhibition, immune regulation, and drug optimisation strategies. The model is available in BioModels (https://www.ebi.ac.uk/biomodels/MODEL2007210001) and the software tool, findCPcli, that implements the computational method is available at https://github.com/findCP/findCPcli

SARS-CoV-2 3D database: Understanding the Coronavirus Proteome and Evaluating Possible Drug Targets.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a rapidly growing infectious disease, widely spread with high mortality rates. Since the release of the SARS-CoV-2 genome sequence in March 2020, there has been an international focus on developing target-based drug discovery, which also requires knowledge of the 3D structure of the proteome. Where there are no experimentally solved structures, our group has created 3D models with coverage of 97.5% and characterised them using state-of-the-art computational approaches. Models of protomers and oligomers, together with predictions of substrate and allosteric binding sites, protein- ligand docking, SARS-CoV-2 protein interactions with human proteins, impacts of mutations, and mapped solved experimental structures are freely available for download. These are imple- mented in SARS CoV-2 3D, a comprehensive and user-friendly database, available at https://sars3d.com/. This provides essential information for drug discovery, both to evaluate targets and design new potential therapeutics.This work is supported and funded by King Abdullah scholarship (Saudi Arabia research coun- cil), and American Leprosy Missions grants (G88726), SET is funded by the Cystic Fibrosis Trust (RG 70975) and Fondation Botnar (RG91317). A.R.J is funded by the Biotechnology and Biological Sciences Research Council (BBSRC) DTP studentship (BB/M011194/1). B.B. is funded by the Cystic Fibrosis Trust and L.C. on a studentship from Ipsen. T.L.B. is funded by a the Wellcome Trust Investigator Award, PHZJ/489 RG83114 (2016-2021

Mycobacterial genomics and structural bioinformatics: opportunities and challenges in drug discovery.

Of the more than 190 distinct species of Mycobacterium genus, many are economically and clinically important pathogens of humans or animals. Among those mycobacteria that infect humans, three species namely Mycobacterium tuberculosis (causative agent of tuberculosis), Mycobacterium leprae (causative agent of leprosy) and Mycobacterium abscessus (causative agent of chronic pulmonary infections) pose concern to global public health. Although antibiotics have been successfully developed to combat each of these, the emergence of drug-resistant strains is an increasing challenge for treatment and drug discovery. Here we describe the impact of the rapid expansion of genome sequencing and genome/pathway annotations that have greatly improved the progress of structure-guided drug discovery. We focus on the applications of comparative genomics, metabolomics, evolutionary bioinformatics and structural proteomics to identify potential drug targets. The opportunities and challenges for the design of drugs for M. tuberculosis, M. leprae and M. abscessus to combat resistance are discussed

L'abordatge de l'assetjament escolar des de la mediació en la Justícia Juvenil de Catalunya

<div><p>The number of paralogs of proteins involved in translation initiation is larger in trypanosomes than in yeasts or many metazoan and includes two poly(A) binding proteins, PABP1 and PABP2, and four eIF4E variants. In many cases, the paralogs are individually essential and are thus unlikely to have redundant functions although, as yet, distinct functions of different isoforms have not been determined. Here, trypanosome PABP1 and PABP2 have been further characterised. PABP1 and PABP2 diverged subsequent to the differentiation of the Kinetoplastae lineage, supporting the existence of specific aspects of translation initiation regulation. PABP1 and PABP2 exhibit major differences in intracellular localization and distribution on polysome fractionation under various conditions that interfere with mRNA metabolism. Most striking are differences in localization to the four known types of inducible RNP granules. Moreover, only PABP2 but not PABP1 can accumulate in the nucleus. Taken together, these observations indicate that PABP1 and PABP2 likely associate with distinct populations of mRNAs. The differences in localization to inducible RNP granules also apply to paralogs of components of the eIF4F complex: eIF4E1 showed similar localization pattern to PABP2, whereas the localisation of eIF4E4 and eIF4G3 resembled that of PABP1. The grouping of translation initiation as either colocalizing with PABP1 or with PABP2 can be used to complement interaction studies to further define the translation initiation complexes in kinetoplastids.</p> </div

Stepwise pathogenic evolution of Mycobacterium abscessus.

Although almost all mycobacterial species are saprophytic environmental organisms, a few, such as Mycobacterium tuberculosis, have evolved to cause transmissible human infection. By analyzing the recent emergence and spread of the environmental organism M. abscessus through the global cystic fibrosis population, we have defined key, generalizable steps involved in the pathogenic evolution of mycobacteria. We show that epigenetic modifiers, acquired through horizontal gene transfer, cause saltational increases in the pathogenic potential of specific environmental clones. Allopatric parallel evolution during chronic lung infection then promotes rapid increases in virulence through mutations in a discrete gene network; these mutations enhance growth within macrophages but impair fomite survival. As a consequence, we observe constrained pathogenic evolution while person-to-person transmission remains indirect, but postulate accelerated pathogenic adaptation once direct transmission is possible, as observed for M. tuberculosis Our findings indicate how key interventions, such as early treatment and cross-infection control, might restrict the spread of existing mycobacterial pathogens and prevent new, emergent ones

Mycobacterial metabolic model development for drug target identification

Antibiotic resistance is increasing at an alarming rate, and three related mycobacteria are sources of widespread infections in humans. According to the World Health Organization, Mycobacterium leprae, which causes leprosy, is still endemic in tropical countries; Mycobacterium tuberculosis is the second leading infectious killer worldwide after COVID-19; and Mycobacteroides abscessus, a group of non-tuberculous mycobacteria, causes lung infections and other healthcare-associated infections in humans. Due to the rise in resistance to common antibacterial drugs, it is critical that we develop alternatives to traditional treatment procedures. Furthermore, an understanding of the biochemical mechanisms underlying pathogenic evolution is important for the treatment and management of these diseases. In this study, metabolic models have been developed for two bacterial pathogens, M. leprae and My. abscessus, and a new computational tool has been used to identify potential drug targets, which are referred to as bottleneck reactions. The genes, reactions, and pathways in each of these organisms have been highlighted; the potential drug targets can be further explored as broad-spectrum antibacterials and the unique drug targets for each pathogen are significant for precision medicine initiatives. The models and associated datasets described in this paper are available in GigaDB, Biomodels, and PatMeDB repositories

Mycobacterial metabolic model development for drug target identification.

Antibiotic resistance is increasing at an alarming rate, and three related mycobacteria are sources of widespread infections in humans. According to the World Health Organization, Mycobacterium leprae, which causes leprosy, is still endemic in tropical countries; Mycobacterium tuberculosis is the second leading infectious killer worldwide after COVID-19; and Mycobacteroides abscessus, a group of non-tuberculous mycobacteria, causes lung infections and other healthcare-associated infections in humans. Due to the rise in resistance to common antibacterial drugs, it is critical that we develop alternatives to traditional treatment procedures. Furthermore, an understanding of the biochemical mechanisms underlying pathogenic evolution is important for the treatment and management of these diseases. In this study, metabolic models have been developed for two bacterial pathogens, M. leprae and My. abscessus, and a new computational tool has been used to identify potential drug targets, which are referred to as bottleneck reactions. The genes, reactions, and pathways in each of these organisms have been highlighted; the potential drug targets can be further explored as broad-spectrum antibacterials and the unique drug targets for each pathogen are significant for precision medicine initiatives. The models and associated datasets described in this paper are available in GigaDB, Biomodels, and PatMeDB repositories

CONTRABASS: exploiting flux constraints in genome-scale models for the detection of vulnerabilities.

Funder: Spanish Ministry of Science and InnovationMOTIVATION: Despite the fact that antimicrobial resistance is an increasing health concern, the pace of production of new drugs is slow due to the high cost and uncertain success of the process. The development of high-throughput technologies has allowed the integration of biological data into detailed genome-scale models of multiple organisms. Such models can be exploited by means of computational methods to identify system vulnerabilities such as chokepoint reactions and essential reactions. These vulnerabilities are appealing drug targets that can lead to novel drug developments. However, the current approach to compute these vulnerabilities is only based on topological data and ignores the dynamic information of the model. This can lead to misidentified drug targets. RESULTS: This work computes flux constraints that are consistent with a certain growth rate of the modelled organism, and integrates the computed flux constraints into the model to improve the detection of vulnerabilities. By exploiting these flux constraints, we are able to obtain a directionality of the reactions of metabolism consistent with a given growth rate of the model, and consequently, a more realistic detection of vulnerabilities can be performed. Several sets of reactions that are system vulnerabilities are defined and the relationships among them are studied. The approach for the detection of these vulnerabilities has been implemented in the Python tool CONTRABASS. Such tool, for which an online web server has also been implemented, computes flux constraints and generates a report with the detected vulnerabilities. AVAILABILITY AND IMPLEMENTATION: CONTRABASS is available as an open source Python package at https://github.com/openCONTRABASS/CONTRABASS under GPL-3.0 License. An online web server is available at http://contrabass.unizar.es. SUPPLEMENTARY INFORMATION: A glossary of terms are available at Bioinformatics online