Relevance of Half-Life in Drug Design - Journal of Medicinal

Nov 7, 2017 - This may pose challenges to achieving optimal efficacy, safety, and patient compliance. If the half-life is too ... Learning Medicinal C...
0 downloads 14 Views 2MB Size
Subscriber access provided by READING UNIV

Perspective

Relevance of Half-life in Drug Design Dennis A Smith, Kevin Beaumont, Tristan S. Maurer, and Li Di J. Med. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jmedchem.7b00969 • Publication Date (Web): 07 Nov 2017 Downloaded from http://pubs.acs.org on November 8, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Medicinal Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Relevance of Half-life in Drug Design

Dennis A. Smith,† Kevin Beaumont,‡ Tristan S. Maurer,‡ and Li Di*§ †4 The Maltings, Walmer, Kent CT14 7AR, U.K. ‡Medicine Design, Pfizer Inc., Cambridge, Massachusetts 02139, United States §Medicine Design, Pfizer Inc., Groton, Connecticut 06340, United States

Invited Miniperspective for Journal of Medicinal Chemistry

1 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Drug half-life has important implications for dosing regimen and peak-to-trough ratio at the steady state. A half-life of 12-48 hours is generally ideal for once daily dosing of oral drugs. If the half-life is too short, it may require more frequent dosing in order to maintain desired exposures and avoid unnecessarily high peak concentrations. This may pose challenges to achieving optimal efficacy, safety and patient compliance. If the half-life it too long, the time over which accumulation and subsequent elimination occur may be prolonged. This may pose problems with managing adverse effects and the design of efficient clinical trials. Half-life is a key parameter for optimization in research and development. Structural modification to affect clearance, and to a lesser extent volume of distribution, is the preferred means of modulating half-life. An effective approach to half-life optimization requires an understanding of the many pitfalls associated with its estimation and interpretation.

Key Words Half-life, pharmacokinetics, drug design, clearance, volume of distribution, structure-property relationship, dosing regimen, peak-to-trough ratio

2 ACS Paragon Plus Environment

Page 2 of 47

Page 3 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Introduction Half-life (t1/2) is a useful parameter to put a rate of decline into an intuitive numerical context. The term was initially coined by Nobel Laureate, Ernest Rutherford, in the study of radioactive decay and this area of application remains the most common in the popular press. From a drug discovery perspective, t1/2 is defined as the time required for the concentration of a drug (typically in blood or plasma) to reduce to half of its initial value when the concentrations of the drug are in simple exponential (log-linear) decline. Despite its relatively simple visualisation and mathematical derivation, the role of t1/2 in drug discovery and development is surprisingly complex.

In our recent perspective on volume of distribution1, the importance of drug

distribution in defining t1/2 was exemplified, as this is likely to be the major determinant of dosing frequency for many drugs. For ease of administration, orally delivered small molecule drugs are often required to be dosed on a once daily basis. Since the duration of pharmacodynamic (PD) effect will be determined, in part, by the duration of exposure, t1/2 is often a critical parameter in drug design. When considering a design effort to optimize t1/2, there are several key questions that need to be understood, namely: •

What is the optimal t1/2 to support once daily administration?



How can t1/2 be modulated?



What are the important considerations in measuring t1/2 during compound optimization?

This manuscript will provide a review and perspective of t1/2 in the context of these generally relevant drug design questions.

3 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 47

Basic concepts, relevance to drug design and optimal space The classical one compartment pharmacokinetic (PK) model provides a simple and generally applicable framework within which to consider the relevance of t1/2 to drug design (Figure 1a). In this model, a dose can be administered directly via IV bolus (Figures 1b and 1c) or absorbed via a first order rate (ka) following oral dosing (Figure 1d and 1e). Drug is eliminated via another first order rate constant (kel) which determines the t1/2 of the drug under conditions in which the rate of absorption is more rapid than the rate of elimination. Under this condition, kel can be estimated from the slope obtained after plotting the Log10 of concentration versus time as indicated in Figures 1c and 1e. The kel parameter determines the t1/2 of a drug according to equation 1.

t/ = ln (2)⁄k 

Eq. 1

T1/2 is an important determinant in almost every aspect of a PK profile. Given the fixed relationship between t1/2 and kel, its relevance can be seen throughout the equations of table 1 (keeping in mind that each instance of kel is exactly equivalent to ln(2) / t1/2). For example, dose requirements to achieve a steady-state maximum or minimum target exposure will be related to the t1/2 through kel. In addition, as a rule of thumb, the time to steady-state exposure will require 3 – 5 t1/2 of multiple dosing. T1/2 is also related to the extent of steady-state accumulation through its relationship to kel. Other determinants include the absorption rate and dosing frequency as depicted in the somewhat complex equation for R in Table 1. An important concept here is that accumulation will increase as t1/2 becomes a larger fraction of the target dose interval (Figure 2A). In contrast, t1/2 is not a determinant of steady-state average concentrations (Cavg,ss 4 ACS Paragon Plus Environment

Page 5 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

See equations in Table 12). However, on a practical level, t1/2 remains an important consideration even when targeting a Cavg,ss. Perhaps the most important aspect to consider in this regard is the relationship of t1/2 to steady-state peak-to-trough exposure ratio (P/Tss). T1/2 values that are short relative to a desired dosing interval (Tau) will potentially incur safety risks due to large swings in exposure around the target steady-state average exposure. As illustrated in Figure 2B, while the Cavg,ss remains constant, P/Tss increases exponentially as t1/2 becomes a smaller fraction of the Tau. This figure is semi-generalized for a one compartment model in that it is universal under the assumption that ka is 1.8 /hr (providing for Tmax values between 1 and 2 hours). As such, a molecule dosed once daily with a t1/2 value of 4 hours (t1/2 / Tau = 0.17), is expected to have a P/Tss ratio greater than 40 (log P/Tss > 1.6). Such large swings in drug exposure around a target C,avg,ss bring the potential for safety-related issues (due to unnecessarily high maximal concentrations). In addition, large swings around a target Cavg,ss can result in suboptimal efficacy. This arises from the fact that target binding and modulation are inherently nonlinear processes, resulting in a disconnect between average exposure and average PD effect as t1/2 decreases. This disconnect is again illustrated in a semi-generalized way (assuming a ka of 1.8/hr and the simple sigmoidal Emax PD model) in Figure 3. In this case, the impact is greater when targeting higher ECxx values. For example, Figure 3 shows that a molecule dosed once daily with a t1/2 value of 4 hours (t1/2 / Tau = 0.17), is expected to achieve an average effect of 17, 39, 65 and 78 % despite maintaining Cavg,ss values equivalent to EC20, EC50, EC80 and EC90, respectively. The disconnect becomes larger as t1/2 continues to decrease to 2 hours, with an average effect being approximately 50% at a Cavg,ss value equivalent to the EC90. As such, regardless of whether the target exposure is Cmin,ss, Cavg,ss or Cmax,ss, t1/2 is an important aspect to consider in drug design as it will influence dose requirements and therapeutic index in the

