Breakthrough: a first

Scientific Reports volume 13, Article number: 3300 (2023) Cite this article

886 Accesses

1 Altmetric

Metrics details

The renin–angiotensin–aldosterone-systems (RAAS) play a central role in the pathophysiology of congestive heart failure (CHF), justifying the use of angiotensin converting enzyme inhibitors (ACEi) in dogs and humans with cardiac diseases. Seminal studies in canine CHF had suggested that the pharmacological action of benazepril was relatively independent of doses greater than 0.25 mg/kg P.O, thereby providing a rationale for the European labeled dose of benazepril in dogs with CHF. However, most of these earlier studies relied on measures of ACE activity, a sub-optimal endpoint to characterize the effect of ACEi on the RAAS. The objectives of this study were (i) to expand on previous mathematical modeling efforts of the dose-exposure–response relationship of benazepril on biomarkers of the RAAS which are relevant to CHF pathophysiology and disease prognosis; and (ii) to develop a software implementation capable of simulating clinical trials in benazepril in dogs bedside dose optimization. Our results suggest that 0.5 mg/kg PO q12h of benazepril produces the most robust reduction in angiotensin II and upregulation of RAAS alternative pathway biomarkers. This model will eventually be expanded to include relevant clinical endpoints, which will be evaluated in an upcoming prospective trial in canine patients with CHF.

Although the exact pathophysiology of the heart diseases underlying congestive heart failure (CHF) differ between man and his best friend, overactivation of the renin–angiotensin–aldosterone system (RAAS) plays a key role in the pathogenesis and development of CHF in both humans and dogs. To reduce RAAS activation, there is a substantial history of using ACEis, such as benazepril, to treat CHF in both species1,2,3. This makes the use of benazepril to treat CHF in canines and humans an excellent case study for applying the One Health Initiative paradigm. This paradigm recognizes that accumulating data on the effect of therapeutics on CHF in canines has the potential to benefit therapeutic management of CHF in humans and vice versa4.

The RAAS is a neurohormonal compensatory system which primarily manages blood volume and pressure by modulating electrolyte transport and vascular tone. The contemporary model of RAAS activation has two main components. The classical RAAS pathway refers to the peptide cascade from angiotensinogen to angiotensin I (AngI), and then from AngI to angiotensin II (AngII). These enzymatic reactions are catalyzed by renin and ACE, respectively, and ultimately lead to increased aldosterone (ALD) production (see Fig. 1;4). Short-term physiologic consequences of classical RAAS activation include vasoconstriction, renal sodium and water retention, and increased blood pressure. Long-term physiological consequences include fluid overload, increased cardiac afterload, and myocardial and vascular fibrosis5,6,7,8,9. Essentially, chronic long-term classical RAAS activation both contributes to, and is stimulated by, the development of CHF10,11, while classical RAAS pathway downregulation has been associated with improved long-term prognosis in CHF9,12,13,14,15. The alternative RAAS pathway acts as a counterregulatory mechanism against classical pathway activation. Activation of the alternative RAAS pathway is characterized by catalysis of AngII to angiotensin (1–7) (i.e. Ang(1–7)) by the enzyme ACE2. In turn, Ang(1–7) activates Mas receptors leading to vasodilatation, diuresis, and natriuresis16. Via this physiological effect, chronic alternative RAAS activation in CHF has been associated with reduced risk of heart failure in patients with reduced and preserved ejection fraction. An ideal therapeutic drug candidate for CHF would therefore modulate both pathways at once, downregulating the activity of the classical RAAS while preserving or upregulating the alternative RAAS pathway13. However, little is known about the effect of benazepril on the alternative RAAS in either humans or dogs.

Biological arms of the RAAS. RAAS activation is thought of as having two main pathways which act as counterregulatory mechanisms for one another. The classical RAAS pathway (in red orange) refers to the peptide cascade from angiotensin I (Ang I) to angiotensin II (Ang II) via ACE. This stimulates aldosterone production which then activates AT1 receptors (AT1R). Physiologic consequences of classical RAAS activation, including vasoconstriction, hypertrophy, and fibrosis, typical worsen congestive heart failure (CHF). Benazepril inhibits ACE, therefore activating the alternative RAAS pathway (in green). Activation of the alternative RAAS pathway is characterized by catalysis of Ang II to Ang1-7 by the enzyme ACE2. In turn, Ang1-7 activates Mas receptors leading to vasodilatation, diuresis, and natriuresis. These effects are protective against CHF. Our goal is to use mathematical modeling to determine a dosage which both reduces classical RAAS pathway activation and stimulates alternative RAAS pathway activation. This hypothetical dosage would maximize CHF-protective effects of benazepril.

Benazepril hydrochloride is a non-sulfhydryl ACEi commonly used for the management of CHF in both humans and dogs. Like other ACEi, benazepril is a prodrug that is rapidly converted through hydrolysis to its active benazeprilat by esterases, mainly in the liver17. Although frequently prescribed, the recommended dosage range of benazepril is quite broad and there is no clear consensus on the ideal dose to be used in patients with CHF. In humans, benazepril is typically prescribed for hypertension at an initial dose of 2.5–10 mg per day and up titrated to 20 or 40 mg per day, administered either once or twice daily (q24h or q12h) which is roughly equivalent to 0.5 mg/kg q12h for a 60 kg adult18. In dogs, the labeled dose of benazepril in the EU is 0.25–1.0 mg/kg PO q24h, whereas ACVIM veterinary consensus statements recommend a dose of 0.5 mg/kg PO q12h19. Pharmacokinetic (PK) and pharmacodynamic (PD) studies comparing various doses of benazepril in healthy dogs have not provided consistent recommendations to date. The study that was used for registration of benazepril in the EU showed that a single PO dose of benazepril effectively suppressed ACE activity for up to 24 h, and that ACE inhibition in plasma was independent of dosage ≥ 0.25 mg/kg2. However, subsequent reanalysis of these data using mathematical modeling suggested that q12h dosing (as opposed to q24h dosing) would achieve greater inhibition of ACE with the same q24h total dose20. Furthermore, a different study of single dose enalapril and benazepril at a dosage of 0.5 mg/kg indicated a much shorter duration of effect, with ACE suppression lasting < 12 h21, and a recent retrospective study in dogs with valvular heart disease suggested improved outcomes with q12h dosage22.

There are several reasons why developing consistent recommendations for dosing of ACEi has proven challenging in veterinary medicine. Historically, ACE activity was used as a surrogate for RAAS activity. Recently, however, ACE activity has been shown to be an inefficient measure of RAAS activation. Numerous studies in humans and dogs have shown a lack of correlation between circulating ACE activity and Ang II concentrations4,23. A second challenge in developing scheduling recommendations is the significant chronobiological modulation of the RAAS. Previous experimental models of RAAS activation failed to consider the chronobiology of the RAAS, while contemporary research has shown that biomarkers of the renin pathway are subject to circadian variations in dogs4,23,24. Finally, existing PKPD studies on the effect of various ACEi have not consistently sampled biomarkers of alternative RAAS activation in addition to biomarkers of classical RAAS activation.

