Gene regulatory network inference and analysis of multidrug-resistant Pseudomonas aeruginosa.

BACKGROUND: Healthcare-associated infections caused by bacteria such as Pseudomonas aeruginosa are a major public health problem worldwide. Gene regulatory networks (GRN) computationally represent interactions among regulatory genes and their targets. They are an important approach to help understand bacterial behaviour and to provide novel ways of overcoming scientific challenges, including the identification of potential therapeutic targets and the development of new drugs. OBJECTIVES: The goal of this study was to reconstruct the multidrug-resistant (MDR) P. aeruginosa GRN and to analyse its topological properties. METHODS: The methodology used in this study was based on gene orthology inference using the reciprocal best hit method. We used the genome of P. aeruginosa CCBH4851 as the basis of the reconstruction process. This MDR strain is representative of the sequence type 277, which was involved in an endemic outbreak in Brazil. FINDINGS: We obtained a network with a larger number of regulatory genes, target genes and interactions as compared to the previously reported network. Topological analysis results are in accordance with the complex network representation of biological processes. MAIN CONCLUSIONS: The properties of the network were consistent with the biological features of P. aeruginosa. To the best of our knowledge, the P. aeruginosa GRN presented here is the most complete version available to date.

Detection of potential safety signals related to the use of remdesivir and tocilizumab in the COVID era during pregnancy, resorting to open data from the FDA adverse event reporting system (FAERS).

Background: The in-hospital treatment for COVID-19 may include medicines from various therapeutic classes, such as antiviral remdesivir and immunosuppressant tocilizumab. Safety data for these medicines are based on controlled clinical trials and case reports, limiting the knowledge about less frequent, rare or unique population adverse events excluded from clinical trials. Objective: This study aims at analyzing the reports of Adverse Drug Events (ADEs) related to these two medicines, focusing on events in pregnant women and foetuses. Methods: Data from the open-access FDA Adverse Event Reporting System (FAERS) from 2020 to 2022 were used to create a dashboard on the Grafana platform to ease querying and analyzing report events. Potential safety signals were generated using the ROR disproportionality measure. Results: Remdesivir was notified as the primary suspect in 7,147 reports and tocilizumab in 19,602. Three hundred and three potential safety signals were identified for remdesivir, of which six were related to pregnant women and foetuses (including abortion and foetal deaths). Tocilizumab accumulated 578 potential safety signals, and three of them were associated with this population (including neonatal death). Discussion: None of the possible signals generated for this population were found in the product labels. According to the NIH and the WHO protocols, both medicines are recommended for pregnant women hospitalized with COVID-19. Conclusion: Despite the known limitations of working with open data from spontaneous reporting systems (e.g., absence of certain clinical data, underreporting, a tendency to report severe events and recent medicines) and disproportionality analysis, the findings suggest concerning associations that need to be confirmed or rejected in subsequent clinical studies.

Optimizing therapeutic targets for breast cancer using boolean network models.

Studying gene regulatory networks associated with cancer provides valuable insights for therapeutic purposes, given that cancer is fundamentally a genetic disease. However, as the number of genes in the system increases, the complexity arising from the interconnections between network components grows exponentially. In this study, using Boolean logic to adjust the existing relationships between network components has facilitated simplifying the modeling process, enabling the generation of attractors that represent cell phenotypes based on breast cancer RNA-seq data. A key therapeutic objective is to guide cells, through targeted interventions, to transition from the current cancer attractor to a physiologically distinct attractor unrelated to cancer. To achieve this, we developed a computational method that identifies network nodes whose inhibition can facilitate the desired transition from one tumor attractor to another associated with apoptosis, leveraging transcriptomic data from cell lines. To validate the model, we utilized previously published in vitro experiments where the downregulation of specific proteins resulted in cell growth arrest and death of a breast cancer cell line. The method proposed in this manuscript combines diverse data sources, conducts structural network analysis, and incorporates relevant biological knowledge on apoptosis in cancer cells. This comprehensive approach aims to identify potential targets of significance for personalized medicine.

Boolean model of the gene regulatory network of Pseudomonas aeruginosa CCBH4851.