5 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

manner discussed above. Although other properties (e.g. intrinsic potency and safety) preclude the establishment of absolute values, the relationships described above support a common design rule of thumb that t1/2 should be no less than 25% of the intended dose interval in order to maintain the desired PD effect and minimize P/Tss. While increasing t1/2 is a common design goal, it is important to note that exceedingly long t1/2 can also have untoward implications. For example, the number of doses required to reach a specified fraction of target steady-state levels increases as t1/2 increases. As such, t1/2 values that are exceedingly long can pose issues to rapidly achieving target concentrations in clinical trials and in the therapeutic setting. Likewise, the time it takes to wash drug out (Twash) by a specified fraction (fel) also increases as t1/2 increases (Table 1). For example, the time it takes to wash out (Twash) 50, 75 and 90% (fel = 0.5, 0.75 and 0.9) of a drug with a 48 hour t1/2 is 48, 96 and 160 hours, respectively. Here, t1/2 values that are exceedingly long can pose problems in managing adverse effects as they arise. For these reasons, it is generally desirable to target half-lives that are no more than 48 hours for once daily dosing. It is also difficult to run cross-over studies when the drug t1/2 is long. In some cases, however, longer terminal t1/2 may result from slow but extensive distribution into and out of particular tissues. The extreme example of prolonged t1/2 of this nature is seen with the oral bisphosphonates. These drugs avidly bind to their therapeutic target which is bone. After oral absorption, the apparent clearance from plasma is rapid by binding to bone quasi-reversibly and excretion by the kidney (20-80%). Alendronate, for example, has an initial volume of distribution excluding bone of 0.4 L/kg and a plasma clearance of 10 ml/min/kg indicative of active secretion in the kidney3. After bone uptake, the return of the drug to the circulation is controlled by the resorption of the bone. This renders the terminal t1/2 of bisphosphonates very long (5-10 years). Drug is not detectable in plasma during this period, but the long t1/2 can be

6 ACS Paragon Plus Environment

Page 6 of 47

Page 7 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

determined by monitoring the excretion of the drug in urine. In PK terms, the system has moderate clearance but the affinity for bone leads to a remarkably high volume of distribution (approximately 70,000 L/kg) and leading to long t1/2. Another similar example was observed during the development of the Cholesteryl Ester Transfer Protein (CETP) inhibitor, anacetrapib, the terminal t1/2 in healthy volunteers and patients increased with increasing length of exposure4. Initially thought to have a clinical t1/2 of around 60 hours, long term studies showed residual concentrations of anacetrapib remain significantly elevated after the cessation of therapy. Twelve weeks after stopping treatment of anacetrapib, plasma concentrations were approximately half the on-treatment levels. In a small subset of patients, anacetrapib was detectable in plasma as late as four years after the last dose. The lipophilicity of the compound prompted further pre-clinical studies which illustrated that white fat retention, rather than brown fat retention, was the probable source of the long t1/25. Not surprisingly simple PK reveals why long term therapy exposes the long terminal t1/2. If twelve weeks is a broad estimate of the terminal t1/2 and it relates to redistribution from white fat, then steady state concentrations of anacetrapib will not be achieved in white fat until approximately one year of dosing (approximately 4-fold the terminal or white fat t1/2). Amiodarone is another drug which accumulates into white fat and exhibits a long terminal t1/26-9. The difficulty of measuring the long t1/2 of multiphasic compounds in short studies (single doses) is illustrated by the variation in t1/2 values quoted. Initially values of 4 and 17 hours were published for the t1/2 following IV and oral administration. Subsequent publications have described the t1/2 as 35 hours before a more realistic value of around 40 days. On long term therapy, it is probably nearer 100 days. A study after chronic dosing in rabbits showed that amiodarone total concentrations were the highest in fat tissue followed by lung, liver, and muscle, with the lowest concentration in serum and only

7 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 47

traces in the brain 10. Interestingly, the active metabolite desethylamiodarone had similar concentrations in liver, lung, and kidney to the parent drug but lower accumulation in fat tissue. The problems with t1/2 determination lead to similar huge variation in volume of distribution with values from 0.9 and 148 L/kg, with the latter probably the most accurate estimate.

Optimization of half-life Volume and clearance modulation Keeping with the one-compartmental model, t1/2 can be appreciated in terms of its determinants (Equation 2).

t/ = ln (2)V ⁄CL

Eq. 2

From this relationship, it is clear that t1/2 will be proportional to volume and inversely proportional to clearance. In turn, Equations 3 and 4 provide conceptually useful information regarding the common physiological determinants of Vss and CL. In these equations, fup, fut, Vp and Vt represent the unbound fraction in plasma, unbound fraction in tissue, volume of the plasma and volume of the tissue, respectively. Further, in the depicted well stirred model of hepatic clearance, Qh, and CLint represent the hepatic blood flow and intrinsic clearance, respectively. As covered in our previous review, Vss is most sensitive to ionization state with acids typically displaying low volumes (high binding to albumin in plasma, low fup) and bases displaying higher volumes (high partitioning to acidic phospholipids in tissues, low fut). This has important implications for design, as the requisite clearance values in order to achieve a 8 ACS Paragon Plus Environment

Page 9 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

particular t1/2 will often depend on ionization state for this reason. For example, in order to achieve a t1/2 of 12 hours, acids (Vss ~0.2 L/kg), neutrals (Vss ~ 1 L/kg) and bases (Vss ~ 4 L/kg) will typically require clearances at least as low as 0.2, 1 and 4 ml/min/kg, respectively (1, 5 and 20 % of human hepatic blood flow, respectively).

V = Vp + ∑fu ⁄fu Vt

Eq. 3

CL = Qh × fu × CL ! "Qh + fu × CL ! 

Eq. 4

Optimization of CLint is often the primary focus in design efforts to improve t1/2. For oxidative metabolism, CLint is relatively easily measured in liver microsomes or hepatocytes in high throughput in vitro systems which facilitate an efficient development of structure-activityrelationships (SAR). In addition, because it is also a determinant of Cavg,ss and hepatic extraction, optimization efforts focused on CLint are also likely to pay additional dividends with respect to dose optimization.

While tempting when looking at Vss or CL in isolation, optimization of t1/2 through changes in plasma protein binding is not generally a productive endeavour. As can be seen in Equations 3 and 4, decreases in fup may improve clearance, but may also reduce Vss. Likewise, increases in fup may improve Vss, but will also tend to increase clearance. Consequently, modulation of plasma protein binding cannot modulate t1/2, unless clearance is reduced or Vss increased

