Pharmacometrics and systems pharmacology for metabolic bone diseases

Application of Quantitative Systems Pharmacology Approaches to Support Pediatric Labeling in Rare Diseases

Quantitative Systems Pharmacology (QSP) models offer a promising approach to extrapolate drug efficacy across different patient populations, particularly in rare diseases. Unlike conventional empirical models, QSP models can provide a mechanistic understanding of disease progression and therapeutic response by incorporating current disease knowledge into the descriptions of biomarkers and clinical endpoints. This allows for a holistic representation of the disease and drug response. The mechanistic nature of QSP models is well suited to pediatric extrapolation concepts, providing a quantitative method to assess disease and drug response similarity between adults and pediatric patients. The application of a QSP-based assessment of the disease and drug similarity in adult and pediatric patients in the clinical development program of olipudase alfa, a treatment for Acid Sphingomyelinase Deficiency (ASMD), illustrates the potential of this approach.

A review of mathematical modeling of bone remodeling from a systems biology perspective

Bone remodeling is an essential, delicately balanced physiological process of coordinated activity of bone cells that remove and deposit new bone tissue in the adult skeleton. Due to the complex nature of this process, many mathematical models of bone remodeling have been developed. Each of these models has unique features, but they have underlying patterns. In this review, the authors highlight the important aspects frequently found in mathematical models for bone remodeling and discuss how and why these aspects are included when considering the physiology of the bone basic multicellular unit, which is the term used for the collection of cells responsible for bone remodeling. The review also emphasizes the view of bone remodeling from a systems biology perspective. Understanding the systemic mechanisms involved in remodeling will help provide information on bone pathology associated with aging, endocrine disorders, cancers, and inflammatory conditions and enhance systems pharmacology. Furthermore, some features of the bone remodeling cycle and interactions with other organ systems that have not yet been modeled mathematically are discussed as promising future directions in the field.

Circular RNA circ-3626 promotes bone formation by modulating the miR-338-3p/Runx2 axis

OBJECTIVE

Disorders of bone homeostasis are the key factors leading to metabolic bone disease, such as senile osteoporosis, which is characterized by age-related bone loss. Bone marrow stromal cells (BMSCs) possess high osteogenic capacity which has been regarded as a practical approach to preventing bone loss. Previous studies have shown that the osteogenic differentiation ability of BMSCs is significantly decreased in senile osteoporosis. Recently, circular RNAs (circRNAs) have been regarded as critical regulators in controlling the osteogenic differentiation of BMSCs by sponging microRNAs (miRNAs). Our study aimed to discover new and critical osteogenesis-related circRNAs that can promote bone formation in senile osteoporosis.

METHODS

We detected the dysregulated circRNAs of BMSCs upon osteogenic differentiation induction and identified the critical osteogenic circRNA (circ-3626). The relationship between circ-3626 and osteoporosis was further verified in clinical bone samples and aged mice by qPCR. Moreover, circ-3626 AAV was constructed to examine the osteogenic effect of circ-3626 on bone formation via using Micro-CT, double calcein labeling, and the three-point bending tests. Bioinformatics analysis, Luciferase report gene assays, FISH, RNA pull-down, qPCR, Western Blots, and alizarin red staining assay explore the effects and mechanisms of circ-3626 on osteogenic differentiation of BMSCs.

RESULTS

Circ-3626 was identified as a pivotal osteogenesis-related circRNA via RNA sequencing. The results of alizarin red staining, Western blots, and qPCR assays suggest that overexpressing circ-3626 dramatically accelerates the osteogenic capability of BMSCs. Furthermore, the bone repair capability of aging mice could be significantly improved by circ-3626 AAV treatment. Micro RNA miR-338-3p was identified as the downstream target of circ-3626. Overexpression of circ-3626 increases the expression of Runx2 by sponging miR-338-3p, thereby promoting the osteogenic differentiation of BMSCs by upregulating the expression of osteogenic genes. In addition, Western blots, and qPCR assays suggest circ-3626 AAV treatment promote the expression of Runx2 and osteogenic marker genes.

CONCLUSION

Thus, we demonstrate that circ-3626 plays a pivotal role in promoting bone formation through the miR-338-3p/Runx2 axis and may provide new strategies for preventing and treating the bone loss of senile osteoporosis.

Pharmacometrics: The Already-Present Future of Precision Pharmacology

The use of mathematical modeling to represent, analyze, make predictions or providing information on data obtained in drug research and development has made pharmacometrics an area of great prominence and importance. The main purpose of pharmacometrics is to provide information relevant to the search for efficacy and safety improvements in pharmacotherapy. Regulatory agencies have adopted pharmacometrics analysis to justify their regulatory decisions, making those decisions more efficient. Demand for specialists trained in the field is therefore growing. In this review, we describe the meaning, history, and development of pharmacometrics, analyzing the challenges faced in the training of professionals. Examples of applications in current use, perspectives for the future, and the importance of pharmacometrics for the development and growth of precision pharmacology are also presented.