Overall, although the effects of ACEi, such as benazepril, on ACE activity have been fairly well characterized, and the benefit of ACE inhibition in CHF has been definitively established in several clinical trials in both humans and dogs (0.25 to 1.0 mg/kg q12h-q24h), little is known about the effect of benazepril on the alternative RAAS pathway in either species. Understanding the dose-dependent effects of benazepril on biomarkers of both the classical and alternative RAAS pathways in dogs would allow exploration of benazepril dosages that produce a downregulation of the classical RAAS while preserving, or upregulating, the alternative RAAS. This would translate into an optimization of the clinical benefit. Accumulating data that inform such a nuanced approach to dose optimization in dogs would provide valuable translational information for similar dose optimization of ACEi in humans. To model and predict the dose-dependent effects of benazepril on the classical and alternative arm of the RAAS, we aimed to build a nonlinear mixed-effects (NLME) model of benazepril PKPD. NLME modeling of benazepril PKPD had already previously been shown to be an efficient method for describing the effect of benazepril on the classical RAAS in canines and is a well-accepted framework for building PKPD models11.

To produce data for this modeling and simulation effort, nine healthy beagles were intensively sampled while administered benazepril at various dosages and frequencies. After producing the data, our objective was to use a quantitative-systems pharmacology (QSP) model to characterize the PKPD relationship of benazepril(at) on biomarkers of the RAAS which are relevant to CHF pathophysiology and associated with morbidity/mortality {angiotensins I, II, III, IV, (1–7)}. QSP modeling is a subgroup of PKPD models which seeks to describe the behavior of a pharmaceutical in terms of the biology of its mechanism of action. After developing and calibrating the model, we further developed a software implementation of the benazeprilat-RAAS QSP model, which is capable of rapidly simulating the effect of benazepril HCL at various doses in a larger population of virtual dogs. By developing an easy-to-use simulation interface for our model, the objective of this work was to make a first prediction of the optimal dose/time of benazepril administration in dogs in support of future investigations in patients with CHF.

All study dogs received all oral doses of benazepril as intended. Dogs were monitored for adverse effects associated with benazepril labeling as well those associated with general animal welfare e.g., vomiting, diarrhea, inappetence, weakness/hypotension, fatigue, incoordination, hypercreatininemia. No adverse effects were observed in the animals during the course of the study, and serial complete blood counts and chemistry panels performed showed no evidence of hematologic or biochemical abnormalities from benazepril dosing.

Data were collated and standardized for mathematical modeling as instructed in the Monolix documentation25. Except for standardization of units as molar amounts and concentrations, raw data were left untransformed. Doses were transformed using the molecular weight of benazepril HCl, while concentrations were transformed using the molecular weight of the active metabolite–benazeprilat. Data were reviewed for bulk trends and data quality both before, during, and after mathematical modeling.

Data below the lower limit of quantification (LLOQ) were modeled by adding to the likelihood function a term describing the probability that the true observation lies between zero and the LLOQ, which is equivalent to the M3 method implemented in the NONMEM (Non-linear Mixed Effects Modeling) software.

The log10 time-course of benazeprilat as well as the relevant RAAS biomarkers are reproduced in Fig. 2. Of note, there was some background experimental noise in pharmacodynamics of some biomarkers that ultimately reduced model prediction quality. The noise was most prominent in the biomarker Angiotensin III (2–8) (i.e. AngIII), where the order of the limit of quantification (2.5 pmol/L) was approximately half of the measurement at the 3rd quartile (5.1 pmol/L). Suspected outliers were flagged and tested as model covariates for statistical significance. However, none of the flagged data points were determined to be significant enough outliers to exclude from model building.

RAAS biomarkers pharmacodynamics. An overview of the plasma time-course of several RAAS biomarkers as well as benazepril's active metabolite, benazeprilat. Each subject's time-course is indicated with a red line and points. The golden curve is the mean time-course value.

Following is a summary of the model building process. The base empirical version of the full model was largely an adaption of the benazeprilat PKPD model reported in Mochel et al.23. In all, over 100 different structural modifications were tested, starting from the empirical base model, to produce our final QSP model. To simplify results reporting, the most important modifications tested are summarized in the following two sections. Despite the division of sections into PK and PD, after building a base model to work from, all model fits were performed on the full PKPD dataset.

The PK portion of the base model was a conventional 2-compartment mammillary model with saturable exchange between the central and peripheral compartments. Building on this initial model, several modifications of the PK structure were evaluated. First, several standard compartment variations were tested i.e., using 1-, 2-, or 3-compartment disposition functions. Overall, a 2-compartment PK model outperformed the other candidate models based on the precision of individual parameters and overall quality of fit. Second, non-specific (low affinity, high capacity) binding of benazeprilat to plasma proteins was represented by a 3rd compartment within the central compartment, i.e., representing the free circulation of benazeprilat. The volume of the non-specific binding compartment (Vns) is a representation of the relative binding capacity of benazeprilat which is distributed in plasma but is not freely circulating or interacting with ACE. Therefore, the total amount of measurable benazeprilat in plasma is a combination of the amount non-specifically bound to plasma proteins (Ins) (low affinity, high capacity), the amount specifically bound to ACE (high affinity, low capacity) and the amount of benazeprilat in free circulation (Ifree)4,20. The variable I was chosen to represent benazeprilat as it inhibits ACE activity.

Zero-, first-, mixed-, and sequential absorption structures were tested to model drug absorption from the depot compartment (i.e., intestinal lumen). A model largely equivalent to sequential absorption, but made to be continuous, was found to outperform other competitive models. This model substructure uses a series of 1-order absorptions but can be seen as a continuous analog to a sequential 0-/1-order absorption.

In this model, the first depot compartment for benazepril after oral administration was called 1abs. First-order absorption either occurred immediately at rate ka1 to the compartment of freely circulating benazeprilat (fr), or absorption was delayed by an absorption rate ka through an intercompartment that was pre-circulatory (pr). As is typical, quantity of benazeprilat passed between compartments is called Im, n, for inhibitor, where indexes m and n represent origin and destination compartments, respectively. Fbio represents the total bioavailability (Eq. 1). Doses are administered in benazepril hcl, but are measured as the metabolite–benazeprilat. To reduce complexity in modeling, but preserve absorption and benazepril to benazeprilat conversion variance, all bioavailable benazepril is treated as benazeprilat in the model. Fbio, or total bioavailability, is estimated in the model purely to preserve this variance and to reduce numerical instability in estimation. However, without IV data, the final estimated Fbio does not have a firm pharmacological interpretation.

In summary, the final mammillary model without ACE binding kinetics (Eq. 2) was a 2-compartment PK model with nonspecific protein binding represented by a 3rd compartment (Ins), and a continuous analog to sequential 0-/1-order absorption from some depot compartment.

Rate of exchange between compartments were governed by rates kf, g where f and g (f ≠ g) were each one of either free circulation (fr), tissue (ts), or non-specifically bound in circulation (ns). Residual error was best modeled using a normal-proportional error function (Eq. 6). The only exception were rates of elimination which were written as kCl, d where clearance represented that the parameter was derived from clearance and d was the compartment of origin.

Benazeprilat primary mechanism of action is inhibiting ACE to prevent the catalysis of AngI into AngII. To account for this mechanism, a logistic saturation model was first implemented. However, the superior model for predicting benazeprilat ACE inhibition was found to be the differential Michaelis–Menten model of catalysis inhibition with ACE being the enzyme (E), benazeprilat the inhibitor (I), AngI being the substrate (S) and AngII being the product (P) (Eq. 3). The distribution of ACE across tissue (ts) and free circulation (fr) was also considered. The nomenclature used throughout Eq. (3) is consistent with previous descriptions of the Michaelis–Menten model26.

Two-compartment mammillary models governed the kinetics of biomarkers AngI, AngII, and Ang(1–7). The amount in these compartments were respectively represented by S (substrate), P (product), and Ang(1–7). The two compartments for these angiotensins were called free circulation (fr) and tissue (ts). Conversion steps in the classical and alternative RAAS pathways were modeled through a series of catalytic steps, as previously described27.

