Skip to main content
Advertisement
  • Loading metrics

Use of Angiotensin-Converting Enzyme Inhibitors and Angiotensin II Receptor Blockers During the COVID-19 Pandemic: A Modeling Analysis

  • Mehrshad Sadria,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada

  • Anita T. Layton

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Validation, Visualization, Methodology, Writing – original draft, Writing – review & editing

    anita.layton@uwaterloo.ca

    Affiliations Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada, Department of Biology, Cheriton School of Computer Science, and School of Pharmacy, University of Waterloo, Waterloo, Ontario, Canada

Abstract

Angiotensin-converting enzyme inhibitors (ACEi) and angiotensin II receptor blockers (ARB) are frequently prescribed for a range of diseases including hypertension, proteinuric chronic kidney disease, and heart failure. There is evidence indicating that these drugs upregulate ACE2, a key component of the renin-angiotensin system (RAS) and is found on the cells of a number of tissues, including the epithelial cells in the lungs. While ACE2 has a beneficial role in many diseases such as hypertension, diabetes, and cardiovascular disease, it also serves as a receptor for both SARS-CoV and SARS-CoV-2 via binding with the spike protein of the virus, thereby allowing it entry into host cells. Thus, it has been suggested that these therapies can theoretically increase the risk of SARS- CoV-2 infection and cause more severe COVID-19. Given the success of ACEi and ARBs in cardiovascular diseases, we seek to gain insights into the implications of these medications in the pathogenesis of COVID-19. To that end, we have developed a mathematical model that represents the RAS, binding of ACE2 with SARS-CoV-2 and the subsequent cell entry, and the host’s acute inflammatory response. The model can simulate different levels of SARS-CoV-2 exposure, and represent the effect of commonly prescribed anti-hypertensive medications, ACEi and ARB, and predict tissue damage. Model simulations indicate that whether the extent of tissue damage may be exacerbated by ACEi or ARB treatment depends on a number of factors, including the level of existing inflammation, dosage, and the effect of the drugs on ACE2 protein abundance. The findings of this study can serve as the first step in the development of appropriate and more comprehensive guidelines for the prescription of ACEi and ARB in the current and future coronavirus pandemics.

Author summary

As we brace for the devastating impact of the COVID-19 pandemic, we must tackle a controversy on how to best minimize the risk of lethal disease among the most vulnerable. Preliminary epidemiological data show an exponential increase in disease severity and mortality among patients with cardiovascular disease and diabetes. The coronavirus enters host cells by binding to a specific enzyme “ACE2” on the cell membrane. ACE2 abundance is increased in patients with cardiovascular disease and diabetes treated with two classes of drugs: angiotensin-converting enzyme inhibitors (ACEi) and angiotensin II receptor blockers (ARB). This has generated controversy regarding the approach for patients taking these drugs during the pandemic, with some advocating for discontinuing these medications, while expert opinions recommended against discontinuation, given the lack of strong evidence. Given the success of ACEi and ARBs in cardiovascular diseases, we aim to help patients and physicians weigh the overall pros and cons. To achieve that goal, we have developed a mathematical model of the invasion of the coronavirus and the host’s immune response. Model simulations indicate how much tissue damage COVID-19 induces in a patient undergoing ACEi or ARB treatment depends on a number of factors, including the level of any existing chronic inflammation, dosage, and certain drugs effects.

Introduction

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has wreaked havoc all over the world. SARS-CoV-2 has a higher transmission rate than the SARS-CoV from 2003, and causes the coronavirus disease 2019 (COVID-19). In some patients COVID-19 causes acute respiratory distress syndrome, which has high morbidity and mortality. There appears to be some indications that COVID-19 tends to be more severe in patients with hypertension and diabetes [1], although it is difficult to assess to what extent that preliminary conclusion can be attributable to bias in age, sex, comorbidities, and existing medication. Nonetheless, concern has emerged that some anti-hypertensive treatments, specifically renin-angiotensin system (RAS) inhibitors, may increase the risk of SARS-CoV-2 infection and lead to more severe COVID-19 owing to the RAS-mediated cell entry mechanism of the virus [2].

Both SARS-CoV and SARS-CoV-2 gain entry into host cells via the binding of its spike protein with the angiotensin converting enzyme 2 (ACE2). ACE2 is a key component of the RAS and is found on the cells of a number of tissues, including the type 2 alveolar epithelial cells in the lungs [3]. The actions of ACE2, and its homologue ACE, are described below. Because ACE2 opposes the actions of the vasoactive angiotensin II (Ang II), ACE2 has a beneficial role in many diseases such as hypertension, diabetes, and cardiovascular disease in which its expression is decreased. Following binding of its spike protein to ACE2 and proteolytic cleavage of ACE2 by transmembrane serine protease 2 (TMPRSS2), the virus enters the cell and replicates.

