Leveraging model-informed approaches for drug discovery and development in the cardiovascular space.

Cardiovascular disease remains a significant global health burden, and development of cardiovascular drugs in the current regulatory environment often demands large and expensive cardiovascular outcome trials. Thus, the use of quantitative pharmacometric approaches which can help enable early Go/No Go decision making, ensure appropriate dose selection, and increase the likelihood of successful clinical trials, have become increasingly important to help reduce the risk of failed cardiovascular outcomes studies. In addition, cardiovascular safety is an important consideration for many drug development programs, whether or not the drug is designed to treat cardiovascular disease; modeling and simulation approaches also have utility in assessing risk in this area. Herein, examples of modeling and simulation applied at various stages of drug development, spanning from the discovery stage through late-stage clinical development, for cardiovascular programs are presented. Examples of how modeling approaches have been utilized in early development programs across various therapeutic areas to help inform strategies to mitigate the risk of cardiovascular-related adverse events, such as QTc prolongation and changes in blood pressure, are also presented. These examples demonstrate how more informed drug development decisions can be enabled by modeling and simulation approaches in the cardiovascular area.

Preclinical and translational evaluation of coagulation factor IXa as a novel therapeutic target.

The benefits of novel oral anticoagulants are hampered by bleeding. Since coagulation factor IX (fIX) lies upstream of fX in the coagulation cascade, and intermediate levels have been associated with reduced incidence of thrombotic events, we evaluated the viability of fIXa as an antithrombotic target. We applied translational pharmacokinetics/pharmacodynamics (PK/PD) principles to predict the therapeutic window (TW) associated with a selective small molecule inhibitor (SMi) of fIXa, compound 1 (CPD1, rat fIXa inhibition constant (Ki, 21 nmol/L) relative to clinically relevant exposures of apixaban (rat fXa Ki 4.3 nmol/L). Concentrations encompassing the minimal clinical plasma concentration (C min) of the 5 mg twice daily (BID) dose of apixaban were tested in rat arteriovenous shunt (AVS/thrombosis) and cuticle bleeding time (CBT) models. An I max and a linear model were used to fit clot weight (CW) and CBT. The following differences in biology were observed: (1) antithrombotic activity and bleeding increased in parallel for apixaban, but to a lesser extent for CPD1 and (2) antithrombotic activity occurred at high (>99%) enzyme occupancy (EO) for fXa or moderate (>65% EO) for fIXa. translational PK/PD analysis indicated that noninferiority was observed for concentrations of CPD1 that provided between 86% and 96% EO and that superior TW existed between 86% and 90% EO. These findings were confirmed in a study comparing short interfering (si)RNA-mediated knockdown (KD) modulation of fIX and fX mRNA. In summary, using principles of translational biology to relate preclinical markers of efficacy and safety to clinical doses of apixaban, we found that modulation of fIXa can be superior to apixaban.

Multiscale prediction of patient-specific platelet function under flow.

During thrombotic or hemostatic episodes, platelets bind collagen and release ADP and thromboxane A(2), recruiting additional platelets to a growing deposit that distorts the flow field. Prediction of clotting function under hemodynamic conditions for a patient's platelet phenotype remains a challenge. A platelet signaling phenotype was obtained for 3 healthy donors using pairwise agonist scanning, in which calcium dye-loaded platelets were exposed to pairwise combinations of ADP, U46619, and convulxin to activate the P2Y(1)/P2Y(12), TP, and GPVI receptors, respectively, with and without the prostacyclin receptor agonist iloprost. A neural network model was trained on each donor's pairwise agonist scanning experiment and then embedded into a multiscale Monte Carlo simulation of donor-specific platelet deposition under flow. The simulations were compared directly with microfluidic experiments of whole blood flowing over collagen at 200 and 1000/s wall shear rate. The simulations predicted the ranked order of drug sensitivity for indomethacin, aspirin, MRS-2179 (a P2Y(1) inhibitor), and iloprost. Consistent with measurement and simulation, one donor displayed larger clots and another presented with indomethacin resistance (revealing a novel heterozygote TP-V241G mutation). In silico representations of a subject's platelet phenotype allowed prediction of blood function under flow, essential for identifying patient-specific risks, drug responses, and novel genotypes.

Systems biology of coagulation initiation: kinetics of thrombin generation in resting and activated human blood.

Blood function defines bleeding and clotting risks and dictates approaches for clinical intervention. Independent of adding exogenous tissue factor (TF), human blood treated in vitro with corn trypsin inhibitor (CTI, to block Factor XIIa) will generate thrombin after an initiation time (T(i)) of 1 to 2 hours (depending on donor), while activation of platelets with the GPVI-activator convulxin reduces T(i) to ∼20 minutes. Since current kinetic models fail to generate thrombin in the absence of added TF, we implemented a Platelet-Plasma ODE model accounting for: the Hockin-Mann protease reaction network, thrombin-dependent display of platelet phosphatidylserine, VIIa function on activated platelets, XIIa and XIa generation and function, competitive thrombin substrates (fluorogenic detector and fibrinogen), and thrombin consumption during fibrin polymerization. The kinetic model consisting of 76 ordinary differential equations (76 species, 57 reactions, 105 kinetic parameters) predicted the clotting of resting and convulxin-activated human blood as well as predicted T(i) of human blood under 50 different initial conditions that titrated increasing levels of TF, Xa, Va, XIa, IXa, and VIIa. Experiments with combined anti-XI and anti-XII antibodies prevented thrombin production, demonstrating that a leak of XIIa past saturating amounts of CTI (and not "blood-borne TF" alone) was responsible for in vitro initiation without added TF. Clotting was not blocked by antibodies used individually against TF, VII/VIIa, P-selectin, GPIb, protein disulfide isomerase, cathepsin G, nor blocked by the ribosome inhibitor puromycin, the Clk1 kinase inhibitor Tg003, or inhibited VIIa (VIIai). This is the first model to predict the observed behavior of CTI-treated human blood, either resting or stimulated with platelet activators. CTI-treated human blood will clot in vitro due to the combined activity of XIIa and XIa, a process enhanced by platelet activators and which proceeds in the absence of any evidence for kinetically significant blood borne tissue factor.

RGS/Gi2alpha interactions modulate platelet accumulation and thrombus formation at sites of vascular injury.

Although much is known about extrinsic regulators of platelet function such as nitric oxide and prostaglandin I(2) (PGI(2)), considerably less is known about intrinsic mechanisms that prevent overly robust platelet activation after vascular injury. Here we provide the first evidence that regulators of G-protein signaling (RGS) proteins serve this role in platelets, using mice with a G184S substitution in G(i2α) that blocks RGS/G(i2) interactions to examine the consequences of lifting constraints on G(i2)-dependent signaling without altering receptor:effector coupling. The results show that the G(i2α)(G184S) allele enhances platelet aggregation in vitro and increases platelet accumulation after vascular injury when expressed either as a global knock-in or limited to hematopoietic cells. Biochemical studies show that these changes occur in concert with an attenuated rise in cyclic adenosine monophosphate levels in response to prostacyclin and a substantial increase in basal Akt activation. In contrast, basal cyclic adenosine monophosphate (cAMP) levels, agonist-stimulated increases in [Ca(++)](i), Rap1 activation, and α-granule secretion were unaffected. Collectively, these observations (1) demonstrate an active role for RGS proteins in regulating platelet responsiveness, (2) show that this occurs in a pathway-selective manner, and (3) suggest that RGS proteins help to prevent unwarranted platelet activation as well as limiting the magnitude of the normal hemostatic response.

Pairwise agonist scanning predicts cellular signaling responses to combinatorial stimuli.

Prediction of cellular response to multiple stimuli is central to evaluating patient-specific clinical status and to basic understanding of cell biology. Cross-talk between signaling pathways cannot be predicted by studying them in isolation and the combinatorial complexity of multiple agonists acting together prohibits an exhaustive exploration of the complete experimental space. Here we describe pairwise agonist scanning (PAS), a strategy that trains a neural network model based on measurements of cellular responses to individual and all pairwise combinations of input signals. We apply PAS to predict calcium signaling responses of human platelets in EDTA-treated plasma to six different agonists (ADP, convulxin, U46619, SFLLRN, AYPGKF and PGE(2)) at three concentrations (0.1, 1 and 10 x EC(50)). The model predicted responses to sequentially added agonists, to ternary combinations of agonists and to 45 different combinations of four to six agonists (R = 0.88). Furthermore, we use PAS to distinguish between the phenotypic responses of platelets from ten donors. Training neural networks with pairs of stimuli across the dose-response regime represents an efficient approach for predicting complex signal integration in a patient-specific disease milieu.

A molecular signaling model of platelet phosphoinositide and calcium regulation during homeostasis and P2Y1 activation.

To quantify how various molecular mechanisms are integrated to maintain platelet homeostasis and allow responsiveness to adenosine diphosphate (ADP), we developed a computational model of the human platelet. Existing kinetic information for 77 reactions, 132 fixed kinetic rate constants, and 70 species was combined with electrochemical calculations, measurements of platelet ultrastructure, novel experimental results, and published single-cell data. The model accurately predicted: (1) steady-state resting concentrations for intracellular calcium, inositol 1,4,5-trisphosphate, diacylglycerol, phosphatidic acid, phosphatidylinositol, phosphatidylinositol phosphate, and phosphatidylinositol 4,5-bisphosphate; (2) transient increases in intracellular calcium, inositol 1,4,5-trisphosphate, and G(q)-GTP in response to ADP; and (3) the volume of the platelet dense tubular system. A more stringent test of the model involved stochastic simulation of individual platelets, which display an asynchronous calcium spiking behavior in response to ADP. Simulations accurately reproduced the broad frequency distribution of measured spiking events and demonstrated that asynchronous spiking was a consequence of stochastic fluctuations resulting from the small volume of the platelet. The model also provided insights into possible mechanisms of negative-feedback signaling, the relative potency of platelet agonists, and cell-to-cell variation across platelet populations. This integrative approach to platelet biology offers a novel and complementary strategy to traditional reductionist methods.

Determination of surface tissue factor thresholds that trigger coagulation at venous and arterial shear rates: amplification of 100 fM circulating tissue factor requires flow.

Protein microarrays presenting spots of collagen and lipidated tissue factor (TF) allowed a determination of the critical surface concentration of TF required to trigger coagulation under flow. Whole blood supplemented with corn trypsin inhibitor (to inhibit factor XIIa) was perfused over microarrays for 5 minutes. Immunofluorescence staining of platelet glycoprotein GPIbalpha and fibrin(ogen) revealed a critical TF concentration (EC50) of 3.6, 8.4, and 10.2 molecules-TF/microm2 at wall shear rates of 100, 500, and 1000 s(-1), respectively. For collagen arrays where only the center lane of spots (in the direction of flow) contained TF, a downstream distance of 14 mm was required for the thrombus to widen enough to reach across a 300-micrometer gap to the adjacent TF-free lanes of collagen spots, in agreement with numerical simulation. To investigate the effect of low levels of circulating TF, whole blood (+/-100 fM added TF) was tested under static and flow conditions. After 5 minutes, the addition of 100 fM TF to whole blood had negligible effect under static conditions, but caused a 2.5-fold increase in fibrin formation under flow. This report defines the threshold concentrations of surface TF required to trigger coagulation under flow.