The rise of multi-drug resistant (MDR) and extensively drug resistant (XDR) tuberculosis around the world, including in industrialized nations, poses a great threat to human health and defines a need to develop new, effective and inexpensive anti-tubercular agents. Previously we developed a chemical systems biology approach to identify off-targets of major pharmaceuticals on a proteome-wide scale. In this paper we further demonstrate the value of this approach through the discovery that existing commercially available drugs, prescribed for the treatment of Parkinson's disease, have the potential to treat MDR and XDR tuberculosis. These drugs, entacapone and tolcapone, are predicted to bind to the enzyme InhA and directly inhibit substrate binding. The prediction is validated by in vitro and InhA kinetic assays using tablets of Comtan, whose active component is entacapone. The minimal inhibition concentration (MIC99) of entacapone for Mycobacterium tuberculosis (M.tuberculosis) is approximately 260.0 µM, well below the toxicity concentration determined by an in vitro cytotoxicity model using a human neuroblastoma cell line. Moreover, kinetic assays indicate that Comtan inhibits InhA activity by 47.0% at an entacapone concentration of approximately 80 µM. Thus the active component in Comtan represents a promising lead compound for developing a new class of anti-tubercular therapeutics with excellent safety profiles. More generally, the protocol described in this paper can be included in a drug discovery pipeline in an effort to discover novel drug leads with desired safety profiles, and therefore accelerate the development of new drugs.

Polypharmacology, which focuses on designing therapeutics to target multiple receptors, has emerged as a new paradigm in drug discovery. Polypharmacological effects are an attribute of most, if not all, drug molecules. The efficacy and toxicity of drugs, whether designed as single- or multitarget therapeutics, result from complex interactions between pharmacodynamic, pharmacokinetic, genetic, epigenetic, and environmental factors. Ultimately, to predict a drug response phenotype, it is necessary to understand the change in information flow through cellular networks resulting from dynamic drug-target interactions and the impact that this has on the complete biological system. Although such is a future objective, we review recent progress and challenges in computational techniques that enable the prediction and analysis of in vitro and in vivo drug-response phenotypes.

Whole-genome sequencing (WGS) of maternal plasma cell-free DNA (cfDNA) can potentially evaluate all 24 chromosomes to identify abnormalities of the placenta, fetus, or pregnant woman. Current bioinformatics algorithms typically only report on chromosomes 21, 18, 13, X, and Y; sequencing results from other chromosomes may be masked. We hypothesized that by systematically analyzing WGS data from all chromosomes, we could identify rare autosomal trisomies (RATs) to improve understanding of feto-placental biology. We analyzed two independent cohorts from clinical laboratories, both of which used a similar quality control parameter, normalized chromosome denominator quality. The entire data set included 89,817 samples. Samples flagged for analysis and classified as abnormal were 328 of 72,932 (0.45%) and 71 of 16,885 (0.42%) in cohorts 1 and 2, respectively. Clinical outcome data were available for 57 of 71 (80%) of abnormal cases in cohort 2. Visual analysis of WGS data demonstrated RATs, copy number variants, and extensive genome-wide imbalances. Trisomies 7, 15, 16, and 22 were the most frequently observed RATs in both cohorts. Cytogenetic or pregnancy outcome data were available in 52 of 60 (87%) of cases with RATs in cohort 2. Cases with RATs detected were associated with miscarriage, true fetal mosaicism, and confirmed or suspected uniparental disomy. Comparing the trisomic fraction with the fetal fraction allowed estimation of possible mosaicism. Analysis and reporting of aneuploidies in all chromosomes can clarify cases in which cfDNA findings on selected "target" chromosomes (21, 18, and 13) are discordant with the fetal karyotype and may identify pregnancies at risk of miscarriage and other complications.

