Pan AMPK Activation Protects Tubules in Rat Ischemic Acute Kidney Injury.

Acute Kidney Injury (AKI) is characterized by an abrupt decline in kidney function and has been associated with excess risks of death, kidney disease progression, and cardiovascular events. The kidney has a high energetic demand with mitochondrial health being essential to renal function and damaged mitochondria has been reported across AKI subtypes. 5' adenosine monophosphate-activated protein kinase (AMPK) activation preserves cellular energetics through improvement of mitochondrial function and biogenesis when ATP levels are low such as under ischemia-induced AKI. We developed a selective potent small molecule pan AMPK activator, compound 1, and tested its ability to increase AMPK activity and preserve kidney function during ischemia/reperfusion injury in rats. A single administration of 1 caused sustained activation of AMPK for at least 24 hours, protected against acute tubular necrosis, and reduced clinical markers of tubular injury such as NephroCheck and Fractional Excretion of Sodium (FENa). Reduction in plasma creatinine and increased Glomerular Filtration Rate (GFR) indicated preservation of kidney function. Surprisingly, we observed a strong diuretic effect of AMPK activation associated with natriuresis both with and without AKI. Our findings demonstrate that activation of AMPK leads to protection of tubular function under hypoxic/ischemic conditions which holds promise as a potential novel therapeutic approach for AKI. Significance Statement No approved pharmacological therapies currently exist for acute kidney injury. We developed Compound 1 which dose-dependently activated AMPK in the kidney and protected kidney function and tubules after ischemic renal injury in the rat. This was accompanied by natriuresis in injured as well as uninjured rats.

Dynamic transition states of ErbB1 phosphorylation predicted by spatial stochastic modeling.

ErbB1 overexpression is strongly linked to carcinogenesis, motivating better understanding of erbB1 dimerization and activation. Recent single-particle-tracking data have provided improved measures of dimer lifetimes and strong evidence that transient receptor coconfinement promotes repeated interactions between erbB1 monomers. Here, spatial stochastic simulations explore the potential impact of these parameters on erbB1 phosphorylation kinetics. This rule-based mathematical model incorporates structural evidence for conformational flux of the erbB1 extracellular domains, as well as asymmetrical orientation of erbB1 cytoplasmic kinase domains during dimerization. The asymmetric dimer model considers the theoretical consequences of restricted transactivation of erbB1 receptors within a dimer, where the N-lobe of one monomer docks with the C-lobe of the second monomer and triggers its catalytic activity. The dynamic nature of the erbB1 phosphorylation state is shown by monitoring activation states of individual monomers as they diffuse, bind, and rebind after ligand addition. The model reveals the complex interplay between interacting liganded and nonliganded species and the influence of their distribution and abundance within features of the membrane landscape.

Quantitative systems pharmacology model of the amyloid pathway in Alzheimer's disease: Insights into the therapeutic mechanisms of clinical candidates.

Despite considerable investment into potential therapeutic approaches for Alzheimer's disease (AD), currently approved treatment options are limited. Predictive modeling using quantitative systems pharmacology (QSP) can be used to guide the design of clinical trials in AD. This study developed a QSP model representing amyloid beta (Aß) pathophysiology in AD. The model included mechanisms of Aß monomer production and aggregation to form insoluble fibrils and plaques; the transport of soluble species between the compartments of brain, cerebrospinal fluid (CSF), and plasma; and the pharmacokinetics, transport, and binding of monoclonal antibodies to targets in the three compartments. Ordinary differential equations were used to describe these processes quantitatively. The model components were calibrated to data from the literature and internal studies, including quantitative data supporting the underlying AD biology and clinical data from clinical trials for anti-Aß monoclonal antibodies (mAbs) aducanumab, crenezumab, gantenerumab, and solanezumab. The model was developed for an apolipoprotein E (APOE) ɛ4 allele carrier and tested for an APOE ɛ4 noncarrier. Results indicate that the model is consistent with data on clinical Aß accumulation in untreated individuals and those treated with monoclonal antibodies, capturing increases in Aß load accurately. This model may be used to investigate additional AD mechanisms and their impact on biomarkers, as well as predict Aß load at different dose levels for mAbs with known targets and binding affinities. This model may facilitate the design of scientifically enriched and efficient clinical trials by enabling a priori prediction of biomarker dynamics in the brain and CSF.

Minimal Physiologically-Based Pharmacokinetic (mPBPK) Metamodeling of Target Engagement in Skin Informs Anti-IL17A Drug Development in Psoriasis.

