Impact Factor 3.394

The world's 3rd most-cited Physiology journal

Original Research ARTICLE

Front. Physiol., 08 February 2013 | https://doi.org/10.3389/fphys.2013.00004

A detailed physiologically based model to simulate the pharmacokinetics and hormonal pharmacodynamics of enalapril on the circulating endocrine renin-angiotensin-aldosterone system

  • 1School of Engineering and Science, Jacobs University Bremen, Bremen, Germany
  • 2Computational Systems Biology, Bayer Technology Services GmbH, Leverkusen, Germany
  • 3Modelling and Simulation, Fraunhofer MEVIS, Bremen, Germany

The renin-angiotensin-aldosterone system (RAAS) plays a key role in the pathogenesis of cardiovascular disorders including hypertension and is one of the most important targets for drugs. A whole body physiologically based pharmacokinetic (wb PBPK) model integrating this hormone circulation system and its inhibition can be used to explore the influence of drugs that interfere with this system, and thus to improve the understanding of interactions between drugs and the target system. In this study, we describe the development of a mechanistic RAAS model and exemplify drug action by a simulation of enalapril administration. Enalapril and its metabolite enalaprilat are potent inhibitors of the angiotensin-converting-enzyme (ACE). To this end, a coupled dynamic parent-metabolite PBPK model was developed and linked with the RAAS model that consists of seven coupled PBPK models for aldosterone, ACE, angiotensin 1, angiotensin 2, angiotensin 2 receptor type 1, renin, and prorenin. The results indicate that the model represents the interactions in the RAAS in response to the pharmacokinetics (PK) and pharmacodynamics (PD) of enalapril and enalaprilat in an accurate manner. The full set of RAAS-hormone profiles and interactions are consistently described at pre- and post-administration steady state as well as during their dynamic transition and show a good agreement with literature data. The model allows a simultaneous representation of the parent-metabolite conversion to the active form as well as the effect of the drug on the hormone levels, offering a detailed mechanistic insight into the hormone cascade and its inhibition. This model constitutes a first major step to establish a PBPK-PD-model including the PK and the mode of action (MoA) of a drug acting on a dynamic RAAS that can be further used to link to clinical endpoints such as blood pressure.

Introduction

LHypertension is a worldwide epidemic, which affects all ages and racial populations. It has a very high incidence and is the leading cause of cardiovascular mortality (Roger et al., 2012). There is still a lack of mechanistic models in the published literature that integrate the available physiological knowledge and relate the pharmacokinetics (PK) of hypertension drugs to their pharmacodynamics (PD) effects. One model addressing among others the RAAS and the blood pressure, is the Guyton’s physiological model (Guyton et al., 1972), it was extended and published recently (Montani and Van Vliet, 2009; Osborn et al., 2009). Although this and other models may contain all the parameters listed in this present study and although they might have been validated against several physiological situations they are not based on a generic whole body physiologically based PK (wb PBPK) background model. Other approaches concerning the hormone cascade, as for example by White et al. (1989), Levitt and Schoemaker (2006), Guillaud and Hannaert (2010), Zhou et al. (2012), have different focuses or approaches.

To overcome the lack of detailed physiologically based models for cardiovascular PD, we aimed to establish a mechanistic model for the RAAS covering all relevant biological processes depicted in Figure 1. Renin is synthesized from the enzyme precursor prorenin, which is produced mainly in the kidney and secreted into the plasma by the granular cells (Krop et al., 2008). The active enzyme renin is stored in the cells of the juxtaglomerular apparatus is released immediately on stimulation of these cells (Danser et al., 1998). After secretion into the plasma, renin converts the hepatically synthesized inactive hormone angiotensinogen (AGT) to angiotensin 1 (Ang1). Ang1 is then converted in the plasma compartment by the membrane bound angiotensin-converting-enzyme (ACE), expressed by the vascular endothelium, to angiotensin 2 (Ang2). Ang2 has an inhibiting effect on renin synthesis and secretion and thus on the plasma levels of Ang1 and Ang2 (Johns et al., 1990). Ang2 is a vasoactive peptide that induces volume (and Na) retaining effects as well as vasoconstricting effects and an increase in blood volume, leading to increased blood pressure. It potentiates aldosterone secretion by interacting with its angiotensin 2 receptor type 1 (AT1). Aldosterone increases the blood pressure as well (Gornall et al., 1960) by sodium and extracellular volume retention. Several cardiovascular disease states are associated with changes in circulating Ang2 (Laragh et al., 1972; Laragh, 1995).

FIGURE 1
www.frontiersin.org

Figure 1. Schematic representation of the renin-angiotensin-aldosterone system (RAAS). Renin converts the hepatically secreted angiotensinogen (AGT) to angiotensin 1 (Ang1). Ang1 is then converted to angiotensin 2 (Ang2) by the membrane bound angiotensin-converting-enzyme (ACE). Ang2 then binds to the angiotensin 2 receptor type 1 (AT1). The most important effects of AT1 binding are depicted either: the aldosterone secretion by the adrenals and the inhibition of renin activation from prorenin. Different processes are depicted by different line styles (see upper right box).