Docking scoring functions are notoriously weak predictors of binding affinity. They typically assign a common set of weights to the individual energy terms that contribute to the overall energy score, however, these weights should be gene family-dependent. In addition, they incorrectly assume that individual interactions contribute towards the total binding affinity in an additive manner. In reality, noncovalent interactions often depend on one another in a nonlinear manner. In this paper we show how the use of support vector machines (SVMs), trained by associating sets of individual energy terms retrieved from molecular docking with the known binding affinity of each compound from high-throughput screening experiments, can be used to improve the correlation between known binding affinities and those predicted by the docking program eHiTS. We construct two prediction models; a regression model trained using IC 50 values from BindingDB, and a classification model trained using active and decoy compounds from the Directory of Useful Decoys (DUD). Moreover, to address the issue of overrepresentation of negative data in highthroughput screening data sets, we have designed a multiple-planar SVM training procedure for the classification model. The increased performance that both SVMs give when compared with the original eHiTS scoring function highlights the potential for using nonlinear methods when deriving overall energy scores from their individual components. We apply the above methodology to train a new scoring function for direct inhibitors of M.tuberculosis (M.tb) InhA. By combining ligand binding site comparison with the new scoring function, we propose that phosphodiesterase inhibitors can potentially be repurposed to target M.tb InhA. Our methodology may be applied to other gene families for which target structures and activity data are available, as demonstrated in the work presented here.

Methods for analyzing complete gene families are becoming of increasing importance to the drug discovery process, because similarities and differences within a family are often the key to understanding functional differences that can be exploited in drug design. We undertake a large-scale structural comparison of protein kinase ATP-binding sites using a geometric hashing method. Subsequently, we propose a relevant classification of the protein kinase family based on the structural similarity of its binding sites. Our classification is not only able to reveal the great diversity of different protein kinases and therefore their different potential for inhibitor selectivity but it is also able to distinguish subtle differences within binding site conformation reflecting the protein activation state. Furthermore, using experimental inhibition profiling, we demonstrate that our classification can be used to identify protein kinase binding sites that are known experimentally to bind the same drug, demonstrating that it has potential as an inverse (protein) virtual screening tool, by identifying which other sites have the potential to bind a given drug. In this way the cross-reactivities of the anticancer drugs Tarceva and Gleevec are rationalized.

While a number of factors trend with fetal fraction across the cohort as a whole, they are not the sole determinants of fetal fraction. In this study, the variability for any one patient does not appear large enough to justify postponing testing to a later gestational age.

How easy is it to reproduce the results found in a typical computational biology paper? Either through experience or intuition the reader will already know that the answer is with difficulty or not at all. In this paper we attempt to quantify this difficulty by reproducing a previously published paper for different classes of users (ranging from users with little expertise to domain experts) and suggest ways in which the situation might be improved. Quantification is achieved by estimating the time required to reproduce each of the steps in the method described in the original paper and make them part of an explicit workflow that reproduces the original results. Reproducing the method took several months of effort, and required using new versions and new software that posed challenges to reconstructing and validating the results. The quantification leads to “reproducibility maps” that reveal that novice researchers would only be able to reproduce a few of the steps in the method, and that only expert researchers with advance knowledge of the domain would be able to reproduce the method in its entirety. The workflow itself is published as an online resource together with supporting software and data. The paper concludes with a brief discussion of the complexities of requiring reproducibility in terms of cost versus benefit, and a desiderata with our observations and guidelines for improving reproducibility. This has implications not only in reproducing the work of others from published papers, but reproducing work from one’s own laboratory.

We report a computational approach that integrates structural bioinformatics, molecular modelling and systems biology to construct a drug-target network on a structural proteome-wide scale. The approach has been applied to the genome of Mycobacterium tuberculosis (M.tb), the causative agent of one of today's most widely spread infectious diseases. The resulting drug-target interaction network for all structurally characterized approved drugs bound to putative M.tb receptors, we refer to as the ‘TB-drugome’. The TB-drugome reveals that approximately one-third of the drugs examined have the potential to be repositioned to treat tuberculosis and that many currently unexploited M.tb receptors may be chemically druggable and could serve as novel anti-tubercular targets. Furthermore, a detailed analysis of the TB-drugome has shed new light on the controversial issues surrounding drug-target networks [1]–[3]. Indeed, our results support the idea that drug-target networks are inherently modular, and further that any observed randomness is mainly caused by biased target coverage. The TB-drugome (http://funsite.sdsc.edu/drugome/TB) has the potential to be a valuable resource in the development of safe and efficient anti-tubercular drugs. More generally the methodology may be applied to other pathogens of interest with results improving as more of their structural proteomes are determined through the continued efforts of structural biology/genomics.