A Mathematical Model of In Vitro Cellular Uptake of Zoledronic Acid and Isopentenyl Pyrophosphate Accumulation

The mevalonate pathway is an attractive target for many areas of research, such as autoimmune disorders, atherosclerosis, Alzheimer’s disease and cancer. Indeed, manipulating this pathway results in the alteration of malignant cell growth with promising therapeutic potential. There are several pharmacological options to block the mevalonate pathway in cancer cells, one of which is zoledronic acid (ZA) (an N-bisphosphonate (N-BP)), which inhibits the farnesyl pyrophosphate (FPP) synthase enzyme, inducing cell cycle arrest, apoptosis, inhibition of protein prenylation, and cholesterol reduction, as well as leading to the accumulation of isopentenyl pyrophosphate (IPP). We extrapolated the data based on two independently published papers that provide numerical data on the uptake of zoledronic acid (ZA) and the accumulation of IPP (Ag) and its isomer over time by using in vitro human cell line models. Two different mathematical models for IPP kinetics are proposed. The first model (Model 1) is a simpler ordinary differential equation (ODE) compartmental system composed of 3 equations with 10 parameters; the second model (Model 2) is a differential algebraic equation (DAE) system with 4 differential equations, 1 algebraic equation and 13 parameters incorporating the formation of the ZA+enzyme+Ag complex. Each of the two models aims to describe two different experimental situations (continuous and pulse experiments) with the same ZA kinetics. Both models fit the collected data very well. With Model 1, we obtained a prevision accumulation of IPP after 24 h of 169.6 pmol/mgprot/h with an IPP decreasing rate per (pmol/mgprot) of ZA (kXGZ) equal to 13.24/h. With Model 2, we have comprehensive kinetics of IPP upon ZA treatment. We calculate that the IPP concentration was equal to 141.6 pmol/mgprot/h with a decreasing rate/percentage of 0.051 (kXGU). The present study is the first to quantify the influence of ZA on the pharmacodynamics of IPP. While still incorporating a small number of parameters, Model 2 better represents the complexity of the biological behaviour for calculating the IPP produced in different situations, such as studies on γδ T cell-based immunotherapy. In the future, additional clinical studies are warranted to further evaluate and fine-tune dosing approaches.

A model-based analysis to guide gonadotropin-releasing hormone receptor antagonist use for management of endometriosis

Aims

To identify linzagolix doses, an oral GnRH receptor antagonist, that effectively lower oestradiol (E2) to relieve endometriosis‐related pelvic pain without compromising bone health.

Methods

Integrated statistical, pharmacokinetic–pharmacodynamic and systems pharmacology models were developed from Phase 1 and 2 clinical trial data in healthy volunteers and patients, receiving linzagolix 25–200 mg daily or placebo, and analysed simultaneously. The main outcome measures were pelvic pain scores for dysmenorrhoea, nonmenstrual pelvic pain (NMPP), uterine bleeding and lumbar spine bone mineral density (BMD).

Results

Linzagolix pharmacokinetics were described by a 2‐compartment model with sequential zero/first‐order absorption process (CL/F: 0.422 L/h). E2 changes over time were well described as a function of linzagolix 24‐hour AUC (AUC₅₀: 1.68 × 10⁵ ng h/mL). For a Caucasian reference patient, a change in E2 from 50–20 pg/mL at 24 weeks increased the odds of relief of dysmenorrhoea 1.33‐fold and NMPP 1.07‐fold (95% CI: 1.22–1.47 and 1.02–1.12, respectively) and decreased bleeding days by 1.55 (95% CI: 1.39–1.72). A previously validated quantitative systems pharmacology BMD model was adjusted to the clinical data. The mean week 24 lumbar spine BMD change from baseline ranged from −0.092% in the 50 mg dose, −1.30% in the 100 mg dose group and −2.67% in the 200 mg dose group.

Discussion

The previously‐reported E2 target range (20–50 pg/mL) to balance efficacy and safety endpoints was confirmed. Linzagolix once daily doses between 75–125 mg daily were expected to meet endometriosis‐associated pain, efficacy, and BMD loss targets in Caucasian patients.

On the pharmacological evaluation of bisphosphonates in humans