9 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 47

independent of the modulation of fraction unbound. It is a common misconception that reducing plasma protein binding can result in longer t1/2. Generally, in order to reduce plasma protein binding, lipophilicity is reduced. It is well known that reducing lipophilicity is likely to reduce intrinsic metabolic clearance with an overall reduction of the clearance of a molecule. It is this reduction in metabolic clearance that drives any t1/2 improvement, not the reduction in plasma protein binding. The impact of plasma protein binding on drug disposition has been thoroughly reviewed previously 11, 12.

Modified release formulations Another option for extending t1/2 is the use of modified release formulations (MR). Such formulations can effectively extend the PK t1/2 of molecules with suboptimal values by providing for sustained release along the intestinal tract. This has been successfully used in a number of marketed drugs13-16. The principles of dose size and frequency of administration as well as the potential value of MR formulations was demonstrated early using the two calcium channel blockers, nifedipine and amlodipine (Table 2)1,

17-20

. Amlodipine exhibits t1/2 value that is

amenable to once daily (QD) administration1. The short t1/2 of nifedipine necessitated a three times daily (TID) dose regimen.

Nifedipine was eventually developed as a slow release

formulation (zero order) termed nifedipine gastrointestinal therapeutic system (GITS) to allow once a day administration with doses of 30 and 60 mg21. A more recent example of the application of modified release to extend the t1/2 and reduce dose frequency is the extended release formulation of tofacitinib22, 23. Tofacitinib is JAK inhibitor with a human t1/2 of 3 hours. It is presented as an immediate release tablet formulation with a dose of 5 mg twice daily. A modified release device using an Extrudable Core System (ECS) has been developed which

10 ACS Paragon Plus Environment

Page 11 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

drives similar AUC to 5 mg twice daily with an 11 mg once daily device. Although a thorough review of the practical considerations for solving t1/2 issues with this technology is beyond the scope of this manuscript, it is important to point out that there are limitations to what this technology can achieve. Generally, the compound must have a degree of aqueous solubility and relative rapid dissolution, since released drug will need to dissolve rapidly under all environments encountered in the human gastrointestinal (GI) tract. Membrane permeability also needs to be sufficient to allow for rapid absorption of dissolved compound across the GI tract wall. Finally, affinity for efflux transporters at the wall of the GI tract needs to be minimized. Generally, immediate release devices present very high drug concentrations to the upper GI tract which are capable of saturating any potential efflux mechanism. Modified release devices steadily present much lower concentration throughout the GI tract, which may be below those required to saturate efflux proteins in the GI tract. Aqueous solubility, permeability and efflux transport are important considerations for MR, since the drug needs to be as well absorbed in the colon as it is in the small intestine. The human colon has less fluid, a lower surface area for absorption and higher expression of efflux transporters (such as P-glycoprotein) than the small intestine. Finally, due to the need for excipients to drive MR, the drug loading (dose) for a MR device needs to be limited. Thus, potency and clearance properties of a candidate for MR need to be optimized to maintain dose requirements.

Prodrug approach Finally, because a number of drugs generate long lived active metabolites and therefore act as (accidental) prodrugs, the idea of adding a pro-moiety to a promising molecule that lacks

11 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 47

duration due to a short t1/2 is regularly raised. This is not likely to be successful for the following reasons: 1) The design of optimal prodrugs from a preclinical perspective is difficult. In vitro assays need to be available in preclinical species and human to examine intestinal, hepatic and extrahepatic hydrolysis of the prodrug moiety to the active drug. Preclinical species can modulate the prodrug to a significantly different extent to humans, making extrapolation and prediction fraught with difficulty. 2) If the short t1/2 of the active drug is due to high unbound intrinsic clearance or total clearance, the total daily dose of the prodrug is likely to be very high, even if there is extensive conversion of the pro-moiety to the active drug. 3) The rate of metabolism of the prodrug to the active drug must be slow and this raises the likelihood of non-productive clearance of the prodrug by metabolism via other routes or other mechanisms leading to the need for higher doses to compensate 24.

Pitfalls in half-life optimization ‘Flip-flop’ kinetics Most drugs are absorbed at a rate (ka) which exceeds the elimination rate (kel).

In this

circumstance, the observed t1/2 is defined as previously described in equations 1 and 2. However, in certain circumstances, the rate of absorption can be slower than the rate of elimination, resulting in a t1/2 which is determined by the rate of absorption rather than the rate of elimination. This phenomenon, known as ‘flip-flop’ kinetics, can occur when dealing with molecules that have particularly poor permeability or solubility. In these cases, the t1/2 following oral administration may exceed that observed following IV dosing due to absorption or

12 ACS Paragon Plus Environment

Page 13 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

dissolution rate-limited absorption.

This phenomenon is not uncommon in pre-clinical PK

studies (particularly in rodents) on molecules with suboptimal PK properties in the course of research and development. There are several pitfalls that arise in this situation. For example, as a compound progresses in development, formulations may improve the dissolution rate and hence the t1/2 value will collapse to that defined by elimination. In addition, accumulation will often be much less than predicted from an oral t1/2 in this circumstance. This results from the fact that there is a discrete window of time over which ka is operational as determined by the physiological intestinal transit times and physiochemical drug properties. Further, this highlights the limitation of empirical compartmental pharmacokinetic models which don’t explicitly account for such physiological limitations and may lead to erroneous expectations when these limitations are not properly considered in experimental design and interpretation.

Finally,

attention to oral t1/2 in this circumstance can lead to a mis-prioritization in favour of less soluble or poorly permeable molecules that will ultimately suffer from absorption challenges in the clinic where formulation options are more limited (e.g. limited absolute bioavailability with an associated high degree of inter and intra-subject variability). A much more appropriate use of ‘flip-flop’ kinetics to extend t1/2 is the use of modified-release formulations.

Often PK data is extrapolated to infinity using the apparent terminal t1/2. This is an exercise in which the concentration in the last time point collected is divided by the apparent kel value (typically determined though regression of last few data points) in order to infer the remainder of the AUC not explicitly measured through infinite collection of samples. When the oral t1/2 exceeds the IV one, there is a dilemma. Extrapolating to infinity with the oral t1/2 is using

13 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 47

potentially the ka and not the kel in preclinical species, which may not be the case in the human situation leading to much shorter t1/2.

