Next Article in Journal
Clinical Amyloid Typing by Proteomics: Performance Evaluation and Data Sharing between Two Centres
Previous Article in Journal
Highly Oxygenated Triterpenoids and Diterpenoids from Fructus Rubi (Rubus chingii Hu) and Their NF-kappa B Inhibitory Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Accelerated Repurposing and Drug Development of Pulmonary Hypertension Therapies for COVID-19 Treatment Using an AI-Integrated Biosimulation Platform

VeriSIM Life, 1 Sansome Street, Suite 3500, San Francisco, CA 94104, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2021, 26(7), 1912; https://doi.org/10.3390/molecules26071912
Submission received: 7 March 2021 / Revised: 25 March 2021 / Accepted: 26 March 2021 / Published: 29 March 2021
(This article belongs to the Section Computational and Theoretical Chemistry)

Abstract

:
The COVID-19 pandemic has reached over 100 million worldwide. Due to the multi-targeted nature of the virus, it is clear that drugs providing anti-COVID-19 effects need to be developed at an accelerated rate, and a combinatorial approach may stand to be more successful than a single drug therapy. Among several targets and pathways that are under investigation, the renin-angiotensin system (RAS) and specifically angiotensin-converting enzyme (ACE), and Ca2+-mediated SARS-CoV-2 cellular entry and replication are noteworthy. A combination of ACE inhibitors and calcium channel blockers (CCBs), a critical line of therapy for pulmonary hypertension, has shown therapeutic relevance in COVID-19 when investigated independently. To that end, we conducted in silico modeling using BIOiSIM, an AI-integrated mechanistic modeling platform by utilizing known preclinical in vitro and in vivo datasets to accurately simulate systemic therapy disposition and site-of-action penetration of the CCBs and ACEi compounds to tissues implicated in COVID-19 pathogenesis.

1. Introduction

The emerging pandemic of COVID-19 presents an extraordinary challenge in identifying effective drugs for prevention and cure [1,2]. Analysis of the cumulative surveillance data have shown progressively growing number of cases, which are now confirmed to be 100 million plus, as demonstrated by the World Health Organization emergency COVID19 informational dashboard [3]. In the current pandemic scenario, the average mortality rate was observed to be almost 3.0% [4]. Although the average mortality rate is declining, 7-day moving averages of daily incidence of COVID-19 cases indicate ongoing community transmission [5]. Outstanding preventive measures enforced in many countries and increased load on the health care system has provoked the highest economic impact in XXI century. Decrease of the weekly economic index in the US during the pandemic period could potentially reach 12%, which was never observed within the recent decades [6]. As prevention and containment of each COVID-19 outbreak is paramount in present situation, contingency measures with experimental therapeutics are being urgently investigated for the immediate unmet medical need.
There are several monotherapies that have been tested for COVID-19. Antiviral drug remdesevir was shown to be effective for the treatment of COVID-19 in adult patients and was officially approved by FDA in the United States [7]. The speed of the COVID 19 vaccine development was [8]. Despite some critical notes [9,10], a few anti-COVID-19 vaccines were approved, and public vaccination procedure have been started. At the same time, various pharmaceuticals compounds against coronavirus infection are still being tested and yet to demonstrate their significantly high efficacy rates over placebo as well as lower mortality. Despite the recent advances in antiviral therapeutic approaches, the current pipeline for drug interventions for COVID-19 consists largely of older antivirals, immunomodulatory agents, and traditional Chinese medicines (Table 1) [11]. Hence, a multitargeted combinatorial approach may have a greater potential to be a more successful therapy than a single drug target due to the multifactorial-polygenic infectious nature of the virus. In this critical scenario, development of a novel combination of antiviral medications is a promising approach and has the potential to be completed in the near future. Currently, the most effective way of the COVID-19 treatment is the use of the antibody drug cocktails such as casirivimab and imdevimab combination (REGN-COV2 or REGEN-COV2) [12,13]. Antibody drug cocktail-based treatments may have potential immunogenicity and hypersensitivity adverse effects in patients with COVID-19 [14], making the use of small molecule pharmaceutical potentially more beneficial.
There is a need for accelerated development of effective interventions as COVID-19 confirmed cases are increasing remarkably at a fast pace. While intensive research and clinical trials to address this critical unmet need are being developed to determine the efficacy of known drugs and identify potential therapeutic targets to develop new drugs for treating COVID-19, results to date on the drug efficacy are inconclusive and inconsistent, and safety profiles are unknown in the context of the disease [15]. Recently, a number of other drugs such as combinations of ACE inhibitors and CCBs were also considered for development [7]. Their anti-inflammatory, anti-fibrotic and vasodilatory roles have been well established in the pathophysiology of pulmonary hypertension (PH), a chronic health condition recognized as a high risk factor in severe COVID-19 disease [16]. ACE inhibitors and CCBs, used in the treatment for PH, have shown therapeutic efficacy in COVID-19 [17,18], when investigated separately. Hence, it can be postulated that CCBs in combination with key proteins such as ACE, pertaining to the renin–angiotensin signaling system (RAS) family, acting indirectly via ACE, can be potential targets to combat viral entry and replication and the post-infection proinflammatory responses known as the “cytokine storm” [19].
The growing knowledge regarding SARS-CoV-2 pathophysiology provides a significant number of potential drug targets, specifically focused on the virion structure and target tissues. SARS-COV-2 is a large (27–32 kb), enveloped, positive-stranded RNA virus. Viral capsid consists of four proteins: membrane protein, envelope protein, spike protein, and nucleocapsid protein [27]. The spike protein facilitates the entry of the virus into host cells [27], and it is a critical factor of viral host range and tissue tropism and is a major inducer of host immune response and disease severity [18,28,29]. Viral spike protein is involved in receptor binding and subsequent viral entry into the host cells [28]. At least initially, it binds to the cellular angiotensin-converting enzyme 2 receptor (ACE2), cellular receptor TMPRSS2, and the calcium channel post entry in nasal secretory goblet cells, lung type II pneumocytes, and gut absorptive enterocytes [30]. Subsequently, these proteins pertaining to the RAS can be considered as potential therapeutic targets. Moreover, angiotensin II type-I receptor blockers (ARBs), as well as thiazolidinediones and ibuprofen have been reported to increase the expression of ACE2, thereby increasing the risk of infection [31]. Therefore, among the components of RAS, ACE, a zinc-metallopeptidase converting Angiotensin (Ang) I to Ang II is considered as a most promising therapeutic target. Ang II mostly exerts its activity via a type 1 and type 2 angiotensin receptor maintaining blood pressure homeostasis, and anti-inflammatory response in addition to salt and fluid balance [32]. ACE2 generates Ang (1–7) from Ang II, and then Ang (1–7) after binding and activating the mitochondrial assembly 1 (MAS) receptor broadly, shifts the balance from vasoconstriction with Ang II to vasodilation, in particularly, in pulmonary vessels [33]. The role of this vasodilatory effect in the pathogenesis of COVID-19 is not studied yet, but some animal data suggest a beneficial effect in lung disorders in PH [33]. Additionally, ACE2 and Ang (1–7) have been found to be protective in several different lung injury models [33,34]. ACE inhibitors may potentially attenuate the COVID-19 associated “cytokine storm” by upregulating ACE2, which converts Ang II to Ang (1–7) and activate MAS receptors producing beneficial vasodilatory and anti-inflammatory effects that were shown to play a potential role regarding post-infection of COVID-19 [35].
The underlying mechanism of action of CCBs on SARS-CoV-2 needs further elucidation as the underlying mechanism is not fully understood. Several pathogenic viruses have been known to induce intracellular calcium influx by hijacking predominantly the voltage-gated Ca2+ channels (VGCC) facilitating viral entry, replication, and prolonged infection period [36,37]. In previous studies with Ebola and similar viruses CCBs displayed inhibition of the replication of viruses after entry [33,34,38]. Case fatality rates were markedly reduced due to the CCB treatment among patients that were infected with Severe Fever with Thrombocytopenia Syndrome Virus (SFTSV) [39]. Similarly to SFTSV patients, CCBs are postulated to interfere with SARS-CoV-2 replication after cellular entry. CCBs may act as host-signaling targeted compounds that reduce the rate of viral mutations and interferes with the replication process via modulation of virus-hijacked host cellular machinery, compared to antivirals which target many viral proteins [40]. As a result, this could be an important factor in the development of compounds against SARS-CoV-2 [40]. The putative mechanism of action of CCBs entails the interference of the intracellular calcium influx instigated by the virus and the blocking of calcium-dependent signaling pathways pivotal to viral replication. The transient receptor potential channel (TRP) is known to be associated with hypersensitivity induced by chemical or thermal stimuli. Infection of human bronchial epithelial cells by respiratory viruses including respiratory syncytial virus (RSV), measles virus (MV) and rhinovirus (RV) was found to increase the expression of TRP channels in human bronchial epithelial cells. The over-expression of TRP proteins provides a favorable environment for propagation of virus [17,41]. Additionally, a key consequence of viral pathogenesis is marked by a strong inflammatory response preceded by a set of sequentially activated signaling pathways such as increased intracellular calcium levels leading to mitochondrial dysfunction and eventually cellular apoptosis [42,43,44]. In prior clinical outcomes, CCBs attenuated markedly the proinflammatory response by modulating the intracellular calcium levels to homeostasis in patients and reduced death rates in septic animal models with the high systemic proinflammatory state [45,46]. Additionally, the global attenuation of proinflammatory cytokines and oxidative stress by CCBs was observed in hypertensive patients compared to baseline [47]. It can be hypothesized that CCBs, besides interfering with viral replication, may attenuate systemic inflammatory responses in patients to impart the clinical benefits synergistically with their antiviral efficacy. These therapies may be particularly relevant to COVID-19 given its association with an extended proinflammatory state in patients [19].
Taken together, combinatorial approaches using ACE inhibitors, acting in an anti-inflammatory fashion, and CCBs countering the virus post-entry-stage, re-establishing Ca2+ homeostasis, and consequently down-regulating the proinflammatory signals, may impart synergistic outcomes in increasing the clinical efficacy via reducing viral load in patients compared to the individual drug treatments (Figure 1). While novel in concept, there is insufficient data on ACE/CCB combinatorial therapies available to verify whether these observations are translatable to humans, as no studies have evaluated the effects of combinations of RAS inhibitors and CCBs in COVID-19 clinical trials. To that end, we propose an in silico based computational approach for repurposing and predictions of optimal dosage and disposition of combination of ACEis and CCBs to potentiate therapeutic efficacies in clinical trials.
In the present study, the BIOiSIM was used for the computational prediction of ACE inhibitor and CCB drug dispositions in the context of tissues related to COVID-19 pathogenesis. Two CCB drugs namely Verapamil and Lacidipine belonging dihydropyridine and phenylalkylamine derivatives were selected as they cover interaction with both T- and L-types of calcium channel in cardiovascular system. Captopril, lisinopril, and captopril were chosen as the typical ACE inhibitors belonging to sulfhydryl- and dicarboxylate-containing agents showing variable interactions with ACE active centers. BIOiSIM is a dynamic, biology-driven platform that provides a scalable computational prediction of in vivo pharmacokinetic-pharmacodynamic (PK-PD) phenomena. A potential treatment strategy for COVID-19 will consist of a two-pronged approach; initially, minimizing the rate and extent of in vivo SARS-CoV-2 infection/replication and second, reducing the systemic inflammatory effect implicated in the “cytokine storm” post-infection. The aim of the study was to use the BIOiSIM platform to conduct in silico modeling of the various CCB/ACE inhibitor compounds with integrated experimental preclinical datasets and prediction of the drug disposition to tissues, which were shown to be a site of virus residence such as nasal epithelium, lungs, and intestine [48]. Investigation of the drug disposition in these tissues would help accelerate the development of the targeted anti-COVID-19 therapy.

