Mathematical modeling of cortisol circadian rhythm and cortisol suppression

Depression diagnostics using a nonlinear mathematical oscillatory model

BACKGROUND AND OBJECTIVES

It is known that long-term stress leads to trauma and very often to depression. Usually, the diagnosis of depression is dealt with by psychiatrists who, based on conversations and questions, diagnose the patient's illness and condition. Unfortunately, this diagnosis is not always reliable. To prevent the development of disease, it is necessary to detect illness in a timely manner. One of the indications of the possibility of the onset of disease is a disturbance in the level of hormones in the body, especially cortisol. The purpose of this study was to develop a mathematical model for cortisol variation resulting from stress which would be useful in making conclusions about depressive states.

METHODS

Rapid changes in cortisol concentration, according to ultradian rhythms, which are much faster than the daily circadian rhythm, is modelled as a truly nonlinear oscillator. The mathematical model contains two coupled first order differential equations. The stress is modeled as a pulsating action, described with a periodic trigonometric function, and cortisol production as a cubic nonlinear one. Three models for cortisol variation are considered: 1) the pure nonlinear model, 2) the periodically excited system, 3) and the chaotic system. The results from the study are supported with experimental measurements.

RESULTS

Without stress, cortisol variation is of an oscillatory type with a constant steady-state amplitude. Intensive stress causes a resonant phenomenon in cortisol oscillatory variation. The occasion is short and is usually without consequences. For long stress periods deterministic chaos occurs which permanently changes the levels of cortisol. This phenomenon is an indicator of depression. Results from the suggested models are compared with experimentally obtained ones and good quantitative agreement is obtained.

CONCLUSIONS

The nonlinear oscillator is a good model for indication of depression. The model provides not only general conclusions, but also individual ones, if personal characteristics are taken into consideration. Response of the model depends not only on the input data related to stress, but also on the system parameters that specify each individual. Findings obtained from this study have implications for the medical diagnosis and treatment of depression.

A quantitative approach to analysing cortisol response in the horse

The cortisol response to glucocorticoid intervention has, in spite of several studies in horses, not been fully characterized with regard to the determinants of onset, intensity and duration of response. Therefore, dexamethasone and cortisol response data were collected in a study applying a constant rate infusion regimen of dexamethasone (0.17, 1.7 and 17 μg/kg) to six Standardbreds. Plasma was analysed for dexamethasone and cortisol concentrations using UHPLC-MS/MS. Dexamethasone displayed linear kinetics within the concentration range studied. A turnover model of oscillatory behaviour accurately mimicked cortisol data. The mean baseline concentration range was 34-57 μg/L, the fractional turnover rate 0.47-1.5 1/h, the amplitude parameter 6.8-24 μg/L, the maximum inhibitory capacity 0.77-0.97, the drug potency 6-65 ng/L and the sigmoidicity factor 0.7-30. This analysis provided a better understanding of the time course of the cortisol response in horses. This includes baseline variability within and between horses and determinants of the equilibrium concentration-response relationship. The analysis also challenged a protocol for a dexamethasone suppression test design and indicated future improvement to increase the predictability of the test.

Quantitative characterization of circadian rhythm of pulmonary function in asthmatic patients treated with inhaled corticosteroids

The aim of this study was to characterize the circadian rhythm observed for forced expiratory volume in 1 s (FEV1) in patients with persistent asthma being treated with inhaled corticosteroids. The database included 3379 FEV1 measurements from 189 patients with mild to moderate asthma. A model using the sum of two Sine functions with periods of 12 and 24 h and a constant component of mean circadian rhythm adequately described the circadian rhythm in FEV1 measurements over time. The model adequateness was evaluated by various approaches including visual predictive check (VPC), prediction-corrected VPC, standardized VPC and normalized prediction distribution error. Covariates tested included age, body weight, height, body mass index, baseline FEV1, and sex. Age and height were found to have significant effects on the mean FEV1 level and no covariate was found to have an effect on the magnitude and timing of circadian rhythm. The model predicted that a minimum FEV1 occurred in the early morning and maximum FEV1 occurred in the early afternoon, with a population mean fluctuation of 170 mL, which is consistent with the finding that asthma symptoms usually exacerbate in the early morning for patients with persistent asthma. This developed model provides the first quantitative approach to describing FEV1 circadian rhythm with ICS background treatment and provided insight in designing future registration trials for asthma drug development.

