Target 2035 - an update on private sector contributions.

Target 2035, an international federation of biomedical scientists from the public and private sectors, is leveraging 'open' principles to develop a pharmacological tool for every human protein. These tools are important reagents for scientists studying human health and disease and will facilitate the development of new medicines. It is therefore not surprising that pharmaceutical companies are joining Target 2035, contributing both knowledge and reagents to study novel proteins. Here, we present a brief progress update on Target 2035 and highlight some of industry's contributions.

Permeability of small molecules through a lipid bilayer: a multiscale simulation study.

The transmembrane permeation of eight small (molecular weight <100) organic molecules across a phospholipid bilayer is investigated by multiscale molecular dynamics simulation. The bilayer and hydrating water are represented by simplified, efficient coarse-grain models, whereas the permeating molecules are described by a standard atomic-level force-field. Permeability properties are obtained through a refined version of the z-constraint algorithm. By constraining each permeant at selected depths inside the bilayer, we have sampled free energy differences and diffusion coefficients across the membrane. These data have been combined, according to the inhomogeneous solubility-diffusion model, to yield the permeability coefficients. The results are generally consistent with previous atomic-level calculations and available experimental data. Computationally, our multiscale approach proves 2 orders of magnitude faster than traditional atomic-level methods.

Cardio-vascular safety beyond hERG: in silico modelling of a guinea pig right atrium assay.

As chemists can easily produce large numbers of new potential drug candidates, there is growing demand for high capacity models that can help in driving the chemistry towards efficacious and safe candidates before progressing towards more complex models. Traditionally, the cardiovascular (CV) safety domain plays an important role in this process, as many preclinical CV biomarkers seem to have high prognostic value for the clinical outcome. Throughout the industry, traditional ion channel binding data are generated to drive the early selection process. Although this assay can generate data at high capacity, it has the disadvantage of producing high numbers of false negatives. Therefore, our company applies the isolated guinea pig right atrium (GPRA) assay early-on in discovery. This functional multi-channel/multi-receptor model seems much more predictive in identifying potential CV liabilities. Unfortunately however, its capacity is limited, and there is no room for full automation. We assessed the correlation between ion channel binding and the GPRA's Rate of Contraction (RC), Contractile Force (CF), and effective refractory frequency (ERF) measures assay using over six thousand different data points. Furthermore, the existing experimental knowledge base was used to develop a set of in silico classification models attempting to mimic the GPRA inhibitory activity. The Naïve Bayesian classifier was used to built several models, using the ion channel binding data or in silico computed properties and structural fingerprints as descriptors. The models were validated on an independent and diverse test set of 200 reference compounds. Performances were assessed on the bases of their overall accuracy, sensitivity and specificity in detecting both active and inactive molecules. Our data show that all in silico models are highly predictive of actual GPRA data, at a level equivalent or superior to the ion channel binding assays. Furthermore, the models were interpreted in terms of the descriptors used to highlight the undesirable areas in the explored chemical space, specifically regions of low polarity, high lipophilicity and high molecular weight. In conclusion, we developed a predictive in silico model of a complex physiological assay based on a large and high quality set of experimental data. This model allows high throughput in silico safety screening based on chemical structure within a given chemical space.

A quantitative coarse-grain model for lipid bilayers.

A simplified particle-based computer model for hydrated phospholipid bilayers has been developed and applied to quantitatively predict the major physical features of fluid-phase biomembranes. Compared with available coarse-grain methods, three novel aspects are introduced. First, the main electrostatic features of the system are incorporated explicitly via charges and dipoles. Second, water is accurately (yet efficiently) described, on an individual level, by the soft sticky dipole model. Third, hydrocarbon tails are modeled using the anisotropic Gay-Berne potential. Simulations are conducted by rigid-body molecular dynamics. Our technique proves 2 orders of magnitude less demanding of computational resources than traditional atomic-level methodology. Self-assembled bilayers quantitatively reproduce experimental observables such as electron density, compressibility moduli, dipole potential, lipid diffusion, and water permeability. The lateral pressure profile has been calculated, along with the elastic curvature constants of the Helfrich expression for the membrane bending energy; results are consistent with experimental estimates and atomic-level simulation data. Several of the results presented have been obtained for the first time using a coarse-grain method. Our model is also directly compatible with atomic-level force fields, allowing mixed systems to be simulated in a multiscale fashion.

