Carrier Mediated Distribution System (CAMDIS): a new approach for the measurement of octanol/water distribution coefficients

From mundane to surprising nonadditivity: drivers and impact on ML models

Nonadditivity (NA) in Structure-Activity and Structure-Property Relationship (SAR) data is a rare but very information rich phenomenon. It can indicate conformational flexibility, structural rearrangements, and errors in assay results and structural assignment. While purely ligand-based conformational causes of NA are rather well understood and mundane, other factors are less so and cause surprising NA that has a huge influence on SAR analysis and ML model performance. We here report a systematic analysis across a wide range of properties (20 on-target biological activities and 4 physicochemical ADME-related properties) to understand the frequency of various different phenomena that may lead to NA. A set of novel descriptors were developed to characterize double transformation cycles and identify trends in NA. Double transformation cycles were classified into "surprising" and "mundane" categories, with the majority being classed as mundane. We also examined commonalities among surprising cycles, finding LogP differences to have the most significant impact on NA. A distinct behavior of NA for on-target sets compared to ADME sets was observed. Finally, we show that machine learning models struggle with highly nonadditive data, indicating that a better understanding of NA is an important future research direction.

A highly potent, orally bioavailable pyrazole-derived cannabinoid CB2 receptor-selective full agonist for in vivo studies

The cannabinoid CB2 receptor (CB2R) is a potential therapeutic target for distinct forms of tissue injury and inflammatory diseases. To thoroughly investigate the role of CB2R in pathophysiological conditions and for target validation in vivo, optimal pharmacological tool compounds are essential. Despite the sizable progress in the generation of potent and selective CB2R ligands, pharmacokinetic parameters are often neglected for in vivo studies. Here, we report the generation and characterization of a tetra-substituted pyrazole CB2R full agonist named RNB-61 with high potency (K i 0.13-1.81 nM, depending on species) and a peripherally restricted action due to P-glycoprotein mediated efflux from the brain. 3H and 14C labelled RNB-61 showed apparent K d values < 4 nM towards human CB2R in both cell and tissue experiments. The >6000-fold selectivity over CB1 receptors and negligible off-targets in vitro, combined with high oral bioavailability and suitable systemic pharmacokinetic (PK) properties, prompted the assessment of RNB-61 in a mouse ischemia-reperfusion model of acute kidney injury (AKI) and in a rat model of chronic kidney injury/inflammation and fibrosis (CKI) induced by unilateral ureteral obstruction. RNB-61 exerted dose-dependent nephroprotective and/or antifibrotic effects in the AKI/CKI models. Thus, RNB-61 is an optimal CB2R tool compound for preclinical in vivo studies with superior biophysical and PK properties over generally used CB2R ligands.

Discovery of Orally Available and Brain Penetrant AEP Inhibitors

Alzheimer's Disease (AD) is the most widespread form of dementia, with one of the pathological hallmarks being the formation of neurofibrillary tangles (NFTs). These tangles consist of phosphorylated Tau fragments. Asparagine endopeptidase (AEP) is a key Tau cleaving enzyme that generates aggregation-prone Tau fragments. Inhibition of AEP to reduce the level of toxic Tau fragment formation could represent a promising therapeutic strategy. Here, we report the first orthosteric, selective, orally bioavailable, and brain penetrant inhibitors with an irreversible binding mode. We outline the development of the series starting from reversible molecules and demonstrate the link between inhibition of AEP and reduction of Tau N368 fragment both in vitro and in vivo.

Drug Design and Success of Prospective Mouse In Vitro-In Vivo Extrapolation (IVIVE) for Predictions of Plasma Clearance (CLp) from Hepatocyte Intrinsic Clearance (CLint)