2. Results

2.1. Sensitivity and Convergence Testing

For cases where parameter values were not available experimentally, BIOiSIM’s optimization algorithms were utilized to obtain parameter values enabling simulation of drug disposition in both the plasma venous compartment as well as the COVID-19-associated sites-of-interest. Figure 2 highlights these outcomes by showing how the objective function (cost) varies over the trials that were tested during the coarse grid search phase of optimization. There is confidence in the optimized parameter values because of the clear absolute minima that can be visualized, looking at the convergence plots for Verapamil, Spirapril, Lisinopril, Lacidipine, and Captopril (convergence point highlighted with an arrow). The final optimized values are summarized in Table 2. Table 3 gives an overview of the experimental conditions in the in vivo datasets used for the study, which were replicated for the simulation outputs.

2.2. Simulation Accuracy

BIOiSIM simulations of plasma concentration over time relative to experimental data are captured in Figure 3 [49,51,55,56,58]. The performance of the core model across the different compounds was assessed using a combination of AFE and AAFE values across the metrics generated from non-compartmental analysis and chi-squared hypothesis testing. Additionally, visual analysis of the simulation fits against the validation data was conducted to assess goodness-of-fit, and served as a key metric in assessing the platform performance in cases where error/variability data was not available from the original source. As highlighted in Table 4, the AAFE for Tmax, Cmax, and AUC0-t for all of the compounds—Lisinopril, Captopril, Spirapril, Lacidipine, and Verapamil—was less than 2, indicating that the simulation results were close to recapitulating the actual in vivo behavior. The average AFE values for the compound metrics (1.03, 0.92, and 1.08 for AUC0–t, Cmax, and Tmax, respectively) are within fold-error of ±0.1 and indicate minimal bias toward systematic over- or under-prediction of output values. Both IV and orally administered compounds were simulated with similar accuracy. Spirapril had slightly higher AAFE values for Cmax and Tmax (1.93, 1.75 respectively); however, this can be attributed to the higher variability in the in vivo datasets, as visualized in Figure 3 and confirmed with the chi-squared statistic (p-value < 0.001).
Visual analysis of the plots is also indicative of overall high accuracy, with the simulation outputs accurately captured the ascending and terminal phase that aligns with the expectation for IV and Oral routes-of-administration. For Lacidipine and Captopril, the two compounds where data was available for multiple doses, the maintained accuracy of simulation outputs (AAFE < 1.3 across all metrics) supports the model’s ability to accurately simulate clearance and absorption after oral administration of the compounds. The chi-squared metrics and associated p-values did not show statistically significant fits; as referenced in Table 2, the majority of publications reviewed for in vivo data did not contain timepoint-specific variability, and assumptions made regarding variability are likely too conservative and result in a failure to reject the null hypothesis of difference in the plots.

2.3. Simulating Distribution to Gut, Lung, and Nasal Epithelium

After successfully validating the systemic plasma PK across compounds, predictive models generated outputs of expected tissue concentration profiles in the three COVID-implicated tissues—gut, lung, and nasal epithelium (Figure 4) with calculated metrics and conditions (Table 5). The conditions that were simulated matched those used during model validation, and were assumed to be in the range of typically prescribed dosing regimens. The simulation values showed that based on the assumed permeability of nasal tissue, the overall exposure of compound in the nasal epithelium (AUC0-t) is effectively the same as in the general tissue; however, the Cmax in the nasal tissue is greater for compounds administered via IV infusion by a factor of 1.4 and 2 for Verapamil and Captopril, respectively. This behavior is likely explained by slower distribution of compound into tissues than it is noted for PO administration allowing two nasal compartments to noticeably equilibrate. Thus, it is more important to assume the permeability to the site-of-action may be an insignificant rate-limiting factor.
One interesting observation from these outcomes is that the effective lag time in Tmax between the nasal epithelium and the other tissues is minimal, which may indicate that permeability to the site-of-action is an insignificant rate-limiting factor. Additionally, the relative Cmax is <2-fold different at the nasal epithelium relative to other tissues, and indeed, the other tissues are also similar to plasma levels specifically for Spirapril, Lisinopril, and Captopril. These compounds may therefore be better candidates for rapid repurposing, given that their systemic plasma profiles are likely to have data around interaction with Ca2+ and ACE in PH patients, and those levels are directly correlated to those in tissue.

3. Discussion