Presence of endogenous prednisolone in human urine

The possibility of an endogenous presence of the glucocorticoid prednisolone has already been demonstrated in bovine and horse urine, with the aim of clarifying its origin in this matrix, which is used by official agencies for the control of illicit treatments. From this point of view, the endogenous nature of prednisolone could be a major topic in doping control of both amateur and professional human athletes. A study was therefore made on 34 human volunteers (13 males and 21 females; aged 22-62) to detect the presence of prednisolone in their urine by HPLC-MS(3). One of the volunteers underwent vernal allergy treatment with betamethasone for two subsequent years. An investigation was carried out with the aim of verifying if the suppression, and the circadian rhythm, of cortisol urinary levels could also apply to prednisolone. The results of the study show that prednisolone was present in the urine of all 34 volunteers, with a concentration very close to 100-times lower that of cortisol, with no dependence on gender. The same ratio (1/100) was observed in the prednisolone and cortisol levels detected during the 24h together with the suppression of prednisolone by betamethasone treatment. These data demonstrate the endogenous nature of low concentrations of prednisolone in human urine, and motivate further studies about the biosynthetic pathways of this corticosteroid and its relationship with stress in humans, as already described in cows.

A pharmacokinetic simulation tool for inhaled corticosteroids

The pharmacokinetic (PK) behavior of inhaled drugs is more complicated than that of other forms of administration. In particular, the effects of certain physiological (mucociliary clearance and differences in membrane properties in central and peripheral (C/P) areas of the lung), formulation (as it relates to drug deposition and particle dissolution rate), and patient-related factors (lung function; effects on C/P deposition ratio) affect the systemic PKs of inhaled drugs. The objectives of this project were (1) to describe a compartmental model that adequately describes the fate of inhaled corticosteroids (ICS) after administration while incorporating variability between and within subjects and (2) based upon the model, to provide a freely available tool for simulation of PK trials after ICS administration. This compartment model allows for mucociliary removal of undissolved particles from the lung, distinguishes between central and peripheral regions of the lung, and models drug entering the systemic circulation via the lung and the gastrointestinal tract. The PK simulation tool is provided as an extension package to the statistical software R ('ICSpkTS'). It allows simulation of PK trials for hypothetical ICS and of four commercially available ICS (budesonide, flunisolide, fluticasone propionate, and triamcinolone acetonide) in a parallel study design. Simulated PK data and parameters agreed well with literature data for all four ICS. The ICSpkTS package is especially suitable to explore the effect of changes in model parameters on PK behavior and can be easily adjusted for other inhaled drugs.

From pharmacokinetics to therapeutics

Whilst pharmacokinetics describe the relationship between dose levels and concentration-time profiles of a drug in the body and pharmacodynamics describe the concentration-response relationships, pharmacokinectics-pharmacodynamics(PK-PD) models link these two items providing a framework for modelling the time course of drug response. In this chapter, PK-PD models, describing the therapeutic effects of drugs used for the therapy of allergic diseases have been reviewed. Emphasis was given also to the description of the receptor occupancy, which is tightly related to the downstream clinical response. PK - PD models describing unwanted effects were also commented. An integrated use of these models allows choosing appropriate dosing regimens and providing an objective evaluation of the benefit/risk balance.

Cortisol suppression as a surrogate marker for inhaled corticosteroid-induced growth retardation in children

Exposure of inhaled corticosteroids (ICSs) in pediatrics results in adrenal suppression and growth inhibition. The objective of this study was to assess the relationship of ICS mediated growth retardation with cortisol suppression in asthmatic children. A meta-analysis approach was performed with 33 published articles. Growth velocity (GV) data were obtained from the literature for evaluation of growth. Cumulative cortisol suppression within 24h (CCS%) was calculated at steady state with a validated algorithm. Consolidated GV and CCS% data were employed for model development. A linear mixed effects model was developed to adequately describe the relationship between GV and CCS%. No impact of tested covariates was observed. Population estimate of the rate of change in GV was -0.06cm/year/CCS% (12.7%, coefficient of variation) for both stadiometry and knemometry methods. However, GV from stadiometry is expected to be approximately three fold higher than that from knemometry when cortisol suppression was not presented. The final model was evaluated with posterior predictive check and pattern check approaches. The results from this study elucidate CCS% as an excellent predictor of ICS mediated growth retardation in asthmatic children.