The RAS consists of an enzymatic cascade that leads to the production of Ang II [4]. The cascade starts with angiotensinogen, which is cleaved by renin, followed by angiotensin converting enzyme (ACE) and neutral endopeptidase (NEP), to yield different forms of angiotensin (see Fig 1) [46]. The first reaction in the RAS is catalyzed by renin, which is often considered the driving force in the cascade. The major products of the RAS include angiotensin (1–7) (Ang (1–7)) and Ang II, which bind to receptors and exert their effects on the brain, heart, kidney, vasculature, and immune system [4]. Ang II binds to two different types of receptors: Ang II type 1 receptor (AT1R) and type 2 receptor (AT2R). When bound to AT1R, Ang II stimulates renal vasoconstriction, promotes sodium reabsorption, inflammation, and fibrosis [5]. In contrast, when bound to AT2R, Ang II induces vasodilation and natriuresis. When bound to the angiotensin type 7 receptor (AT7R), also known as the MAS receptor, Ang (1-7) induces vasodilation and natriuresis, elevates nitric oxide bioavailability, and attenuates inflammation and fibrosis [5, 6]. The conversion of Ang II to Ang (1-7) is promoted by ACE2. See Fig 1.

thumbnail
Fig 1. Schematic model diagram.

Orange boxes, inflammatory response; blue boxes, angiotensin converting enzyme 2 (ACE2) dynamics; purple boxes, renin-angiotensin system; white boxes, model parameters. Arrows represent conversion; square end, activation; flat end, inhibition.

https://doi.org/10.1371/journal.pcbi.1008235.g001

Due to the vasoactive and renal effects of AT1R-bound Ang II and AT2R-bound Ang II, the RAS is a target for blood pressure control. Indeed, some of the RAS inhibitors, such as ACE inhibitors (ACEi) and angiotensin II receptor blockers (ARBs), are prescribed to millions of patients worldwide. Because ACEi and ARBs are known to upregulate ACE2 mRNA level [7], it has been suggested that these therapies can theoretically increase the risk of SARS- CoV-2 infection and cause more severe COVID-19 owing to the role of ACE2 as the viral binding site [2]. However, given the success of ACEi and ARBs in cardiovascular diseases, the decision to discontinue their treatment should not be made lightly. Hence, the goal of this study is to develop a mathematical model to gain insights into the effect of these drugs on the pathogenesis of COVID-19.

Results

Exposure level

We first consider the effect of differing levels of exposure on disease severity, i.e., quickly passing by an asymptomatic patient versus working alongside a symptomatic patient with poor personal hygiene. To simulate low versus high level of exposure, we initialize plasma SARS-CoV-2 to the corresponding level (1 versus 10 pathogen units). Key simulations are summarized in Fig 2, panels A1-A3 and B1-B3.

thumbnail
Fig 2. Model response to low level initial exposure to covid-19 (panels A1-A3) versus high level (B1-B3), and effect of a stronger immune response with high exposure level (C1-C3).

Panels A1-C1, tissue damage (D). Recovery is possible at low exposure level (A1). High exposure level yields irrecoverable tissue damage (possibly death; B1), but recovery is possible if immune response is stronger (C1). Panels A2-C2, cellular SARS-CoV-2 level. Viral clearance is possible at low exposure level, or high level with a sufficiently strong immune response (A2 and C2), but clearance fails at high exposure level (B2). Panels A3-C3, inflammatory agents (e.g. activated phagocytes). Note different ordinate range for the B panels (nondimensional units).

https://doi.org/10.1371/journal.pcbi.1008235.g002

Plasma SARS-CoV-2 combines with membrane ACE2 to gain cell entry, increasing cellular SARAS-CoV-2 population (Fig 2A2 and 2B2) and causing tissue damage (Fig 2A1 and 2B1). When internalized SARS-CoV-2 exceeds a threshold (Pc), an immune response is triggered (Fig 2A3 and 2B3). At low exposure level, the immune response can suppress the viral population to a sufficiently low level for recovery (Fig 2A1). At sufficiently high exposure, however, tissue damage becomes severe and tissue recovery is not achieved (Fig 2B1). For a given set of parameters, the model prediction is interpreted as “tissue recovery” if the long-term tissue damage (D) approaches its initial small value or as “sustained tissue damage,” when a sufficiently large steady-state D value is predicted.

Strength of immune response

We conduct additional sets of simulation in which we either (i) assume that an effective immune response is triggered earlier (i.e., at half the baseline internalized SARS-CoV-2 threshold Pc) or (ii) the immune response, once triggered, is more effective in reducing the viral load (specifically, we double kpm and kpn in Eq 10 in methods; kpm and kpn characterize the strength of the local immune response and the inhibitory action of the activated phagocytes). In both cases, tissue recovery is predicted. Key results for the case with the stronger immune response is shown in Fig 2C1–2C3; results for the earlier trigger case are similar. These results highlight the importance of individual difference in immune response in determining disease outcome.

Hypertension, ACEi, and ARBs

In the next set of simulations, we seek to understand the complex story around hypertension and some of its treatments in the context of COVID-19. Hypertension is a complex and multifactorial disease. We focus on the observation that low-grade chronic inflammation is often involved in the development and pathophysiology of hypertension [8]. Inflammatory status can be altered by a number of factors. For simplicity, we simulate an inflammatory state by increasing the generate rate of activated neutrophils (snr in Eq 11 in methods).