The COVID-19 pandemic has been converted to a global crisis on an unprecedented scale. With synchronous and monumental efforts among researchers worldwide, several novel and old therapeutic treatments have been investigated and promoted without definitive or explicit protocols. Few therapeutic regimens may have developed at a risk; fortuitously, some demonstrate initial hope and potential efficacy. Albeit, current published results of exhaustive clinical trials are yet to be interpreted and evaluated for therapeutic impacts in patients. Several clinical studies that are focusing on the treatments of patients with COVID-19 have demonstrated an urgent unmet medical need to determine and establish the optimal treatment for therapeutic approaches currently being tested in a clinical context. To that end, we developed an in silico modeling approach in repurposing of CCBs and ACEi compounds, using an AI-integrated mechanistic modeling platform by utilizing known preclinical in vitro and in vivo datasets to accurately simulate systemic therapy disposition and site-of-action penetration of the CCBs and ACEi compounds to tissues implicated in COVID-19 pathogenesis. Numerous studies have been carried out so far confirming increased ACE expression in response to the ACEi administration. Several hypotheses exist about how increased tissue ACE2 expression may be protective rather than harmful during SARS-CoV-2 infection. For example, increased ACE2 expression may lead to enhanced sequestration of SARS-CoV-2. Moreover, ACEis’ lead to competition with Ang II for AT1R, resulting in increased Ang II to be processed by ACE2. This increases Ang (1–7) levels, which results in vasodilating and anti-fibrotic effects, providing crucial protection during coronavirus infections. Furthermore, increased binding of ACE2 to circulating Ang II could induce a conformational change resulting in less favorable binding of SARSCoV-2 to its receptor and decreased internalization of the virus when bound to ACE2 [59]. Subsequently, numerous clinical studies have shown no evidence for deleterious effects of ACEi in COVID-19 patient. In fact, discontinuing these life-saving medications potentially can have adverse effects in these groups of patients [60].
In the present study, we have demonstrated that representatives of two drug classes, namely ACE2 inhibitors and CCB, which were recently shown to provide healing effects in COVIS-19 patients [15,57], possess beneficial pharmacokinetic profile with good accumulation rate in lung and nasal epithelium tissues. Despite some recent publication confirm significant improvement of the COVID-19 patients’ condition due to monotherapy with these drugs or their combination, there are neither computation data no experimental results published confirming their suitable pharmacokinetics in the site of coronavirus potential residence i.e., lung, nala epithelium, and intestine. The simulation results presented here may serve as potential baselines for metrics designed to repurpose and prioritize compound candidates based on their PK disposition. The simulation of compound penetration to lung, gut, and nasal tissues was based on direct optimization of partition coefficients to simulate the systemic compound distribution and/or calculation using the Rodgers-Rowland method (Table 2). Thus, fine-tuned dosing optimization of drug exposure to these specific tissues could vary relative to actual clinical outcomes and the kinetic properties (e.g., IC50) of the compounds themselves. However, as the PK profile simulation accuracy is high the effective distribution and overall potency/penetration of these compounds at a certain dose can be assumed to be relatively accurate and may have utility for prioritizing compounds based on potential for reaching a desired site-of-action. Given the potential variability inherent in the methods for compound partitioning, further experimental validation in vitro or in vivo would be required to gain full confidence in the distribution simulated to these sites.
The approach taken serves as initial validation of computational tools for rapid repurposing of therapeutics with an established mechanism-of-action. Specifically, this lays the groundwork for a potential repurposing pipeline outside of COVID-19 and inclusive of unmet medical needs where existing drugs may target a pathway relevant to the condition. The general workflow of scanning a database of compounds, identifying accurately simulated datasets, and using proprietary optimization algorithms filled in knowledge gaps by optimizing missing experimental information within the bounds of physiological expectations enables the extension of simulation results beyond experimentally measured data and into case-specific prediction (Figure 5). Thus, centralized model-integrated databases well-curated datasets of in vivo compound PK are also critical for accelerated prioritization of candidates for emerging diseases. There is value in augmenting the data as selected candidates undergo further in vivo preclinical or clinical testing, specifically for building confidence in model simulation results.
Further extension and application of these results may enable greater insight into the potential therapeutic effects of ACEi and CCB combinations, as highlighted in Figure 6. One potential includes the development of PD models relating the extent of ACE inhibition to Ang II levels (postulated to play a role in the “cytokine storm” secondary to COVID-19) that could drive hypothesis-driven modeling and simulation of patient outcomes. Additionally, optimal dosing regimens of the compounds can be derived by exploring the expression levels of Ca2+ channels and ACE across the relevant tissues and extrapolating from in vitro studies of channel blocking and inhibition to assess the likelihood that a sufficient concentration of drug can reach a site-of-action to either minimize viral replication or reduce inflammatory effects. The levels of such endogenous mediators can be linked to the expected drug concentrations in plasma or active site to gain a deeper understanding of the interaction of drug with the active site components and their downstream pharmacological effects—integration of these PK-PD parameters with BIOiSIM platform can extend the prediction of drug disposition and safety and efficacy applicable to other drug therapies in healthy and diseased populations.
As more data is collected, and the pathophysiology of COVID-19 in healthy and diseased populations is understood, more customized prescriptions of specific CCB/ACEi combinations can be derived by looking at susceptibility to drug–drug interaction and variability in relevant gene expression levels between subjects. Overall, these results give additional confidence in the BIOiSIM platform’s ability to rapidly identify and simulate drug disposition for compounds that may be efficacious therapies for rapidly emerging, deadly conditions. Considering the universality ingrained in the core of the BIOiSIM platform, its predictive potential could be implemented across any modalities for repurposing and consequent treatment of infectious diseases as well as other clinical conditions.

4. Materials and Methods

4.1. Overview of BIOiSIM Platform

The core functionality of the in silico simulation platform described here was outlined in a recent work around modeling penetration of transdermal formulations to systemic circulation [47]. Briefly, the platform is comprised of a 16-organ model of compound PK, validated across an internal database of small molecule compounds. Auxiliary models (specific PK or PD) can be integrated into the centralized framework to expand the ability to make specific predictions of compound disposition. Artificial intelligence-Machine learning (AI-ML) algorithms are utilized to either optimize missing parameters in the case of insufficient experimental data, or as predictive solutions to train on existing in vivo/in vitro datasets. The software systems are hosted on Amazon Web Services (AWS) cloud, enable high throughput of parallelized PK simulations (PK simulator performance is ~0.08 s to simulate 1 hr of drug exposure at 0.36 s resolution) which is available in the VeriSIM Life customer portal.
The drug-dependent parameters used in the model are either experimentally determined or predicted/optimized using a combination of ML cost minimization algorithms; these include integrating iterative optimization algorithms with random walk methods (similar to Markov Chain/Monte Carlo) to converge on a global minimum. Details of the approach are discussed in our recent publication [61].

4.2. Modeling Compound Disposition

Computing the mass balance between the compartments follows the general form:
V o r g a n d C o r g a n d t = Q o r g a n ( C b l o o d , i n C o r g a n K T : P , u n b o u n d f u n b o u n d , p l a s m a B : P )
For metabolizing and eliminating organs (e.g., liver, kidney) the general relationship in the model is defined as:
V o r g a n d C o r g a n d t = Q o r g a n ( C b l o o d , i n C o r g a n K T : P , u n b o u n d f u n b o u n d , p l a s m a B : P C L o r g a n C o r g a n )
where tissue-dependent parameters are expressed as V (organ volume), Q (flow rate), CL (organ-level clearance), and K (unbound tissue:plasma partition coefficient); B:P represents whole blood to plasma partition ratio, and C is the drug concentration in the specific compartment.
This global approach has been previously described in our manuscript [61].
Given the identified relevant tissues for SARS-CoV-2 infection (nasal, lung, and intestine), a critical aspect of the approach was enabling simulation of compound disposition to those relevant sites-of-action. Generally, the tissues were assumed to be well-mixed; thus, distribution was characterized using a partition coefficient:
K T i s s u e : P l a s m a = C t i s s u e C p l a s m a
where Ctissue is concentration of the compound-of-interest in a particular organ/tissue, and KTissue:Plasma is the drug/subject-specific partition coefficient in the tissue. The methodology utilized to determine the partition coefficient was derived from the simulation accuracy of the plasma-concentration curve. There are multiple different methodologies that have been utilized in an attempt to predict these partition coefficients accurately, including industry standards such as the Rodgers-Rowland and Schmitt equations. However, even these relationships are susceptible to significant variability [62,63,64]. These models predict compound partitioning as a function of logP (octanol-water partition coefficient), pKa, and plasma protein binding (fu,p); however, for specific groups of compounds, they are found to underperform in their predictive capabilities based on the simplifying assumptions used [62,63,64,65,66]. Therefore, we utilized the standard Rodgers-Rowland equation, optimization of the octanol-water partition coefficient, and direct optimization of an average Kp to increase the simulation accuracy for in vivo PK. Partition coefficient values were directly optimized for Captopril, Lacidipine, and Verapamil, while octanol-water partition values (logP) were optimized for Spirapril.
The nasal epithelium mechanism for compound distribution involved integration with a permeation across a barrier. Given the location of nasal goblet cells in vivo, compound into the epithelium was modeled with assumptions of 1D flux across the barrier derived from Fick’s Law. To approximate compound buildup in tissue and effective reduction of the concentration gradient without incorporating semi-infinite sink assumptions, the epithelium thickness was scaled down by a factor of two. To model penetration to the nasal epithelium, a model of small molecule penetration across a permeable barrier was utilized. The developed model was identified and adapted from a combination of approaches utilized for predicting penetration across the barrier:
d C N a d t = Q N a ( C a   B P C N a k p n a   ) P n a s a l A   ( C N a K p N a   C n a s a l _ e p i f n a s a l   )
d C n a s a l _ e p i d t = P n a s a l A s u r f ( C N a K p , N a   C n a s a l _ e p i f n a s a l   )
where Na refers to the nasal tissue compartment, Cnasal_epi the concentration in the nasal epithelium, fnasal the fraction of compound bound to the tissue, Pnasal the permeability across the barrier, and A the area across which compound diffusion can occur. Physiological parameters used in the model were aggregated from a diversity of sources [65,66,67]. For the purpose of these simulations, permeability was assumed to be 0.000101 cm/s for all of the compounds studied, and the partition coefficient into the nasal tissue was assumed to be equivalent to that of lung tissue, as reported in Table 2 [68,69,70].