Assessment of the impact of dosing time on the pharmacokinetics/pharmacodynamics of prednisolone

Prednisolone is widely used for the treatment of inflammation and auto-immune diseases. It exhibits nonlinear pharmacokinetics (PK); and its induced systemic effects (pharmacodynamics (PD)) are commonly evaluated with two biomarkers, cortisol and blood lymphocytes in plasma. Circadian patterns are observed in both biomarkers. Furthermore, the disease itself may show a circadian pattern. For example, in rheumatoid arthritis patients, better therapeutic outcomes have been reported when prednisolone was administered in the very early morning. The aim of this study is to evaluate the impact of dosing time on the PK/PD of prednisolone with a simulation approach using an interactive algorithm. A series of simulations were performed with either intravenous or oral administration of prednisolone or prednisone. The results showed that the initial or maximum concentration and trough concentration of total prednisolone were lower when the drug was administered in the early morning around 6 AM: . Oscillation patterns were observed in cumulative cortisol suppression (CCS) and alteration of total lymphocyte trafficking in blood. When the drug was given in the morning within the therapeutic dose range, or around 6 PM: for a small dose amount (<1 mg), the minimum CCS and maximum effect on lymphocytes were observed. These results indicated that the PK/PD of prednisolone are time- and dose-dependent, and suggested that it is necessary to consider the application of chronotherapy to achieve better clinical outcomes with fewer side effects of prednisolone, and a PK/PD simulation approach could provide a valuable tool to evaluate and predict time-dependency in the system.

A Pharmacokinetic/Pharmacodynamic Approach to Predict Total Prednisolone Concentrations in Human Plasma

Prednisolone and prednisone are two widely used corticosteroids for various inflammatory and immune diseases. Prednisolone is the active form of prednisone in vivo. Total prednisolone in plasma exhibits nonlinear pharmacokinetics mainly due to its nonlinear protein binding. Other factors such as reversible metabolism (or interconversion between prednisolone and prednisone), competitive protein binding from endogenous cortisol, cortisol circadian rhythm, and prednisolone mediated cortisol suppression complicate prednisolone pharmacokinetics. This study was aimed to develop a new approach to describe the nonlinear pharmacokinetics of total prednisolone and predict total prednisolone concentrations in plasma. Based on literature datasets, a linear two-compartment pharmacokinetic model was developed to adequately describe the reversible metabolism between free prednisone and prednisolone. Cortisol and prednisolone protein binding were described via the sum of a Langmuir and linear type binding. The endogenous cortisol circadian rhythm and cortisol suppression during prednisone or prednisolone exposure were described with a previously reported linear release rate pharmacokinetic/pharmacodynamic (PK/PD) model. By combining the pharmacokinetic models for free prednisone and prednisolone, the linear release rate model for cortisol suppression, and competitive protein binding between cortisol and prednisolone, we were able to predict total prednisolone concentrations in plasma. The predicted total prednisolone concentrations in plasma were in good agreement with the literature reported data. Thus, this PK/PD approach shows that the combination of nonlinear protein binding, cortisol circadian rhythm, and cortisol suppression could account for the nonlinearity of total prednisolone. In addition, it also allows a valid prediction of total prednisolone in plasma after either prednisone or prednisolone administration.

Quantitative structure-pharmacokinetic/pharmacodynamic relationships