Existing inflammation may increase the severity of a COVID-19 infection.

We simulate the immune response of patients with an existing mild inflammation (snr × 2) and with a more severe inflammation (snr × 5). A low initial exposure level is assumed, so that the control (normotensive, without an existing inflammation) achieves tissue recovery (see Fig 2 and Table 1, “Normotensive” without drug treatment). The mild inflammation case attains tissue recovery as well (Table 1, “HTN mild inflame,” without drug treatment). In contrast, the more severe inflammation case fails to reach tissue recovery, even with only low initial exposure. This result indicates that a sufficiently severe existing inflammation may increase the likelihood of serious complications following a COVID-19 infection.

thumbnail
Table 1. Simulated COVID-19 outcome for normotensive patients, hypertensive (HTN) patients with an existing mild inflammation, and HTN patients with a more severe inflammation.

ACEi, angiotensin-converting enzyme inhibitor ARB, angiotensin II receptor blocker. In the “ACEi” column, “Yes” corresponds to full effect of ACEi, i.e., reduced conversion rate from Ang I to Ang II and upregulation of ACE2; “w/o ACE2 ↑” and “ACE2 ↑ only 50%” represent only the reduced Ang I-to-Ang II conversion rate, with ACE2 level kept at baseline or 2.5-fold increase, respectively; “Half dose” represents half the ACEi effect, in both Ang I-to-Ang II conversion and ACE2 upregulation. Labels under the “ARB” column are defined analogously.

https://doi.org/10.1371/journal.pcbi.1008235.t001

ACEi, ARBs and COVID-19.

As previously noted, concerns have been expressed regarding the possibility that ACEi may aggravate the severity of a SARS-CoV-2 infection, due to the upregulation of ACE2 [7], which facilitates SARS-CoV-2 cell entry, following the treatment [9]. The same concern may apply to ARBs, which have also been reported to upregulate ACE2 mRNA [7]. But given their major cardiovascular benefits, a premature rejection of these therapies might do more harm than good. Thus, we seek to glean insight into the complex story around ACEi, ARBs, and COVID-19. Specifically, we conduct simulations to determine whether the administration of ACEi and ARBs in different patient subpopulations changes tissue recovery outcomes. Our findings are summarized in Table 1 and Fig 3.

thumbnail
Fig 3.

Schematic diagram summarizing predicted effects, and the underlying mechanisms, of ACEi (a and c) and ARB (b and d) under conditions of no or mild inflammation (a and b) or more severe inflammation (c and d).

https://doi.org/10.1371/journal.pcbi.1008235.g003

Following our previous work [10], ACE inhibition is modelled by reducing the conversion of Ang I to Ang II by 78%. Thus, ACE inhibition lowers Ang II, as well as AT1R-bound Ang II (pro-inflammatory) and AT2R-bound Ang II (anti-inflammatory) downstream. ACE inhibition is known to upregulate ACE2 mRNA level [7]. The higher abundance of ACE2 would increase (i) the rate of conversion of Ang II to Ang (1-7) which is anti-inflammatory, and (ii) the rate of SARS-CoV-2 cell entry. Both rates are increased by 5 folds based on the change in ACE2 mRNA level reported in Ref. [7], with the assumption that these rates are directly proportional to ACE2 mRNA level. Similarly, ARB is modelled by reducing the conversion of Ang II to AT1R-bound Ang II by 67% [10]. Taken in isolation, ARB does not affect Ang II, but lowers the amount of AT1R-bound Ang II (pro-inflammatory) and increases the amount of AT2R-bound Ang II (anti-inflammatory). ACE2 abundance is assumed to increase 3 folds by ARB, again based on the change in ACE2 mRNA level [7].

Normotension.

Consider a normotensive individual, without a significant existing inflammation. Following a low initial exposure to SARS-CoV-2, the model predicts tissue recovery in the absence of ACEi or ARBs. However, when ongoing ACEi treatment is simulated, tissue recovery is no longer possible following the same SARS-CoV-2 exposure. Why does ACEi impede tissue recovery? To understand the underlying mechanism, we eliminate the effect of ACEi on ACE2 abundance. That is, we only reduce the conversion of Ang I to Ang II by the same degree, but keep ACE2 abundance at baseline (specifically, we keep the rate of SARS-CoV-2 cell entry at baseline values). With this configuration, the mode predicts tissue recovery following SARS-CoV-2 exposure. This result suggests that, in this case, the upregulated ACE2 abundance, and the resulting higher SARS-CoV-2 cell entry rate, may be responsible, at least in part, for the potential negative impact of ACEi on tissue recovery in COVID-19. These results and the underlying mechanisms are summarized in Table 1 and Fig 3(A).