4.3. Rapid Repurposing Workflow

4.3.1. Test Dataset

An internal database largest proprietary curated database consisting of structure-related data for > 1 M compounds, >3700+ unique in vivo plasma concentration-time validation datasets from public and proprietary sources (signed data sharing partnerships) representing ~2000 unique compounds and 83 different subject populations (different species, gender, strain, sub-strain).
The dataset was compiled through a combination of automated data scraping and manual data curation, resulting in the development of over 41 coded/automated consistency checks that detect outliers and corrupt data coming from either published literature on FDA compounds or from ongoing data partnerships. This innovative data curation approach was applied to CCB/ACE data curation, import and validation.
Screening through the internal database of compounds, comprising >2000 small molecules PK parameters & in vivo datasets identified five compounds that targeted the pathway-of-interest and had sufficient validation data available for confident & rapid simulation of outcomes. To establish the combinatorial effect of ACEis/CCBs combinations on COVID-19 infection, an accurate prediction of drug concentration at a site-of-action is required. Figure 2 highlights the approach taken to identify the 5 drugs (3 ACEis, 2 CCBs) that were candidates for the study as a result of the available in vivo clinical data and preclinical PK parameters that were already established. to validate the ability to predict plasma concentration accurately (as the standard metric for PK).

4.3.2. Statistics and Tools

The statistical methodologies utilized for the data analysis have been detailed in a previous work [61]. Briefly, in vivo plasma concentration datasets and associated error bars, when available, were manually digitized from source publications using “WebPlotDigitizer” version 4.2.34. Model development and validation were done using the in-house platform in Python with Cython integration; matplotlib (v2.0.2) and Numpy (v1.14.2) were auxiliary packages used in simulation deployment and analysis. Model validation and analysis of model goodness-of-fit/accuracy were conducted using three quantitative metrics: absolute average fold error (AAFE), average fold-error (AFE), and chi-squared statistic (Χ2) with associated p-value (null hypothesis defined as no difference in predicted vs. experimental measurements). Sensitivity of the model was evaluated using convergence plots generated during optimization of the parameters. Non-compartmental calculations were utilized to compare the accuracy of the simulations to the experimental data; AAFE and AFE were utilized to evaluate the accuracy of the PK outputs AUC0–t, Cmax, and tmax using the general equations:
A F E = A v e r a g e   f o l d   e r r o r = 10 1 n i = 1 n log ( p r e d i c t e d i o b s e r v e d i )
A A F E = A b s o l u t e   a v e r a g e   f o l d   e r r o r = 10 1 n i = 1 n | log ( p r e d i c t e d i o b s e r v e d i ) |
where n is the total number of compounds used in the analysis and Predictedi/Observedi correspond to predicted and observed values of PK parameters, respectively. χ2 was calculated using the relationship:
χ 2 = 1 n i = 1 n ( p r e d i c t e d i o b s e r v e d i o b s e r v e d e r r o r ,   i ) 2  
where observed error is the standard error/deviation in the measurements in the individual experimental data timepoints, obtained by digitizing the error bars from the respective publications. Optimization convergence was driven and measured by AAFE of AUC, Cmax, and tmax. Statistical calculations of AFE, AAFE, and visual plot analysis were done in GraphPad Prism version 8.4.1 (GraphPad Software, San Diego, CA, USA) and Microsoft Excel (2016).

4.3.3. Subjects

Species-specific parameters used in the simulations were adapted from literature sources and included parameters such as organ flow rates, composition, volumes, and protein levels for the physiological compartments in the BIOiSIM model [71,72,73,74,75]. Parameters for the model of small molecule diffusion across the nasal tissue and epithelium were obtained from a thorough literature review [65,66,76,77].

5. Conclusions

As cumulative datasets are collected, and the pathophysiology of COVID-19 in healthy and diseased populations is understood, more customized prescriptions of specific CCB/ACEi combinations can be derived by looking at susceptibility to drug–drug interaction and the associated variability in relevant gene expression levels between subjects. The findings of the present study confirm the potential use of a combinatorial approach of CCB and ACEi as therapeutic agents against COVID-19 infection due to their favorable tissue distribution with sufficient level of accumulation in the gut, lung, and airway epithelium. These organs have been demonstrated to be sites of virus infection, residence, and replication. Therefore, the combinatorial therapy using these drug classes can provide therapeutic efficacy against COVD-19. Overall, these results evaluate a new paradigm in using AI/ML-driven computational modeling for repurposing and accelerating the drug development process for swiftly identifying and simulating drug disposition for compounds that may be effective therapies for rapidly emerging, deadly clinical conditions.

Author Contributions

