Next Article in Journal
Binding of CCCTC-Binding Factor (CTCF) to the Minute Virus of Mice Genome Is Important for Proper Processing of Viral P4-Generated Pre-mRNAs
Next Article in Special Issue
Quantification of the Tradeoff between Test Sensitivity and Test Frequency in a COVID-19 Epidemic—A Multi-Scale Modeling Approach
Previous Article in Journal
Special Issue “Hepatitis B Virus Infection: From Diagnostics to Treatments”
Previous Article in Special Issue
Convolutional Neural Network Based Approach to In Silico Non-Anticipating Prediction of Antigenic Distance for Influenza Virus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Molecular Impact of SARS-CoV-2 Infection on the Renin-Angiotensin System

1
Computational Biology and Bioinformatics, Université Libre de Bruxelles, CP 165/61, Roosevelt Ave. 50, 1050 Brussels, Belgium
2
Biosystems Modeling and Control, Université Libre de Bruxelles, CP 165/61, Roosevelt Ave. 50, 1050 Brussels, Belgium
*
Author to whom correspondence should be addressed.
Submission received: 15 October 2020 / Revised: 24 November 2020 / Accepted: 25 November 2020 / Published: 30 November 2020
(This article belongs to the Special Issue Mathematical Modeling of Viral Infection)

Abstract

:
SARS-CoV-2 infection is mediated by the binding of its spike protein to the angiotensin-converting enzyme 2 (ACE2), which plays a pivotal role in the renin-angiotensin system (RAS). The study of RAS dysregulation due to SARS-CoV-2 infection is fundamentally important for a better understanding of the pathogenic mechanisms and risk factors associated with COVID-19 coronavirus disease and to design effective therapeutic strategies. In this context, we developed a mathematical model of RAS based on data regarding protein and peptide concentrations; the model was tested on clinical data from healthy normotensive and hypertensive individuals. We used our model to analyze the impact of SARS-CoV-2 infection on RAS, which we modeled through a downregulation of ACE2 as a function of viral load. We also used it to predict the effect of RAS-targeting drugs, such as RAS-blockers, human recombinant ACE2, and angiotensin 1–7 peptide, on COVID-19 patients; the model predicted an improvement of the clinical outcome for some drugs and a worsening for others. Our model and its predictions constitute a valuable framework for in silico testing of hypotheses about the COVID-19 pathogenic mechanisms and the effect of drugs aiming to restore RAS functionality.

1. Introduction

Since December 2019, the world has been facing a global viral pandemic of the novel severe acute respiratory syndrome coronavirus 2, “SARS-CoV-2”; this pandemic has, to date, caused millions of people to be infected and hundreds of thousands to die [1]. First detected in the city of Wuhan (China) [2,3,4,5], SARS-CoV-2 spread rapidly throughout the world. The coronavirus family, to which SARS-CoV-2 belongs, includes a number of viruses, such as SARS-CoV and MERS-CoV, which have been implicated in serious epidemics that cause acute respiratory distress syndrome (ARDS). There is not yet consensus on the origin of SARS-CoV-2 [6,7,8,9], but the primary hypothesis is that it originated from bat (Rhinolophus affisor) or pangolin (Manis javanica), since the genomes of these two viral species share high sequence identity with SARS-CoV-2.
Coronaviral genomes encode a series of structural proteins, one of which is the spike glycoprotein or S-protein that protrudes from the membrane surface [9]. Similar to the SARS-CoV virus that was identified in 2003, the S-protein of SARS-CoV-2 has been shown to bind to the angiotensin-converting enzyme 2 (ACE2), so that it can be used as an entry receptor to the cell [9,10,11,12,13]. This protein plays a pivotal role in the renin-angiotensin system (RAS) signaling pathway [14] by cleaving angiotensin I and II peptides to generate angiotensin 1–9 and the biologically active peptide angiotensin 1–7, respectively [15,16]. ACE2 is highly expressed in type II alveolar cells of lung, epithelial cells of oral mucosa, colon enterocytes, myocardial cells, and kidney proximal tubule cells. The protective role of ACE2 in severe ARDS is also widely recognized [17,18]. Indeed, it has been shown, both in in vitro and in vivo mouse models, that a loss of ACE2 expression causes increased production of angiotensin II and that this contributes to lung failure [18].
It has already been established years ago that the SARS-CoV spike protein interferes with RAS due to its binding to ACE2 [19], thus causing ACE2 downregulation; this has opened up a number of interesting means of tackling SARS-CoV infection through RAS modulation. Indeed, injection of a soluble form of recombinant human ACE2 (rhACE2, GSK2586881) into mice infected with SARS-CoV appears to have a double role [18]: it slows the viral infection by binding to the S-protein and rescues ACE2 activity, thus causing angiotensin II reduction and protecting lung from severe failure.
rhACE2 has been tested in phase II trials for its ability to ameliorate ARDS [20]. Although rhACE2 treatment is well tolerated by patients and offers a significant reduction in the angiotensin II level, clinical distress severity was not reduced in a recent pilot study [20]. Further studies are needed to understand the biological differences between the responses of animal models and humans.
Since SARS-CoV-2 also targets ACE2 receptors when it infects cells, it is logical to hypothesize that rhACE2 might help reduce the severity of COVID-19 disease [21]. Indeed, it has been shown that rhACE2 inhibits SARS-CoV-2 infection in vitro and that this inhibition depends both on the initial quantity of the virus and on the rhACE2 concentration [22]. Following these exciting results, a clinical trial with exogenous submission of rhACE2 recently started [23]. A number of other clinical trials are also underway that target the dysregulated RAS to restore its functionality [24,25,26,27,28].
Hypertension and cardiovascular disease have been shown to be risk factors in cases of SARS-CoV-2 infection. This brings into question what might be the potential effects on the COVID-19 development of the RAS-targeting drugs that are used to treat hypertension and cardiovascular disease. RAS-targeting drugs fall into three categories: (i) angiotensin-converting enzyme inhibitors (ACE-I), (ii) angiotensin receptor blockers (ARBs), and (iii) direct renin inhibitors (DRIs) (Figure 1). Several recent studies on large patient cohorts [29,30,31] concluded that there is only a weak correlation between treatment with drugs from these categories and any substantial increase in the risk of COVID-19.
Despite these interesting findings, there is not yet a detailed understanding of how SARS-CoV-2 infection leads to a dysregulation of RAS and, in severe cases, to ARDS. It is of fundamental importance that we gain better insights into the perturbed RAS in order to properly elucidate the pathogenic mechanisms and associated risk factors of SARS-CoV-2 infection; this, in turn, will enable novel therapeutic strategies to be designed and tested so that disease progression can be inhibited.

2. Methods

2.1. Modeling the Renin-Angiotensin System

RAS has been widely studied both experimentally [32,33,34] and computationally [35,36,37,38]. It plays a key role in the regulation of a large series of physiological systems including the renal, lung, and cardiovascular systems. Consequently, its dysregulation is related to multiple pathological conditions such as hypertension and ARDS, just to mention a few [39,40,41,42,43].
There are two different types of RAS: the circulating RAS that is localized in the plasma and is involved in the regulation of the cardiovascular system and the tissue-localized systems that act intracellularly or interstitially within different organs in association with the systemic RAS or independently of it. Here, we focus on the local RAS within the pulmonary circulation and model its network of biochemical reactions, as schematically depicted in Figure 1.
When the blood pressure decreases, the juxtaglomerular kidney cells that sense changes in renal perfusion pressure secret an aspartic protease protein called renin (RE, EC 3.4.23.15). The activity of this enzyme, called plasma renin activity ( P R A ), is the common measure used in clinical practice to set up the diagnosis and treatment design of essential hypertension.
The dynamics of the renin concentration can be modeled as:
d [ RE ] d t = β Log 2 h r e [ RE ]
where h r e is renin’s half-life and β its production rate. The latter is not constant, but depends on other elements of RAS, which we will discuss later in the section. The role of renin is to cleave the N-terminus of a protein from the serine protease inhibitor family called angiotensinogen (AGT) to form the decapeptide hormone angiotensin I (AngI). The dynamics of the angiotensinogen can be written as:
d [ AGT ] d t = k a g t c r e [ RE ] Log 2 h a g t [ AGT ]
where the reaction rate c r e relates the renin concentration to its activity, k a g t is AGT’s production rate, and h a g t its half-life.
The AngI peptide is further cleaved by different enzymes:
  • The angiotensin-converting enzyme (ACE, EC3.4.15.1) is a zinc metalloproteinase located mainly in the capillaries of lung and in the endothelial cells. It catalyzes the transformation of AngI into the octapeptide angiotensin II (AngII).
  • Chymase (CHY, EC 3.4.21.39), a serine protease that is mainly localized in blood vessels and heart, also catalyzes the transformation of AngI into AngII.
  • Neprilysin (NEP, EC3.4.24.11), another zinc metalloproteinase that is expressed in a wide variety of tissues, catalyzes the transformation of AngI into the heptapeptide hormone angiotensin-(1-7) (Ang1-7).