Hepatocyte intrinsic clearance (CL_int) and methods of in vitro-in vivo extrapolation (IVIVE) are often used to predict plasma clearance (CL_p) in drug discovery. While the prediction success of this approach is dependent on the chemotype, specific molecular properties and drug design features that govern these outcomes are poorly understood. To address this challenge, we investigated the success of prospective mouse CL_p IVIVE across 2142 chemically diverse compounds. Dilution scaling, which assumes that the free fraction in hepatocyte incubations (f_u,inc) is governed by binding to the 10% of serum in the incubation medium, was used as our default CL_p IVIVE approach. Results show that predictions of CL_p are better for smaller (molecular weight (MW) < 500 Da), less polar (total polar surface area (TPSA) < 100 Ų, hydrogen bond donor (HBD) ≤1, hydrogen bond acceptor (HBA) ≤ 6), lipophilic (log D > 3), and neutral compounds, with low HBD count playing the key role. If compounds are classified according to their chemical space, predictions were good for compounds resembling central nervous system (CNS) drugs [average absolute fold error (AAFE) of 2.05, average fold error (AFE) of 0.90], moderate for classical druglike compounds (according to Lipinski, Veber, and Ghose guidelines; AAFE of 2.55; AFE of 0.68), and poor for nonclassical "beyond the rule of 5" compounds (AAFE of 3.31; AFE of 0.41). From the perspective of measured druglike properties, predictions of CL_p were better for compounds with moderate-to-high hepatocyte CL_int (>10 μL/min/10⁶ cells), high passive cellular permeability (P_app > 100 nm/s), and moderate observed CL_p (5-50 mL/min/kg). Influences of plasma protein binding (f_u,p) and P-glycoprotein (Pgp) apical efflux ratio (AP-ER) were less pronounced. If the extended clearance classification system (ECCS) is applied, predictions were good for class 2 (P_app > 50 nm/s; neutral or basic; AAFE of 2.35; AFE of 0.70) and acceptable for class 1A compounds (AAFE of 2.98; AFE of 0.70). Classes 1B, 3 A/B, and 4 showed poor outcomes (AAFE > 3.80; AFE < 0.60). Functional groups trending toward weaker CL_p IVIVE were esters, carbamates, sulfonamides, carboxylic acids, ketones, primary and secondary amines, primary alcohols, oxetanes, and compounds liable to aldehyde oxidase metabolism, likely due to multifactorial reasons. Multivariate analysis showed that multiple properties are relevant, combining together to define the overall success of CL_p IVIVE. Our results indicate that the current practice of prospective CL_p IVIVE is suitable only for CNS-like compounds and well-behaved classical druglike space (e.g., high permeability or ECCS class 2) without challenging functional groups. Unfortunately, based on existing mouse data, prospective CL_p IVIVE for complex and nonclassical chemotypes is poor and hardly better than random guessing. This is likely due to complexities such as extrahepatic metabolism and transporter-mediated disposition which are poorly captured by this methodology. With small-molecule drug discovery increasingly evolving toward nonclassical and complex chemotypes, existing CL_p IVIVE methodology will require improvement. While empirical correction factors may bridge the gap in the near future, improved and new in vitro assays, data integration models, and machine learning (ML) methods are increasingly needed to address this challenge and reduce the number of nonclinical pharmacokinetic (PK) studies.

Sensitive and valid assay for reliable evaluation of drug interactions mediated by human organic anion transporter 1 and 3 using 5-carboxyfluorescein

Drug interactions can induce significant clinical impacts, either by increasing adverse effects or by decreasing the therapeutic effect of drugs, and thus, need to be explored thoroughly. Clinically significant drug interactions can be induced by organic anion transporter 1 (OAT1) and OAT3 when concomitant medications competitively interact with the transporters. The purposes of this study were to develop and validate a sensitive and selective analytical method for 5-carboxyfluorescein (5-CF) and optimize the experimental conditions for interaction studies. An analytical method using high-performance liquid chromatography (HPLC) equipped with a fluorescence detector was validated for accuracy, precision, matrix effect, recovery, stability, dilutional integrity, and carry-over effect. In addition, the 5-CF concentration, incubation period, and washing conditions for interaction study were optimized. Using a valid analytical method and optimized conditions, we performed an interaction study for OAT1 and OAT3 using 26 test articles. Some of the test articles showed strong inhibitory potency for the transporters, with IC₅₀ values close to or less than 10 μM. The valid analysis method and optimized systems developed in this study can be utilized to improve the predictability of drug interactions in humans and consequently aid in successful disease treatment by maintaining appropriate systemic exposures.