Consider the same normotensive individual, now under ARB treatment. The model again predicts sustained tissue damage. This time tissue recovery is not possible even with ACE2 abundance kept at baseline. By inhibiting the binding of Ang II to AT1R, ARB reduces the pro-inflammatory AT1R-bound Ang II (smaller N*) while increasing the anti-inflammatory AT2R-bound Ang II (larger CA). While these changes attenuate tissue damages the (the kdnfs(f(N*)) term in Eq 12, methods), they also impede the defense against pathogens (the kpnf(N*)P term in Eq 10). With the invasion of a new virus, attenuation of the inflammatory defense can be costly. (See Table 1 and Fig 3(A)). To confirm this line of reasoning, we conduct the same simulation, with the pro-inflammatory effect of AT1R-bound Ang II on N* enhanced (by increasing kAT1R 10 folds) and the anti-inflammatory effect of AT2R-bound Ang II on CA reduced (by decreasing kAT2R 10 folds). With these changes, the model predicts eventual recovery from COVID-19 induced tissue damage even under ARB treatment.

Hypertension with mild inflammation.

Next we consider a hypertensive patient, with an existing mild chronic inflammation. Following a low initial exposure to SARS-CoV-2, the model predicts tissue recovery. However, that tissue recovery is lost when combined with ACEi treatment. Again, why does ACEi impede tissue recovery? As in the normotensive case, the explanation lies with the ACEi-induced upregulation of ACE2. It is noteworthy that we assume 5-fold increase in ACE2, based on a study in rat cardiac tissues [7]. It is not clear to what extent the same increase should be found in human lung tissues. If the ACEi-induced ACE2 increase is lowered to 50%, then the model predicts tissue recovery. See Table 1 and Fig 3(a). This result suggests that the concern about the use of ACEi in patients with mild inflammation may depend on the drug’s effect on ACE2 abundance, which has yet to be sufficiently well characterized in humans.

The model predicts that with ARB treatment, with or without ACE2 upregulation, tissue recovery is not possible (Table 1 and Fig 3(b)). The explanation is similar to the normotensive case above. It is noteworthy that the effects of ACEi and ARBs on the model’s pro- and anti-inflammatory signals are not well quantified. Furthermore, our simulations indicate that when the effect of ACEi or ARB is halved (conversion of Ang I to Ang II reduced by 39% in ACEi and binding of Ang II to AT1R reduced by 33.5%), tissue recovery is attained. The attenuation of the effect of ACEi or ARB can be interpreted as a lower dosage. Thus, this result highlights the importance of dosage in determining the effects of ACEi and ARBs in the content of COVID-19.

Hypertension with more severe inflammation.

Finally, we consider a hypertensive patient with an existing more severe chronic inflammation. Unlike the previous cases, without ACEi or ARB treatment, sustained tissue damage is predicted following a low initial exposure to SARS-CoV-2. Interestingly, the model predicts that with the ACEi or ARB treatment, tissue recovery becomes possible. See Table 1 and (Fig 3(c) and 3(d)). It is important to note that this result is qualitatively different from the normotensive or mild inflammation cases. With an already elevated level of inflammation and tissue damage, the anti-inflammatory benefits of these drugs outweigh other factors, such as the increased ACE2-mediated cell entry by SARS-CoV-2. Thus, ACEi and ARBs may facilitate recovery from COVID-19 induced tissue damage in patients with an existing more severe chronic inflammation.

Sensitivity analysis

Many of the model parameters have not been characterized. To assess the sensitivity of model prediction on variations in those parameters, we conduct a sensitivity analysis (Table 2) in all model parameters listed in Tables 3 and 4, except for the half-life of ACE2 which has been measured experimentally [11]. For parameters whose values have not been directly measured but were inferred from experimental findings (kAT1R, kACE2p, kAng(1-7), and kAT2R), we vary their baseline values by ±10%; for other parameters which are largely uncharacterized, we multiply and divide their baseline values by 2. For each of these variations, we conduct simulation of low-level SARS-CoV-2 exposure, and determine whether tissue recovery is attained, and if so, compare the peak viral load (P) and tissue damage (D) to baseline values.

thumbnail
Table 2. Parameter sensitivity study.

Percentage changes in peak viral load (P) and tissue damage (D) are reported only for cases that achieve tissue recovery.

https://doi.org/10.1371/journal.pcbi.1008235.t002

thumbnail
Table 3. Predicted variables and parameters for ACE2 dynamics.

Baseline values are given for model parameters.

https://doi.org/10.1371/journal.pcbi.1008235.t003

thumbnail
Table 4. Predicted variables and parameters for immune response.

Baseline values are given for model parameters.

https://doi.org/10.1371/journal.pcbi.1008235.t004

Results of the sensitivity analysis are shown in Table 2. These results indicate that model results are largely insensitive to kburst, the rate constant for release of cellular virus into plasma, but for all other parameters, a sufficiently large deviation from baseline values in at least one direction would render tissue recovery impossible. In particular, this analysis indicates that model results are particularly sensitive to the pro-inflammatory effect of AT1R-bound Ang II (kAT1R) and to the anti-inflammatory effect of ACE2 (kACE2p). Interestingly, a sufficiently large deviation in these parameters from their baseline values in either direction would result in sustained tissue damage. This result highlights the role of inflammation as a double-edged sword. A weakened in the inflammatory response simultaneously attenuates tissue damage and impedes the system’s defense against pathogens. With the invasion of a new virus, attenuation of the inflammatory defense may yield sustained tissue damage.

