A computational analysis of the long-term regulation of arterial pressure

The asserted dominant role of the kidneys in the chronic regulation of blood pressure and in the etiology of hypertension has been debated since the 1970s. At the center of the theory is the observation that the acute relationships between arterial pressure and urine production—the acute pressure-diuresis and pressure-natriuresis curves—physiologically adapt to perturbations in pressure and/or changes in the rate of salt and volume intake. These adaptations, modulated by various interacting neurohumoral mechanisms, result in chronic relationships between water and salt excretion and pressure that are much steeper than the acute relationships. While the view that renal function is the dominant controller of arterial pressure has been supported by computer models of the cardiovascular system known as the “Guyton-Coleman model”, no unambiguous description of a computer model capturing chronic adaptation of acute renal function in blood pressure control has been presented. Here, such a model is developed with the goals of: 1. representing the relevant mechanisms in an identifiable mathematical model; 2. identifying model parameters using appropriate data; 3. validating model predictions in comparison to data; and 4. probing hypotheses regarding the long-term control of arterial pressure and the etiology of primary hypertension. The developed model reveals: long-term control of arterial blood pressure is primarily through the baroreflex arc and the renin-angiotensin system; and arterial stiffening provides a sufficient explanation for the etiology of primary hypertension associated with ageing. Furthermore, the model provides the first consistent explanation of the physiological response to chronic stimulation of the baroreflex.

Theoretical analysis and observations of the control of blood volume and salt content by the kidneys has led to the hypothesis that arterial pressure is determined in the long-term (over time scales of days or more) by the balance between the level of salt intake and the acute relationship between pressure and salt excretion by the kidneys. In fact, it is stated that the renal pressure/volume control system adjusts arterial pressure with "infinite gain" and that the renal function curve and rate of salt and water intake are the "two sole determinants of the long" 1 . Yet while it is certain that at any steady level of arterial pressure the rates of salt and water intake and excretion are balanced, there is considerable debate over how this balance is achieved, and thus what are the long-term determinants of arterial pressure [2][3][4][5] .
One viewpoint is that the acute relationship between pressure and salt excretion-the acute pressure-natriuresis mechanism-represents a physiological input-output relationship and that alterations to this mechanism underlie most (if not all) chronic changes in pressure 1,6 . While competing ideas identify the nervous system and its influence on the heart and vasculature as the primary long-term controllers of arterial pressure 3,4 , the hypothesis that long-term control of blood pressure is achieved through the renal pressure-diuresis represents the dominant thinking in the field 7 .
The most widely recognized model of long-term blood pressure control is the model developed over several decades by Guyton and colleagues. The 1972 realization of the model 6 invokes approximately 160 variables, and several hundred adjustable parameters 8 . More recent versions of the model involve many thousands of variables and tens of thousands of parameters 9,10 . While several realizations of the Guyton model are or have been disseminated as computer programs, human readable expositions of the various versions of the model are lacking. Moreover, there is no published report of any of these models in which the governing equations are defined, the parameter values are provided, the data that were used for model identification are reported, and the model identification process and results are reported. Indeed, it is unlikely if not impossible that even the relatively simple 1972 version of the model was formally identified. These facts unfortunately make this seminal work largely impenetrable and possibly irreproducible. Most critically, the Guyton-Coleman model is based on the assertion that arterial blood pressure is, under all circumstances, controlled primarily by the kidney 11 . Therefore the Guyton-Coleman model cannot be used to explore alternative hypotheses.
In the Systems Approach for PHysiological Integration of Renal, cardiac and respiratory functions (SAPHIR) project, Thomas et al. 8 aimed to develop and disseminate a core "model of human physiology targeting the short-and long-term regulation of blood pressure, body fluids and homeostasis of the major solutes" including "the main regulatory sensors (baro-and chemoreceptors) and nervous (autonomic control) and hormonal regulators (antidiuretic hormone, aldosterone and angiotensin)". The project goal is to, for the first time, provide an open-source transparent model integrating these systems. In 2008, the SAPHIR group published a core model of circulation and volume exchange across fluid compartments that is derived from the Guyton models, consisting of 20 equations and perhaps less than 100 parameters 8 . This model does not account for adaptation of the acute pressure-diuresis/natriuresis function to chronic changes in arterial pressure or volume/salt loading. Thus, it does not capture the concept of the chronic renal function curve, which is thought to be of fundamental importance to the long-term regulation of arterial pressure. Recently Averina et al. 4 reported a mathematical model of long-term control of arterial pressure that captures the concept of chronic adaptation of pressure-natriuresis and, to our knowledge, represents the only published model incorporating the concept of adaptation of the acute pressure-natriuresis relationship in response to changes in salt/volume loading. Thus, the model of Averina et al. represents the model of record of regulation of blood pressure accounting for the chronic renal function curve. Yet, this model was developed to illustrate that the chronic renal function curve need not represent a physiological input-output mechanism, and to introduce an alternative description of the long-term regulation of arterial pressure and the etiology of salt-and angiotensindependent hypertension.
Similar to the goals of the SAPHIR project, we developed a mathematical model of the Guyton concept of long-term control of arterial pressure, analysed its behavior and compared its predictions to experimental observations. Similarly, an open-source transparent model integrating these systems is provided. The approach pursued here differs from that of Thomas and colleagues in that while the SAPHIR project explicitly centers on translating the components of the Guyton model into an explicitly defined core model, our goal is

Changes from Version 1
We are grateful to the reviewers for their insightful commentaries and constructive criticism. An overriding issue in the reviews is the "phenomenological" nature of the model equations used in this study. Although the phenomenological nature of the model was not judged to be a shortcoming, the issue was raised repeatedly. Yet all models are "phenomenological" at some level. Indeed, even the models pointed out by the reviewers incorporate exactly the same sort of phenomenological representations that we have invoked. Thus, rather than the "phenomenological" nature of our model, the more specific issue may be our model's simplicity. Indeed, certain aspects of our model are simple compared to certain aspects of certain other models. That simplicity could be judged to be a strength or a weakness, depending on the context or application. Thus, it is our opinion the important questions are: Does the model capture the phenomena that we set out to capture? Is the model validated against experimental data that were not used to identify the model? Does the model facilitate any insights into physiological phenomena? In other words, we invoke the often repeated words of George Box, "All models are wrong, but some are useful". So, how is our model useful? The reviews have helped us a great deal to spell out clear answers to that question. In the revised text we haveclarified the goals and stated the conclusions so that they are clearly justified by the analysis presented in the paper. In addition, a basic parameter sensitivity has been performed with results presented in Figure 1. This analysis revealed a slightly improved parameterization of the model, reflected in the updated Table 1 and Figure 6- Figure 10. Finally, we have added a panel to Figure 4 illustrating the overall organization of the model.

See referee reports
REVISED text we have clarified the goals and stated the conclusions so that to develop a model from scratch in which all model components and associated data are formally identified based on experimental data. Furthermore, given the demands of model identification and validation, our goals are necessarily more narrowly constrained than those of the SAPHIR project and of the original Guyton-Coleman models. In our study, a model is constructed in the spirit of the Guyton-Coleman models, adopting a practical approach to datadriven phenomenological representations of physiological systems, rather than aiming for physical and mechanistic realism.
Specifically, we have developed a model of the long-term control of arterial pressure that captures the Guyton concept of pressurediuresis/natriuresis as physiological input-output relationships. The model is designed to meet the following criteria: 1. The model accounts for the effects of the baroreflex on the circulation, the heart, and the kidney, and the effects of the renin-angiotensin system on the circulation, the heart, and the kidney; 2. The model is explicitly documented for definition and reproducibility: all model equations are reported herein and justified based on experimental data; 3. Parameter values are reported and justified: All parameter values are estimated by comparing model simulations to measured data. The identified model is used to probe physiological mechanisms underlying the chronic renal function curve, to provide insight into how chronic stimulation of the baroreflex leads to chronic reductions in arterial pressure, and to generate hypotheses regarding the etiology of primary hypertension.