Structural studies of Salmonella typhimurium ArnB (PmrH) aminotransferase: a 4-amino-4-deoxy-L-arabinose lipopolysaccharide-modifying enzyme.

Lipid A modification with 4-amino-4-deoxy-L-arabinose confers on certain pathogenic bacteria, such as Salmonella, resistance to cationic antimicrobial peptides, including those derived from the innate immune system. ArnB catalysis of amino group transfer from glutamic acid to the 4"-position of a UDP-linked ketopyranose molecule to form UDP-4-amino-4-deoxy-L-arabinose represents a key step in the lipid A modification pathway. Structural and functional studies of the ArnB aminotransferase were undertaken by combining X-ray crystallography with biochemical analyses. High-resolution crystal structures were solved for two native forms and one covalently inhibited form of S. typhimurium ArnB. These structures permitted identification of key residues involved in substrate binding and catalysis, including a rarely observed nonprolyl cis peptide bond in the active site.

Structure-Based Design of a Potent, Selective, and Brain Penetrating PDE2 Inhibitor with Demonstrated Target Engagement.

Structure-guided design led to the identification of the novel, potent, and selective phosphodiesterase 2 (PDE2) inhibitor 12. Compound 12 demonstrated a >210-fold selectivity versus PDE10 and PDE11 and was inactive against all other PDE family members up to 10 µM. In vivo evaluation of 12 provided evidence that it is able to engage the target and to increase cGMP levels in relevant brain regions. Hence, 12 is a valuable tool compound for the better understanding of the role of PDE2 in cognitive impairment and other central nervous system related disorders.

Discovery of a potent, selective, and orally active phosphodiesterase 10A inhibitor for the potential treatment of schizophrenia.

We report the discovery of a series of imidazo[1,2-a]pyrazine derivatives as novel inhibitors of phosphodiesterase 10A (PDE10A). In a high-throughput screening campaign we identified the imidazopyrazine derivative 1, a PDE10A inhibitor with limited selectivity versus the other phosphodiesterases (PDEs). Subsequent investigation of 1 and replacement of the trimethoxyphenyl group by a (methoxyethyl)pyrazole moiety maintained PDE10A inhibition but enhanced selectivity against the other PDEs. Systematic examination and analysis of structure-activity and structure-property relationships resulted in the discovery of 2, an in vitro potent and selective inhibitor of PDE10A with high striatal occupancy of PDE10A, promising in vivo efficacy in different rodent behavioral models of schizophrenia, and a good pharmacokinetic profile in rats.

Towards a better prediction of peak concentration, volume of distribution and half-life after oral drug administration in man, using allometry.

BACKGROUND: It is imperative that new drugs demonstrate adequate pharmacokinetic properties, allowing an optimal safety margin and convenient dosing regimens in clinical practice, which then lead to better patient compliance. Such pharmacokinetic properties include suitable peak (maximum) plasma drug concentration (C(max)), area under the plasma concentration-time curve (AUC) and a suitable half-life (t(½)). The C(max) and t(½) following oral drug administration are functions of the oral clearance (CL/F) and apparent volume of distribution during the terminal phase by the oral route (V(z)/F), each of which may be predicted and combined to estimate C(max) and t(½). Allometric scaling is a widely used methodology in the pharmaceutical industry to predict human pharmacokinetic parameters such as clearance and volume of distribution. In our previous published work, we have evaluated the use of allometry for prediction of CL/F and AUC. In this paper we describe the evaluation of different allometric scaling approaches for the prediction of C(max), V(z)/F and t(½) after oral drug administration in man. METHODS: Twenty-nine compounds developed at Janssen Research and Development (a division of Janssen Pharmaceutica NV), covering a wide range of physicochemical and pharmacokinetic properties, were selected. The C(max) following oral dosing of a compound was predicted using (i) simple allometry alone; (ii) simple allometry along with correction factors such as plasma protein binding (PPB), maximum life-span potential or brain weight (reverse rule of exponents, unbound C(max) approach); and (iii) an indirect approach using allometrically predicted CL/F and V(z)/F and absorption rate constant (k(a)). The k(a) was estimated from (i) in vivo pharmacokinetic experiments in preclinical species; and (ii) predicted effective permeability in man (P(eff)), using a Caco-2 permeability assay. The V(z)/F was predicted using allometric scaling with or without PPB correction. The t(½) was estimated from the allometrically predicted parameters CL/F and V(z)/F. Predictions were deemed adequate when errors were within a 2-fold range. RESULTS: C(max) and t(½) could be predicted within a 2-fold error range for 59% and 66% of the tested compounds, respectively, using allometrically predicted CL/F and V(z)/F. The best predictions for C(max) were obtained when k(a) values were calculated from the Caco-2 permeability assay. The V(z)/F was predicted within a 2-fold error range for 72% of compounds when PPB correction was applied as the correction factor for scaling. CONCLUSIONS: We conclude that (i) C(max) and t(½) are best predicted by indirect scaling approaches (using allometrically predicted CL/F and V(z)/F and accounting for k(a) derived from permeability assay); and (ii) the PPB is an important correction factor for the prediction of V(z)/F by using allometric scaling. Furthermore, additional work is warranted to understand the mechanisms governing the processes underlying determination of C(max) so that the empirical approaches can be fine-tuned further.