The pharmacologic effect(s) of biotherapeutics directed against soluble targets are driven by the magnitude and duration of free target suppression at the tissue site(s) of action. Interleukin (IL)-17A is an inflammatory cytokine that plays a key role in the pathogenesis of psoriasis. In this work, clinical trial data from two monoclonal antibodies (mAbs) targeting IL-17A for treatment of psoriasis (secukinumab and ixekizumab) were analyzed simultaneously to quantitatively predict their target engagement (TE) profiles in psoriatic skin. First, a model-based meta-analysis (MBMA) for clinical responses was conducted separately for each drug based on dose. Next, a minimal physiologically-based pharmacokinetic (mPBPK) model was built to assess skin site IL-17A target engagement for ixekizumab and secukinumab simultaneously. The mPBPK model captured the observed drug PK, serum total IL-17A, and skin drug concentration-time profiles reasonably well across the different dosage regimens investigated. The developed mPBPK model was then used to predict the average TE (i.e., free IL-17A suppression) in skin achieved over a 12-weeks treatment period for each drug following their respective regimens and subsequently assess the TE-efficacy response relationship. It was predicted that secukinumab achieved 98.6% average TE in the skin at 300 mg q4w SC while ixekizumab achieved 99.9% average TE under 160 mg (loading) followed by 80 mg q2w SC. While direct quantification of free IL-17A levels at the site of action is technically challenging, integrated mPBPK-MBMA approaches offer quantitative predictions of free IL-17A levels at the site of action to facilitate future drug development via IL-17A suppression in psoriasis.

Higher naloxone dosing in a quantitative systems pharmacology model that predicts naloxone-fentanyl competition at the opioid mu receptor level.

Rapid resuscitation of an opioid overdose with naloxone, an opioid antagonist, is critical. We developed an opioid receptor quantitative systems pharmacology (QSP) model for evaluation of naloxone dosing. In this model we examined three opioid exposure levels that have been reported in the literature (25 ng/ml, 50 ng/ml, and 75 ng/ml of fentanyl). The model predicted naloxone-fentanyl interaction at the mu opioid receptor over a range of three naloxone doses. For a 2 mg intramuscular (IM) dose of naloxone at lower fentanyl exposure levels (25 ng/ml and 50 ng/ml), the time to decreasing mu receptor occupancy by fentanyl to 50% was 3 and 10 minutes, respectively. However, at a higher fentanyl exposure level (75 ng/ml), a dose of 2 mg IM of the naloxone failed to reduce mu receptor occupancy by fentanyl to 50%. In contrast, naloxone doses of 5 mg and 10 mg IM reduced mu receptor occupancy by fentanyl to 50% in 5.5 and 4 minutes respectively. These results suggest that the current doses of naloxone (2 mg IM or 4 mg intranasal (IN)) may be inadequate for rapid reversal of toxicity due to fentanyl exposure and that increasing the dose of naloxone is likely to improve outcomes.

Orchestration of ErbB3 signaling through heterointeractions and homointeractions.

Members of the ErbB family of receptor tyrosine kinases are capable of both homointeractions and heterointeractions. Because each receptor has a unique set of binding sites for downstream signaling partners and differential catalytic activity, subtle shifts in their combinatorial interplay may have a large effect on signaling outcomes. The overexpression and mutation of ErbB family members are common in numerous human cancers and shift the balance of activation within the signaling network. Here we report the development of a spatial stochastic model that addresses the dynamics of ErbB3 homodimerization and heterodimerization with ErbB2. The model is based on experimental measures for diffusion, dimer off-rates, kinase activity, and dephosphorylation. We also report computational analysis of ErbB3 mutations, generating the prediction that activating mutations in the intracellular and extracellular domains may be subdivided into classes with distinct underlying mechanisms. We show experimental evidence for an ErbB3 gain-of-function point mutation located in the C-lobe asymmetric dimerization interface, which shows enhanced phosphorylation at low ligand dose associated with increased kinase activity.

Analytical solution of steady-state equations for chemical reaction networks with bilinear rate laws.

True steady states are a rare occurrence in living organisms, yet their knowledge is essential for quasi-steady-state approximations, multistability analysis, and other important tools in the investigation of chemical reaction networks (CRN) used to describe molecular processes on the cellular level. Here, we present an approach that can provide closed form steady-state solutions to complex systems, resulting from CRN with binary reactions and mass-action rate laws. We map the nonlinear algebraic problem of finding steady states onto a linear problem in a higher-dimensional space. We show that the linearized version of the steady-state equations obeys the linear conservation laws of the original CRN. We identify two classes of problems for which complete, minimally parameterized solutions may be obtained using only the machinery of linear systems and a judicious choice of the variables used as free parameters. We exemplify our method, providing explicit formulae, on CRN describing signal initiation of two important types of RTK receptor-ligand systems, VEGF and EGF-ErbB1.