The dynamics of AngI can thus be described as:
d [ AngI ] d t = c r e [ RE ] c a c e + c c h y + c n e p [ AngI ] Log 2 h a n g I [ AngI ]
where c a c e , c c h y , and c n e p are the reaction rates associated with the corresponding enzymatic reactions.
The role of AngII in RAS is central since it has a vasoconstriction effect, enhances blood pressure, and triggers inflammatory processes and fibrosis. In lung, the capillary blood vessels are among the sites that have the highest ACE expression and production of AngII. Its dysregulation has frequently been related to a wide series of chronic and acute diseases such as pulmonary fibrosis and ARDS.
AngII’s effects are mediated by two G-protein coupled receptors (GPCR) called angiotensin II type 1 (AT1R) and type 2 (AT2R). In addition, it can be cleaved by different enzymes. For example, ACE2 generates Ang1-7 peptides, and aminopeptidase A (APA, EC 3.4. 11.7) generates other peptides such as angiotensin III (AngIII), which is further cleaved to AngIV. In our model, we skipped all the details about the enzymatic reactions AngII-AngIII-AngIV and kept only a single equation for their transformation. The dynamics of AngII and AngIV can thus be written as:
d [ AngII ] d t = ( c a c e + c c h y ) [ AngI ] c a c e 2 + c a n g I V + c a t 1 r + c a t 2 r [ AngII ] Log 2 h a n g I I [ AngII ]
d [ AngIV ] d t = c a n g I V [ AngII ] Log 2 h a n g I V [ AngIV ]
where h a n g I I and h a n g I V are the half-lives of the peptides and c a c e 2 , c a n g I V , c a t 1 r , and c a t 2 r the rates of the enzymatic reactions.
The dynamics of the peptide-bound form of the GPCRs are expressed as:
d [ AT 1 R AngII ] d t = c a t 1 r [ AngII ] Log 2 h a t 1 r [ AT 1 R AngII ]
d [ AT 2 R AngII ] d t = c a t 2 r [ AngII ] Log 2 h a t 2 r [ AT 2 R AngII ]
where [AT1R-AngII] and [AT1R-AngII] are the concentrations of the bound forms of the receptors and h a t 1 r and h a t 2 r their half-lives.
Until now, we have modeled the ACE/AngII/AT1R regulatory axis of RAS. Since the last two decades, it became clear that there is another RAS axis that acts as a counterregulator of the first axis [44]. The key role of this second axis is played by the Ang1-7 peptide that is mainly produced from AngII by the ACE2 enzyme and binds to the transmembrane GPCR called MAS. However, Ang1-7 can also be obtained as an enzymatic product from AngI via the catalytic activity of NEP and, to a lesser extent, from Ang1-9 via ACE and NEP. We overlooked the Ang1-9-related enzymatic reactions in our model, as they contribute less to Ang1-7 formation [33,34]. The dynamical equations for the Ang1-7 peptide and the MAS-bound receptor are as follows:
d [ Ang 1 7 ] d t = c n e p [ AngI ] + c a c e 2 [ AngII ] c m a s [ Ang 1 7 ] Log 2 h a n g 1 7 [ Ang 1 7 ]
d [ MAS Ang 1 7 ] d t = c m a s [ Ang 1 7 ] Log 2 h m a s [ MAS Ang 1 7 ]
Let us now go back to Equation (1) in which we simply expressed the renin production as a baseline term β . To describe the autoregulatory nature of RAS, this term has to depend on the production of other species, thus introducing a feedback regulation. It is known that this feedback depends on AT1R bound to AngII. Following other models [37,38], we express β as::
β = β 0 + [ AT 1 R AngII ] 0 N [ AT 1 R AngII ] δ 1
where β 0 is a constant parameter to be identified and [ AT 1 R AngII ] 0 N the equilibrium concentration for healthy normotensive humans. δ is a positive number that we fixed to 0.8 [37].
Technical details on the procedure used to solve the model and on model stability are given in Section 2.6 and Section 2.7.

2.2. Modeling Blood Pressure

Blood pressure is well known to be increased by the concentration of AngII bound to AT1R. It has also been described to be decreased by the concentration of MAS bound to Ang1-7 and of AT2R bound to AngII, but the precise mechanism is not yet known [45,46,47]. Therefore, we did not introduce in our model a feedback between these concentrations and renin production, as we did for AT1R-AngII, and modeled the diastolic blood pressure (DBP) simply from the AT1R-AngII concentration:
D B P = P 0 + P 1 [ AT 1 R AngII ]
We identified the two parameters P 0 and P 1 by fixing DBP equal to 80 mmHg for normotensive individuals and to 110 mmHg for hypertensive individuals. Hence, P 0 + P 1 [ AT 1 R AngII ] 0 N = 80 mmHg and P 0 + P 1 [ AT 1 R AngII ] 0 H = 110 mmHg, where the N and H superscripts denote the concentration in normotensive and hypertensive individuals and the 0 subscript the equilibrium concentrations.

2.3. Modeling RAS-Blocker Effects

Since dysregulated RAS with high levels of AngII is related to essential hypertension, a wide range of RAS-targeting drugs have been developed in the last forty years [48]. They can be classified into three different categories based on their pharmacological target [49]:
  • Angiotensin-converting enzyme inhibitors (ACE-I) that bind to ACE and thus inhibit the formation of angiotensin II and the associated vasoconstriction and inflammatory cascades. Examples of this type of drug are enalapril, lisinopril, and captopril.
  • Angiotensin receptor blockers (ARB) that block the binding of AngII to AT1R and thus act in antagonism with AngII. Examples are candesartan, losartan, and valsartan.
  • Direct renin inhibitors (DRI) that act on renin and thus inhibit the conversion of AGT to AngI. Examples are aliskiren, enalkiren, and remikiren.
We modeled the action of these three types of drugs by modifying the reaction rates associated with their targets as:
c a c e c a c e × 1 γ ACE I c a t 1 r c a t 1 r × 1 γ ARB c r e c r e × 1 γ DRI
where γ ACE I , γ ARB , and γ DRI are parameters describing the drug activity.

2.4. Modeling SARS-CoV-2 Infection

Since ACE2 is the entry point of SARS-CoV-2 [19], it is downregulated upon infection, and this impacts substantially the local and systemic RASs. In order to model the downregulation effect due to the virus, we modified the ACE2 reaction rate with the function γ CoV :
c a c e 2 c a c e 2 × 1 γ CoV ( C t )
We chose γ CoV to be a sigmoid function of the cycle threshold value C t , which is inversely related to the viral load [50]:
γ CoV = 1 1 + e a C t b
where a and b are positive real numbers. C t values of 31.5, 27.6, and 23.8 correspond to mild, moderate, and severe disease, respectively, and C t > 40 to undetected viral infection [51]. We thus chose the inflection point of the sigmoid at C t = 31.5 and imposed γ CoV > 0.99 for C t > 40 . Using these relations, we identified the values of the parameters a and b. They are reported in Table 1.

2.5. Modeling ARDS Severity

To model ARDS severity and how the lungs of SARS-CoV-2 patients evolve in response to RAS dysregulation, we introduced a phenomenological relation to estimate the PaO2/FiO2 ratio, defined as the ratio between the partial pressure of arterial oxygen (PaO2) and the fraction of inspired oxygen (FiO2). This quantity plays a key role in the assessment of ARDS patients [52,53]. The normal range of PaO2/FiO2 is between 400 and 500 mmHg. Mild and moderate ARDS are characterized by PaO2/FiO2 values in the range [200–300] mmHg and [100–200] mmHg, respectively. ARDS is severe for values below 100 mmHg.
We predicted the PaO2/FiO2 ratio as a function of the AngII and Ang1-7 concentrations:
P A O 2 / F I O 2 = A 0 + A 1 [ AngII ] [ AngII ] 0 + [ Ang 1 7 ] [ Ang 1 7 ] 0
where A 0 and A 1 are two parameters that we identified on the basis of our model by comparing the baseline RAS with the same system in which ACE2 is knocked out. In the former case, we fixed PaO2/FiO2 = 450 mmHg and, in the latter, PaO2/FiO2 = 50 mmHg.

2.6. Solving the RAS Model

The mathematical model of RAS described in Equations (1)–(10) is a system of ordinary differential equations (ODEs), which are linear except for the feedback loop of Equation (10).
We collected from the literature the values of the equilibrium concentrations of all proteins and peptides except renin and MAS bound to Ang1-7, for normotensive and hypertensive humans (Table 2). From these values, we fixed the parameters that appear in the phenomenological relations (11) and (15) for DBP and PaO2/FiO2 (Table 1). We also collected the values of the half-life of all proteins and peptides but MAS; we assumed the latter to be equal to that of the other membrane receptors (Table 1). Moreover, we estimated the value of reaction rate c r e from [36,54].
Using these concentration and parameter values, we solved the system of nine ODEs (Equations (1)–(9)) in the stationary state to identify the unknown parameters and concentrations. However, these equations have 12 unknowns: k a g t , β 0 , c a c e , c a c e 2 , c a n g I V , c a t 1 r , c a t 2 r , c m a s , c c h y , c n e p , [RE] 0 , and [MAS-Ang1-7] 0 . We had thus to assume three additional relations, which are:
c m a s = c a t 2 r
c c h y = 0
c n e p = 0
Since no quantitative data related to the MAS receptor can be found in the literature, we hypothesized the first relation assuming MAS and AT2R to be equally expressed and the affinity of Ang1-7 for MAS to be similar to the affinity of AngII for AT2R [46]. Moreover, we assumed c c h y = 0 and c n e p = 0 , but discuss the effect of non-vanishing values in Section 4.
By imposing these three additional relations, we solved the system of nine ODEs in the stationary state. The values obtained for [RE] 0 , [MAS-Ang1-7] 0 , k a g t , β 0 , c a c e , c a c e 2 , c a n g I V , c a t 1 r , and c a t 2 r for normotensive and hypertensive humans, are given in Table 2.

2.7. Stability of the RAS Model