Inhibitors of the RAAS can be classified by their mode of action (MoA; Cagnoni et al., 2010; Eisenberg, 2011; Aronson and Krum, 2012) into five categories:

ACE inhibitors (ACEIs),

Direct renin inhibitors (DRIs),

Angiotensin-2-receptor-blockers (ARBs),

Aldosterone synthase inhibitors (ASIs), and

Potassium sparing diuretics.

A common drug of the first type of inhibitors is enalaprilat (Enaat), a potent and reversible ACE inhibitor (MacFadyen et al., 1993; Dhareshwar, 2007). Enaat is poorly absorbed from the gastrointestinal tract. Therefore, an ethyl ester promoiety was incorporated within the administered parent drug enalapril (Ena) to increase its oral bioavailability. After absorption of inactive Ena by the peptide carrier-mediated transport system in the intestine, the prodrug Ena is converted to Enaat by the carboxylesterase family of enzymes activating it for ACE inhibition (MacFadyen et al., 1993; Dhareshwar, 2007).

Materials and Methods

A model of the RAAS based on available information was established (Eissing et al., 2011). To evaluate the RAAS model, a generic PBPK model of Ena and Enaat was developed and validated with clinical data that was then coupled to the RAAS model.

Software

The RAAS model was created using the Systems Biology software suite that consists of the wb PBPK software tool PK-Sim® (Version 4.2, Bayer Technology Services, Leverkusen, Germany) and the general purpose modeling software MoBi® (Version 2.3, Bayer Technology Services, Leverkusen, Germany). This software has successfully been applied for diverse modeling studies in laboratory animals and humans (Willmann et al., 2003, 2004, 2007, 2010; Edginton et al., 2006a,b,c, 2008; Brochot et al., 2010; Eissing et al., 2011). The parameter identifications for MoBi models were done using MATLAB® (R2010b, The MathWorks, Inc., Natick, MA, USA) and the MoBi Toolbox for MATLAB. Identification algorithm for Ena and Enaat models was Monte Carlo, while RAAS identification algorithm was fminsearch (fminsearch uses the Nelder–Mead Simplex method). The identification error function was in both cases error root mean square. For additional information, see the software manual or http://www.systems-biology.com/products.

Development of the RAAS-PBPK Model

Physiological data presented in Table 1 contains the information gained from the literature or, if not available by optimization based on data for Ang1, Ang2, aldosterone, renin, prorenin, ACE, AGT, and the AT1. Optimization was done with the overall coupled model and with all RAAS parameters (Table 1) simultaneously. For each of them an individual PBPK model was developed which contains many physiological background parameters (Eissing et al., 2011) that are derived from the few input parameters presented in Table 1. Mass transfer within the plasma (r1) between organs is calculated as organ specific blood flow (Qorgan) times the drug concentration within the plasma of the respective organ (Cpls_organ) times the overall plasma fraction (1 − HCT) which is one subtracted by the hematocrit (HCT) value. Passive transport into the interstitial and intracellular space is realized by multiplying organ specific partition coefficients Korg.

r1=Qorgan1-HCTCpls_organ
TABLE 1
www.frontiersin.org

Table 1. Parameters and variables used for the different actors in the RAAS.

The interactions in the coupled model were described as depicted in Figure 1 reflecting the biological processes as detailed in the following. AGT, prorenin, renin, and aldosterone are the four hormones and enzymes synthesized in cells and then secreted endocrinally. It was assumed that these hormones and enzymes are produced intracellular and secreted into the plasma by a first order rate (r2 in μmol· min−1) proportional to the intracellular concentration (CCellularHormone) in the respective organ with the rate constant kSecretion (in L· min−1). The renin secretion r3 (μmol· min−1) is described by a “competitive” inhibition equation including the conversion of prorenin to renin influenced by Ang2 (see Figure 1), describing the conversion from prorenin to renin, and simultaneously describing the secretion from intracellular to plasma compartment of the kidney:

r3=VMaxKidVKidneyCProrenintCProrenint+KmKid1+CAng2tKinhibition

KmKid and VMaxKid are the respective enzymatic activity constants, Vkidney is the intracellular volume of the kidney, CProrenin is the renal intracellular prorenin concentration (parameters under prorenin in Table 1). Kinhibition represents the inhibitory constant for Ang2 and CAng2 is the renal plasma concentration of Ang2. The degradation rate r4 (μmol· min−1· L−1) is described by Guillaud and Hannaert (2010):

r4=Ctln2t12

