Pharmacodynamic modelling of reversible gastric acid pump inhibition in dog and man

Pharmacodynamic model of slow reversible binding and its applications in pharmacokinetic/pharmacodynamic modeling: review and tutorial

Therapeutic responses of most drugs are initiated by the rate and degree of binding to their receptors or targets. The law of mass action describes the rate of drug-receptor complex association (k_on) and dissociation (k_off) where the ratio k_off/k_on is the equilibrium dissociation constant (K_d). Drugs with slow reversible binding (SRB) often demonstrate delayed onset and prolonged pharmacodynamic effects. This report reviews evidence for drugs with SRB features, describes previous pharmacokinetic/pharmacodynamic (PK/PD) modeling efforts of several such drugs, provides a tutorial on the mathematics and properties of SRB models, demonstrates applications of SRB models to additional compounds, and compares PK/PD fittings of SRB with other mechanistic models. We identified and summarized 52 drugs with in vitro-confirmed SRB from a PubMed literature search. Simulations with a SRB model and observed PK/PD profiles showed delayed and prolonged responses and that increasing doses/k_on or decreasing k_off led to greater expected maximum effects and a longer duration of effects. Recession slopes for return of responses to baseline after single doses were nearly linear with an inflection point that approaches a limiting value at larger doses. The SRB model newly captured literature data for the antihypertensive effects of candesartan and antiallergic effects of noberastine. Their PD profiles could also be fitted with indirect response and biophase models with minimal differences. The applicability of SRB models is probably commonplace, but underappreciated, owing to the need for in vitro confirmation of binding kinetics and the similarity of PK/PD profiles to models with other mechanistic determinants.

PREDICTING EFFICACY OF PROTON PUMP INHIBITORS IN REGULATING GASTRIC ACID SECRETION

Developing drugs to treat gastric acid related illnesses such as ulcers and acid reflux disease is the leading focus of pharmaceutical companies. In fact, expenditure for treating these disorders is highest among all illnesses in the US. Over the last few decades, a class of drugs known as a proton pump inhibitors (PPIs) appeared on the market and are highly effective at abating gastric illnesses by raising stomach pH (reducing gastric acid levels). While much is known about the action of PPIs , there are still open questions regarding their efficacy, dosing and long-term effects. Here we extend a previous gastric acid secretion model developed by our group to incorporate a pharmacodynamic/pharmacokinetic model to study proton pump inhibitor (PPI) action. Model-relevant parameters for specific drugs such as omeprazole (OPZ) , lansoprazole (LPZ) and pantoprazole (PPZ) were used from published data, and we conducted simulations to study various aspects of PPI treatment. Clinical data suggests that duration of acid suppression is dependent on proton pump turnover rates and this is supported by our model. We found the order of efficacy of the different PPIs to be OPZ>PPZ>LPZ for clinically recommended dose values, and OPZ>PPZ=LPZ for equal doses. Our results indicate that a breakfast dose for once-daily dosing regimens and a breakfast-lunch dose for twice-daily dosing regimens is recommended. Simulation of other gastric disorders using our model provides atypical applications for the study of drug treatment on homeostatic systems and identification of potential side-effects.

Long-lasting target binding and rebinding as mechanisms to prolong in vivo drug action

An increasing number of examples in the literature suggest that the in vivo duration of drug action not only depends on macroscopic pharmacokinetic properties like plasma half-life and the time needed to equilibrate between the plasma and the effect compartments, but is also influenced by long-lasting target binding and rebinding. The present review combines information from different research areas and simulations to explore the nature of these mechanisms and the conditions in which they are most prevalent. Simulations reveal that these latter phenomena become especially influential when there is no longer sufficient free drug around to maintain high levels of receptor occupancy. There is not always a direct link between slow dissociation and long-lasting in vivo target protection, as the rate of free drug elimination from the effect compartment is also a key influencing factor. Local phenomena that hinder the diffusion of free drug molecules away from their target may allow them to consecutively bind to the same target and/or targets nearby (denoted as 'rebinding') even when their concentration in the bulk phase has already dropped to insignificant levels. The micro-anatomic properties of many effect compartments are likely to intensify this phenomenon. By mimicking the complexity of tissues, intact cells offer the opportunity to investigate both mechanisms under the same, physiologically relevant conditions.

Applicability of the dopamine and rate hypotheses in explaining the differences in behavioral pharmacology of the chloro-benztropine analogs: studies conducted using intracerebral microdialysis and population pharmacodynamic modeling