The system of nine ODEs (Equations (1)–(9)) can be summarized in the form:
d x ( t ) d t = f ( x ( t ) , θ )
where x ( t ) is the vector containing the nine state variables, i.e., the concentrations of all proteins and peptides at time t, θ is the vector with all the production, kinetic, and half-life parameters, and f represents the vector that corresponds to the right-hand sides of Equations (1)–(9). In order to analyze the stability of the two steady states x 0 N and x 0 H for normotensive and hypertensive individuals, respectively, we computed the eigenvalues of the Jacobian matrix:
J ( x 0 ) = f ( x , θ ) x | x = x 0
where x 0 stands for either x 0 N or x 0 H .
In both the normotensive and hypertensive cases, seven strictly negative real values were obtained, together with two complex conjugate eigenvalues with strictly negative real parts. Both steady states x 0 N and x 0 H are therefore stable. The nonzero imaginary parts of the two complex conjugate eigenvalues are responsible for some damped oscillations in transient responses to parameter changes, but the overshoots are limited. It is interesting to note that the imaginary part is more than three times lower in the hypertensive case, hence leading to more damped responses in comparison with the normotensive case.
To quantify the state variable transients and the aforementioned overshoots, we simulated step responses corresponding to a 10% increase in the normal baseline for renin production β 0 . We observed some damped oscillations during the transient phase of the normotensive case, with very limited overshoots, e.g., 1.3% for the RE concentration. In the hypertensive case, the imaginary part of the complex conjugate eigenvalues is so low that the overshoots become almost undetectable (0.025%).

3. Results

The main objective of this paper is to investigate the effect of RAS-targeting drugs and SARS-CoV-2 infection, both individually and in combination, on the RAS of normotensive and hypertensive individuals. The robustness and predictive power of our model was first assessed by investigating the effects on RAS of three types of antihypertensive drugs: (i) ACE-I, (ii) ARB, and (iii) DRI (described in Section 2.3). This assessment included a comparison of model simulations with patient clinical data. Following the confirmation of model robustness and accuracy, ACE2 downregulation due to viral infection was introduced into the model to quantitatively predict how RAS is perturbed in COVID-19.

3.1. Model Predictions and Clinical Data on RAS-Blocker Drugs

The effect of enalapril, an ACE-I type drug, on plasma ACE activity and on plasma levels of AngI and AngII has been measured in normotensive individuals who received a single oral dose of 20 mg [60]. To compare these data with model predictions, we first fitted the γ ACE I parameter introduced in Equation (12) to the ACE activity values during enalapril administration divided by the pre-treatment activity (measured by an antibody-trapping assay). Once γ ACE I was set, we used our model to predict the dynamical response of RAS to this inhibitor drug. The time-dependent values of the AngI and AngII concentrations, normalized by their concentration at time 0, are shown in Figure 2a,b, both for our model predictions and experimental enalapril data; there is very good agreement between the two curves, without any further parameter fitting. The excellent correspondence between model prediction and experimental data is also clear from the root mean squared deviation (rmsd) between model prediction and experimental data on all time points following drug administration, as shown in Table 3.
Our model, thus, captures the known dynamics of ACE inhibition, (i.e., increased AngI levels and decreased AngII levels); this has the effect of lowering the concentration of AngII bound to AT1R and, thus, also lowers the blood pressure (Equation (11)).
To study the effect of ARB antihypertensive drugs on RAS, we considered data from [61], which measured the effects of different types of AT1R blocking molecules on the plasma levels of AngII in normotensive individuals. Specifically, the study participants received a single 50 mg dose of losartan, 80 mg of valsartan, or 150 mg of irbesartan. First, we fitted the γ ARB parameter (defined in Equation (12)) to the in vitro ability of the administered drug to induce the AngII receptor blockade, as measured by an AT1R radioreceptor binding assay [61]. We then used our model to predict the time-dependent AngI level, which was normalized by its concentration prior to drug administration. The results were evaluated through the rmsd between experimental and predicted values of AngI/AngI 0 at different time points after drug administration. The results, which are detailed in Table 3, clearly show that our model accurately predicts the RAS response to ARBs.
We also studied the effect of DRI-type drugs using experimental data that describe PRA activity and RE, AngI, and AngII concentrations, when different doses of aliskiren were administered orally to normotensive individuals [62]. We used the PRA activity data to fit the γ DRI parameter (introduced in Equation (12)), and we used our model to calculate the normalized AngI and AngII levels as a function of time. Here also, the results from our model and the experimental concentration data agree very well, as shown in Table 3.
In summary, the rmsd between predicted and experimental values of normalized AngI and AngII levels, averaged over all tested drugs, dosages, and a total of 38 time points, is 0.57 and 0.18, respectively (Table 3). These values should be compared with average experimental values of 1.7 and 0.5, respectively, demonstrating excellent agreement between experimental data and model predictions.
It should be noted that all reported experimental data were obtained after administration of single doses of RAS-targeting drugs. However, for hypertensive patients receiving long-term treatment, the expression of some enzymes involved in RAS could be either up- or down-regulated; we will return to this point in Section 4.
Finally, we compared model predictions against clinical data from large cohorts of patients describing the effect of ACE-I and ARB drug administration on blood pressure [63,64]. We first analyzed the response to ACE-I drugs alone. We plot measured DBP values averaged over more than ten ACE-I drug types as a function of the normalized dosage ξ ACE I [63] in Figure 2c, as well as predicted DBP values. We first considered a linear relation between γ ACE I and the normalized drug dosage ( γ ACE I = 0.5 ξ ACE I ). Despite this simplification, chosen to limit the number of parameters and thus overfitting issues, the curve reproduces the experimental data reasonably well. We also defined a non-linear relationship between these two quantities by introducing additional parameters: γ ACE I = 0.4 ξ ACE I 1 / 4 . We thus obtained a better fit as shown in Figure 2c.
We then studied the effect of the combined administration of the two drugs, ARB and ACE-I, on blood pressure, plotting the predicted DBP values as a function of both ξ ACE I and ξ ARB (see Figure 2d). We found that combined administration of ARB and ACE-I reduces DBP by 4 mmHg when compared with ARB monotherapy and by 12 mmHg when compared with ACE-I monotherapy. These predictions should be compared with clinical DBP values of 3 mmHg for combined administration compared to either monotherapy [64]. Thus, our model again provides an excellent prediction of experimental clinical data; further improvements to the model’s predictive strength are possible by fixing the γ ARB value at the maximum dose to be slightly lower than the corresponding γ ACE I value.

3.2. RAS in COVID-19

It is known that ACE2 is the cellular receptor of the spike glycoprotein of SARS-CoV-2 [9,10,11,12,13] and that it triggers the entry of SARS-COV-2 into the host cell. Although ACE2 is expressed in a variety of tissues [65,66,67], it is expressed mainly in the alveolar epithelial cells of lung, in the gastrointestinal tract, and in the kidney proximal tubular cells.
Here, we used our model to predict how RAS is perturbed by the SARS-CoV-2 virus. Simulation results of AngII and Ang1-7 concentrations and of the physiological value of PaO2/FiO2 as a function of SARS-CoV-2 viral load are presented in Figure 3 and in Table 4.
We observe that the AngII level increases with increasing viral load, with a much stronger increase for hypertensive than for normotensive patients. The AngII level is predicted to increase by approximately 15% for patients with moderate and severe COVID-19 (Table 4); this prediction is in very good agreement with the experimental value of 16% found in [68], but in poorer agreement with the value of 35% resulting from a study of only 12 patients [69].
We also observe that our model predicts a severe reduction of the Ang1-7 level, due to ACE2 downregulation; this reduction is the same for hypertensive and normotensive patients.
The overall result of the model is that RAS becomes imbalanced upon SARS-CoV-2 infection, with the harmful AngII axis upregulated and the counteracting Ang1-7 axis severely downregulated. This imbalance can be related to multiple clinical manifestations of COVID-19. More specifically, increased AngII levels cause hyperinflammation, which, in turn, increases plasma proinflammatory cytokine levels (in particular, IL-6) [70,71]. In addition, thrombotic events are observed, since AngII promotes the expression of plasminogen activator inhibitor-1 (PAI-1) and tissue-factors (TFs) [72,73]. Ang1-7, which normally counteracts these various effects [44], is downregulated by SARS-CoV-2 infection, such that COVID-19 clinical manifestations become increasingly severe as the disease develops.
Moreover, our model predicts severe ARDS with PaO2/FiO2 < 100 mmHg for normotensive and hypertensive patients whose C t values are smaller than 24.1 and 27.0, respectively. Our model predicts moderate ARDS, characterized by a PaO2/FiO2 ratio in the range of 100-200 mmHg, for normotensive and hypertensive patients having 24.1 < C t < 29.3 and 27.0 < C t < 29.7 , respectively, and mild ARDS, characterized by a PaO2/FiO2 ratio in the range of 200–300 mmHg for normotensive and hypertensive patients having 29.3 < C t < 31.4 and 29.7 < C t < 31.6 , respectively.
Our modeling approach suggests a weak relationship between hypertension and ARDS severity resulting from SARS-CoV-2 infection. The mean value of the P A O 2 / F I O 2 ratio over the entire C t range is approximately 20 mmHg lower for hypertensive than for normotensive patients. Indeed, the large difference in AngII levels between normotensive and hypertensive patients is partially compensated by the absence of any difference in Ang1-7 levels.

3.3. Impact of RAS-Modulating Drugs on COVID-19 Severity