t1/2 determines the half-life and C(t) being the respective plasma concentration over time. These half-lives were optimized for the hormones and enzymes (see Table 1). The conversion of AGT to Ang1 and from Ang1 to Ang2 take place in the plasma compartments of all organs (Eissing et al., 2011) and is an enzymatically regulated process described by a Michaelis–Menten equation (Gould et al., 1980; Ehlers and Kirsch, 1988). The change of hormone mass over time (r5 in μmol· min−1) then depends on the enzyme concentration (C0) and the prohormone (CProhormone) concentration in the respective organs, as well as on the enzyme activity represented by the Km and kcat of the enzyme (Meyer et al., 2012). This can be expressed by:

r5=kcatV0C0tCProhormonetCProhormonet+Km

The calculation refers to the respective organ plasma volumes V0 in L (taken from PK-Sim® database). ACE is most notably bound to the endothelium (Dzau et al., 2001) and thus not circulating in the model. The concentration of ACE and AT1 in each organ (C0) is calculated via a reference concentration (CRef, see Table 3) times the percentage of reference expression (ERef) from literature (see Table 4, factor ERef = 100 for the highest expression value):

C0=CRefERef100

Ang2 binds to the AT1 receptor, which influences the aldosterone secretion. AT1 binding rate (r7 in μmol· min−1· L−1) is described by a typical equation with k+ and k binding constants for association and dissociation and the respective concentrations for AT1 (C0), Ang2 (C2), and the AT1-Ang2 complex (CComplex):

r7=k+C0tC2t-k_Ccomplext

The AT1 concentrations were calculated as explained in Eq. 6, the respective relative expression values (ERef) are shown in Table 4, and the reference concentration (CRef) is presented in Table 1. As shown in Figure 1 AT1 is membrane bound too and thereby does not circulate in the plasma. Aldosterone synthesis (r8 in μmol· min−1) in our model was directly linked to the complex formation of Ang2 bound to the AT1 (CComplex), which is modeled as

r8=ksecretionCcomplext+kproduction

The secretion is thereby described to be dependent on intracellular concentration of Ang2 receptor complex multiplied by the rate constant kSecretion (in L· min−1) and additionally represents the continuous secretion as kProduction (Veldhuis et al., 2008).

Development of the Enalapril/Enalaprilat PBPK Model

A PBPK model for oral Ena including the metabolism to active Enaat was developed first and used as test case for the dynamic behavior of the RAAS model. Based on intravenous (IV) data, parameters for Enaat were determined and then used in the parent-metabolite model. The required physicochemical input properties and clearance processes for the model are presented in Table 2. The model simulations for oral Ena were compared to different experimental data (Biollaz et al., 1982; Noormohamed et al., 1990; Nussberger et al., 2002; Lee et al., 2003; Najib et al., 2003; Gu et al., 2004; Lu et al., 2009). The six different subject and dataset properties (all included to the respective simulations) are shown in Table 3. Only the five parameters presented in Table 3 were optimized with MATLAB®. The respective IV Enaat simulations were compared to experimental data of Enaat plasma concentrations after IV Enaat (Hockings et al., 1986) administration. The coupling of parent and metabolite was realized with a hydrolytic reaction in the hepatocytes (Dhareshwar, 2007). This reaction is mediated by different carboxylesterase enzymes. Km values for metabolism to Enaat in Table 1 were taken from measurements in rat hepatocytes (Tabata et al., 1990; Abu-Zahra and Pang, 2000). The conversion rate (r9 in μmol· min−1) from Ena into Enaat in the hepatocytes with the respective intracellular volume VLiver (see Table 3), the concentrations of Ena (CEna) and Enaat (CEnaat) is modeled by Michaelis–Menten kinetics:

r9=VmaxLivVLiverCEnatCEnat+KmLiv
TABLE 2
www.frontiersin.org

Table 2. Physicochemical and PK properties of enalapril and enalaprilat.

TABLE 3
www.frontiersin.org

Table 3. Population specific information about subjects, optimized parameters, and model specific parameters.

The respective parameters for Michaelis–Menten kinetics are VMaxLiv and KmLiv.

Assembling of the Models to a Full Dynamic PBPK-RAAS Model

To couple the drug model to the full RAAS model the MoA was included. Enaat binds to ACE in the RAAS (see Table 1, ACE parameters) generating an inactive complex. The formation rate of Enaat-ACE complex (CComplex) is realized as presented in equation 6 for all plasma compartments. ACE concentration is based on the relative expression value in Table 4 which was obtained from Unigene (2010) and processed as described in Meyer et al. (2012). The inhibition of ACE occurs by ACE binding to Enaat, forming an ACE-Enaat- complex, thereby blocking ACE for the interaction with Ang1. The binding is implemented in the plasma compartment of each organ.

TABLE 4
www.frontiersin.org

Table 4. Reference expression values of ACE per organ as included in PK-Sim 5® expression database (Meyer et al., 2012) from Unigene (2010).

Results

Oral Administration of Enalapril