Model Components 1: Aorta/large-artery mechanics
Based on the simple approximation of a thin-walled cylinder, the strain ε in the aorta is computed as a function of volume V Ao where V 0 is a parameter representing the unstressed volume, and d A /d 0 is the ratio of diameter to unstressed diameter. (Here ε is defined to be equal to 1 in the unstressed state when d A = d 0 ). The aortic pressure-volume relationship is assumed to be governed by where C Ao represents an acute compliance and V sAo (t) accounts for creep mechanics of aortic wall, simulated according to where τ cAo is the time constant of stress relaxation and C Ao /C ∞ is the ratio of acute to effective chronic compliance of the vessel.
where the measured aortic pressure waveform ( Figure 1A) is used to numerically approximate dP A0 (t)/dt in integrating Equation (4). The parameters γ Ao , C Ao , and τ cA , were adjusted to match the data in Figure 1B; the values of V 0 and d 0 were set arbitrarily by assuming a cylindrical vessel of length 30 mm. All parameter values are listed in Table 1.

Model Components 2: Kinetics of baroreflex afferent firing
The baroreceptor afferent firing rate is assumed to be governed by the rate of change of strain in the vessel wall. The model invokes a moving average strain value ε -(t), which is computed The time constant τ s is an adjustable parameter. The baroreceptor firing rate is assumed proportional to δ ε = max(ε-ε -,0) via the satu- where s(t) represents the fraction of baroreceptor afferents that are in an active/permissible state, and δ 0 and f 0 are additional adjustable parameters. Equation (6) is a static nonlinearity 13 that enforces a saturating response. It is assumed that the baroreceptors within the population transition from an active to inactive state at a rate proportional to the firing rate, and transition to the active state at a constant rate: The adjustable parameters in the baroreflex afferent model (τ s , δ 0 , f 0 , a, and b) are identified based on measurements following step changes in non-pulsatile carotid sinuspressure 14 and ramps in in vivo pulsatile aortic pressure 12 .
The data in Figure 2A, obtained from Chapleau et al. 14 , correspond to the multi-fiber response of the carotid sinus baroreceptors following a step change in pressure from 40 mmHg to 93 mmHg. To simulate this experiment, the pressure is assumed to follow a time course described by where pressure is expressed in units of mmHg and time in seconds. Thus it is assumed that pressure is increased by +53 mmHg over an interval from -1 to 0 seconds. Since the overall model does not differentiate between aortic and carotid baroreflex signals, data from the carotid pressure step are matched to the model of Equation (4)-Equation (7) in Figure 2A. Following the step pressure increase the afferent firing rate f BR rapidly increases to a maximum value before decaying over a time scale of several hundred seconds. This decay is determined in the model by the combined action of mechanical relaxation governed by Equation (5) and inactivation governed by Equation (7). Responses to ramps of pressure in vivo are compared to data of Coleridge et al. in Figure 3. Coleridge et al. adjusted in vivo pressure in the aortic sinus by placing hydraulic occlusion at various positions along the aorta. Model simulations are based on applying a constant increase/decrease to measured P Ao (t), resulting in the pressure time series plotted in Figure 3A and 3C. It is assumed that mean pressure is adjusted at a rate of ±10 mmHg sec -1 , corresponding to the experimental measurements. The top panel in the figure corresponds to simulations and data associated with the baseline state, where initial mean pressure is 100 mmHg. Data on afferent firing rate for this experiment are plotted as open circles in Figure 3B. The solid line in Figure 3B represent simulated f BR (t) versus P Ao (t). The red line   . Baroreflex response to in vivo pressure ramps. A. Applied pressure transients for increasing and decreasing pressure starting at a baseline mean pressure of 100 mmHg. Pressure is increased/decreased so that the mean pressure changes at the rate of ±10 mmHg sec -1 . B. Simulated baroreceptor firing response to pressure ramps from A is compared to data from Coleridge et al. 11 . C. Applied pressure transients for increasing and decreasing pressure starting at a baseline mean pressure of 125 mmHg. D. Simulated baroreceptor firing response to pressure ramps from C is compared to data from Coleridge et al. 11 . See text for details on simulation protocol. In B and D, experimental data are plotted as circles; black line represents the simulation prediction; red line represents simulation predictions averaged over each heart beat.
represents f BR (t) and P Ao (t) averaged over each beat. The lower panel in Figure 3 corresponds to a hypertensive state where the aortic pressure was held at an elevated level with mean of 125 mmHg for 20 minutes prior to the pressure ramp experiment. Thus the system has had 20 minutes to adapt or "reset" to a new mean pressure, resulting in a right shift of the baroreflex response f BR versus P Ao .
The values of the parameters τ s , δ 0 , f 0 , a, and b were adjusted to match the data in Figure 2A and Figure 3.

Model Components 3: Mechanics of the heart and circulation
The circulation is modeled as a closed-loop lumped-parameter circuit illustrated in Figure 4A, which ignores the pulmonary circulation and treats the heart as time-varying elastance representing the left ventricle 15 . The left-ventricular pressure is described by where E LV (t) is the left-ventricular elastance and V LV (t) is the volume of blood in the ventricle. The elastance is simulated using a smooth function that increases in systole and decreases as the heart relaxes: where θ ∈ (0,1) is the fraction of a heart beat that has elapsed at a given time. The parameters E max and E min represent minimum effective elastance in the ventricle during systole and diastole. The factor (0.75 + φ SN ) multiplying E max accounts for the effect of sympathetic tone on contractility, where φ SN (t) ∈ (0,1) is a model variable (see below) representing the sympathetic drive. Under baseline conditions φ SN ≈ 0.25 and therefore, under maximal sympathetic stimulation, cardiac contractility is approximately 175% of baseline.
The variable θ is simulated via where H = H 0 + H 1 (φ SN -0.25) is the heart rate, and θ(t) is reset to 0 each time the value reaches 1. The parameters H 0 and H 1 are set to give a baseline heart rate of 75 beats min -1 and a heart rate of 150 beats min -1 under maximal sympathetic tone.
The circuit model is simulated based on equations for the state vari- , and V sV (t). The governing equations for the six volume variables are dV dt where P LV (t) is determined by Equation (9).
The flows Q in and Q urine represent the rates of volume uptake/infusion and urine production, described below.
The max(0,•) terms in Equation (12) account for the valves, which permit flow only in the direction indicated in Figure 4A.
The neurohumoral control mechanisms captured by the model are diagrammed in Figure 4B. The R Ao and R V resistances are set to constant values, while other resistances and capacitances vary with sympathetic tone and angiotensin II level. Specifically, where C A0 , C V0 , and R A0 are constants, and α 1 , α 2 , α 3 , and α 4 are constants that determine the magnitude of the effects of vasoconstriction via sympathetic tone and angiotensin II. The variable φ A2 represents the plasma angiotensin II activity. The governing equations for φ A2 are described below.
The R A resistance is further governed by a whole-body autoregulation phenomenon, as incorporated into the Guyton-Coleman models 6,8,9 . The function r AR (t) ∈ (0,1) accounts for autoregulatory effects on systemic arterial conductivity based on where F -(t) is a variable that averages the mean cardiac output with a moving average defined by a first-order process: The parameters F 0 and F 1 determine the autoregulatory response of resistance to changes in flow. The pressures are computed from the relationships in which venous compliance is simulated using a linear formulation of stress relaxation similar to that used for the aorta. Specifically, venous stress relaxation kinetics are governed by The constant V V01 in Equation (17) represents an unstressed volume for the overall cardiovascular system.
There are a total of 23 parameters associated with this component of the model. Assignment of the values listed in Table 1 was guided by a variety of data sets and previous computational models. The values of E max and E max were obtained by scaling the model of 16 to provide reasonable pressure under baseline conditions at ventricular volumes appropriate for dog. The cardiac cycle timing parameters T M and T R were set to the values used by Beard 16 . Heart rate parameters H 0 and H 1 were determined as described above. The baseline resistance and compliance parameters R out , R Ao , R A0 , R V , C A0 , C V0 , and V V01 were chosen so that under baseline conditions ( Figure 5), the mean pressure is 100 mmHg, the diastolic and systolic pressures are 85 and 115 mmHg, and the ejection fraction is 0.58. (The values for these seven parameters do not represent a unique set that gives these outputs under baseline conditions). The parameters γ V and τ cV were set to match measurement of stress relaxation in the canine jugular 17 .
The remaining seven parameters in the cardiac and circulation model component (α 1 , α 2 , α 3 , α 4 , F 0 , F 1 , and τ AR ) were identified by simulating the responses of the system to volume infusion and hemorrhage, as detailed below in Results.