Introduction: Pseudomonas aeruginosa infections are one of the leading causes of death in immunocompromised patients with cystic fibrosis, diabetes, and lung diseases such as pneumonia and bronchiectasis. Furthermore, P. aeruginosa is one of the main multidrug-resistant bacteria responsible for nosocomial infections worldwide, including the multidrug-resistant CCBH4851 strain isolated in Brazil. Methods: One way to analyze their dynamic cellular behavior is through computational modeling of the gene regulatory network, which represents interactions between regulatory genes and their targets. For this purpose, Boolean models are important predictive tools to analyze these interactions. They are one of the most commonly used methods for studying complex dynamic behavior in biological systems. Results and discussion: Therefore, this research consists of building a Boolean model of the gene regulatory network of P. aeruginosa CCBH4851 using data from RNA-seq experiments. Next, the basins of attraction are estimated, as these regions and the transitions between them can help identify the attractors, representing long-term behavior in the Boolean model. The essential genes of the basins were associated with the phenotypes of the bacteria for two conditions: biofilm formation and polymyxin B treatment. Overall, the Boolean model and the analysis method proposed in this work can identify promising control actions and indicate potential therapeutic targets, which can help pinpoint new drugs and intervention strategies.

A Data Science Approach for the Identification of Molecular Signatures of Aggressive Cancers.

The main hallmarks of cancer include sustaining proliferative signaling and resisting cell death. We analyzed the genes of the WNT pathway and seven cross-linked pathways that may explain the differences in aggressiveness among cancer types. We divided six cancer types (liver, lung, stomach, kidney, prostate, and thyroid) into classes of high (H) and low (L) aggressiveness considering the TCGA data, and their correlations between Shannon entropy and 5-year overall survival (OS). Then, we used principal component analysis (PCA), a random forest classifier (RFC), and protein-protein interactions (PPI) to find the genes that correlated with aggressiveness. Using PCA, we found GRB2, CTNNB1, SKP1, CSNK2A1, PRKDC, HDAC1, YWHAZ, YWHAB, and PSMD2. Except for PSMD2, the RFC analysis showed a different list, which was CAD, PSMD14, APH1A, PSMD2, SHC1, TMEFF2, PSMD11, H2AFZ, PSMB5, and NOTCH1. Both methods use different algorithmic approaches and have different purposes, which explains the discrepancy between the two gene lists. The key genes of aggressiveness found by PCA were those that maximized the separation of H and L classes according to its third component, which represented 19% of the total variance. By contrast, RFC classified whether the RNA-seq of a tumor sample was of the H or L type. Interestingly, PPIs showed that the genes of PCA and RFC lists were connected neighbors in the PPI signaling network of WNT and cross-linked pathways.

Potential safety signals for antibacterial agents from the Brazilian national pharmacovigilance database (Vigimed/VigiFlow).

Antibacterial drugs are a widely used drug class due to the frequency of infectious diseases globally. Risks knowledge should ground these medicines' selection. Data mining in large databases is essential to identify early safety signals and to support pharmacovigilance systems. We conducted a cross-sectional study to assess adverse drug events related to antibiotics reporting between December 2018 and December 2021 in the Brazilian database (Vigimed/VigiFlow). We used the Reporting Odds Ratio (ROR) disproportionality analysis method to identify disproportionate reporting signals (SDR), referring to statistical combinations between drugs and adverse events. Vancomycin was the most reported antibiotic (n = 1,733), followed by ceftriaxone (n = 1,277) and piperacillin and tazobactam (n = 1,024). We detected 294 safety signals related to antibacterials. We identified azithromycin leading in the number of safety signals (n = 49), followed by polymyxin B (n = 25). Of these, 95 were not provided for in the drug label and had little or no reports in the medical literature. Three serious events are associated with ceftazidime and avibactam, a new drug in the Brazilian market. We also found suicide attempts as a sign associated with amoxicillin/clavulanate. Gait disturbance, a worrying event, especially in the elderly, was associated with azithromycin. Our findings may help guide further pharmacoepidemiologic studies and monitoring safety signals in pharmacovigilance.

CABGen: A Web Application for the Bioinformatic Analysis of Bacterial Genomes.

