One small step in time, one giant leap for DMPK kind - a CRO perspective of the evolving core discipline of drug development.

As the Space Race or Formula 1 drives innovation, efficiency and progress in home technology and home car markets, Drug Metabolism and Pharmacokinetics (DMPK) drives scientific innovation and value for drug development companies. Stand still and fall behind as the saying goes, and these analogies are true as much in the design and conduct of DMPK studies as they are in the technology and manufacturing sectors.This short review showcases the impact that DMPK has had on drug development and how it has changed in the last 10 years, illustrating the value added scientific benefit, cost and time saving, that innovative DMPK program design and study conduct have. Examples and case studies spanning novel in vitro alternatives such as organ-on-a-chip (OOAC) developments; use of in vivo microsampling across small and large animal species; challenging historical paradigms in Absorption, Distribution, Metabolism and Excretion (ADME) studies; and embracing new technologies to address regulatory concerns, are presented.The continual pace of change has kept DMPK at the core of pharmaceutical, crop and chemical evaluation, and this is set to continue as regulators use this discipline to inform decision-making. With new modalities and new scientific questions, DMPK will continue to evolve, with the likes of new in vitro, in vivo and in silico models becoming central to candidate selection and progression.

An OECD TG 428 study ring trial with 14C-Caffeine demonstrating repeatability and robustness of the dermal absorption in vitro method.

The dermal absorption potential of 14C-Caffeine applied as a 4 mg/mL concentration (10 µL/cm2 finite dose) was investigated in six laboratories under Good Laboratory Practice conditions using an OECD TG 428-compliant in vitro assay with flow-through cells and split-thickness human skin. Potential sources of variation were reduced by a standardized protocol, test item and skin source. Particularly, skin samples from same donors were distributed over two repeats and between labs in a non-random, stratified design. Very similar recovery was achieved in the various assay compartments between laboratories, repeats and donors, demonstrating that the assay can be robustly and reliably performed. The absorption in one laboratory was 5-fold higher than in the others. This did not clearly correlate with skin integrity parameters but might be associated with an accidental COVID-19 pandemic-related interruption in sample shipment. It is possible that other factors may affect dermal absorption variation not routinely assessed or considered in the current method. The mean receptor fluid recovery, potential absorption (recovery in receptor fluid and skin except tape strips 1 and 2) and mass balance of caffeine was 6.99%, 7.14% and 99.13%, respectively, across all and 3.87%, 3.96% and 99.00% in the subset of five laboratories.

PDielec: The calculation of infrared and terahertz absorption for powdered crystals.

The Python package PDielec is described, which calculates the infrared absorption characteristics of a crystalline material supported in a non-absorbing medium. PDielec post processes solid-state quantum mechanical and molecular mechanical calculations of the phonons and dielectric response of the crystalline material. Using an effective medium method, the package calculates the internal electric field arising from different particle morphologies and calculates the resulting shift in absorption frequency and intensity arising from the coupling between a phonon and the internal field. The theory of the approach is described, followed by a description of the implementation within PDielec. Finally, a section providing several examples of its application is given. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Effect of molecular size and particle shape on the terahertz absorption of a homologous series of tetraalkylammonium salts.

The absorption coefficient and refractive index have been measured for a homologous series of tetraalkylammonium bromides over the frequency range 0.3-5.5 THz. Spectral features are found to shift to lower frequencies as the molecular mass is increased, as expected. However, to understand the detailed structure of the observed spectral features, density functional perturbation theory calculations have been performed on the first four crystalline compounds in the series. From these calculations, we find that each spectrum is dominated by three translatory modes involving asymmetric motion of the ammonium cation and bromine counterion, although the overall number of active modes increases with increasing molecular size. The experimentally observed absorption is not completely described by the infrared active phonon modes alone. We show that it is also necessary to include the coupling of the phonon modes with the macroscopic field generated by the collective displacement of the vibrating ions, and we have applied an effective medium theory, which accounts for particle shape to allow for this effect in the calculation of the terahertz spectra.

Exploring the Stability and Disorder in the Polymorphs of L-Cysteine through Density Functional Theory and Vibrational Spectroscopy.