Model Components 4: Autonomic system
The whole-body sympathetic tone is represented in the model by the variable φ SN (t) ∈ (0,1) and is determined by the baroreflex arc: Thus, in the absence of baroreflex firing, the sympathetic tone will approach the maximum value of 1. The constant parameter f SN is set so that under baseline conditions φ SN (t) = 0.25. Thus, following a rapid severe drop in pressure φ SN (t) will approach a value that is four times the baseline value.

Model Components 5: Renin-angiotensin system
The state of the renin-angiotensin system is captured by two variables: φ R (t) ∈ (0,1) and φ A2 (t) ∈ (0,1), which represents plasma renin and angiotensin II activity governed by the combined action of sympathetic tone and pressure on renin release.
The plasma renin variable φ R (t) is governed by where φ R ∞ decreases with increasing time-averaged arterial pressure P -. The parameters P 1 and P 2 determine the shape of the steady state relationship between renin production and pressure. This formulation assumes that the relationship between steady-state φ R and pressure P is shifted by the sympathetic tone. Based on observations showing that plasma renin activity and angiotensin II levels are nearly perfectly linearly related in vivo 18 , we assumed that φ A2 (t) follows φ R (t) according to These expressions simplify the known (and unknown) mechanisms governing renin production into a simple phenomenological relationship between pressure, sympathetic tone, and the kinetics of the renin-angiotensin system. Here the kinetics of the renin-angiotensin system are treated with much less detail than other models, which account for whole-body pharmacokinetics 19 and saturable enzyme kinetics 20 . The five parameters invoked in this model component (τ R , τ A2 , g, P 1 , and P 2 ) are identified based on comparing simulations to measurements of pressure, heart rate, and plasma renin activity in rabbit during graded hemorrhage, as detailed below in Results. (The simulation parameter τ P was arbitrarily set at 15 seconds).

Model Components 6: Neurohumoral control of pressurediuresis/natriuresis
Regulation of body-fluid volume is assumed governed by a linear pressure-diuresis relationship: where Q urine is rate of volume output via the kidneys and k 1 is the slope of the relationship between Q urine and pressure. The model variable P s (t) is the variable offset of the pressure-diuresis relationship, which is controlled by sympathetic tone and angiotensin II level: The constant parameters P s,min , P s,max , φ 0 , and φ 1 determine how the acute pressure-diuresis relationship shifts in response to changes in the tone variables φ SN and φ A2 . Thus it is assumed that φ SN and φ A2 have equal and additive effects on renal function. Equation (22) assumes that changes in P s governed by sympathetic tone and the renin-angiotensin system occur with a time constant τ k . Furthermore, the maximal rate of urine production is set to 10 ml min -1 .
The parameters k 1 , P s,max , P s,min , φ 0 , and φ 1 , were set by matching model predictions to data on the pressure diuresis relationship under physiological conditions, with angiotensin II infusion, and with administration of an angiotensin converting enzyme (ACE) inhibitor. To compare predictions of Equation (21) and Equation (22) to renal output measurements under ACE inhibition, we assumed that for the chronic measurements, sympathetic tone was maintained at its baseline value and φ A2 = 0. Fitting data from Hall et al. 21 , it is estimated that P s ∞ (φ SN = 0.25, φ A2 = 0) = 75 mmHg. Similarly, with angiotensin II infusion, we assumed that sympathetic tone was maintained at its baseline value and φ A2 = 1. The data of Hall et al. yield an estimate of P s ∞ (φ SN = 0.25, φ A2 = 1.0) = 125 mmHg for these conditions. Finally, under baseline conditions, with P -= 100 mmHg, urine output was estimated to be Q urine = 0.5835 ml min -1 , based on volume infusion experiments described in the Results. With k 1 = 0.125 ml sec -1 mmHg -1 (see Results), and using the baseline φ A2 value of 0.1864 (which is determined from blood withdrawal experiments; see Results), it follows that P s ∞ (φ SN = 0.25, φ A2 = 0.1864) = 95.3 mmHg. These estimates of P s ∞ at three different values of φ SN + φ A2 were used to estimate the values of P s,max , P s,min , φ 0 , and φ 1 .
The time constant τ k was determined from data on the rate of urine production following infusion of blood leading to an acute increase in pressure and drop in φ SN and φ A2 . Analysis of data from these experiments is detailed below in the Results together with a direct comparison between the data of Hall et al. 21 and model simulations.

Results 1: Response to volume infusion
Guyton and colleagues conducted experiments in which a large amount of blood was infused into an anesthetized dog, resulting in a rapid increase in total blood volume of 45% compared to initial baseline value. The model of Equation (1)-Equation (22) was simulated and compared to data from this experiment to identify adjustable model parameters and to probe the predicted response of unmeasured model variables to this protocol. Figure 6 shows the predicted effects of infusing 45% of initial baseline blood volume in a normal animal on arterial pressure, cardiac output, rate of urine formation, blood volume, sympathetic tone, and renin-angiotensin level together with experimental data obtained from Dobbs et al. 22  Since all other data sets analyzed here have a baseline pressure of 100 mmHg, 15 mmHg was subtracted from the data from these experiments to achieve consistency). The initial steady state of the model was obtained based on a constant infusion of Q in = 0.5835 ml min -1 , which matches the mean reported rate of volume output for the baseline data. For times 0 < t < 5 min, Q in was increased to 0.5835 + 126.82 ml min -1 , to result in a total excess volume of 634 ml. After the five-minute infusion, Q in was returned to the baseline value of 0.5835 ml min -1 . The model is able to effectively reproduce the trends in the experimentally measured variables, and predicts that φ SN and φ R drop to less than one half of their resting values following the volume infusion in response to the transient increase in pressure. These reductions in φ SN and φ R confer increases in vascular compliances helping to accommodate the substantial volume increase. As volume is removed from the system φ SN and φ R return toward the baseline values, with φ SN responding substantially faster than φ R .
The data plotted in Figure 6, along with data observations on graded hemorrhage, were used in identifying several adjustable parameters in the model. (See following section for details).