Due to recent developments in NGS technologies, genome sequencing is generating large volumes of new data containing a wealth of biological information. Understanding sequenced genomes in a biologically meaningful way and delineating their functional and metabolic landscapes is a first-level challenge. Considering the global antimicrobial resistance (AMR) problem, investments to expand surveillance and improve existing genome analysis technologies are pressing. In addition, the speed at which new genomic data is generated surpasses our capacity to analyze it with available bioinformatics methods, thus creating a need to develop new, user-friendly and comprehensive analytical tools. To this end, we propose a new web application, CABGen, developed with open-source software. CABGen allows storing, organizing, analyzing, and interpreting bioinformatics data in a friendly, scalable, easy-to-use environment and can process data from bacterial isolates of different species and origins. CABGen has three modules: Upload Sequences, Analyze Sequences, and Verify Results. Functionalities include coverage estimation, species identification, de novo genome assembly, and assembly quality, genome annotation, MLST mapping, searches for genes related to AMR, virulence, and plasmids, and detection of point mutations in specific AMR genes. Visualization tools are also available, greatly facilitating the handling of biological data. The reports include those results that are clinically relevant. To illustrate the use of CABGen, whole-genome shotgun data from 181 bacterial isolates of different species collected in 5 Brazilian regions between 2018 and 2020 were uploaded and submitted to the platform's modules.

Pulmonary Mesenchymal Stem Cells in Mild Cases of COVID-19 Are Dedicated to Proliferation; In Severe Cases, They Control Inflammation, Make Cell Dispersion, and Tissue Regeneration.

Mesenchymal stem cells (MSCs) are multipotent adult stem cells present in virtually all tissues; they have potent self-renewal capacity and differentiate into multiple cell types. For many reasons, these cells are a promising therapeutic alternative to treat patients with severe COVID-19 and pulmonary post-COVID sequelae. These cells are not only essential for tissue regeneration; they can also alter the pulmonary environment through the paracrine secretion of several mediators. They can control or promote inflammation, induce other stem cells differentiation, restrain the virus load, and much more. In this work, we performed single-cell RNA-seq data analysis of MSCs in bronchoalveolar lavage samples from control individuals and COVID-19 patients with mild and severe clinical conditions. When we compared samples from mild cases with control individuals, most genes transcriptionally upregulated in COVID-19 were involved in cell proliferation. However, a new set of genes with distinct biological functions was upregulated when we compared severely affected with mild COVID-19 patients. In this analysis, the cells upregulated genes related to cell dispersion/migration and induced the γ-activated sequence (GAS) genes, probably triggered by IFNGR1 and IFNGR2. Then, IRF-1 was upregulated, one of the GAS target genes, leading to the interferon-stimulated response (ISR) and the overexpression of many signature target genes. The MSCs also upregulated genes involved in the mesenchymal-epithelial transition, virus control, cell chemotaxis, and used the cytoplasmic RNA danger sensors RIG-1, MDA5, and PKR. In a non-comparative analysis, we observed that MSCs from severe cases do not express many NF-κB upstream receptors, such as Toll-like (TLRs) TLR-3, -7, and -8; tumor necrosis factor (TNFR1 or TNFR2), RANK, CD40, and IL-1R1. Indeed, many NF-κB inhibitors were upregulated, including PPP2CB, OPTN, NFKBIA, and FHL2, suggesting that MSCs do not play a role in the "cytokine storm" observed. Therefore, lung MSCs in COVID-19 sense immune danger and act protectively in concert with the pulmonary environment, confirming their therapeutic potential in cell-based therapy for COVID-19. The transcription of MSCs senescence markers is discussed.

Galaxy and MEAN Stack to Create a User-Friendly Workflow for the Rational Optimization of Cancer Chemotherapy.

One aspect of personalized medicine is aiming at identifying specific targets for therapy considering the gene expression profile of each patient individually. The real-world implementation of this approach is better achieved by user-friendly bioinformatics systems for healthcare professionals. In this report, we present an online platform that endows users with an interface designed using MEAN stack supported by a Galaxy pipeline. This pipeline targets connection hubs in the subnetworks formed by the interactions between the proteins of genes that are up-regulated in tumors. This strategy has been proved to be suitable for the inhibition of tumor growth and metastasis in vitro. Therefore, Perl and Python scripts were enclosed in Galaxy for translating RNA-seq data into protein targets suitable for the chemotherapy of solid tumors. Consequently, we validated the process of target diagnosis by (i) reference to subnetwork entropy, (ii) the critical value of density probability of differential gene expression, and (iii) the inhibition of the most relevant targets according to TCGA and GDC data. Finally, the most relevant targets identified by the pipeline are stored in MongoDB and can be accessed through the aforementioned internet portal designed to be compatible with mobile or small devices through Angular libraries.