A PBPK model of IV Enaat has been developed with the data indicated in Table 2. In Figure 2 the resulting PK-simulations for IV Enaat are shown and the goodness-of-fit of the model based on data is depicted (Hockings et al., 1986). This model was used to determine the parameters related to Enaat, which are log P, Molecular weight, pKa, plasma protein binding and clearance (see Table 2). Additionally the constants of ACE binding were identified (see Table 1). The constant for association is much higher than the dissociation constant, as stated by Reynolds (1984). This parameterization was used within the coupled model without changes.

FIGURE 2
www.frontiersin.org

Figure 2. (A) Plasma concentration-time profiles of Enaat after IV Enaat administration. The solid blue line is the simulation result compared to the data of Hockings et al. (1986) which are indicated as points (adolescents) and triangles (elderly). (B) The simulated versus observed plot for the same simulation with the Hockings data are shown, indicated as points (adolescents) and triangles (elderly). The solid black regression line denotes the line of identity (observed outcome is equal to simulated outcome).

Resulting plasma concentrations of ACE based on the explained calculation (peripheral venous plasma concentration in the model validated by data from Gu et al. (2004): 7.17E−4 μmol· L−1) correspond nicely with literature data of 2.52E−3 μmol· L−1 by Brice et al. (1995). After integrating the metabolite model into the oral model, the coupled model exhibits the behavior presented in Figure 3. Including the population dependent parameters for the different publications (see Table 3), the six model simulations show excellent agreement with the data. The serum concentration profile of Enaat exhibits a prolonged terminal phase, apparently representing the small fraction of the administered dose that has been bound to ACE. As obvious from Figure 3, the model Ena simulations and Enaat predictions accurately fit the observed data.

FIGURE 3
www.frontiersin.org

Figure 3. (A) Exemplary plasma concentration-time curves of Ena and Enaat from the coupled Ena-Enaat model after oral Ena application compared to experimental data (Gu et al., 2004). The red solid line is the simulated Ena plasma concentration, while the blue solid line is the simulation result of Enaat. The respective data of Ena are indicated as circles, of Enaat as triangles. (B) Simulated versus observed plot for Enalapril from the coupled model, compared to the data by Biollaz et al. (1982), Noormohamed et al. (1990), Lee et al. (2003), Najib et al. (2003), Gu et al. (2004), Lu et al. (2009), with SD. The solid black regression line denotes the line of identity (observed outcome is equal to simulated outcome). (C) Predicted versus observed plot for Enalaprilat from the coupled Ena-Enaat model, compared to the data by Biollaz et al. (1982), Noormohamed et al. (1990), Lee et al. (2003), Najib et al. (2003), Gu et al. (2004), Lu et al. (2009). The solid black regression line denotes the line of identity (observed outcome is equal to simulated outcome).

RAAS Inhibition by Enalaprilat

The PBPK model for the circulating, endocrine RAAS was developed to describe the steady state described in literature. It corresponds to the pre-administration steady state discussed further in the following. The model for Ena was then coupled to the RAAS model to explore the dynamical behavior of the system under inhibition. In Figure 4 the results for the full RAAS model under administration of 10 mg Ena are shown at t = 1500 min (administration at t = 1500 indicated by the red arrow). All parameters of the Ena-Enaat model were used without changes. The pre- and post-administration steady states as well as the dynamical transition after administration at t = 1500 min (red arrow) show excellent agreement to literature data (Table 1, Figure 4A). In detail, Figure 4 shows the MoA of Ena. After the administration of Ena, ACE is inhibited. As can be seen, the conversion of Ang1 to Ang2 is inhibited, lowering the Ang2 and aldosterone plasma levels. The negative feedback loop from the Ang2 to the renin synthesis is inhibited, leading to an increase of the renin plasma concentration, and thereby the Ang1 levels. The comparison of these induced changes in plasma levels by a scatter plot (Figure 4B) to study data from Nussberger et al. (2002) shows good agreement and corresponds also to older experimental data (Given et al., 1984; Johnston et al., 1984; Juillerat et al., 1990; Azizi et al., 1997).

FIGURE 4
www.frontiersin.org

Figure 4. (A) Plasma concentration-time curves of hormones included into the RAAS. Hormones are in a steady state until administration of oral Ena at t = 1500 min indicated by the red arrow. The respective data are indicated by different symbols (Nussberger et al., 2002). (B) Simulated versus observed plasma concentrations of different actors of the RAAS before and after inhibition of ACE by Ena. The different hormones are depicted by different symbols (Nussberger et al., 2002) with their SD. The solid black regression line denotes the line of identity (observed outcome is equal to simulated outcome).

The simulation of exemplary Ena administration shows, that the full dynamical behavior of inhibition and the other RAAS-hormone levels can be represented accurately (see Figure 4).

Discussion