Results 2: Response to hemorrhage
Data from Quail et al. 24 on heart rate, pressure, and plasma renin activity following graded blood withdrawal in rabbits were used to identify parameters associated with the cardiovascular mechanics and the renin-angiotensin system. Figure 7 shows data measured in Figure 6. Simulation of arterial pressure, cardiac output, rate of urine formation, blood volume, φ SN , and φ R following infusion of 45% of initial baseline blood volume in a normal animal. Data on arterial pressure, cardiac output, rate of urine formation, and blood volume are obtained from Dobbs et al. 18 , Guyton et al. 1 , and Prather et al. 19 . The initial steady state of the model was obtained based on a constant infusion of Q in = 0.5835 ml min -1 . For times 0 < t < 5 min, Q in was set to 0.5835 + 126.82 ml min -1 , to result in a total excess volume of 634 ml. The initial condition for the simulation was obtained by setting Q in = 0.5835 ml min -1 and running the model to obtain the steady state.
normal rabbits, where the heart rate measurements have been scaled to a baseline value of 75 beats min -1 , renin activity is scaled to maximum value of 1, and mean pressure is scaled to an initial baseline value of 100 mmHg. In the experiments, blood was withdrawn from the animals at a rate of 2% blood volume per minute. In model simulations, blood is withdrawn at the experimental rate starting at time 0 and with withdrawal stopped at time 17.5 minutes. For this experiment, the initial condition was identical to that used for the volume infusion experiments of Figure 6, and volume infusion and urine output were set to zero for the time course.
Model predictions of H, P -, and φ R compare favorably to experimental data, with mean arterial pressure dropping from the initial value of 100 mmHg to approximately 93 mmHg at the end of the blood withdrawal. The observed increase in heart rate is associated with an increase in φ SN from baseline level of 0.25 to 0.81, and an increase of φ R from baseline level of 0.186 to 0.668, at the termination of the withdrawal period. During this period, arterial blood pressure is protected from a more severe reduction in large part by an increase in vascular tone mediated by increases in sympathetic tone and angiotensin II level. Simulated cardiac output (not shown) also drops during blood withdrawal, reaching a minimum value of 690 ml min -1 (approximately 46% of baseline flow) at 17.5 minutes.
The model is not able to capture the sharp, almost step-like, observed increase in renin that occurs following the 10-minute time point (corresponding to 20% of blood withdrawn). Specifically, the model predicts a more graded response, without the delay observed in the data. A better match to the data might be achieved by incorporating a delay into the governing equations, or by simulating the renin-angiotensin system using a higher-order system of differential equations than Equation (19) and Equation (20). The model predicts that φ SN and H peak at the end of the withdrawal period and partially recover toward baseline levels within a few minutes after the end of blood withdrawal. The renin-angiotensin tone is predicted to remain at an elevated level after the withdrawal. The partial drops in φ SN and H are associated with a graded decrease in pressure continuing over the entire simulated time course, with pressure reaching approximately 85 mmHg at 30 minutes. Pressure continues to drop after the end of the withdrawal period due in part to the action of the whole-body autoregulatory mechanism governed by Equation (16) and Equation (17). Via this mechanism, the drop in flow associated with hemorrhage causes a reduction in arterial resistance, resulting in a recovery in flow and a continuing drop in pressure.
The model predictions in Figure 6 and Figure 7 are most sensitive to the values of 13 parameters from the model components of circulatory mechanics, the renin-angiotensin system, and the kinetics of pressure-diuresis: α 1 , α 2 , α 3 , α 4 , F 0 , F 1 , τ AR , τ R , τ A2 , g, P 1 , P 2 , and t K . The values of these parameters were adjusted to match the data plotted in Figure 6 and Figure 7.  Figure 6. Data from heart rate, pressure, and plasma renin activity following graded blood withdrawal are from Quail et al. 18 with pressure and heart rate scaled as described in the text.

Results 3: Chronic renal function curves
The steady-state relationship between arterial pressure and rate of urine output is obtained from model simulations by varying Q in , the rate of volume infusion, and attaining steady-state model predictions of steady state pressure where Q in = Q urine . Figure 8A plots the predicted relationship between rate of urine output and mean arterial pressure under three different conditions: (1) the normal physiological state; (2) φ R clamped at 0, representing complete block of angiotensin converting enzyme (ACE); and (3) φ R clamped at 1, representing infusion of saturating levels of angiotensin II. Note that the pressure-diuresis relationship is plotted in the traditional manner with arterial pressure on the abscissa and rate of volume excretion on the ordinate, even though in these simulations Q in is varied and steady-state P is computed as a function of Q in .
Model predictions are compared to data from Hall et al. 21 on steady-state pressure at different levels of salt/volume loading in normal dogs, dogs infused with ACE inhibitor, and dogs infused with angiotensin II. Model simulations effectively match experimental data, validating the assignment of parameter values described above under "Neurohumoral Control of Pressure-Diuresis/Natriuresis". Figure 8B illustrates model-predicted steady-state renin and angiotensin II activities (φ R and φ A2 ) and sympathetic tone (φ SN ) as functions of Q in for the normal case. (Recall that from Equation (20) angiotensin II activity φ A2 is equal to φ R in the steady state). As Q in is increased φ R and φ A2 decrease to maintain pressure at nearly a constant level of the simulated range of volume/salt loading. The sympathetic tone on the other hand remains nearly constant around the baseline level of 0.25. Sympathetic tone remains constant because in the model the only determinant of φ SN is the baroreceptor afferent firing rate, which effectively adapts to the small changes in pressure that are associated with the simulated range of Q in . (Thus, the model does not capture suppression of sympathetic tone typically observed with chronic salt/volume loading).