A Systematic Strategy to Find Potential Therapeutic Targets for Pseudomonas aeruginosa Using Integrated Computational Models.

Pseudomonas aeruginosa is an opportunistic human pathogen that has been a constant global health problem due to its ability to cause infection at different body sites and its resistance to a broad spectrum of clinically available antibiotics. The World Health Organization classified multidrug-resistant Pseudomonas aeruginosa among the top-ranked organisms that require urgent research and development of effective therapeutic options. Several approaches have been taken to achieve these goals, but they all depend on discovering potential drug targets. The large amount of data obtained from sequencing technologies has been used to create computational models of organisms, which provide a powerful tool for better understanding their biological behavior. In the present work, we applied a method to integrate transcriptome data with genome-scale metabolic networks of Pseudomonas aeruginosa. We submitted both metabolic and integrated models to dynamic simulations and compared their performance with published in vitro growth curves. In addition, we used these models to identify potential therapeutic targets and compared the results to analyze the assumption that computational models enriched with biological measurements can provide more selective and (or) specific predictions. Our results demonstrate that dynamic simulations from integrated models result in more accurate growth curves and flux distribution more coherent with biological observations. Moreover, identifying drug targets from integrated models is more selective as the predicted genes were a subset of those found in the metabolic models. Our analysis resulted in the identification of 26 non-host homologous targets. Among them, we highlighted five top-ranked genes based on lesser conservation with the human microbiome. Overall, some of the genes identified in this work have already been proposed by different approaches and (or) are already investigated as targets to antimicrobial compounds, reinforcing the benefit of using integrated models as a starting point to selecting biologically relevant therapeutic targets.

Genetic Basis of Antimicrobial Resistant Gram-Negative Bacteria Isolated From Bloodstream in Brazil.

Multidrug-resistant microorganisms are a well-known global problem, and gram-negative bacilli are top-ranking. When these pathogens are associated with bloodstream infections (BSI), outcomes become even worse. Here we applied whole-genome sequencing to access information about clonal distribution, resistance mechanism diversity and other molecular aspects of gram-negative bacilli (GNB) isolated from bloodstream infections in Brazil. It was possible to highlight international high-risk clones circulating in the Brazilian territory, such as CC258 for Klebsiella pneumoniae, ST79 for Acinetobacter baumannii and ST233 for Pseudomonas aeruginosa. Important associations can be made such as a negative correlation between CRISPR-Cas and K. pneumoniae CC258, while the genes bla TEM, bla KPC and bla CTX-M are highly associated with this clone. Specific relationships between A. baumannii clones and bla OXA-51 variants were also observed. All P. aeruginosa ST233 isolates showed the genes bla VIM and bla OXA486. In addition, some trends could be identified, where a new P. aeruginosa MDR clone (ST3079), a novel A. baumannii clonal profile circulating in Brazil (ST848), and important resistance associations in the form of bla VIM-2 and bla IMP-56 being found together in one ST233 strain, stand out. Such findings may help to develop approaches to deal with BSI and even other nosocomial infections caused by these important GNB.

Multiparametric Color Tendency Analysis (MCTA): A Method to Analyze Several Flow Cytometry Labelings Simultaneously.

