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

How to Personalize General Anesthesia-A Prospective Theoretical Approach to Conformational Changes of Halogenated Anesthetics in Fire Smoke Poisoning

Smoke intoxication is a central event in mass burn incidents, and toxic smoke acts at different levels of the body, blocking breathing and oxygenation. The majority of these patients require early induction of anesthesia to preserve vital functions. We studied the influence of hemoglobin (HMG) and myoglobin (MGB) blockade by hydrochloric acid (HCl) in an interaction model with gaseous anesthetics using molecular docking techniques. In the next part of the study, molecular dynamics (MD) simulations were performed on the top-scoring ligand-receptor complexes to investigate the stability of the ligand-receptor complexes and the interactions between ligands and receptors in more detail. Through docking analysis, we observed that hemoglobin creates more stable complexes with anesthetic gases than myoglobin. Intoxication with gaseous hydrochloric acid produces conformational and binding energy changes of anesthetic gases to the substrate (both the pathway and the binding site), the most significant being recorded in the case of desflurane and sevoflurane, while for halothane and isoflurane, they remain unchanged. According to our theoretical model, the selection of anesthetic agents for patients affected by fire smoke containing hydrochloric acid is critical to ensure optimal anesthetic effects. In this regard, our model suggests that halothane and isoflurane are the most suitable choices for predicting the anesthetic effects in such patients when compared to sevoflurane and desflurane.

Beyond Rule of Five and PROTACs in Modern Drug Discovery: Polarity Reducers, Chameleonicity, and the Evolving Physicochemical Landscape

Developing orally bioavailable drugs demands an understanding of absorption in early drug development. Traditional methods and physicochemical properties optimize absorption for rule of five (Ro5) compounds; beyond rule of five (bRo5) drugs necessitate advanced tools like the experimental measure of exposed polarity (EPSA) and the AbbVie multiparametric score (AB-MPS). Analyzing AB-MPS and EPSA against ∼1000 compounds with human absorption data and ∼10,000 AbbVie tool compounds (∼1000 proteolysis targeting chimeras or PROTACs, ∼7000 Ro5s, and ∼2000 bRo5s) revealed new patterns of physicochemical trends. We introduced a high-throughput "polarity reduction" descriptor: ETR, the EPSA-to-topological polar surface area (TPSA) ratio, highlights unique bRo5 and PROTAC subsets for specialized drug design strategies for effective absorption. Our methods and guidelines refine drug design by providing innovative in vitro approaches, enhancing physicochemical property optimization, and enabling accurate predictions of intestinal absorption in the complex bRo5 domain.

Leveraging nitrogen occurrence in approved drugs to identify structural patterns

BACKGROUND

The process of drug development and discovery is costly and slow. Although an understanding of molecular design principles and biochemical processes has progressed, it is essential to minimize synthesis-testing cycles. An effective approach is to analyze key heteroatoms, including oxygen and nitrogen. Herein, we present an analysis focusing on the utilization of nitrogen atoms in approved drugs.

RESEARCH DESIGN AND METHODS

The present work examines the frequency, distribution, prevalence, and diversity of nitrogen atoms in a dataset comprising 2,049 small molecules approved by different regulatory agencies (FDA and others). Various types of nitrogen atoms, such as sp3-, sp2-, sp-hybridized, planar, ring, and non-ring are included in this investigation.

RESULTS

The results unveil both previously reported and newly discovered patterns of nitrogen atom distribution around the center of mass in the majority of drug molecules.

CONCLUSIONS

This study has highlighted intriguing trends in the role of nitrogen atoms in drug design and development. The majority of drugs contain 1-3 nitrogen atoms within 5Å from the center of mass (COM) of a molecule, with a higher preference for the ring and planar nitrogen atoms. The results offer invaluable guidance for the multiparameter optimization process, thus significantly contributing toward the conversion of lead compounds into potential drug candidates.

Molecular chameleons in drug discovery