We analyzed the effect of administering a selection of drugs to normotensive and hypertensive patients who were infected with SARS-CoV-2. More specifically, we considered RAS-blocking drugs that are already commonly used to treat hypertension, as well as drugs that are currently undergoing clinical trials in the context of COVID-19, such as rhACE2 and Ang1-7.
  • Antihypertensive RAS-blocking drugs: We combined the effect of each of the three RAS-blocking ACE-I, ARB, and DRI drugs, which were modeled by the enzyme-inhibiting γ functions (introduced in Equation (12)), with the ACE2-inhibiting C t -dependent γ CoV function (defined in Equation (14)), which mimics SARS-CoV-2 infection. the PaO2/FiO2 values predicted by our model are presented in Figure 4.
Our model predicts that administration of ACE-I and DRI drugs protect from the adverse effects of ARDS, especially for hypertensive patients, while ARB drugs are predicted to worsen ARDS severity, especially for normotensive patients.
Model predictions for ACE inhibitors are in agreement with clinical data, which indicate that treatment with ACE inhibitors is associated with better survival among COVID-19 patients [31,74]. Indeed, only 3% of non-surviving COVID-19 patients that were monitored were treated with ACE-I drugs compared to 9% of surviving COVID-19 patients [31]. Moreover, in a meta-analysis [74], hypertensive patients treated with ACE-I drugs were associated with a reduced mortality of 35% when compared to patients who were not treated with ACE-I drugs. in another clinical analysis [75], older patients who were treated with ACE-I drugs had a 40% lower risk of hospitalization than those who were not treated with ACE-I drugs.
No data are currently available to validate our model prediction that COVID-19 attenuation due to ACE-I drug treatment is stronger in hypertensive than in normotensive patients. Furthermore, no data are currently available to validate our model prediction that DRI and ACE-I drug treatments cause similar levels of COVID-19 disease attenuation.
In contrast to DRI and ACE-I drugs, our model predicts that ARB drug treatment worsens COVID-19 severity, with the effect being stronger for normotensive compared to hypertensive patients. Here, the agreement between model predictions and clinical data is less clear, with some clinical data in agreement with our model prediction [31,75], while other clinical data suggest that ARB drug treatment does not affect hospitalization risk [75] or mortality [74,76]. This lack of agreement must be further investigated with additional clinical data.
Moreover, we performed a quantitative prediction of the drug effects on disease severity by calculating the RAS peptide concentrations, PaO2/FiO2 values, and DPB for moderate COVID-19 patients. The results are presented in Table 5.
Administration of ACE-I drugs, modeled by γ A C E I = 0.5 , increases the P A O 2 / F I O 2 value by approximately 50 and 70 mmHg for normotensive and hypertensive patients, respectively. an equivalent administration of DRI drugs increases this ratio even more, by 70 and 150 mmHg, while ARB administration decreases it by 140 and 30 mmHg for normotensive and hypertensive patients, respectively.
The opposite effect of ARBs administration compared to ACE-I and DRI drugs can be attributed to the substantial increase in AngII concentration, which is only partially balanced by a relatively small increase in Ang1-7 concentration, given that ACE2 is downregulated in SARS-CoV-2 infection.
Note that a number of ARB drugs, including valsartan and losartan, are currently being tested in clinical trials, with the hope that they will rescue RAS in COVID-19 patients [25,26,27]. Our model predicts that this will not be the case.
Finally, as shown in Table 5, the blood pressure is predicted to be unaffected by the administration of either ACE-I, ARB, or DRI to normotensive COVID-19 patients, but to be reduced by approximately 10–20 mmHg by administration to hypertensive patients.
  • Other RAS-targeting drugs: We used our model to test the potential of other drugs that are currently in clinical trials to restore the functional activity of the perturbed RAS upon viral infection. First, we modeled how the administration of an exogenous supplement of rhACE2 (GSK2586881) affects RAS by modifying the reaction rate c a c e 2 defined in Equation (13). This rate already includes the function γ CoV that mimics SARS-CoV-2 infection, and we simply added a second function γ rhACE 2 associated with the effects of rhACE2 administration:
    c a c e 2 c a c e 2 × 1 + γ rhACE 2 γ CoV ( C t )
Our model predicts an increase in P A O 2 / F I O 2 following the administration of exogenous rhACE2, thus predicting an alleviation of disease severity, as shown in Figure 5 and Table 5. Specifically, P A O 2 / F I O 2 is predicted to increase by approximately 200 mmHg when γ rhACE 2 is fixed at 0.5. Our model also predicts, as expected, a reduction in AngII level and an increase in Ang1-7 level.
These predictions are in agreement with both animal and in vitro studies [18,22], whereby rhACE2 is observed to alleviate virus-related ARDS severity through a double action. First, by rhACE2 binding to the virus spike protein, interaction with endogenous ACE2 is prevented, and infection is slowed down. Second, rhACE2 administration increases ACE2 activity, thus causing a reduction in AngII level and an increase in Ang1-7 level; this protects lung against severe failure.
Current clinical trial data concerning the administration of different doses of rhACE2 (0.1, 0.2, 0.4, and 0.8 mg/kg) to SARS-CoV patients at different time intervals (2, 4, and 18 h) are only in partial agreement with our model predictions [20]. Specifically, while clinical data followed the predicted decrease in [AngII] and the predicted increase in [Ang1-7], there was no sustained increase in P A O 2 / F I O 2 compared with placebo. It has been suggested that the drug concentrations used in these clinical trials were too low to have a measurable effect on the respiratory system and that drug administration via infusion would have been more sustained [20]. More experimental and clinical data are clearly needed to further investigate the effect of rhACE2 on coronavirus-related ARDS.
Another method of boosting the second RAS axis, ACE2/Ang1-7/MAS, which is downregulated by SARS-CoV-2 infection, is to administer Ang1-7 peptides as a means of triggering anti-inflammatory and anti-fibrotic mechanisms. We modeled Ang1-7 peptide administration by introducing a new parameter, the production rate η A n g 17 , to the dynamical Equation (8) of [Ang1-7]; this allows the model to describe the exogenous Ang1-7 level, which is added to the endogenous Ang1-7 baseline. as shown in Figure 5b and Table 5, our model predicts a clear alleviation of COVID-19 severity, with P A O 2 / F I O 2 increasing by 50 and 130 mmHg for hypertensive and normotensive patients, respectively, upon administration of η A n g 17 = 25 fmol/(mL min) Ang1-7 in infusion. Note that COVID-19 alleviation is significantly stronger in normotensive compared to hypertensive patients for the same drug concentrations; a slightly stronger concentration of Ang1-7 must be administered to hypertensive patients for an equivalent effect.
Our model predicts a quantitative reduction in ARDS severity in COVID-19 patients, in agreement with the known anti-inflammation and anti-fibrosis nature of Ang1-7. Model predictions nicely agree with data from animal studies without the need for any additional fitting. For example, administration of Ang1-7 by infusion to acid-injured rats suffering from ARDS increases the baseline Ang1-7 level by a factor 2.5, leading to an increase in P A O 2 / F I O 2 of approximately 70 mmHg [77]. However, while the P A O 2 / F I O 2 value increases linearly in our model as a function of Ang1-7 concentration, it reaches a plateau in rats; this suggests that our model is probably oversimplified, since P A O 2 / F i O 2 is not a linear function of Ang1-7 concentration. Further work on this aspect of our model will be possible when more data become available.

4. Discussion

The spike protein of SARS-CoV-2 interferes with RAS by binding to the ACE2 receptor, a key element of RAS. Despite recent progress in understanding COVID-perturbed RAS and how its functionality can be restored, more work is urgently needed in the context of the current COVID-19 pandemic.
We here present a simple computational approach to modeling RAS evolution in the context of SARS-CoV-2 infection. Inspired by a number of existing RAS models [35,36,37,38], we searched the literature for measured half-lives and concentrations of angiotensin peptides and their receptors in healthy normotensive and hypertensive individuals and then identified the unknown production and reaction rate parameters from the model. As an initial test, we compared our model predictions of how the administration of RAS-blocking drugs would affect Ang peptide concentrations and blood pressure with relevant experimental data; we found good quantitative agreement between our model and experimental data, without the need for further parameter fitting. We then modeled the effect of SARS-CoV-2 infection on RAS through the downregulation of ACE2, which we related to the SARS-CoV-2 viral load.
A focal point of our work was to investigate how a series of RAS-targeting drugs affected COVID-19 patients. We found that the administration of two antihypertensive drugs, ACE-I and DRI, tended to reduce the severity of COVID-19, while ARB drugs worsened it. Clinical data generally support the model’s predictions for the administration of ACE-I drugs, but they are either absent or partially contradict the model’s prediction for DRI and ARB administration. Additionally, we modeled a potential treatment that is currently under clinical trial in COVID-19 patients: administration of rhACE2 or Ang1-7 by drug infusion. Our model predicts improved clinical outcomes in these cases, in agreement with a series of experimental data on animal models.
It is important to note that, despite its simplicity, our model has excellent accuracy in reproducing clinical and experimental data on the perturbed RAS. Furthermore, the model’s predictions of changes in COVID-19 severity due to drug administration are blind predictions, without the fitting of any additional parameters.
Many challenges remain in our current understanding of RAS perturbation in COVID-19 patients. Importantly, more data regarding angiotensin peptide concentrations upon SARS-CoV-2 infection are urgently needed, since currently available data are often inconsistent or conflicting so that reliable comparisons between model predictions and experimental data cannot be made. Even in healthy individuals, angiotensin peptide levels can vary substantially due to their low circulating concentrations, the experimental techniques used to measure them, and inter-patient variability.
When developing our model, we chose not to consider two enzymes that are active in RAS through the cancellation of their reaction rates: CHY and NEP (see Equations (17) and (18)). The CHY enzyme is expressed in mast cells present in interstitial lung connective tissues, and it cleaves AngI to form AngII. the addition of this enzymatic reaction in the model would not really influence the predictions since it would essentially be a reparametrization of ACE activity and of ACE-I action. It might, nevertheless, be interesting to add the CHY enzymatic reaction, which yields ACE-independent synthesis of AngII and has been suggested (although debated) to be upregulated in the case of long-term ACE-I administration [78]; this would enable an explanation of why ACE-I fails to inhibit AngII formation after some time [78,79].
The NEP enzyme is expressed in a wide range of tissues, being particularly abundant in kidney, and it cleaves AngI to form Ang1-7. It influences the counterregulatory RAS axis through its connection to Ang1-7 levels, thus affecting COVID-19 severity. However, NEP’s role is far from clear, and the literature contains contradictory findings. Experimental data from rats with ARDS suggest that NEP is severely downregulated in both plasma and lung tissues [80]. Note that NEP also cleaves natriuretic peptides, which have both anti-inflammatory and anti-fibrotic effects [81]. Therefore, the combined administration of NEP-inhibiting and ARB drugs has been suggested to treat SARS-CoV-2 patients [82].
Our future work will include building more complexity into our model by explicitly considering the communication between local and systemic RASs [33,34] and by including the interaction between RAS and the immune system [83]. This model extension is necessary for an improved quantitative understanding of RAS dysregulation upon a variety of perturbations, including SARS-CoV-2 infection.
Our model and its predictions provide a valuable and robust framework for in silico testing of hypotheses regarding COVID-19’s pathogenic mechanisms and the effect of drugs that are aimed at restoring RAS functionality. Our work also opens a broader discussion on the role of the full RAS in COVID-19, a topic that has received little attention to date, perhaps due to the current focus on the ACE2 enzyme, which, although very important as it is directly targeted by the virus, constitutes only one part of a much more complex system.

