A deep neural network: mechanistic hybrid model to predict pharmacokinetics in rat.

An important aspect in the development of small molecules as drugs or agrochemicals is their systemic availability after intravenous and oral administration. The prediction of the systemic availability from the chemical structure of a potential candidate is highly desirable, as it allows to focus the drug or agrochemical development on compounds with a favorable kinetic profile. However, such predictions are challenging as the availability is the result of the complex interplay between molecular properties, biology and physiology and training data is rare. In this work we improve the hybrid model developed earlier (Schneckener in J Chem Inf Model 59:4893-4905, 2019). We reduce the median fold change error for the total oral exposure from 2.85 to 2.35 and for intravenous administration from 1.95 to 1.62. This is achieved by training on a larger data set, improving the neural network architecture as well as the parametrization of mechanistic model. Further, we extend our approach to predict additional endpoints and to handle different covariates, like sex and dosage form. In contrast to a pure machine learning model, our model is able to predict new end points on which it has not been trained. We demonstrate this feature by predicting the exposure over the first 24 h, while the model has only been trained on the total exposure.

Identification of the Highly Active, Species Cross-Reactive Complex I Inhibitor BAY-179.

Mitochondria are key regulators of energy supply and cell death. Generation of ATP within mitochondria occurs through oxidative phosphorylation (OXPHOS), a process which utilizes the four complexes (complex I-IV) of the electron transport chain and ATP synthase. Certain oncogenic mutations (e.g., LKB1 or mIDH) can further enhance the reliance of cancer cells on OXPHOS for their energetic requirements, rendering cells sensitive to complex I inhibition and highlighting the potential value of complex I as a therapeutic target. Herein, we describe the discovery of a potent, selective, and species cross-reactive complex I inhibitor. A high-throughput screen of the Bayer compound library followed by hit triaging and initial hit-to-lead activities led to a lead structure which was further optimized in a comprehensive lead optimization campaign. Focusing on balancing potency and metabolic stability, this program resulted in the identification of BAY-179, an excellent in vivo suitable tool with which to probe the biological relevance of complex I inhibition in cancer indications.

Prediction of Oral Bioavailability in Rats: Transferring Insights from in Vitro Correlations to (Deep) Machine Learning Models Using in Silico Model Outputs and Chemical Structure Parameters.

Oral administration of drug products is a strict requirement in many medical indications. Therefore, bioavailability prediction models are of high importance for prioritization of compound candidates in the drug discovery process. However, oral exposure and bioavailability are difficult to predict, as they are the result of various highly complex factors and/or processes influenced by the physicochemical properties of a compound, such as solubility, lipophilicity, or charge state, as well as by interactions with the organism, for instance, metabolism or membrane permeation. In this study, we assess whether it is possible to predict intravenous (iv) or oral drug exposure and oral bioavailability in rats. As input parameters, we use (i) six experimentally determined in vitro and physicochemical endpoints, namely, membrane permeation, free fraction, metabolic stability, solubility, pKa value, and lipophilicity; (ii) the outputs of six in silico absorption, distribution, metabolism, and excretion models trained on the same endpoints, or (iii) the chemical structure encoded as fingerprints or simplified molecular input line entry system strings. The underlying data set for the models is an unprecedented collection of almost 1900 data points with high-quality in vivo experiments performed in rats. We find that drug exposure after iv administration can be predicted similarly well using hybrid models with in vitro- or in silico-predicted endpoints as inputs, with fold change errors (FCE) of 2.28 and 2.08, respectively. The FCEs for exposure after oral administration are higher, and here, the prediction from in vitro inputs performs significantly better in comparison to in silico-based models with FCEs of 3.49 and 2.40, respectively, most probably reflecting the higher complexity of oral bioavailability. Simplifying the prediction task to a binary alert for low oral bioavailability, based only on chemical structure, we achieve accuracy and precision close to 70%.

Binding Kinetics Survey of the Drugged Kinome.

Target residence time is emerging as an important optimization parameter in drug discovery, yet target and off-target engagement dynamics have not been clearly linked to the clinical performance of drugs. Here we developed high-throughput binding kinetics assays to characterize the interactions of 270 protein kinase inhibitors with 40 clinically relevant targets. Analysis of the results revealed that on-rates are better correlated with affinity than off-rates and that the fraction of slowly dissociating drug-target complexes increases from early/preclinical to late stage and FDA-approved compounds, suggesting distinct contributions by each parameter to clinical success. Combining binding parameters with PK/ADME properties, we illustrate in silico and in cells how kinetic selectivity could be exploited as an optimization strategy. Furthermore, using bio- and chemoinformatics we uncovered structural features influencing rate constants. Our results underscore the value of binding kinetics information in rational drug design and provide a resource for future studies on this subject.

Quantifying the degree of persistence in random amoeboid motion based on the Hurst exponent of fractional Brownian motion.

