Pharmacokinetic-pharmacodynamic modelling of acute N-terminal pro B-type natriuretic peptide after doxorubicin infusion in breast cancer

Mathematical modeling of cardio-oncology: Modeling the systemic effects of cancer therapeutics on the cardiovascular system

Cardiotoxicity is a common side-effect of many cancer therapeutics; however, to-date there has been very little push to understand the mechanisms underlying this group of pathologies. This has led to the emergence of cardio-oncology, a field of medicine focused on understanding the effects of cancer and its treatment on the human heart. Here, we describe how mechanistic modeling approaches have been applied to study open questions in the cardiovascular system and how these approaches are being increasingly applied to advance knowledge of the underlying effects of cancer treatments on the human heart. A variety of mechanistic, mathematical modeling techniques have been applied to explore the link between common cancer treatments, such as chemotherapy, radiation, targeted therapy, and immunotherapy, and cardiotoxicity, nevertheless there is limited coverage in the different types of cardiac dysfunction that may be associated with these treatments. Moreover, cardiac modeling has a rich heritage of mathematical modeling and is well suited for the further development of novel approaches for understanding the cardiotoxicities associated with cancer therapeutics. There are many opportunities to combine mechanistic, bottom-up approaches with data-driven, top-down approaches to improve personalized, precision oncology to better understand, and ultimately mitigate, cardiac dysfunction in cancer patients.

New Insights in the Era of Clinical Biomarkers as Potential Predictors of Systemic Therapy-Induced Cardiotoxicity in Women with Breast Cancer: A Systematic Review

Cardiotoxicity induced by breast cancer therapies is a potentially serious complication associated with the use of various breast cancer therapies. Prediction and better management of cardiotoxicity in patients receiving chemotherapy is of critical importance. However, the management of cancer therapy-related cardiac dysfunction (CTRCD) lacks clinical evidence and is based on limited clinical studies.

AIM

To provide an overview of existing and potentially novel biomarkers that possess a promising predictive value for the early and late onset of CTRCD in the clinical setting.

METHODS

A systematic review of published studies searching for promising biomarkers for the prediction of CTRCD in patients with breast cancer was undertaken according to PRISMA guidelines. A search strategy was performed using PubMed, Google Scholar, and Scopus for the period 2013-2023. All subjects were >18 years old, diagnosed with breast cancer, and received breast cancer therapies.

RESULTS

The most promising biomarkers that can be used for the development of an alternative risk cardiac stratification plan for the prediction and/or early detection of CTRCD in patients with breast cancer were identified.

CONCLUSIONS

We highlighted the new insights associated with the use of currently available biomarkers as a standard of care for the management of CTRCD and identified potentially novel clinical biomarkers that could be further investigated as promising predictors of CTRCD.

Population pharmacokinetics of doxorubicin: A systematic review

Because of the high interindividual pharmacokinetic variability, several population pharmacokinetic (PopPK) models of doxorubicin (DOX) were developed to characterize factors influencing such variability. However, significant predictors for DOX pharmacokinetics identified using PopPK models varied across studies. Thus, this review aims to summarize PopPK models of DOX and its metabolites (if any) as well as significant covariates influencing DOX (and its metabolites) pharmacokinetic variability. A systematic search from PubMed, CINAHL Complete, Science Direct, and SCOPUS databases identified 503 studies. Of these, 16 studies met the inclusion criteria and were included in this review. DOX pharmacokinetics was described with two- or three-compartment models. Most studies found a significant increase in DOX clearance with an increase in body surface area from the median value of 1.8 m² . Moreover, this review identified that while a 10-year increase in patient age resulted in a decrease in DOX clearance in adults and the elderly, younger children had lower DOX clearance compared to older children. Further, low DOX exposure was observed in pregnant women, and thus dosage adjustment is required. Concerning model applicability, predictive performance assessment of these published models should be performed before implementing such models in clinical practice.

Population Pharmacokinetics and Pharmacodynamics of Vericiguat in Patients with Heart Failure and Reduced Ejection Fraction

Background

Vericiguat, a stimulator of soluble guanylate cyclase, has been developed as a first-in-class therapy for worsening chronic heart failure in adults with left ventricular ejection fraction < 45%.

Objective

The objective of this article was to characterize the pharmacokinetics and pharmacokinetic variability of vericiguat combined with guideline-directed medical therapy (standard of care), and identify exposure–response relationships for safety (hemodynamics) and pharmacodynamic markers of efficacy (N-terminal pro-B-type natriuretic peptide concentration [NT-proBNP]) in patients with heart failure and left ventricular ejection fraction < 45% in the SOCRATES-REDUCED study (NCT01951625).

Methods

Vericiguat and NT-proBNP plasma concentrations in 454 and 432 patients in SOCRATES-REDUCED, respectively, were analyzed using nonlinear mixed-effects modeling.