Quantitative structure-activity relationships have long been considered a vital component of drug discovery and development, providing insight into the role of molecular properties in the biological activity of similar and unrelated compounds. Recognition that in vitro bioassay and/or pre-clinical activity are insufficient for anticipating which compounds are suitable leads for further development has shifted the focus toward integrated pharmacokinetic (PK) and pharmacodynamic (PD) processes. Over the last decade, considerable progress has been made in constructing empirical and mechanistic quantitative structure-PK relationships (QSPKR), as well as diverse mechanism-based pharmacodynamic models of drug effects. In this review, traditional and contemporary approaches to developing QSPKR models are discussed, along with selected examples of attempts to couple QSPKR and pharmacodynamic models to anticipate the intensity and time-course of the pharmacological effects of new or related compounds, or quantitative structure-pharmacodynamic relationships modeling. Such models are in accordance with the goals of systems biology and the ideal of designing drugs and delivery systems from first principles.

Model-based covariate pharmacokinetic analysis and lack of cortisol suppression by the new inhaled corticosteroid ciclesonide using a novel cortisol release model

Ciclesonide is a novel and effective inhaled corticosteroid for the treatment of asthma that is converted into the active metabolite C21-desisobutyryl-ciclesonide (des-CIC) in the lung. The objectives of this analysis were to characterize covariate effects on des-CIC pharmacokinetics (PK) and circadian cortisol release. In addition, the effect of systemic des-CIC exposure on circadian cortisol release was also evaluated. Data from 12 phase I and 3 phase III studies in adults and 2 phase III studies in children (ex-actuator dose of ciclesonide: 40-2880 microg) were pooled. There were 635 subjects in the analysis with 5238 and 4470 observations recorded for concentrations of des-CIC and cortisol, respectively. The des-CIC PK was described using a 1-compartment model with first-order absorption and no clinically relevant differences in des-CIC PK due to body weight, race, asthma severity, or gender. The mean PK parameters estimates of apparent clearance and apparent volume of distribution were 302 L/h and 1310 L, respectively. A 1-compartment model with first-order absorption, an endogenous "predose" cortisol concentration at dose interval and a lag time based on a fixed, hypothetical cortisol dosing time of 10:00 PM could adequately characterize the circadian rhythm of endogenous cortisol release. The potential effect of des-CIC on the circadian rhythm of endogenous cortisol release was evaluated using dose and AUC as covariates and applying the Emax model. Less than 1% of all observed des-CIC concentrations were higher than the EC50 for cortisol suppression, indicating negligible changes in cortisol concentrations at therapeutically relevant doses.

Pharmacokinetics and pharmacodynamics of systemically administered glucocorticoids

Glucocorticoids have pleiotropic effects that are used to treat diverse diseases such as asthma, rheumatoid arthritis, systemic lupus erythematosus and acute kidney transplant rejection. The most commonly used systemic glucocorticoids are hydrocortisone, prednisolone, methylprednisolone and dexamethasone. These glucocorticoids have good oral bioavailability and are eliminated mainly by hepatic metabolism and renal excretion of the metabolites. Plasma concentrations follow a biexponential pattern. Two-compartment models are used after intravenous administration, but one-compartment models are sufficient after oral administration.The effects of glucocorticoids are mediated by genomic and possibly nongenomic mechanisms. Genomic mechanisms include activation of the cytosolic glucocorticoid receptor that leads to activation or repression of protein synthesis, including cytokines, chemokines, inflammatory enzymes and adhesion molecules. Thus, inflammation and immune response mechanisms may be modified. Nongenomic mechanisms might play an additional role in glucocorticoid pulse therapy. Clinical efficacy depends on glucocorticoid pharmacokinetics and pharmacodynamics. Pharmacokinetic parameters such as the elimination half-life, and pharmacodynamic parameters such as the concentration producing the half-maximal effect, determine the duration and intensity of glucocorticoid effects. The special contribution of either of these can be distinguished with pharmacokinetic/pharmacodynamic analysis. We performed simulations with a pharmacokinetic/pharmacodynamic model using T helper cell counts and endogenous cortisol as biomarkers for the effects of methylprednisolone. These simulations suggest that the clinical efficacy of low-dose glucocorticoid regimens might be increased with twice-daily glucocorticoid administration.

Population pharmacokinetics and pharmacodynamics of ciclesonide