Static and dynamic density functional calculations are reported for the four known polymorphs of l-cysteine. Static calculations are used to explore the relative free energies (within the harmonic approximation) of the polymorphs as a function of pressure. An important feature of the structural differences between the polymorphs is shown to be the dihedral angle of the C-C-S-H bond. It is shown that, by varying this angle, it is possible to move between hydrogen bonding motifs S-H···S and S-H···O in all four polymorphs. The energetics for dihedral angle rotation are explored, and the barriers for rotation between the hydrogen bonding motifs have been calculated for each polymorph. Two possible models for the experimental disorder observed in Form I at room temperature are explored using both static and dynamic methods; a domain disorder model, where the disorder is localized, and a dispersed disorder model, where the disorder is randomly distributed throughout the crystal. Molecular dynamics calculations show transitions between the two hydrogen bonding motifs occurring in the dispersed disorder model at 300 and 350 K. In addition, molecular dynamics calculations of Form IV also showed the onset of hydrogen bond disorder at 300 K. Calculations of the predicted infrared and terahertz absorption are performed for both the static and dynamic simulations, and the results are compared with experimental results to understand the influence of disorder on the observed spectra.

The Changing Landscape for Human Absorption, Metabolism, and Excretion: Practical Experiences From a Data Analysis of 500 Studies.

Coincidental with the intensified regulatory and industry focus on the design and conduct of human absorption, metabolism, and excretion (hAME) studies in the past 12 months, we have recently completed our 500th cohort involving radiolabeled test item administration to humans. Here, we build upon a recent industry white paper in this journal1 and share some of our own experiences as a Contract Research Organization based upon collaborations with numerous pharma companies and their differing approaches to design and timing, to add further context to the discussion regarding hAME studies and the pivotal role that drug metabolism and pharmacokinetics plays. In this article, we explore how both changing relationships within the industry and shifting regulatory guidelines are impacting strategies, and compare EU and US pre-study approval requirements, before evaluating the trends from over 500 studies conducted at our global facilities conducted over more than 30 years. We conclude with a review of how improved technical capabilities and strategies are influencing the design and conduct of hAME studies, before speculating on some of the driving factors which may shape the direction they take in the future.

Empirical van der Waals corrections to solid-state density functional theory: Iodine and phosphorous containing molecular crystals.

Parameters are derived for a molecular mechanics type dispersive correction to solid-state density functional theory calculations on molecular crystals containing iodine and phosphorous. The molecular C(6) coefficients are derived from photoabsorption differential oscillator strength spectra determined from accurate (e,e) dipole spectra. The cross-over parameters, which ensure correct behavior at short internuclear distances, are obtained by fitting predicted crystal lattice parameters to experimental data. The accuracy of the parameterization is assessed by optimizing the experimental structures of several additional phosphorous and iodine containing molecular crystals and by examining the relative stabilities of the known polymorphs of phosphorous pentoxide and the stabilities of different packings of an iodine containing molecular crystal, 2,9-bis(iodo)anthanthrone, which has been the subject of a crystal structure prediction study. Optimizations of the experimental crystal structures did not lead to significant geometric deviations. The optimized experimental structure of 2,9-bis(iodo)anthanthrone is the lowest energy packing found, indicating a satisfactory description of both energy and structure for these molecular crystals.

Progress in crystal structure prediction.

The results of the application of a density functional theory method incorporating dispersive corrections in the 2010 crystal structure prediction blind test are reported. The method correctly predicted four out of the six experimental structures. Three of the four correct predictions were found to have the lowest lattice energy of any crystal structure for that molecule. The experimental crystal structures for all six compounds were found during the structure generation phase of the simulations, indicating that the tailor-made force fields used for screening structures were valid and that the structure generation engine, which combines a Monte Carlo parallel tempering algorithm with an efficient lattice energy minimiser, was working effectively. For the three compounds for which the experimental crystal structures did not correspond to the lowest energy structures found, the method for calculating the lattice energy needs to be further refined or there may be other polymorphs that have not yet been found experimentally.

Predictability of the polymorphs of small organic compounds: crystal structure predictions of four benchmark blind test molecules.