Discussion

The principal goal of this study is to gain insights into the implications of hypertension and some of its medications in the pathogenesis of COVID-19. To that end, we have developed a mathematical model that represents the RAS, the binding of a key RAS component, ACE2, with SARS-CoV-2 and the subsequent cell entry, and the host’s inflammatory response. The model can simulate different levels of SARS-CoV-2 exposure, and the effect of commonly prescribed anti-hypertensive medications, ACEi and ARBs, and predict tissue damage and recovery.

What are the factors that determine the severity of COVID-19? Our simulation results indicate that, not surprisingly, the level of exposure to SARS-CoV-2 is a key factor, as is the robustness of one’s immune system (see Fig 2). Aspects of the immune system is known to be impaired in the aged population [12], which may explain, in part, the higher COVID-19 mortality in that subpopulation. Pre-existing conditions such as hypertension are also known play a role. Hypertension is a highly complex disease that is difficult to simulate. Instead, we focus on the low-grade inflammation that is frequently associated with hypertension. Not surprisingly, our simulations indicate that a pre-existing chronic inflammation may increase the difficulty of recovering from tissue damaged sustained in COVID-19 (Table 1 and Fig 3). Chronic inflammation is also found in diabetes, cardiovascular diseases, and allergies. In a recent study, elevated extracellular traps from neutrophils have been observed in COVID-19 patients and correlates with worse outcomes [13].

SARS-CoV-2 gains entry into host cells by first binding with a membrane receptor ACE2, an essential component of the RAS. This mechanism may have implications on the treatment of hypertension. Although neither ACEi nor ARBs directly affects ACE2 activities, studies in rodents have shown that these agents upregulate the expression and activity of ACE2 in cardiac tissue [7]. This observation has prompted concerns about the potential enhanced susceptibility of patients receiving these drugs to SARS-CoV-2 infection and COVID-19 severity. How might ACEi and ARBs impact recovery from tissue damage in COVID-19? The answer is, it depends, as summarized in Table 1 and Fig 3:

  • For individuals with only mild chronic inflammation or none at all, the administration of ACEi may indeed exacerbate tissue damage. The elevated viral entry outweighs the pro-inflammatory benefits of the drug. The caveat is that this balance depends on our baseline assumption of ACEI-mediated ACE2 upregulation: a 5-fold increase, based on a study in the rat cardiac muscle [7]. If ACEi increases ACE2 expression in human pulmonary tissue to a lesser extent or not at all, then ACEi may not impede tissue recovery.
  • For the same cohort (mild or no chronic inflammation), administration of ARBs inappropriately disrupts the body’s inflammatory and immune response, and is predicted to worsen tissue damage. A caveat is that how much ARBs affect the pro- and anti-inflammatory signals is not well quantified.
  • For individuals with more severe inflammation, ACEi and ARBs may facilitate tissue recovery. This is attributable to the anti-inflammatory benefits of these drugs.

ACEi and ARBs are prescribed for a range of diseases including hypertension, proteinuric chronic kidney disease, and heart failure. Although other drug classes may be available for the management of hypertension, patients with proteinuric chronic kidney disease or heart failure with impaired left ventricular function have limited alternatives. Given the renal and cardiovascular protective effects of these drugs, discontinuing their use will likely have significant negative health impacts. Indeed, while some have advocated for discontinuing these medications during the pandemic, expert opinions recommended against discontinuation, given the lack of strong evidence, or evidence indicating the lack of association of ACEi or ARB treatment with increased risk for acquiring COVID-19 or for poorer outcomes [14]. Therefore, a decision regardless the prescription and usage of ACEi and ARB should not be made lightly. It is our hope that the findings of this study serve as a step in the development of appropriate and more comprehensive guidelines for the prescription of these drugs in the current and future coronavirus pandemics.

This study presents the first mechanistic model that simulates the interactions between the two of the primary players in SARS-CoV-2 infection: the RAS and the acute inflammatory response. The model allows us to weigh the anti-inflammatory benefit of ACE2 against its role in facilitating viral entry into host cells, which is an important consideration in COVID-19 pathogenesis. Despite its contributions, major limitations of the model must be acknowledged. Compared to our published RAS model [10], the present model contains a simplified representation of the RAS that assumes a constant plasma renin activity (PRA), which implies that the feedback from AT1R-bound Ang II to PRA is not represented. Thus, the model does not account for the feedback effects in which a drop in AT1R-bound Ang II, following the administration ACEi and ARBs, increases renin and Ang I, and subsequently increases Ang (1-7).

Additionally, the model only represents the acute inflammatory response; the adaptive or chronic immune responses are not included. During the incubation and non-severe stages, a specific adaptive immune response is required to eliminate the virus and to preclude disease progression to severe stages. Therefore, modeling the acute immune response in COVID-19 (e.g., Ref. [15]) may accelerate the development of strategies for boosting immune responses (anti-sera or pegylated IFNα) at this stage. Whether SARS-CoV-2 may induce a chronic immune response, as seen in HIV infection, has yet to be determined.

