Next Article in Journal
Neonatal Pharmacokinetics and Biodistribution of Polymeric Nanoparticles and Effect of Surfactant
Next Article in Special Issue
Donepezil Brain and Blood Pharmacokinetic Modeling after Nasal Film and Oral Solution Administration in Mice
Previous Article in Journal
Combined Nanofibrous Face Mask: Co-Formulation of Lipases and Antibiotic Agent by Electrospinning Technique
Previous Article in Special Issue
Understanding the Mechanisms and Treatment of Heart Failure: Quantitative Systems Pharmacology Models with a Focus on SGLT2 Inhibitors and Sex-Specific Differences
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of a Cardiovascular Systems Model for Design and Analysis of Hemodynamic Safety Studies

1
Leiden Academic Centre for Drug Research, Leiden University, Einsteinweg 55, 2333 CC Leiden, The Netherlands
2
LAP&P Consultants BV, Archimedesweg 31, 2333 CM Leiden, The Netherlands
3
Certara QSP, Canterbury CT2 7FG, UK
*
Author to whom correspondence should be addressed.
Pharmaceutics 2023, 15(4), 1175; https://doi.org/10.3390/pharmaceutics15041175
Submission received: 11 January 2023 / Revised: 20 March 2023 / Accepted: 5 April 2023 / Published: 7 April 2023

Abstract

:
Early prediction, quantification and translation of cardiovascular hemodynamic drug effects is essential in pre-clinical drug development. In this study, a novel hemodynamic cardiovascular systems (CVS) model was developed to support these goals. The model consisted of distinct system- and drug-specific parameter, and uses data for heart rate (HR), cardiac output (CO), and mean atrial pressure (MAP) to infer drug mode-of-action (MoA). To support further application of this model in drug development, we conducted a systematic analysis of the estimation performance of the CVS model to infer drug- and system-specific parameters. Specifically, we focused on the impact on model estimation performance when considering differences in available readouts and the impact of study design choices. To this end, a practical identifiability analysis was performed, evaluating model estimation performance for different combinations of hemodynamic endpoints, drug effect sizes, and study design characteristics. The practical identifiability analysis showed that MoA of drug effect could be identified for different drug effect magnitudes and both system- and drug-specific parameters can be estimated precisely with minimal bias. Study designs which exclude measurement of CO or use a reduced measurement duration still allow the identification and quantification of MoA with acceptable performance. In conclusion, the CVS model can be used to support the design and inference of MoA in pre-clinical CVS experiments, with a future potential for applying the uniquely identifiable systems parameters to support inter-species scaling.

1. Introduction

Cardiovascular toxicity remains one of the most significant causes of market withdrawal [1,2,3]. In this context, pre-clinical and clinical cardiovascular safety studies and associated quantitative analyses have primarily focused on early identification of QT prolongation effects, associated with an increased risk of torsade de pointes. However, drugs may also be associated with other non-QT cardiovascular hemodynamic drug effects; such effects are now increasingly being considered in cardiovascular safety studies [4]. To assess these hemodynamic drug effects, studies in animal models, such as dog and rat, are commonly conducted while measuring readouts, such as heart rate (HR), mean arterial pressure (MAP), and cardiac output (CO) or cardiac contractility. One key purpose of such pre-clinical studies is identifying the specifical hemodynamic endpoint on which the drug is acting, i.e., the mode-of-action (MoA). Due to the complex regulatory system which drives the cardiovascular hemodynamic response and the subsequent relationship between observed changes in hemodynamic readouts, identifying the MoA in pre-clinical cardiovascular studies is not simple [5]. In addition, data from such pre-clinical studies should ideally be used to predict hemodynamic drug effects in other species, including humans; however, it is hampered by potential inter-species differences.
Quantitative mathematical models may have potential to support the optimization of experimental designs of pre-clinical and clinical hemodynamic studies. In past studies, a quantitative hemodynamic cardiovascular systems model (CVS) was developed using experimental hemodynamic study data for multiple drugs in rats, including measurements of HR, CO, and MAP [6,7]. This CVS model characterizes short-term drug effects on the inter-relationship between five hemodynamic variables, i.e., HR, stroke volume (SV), total peripheral resistance (TPR), CO, and MAP. The model allows incorporation of drug effects associated with these variables to evaluate and quantify the likely drug mode-of-action (MoA). The use of distinct drug- and system-specific parameters in the model may also enable inter-species scaling of drug effects towards other species, such as dogs [8] and humans [9]. As such, the CVS model has potential value to support classification of the drug MoA and the quantification of its effect size, supporting pre-clinical study design optimization; it could also produce translational cross-species predictions, such as supporting dose selection in first-time-in-human studies. The CVS model was developed based on experimental data for a selected set of compounds with specific characteristics. Broader application of the CVS model is currently limited by a lack of general quantitative characterization of the model’s performance in terms of correct identification and unbiased estimation of hemodynamic drug effects; any broad application must take into account study design considerations, such as the effect of differences in the type of available hemodynamic readouts, sample size, drug doses, and the magnitude of the effect size.
In this study, we describe a systematic analysis of the estimation performance of the CVS model to infer drug- and system-specific parameters, evaluating its utility to support the design and analysis of pre-clinical cardiovascular studies. Specifically, we aimed to evaluate the following elements for the CVS model: (1) the practical identifiability of drug- and system-specific parameters; and (2) the estimation performance in terms of correct identification and parameter estimation precision in relation to drug effects, available hemodynamic readouts, and experimental study design characteristics.