Ciclesonide is a novel glucocorticoid that is converted into ciclesonide--active principle (CIC-AP) in the lung. The study objectives were to identify a structural model for population pharmacokinetic (PK) analysis of CIC-AP using nonlinear mixed-effects modeling, assess the influence of select covariates on PK and/or pharmacodynamic (PD) parameters, and investigate the effects of CIC-AP on endogenous cortisol. Pooled concentration data from nine phase I studies (dose: 400-3600 micrograms) involving healthy and asthmatic patients were included in the PK analysis. There were 151 subjects (3300 observations) for the CIC-AP population PK analysis. Various models examined inter- and intrasubject variability for the PK parameters. Population estimates of the PK parameters of clearance and volume of distribution were 396 L/h (64.8% co-efficient of variation [CV]) and 1190 L (41.2% CV), respectively. Pharmacodynamic population estimates included maximum cortisol release rate, 3140 ng/h (5.4% CV). The EC50 of CIC-AP was 0.88 ng/mL. Ciclesonide is a safe corticosteroid that causes negligible cortisol suppression. The disposition and effect of CIC-AP can be described using mixed-effect modeling. The estimated EC50 is similar to mean Cmax from an 800-micrograms dose, further suggesting CIC-AP has little effect on cortisol suppression.

Quantitative structure-pharmacokinetic/pharmacodynamic relationships of corticosteroids in man

The purpose of this study was to develop quantitative structure-activity/pharmacokinetic relationships (QSAR/QSPKR) for 11 selected corticosteroids in man. Multiple linear regression analysis with an automatic forward step-wise inclusion algorithm was used to construct QSAR/QSPKR models from molecular and submolecular descriptors that were generated using the SYBYL and KowWin computer programs. The final equations describing steroid relative receptor affinity, systemic clearance, volume of distribution, fraction unbound in plasma, and percent of oral absorption, all showed significant correlations (R(2) range 0.841 to 0.977). These relationships were crossvalidated using the leave-one-out method, and appeared to have good predictive performance (Q(2) range 0.715 to 0.912). In addition, a general method for integrating QSAR/QSPKR data to predict the time course of pharmacologic effects is presented. This approach, termed quantitative structure-pharmacodynamic relationships modeling, was successfully applied to predict the rapid cortisol suppressive effects of triamcinolone acetonide after a 2-mg intravenous bolus dose in healthy volunteers.

Measurement of cortisol secretion rate by stable isotope methodology following oral administration of 13C-labeled cortisol to a human subject

Plasma concentration measurements of 13C-labeled cortisol ([1,2,4,19-13C(4)]cortisol, cortisol-13C(4)) and its metabolite cortisone-13C(4) were made simultaneously with measurements of endogenous cortisol and cortisone by gas chromatography-mass spectrometry (GC-MS). After administering a small amount (3mg) of cortisol-13C(4) to a human subject, changes in cortisol secretion rates were estimated by deconvolution techniques from the measured plasma cortisol and cortisone levels and the rates of elimination and interconversion of cortisol and cortisone were obtained from the plasma concentration-time data of cortisol-13C(4) and cortisone-13C(4). The objective of this study was to look for a novel approach to quantitate rates of minute-to-minute cortisol secretion in man by taking into account the interconversion of cortisol and cortisone by 11beta-hydroxysteroid dehydrogenase (11beta-HSD).

A new solution-based intranasal triamcinolone acetonide formulation in patients with perennial allergic rhinitis: how does the pharmacokinetic/pharmacodynamic profile for cortisol suppression compare with an aqueous suspension-based formulation?

