The Genotype and Phenotypes in Families (GPF) platform manages the large and complex data at SFARI

The exploration of genotypic variants impacting phenotypes is a cornerstone in genetics research. The emergence of vast collections containing deeply genotyped and phenotyped families has made it possible to pursue the search for variants associated with complex diseases. However, managing these large-scale datasets requires specialized computational tools tailored to organize and analyze the extensive data. GPF (Genotypes and Phenotypes in Families) is an open-source platform ( https://github.com/iossifovlab/gpf ) that manages genotypes and phenotypes derived from collections of families. The GPF interface allows interactive exploration of genetic variants, enrichment analysis for de novo mutations, and phenotype/genotype association tools. In addition, GPF allows researchers to share their data securely with the broader scientific community. GPF is used to disseminate two large-scale family collection datasets (SSC, SPARK) for the study of autism funded by the SFARI foundation. However, GPF is versatile and can manage genotypic data from other small or large family collections. Our GPF-SFARI GPF instance ( https://gpf.sfari.org/ ) provides protected access to comprehensive genotypic and phenotypic data for the SSC and SPARK. In addition, GPF-SFARI provides public access to an extensive collection of de novo mutations identified in individuals with autism and related disorders and to gene-level statistics of the protected datasets characterizing the genes' roles in autism. Here, we highlight the primary features of GPF within the context of GPF-SFARI.

BET inhibition induces vulnerability to MCL1 targeting through upregulation of fatty acid synthesis pathway in breast cancer

Therapeutic options for treatment of basal-like breast cancers remain limited. Here, we demonstrate that bromodomain and extra-terminal (BET) inhibition induces an adaptive response leading to MCL1 protein-driven evasion of apoptosis in breast cancer cells. Consequently, co-targeting MCL1 and BET is highly synergistic in breast cancer models. The mechanism of adaptive response to BET inhibition involves the upregulation of lipid synthesis enzymes including the rate-limiting stearoyl-coenzyme A (CoA) desaturase. Changes in lipid synthesis pathway are associated with increases in cell motility and membrane fluidity as well as re-localization and activation of HER2/EGFR. In turn, the HER2/EGFR signaling results in the accumulation of and vulnerability to the inhibition of MCL1. Drug response and genomics analyses reveal that MCL1 copy-number alterations are associated with effective BET and MCL1 co-targeting. The high frequency of MCL1 chromosomal amplifications (>30%) in basal-like breast cancers suggests that BET and MCL1 co-targeting may have therapeutic utility in this aggressive subtype of breast cancer.

Yan et al. show that pharmacological co-targeting of MCL1 and BET is highly effective in breast cancer cells. The proposed combination therapy may be effective for treatment of patients with aggressive subtypes of breast cancers whose tumors carry genetic aberrations associated with cell-death evasion.

A novel, multitargeted endogenous metabolic modulator composition impacts metabolism, inflammation, and fibrosis in nonalcoholic steatohepatitis-relevant primary human cell models

Nonalcoholic steatohepatitis (NASH) is a complex metabolic disease of heterogeneous and multifactorial pathogenesis that may benefit from coordinated multitargeted interventions. Endogenous metabolic modulators (EMMs) encompass a broad set of molecular families, including amino acids and related metabolites and precursors. EMMs often serve as master regulators and signaling agents for metabolic pathways throughout the body and hold the potential to impact a complex metabolic disease like NASH by targeting a multitude of pathologically relevant biologies. Here, we describe a study of a novel EMM composition comprising five amino acids and an amino acid derivative (Leucine, Isoleucine, Valine, Arginine, Glutamine, and N-acetylcysteine [LIVRQNac]) and its systematic evaluation across multiple NASH-relevant primary human cell model systems, including hepatocytes, macrophages, and stellate cells. In these model systems, LIVRQNac consistently and simultaneously impacted biology associated with all three core pathophysiological features of NASH—metabolic, inflammatory, and fibrotic. Importantly, it was observed that while the individual constituent amino acids in LIVRQNac can impact specific NASH-related phenotypes in select cell systems, the complete combination was necessary to impact the range of disease-associated drivers examined. These findings highlight the potential of specific and potent multitargeted amino acid combinations for the treatment of NASH.

Antibiotic susceptibility signatures identify potential antimicrobial targets in the Acinetobacter baumannii cell envelope

A unique, protective cell envelope contributes to the broad drug resistance of the nosocomial pathogen Acinetobacter baumannii. Here we use transposon insertion sequencing to identify A. baumannii mutants displaying altered susceptibility to a panel of diverse antibiotics. By examining mutants with antibiotic susceptibility profiles that parallel mutations in characterized genes, we infer the function of multiple uncharacterized envelope proteins, some of which have roles in cell division or cell elongation. Remarkably, mutations affecting a predicted cell wall hydrolase lead to alterations in lipooligosaccharide synthesis. In addition, the analysis of altered susceptibility signatures and antibiotic-induced morphology patterns allows us to predict drug synergies; for example, certain beta-lactams appear to work cooperatively due to their preferential targeting of specific cell wall assembly machineries. Our results indicate that the pathogen may be effectively inhibited by the combined targeting of multiple pathways critical for envelope growth.

Guided screen for synergistic three-drug combinations

Combinations of three or more drugs are routinely used in various medical fields such as clinical oncology and infectious diseases to prevent resistance or to achieve synergistic therapeutic benefits. The very large number of possible high-order drug combinations presents a formidable challenge for discovering synergistic drug combinations. Here, we establish a guided screen to discover synergistic three-drug combinations. Using traditional checkerboard and recently developed diagonal methods, we experimentally measured all pairwise interactions among eight compounds in Erwinia amylovora, the causative agent of fire blight. Showing that synergy measurements of these two methods agree, we predicted synergy/antagonism scores for all possible three-drug combinations by averaging the synergy scores of pairwise interactions. We validated these predictions by experimentally measuring 35 three-drug interactions. Therefore, our guided screen for discovering three-drug synergies is (i) experimental screen of all pairwise interactions using diagonal method, (ii) averaging pairwise scores among components to predict three-drug interaction scores, (iii) experimental testing of top predictions. In our study, this strategy resulted in a five-fold reduction in screen size to find the most synergistic three-drug combinations.

Design of high-order antibiotic combinations against M. tuberculosis by ranking and exclusion

Combinations of more than two drugs are routinely used for the treatment of pathogens and tumors. High-order combinations may be chosen due to their non-overlapping resistance mechanisms or for favorable drug interactions. Synergistic/antagonistic interactions occur when the combination has a higher/lower effect than the sum of individual drug effects. The standard treatment of Mycobacterium tuberculosis (Mtb) is an additive cocktail of three drugs which have different targets. Herein, we experimentally measured all 190 pairwise interactions among 20 antibiotics against Mtb growth. We used the pairwise interaction data to rank all possible high-order combinations by strength of synergy/antagonism. We used drug interaction profile correlation as a proxy for drug similarity to establish exclusion criteria for ideal combination therapies. Using this ranking and exclusion design (R/ED) framework, we modeled ways to improve the standard 3-drug combination with the addition of new drugs. We applied this framework to find the best 4-drug combinations against drug-resistant Mtb by adding new exclusion criteria to R/ED. Finally, we modeled alternating 2-order combinations as a cycling treatment and found optimized regimens significantly reduced the overall effective dose. R/ED provides an adaptable framework for the design of high-order drug combinations against any pathogen or tumor.

Miniaturized Checkerboard Assays to Measure Antibiotic Interactions

Drugs may have synergistic or antagonistic interactions when combined. Checkerboard assays, where two drugs are combined in many doses, allow sensitive measurement of drug interactions. Here, we describe a protocol to measure the pairwise interactions among three antibiotics, in duplicate, in 5 days, using only two 96-well microplates and standard laboratory equipment.

Characterizing ABC-Transporter Substrate-Likeness Using a Clean-Slate Genetic Background

Mutations in ATP Binding Cassette (ABC)-transporter genes can have major effects on the bioavailability and toxicity of the drugs that are ABC-transporter substrates. Consequently, methods to predict if a drug is an ABC-transporter substrate are useful for drug development. Such methods traditionally relied on literature curated collections of ABC-transporter dependent membrane transfer assays. Here, we used a single large-scale dataset of 376 drugs with relative efficacy on an engineered yeast strain with all ABC-transporter genes deleted (ABC-16), to explore the relationship between a drug’s chemical structure and ABC-transporter substrate-likeness. We represented a drug’s chemical structure by an array of substructure keys and explored several machine learning methods to predict the drug’s efficacy in an ABC-16 yeast strain. Gradient-Boosted Random Forest models outperformed all other methods with an AUC of 0.723. We prospectively validated the model using new experimental data and found significant agreement with predictions. Our analysis expands the previously reported chemical substructures associated with ABC-transporter substrates and provides an alternative means to investigate ABC-transporter substrate-likeness.

Prediction of ultra-high-order antibiotic combinations based on pairwise interactions

Drug combinations are a promising approach to achieve high efficacy at low doses and to overcome resistance. Drug combinations are especially useful when drugs cannot achieve effectiveness at tolerable doses, as occurs in cancer and tuberculosis (TB). However, discovery of effective drug combinations faces the challenge of combinatorial explosion, in which the number of possible combinations increases exponentially with the number of drugs and doses. A recent advance, called the dose model, uses a mathematical formula to overcome combinatorial explosion by reducing the problem to a feasible quadratic one: using data on drug pairs at a few doses, the dose model accurately predicts the effect of combinations of three and four drugs at all doses. The dose model has not yet been tested on higher-order combinations beyond four drugs. To address this, we measured the effect of combinations of up to ten antibiotics on E. coli growth, and of up to five tuberculosis (TB) drugs on the growth of M. tuberculosis. We find that the dose model accurately predicts the effect of these higher-order combinations, including cases of strong synergy and antagonism. This study supports the view that the interactions between drug pairs carries key information that largely determines higher-order interactions. Therefore, systematic study of pairwise drug interactions is a compelling strategy to prioritize drug regimens in high-dimensional spaces.

Chemogenomic model identifies synergistic drug combinations robust to the pathogen microenvironment

Antibiotics need to be effective in diverse environments in vivo. However, the pathogen microenvironment can have a significant impact on antibiotic potency. Further, antibiotics are increasingly used in combinations to combat resistance, yet, the effect of microenvironments on drug-combination efficacy is unknown. To exhaustively explore the impact of diverse microenvironments on drug-combinations, here we develop a computational framework—Metabolism And GENomics-based Tailoring of Antibiotic regimens (MAGENTA). MAGENTA uses chemogenomic profiles of individual drugs and metabolic perturbations to predict synergistic or antagonistic drug-interactions in different microenvironments. We uncovered antibiotic combinations with robust synergy across nine distinct environments against both E. coli and A. baumannii by searching through 2556 drug-combinations of 72 drugs. MAGENTA also accurately predicted the change in efficacy of bacteriostatic and bactericidal drug-combinations during growth in glycerol media, which we confirmed experimentally in both microbes. Our approach identified genes in glycolysis and glyoxylate pathway as top predictors of synergy and antagonism respectively. Our systems approach enables tailoring of antibiotic therapies based on the pathogen microenvironment.

The antibiotic resistance epidemic has created a pressing need to understand factors that influence antibiotic efficacy. An often-overlooked factor in the search for new treatments is the pathogen environment. Understanding the differences in pathogen sensitivity to antibiotics in lab conditions versus inside the host is necessary for translating new discoveries into the clinic. Hence, we experimentally measured the sensitivity of E. coli to drugs and drug combinations in different metabolic conditions. Our data revealed that the environment dramatically changes treatment potency. Each antibiotic class was affected uniquely by each metabolic condition. The large number of metabolic conditions inside the host greatly complicates the identification of effective therapies. To address this challenge, we present a computational approach called MAGENTA that accurately predicted efficacy of antibiotic regimens in different conditions, which we confirmed experimentally. Furthermore, we show that MAGENTA can be applied to other bacterial pathogens such as A. baumannii and M. tuberculosis without the need for generating expensive data in each organism. MAGENTA accurately predicted efficacy in the pathogen A. baumannii using data from E. coli by identifying genes that are common between the two bacteria. Our study revealed the significant yet predictable impact of environment on drug combination potency.

Modeling the impact of drug interactions on therapeutic selectivity

Combination therapies that produce synergistic growth inhibition are widely sought after for the pharmacotherapy of many pathological conditions. Therapeutic selectivity, however, depends on the difference between potency on disease-causing cells and potency on non-target cell types that cause toxic side effects. Here, we examine a model system of antimicrobial compound combinations applied to two highly diverged yeast species. We find that even though the drug interactions correlate between the two species, cell-type-specific differences in drug interactions are common and can dramatically alter the selectivity of compounds when applied in combination vs. single-drug activity-enhancing, diminishing, or inverting therapeutic windows. This study identifies drug combinations with enhanced cell-type-selectivity with a range of interaction types, which we experimentally validate using multiplexed drug-interaction assays for heterogeneous cell cultures. This analysis presents a model framework for evaluating drug combinations with increased efficacy and selectivity against pathogens or tumors.

Diagonal Method to Measure Synergy Among Any Number of Drugs

A synergistic drug combination has a higher efficacy compared to the effects of individual drugs. Checkerboard assays, where drugs are combined in many doses, allow sensitive measurement of drug interactions. However, these assays are costly and do not scale well for measuring interaction among many drugs. Several recent studies have reported drug interaction measurements using a diagonal sampling of the traditional checkerboard assay. This alternative methodology greatly decreases the cost of drug interaction experiments and allows interaction measurement for combinations with many drugs. Here, we describe a protocol to measure the three pairwise interactions and one three-way interaction among three antibiotics in duplicate, in five days, using only three 96-well microplates and standard laboratory equipment. We present representative results showing that the three-antibiotic combination of Levofloxacin + Nalidixic Acid + Penicillin G is synergistic. Our protocol scales up to measure interactions among many drugs and in other biological contexts, allowing for efficient screens for multi-drug synergies against pathogens and tumors.

A Parallel Adder Coordinates Mycobacterial Cell-Cycle Progression and Cell-Size Homeostasis in the Context of Asymmetric Growth and Organization

In model bacteria, such as E. coli and B. subtilis, regulation of cell-cycle progression and cellular organization achieves consistency in cell size, replication dynamics, and chromosome positioning [1-3]. Mycobacteria elongate and divide asymmetrically, giving rise to significant variation in cell size and elongation rate among closely related cells [4, 5]. Given the physical asymmetry of mycobacteria, the models that describe coordination of cellular organization and cell-cycle progression in model bacteria are not directly translatable [1, 2, 6-8]. Here, we used time-lapse microscopy and fluorescent reporters of DNA replication and chromosome positioning to examine the coordination of growth, division, and chromosome dynamics at a single-cell level in Mycobacterium smegmatis (M. smegmatis) and Mycobacterium bovis Bacillus Calmette-Guérin (BCG). By analyzing chromosome and replisome localization, we demonstrated that chromosome positioning is asymmetric and proportional to cell size. Furthermore, we found that cellular asymmetry is maintained throughout the cell cycle and is not established at division. Using measurements and stochastic modeling of mycobacterial cell size and cell-cycle timing in both slow and fast growth conditions, we found that well-studied models of cell-size control are insufficient to explain the mycobacterial cell cycle. Instead, we showed that mycobacterial cell-cycle progression is regulated by an unprecedented mechanism involving parallel adders (i.e., constant growth increments) that start at replication initiation. Together, these adders enable mycobacterial populations to regulate cell size, growth, and heterogeneity in the face of varying environments.

Efficient measurement and factorization of high-order drug interactions in Mycobacterium tuberculosis

Combinations of three or more drugs are used to treat many diseases, including tuberculosis. Thus, it is important to understand how synergistic or antagonistic drug interactions affect the efficacy of combination therapies. However, our understanding of high-order drug interactions is limited because of the lack of both efficient measurement methods and theoretical framework for analysis and interpretation. We developed an efficient experimental sampling and scoring method [diagonal measurement of n-way drug interactions (DiaMOND)] to measure drug interactions for combinations of any number of drugs. DiaMOND provides an efficient alternative to checkerboard assays, which are commonly used to measure drug interactions. We established a geometric framework to factorize high-order drug interactions into lower-order components, thereby establishing a road map of how to use lower-order measurements to predict high-order interactions. Our framework is a generalized Loewe additivity model for high-order drug interactions. Using DiaMOND, we identified and analyzed synergistic and antagonistic antibiotic combinations against Mycobacteriumtuberculosis. Efficient measurement and factorization of high-order drug interactions by DiaMOND are broadly applicable to other cell types and disease models.

Synthesis and characterization of amino acid-functionalized calcium phosphate nanoparticles for siRNA delivery

Small interfering RNAs (siRNA) are short nucleic acid fragments of about 20-27 nucleotides, which can inhibit the expression of specific genes. siRNA based RNAi technology has emerged as a promising method for the treatment of a variety of diseases. However, a major limitation in the therapeutic use of siRNA is its rapid degradation in plasma and cellular cytoplasm, resulting in short half-life. In addition, as siRNA molecules cannot penetrate into the cell efficiently, it is required to use a carrier system for its delivery. In this work, chemically and morphologically different calcium phosphate (CaP) nanoparticles, including spherical-like hydroxyapatite (HA-s), needle-like hydroxyapatite (HA-n) and calcium deficient hydroxyapatite (CDHA) nanoparticles were synthesized by the sol-gel technique and the effects of particle characteristics on the binding capacity of siRNA were investigated. In order to enhance the gene loading efficiency, the nanoparticles were functionalized with arginine and the morphological and their structural characteristics were analyzed. The addition of arginine did not significantly change the particle sizes; however, it provided a significantly increased binding of siRNA for all types of CaP nanoparticles, as revealed by spectrophotometric measurements analysis. Arginine functionalized HA-n nanoparticles showed the best binding behavior with siRNA among the other nanoparticles due to its high, positive zeta potential (+18.8mV) and high surface area of Ca⁺⁺ rich "c" plane. MTT cytotoxicity assays demonstrated that all the nanoparticles tested herein were biocompatible. Our results suggest that high siRNA entrapment in each of the three modified non-toxic CaP nanoparticles make them promising candidates as a non-viral vector for delivering therapeutic siRNA molecules to treat cancer.

Prediction of Antibiotic Interactions Using Descriptors Derived from Molecular Structure

Combination antibiotic therapies are clinically important in the fight against bacterial infections. However, the search space of drug combinations is large, making the identification of effective combinations a challenging task. Here, we present a computational framework that uses substructure profiles derived from the molecular structures of drugs and predicts antibiotic interactions. Using a previously published data set of 153 drug pairs, we showed that substructure profiles are useful in predicting synergy. We experimentally measured the interaction of 123 new drug pairs, as a prospective validation set for our approach, and identified 37 new synergistic pairs. Of the 12 pairs predicted to be synergistic, 10 were experimentally validated, corresponding to a 2.8-fold enrichment. Having thus validated our methodology, we produced a compendium of interaction predictions for all pairwise combinations among 100 antibiotics. Our methodology can make reliable antibiotic interaction predictions for any antibiotic pair within the applicability domain of the model since it solely requires chemical structures as an input.

Synergistic activity of vorinostat combined with gefitinib but not with sorafenib in mutant KRAS human non-small cell lung cancers and hepatocarcinoma

Development of drug resistance limits the efficacy of targeted therapies. Alternative approaches using different combinations of therapeutic agents to inhibit several pathways could be a more effective strategy for treating cancer. The effects of the approved epidermal growth factor receptor (EGFR)-tyrosine kinase inhibitor (gefitinib) or a multi-targeted kinase inhibitor (sorafenib) in combination with a histone deacetylase inhibitor (vorinostat) on cell proliferation, cell cycle distribution, apoptosis, and signaling pathway activation in human lung adenocarcinoma and hepatocarcinoma cells with wild-type EGFR and mutant KRAS were investigated. The effects of the synergistic drug combinations were also studied in human lung adenocarcinoma and hepatocarcinoma cells in vivo. The combination of gefitinib and vorinostat synergistically reduced cell growth and strongly induced apoptosis through inhibition of the insulin-like growth factor-1 receptor/protein kinase B (IGF-1R/AKT)-dependent signaling pathway. Moreover, the gefitinib and vorinostat combination strongly inhibited tumor growth in mice with lung adenocarcinoma or hepatocarcinoma tumor xenografts. In contrast, the combination of sorafenib and vorinostat did not inhibit cell proliferation compared to a single treatment and induced G₂/M cell cycle arrest without apoptosis. The sorafenib and vorinostat combination sustained the IGF-1R-, AKT-, and mitogen-activated protein kinase-dependent signaling pathways. These results showed that there was synergistic cytotoxicity when vorinostat was combined with gefitinib for both lung adenocarcinoma and hepatocarcinoma with mutant KRAS in vitro and in vivo but that the combination of vorinostat with sorafenib did not show any benefit. These findings highlight the important role of the IGF-1R/AKT pathway in the resistance to targeted therapies and support the use of histone deacetylase inhibitors in combination with EGFR-tyrosine kinase inhibitors, especially for treating patients with mutant KRAS resistant to other treatments.

Chemogenomics and orthology-based design of antibiotic combination therapies

Combination antibiotic therapies are being increasingly used in the clinic to enhance potency and counter drug resistance. However, the large search space of candidate drugs and dosage regimes makes the identification of effective combinations highly challenging. Here, we present a computational approach called INDIGO, which uses chemogenomics data to predict antibiotic combinations that interact synergistically or antagonistically in inhibiting bacterial growth. INDIGO quantifies the influence of individual chemical–genetic interactions on synergy and antagonism and significantly outperforms existing approaches based on experimental evaluation of novel predictions in Escherichia coli. Our analysis revealed a core set of genes and pathways (e.g. central metabolism) that are predictive of antibiotic interactions. By identifying the interactions that are associated with orthologous genes, we successfully estimated drug‐interaction outcomes in the bacterial pathogens Mycobacterium tuberculosis and Staphylococcus aureus, using the E. coli INDIGO model. INDIGO thus enables the discovery of effective combination therapies in less‐studied pathogens by leveraging chemogenomics data in model organisms.

Harnessing Connectivity in a Large-Scale Small-Molecule Sensitivity Dataset

Identifying genetic alterations that prime a cancer cell to respond to a particular therapeutic agent can facilitate the development of precision cancer medicines. Cancer cell-line (CCL) profiling of small-molecule sensitivity has emerged as an unbiased method to assess the relationships between genetic or cellular features of CCLs and small-molecule response. Here, we developed annotated cluster multidimensional enrichment analysis to explore the associations between groups of small molecules and groups of CCLs in a new, quantitative sensitivity dataset. This analysis reveals insights into small-molecule mechanisms of action, and genomic features that associate with CCL response to small-molecule treatment. We are able to recapitulate known relationships between FDA-approved therapies and cancer dependencies and to uncover new relationships, including for KRAS-mutant cancers and neuroblastoma. To enable the cancer community to explore these data, and to generate novel hypotheses, we created an updated version of the Cancer Therapeutic Response Portal (CTRP v2).

SIGNIFICANCE

We present the largest CCL sensitivity dataset yet available, and an analysis method integrating information from multiple CCLs and multiple small molecules to identify CCL response predictors robustly. We updated the CTRP to enable the cancer research community to leverage these data and analyses.

Target-independent prediction of drug synergies using only drug lipophilicity

Physicochemical properties of compounds have been instrumental in selecting lead compounds with increased drug-likeness. However, the relationship between physicochemical properties of constituent drugs and the tendency to exhibit drug interaction has not been systematically studied. We assembled physicochemical descriptors for a set of antifungal compounds (“drugs”) previously examined for interaction. Analyzing the relationship between molecular weight, lipophilicity, H-bond donor, and H-bond acceptor values for drugs and their propensity to show pairwise antifungal drug synergy, we found that combinations of two lipophilic drugs had a greater tendency to show drug synergy. We developed a more refined decision tree model that successfully predicted drug synergy in stringent cross-validation tests based on only lipophilicity of drugs. Our predictions achieved a precision of 63% and allowed successful prediction for 58% of synergistic drug pairs, suggesting that this phenomenon can extend our understanding for a substantial fraction of synergistic drug interactions. We also generated and analyzed a large-scale synergistic human toxicity network, in which we observed that combinations of lipophilic compounds show a tendency for increased toxicity. Thus, lipophilicity, a simple and easily determined molecular descriptor, is a powerful predictor of drug synergy. It is well established that lipophilic compounds (i) are promiscuous, having many targets in the cell, and (ii) often penetrate into the cell via the cellular membrane by passive diffusion. We discuss the positive relationship between drug lipophilicity and drug synergy in the context of potential drug synergy mechanisms.

Strength of selection pressure is an important parameter contributing to the complexity of antibiotic resistance evolution

Revealing the genetic changes responsible for antibiotic resistance can be critical for developing novel antibiotic therapies. However, systematic studies correlating genotype to phenotype in the context of antibiotic resistance have been missing. In order to fill in this gap, we evolved 88 isogenic Escherichia coli populations against 22 antibiotics for 3 weeks. For every drug, two populations were evolved under strong selection and two populations were evolved under mild selection. By quantifying evolved populations’ resistances against all 22 drugs, we constructed two separate cross-resistance networks for strongly and mildly selected populations. Subsequently, we sequenced representative colonies isolated from evolved populations for revealing the genetic basis for novel phenotypes. Bacterial populations that evolved resistance against antibiotics under strong selection acquired high levels of cross-resistance against several antibiotics, whereas other bacterial populations evolved under milder selection acquired relatively weaker cross-resistance. In addition, we found that strongly selected strains against aminoglycosides became more susceptible to five other drug classes compared with their wild-type ancestor as a result of a point mutation on TrkH, an ion transporter protein. Our findings suggest that selection strength is an important parameter contributing to the complexity of antibiotic resistance problem and use of high doses of antibiotics to clear infections has the potential to promote increase of cross-resistance in clinics.

A drug similarity network for understanding drug mechanism of action

Chemogenomic experiments, where genetic and chemical perturbations are combined, provide data for discovering the relationships between genotype and phenotype. Traditionally, analysis of chemogenomic datasets has been done considering the sensitivity of the deletion strains to chemicals, and this has shed light on drug mechanism of action and detecting drug targets. Here, we computationally analyzed a large chemogenomic dataset, which combines more than 300 chemicals with virtually all gene deletion strains in the yeast S. cerevisiae. In addition to sensitivity relation between deletion strains and chemicals, we also considered the deletion strains that are resistant to chemicals. We found a small set of genes whose deletion makes the cell resistant to many chemicals. Curiously, these genes were enriched for functions related to RNA metabolism. Our approach allowed us to generate a network of drugs and genes that are connected with resistance or sensitivity relationships. As a quality assessment, we showed that the higher order motifs found in this network are consistent with biological expectations. Finally, we constructed a biologically relevant network projection pertaining to drug similarities, and analyzed this network projection in detail. We propose this drug similarity network as a useful tool for understanding drug mechanism of action.

Large-scale identification and analysis of suppressive drug interactions

One drug may suppress the effects of another. Although knowledge of drug suppression is vital to avoid efficacy-reducing drug interactions or discover countermeasures for chemical toxins, drug-drug suppression relationships have not been systematically mapped. Here, we analyze the growth response of Saccharomyces cerevisiae to anti-fungal compound ("drug") pairs. Among 440 ordered drug pairs, we identified 94 suppressive drug interactions. Using only pairs not selected on the basis of their suppression behavior, we provide an estimate of the prevalence of suppressive interactions between anti-fungal compounds as 17%. Analysis of the drug suppression network suggested that Bromopyruvate is a frequently suppressive drug and Staurosporine is a frequently suppressed drug. We investigated potential explanations for suppressive drug interactions, including chemogenomic analysis, coaggregation, and pH effects, allowing us to explain the interaction tendencies of Bromopyruvate.

Oxidative stress is a mediator for increased lipid accumulation in a newly isolated Dunaliella salina strain

Green algae offer sustainable, clean and eco-friendly energy resource. However, production efficiency needs to be improved. Increasing cellular lipid levels by nitrogen depletion is one of the most studied strategies. Despite this, the underlying physiological and biochemical mechanisms of this response have not been well defined. Algae species adapted to hypersaline conditions can be cultivated in salty waters which are not useful for agriculture or consumption. Due to their inherent extreme cultivation conditions, use of hypersaline algae species is better suited for avoiding culture contamination issues. In this study, we identified a new halophilic Dunaliella salina strain by using 18S ribosomal RNA gene sequencing. We found that growth and biomass productivities of this strain were directly related to nitrogen levels, as the highest biomass concentration under 0.05 mM or 5 mM nitrogen regimes were 495 mg/l and 1409 mg/l, respectively. We also confirmed that nitrogen limitation increased cellular lipid content up to 35% under 0.05 mM nitrogen concentration. In order to gain insight into the mechanisms of this phenomenon, we applied fluorometric, flow cytometric and spectrophotometric methods to measure oxidative stress and enzymatic defence mechanisms. Under nitrogen depleted cultivation conditions, we observed increased lipid peroxidation by measuring an important oxidative stress marker, malondialdehyde and enhanced activation of catalase, ascorbate peroxidase and superoxide dismutase antioxidant enzymes. These observations indicated that oxidative stress is accompanied by increased lipid content in the green alga. In addition, we also showed that at optimum cultivation conditions, inducing oxidative stress by application of exogenous H₂O₂ leads to increased cellular lipid content up to 44% when compared with non-treated control groups. Our results support that oxidative stress and lipid overproduction are linked. Importantly, these results also suggest that oxidative stress mediates lipid accumulation. Understanding such relationships may provide guidance for efficient production of algal biodiesels.

Drugs and their interactions

Small molecules with selective efficacy can be used as drugs. Drugs typically have a therapeutic dose range where the benefit from intended effects outweighs the cost of unintended (side) effects. Herein, I use case scenarios to illustrate potential advantages and pitfalls of treatment with drugs alone or in combination. Combinations of drugs may show surprising effects given the effects of individual drugs, in phenomena known as drug interactions. Drug interactions can be classified as synergistic or antagonistic, if the effect of the combination is higher or lower than expected, respectively. A better understanding of drug interactions and their relationship to phenotypes offers hope in finding drug combinations that have high therapeutic values.

Novel anti-HER2 monoclonal antibodies: synergy and antagonism with tumor necrosis factor-α

Background

One-third of breast cancers display amplifications of the ERBB2 gene encoding the HER2 kinase receptor. Trastuzumab, a humanized antibody directed against an epitope on subdomain IV of the extracellular domain of HER2 is used for therapy of HER2-overexpressing mammary tumors. However, many tumors are either natively resistant or acquire resistance against Trastuzumab. Antibodies directed to different epitopes on the extracellular domain of HER2 are promising candidates for replacement or combinatorial therapy. For example, Pertuzumab that binds to subdomain II of HER2 extracellular domain and inhibits receptor dimerization is under clinical trial. Alternative antibodies directed to novel HER2 epitopes may serve as additional tools for breast cancer therapy. Our aim was to generate novel anti-HER2 monoclonal antibodies inhibiting the growth of breast cancer cells, either alone or in combination with tumor necrosis factor-α (TNF-α).

Methods

Mice were immunized against SK-BR-3 cells and recombinant HER2 extracellular domain protein to produce monoclonal antibodies. Anti-HER2 antibodies were characterized with breast cancer cell lines using immunofluorescence, flow cytometry, immunoprecipitation, western blot techniques. Antibody epitopes were localized using plasmids encoding recombinant HER2 protein variants. Antibodies, either alone or in combination with TNF-α, were tested for their effects on breast cancer cell proliferation.

Results

We produced five new anti-HER2 monoclonal antibodies, all directed against conformational epitope or epitopes restricted to the native form of the extracellular domain. When tested alone, some antibodies inhibited modestly but significantly the growth of SK-BR-3, BT-474 and MDA-MB-361 cells displaying ERBB2 amplification. They had no detectable effect on MCF-7 and T47D cells lacking ERBB2 amplification. When tested in combination with TNF-α, antibodies acted synergistically on SK-BR-3 cells, but antagonistically on BT-474 cells. A representative anti-HER2 antibody inhibited Akt and ERK1/2 phosphorylation leading to cyclin D1 accumulation and growth arrest in SK-BR-3 cells, independently from TNF-α.

Conclusions

Novel antibodies against extracellular domain of HER2 may serve as potent anti-cancer bioactive molecules. Cell-dependent synergy and antagonism between anti-HER2 antibodies and TNF-α provide evidence for a complex interplay between HER2 and TNF-α signaling pathways. Such complexity may drastically affect the outcome of HER2-directed therapeutic interventions.

Systematic exploration of synergistic drug pairs

Drug synergy allows a therapeutic effect to be achieved with lower doses of component drugs. Drug synergy can result when drugs target the products of genes that act in parallel pathways ('specific synergy'). Such cases of drug synergy should tend to correspond to synergistic genetic interaction between the corresponding target genes. Alternatively, 'promiscuous synergy' can arise when one drug non-specifically increases the effects of many other drugs, for example, by increased bioavailability. To assess the relative abundance of these drug synergy types, we examined 200 pairs of antifungal drugs in S. cerevisiae. We found 38 antifungal synergies, 37 of which were novel. While 14 cases of drug synergy corresponded to genetic interaction, 92% of the synergies we discovered involved only six frequently synergistic drugs. Although promiscuity of four drugs can be explained under the bioavailability model, the promiscuity of Tacrolimus and Pentamidine was completely unexpected. While many drug synergies correspond to genetic interactions, the majority of drug synergies appear to result from non-specific promiscuous synergy.

The genetic landscape of a cell

A genome-scale genetic interaction map was constructed by examining 5.4 million gene-gene pairs for synthetic genetic interactions, generating quantitative genetic interaction profiles for approximately 75% of all genes in the budding yeast, Saccharomyces cerevisiae. A network based on genetic interaction profiles reveals a functional map of the cell in which genes of similar biological processes cluster together in coherent subsets, and highly correlated profiles delineate specific pathways to define gene function. The global network identifies functional cross-connections between all bioprocesses, mapping a cellular wiring diagram of pleiotropy. Genetic interaction degree correlated with a number of different gene attributes, which may be informative about genetic network hubs in other organisms. We also demonstrate that extensive and unbiased mapping of the genetic landscape provides a key for interpretation of chemical-genetic interactions and drug target identification.

A protein domain-based interactome network for C. elegans early embryogenesis

Many protein-protein interactions are mediated through independently folding modular domains. Proteome-wide efforts to model protein-protein interaction or "interactome" networks have largely ignored this modular organization of proteins. We developed an experimental strategy to efficiently identify interaction domains and generated a domain-based interactome network for proteins involved in C. elegans early-embryonic cell divisions. Minimal interacting regions were identified for over 200 proteins, providing important information on their domain organization. Furthermore, our approach increased the sensitivity of the two-hybrid system, resulting in a more complete interactome network. This interactome modeling strategy revealed insights into C. elegans centrosome function and is applicable to other biological processes in this and other organisms.

Visualizing evolution and impact of biomedical fields

We describe a new tool for visualization of biomedical scientific trends. The method captures variations in scientific impact over time to allow for a comparison of relative significance and evolution of fields similar to a financial market scorecard. The tool is available at SciTrends (http://www.scitrends.net), depicting the evolution of almost 200 thousand biomedical fields in time. With millions of articles on thousands of topics published in biomedicine, we envision that only with such large-scale tools researchers can objectively understand the ever-changing interests in the biomedical sciences and make more informed decisions.

A recipe for high impact

Our analysis highlights common statistical features of high-impact articles; we also show how information flows among various publication types.

Finding nuclear localization signals

A variety of nuclear localization signals (NLSs) are experimentally known although only one motif was available for database searches through PROSITE. We initially collected a set of 91 experimentally verified NLSs from the literature. Through iterated 'in silico mutagenesis' we then extended the set to 214 potential NLSs. This final set matched in 43% of all known nuclear proteins and in no known non-nuclear protein. We estimated that >17% of all eukaryotic proteins may be imported into the nucleus. Finally, we found an overlap between the NLS and DNA-binding region for 90% of the proteins for which both the NLS and DNA-binding regions were known. Thus, evolution seemed to have used part of the existing DNA-binding mechanism when compartmentalizing DNA-binding proteins into the nucleus. However, only 56 of our 214 NLS motifs overlapped with DNA-binding regions. These 56 NLSs enabled a de novo prediction of partial DNA-binding regions for approximately 800 proteins in human, fly, worm and yeast.