A full mechanistic model of the circulating RAAS was developed. The model represents the steady state and dynamic behavior of the RAAS in an accurate manner. This was demonstrated by simulating the effects of representative drug acting on the system. The Ang1 concentration is influenced by different reactions, Ang1 generation, Ang2 generation, degradation of Ang1, which lead to the pre-administration steady state. So it is rather difficult to evaluate the condition of substrate excess for the second, since the new generation of Ang1 is fast. In addition the reference concentrations C0 are artificial constructs in a product with kcat, where both and Km are fitted parameters. Thus Michaelis–Menten equations were used, although the C0 is much higher than the Ang1 concentration for the steady state, as it was used for these enzymes to catalyze this reaction in the described way (Ehlers and Kirsch, 1988; Guillaud and Hannaert, 2010). Although we find an excellent correlation of plasma levels and a good correlation of parameters to the articles cited before, it has not been proven in this article, that the optimized parameter values are unique to reach the reference data, since not for all parameters measured values are available. The Ena model was validated independently from the parameterization of the RAAS and can thus also serve as a validation for the RAAS. Following the development of the Ena model based on literature data, subsequently only population dependent parameters were changed (listed in Table 3) to represent physiological differences in the studies. After coupling to the RAAS model the overall model was able to represent detailed information of all main components and predict circulating plasma concentrations of all included hormones by one coupled PBPK/PD model. The excellent representation of the hormone levels at steady state and during inhibition, implicate a good representation of the underlying physiological processes in the RAAS.

The detailed physiological representation of relevant processes spanning from wb distribution to molecular interactions allows for, e.g., a mechanistic implementation of independent knowledge on drug-target interactions, and thereby an integration and translation of available knowledge.

Thus the model constitutes a first major step to establish a PBPK-PD-model for drugs acting on blood pressure and heart rate. Generally blood pressure and heart rate represent the state of the art endpoint in clinical studies. Currently blood pressure is not yet included in the model. The inclusion of inhibitors from other classes of MoA’s can be used to explore pharmacodynamic effects in the same model framework. Further extensions of this model and an inclusion of the dynamic behavior of the blood pressure in response to changes in RAAS would allow a full mechanistic representation of the PBPK-PD relationship in cardiovascular diseases. Such a model could thus rationalize treatment decisions in the area of cardiovascular diseases. By inclusion of a full process description from administration of a drug to its effect, i.e., the lowering of blood pressure, such a model allows to support cardiovascular drug research from target identification and validation to dosing regimen decision.

Conflict of Interest Statement

Karina Claassen is funded by and Thomas Eissing, Michael Block, and Stefan Willmann are employes of Bayer Technology Services GmbH, the company that owns and commercializes the software platform used for the simulations described in the manuscript (PK-Sim® and MoBi®), as well as parent company stock owners.

References

Abu-Zahra, T. N., and Pang, K. S. (2000). Effect of zonal transport and metabolism on hepatic removal: enalapril hydrolysis in zonal, isolated rat hepatocytes in vitro and correlation with perfusion data. Drug Metab. Dispos. 28, 807–813.

Pubmed Abstract | Pubmed Full Text

Aronson, D., and Krum, H. (2012). Novel therapies in acute and chronic heart failure. Pharmacol. Ther. 23, 23.

Azizi, M., Guyene, T. T., Chatellier, G., Wargon, M., and Menard, J. (1997). Additive effects of losartan and enalapril on blood pressure and plasma active renin. Hypertension 29, 634–640.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Biollaz, J., Schelling, J. L., Jacot Des Combes, B., Brunner, D. B., Desponds, G., Brunner, H. R., et al. (1982). Enalapril maleate and a lysine analogue (MK-521) in normal volunteers; relationship between plasma drug levels and the renin angiotensin system. Br. J. Clin. Pharmacol. 14, 363–368.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Boyd, G. W., Landon, J., and Peart, W. S. (1967). Radioimmunoassay for determining plasms-levels of angiotensin II in man. Lancet 2, 1002–1005.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Brice, E. A., Friedlander, W., Bateman, E. D., and Kirsch, R. E. (1995). Serum angiotensin-converting enzyme activity, concentration, and specific activity in granulomatous interstitial lung disease, tuberculosis, and COPD. Chest 107, 706–710.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Brochot, A., Zamacona, M., and Stockis, A. (2010). Physiologically based pharmacokinetic/pharmacodynamic animal-to-man prediction of therapeutic dose in a model of epilepsy. Basic Clin. Pharmacol. Toxicol. 106, 256–262.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cagnoni, F., Njwe, C. A., Zaninelli, A., Ricci, A. R., Daffra, D., D’Ospina, A., et al. (2010). Blocking the RAAS at different levels: an update on the use of the direct renin inhibitors alone and in combination. Vasc. Health Risk Manag. 6, 549–559.

Pubmed Abstract | Pubmed Full Text

Danser, A. H., Derkx, F. H., Schalekamp, M. A., Hense, H. W., Riegger, G. A., and Schunkert, H. (1998). Determinants of interindividual variation of renin and prorenin concentrations: evidence for a sexual dimorphism of (pro)renin levels in humans. J. Hypertens. 16, 853–862.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Dhareshwar, S. S. (2007). “Case study: enalapril: a prodrug of enalaprilat,” in Prodrugs, eds. V. J. Stella, R. T. Borchardt, M. J. Hageman, R. Oliyai, H. Maag, and J. W. Tilley (New York: Springer), 1221–1229.