Author Contributions

Conceptualization, F.P., P.B. and M.R.; formal analysis and investigation F.P. and M.R.; methodology and validation, F.P., P.B. and M.R.; writing–original draft preparation, F.P. and M.R.; writing–review and editing F.P., P.B. and M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by the F.R.S.-FNRS Fund for Scientific Research through a COVID—Exceptional Research Project.

Acknowledgments

We thank Filippo Annoni and Fabio Taccone for enlightening discussions. F.P. and M.R. are Scientific Collaborator and Research Director, respectively, at the F.R.S.-FNRS Fund for Scientific Research.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Data Availability

The code used to generate all the results of this paper is freely available on the GitHub repository (https://github.com/3BioCompBio/RASinCOVID).

References

  1. Dong, E.; Du, H.; Gardner, L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect. Dis. 2020, 20, 533–534. [Google Scholar] [CrossRef]
  2. Chan, J.F.W.; Yuan, S.; Kok, K.H.; To, K.K.W.; Chu, H.; Yang, J.; Xing, F.; Liu, J.; Yip, C.C.Y.; Poon, R.W.S.; et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster. Lancet 2020, 395, 514–523. [Google Scholar] [CrossRef] [Green Version]
  3. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [Green Version]
  4. Chen, N.; Zhou, M.; Dong, X.; Qu, J.; Gong, F.; Han, Y.; Qiu, Y.; Wang, J.; Liu, Y.; Wei, Y.; et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study. Lancet 2020, 395, 507–513. [Google Scholar] [CrossRef] [Green Version]
  5. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.M.; Wang, W.; Song, Z.G.; Hu, Y.; Tao, Z.W.; Tian, J.H.; Pei, Y.Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [Green Version]
  6. Andersen, K.G.; Rambaut, A.; Lipkin, W.I.; Holmes, E.C.; Garry, R.F. The proximal origin of SARS-CoV-2. Nat. Med. 2020, 26, 450–452. [Google Scholar] [CrossRef] [Green Version]
  7. Benvenuto, D.; Giovanetti, M.; Ciccozzi, A.; Spoto, S.; Angeletti, S.; Ciccozzi, M. The 2019-new coronavirus epidemic: Evidence for virus evolution. J. Med. Virol. 2020, 92, 455–459. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, T.; Wu, Q.; Zhang, Z. Probable pangolin origin of SARS-CoV-2 associated with the COVID-19 outbreak. Curr. Biol. 2020, 30, 1346–1351.e2. [Google Scholar] [CrossRef]
  9. 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] [Green Version]
  10. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, 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]
  11. Zhang, H.; Penninger, J.M.; Li, Y.; Zhong, N.; Slutsky, A.S. Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: Molecular mechanisms and potential therapeutic target. Intensive Care Med. 2020, 46, 586–590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Cao, Y.; Li, L.; Feng, Z.; Wan, S.; Huang, P.; Sun, X.; Wen, F.; Huang, X.; Ning, G.; Wang, W. Comparative genetic analysis of the novel coronavirus (2019-nCoV/SARS-CoV-2) receptor ACE2 in different populations. Cell Discov. 2020, 6, 1–4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lu, R.; Zhao, X.; Li, J.; Niu, P.; Yang, B.; Wu, H.; Wang, W.; Song, H.; Huang, B.; Zhu, N.; et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding. Lancet 2020, 395, 565–574. [Google Scholar] [CrossRef] [Green Version]
  14. Burrell, L.M.; Johnston, C.I.; Tikellis, C.; Cooper, M.E. ACE2, a new regulator of the renin–angiotensin system. Trends Endocrinol. Metab. 2004, 15, 166–169. [Google Scholar] [CrossRef]
  15. Donoghue, M.; Hsieh, F.; Baronas, E.; Godbout, K.; Gosselin, M.; Stagliano, N.; Donovan, M.; Woolf, B.; Robison, K.; Jeyaseelan, R.; et al. A novel angiotensin-converting enzyme–related carboxypeptidase (ACE2) converts angiotensin I to angiotensin 1-9. Circ. Res. 2000, 87, e1–e9. [Google Scholar] [CrossRef]
  16. Tipnis, S.R.; Hooper, N.M.; Hyde, R.; Karran, E.; Christie, G.; Turner, A.J. A human homolog of angiotensin- converting enzyme cloning and functional expression as a captopril-insensitive carboxypeptidase. J. Biol. Chem. 2000, 275, 33238–33243. [Google Scholar] [CrossRef] [Green Version]
  17. Imai, Y.; Kuba, K.; Rao, S.; Huan, Y.; Guo, F.; Guan, B.; Yang, P.; Sarao, R.; Wada, T.; Leong-Poi, H.; et al. Angiotensin-converting enzyme 2 protects from severe acute lung failure. Nature 2005, 436, 112–116. [Google Scholar] [CrossRef]
  18. 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 angiotensin converting enzyme 2 (ACE2) in SARS coronavirus–induced lung injury. Nat. Med. 2005, 11, 875–879. [Google Scholar] [CrossRef]
  19. Li, W.; Moore, M.J.; Vasilieva, N.; Sui, J.; Wong, S.K.; Berne, M.A.; Somasundaran, M.; Sullivan, J.L.; Luzuriaga, K.; Greenough, T.C.; et al. Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature 2003, 426, 450–454. [Google Scholar] [CrossRef] [Green Version]
  20. Khan, A.; Benthin, C.; Zeno, B.; Albertson, T.E.; Boyd, J.; Christie, J.D.; Hall, R.; Poirier, G.; Ronco, J.J.; Tidswell, M.; et al. A pilot clinical trial of recombinant human angiotensin-converting enzyme 2 in acute respiratory distress syndrome. Crit. Care 2017, 21, 1–9. [Google Scholar] [CrossRef] [Green Version]
  21. Gheblawi, M.; Wang, K.; Viveiros, A.; Nguyen, Q.; Zhong, J.C.; Turner, A.J.; Raizada, M.K.; Grant, M.B.; Oudit, G.Y. Angiotensin-converting enzyme 2: SARS-CoV-2 receptor and regulator of the renin-angiotensin system: Celebrating the 20th anniversary of the discovery of ACE2. Circ. Res. 2020, 126, 1456–1474. [Google Scholar] [CrossRef] [PubMed]
  22. Monteil, V.; Kwon, H.; Prado, P.; Hagelkrüys, A.; Wimmer, R.A.; Stahl, M.; Leopoldi, A.; Garreta, E.; Del Pozo, C.H.; Prosper, F.; et al. Inhibition of SARS-CoV-2 infections in engineered human tissues using clinical-grade soluble human ACE2. Cell 2020, 181, 905–913.e7. [Google Scholar] [CrossRef] [PubMed]
  23. NCT04287686. Recombinant Human Angiotensin-Converting Enzyme 2 (rhACE2) as a Treatment for Patients with COVID-19. Available online: https://clinicaltrials.gov/ct2/show/NCT04287686 (accessed on 27 February 2020).
  24. NCT04332666. Angiotensin-(1,7) Treatment in COVID-19: The ATCO Trial (ATCO). Available online: https://www.clinicaltrials.gov/ct2/show/NCT04332666 (accessed on 3 April 2020).
  25. NCT04335786. Valsartan for Prevention of Acute Respiratory Distress Syndrome in Hospitalized Patients With SARS-COV-2 (COVID-19) Infection Disease. Available online: https://clinicaltrials.gov/ct2/show/NCT04335786 (accessed on 6 April 2020).
  26. NCT04312009. Losartan for Patients with COVID-19 Requiring Hospitalization. Available online: https://clinicaltrials.gov/ct2/show/NCT04312009 (accessed on 17 March 2020).
  27. NCT04311177. Losartan for Patients with COVID-19 Not Requiring Hospitalization. Available online: https://clinicaltrials.gov/ct2/show/NCT04311177 (accessed on 17 March 2020).
  28. NCT04318418. ACE Inhibitors, Angiotensin II Type-I Receptor Blockers and Severity of COVID-19 (CODIV-ACE). Available online: https://clinicaltrials.gov/ct2/show/NCT04318418 (accessed on 24 March 2020).
  29. Reynolds, H.R.; Adhikari, S.; Pulgarin, C.; Troxel, A.B.; Iturrate, E.; Johnson, S.B.; Hausvater, A.; Newman, J.D.; Berger, J.S.; Bangalore, S.; et al. Renin–angiotensin–aldosterone system inhibitors and risk of Covid-19. N. Engl. J. Med. 2020, 382, 2441–2448. [Google Scholar] [CrossRef] [PubMed]
  30. Mancia, G.; Rea, F.; Ludergnani, M.; Apolone, G.; Corrao, G. Renin–angiotensin–aldosterone system blockers and the risk of Covid-19. N. Engl. J. Med. 2020, 382, 2431–2440. [Google Scholar] [CrossRef] [PubMed]
  31. Mehra, M.R.; Desai, S.S.; Kuy, S.; Henry, T.D.; Patel, A.N. Cardiovascular disease, drug therapy, and mortality in COVID-19. N. Engl. J. Med. 2020, 382, e102. [Google Scholar] [CrossRef]
  32. Paul, M.; Poyan Mehr, A.; Kreutz, R. Physiology of local renin-angiotensin systems. Physiol. Rev. 2006, 86, 747–803. [Google Scholar] [CrossRef]
  33. Raizada, M.K.; Phillips, M.I.; Sumners, C. Cellular and Molecular Biology of the Renin-Angiotensin System; CRC Press: Boca Raton, FL, USA, 1993. [Google Scholar]
  34. Casarini, D.E.; Arita, D.Y.; Cunha, T.S.; Colucci, J.A. New Aspects of the Renin Angiotensin System in Cardiovascular and Renal Diseases; Bentham: Sharjah, UAE, 2016. [Google Scholar]
  35. Hallow, K.M.; Lo, A.; Beh, J.; Rodrigo, M.; Ermakov, S.; Friedman, S.; de Leon, H.; Sarkar, A.; Xiong, Y.; Sarangapani, R.; et al. A model-based approach to investigating the pathophysiological mechanisms of hypertension and response to antihypertensive therapies: Extending the Guyton model. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2014, 306, R647–R662. [Google Scholar] [CrossRef]
  36. Versypt, A.N.F.; Harrell, G.K.; McPeak, A.N. A pharmacokinetic/pharmacodynamic model of ACE inhibition of the renin-angiotensin system for normal and impaired renal function. Comput. Chem. Eng. 2017, 104, 311–322. [Google Scholar] [CrossRef]
  37. Leete, J.; Gurley, S.; Layton, A.T. Modeling sex differences in the renin angiotensin system and the efficacy of antihypertensive therapies. Comput. Chem. Eng. 2018, 112, 253–264. [Google Scholar] [CrossRef]
  38. Leete, J.; Layton, A.T. Sex-specific long-term blood pressure regulation: Modeling and analysis. Comput. Biol. Med. 2019, 104, 139–148. [Google Scholar] [CrossRef] [Green Version]
  39. Ruiz-Ortega, M.; Lorenzo, O.; Ruperez, M.; Esteban, V.; Suzuki, Y.; Mezzano, S.; Plaza, J.; Egido, J. Role of the renin-angiotensin system in vascular diseases: Expanding the field. Hypertension 2001, 38, 1382–1387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. De Man, F.S.; Tu, L.; Handoko, M.L.; Rain, S.; Ruiter, G.; François, C.; Schalij, I.; Dorfmüller, P.; Simonneau, G.; Fadel, E.; et al. Dysregulated renin–angiotensin–aldosterone system contributes to pulmonary arterial hypertension. Am. J. Respir. Crit. Care Med. 2012, 186, 780–789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Jia, G.; Aroor, A.R.; Hill, M.A.; Sowers, J.R. Role of renin-angiotensin-aldosterone system activation in promoting cardiovascular fibrosis and stiffness. Hypertension 2018, 72, 537–548. [Google Scholar] [CrossRef] [PubMed]
  42. Te Riet, L.; van Esch, J.H.; Roks, A.J.; van den Meiracker, A.H.; Danser, A.J. Hypertension: Renin–angiotensin– aldosterone system alterations. Circ. Res. 2015, 116, 960–975. [Google Scholar] [CrossRef] [PubMed]
  43. Kobori, H.; Nangaku, M.; Navar, L.G.; Nishiyama, A. The intrarenal renin-angiotensin system: From physiology to the pathobiology of hypertension and kidney disease. Pharmacol. Rev. 2007, 59, 251–287. [Google Scholar] [CrossRef]
  44. Simões e Silva, A.; Silveira, K.; Ferreira, A.; Teixeira, M. ACE2, angiotensin-(1-7) and M as receptor axis in inflammation and fibrosis. Br. J. Pharmacol. 2013, 169, 477–492. [Google Scholar] [CrossRef] [Green Version]
  45. Povlsen, A.L.; Grimm, D.; Wehland, M.; Infanger, M.; Krüger, M. The vasoactive Mas receptor in essential hypertension. J. Clin. Med. 2020, 9, 267. [Google Scholar] [CrossRef] [Green Version]
  46. Santos, R.A.; e Silva, A.C.S.; Maric, C.; Silva, D.M.; Machado, R.P.; de Buhr, I.; Heringer-Walther, S.; Pinheiro, S.V.B.; Lopes, M.T.; Bader, M.; et al. Angiotensin-(1–7) is an endogenous ligand for the G protein-coupled receptor Mas. Proc. Natl. Acad. Sci. USA 2003, 100, 8258–8263. [Google Scholar] [CrossRef] [Green Version]
  47. Carey, R.M. AT2 receptors: Potential therapeutic targets for hypertension. Am. J. Hypertens. 2017, 30, 339–347. [Google Scholar]
  48. Zaman, M.A.; Oparil, S.; Calhoun, D.A. Drugs targeting the renin–angiotensin–aldosterone system. Nat. Rev. Drug Discov. 2002, 1, 621–636. [Google Scholar] [CrossRef]
  49. Williams, B. Drug discovery in renin–angiotensin system intervention: Past and future. Ther. Adv. Cardiovasc. Dis. 2016, 10, 118–125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Borg, I.; Rohde, G.; Löseke, S.; Bittscheidt, J.; Schultze-Werninghaus, G.; Stephan, V.; Bufe, A. Evaluation of a quantitative real-time PCR for the detection of respiratory syncytial virus in pulmonary diseases. Eur. Respir. J. 2003, 21, 944–951. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Zheng, S.; Fan, J.; Yu, F.; Feng, B.; Lou, B.; Zou, Q.; Xie, G.; Lin, S.; Wang, R.; Yang, X.; et al. Viral load dynamics and disease severity in patients infected with SARS-CoV-2 in Zhejiang province, China, January-March 2020: Retrospective cohort study. BMJ 2020, 369, m1443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Villar, J.; Pérez-Méndez, L.; Blanco, J.; Añón, J.M.; Blanch, L.; Belda, J.; Santos-Bouza, A.; Fernández, R.L.; Kacmarek, R.M. A universal definition of ARDS: The PaO2/FiO2 ratio under a standard ventilatory setting—a prospective, multicenter validation study. Intensive Care Med. 2013, 39, 583–592. [Google Scholar] [CrossRef]
  53. Ware, L.B.; Matthay, M.A. The acute respiratory distress syndrome. New Engl. J. Med. 2000, 342, 1334–1349. [Google Scholar] [CrossRef]
  54. Streatfeild-James, R.M.; Williamson, D.; Pike, R.N.; Tewksbury, D.; Carrell, R.W.; Coughlin, P.B. Angiotensinogen cleavage by renin: Importance of a structurally constrained N-terminus. FEBS Lett. 1998, 436, 267–270. [Google Scholar] [CrossRef] [Green Version]
  55. Katsurada, A.; Hagiwara, Y.; Miyashita, K.; Satou, R.; Miyata, K.; Ohashi, N.; Navar, L.G.; Kobori, H. Novel sandwich ELISA for human angiotensinogen. Am. J. Physiol. Ren. Physiol. 2007, 293, F956–F960. [Google Scholar] [CrossRef] [Green Version]
  56. Chappell, M.C. Biochemical evaluation of the renin-angiotensin system: The good, bad, and absolute? Am. J. Physiol. Heart Circ. Physiol. 2016, 310, H137–H152. [Google Scholar] [CrossRef] [Green Version]
  57. Pendergrass, K.D.; Pirro, N.T.; Westwood, B.M.; Ferrario, C.M.; Brosnihan, K.B.; Chappell, M.C. Sex differences in circulating and renal angiotensins of hypertensive mRen. Lewis but not normotensive Lewis rats. Am. J. Physiol. Heart Circ. Physiol. 2008, 295, H10–H20. [Google Scholar] [CrossRef] [Green Version]
  58. Sullivan, J.C.; Rodriguez-Miguelez, P.; Zimmerman, M.A.; Harris, R.A. Differences in angiotensin (1–7) between men and women. Am. J. Physiol.-Heart Circ. Physiol. 2015, 308, H1171–H1176. [Google Scholar] [CrossRef]
  59. Nussberger, J.; Brunner, D.B.; Waeber, B.; Brunner, H.R. Specific measurement of angiotensin metabolites and in vitro generated angiotensin II in plasma. Hypertension 1986, 8, 476–482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Nussberger, J.; Brunner, D.; Keller, I.; Brunner, H.R. Measurement of converting enzyme activity by antibody-trapping of generated angiotensin II: Comparison with two other methods. Am. J. Hypertens. 1992, 5, 393–398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Mazzolai, L.; Maillard, M.; Rossat, J.; Nussberger, J.; Brunner, H.R.; Burnier, M. Angiotensin II receptor blockade in normotensive subjects: A direct comparison of three AT1 receptor antagonists. Hypertension 1999, 33, 850–855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Nussberger, J.; Wuerzner, G.; Jensen, C.; Brunner, H.R. Angiotensin II Suppression in humans by the orally active renin inhibitor aliskiren (SPP100) comparison with enalapril. Hypertension 2002, 39, e1–e8. [Google Scholar] [CrossRef] [Green Version]
  63. Heran, B.S.; Wong, M.M.; Heran, I.K.; Wright, J.M. Blood pressure lowering efficacy of angiotensin converting enzyme (ACE) inhibitors for primary hypertension. Cochrane Database Syst. Rev. 2008. [Google Scholar] [CrossRef]
  64. Doulton, T.W.; He, F.J.; MacGregor, G.A. Systematic review of combined angiotensin-converting enzyme inhibition and angiotensin receptor blockade in hypertension. Hypertension 2005, 45, 880–886. [Google Scholar] [CrossRef] [Green Version]
  65. Chen, L.; Li, X.; Chen, M.; Feng, Y.; Xiong, C. The ACE2 expression in human heart indicates new potential mechanism of heart injury among patients infected with SARS-CoV-2. Cardiovasc. Res. 2020, 116, 1097–1100. [Google Scholar] [CrossRef] [Green Version]
  66. Xu, H.; Zhong, L.; Deng, J.; Peng, J.; Dan, H.; Zeng, X.; Li, T.; Chen, Q. High expression of ACE2 receptor of 2019-nCoV on the epithelial cells of oral mucosa. Int. J. Oral Sci. 2020, 12, 1–5. [Google Scholar] [CrossRef]
  67. Sungnak, W.; Huang, N.; Bécavin, C.; Berg, M.; Queen, R.; Litvinukova, M.; Talavera-López, C.; Maatz, H.; Reichart, D.; Sampaziotis, F.; et al. SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes. Nat. Med. 2020, 26, 681–687. [Google Scholar] [CrossRef] [Green Version]
  68. Liu, N.; Hong, Y.; Chen, R.G.; Zhu, H.M. High rate of increased level of plasma Angiotensin II and its gender difference in COVID-19: An analysis of 55 hospitalized patients with COVID-19 in a single hospital, WuHan, China. medRxiv 2020. [Google Scholar] [CrossRef]
  69. Liu, Y.; Yang, Y.; Zhang, C.; Huang, F.; Wang, F.; Yuan, J.; Wang, Z.; Li, J.; Li, J.; Feng, C.; et al. Clinical and biochemical indexes from 2019-nCoV infected patients linked to viral loads and lung injury. Sci. China Life Sci. 2020, 63, 364–374. [Google Scholar] [CrossRef] [Green Version]
  70. Merad, M.; Martin, J.C. Pathological inflammation in patients with COVID-19: A key role for monocytes and macrophages. Nat. Rev. Immunol. 2020, 20, 355–362. [Google Scholar] [CrossRef] [PubMed]
  71. Satou, R.; Penrose, H.; Navar, L.G. Inflammation as a regulator of the renin-angiotensin system and blood pressure. Curr. Hypertens. Rep. 2018, 20, 100. [Google Scholar] [CrossRef] [PubMed]
  72. Vaughan, D.E.; Lazos, S.A.; Tong, K. Angiotensin II regulates the expression of plasminogen activator inhibitor-1 in cultured endothelial cells. A potential link between the renin-angiotensin system and thrombosis. J. Clin. Investig. 1995, 95, 995–1001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Vaughan, D.E. PAI-1 and atherothrombosis. J. Thromb. Haemost. 2005, 3, 1879–1883. [Google Scholar] [CrossRef]
  74. Ssentongo, A.; Ssentongo, P.; Heilbrunn, E.S.; Lekoubou, A.; Du, P.; Liao, D.; Oh, J.S.; Chinchilli, V.M. Renin-angiotensin-aldosterone system inhibitors and mortality in patients with hypertension hospitalized for COVID-19: A systematic review and meta-analysis. medRxiv 2020. [Google Scholar] [CrossRef]
  75. Khera, R.; Clark, C.; Lu, Y.; Guo, Y.; Ren, S.; Truax, B.; Spatz, E.S.; Murugiah, K.; Lin, Z.; Omer, S.B.; et al. Association of Angiotensin-Converting Enzyme Inhibitors and Angiotensin Receptor Blockers with the Risk of Hospitalization and Death in Hypertensive Patients with Coronavirus Disease-19. medRxiv 2020. [Google Scholar] [CrossRef]
  76. Baral, R.; White, M.; Vassiliou, V.S. Impact of hospitalised patients with COVID-19 taking Renin-Angiotensin-Aldosterone System inhibitors: A systematic review and meta-analysis. medRxiv 2020. [Google Scholar] [CrossRef]
  77. Zambelli, V.; Bellani, G.; Borsa, R.; Pozzi, F.; Grassi, A.; Scanziani, M.; Castiglioni, V.; Masson, S.; Decio, A.; Laffey, J.G.; et al. Angiotensin-(1-7) improves oxygenation, while reducing cellular infiltrate and fibrosis in experimental Acute Respiratory Distress Syndrome. Intensive Care Med. Exp. 2015, 3, 1–17. [Google Scholar] [CrossRef] [Green Version]
  78. Chester, A.; Borland, J. Chymase-dependent angiotensin II formation in human blood vessels. J. Hum. Hypertens. 2000, 14, 373–376. [Google Scholar] [CrossRef] [Green Version]
  79. Athyros, V.G.; Mikhailidis, D.P.; Kakafika, A.I.; Tziomalos, K.; Karagiannis, A. Angiotensin II reactivation and aldosterone escape phenomena in renin–angiotensin–aldosterone system blockade: Is oral renin inhibition the solution? Expert Opin. Pharmacother. 2007, 8, 529–535. [Google Scholar] [CrossRef] [PubMed]
  80. Hashimoto, S.; Amaya, F.; Oh-hashi, K.; Kiuchi, K.; Hashimoto, S. Expression of neutral endopeptidase activity during clinical and experimental acute lung injury. Respir. Res. 2010, 11, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Bayes-Genis, A.; Morant-Talamante, N.; Lupón, J. Neprilysin and natriuretic peptide regulation in heart failure. Curr. Heart Fail. Rep. 2016, 13, 151–157. [Google Scholar] [CrossRef] [PubMed]
  82. Acanfora, D.; Ciccone, M.M.; Scicchitano, P.; Acanfora, C.; Casucci, G. Neprilysin inhibitor–angiotensin II receptor blocker combination (sacubitril/valsartan): Rationale for adoption in SARS-CoV-2 patients. Eur. Heart-J.–Cardiovasc. Pharmacother. 2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Crowley, S.D.; Rudemiller, N.P. Immunologic effects of the renin-angiotensin system. J. Am. Soc. Nephrol. 2017, 28, 1350–1361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic representation of RAS. In the unperturbed system, soluble proteins that are explicitly considered in the model are in blue grey, the peptides in light blue, and the peptide-bound membrane proteins in medium blue. The activities and enzymes considered only through reaction rates are in green. The feedback loop is indicated in blue. In the perturbed system, the drugs are in orange and SARS-CoV-2 in dark red.