Multiorgan Crystal Deposition of an Amphoteric Drug in Rats Due to Lysosomal Accumulation and Conversion to a Poorly Soluble Hydrochloride Salt

Poor solubility of drug candidates mainly affects bioavailability, but poor solubility of drugs and metabolites can also lead to precipitation within tissues, particularly when high doses are tested. RO0728617 is an amphoteric compound bearing basic and acidic moieties that has previously demonstrated good solubility at physiological pH but underwent widespread crystal deposition in multiple tissues in rat toxicity studies. The aim of our investigation was to better characterize these findings and their underlying mechanism(s), and to identify possible screening methods in the drug development process. Main microscopic features observed in rat RO0728617 toxicity studies were extensive infiltrates of crystal-containing macrophages in multiple organs. Matrix-assisted laser desorption/ionization Fourier transform ion cyclotron resonance mass spectrometry revealed that these crystals contained the orally administered parent compound, and locality was confirmed to be intracytoplasmic and partly intralysosomal by electron microscopic examination. Crystal formation was explained by lysosomal accumulation of the compound followed by precipitation of the hydrochloride salt under physiological conditions in the lysosomes, which have a lower pH and higher chloride concentration in comparison to the cytosol. This study demonstrates that risk of drug precipitation can be assessed by comparing the estimated lysosomal drug concentration at a given dose with the solubility of the compound at lysosomal conditions.

Modulation of γ-Secretase Activity by a Carborane-Based Flurbiprofen Analogue

All over the world, societies are facing rapidly aging populations combined with a growing number of patients suffering from Alzheimer's disease (AD). One focus in pharmaceutical research to address this issue is on the reduction of the longer amyloid-β (Aβ) fragments in the brain by modulation of γ-secretase, a membrane-bound protease. R-Flurbiprofen (tarenflurbil) was studied in this regard but failed to show significant improvement in AD patients in a phase 3 clinical trial. This was mainly attributed to its low ability to cross the blood-brain barrier (BBB). Here, we present the synthesis and in vitro evaluation of a racemic meta-carborane analogue of flurbiprofen. By introducing the carborane moiety, the hydrophobicity could be shifted into a more favourable range for the penetration of the blood-brain barrier, evident by a logD7.4 value of 2.0. Furthermore, our analogue retained γ-secretase modulator activity in comparison to racemic flurbiprofen in a cell-based assay. These findings demonstrate the potential of carboranes as phenyl mimetics also in AD research.

Identification and Preclinical Development of a 2,5,6-Trisubstituted Fluorinated Pyridine Derivative as a Radioligand for the Positron Emission Tomography Imaging of Cannabinoid Type 2 Receptors

Despite the broad implications of the cannabinoid type 2 receptor (CB2) in neuroinflammatory processes, a suitable CB2-targeted probe is currently lacking in clinical routine. In this work, we synthesized 15 fluorinated pyridine derivatives and tested their binding affinities toward CB2 and CB1. With a sub-nanomolar affinity (K_i for CB2) of 0.8 nM and a remarkable selectivity factor of >12,000 over CB1, RoSMA-18-d₆ exhibited outstanding in vitro performance characteristics and was radiofluorinated with an average radiochemical yield of 10.6 ± 3.8% (n = 16) and molar activities ranging from 52 to 65 GBq/μmol (radiochemical purity > 99%). [¹⁸F]RoSMA-18-d₆ showed exceptional CB2 attributes as demonstrated by in vitro autoradiography, ex vivo biodistribution, and positron emission tomography (PET). Further, [¹⁸F]RoSMA-18-d₆ was used to detect CB2 upregulation on postmortem human ALS spinal cord tissues. Overall, these results suggest that [¹⁸F]RoSMA-18-d₆ is a promising CB2 PET radioligand for clinical translation.