Molecular chameleons possess a flexibility that allows them to dynamically shield or expose polar functionalities in response to the properties of the environment. Although the concept of molecular chameleons was introduced already in 1970, interest in them has grown considerably since the 2010s, when drug discovery has focused to an increased extent on new chemical modalities. Such modalities include cyclic peptides, macrocycles and proteolysis-targeting chimeras, all of which reside in a chemical space far from that of traditional small-molecule drugs. Both cell permeability and aqueous solubility are required for the oral absorption of drugs. Engineering these properties, and potent target binding, into the larger new modalities is a more daunting task than for traditional small-molecule drugs. The ability of chameleons to adapt to different environments may be essential for success. In this Review, we provide both general and theoretical insights into the realm of molecular chameleons. We discuss why chameleons have come into fashion and provide a do-it-yourself toolbox for their design; we then provide a glimpse of how advanced in silico methods can support molecular chameleon design.

Controlling the Chameleonic Behavior and Membrane Permeability of Cyclosporine Derivatives via Backbone and Side Chain Modifications

Some macrocycles exhibit enhanced membrane permeability through conformational switching in different environmental polarities, a trait known as chameleonic behavior. In this study, we demonstrate specific backbone and side chain modifications that can control chameleonic behavior and passive membrane permeability using a cyclosporin O (CsO) scaffold. To quantify chameleonic behavior, we used a ratio of the population of the closed conformation obtained in polar solvent and nonpolar solvent for each CsO derivative. We found that β-hydroxylation at position 1 (1 and 3) can encode chameleonicity and improve permeability. However, the conformational stabilization induced by adding an additional transannular H-bond (2 and 5) leads to a much slower rate of membrane permeation. Our CsO scaffold provides a platform for the systematic study of the relationship among conformation, membrane permeability, solubility, and protein binding. This knowledge contributes to the discovery of potent beyond the rule of five (bRo5) macrocycles capable of targeting undruggable targets.

Targeted Protein Degradation: Advances, Challenges, and Prospects for Computational Methods

The therapeutic approach of targeted protein degradation (TPD) is gaining momentum due to its potentially superior effects compared with protein inhibition. Recent advancements in the biotech and pharmaceutical sectors have led to the development of compounds that are currently in human trials, with some showing promising clinical results. However, the use of computational tools in TPD is still limited, as it has distinct characteristics compared with traditional computational drug design methods. TPD involves creating a ternary structure (protein-degrader-ligase) responsible for the biological function, such as ubiquitination and subsequent proteasomal degradation, which depends on the spatial orientation of the protein of interest (POI) relative to E2-loaded ubiquitin. Modeling this structure necessitates a unique blend of tools initially developed for small molecules (e.g., docking) and biologics (e.g., protein-protein interaction modeling). Additionally, degrader molecules, particularly heterobifunctional degraders, are generally larger than conventional small molecule drugs, leading to challenges in determining drug-like properties like solubility and permeability. Furthermore, the catalytic nature of TPD makes occupancy-based modeling insufficient. TPD consists of multiple interconnected yet distinct steps, such as POI binding, E3 ligase binding, ternary structure interactions, ubiquitination, and degradation, along with traditional small molecule properties. A comprehensive set of tools is needed to address the dynamic nature of the induced proximity ternary complex and its implications for ubiquitination. In this Perspective, we discuss the current state of computational tools for TPD. We start by describing the series of steps involved in the degradation process and the experimental methods used to characterize them. Then, we delve into a detailed analysis of the computational tools employed in TPD. We also present an integrative approach that has proven successful for degrader design and its impact on project decisions. Finally, we examine the future prospects of computational methods in TPD and the areas with the greatest potential for impact.

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.

Bifunctional degraders of cyclin dependent kinase 9 (CDK9): Probing the relationship between linker length, properties, and selective protein degradation

