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.

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.

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.

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.

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.

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.

Maximizing the success of bile duct cannulation studies in rats: recommendations for best practice.

Bile duct cannulation (BDC) studies are usually carried out in the rat to support the absorption, distribution, metabolism and excretion profiling of novel agrochemicals and pharmaceuticals. The different aspects of these studies (e.g. surgical preparation, dosing and collection of bile) can be intricate and/or technically complex. The animals are often kept singly housed for the duration of the studies following surgical implantation of the cannulas. The generation of insufficient data to meet the study objectives, for example due to failure in cannula patency, can result in the need to repeat these studies. A working group of contract research organizations that routinely carry out BDC studies was brought together by the National Centre for the Replacement, Refinement and Reduction of Animals in Research (NC3Rs) to share their experiences, and to establish the key factors necessary to ensure routinely high success rates. Through these discussions the group has identified opportunities for best practice across various aspects of the studies. The aim of these recommendations is to support all staff involved in conducting BDC studies to maximize the amount of useful data generated using the fewest animals possible, while ensuring the highest possible standards of animal welfare.

Minipig and Human Metabolism of Aldehyde Oxidase Substrates: In Vitro-In Vivo Comparisons.

The importance of aldehyde oxidase (AOX) is becoming increasingly recognized in the prediction of human pharmacokinetic parameters from animal data. The objectives of these studies were to ascertain whether an in vitro-in vivo correlation existed in the clearance and metabolic pathways of AOX substrates and to establish whether the minipig represented an appropriate non-rodent model for man in the pre-clinical development of drugs metabolized by AOX. Using the AOX substrates carbazeran, 6-deoxypenciclovir and zaleplon, clearance was estimated from in vitro depletion experiments with minipig and human liver cytosol and microsomes and scaled before comparison with data generated in parallel in vivo studies in minipigs. In vitro and in vivo metabolic pathways were characterized by LC-MS/MS. Scaling of in vitro metabolism data to predict in vivo clearance underestimated in vivo values, although the rank order of clearance for the three compounds was preserved. Prediction of human in vivo clearance from scaled minipig in vivo data produced results which correlated well with published clinical values. Overall, this study is the first to compare minipig in vitro metabolism data with in vivo pharmacokinetic data for compounds metabolized by AOX and provides a scientific rationale for the selection of this species as a model for humans in the development of drugs which are substrates of AOX.

Report on the sixth blind test of organic crystal structure prediction methods.

The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations. All of the targets, apart from a single potentially disordered Z' = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.

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.

A mechanophysical phase transition provides a dramatic example of colour polymorphism: the tribochromism of a substituted tri(methylene)tetrahydrofuran-2-one.

BACKGROUND: Derivatives of fulgides have been shown to have interesting photochromic properties. We have synthesised a number of such derivatives and have found, in some cases, that crystals can be made to change colour on crushing, a phenomenon we have termed "tribochromism". We have studied a number of derivatives by X-ray crystallography, to see if the colour is linked to molecular structure or crystal packing, or both, and our structural results have been supported by calculation of molecular and lattice energies. RESULTS: A number of 5-dicyanomethylene-4-diphenylmethylene-3-disubstitutedmethylene-tetrahydrofuran-2-one compounds have been prepared and structurally characterised. The compounds are obtained as yellow or dark red crystals, or, in one case, both. In two cases where yellow crystals were obtained, we found that crushing the crystals gave a deep red powder. Structure determinations, including those of the one compound which gave both coloured forms, depending on crystallisation conditions, showed that the yellow crystals contained molecules in which the structure comprised a folded conformation at the diphenylmethylene site, whilst the red crystals contained molecules in a twisted conformation at this site. Lattice energy and molecular conformation energies were calculated for all molecules, and showed that the conformational energy of the molecule in structure IIIa (yellow) is marginally higher, and the conformation thus less stable, than that of the molecule in structure IIIb (red). However, the van der Waals energy for crystal structure IIIa, is slightly stronger than that of structure IIIb - which may be viewed as a hint of a metastable packing preference for IIIa, overcome by the contribution of a more stabilising Coulomb energy to the overall more favourable lattice energy of structure IIIb. CONCLUSIONS: Our studies have shown that the crystal colour is correlated with one of two molecular conformations which are different in energy, but that the less stable conformation can be stabilised by its host crystal lattice. Graphical abstractGraphical representation of the structural and colour change in the tribochromic compound (III).

Co-crystallization with nicotinamide in two conformations lowers energy but expands volume.

A dispersion-corrected density functional theory method has been used to study the formation energies and volumes of cocrystals. For four cocrystals of nicotinamide (NIC) and (R)-mandelic acid, a broad agreement is found between experimental and computed values. We report that cocrystals containing NIC are anomalous as their formation generally decreases energy but expands volume. In this respect, the formation of NIC cocrystals is in contrast to most physical processes, but similar to water freezing. As in the case of water freezing, the cocrystallization with NIC leads to stronger hydrogen bonds and looser molecular packing, a combination that is likely responsible for the anticorrelation between energy and volume. NIC has two conformers 4 kJ/mol apart in energy and both can form cocrystals, with the resulting structures having comparable formation energies and volumes. These results are relevant for understanding the nature of cocrystallization and why NIC is a prolific cocrystal former.

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.

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.

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.

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.

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.