A major limitation of the model is the lack of available data, both for estimating model parameters and for validating model predictions. This is unfortunately often the case with diseases without an animal model. The sensitivity analysis (Table 2) suggests that a better characterization of the pro-inflammatory effect of AT1R-bound Ang II and of the anti-inflammatory effect of ACE2 would improve the model’s prediction on the extent of tissue damage resulting from COVID-19. The simulations of ACEi and ARB indicate that a better understanding of the effect of these drugs on ACE2 abundance and viral entry would yield a clearer indication of these drugs on the pathogenesis of COVID-19. Carefully measured time-courses of viral titer and the concentrations of various immune agents would be useful in the development of a model with a more comprehensive immune response and in validating its predictions.

An important message of this study is that there is no one-size-fits-all approach in drug prescription. Whether ACEi and ARBs might improve or impede tissue recovery in COVID-19 depends on the conditions of the patient. It is noteworthy that significant sex differences have been reported in the RAS, and the current RAS model parameters are taken for male [10]. We conducted simulates using parameters fitted for female, and model predictions are qualitatively similar. Many other individual differences exist in hypertension besides sex [1620], including race [21], comorbidities, as well as other medications, including metformin and SGLT2 inhibitors [22, 23] for diabetes, and these factors should be taken into account.

Methods

Mathematical model

The model represents the RAS and the key role of one of its components, ACE2, in the host cell entry of SARS-CoV-2, as well as the inflammatory response that it triggers. For a susceptible individual, following an initial viral exposure, the model predicts the population dynamics of SARS-CoV-2, inflammatory agents, anti-inflammatory mediators, tissue damage, ACE2, and other key RAS components. The long-term model solution predicts the severity of tissue damage.

RAS components.

Recall that ACE2, which serves as the entry receptor for SARS-CoV and SARS-CoV-2, promotes the conversion of Ang II to Ang (1-7). The RAS model is based on our published model [10] and includes key parts of the RAS cascade; see Fig 1 (purple boxes). Ang I is converted into other forms through ACE and NEP; the respective reaction rate constants are denoted cACE and cNEP. The rate of change of [Ang I] is given by (1) where hAng I denotes the half-life of Ang I.

Ang II is converted from Ang I through ACE and then converted into Ang (1–7) through ACE2 with the rate constant cACE2. Ang II also binds to the AT1R and AT2R with rate constants cAT1R and cAT2R, respectively. We assume that Ang II has a half-life of hAng II. With this notation the rate of change of [Ang II] is given by (2)

Ang II binds to AT1R with rate constant cAT1R to produce AT1R- bound Ang II. AT1R-bound Ang II decays with a half-life of hAT1R, giving the rate of change of [AT1R-bound Ang II] as (3)

Ang II binds to the AT2R with rate constant cAT2R to produce AT2R-bound Ang II. AT2R-bound Ang II decays with a half-life of hAT2R, giving the rate of change of [AT2R-bound Ang II] as (4)

Ang (1-7) is converted by NEP from Ang I and by ACE2 from Ang II, and it decays with a half life of hAng(1-7). Thus, the rate of change of [Ang (1-7)] is (5)

ACE2 dynamics.

ACE2 is a type 1 transmembrane protein (N-terminus outside, C-terminus intracellular), predominantly localized on endothelial cells where its catalytic site is exposed (the so-called “ectoenzyme”) to circulating vasoactive peptides [24]. The expression level of ACE2 can be modulated through its cleavage from the cell membrane. This cleavage or shedding releases the catalytically active ectodomain and is mediated by a disintegrin and metalloprotease (ADAM 17) [25]. The ACE2 that has shedded into plasma may bind to SARS-Co-V-2 and essentially deactivate it.

Membrane-bound ACE2 (denoted ACE2m) is generated at a rate of SACE2m. It is cleaved by ADAM 17 and sheds with rate constant cADAM17 and internalized by AT1R-bound Ang II with rate constant cAT1R. Additionally, ACE2m binds to SARS-CoV-2 with rate constant cCov19. Thus, the rate of change of ACE2m is given by (6) where hACE2 denotes the half-life of ACE2.

The shedding of ACE2m generates plasma ACE2 (denoted ACE2p).

(7)

The AT1R-bound Ang II-mediated internalization of ACE2m and generates ACE2i, as does the entry of the ACE2-bound SARS-CoV-2: (8)

We also track plasma SARS-CoV-2. Its population (Pplasma) is determined by the release of internalized SARS-CoV-2 (P in Eq 10) into plasma, characterized by rate constant kburst, and the binding of plasma SARS-CoV-2 with ACE2m and ACEp, the latter characterized by rate constant kACE2-P: (9)

Inflammatory response components.