Cyclin-dependent kinase 9 (CDK9) is a promising therapeutic target in multiple cancer types, including acute myeloid leukemia (AML). Protein degraders, also known as proteolysis targeting chimeras (PROTACs), have emerged as tools for the selective degradation of cancer targets, including CDK9, complementing the activity of traditional small-molecule inhibitors. These compounds typically incorporate previously reported inhibitors and a known E3 ligase ligand to induce ubiquitination and subsequent degradation of the target protein. Although many protein degraders have been reported in the literature, the properties of the linker necessary for efficient degradation still require special attention. In this study, a series of protein degraders was developed, employing the clinically tested CDK inhibitor AT7519. The purpose of this study was to examine the effect that linker composition, specifically chain length, would have on potency. In addition to establishing a baseline of activity for various linker compositions, two distinct homologous series, a fully alkyl series and an amide-containing series, were prepared, demonstrating the dependence of degrader potency in these series on linker length and the correlation with predicted physicochemical properties.

Delivering on the promise of protein degraders

Over the past 3 years, the first bivalent protein degraders intentionally designed for targeted protein degradation (TPD) have advanced to clinical trials, with an initial focus on established targets. Most of these clinical candidates are designed for oral administration, and many discovery efforts appear to be similarly focused. As we look towards the future, we propose that an oral-centric discovery paradigm will overly constrain the chemical designs that are considered and limit the potential to drug novel targets. In this Perspective, we summarize the current state of the bivalent degrader modality and propose three categories of degrader designs, based on their likely route of administration and requirement for drug delivery technologies. We then describe a vision for how parenteral drug delivery, implemented early in research and supported by pharmacokinetic-pharmacodynamic modelling, can enable exploration of a broader drug design space, expand the scope of accessible targets and deliver on the promise of protein degraders as a therapeutic modality.

It's ok to be outnumbered - sub-stoichiometric modulation of homomeric protein complexes

An arsenal of molecular tools with increasingly diversified mechanisms of action is being developed by the scientific community to enable biological interrogation and pharmaceutical modulation of targets and pathways of ever increasing complexity. While most small molecules interact with the target of interest in a 1 : 1 relationship, a noteworthy number of recent examples were reported to bind in a sub-stoichiometric manner to a homomeric protein complex. This approach requires molecular understanding of the physiologically relevant protein assemblies and in-depth characterization of the compound's mechanism of action. The recent literature examples summarized here were selected to illustrate methods used to identify and characterize molecules with such mechanisms. The concept of one small molecule targeting a homomeric protein assembly is not new but the subject deserves renewed inspection in light of emerging technologies and increasingly diverse target biology, to ensure relevant in vitro systems are used and valuable compounds with potentially novel sub-stoichiometric mechanisms of action aren't overlooked.

Conformational Sampling Deciphers the Chameleonic Properties of a VHL-Based Degrader

Chameleonicity (the capacity of a molecule to adapt its conformations to the environment) may help to identify orally bioavailable drugs in the beyond-Rule-of-5 chemical space. Computational methods to predict the chameleonic behaviour of degraders have not yet been reported and the identification of molecular chameleons still relies on experimental evidence. Therefore, there is a need to tune predictions with experimental data. Here, we employ PROTAC-1 (a passively cell-permeable degrader), for which NMR and physicochemical data prove the chameleonic behaviour, to benchmark the capacity of two conformational sampling algorithms and selection schemes. To characterize the conformational ensembles in both polar and nonpolar environments, we compute three molecular properties proven to be essential for cell permeability: conformer shape (radius of gyration), polarity (3D PSA), and the number of intramolecular hydrogen bonds. Energetic criteria were also considered. Infographics monitored the simultaneous variation of those properties in computed and NMR conformers. Overall, we provide key points for tuning conformational sampling tools to reproduce PROTAC-1 chameleonicity according to NMR evidence. This study is expected to improve the design of PROTAC drugs and the development of computational sustainable strategies to exploit the potential of new modalities in drug discovery.

Permeability of New Antifungal Fluconazole Derivatives through a Lipophilic Membrane: Experiment and Modeling

Relationships between the structures of molecules and their properties form the basis of modern chemistry and lay the foundation for structure-based drug design. Being the main two determinants of bioavailability, solubility and permeability of drugs are widely investigated experimentally and predicted from physicochemical parameters and structural descriptors. In the present study, we measure the passive diffusion permeability of a series of new fluconazole derivatives with triazole and thiazolo-pyrimidine moieties connected by different linker bridges through the PermeaPad barrier-a relatively new biomimetic lipophilic membrane that has been increasingly used in recent years. The permeability coefficients of new derivatives are shown to be dependent both on the structure of the linker fragment and on the substituent in the phenyl ring of the thiazolo-pyrimidine moiety. The impact of the compound ionization state on the permeability is revealed. Reliable correlations of the permeability with the antifungal activity and distribution coefficient are found. In addition, the solubility-diffusion approach is shown to be able to successfully predict the permeability of the studied derivatives. The obtained results can be considered another step in the development of permeability databases and design of schemes for in vitro permeability prediction.

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.

Key Considerations in Targeted Protein Degradation Drug Discovery and Development

Targeting proteins’ enzymatic functions with small molecule inhibitors, as well as functions of receptor proteins with small-molecule agonists and antagonists, were the major forms of small-molecule drug development. These small-molecule modulators are based on a conventional occupancy-driven pharmacological approach. For proteome space traditionally considered undruggable by small-molecule modulators, such as enzymes with scaffolding functions, transcription factors, and proteins that lack well-defined binding pockets for small molecules, targeted protein degraders offer the opportunity to drug the proteome with an event-driven pharmacological approach. A degrader molecule, either PROTAC or molecular glue, brings the protein of interest (POI) and E3 ubiquitin ligase in close proximity and engages the ubiquitin-proteasome system (UPS), the cellular waste disposal system for the degradation of the POI. For the development of targeted protein degraders to meet therapeutic needs, several aspects will be considered, namely, the selective degradation of disease-causing proteins, the oral bioavailability of degraders beyond Lipinski’s rule of five (bRo5) scope, demands of new E3 ubiquitin ligases and molecular glue degraders, and drug resistance of the new drug modality. This review will illustrate several under-discussed key considerations in targeted protein degradation drug discovery and development: 1) the contributing factors for the selectivity of PROTAC molecules and the design of PROTACs to selectively degrade synergistic pathological proteins; 2) assay development in combination with a multi-omics approach for the identification of new E3 ligases and their corresponding ligands, as well as molecular glue degraders; 3) a molecular design to improve the oral bioavailability of bRo5 PROTACs, and 4) drug resistance of degraders.

Designing Soluble PROTACs: Strategies and Preliminary Guidelines

Solubility optimization is a crucial step to obtaining oral PROTACs. Here we measured the thermodynamic solubilities (log S) of 21 commercial PROTACs. Next, we measured BRlogD and log k_w^IAM (lipophilicity), EPSA, and Δ log k_w^IAM (polarity) and showed that lipophilicity plays a major role in governing log S, but a contribution of polarity cannot be neglected. Two-/three-dimensional descriptors calculated on conformers arising from conformational sampling and steered molecular dynamics failed in modeling solubility. Infographic tools were used to identify a privileged region of soluble PROTACs in a chemical space defined by BRlogD, log k_w^IAM and topological polar surface area, while machine learning provided a log S classification model. Finally, for three pairs of PROTACs we measured the solubility, lipophilicity, and polarity of the building blocks and identified the limits of estimating PROTAC solubility from the synthetic components. Overall, this paper provides promising guidelines for optimizing PROTAC solubility in early drug discovery programs.

Predicting Solubility of Newly-Approved Drugs (2016–2020) with a Simple ABSOLV and GSE(Flexible-Acceptor) Consensus Model Outperforming Random Forest Regression

This study applies the ‘Flexible-Acceptor’ variant of the General Solubility Equation, GSE(Φ,B), to the prediction of the aqueous intrinsic solubility, log₁₀S₀, of FDA recently-approved (2016–2020) ‘small-molecule’ new molecular entities (NMEs). The novel equation had been shown to predict the solubility of drugs beyond Lipinski’s ‘Rule of 5’ chemical space (bRo5) to a precision nearly matching that of the Random Forest Regression (RFR) machine learning method. Since then, it was found that the GSE(Φ,B) appears to work well not only for bRo5 NMEs, but also for Ro5 drugs. To put context to GSE(Φ,B), Yalkowsky’s GSE(classic), Abraham’s ABSOLV, and Breiman’s RFR models were also applied to predict log₁₀ S₀ of 72 newly-approve NMEs, for which useable reported solubility values could be accessed (nearly 60% from FDA New Drug Application published reports). Except for GSE (classic), the prediction models were retrained with an enlarged version of the Wiki-pS₀ database (nearly 400 added log₁₀ S₀ entries since our recent previous study). Thus, these four models were further validated by the additional independent solubility measurements which the newly-approved drugs introduced. The prediction methods ranked RFR ~ GSE (Φ,B) > ABSOLV > GSE (classic) in performance. It was further demonstrated that the biases generated in the four separate models could be nearly eliminated in a consensus model based on the average of just two of the methods: GSE (Φ,B) and ABSOLV. The resulting consensus prediction equation is simple in form and can be easily incorporated into spreadsheet calculations. Even more significant, it slightly outperformed the RFR method.

Managing Experimental 3D Structures in the Beyond-Rule-of-5 Chemical Space: The Case of Rifampicin

The beyond-Rule-of-5 (bRo5) chemical space is a source of new oral drugs and includes large and flexible compounds. Because of their size and conformational variability, bRo5 molecules assume different privileged conformations in the compartments of human body, i. e., they can exhibit chameleonic properties. The elucidation of the ensemble of 3D structures explored by such molecules under different conditions is therefore critical to check the role played by chameleonicity to modulate cell permeability. Here we characterized the conformational ensembles of rifampicin, a bRo5 drug, in polar and nonpolar solvents and in the solid state. We performed NMR experiments, analyzed their results with a novel algorithm and set-up a pool of ad hoc in silico strategies to investigate crystallographic structures retrieved from the CSD. Moreover, a polarity descriptor often related to permeability (SA-3D-PSA) was calculated for all the conformers and its variation with the environment analyzed. Results showed that the conformational behavior of rifampicin in solution and in the solid state is not superposable. The identification of dynamic intramolecular hydrogen bonds can be assessed by NMR spectroscopy but not by X-ray structures. Moreover, SA-3D-PSA revealed that dynamic IMHBs do not provide rifampicin with chameleonic properties. Overall, this study highlights that the peculiarity of rifampicin, which is cell permeable probably because of the presence of static IMHBs but is devoid of any chameleonic behavior, can be assessed by a proper analysis of experimental 3D structures.

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.

Interplay among Conformation, Intramolecular Hydrogen Bonds, and Chameleonicity in the Membrane Permeability and Cyclophilin A Binding of Macrocyclic Peptide Cyclosporin O Derivatives

A macrocyclic peptide scaffold with well-established structure-property relationship is desirable for tackling undruggable targets. Here, we adopted a natural macrocycle, cyclosporin O (CsO) and its derivatives (CP1-3), and evaluated the impact of conformation on membrane permeability, cyclophilin A (CypA) binding, and the pharmacokinetic (PK) profile. In nonpolar media, CsO showed a similar conformation to cyclosporin A (CsA), a well-known chameleonic macrocycle, but less chameleonic behavior in a polar environment. The weak chameleonicity of CsO resulted in decreased membrane permeability; however, the more rigid conformation of CsO was not detrimental to its PK profile. CsO exhibited a higher plasma concentration than CsA, which resulted from minimal CypA binding and lower accumulation in red blood cells and moderate oral bioavailability (F = 12%). Our study aids understanding of CsO, a macrocyclic peptide that is less explored than CsA but with greater potential for diversity generation and rational design.

Target-Based Evaluation of "Drug-Like" Properties and Ligand Efficiencies

Physicochemical descriptors commonly used to define "drug-likeness" and ligand efficiency measures are assessed for their ability to differentiate marketed drugs from compounds reported to bind to their efficacious target or targets. Using ChEMBL version 26, a data set of 643 drugs acting on 271 targets was assembled, comprising 1104 drug-target pairs having ≥100 published compounds per target. Taking into account changes in their physicochemical properties over time, drugs are analyzed according to their target class, therapy area, and route of administration. Recent drugs, approved in 2010-2020, display no overall differences in molecular weight, lipophilicity, hydrogen bonding, or polar surface area from their target comparator compounds. Drugs are differentiated from target comparators by higher potency, ligand efficiency (LE), lipophilic ligand efficiency (LLE), and lower carboaromaticity. Overall, 96% of drugs have LE or LLE values, or both, greater than the median values of their target comparator compounds.

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.

Steering New Drug Discovery Campaigns: Permeability, Solubility, and Physicochemical Properties in the bRo5 Chemical Space

An increasing number of drug discovery programs concern compounds in the beyond rule of 5 (bRo5) chemical space, such as cyclic peptides, macrocycles, and degraders. Recent results show that common paradigms of property-based drug design need revision to be applied to larger and more flexible compounds. A virtual event entitled "Solubility, permeability and physico-chemical properties in the bRo5 chemical space" was organized to provide preliminary guidance on how to make the discovery of oral drugs in the bRo5 space more effective. The four speakers emphasized the importance of the bRo5 space as a source of new oral drugs and provided examples of experimental and computational methods specifically tailored for design and optimization in this chemical space.

"Flexible-Acceptor" General Solubility Equation for beyond Rule of 5 Drugs

This study describes a novel nonlinear variant of the well-known Yalkowsky general solubility equation (GSE). The modified equation can be trained with small molecules, mostly from the Lipinski Rule of 5 (Ro5) chemical space, to predict the intrinsic aqueous solubility, S₀, of large molecules (MW > 800 Da) from beyond the rule of 5 (bRo5) space, to an accuracy almost equal to that of a recently described random forest regression (RFR) machine learning analysis. The new approach replaces the GSE constant factors in the intercept (0.5), the octanol-water log P (-1.0), and melting point, mp (-0.01) terms with simple exponential functions incorporating the sum descriptor, Φ+B (Kier Φ molecular flexibility and Abraham H-bond acceptor potential). The constants in the modified three-variable (log P, mp, Φ+B) equation were determined by partial least-squares (PLS) refinement using a small-molecule log S₀ training set (n = 6541) of mostly druglike molecules. In this "flexible-acceptor" GSE(Φ,B) model, the coefficient of log P (normally fixed at -1.0) varies smoothly from -1.1 for rigid nonionizable molecules (Φ+B = 0) to -0.39 for typically flexible (Φ ∼ 20, B ∼ 6) large molecules. The intercept (traditionally fixed at +0.5) varies smoothly from +1.9 for completely inflexible small molecules to -2.2 for typically flexible large molecules. The mp coefficient (-0.007) remains practically constant, near the traditional value (-0.01) for most molecules, which suggests that the small-to-large molecule continuum is mainly solvation responsive, apparently with only minor changes in the crystal lattice contributions. For a test set of 32 large molecules (e.g., cyclosporine A, gramicidin A, leuprolide, nafarelin, oxytocin, vancomycin, and mostly natural-product-derived therapeutics used in infectious/viral diseases, in immunosuppression, and in oncology) the modified equation predicted the intrinsic solubility with a root-mean-square error of 1.10 log unit, compared to 3.0 by the traditional GSE, and 1.07 by RFR.

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.