Here, it is worth pointing out that the t1/2 of metabolite can be equally confounded by formation –rate limited disposition (analogous to ‘flip-flop’ kinetics) when the rate of formation from parent is slower than the rate of metabolite elimination25. Thus, a metabolite may, if administered discretely show a shorter t1/2 than parent, but when the parent is administered, the decline of the metabolite will parallel that of the parent. Consequently, the t1/2 of a metabolite when derived from administration of a parent compound can never be shorter than that of the parent. Of course, some metabolites can possess longer t1/2 values than the parent, normally due to their being more stable than the parent to metabolic processes or they exhibit Vss greater than the parent. Prime examples are the N-dealkylation reactions of tertiary amines where the generated metabolites exhibit lower intrinsic clearance by subsequent metabolism combined with potentially higher Vss values.

Dose-dependent half-lives Dose-dependent PK usually occurs due to saturation of absorption or clearance processes as doses are increased. The most common dose-dependent processes observed are saturation of clearance in pre-clinical toxicity studies, where doses are high. In many low dose PK studies, the unbound concentrations of drug are likely to be below the Km (e.g. substrate concentration that gives half maximal velocity of transport or metabolism) for the clearance process. Under these conditions, the parameters of clearance and t1/2 will be dose independent (i.e., will not change with increasing doses). However, if doses are used that drive to unbound concentrations

14 ACS Paragon Plus Environment

Page 15 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

that exceed the Km, then the PK will change from first order to zero order processes. If the dose is high enough for zero order kinetics to be observed, the plot of Cp versus time will be linear and a Log10Cp against time will be curvilinear. Once concentrations fall below the Km, the normal linear decline is observed for the Log10Cp against time curve. Therefore, clearance will decrease and t1/2 will increase with increasing dose (Figure 4). Consequently, it is important that PK studies which are extrapolated to the clinical situation are completed at doses that generate unbound plasma concentrations below the Km of the rate-determining clearing process. In rare instances, nonlinear pharmacokinetics may result from saturation of plasma protein binding. In these cases, concentrations exceeding the binding capacity of plasma proteins lead to higher free fractions which, in turn, will increase clearance and increase volume of distribution (all else being equal). The potential for this phenomenon should be considered in circumstances where plasma drug concentrations approach physiological levels of the major binding proteins in plasma (e.g. ~ 15 µM and 600 µM for α1-acid glycoprotein and albumin, respectively). Valproic acid is an example of a drug associated with therapeutic exposures in this range (total plasma concentrations ~ 350 – 700 µM) and which also exhibits saturable protein binding and nonlinear clearance. In other, more complex situations, saturable tissue binding may occur. In these cases, a decrease in volume of distribution may be seen at concentrations which begin to saturate tissue binding.

Multicompartmental PK profiles In the above discussion, a one compartmental model, which provides for a monoexponential decline in drug concentration, was assumed in order to discuss basic concepts related to the t1/2. This one compartment situation requires essentially instantaneous (or very rapid) distribution of

15 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 47

a compound into tissues. In the case of most acidic drugs, the majority of total drug is confined to plasma and so a single compartment model is often observed. However, for lipophilic neutral and basic compounds, this is seldom the case. The decline of plasma concentration – time course, particularly after IV injection, typically shows at least two or more discrete phases when viewed on a log10Cp against time plot (Figure 5). This is due to the rate of entry into compartment 2 from the first compartment (k12) and from compartment 2 back to the first compartment (k21). Consequently, t1/2 can no longer be described as the time taken for the plasma concentrations to fall by half, as there are at least two distinct phases with different rates of decline. The initial decline is often viewed as a distribution phase followed by a terminal elimination rate (β). In this model, Vss, CL, β and terminal t1/2 are defined by equations 5, 6, 7 and 8 respectively. V = V1 $1 + CL = V1 × k 

k "k % 

Eq. 5 Eq. 6

β = 0.5*(k + k  + k  ) − ,(k + k  + k  ) − 4 × k  × k  .

terminal t/ = ln (2)⁄β

Eq. 7 Eq.8

Although the relationships appear to be much more complicated, it is important to note that, β is closely approximated by CL / Vss (i.e. another way to define kel in the one compartment model) so long as the distributional rates exceed that of the elimination rate. As such, under most circumstances, the basic concepts for modulating terminal half-life according to Vss and CL remain relevant with this more complicated profile. In addition, it is not uncommon for the multi-exponential decline to be lost or attenuated following oral administration due to the relative rates of absorption versus distribution. However, in a manner analogous to ‘flip-flop’

16 ACS Paragon Plus Environment

Page 17 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

kinetics, these rates may be difficult to interpret once the distribution rates become slower than the elimination rate. In some cases, a prolonged terminal half-life may be observed after concentrations have declined significantly providing for a portion of the concentration-time curve containing very little of the overall area under the curve (AUC). In this case, while the terminal elimination phase t1/2 may well approach a value that would be associated with once daily administration in humans, it would not lead to significant accumulation of drug exposure at steady state owing to the fact that it contributes very little to the overall AUC and clearance of the drug. Such a drug would still need to be administered more than once daily to achieve appropriate steady state exposure. To account for this complexity in a practical way, the concept of an “effective” or “operational” half-life has been developed which describes half-life in relation to the observed or predicted degree of accumulation and shape of the concentration-time profile upon steady-state dosing26, 27. This is exemplified by the NK2 antagonist UK-224671 [(S)-(1)-4-(1-h2-[1-(cyclopropylmethyl)-3-(3,4-dichlorophenyl)-6-oxo-3piperidyl]ethyljazetidin-3-yl)-1-piperazine sulphonamide] 28. The PK of UK-224671 in human following IV administration is characterized by a multiple compartmental profile (at least 3 phases, Figure 6). The terminal elimination phase t1/2 is calculated at 28 hours. However, this terminal phase describes approximately 20% of the AUC. Consequently, a once daily regimen for this compound would not accumulate to an appreciable steady state trough concentration. Paradoxically, this compound would be best suited to three times daily administration, despite an apparently once daily t1/2.

Estimation of half-life

17 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 47

The normally quoted (terminal) t1/2 is commonly determined by estimating the final slope of the plasma concentration-time curve. This now brings error in data analysis. After an IV dose, time points will often be taken with three sampling in the first hour and then hourly up to 4 hours followed by samplings at perhaps 6, 8 and 24 hours. The sampling time of the later time points is mostly governed by practical reasons. The t1/2 calculation will probably use the final three or so time points that have measurable concentrations. In the estimation of t1/2, at least three points should be used. This is an inadequate number to gain a robust estimate of t1/2. Having any drug detectable at 24 hours post dose usually significantly prolongs the apparent t1/2, which means assay sensitivity between compounds may be a major factor in selecting compounds rather than real PK differences. The assay differences may be compounded across species. For instance rodent studies usually have a lower number of samples of the different time points with a smaller volume of plasma samples than larger species like dog. Therefore, the larger species has the potential to generate a significantly longer t1/2 than rodents simply because of analytical sensitivity. The same analytical sensitivity issue will be seen in the area under curve (AUC) and area under the first moment curve (AUMC) so that clearance and volume will also be defined by the analytical sensitivity issue. Clearly, these are aspects of the PK study design that may complicate the estimation of the true t1/2 of the compound being tests and affect the extrapolation to human PK.