Results 4: Response to chronic baroreflex stimulation
Numerous studies have demonstrated that chronic stimulation of the carotid baroreceptor afferent nerve with implantable devices results in a prolonged decrease in mean arterial pressure 25,26 . Figure 9 illustrates a representative data set from Lohmeier et al. 27,28 for which baroreflex stimulation was applied in dogs for a seven day period, followed by several days of recovery. During the stimulation period, mean arterial pressure drops by approximately 20%, renin activity by approximately 35%, and sympathetic tone (experimentally assayed indirectly by plasma norepinephrine) drops by approximately 55%.
To simulate this experiment, a constant was added to the modelpredicted baroreceptor firing rate in the equation for sympathetic tone. Thus Equation (18) was replaced by where f stim is a parameter adjusted to match the experimental data.
Using f stim = 6.5 sec -1 yields predictions that effectively match the observations of Lohmeier et al. 27,28 (Figure 9). Indeed, given that only one parameter was adjusted, the comparison of model predictions to the five variables plotted in the figure represents a strong validation test for the model. In particular, the observed reduction in renin activity with baroreflex stimulation represents a phenomenon that has not been captured by previous modeling efforts 26,28 . Specifically, the model of Illiescu and Lohmeier predicts a substantial increase in renin production following initiation of stimulation followed by a return to the initial baseline level. In Illiescu-Lohmeier model, the transient increase in renin is a result of the drop in pressure causing an initial increase in renin production. The predicted initial increase is later offset by decreased sympathetic nerve activity 24 . In our model the observed phenomenon of reduction in plasma renin activity emerges as a property of the integrated computational model: stimulation of the baroreflex and associated drop in sympathetic tone is predicted to result in a drop in renin production even with the long-term drop in mean pressure.
Thus our model prediction differs from that of the Illiescu-Lohmeier model in two critical ways. First, our model captures the observed long-term response of renin while the Illiescu-Lohmeier model does not. Second, our model predicts no increase in renin during the first day of baroreflex stimulation, while the Illiescu-Lohmeier model does. Illiescu and Lohmeier predict that an increase in renin "was not measured experimentally because initial blood sampling was not made until 24 h after initiating baroreflex activation" 24 . In contrast, our model predicts that the increase in renin was not measured experimentally because it did not happen.
Results 5: Exploration of the etiology of primary hypertension Pettersen et al. 29 demonstrate that impairment of the baroreflex caused by stiffening of arterial vessels represents a viable hypothesis for the etiology of primary hypertension. Specifically, using a closed-loop cardiovascular mechanics model coupled to a detailed model of the baroreflex arc 30 , Pettersen et al. 29 show that when the stiffness of the aorta is increased to represent changes in mechanical properties associated with ageing, the model based on their mechanogenic hypothesis predicts a substantial increase in mean arterial pressure with age. By assuming a constant blood volume for all age groups, the model explicitly did not account for the regulation of plasma volume and salt through the kidney and the renin-angiotensin system following from the hypothesized shift in the renal pressure-diuresis/natriuresis function curve. The rationale for not including adaptive mechanisms likely to partially ameliorate the effects of arterial stiffening on blood pressure was to test the explanatory sufficiency of the mechanogenic hypothesis by showing that it would predict a stronger relation between blood pressure increase and age than empirically observed. As a first step towards a complete merge of physiological renal function with the model developed by Pettersen et al. to expand the prediction space of the mechanogenic hypothesis, we studied how the current model responded to changes in the aortic compliance.
Using the current model, which does account for blood volume regulation by the kidneys, the hypothesis of Pettersen et al. may be further analyzed by determining how simulations of the current model respond to changes in the aortic compliance. Figure 10A (solid line) plots predicted steady-state mean arterial pressure as a function of relative aortic stiffness, C Ao 0 /C Ao , where C Ao 0 is the baseline normal value of aortic compliance, and C Ao is the value used to obtain the pressures reported in the Figure. The maximum simulated relative stiffness, C Ao 0 /C Ao = 4, is approximately the average relative stiffness for the 75-year-old population simulated in Pettersen et al. 29 . This 4-fold increase in aortic stiffness increases the predicted mean pressure from 100 to 128 mmHg, with systolic/diastolic ratio increasing from 115/90 to 139/117. Furthermore, over the simulated range of stiffness model-predicted φ SN increases from 0.25 to 0.45 and φ R increases from 0.186 to 0.59 ( Figure 10B) with a heart rate increase of approximately 25% as mean pressure increases from 100 to 128 mmHg.
As expected, this predicted pressure increase is less than that predicted for a model that does not account for blood volume regulation. (The model of Pettersen et al. predicts a mean arterial pressure 150 mmHg for 75-year-old group). The difference between the predictions of the current model and that of Pettersen et al. is explained primarily by blood volume regulation: with a 4-fold increase in aortic stiffness the current model predicts the active blood volume to decrease from 1364 to 1050 ml ( Figure 10C). This result supports the view that normal functioning kidneys are able to partially compensate for age-related increases in arterial stiffness 29 .
Therefore this model demonstrates that it is possible to capture the age-related phenotype of primary hypertension as emerging solely from changes to the mechanical properties of the aorta and carotid arteries. This may seem surprising since the simulation results plotted in Figure 10 are obtained with mechanical properties of all components of the vasculature other than large vessels associated with the baroreflex remaining normal. Indeed, it is expected that vascular stiffening would directly impact other physiological processes. For example, Beard and Mescam demonstrated that stiffening of renal arteries can explain the shift in the acute pressure-diuresis curve associated with changes in mean pressure in Dahl S rats 31 . For the model results illustrated in Figure 10 (representing normal kidney function), the predicted increase in mean pressure associated with increased vascular stiffness is partly due to an increase in the pressure that is required to elicit a given baroreceptor firing rate. As a result, the relationship between baroreceptor firing and mean arterial pressure is shifted to higher pressure in hypertension, as has been observed clinically and in animal models 7 . This phenomenon has been called "adaptation" or "resetting". Yet, even though model predictions are in agreement with the observed adaptation of the baroreflex, model predictions are not in agreement with the interpretation that, since baroreceptors adapt, "they cannot participate in the long-term control of arterial pressure" 7 . By contrast, the observed adaptation is a crucial component of the mechanogenic mechanism simulated here. Because stiffening causes a long-term adaptation of the relationship between pressure and baroreceptor firing rate, vascular stiffening can cause a chronic increase in stable mean pressure. In other words, since baroreceptors adapt, they can contribute to the etiology of hypertension.
These model predictions are in agreement with experimental observations of Thrasher, who has shown that mechanically suppressing baroreceptor sensitivity leads to sustained elevation in arterial blood pressure in a canine model 32,33 . Our model predictions are consistent with Thrasher's conclusion that mechanical unloading of baroreceptors leads to sustained increases in arterial pressure 34 .

Parameter sensitivity
Sensitivity values of parameter estimates were estimated by computing the relative change in the mean-square error difference between data and model simulation associated with a change in the parameter value. For example, as explained above, parameters τ s , δ 0 , f 0 , a, and b were estimated by matching model predictions to data in Figure 2A and Figure 3. If E is the minimum normalized mean-square error difference between model prediction and data, then the sensitivity associated with parameter τ s is estimated as |E(τ s + 0.1τ s ) -E(τ s ) | /0.1. Thus a sensitivity value of 1 for a given parameter would mean that a 10% change in the parameter is associated with a 10% in the error. Sensitivities for all parameters identified based on the data in Figure Table 1. Computed sensitivity values range from 0.05 (for α 2 ) to greater than 100 (for C Ao ).
The low sensitivity of results to α 2 can be understood by examining the expression of R A (t) in Equation (13). Since the estimated value of α 2 is greater than 10, the expression is effectively approximated by R A (t) = R A0 α 2 φ SN (1 + α 4 φ A2 ), where the constants R A0 and α 2 cannot be separately identified. Since the baseline resistance and compliance parameters (including R A0 ) are adjusted to (for given values of α 1 , α 2 , α 3 , and α 4 ) are adjusted so that the model generates the reported baseline pulse pressure, cardiac output, and ejection fraction, the value of any of α 1 , α 2 , α 3 , and/or α 4 becomes arbitrary for values significantly greater than 10.

