Simulation of the Impact of Atropisomer Interconversion on Plasma Exposure of Atropisomers of an Endothelin Receptor Antagonist

PBPK Modeling of PAXLOVIDTM: Incorporating Rotamer Conversion Kinetics to Advanced Dissolution and Absorption Model

PAXLOVID™ is a combination medicine of nirmatrelvir tablets co-packaged with ritonavir tablets. Nirmatrelvir is a peptidomimetic inhibitor of SARS-CoV2 main protease (Mpro), developed for the treatment of COVID-19. Ritonavir is co-administered as a pharmacokinetics (PK) enhancer to inhibit CYP3A mediated metabolism increasing exposures of nirmatrelvir. In the solid form, nirmatrelvir exists in a stable single conformational state (ANTI form). However, nirmatrelvir exhibits atropisomerism in solution whereby upon dissolution the ANTI rotational isomer reversibly converts to another conformation state (SYN form). Nirmatrelvir rotamer conversion follows pseudo first order kinetics with a conversion half-life of approximately 15 min in aqueous solutions, which is on a similar time scale of diffusion mediated dissolution from the solid form. In vitro dissolution studies further indicated that rotamer conversion is one of the processes controlling nirmatrelvir dissolution. It was hypothesized that rotamer conversion kinetics would affect oral absorption of nirmatrelvir in vivo. Consequently, a physiologically based pharmacokinetic (PBPK) model for Paxlovid was developed in Simcyp™ using the advanced dissolution, absorption, and metabolism model (ADAM) by incorporating rotamer conversion kinetics to achieve a more mechanistic description of nirmatrelvir oral absorption. The results demonstrate that the established absorption model with rotamer kinetics adequately described observed clinical data from various nirmatrelvir doses, dosage forms, and dosing regimens. The predicted vs. observed AUCinf and Cmax ratios were within 2-fold. The model has been internally used to inform clinical studies and dose recommendations for pediatrics.

Review and Modeling of Crystal Growth of Atropisomers from Solutions

In this paper, theories on anisotropic crystal growth and crystallization of atropisomers are reviewed and a model for anisotropic crystal growth from solution containing slow inter-converting conformers is presented. The model applies to systems with growth-dominated crystallization from solutions and assumes that only one conformation participates in the solute integration step and is present in the crystal lattice. Other conformers, defined as the wrong conformers, must convert to the right conformer before they can assemble to the crystal lattice. The model presents a simple implicit method for evaluating the growth inhibition effect by the wrong conformers. The crystal growth model applies to anisotropic growth in two main directions, namely a slow-growing face and a fast-growing face and requires the knowledge of solute crystal face integration coefficients in both directions. A parameter estimation algorithm was derived to extract those coefficients from data about temporal concentration and crystal size during crystallization and was designed to have a short run time, while providing a high-resolution estimation. The model predicts a size-dependent growth rate and simulations indicated that for a given seed size and solvent system and for an isothermal anti-solvent addition crystallization, the seed loading and the supersaturation at seeding are the main factors impacting the final aspect ratio. The model predicts a decrease of the growth inhibition effect by the wrong conformer with increasing temperature, likely due to faster equilibration between conformers and/or a decrease of the population of the wrong conformer, if of low energy, at elevated temperatures. Finally, the model predicts that solute surface integration becomes the rate-limiting mechanism for high solute integration activation energies, resulting in no impact of the WC on the overall crystal growth process.

Characterization of Stereoselective Metabolism, Inhibitory Effect on Uric Acid Uptake Transporters, and Pharmacokinetics of Lesinurad Atropisomers