At last, the function fCT(t) governs the effect of chronobiology on the production rate of the substrate (rs). fct(t) is a scaled cosine function where the wavelength (or period) is matched to 24 h, the relative maximum amplitude is the scalar PRA (peak renin amplitude), and the scale of that amplitude is governed by δ24hr. Chronobiology is herein only modeled relative to AngI production (Eq. 4).

The catalyses of AngII to AngIII, and AngIII to AngIV were modeled via a series of catalytic conversion models (Eq. 5). Cleavages of AngII to AngIII, and AngIII to IV, are primarily performed by renally-bound aminopeptidases A and N, respectively28,29,30. Vfree was subdivided into two circulatory system volumes of distribution; a small renal volume (Vrn) and a larger plasma volume (Vpl). All catabolism of AngII to AngIII, and AngIII to AngIV were linked to the renal volume as this is where aminopeptidases A and N are physiologically located.

All analytes residuals were best described by proportional error models (Eq. 6), with the concentration of a given biomarker scaled by ε. ε is a normal distribution distributed with standard deviation b, i.e. ε ~ N(0, b).

ANOVA tests on covariates indicated that model performance would not be significantly improved by the inclusion of any covariate effects. The full model written in Mlxtran is available in the supplemental files, and a model diagram detailing the full structure is reproduced in Fig. 3. In S1 Table, the reader can find a detailed description of all mathematical symbols defined in the "Results" section.

Detailed model diagram. Detailed diagram of the final model structure. Benazeprilat pharmacokinetics were modeled using a 2-compartment model with a mix of 1-order and 1-order delayed by 1-order transfer absorption from the depot compartment. Both volumes of distribution, free and tissue, were modeled with a fixed amount of ACE with which Benazeprilat could act on. Non-specific binding affected the free circulation compartment. A series of direct response models were used to describe the transformation of angiotensin I into its various metabolites. The free volume of distribution was subdivided into plasma and kidney volumes for Ang III and Ang IV. A Michaelis–Menten kinetic model of inhibitor, substrate, and enzyme interaction was used to describe the competitive inhibition of ACE by benazeprilat. k3 and k-3 were the parameters governing ACE-benazeprilat (enzyme-inhibitor) association and dissociation, while k1 and k-1 determined the rate of ACE-angiotensin (enzyme–substrate) association and dissociation. k2 controlled the production rate of angiotensin II from angiotensin I via ACE. An independent clearance for each metabolite as well as benazeprilat controlled the rate of removal of various molecules from the plasma.

Inspection of the SAEM search and a sensitivity analysis on initial parameter values revealed a stable and precise search for all parameter estimates. The final selected model had high precision in parameter estimates as evaluated via RSE (majority of estimates < 35%). A summary of model parameter estimates, including typical value, RSE (%) and inter-individual variability (IIV) can be found in Table 1.

Inspection of goodness-of-fit summary plots (Figs. 4, 5, 6) indicate that benazeprilat predictions from the model are largely in line with experimental measurements. Importantly, the final PKPD model, which enabled the simultaneous fit of all angiotensins, was found to characterize the time-varying changes of the both the classical and alternative arm of the RAAS satisfactorily, as shown by the standard goodness-of-fit diagnostics of observations vs predictions (Fig. 4), the individual predictions (Fig. 5), and the simulation-based validation diagnostics (i.e., NPDEs, Fig. 6).

Observations vs predictions. The observations plotted against the predictions for all metabolites and drug concentration data. This gives a complete picture of model performance. The golden line is the LOESS curve showing the correlation between observations and predictions. The black line plotted diagonally represents an ideal model performance with no misspecification. The general agreement between LOESS and idealized curve indicates that there is little misspecification in model structure.

Sample of individual predictions. A sample of individual observations vs predictions randomly sampled from the concentration and metabolite data. The general agreement between plasma concentration time-course and individual predictions indicates that the model reproduces the observations with high accuracy.

Normalized prediction distribution errors. Normalized prediction distribution errors (NPDEs) are an analog to residuals used for diagnosing both model structural misspecifications as well as the performance of the residual error model. The distribution of a well-specified model is normal, ideally. Bands represent the 90% prediction band for the 95th, 50th, and 5th percentiles, respectively. Curves are the observed percentiles for the 95th, 50th, and 5th percentiles, respectively. Data are binned at regular intervals to derive these average trends. For each binning range, if the structural model fits the data well, observed percentiles will be symmetrically distributed across a 50th percentile curve which falls within the 50th percentile band. If the error model is well specified, observed percentiles will fall within prediction bands. Any misspecification is ideally random. The model appears to underpredict angiotensin (1–7) for small measured values slightly, but otherwise there is high agreement between model and data.

There are three primary views in this application (Fig. 7). In all views, the time of first dose of benazepril is specified in a 24-h clock format. On the left-hand side of the application is a menu for specifying dosing, parameters of the simulation, and modalities for calculating the area under the effect curve (AUEC) that quantifies the effect of the active benazeprilat on the RAAS at various doses vs. placebo control. Note that the application menu can be hidden to increase the size of the plotting area.

Application views. Our application provides a user-friendly way to apply our model of RAAS response to various administration schedules of benazepril. A collapsible left-hand widget allows the user to specify the simulation. The application has 3 separate panels. In the first panel, a single dosage scheme can be applied to a large simulated population of animals. Then the prediction distribution of simulated patient responses is plotted for study. In the second panel, the user can make a comparison between several proposed administration schedules. The plots produced in this panel depict the median time-course of the various metabolites in response to the proposed schedules. The user also has access to an x-axis zoom and AUC comparison summaries on the right-hand side. The final panel is simply documentation of the simulation engine code and tips for usage.

The left-hand menu is split into three tabs which allow the user to define parameters of the simulation. The dosing tab permits the user to define the dosing schedule in terms of time of first dose, number of doses, size of dosage, and interdose interval. The simulation parameters tab gives the user access to the timescale of the simulation, the fineness of the grid used for simulation, and the sample size used to calculate the median and prediction intervals of the simulated PKPD. Finally, the AUEC tab provides a means to compare pharmacodynamic effects between competing dosing scenarios by defining a time period for which to calculate AUEC estimates.

The first tab gives the user tools for analyzing the distribution of responses after a single schedule of benazepril. The distribution is specified in terms of median effect (blue line), median effect of placebo using same simulated individuals (dashed black line), and 90% prediction interval (blue bands) in steps of 10% i.e., 5% to 95%, 15% to 85%, etc. The AUEC of treatment vs. placebo can be compared for the timespan between the dashed vertical lines. The percent difference between those two AUECs is documented in the hovering label.

The second panel allows the user to compare up to four competitive dosing schedules to placebo. In this panel, the user can see the median response (key at bottom) and the placebo effect (dashed black line), but not the distribution of responses. On the right-hand side of the page, the user can compare the percent difference from placebo in the RAAS components modeled in this study by paging through the various data tables. These comparisons are percentages relative to placebo.

The last view is simply documentation on the model that powers the simulation engine and general recommendations for using the software. It also provides a brief summary of the design of the application and gives a full reproduction of the R code that makes up the model. In this view, the application menu also gives a brief summary of user warnings.