Summary and conclusions
A computational model (Source code is permanently available on: http://zenodo.org/record/7126 (10.5281/zenodo.7126)) was developed and identified to serve as a phenomenological representation of the major physiological processes controlling arterial blood pressure. The design of the model was based on balancing the compromise between complexity/physiological fidelity and identifiability, yielding a model of 16 state variables and 46 parameters. Of the 46 parameters, 23 were set to previously established values or arbitrary values as explained above, and 23 were adjusted to match experimental data illustrated in Figure 1-Figure 3, Figure 6-Figure 8. The model captures physiological phenomena occurring on times scales ranging from milliseconds (e.g. response of baroreceptor firing to arterial pressure) to days (e.g. response of mean pressure to chronic stimulation of baroreflex). It also effectively captures the response of the cardiovascular system to drastic perturbations, including infusion of 45% of blood volume, withdrawal of 35% of blood volume, electrical stimulation of the baroreflex, and the partial compensatory response of the kidneys to blood pressure increase accompanying arterial stiffening.
The major conclusions from this work are: 1. The observed adapting/resetting of the baroreflex to long-term changes in pressure does not indicate that changes in baroreflex function cannot play a role in long-term changes in arterial pressure. Model simulation reveals how the baroreflex arc and the reninangiotensin system may both contribute in a coordinated manner to affect long terms changes in arterial pressure.
2. It is demonstrated how renal function is not the central or primary determinant of long-term arterial pressure in a computer model that captures the observed behavior of the renin-angiotensin system, and chronic adaptation of the acute pressure-diuresis/natriuresis phenomena. This result contradicts the classical 'renocentric' view of the long-term control of blood pressure even though the model is built on the core assumption of the renocentric view that the pressure-diuresis relationship acts as a physiological input-output relationship in which pressure determines renal output. As a result, it is shown that hypertension does not necessarily indicate any degree of renal dysfunction.
3. The physiological response to chronic stimulation of the baroreflex can be explained by a mechanism in which arterial blood pressure is regulated primarily through the baroreflex arc and the interoperation of the baroreflex arc and the renin-angiotensin system. 4. Model simulations are consistent with the mechanogenic hypothesis 29 that arterial stiffening represents a contributing factor causing changes in pressure, sympathetic tone, and renal function associated with ageing. Through pressure-diuresis, normal kidney function is predicted to be able to partially ameliorate the effects of arterial stiffening on blood pressure.
Regarding the forth conclusion, model simulations predict that stiffening of the large arterial vessels that are associated with baroreceptors can contribute substantially to long-term changes in pressure. This result contradicts the 'renocentric' theory of blood pressure control even though the core assumption of the Guyton 'renocentric' theory of blood pressure control (that the pressure-diuresis relationship acts as a physiological input-output relationship in which pressure determines renal output) is a core assumption of this model. Although it has been argued that because its effective gain resets to a given mean arterial pressure the baroreflex cannot play an important role in long-term pressure regulation, this analysis shows how a dynamic system representing the baroreflex and capturing the observed resetting phenomenon can play a key role in the long-term control of pressure and the etiology of hypertension. Model simulations reveal that vascular stiffening increases sympathetic tone and shifts the effective baroreflex response to an increased pressure baseline. Because these changes cause an increase in angiotensin II activity and associated shift in the acute pressurediuresis relationship, the renocentric 7 and mechanogenic 29 hypotheses for the etiology of primary hypertension could potentially be interpreted as compatible. However, as changes to renal function arise as downstream consequences of the mechanical remodeling, the mechanogenic explanation subsumes more of the biology involved in the etiology of hypertension compared with a renocentric explanation. This is a major conceptual advance that also provides a new interpretational framework for available experimental and clinical data.
The model developed here represents a simplification and several components are represented largely in phenomenological terms. A more detailed model may be constructed through integration of more biophysically based and detailed component models 30,31,35 into the framework developed here. Regardless, the simplifications invoked in the model do not impact the major conclusions. Perhaps as important as the simplified modeling framework is the fact that a model identified from data obtained from experiments on dogs and rabbits, is used to draw insight into human physiology and pathophysiology. The implicit assumption here is that arterial pressure regulation processes are at least qualitatively similar in dog, rabbit and man.

F1000Research
No competing interests were disclosed.

Strong points: Sound "horizontal" (multi-organ) physiological modelling & simulation (following a systems/engineering analysis and conception). Modular, integrative, multi-dynamic approach.
Targeted, specific hypothesis -proposes one possible/sufficient etiogenic mechanism in the context of age-related hypertension.
Weak points: Phenomenological equations. Some conclusions drawn are too assertive (especially with respect to such a complex, multi-factorial, and evolving, pathology: hypertension); the likeliness of the conclusion appears strongly dependent on the limited physiology (and pathology) that is implemented.

GENERAL
In bio-medical research, it has become understood that mathematical modelling and simulation (MS) can provide a complementary approach to experiments and are a necessary response to extreme biological complexity. Indeed, biochemical and genetic circuits modelling are becoming mature and productive. On the other hand, there remains a very serious "modelling gap" to be filled between the ("low-level") biological/physiological processes and ("higher-level") pathophysiology.
Simply put, the "CV-circulatory-pressive" construct by Guyton  not all modelling endeavours refer to it in one way or another. This is because it was an exhaustive attempt to describe the CV circulatory system in its physiological entirety, including high-level interactions such as an integrated CV and renal system. However, as mentioned by the Authors, for various reasons it is difficult or impossible to rationally understand and exploit those models (e.g. "too large" in the case of HumMod, or "not identified" as in the case of Guyton's models of 1972 or 1992).

Beard
propose here a "top-down", modular, "engineering" high-quality MS work, as a rationally et al. exploitable (numerically identifiable) alternative to Guyton's model, in order to explore hypertension etiogenic hypothesis related to arterial stiffness, baroreflex and the renin-angiotensin system (RAS). The main incentive relates to the controversed Guytonian view that long-term regulation of arterial pressure is driven by the kidney. In that context the Authors try to demonstrate that arterial stiffening (as in ageing) and/or the central nervous system (here simplified as the "baroreflex") could provide a sufficient explanation(s) of hypertension, through interactions with the kidney pressure-natriuresis function and the RAS.
The significance and relevance of this work for CV regulations and hypertension is obvious, although (i) it is limited by the model itself (i.e. effectively implemented features, simplifications, phenomenological approach for equations, etc.) and (ii) some of the conclusions/assertions should be tempered (among other reasons: given the multifactorial/multi-dynamic nature of HT, that it has been resistant to investigation for more than a century now, and given the simplicity of the proposed construct).

MAJOR COMMENTS
In the targeted domain of the pathophysiology of blood pressure regulation and primary hypertension, this work proposes a rigorously constructed and parameterized mathematical model of the "long-term control (but see below). The authors should emphasize however that the model numerical of arterial pressure" instantiation centers on dog (or rabbit). Also of important note (this should be addressed in the revised version), no mention is made of the parameterization procedure and tools or sensitivity analysis.
Moreover, the parameterization is done under the (quite frequent) implicit hypothesis that arterial pressure regulation and other CV-related processes only differ quantitatively (e.g. in parameter values) and are otherwise qualitatively similar e.g. in dog, rabbit and man. This very point is rarely, if ever, addressed. We do feel that in the present turmoil of "physiological" models targeting *human* pathology, this should be, at minimum, mentioned and preferably justified. Given that the Authors target hypertension, and human especially given the phenomenological approach chosen, why did they not directly parameterize their model as a "human" CVR model?
As a matter of fact, the anatomy and mechanical/physiological/biochemical properties and functions of the human CV system, arteries, heart, kidney and RAS are, to the best of our knowledge, now understood and numerically characterized enough to allow for such a human quantitative approach to be adopted (especially with this articles chosen phenomenological "input-output" approach). This would have given some more resonance and relevance to the work. In particular, even if all CV-pressive, renal and hormonal processes proved to be the same between all these species, the possibility would remain that their relative contributions (and dynamics) to arterial pressure (dys-) regulations are different.
As their objectives, the Authors aimed to: develop a "computer model capturing chronic adaptation of acute renal function in blood pressure as a more rigorous alternative to Guyton These are very strong and ambitious objectives, and possibly flawed by the phenomenological ab initio (but otherwise legitimate) approach chosen.
Further, that this is "a model capturing chronic adaptation of acute renal function in blood pressure" cannot be claimed with such a simple kidney function, for example: the pressure-natriuresis ignoring actions of the ANS (autonomic nervous system) upon renal hemodynamics and tubular function and the angiotensin action on glomerular and tubular function (for instance Glomerular Tubulo Balance and Na/water reabsorption). Kidney nerve are also ignored. afferents What is really addressed in this work are the interactions (possibly etiogenic) between the physiological (input/output) modules that are implemented.
Finally, after more than a century of physiological and clinical research, the etiology of primary hypertension remains unresolved; although it is established that (numerous), genetic, environmental, age-related and adaptive factors do intervene. Moreover, when in the realm of age-related hypertension the Authors do select putative etiogenic processes (arterial stiffening, reduced BR receptors sensitivity) and explore potential interactions with RAS and kidney function (phenomenological pressure-natriuresis) in a simplified CV construct. Driven by the objective (age related HT/arterial stiffening), the proposed model features six central physiological functionalities, as modules (model components) implementing: MC1. "Aorta-large artery mechanics" MC2. "Kinetics of baroreflex afferent firing" MC3. "Mechanics of the heart and circulation" MC4. "Autonomic system" MC5. "Neuro humoral control of pressure diuresis/natriuresis" MC6. "Renin-angiotensin system" For each one of these model components (physiological organ/process), the Authors should give, where possible, at least one published "first-principle" (physiologically-based) model as a reference, e.g. amongst others (i) for the kidney/pressure-natriuresis under hormonal influence, Moss and Thomas, 2013 (ii) and/or for the RAS. Secondly, the model is validated (long-term = days to weeks) using published steady-state pressure-natriuresis kidney function data. Although the Authors target long-term processes, 2 weeks (maximal simulated time here) is far too insufficient to refer to the CV, renal and central physiological mechanisms truly involved (months-to-years in dogs, years and decades in man).

OTHER COMMENTS
Is one of the input arrows missing in the model figure (fig 4)? A model diagram (rather than the electrical analog presented), as well as for the main hypothesis (vs alternative(s)) would be welcome.
To what extent does this work differ from the previous work by the same group (cited, see

F1000Research
To what extent does this work differ from the previous work by the same group (cited, see )? Pettersen 2013 et al., The Matlab(C) model code and scripts proposed ( ) http://zenodo.org/record/7126#.UpOLNOJ21pg have been tested and verified by Dr F. Guillaud (ModTeam/Inserm U1082). The scripts could be run, the model could be simulated and results that were visually identical to those in the paper were obtained.
We have read this submission. We believe that we have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however we have significant reservations, as outlined above.
No competing interests were disclosed. The reviewers have provided an extensive and insightful commentary on our paper. We agree with most of the points raised. We have extracted specific criticisms and responded to them below.

Weak points: Phenomenological equations;
F1000Research models are built based on data from multiple sources, multiple species, different conditions (in vitro, in vivo, etc.) In the end, all of our models are chimeras. We have taken the liberty of adapting the reviewers' language at the end of the discussion section of the revised paper. (We do point out that our implicit assumption is the central assumption of all biomedical research using animal models. It might not be a perfect assumption, but it does give us a starting point.) Also of important note... no mention is made of the parameterization procedure and tools or sensitivity analysis.
Please see our detailed response to this point in the response to the first review.
For each one of these model components (physiological organ/process), the Authors should give, where possible, at least one published "first-principle" (physiologically-based) model as a reference This is an excellent idea in principle. However, it is not clear if it is possible. For example, we are unclear as to how the Moss and Thomas paper could be characterized as representing the hormonal influence on pressure-natriuresis from first principles. Specifically, the Moss and Thomas paper represents the action of angiotensin-II "via the parameter FrANG, which adjusts the " The parameter FrANG appears in the model response of the pressure-natriuresis mechanism. equations as a factor multiplying a phenomenological expression for the fractional reabsorption, and the effect of the hormone is simulated by varying FrANG between "a maximal value of 1 ". Thus, the model is (minimal angiotensin II) and a minimal value of 1/6 (maximal angiotensin II) every bit as "phenomenological" as the one presented here. (This point is not to be taken as a criticism of the excellent study of Moss and Thomas.) The papers by Guillaud and Hannaert, and Claassen et al. do account for the RAS system with more physiological fidelity/detail than our model. We do disagree with the assertion that any of the models discussed here are built from "first" principles. All of these papers invoke the same sort of phenomenological expressions for certain components.
Although the Authors target long-term processes, 2 weeks (maximal simulated time here) is far too insufficient to refer to the CV, renal and central physiological mechanisms truly involved...
The reviewers may be correct that our model is far too insufficient to provide insight into a pathophysiological mechanism that emerges over decades. Yet, we refer the reviewers (and the readers) to the stated conclusions in the paper. figure (fig 4)?

Is one of the input arrows missing in the model
The net volume output rate is Q_urine -Q_in. The figure has been revised to clarify this.
A model diagram (rather than the electrical analog presented), as well as for the main hypothesis (vs alternative(s)) would be welcome.
This has been added. (See Figure 4B.) The main hypothesis here is that pressure-diuresis acts as a physiological input-out relationship, and that renal function is controlled through the RAS and the autonomic system. Thus, fundamentally, the model is constructed to capture the Guyton  Thank you for the attention to our paper and the constructive criticisms.
Most model equations are in some sense ad-hoc choices (or phenomenological, as the authors mention at the very end). These equations are not based on first principles and hence they do not represent the underlying mechanisms... In addition, the authors make a lot of self-references...
We agree with the reviewer's assessment. To reiterate a statement from the paper summary, the model represents a balance complexity/physiological fidelity and identifiability. The sensitivity analysis added to the revised paper indirectly addresses the issue of the phenomenological nature of model components. Finally, we have revised the text to clearly state the nature of the model earlier than the last paragraph of the paper. See third paragraph on page two, section on "Model Components 5", and first paragraph of the summary.
Also, we point out that while we have cited some previous modelling works, the majority of components used here have been constructed , with the overall goal in mind. The varying de novo elastance heart model is adapted from the cited reference. And of course there is nothing new about lumped-parameter models of the circulation. But most of the other (phenomenological) formulae are constructed to match the observed dynamics of responses of these systems to the experiments used to identify the model.

Formula (3): Should the fraction be / instead of /? C Ao C Ao
No, it is correct in the paper. One way to see this is to think about the limit where goes to C Ao infinity. In this case the vessel becomes infinitely compliant if we wait long enough and thus does not support any finite pressure.

Formula (4): An index o is missing on V sA
F1000Research

Formula (4): An index o is missing on V sA
Fixed.
On page 4, column 2, line 2: should be d 0

Fixed
On page 4, column 2, line 5: ' '… the population transition from an active state at a rate … Please state reference (or state that this is a hypothetical assumption)!
The assumption of Equation (7) is stated: "It is assumed that the baroreceptors within the ". We do population transition from an active to inactive state at a rate proportional to the firing rate not know what factors are responsible for the slow "resetting" of the baroreflex response. But, it is compelling that this simple equation is able to capture the nonlinear phenomena of This paper by Beard makes three major claims: 1) to provide a mathematical description for Guyton's et al existing model; 2) to show how stimulation of the baroreflex arc can explain long-term changes in blood pressure via the renal function curve employed in Guyton's model; and 3) to show how arterial stiffening can lead to hypertension without renal injury. I commend the authors for making a much needed mathematical contribution to the field of computational physiology, where one of the most widely recognized and used models of long-term control of arterial pressure has only very rudimentary mathematical description and analysis published. Moreover, the authors are bridging the gap between the renocentric view of hypertension etiology and the emerging experimental data that points to an independent long-term role of neural regulation. For these reasons the paper is a valuable addition to the existing body of knowledge.
There are several issues however, that need to be addressed in a revised version of the paper. The overarching recommendation is to correct the paper's focus and claims to make them more consistent with the model design, as well as to make the model's description easier to understand. In particular, I believe the model should not be presented as a generic model with the "representation of the major but instead be positioned as a conceptual physiological processes controlling arterial blood pressure" framework which aims to illustrate how the baroreflex arc be involved in long-term blood pressure may control. Phenomenological mathematical representation, a small amount of detail, limited analysis, and quality of shown simulations are appropriate for a conceptual model but are insufficient for a model from which one could be extracting conclusions on the particular physiological mechanisms to be implicated in hypertension development.

2.
3. which one could be extracting conclusions on the particular physiological mechanisms to be implicated in hypertension development.

Suggested revisions:
There are several versions of Guyton's model ranging from very simplistic to very complex . Some 1 models are extremely detailed but have a closed architecture (e.g. ). The model created by hummod.org the authors is not sophisticated enough to correspond to most of these versions. Please state clearly which version of Guyton's model is being considered here. In my opinion, the model proposed here seems of approximately the same level of complexity to the one presented by Guyton and Coleman in the text 'Long-term regulation of the circulation: interrelationships with body fluid volumes' Physical . (1967) bases of circulatory transport: Regulation and exchange 2 Focus the goals according to what was accomplished in the paper: "capturing the relevant mechanisms in an identifiable mathematical model" is too over-reaching since the three key equations (18, 19, and 21) do not capture mechanisms but are of phenomenological nature. "probing hypotheses regarding the long-term control of arterial pressure" -the model is both too limited and phenomenological to allow for any of various hypotheses. It is, testing however, fit for providing for two hypotheses: roles of the baroreflex a possible explanation involvement and arterial stiffening. Otherwise, for testing, one would require a more mechanistic description, inclusion of other key organs and regulatory mechanisms, and a more rigorous mathematical analyses. The authors provided an explanation of how the baroreflex may be involved in long-term control of blood pressure. This is a novel contribution to the field and should be the focus of the discussion: Add the discussion and references of experimental knowledge of baroreflex stimulation, inhibition, and especially long-term resetting. Describe how Guyton's models incorporate neural control and explain the above finding, especially the effects of chronic stimulation and how it is different from the proposed model (e.g. ) 3 Discuss existing knowledge or a possible physiological explanation for the finding that renin activity will be diminished with baroreflex stimulation. Discuss what other important physiological mechanisms are not reflected in the model and how they may play a role. In particular, explain why macula densa sodium delivery is thought not to be important enough to include in the model, despite the fact that it is known to be one of the key regulators of renin activity.
Do not confuse the model's predictions with its assumptions when discussing simulation results. The model assumes and is explicitly designed to tie the baroreflex, renin, and SNA with long-term pressure control. An explicitly designed assumption is not the same as prediction. Thus, statements like Conclusion 1 in "Summary and Conclusions" are inappropriate. The fact that the model's assumptions allow for physiologically plausible simulations merely supports the notion that baroreflex be important in the way the may authors have described it. Improve the clarity of the mathematical description: Provide motivation for the chosen elements of the model and outline the model in its entirety first. For example, a flowchart or a diagram of the model may help readers to grasp the overall idea. Add at least a basic analysis of parameter sensitivity, with clear description in the Methods section.
In Table 1 please provide physiologically-based description of parameters rather than

F1000Research
In Table 1 please provide physiologically-based description of parameters rather than referring to the equations where they are used (e.g. Emax is the maximal left-ventricular elastance, rather than ). Otherwise, the "Varying elastance heart model, " Equation (10) definitions are circular: parameters are used in the equations without an explanation of their meaning in either the text or the (2). | 37 PubMed Abstract Publisher Full Text We have read this submission. We believe that we have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however we have significant reservations, as outlined above.
No competing interests were disclosed. We are grateful to the reviewers for paying such serious attention to our paper and offering many constructive criticisms. While we have revised the paper to address the concerns and clarify the presentation, we disagree with the general interpretation of our model as "a conceptual framework which aims to illustrate how the baroreflex arc may be involved in long-term blood pressure control. " As stated in the summary, the model "was developed and identified to serve as a representation ." For example, there is no of the major physiological processes controlling arterial blood pressure equation or parameter value artificially injected into the model to require pressure-diuresis to not be able to return the pressure to a baseline value during baroreflex stimulation. In fact, we could parameterize the model so that the chronic pressure-diuresis phenomenon acts as an "infinite-gain" control system, always returning pressure to its set point. But in doing so, the model would not be able to match the data used to identify it, or to predict the data used to validate it.
Regarding the level of detail captured by the model, the reviewers are correct that at some level, all components of the model are phenomenological in nature. (This is true of all mathematical models.) Yet regarding the question around the level of detail and complexity of the model, we stand by the assertion that this model represents the most detailed model of long-term control of arterial pressure that is identified in the literature. (Here we are invoking the word "identified" to indicate that the model equations and associated parameter values are reported along with comparisons of model simulations of data used to identify the parameters.) F1000Research comparisons of model simulations of data used to identify the parameters.) In fact, one of the unstated goals of this paper is to inject some clarity and rigor into this field with precise and concrete language. Therefore we do not know how to respond to the statement that our model " " the versions of the Guyton model is not sophisticated enough to correspond to most of because it is not clear, precisely, what is meant by "sophisticated" or "correspond". In discussing the Guyton models, we clearly state that: "The 1972 realization of the model invokes approximately . 160 variables, and several hundred adjustable parameters More recent versions of the model ." Although it is true that involve many thousands of variables and tens of thousands of parameters our model has less variables and parameters than any version of the Guyton model that we are aware of, is the sophistication of a model measured purely by the number of variables or adjustable parameters? Since no version of the Guyton model captures pulsatile mechanics, no version of the model can capture the baroreceptor resetting phenomena illustrated in Figures 2 and 3 of our paper. Nor does the HumMod model capture the system response to baroreflex stimulation illustrated in Figure 9. Therefore one could argue that our 16-variable model is more sophisticated than the 5,000-variable Guyton model.