Lesinurad [Zurampic; 2-(5-bromo-4-(4-cyclopropylnaphthalen-1-yl)-4H-1,2,4-triazol-3-ylthio)], a selective inhibitor of uric acid reabsorption transporters approved for the treatment of gout, is a racemate of two atropisomers. The objective of this investigation was to evaluate the stereoselectivity of metabolism, the inhibitory potency on kidney uric acid reabsorption transporters (URAT1 and OAT4), and the clinical pharmacokinetics of the lesinurad atropisomers. Incubations with human liver microsomes (HLM), recombinant CYP2C9, and recombinant CYP3A4 were carried out to characterize the stereoselective formation of three metabolites: M3 (hydroxylation), M4 (a dihydrodiol metabolite), and M6 (S-dealkylation). The formation of M3 in HLM with atropisomer 1 was approximately twice as much as that with atropisomer 2, whereas formation of M4 with atropisomer 1 was 8- to 12-fold greater than that with atropisomer 2. There were no significant differences in the plasma protein binding among lesinurad and the atropisomers. Following oral administration of 400 mg lesinurad once daily for 14 days to healthy human volunteers, the systemic exposure (C _max at steady state and area under the concentration-time curve from time zero to the time of dosing interval) of atropisomer 1 was approximately 30% lower than that of atropisomer 2, whereas renal clearance was similar. In vitro cell-based assays using HEK293 stable cells expressing URAT1 and OAT4 demonstrated that atropisomer 2 was approximately 4-fold more potent against URAT1 than atropisomer 1 and equally active against OAT4. In conclusion, lesinurad atropisomers showed stereoselectivity in clinical pharmacokinetics, metabolism, and inhibitory potency against URAT1.

Biological and Structural Characterization of Rotamers of C-C Chemokine Receptor Type 5 (CCR5) Inhibitor GSK214096

We recently reported the discovery of preclinical CCR5 inhibitor GSK214096, 1 (J. Med. Chem. 2011, 54, 756). Detailed characterization of 1 revealed that it exists as a mixture of four separable atropisomers A-D. The two slow-interconverting pairs of rotamers A + B and C + D were separated and further characterized. HIV and CCR5-mediated chemotaxis data strongly suggest that the antiviral potency of 1 is due to rotamers A + B and not C + D. Furthermore, integrated UV, vibrational circular dichroism VCD and computational approach allowed to determine the M chirality in C + D (and P chirality in A + B). These findings imply additional avenues to be pursued toward new CCR5 antagonists.

Vibrational circular dichroism (VCD) chiral assignment of atropisomers: application to γ-aminobutyric acid (GABA) modulators designed as potential anxiolytic drugs

Atropisomers exist when axial chirality is present as a result of conformationally restricted rotation around a single bond. The interconversion rate of the individual atropisomers is critical to the assessment of chiral stability of a drug throughout scale-up, development, production, and storage as well as in vivo pharmacokinetics. We describe the application of vibrational circular dichroism spectroscopy coupled with quantum mechanics simulations to assign the absolute axial chirality and measure the racemization half-life of a series of potential anxiolytic drugs that act as γ-aminobutyric acid modulators.

Atropisomeric property of 1-(2,6-difluorobenzyl)-3-[(2R)-amino-2-phenethyl]-5-(2-fluoro-3-methoxyphenyl)-6-methyluracil

1-(2,6-Difluorobenzyl)-3-[(2R)-amino-2-phenethyl]-5-(2-fluoro-3-methoxyphenyl)-6-methyluracil (6), a potent and orally active antagonist of the human gonadotropin-releasing hormone receptor, exists as a pair of atropisomers in solution, which was detected by NMR spectroscopy, and separable by HPLC. In addition to a (R)-configured benzylamine, there is a second stereogenic element due to the presence of a chiral axis between the substituted 5-phenyl group and the uracil core. The rate constant of the interconversion (k = 5.07 x 10(-5) s(-1)) of these two atropisomers was determined by proton NMR analysis of a diastereoisomer-enriched sample in aqueous solution at 25 degrees C, and the corresponding Gibbs free energy DeltaG(#) of rotation barrier (97.4 kJ mol(-1)) was calculated using the Eyring equation. The diastereoisomer half-life at physiological temperature (37 degrees C) in aqueous media was estimated to be about 46 min.