2. Methods

2.1. Model Structure

The hemodynamic systems model used in this analysis was developed by Snelder et al. [7]; it will hereafter be referred to as the CVS model. In summary, the CVS model includes the inter-relationship of five hemodynamic variables, including HR, SV, CO, TPR, and MAP (Figure 1). CO is defined as the product of HR and SV, while MAP is the product of CO and TPR. The drug effect is described with an indirect response model reflecting the impact of drug on the turnover of primary hemodynamic biomarkers HR, SV, or TPR, while CO, MAP, and HR data were used as readouts for model estimation. The previously described handling effect was excluded from this analysis since its magnitude is much smaller than the magnitude of circadian rhythm or drug effect. The differential equations of the model are as follows:
d H R d t = K i n _ H R · ( 1 + C R H R ) · 1 F B M A P · ( 1 + E F F ) K o u t _ H R · H R
d S V * d t = K i n _ S V · 1 F B M A P · ( 1 + E F F ) K o u t _ S V · S V *
d T P R d t = K i n _ T P R · ( 1 + C R T P R ) · 1 F B M A P · ( 1 + E F F ) K o u t _ T P R · T P R
S V = S V * × [ 1 H R S V × L N H R B S L H R ]
C O = H R · S V
M A P = C O · T P R
C R H R = a m p H R · c o s ( 2 π · t + h o r H R 24 )
C R T P R = a m p T P R · c o s ( 2 π · t + h o r T P R 24 )
E F F = E m a x · C E C 50 + C
The model equations consist of system-specific parameters (i.e., feedback-, turnover-, and circadian rhythm-related parameters) and drug-specific parameters (i.e., EC50 and Emax). The FB constant represents the negative feedback of MAP on the turnover of HR, SV, and TPR. The parameters Kin_HR, Kin_SV, and Kin_TPR are the zero-order production constants for HR, SV and TPR, while Kout_HR, Kout_SV, and Kout_TPR are the first-order dissipation constant for HR, SV and TPR, respectively. BSL is the baseline value of each hemodynamic variable. Circadian rhythms (CR) for HR and TPR were modelled with a cosine function, as previously described. For the CR equations, amp and hor represent amplitude and horizontal displacement, respectively. HRSV is a constant to describe the magnitude of direct effect of HR and SV. Regarding pharmacokinetics (PK), we assumed a hypothetical drug with the PK described according to a one-compartment model. The elimination rate constant (K) of the hypothetical drug was set as 0.173 h−1 corresponding to a realistic half-life of 4 h. A typical Emax model was assumed to link the drug exposure to the drug effect (EFF) on either HR, SV, or TPR. C represents the drug concentration at the central compartment. Drug-specific parameters EC50 and Emax represent the concentration at which half of the maximum and maximum effects are achieved, respectively. Inter-individual variability was included on baseline parameters (BSL) for HR, CO and MAP. The full model code of the model is included in Supplemental Materials S1.

2.2. Simulation Scenarios

All simulation scenarios, including study design characteristics and model parameters used, are summarized in detail in Table S1; the scenarios included one scenario for the quantification of system- and drug-specific parameters and three scenarios for identification of the MoA with different drug effects, different observation duration, or different numbers of animals. We briefly simulated dosing regimens for animals receiving either (i) three ascending doses: 0.1 mg/kg, 1 mg/kg, and 10 mg/kg i.v. bolus on day 1, 2, and 3; or (ii) five ascending doses: 0.1 mg/kg, 0.3 mg/kg, 1 mg/kg, 3 mg/kg, and 10 mg/kg i.v. bolus on days 1 to 5. These dosing regimens follow commonly applied study designs. We adopted a dense sampling design for the simulations, taking hourly measurements after drug administration of HR, CO, and MAP as, in these studies, rich sampling using telemetry is common practice. For all the scenarios focusing on identification of the MoA, drug-specific parameters (Emax and EC50), interindividual variability (IIV), and residual errors of each type of measurement were estimated; the system-specific parameters were fixed to the parameter values from the published CVS model [7].

