The idea of synergistic interactions between drugs and chemicals has been an important issue in the biomedical world for over a century. As complex diseases, especially cancer, are being treated with ...various drug cocktails, understanding the interactions among these drugs is increasingly vital to ensuring successful treatment regimens. However, the idea of synergy is not limited to only the biomedical realm and these ideas have developed across many different disciplines, as well. In this review, we first discuss the various terminology surrounding the idea of synergy, providing a comprehensive list of terms defined across numerous disciplines. We then review the most common methodology for detection and quantification of synergy, including the two most prominent reference models for describing additive interactions: Loewe Additivity and Bliss Independence. We also discuss advantages and limitations to each method, with a focus on the Chou-Talalay Combination Index method. Finally, we describe how methods development and terminology have developed among disciplines outside of biomedicine and pharmacology, to synthesize the literature for readers.
A GROMOS force-field parameter set 54A8 is developed, which is based on the latest 54A7 set Schmid et al. Eur. Biophys. J. 2011, 40, 843–856 and involves a recalibration of the nonbonded interaction ...parameters for the charged amino acid side chains, based on ionic side chain analogs. After a thorough analysis of the available experimental data, conventional hydration free energies for the ammonium; mono-, di-, tri-, and tetramethyl-ammonium; formate; acetate; propanoate; imidazolium; and guanidinium ions are combined with a standard absolute intrinsic proton hydration free energy ΔG hyd ⊖Hg + = −1100 kJ·mol–1 to yield absolute intrinsic single-ion hydration free energies serving as experimental target data. The raw hydration free energies calculated from atomistic simulations are affected by electrostatic and finite-size artifacts, and corrections are applied to reach methodological independence prior to comparison with these experimental values. Except for monomethyl-ammonium, ions with parameters derived directly from the 54A7 force field considerably underestimate (ammonium, formate, acetate, propanoate, guanidinium) or overestimate (di-, tri-, and tetramethyl-ammonium; imidazolium) the magnitude of the intrinsic hydration free energy, the largest deviation affecting the acetate ion (40.0 kJ·mol–1). After reparameterization into 54A8, the mean and maximal absolute deviations between simulated and experimental data over the set of 10 ions are reduced from 23.1 and 40.0 kJ·mol–1, respectively, to 1.8 and 6.3 kJ·mol–1, respectively. Although the 54A7 and 54A8 parameter sets differ significantly in terms of the hydration free energies of the ions considered, other properties such as ion–water radial distribution functions and ion–ion potentials of mean force appear to be only moderately sensitive to this change. These properties are similar for the two sets and, in the case of the ion–water radial distribution functions, in good agreement with available experimental data.
Cancer patients exhibit a broad range of inter-individual variability in response and toxicity to widely used anticancer drugs, and genetic variation is a major contributor to this variability. To ...identify new genes that influence the response of 44 FDA-approved anticancer drug treatments widely used to treat various types of cancer, we conducted high-throughput screening and genome-wide association mapping using 680 lymphoblastoid cell lines from the 1000 Genomes Project. The drug treatments considered in this study represent nine drug classes widely used in the treatment of cancer in addition to the paclitaxel + epirubicin combination therapy commonly used for breast cancer patients. Our genome-wide association study (GWAS) found several significant and suggestive associations. We prioritized consistent associations for functional follow-up using gene-expression analyses. The NAD(P)H quinone dehydrogenase 1 (NQO1) gene was found to be associated with the dose-response of arsenic trioxide, erlotinib, trametinib, and a combination treatment of paclitaxel + epirubicin. NQO1 has previously been shown as a biomarker of epirubicin response, but our results reveal novel associations with these additional treatments. Baseline gene expression of NQO1 was positively correlated with response for 43 of the 44 treatments surveyed. By interrogating the functional mechanisms of this association, the results demonstrate differences in both baseline and drug-exposed induction.
Drawing integrated conclusions from diverse source data requires synthesis across multiple types of information. The ToxPi (Toxicological Prioritization Index) is an analytical framework that was ...developed to enable integration of multiple sources of evidence by transforming data into integrated, visual profiles. Methodological improvements have advanced ToxPi and expanded its applicability, necessitating a new, consolidated software platform to provide functionality, while preserving flexibility for future updates.
We detail the implementation of a new graphical user interface for ToxPi (Toxicological Prioritization Index) that provides interactive visualization, analysis, reporting, and portability. The interface is deployed as a stand-alone, platform-independent Java application, with a modular design to accommodate inclusion of future analytics. The new ToxPi interface introduces several features, from flexible data import formats (including legacy formats that permit backward compatibility) to similarity-based clustering to options for high-resolution graphical output.
We present the new ToxPi interface for dynamic exploration, visualization, and sharing of integrated data models. The ToxPi interface is freely-available as a single compressed download that includes the main Java executable, all libraries, example data files, and a complete user manual from http://toxpi.org .
There are currently 85,000 chemicals registered with the Environmental Protection Agency (EPA) under the Toxic Substances Control Act, but only a small fraction have measured toxicological data. To ...address this gap, high-throughput screening (HTS) and computational methods are vital. As part of one such HTS effort, embryonic zebrafish were used to examine a suite of morphological and mortality endpoints at six concentrations from over 1,000 unique chemicals found in the ToxCast library (phase 1 and 2). We hypothesized that by using a conditional generative adversarial network (cGAN) or deep neural networks (DNN), and leveraging this large set of toxicity data we could efficiently predict toxic outcomes of untested chemicals. Utilizing a novel method in this space, we converted the 3D structural information into a weighted set of points while retaining all information about the structure. In vivo toxicity and chemical data were used to train two neural network generators. The first was a DNN (Go-ZT) while the second utilized cGAN architecture (GAN-ZT) to train generators to produce toxicity data. Our results showed that Go-ZT significantly outperformed the cGAN, support vector machine, random forest and multilayer perceptron models in cross-validation, and when tested against an external test dataset. By combining both Go-ZT and GAN-ZT, our consensus model improved the SE, SP, PPV, and Kappa, to 71.4%, 95.9%, 71.4% and 0.673, respectively, resulting in an area under the receiver operating characteristic (AUROC) of 0.837. Considering their potential use as prescreening tools, these models could provide in vivo toxicity predictions and insight into the hundreds of thousands of untested chemicals to prioritize compounds for HT testing.
Background: Thirty years of pesticide registration toxicity data have been historically stored as hardcopy and scanned documents by the U.S. Environmental Protection Agency (EPA). A significant ...portion of these data have now been processed into standardized and structured toxicity data within the EPA's Toxicity Reference Database (ToxRefDB), including chronic, cancer, developmental, and reproductive studies from laboratory animals. These data are now accessible and mineable within ToxRefDB and are serving as a primary source of validation for U.S. EPA's ToxCast research program in predictive toxicology. Objectives: We profiled in vivo toxicities across 310 chemicals as a model application of ToxRefDB, meeting the need for detailed anchoring end points for development of ToxCast predictive signatures. Methods: Using query and structured data-mining approaches, we generated toxicity profiles from ToxRefDB based on long-term rodent bioassays. These chronic/cancer data were analyzed for suitability as anchoring end points based on incidence, target organ, severity, potency, and significance. Results: Under conditions of the bioassays, we observed pathologies for 273 of 310 chemicals, with greater preponderance (> 90%) occurring in the liver, kidney, thyroid, lung, testis, and spleen. We observed proliferative lesions for 225 chemicals, and 167 chemicals caused progression to cancer-related pathologies. Conclusions: Based on incidence, severity, and potency, we selected 26 primarily tissue-specific pathology end points to uniformly classify the 310 chemicals. The resulting toxicity profile classifications demonstrate the utility of structuring legacy toxicity information and facilitating the computation of these data within ToxRefDB for ToxCast and other applications.
The potential for cardiotoxicity is carefully evaluated for pharmaceuticals, as it is a major safety liability. However, environmental chemicals are seldom tested for their cardiotoxic potential. ...Moreover, there is a large variability in both baseline and drug-induced cardiovascular risk in humans, but data are lacking on the degree to which susceptibility to chemically-induced cardiotoxicity may also vary. Human induced pluripotent stem cell (iPSC)-derived cardiomyocytes have become an important in vitro model for drug screening. Thus, we hypothesized that a population-based model of iPSC-derived cardiomyocytes from a diverse set of individuals can be used to assess potential hazard and inter-individual variability in chemical effects on these cells. We conducted concentration-response screening of 134 chemicals (pharmaceuticals, industrial and environmental chemicals and food constituents) in iPSC-derived cardiomyocytes from 43 individuals, comprising both sexes and diverse ancestry. We measured kinetic calcium flux and conducted high-content imaging following chemical exposure, and utilized a panel of functional and cytotoxicity parameters in concentration-response for each chemical and donor. We show reproducible inter-individual variability in both baseline and chemical-induced effects on iPSC-derived cardiomyocytes. Further, chemical-specific variability in potency and degree of population variability were quantified. This study shows the feasibility of using an organotypic population-based human in vitro model to quantitatively assess chemicals for which little cardiotoxicity information is available. Ultimately, these results advance in vitro toxicity testing methodologies by providing an innovative tool for population-based cardiotoxicity screening, contributing to the paradigm shift from traditional animal models of toxicity to in vitro toxicity testing methods.
•Cardiotoxicity information is lacking for most environmental chemicals.•Cardiovascular disease risks vary greatly across the population.•We demonstrated that hiPSC-derived cardiomyocytes can address these gaps.•Potency and population variability in cardiotoxicity varied across chemicals.•This population-based model can substantially improve chemical safety evaluations.
The asymmetric solvation of ions can be defined as the tendency of a solvent to preferentially solvate anions over cations or cations over anions, at identical ionic charge magnitudes and effective ...sizes. Taking water as a reference, these effects are quantified experimentally for many solvents by the relative acity (A) and basity (B) parameters of the Swain scale. The goal of the present study is to investigate the asymmetric solvation of ions using molecular dynamics simulations, and to connect the results to this empirical scale. To this purpose, the charging free energies of alkali and halide ions, and of their hypothetical oppositely charged counterparts, are calculated in a variety of solvents. In a first set of calculations, artificial solvent models are considered that present either a charge or a shape asymmetry at the molecular level. The solvation asymmetry, probed by the difference in charging free energy between the two oppositely charged ions, is found to encompass a term quadratic in the ion charge, related to the different solvation structures around the anion and cation, and a term linear in the ion charge, related to the solvation structure around the uncharged ion-sized cavity. For these simple solvent models, the two terms are systematically counteracting each other, and it is argued that only the quadratic term should be retained when comparing the results of simulations involving physical solvents to experimental data. In a second set of calculations, 16 physical solvents are considered. The theoretical estimates for the acity A are found to correlate very well with the Swain parameters, whereas the correlation for B is very poor. Based on this observation, the Swain scale is reformulated into a new scale involving an asymmetry parameter Σ, positive for acitic solvents and negative for basitic ones, and a polarity parameter Π. This revised scale has the same predictive power as the original scale, but it characterizes asymmetry in an absolute sense, the atomistic simulations playing the role of an extra-thermodynamic assumption, and is optimally compatible with the simulation results. Considering the 55 solvents in the Swain set, it is observed that a moderate basity (Σ between −0.9 and −0.3, related to electronic polarization) represents the baseline for most solvents, while a highly variable acity (Σ between 0.0 and 3.0, related to hydrogen-bond donor capacity modulated by inductive effects) represents a landmark of protic solvents.