As a final consideration in our study, we directly compared four dosage scenarios in our simulation engine: 0.25 mg/kg q24h in the AM, 0.25 mg/kg q24h in the PM, 0.25 mg/kg q12h, and 0.5 mg/kg q12h. In our simulation application, we set the engine to compare the median AUEC of 500 dogs (matched between virtual trials), at a sampling rate of 500 times over a period of 25 virtual days. The median AUEC comparison was made on day 20 over a period of 24 h. The long virtual time of simulation assured the simulated dogs reached steady state PD of benazeprilat. For 0.25 mg/kg q24h, independent of the time of dosing, we saw an approximate 55% decrease compared to placebo for AngII and a 95% increase in Ang(1–7). With a schedule of 0.5 mg/kg q12h, we saw an approximate 80% decrease versus placebo for AngII and 135% increase in Ang(1–7). Overall, a greater daily biomarker variance was observed with q24h (vs. q12h) dosing. Summary results are tabulated in Table 2 while median time-courses are plotted in Fig. 8.

Simulation summary. Four administration schedules are compared in this simulation scenario: (1) 0.25 mg/kg every 24 h at 8 a.m.; (2) 0.25 mg/kg every 24 h at 8 p.m.; (3) 0.25 mg/kg twice a day at 8 a.m. and 8 p.m.; (4) 0.5 mg/kg twice a day at 8 a.m. and 8 p.m. 500 individuals were simulated for each scenario (for a total of 2500 individuals, with placebo). Curves are the median time-course of the molecule from this simulated population. Comparisons between schedules are made by calculating the percent difference between the median AUC and placebo (of each schedule). Summary results are tabulated in Table 2.

In both canines and humans, classical RAAS overactivation plays a key role in the pathogenesis and development of CHF4. To modulate classical RAAS overactivation in CHF, there is a substantial history of using ACEi such as benazepril in both species1,2,3. Improving our understanding of the effect of ACEi on CHF in canines has the potential to benefit therapeutic management of CHF in humans and vice versa4. Activation of the alternative RAAS pathway is characterized by catalysis of AngII to Ang(1–7) by the enzyme ACE2. In turn, Ang(1–7) activates Mas receptors (Esteban PloS One 2009). In direct contrast to the effect of classical pathway overactivation on CHF, activation of the alternative pathway is associated with improved clinical outcomes via reduced and preserved ejection fraction and reduced risk of heart failure.

An ideal therapeutic for CHF would modulate both the classical and alternative RAAS pathways, downregulating the activity of the classical RAAS pathway while preserving or upregulating the alternative RAAS pathway13. Although the effect of ACEi on ACE activity and AngII has been fairly well-characterized both in the veterinary and human literature, little is known about the effect of ACEi on the alternative arm of the RAAS. Therefore, any further dose optimization of ACEi is dependent on studying the effect of this therapeutic drug on the alternative pathway. Characterizing this effect in dogs essentially furthers our understanding of how to optimize therapeutic management of canine CHF while providing valuable insight into the molecular effects of ACEi for translation to human CHF31.

Additionally, in humans and in dogs, the recommended dosage range for benazepril is quite broad and there is no clear consensus on the ideal dose in patients with CHF. PKPD studies comparing various doses of benazepril in healthy dogs have not provided consistent recommendations. In King et al. a single PO dose of benazepril effectively suppressed ACE activity for up to 24 h, and that ACE inhibition in plasma was independent of dosage ≥ 0.25 mg/kg17. Subsequent reanalysis of these data using PK modeling suggested that q12h dosing would achieve greater inhibition of ACE with the same total q24h dose20. Later, Hamlin and Nakayama found a single dose of benazepril at 0.5 mg/kg suppressed ACE for < 12 h21. Lastly, a recent retrospective study in dogs with valvular heart disease suggested improved outcomes with q12h dosage22.

In this study, we have attempted to address knowledge gaps in optimal benazepril dosing by describing the dose-dependent effects of benazepril on biomarkers of both the classical and alternative RAAS pathways in dogs. The solution we implemented to address these gaps in knowledge was to build a modeling and simulation platform of the effect of benazepril on both arms of the RAAS. This simulation engine allows for exploration of benazepril dosages that produce both a substantial downregulation of the classical RAAS while preserving or upregulating the alternative RAAS. This allows for side-by-side comparison of several dosage schemes using virtual clinical trials and ultimately optimization of clinical benefits.

Several metrics were used in assessing model performance. Population-level goodness-of-fit diagnostic plots like observations vs. predictions and NPDEs indicated that structural misspecification was very low. Focusing on individual predictions, we see a dynamic model capable of fitting complex individual variations without being overly fit to noise and spurious trends in the data. Importantly, for such a large set of parameters, our precision of parameter estimates was very high. Some parameters had to be fixed to exploratory values to achieve that precision, but this is an expected outcome when using enzyme kinematic models. Of note, although the precision of these parameters was high, they should not be overly interpreted without experimental verification.

Using our simulation engine, we chose to compare several reasonable dosing schedules, including 0.25 mg/kg PO q24h and 0.5 mg/kg PO q12h—covering the range of dosing schedules most commonly used in the EU and US. Comparisons between dosing schedules were made based on area under the effect curve of the biomarker response (AUEC) relative to placebo.

Chronobiology played a modest role in the scheduling of benazepril in this study. While using the simulation engine to explore various dosages, evening dosing appears to produce the lowest variance in the classical RAAS, morning dosing appears to produce the lowest variance in AngIII and AngIV PD. Morning and evening administration both produce the same relative AUC improvement over placebo across the biomarker panel.

Of the schedules compared, the most robust downregulation of classical RAAS biomarkers and upregulation of alternative RAAS biomarkers was observed with the 0.5 mg/kg PO q12h dose of benazepril. However, this should be put in context because the improvement of 0.5 mg/kg q12h over 0.25 mg/kg q12h is 12.5%, 13.3%, 14.2%, 13.1% and 14.3% for Ang(1–7), AngI(1–10), AngII(2–8), AngIII, and AngIV, respectively. Those respective changes in classical and alternative RAAS biomarker concentrations will need to be linked to clinical outcomes in order to inform a clinical dosage selection.

An interesting finding from our simulations is that there is generally high agreement between q24h and q12h dosing of benazepril, as long as the total dosage per day is kept constant. For example, 0.5 mg/kg q24h and 0.25 mg/kg q12h produce similar pharmacodynamic effects on the RAAS, although q12h dosing led to fewer fluctuations of the angiotensins in plasma compared to q24h dosing. Noteworthily, although gender-related differences in heart failure biomarkers have been reported in human patients32, our experimental study failed to demonstrate significant differences in RAAS pharmacodynamics between male and female dogs. However, this should be interpreted cautiously in light of the small sample size in our study.

To the best of the authors’ knowledge, this is the first description of a simulation engine designed to optimize dosages of therapeutic drugs in veterinary medicine. Although experimental validation is necessary for bedside application, the model structure readily lends itself to such validation. Many parameters are directly measurable because they have meaningful pharmacological interpretations. For example, the Michaelis–Menten kinetic parameters which govern ACE inhibition could potentially be measured independently—possibly in vitro. Additionally, the kinetics of each individual molecule can be measured separately in IV-PK time-course studies in vivo. Furthermore, by using a combination of in vitro studies for defining enzyme kinematic parameters, in vivo studies of individual metabolite kinetics, literature values for baseline parameters, and allometric scaling of mammillary compartment parameters, this model can be easily adapted into a bedside tool in both canines and humans. In doing so, other physiological variables such as gender and ethnicity, which are significant covariates of treatment response to ACEi in human patients with CHF would need to be considered33,34. Yet, another important variable that needs including in an attempt to scale our model to humans, is the fact that many endogenous peptides as well as albumin are known to suppress ACE activity in vivo35,36. Indeed, early reports from Ryan et al.37 suggest that many endogenous small molecules (< 10 kDa) would already contribute significantly to ACE inhibition. Whether these observations hold true in dogs remains to be shown. Interestingly though, endogenous ACE inhibition has also been reported in donkey, goat and bovine sera35.