Conceptualization and writing original draft, K.C., V.G.A., J.V.; Writing review, editing and data curation, M.K., N.M., H.H., A.J.; review and editing software, N.G.; validation and editing, Y.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors are grateful to Alexander Harmsen for useful discussions, and his help in reviewing and managing the writing process.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Nicola, M.; O’Neill, N.; Sohrabi, C.; Khan, M.; Agha, M.; Agha, R. Evidence based management guideline for the COVID-19 pandemic—Review article. Int. J. Surg. 2020, 77, 206–216. [Google Scholar] [CrossRef] [PubMed]
  2. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef] [PubMed]
  3. WHO. WHO Coronavirus Disease (COVID-19) Dashboard; WHO: Geneva, Switzerland, 2021. [Google Scholar]
  4. Rehman, M.F.U.; Fariha, C.; Anwar, A.; Shahzad, N.; Ahmad, M.; Mukhtar, S.; Ul Haque, F.M. Novel coronavirus disease (COVID-19) pandemic: A recent mini review. Comput. Struct. Biotechnol. J. 2021, 19, 612–623. [Google Scholar] [CrossRef] [PubMed]
  5. Ledford, H. Why do COVID death rates seem to be falling? Nature 2020, 587, 190–192. [Google Scholar] [CrossRef] [PubMed]
  6. Mosser, P.C. Central bank responses to COVID-19. Bus. Econ. 2020, 55, 191–201. [Google Scholar] [CrossRef] [PubMed]
  7. Guy, R.K.; DiPaola, R.S.; Romanelli, F.; Dutch, R.E. Rapid repurposing of drugs for COVID-19. Science 2020, 368, 829–830. [Google Scholar] [CrossRef] [PubMed]
  8. Heaton, P.M. The Covid-19 Vaccine-Development Multiverse. N. Engl. J. Med. 2020, 383, 1986–1988. [Google Scholar] [CrossRef] [PubMed]
  9. Lebrasseur, A.; Fortin-Bédard, N.; Lettre, J.; Bussières, E.-L.; Best, K.; Boucher, N.; Hotton, M.; Beaulieu-Bonneau, S.; Mercier, C.; Lamontagne, M.-E.; et al. Impact of COVID-19 on people with physical disabilities: A rapid review. Disabil. Health J. 2020, 14, 101014. [Google Scholar] [CrossRef]
  10. McMasters, M.; Blair, B.M.; Lazarus, H.M.; Alonso, C.D. Casting a wider protective net: Anti-infective vaccine strategies for patients with hematologic malignancy and blood and marrow transplantation. Blood Rev. 2020, 100779. [Google Scholar] [CrossRef]
  11. Zhao, M.; Zhang, J.; Li, H.; Luo, Z.; Ye, J.; Xu, Y.; Wang, Z.; Ye, D.; Liu, J.; Li, D.; et al. Recent progress of antiviral therapy for coronavirus disease 2019. Eur. J. Pharmacol. 2020, 890, 173646. [Google Scholar] [CrossRef]
  12. Baum, A.; Ajithdoss, D.; Copin, R.; Zhou, A.; Lanza, K.; Negron, N.; Ni, M.; Wei, Y.; Mohammadi, K.; Musser, B.; et al. REGN-COV2 antibodies prevent and treat SARS-CoV-2 infection in rhesus macaques and hamsters. Science 2020, 370, 1110–1115. [Google Scholar] [CrossRef]
  13. Weinreich, D.M.; Sivapalasingam, S.; Norton, T.; Ali, S.; Gao, H.; Bhore, R.; Musser, B.J.; Soo, Y.; Rofail, D.; Im, J.; et al. REGN-COV2, a Neutralizing Antibody Cocktail, in Outpatients with Covid-19. N. Engl. J. Med. 2021, 384, 238–251. [Google Scholar] [CrossRef]
  14. DeFrancesco, L. COVID-19 antibodies on trial. Nat. Biotechnol. 2020, 38, 1242–1252. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, J.; Xie, B.; Hashimoto, K. Current status of potential therapeutic candidates for the COVID-19 crisis. Brain Behav. Immun. 2020, 87, 59–73. [Google Scholar] [CrossRef] [PubMed]
  16. Badagliacca, R.; Sciomer, S.; Petrosillo, N. Endothelin receptor antagonists for pulmonary arterial hypertension and COVID-19: Friend or foe? J. Heart Lung Transplant. 2020, 39, 729–730. [Google Scholar] [CrossRef]
  17. Chen, X.; Cao, R.; Zhong, W. Host Calcium Channels and Pumps in Viral Infections. Cells 2019, 9, 94. [Google Scholar] [CrossRef] [Green Version]
  18. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Kruger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280.e8. [Google Scholar] [CrossRef] [PubMed]
  19. Ragab, D.; Salah Eldin, H.; Taeimah, M.; Khattab, R.; Salem, R. The COVID-19 Cytokine Storm; What We Know So Far. Front Immunol. 2020, 11, 1446. [Google Scholar] [CrossRef]
  20. Babaei, F.; Mirzababaei, M.; Nassiri-Asl, M.; Hosseinzadeh, H. Review of registered clinical trials for the treatment of COVID-19. Drug Dev. Res. 2020. [Google Scholar] [CrossRef] [PubMed]
  21. Luo, H.; Yang, M.; Tang, Q.L.; Hu, X.Y.; Willcox, M.L.; Liu, J.P. Characteristics of registered clinical trials on traditional Chinese medicine for coronavirus disease 2019 (COVID-19): A scoping review. Eur. J. Integr. Med. 2020, 41, 101251. [Google Scholar] [CrossRef]
  22. Bakovic, A.; Risner, K.; Bhalla, N.; Alem, F.; Chang, T.L.; Weston, W.; Harness, J.A.; Narayanan, A. Brilacidin, a COVID-19 Drug Candidate, Exhibits Potent In Vitro Antiviral Activity Against SARS-CoV-2. bioRxiv 2020. [Google Scholar] [CrossRef]
  23. Omrani, A.S.; Pathan, S.A.; Thomas, S.A.; Harris, T.R.E.; Coyle, P.V.; Thomas, C.E.; Qureshi, I.; Bhutta, Z.A.; Mawlawi, N.A.; Kahlout, R.A.; et al. Randomized double-blinded placebo-controlled trial of hydroxychloroquine with or without azithromycin for virologic cure of non-severe Covid-19. EClinicalMedicine 2020, 29. [Google Scholar] [CrossRef]
  24. Breining, P.; Frolund, A.L.; Hojen, J.F.; Gunst, J.D.; Staerke, N.B.; Saedder, E.; Cases-Thomas, M.; Little, P.; Nielsen, L.P.; Sogaard, O.S.; et al. Camostat mesylate against SARS-CoV-2 and COVID-19-Rationale, dosing and safety. Basic Clin. Pharmacol. Toxicol. 2020, 128, 204–212. [Google Scholar] [CrossRef] [PubMed]
  25. Takahashi, W.; Yoneda, T.; Koba, H.; Ueda, T.; Tsuji, N.; Ogawa, H.; Asakura, H. Potential mechanisms of nafamostat therapy for severe COVID-19 pneumonia with disseminated intravascular coagulation. Int. J. Infect. Dis. 2021, 102, 529–531. [Google Scholar] [CrossRef] [PubMed]
  26. Dabbous, H.M.; Abd-Elsalam, S.; El-Sayed, M.H.; Sherief, A.F.; Ebeid, F.F.S.; El Ghafar, M.S.A.; Soliman, S.; Elbahnasawy, M.; Badawi, R.; Tageldin, M.A. Efficacy of favipiravir in COVID-19 treatment: A multi-center randomized study. Arch. Virol. 2021, 166, 949–954. [Google Scholar] [CrossRef] [PubMed]
  27. Wrapp, D.; Wang, N.; Corbett, K.S.; Goldsmith, J.A.; Hsieh, C.L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 2020, 367, 1260–1263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Li, F. Structure, Function, and Evolution of Coronavirus Spike Proteins. Annu. Rev. Virol. 2016, 3, 237–261. [Google Scholar] [CrossRef] [Green Version]
  29. Zhou, P.; Yang, X.L.; Wang, X.G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.R.; Zhu, Y.; Li, B.; Huang, C.L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Ziegler, C.G.K.; Allon, S.J.; Nyquist, S.K.; Mbano, I.M.; Miao, V.N.; Tzouanas, C.N.; Cao, Y.; Yousif, A.S.; Bals, J.; Hauser, B.M.; et al. SARS-CoV-2 Receptor ACE2 Is an Interferon-Stimulated Gene in Human Airway Epithelial Cells and Is Detected in Specific Cell Subsets across Tissues. Cell 2020, 181, 1016–1035. [Google Scholar] [CrossRef] [PubMed]
  31. Kai, H.; Kai, M. Interactions of coronaviruses with ACE2, angiotensin II, and RAS inhibitors—Lessons from available evidence and insights into COVID-19. Hypertens. Res. 2020, 43, 648–654. [Google Scholar] [CrossRef]
  32. Baudin, B. New aspects on angiotensin-converting enzyme: From gene to disease. Clin. Chem. Lab. Med. 2002, 40, 256–265. [Google Scholar] [CrossRef]
  33. Rice, G.I.; Thomas, D.A.; Grant, P.J.; Turner, A.J.; Hooper, N.M. Evaluation of angiotensin-converting enzyme (ACE), its homologue ACE2 and neprilysin in angiotensin peptide metabolism. Biochem. J. 2004, 383, 45–51. [Google Scholar] [CrossRef] [PubMed]
  34. Kuba, K.; Imai, Y.; Rao, S.; Gao, H.; Guo, F.; Guan, B.; Huan, Y.; Yang, P.; Zhang, Y.; Deng, W.; et al. A crucial role of angio-tensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury. Nat. Med. 2005, 11, 875–879. [Google Scholar] [CrossRef] [PubMed]
  35. Patel, A.B.; Verma, A. COVID-19 and Angiotensin-Converting Enzyme Inhibitors and Angiotensin Receptor Blockers: What Is the Evidence? JAMA 2020, 323, 1769–1770. [Google Scholar] [CrossRef] [PubMed]
  36. Donate-Macian, P.; Jungfleisch, J.; Perez-Vilaro, G.; Rubio-Moscardo, F.; Peralvarez-Marin, A.; Diez, J.; Valverde, M.A. The TRPV4 channel links calcium influx to DDX3X activity and viral infectivity. Nat. Commun. 2018, 9, 2307. [Google Scholar] [CrossRef]
  37. Ueda, M.; Daidoji, T.; Du, A.; Yang, C.S.; Ibrahim, M.S.; Ikuta, K.; Nakaya, T. Highly pathogenic H5N1 avian influenza virus induces extracellular Ca2+ influx, leading to apoptosis in avian cells. J. Virol. 2010, 84, 3068–3078. [Google Scholar] [CrossRef] [Green Version]
  38. Nathan, L.; Lai, A.L.; Millet, J.K.; Straus, M.R.; Freed, J.H.; Whittaker, G.R.; Daniel, S. Calcium Ions Directly Interact with the Ebola Virus Fusion Peptide To Promote Structure–Function Changes That Enhance Infection. ACS Infect. Dis. 2020, 6, 250–260. [Google Scholar] [CrossRef] [PubMed]
  39. Li, H.; Zhang, L.K.; Li, S.F.; Zhang, S.F.; Wan, W.W.; Zhang, Y.L.; Xin, Q.L.; Dai, K.; Hu, Y.Y.; Wang, Z.B.; et al. Calcium channel blockers reduce severe fever with thrombocytopenia syndrome virus (SFTSV) related fatality. Cell Res. 2019, 29, 739–753. [Google Scholar] [CrossRef] [PubMed]
  40. Lavanya, M.; Cuevas, C.D.; Thomas, M.; Cherry, S.; Ross, S.R. siRNA screen for genes that affect Junin virus entry uncovers voltage-gated calcium channels as a therapeutic target. Sci. Transl. Med. 2013, 5, 204ra131. [Google Scholar] [CrossRef] [Green Version]
  41. Omar, S.; Clarke, R.; Abdullah, H.; Brady, C.; Corry, J.; Winter, H.; Touzelet, O.; Power, U.F.; Lundy, F.; McGarvey, L.P.A.; et al. Respiratory virus infection up-regulates TRPV1, TRPA1 and ASICS3 receptors on airway cells. PLoS ONE 2017, 12, e0171681. [Google Scholar] [CrossRef] [Green Version]
  42. Horng, T. Calcium signaling and mitochondrial destabilization in the triggering of the NLRP3 inflammasome. Trends Immunol. 2014, 35, 253–261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Maher, P.; van Leyen, K.; Dey, P.N.; Honrath, B.; Dolga, A.; Methner, A. The role of Ca2+ in cell death caused by oxidative glutamate toxicity and ferroptosis. Cell Calcium 2018, 70, 47–55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Yaron, J.R.; Gangaraju, S.; Rao, M.Y.; Kong, X.; Zhang, L.; Su, F.; Tian, Y.; Glenn, H.L.; Meldrum, D.R. K+ regulates Ca2+ to drive inflammasome signaling: Dynamic visualization of ion flux in live cells. Cell Death Dis. 2015, 6, e1954. [Google Scholar] [CrossRef] [Green Version]
  45. D’Elia, J.A.; Weinrauch, L.A. Calcium Ion Channels: Roles in Infection and Sepsis Mechanisms of Calcium Channel Blocker Benefits in Immunocompromised Patients at Risk for Infection. Int. J. Mol. Sci. 2018, 19, 2465. [Google Scholar] [CrossRef] [Green Version]
  46. Silva, I.V.G.; de Figueiredo, R.C.; Rios, D.R.A. Effect of Different Classes of Antihypertensive Drugs on Endothelial Function and Inflammation. Int. J. Mol. Sci. 2019, 20, 3458. [Google Scholar] [CrossRef] [Green Version]
  47. Kim, H.J.; Han, S.J.; Kim, D.J.; Jang, H.C.; Lim, S.; Choi, S.H.; Kim, Y.H.; Shin, D.H.; Kim, S.H.; Kim, T.H.; et al. Effects of valsartan and amlodipine on oxidative stress in type 2 diabetic patients with hypertension: A randomized, multicenter study. Korean J. Intern. Med. 2017, 32, 497–504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Tang, H.; Lu, X.; Qie, S.; Xi, J. Thoughts on detecting tissue distribution of potential COVID-19 receptors. Future Virol. 2020, 15, 489–496. [Google Scholar] [CrossRef]
  49. Lombardo, F.; Berellini, G.; Obach, R.S. Trend Analysis of a Database of Intravenous Pharmacokinetic Parameters in Humans for 1352 Drug Compounds. Drug Metab. Dispos. 2018, 46, 1466–1477. [Google Scholar] [CrossRef] [Green Version]
  50. Benet, L.Z.; Broccatelli, F.; Oprea, T.I. BDDCS applied to over 900 drugs. AAPS J. 2011, 13, 519–547. [Google Scholar] [CrossRef] [Green Version]
  51. Krahenbuhl, S.; Grass, P.; Surve, A.; Kutz, K.; Reichen, J. Pharmacokinetics and haemodynamic effects of a single oral dose of the novel ACE inhibitor spirapril in patients with chronic liver disease. Eur. J. Clin. Pharmacol. 1993, 45, 247–253. [Google Scholar] [CrossRef]
  52. Yang, B.; Wu, C.; Ji, B.; Wu, M.; He, Z.; Shang, L.; Sun, J. Virtual population pharmacokinetic using physiologically based pharmacokinetic model for evaluating bioequivalence of oral lacidipine formulations in dogs. Asian J. Pharm. Sci. 2017, 12, 98–104. [Google Scholar] [CrossRef] [Green Version]
  53. Lee, C.R.; Bryson, H.M. Lacidipine. A review of its pharmacodynamic and pharmacokinetic properties and therapeutic po-tential in the treatment of hypertension. Drugs 1994, 48, 274–296. [Google Scholar] [CrossRef] [PubMed]
  54. Ye, M.; Nagar, S.; Korzekwa, K. A physiologically based pharmacokinetic model to predict the pharmacokinetics of highly protein-bound drugs and the impact of errors in plasma protein binding. Biopharm. Drug Dispos. 2016, 37, 123–141. [Google Scholar] [CrossRef]
  55. Eichelbaum, M.; Mikus, G.; Vogelgesang, B. Pharmacokinetics of (+)-, (-)- and (+/-)-verapamil after intravenous administration. Br. J. Clin. Pharmacol. 1984, 17, 453–458. [Google Scholar] [CrossRef] [Green Version]
  56. Beermann, B. Pharmacokinetics of lisinopril. Am. J. Med. 1988, 85, 25–30. [Google Scholar] [CrossRef]
  57. Albashir, A.A.D. Renin-Angiotensin-Aldosterone System (RAAS) Inhibitors and Coronavirus Disease 2019 (COVID-19). South Med. J. 2021, 114, 51–56. [Google Scholar] [CrossRef] [PubMed]
  58. Da Ros, L.; Squassante, L.; Milleri, S. Dose linearity of lacidipine pharmacokinetics after single and repeated oral doses in healthy volunteers. Clin. Pharmacokinet. 2003, 42, 99–106. [Google Scholar] [CrossRef]
  59. Bourgonje, A.R.; Abdulle, A.E.; Timens, W.; Hillebrands, J.-L.; Navis, G.J.; Gordijn, S.J.; Bolling, M.C.; Dijkstra, G.; Voors, A.A.; Osterhaus, A.D.; et al. Angiotensin-converting enzyme 2 (ACE2), SARS-CoV-2 and the pathophysiology of coronavirus disease 2019 (COVID-19). J. Pathol. 2020, 251, 228–248. [Google Scholar] [CrossRef] [PubMed]
  60. Rossi, G.P.; Sanga, V.; Barton, M. Potential harmful effects of discontinuing ACE-inhibitors and ARBs in COVID-19 patients. eLife 2020, 9, e57278. [Google Scholar] [CrossRef] [Green Version]
  61. Maharao, N.; Antontsev, V.; Hou, H.; Walsh, J.; Varshney, J. Scalable in silico Simulation of Transdermal Drug Permeability: Application of BIOiSIM Platform. Drug Des. Devel. Ther. 2020, 14, 2307–2317. [Google Scholar] [CrossRef] [PubMed]
  62. Rodgers, T.; Leahy, D.; Rowland, M. Physiologically based pharmacokinetic modeling 1: Predicting the tissue distribution of moderate-to-strong bases. J. Pharm. Sci. 2005, 94, 1259–1276. [Google Scholar] [CrossRef] [PubMed]
  63. Rodgers, T.; Rowland, M. Physiologically based pharmacokinetic modelling 2: Predicting the tissue distribution of acids, very weak bases, neutrals and zwitterions. J. Pharm. Sci. 2006, 95, 1238–1257. [Google Scholar] [CrossRef]
  64. Schmitt, W. General approach for the calculation of tissue to plasma partition coefficients. Toxicol. In Vitro 2008, 22, 457–467. [Google Scholar] [CrossRef] [PubMed]
  65. Andersen, M.; Sarangapani, R.; Gentry, R.; Clewell, H.; Covington, T.; Frederick, C.B. Application of a hybrid CFD-PBPK nasal dosimetry model in an inhalation risk assessment: An example with acrylic acid. Toxicol. Sci. 2000, 57, 312–325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Boger, E.; Evans, N.; Chappell, M.; Lundqvist, A.; Ewing, P.; Wigenborg, A.; Friden, M. Systems Pharmacology Approach for Prediction of Pulmonary and Systemic Pharmacokinetics and Receptor Occupancy of Inhaled Drugs. CPT Pharmacomet. Syst. Pharmacol. 2016, 5, 201–210. [Google Scholar] [CrossRef] [PubMed]
  67. Frederick, C.B.; Bush, M.L.; Lomax, L.G.; Black, K.A.; Finch, L.; Kimbell, J.S.; Morgan, K.T.; Subramaniam, R.P.; Morris, J.B.; Ultman, J.S. Application of a hybrid computational fluid dynamics and physiologically based inhalation model for interspecies dosimetry extrapolation of acidic vapors in the upper airways. Toxicol. Appl. Pharmacol. 1998, 152, 211–231. [Google Scholar] [CrossRef]
  68. Jia, L.; Wong, H. In vitro and in vivo assessment of cellular permeability and pharmacodynamics of S-nitrosylated captopril, a nitric oxide donor. Br. J. Pharmacol. 2001, 134, 1697–1704. [Google Scholar] [CrossRef]
  69. O’Hagan, S.; Kell, D.B. The apparent permeabilities of Caco-2 cells to marketed drugs: Magnitude, and independence from both biophysical properties and endogenite similarities. PeerJ 2015, 3, e1405. [Google Scholar] [CrossRef] [Green Version]
  70. Wu, C.; Kou, L.; Ma, P.; Gao, L.; Li, B.; Li, R.; Luo, C.; Shentu, J.; He, Z.; Sun, J. Interspecies prediction of oral pharmacokinetics of different lacidipine formulations from dogs to human: Physiologically based pharmacokinetic modelling combined with bi-orelevant dissolution. RSC Adv. 2015, 5, 19844–19852. [Google Scholar] [CrossRef]
  71. Integrated Physiology Database—PK-Sim®; Bayer Technology Services GmbH: Leverkusen, Germany, 2012; Available online: http://www.systems-biology.com/products/pk-sim.html (accessed on 10 October 2020).
  72. US Environmental Protection Agency Office of Health and Environmental Assessment. Physiological Parameter Values for PBPK Models: A Report Prepared by the International Life Sciences Institute and Risk Sciences Institute; International Life Sciences Institute and Risk Sciences Institute: Washington, DC, USA, 1994.
  73. ICRP. Basic Anatomical and Physiological Data for Use in Radiological Protection Reference Values. Ann. ICRP 2002, 32, 3–4. [Google Scholar]
  74. Maharao, N.; Antontsev, V.; Wright, M.; Varshney, J. Entering the era of computationally driven drug development. Drug Metab. Rev. 2020, 52, 283–298. [Google Scholar] [CrossRef] [PubMed]
  75. Peters, S. Physiologically-Based Pharmacokinetic (PBPK) Modeling and Simulations: Principles, Methods, and Applications in the Pharmaceutical Industry; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012. [Google Scholar] [CrossRef]
  76. Campbell, J.L.; Andersen, M.E.; Clewell, H.J. A hybrid CFD-PBPK model for naphthalene in rat and human with IVIVE for nasal tissue metabolism and cross-species dosimetry. Inhal. Toxicol. 2014, 26, 333–344. [Google Scholar] [CrossRef] [PubMed]
  77. Sweeney, L.M.; Kirman, C.R.; Gargas, M.L.; Carson, M.L.; Tardiff, R.G. Development of a physiologically-based toxicokinetic model of acrylamide and glycidamide in rats and humans. Food Chem. Toxicol. 2010, 48, 668–685. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Potential model for the combinatorial actions of angiotensin-converting enzyme inhibitors (ACEi) and calcium channel blockers (CCBs) on SARS-CoV-2 infection, replication, and proinflammatory response. ACEi acting indirectly through ACE2 and the mitochondrial assembly 1 (MAS), ATR2 axis imparts an anti-inflammatory and anti-fibrotic response, while CCBs acts at various steps, restores intracellular Ca2+ flux, consequently inhibiting post-infection virus internalization and genome replication. Both ACEi and CCBs potentially can induce a synergistic anti-inflammatory and anti-fibrotic response in the attenuation of the “cytokine storm”.