Predicting the crystal structure of an organic molecule from first principles has been a major challenge in physical chemistry. Recently, the application of Density Functional Theory including a dispersive energy correction (the DFT(d) method) has been shown to be a reliable method for predicting experimental structures based purely on their ranking according to lattice energy. Further validation results of the application of the DFT(d) method to four organic molecules are presented here. The compounds were targets (labelled molecule II, VI, VII and XI) in previous blind tests of crystal structure prediction, and their structures proved difficult to predict. However, this study shows that the DFT(d) approach is capable of predicting the solid state structures of these small molecules. For molecule VII, the most stable (rank 1) predicted crystal structure corresponds to the experimentally observed structure. For molecule VI, the rank 1, 2 and 3 predicted structures correspond to the three experimental polymorphs, forms I, III and II, respectively. For molecules II and XI, their rank 1 predicted structures are energetically more stable than those corresponding to the experimental crystal structures, and were not found amongst the structures submitted by the participants in the blind tests. The rank 1 structure of molecule II is predicted to exist under high pressure, whilst the rank 1 structure predicted for molecule XI has the same space group and hydrogen bonding pattern as observed in the crystal of 1-amino-1-methyl-cyclopropane, which is structurally related to molecule XI. The experimental crystal structure of molecule II corresponds to the rank 4 prediction, 0.8 kJ mol(-1) above the global minimum structure, and the experimental structure of molecule XI corresponds to the rank 2 prediction, 0.4 kJ mol(-1) above the global minimum.

Towards crystal structure prediction of complex organic compounds--a report on the fifth blind test.

Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories - a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.

A new group housing approach for non-human primate metabolism studies.

Understanding the absorption, distribution, metabolism and excretion (ADME) of candidate drugs in preclinical species is an integral part of the safety and efficacy evaluation in drug development. For this purpose, the housing of single animals in metabolism cages has historically been common practice for ADME studies. Whilst mini-pigs and dogs are selected wherever possible, non-human primates (NHPs) are used where there is no suitable scientific alternative. Having undergone only minimal revisions over the past 30 years, the traditional single-housing metabolism cage design for NHPs significantly limits normal vertical movement and social behaviours in primates. Minimising animal suffering and improving welfare is an important aspect of working with animals in research and Novo Nordisk A/S, together with collaborators, has focused on this area for many years. A novel metabolism cage for group housing of NHPs has been designed in a joint collaboration between Novo Nordisk A/S and Covance Inc. The advantages of this novel cage are extensive, including a significantly increased cage volume and ability for socialisation, as well as improvements to alleviate stress and boredom. The excretion balance data from six male NHPs housed in single or group metabolism cages were compared using the radiolabelled test compound [14C]-quetiapine. Welfare, in terms of stress and behaviour, when animals were single or group housed was also assessed. Mean recoveries of radioactivity were shown to be comparable irrespective of housing design (83.2% for group-housed animals vs. 87.1% for single-housed animals), supporting the potential suitability of NHP group housing for future metabolism ADME studies.

Crystal structure prediction and isostructurality of three small molecules.

A crystal structure prediction (CSP) study of three small, rigid and structurally related organic compounds (differing only in the position and number of methyl groups) is presented. A tailor-made force field (TMFF; a non-transferable force field specific for each molecule) was constructed with the aid of a dispersion-corrected density functional theory method (the hybrid method). Parameters for all energy terms in each TMFF were fitted to reference data generated by the hybrid method. Each force field was then employed during structure generation. The experimentally observed crystal structures of two of the three molecules were found as the most stable crystal packings in the lists of their force-field-optimised structures. A number of the most stable crystal structures were re-optimised with the hybrid method. One experimental crystal structure was still calculated to be the most stable structure, whereas for another compound the experimental structure became the third most stable structure according to the hybrid method. For the third molecule, the experimentally observed polymorph, which was found to be the fourth most stable form using its TMFF, became the second most stable form. Good geometrical agreements were observed between the experimental structures and those calculated by both methods. The average structural deviation achieved by the TMFFs was almost twice that obtained with the hybrid method. The TMFF approach was extended by exploring the accuracy of a more general TMFF (GTMFF), which involved fitting the force-field parameters to the reference data for all three molecules simultaneously. This GTMFF was slightly less accurate than the individual TMFFs but still of sufficient accuracy to be used in CSP. A study of the isostructural relationships between these molecules and their crystal lattices revealed a potential polymorph of one of the compounds that has not been observed experimentally and that may be accessible in a thorough polymorph screen, through seeding, or through the use of a suitable tailor-made additive.