Results

Vericiguat pharmacokinetics were well described by a one-compartment model with apparent clearance, apparent volume of distribution, and absorption rate constant. Age, bodyweight, plasma bilirubin, and creatinine clearance were identified as significant covariates on apparent clearance; sex and bodyweight on apparent volume of distribution; and bodyweight and plasma albumin level on absorption rate constant. Pharmacokinetic/pharmacodynamic analysis showed initial minor and transient effects of vericiguat on blood pressure with low clinical impact. There were no changes in heart rate following initial or repeated vericiguat administration. An exposure-dependent and time-dependent turnover pharmacokinetic/pharmacodynamic model for NT-proBNP described production and elimination rates and an demonstrated exposure-dependent reduction in [NT-proBNP] by vericiguat plus standard of care compared with placebo plus standard of care. This effect was dependent on baseline [NT-proBNP].

Conclusions

Vericiguat has predictable pharmacokinetics, with no long-term effects on blood pressure in patients with heart failure and left ventricular ejection fraction < 45%. A pharmacokinetic/pharmacodynamic model described a vericiguat exposure-dependent reduction of NT-proBNP.

Clinical Trial Identifier

NCT01951625.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-021-01024-y.

Drug Exposure to Establish Pharmacokinetic-Response Relationships in Oncology

In the oncology field, understanding the relationship between the dose administered and the exerted effect is particularly important because of the narrow therapeutic index associated with anti-cancer drugs and the high interpatient variability. Therefore, in this review, we provide a critical perspective of the different methods of characterising treatment exposure in the oncology setting. The increasing number of modelling applications in oncology reflects the applicability and the impact of pharmacometrics on all phases of the drug development process and patient management as well. Pharmacometric modelling is a worthy component within the current paradigm of model-based drug development, but pharmacometric modelling techniques are also accessible for the clinician in the optimisation of current oncology therapies. Consequently, the application of population models in a hospital setting by generating close collaborations between physicians and pharmacometricians is highly recommended, providing a systematic means of developing and assessing model-based metrics as 'drivers' for various responses to treatments, which can then be evaluated as predictors for treatment success. Characterising the key determinants of variability in exposure is of particular importance for anticancer agents, as efficacy and toxicity are associated with exposure. We present the different strategies to describe and predict drug exposure that can be applied depending on the data available, with the objective of obtaining the most useful information in the patients' favour throughout the full drug cycle. Therefore, the objective of the present article is to review the different approaches used to characterise a patient's exposure to oncology drugs, which will result in a better understanding of the time course of the response and the magnitude of interpatient variability.

In vitro to Clinical Translation of Combinatorial Effects of Doxorubicin and Abemaciclib in Rb-Positive Triple Negative Breast Cancer: A Systems-Based Pharmacokinetic/Pharmacodynamic Modeling Approach

Background

Doxorubicin (DOX) and its pegylated liposomal formulation (L_DOX) are the standard of care for triple-negative breast cancer (TNBC). However, resistance to DOX often occurs, motivating the search for alternative treatment approaches. The retinoblastoma protein (Rb) is a potential pharmacological target for TNBC treatment since its expression has been associated with resistance to DOX-based therapy.

Methods

DOX (0.01–20 μM) combination with abemaciclib (ABE, 1–6 μM) was evaluated over 72 hours on Rb-positive (MDA-MB-231) and Rb-negative (MDA-MB-468) TNBC cells. Combination indices (CI) for DOX+ABE were calculated using Compusyn software. The TNBC cell viability time-course and fold-change from the control of phosphorylated-Rb (pRb) protein expression were measured with CCK8-kit and enzyme-linked immunosorbent assay. A cell-based pharmacodynamic (PD) model was developed, where pRb protein dynamics drove cell viability response. Clinical pharmacokinetic (PK) models for DOX, L_DOX, and ABE were developed using data extracted from the literature. After scaling cancer cell growth to clinical TNBC tumor growth, the time-to-tumor progression (TTP) was predicted for human dosing regimens of DOX, ABE, and DOX+ABE.

Results

DOX and ABE combinations were synergistic (CI<1) in MDA-MB-231 and antagonistic (CI>1) in MDA-MB-468. The maximum inhibitory effects (Imax) for both drugs were set to one. The drug concentrations producing 50% of Imax for DOX and ABE were 0.565 and 2.31 μM (MDA-MB-231) and 0.121 and 1.61 μM (MDA-MB-468). The first-orders rate constants of abemaciclib absorption (k_a) and doxorubicin release from L_DOX (k_Rel) were estimated at 0.31 and 0.013 h⁻¹. Their linear clearances were 21.7 (ABE) and 32.1 L/h (DOX). The estimated TTP for intravenous DOX (75 mg/m2 every 21 days), intravenous L_DOX (50 mg/m² every 28 days), and oral ABE (200 mg twice a day) were 125, 31.2, and 8.6 days shorter than drug-free control. The TTP for DOX+ABE and L_DOX+ABE were 312 days and 47.5 days shorter than control, both larger than single-agent DOX, suggesting improved activity with the DOX+ABE combination.