Regarding the specific suggestions:
"capturing the relevant mechanisms in an identifiable mathematical model" is too over-reaching since the three key equations (18, 19, and 21) do not capture mechanisms but are of phenomenological nature.
More specifically, we stated that one of the "goal" of the study is "capturing the relevant ." We believe we have done that. However, in mechanisms in an identifiable mathematical model the revised paper we have changed the goal to "representing the relevant mechanisms in an ." identifiable mathematical model "probing hypotheses regarding the long-term control of arterial pressure" -the model is both too limited and phenomenological to allow for any testing of various hypotheses.
Again, we counter that we have met the stated goal. We are not claiming to have proved any hypothesis. Yet we have introduced some important disproofs, as outlined in some of our major conclusions. Specifically, it is stated in the literature as 'scientific facts' that; the renal function curve represents the dominant long-term controller of blood pressure and, that since the baroreflex "resets" it cannot play a role in chronic changes in mean pressure. Our model analysis has cast significant doubt upon the validity of both of these hypotheses. Admittedly, these findings represent somewhat of a double negative i.e., a model capturing baroreflex resetting is used to demonstrate how a baroreflex that resets could play a role in long-term changes in pressure. Therefore, we have shown that the fact that the baroflex resets, does not mean that the baroreflex cannot play a role in long-term pressure control. Therefore we are confident that it is appropriate to state that we have constructed a model with the goal of "probing hypotheses regarding the " long-term control of arterial pressure.
Add the discussion and references of experimental knowledge of baroreflex stimulation, inhibition, and especially long-term resetting. Figures 2 and 3 represent direct experimental interrogation of baroreflex stimulation, inhibition, and long-term resetting. The resetting phenomenon is also indirectly probed in Figure 6 where, over the time course of the experiment, the baroreflex response to the acute pressure F1000Research where, over the time course of the experiment, the baroreflex response to the acute pressure increase gradually diminishes. The phenomenon and experimental observations are discussed in association with our simulations of these experiments. These are also discussed in the section on arterial stiffening and the etiology of hypertension. Finally, we have added a couple of sentences (end of Results Section 5) comparing our model predictions with the long-term observations of Thrasher on baroreceptor unloading.