Discovery of Pyrido[2,3-b]indole Derivatives with Gram-Negative Activity Targeting Both DNA Gyrase and Topoisomerase IV

The rise of multidrug resistant (MDR) Gram-negative (GN) pathogens and the decline of available antibiotics that can effectively treat these severe infections are a major threat to modern medicine. Developing novel antibiotics against MDR GN pathogens is particularly difficult as compounds have to permeate the GN double membrane, which has very different physicochemical properties, and have to circumvent a plethora of resistance mechanisms such as multiple efflux pumps and target modifications. The bacterial type II topoisomerases DNA gyrase (GyrA₂B₂) and Topoisomerase IV (ParC₂E₂) are highly conserved targets across all bacterial species and validated in the clinic by the fluoroquinolones. Dual inhibitors targeting the ATPase domains (GyrB/ParE) of type II topoisomerases can overcome target-based fluoroquinolone resistance. However, few ATPase inhibitors are active against GN pathogens. In this study, we demonstrated a successful strategy to convert a 2-carboxamide substituted azaindole chemical scaffold with only Gram-positive (GP) activity into a novel series with also potent activity against a range of MDR GN pathogens. By systematically fine-tuning the many physicochemical properties, we identified lead compounds such as 17r with a balanced profile showing potent GN activity, high aqueous solubility, and desirable PK features. Moreover, we showed the bactericidal efficacy of 17r using a neutropenic mouse thigh infection model.

Matched Molecular Series Analysis for ADME Property Prediction

Generation and prioritization of new molecules are the most central part of the drug design process. Matched molecular series analysis (MMSA) has recently been proposed as a formal approach that captures both of these key elements of design. In order to better understand the power of MMSA and its specific limitations, we here evaluate its performance as an ADME property prediction tool. We use four large and diverse inhouse data sets, logD, microsomal clearance, CYP2C9, and CYP3A4 inhibition. MMSA follows the concept of parallel structure-activity relationship (SAR), where if two identical substituent series on different scaffolds show similarity in their property profiles, SAR from one series can be transferred to the other series. We test four different similarity metrics to identify pairs of molecular series where information can be transferred. We find that the best prediction performance is achieved by a combination of centered root-mean-square deviation (cRMSD) and a network score approach previously published by Keefer et al. However, cRMSD alone strikes the best balance between accuracy and the number of predictions that can be made. We identify statistical metrics that allow estimating when MMSA predictions will work, similar to the well-known applicability domain concept in machine learning. MMSA achieves a prediction accuracy that is comparable to a standard machine-learning model and matched molecular pair analysis. In contrast to machine learning, however, it is very easy to understand where MMSA predictions are coming from. Finally, to prospectively test the power of MMSA, we retested compounds that were strong outliers in the initial predictions and show how the MMSA model can help to identify erroneous data points.

A Comparative Study of in vitro Assays for Predicting the Nonspecific Binding of PET Imaging Agents in vivo

Nonspecific binding (NSB) is a key parameter in optimizing PET imaging tracers. We compared the ability to predict NSB of three available methods: LIMBA, rat f_u,brain , and CHI(IAM). Even though NSB is often associated with lipophilicity, we observed that logD does not correlate with any of these assays, clearly indicating that lipophilicity, while influencing NSB, is insufficient to predict it. A cross-comparison of the methods showed that all three correlate and are useful predictors of NSB. The three assays, however, rank the molecules slightly differently, illustrating the challenge of comparing molecules within a narrow chemical space. We also noted that CHI(IAM) values more effectively predict V_NS , a measure of in vivo NSB in the human brain. CHI(IAM) measurements might be a closer model of the actual physicochemical interaction between PET tracer candidates and cell membranes, and seems to be the method of choice for the optimization of in vivo NSB.