Calculation and measurement of terahertz active normal modes in crystalline PETN.

The terahertz frequency spectrum of pentaerythritol tetranitrate (PETN) is calculated using Discover with the COMPASS force field, CASTEP and PWscf. The calculations are compared to each other and to terahertz spectra (0.3-3 THz) of crystalline PETN recorded at 4 K. A number of analysis methods are used to characterise the calculated normal modes.

Crystal structure prediction and isostructurality of three small organic halogen compounds.

A theoretical investigation of the packing stabilities of three small organic halogen compounds is presented based on a crystal structure prediction (CSP) study. Each compound has four identical halogen atoms (fluorine, chlorine, and bromine) and a four-membered ring consisting of carbon and sulfur atoms arranged alternately. Two halogen atoms are attached to each carbon and two oxygen atoms are attached to each sulfur forming SO(2) functional groups. The crystal structures of these compounds have been determined experimentally and show distinct packing arrangements. Utilising the computational approaches implemented in the GRACE software package, each compound is subjected to a full CSP study using a force field specific for each molecule (called the tailor-made force field or TMFF) and a dispersion corrected solid-state density functional method (or DFT(d) method). Energetically feasible crystal structures are generated in all 230 space groups restricted to a single molecule in the crystallographic asymmetric unit (Z' = 1) using the TMFF of each molecule. Next, a selection of structures with low TMFF lattice energies are further refined with the DFT(d) method. The CSP results show that the experimental crystal structures of the molecules containing fluorine and chlorine are well described energetically and geometrically by their TMFFs and the DFT(d) method. Both approaches locate their experimental lattices as the most stable structures. For the molecule containing bromine, a crystal structure corresponding to the force field optimised experimental structure is located as the second structure in the list of force field predicted structures, ranked by calculated lattice energy. Despite the structural similarity of the predicted and experimental structures, close examination of the DFT(d) optimisation results of the experimental structure reveals a slightly lower energy structure than that found by the CSP simulations. Furthermore, minimisation of the force field optimised structure using the DFT(d) method does not lead to the same minimum as the DFT(d) optimised experimental structure. Based on the CSP results and isostructurality among these three compounds, two new potential polymorphs for the molecules containing chlorine and bromine are proposed. These polymorphs might be obtained experimentally under the right crystallisation conditions.

A novel welfare and scientific approach to conducting dog metabolism studies allowing dogs to be pair housed.

Metabolism cages are designed to conduct absorption, distribution, metabolism and excretion (ADME) studies, enabling an 'excretion balance' scientific objective to be met. Historically, the design of dog metabolism cages has involved single housing. This type of housing has limitations for normal social behaviours and has been largely unchanged for 25-30 years. Improving animal welfare is a focus area for the authorities as well as the industry throughout the European Union. A collaboration was developed between Novo Nordisk and Covance to enhance the design of metabolism cages, allowing dogs to be pair housed. The purpose of the study was to compare excretion balance data from pair-housed and singly housed dogs in order to demonstrate that conducting excretion balance studies with a pair-housing design improves animal welfare without compromising the scientific integrity of the study. A radiolabelled test compound, [14C]-Quetiapine, was selected for this investigation based on its excretion profile. The assessment of the dogs' stress levels was investigated by measuring the levels of serum cortisol as an indicative biomarker. Results were inconclusive due to large variations in cortisol levels. However, dogs appeared calmer in the pair-housing setting. The overall mean recovery (±standard deviation) for pair-housed animals (94.0 ± 0.66% of the dose) was equivalent to that from singly housed dogs (93.0 ± 2.29%). Based on these data, we conclude that pair housing of dogs for future metabolism ADME studies does not compromise the scientific integrity, and therefore is a major progression in the design of these studies, enhancing welfare.

Significant progress in predicting the crystal structures of small organic molecules--a report on the fourth blind test.

We report on the organization and outcome of the fourth blind test of crystal structure prediction, an international collaborative project organized to evaluate the present state in computational methods of predicting the crystal structures of small organic molecules. There were 14 research groups which took part, using a variety of methods to generate and rank the most likely crystal structures for four target systems: three single-component crystal structures and a 1:1 cocrystal. Participants were challenged to predict the crystal structures of the four systems, given only their molecular diagrams, while the recently determined but as-yet unpublished crystal structures were withheld by an independent referee. Three predictions were allowed for each system. The results demonstrate a dramatic improvement in rates of success over previous blind tests; in total, there were 13 successful predictions and, for each of the four targets, at least two groups correctly predicted the observed crystal structure. The successes include one participating group who correctly predicted all four crystal structures as their first ranked choice, albeit at a considerable computational expense. The results reflect important improvements in modelling methods and suggest that, at least for the small and fairly rigid types of molecules included in this blind test, such calculations can be constructively applied to help understand crystallization and polymorphism of organic molecules.

Rapid Ring-Opening Metathesis Polymerization of Monomers Obtained from Biomass-Derived Furfuryl Amines and Maleic Anhydride.

Well-controlled and extremely rapid ring-opening metathesis polymerization of unusual oxanorbornene lactam esters by Grubbs third-generation catalyst is used to prepare a range of bio-based homo- and copolymers. Bio-derived oxanorbornene lactam monomers were prepared at room temperature from maleic anhydride and secondary furfuryl amines by using a 100 % atom economical, tandem Diels-Alder lactamization reaction, followed by esterification. Several of the resulting homo- and copolymers show good control over polymer molecular weight and have narrow molecular weight distributions.

Pharmacokinetic investigation of imatinib using accelerator mass spectrometry in patients with chronic myeloid leukemia.

PURPOSE: To investigate the potential use of accelerator mass spectrometry (AMS) in the study of the clinical pharmacology of imatinib. EXPERIMENTAL DESIGN: Six patients who were receiving imatinib (400 mg/d) as part of their ongoing treatment for chronic myeloid leukemia (CML) received a dose containing a trace quantity (13.6 kBq) of (14)C-imatinib. Blood samples were collected from patients before and at various times up to 72 h after administration of the test dose and were processed to provide samples of plasma and peripheral blood lymphocytes (PBL). Samples were analyzed by AMS, with chromatographic separation of parent compound from metabolites. In addition, plasma samples were analyzed by liquid chromatography/mass spectrometry (LCMS). RESULTS: Analysis of the AMS data indicated that imatinib was rapidly absorbed and could be detected in plasma up to 72 h after administration. Imatinib was also detectable in PBL at 24 h after administration of the (14)C-labeled dose. Comparison of plasma concentrations determined by AMS with those derived by LCMS analysis gave similar average estimates of area under plasma concentration time curve (26 +/- 3 versus 27 +/- 11 microg/mL.h), but with some variation within each individual. CONCLUSIONS: Using this technique, data were obtained in a small number of patients on the pharmacokinetics of a single dose of imatinib in the context of chronic dosing, which could shed light on possible pharmacologic causes of resistance to imatinib in CML.

Modeling the interplay of inter- and intramolecular hydrogen bonding in conformational polymorphs.

The predicted stability differences of the conformational polymorphs of oxalyl dihydrazide and ortho-acetamidobenzamide are unrealistically large when the modeling of intermolecular energies is solely based on the isolated-molecule charge density, neglecting charge density polarization. Ab initio calculated crystal electron densities showed qualitative differences depending on the spatial arrangement of molecules in the lattice with the greatest variations observed for polymorphs that differ in the extent of inter- and intramolecular hydrogen bonding. We show that accounting for induction dramatically alters the calculated stability order of the polymorphs and reduces their predicted stability differences to be in better agreement with experiment. Given the challenges in modeling conformational polymorphs with marked differences in hydrogen bonding geometries, we performed an extensive periodic density functional study with a range of exchange-correlation functionals using both atomic and plane wave basis sets. Although such electronic structure methods model the electrostatic and polarization contributions well, the underestimation of dispersion interactions by current exchange-correlation functionals limits their applicability. The use of an empirical dispersion-corrected density functional method consistently reduces the structural deviations between the experimental and energy minimized crystal structures and achieves plausible stability differences. Thus, we have established which types of models may give worthwhile relative energies for crystal structures and other condensed phases of flexible molecules with intra- and intermolecular hydrogen bonding capabilities, advancing the possibility of simulation studies on polymorphic pharmaceuticals.