Often oral studies use higher doses than IV studies for safety reasons (oral exposure is usually less that for the same intravenous dose due to oral bioavailability). This adds to the complexity of t1/2 determination since the oral study concentrations will be higher than IV dosing in the terminal phase, assuming the same analytical method is used. Consequently, the ‘true’ terminal

18 ACS Paragon Plus Environment

Page 19 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

t1/2 of the molecule could be exposed with the chance of therefore a significantly longer apparent t1/2 after oral administration. The difficulty in interpretation will then be whether this is a fundamental parameter of the compound, or driven by slow dissolution in the GI tract (as discussed previously). The problems of detection of drug in later samples or not sampling long enough also apply to the estimation of CL and Vss. Although model free approaches are used largely in drug discovery there is still a need for determination of the final elimination phase to extrapolate AUC values to infinity. Underestimating the terminal t1/2 therefore results in an underestimation of Vss and an overestimation of CL. When this is exaggerated by the low intravenous dose compared to the oral dose in a crossover study overestimates (including >100%) will be made for bioavailability.

Population Variability of Half-Life Variability of the t1/2 across populations can be of major concern. Variability in t1/2 should be a consideration in any patient population for which a high degree of variability in CL is expected due to the known mechanisms of clearance (hepatic impairment, renal impairment, genotype, ethnicity, age, etc.).

In other cases, variability can arise due to much more complicated

mechanisms. For example, high variation in t1/2 is accentuated when the drug’s t1/2 is partially dependent on reversible metabolism (metabolic conversion back to parent). This is illustrated by benoxaprofen, a nonselective COX inhibitor (Figure 7). Most of the drugs in this class contain a carboxylic acid moiety. It is difficult to obtain a long t1/2 with such agents due to the low Vss. Benoxaprofen was developed because the plasma elimination t1/2 of the drug was in the range of 19-26 h29, almost ideal for a once a day administration. However, the development studies were conducted only in younger volunteers and patients. At the time a once a day drug with a simple

19 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 47

single dose size was commercially attractive against the many other drugs in the class. Very soon after launch, toxicity problems were noted affecting the skin and liver. With both skin toxicity and hepatotoxicity, the major proportion of effects were observed in elderly patients30. When PK of these patients were studied, t1/2 extended out to 150 hours which coincided with reduced creatinine clearance31.

Acyl glucuronidation was the major clearance pathway of

benoxaprofen32. Excretion of this metabolite could occur renally and hepatobiliary. Enterohepatic recirculation of this metabolite after hepatobiliary clearance and subsequent hydrolysis in the gut back to benoxaprofen was likely to be an important contributor to the long t1/2. Gastrointestinal stasis leading to higher recirculation and lowered kidney function lowering the renal clearance of the acyl glucuronide contributed to accumulation of the drug in the elderly. Today, studies in elderly would have been conducted during drug development in addition to more detailed studies to understand the clearance and metabolism. These would have probably rendered the drug much less commercially attractive. The drug illustrates that thinking and basing decisions on a single average value for t1/2 is problematic. Drugs can exhibit quite different t1/2 in different patients and populations leading to the need for different dosing regimens.

Perspective The simplest PK definition of t1/2 as the time taken for plasma concentrations to fall by half only applies under one compartment PK following IV administration. Under other conditions, such as oral administration and multi-compartmental PK, the effect of rate constants other than kel need to be taken into consideration, however, Vss and CL remain important factors to consider in any approach to optimize t1/2. 20 ACS Paragon Plus Environment

Page 21 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

We introduced this perspective by posing several important questions with respect to t1/2 during a drug discovery project. •

What t1/2 is optimal to support once daily administration?

Deriving dose regimen from t1/2 should drive to minimize peak-to-trough ratio, which will limit the potential for Cmax driven side-effects and lack of efficacy at Cmin. An acceptable peak-totrough ratio is below 4, which leads to a t1/2 greater than 12 hours (if the PK is one compartment). T1/2 of 24 to 48 hours would further minimize peak-to-trough ratio and potentially allow ‘forgiveness’ of a single missed dose in minimizing loss of steady state exposure. Compounds with t1/2 values greater than this can lead to potential issues in terms of time to achieve steady-state, washout in the event of adverse effects and longer drug development time.



How can t1/2 be modulated?

T1/2 is determined by the plasma CL and Vss of a compound. To improve t1/2, medicinal chemistry strategy should mainly focus on reducing plasma CL (e.g., increase metabolic stability, reduce renal clearance, reduce hepatic uptake by active transporters and biliary clearance) and/or secondly increasing Vss1 . Generally, modulation of plasma protein binding in the absence of changing CL or Vss will not affect t1/2. This is a consequence of the relationship of t1/2 to CL and Vss (Equation 2). Changing fraction unbound in plasma will often change CL and Vss in equal and opposite directions. Should a clinical agent not possess the appropriate t1/2 to suit its proposed dose regimen, it is possible to improve the t1/2 by modifying the delivery of the compound. Subject to appropriate properties of aqueous solubility sufficient to support dose

21 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

requirements, high membrane permeability and low efflux transport, a compound can be formulated as MR to enhance t1/2. This formulation will constantly release the compound throughout its passage through the GI tract. Effectively, the MR formulation gives zero order delivery of the compound and results in ‘flip-flop’ PK, meaning t1/2 in plasma will reflect the rate of delivery.



What are the important considerations for measuring t1/2 during compound optimization?

Complexities related to the multiple rates (absorption, distribution and elimination) can confound efforts to optimize t1/2 by obscuring the interpretation of preclinical and clinical data. A careful analysis of PK profiles under different dose regimens and formulations (particularly IV and oral administration) are required in order to put an observed t1/2 into a context useful to design. In addition, nonlinearities in absorption and elimination can equally confound the interpretation of t1/2. PK studies employing a range of doses and exposure levels are often required to detect these issues. Finally, population variability in t1/2 (particularly in the clinic) can confound the interpretation of t1/2. Here an understanding of the mechanisms of clearance and distribution specific to a molecule and how those mechanisms vary across potential patient populations is important to ensure the appropriate interpretation of observed t1/2.

Acknowledgements The authors greatly appreciate the technical assistance from David Tess.

22 ACS Paragon Plus Environment

Page 22 of 47

Page 23 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Corresponding Author Information: Li Di Tel. 860-715-6172 Email: [email protected]