The present study was undertaken to describe the pharmacokinetics of a new solution-based intranasal triamcinolone acetonideformulation (Tri-Nasal) in patients with perennial allergic rhinitis and to use a pharmacokinetic/pharmacodynamic (PK/PD) simulation approach to compare the potential effects on plasma cortisol with that of an aqueous suspension-based nasal triamcinolone acetonide formulation (Nasacort AQ). Data from an open-label, randomized, three-way crossover study in patients with perennial allergic rhinitis receiving three doses (100, 200, and 400 microg) of a nasal solution-based triamcinolone acetonide formulation (Tri-Nasal) over 7 days were used to describe the pharmacokinetics of this formulation. Available literature data for a suspension-based aqueous triamcinolone acetonide formulation (Nasacort AQ) were used to describe its pharmacokinetic profile after similar single doses of 110, 220, and 440 microg. A PK/PD simulation approach was used to predict the anticipated cumulative cortisol suppression (CCS) of these two formulations. These simulations suggested a cortisol suppression of 8% to 16% for the single and steady-state doses of the solution-based product. Similar CCS estimates were predicted for equivalent doses of the aqueous suspension-based triamcinolone acetonide formulation with no difference between both formulations. Post hoc power analysis suggested that the predicted cortisol suppression is not likely to be significant for either preparation, including the clinically recommended doses of 200 and 220 microg of the solution-based and suspension-based formulations, respectively. In summary, based on the results of this PK/PD simulation, the plasma levels observed afternasal administration of the solution or the aqueous suspension are unlikely to induce a clinically relevant cortisol suppression, especially for the recommended dosing regimens of 200 and 220 microg/day.

An interactive algorithm for the assessment of cumulative cortisol suppression during inhaled corticosteroid therapy

The objective of the study was to develop an algorithm based on a pharmacokinetic-pharmacodynamic (PK/PD) modeling approach to quantify and predict cumulative cortisol suppression (CCS) as a surrogate marker for the systemic activity of inhaled corticosteroid therapy. Two Excel spreadsheets, one for single dose and another for steady-state multiple doses of inhaled steroids, were developed for predicting CCS. Four of the commonly used inhaled steroids were chosen for the purposes of simulation: fluticasone propionate (FP), budesonide (BUD), flunisolide (FLU), and triamcinolone acetonide (TAA). Drug-specific PK and PD parameters were obtained from previous single- and multiple-dose studies. In cases in which multiple-dose data were not available, the single-dose data were extrapolated. The algorithm was designed to calculate CCS based on 5 input parameters: name of drug, dose, dosing interval, time(s) of dosing, and type of inhaler device. In addition, a generalized algorithm was set up to calculate CCS based on clearance, volume of distribution, absorption rate, protein binding, pulmonary deposition, oral bioavailability, and unbound EC50 of the corticosteroid of interest. The spreadsheet allowed predictions of CCS for single doses as well as steady-state conditions. A simple method has been developed that facilitates comparisons between various drugs and dosing regimens and has the potential to significantly reduce the number of comparative clinical trials to be performed for evaluating the short-term systemic activity of inhaled corticosteroids.

Algorithm for application of Fourier analysis for biorhythmic baselines of pharmacodynamic indirect response models

The change of an indirect pharmacological response R(t) can be described by a periodic time-dependent production rate kin(t) and a first-order loss constant kout. If kin(t) follows some biological rhythm (e.g., circadian), then the response R(t) also displays a periodic behavior. A new approach for describing the input function in indirect response models with biorhythmic baselines of physiologic substances is introduced. The present approach uses the baseline (placebo) response Rb(t) to recover the equation for kin(t). Fourier analysis provides an approximate equation for Rb(t) that consists of terms (usually two or three) of the Fourier series (harmonics) that contribute most to the overall sum. The model differential equation is solved backward for kin(t), yielding the equation involving Rb(t). A computer program was developed to perform the square L2-norm approximation technique. Fourier analysis was also performed based on nonlinear regression. Cortisol suppression after inhalation of fluticasone propionate (FP) was modeled based on the inhibition of the secretion rate kin(t) using ADAPT II. The pharmacodynamic parameters kout and IC50 were estimated from the model equation with kin(t) derived by the new approach. The proposed method of describing the input function needs no assumption about the behavior of kin(t), is as efficient as methods used previously, and is more flexible in describing the baseline data than the nonlinear regression method.

Pharmacokinetic and adrenal interactions of IL-10 and prednisone in healthy volunteers

