Human physiologically based pharmacokinetic model for ACE inhibitors: ramipril and ramiprilat

BACKGROUND: The angiotensin-converting enzyme (ACE) inhibitors have complicated and poorly characterized pharmacokinetics. There are two binding sites per ACE (high affinity "C", lower affinity "N") that have sub-nanomolar affinities and dissociation rates of hours. Most inhibitors are given orally in a prodrug form that is systemically converted to the active form. This paper describes the first human physiologically based pharmacokinetic (PBPK) model of this drug class. METHODS: The model was applied to the experimental data of van Griensven et. al for the pharmacokinetics of ramiprilat and its prodrug ramipril. It describes the time course of the inhibition of the N and C ACE sites in plasma and the different tissues. The model includes: 1) two independent ACE binding sites; 2) non-equilibrium time dependent binding; 3) liver and kidney ramipril intracellular uptake, conversion to ramiprilat and extrusion from the cell; 4) intestinal ramipril absorption. The experimental in vitro ramiprilat/ACE binding kinetics at 4°C and 300 mM NaCl were assumed for most of the PBPK calculations. The model was incorporated into the freely distributed PBPK program PKQuest. RESULTS: The PBPK model provides an accurate description of the individual variation of the plasma ramipril and ramiprilat and the ramiprilat renal clearance following IV ramiprilat and IV and oral ramipril. Summary of model features: Less than 2% of total body ACE is in plasma; 35% of the oral dose is absorbed; 75% of the ramipril metabolism is hepatic and 25% of this is converted to systemic ramiprilat; 100% of renal ramipril metabolism is converted to systemic ramiprilat. The inhibition was long lasting, with 80% of the C site and 33% of the N site inhibited 24 hours following a 2.5 mg oral ramipril dose. The plasma ACE inhibition determined by the standard assay is significantly less than the true in vivo inhibition because of assay dilution. CONCLUSION: If the in vitro plasma binding kinetics of the ACE inhibitor for the two binding sites are known, a unique PBPK model description of the Griensven et. al. experimental data can be obtained

Estimating impossible curves using NONMEM

1On fitting model equations to experimental data, the situation may arise that individual subjects provide insufficient information to obtain adequate parameter estimates due to the fact that not all aspects are exhibited by all subjects or that the models are simply too complex. This may be solved by applying nonlinear mixed effect modelling to the data, which integrates the information provided by different subjects

Relationship between airway responsiveness to neurokinin A and methacholine in asthma

Non-adrenergic non-cholinergic (NANC) nerves release bronchoactive tachykinins such as substance P (SP) and neurokinin A (NKA) that can induce features of asthma. The airway response to NKA in humans closely resembles that of methacholine (M). Hence, we investigated the relationship between airway responsiveness to NKA and M in subjects with asthma. To this end, we analyzed baseline data of 27 subjects with mild persistent asthma (20F/7M) 19-46 y; FEV1 81-136% pred.; PC20FEV1 (M) or = 20% fall from baseline (PC20FEV1). Twenty-two subjects reached a PC20FEV1 on both occasions. The PC20FEV1 values of both agonists correlated significantly (Spearman's r=-0.721; p=0.0002), and the relationship was given by 10log(PC20FEV1(NKA))= -1.36 + (0.60 x 10log(PC20FEV1(M)). We have demonstrated a significant relationship between airway responsiveness to NKA and methacholine in asthma. This suggests that both agonists may share common final pathways in causing bronchoconstriction in patients with mild persistent asthma. Based on our data and previous studies in asthma, it can be hypothesized that this direct NKA-induced bronchoconstrictor response may be mediated by predominant stimulation of the tachykinin NK-2 receptors on airway smooth muscle cell

Estimation of growth hormone secretion rate: Impact of kinetic assumptions intrinsic to the analytical approach

We compared four common mathematical techniques to determine daily endogenous growth hormone (GH) secretion rates from diurnal plasma GH concentration profiles in 24 women (16 upper- or lower-body obese and 8 normal-weight individuals). Two forms of deconvolution analysis and two techniques based on a priori determined GH clearance estimates were employed. Deconvolution analyses revealed significant differences in the 24-h GH secretion rate between normal-weight and upper-body obese women, whereas the other two techniques did not. Moreover, deconvolution analyses predicted that the reduction in mean plasma GH concentrations in upper-body obese women was accounted for by impaired GH secretion, whereas the other methods suggested that obesity increases GH metabolic clearance. Thus we infer that disparate conclusions concerning GH secretion can be drawn from the same primary data set. The different inferences likely reflect dissimilar kinetic assumptions and the particular limitations intrinsic to each analytical approach. Accordingly, we urge caution in the facile comparison of calculated GH secretion data in humans, especially when kinetic and secretion measurements are performed under different conditions. The most appropriate way to determine the GH secretion rate in humans must be balanced by the exact intent of the experiment and the acceptability of different assumptions in that context

Circadian rhythm of plasma leptin levels in upper and lower body obese women: Influence of body fat distribution and weight loss

Plasma leptin concentrations were measured every 20 min for 24 h in eight normal weight women and in eight upper body and eight lower body obese women matched for body mass index. The circadian rhythm of leptin, which could mathematically be described by a cosine, was characterised by an acrophase just after midnight in all subjects. The amplitude of a cosine fit as well as the average 24-h leptin concentration were increased by 280% and 420%, respectively, in obese compared to normal weight women. All characteristics of leptin concentration profiles were similar in upper body and lower body obese women, except for a significantly higher amplitude in the lower body obese group. Visceral and sc body fat depots were measured using magnetic resonance imaging in all three groups. Average 24-h leptin concentrations were strongly correlated with sc fat (r = 0.84), whereas visceral fat was not an independent predictor of the plasma leptin level. A loss of 50% of the overweight was associated with a 55% decrease in the average 24-h leptin concentrations in obese women (95% confidence interval, 12.3, 26.6), whereas the characteristics of the circadian rhythm of leptin remained unchanged. Finally, it was observed that a fasting plasma leptin concentration is not an acceptable indicator of the average leptin concentration over 24 h

Influence of obesity and body fat distribution on growth hormone kinetics in humans

We studied the kinetics of exogenous recombinant 22-kDa human growth hormone (rhGH) in premenopausal women with upper body obesity (UBO), lower body obesity (LBO), or normal body weight. A bolus of 100 mU rhGH was administered during a continuous infusion of somatostatin to suppress endogenous GH secretion. GH kinetics were investigated with noncompartmental analysis of plasma GH curves. GH peak values in response to GH infusion and plasma half-life of GH were not significantly different between normal weight and obese subjects. In contrast, GH clearance was 33% higher in LBO women and 51% higher in UBO women compared with clearance in normal weight controls. The difference in clearance between LBO and UBO was not statistically significant. Altered GH clearance characteristics contribute to low circulating GH levels in obese humans. Body fat distribution does not appear to affect GH kinetics