Establishment of an In Vitro-In Vivo Correlation for Melanin Binding and the Extension of the Ocular Half-Life of Small-Molecule Drugs

A large variety of drugs bind effectively to melanin, and this binding influences their ocular pharmacokinetic and distribution profiles. We aimed to establish a correlation between in vitro melanin binding and in vivo ocular pharmacokinetics (PK). The extent of melanin binding in vitro was determined for a set of model drugs; binding kinetics and binding isotherms were generated and fitted to a mechanistic model to derive the drug-melanin binding parameters (B_max, K_D, k_on, and k_off). In addition, in vitro ADME properties such as cellular permeability, P-glycoprotein-mediated efflux, plasma protein binding, and octanol partition coefficients were determined. Moreover, cellular uptake was measured in the nonpigmented ARPE-19 cells and in lightly pigmented human epidermal melanocytes. Finally, in vivo ocular PK studies were performed in albino and pigmented rats using intravenous injections. Substantial drug enrichment accompanied by a very long residence time was observed in pigmented ocular tissues, which could be linked to the melanin binding determined in vitro and to the intracellular drug uptake into the pigmented cells. The resulting ocular PK profile is shown to be a consequence of the interplay of melanin binding with concurrent processes such as systemic clearance, plasma protein binding, cellular permeation, P-glycoprotein efflux, pH partitioning, and tissue binding. Understanding this interplay at a mechanistic level could help in the rational design and development of new small-molecule drug candidates with the desired PK/pharmacodynamic profile to target the back of the eye.

Kinetics of lipid bilayer permeation of a series of ionisable drugs and their correlation with human transporter-independent intestinal permeability

For low molecular weight drugs, lipid bilayer permeation is considered the major route for in vivo cell barrier passage. We recently introduced a fluorescence assay with liposomes to determine permeation kinetics of ionisable compounds across the lipid bilayer by monitoring drug-induced pH changes inside the liposomes. Here, we determined the permeability coefficients (P_FLipP, FLipP for "Fluorescence Liposomal Permeability") across 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers of 35 ionisable drugs at pH6.0 and compared them to available in vivo human jejunal permeability (P_eff) data. P_FLipP values were furthermore compared with published Caco-2 cell permeability coefficients (P_Caco-2), permeability coefficients determined with the parallel artificial membrane permeability assay (PAMPA) and with log D (pH6.0). The log P_FLipP, corrected for predicted para-cellular diffusion, and log P_Caco-2 correlated best with log P_eff, with similar adjusted R² (0.75 and 0.74, n=12). Our results suggest that transporter-independent intestinal drug absorption is predictable from liposomal permeability.

31P and 1H NMR Studies of the Molecular Organization of Lipids in the Parallel Artificial Membrane Permeability Assay

The parallel artificial membrane permeability assay (PAMPA) has emerged as a widely used primary in vitro screen for passive permeability of potential drug candidates. However, the molecular structure of the permeation barrier (consisting of a filter-supported dodecane-egg lecithin mixture) has never been characterized. Here, we investigated the long-range order of phospholipids in the PAMPA barrier by means of ³¹P static solid-state NMR. Diffusion constants of PAMPA membrane components were derived from liquid state NMR and, in addition, drug distribution between the PAMPA lipid phase and buffer (log D_PAMPA at pH 7.4) was systematically investigated. Increasing concentration of n-dodecane to the system egg lecithin-water (lamellar phase, L_α) induces formation of inverted hexagonal (H_ii) and isotropic phases. At n-dodecane concentrations matching those used in PAMPA (9%, w/v) a purely "isotropic" phase was observed corresponding to lipid aggregates with a diameter in the range 4-7 nm. Drug distribution studies indicate that these reverse micelles facilitate the binding to, and in turn the permeation across, the PAMPA dodecane barrier, in particular for amphiphilic solutes. The proposed model for the molecular architecture and function of the PAMPA barrier provides a fundamental, hitherto missing framework to evaluate the scope but also limitations of PAMPA for the prediction of in vivo membrane permeability.