Dzau, V. J., Bernstein, K., Celermajer, D., Cohen, J., Dahlof, B., Deanfield, J., et al. (2001). The relevance of tissue angiotensin-converting enzyme: manifestations in mechanistic and endpoint data. Am. J. Cardiol. 88, 1L–20L.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Edginton, A. N., Schmitt, W., Voith, B., and Willmann, S. (2006a). A mechanistic approach for the scaling of clearance in children. Clin. Pharmacokinet. 45, 683–704.

CrossRef Full Text

Edginton, A. N., Schmitt, W., and Willmann, S. (2006b). Application of physiology-based pharmacokinetic and pharmacodynamic modeling to individualized target-controlled propofol infusions. Adv. Ther. 23, 143–158.

CrossRef Full Text

Edginton, A. N., Schmitt, W., and Willmann, S. (2006c). Development and evaluation of a generic physiologically based pharmacokinetic model for children. Clin. Pharmacokinet. 45, 1013–1034.

CrossRef Full Text

Edginton, A. N., Theil, F. P., Schmitt, W., and Willmann, S. (2008). Whole body physiologically-based pharmacokinetic models: their use in clinical drug development. Expert Opin. Drug Metab. Toxicol. 4, 1143–1152.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ehlers, M. R., and Kirsch, R. E. (1988). Catalysis of angiotensin I hydrolysis by human angiotensin-converting enzyme: effect of chloride and pH. Biochemistry 27, 5538–5544.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Eisenberg, J. M. (2011). ACEIs, ARBs, or DRI for Adults with Hypertension. Rockville, MD: Agency for Healthcare Research and Quality.

Eissing, T., Kuepfer, L., Becker, C., Block, M., Coboeken, K., Gaub, T., et al. (2011). A computational systems biology software platform for multiscale modeling and simulation: integrating whole-body physiology, disease biology, and molecular reaction networks. Front. Physiol. 2:4. doi:10.3389/fphys.2011.00004

CrossRef Full Text

Given, B. D., Taylor, T., Hollenberg, N. K., and Williams, G. H. (1984). Duration of action and short-term hormonal responses to enalapril (MK 421) in normal subjects. J. Cardiovasc. Pharmacol. 6, 436–441.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gornall, A. G., Grundy, H. M., and Koladich, C. J. (1960). Studies on aldosterone. III. Chronic effect on the blood pressure of rats. Can. J. Biochem. Physiol. 38, 43–50.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gould, A. B., Dewolf, R., Goodman, S., Onesti, G., and Swartz, C. (1980). Kinetic studies of the human renin and human substrate reaction. Biochem. Med. 24, 321–326.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gu, Q., Chen, X., Zhong, D., and Wang, Y. (2004). Simultaneous determination of enalapril and enalaprilat in human plasma by liquid chromatography-tandem mass spectrometry. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 813, 337–342.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Guillaud, F., and Hannaert, P. (2010). A computational model of the circulating renin-angiotensin system and blood pressure regulation. Acta Biotheor. 58, 143–170.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Guyton, A. C., Coleman, T. G., and Granger, H. J. (1972). Circulation: overall regulation. Annu. Rev. Physiol. 34, 13–46.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hockings, N., Ajayi, A. A., and Reid, J. L. (1986). Age and the pharmacokinetics of angiotensin converting enzyme inhibitors enalapril and enalaprilat. Br. J. Clin. Pharmacol. 21, 341–348.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Johns, D. W., Peach, M. J., Gomez, R. A., Inagami, T., and Carey, R. M. (1990). Angiotensin II regulates renin gene expression. Am. J. Physiol. 259, F882–F887.

Pubmed Abstract | Pubmed Full Text

Johnston, C. I., Jackson, B. J., Larmour, I., Cubella, R., and Casley, D. (1984). Plasma enalapril levels and hormonal effects after short- and long-term administration in essential hypertension. Br. J. Clin. Pharmacol. 18(Suppl. 2), 233S–239S, 241S.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Juillerat, L., Nussberger, J., Menard, J., Mooser, V., Christen, Y., Waeber, B., et al. (1990). Determinants of angiotensin II generation during converting enzyme inhibition. Hypertension 16, 564–572.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Katsurada, A., Hagiwara, Y., Miyashita, K., Satou, R., Miyata, K., Ohashi, N., et al. (2007). Novel sandwich ELISA for human angiotensinogen. Am. J. Physiol. Renal Physiol. 293, 6.

CrossRef Full Text

Knox, C., Law, V., Jewison, T., Liu, P., Ly, S., Frolkis, A., et al. (2011). DrugBank 3.0: a comprehensive resource for “omics” research on drugs. Nucleic Acids Res. 39, 8.

CrossRef Full Text