2.3. Practical Identifiability Analysis

2.3.1. Stochastic Stimulation and Re-Estimation

The practical identifiability analyses of the CVS model for identification of drug- and systems parameters, and the evaluation of drug effect, available readouts, and study designs, were all performed using stochastic simulations and re-estimation (SSE) analyses. SSEs were performed using Perl-speaks-NONMEM toolkit (PsN version 4.8.1, Uppsala University, Uppsala, Sweden) in conjunction with NONMEM (version 7.4.3, Icon Development Solutions, Ellicott, MD, USA). The SSE is a commonly used tool for model comparison and hypothesis testing. In SSE analyses, a number of simulated datasets were generated using the input model (original model). The simulated datasets were then used to fit the original model and a set of alternative models. Ultimately, bias and precision of parameter estimates, as well as statistical significance of tested models, can be computed. A classical first-order conditional estimation with interaction (FOCE-I) method was used for model estimation.

2.3.2. Identifiability of Drug- and System-Specific Parameters

We evaluated whether the CVS model could be used to identify drug- and system-specific parameters based on common hemodynamic study designs. For this assessment, we implemented CVS model structures that included drug effect parameters on one of the three primary hemodynamic variables (HR, TPR, and SV) and a model which did not include a drug effect on hemodynamics as the alternative model. We first simulated a number of datasets using the models with drug effects on either HR, SV, or TPR, respectively (original models). We then re-fitted the simulated datasets using the original (i.e., true) models, as well as alternative models with either drug effect on the other hemodynamic variables (i.e., incorrect models) or no drug effect included at all. For each scenario, we performed 100 simulations and used them to fit both original models and alternative models (re-estimations). For each re-estimated model, we evaluated the precision of parameter estimates using the relative prediction error (RPE), calculated as follows:
R P E % = P e s t P s i m P s i m × 100 %
with Pest being the re-estimated parameter value and Psim the original, simulated value.
To assess the ability of each model to identify the correct drug effect, we calculated the decrease in objective function value (ΔOFV) between the base and the alternative models for each simulation scenario (Figure 2). We classified a significant drug effect as a decrease in OFV > 5.99 (p < 0.05, degree of freedom = 2 in chi square test), compared to the model with no drug effect.

2.3.3. Drug Effect Magnitudes

We performed SSEs with either an inhibitory or stimulatory drug effect on HR, SV, or TPR to confirm if the MoA could be identified correctly for different drug effect magnitudes and directionalities. The Emax in the original models was set to 1 for inhibitory effects so that the production rate term of the primary biomarkers was positive. The value of Emax was set as 1 or 10 for stimulatory effects. EC50 was tested for 100 ng/mL or 1000 ng/mL to simulate different magnitudes of drug effects.

2.3.4. Observation Duration

We performed SSEs for study durations of 3, 6, 12, and 24 h after administration of three ascending doses to compare the identification performance of MoA and the precision of drug-specific parameters with different observation durations. Both a large drug effect (EC50 = 100 ng/mL) and a small drug effect (EC50 = 100 µg/mL) were tested.

2.3.5. Number of Animals

To determine how the number of animals used in the hemodynamic studies influences the identification of MoA and precision of parameter estimation, SSEs using three, four, or five animals were conducted. Both a large drug effect (EC50 = 100 ng/mL) and a small drug effect (EC50 = 100 µg/mL) were tested.

3. Results

3.1. Identifiability of Drug- and System-Specific Parameters

The mean relative prediction errors of drug- and system-specific parameters for both studied dosing regimens showed minimal bias (relative prediction error < 20%), suggesting that all the system- and drug-specific parameters can be adequately estimated (Figure 3). Depending on the action of the drug effect, Kout parameters, CO0, and HRSV frequently showed an increased prediction error of up to 20%, in comparison to other parameters’ prediction errors of <5%. Inclusion of CO data, besides HR and MAP data, decreases the prediction errors, albeit not substantially. Not all test models could minimize successfully. The model refitted to the data including CO observations showed superior overall stability of estimations with a significantly (p < 0.05) higher percentage of successful minimizations compared to simulations that excluded CO observations. It was found that the median of the relative prediction error of parameters in SSE with five ascending doses was closer to 0 compared to the ones in SSE with three ascending doses (results on file); this indicates that the study design of five ascending doses could help to improve the accuracy of parameter estimation. Finally, we compared the final estimates of EC50 of models with unfixed and fixed system-specific parameters. The relative prediction errors of EC50 were found within ±10%, which indicates that fixing system-specific parameters has no significant influence on the estimation of EC50. (Figure S1); this finding may be considered when applying the model for analysis of hemodynamic studies.