Pathogens induce an innate immune response, through which the organism aims to eliminate threats to survival and to promote tissue repair and healing. We have extended the model in [26] to represent its interactions with the RAS and the internalization of SARS-CoV-2; see Fig 1 (orange boxes). It represents the interactions among internalized (i.e., cellular) SARS-CoV-2 (denoted P), inflammatory agents (N*), tissue damage (D), and anti-inflammatory mediators (CA). The internalized SARS-CoV-2 population depends on its cell entry (first time in Eq 10), with rate constant cCov19 and scaling constant kCov19, a logistic term with growth rate kpg and carrying capacity P (second term), local immune response (third term, see derivation in [26]), inhibitory action of the activated phagocytes (fourth term), and the release of internalized viral content into plasma following cell eruption (five term): (10) where describes the inhibitory action of the anti-inflammatory agents. To simulate the delay of immune response, kpm and kpn are initialized to 10% of their baseline value while P < Pc, where the critical Pc is taken to be 2 pathogen units.

Inflammatory agents such as phagocytes are activated in response to damaged tissue and pathogens, as well as its own positive feedback loop (first term in Eq 11). Additionally, AT1R-bound Ang II promotes inflammation (second term). N* decays with a rate constant μn. (11) where

Tissue damage is enhanced by the activated phagocytes but attenuated by the anti-inflammatory agents; it also decays with a rate constant μd. (12) where

The CA equation contains a source of CA, denoted sc, and a term modeling the production of anti-inflammatory mediator from damage and activated phagocytes, which takes the form and is inhibited by CA. This expression is a Michaelis-Menten-type term, in which kcnd controls the effectiveness of damage, relative to activated phagocytes, in producing CA. The model also represents the effect of anti-inflammatory agents: plasma ACE2, Ang (1-7) and AT2R-bound Ang II, as follows (13) where μc denotes the decay rate constant and

Model parameters

Parameters for the RAS are taken from the male model in [10]. Parameters describing ACE2 dynamics are given in Table 3. SACE2m, cADAM17, and kAT1R are estimated by setting Eqs 69 to steady state, assuming a steady-state distribution of ACE2m : ACE2p : ACE2i of 0.94:0.045:0.015. cCov19 is chosen to fit the viral growth rate reported in [27]. Parameters for the innate immune response are taken from Ref. [26]. Additional parameters are shown in Table 4.

Initial conditions.

RAS variables are initialized to their steady-state values in the corresponding standalone RAS model, i.e., with fixed ACE2 expression. Similarly, the immune response variables are initialized to their steady-state values in a standalone immune response model with no pathogen (i.e., P = 0).