Amoebae explore their environment in a random way, unless external cues like, e.g., nutrients, bias their motion. Even in the absence of cues, however, experimental cell tracks show some degree of persistence. In this paper, we analyzed individual cell tracks in the framework of a linear mixed effects model, where each track is modeled by a fractional Brownian motion, i.e., a Gaussian process exhibiting a long-term correlation structure superposed on a linear trend. The degree of persistence was quantified by the Hurst exponent of fractional Brownian motion. Our analysis of experimental cell tracks of the amoeba Dictyostelium discoideum showed a persistent movement for the majority of tracks. Employing a sliding window approach, we estimated the variations of the Hurst exponent over time, which allowed us to identify points in time, where the correlation structure was distorted ("outliers"). Coarse graining of track data via down-sampling allowed us to identify the dependence of persistence on the spatial scale. While one would expect the (mode of the) Hurst exponent to be constant on different temporal scales due to the self-similarity property of fractional Brownian motion, we observed a trend towards stronger persistence for the down-sampled cell tracks indicating stronger persistence on larger time scales.

Estimating HIV-1 fitness characteristics from cross-sectional genotype data.

Despite the success of highly active antiretroviral therapy (HAART) in the management of human immunodeficiency virus (HIV)-1 infection, virological failure due to drug resistance development remains a major challenge. Resistant mutants display reduced drug susceptibilities, but in the absence of drug, they generally have a lower fitness than the wild type, owing to a mutation-incurred cost. The interaction between these fitness costs and drug resistance dictates the appearance of mutants and influences viral suppression and therapeutic success. Assessing in vivo viral fitness is a challenging task and yet one that has significant clinical relevance. Here, we present a new computational modelling approach for estimating viral fitness that relies on common sparse cross-sectional clinical data by combining statistical approaches to learn drug-specific mutational pathways and resistance factors with viral dynamics models to represent the host-virus interaction and actions of drug mechanistically. We estimate in vivo fitness characteristics of mutant genotypes for two antiretroviral drugs, the reverse transcriptase inhibitor zidovudine (ZDV) and the protease inhibitor indinavir (IDV). Well-known features of HIV-1 fitness landscapes are recovered, both in the absence and presence of drugs. We quantify the complex interplay between fitness costs and resistance by computing selective advantages for different mutants. Our approach extends naturally to multiple drugs and we illustrate this by simulating a dual therapy with ZDV and IDV to assess therapy failure. The combined statistical and dynamical modelling approach may help in dissecting the effects of fitness costs and resistance with the ultimate aim of assisting the choice of salvage therapies after treatment failure.

HIV quasispecies dynamics during pro-active treatment switching: impact on multi-drug resistance and resistance archiving in latent reservoirs.

The human immunodeficiency virus (HIV) can be suppressed by highly active anti-retroviral therapy (HAART) in the majority of infected patients. Nevertheless, treatment interruptions inevitably result in viral rebounds from persistent, latently infected cells, necessitating lifelong treatment. Virological failure due to resistance development is a frequent event and the major threat to treatment success. Currently, it is recommended to change treatment after the confirmation of virological failure. However, at the moment virological failure is detected, drug resistant mutants already replicate in great numbers. They infect numerous cells, many of which will turn into latently infected cells. This pool of cells represents an archive of resistance, which has the potential of limiting future treatment options. The objective of this study was to design a treatment strategy for treatment-naive patients that decreases the likelihood of early treatment failure and preserves future treatment options. We propose to apply a single, pro-active treatment switch, following a period of treatment with an induction regimen. The main goal of the induction regimen is to decrease the abundance of randomly generated mutants that confer resistance to the maintenance regimen, thereby increasing subsequent treatment success. Treatment is switched before the overgrowth and archiving of mutant strains that carry resistance against the induction regimen and would limit its future re-use. In silico modelling shows that an optimal trade-off is achieved by switching treatment at days after the initiation of antiviral therapy. Evaluation of the proposed treatment strategy demonstrated significant improvements in terms of resistance archiving and virological response, as compared to conventional HAART. While continuous pro-active treatment alternation improved the clinical outcome in a randomized trial, our results indicate that a similar improvement might also be reached after a single pro-active treatment switch. The clinical validity of this finding, however, remains to be shown by a corresponding trial.

Drug-class specific impact of antivirals on the reproductive capacity of HIV.

Predictive markers linking drug efficacy to clinical outcome are a key component in the drug discovery and development process. In HIV infection, two different measures, viral load decay and phenotypic assays, are used to assess drug efficacy in vivo and in vitro. For the newly introduced class of integrase inhibitors, a huge discrepancy between these two measures of efficacy was observed. Hence, a thorough understanding of the relation between these two measures of drug efficacy is imperative for guiding future drug discovery and development activities in HIV. In this article, we developed a novel viral dynamics model, which allows for a mechanistic integration of the mode of action of all approved drugs and drugs in late clinical trials. Subsequently, we established a link between in vivo and in vitro measures of drug efficacy, and extract important determinants of drug efficacy in vivo. The analysis is based on a new quantity-the reproductive capacity-that represents in mathematical terms the in vivo analog of the read-out of a phenotypic assay. Our results suggest a drug-class specific impact of antivirals on the total amount of viral replication. Moreover, we showed that the (drug-)target half life, dominated by immune-system related clearance processes, is a key characteristic that affects both the emergence of resistance as well as the in vitro-in vivo correlation of efficacy measures in HIV treatment. We found that protease- and maturation inhibitors, due to their target half-life, decrease the total amount of viral replication and the emergence of resistance most efficiently.