One of the key parameters for a successful treatment with any drug is the use of an optimal dose regimen. Bisphosphonates (BPs) have been in clinical use for over five decades and during this period clinical pharmacokinetic (PK) and pharmacodynamic (PD) evaluations have been instrumental for the identification of optimal dose regimens in patients. Ideal clinical PK and PD studies help drug developers explain variability in responses and enable the identification of a dose regimen with an optimal effect. PK and PD studies of the unique and rather complex pharmacological properties of BPs also help determine to a significant extent ideal dosing for these drugs. Clinical PK and PD evaluations of BPs preferably use study designs and assays that enable the assessment of both short- (days) and long-term (years) presence and effect of these drugs in patients. BPs are mainly used for metabolic bone diseases because they inhibit osteoclast-mediated bone resorption and the best way to quantify their effects in humans is therefore by measuring biochemical markers of bone resorption in serum and urine. In these very same samples BP concentrations can also be measured. Short-term serum and urine data after both intravenous (IV) and oral administration enable the assessment of oral bioavailability as well as the amount of BP delivered to the skeleton. Longer-term data provide information on the anti-resorptive effect as well as the elimination of the BP from the skeleton. Using PK-PD models to mathematically link the anti-resorptive action of the BPs to the amount of BP at the skeleton provides a mechanism-based explanation of the pattern of bone resorption during treatment. These models have been used successfully during the clinical development of BPs. Newer versions of such models, which include systems pharmacology and disease progression models, are more comprehensive and include additional PD parameters such as BMD and fracture risk. Clinical PK and PD studies of BPs have been useful for the identification of optimal dose regimens for metabolic bone diseases. These analyses will also continue to be important for newer research directions, such as BP use in the delivery of other drugs to the bone to better treat bone metastases and bone infections, as well as the potential benefit of BPs at non-skeletal targets for the prevention and treatments of soft tissue cancers, various fibroses, and other cardiovascular and neurodegenerative diseases, and reduction in mortality and extension of lifespan.

On precision dosing of oral small molecule drugs in oncology

Personalization of oral small molecule anticancer drug doses based on individual patient blood drug levels, also known as therapeutic drug monitoring (TDM), has the potential to significantly improve the effectiveness of treatment by maximizing drug efficacy and minimize toxicity. However, this option has not yet been widely embraced by the oncology community. Some reasons for this include increased logistical complexity of dose individualization, the lack of clinical laboratories that measure small molecule drug concentrations in support of patient care, and the lack of reimbursement of costs. However, the main obstacle may be the lack of studies clearly demonstrating that monitoring of oral small molecule anticancer drug levels actually improves clinical outcomes. Without unequivocal evidence in support of TDM-guided dose individualization, especially demonstration of improved survival with TDM in randomized controlled trials, wide acceptance of this approach by oncologists and reimbursement by insurance companies is unlikely, and patients may continue to suffer as a result of receiving incorrect drug doses. This article reviews the current status of TDM of oral small molecule drugs in oncology and intends to provide strategic insights into the design of studies for evaluating the utility of TDM in this clinical context.

The Role of Osteoprotegerin in Vascular Calcification and Bone Metabolism: The Basis for Developing New Therapeutics

Osteoporosis (OP) and cardiovascular diseases (CVD) are both important causes of mortality and morbidity in aging patients. There are common mechanisms underlying the regulation of bone remodeling and the development of smooth muscle calcification; a temporal relationship exists between osteoporosis and the imbalance of mineral metabolism in the vessels. Vascular calcification appears regulated by mechanisms that include both inductive and inhibitory processes. Multiple factors are implicated in both bone and vascular metabolism. Among these factors, the superfamily of tumor necrosis factor (TNF) receptors including osteoprotegerin (OPG) and its ligands has been established. OPG is a soluble decoy receptor for receptor activator of nuclear factor-kB ligand (RANKL) and TNF-related apoptosis-inducing ligand (TRAIL). OPG binds to RANKL and TRAIL, and inhibits the association with their receptors, which have been labeled as the receptor activator of NF-kB (RANK). Sustained release of OPG from vascular endothelial cells (ECs) has been demonstrated in response to inflammatory proteins and cytokines, suggesting that OPG/RANKL/RANK system plays a modulatory role in vascular injury and inflammation. For the development of potential therapeutic strategies targeting vascular calcification, critical consideration of the implications for bone metabolism must be taken into account to prevent potentially detrimental effects to bone metabolism.

Pharmacometrics and systems pharmacology for metabolic bone diseases

Mathematical modelling and simulation (M&S) of drug concentrations, pharmacologic effects and the (patho)physiologic systems within which they interact can be powerful tools for the preclinical, translational and clinical development of drugs. Indeed, the Prescription Drug User Fee Act (PDUFA VI), incorporated as part of the FDA Reauthorization Act of 2017 (FDARA), highlights the goal of advancing model-informed drug development (MIDD). MIDD can benefit development across many drug classes, including for metabolic bone diseases such as osteoporosis, cancer-related and numerous rare metabolic bone diseases; conditions characterized by significant morbidity and mortality. A drought looms in terms of the availability of new drugs to better treat these devastating diseases. This review provides an overview of several M&S approaches ranging from simple pharmacokinetic to integrated pharmacometric and systems pharmacology modelling. Examples are included to illustrate the use of these approaches during the development of several drugs for metabolic bone diseases such as bisphosphonates, denosumab, teriparatide and sclerostin inhibitors (romosozumab and blosozumab).