Despite the remarkable evolution of flow cytometers, fluorescent probes, and flow cytometry analysis software, most users still follow the same ways for data analysis. Conventional flow cytometry analysis relies on the creation of dot plot sequences, based on two fluorescence parameters at a time, to evidence phenotypically distinct populations. Thus, reaching conclusions about the biological characteristics of the samples is a fragmented and challenging process. We present here the MCTA (Multiparametric Color Tendency Analysis), a method for data analysis that considers multiple labelings simultaneously, extending and complementing conventional analysis. The MCTA method executes the background fluorescence exclusion, spillover compensation, and a user-defined gating strategy for subpopulation analysis. The results are then presented in conventional FSC x SSC dot plots with statistical data. For each event, the method converts each of the multiple fluorescence colors under analysis into a vector, with longer vectors being attributed to more intense labelings. Then, the MCTA generates a resultant vector, which is therefore mostly influenced by predominant labelings. The radial position of this resultant vector corresponds to a resultant color, making it easy to visualize phenotypic modulations among cellular subpopulations. Besides, it is a deterministic method that quickly assigns a resulting color to all events that obey the gating strategy, with no polymeric regions defined by the user or downsampling. The MCTA application generates a single dot plot showing all events in the FCS file, but a resultant color is attributed only to those that obey the gating strategy. Therefore, it can also help to evidence rare events or unpredicted subpopulations naturally excluded from the regions defined by the user. We believe that the MCTA method adds a new perspective over multiparametric flow cytometry analysis while evidencing modulations of molecular labeling profiles based on multiple fluorescences. Availability and implementation: The instructions for the MCTA application is freely available at https://github.com/flowcytometry/MCTA.

Modeling Basins of Attraction for Breast Cancer Using Hopfield Networks.

Cancer is a genetic disease for which traditional treatments cause harmful side effects. After two decades of genomics technological breakthroughs, personalized medicine is being used to improve treatment outcomes and mitigate side effects. In mathematical modeling, it has been proposed that cancer matches an attractor in Waddington's epigenetic landscape. The use of Hopfield networks is an attractive modeling approach because it requires neither previous biological knowledge about protein-protein interactions nor kinetic parameters. In this report, Hopfield network modeling was used to analyze bulk RNA-Seq data of paired breast tumor and control samples from 70 patients. We characterized the control and tumor attractors with respect to their size and potential energy and correlated the Euclidean distances between the tumor samples and the control attractor with their corresponding clinical data. In addition, we developed a protocol that outlines the key genes involved in tumor state stability. We found that the tumor basin of attraction is larger than that of the control and that tumor samples are associated with a more substantial negative energy than control samples, which is in agreement with previous reports. Moreover, we found a negative correlation between the Euclidean distances from tumor samples to the control attractor and patient overall survival. The ascending order of each node's density in the weight matrix and the descending order of the number of patients that have the target active only in the tumor sample were the parameters that withdrew more tumor samples from the tumor basin of attraction with fewer gene inhibitions. The combinations of therapeutic targets were specific to each patient. We performed an initial validation through simulation of trastuzumab treatment effects in HER2+ breast cancer samples. For that, we built an energy landscape composed of single-cell and bulk RNA-Seq data from trastuzumab-treated and non-treated HER2+ samples. The trajectory from the non-treated bulk sample toward the treated bulk sample was inferred through the perturbation of differentially expressed genes between these samples. Among them, we characterized key genes involved in the trastuzumab response according to the literature.

Serious Adverse Drug Reactions and Safety Signals in Children: A Nationwide Database Study.

Children are more exposed to inappropriate medicine use and its consequent harms. Spontaneous reporting of suspected Serious Adverse Drug Reactions (SADR) increases knowledge and prevention of pharmacotherapy risk. Disproportionality measures are useful to quantify unexpected safety issues associated with a given drug-event pair (signals of disproportionality). This cross-sectional study aimed to assess SADR reporting and safety signals for Brazilian children from 0-12 years old, notified between January 2008 and December 2013 from the Brazilian Surveillance Agency (Notivisa). Information from serious reports (gender and age of the patient, event description, suspected drug) was included. Disproportionality analysis based on Reporting Odds Ratios with a confidence interval of 95% was conducted to identify possible signals of disproportionate reporting (SDR). Almost 30% of 1,977 suspected SADR was related to babies (0-1-year-old). 69% of reports happened with intravenous dosage forms, and 35% of suspected SADR involved off label use according to age. Laronidase, miglustat, imipenem/cilastatin, and clofarabine were involved in six or more suspected deaths among 75 deaths reported. There were 107 SDRs, of which 16 events (15%) were not described in the product labels. There was a relatively higher number of SADRs in Brazilian children compared with studies from other countries. SDRs found, (especially drug-event pairs 'imipenen/cilastatin-pneumonia' and 'laronidase-respiratory insufficiency') should be investigated more. The reports of SADR with IV dosage forms and OL drug use suggest the need for drug research and the use of better dosage forms for children in Brazil.