The pharmacokinetic and adrenal interactions of recombinant human interleukin-10 and prednisolone were examined in this open-label, randomized, four-way crossover study in 12 healthy adult male volunteers. Single doses of IL-10 (8 micrograms/kg s.c.), IL-10 with prednisone (15 mg p.o.), placebo with prednisone, or placebo were administered on four separate occasions with at least 3-week interceding washout periods. Measurements included plasma prednisone, prednisolone and cortisol, unbound prednisolone, and serum IL-10 concentrations. Pharmacokinetic parameters were determined using noncompartmental and model-fitting analysis, while area analysis and an indirect response model were used to assess cortisol dynamics. IL-10 exhibited prolonged serum concentrations owing to dual-absorption processes that were largely unaffected by prednisone. The Cmax values were about 3 ng/mL, while the tmax occurred at 7 to 9 hours. Prednisolone exhibited rapid systemic kinetics with a Cmax of 235 ng/mL, tmax at 1.11 hours, and t1/2 of 2.54 hours with no significant alterations owing to IL-10. Both prednisolone and prednisolone/IL-10 caused marked suppression of cortisol concentrations with similar magnitude and IC50 values; however, IL-10 alone significantly increased the 24-hour AUC of cortisol by 20%. Thus, IL-10 and prednisolone do not interact in disposition or adrenal suppression to a clinically significant degree.

Pharmacokinetics and pharmacodynamics of inhaled corticosteroids

There are significant differences in the pharmacokinetic properties of inhaled corticosteroids currently used in medical practice. All are rapidly cleared from the body but they show varying levels of oral bioavailability and more importantly variation in the rate of absorption after inhalation. Oral bioavailability is lowest for fluticasone propionate, indicating a low potential for unwanted systemic corticosteroid effects. Mathematical modeling has shown pulmonary residence times to be longest for fluticasone propionate and triamcinolone acetonide but shortest for budesonide and flunisolide. These properties appear to relate to pulmonary solubility, which appears to be the rate-limiting step in the absorption process.

Characteristics of indirect pharmacodynamic models and applications to clinical drug responses

This review describes four basic physiologic indirect pharmacodynamic response (IDR) models which have been proposed to characterize the pharmacodynamics of drugs that act by indirect mechanisms such as inhibition or stimulation of the production or dissipation of factors controlling the measured effect. The principles underlying IDR models and their response patterns are described. The applicability of these basic IDR models to characterize pharmacodynamic responses of diverse drugs such as inhibition of gastric acid secretion by nizatidine and stimulation of MX protein synthesis by interferon alpha-2a is demonstrated. A list of other uses of these models is provided. These models can be readily extended to accommodate additional complexities such as nonstationary or circadian baselines, equilibration delay, depletion or accumulation of a precursor pool, sigmoidicity, or other mechanisms. Indirect response models which have a logical mechanistic basis account for time-delays in many responses and are widely applicable in clinical pharmacology.

Pharmacokinetic and pharmacodynamic properties of inhaled corticosteroids in relation to efficacy and safety

There are significant differences in the pharmacokinetic properties of inhaled corticosteroids currently available for use in treatment of asthma and this can result in differences in pharmacodynamic activity. All currently used inhaled corticosteroids are rapidly cleared from the body, but show varying levels of oral bioavailability, with fluticasone propionate having the lowest. Following inhalation, there is also considerable variability in the rate of absorption from the lung, and pulmonary residence times are greatest for fluticasone propionate and triamcinolone acetonide, and shortest for budesonide and flunisolide. Cortisol suppression is frequently used as a surrogate marker of systemic corticosteroid activity. Cortisol release displays a circadian rhythm, which can be mathematically modelled and the effects of exogenous corticosteroids on cortisol suppression established. However, when interpreting the effects of inhaled corticosteroids on cumulative cortisol suppression, it is important to take into consideration the pharmacokinetic properties of each particular drug, together with the study design and the time of administration.

Pharmacokinetic/pharmacodynamic evaluation of systemic effects of flunisolide after inhalation