Application of PAMPA-models to predict BBB permeability including efflux ratio, plasma protein binding and physicochemical parameters.

This study examines whether algorithms to predict brain penetration of 88 drug candidates could benefit from inclusion of PAMPA data such as P(eff), flux and membrane retention. Specifically the ability to fit experimentally derived LogBB data with PAMPA information and compound related physicochemical and structural parameters was assessed. Collected data were analyzed by partial least square analysis and various regression models for LogBB. Four PAMPA methodologies were evaluated in this study including: (1) a PAMPA-BLM (black lipid membrane) model, (2) a PAMPA-DS (double sink) model, (3) a PAMPA-BBB (blood-brain barrier) model and (4) a PAMPA-BBB-UWL (unstirred water layer). Additionally, plasma protein binding (PPB) experiments and a Caco-2 assay were performed to determine the unbound fraction in plasma and the efflux ratio, respectively, for subsets of the selected compounds. This information was combined with the obtained PAMPA data in an effort to improve the predictions of LogBB. Taken in aggregate, the results presented, suggest that the PAMPA-BLM parameters are the most important contributors to predict the LogBB. The optimized multiple linear regression (MLR) relationship including the PAMPA-BLM properties demonstrated a slightly improved prediction compared to the model without the PAMPA-BLM parameters. Including the plasma protein binding of 15 compounds resulted in a significantly improved PAMPA-BLM prediction of LogBB, while integrating the efflux ratio with PAMPA-BLM or PAMPA-BBB P(eff) values, resulted in improved classification of brain permeable [BBB+(LogBB>/=0)] and impermeable [BBB-(LogBB<0)] compounds.

Application of PAMPA-models to predict BBB permeability including efflux ratio, plasma protein binding and physicochemical parameters.

This study examines whether algorithms to predict brain penetration of 88 drug candidates could benefit from inclusion of PAMPA data such as Peff, flux and membrane retention. Specifically the ability to fit experimentally derived LogBB data with PAMPA information and compound related physicochemical and structural parameters was assessed. Collected data were analyzed by partial least square analysis and various regression models for LogBB. Four PAMPA methodologies were evaluated in this study including: (1) a PAMPA-BLM (black lipid membrane) model, (2) a PAMPA-DS (double sink) model, (3) a PAMPA-BBB (blood-brain barrier) model and (4) a PAMPA-BBB-UWL (unstirred water layer). Additionally, plasma protein binding (PPB) experiments and a Caco-2 assay were performed to determine the unbound fraction in plasma and the efflux ratio, respectively, for subsets of the selected compounds. This information was combined with the obtained PAMPA data in an effort to improve the predictions of LogBB. Taken in aggregate, the results presented, suggest that the PAMPA-BLM parameters are the most important contributors to predict the LogBB. The optimized multiple linear regression (MLR) relationship including the PAMPA-BLM properties demonstrated a slightly improved prediction compared to the model without the PAMPA-BLM parameters. Including the plasma protein binding of 15 compounds resulted in a significantly improved PAMPA-BLM prediction of LogBB, while integrating the efflux ratio with PAMPA-BLM or PAMPA-BBB Peff values, resulted in improved classification of brain permeable [BBB + (LogBB >or= 0)] and impermeable [BBB--(LogBB < 0)] compounds.