References

  1. 1. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet, 2020. 395(10229):1054–1062. pmid:32171076
  2. 2. South AM, Tomlinson L, Edmonston D, Hiremath S, Sparks MA. Controversies of renin-angiotensin system inhibition during the COVID-19 pandemic. Nat Rev Nephrol. 2020;16(6):305–307. pmid:32246101
  3. 3. Clarke NE, Turner AJ. Angiotensin-converting enzyme 2: the first decade. Int J Hypertens. 2012;2012:307315. pmid:22121476
  4. 4. Sparks MA, Crowley SD, Gurley SB, Mirotsou M, Coffman TM. Classical Renin-Angiotensin system in kidney physiology. Compr Physiol. 2014;4(3):1201–28. pmid:24944035
  5. 5. Chappell MC. Nonclassical renin-angiotensin system and renal function. Compr Physiol. 2012;2(4):2733–52. pmid:23720263
  6. 6. Chappell MC, Marshall AC, Alzayadneh EM, Shaltout HA, Diz DI. Update on the Angiotensin converting enzyme 2-Angiotensin (1-7)-MAS receptor axis: fetal programing, sex differences, and intracellular pathways. Front Endocrinol (Lausanne). 2014;4:201.
  7. 7. Ferrario CM, Jessup J, Chappell MC, Averill DB, Brosnihan KB, Tallant EA, et al. Effect of angiotensin-converting enzyme inhibition and angiotensin II receptor blockers on cardiac angiotensin-converting enzyme 2. Circulation. 2005;111(20):2605–10. pmid:15897343
  8. 8. Savoia C, Schiffrin EL. Inflammation in hypertension. Curr Opin Nephrol Hypertens. 2006;15(2): 152–8. pmid:16481882
  9. 9. Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S.et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell. 2020;181(2):271–280.e8. pmid:32142651
  10. 10. Leete J, Gurley S, Layton A. Modeling Sex Differences in the Renin Angiotensin System and the Efficacy of Antihypertensive Therapies. Computers & Chemical Engineering. 2018;112:253–264.
  11. 11. Wysocki J., et al., Targeting the Degradation of Angiotensin II With Recombinant Angiotensin-Converting Enzyme 2: Prevention of Angiotensin II–Dependent Hypertension. Hypertension, 2010. 55(1): p. 90–98. pmid:19948988
  12. 12. Gomez C.R., Boehmer E.D., and Kovacs E.J., The aging innate immune system. Curr Opin Immunol, 2005. 17(5): p. 457–62. pmid:16084711
  13. 13. Zuo Y, Yalavarthi S, Shi H, Gockman K, Zuo M, Madison JA, et al. Neutrophil extracellular traps (NETs) as markers of disease severity in COVID-19. medRxiv, 2020. https://doi.org/10.1101/2020.04.09.20059626
  14. 14. Reynolds HR, Adhikari S, Pulgarin C, Troxel AB, Iturrate E, Johnson SB. et al., Renin-Angiotensin-Aldosterone System Inhibitors and Risk of Covid-19. N Engl J Med. 2020;382:2441–2448. pmid:32356628
  15. 15. Sadria M, Layton AT. Modeling Within-host SARS-CoV-2 Infection Dynamics and Potential Treatments. Scientific reports, 2020.
  16. 16. Leete J, Layton AT. Sex-specific long-term blood pressure regulation: Modeling and analysis. Comput Biol Med. 2019;104: 139–148. pmid:30472496
  17. 17. Chen Y, Sullivan JC, Edwards A, Layton AT. Sex-specific computational models of the spontaneously hypertensive rat kidneys: factors affecting nitric oxide bioavailability. Am J Physiol Renal Physiol. 2017;313(2):F174–F183. pmid:28356289
  18. 18. Ahmed S, Layton AT. Sex-specific computational models for blood pressure regulation in the rat. Am J Physiol Renal Physiol. 2019.
  19. 19. Hu R, McDonough AA, Layton AT. Functional implications of the sex differences in transporter abundance along the rat nephron: modeling and analysis. Am J Physiol Renal Physiol. 2019;317(6):F1462–F1474. pmid:31566436
  20. 20. Hu R, McDonough AA, Layton AT. Sex-Differences in Solute Transport Along the Nephrons: Effects of Na+ Transport Inhibition. American Journal of Physiology-Renal Physiology. 2020.
  21. 21. Gillum RF. Pathophysiology of hypertension in blacks and whites. A review of the basis of racial blood pressure differences. Hypertension. 1979;1(5):468–75. pmid:541040
  22. 22. Layton AT, Vallon V. SGLT2 Inhibition in a Kidney with Reduced Nephron Number: Modeling and Analysis of Solute Transport and Metabolism. Am J Physiol Renal Physiol. 2018;314(5):F969–F984. pmid:29361669
  23. 23. Layton AT, Vallon V. Cardiovascular benefits of SGLT2 inhibition in diabetes and chronic kidney diseases. Acta Physiol (Oxf). 2018;222(4):e13050.
  24. 24. Warner FJ, Lew RA, Smith AI, Lambert DW, Hooper NM, Turner AJ. Angiotensin-converting enzyme 2 (ACE2), but not ACE, is preferentially localized to the apical surface of polarized kidney cells. J Biol Chem. 2005:280(47):39353–62. pmid:16166094
  25. 25. Lambert W, Yarski M, Warner FJ, Thornhill P, Parkin ET , Smith AI,et al. Tumor Necrosis Factor- Convertase (ADAM17) Mediates Regulated Ectodomain Shedding of the Severe-acute Respiratory Syndrome-Coronavirus (SARS-CoV) Receptor, Angiotensin-converting Enzyme-2 (ACE2). Journal of Biological Chemistry. 2005;280:0113–30119.
  26. 26. Reynolds A, Rubin J, Clermont G, Day J, Vodovotz Y, Ermentrout GB. A reduced mathematical model of the acute inflammatory response: I. Derivation of model and analysis of anti-inflammation. J Theor Biol. 2006;242(1):220–36. pmid:16584750
  27. 27. Pan Y, Zhang D, Yang P, Poon LLM, Wang Q. Viral load of SARS-CoV-2 in clinical samples. Lancet Infect Dis. 2020;20(4):411–412. pmid:32105638
  28. 28. Jiang JS, Lang YD, Chou HC, Shih CM, Wu MY, Chen CM, et al. Activation of the renin-angiotensin system in hyperoxia-induced lung fibrosis in neonatal rats. Neonatology. 2012;101(1):47–54. pmid:21791939
  29. 29. Ferreira AJ, Shenoy V, Yamazato Y, Sriramula S, Francis J, Yuan L, et al. Evidence for Angiotensin-converting Enzyme 2 as a Therapeutic Target for the Prevention of Pulmonary Hypertension. Am J Res Crit Care Med. 2009;179(11): 1048–1054.
  30. 30. Shenoy V., Ferreira AJ, Qi Y, Fraga-Silva RA, Díez-Freire C, Dooies A, et al., The angiotensin-converting enzyme 2/angiogenesis-(1–7)/Mas axis confers cardiopulmonary protection against lung fibrosis and pulmonary hypertension. Am J Res Crit Care Med. 2010;182(8):1065–1072.
  31. 31. Wagenaar GT, Laghmani EH, Fidder M, Sengers RMA, de Visser YP, de Vries L, Rink R, et al. Agonists of MAS oncogene and angiotensin II type 2 receptors attenuate cardiopulmonary disease in rats with neonatal hyperoxia-induced lung injury. Am J Physiol Lung Cell Mol Physiol. 2013;305(5):L341–51. pmid:23812633
  32. 32. Wang H, Bloom O, Zhang M, Vishnubhakat JM, Ombrellino M, Che J, et al. HMG-1 as a late mediator of endotoxin lethality in mice. Science. 1999;285(5425):248–251. pmid:10398600