23 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Biographical Sketch Dr. Dennis A. Smith worked in the pharmaceutical industry for 32 years after gaining his Ph.D from the University of Manchester. During this period he directly helped in the Discovery and Development of eight marketed NCEs. More recently his roles include advisory boards or expert panels with many drug research based organisations, both industrial and academic, including Medicines for Malaria Venture and Cancer Research UK. His research interests and publications span all aspects of Drug Discovery and Development particularly where drug metabolism knowledge can impact on the design of more efficacious and safer drugs.

Mr. Kevin Beaumont has worked extensively in the discovery and development drug metabolism field, throughout his 30 years in the pharmaceutical industry. His major area of expertise is in the modulation of physicochemistry to affect drug disposition and prediction of human PK. He is author on over 40 peer reviewed publications. Overall, Kevin has worked on many drug discovery and development projects throughout his career. He has been responsible for the DMPK input to at least 30 FIH studies as well as 10 Phase II compounds, including 1 marketed agent. Kevin now provides DMPK input to the Inflammation and Immunology and Rare Disease Research Units, based in Cambridge Massachusetts.

Dr. Tristan S. Maurer received his PharmD from the University of Georgia in 1993 and his PhD from the University of Buffalo, State University of New York in 1999. During his 18 year tenure with Pfizer, his work has focused on the development and application of quantitatively rigorous, biologically-based methods to predict human pharmacokinetics & pharmacodynamics 24 ACS Paragon Plus Environment

Page 24 of 47

Page 25 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

from preclinical data. He has co-authored over 60 manuscripts illustrating the utility of these methods to drug design and early clinical development. Currently, Dr. Maurer sits on the Medicine Design leadership team responsible for scientific and operational strategies spanning from idea to loss of exclusivity. He also leads a modelling and simulation group responsible for both computational chemistry and quantitative translational pharmacology across Pfizer’s small molecule portfolio.

Dr. Li Di has about 20 years of experience in the pharmaceutical industry including Pfizer, Wyeth and Syntex. She is currently a research fellow at Medicine Design Department, Pfizer Global Research and Development, Groton, CT. Her research interests include the areas of drug metabolism, absorption, transporters, pharmacokinetics, blood–brain barrier, and drug-drug interactions. She has over 120 publications including two books and presented more than 80 invited lectures. She is a recipient of the Thomas Alva Edison Patent Award, the New Jersey Association for Biomedical Research Outstanding Woman in Science Award, the Wyeth President’s Award, Peer Award for Excellence and Publication Award.

25 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abbreviations AUC = area under the curve C0 = initial drug concentration Cavg,ss = average steady state drug concentration Cavg,ss, u = average unbound steady state drug concentration Cavg,ss = steady state average concentration Cmin, ss = steady state minimum concentration Cmax,ss = steady state maximum concentration CL = clearance CLint = intrinsic clearance CLint,u = unbound intrinsic clearance Cmax = maximum drug concentration Cmin = minimum drug concentration Cp = drug concentration in plasma ECxx = Effective concentration at certain “xx” percentage of maximal, e.g., EC50 is the concentration that causes 50% of maximum effect fa = fraction absorbed fel = fraction of drug being washed out fg = fraction escaped from gut 26 ACS Paragon Plus Environment

Page 26 of 47

Page 27 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

fss = fraction of target steady-state level fup = unbound fraction in plasma fut = unbound fraction in tissue GITS = gastrointestinal therapeutic system IV = intravenous ka = first-order absorption rate constant kel = first-order elimination rate constant Km = substrate concentration that gives half maximal velocity of an enzymatic reaction MR = modified release formulation MRT = mean residence time PBPK = physiologically based pharmacokinetic modelling PD = pharmacodynamics PK = pharmacokinetics P/Tss = steady-state peak-to-trough concentration ratio QD = once daily Qh = hepatic plasma flow t = time t1/2 = half-life

27 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Tau = dosing interval TID = three times daily Tmax = time at maximum plasma drug concentration Twash = the time it takes to wash drug out by a specified fraction Vd = volume of distribution Vss = steady state volume of distribution Vp = volume of the plasma Vt = volume of the tissue

28 ACS Paragon Plus Environment

Page 28 of 47

Page 29 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

References 1.

Smith, D. A.; Beaumont, K.; Maurer, T. S.; Di, L. Volume of Distribution in Drug Design. J. Med.

Chem. 2015, 58, 5691-5698. 2.

Wagner, J. G.; Northam, J. I.; Alway, C. D.; Carpenter, O. S. Blood levels of drug at the

equilibrium state after multiple dosing. Nature 1965, 207, 1301-1302. 3.

Lin, J. H. Bisphosphonates: A review of their pharmacokinetic properties. Bone (N. Y.) 1996, 18,

75-85. 4.

Mohammadpour, A. H.; Akhlaghi, F. Future of Cholesteryl Ester Transfer Protein (CETP)

Inhibitors: A Pharmacological Perspective. Clin. Pharmacokinet. 2013, 52, 615-626. 5.

Hartmann, G.; Kumar, S.; Johns, D.; Gheyas, F.; Gutstein, D.; Shen, X.; Burton, A.; Lederman, H.;

Lutz, R.; Jackson, T.; Chavez-Eng, C.; Mitra, K. Disposition into Adipose Tissue Determines Accumulation and Elimination Kinetics of the Cholesteryl Ester Transfer Protein Inhibitor Anacetrapib in Mice. Drug Metab. Dispos. 2016, 44, 428-434. 6.

Holt, D. W.; Tucker, G. T.; Jackson, P. R.; Storey, G. C. Amiodarone pharmacokinetics. Am. Heart

J. 1983, 106, 840-847. 7.

Plomp, T. A.; Van Rossum, J. M.; Robles de Medina, E. O.; Van Lier, T.; Maes, R. A. A.

Pharmacokinetics and body distribution of amiodarone in man. Arzneim.-Forsch. 1984, 34, 513-520. 8.

Anastasiou-Nana, M.; Levis, G. M.; Moulopoulos, S. Pharmacokinetics of amiodarone after

intravenous and oral administration. Int. J. Clin. Pharmacol. Ther. Toxicol. 1982, 20, 524-529. 9.

Latini, R.; Tognoni, G.; Kates, R. E. Clinical pharmacokinetics of amiodarone. Clin. Pharmacokinet.

1984, 9, 136-156. 10.

Kannan, R.; Miller, S.; Singh, B. N. Tissue uptake and metabolism of amiodarone after chronic

administration in rabbits. Drug Metab. Dispos. 1985, 13, 646-650.

29 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

11.

Smith, D. A.; Di, L.; Kerns, E. H. The effect of plasma protein binding on in vivo efficacy:

misconceptions in drug discovery. Nat. Rev. Drug Discovery 2010, 9, 929-939. 12.