Previous studies indicated that the chloro-benztropine analogs differed significantly in their cocaine-like activity, which was not expected based on the similarity in their in vitro binding affinity and functional potency at the dopamine transporter (DAT). The present study was designed to extend the understanding of the involvement of both pharmacokinetic and pharmacodynamic factors in mediating the behavioral differences among these analogs. The pharmacokinetics of 3'-chloro-3alpha-(diphenylmethoxy)tropane (3'-Cl BZT), the analog showing a cocaine-like behavioral profile in rodents, was compared with previously reported pharmacokinetic characteristics of cocaine and 4',4''-dichloro-3alpha-(diphenylmethoxy)tropane (4',4''-diCl BZT), an analog totally devoid of cocaine-like actions. Microdialysis studies in rats were conducted to determine whether 3'-Cl and 4',4''-diCl BZT differed significantly in their effect on nucleus accumbens extracellular dopamine levels, with cocaine serving as a reference. A mechanistic model based on DAT association/dissociation kinetics was used to describe the time delay between the plasma concentrations of the chloro-analogs and their dopaminergic effects. 3'-Cl BZT had plasma elimination half-life of 1.9 h versus 0.5 and 21.1 h for cocaine and 4',4''-diCl BZT, respectively. 4',4''-diCl BZT increased the DA levels at a slower rate and to a significantly lower extent relative to 3'-Cl BZT that were, in turn, lower than cocaine. The duration of dopamine elevation was as follows: 4',4''-diCl BZT > 3'-Cl BZT > cocaine. The model indicated faster association and dissociation with DAT for 3'-Cl BZT relative to 4',4''-diCl BZT. The present results indicate that behavioral differences among the chloro-analogs may be explainable based on both the dopamine and rate hypotheses of drug abuse.

Application of a combined effect compartment and binding model for gastric acid inhibition of AR-HO47108: a potassium competitive acid blocker, and its active metabolite AR-HO47116 in the dog

The effect of AR-HO47108, a potassium competitive acid blocker, and its active metabolite AR-HO47116 was studied in Heidenhain pouch dogs following administration of single oral and intravenous doses of the two compounds. The histamine-stimulated acid secretion was measured in different periods after dose up to 24h. All data obtained in the different studies was pooled and analyzed by non-linear mixed effects modelling. It was found that there is a delay between the plasma concentration-time peak and the maximum inhibitory effect and that the effect persisted longer than anticipated from the plasma concentration half-lives of the compounds. In addition, it was found that the peak effect was reached earlier at higher doses. The effect data was well described by a combined effect compartment and binding model and both distribution to the biophase i.e. the canaliculus of the parietal cell and a rate limiting binding interaction between drug and enzyme appear to contribute to the observed delay. In addition, a secretion rate dependent washout from the biophase may contribute. Furthermore, because the parent compound and metabolite bind to the same enzyme, the effect is determined by competition between the two for the same enzyme. The metabolite was found to be less potent than the parent compound, with Kd values of the combined model of 125 and 11.2 nM for the metabolite and parent compound, respectively. However, the metabolite is generated in high concentrations that rapidly exceed the concentration of parent compound after oral administration of parent compound, and this together with its longer plasma half-life will make its contribution to the overall effect increase with time and slightly prolong the duration of the effect.

Assessment of the relative in vivo potency of the hydroxylated metabolite of darifenacin in its ability to decrease salivary flow using pooled population pharmacokinetic-pharmacodynamic data

AIMS

To describe the population pharmacokinetic-pharmacodynamic relationship between darifenacin (UK-88,525) and its hydroxylated metabolite (UK-148,993), and the reduction in salivary flow (SF, a M3-mediated response). This enabled an estimation of the in vivo potency of the metabolite to decrease SF relative to that of the parent drug.

METHODS

A total of 262 individuals were pooled from 11 Phase 1 studies and one Phase 2 study. A comparison was made between a series of pharmacodynamic models (direct-effect, indirect-effect, link and binding model) using NONMEM.

RESULTS

The binding model yielded the best description of the decrease in SF by fully accounting for the time course of the pharmacodynamic effect. An internal validation exercise demonstrated the robustness of this model. Covariate analysis identified a circadian rhythm in SF. This model, with confidence intervals (CI) determined by likelihood profiling, indicated that the relative potency of the metabolite to darifenacin to reduce SF was 11.1% (95% CI 3.8, 19.6). This implied that the metabolite was ninefold less potent than darifenacin in vivo. Accounting for the unbound fraction of darifenacin (2%) and its metabolite (13%), the in vivo protein binding-corrected relative potency was estimated to be 2.1%, indicating that the metabolite was 50-fold less potent than the parent drug. The model supported the assumption that no other metabolites contributing to the impairment of the SF were formed during first-pass, and that the development of sensitization or tolerance was not evident over time. The validation process indicated that the i.v.-oral crossover study was necessary for the estimation of the relative potency.

CONCLUSIONS

Population modelling of darifenacin and its hydroxylated metabolite yielded individual pharmacokinetic predictions that could be used to assess the in vivo potency of the metabolite to decrease SF relative to that of the parent drug. The metabolite had a negligible effect on SF.