Data in
Describe how Guyton's models incorporate neural control and explain the above finding, especially the effects of chronic stimulation and how it is different from the proposed model...
A discussion on this point has been added to the section "Results 4". Discuss existing knowledge or a possible physiological explanation for the finding that renin activity will be diminished with baroreflex stimulation.
It is not our finding that renin activity will diminish with baroreflex stimulation. That is the prior experimental observation. The mechanism in the model is through the observed dependence of renin release on sympathetic tone. At a deeper level, the mechanism may simply be explained as follows: decreased sympathetic tone leads to decreased renal arterial tone, leading to increased glomerular filtration and/or medullary blood flow, leading to increased flux of sodium to the macula densa, leading to decreased renin production. While we would rather leave such speculation out of the manuscript, it is published here as part of the open review process.
Discuss what other important physiological mechanisms are not reflected in the model and how they may play a role... explain why macula densa sodium delivery is thought not to be important enough to include in the model, despite the fact that it is known to be one of the key regulators of renin activity.
It is not that we believe that "macula densa sodium delivery is…not...important enough to include in ." Rather the phenomenon is captured in a simplified way. To clarify this point, we have the model added the following text to the section on "Neurohumoral control of pressure-diuresis/natriuresis": "These expressions simplify the known (and unknown) mechanisms governing renin production into a simple phenomenological relationship between pressure, sympathetic tone, and the kinetics ." of the renin-angiotensin system A more biophysically detailed model would be required to represent the direct link between sodium transport and renin production. At the bare minimum, such a model would represent the biophysical processes linking pressure and solute transport in the kidney, including effects of sympathetic innervation on arterial smooth muscle stress development, mechanics of blood and filtrate flow, solute advection, diffusion, permeation, exchange, and pumping.
Do not confuse the model's predictions with its assumptions when discussing simulation results... Conclusion 1 in "Summary and Conclusions" is inappropriate. The fact that the model's assumptions allow for physiologically plausible simulations merely supports the notion that baroreflex may be important in the way the authors have described it.
We agree that this first stated conclusion does not represent a conclusion and what we are trying to say is now better captured in the revised text for Conclusion 1.