Liu, X.; Wright, M.; Hop, C. E. C. A. Rational Use of Plasma Protein and Tissue Binding Data in

Drug Design. J. Med. Chem. 2014, 57, 8238-8248. 13.

Ding, H. Modified-Release Drug Products and Drug Devices In Applied Biopharmaceutics and

Pharmacokinetics, 7th ed.; Shargel, L.; Yu, A. B. C., Eds. McGraw-Hill Education: New York, NY, 2016; pp 567-613. 14.

Ummadi, S.; Shravani, B.; Rao, N. G. R.; Reddy, M. S.; Nayak, B. S. Overview on Controlled

Release Dosage Form. International Journal of Pharma Sciences 2013, 3, 258-269. 15.

Gupta, B. P.; Thakur, N.; Jain, N. P.; Banweer, J.; Jain, S. Osmotically controlled drug delivery

system with associated drugs. J. Pharm. Pharm. Sci. 2010, 13, 571-588. 16.

Keraliya, R. A.; Patel, C.; Patel, P.; Keraliya, V.; Soni, T. G.; Patel, R. C.; Patel, M. M. Osmotic drug

delivery system as a part of modified release dosage form. ISRN Pharm. 2012, 528079, 9 pp. 17.

Faulkner, J. K.; McGibney, D.; Chasseaud, L. F.; Perry, J. L.; Taylor, I. W. The pharmacokinetics of

amlodipine in healthy volunteers after single intravenous and oral doses and after 14 repeated oral doses given once daily. Br. J. Clin. Pharmacol. 1986, 22, 21-5. 18.

Abernethy, D. R.; Schwartz, J. B. Pharmacokinetics of calcium antagonists under development.

Clin. Pharmacokinet. 1988, 15, 1-14. 19.

Hamann, S. R.; Piascik, M. T.; McAllister, R. G., Jr. Aspects of the clinical pharmacology of

nifedipine, a dihydropyridine calcium-entry antagonist. Biopharm. Drug Dispos. 1986, 7, 1-10. 20.

Obach, R. S.; Lombardo, F.; Waters, N. J. Trend analysis of a database of intravenous

pharmacokinetic parameters in humans for 670 drug compounds. Drug Metab. Dispos. 2008, 36, 13851405.

30 ACS Paragon Plus Environment

Page 30 of 47

Page 31 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

21.

Grundy, J. S.; Foster, R. T. The nifedipine gastrointestinal therapeutic system (GITS). Evaluation

of pharmaceutical, pharmacokinetic and pharmacological properties. Clin. Pharmacokinet. 1996, 30, 2851. 22.

Lamba, M.; Wang, R.; Fletcher, T.; Alvey, C.; Kushner, J. I. V.; Stock, T. C. Extended-Release Once-

Daily Formulation of Tofacitinib: Evaluation of Pharmacokinetics Compared With Immediate-Release Tofacitinib and Impact of Food. J. Clin. Pharmacol. 2016, 56, 1362-1371. 23.

Lamba, M.; Hutmacher, M. M.; Furst, D. E.; Dikranian, A.; Dowty, M. E.; Conrado, D.; Stock, T.;

Nduaka, C.; Cook, J.; Krishnaswami, S. Model-Informed Development and Registration of a Once-Daily Regimen of Extended-Release Tofacitinib. Clin. Pharmacol. Ther. (Hoboken, NJ, U. S.) 2017, 101, 745-753. 24.

Beaumont, K.; Webster, R.; Gardner, I.; Dack, K. Design of ester prodrugs to enhance oral

absorption of poorly permeable compounds: Challenges to the discovery scientist. Curr. Drug Metab. 2003, 4, 461-485. 25.

Houston, J. B.; Taylor, G. Drug metabolite concentration-time profiles: influence of route of drug

administration. Br. J. Clin. Pharmacol. 1984, 17, 385-94. 26.

Sahin, S.; Benet, L. Z. The Operational Multiple Dosing Half-life: A Key to Defining Drug

Accumulation in Patients and to Designing Extended Release Dosage Forms. Pharm. Res. 2008, 25, 28692877. 27.

Boxenbaum, H.; Battle, M. Effective half-life in clinical pharmacology. J. Clin. Pharmacol. 1995,

35, 763-766. 28.

Beaumont, K.; Harper, A.; Smith, D. A.; Abel, S. Pharmacokinetics and metabolism of a

sulphamide NK2 antagonist in rat, dog and human. Xenobiotica 2000, 30, 627-642. 29.

Nash, J. F.; Carmichael, R. H.; Ridolfo, A. S.; Spradlin, C. T. Pharmacokinetic studies of

benoxaprofen after therapeutic doses with a review of related pharmacokinetic and metabolic studies. J. Rheumatol., Suppl. 1980, 6, 12-19.

31 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

30.

Halsey, J. P.; Cardoe, N. Benoxaprofen: side-effect profile in 300 patients. Br. Med. J. (Clin. Res.

Ed.) 1982, 284, 1365-1368. 31.

Hamdy, R. C.; Murnane, B.; Perera, N.; Woodcock, K.; Koch, I. M. The pharmacokinetics of

benoxaprofen in elderly subjects. Eur. J. Rheumatol. Inflammation 1982, 5, 69-75. 32.

Chatfield, D. H.; Green, J. N. Disposition and metabolism of benoxaprofen in laboratory animals

and man. Xenobiotica 1978, 8, 133-144.

32 ACS Paragon Plus Environment

Page 32 of 47

Page 33 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

List of Tables and Figures (All are original) Table 1. Key equations which characterize the one-compartment PK model Table 2. Comparison of nifedipine and the long t1/2 calcium channel blocker amlodipine. Daily dose size reflects the potency (IC50) and CLint,u of the drugs. Dose frequency reflects the t1/2. Figure 1. Basic characteristics of a one compartment PK model. Panel A illustrates the model structure shown with first order absorption and elimination rates. Panels b and c illustrate the characteristic exponential decline following IV bolus administrations. Panels d and e illustrate the absorptive phase followed by monoexponential decline that is typical following oral administration. Figure 2. Semi-generalized relationships of the extent of steady-state accumulation (A) and peakto-trough concentration ratio (B) (plotted a log P/Tss values) to half-life-to-dose interval ratio (t1/2 / Tau). Relationships were simulated using the corresponding R(t) and P/Tss equations in Table 1 and assuming an absorption rate of 1.8/hr. A: Accumulation (as calculated as the ratio of AUC at steady-state to AUC following a single dose) increases as the as t1/2 becomes a smaller fraction of the target dose interval. B: While the steady-state average concentration is constant, P/Tss increases exponentially as t1/2 becomes a smaller fraction of the target dose interval. Figure 3. Semi-generalized relationship between average pharmacodynamics effect (expressed as a fraction of maximal) and half-life-to-dose interval ratio (t1/2 / Tau). Relationship was simulated using the C(t)ss equation of table 1 assuming a ka of 1.8/hour. Simulations depict 4 dosing levels designed to deliver Cavg,ss values equivalent to arbitrarily selected EC20, EC50, EC80 and EC90 values. Actual average effect was calculated as the dose-interval normalized AUC under an effect curve defined by (100 x C(t)ss) / (C(t)ss + EC50). Volume was varied in order to affect t1/2 / Tau while maintaining a constant Cavg,ss. 33 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. The plasma concentration versus time curve for a compound following escalating oral administration exhibiting dose-dependent (non-linear) PK due to saturation of clearance with increasing dose. Figure 5. Basic characteristics of a two compartment PK model. Panel A illustrates the model structure shown with first-order rates of absorption, elimination and distribution. Panel b illustrates the characteristic biexponential decline in plasma concentrations following IV bolus administration. Figure 6. The PK of UK-224,671 in human following single intravenous administration of 0.1 mg/kg. Panel a: Overall PK profile. Panel b: Isolated initial phase profile. Panel c: Isolated terminal phase profile