3.2. Drug Effect Magnitudes

We found that for all simulated drug effect magnitudes and effect sites (i.e., HR, SV, TPR), the correct MoA could always be identified using a standard study design, i.e., for all re-estimated runs, substantially lower ΔOFV values were consistently observed for the correct model in comparison to alternative models with an incorrect MoA (results on file). These results confirm the utility of the CVS model to correctly identify drug MoA using standard study designs.

3.3. Observation Duration

A variable impact was observed of the duration of experimental observations on the correct identification of drug effect sites for high EC50 values (100 ug/mL), i.e., when the drug effect is close to the circadian variation (Figure 4). In contrast, for effects on HR an observation duration of 6 h is sufficient to identify the correct MoA; identification of the correct MoA for TPR and SV is only possible with an intermediate or low percentage of runs. Here, addition of CO data substantially increases the models’ ability to identify the correct MoA. When considering the relative prediction error for EC50 for these scenarios, we found that the extent of bias was correlated with the ability to identify the correct drug effect, with substantial bias in the EC50 for TPR and SV (50–100%) for shorter study durations; there was also a positive effect in reducing bias (Figure 5). The added value of including CO in terms of bias was less distinct. In contrast, when a lower EC50 value is used (Figure 6), the model can identify the correct mode-of-action for all drug MoAs; this is true even for a minimal study duration of 3 h, with only minimal bias in the EC50 of <5% (Figure S2). These findings highlight the ability of this model to infer hemodynamic MoA for short observation times and even to detect modest drug effects.

3.4. Number of Animals

A distinct effect of sample size (ranging between 3 to 5 animals) was found for different MoA for high EC50 values (100 ug/mL). For drug effect on HR, for all sample sizes the correct MoA could always be identified, whereas for SV and TPR the ability to identify a correct drug effect was much lower (Figure 7). Increasing both the number of animals and the additional measurement of CO substantially increased the model’s ability to identify the MoA for SV and TPR. Similar trends were observed for the relative prediction error distributions (Figure 8). When a lower EC50 value is used (Figure 6), the model can identify the correct mode-of-action for all drug MoAs for all sample sizes studied, with only minimal bias in the EC50 (Figure S3). These findings highlight the ability of this model to infer hemodynamic MoA for short observation times and even to detect modest drug effects.

4. Discussion