Conclusion

The developed translational systems-based PK/PD model provides an in vitro-to-clinic modeling platform for DOX+ABE in TNBC. Although model-based simulations suggest improved outcomes with combination over monotherapy, tumor relapse was not prevented with the combination. Hence, DOX+ABE may not be an effective treatment combination for TNBC.

Drug-Eluting Hydrophilic Coating Modification of Intraocular Lens via Facile Dopamine Self-Polymerization for Posterior Capsular Opacification Prevention

Posterior capsular opacification (PCO) is the most important complication in cataract phacoemulsification and intraocular lens (IOL) implantation surgery, mainly stemming from the adhesion, proliferation, and transdifferentiation of the postsurgically residual lens epithelial cells (LECs). Previous investigations mainly focused on the hydrophilic surface modification of the IOLs for PCO prevention, such as heparinization. However, the long-term clinical investigations show that there is no significant difference between pristine and heparinized IOLs. In the present study, a synergetic coating with properties of drug-eluting and hydrophilicity was designed and modified onto the IOL surface via facile dopamine self-polymerization. The antiproliferative drug doxorubicin (DOX) was loaded when a polydopamine (PDA) coating was formed on the IOL surface. The hydrophilic 2-methacryloyloxyethyl phosphorylcholine (MPC) could be subsequently grafted onto the drug-loaded PDA coating surface easily. The hydrophilic outer layer could slow down drug-eluting from underneath the drug-incorporated coating. In vitro and in vivo investigations demonstrated that such multifunctionalized coating-modified IOLs could not only thoroughly and effectively prevent PCO development by induced cell apoptosis but also render safety and biocompatibility to the surrounding tissues.

Multiscale and Translational Quantitative Systems Toxicology, Pharmacokinetic-Toxicodynamic Modeling Analysis for Assessment of Doxorubicin-Induced Cardiotoxicity

Dose-dependent life-threatening doxorubicin-induced cardiotoxicity (DIC) is a major clinical challenge that needs to be addressed. Here, we developed an integrated multiscale and translational quantitative systems toxicology and pharmacokinetic-toxicodynamic (QST-PK/TD) model for optimization of doxorubicin dosing regimens for early monitoring and minimization of DIC. A QST model was established by exposing human cardiomyocytes, AC16 cells, to doxorubicin over a time course, and measuring the dynamics of intracellular signaling proteins, AC16 cell viability and released biomarkers of cardiomyocyte injury such as the B-type natriuretic peptide (BNP). Experiments were scaled up to a three-dimensional and dynamic (3DD) cell culture system to evaluate DIC under various dosing regimens. The PK determinants of doxorubicin influencing DIC were identified in vitro and then translated to the in vivo setting through hybrid physiologically based PK (PBPK)/TD models using preclinical- and clinical-level data extracted from literature. The developed cellular-level QST model captured well the observed dynamics of intracellular proteins, AC16 cell viability and BNP kinetics. In the 3DD setting, dose fractionation of doxorubicin displayed a significant reduction in cardiotoxicity compared to single intravenous doses with equal exposure, implying doxorubicin peak concentrations as the PK determinant for DIC. The in vivo hybrid PBPK/TD models captured well doxorubicin PK and DIC. Peak doxorubicin concentrations correlated well with acute DIC for dose-fractionated regimens, while maximum 48-h moving average concentrations correlated with DIC for dose-fractionated and long-term infusion regimens in vivo. The developed multiscale and translational QST-PK/TD modeling platform may serve as an in silico tool for assessment of early toxicity and/or efficacy of developmental drugs in vitro.

Is there scope for better individualisation of anthracycline cancer chemotherapy?