Characterization of a SPM-1 metallo-beta-lactamase-producing Pseudomonas aeruginosa by comparative genomics and phenotypic analysis.

Pseudomonas aeruginosa is one of the most common pathogens related to healthcare-associated infections. The Brazilian isolate, named CCBH4851, is a multidrug-resistant clone belonging to the sequence type 277. The antimicrobial resistance mechanisms of the CCBH4851 strain are associated with the presence of the bla[Formula: see text] gene, encoding a metallo-beta-lactamase, in combination with other exogenously acquired genes. Whole-genome sequencing studies focusing on emerging pathogens are essential to identify key features of their physiology that may lead to the identification of new targets for therapy. Using both Illumina and PacBio sequencing data, we obtained a single contig representing the CCBH4851 genome with annotated features that were consistent with data reported for the species. However, comparative analysis with other Pseudomonas aeruginosa strains revealed genomic differences regarding virulence factors and regulatory proteins. In addition, we performed phenotypic assays that revealed CCBH4851 is impaired in bacterial motilities and biofilm formation. On the other hand, CCBH4851 genome contained acquired genomic islands that carry transcriptional factors, virulence and antimicrobial resistance-related genes. Presence of single nucleotide polymorphisms in the core genome, mainly those located in resistance-associated genes, suggests that these mutations may also influence the multidrug-resistant behavior of CCBH4851. Overall, characterization of Pseudomonas aeruginosa CCBH4851 complete genome revealed the presence of features that strongly relates to the virulence and antibiotic resistance profile of this important infectious agent.

Signaling Complexity Measured by Shannon Entropy and Its Application in Personalized Medicine.

Traditional approaches to cancer therapy seek common molecular targets in tumors from different patients. However, molecular profiles differ between patients, and most tumors exhibit inherent heterogeneity. Hence, imprecise targeting commonly results in side effects, reduced efficacy, and drug resistance. By contrast, personalized medicine aims to establish a molecular diagnosis specific to each patient, which is currently feasible due to the progress achieved with high-throughput technologies. In this report, we explored data from human RNA-seq and protein-protein interaction (PPI) networks using bioinformatics to investigate the relationship between tumor entropy and aggressiveness. To compare PPI subnetworks of different sizes, we calculated the Shannon entropy associated with vertex connections of differentially expressed genes comparing tumor samples with their paired control tissues. We found that the inhibition of up-regulated connectivity hubs led to a higher reduction of subnetwork entropy compared to that obtained with the inhibition of targets selected at random. Furthermore, these hubs were described to be participating in tumor processes. We also found a significant negative correlation between subnetwork entropies of tumors and the respective 5-year survival rates of the corresponding cancer types. This correlation was also observed considering patients with lung squamous cell carcinoma (LUSC) and lung adenocarcinoma (LUAD) based on the clinical data from The Cancer Genome Atlas database (TCGA). Thus, network entropy increases in parallel with tumor aggressiveness but does not correlate with PPI subnetwork size. This correlation is consistent with previous reports and allowed us to assess the number of hubs to be inhibited for therapy to be effective, in the context of precision medicine, by reference to the 100% patient survival rate 5 years after diagnosis. Large standard deviations of subnetwork entropies and variations in target numbers per patient among tumor types characterize tumor heterogeneity.

Suspected Adverse Drug Reactions Related to Breast Cancer Chemotherapy: Disproportionality Analysis of the Brazilian Spontaneous Reporting System.