In this analysis, we demonstrated the utility of the CVS model framework for inference of drug- and system-specific parameters and identifying the correct MoA for different drug effect magnitudes and study design choices through systematic simulation and re-estimation-based practical identifiability analyses.
In our simulation analyses, we found that for commonly used experimental study designs, system- and drug-specific parameters could be identified with acceptable parameter estimation precision and bias, based on single-drug studies. We also found that the availability of additional dose levels (i.e., 3 versus 5 dose levels as tested) improved the stability and precision of the estimation procedure. Alternatively, as described previously by Snelder et al. [6,7], fitting datasets across multiple drugs may aid in more reliable estimation of system-specific parameters. Once such systems parameters are established, e.g., for a specific experimental setting, the systems parameters may be fixed to make the model more robust in identifying MoA and drug effect parameters, when applying the model to other drugs or drug candidates
We demonstrated how experimental study designs, such as duration of experimental observations and the sample size, can substantially influence the model’s ability to correctly identify drug MoA and estimate the associated drug effect parameters. In particular, when drug effects have a more modest magnitude, there is a substantial impact from the duration of observations and sample size. Importantly, for more modest drug effects, we showed that a model-based approach can still characterize the hemodynamic drug effects effectively. From our analyses, it is also clear that the power to identify drug effects at specific effect sites, i.e., HR, SV, and TPR, also varies substantially. Because the circadian rhythm magnitudes and residual errors of these variables are similar, it is likely this difference from the substantially different kout values for these parameters is related to the magnitude of the circadian rhythm effects for these parameters relative to the drug effects (Figure 6). In this context, the identification of drug effects on SV was found to be most challenging, while inclusion of CO measurements dramatically improves identification of SV-associated drug effect. In addition, extension of the observation duration (Figure 5) or increasing the amount of animals was shown to enhance ability to identify SV-associated drug effects.
Our results show how design considerations, such as dose, sample size, and observation duration, all impact the estimation performance of hemodynamic drug effects and other model parameters. To this end, model-informed optimization of experimental designs is a relevant approach to ensure that both optimally informative experiments are performed and the amount of animals required can be minimized for ethical considerations [10,11,12]. Since intensive sampling of hemodynamic variables with telemetry is common in pre-clinical experiments, especially for measurements of HR and MAP, comparison of model performance with intensive and sparse sampling designs was not considered in our analysis. However, achieving a dense sampling frequency of hemodynamic variables can be more challenging to achieve, in particular in the clinic. To further explore the impact of sampling frequency in clinical contexts on parameter estimation and identification of mode-of-action, potential study designs can be evaluated using stochastic simulation–estimation strategies using the CVS model in the future, similar to the workflow described in the current report. Alternatively, the use of optimal design methodologies, such as D-optimal design, could be a relevant and more efficient strategy for evaluation and optimization of pre-clinical and clinical hemodynamic studies.
We found that CO data may not be required for the identification of MoA. This is important because CO is often measured invasively using arterial cannulation and a pulmonary artery catheter placement, which makes CO more complex to measure in comparison to HR and MAP data. The CVS model can, thus, be practically convenient for studying MoA in drug development. However, if available, CO data may be of value since we also found that the inclusion of CO data improved the precision of the parameter estimates and increased the rate of successful minimization in the second SSE analyses. In recent years, impedance cardiography has been used as an accurate non-invasive and lower-cost approach for obtaining frequent serial measurements of CO in clinical contexts. This may potentially benefit the use of the models in hemodynamic studies when the new technology becomes more routinely available in the future [13,14].
Although our analysis is fully based on simulation studies, we expect that our conclusions will also be applicable to experimental hemodynamic studies. Importantly, the CVS model used in our analyses has been developed based on extensive and densely sampled hemodynamic experiments for multiple drugs as previously described [7]; these experiments were able to accurately describe the observed hemodynamic responses as well as the circadian rhythms of cardiovascular variables, thereby forming the basis for the simulations in the current analysis. One factor not included in the current analysis is the potential situation when compounds have a dual mode-of-action, i.e., which affect multiple hemodynamic endpoints at the same time. In addition, we focused on drugs with a typical half-life of 4 h. For drugs with for example much longer half-lives, we expect that the estimation of kout might be affected and could require extended study durations to achieve the same power to estimate the kout and/or drug effect parameters. The parameter values used in this analysis were based on rat hemodynamic studies. However, we expect that our main conclusions can also be applied to other species if a species-specific model has previously been derived. Indeed, the hemodynamic model has recently been applied to dog hemodynamic studies [8].
In conclusion, we demonstrated the utility of the CVS model for identifying and quantifying drug MoA in pre-clinical cardiovascular studies. The MoA can be identified using the CVS model based on MAP and HR data, while the inclusion of CO data could further improve the model estimation. The findings in this study can also contribute to optimization of experimental study design and reduce the number of animals used in pre-clinical studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pharmaceutics15041175/s1, Supplemental Material S1: Model code (NONMEM); Table S1: Simulation scenarios of the practical identifiability analysis a; Figure S1: Comparison of final estimates of EC50 in models with unfixed and fixed system-specific parameters; Figure S2: Relative prediction error of EC50 in SSE analyses with observations of HR and MAP or observations of HR, CO and MAP within 3 h, 6 h, 12 h and 24 h, while Emax was fixed as −1 and EC50 was fixed as 100 ng/mL in the original model. The relative prediction error is defined as the ratio of the difference of the parameters’ estimates and the parameters’ true values over the parameters’ true values; Figure S3: Relative prediction error of EC50 in SSE analyses using data for 3, 4, or 5 animals with observations of HR and MAP or observations of HR, CO and MAP, while Emax was fixed as −1 and EC50 was fixed as 100 ng/mL in the original model. The relative prediction error is defined as the ratio of the difference of the parameters’ estimates and the parameters’ true values over the parameters’ true values.

Author Contributions

Conceptualization, Y.F., N.S. and J.G.C.v.H.; methodology, Y.F. and N.S.; formal analysis, Y.F.; writing—original draft preparation, Y.F. and J.G.C.v.H.; writing—review and editing, Y.F., N.S., T.G., P.H.v.d.G. and J.G.C.v.H.; supervision, J.G.C.v.H., N.S. and P.H.v.d.G. All authors have read and agreed to the published version of the manuscript.

Funding

This project has received funding from the Innovative Medicines Initiative 2 Joint Undertaking under grant agreement No. 116030. This Joint Undertaking receives support from the European Union’s Horizon 2020 research and innovation program and EFPIA. Y.F. received funding from the China Scholarship Council.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw data is available upon request.

Conflicts of Interest

The authors declared no competing interests for this work.