Figure 1. Schematic representation of RAS. In the unperturbed system, soluble proteins that are explicitly considered in the model are in blue grey, the peptides in light blue, and the peptide-bound membrane proteins in medium blue. The activities and enzymes considered only through reaction rates are in green. The feedback loop is indicated in blue. In the perturbed system, the drugs are in orange and SARS-CoV-2 in dark red.
Viruses 12 01367 g001
Figure 2. Dynamical response of RAS to ACE-I (enalapril) administration. Comparison between the computational prediction (green) and the experimental data (brown) of normalized AngI (a) and AngII (b) as a function of time (in hours) after the single dose administration. Continuous lines are obtained through data interpolation. (c) Measured DBP averaged over more than ten ACE-I types as a function of the normalized dosage ξ ACE I (dosage divided by maximal dosage) (brown points) and predicted DBP as a function of ξ ACE I values considering γ ACE I = 0.5 ξ ACE I (continuous green line) and γ ACE I = 0.4 ξ ACE I 1 / 4 (dashed green line). (d) Predicted effect of the combination of ACE-I and ARB on DPB values as a function of the normalized drug dosages ξ ACE I and ξ ARB , considering γ ACE I = 0.5 ξ ACE I and γ ARB = 0.5 ξ ARB .
Figure 2. Dynamical response of RAS to ACE-I (enalapril) administration. Comparison between the computational prediction (green) and the experimental data (brown) of normalized AngI (a) and AngII (b) as a function of time (in hours) after the single dose administration. Continuous lines are obtained through data interpolation. (c) Measured DBP averaged over more than ten ACE-I types as a function of the normalized dosage ξ ACE I (dosage divided by maximal dosage) (brown points) and predicted DBP as a function of ξ ACE I values considering γ ACE I = 0.5 ξ ACE I (continuous green line) and γ ACE I = 0.4 ξ ACE I 1 / 4 (dashed green line). (d) Predicted effect of the combination of ACE-I and ARB on DPB values as a function of the normalized drug dosages ξ ACE I and ξ ARB , considering γ ACE I = 0.5 ξ ACE I and γ ARB = 0.5 ξ ARB .
Viruses 12 01367 g002
Figure 3. Simulated response of RAS to viral infection. (a) The γ CoV function used to model the effect of the infection as a function of C t , the cycle threshold of the virus. (bd) Predictions obtained from our model for the normalized levels of AngII and Ang1-7 and for the physiological PaO2/FiO2 value, as a function of C t , for normotensive (blue) and hypertensive (red) individuals.
Figure 3. Simulated response of RAS to viral infection. (a) The γ CoV function used to model the effect of the infection as a function of C t , the cycle threshold of the virus. (bd) Predictions obtained from our model for the normalized levels of AngII and Ang1-7 and for the physiological PaO2/FiO2 value, as a function of C t , for normotensive (blue) and hypertensive (red) individuals.
Viruses 12 01367 g003
Figure 4. Impact of different RAS-blocking drugs in normotensive (blue) and hypertensive (red) SARS-CoV-2 patients. Predicted P A O 2 / F I O 2 value as a function of the cycle threshold value C t and (a) γ A C E I , (b) γ A R B , and (c) γ D R I functions that model the administration of the corresponding drugs.
Figure 4. Impact of different RAS-blocking drugs in normotensive (blue) and hypertensive (red) SARS-CoV-2 patients. Predicted P A O 2 / F I O 2 value as a function of the cycle threshold value C t and (a) γ A C E I , (b) γ A R B , and (c) γ D R I functions that model the administration of the corresponding drugs.
Viruses 12 01367 g004
Figure 5. Impact on the P A O 2 / F I O 2 value of the administration of rhACE2 and Ang1-7 in normotensive (blue) and hypertensive (red) SARS-CoV-2 patients. (a) Predicted P A O 2 / F I O 2 values as a function of C t and γ rhACE 2 . (b) Predicted P A O 2 / F I O 2 values as a function of C t and η A n g 17 .
Figure 5. Impact on the P A O 2 / F I O 2 value of the administration of rhACE2 and Ang1-7 in normotensive (blue) and hypertensive (red) SARS-CoV-2 patients. (a) Predicted P A O 2 / F I O 2 values as a function of C t and γ rhACE 2 . (b) Predicted P A O 2 / F I O 2 values as a function of C t and η A n g 17 .
Viruses 12 01367 g005
Table 1. Half-lives of the species involved in RAS and other parameters of the model. “Fitted” means fitted on experimental data.
Table 1. Half-lives of the species involved in RAS and other parameters of the model. “Fitted” means fitted on experimental data.
ParameterUnitValuesReference
h a g t min600[35]
h a n g 1 - - 7 min0.5[35]
h a n g I min0.5[35]
h a n g I I min0.5[35]
h a n g I V min0.5[35]
h a t 1 r min12[35]
h a t 2 r min12[35]
h r e min12[35]
h m a s min12-
c r e 1/min20[36,54]
A 0 mmHg450Fitted
A 1 mmHg267Fitted
P 0 mmHg73.6Fitted
P 1 mmHg mL/fmol0.43Fitted
a-0.53Fitted
b-16.7Fitted
Table 2. Equilibrium concentrations of the species involved in RAS and production and reaction rate parameters, for healthy normotensive and hypertensive humans. “Solved” means solved from the model.
Table 2. Equilibrium concentrations of the species involved in RAS and production and reaction rate parameters, for healthy normotensive and hypertensive humans. “Solved” means solved from the model.
ParameterUnitNormotensiveHypertensiveReference
[AGT] 0 fmol/mL6 × 10 5 6 × 10 5 [55]
[AngI] 0 fmol/mL70110[56,57]
[AngII] 0 fmol/mL28156[56,57]
[Ang1-7] 0 fmol/mL3692[56,57,58]
[AngIV] 0 fmol/mL11[59]
[AT1R-AngII] 0 fmol/mL1585[37]
[AT2R-AngII] 0 fmol/mL527[37]
[RE] 0 fmol/mL9.4325.25Solved
[MAS-Ang1-7] 0 fmol/mL6.4315.92Solved
k a g t fmol/(mL min)881.821198.22Solved
β 0 fmol/(mL min)0.542.21Solved
c a c e 1/min1.313.21Solved
c a c e 2 1/min1.800.82Solved
c a n g I V 1/min0.050.01Solved
c a t 1 r 1/min0.030.03Solved
c a t 2 r 1/min0.010.01Solved
Table 3. Comparison between model predictions and experimental values of AngI and AngII levels normalized by their value before the administration of the drugs. Range is the interval of experimental values, and rmsd is the root mean square deviation between experimental and predicted values, computed over all time points; Np is the number of time points. ACE-I, angiotensin-converting enzyme inhibitor; ARB, angiotensin receptor blocker; DRI, direct renin inhibitor.
Table 3. Comparison between model predictions and experimental values of AngI and AngII levels normalized by their value before the administration of the drugs. Range is the interval of experimental values, and rmsd is the root mean square deviation between experimental and predicted values, computed over all time points; Np is the number of time points. ACE-I, angiotensin-converting enzyme inhibitor; ARB, angiotensin receptor blocker; DRI, direct renin inhibitor.
DrugsClassDose[AngI](t)/[AngI] 0 [AngII](t)/[AngII] 0 NpRef.
(mg)rmsd (Range)rmsd (Range)
EnalaprilACE-I201.31 [1.0–9.2]0.09 [0.2–1.0]5[60]
LosartanARB500.61 [1.0–2.1]-3[61]
ValsartanARB8500.83 [1.0–2.2]-3[61]
IrbesartanARB1500.97 [1.0–4.4]-3[61]
AliskirenDRI400.13 [0.4–1.1]0.14 [0.5–1.0]6[62]
AliskirenDRI800.15 [0.4–1.0]0.16 [0.4–1.0]6[62]
AliskirenDRI1600.26 [0.2–1.0]0.20 [0.3–1.0]6[62]
AliskirenDRI6400.29 [0.1–1.0]0.29 [0.1–1.0]6[62]
Mean     0.57       0.18
Table 4. Prediction of biochemical and clinical features of SARS-CoV-2 patients.
Table 4. Prediction of biochemical and clinical features of SARS-CoV-2 patients.
UninfectedMildModerateSevere
C t 40.031.527.623.8
                  Normotensive