Spontaneous reporting systems may generate a large volume of information in real world conditions with a relatively low cost. Disproportionality measures are useful to indicate and quantify unexpected safety issues associated with a given drug-event pair (signals of disproportionality), based upon differences compared to the background reporting frequency. This cross-sectional study (2008 to 2013) aimed to analyse the feasibility of detecting such signals in the Brazilian Pharmacovigilance Database comprising suspected adverse drug reactions related to the use of doxorubicin, cyclophosphamide, carboplatin, trastuzumab, docetaxel, and paclitaxel for breast cancer chemotherapy. We first accessed overall database features (patient information and suspected adverse drug reactions) and further conducted a disproportionality analysis based on Reporting Odds Ratios with a confidence interval of 95% in order to identify possible signals of disproportionate reporting, only among serious suspected adverse drug reactions. Of all data reports of adverse reactions (n = 2603), 83% were classified as serious, with the highest prevalence with docetaxel (78.1%). The final analysis was performed using 1,309 reports with 3,139 drug-reaction pairs. The following signals of disproportionate reporting, some rare or not mentioned on labels, were observed: tachypnea with docetaxel; bronchospasm, syncope, cyanosis, and anaphylactic reaction with paclitaxel; and anaphylactic shock with trastuzumab. Structured management of spontaneous adverse drug reaction reporting is essential for monitoring the safe use of drugs and detecting early safety signals. Disproportionality signal analysis represents a viable and applicable strategy for oncology signal screening in the Brazilian Pharmacovigilance Database.

FindTargetsWEB: A User-Friendly Tool for Identification of Potential Therapeutic Targets in Metabolic Networks of Bacteria.

Background: Healthcare-associated infections (HAIs) are a serious public health problem. They can be associated with morbidity and mortality and are responsible for the increase in patient hospitalization. Antimicrobial resistance among pathogens causing HAI has increased at alarming levels. In this paper, a robust method for analyzing genome-scale metabolic networks of bacteria is proposed in order to identify potential therapeutic targets, along with its corresponding web implementation, dubbed FindTargetsWEB. The proposed method assumes that every metabolic network presents fragile genes whose blockade will impair one or more metabolic functions, such as biomass accumulation. FindTargetsWEB automates the process of identification of such fragile genes using flux balance analysis (FBA), flux variability analysis (FVA), extended Systems Biology Markup Language (SBML) file parsing, and queries to three public repositories, i.e., KEGG, UniProt, and DrugBank. The web application was developed in Python using COBRApy and Django. Results: The proposed method was demonstrated to be robust enough to process even non-curated, incomplete, or imprecise metabolic networks, in addition to integrated host-pathogen models. A list of potential therapeutic targets and their putative inhibitors was generated as a result of the analysis of Pseudomonas aeruginosa metabolic networks available in the literature and a curated version of the metabolic network of a multidrug-resistant P. aeruginosa strain belonging to a clone endemic in Brazil (P. aeruginosa ST277). Genome-scale metabolic networks of other gram-positive and gram-negative bacteria, such as Staphylococcus aureus, Klebsiella pneumoniae, and Haemophilus influenzae, were also analyzed using FindTargetsWEB. Multiple potential targets have been found using the proposed method in all metabolic networks, including some overlapping between two or more pathogens. Among the potential targets, several have been previously reported in the literature as targets for antimicrobial development, and many targets have approved drugs. Despite similarities in the metabolic network structure for closely related bacteria, we show that the method is able to selectively identify targets in pathogenic versus non-pathogenic organisms. Conclusions: This new computational system can give insights into the identification of new candidate therapeutic targets for pathogenic bacteria and discovery of new antimicrobial drugs through genome-scale metabolic network analysis and heterogeneous data integration, even for non-curated or incomplete networks.

Gene regulatory network inference and analysis of multidrug-resistant Pseudomonas aeruginosa

BACKGROUND Healthcare-associated infections caused by bacteria such as Pseudomonas aeruginosa are a major public health problem worldwide. Gene regulatory networks (GRN) computationally represent interactions among regulatory genes and their targets. They are an important approach to help understand bacterial behaviour and to provide novel ways of overcoming scientific challenges, including the identification of potential therapeutic targets and the development of new drugs. OBJECTIVES The goal of this study was to reconstruct the multidrug-resistant (MDR) P. aeruginosa GRN and to analyse its topological properties. METHODS The methodology used in this study was based on gene orthology inference using the reciprocal best hit method. We used the genome of P. aeruginosa CCBH4851 as the basis of the reconstruction process. This MDR strain is representative of the sequence type 277, which was involved in an endemic outbreak in Brazil. FINDINGS We obtained a network with a larger number of regulatory genes, target genes and interactions as compared to the previously reported network. Topological analysis results are in accordance with the complex network representation of biological processes. MAIN CONCLUSIONS The properties of the network were consistent with the biological features of P. aeruginosa. To the best of our knowledge, the P. aeruginosa GRN presented here is the most complete version available to date.