Anthracyclines are used to treat solid and haematological cancers, particularly breast cancers, lymphomas and childhood cancers. Myelosuppression and cardiotoxicity are the primary toxicities that limit treatment duration and/or intensity. Cardiotoxicity, particularly heart failure, is a leading cause of morbidity and mortality in cancer survivors. Cumulative anthracycline dose is a significant predictor of cardiotoxicity risk, suggesting a role for anthracycline pharmacokinetic variability. Population pharmacokinetic modelling in children has shown that doxorubicin clearance in the very young is significantly lower than in older children, potentially contributing to their higher risk of cardiotoxicity. A model of doxorubicin clearance based on body surface area and age offers a patient-centred dose-adjustment strategy that may replace the current disparate initial-dose selection tools, providing a rational way to compensate for pharmacokinetic variability in children aged <7 years. Population pharmacokinetic models in adults have not adequately addressed older ages, obesity, hepatic and renal dysfunction, and potential drug-drug interactions to enable clinical application. Although candidate gene and genome-wide association studies have investigated relationships between genetic variability and anthracycline pharmacokinetics or clinical outcomes, there have been few clinically significant reproducible associations. Precision-dosing of anthracyclines is currently hindered by lack of clinically useful pharmacokinetic targets and models that predict cumulative anthracycline exposures. Combined with known risk factors for cardiotoxicity, the use of advanced echocardiography and biomarkers, future validated pharmacokinetic targets and predictive models could facilitate anthracycline precision dosing that truly maximises efficacy and provides individualised early intervention with cardioprotective therapies in patients at risk of cardiotoxicity.

Fisetin, a plant flavonoid ameliorates doxorubicin-induced cardiotoxicity in experimental rats: the decisive role of caspase-3, COX-II, cTn-I, iNOs and TNF-α

Doxorubicin (DOX) is a widely used anthracycline antibiotic for the management of carcinoma. However, it is associated with cardiotoxicity. Fisetin is a plant flavonoid reported to have anti-inflammatory and antiapoptotic potential. To evaluate the cardioprotective potential of fisetin in DOX-induced cardiotoxicity in experimental rats. Sprague-Dawley rats were pre-treated with either fisetin (10, 20 and 40 mg/kg) or sitagliptin (10 mg/kg, p.o.) for 7 days. Cardiac toxicity was induced in rats (except the normal group) by doxorubicin (15 mg/kg i.p.) on 8th day. Various behavioral, biochemical, molecular and histological parameters were assessed in cardiac tissue. DOX-induced alterations in electrocardiographic, hemodynamic and left ventricular function were significantly (p < 0.05) inhibited by fisetin (20 and 40 mg/kg) treatment. Fisetin significantly decrease (p < 0.05) DOX-induced elevated serum CK-MB, LDH, AST, ALT and ALP levels. DOX-induced elevated cardiac oxido-nitrosative (SOD, GSH, MDA and NO) was significantly inhibited (p < 0.05) by fisetin. Up-regulated cardiac caspase-3, COX-II, cTn-I, iNOs, TNF-α, and IL-1β mRNA, as well as protein expressions were significantly decreased (p < 0.05) by fisetin treatment. It also significantly (p < 0.05) attenuated DOX-induced histopathological alterations in cardiac tissue. In conclusion, the fisetin exerts its cardioprotective potential against DOX-induced toxicity via inhibition of multiple pathways including oxidative stress (SOD, GSH, MDA and NO), inflammation (COX-II, TNF-α, and IL-1β), and apoptosis (Caspase-3). Therefore, fisetin can be considered as a potential cardioprotective agent during the management of carcinoma using doxorubicin anthracyclines.

Pharmacokinetic-pharmacodynamic modelling of acute N-terminal pro B-type natriuretic peptide after doxorubicin infusion in breast cancer

AIMS

The aim of the present study was to develop a pharmacokinetic-pharmacodynamic (PK-PD) model to characterize the relationship between plasma doxorubicin and N-terminal pro B-type natriuretic peptide (NT-proBNP) concentrations within 48 h of doxorubicin treatment.

METHODS

The study enrolled 17 female patients with stages 1-3 breast cancer and receiving adjuvant doxorubicin (60 mg m(-2) ) and cyclophosphamide (600 mg m(-2) ) every 14 days for four cycles. In two consecutive cycles, plasma concentrations of doxorubicin, doxorubicinol, troponin and NT-proBNP were collected before infusion, and up to 48 h after the end of doxorubicin infusion. Nonlinear mixed-effects modelling was used to describe the PK-PD relationship of doxorubicin and NT-proBNP.

RESULTS

A three-compartment parent drug with a one-compartment metabolite model best described the PK of doxorubicin and doxorubicinol. Troponin concentrations remained similar to baseline. An indirect PD model with transit compartments best described the relationship of doxorubicin exposure and acute NT-proBNP response. Estimated PD parameters were associated with large between-subject variability (total assay variability 38.8-73.9%). Patient clinical factors, including the use of enalapril, were not observed to be significantly associated with doxorubicin PK or NT-proBNP PD variability.

CONCLUSION

The relationship between doxorubicin concentration and the acute NT-proBNP response was successfully described with a population PK-PD model. This model will serve as a valuable framework for future studies to identify clinical factors associated with the acute response to doxorubicin. Future studies are warranted to examine the relationship between this acute response and subsequent heart failure. Should such a relationship be established, this model could provide useful information on patients' susceptibility to doxorubicin-induced long-term cardiotoxicity.