[AngII] (fmol/mL)28323638
[Ang1-7] (fmol/mL)362151
PaO2/FiO2 (mmHg)45030014598
DBP (mmHg)80818282
                  Hypertensive
[AngII] (fmol/mL)156186221231
[Ang1-7] (fmol/mL)9255152
PaO2/FiO2 (mmHg)45029211560
DBP (mmHg)110117125128
Table 5. Predicted effects on AngII and Ang1-7 levels, PaO2/FiO2, and DBP upon drug administration to normotensive and hypertensive COVID-19 patients. the drug administrations are modeled by γ A C E I , γ ARB , γ DRI , γ r h A C E 2 = 0.5 , and η A n g 17 = 25 fmol/(mL min) and moderate SARS-CoV-2 infection by γ C o V = 27.6 .
Table 5. Predicted effects on AngII and Ang1-7 levels, PaO2/FiO2, and DBP upon drug administration to normotensive and hypertensive COVID-19 patients. the drug administrations are modeled by γ A C E I , γ ARB , γ DRI , γ r h A C E 2 = 0.5 , and η A n g 17 = 25 fmol/(mL min) and moderate SARS-CoV-2 infection by γ C o V = 27.6 .
DrugsNo DrugsACE-IARBDRIrhACE2Ang1–7
Normotensive—Moderate Infection
[AngII]/[AngII] 0 1.291.101.980.991.101.29
[Ang1-7]/[Ang1-7] 0 0.150.130.230.110.680.64
PaO2/FiO2 (mmHg)1451880216337278
DBP (mmHg)828180808182
Hypertensive—Moderate Infection
[AngII]/[AngII] 0 1.421.121.550.771.141.42
[Ang1-7]/[Ang1-7] 0 0.160.130.180.090.700.36
PaO2/FiO2 (mmHg)11518583268332167
DBP (mmHg)125114101102115125
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pucci, F.; Bogaerts, P.; Rooman, M. Modeling the Molecular Impact of SARS-CoV-2 Infection on the Renin-Angiotensin System. Viruses 2020, 12, 1367. https://0-doi-org.brum.beds.ac.uk/10.3390/v12121367

AMA Style

Pucci F, Bogaerts P, Rooman M. Modeling the Molecular Impact of SARS-CoV-2 Infection on the Renin-Angiotensin System. Viruses. 2020; 12(12):1367. https://0-doi-org.brum.beds.ac.uk/10.3390/v12121367

Chicago/Turabian Style

Pucci, Fabrizio, Philippe Bogaerts, and Marianne Rooman. 2020. "Modeling the Molecular Impact of SARS-CoV-2 Infection on the Renin-Angiotensin System" Viruses 12, no. 12: 1367. https://0-doi-org.brum.beds.ac.uk/10.3390/v12121367

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop