A Fast Chromatographic Method for Estimating Lipophilicity and Ionization in Nonpolar Membrane-Like Environment

Triazoles in Medicinal Chemistry: Physicochemical Properties, Bioisosterism, and Application

Triazole demonstrates distinctive physicochemical properties, characterized by weak basicity, various dipole moments, and significant dual hydrogen bond acceptor and donor capabilities. These features are poised to play a pivotal role in drug-target interactions. The inherent polarity of triazole contributes to its lower logP, suggesting the potential improvement in water solubility. The metabolic stability of triazole adds additional value to drug discovery. Moreover, the metal-binding capacity of the nitrogen atom lone pair electrons of triazole has broad applications in the development of metal chelators and antifungal agents. This Perspective aims to underscore the unique physicochemical attributes of triazole and its application. A comparative analysis involving triazole isomers and other heterocycles provides guiding insights for the subsequent design of triazoles, with the hope of offering valuable considerations for designing other heterocycles in medicinal chemistry.

Chamelogk: A Chromatographic Chameleonicity Quantifier to Design Orally Bioavailable Beyond-Rule-of-5 Drugs

New chemical modalities in drug discovery include molecules belonging to the bRo5 chemical space. Because of their complex and flexible structure, bRo5 compounds often suffer from a poor solubility/permeability profile. Chameleonicity describes the capacity of a molecule to adapt to the environment through conformational changes; the design of molecular chameleons is a medicinal chemistry strategy simultaneously optimizing solubility and permeability. A default method to quantify chameleonicity in early drug discovery is still missing. Here we introduce Chamelogk, an automated, fast, and cheap chromatographic descriptor of chameleonicity. Moreover, we report measurements for 55 Ro5 and bRo5 compounds and validate our method with literature data. Then, selected case studies (macrocycles, nonmacrocyclic compounds, and PROTACs) are used to illustrate the application of Chamelogk in combination with lipophilicity (BRlogD) and polarity (Δ log kwIAM) descriptors. Overall, we show how Chamelogk deserves being included in property-based drug discovery strategies to design oral bioavailable bRo5 compounds.

Examination of the Impact of Triazole Position within Linkers on Solubility and Lipophilicity of a CDK9 Degrader Series

Optimization of degrader properties is often a challenge due to their beyond-rule-of-5 nature. Given the paucity of known E3 ligases and the often-limited choice of ligands with varied chemical structures for a given protein target, degrader linkers represent the best position within the chimeric molecules to modify their overall physicochemical properties. In this work, a series of AT7519-based CDK9 degraders was assembled using click chemistry, facilitating the tuning of aqueous solubility and lipophilicity while retaining their linker type and molecular weight. Using chromatographic logD and kinetic solubility experiments, we show that degraders with similar chemical constitution but varied position of the embedded triazole demonstrate different lipophilicity and aqueous solubility properties. Overall, this work highlights the impact of triazole placement on linker composition through application of click chemistry for degrader synthesis and its ability to be used to promote the achievement of favorable physicochemical properties.

Today's drug discovery and the shadow of the rule of 5

INTRODUCTION

The rule of 5 developed by Lipinski et al., a landmark and prescient piece of scholarship, focused the minds of drug hunters by systematically characterizing the physical make-up of drug molecules for the first time, noting many sub-optimal compounds identified by high-throughput screening practices. Its profound influence on thinking and practices, whilst providing benefit, perhaps etched the guidelines too strongly in the minds of some drug hunters who applied the bounds too literally without understanding the implications of the underlying statistics.

AREAS COVERED

This opinion is based on recent key developments that take thinking, measurements, and standards beyond those first set out, particularly the influences of molecular weight and the understanding, measurement, and calculation of lipophilicity.

EXPERT OPINION

Techniques and technologies for physicochemical estimations set new standards. It is timely to celebrate the significance and influence of the rule of 5, whilst taking thinking to new levels with better characterizations. The shadow of the rule of 5 may be long, but it is not dark, as new measurements, predictions and principles emerge as guiding lights in the design and prioritization of higher-quality molecules redefining the meaning of beyond the rule of 5.

Quantum Mechanical Prediction of Dissociation Constants for Thiazol-2-imine Derivatives