The pharmacokinetics and pharmacodynamics of flunisolide were studied in healthy volunteers after inhalation. In the morning on the day the study began, volunteers inhaled 0.5 mg of flunisolide with and without oral administration of charcoal, or 1 mg, 2 mg, and 3 mg of flunisolide with concomitant administration of charcoal. A placebo group was used to assess the endogenous cortisol, granulocyte, and lymphocyte baseline levels. Flunisolide plasma levels were determined by high-performance liquid chromatography using a tandem mass spectrometer as detector (HPLC/MS/MS). Cortisol plasma levels and differential white blood cell counts were obtained over 12 hours. An integrated pharmacokinetic/pharmacodynamic (PK/PD) model was applied to link the flunisolide plasma concentrations with the effects on lymphocytes, granulocytes, and cortisol. Maximum concentration levels of 3 to 9 ng/mL of flunisolide were observed after 0.2 to 0.3 hours for all of the investigated doses. The terminal half-life ranged from 1.3 to 1.7 hours. There was no statistical difference between treatments in the presence or absence of orally administered charcoal. The pharmacokinetic/pharmacodynamic (PK/PD) models satisfactorily described the time-courses of the effects on granulocytes, lymphocytes, and cortisol suppression. The resulting E50-values (concentrations to induce 50% of the maximum effect) concurred with the reported values of in vitro receptor binding affinities. The duration of the systemic effects were short because of the short half-life of the drug. Cumulative cortisol suppression increased with dose administration and ranged from 20% to 36%. The PK/PD simulations resulted in a smaller degree of cortisol suppression for the drug administered at 10 PM. The cumulative change from baseline was slightly smaller for the effects on granulocytes and lymphocytes than those on cortisol. This information promotes the comparison with other inhaled glucocorticoids.

Dependency of cortisol suppression on the administration time of inhaled corticosteroids

Endogenous cortisol suppression is one of the major systemic side effects of inhaled corticosteroids in the treatment of asthma. A previously developed pharmacokinetic/ pharmacodynamic approach was used to evaluate the influence of administration time on the cumulative cortisol suppression (CCS) after single doses of the inhaled corticosteroids flunisolide and fluticasone propionate. Administration time-dependent simulations of CCS were performed with drug-specific pharmacokinetic and pharmacodynamic parameters obtained from previous clinical trials. Both drugs showed similar diurnal variation in CCS, dependent on the administration time, with maximum suppression when administered in the early morning at approximately 3 AM. The optimum administration time for minimized CCS was in the afternoon but was shifted from 3 PM for fluticasone propionate to later time points around 7 PM for flunisolide, probably because of the shorter terminal elimination half-life of flunisolide. Regarding peak to trough fluctuation, however, CCS after fluticasone propionate showed only half the administration time dependency as after flunisolide. Therefore, the ratio between CCS after flunisolide and after fluticasone propionate also followed administration time-dependent variations. This led to the conclusion that administration time has to be considered as a pivotal influential factor in clinical studies comparing CCS among different inhaled corticosteroids.

Fluticasone propionate--an update on preclinical and clinical experience

Fluticasone propionate (FP) is a novel androstane glucocorticoid with potent anti-inflammatory activity which has been effectively used, intranasally, as therapy for seasonal and allergic perennial rhinitis. When taken by the inhaled route, FP has shown significant therapeutic efficacy in the management of asthma. Fluticasone propionate is a highly lipophilic molecule with good uptake, binding and retention characteristics in human lung tissue. Fluticasone propionate has high glucocorticoid receptor selectivity and affinity, demonstrating rapid receptor association and slow receptor dissociation. In vitro, FP has been shown to potently inhibit T lymphocyte proliferation, cytokine generation, tumour necrosis factor alpha (TNF-alpha)-induced adhesion molecule expression, interleukin-5-induced eosinophilia, mucosal oedema and toluene 2,4-diisocyanate-induced mast cell proliferation, while promoting secretory leucocyte protease inhibitor production and eosinophil apoptosis. In human studies, FP has demonstrated marked vasoconstrictor potency in normal subjects and inhibited antigen-induced mucosal platelet activating factor/eicosanoid production, T lymphocytes and CD25+ cells in patients with rhinitis. Biopsy data from mild asthmatics demonstrate FP-associated reduction in CD3, CD4, CD8 and CD25 cells, with an accompanying reduction in eosinophil and mast cell markers. Clinical studies have evaluated lung function, bronchial reactivity, exacerbation rates and oral corticosteroid-sparing effect. Results show that FP has at least twice the clinical potency of beclomethasone dipropionate and budesonide. This appears to be achieved without an accompanying increase in systemic effects, suggesting a therapeutic index which may be higher than other currently available inhaled corticosteroids.