There are several practical limitations in this study worth recognizing to guide future model refinement. First, our results were derived from an experimental model of RAAS activation in response to a single dose of benazepril, rather than repeated dosing in a clinical trial with canine CHF patients. Also, the sample size of the study was quite limited compared to a mature clinical trial design. Consequently, our simulation engine does not consider multiple physiological factors that could potentially influence the effect of benazepril in a patient population, such as the impact of age, breed, and other co-medications that are known to modulate the RAAS, such as furosemide38 and spironolactone39. This information will be provided in a prospective clinical trial in canine patients with CHF to adjust our model-based predictions. Lastly, the therapeutic window was not considered when making comparisons between dosing schedules. Specifically, there are high doses of benazepril the user can specify in our simulation engine (max. 2 mg/kg q6h) that have not been tested experimentally; although safety up to 1 mg/kg q24h has been established in previous studies22.

This extensive QSP model of the RAAS in response to benazepril and the development of that model into a tool for bedside optimization of benazepril for CHF and similar diseases is highly novel. By developing an easy-to-use simulation interface for our model, we are now able to make a first prediction of the optimal dose/time of benazepril administration in dogs in support of future investigations in patients with CHF. Beyond the research presented in this manuscript with our tool, simulation tools can continuously expand the impact of scientific research by being used to test new hypotheses surrounding dose optimization and being improved when new data becomes available to refine parameter estimates. The data we measured in an experimental model of RAAS activation has yet to be linked directly to clinical outcomes in CHF, so there is room to expand the model into disease outcomes with a link function. This will be necessary for ultimate application in dosage selection. The most important extension is to experimentally validate a relevant selection of simulations; this option is currently being explored. This model-based approach is now supporting the design of an upcoming prospective multicenter clinical trial in canine patients with CHF to confirm findings from our simulator and refine our model-based predictions with actual clinical information. This clinical trial will help to confirm if the difference observed in PKPD between different dosages is translated into clinical benefits in dogs with naturally occurring CHF.

Experimental procedures were performed in accordance with relevant guidelines and regulations at Iowa State University. The study was approved by the Institutional Animal Care and Use Committee at Iowa State University under Protocol 19-344. Authors complied with ARRIVE guidelines in the completion of this study.

Nine purpose-bred laboratory beagles (5 castrated males and 4 spayed females), 40–42 months old, weighing 9.0–13.5 kg were randomized based on body weight and sex into three oral dosing groups of benazepril. Systemic and cardiovascular health of all dogs was confirmed prior to the study with physical examination, routine laboratory screening (complete blood count, serum biochemical analysis), blood pressure measurement, and echocardiography.

Study dogs were housed in the Laboratory Animal Resources unit at the Iowa State University College of Veterinary Medicine. Dogs were acclimatized to the facility for > 1 month prior to the experiment. Dogs were pair-housed in adjoining pens (approximately 2 m2 per dog or 4 m2 per pair) on elevated rubber-coated grated flooring. Housing conditions were standardized with ambient temperatures of 18 °C, a 12-h light cycle (07:00 to 19:00), and access to water ad libitum. During intensive sampling days only (D1, D18, and D35), dogs were separated into single housing units and water consumption was quantified every 8 h for the 24-h period.