Kobori, H., Alper, A. B. Jr., Shenava, R., Katsurada, A., Saito, T., Ohashi, N., et al. (2009). Urinary angiotensinogen as a novel biomarker of the intrarenal renin-angiotensin system status in hypertensive patients. Hypertension 53, 344–350.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Krop, M., De Bruyn, J. H., Derkx, F. H., and Danser, A. H. (2008). Renin and prorenin disappearance in humans post-nephrectomy: evidence for binding? Front. Biosci. 13, 3931–3939.

Pubmed Abstract | Pubmed Full Text

Laragh, J. H. (1995). Renin-angiotensin-aldosterone system for blood pressure and electrolyte homeostasis and its involvement in hypertension, in congestive heart failure and in associated cardiovascular damage (myocardial infarction and stroke). J. Hum. Hypertens. 9, 385–390.

Pubmed Abstract | Pubmed Full Text

Laragh, J. H., Baer, L., Brunner, H. R., Buhler, F. R., Sealey, J. E., and Vaughan, E. D. Jr. (1972). Renin, angiotensin and aldosterone system in pathogenesis and management of hypertensive vascular disease. Am. J. Med. 52, 633–652.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lee, J., Son, J., Lee, M., Lee, K. T., and Kim, D. H. (2003). Simultaneous quantitation of enalapril and enalaprilat in human plasma by 96-well solid-phase extraction and liquid chromatography/tandem mass spectrometry. Rapid Commun. Mass Spectrom. 17, 1157–1162.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Levitt, D. G., and Schoemaker, R. C. (2006). Human physiologically based pharmacokinetic model for ACE inhibitors: ramipril and ramiprilat. BMC Clin. Pharmacol. 6:1. doi:10.1186/1472-6904-6-1

CrossRef Full Text

Lu, S., Jiang, K., Qin, F., Lu, X., and Li, F. (2009). Simultaneous quantification of enalapril and enalaprilat in human plasma by high-performance liquid chromatography-tandem mass spectrometry and its application in a pharmacokinetic study. J. Pharm. Biomed. Anal. 49, 163–167.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

MacFadyen, R. J., Meredith, P. A., and Elliott, H. L. (1993). Enalapril clinical pharmacokinetics and pharmacokinetic-pharmacodynamic relationships. An overview. Clin. Pharmacokinet. 25, 274–282.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Meyer, M., Schneckener, S., Ludewig, B., Kuepfer, L., and Lippert, J. (2012). Using expression data for quantification of active processes in physiologically-based pharmacokinetic modeling. Drug Metab. Dispos. 40, 892–901.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Montani, J. P., and Van Vliet, B. N. (2009). Understanding the contribution of Guyton’s large circulatory model to long-term control of arterial pressure. Exp. Physiol. 94, 382–388.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Najib, N. M., Idkaidek, N., Adel, A., Admour, I., Astigarraga, R. E., De Nucci, G., et al. (2003). Bioequivalence evaluation of two brands of enalapril 20 mg tablets (Narapril and Renitec) in healthy human volunteers. Biopharm. Drug Dispos. 24, 315–320.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Noda, T., Yaginuma, T., O’Rourke, M. F., and Hosoda, S. (2006). Effects of nifedipine on systemic and pulmonary vascular impedance in subjects undergoing cardiac catheterization. Hypertens. Res. 29, 505–513.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Noormohamed, F. H., McNabb, W. R., and Lant, A. F. (1990). Pharmacokinetic and pharmacodynamic actions of enalapril in humans: effect of probenecid pretreatment. J. Pharmacol. Exp. Ther. 253, 362–368.

Pubmed Abstract | Pubmed Full Text

Nussberger, J., Wuerzner, G., Jensen, C., and Brunner, H. R. (2002). Angiotensin II suppression in humans by the orally active renin inhibitor Aliskiren (SPP100): comparison with enalapril. Hypertension 39, E1–E8.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Osborn, J. W., Averina, V. A., and Fink, G. D. (2009). Current computational models do not reveal the importance of the nervous system in long-term control of arterial pressure. Exp. Physiol. 94, 389–396.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pardridge, W. M. (2011). Targeted delivery of hormones to tissues by plasma proteins. Compr. Physiol. 335–382.

Ranadive, S. A., Chen, A. X., and Serajuddin, A. T. (1992). Relative lipophilicities and structural-pharmacological considerations of various angiotensin-converting enzyme (ACE) inhibitors. Pharm. Res. 9, 1480–1486.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Remko, M. (2007). Acidity, lipophilicity, solubility, absorption, and polar surface area of some ACE inhibitors. Chem. Pap. 61, 133–141.

CrossRef Full Text