34 ACS Paragon Plus Environment

Page 34 of 47

Page 35 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Table 1: Key equations which characterize the one-compartment PK model2 Parameter C(t)ss

Equation C(t) =

F5 D5 k8 e:;

CA8B, =

F5 D5 k8 e:;

CA !, =

F5 D5 k8 1 1 ∙ ∙G − H V k 8 − k  1 − e:; 1 − e:;? 8>

Steady-state concentration at time =t Cmax,ss Steady-state maximum concentration

Cmin,ss Steady-state minimum concentration

C8IJ,KK =

Cavg,ss

F5 D5 CL tau

Steady-state average concentration t A8B, =

tmax,ss

1 k 8 1 − e:;  ∙ ln 9 @ k 8 − k  k  (1 − e:;? 8> )

Steady-state post dosing time at which maximum concentration occurs

R

R=9

e:; − e;

e;? 8> − 1 ∙ e; − 1 ∙ G

trough concentration ratio

Twash

TP8Q = −

ln(1 − f  ) ln (2)⁄R/

The time it takes to “wash out” a specified fraction (fel) of the drug concentration after sessation of dosing All relationships correspond to a one compartment pharmacokinetic model. Fpo, Dpo, Vss, CL and tau represent oral bioavailability, oral dose, steady-state volume of distribution, clearance rate and dosing interval, respectively. kel and ka represent the first order rates of elimination and absorption.

36 ACS Paragon Plus Environment

Page 36 of 47

Page 37 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Table 2: Comparison of nifedipine and the long t1/2 calcium channel blocker amlodipine. Daily dose size reflects the potency (IC50) and CLint,u of the drugs. Dose frequency reflects the t1/2.

IC50 (nM)

Cav,u (nM)

t1/2 (h)

CLint, u

Daily Dose

(ml/min/kg)

(mg)

Nifedipine

4

6

2

175

10 TID

Amlodipine

2

2

40

85

5-10 QD

37 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Graph

Just Right

Half Life Too

Too Short

Long Structural Modification

38 ACS Paragon Plus Environment

Page 38 of 47

Page 39 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Figure 1. Basic characteristics of a one compartment PK model. Panel A illustrates the model structure shown with first order absorption and elimination rates. Panels b and c illustrate the characteristic exponential decline following IV bolus administrations. Panels d and e illustrate the absorptive phase followed by monoexponential decline that is typical following oral administration.

39 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

40 ACS Paragon Plus Environment

Page 40 of 47

Page 41 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

41 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 47

Figure 2. Semi-generalized relationships of the extent of steady-state accumulation (A) and peak-to-trough concentration ratio (B) (plotted a log P/Tss values) to half-life-to-dose interval ratio (t1/2 / Tau). Relationships were simulated using the corresponding R(t) and P/Tss equations in Table 1 and assuming an absorption rate of 1.8/hr. A: Accumulation (as calculated as the ratio of AUC at steady-state to AUC following a single dose) increases as the as t1/2 becomes a smaller fraction of the target dose interval. B: While the steady-state average concentration is constant, P/Tss increases exponentially as t1/2 becomes a smaller fraction of the target dose interval.

42 ACS Paragon Plus Environment

Page 43 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Figure 3. Semi-generalized relationship between average pharmacodynamics effect (expressed as a fraction of maximal) and half-life-to-dose interval ratio (t1/2 / Tau). Relationship was simulated using the C(t)ss equation of table 1 assuming a ka of 1.8/hour. Simulations depict 4 dosing levels designed to deliver Cavg,ss values equivalent to arbitrarily selected EC20, EC50, EC80 and EC90 values. Actual average effect was calculated as the dose-interval normalized AUC under an effect curve defined by (100 x C(t)ss) / (C(t)ss + EC50). Volume was varied in order to affect t1/2 / Tau while maintaining a constant Cavg,ss.

43 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

Figure 4. The plasma concentration versus time curve for a compound following escalating oral administration exhibiting dose-dependent (non-linear) PK due to saturation of clearance with increasing dose.

Plasma Concentration (ng/ml)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 47

10000 1000 100

10mg/kg 30mg/kg

10

100mg/kg

1 0

20 40 Time (hours)

44 ACS Paragon Plus Environment

60

Page 45 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Figure 5. Basic characteristics of a two compartment PK model. Panel A illustrates the model structure shown with first-order rates of absorption, elimination and distribution. Panel b illustrates the characteristic biexponential decline in plasma concentrations following IV bolus administration.

45 ACS Paragon Plus Environment

Journal of Medicinal Chemistry

Figure 6. The PK of UK-224,671 in human following single intravenous administration of 0.1 mg/kg. Panel a: Overall PK profile. Panel b: Isolated initial phase profile. Panel c: Isolated

1000

Blood Concentration (ng/ml)

Blood Concentration (ng/ml)

terminal phase profile.

Initial Phase 100

a

Terminal Phase

10 1

1000 0

20

40 60 Time (hours)

80

b

100 10

Initial Phase AUC = 676ng.h/ml 1 0

Blood Concentration (ng/ml)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 47

20

40 60 Time (hours)

1000

80

c

100

Terminal Phase AUC = 171ng.h/ml

10 1 0

20

40 60 Time (hours)

46 ACS Paragon Plus Environment

80

Page 47 of 47

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Medicinal Chemistry

Figure 7. Structure of benoxaprofen and acylglucuronidation metabolite. Oxidative metabolism is likely to be very low due the chlorine substituent and oxidazole ring leading to acylglucuronidation as the major clearance pathway32.

47 ACS Paragon Plus Environment