After sample collection on baseline sampling days (D-5, D12, and D29), dogs were offered a low-sodium diet (Hill's Prescription Diet h/d, 17 mg sodium per 100 kcal) at 23:00 q24h for 5 days to attain a steady activation of RAAS4,23,24. After data collection on D2, D35, and D36, dogs began a 10-day wash-out period between cycles during which they were offered their standard diet (Royal Canin Beagle Adult, 110 mg sodium per 100 kcal) q24h at 09:00. Volume of low-sodium diet was calculated so that the dogs received the same caloric intake throughout the study.

This 35-day prospective study was divided into three periods with three different benazepril dosing groups: (A) 0.125 mg/kg q12hr PO, (B) 0.25 mg/kg q12hr PO, and (C) 0.5 mg/kg q24hr PO. All dogs received all treatments using a partial crossover (ABC/BCA/CAB) design. Dogs were sampled in the same order at each time point and exact time of sampling was recorded. Blood sample collection was divided into baseline sampling days (D-5, D12, and D29), sparse sampling days (days 0, 17, and 34), and intensive sampling days (D1, D18, and D35). Baseline and sparse sampling occurred at 07:00.

On intensive sampling days (D1, D18, and D35), blood was collected starting at 07:00 (0 h, immediately before oral benazepril dosing) and repeated at + 0.5, 1, 2, 4, 8, 12, 12.5, 13, 14, 16, 20, and 24 h post-dosing. Benazepril (NELIO® 5 mg chewable tablets, Ceva Sante Animale) was administered on intensive sampling days following the 0-h blood sampling (all dose groups) and the 12-h sampling (q12hr dose groups only). The benazepril dose was calculated to the nearest 1.25 mg increment.

Venous blood samples were collected from an external jugular or cephalic vein with a 1 inch, 20-gauge or 22-gauge needle attached to a 6 mL syringe. Dogs were kept and maintained in the same position (seated with neck extended) during blood collection. On intensive sampling days, approximately 4 mL of whole blood was collected at each time point with 2 mL transferred to an additive-free collection tube and 2 mL transferred to a lithium heparin tube containing 11.2 µL of dichlorvos prepared as a 6 mg/mL solution in acetonitrile. On baseline days, approximately 6 mL of whole blood was collected with 2 mL placed in an additive-free tube for RAAS analysis, 2 mL placed in an EDTA tube for complete blood count, and 2 mL placed in an additive-free tube for serum chemistry panel. On sparse sampling days, 2 mL of blood was collected and placed in an additive-free tube. All samples intended for pharmacokinetic or RAAS analysis were centrifuged for 15 min, after which plasma or serum was transferred into cryovials that were then stored at − 80 °C for later analysis. Samples for complete blood counts or serum chemistry panels were analyzed by the Iowa State Clinical Pathology Laboratory.

Plasma benazeprilat analysis was performed by the Iowa State University Analytical Chemistry Laboratory. Benazeprilat and benazeprilat-d5 analytical standards were purchased from Toronto Research Chemicals (Ontario, Canada). Benazepril and benazepril-d5 analytical standards were purchased from Cayman Chemical (Ann Arbor, Michigan, USA). Benazeprilat and benazeprilat-d5 stock standard solutions were prepared at 0.25 mg/mL in 2:1:1 acetonitrile:water:DMSO. The benazepril and benazepril-d5 stock solutions were prepared at 1 mg/mL in acetonitrile. Control beagle plasma was purchased from Equitech Bio (Kerrville, TX, USA). All solvents used for sample preparation and the chromatography portion of the analytical method were purchased from Fisher Scientific (Waltham MA, USA).

A sample volume of 150 µL was fortified with 15 µL of a benazeprilat-d5 solution at 0.1 ppm. Plasma samples were precipitated with 600 µL of acetonitrile containing 0.5% formic acid and vortexed by hand for several seconds. All samples were centrifuged at 10,000 rpm for 5 min. A 600 µL volume of each sample was transferred to a clean 2 mL flip-top tube. All flip-top tubes were placed in the CentriVap Concentrator system (Labconco Corp., Kansas City, MO, USA) and concentrated to dryness. All samples were reconstituted in 100 µL of 50:50 acetonitrile:water and centrifuged at 10,000 rpm for 5 min prior to LC–MS/MS analysis. All samples were analyzed using an injection volume of 2 µL.

A Vanquish Flex LC pump interfaced with a TSQ Altis mass spectrometer (Thermo Fisher Scientific, San Jose, CA, USA) were used for the analysis. The source conditions were as follows: spray voltage—3500 V, sheath gas—40.6 Arb, auxiliary gas—23 Arb, sweep gas—0.4 Arb, ion transfer tube temperature—325 °C, and vaporizer temperature—350 °C. The total run time of the method was 3 min. The resolution of Q1 and Q3 was 0.7 FWHM. The CID gas was set to 2 mTorr. The chromatographic peak width was 2 s. and the cycle time was 0.2 s. The mass spectrometer was operated in positive ion electrospray ionization mode. Data was acquired using a multiple reaction monitoring (MRM) method that selected for the benazeprilat ([M+H]+ 397.2) and benazeprilat-d5 ([M+H]+ 402.2) precursor ions.

The column used for the analysis was Hypersilgold Aq 50 × 2.1 mm, 1.9 µm (Thermo Fisher Scientific, Waltham, MA, USA). Mobile Phase A was water + 0.1% formic acid and Mobile Phase B was acetonitrile + 0.1% formic acid. The column oven temperature was set to 35 °C. The chromatography gradient was as follows: Start at 0% B and linear ramp to 100%B in 2.0 min, hold at 100% B for 0.4 min, drop to 0% B in 0.01 min, and hold at 0% B for 0.59 min. The flow rate of the method was 0.4 mL/min.

Equilibrium concentrations of Ang I, Ang II, Ang III, Ang IV, Ang 1–9, Ang 1–7, and Ang 1–5 were quantified in serum samples by liquid chromatography-mass spectrometry/mass spectroscopy performed at a commercial laboratory (Attoquant Diagnostics, Vienna, Austria), as previously described4,40. Briefly, samples were spiked with a stable isotope labelled internal standard for each angiotensin after ex vivo equilibration and analytes were extracted using C18-based solid- phase extraction. Extracts samples were analyzed using mass spectrometry analysis using a reversed- analytical column (Acquity UPLC C18, Waters) operating in line with a XEVO TQ-S triple quadrupole mass spectrometer (Waters Xevo TQ/S, Milford, MA) in multiple reaction monitoring mode. Internal standards were used to correct for analyte recovery across the sample preparation procedure in each individual sample. Analyte concentrations were calculated from integrated chromatograms considering the corresponding response factors determined in appropriate calibration curves in serum matrix, when integrated signals exceeded a signal-to-noise ratio of 10. The lower limits of quantification for the analytes in canine serum were 3 pmol/L (Ang I), 2 pmol/L (Ang II), 2,5 pmol/L (Ang III), 2 pmol/L (Ang IV), 2,5 pmol/L (Ang 1–7), and 2 pmol/L (Ang 1–5), respectively41,42.

Pharmacokinetic and pharmacodynamic data were imported into R 4.0.243 for model preparation. To maintain consistency across administration units and measurement units, all concentration data were converted to micromoles per liter. Benazepril HCl and benazeprilat molecular weights were obtained from PubChem44,45.

Historical control data from a previous study was used in lieu of a placebo group23. After referencing parameters to this baseline during fit, all placebo simulation response generated was equivalent to the baseline chronobiology function with parameters estimated from the simultaneous modeling of all biomarkers response.

The recorded data (yij) were imported into Monolix 20120 R1 (Lixoft, France) and used to estimate population parameters (μ) and variance via the stochastic approximation expectation maximization algorithm (SAEM)46. Individual parameters (ϕi) were determined via the modes of the individual posterior distributions—which were estimated using a Markov-Chain Monte-Carlo (MCMC) procedure. NLME models were written as previously described (Eq. 7)47,48.

Model predictions (F(ϕi, βi, tij)) for the ith individual at the jth timepoint—time (tij)—were parameterized using individual parameters and individual covariates (βi). The residuals were modeled as G(ϕi, tij) · εij, where G is an arithmetic combination of proportional and additive error distributions.

Individual parameters were modeled as a function of the population parameters, interindividual variability (ηi), and individual covariates via the interindividual variation function h(μ, ηi, βi). Individual variability was modeled with a normal distribution of mean 0, variance–covariance matrix Ω, and variance ω2. Typically, h(μ, ηi, βi) is a log-normal link function (Eq. 8) or a logit-normal link function (Eq. 9) in those cases where ϕi is bounded to be between 0 and 1.

The systems pharmacology model was built in two consecutive phases. First, a largely empirical PKPD model with minimum parameterization was built to capture basic biological variations in the data. In practice, this entailed fitting a basic 1-, 2-, or 3-compartment model to the pharmacokinetics of benazeprilat, and then linking the PK to the RAAS biomarkers concentrations via various indirect and direct response models. Whenever possible we opted for direct over indirect effects models, and fewer compartments, to reduce the number of total estimated parameters during model fit.

Then, in an iterative fashion, model components were replaced with more mechanistic structures. To do this, we modeled the cascade of peptides which define the alternative and classical RAAS pathways. We also tested whether we could expand the mathematical model to include important biological systems such as the clearance of angiotensins via the liver vs. kidneys, non-specific plasma binding, first-pass metabolism, and site-specific metabolism. Some components and parameters of the model structure were arbitrarily fixed to literature or exploratory values to preserve fidelity to relevant biological systems. For example, our model equations were rewritten so that the production of AngII was always one-to-one proportional with catalysis of AngI via the angiotensin converting enzyme. The final model was refined through various arithmetic simplifications and parameter search optimizations to improve precision of parameter estimates as much as possible without comprising fit to experimental data. The significance of bodyweight, sex, sodium intake, and benazepril dose on parameters estimates was further evaluated using the automated Pearson's correlation test and ANOVA method as implemented in Monolix 2020 R1.

Experimental data were collated and imported into the 2020R1 Monolix Suite for data exploration, model development and evaluation. Quality of fit was evaluated using standard goodness-of-fit diagnostics (e.g., observed vs. predictions, scatter plot of residuals), as well as numerical summaries of fit, such as the corrected Bayesian information criteria (BICc). Precision of parameter estimates were determined using residual standard error (RSE%).

Normalized prediction distribution errors (NPDEs) are recommended for evaluating model misspecification when study design is heterogenous with respect to dosing groups49. In this study, we had several different study groups under varying dosing schedules; therefore, NPDEs were chosen to determine quality of fit over more conventional visual predictive checks. In practice, NPDEs evaluate the percentage of the distribution of predictions at each mean prediction under the observation, thus forming a heterogenous uniform distribution. Therefore, an inverse cumulative distribution function is applied to each value to obtain a probability density function vs. population prediction.

All programming relevant to the simulation application was written in R v4.0.243. Models were translated from Mlxtran (Monolix Suite 2020R1) to the domain specific language, and R package, Odin v1.2.150 to simulate clinical trials. Odin provides an interface for the ODE solvers, and R packages, deSolve v1.3051 and dde v1.0.152. Odin was used because of its superior ability to solve large systems of ODEs over other R packages.

Finally, the application for simulating clinical trials was built in Shiny v1.6.0. Shiny is an R package that automatically generates HTML applications from R code. All Shiny applications are designed to (1) generate an HTML-based graphical user interface (GUI) that allows users to interact with the computer hosting the Shiny application, called the server, and (2) execute R code on the server based on the user's interactions with the GUI. To facilitate the use of our simulator, a user-friendly GUI was developed that allows to specify modalities for a clinical trial simulation in R (i.e., defining parameters such as dosage, dosing interval, size of trial) on a website server.

The data generated during this study are available from the corresponding author on reasonable request.

BENCH (BENazepril in Canine Heart disease) Study Group. The effect of benazepril on survival times and clinical signs of dogs with congestive heart failure: Results of a multicenter, prospective, randomized, double-blinded, placebo-controlled, long-term clinical trial. J. Vet. Cardiol. 1(1), 7–18 (1999).

Article Google Scholar

Woodfield, J. A. Controlled clinical evaluation of enalapril in dogs with heart failure: Results of the Cooperative Veterinary Enalapril Study Group The COVE Study Group. J. Vet. Intern. Med. 9(4), 243–252 (1995).

Article Google Scholar

Swedberg, K. Long-term survival in severe heart failure in patients treated with enalapril; ten year follow-up of CONSENSUS I. Eur. Heart J. 20(2), 136–139 (1999).

Article CAS PubMed Google Scholar

Mochel, J. P. et al. Pharmacokinetic/pharmacodynamic modeling of renin-angiotensin aldosterone biomarkers following angiotensin-converting enzyme (ACE) inhibition therapy with benazepril in dogs. Pharm. Res. 32(6), 1931–1946 (2015).

Article CAS PubMed Google Scholar

Seifarth, C., Trenkel, S., Schobel, H., Hahn, E. G. & Hensen, J. Influence of antihypertensive medication on aldosterone and renin concentration in the differential diagnosis of essential hypertension and primary aldosteronism: Antihypertensive medication and primary aldosteronism. Clin. Endocrinol. 57(4), 457–465 (2002).

Article CAS Google Scholar

Chobanian, A. V. et al. Seventh report of the joint national committee on prevention, detection, evaluation, and treatment of high blood pressure. Hypertension 42(6), 1206–1252 (2003).

Article CAS PubMed Google Scholar

Pacurari, M., Kafoury, R., Tchounwou, P. B. & Ndebele, K. The renin-angiotensin-aldosterone system in vascular inflammation and remodeling. Int. J. Inflamm. 2014, 1–13 (2014).

Article Google Scholar

Tummala, P. E. et al. Angiotensin II induces vascular cell adhesion molecule-1 expression in rat vasculature: A potential link between the renin-angiotensin system and atherosclerosis. Circulation 100(11), 1223–1229 (1999).

Article CAS PubMed Google Scholar

Qi, Y. et al. Moderate cardiac-selective overexpression of angiotensin II type 2 receptor protects cardiac functions from ischaemic injury: Cardioprotection by overexpression of AT 2 receptors. Exp. Physiol. 97(1), 89–101 (2012).

Article CAS PubMed Google Scholar

Regulski, M. et al. Chemistry and pharmacology of angiotensin-converting enzyme inhibitors. CPD. 21(13), 1764–1775 (2015).

Article CAS Google Scholar

Toutain, P. L. & Lefebvre, H. P. Pharmacokinetics and pharmacokinetic/pharmacodynamic relationships for angiotensin-converting enzyme inhibitors. J. Vet. Pharmacol. Ther. 27(6), 515–525 (2004).

Article CAS PubMed Google Scholar

Ames, M. K., Atkins, C. E. & Pitt, B. The renin-angiotensin-aldosterone system and its suppression. J. Vet. Intern. Med. 33(2), 363–382 (2019).

Article PubMed PubMed Central Google Scholar

Larouche-Lebel, É., Loughran, K. A., Huh, T. & Oyama, M. A. Effect of angiotensin receptor blockers and angiotensin converting enzyme 2 on plasma equilibrium angiotensin peptide concentrations in dogs with heart disease. J. Vet. Intern. Med. 35(1), 22–32 (2021).

Article PubMed Google Scholar

Patel, V. B., Zhong, J. C., Grant, M. B. & Oudit, G. Y. Role of the ACE2/angiotensin 1–7 axis of the renin-angiotensin system in heart failure. Circ. Res. 118(8), 1313–1326 (2016).

Article CAS PubMed PubMed Central Google Scholar

Eisenberg, M. J. & Gioia, L. C. Angiotensin II receptor blockers in congestive heart failure: Cardiology in review. Cardiol. Rev. 14(1), 26–34 (2006).

Article PubMed Google Scholar

Silveira, K. D. et al. Beneficial effects of the activation of the angiotensin-(1–7) Mas receptor in a murine model of adriamycin-induced nephropathy. PLoS ONE 8(6), e66082 (2013).

Article ADS CAS PubMed PubMed Central Google Scholar

King, J. N., Mauron, C. & Kaiser, G. Pharmacokinetics of the active metabolite of benazepril, benazeprilat, and inhibition of plasma angiotensin-converting enzyme activity after single and repeated administrations to dogs. Am. J. Vet. Res. 56(12), 1620–1628 (1995).

CAS PubMed Google Scholar

Authors/Task Force Members et al. ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure 2008: The Task Force for the Diagnosis and Treatment of Acute and Chronic Heart Failure 2008 of the European Society of Cardiology. Developed in collaboration with the Heart Failure Association of the ESC (HFA) and endorsed by the European Society of Intensive Care Medicine (ESICM). Eur. Heart J. 29(19), 2388–2442 (2008).

Article Google Scholar

Keene, B. W. et al. ACVIM consensus guidelines for the diagnosis and treatment of myxomatous mitral valve disease in dogs. J. Vet. Intern. Med. 33(3), 1127–1140 (2019).

Article PubMed PubMed Central Google Scholar

Toutain, P. L., Lefebvre, H. P. & King, J. N. Benazeprilat disposition and effect in dogs revisited with a pharmacokinetic/pharmacodynamic modeling approach. J. Pharmacol. Exp. Ther. 292(3), 1087–1093 (2000).

CAS PubMed Google Scholar

Hamlin, R. L. & Nakayama, T. Comparison of some pharmacokinetic parameters of 5 angiotensin-converting enzyme inhibitors in normal beagles. J. Vet. Intern. Med. 12(2), 93–95 (1998).

Article CAS PubMed Google Scholar

Ward, J. L., Chou, Y., Yuan, L., Dorman, K. S. & Mochel, J. P. Retrospective evaluation of a dose-dependent effect of angiotensin-converting enzyme inhibitors on long-term outcome in dogs with cardiac disease. J. Vet. Intern. Med. 35(5), 2102–2111 (2021).

Article PubMed PubMed Central Google Scholar

Mochel, J. P. et al. Chronobiology of the renin-angiotensin-aldosterone system in dogs: Relation to blood pressure and renal physiology. Chronobiol. Int. 30(9), 1144–1159 (2013).

Article CAS PubMed Google Scholar

Mochel, J. P. et al. Influence of feeding schedules on the chronobiology of renin activity, urinary electrolytes and blood pressure in dogs. Chronobiol. Int. 31(5), 715–730 (2014).

Article CAS PubMed Google Scholar

Defining a data set for Monolix. Monolix 2017. https://monolix.lixoft.com/data-and-models/creating-data-set/ (Accessed 15 Dec 2021).

Murray, J. D. 6.1 enzyme kinetics: basic enzyme reaction. In Mathematical biology I An introduction. 3rd ed. 175–8. (Springer, 17).

Felmlee, M. A., Morris, M. E. & Mager, D. E. Mechanism-based pharmacodynamic modeling. Methods Mol. Biol. 929, 583–600 (2012).

Article CAS PubMed PubMed Central Google Scholar

KEGG PATHWAY: Renin-angiotensin system—Homo sapiens (human). https://www.genome.jp/pathway/hsa04614 (Accessed 31 Jan 2022).

Bodineau, L., Frugière, A., Marc, Y., Claperon, C. & Llorens-Cortes, C. Aminopeptidase A inhibitors as centrally acting antihypertensive agents. Heart Fail. Rev. 13(3), 311–319 (2008).

Article CAS PubMed Google Scholar

Danziger, R. S. Aminopeptidase N in arterial hypertension. Heart Fail. Rev. 13(3), 293–298 (2008).

Article CAS PubMed Google Scholar

Schneider, B. et al. Model-based reverse translation between veterinary and human medicine: The one health initiative: Model-based reverse translational pharmacology. CPT Pharmacometrics Syst. Pharmacol. 7(2), 65–68 (2018).

Article CAS PubMed Google Scholar

Cediel, G. et al. Gender-related differences in heart failure biomarkers. Front. Cardiovasc. Med. 5(7), 617705 (2021).

Article Google Scholar

Falconnet, C., Bochud, M., Bovet, P., Maillard, M. & Burnier, M. Gender difference in the response to an angiotensin-converting enzyme inhibitor and a diuretic in hypertensive patients of African descent. J. Hypertens. 22(6), 1213–1220 (2004).

Article CAS PubMed Google Scholar

Li, J. S. et al. Racial differences in blood pressure response to angiotensin-converting enzyme inhibitors in children: A meta-analysis. Clin. Pharmacol. Ther. 84(3), 315–319 (2008).

Article CAS PubMed Google Scholar

Fagyas, M. et al. New perspectives in the renin-angiotensin-aldosterone system (RAAS) I: Endogenous angiotensin converting enzyme (ACE) inhibition. PLoS ONE 9(4), e87843 (2014).

Article ADS PubMed PubMed Central Google Scholar

Fagyas, M. et al. New perspectives in the renin-angiotensin-aldosterone system (RAAS) II: Albumin suppresses angiotensin converting enzyme (ACE) activity in human. PLoS ONE 9(4), e87844 (2014).

Article ADS PubMed PubMed Central Google Scholar

Ryan, J. W., Martin, L. C., Chung, A. & Pena, G. A. Mammalian inhibitors of angiotensin converting enzyme (kininase II). Adv. Exp. Med. Biol. 120B, 599–606 (1979).

CAS PubMed Google Scholar

Bieth, B., Bornkamp, B., Toutain, C., Garcia, R. & Mochel, J. P. Multiple comparison procedure and modeling: A versatile tool for evaluating dose-response relationships in veterinary pharmacology—A case study with furosemide. J. Vet. Pharmacol. Ther. 39(6), 539–546 (2016).

Article CAS PubMed Google Scholar

Pitt, B. et al. The effect of spironolactone on morbidity and mortality in patients with severe heart failure. Randomized Aldactone Evaluation Study Investigators. N. Engl. J. Med. 341(10), 709–717 (1999).

Article CAS PubMed Google Scholar

Adin, D. et al. Renin-angiotensin aldosterone profile before and after angiotensin-converting enzyme-inhibitor administration in dogs with angiotensin-converting enzyme gene polymorphism. J. Vet. Intern. Med. 34(2), 600–606 (2020).

Article PubMed PubMed Central Google Scholar

Guo, Z. et al. Effects of ramipril on the aldosterone/renin ratio and the aldosterone/angiotensin II ratio in patients with primary aldosteronism. Hypertension 76(2), 488–496 (2020).

Article CAS PubMed Google Scholar

Guo, Z. et al. Measurement of equilibrium angiotensin II in the diagnosis of primary aldosteronism. Clin. Chem. 66(3), 483–492 (2020).

Article PubMed Google Scholar

Index of /bin/windows/base/old/4.0.2. https://cran.r-project.org/bin/windows/base/old/4.0.2/ (Accessed 10 Dec 2021).

PubChem. Benazepril hydrochloride. https://pubchem.ncbi.nlm.nih.gov/compound/5362123 (Accessed 31 Jan 2022).

PubChem. Benazeprilat. https://pubchem.ncbi.nlm.nih.gov/compound/5463984 (Accessed 31 Jan 2022).

Monolix Methodology, a software for the analysis of nonlinear mixed effects models. (2014) http://lixoft.com/wp-content/uploads/2016/03/monolixMethodology.pdf (Accessed 10 May 2020).

Pelligand, L., Soubret, A., King, J. N., Elliott, J. & Mochel, J. P. Modeling of large pharmacokinetic data using nonlinear mixed-effects: A paradigm shift in veterinary pharmacology. A case study with robenacoxib in cats. CPT Pharmacometrics Syst. Pharmacol. 5(11), 625–635 (2016).

Article CAS PubMed PubMed Central Google Scholar

Sheiner, L. B. & Ludden, T. M. Population pharmacokinetics/dynamics*. Annu. Rev. Pharmacol. Toxicol. 32(1), 185–209 (1992).

Article CAS PubMed Google Scholar

Mentré, F. & Escolano, S. Prediction discrepancies for the evaluation of nonlinear mixed-effects models. J. Pharmacokinet. Pharmacodyn. 33(3), 345–367 (2006).

Article PubMed Google Scholar

FitzJohn, R., Jombart, T., Science IC of, Medicine T and. odin: ODE Generation and Integration. (2021) https://CRAN.R-project.org/package=odin (Accessed 10 Dec 2021).

Soetaert K, Petzoldt T, Setzer RW, ddaspk.f PNB (files, dvode.f, zvode.f), et al. deSolve: Solvers for Initial Value Problems of Differential Equations ("ODE", "DAE", ’DDE’) (2021) https://CRAN.R-project.org/package=deSolve (Accessed 10 Dec 2021).

FitzJohn, R., Hinsley, W., Science IC of, Medicine T and. dde: Solve Delay Differential Equations. (2020) https://CRAN.R-project.org/package=dde (Accessed 10 Dec 2021).

Download references

The authors would like to thank Dr. Bourgois-Mochel for her logistical support of the study and Dr. Oliver Domenig for the quantitative analysis of the pharmacodynamic data.

Funding for this study was provided by Ceva Sante Animale.

SMART Pharmacology, Iowa State University College of Vet. Medicine, 2448 Lloyd, 1809 S Riverside Dr., Ames, IA, 50011-1250, USA

Benjamin K. Schneider & Jonathan P. Mochel

Veterinary Clinical Sciences, Iowa State University, Ames, IA, 50011-1250, USA

Jessica Ward & Samantha Sotillo

Ceva Santé Animale, 33 500, Libourne, France

Catherine Garelli-Paar, Emilie Guillot & Marc Prikazsky

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

B.K.S.: methodology, investigation, formal analysis, software, validation, visualization, writing—original draft preparation; J.W.: conceptualization, data curation, investigation, writing—original draft preparation; S.S.: data curation, investigation, writing—review and editing; C.G.P., M.P.: conceptualization, investigation, writing—review and editing; E.G.: conceptualization, investigation, writing—original draft preparation; J.P.M.: supervision, project administration, conceptualization, methodology, investigation, funding acquisition, resources, formal analysis, validation, writing—original draft preparation.

Correspondence to Jonathan P. Mochel.

Authors JPM and JW have served as consultants for Ceva Sante Animale and have received reimbursement and honoraria for consulting, expert testimony, travel, and service as key opinion leaders (KOLs). Authors Garelli-Paar, Guillot and Prikazsky are employees of Ceva Sante Animale. Authors SS and BKS do not have any conflict of interest.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Schneider, B.K., Ward, J., Sotillo, S. et al. Breakthrough: a first-in-class virtual simulator for dose optimization of ACE inhibitors in translational cardiovascular medicine. Sci Rep 13, 3300 (2023). https://doi.org/10.1038/s41598-023-30453-x

Download citation

Received: 21 November 2022

Accepted: 23 February 2023

Published: 26 February 2023

DOI: https://doi.org/10.1038/s41598-023-30453-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.