Figure 1. Potential model for the combinatorial actions of angiotensin-converting enzyme inhibitors (ACEi) and calcium channel blockers (CCBs) on SARS-CoV-2 infection, replication, and proinflammatory response. ACEi acting indirectly through ACE2 and the mitochondrial assembly 1 (MAS), ATR2 axis imparts an anti-inflammatory and anti-fibrotic response, while CCBs acts at various steps, restores intracellular Ca2+ flux, consequently inhibiting post-infection virus internalization and genome replication. Both ACEi and CCBs potentially can induce a synergistic anti-inflammatory and anti-fibrotic response in the attenuation of the “cytokine storm”.
Molecules 26 01912 g001
Figure 2. Convergence onto discrete parameter combinations during coarse optimization for (A) Verapamil, (B) Spirapril, (C) Lisinopril, (D) Lacidipine, and (E) Captopril. Relative cost on the y-axis is calculated as: Cost_relative = AAFECmax*AAFEtmax*AAFEAUC where AAFECmax, AAFEtmax, AAFEAUC—absolute average fold error for Cmax, Tmax and AUC predictions.
Figure 2. Convergence onto discrete parameter combinations during coarse optimization for (A) Verapamil, (B) Spirapril, (C) Lisinopril, (D) Lacidipine, and (E) Captopril. Relative cost on the y-axis is calculated as: Cost_relative = AAFECmax*AAFEtmax*AAFEAUC where AAFECmax, AAFEtmax, AAFEAUC—absolute average fold error for Cmax, Tmax and AUC predictions.
Molecules 26 01912 g002
Figure 3. Dose and route-of-administration dependent prediction of compound plasma concentration for (A) Verapamil (i.v.), (B) Spirapril (p.o.), (C) Lisinopril (p.o.), (D) Lacidipine (p.o.), and (E) Captopril (i.v.). Dashed lines correspond to BIOiSIM simulation outputs. Error bars and individual data points were digitized from the original publications (49,51,55,56,58) and correspond to standard error/standard deviation, as presented in the original work.
Figure 3. Dose and route-of-administration dependent prediction of compound plasma concentration for (A) Verapamil (i.v.), (B) Spirapril (p.o.), (C) Lisinopril (p.o.), (D) Lacidipine (p.o.), and (E) Captopril (i.v.). Dashed lines correspond to BIOiSIM simulation outputs. Error bars and individual data points were digitized from the original publications (49,51,55,56,58) and correspond to standard error/standard deviation, as presented in the original work.
Molecules 26 01912 g003
Figure 4. Prediction of the site-specific distribution of (A) Verapamil (i.v.), (B) Spirapril (p.o.), (C) Lisinopril (p.o.), (D) Lacidipine (p.o.), and (E) Captopril (i.v.) to tissues implicated in COVID-19 replication and pathogenesis. Subject conditions are described in Table 3.
Figure 4. Prediction of the site-specific distribution of (A) Verapamil (i.v.), (B) Spirapril (p.o.), (C) Lisinopril (p.o.), (D) Lacidipine (p.o.), and (E) Captopril (i.v.) to tissues implicated in COVID-19 replication and pathogenesis. Subject conditions are described in Table 3.
Molecules 26 01912 g004
Figure 5. Metrics for candidate compound identification unitilizing internal database of small molecule PK. The numbers in blue boxes correspond to the number of compounds identified as acceptable during each stage of curatrion.
Figure 5. Metrics for candidate compound identification unitilizing internal database of small molecule PK. The numbers in blue boxes correspond to the number of compounds identified as acceptable during each stage of curatrion.
Molecules 26 01912 g005
Figure 6. Diagram of modeling approach taken that extends beyond pharmacokinetic modeling of compound disposition and can explore more complex viral and combination effects.
Figure 6. Diagram of modeling approach taken that extends beyond pharmacokinetic modeling of compound disposition and can explore more complex viral and combination effects.
Molecules 26 01912 g006
Table 1. Select compounds currently in clinical trials for treating COVID-19.
Table 1. Select compounds currently in clinical trials for treating COVID-19.
Compound NameDevelopmentOverviewSource
RemdesivirGilead, Approved by FDAAntiviral, host factor-targeted. RNA-dependent/RNA polymerase-targeted[20]
APN01APEIRON Biologics, Phase IPilot trial ongoing in China[21]
BrilacidinInnovation Pharmaceuticals, Phase IIDefensin mimetic drug candidate. Has shown antibacterial, anti-inflammatory, and immunomodulatory properties in several clinical trials[22]
HydroxychloroquineRepurposed, RejectedHost factor-targeted. Antimalarial drug that affects endosomal function and blocks autophagosome-lysosome fusion[23]
AzithromycinRepurposedHost factor-targeted. Broad-spectrum antibiotic, blocks autophagosome clearance in human cells[23]
CamostatRepurposedHost factor-targeted. TMPRSS2 inhibitor[24]
NafamostatRepurposedHost factor-targeted. TMPRSS2 inhibitor[25]
FavipiravirRepurposed, Approved in Russia, JapanHost factor-targeted. RNA-dependent/RNA polymerase-targeted[26]
Table 2. Classification of identified ACEi, CCB compounds, and utilized physicochemical and pharmacokinetic (PK) parameters values obtained from literature, default value approximation, or machine learning (ML) optimization.
Table 2. Classification of identified ACEi, CCB compounds, and utilized physicochemical and pharmacokinetic (PK) parameters values obtained from literature, default value approximation, or machine learning (ML) optimization.
Drug NameClasskaLogPpKaFu,pB:PClearance (L/h/kg)KplungKpgut
LisinoprilACEi0.17 *−1.115 [49]3.17 (acid), 10.21 (base) [49]0.99 [49]0.71 *0.072 [49]0.570.50
CaptoprilACEiN/A0.34 [50]4.01 (acid),
−1.2 [49]
0.73 [49]0.45 *0.72 [49]0.15 *0.15 *
SpiraprilACEi0.53 [51]0 *3.62 (acid), 5.2 (base) [49]0.314678 *0.74 **0.43 [51]0.21 *0.16 *
LacidipineCCB1.7843 *5.51 [52]19.47 (acid),
−6.4 (base) [49]
0.05 [53]0.70 *1.23 [49]11.72 *11.72 *
VerapamilCCBN/A3.795 [54]9.68 (base) [49]0.064 [55]0.88 [36,54]0.84 [49,54,55]3.69 *3.69 *
*—optimized values, **—default values.
Table 3. Background on datasets used for systemic plasma-venous compartment disposition simulation and optimization of missing PK parameters.
Table 3. Background on datasets used for systemic plasma-venous compartment disposition simulation and optimization of missing PK parameters.
Drug NameFormulationExperimental SetupReference
Lisinopril20 mg, oral dose20 mg of Lisinopril was given orally for 10 consecutive days. 8 subjects in the study.[56]
Captopril2.78 mg, 5.67 mg, 11.4 mg, IV dose1 mL of intravenous injection at three different dosage levels was administered to 7 healthy subjects.[57]
Spirapril25 mg, oral dose25 mg spirapril p.o. prepared by dissolving 25 mg of lyophilized spirapril in 50 mL tap water was given to the subjects. 16 subjects.[51]
Lacidipine2 mg, 4 mg, oral doseSingle dose of 2 mg and 4 mg of Lacidipine was administered. The study has a total of 24 subjects (12 male, 12 female)[58]
Verapamil50 mg, IV dose5 subjects received 5 mg verapamil dissolved in 30 mL of saline infused over 5 min.[55]
Table 4. Comparison of plasma venous PK metrics between model simulation outputs and experimentally-derived measurements.
Table 4. Comparison of plasma venous PK metrics between model simulation outputs and experimentally-derived measurements.
CompoundsLisinoprilCaptoprilSpiraprilLacidipineVerapamil
Output metricsROAOralIVIVIVOralOralOralIV
Dose, mg202.785.6711.425245
AUC0-t, µg·h/LObservedN/AN/AN/AN/AN/A3.667.66N/A
Calculated752.0042.9793.93215.69991.833.126.80703.03
Predicted823.5049.1692.55212.12977.603.278.80494.26
AAFE1.101.141.011.021.011.051.291.42
AFE1.101.140.990.980.991.051.290.70
Cmax, µg/LObservedN/AN/AN/AN/A430.001.243.09N/A
Calculated57.40104.64234.55454.71378.001.172.871176.82
Predicted53.93105.81183.46497.25196.171.002.011696.96
AAFE1.061.011.281.091.931.161.431.44
AFE0.941.010.781.090.520.860.701.44
Tmax, hObservedN/AN/AN/AN/A0.901.131.25N/A
Calculated6.040.150.190.131.001.231.050.09
Predicted6.030.150.190.131.751.051.050.09
AAFE1.001.001.001.001.751.171.001.00
AFE1.001.001.001.001.750.861.001.00
StatisticsChi-squared2803.13 *22.5623.0126.852.70107.9 *182.15 *11.57 *
p-values>0.50>0.50>0.50>0.50<0.001>0.50>0.500.36
Note: outputs marked as “observed” were extracted directly from the source manuscripts. “Calculated” corresponds to recalculation of the output values using internal non-compartmental methods. Predicted Cmax values correspond to the maximum sampling concentration within the time range of observed timepoints. For observed datapoints with poorly visible error bars, an effective error was calculated using the variability from other timepoints. Chi-squared values and associated p-values were calculated assuming a standard deviation of 0.1× in observed values for manuscript without reported measurement error. * differences considered statistically significant.
Table 5. PK metrics calculated from simulations of drug distribution to ACE-expressing organs and tissues.
Table 5. PK metrics calculated from simulations of drug distribution to ACE-expressing organs and tissues.
MetricsTissueLisinopril,
20 mg, oral
Captopril,
2.78 mg, IV
Spirapril,
25 mg, oral
Lacidipine,
2 mg, oral
Verapamil,
50 mg, IV
AUC0–t,
µg·h/L
Lung341.2452.362315.6413.86344.12
Gut354.03196.742310.3113.86354.74
Nasal tissue196.911.48207.7431.63360.48
Nasal epithelium196.341.48207.5831.63360.48
Cmax, µg/LLung34.0715.267514.65140.2465.26
Gut33.68103.163213.0092.7770.37
Nasal tissue19.670.43657.99256.3668.36
Nasal epithelium19.670.43478.23128.7868.35
Tmax, hLung6.651.130.080.021.84
Gut5.680.570.130.031.30
Nasal tissue6.661.130.090.021.85
Nasal epithelium6.691.170.110.051.89
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chakravarty, K.; Antontsev, V.G.; Khotimchenko, M.; Gupta, N.; Jagarapu, A.; Bundey, Y.; Hou, H.; Maharao, N.; Varshney, J. Accelerated Repurposing and Drug Development of Pulmonary Hypertension Therapies for COVID-19 Treatment Using an AI-Integrated Biosimulation Platform. Molecules 2021, 26, 1912. https://doi.org/10.3390/molecules26071912

AMA Style

Chakravarty K, Antontsev VG, Khotimchenko M, Gupta N, Jagarapu A, Bundey Y, Hou H, Maharao N, Varshney J. Accelerated Repurposing and Drug Development of Pulmonary Hypertension Therapies for COVID-19 Treatment Using an AI-Integrated Biosimulation Platform. Molecules. 2021; 26(7):1912. https://doi.org/10.3390/molecules26071912

Chicago/Turabian Style

Chakravarty, Kaushik, Victor G. Antontsev, Maksim Khotimchenko, Nilesh Gupta, Aditya Jagarapu, Yogesh Bundey, Hypatia Hou, Neha Maharao, and Jyotika Varshney. 2021. "Accelerated Repurposing and Drug Development of Pulmonary Hypertension Therapies for COVID-19 Treatment Using an AI-Integrated Biosimulation Platform" Molecules 26, no. 7: 1912. https://doi.org/10.3390/molecules26071912

Article Metrics

Back to TopTop