Reynolds, C. H. (1984). Kinetics of inhibition of angiotensin converting enzyme by captopril and by enalapril diacid. Biochem. Pharmacol. 33, 1273–1276.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Roger, V. L., Go, A. S., Lloyd-Jones, D. M., Benjamin, E. J., Berry, J. D., Borden, W. B., et al. (2012). Heart disease and stroke statistics–2012 update: a report from the American Heart Association. Circulation 125, e2–e220.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sealey, J. E., Atlas, S. A., and Laragh, J. H. (1980). Prorenin and other large molecular weight forms of renin. Endocr. Rev. 1, 365–391.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sirianni, G. L., and Pang, K. S. (1998). Intracellular and not intraluminal esterolysis of enalapril in kidney. Studies with the single pass perfused nonfiltering rat kidney. Drug Metab. Dispos. 26, 324–331.

Pubmed Abstract | Pubmed Full Text

Tabata, S., Yamazaki, H., Ohtake, Y., and Hayashi, S. (1990). Rapid conversion of the new angiotensin converting enzyme inhibitor ramipril to its active metabolite in rats. Arzneimittelforschung 40, 865–867.

Pubmed Abstract | Pubmed Full Text

Thevananther, S., and Brecher, A. S. (1999). Isolation of angiotensin converting enzyme (ACE) binding protein from human serum with an ACE affinity column. Can. J. Physiol. Pharmacol. 77, 216–223.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tu, W., Eckert, G. J., Pratt, J. H., and Jan Danser, A. H. (2012). Plasma levels of prorenin and renin in blacks and whites: their relative abundance and associations with plasma aldosterone concentration. Am. J. Hypertens. 25, 1030–1034.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Unigene, (2010). National Center for Biotechnology Information (NCBI). Available at: www.ncbi.nlm.nih.gov/unigene

UniProt Consortium, S. (2012). Reorganizing the protein space at the Universal Protein Resource (UniProt). Nucleic Acids Res. 40, 18.

CrossRef Full Text

Veldhuis, J. D., Keenan, D. M., and Pincus, S. M. (2008). Motivations and methods for analyzing pulsatile hormone secretion. Endocr. Rev. 29, 823–864.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wang, Y., Xiao, J., Suzek, T. O., Zhang, J., Wang, J., Zhou, Z., et al. (2012). PubChem’s BioAssay Database. Nucleic Acids Res. 40, 2.

CrossRef Full Text

White, K. P. Jr., Radke-Sharpe, N. F., Kaiser, D. L., and Owens, G. K. (1989). A dynamic model of angiotensin II infusion experiments. J. Biomed. Eng. 11, 63–71.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Willmann, S., Hohn, K., Edginton, A., Sevestre, M., Solodenko, J., Weiss, W., et al. (2007). Development of a physiology-based whole-body population model for assessing the influence of individual variability on the pharmacokinetics of drugs. J. Pharmacokinet. Pharmacodyn. 34, 401–431.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Willmann, S., Schmitt, W., Keldenich, J., and Dressman, J. B. (2003). A physiologic model for simulating gastrointestinal flow and drug absorption in rats. Pharm. Res. 20, 1766–1771.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Willmann, S., Schmitt, W., Keldenich, J., Lippert, J., and Dressman, J. B. (2004). A physiological model for the estimation of the fraction dose absorbed in humans. J. Med. Chem. 47, 4022–4031.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Willmann, S., Thelen, K., Becker, C., Dressman, J. B., and Lippert, J. (2010). Mechanism-based prediction of particle size-dependent dissolution and absorption: cilostazol pharmacokinetics in dogs. Eur. J. Pharm. Biopharm. 76, 83–94.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wishart, D. S., Knox, C., Guo, A. C., Eisner, R., Young, N., Gautam, B., et al. (2009). HMDB: a knowledgebase for the human metabolome. Nucleic Acids Res. 37, 25.

CrossRef Full Text

Zhou, X., Shang, D., Zhang, T., Li, L., Zhou, T., and Lu, W. (2012). Modeling of angiotensin II-angiotensin-(1-7) counterbalance in disease progression in spontaneously hypertensive rats treated with/without perindopril. Pharmacol. Res. 9, 9.

Keywords: physiologically based pharmacokinetic model, cardiovascular, renin-angiotensin-aldosterone system, enalapril, enalaprilat

Citation: Claassen K, Willmann S, Eissing T, Preusser T and Block M (2013) A detailed physiologically based model to simulate the pharmacokinetics and hormonal pharmacodynamics of enalapril on the circulating endocrine renin-angiotensin-aldosterone system. Front. Physio. 4:4. doi: 10.3389/fphys.2013.00004

Received: 29 July 2012; Accepted: 08 January 2013;
Published online: 08 February 2013.

Edited by:

Robert Hester, University of Mississippi, USA

Reviewed by:

Radu Iliescu, University of Medicine and Pharmacy Gr. T. Popa “Iasi,” Romania
Marek Matejak, Charles University, Czech Republic
Patrick Hannaert, Institut National de la Santé et de la Recherche Médicale, France

Copyright: © 2013 Claassen, Willmann, Eissing, Preusser and Block. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.

*Correspondence: Michael Block, Bayer Technology Services GmbH, TD-ET-SB, Building 9115, 51368 Leverkusen, Germany. e-mail: michael.block@bayer.com