References

  1. Laverty, H.G.; Benson, C.; Cartwright, E.; Cross, M.; Garland, C.; Hammond, T.; Holloway, C.; McMahon, N.; Milligan, J.; Park, B.; et al. How can we improve our understanding of cardiovascular safety liabilities to develop safer medicines? Br. J. Pharmacol. 2011, 163, 675–693. [Google Scholar] [CrossRef] [PubMed]
  2. Munos, B. Lessons from 60 years of pharmaceutical innovation. Nat. Rev. Drug Discov. 2009, 8, 959–968. [Google Scholar] [CrossRef] [PubMed]
  3. Cook, D.; Brown, D.; Alexander, R.; March, R.; Morgan, P.; Satterthwaite, G.; Pangalos, M.N. Lessons learned from the fate of AstraZeneca’s drug pipeline: A five-dimensional framework. Nat. Rev. Drug Discov. 2014, 13, 419–431. [Google Scholar] [CrossRef] [PubMed]
  4. FDA Assessment of Pressor Effects of Drugs Guidance for Industry DRAFT GUIDANCE. 2022. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/assessment-pressor-effects-drugs-guidance-industry (accessed on 1 January 2023).
  5. Oparil, S.; Zaman, M.A.; Calhoun, D.A. Pathogenesis of Hypertension. Ann. Intern. Med. 2003, 139, 761. [Google Scholar] [CrossRef] [PubMed]
  6. Snelder, N.; Ploeger, A.B.; Luttringer, O.; Rigel, D.F.; Webb, R.L.; Feldman, D.; Fu, F.; Beil, M.; Jin, L.; Stanski, D.R.; et al. PKPD modelling of the interrelationship between mean arterial BP, cardiac output and total peripheral resistance in conscious rats. Br. J. Pharmacol. 2013, 169, 1510–1524. [Google Scholar] [CrossRef] [PubMed]
  7. Snelder, N.; Ploeger, A.B.; Luttringer, O.; Rigel, D.F.; Fu, F.; Beil, M.; Stanski, D.R.; Danhof, M. Drug effects on the CVS in conscious rats: Separating cardiac output into heart rate and stroke volume using PKPD modelling. Br. J. Pharmacol. 2014, 171, 5076–5092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Venkatasubramanian, R.; Collins, T.A.; Lesko, L.J.; Mettetal, J.T.; Trame, M.N. Semi-mechanistic modelling platform to assess cardiac contractility and haemodynamics in preclinical cardiovascular safety profiling of new molecular entities. Br. J. Pharmacol. 2020, 177, 3568–3590. [Google Scholar] [CrossRef] [PubMed]
  9. Bahnasawy, S.; Al-Sallami, H.; Duffull, S. A minimal model to describe short-term haemodynamic changes of the cardiovascular system. Br. J. Clin. Pharmacol. 2020, 87, 1411–1421. [Google Scholar] [CrossRef] [PubMed]
  10. Huang, W.; du Sert, N.P.; Vollert, J.; Rice, A.S.C. General Principles of Preclinical Study Design. Good Res. Pract. Non Clin. Pharmacol. Biomed. 2019, 257, 55–69. [Google Scholar]
  11. Landis, S.C.; Amara, S.G.; Asadullah, K.; Austin, C.P.; Blumenstein, R.; Bradley, E.W.; Crystal, R.G.; Darnell, R.B.; Ferrante, R.J.; Fillit, H.; et al. A call for transparent reporting to optimize the predictive value of preclinical research. Nature 2012, 490, 187–191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Curtis, M.J.; Hancox, J.C.; Farkas, A.; Wainwright, C.L.; Stables, C.L.; Saint, D.A.; Clements-Jewery, H.; Lambiase, P.D.; Billman, G.E.; Janse, M.J.; et al. The Lambeth Conventions (II): Guidelines for the study of animal and human ventricular and supraventricular arrhythmias. Pharmacol. Ther. 2013, 139, 213–248. [Google Scholar] [CrossRef] [PubMed]
  13. Ventura, H.O.; Taler, S.J.; Strobeck, J.E. Hypertension as a hemodynamic disease: The role of impedance cardiography in diagnostic, prognostic, and therapeutic decision making. Am. J. Hypertens. 2005, 18, 26–43. [Google Scholar] [CrossRef] [PubMed]
  14. Koobi, T.; Kaukinen, N.A.; Turjanmaa, V.M.H.; Uusitalo, A.J. Whole-body impedance cardiography in the measurement of cardiac output. Crit. Care Med. 1997, 25, 779–785. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Model structure of the CVS model consisting of five hemodynamic variables of HR, SV, TPR, CO, and MAP and a one-compartment pharmacokinetic model of the hypothetical drug [7]. CO is defined as the product of HR and SV, while MAP is the product of CO and TPR. The drug affects the production rate of HR, SV, or TPR. C: drug concentration (ng/mL); HR: heart rate (beats/min); SV: stroke volume (mL/beat); TPR: total peripheral resistance (mmHg∙min/mL); CO: cardiac output (mL/min); MAP: mean arterial pressure (mmHg); FB: feedback effect of MAP (mmHg−1); Edrug: drug effect on three primary effect sites. Kin: production rate constant (Kin_HR: beats*h−1/min, Kin_SV: mL*h−1/beat, Kin_TPR: mmHg*min*h−1/mL); Kout: dissipation rate constant (h−1); k: elimination rate constant of the hypothetical drug (h−1).