As weak acids or bases, in solution, drug molecules are in either their ionized or nonionized states. A high degree of ionization is essential for good water solubility of a drug molecule and is required for drug-receptor interactions, whereas the nonionized form improves a drug's lipophilicity, allowing the ligand to cross the cell membrane. The penetration of a drug ligand through cell membranes is mainly governed by the pKa of the drug molecule and the membrane environment. In this study, with the aim of predicting the acetonitrile pKa's (pKa(MeCN)) of eight drug-like thiazol-2-imine derivatives, we propose a very accurate and computationally affordable protocol by using several quantum mechanical approaches. Benchmark studies were conducted on a set of training molecules, which were selected from the literature with known pKa(water) and pKa(MeCN). Highly well-correlated pKa values were obtained when the calculations were performed with the isodesmic method at the M062X/6-31G** level of theory in conjunction with SMD solvation model for nitrogen-containing heterocycles. Finally, experimentally unknown pKa(MeCN) values of eight thiazol-2-imine structures, which were previously synthesized by some of us, are proposed.

Ionization and lipophilicity in nonpolar media mimicking the cell membrane interior

Ionization and lipophilicity may vary with the environment. Therefore, in this study we provide some insight in the performances of different experimental techniques (potentiometry, UV-vis, shake-flask and chromatography) to determine ionization and lipophilicity in more nonpolar systems than those commonly used in drug discovery. To this purpose a pool of 11 compounds of pharmaceutical interest was firstly submitted to a few experimental techniques to measure pK_a in water, water/acetonitrile mixtures and pure acetonitrile. Then we measured logP/logD with shake-flask and potentiometry in octanol/water and toluene/water and also determined a chromatographic lipophilicity index (log k'80 PLRP-S) in a nonpolar system. Results show that ionization decreases for both acids and bases in a coherent, significant but not dramatical extent when water is present in the system, but the picture is completely different in pure acetonitrile. Lipophilicity may vary or not with the environment according to the chemical structure of the investigated compounds as also revealed by electrostatic potential maps. Since the internal core of cell membranes is largely nonpolar, our results support the need of extending the pool of physicochemical descriptors to be determined in the various stages of drug discovery programs and indicate some experimental strategies for their determination.

Refinement of Computational Access to Molecular Physicochemical Properties: From Ro5 to bRo5

There is a need of computational tools to rank bRo5 drug candidates in the very early phases of drug discovery when chemical matter is unavailable. In this study, we selected three compounds: (a) a Ro5 drug (Pomalidomide), (b) a bRo5 orally available drug (Saquinavir), and (c) a polar PROTAC (CMP 98) to focus on computational access to physicochemical properties. To provide a benchmark, the three compounds were first experimentally characterized for their lipophilicity, polarity, IMHBs, and chameleonicity. To reproduce the experimental information content, we generated conformer ensembles with conformational sampling and molecular dynamics in both water and nonpolar solvents. Then we calculated Rgyr, 3D PSA, and IMHB number. An innovative pool of strategies for data analysis was then provided. Overall, we report a contribution to close the gap between experimental and computational methods for characterizing bRo5 physicochemical properties.

Rational Control of Molecular Properties Is Mandatory to Exploit the Potential of PROTACs as Oral Drugs

To obtain new oral drugs in the beyond rule of five space, PROTACs among others, molecular properties should be optimized in early drug discovery. Degraders call for design strategies which focus on intramolecular interaction and chameleonicity. In parallel, tailored revalidation of permeability assessment and prediction methods becomes fundamental in this innovative chemical space.

Rifampicin as an example of beyond-rule-of-5 compound: Ionization beyond water and lipophilicity beyond octanol/water

Ionization and lipophilicity in early drug discovery are commonly characterized in water and octanol/water, respectively and thus do not consider the non-polar features of the biomembrane core. This is particularly limiting for bRo5 compounds which may adapt their properties (e.g. ionization and lipophilicity) to the environment. In this paper we used experimental methods to characterize rifampicin for its ionization properties in various water/cosolvent mixtures and in pure MeCN and its lipophilicity in octanol/water and toluene/water systems. Moreover, we also measured log k'80 PLRP-S, a chromatographic index of lipophilicity in non-polar media. Results show that the existence domain of neutral rifampicin is limited compared to the zwitterion, but the lipophilic cationic species is extremely relevant in non-polar environments.

Degraders early developability assessment: face-to-face with molecular properties

Pharmaceutical scientists have huge expectations from heterobifunctional small molecule degraders to treat diseases with an unmet medical need. However, degraders are large and flexible and pose significant challenges in terms of cellular uptake and bioavailability. An efficient property-based design is therefore required to discover new oral degrader medicines. Here, we show the non‑transferability to degraders of in silico tools routinely implemented in small molecule drug discovery programs; and provide ionization, lipophilicity, polarity and chameleonicity data for a series of seven degraders. We also reveal that permeability can be modeled by Δlog k_W^IAM - an experimental polarity descriptor. Overall, the paper is a proof-of-concept that shows to discover new oral degrader drugs ad hoc property-based design strategies are required.

Pharmaceuticals' removal by constructed wetlands: a critical evaluation and meta-analysis on performance, risk reduction, and role of physicochemical properties on removal mechanisms

This paper presents a comprehensive and critical analysis of the removal of pharmaceuticals (PhCs), the governing physicochemical properties, and removal mechanisms in constructed wetlands (CWs). The average removal efficiency of the most widely studied 34 PhCs ranges from 21% to 93%, with the exception of one PhC that exhibited negative removal. Moreover, CWs are effective in significantly reducing the environmental risk caused by many PhCs. Based on risk assessment, 12 PhCs were classified under high risk category (oxytetracycline > ofloxacin > sulfamethoxazole > erythromycin > sulfadiazine > gemfibrozil > ibuprofen > acetaminophen > salicylic acid > sulfamethazine > naproxen > clarithromycin), which could be considered for regular monitoring, water quality standard formulation and control purposes. Biodegradation (aerobic and anaerobic) is responsible for the removal of the majority of PhCs, often in conjunction with other mechanisms (e.g., adsorption/sorption, plant uptake, and photodegradation). The physicochemical properties of molecules play a pivotal role in the elimination processes, and could serve as important predictors of removal. The correlation and multiple linear regression analysis suggest that organic carbon sorption coefficient (Log Koc), octanol-water distribution coefficient (Log Dow), and molecular weight form a good predictive linear regression model for the removal efficiency of PhCs (R² = 0.65, P-value <0.05).

Facts, Patterns, and Principles in Drug Discovery: Appraising the Rule of 5 with Measured Physicochemical Data

The rule of 5 was designed to estimate the likelihood of poor absorption or permeation, noting the impact of poor solubility. This Perspective explores the impact of various physicochemical descriptors and contemporary lipophilicity measurements on permeability and solubility, showing that the distribution coefficient log D_7.4 (rather than log P) is the most impactful parameter. Molecular weight, almost invariably the defining characteristic of "beyond the rule of 5" compounds, has little impact on solubility when log D_7.4 measurements and aromaticity are considered. Predicting permeation is more complex, given passive and carrier transport mechanisms; however, notable patterns of behavior are apparent, giving insight even "beyond the rule of 5". Recommended best practices should involve using the facts (measurements) and the patterns they reveal to establish informative principles rather than fastidious rules.

Updating the portfolio of physicochemical descriptors related to permeability in the beyond the rule of 5 chemical space

Beyond rule of 5 (bRo5) molecules are attracting significant interest in modern drug discovery mostly because many novel targets require large and more flexible structures. The main aim of this paper is the identification of ad hoc bRo5 physicochemical descriptors of ionization, lipophilicity, polarity and chameleonicity and their measurement. We used different methods to collect ionization (pK_a measures and log k'80 PLRP-S trends), lipophilicity (in octanol/water, in apolar systems and in biomimetic environments), polarity (Δlog P_oct-tol, EPSA and Δlog K_W^IAM) and chameleonicity (ChameLogD) descriptors for 26 bRo5 drugs. A second aim was to check the relationship between physicochemical descriptors and permeability for a subset of compounds for which solid permeability values are reported in the literature. Results showed that the physicochemical profile in the bRo5 chemical space is often experimentally accessible, albeit more tools are required to overcome limitations of individual methods. For the investigated compounds, permeability is governed by Δlog P_oct-tol and preliminary data support that chameleonicity could also have an impact.

Experimental lipophilicity for beyond Rule of 5 compounds

Aim: To set up a chromatographic strategy for the determination of log P for beyond Rule of 5 (bRo5) drugs. Materials & methods: Capacity factors measured by reverse phase-HPLC. Balance of intermolecular interactions governing systems assessed by partial least squares regression (PLSR) coupled with block relevance anaysis (PLSR-BR) and multiblock PLSR (MBPLSR). Determination of virtual log P obtained through conformational sampling. Results: log k′60 is highly correlated with log P for a dataset of 36 Ro5 compliant compounds (R2 = 0.93, Q2 = 0.90). We refer to the value generated via this method as BRlogP. The balance of intermolecular forces controlling BRlogP and log P are very similar. The ElogPs measured for the bRo5 dataset are significantly higher than corresponding BRlogP. Conclusion: The combination of BRlogP and ElogP provides an experimental lipophilicity range for bRo5 compounds.

Permeability prediction for zwitterions via chromatographic indexes and classification into 'certain' and 'uncertain'

Background: The development of zwitterions to a drug is likely to be more challenging than compounds of other charge types. Results: Two chromatographic indexes (log k'80 PLRP-S and log K_W^IAM) can be successfully used as permeability classifiers of ampholytes. Moreover, a pragmatic classification into ordinary ampholytes; zwitterions 'certain' (i.e., the zwitterionic species is dominant in the physiological pH range); and zwitterions 'uncertain' (multiple species are present in the physiological pH range) enables to study the permeability of ampholytic compounds in relation to species distribution. Methodology: Potentiometry (pK_a), reversed-phase (RP)-chromatography, tri-layer parallel artificial membrane permeability assays, quantitative structure-property relationships (QSPR) and block relevance (BR) analysis, receiver operating characteristic (ROC) curves. Conclusion: Structures considered as poorly permeable like zwitterions can be integrated in drug discovery programs by applying ad hoc experimental and computational tools.

Intramolecular hydrogen bonding: An opportunity for improved design in medicinal chemistry

Recent literature shows that intramolecular hydrogen bond (IMHB) formation can positively impact upon the triad of permeability, solubility, and potency of drugs and candidates. IMHB modulation can be applied to compounds in any chemical space as a means for discovering drug candidates with both acceptable potency and absorption, distribution, metabolism, and excretion-Tox profiles. Integrating IMHB formation in design of drugs is, therefore, an exciting and timely challenge for modern medicinal chemistry. In this review, we first provide some background about IMHBs from the medicinal chemist's point of view and highlight some IMHB-associated misconceptions. Second, we propose a classification of IMHBs for drug discovery purposes, review the most common in silico tactics to include IMHBs in lead optimization and list some experimental physicochemical descriptors, which quantify the propensity of compounds to form IMHBs. By focusing on the compounds size and the number of IMHBs that can potentially be formed, we also outline the major difficulties encountered when designing compounds based on the inclusion of IMHBs. Finally, we discuss recent case studies illustrating the application of IMHB to optimize cell permeability and physicochemical properties of small molecules, cyclic peptides and macrocycles.

Log P as a tool in intramolecular hydrogen bond considerations

Intramolecular hydrogen bonding (IMHB) considerations are gaining relevance in drug discovery and a molecular descriptor which can predict very early the capacity of a compound to form IMHB is needed to speed up the optimization process of drug candidates. Although log P_oct is largely used for optimization purposes, in this paper we firstly use the Block Relevance (BR) analysis to theoretically show how log P_oct is not a convenient choice to assess IMHB properties of candidates. Then we discuss the limits of log P_oct and introduce Δlog P_oct-tol, i.e. the difference between log P_oct and log P_tol (the logarithm of the partition coefficient in the toluene/water system). Finally, we provided some examples also including bRo5 protease inhibitors, to clarify how to interpret Δlog P_oct-tol values.

Learning how to use IAM chromatography for predicting permeability

The interest for IAM (Immobilized Artificial Membranes) chromatography in the prediction of drug permeability is increasing. Here we firstly set-up a dataset of 253 molecules including neutral and ionized drugs and few organic compounds for which we either measured or retrieved from the literature IAM.PC.DD2 log K_w^IAM data. Then we applied block relevance (BR) analysis to extract from PLS models the relative contribution of intermolecular forces governing log K_w^IAM and Δlog K_w^IAM (a combined descriptor calculated from log K_w^IAM). Finally, the relationship between log K_w^IAM, Δlog K_w^IAM and passive permeability determined in both PAMPA and MDCK-LE systems was looked for. Models provided the basis for a rational application of IAM chromatography in permeability prediction.

High throughput methods to measure the propensity of compounds to form intramolecular hydrogen bonding

Implementation of IMHB considerations in drug discovery needs robust and validated descriptors to experimentally verify the propensity of compounds to exhibit IMHBs. The first part of the paper presents an overview of the most common techniques to measure the propensity of compounds to form IMHBs. Then we review and discuss recently proposed high throughput (HT) physicochemical descriptors (i.e. Δlog Poct-tol, EPSA and log k'80 PLRP-S) which provide the same information. Analysis of the available data enabled us to extract guidelines for the application of these descriptors in drug discovery programs.

Updating molecular properties during early drug discovery

Current multiparameter optimization (MPO) strategies make use of few experimental physicochemical descriptors (i.e., solubility at physiological pH and lipophilicity in the octanol/water system). Here, we show how new trends in drug discovery (i.e., large and flexible molecules for 'difficult' targets) call for the integration of ad hoc descriptors in MPO approaches. In particular, to rank, select, and optimize drug candidates, it could be relevant to have experimental data relating to the acid-base properties and the folding of the molecule to mask polar groups (so-called 'chameleonic' properties). We propose two strategies to quantify ionization and chameleonic properties and discuss their practical integration in property criteria profiles.