Figure 1. Model structure of the CVS model consisting of five hemodynamic variables of HR, SV, TPR, CO, and MAP and a one-compartment pharmacokinetic model of the hypothetical drug [7]. CO is defined as the product of HR and SV, while MAP is the product of CO and TPR. The drug affects the production rate of HR, SV, or TPR. C: drug concentration (ng/mL); HR: heart rate (beats/min); SV: stroke volume (mL/beat); TPR: total peripheral resistance (mmHg∙min/mL); CO: cardiac output (mL/min); MAP: mean arterial pressure (mmHg); FB: feedback effect of MAP (mmHg−1); Edrug: drug effect on three primary effect sites. Kin: production rate constant (Kin_HR: beats*h−1/min, Kin_SV: mL*h−1/beat, Kin_TPR: mmHg*min*h−1/mL); Kout: dissipation rate constant (h−1); k: elimination rate constant of the hypothetical drug (h−1).
Pharmaceutics 15 01175 g001
Figure 2. Overview identifiability analyses. SSE: stochastic simulation and re-estimation; PsN: Perl-speaks-NONMEM toolkit; HR: heart rate; SV: stroke volume; TPR: total peripheral resistance; MoA: mode-of-action; OFV: objective function value.
Figure 2. Overview identifiability analyses. SSE: stochastic simulation and re-estimation; PsN: Perl-speaks-NONMEM toolkit; HR: heart rate; SV: stroke volume; TPR: total peripheral resistance; MoA: mode-of-action; OFV: objective function value.
Pharmaceutics 15 01175 g002
Figure 3. Identifiability of drug- and system-specific parameters. HR: heart rate; CO: cardiac output; SV: stroke volume; TPR: total peripheral resistance; MAP: mean arterial pressure. Relative prediction errors after simulation and re-estimation (SSE) analysis using observations for HR and MAP (top) or observations of HR, CO, and MAP (bottom) following three ascending doses. EC50: the concentration at which half of maximum effect is achieved; FB: regulatory feedback effect from MAP to HR, SV, and TPR; HR0: baseline value of HR; MAP0: baseline value of MAP; CO0: baseline value of CO; KoutHR: dissipation rate of HR; KoutSV: dissipation rate of SV; KoutTPR: dissipation rate of TPR; HRSV: direct effect constant between HR and SV; HorHR: horizontal displacement of circadian rhythm of HR; AmpHR: amplitude of circadian rhythm of HR; HorTPR: horizontal displacement of circadian rhythm of TPR.
Figure 3. Identifiability of drug- and system-specific parameters. HR: heart rate; CO: cardiac output; SV: stroke volume; TPR: total peripheral resistance; MAP: mean arterial pressure. Relative prediction errors after simulation and re-estimation (SSE) analysis using observations for HR and MAP (top) or observations of HR, CO, and MAP (bottom) following three ascending doses. EC50: the concentration at which half of maximum effect is achieved; FB: regulatory feedback effect from MAP to HR, SV, and TPR; HR0: baseline value of HR; MAP0: baseline value of MAP; CO0: baseline value of CO; KoutHR: dissipation rate of HR; KoutSV: dissipation rate of SV; KoutTPR: dissipation rate of TPR; HRSV: direct effect constant between HR and SV; HorHR: horizontal displacement of circadian rhythm of HR; AmpHR: amplitude of circadian rhythm of HR; HorTPR: horizontal displacement of circadian rhythm of TPR.
Pharmaceutics 15 01175 g003
Figure 4. Impact of duration of experimental observations on the percentage of runs with correct identification of mechanisms of action. Results from simulation and re-estimation analyses using either observations of HR and MAP or observations of HR, CO, and MAP within 3 h, 6 h, 12 h, and 24 h. Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in the original model. Emax was fixed to −1 and EC50 was fixed to 100 ug/mL in the original model. Correct identification is based on applying a likelihood ratio test compared to the base (no drug) model.
Figure 4. Impact of duration of experimental observations on the percentage of runs with correct identification of mechanisms of action. Results from simulation and re-estimation analyses using either observations of HR and MAP or observations of HR, CO, and MAP within 3 h, 6 h, 12 h, and 24 h. Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in the original model. Emax was fixed to −1 and EC50 was fixed to 100 ug/mL in the original model. Correct identification is based on applying a likelihood ratio test compared to the base (no drug) model.
Pharmaceutics 15 01175 g004
Figure 5. Impact of duration of experimental observations on relative prediction error of drug effect estimate EC50. Results from simulation and re-estimation analyses using either observations of HR and MAP or observations of HR, CO, and MAP within 3 h, 6 h, 12 h, and 24 h, while Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in the original model.
Figure 5. Impact of duration of experimental observations on relative prediction error of drug effect estimate EC50. Results from simulation and re-estimation analyses using either observations of HR and MAP or observations of HR, CO, and MAP within 3 h, 6 h, 12 h, and 24 h, while Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in the original model.
Pharmaceutics 15 01175 g005
Figure 6. Comparison of magnitudes of circadian rhythm on HR or TPR and magnitudes of drug effect. Comparison of EC50 as 100 ng/mL (left), 1000 ng/mL (middle), and 100 ug/mL (right). Dark blue solid line: drug effect following a dose of 0.1 mg/kg; orange solid line: drug effect following a dose of 1 mg/kg; green solid line: drug effect following a dose of 10 mg/kg; black dashed line: magnitude of circadian rhythm effect of HR; gray dashed line: circadian rhythm effect of TPR.
Figure 6. Comparison of magnitudes of circadian rhythm on HR or TPR and magnitudes of drug effect. Comparison of EC50 as 100 ng/mL (left), 1000 ng/mL (middle), and 100 ug/mL (right). Dark blue solid line: drug effect following a dose of 0.1 mg/kg; orange solid line: drug effect following a dose of 1 mg/kg; green solid line: drug effect following a dose of 10 mg/kg; black dashed line: magnitude of circadian rhythm effect of HR; gray dashed line: circadian rhythm effect of TPR.
Pharmaceutics 15 01175 g006
Figure 7. Impact of number of animals used in experiments on percentage of runs with correct identification of mechanism of action. Derived using SSE analyses using data for 3, 4, or 5 animals with observations of HR and MAP or observations of HR, CO, and MAP, while Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in original model. Correct identification is based on applying a likelihood ratio test compared to the base (no drug) model.
Figure 7. Impact of number of animals used in experiments on percentage of runs with correct identification of mechanism of action. Derived using SSE analyses using data for 3, 4, or 5 animals with observations of HR and MAP or observations of HR, CO, and MAP, while Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in original model. Correct identification is based on applying a likelihood ratio test compared to the base (no drug) model.
Pharmaceutics 15 01175 g007
Figure 8. Impact of number of animals used in experiments on relative prediction error of drug effect estimate EC50. Derived using SSE analyses using data for 3, 4, or 5 animals with observations of HR and MAP or observations of HR, CO, and MAP, while Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in the original model.
Figure 8. Impact of number of animals used in experiments on relative prediction error of drug effect estimate EC50. Derived using SSE analyses using data for 3, 4, or 5 animals with observations of HR and MAP or observations of HR, CO, and MAP, while Emax was fixed as −1 and EC50 was fixed as 100 ug/mL in the original model.
Pharmaceutics 15 01175 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fu, Y.; Snelder, N.; Guo, T.; van der Graaf, P.H.; van Hasselt, J.G.C. Evaluation of a Cardiovascular Systems Model for Design and Analysis of Hemodynamic Safety Studies. Pharmaceutics 2023, 15, 1175. https://doi.org/10.3390/pharmaceutics15041175

AMA Style

Fu Y, Snelder N, Guo T, van der Graaf PH, van Hasselt JGC. Evaluation of a Cardiovascular Systems Model for Design and Analysis of Hemodynamic Safety Studies. Pharmaceutics. 2023; 15(4):1175. https://doi.org/10.3390/pharmaceutics15041175

Chicago/Turabian Style

Fu, Yu, Nelleke Snelder, Tingjie Guo, Piet H. van der Graaf, and Johan. G. C. van Hasselt. 2023. "Evaluation of a Cardiovascular Systems Model for Design and Analysis of Hemodynamic Safety Studies" Pharmaceutics 15, no. 4: 1175. https://doi.org/10.3390/pharmaceutics15041175

APA Style

Fu, Y., Snelder, N., Guo, T., van der Graaf, P. H., & van Hasselt, J. G. C. (2023). Evaluation of a Cardiovascular Systems Model for Design and Analysis of Hemodynamic Safety Studies. Pharmaceutics, 15(4), 1175. https://doi.org/10.3390/pharmaceutics15041175

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

Article Metrics

Back to TopTop