Exploring the chemical space of orally bioavailable PROTACs.

A principal challenge in the discovery of proteolysis targeting chimeras (PROTACs) as oral medications is their bioavailability. To facilitate drug design, it is therefore essential to identify the chemical space where orally bioavailable PROTACs are more likely to be situated. To this aim, we extracted structure-bioavailability insights from published data using traditional 2D descriptors, thereby shedding light on their potential and limitations as drug design tools. Subsequently, we describe cutting-edge experimental, computational and hybrid design strategies based on 3D descriptors, which show promise for enhancing the probability of discovering PROTACs with high oral bioavailability.

DegraderTCM: A Computationally Sparing Approach for Predicting Ternary Degradation Complexes.

Proteolysis targeting chimeras (PROTACs or degraders) represent a novel therapeutic modality that has raised interest thanks to promising results and currently undergoing clinical testing. PROTACs induce the selective proteasomal degradation of undesired proteins by the formation of ternary complexes (TCs). Having knowledge of the 3D structure of TCs is crucial for the design of PROTAC drugs. Here, we describe DegraderTCM, a new computational method for modeling PROTAC-mediated TCs that requires low computational power and provides sound results in a short time span. We validated DegraderTCM against a selected set of experimentally determined structures and defined a method to predict the PROTAC degradation activity based on the computed TC structure. Finally, we modeled TCs of known degraders holding significance for defining the method's applicability domain. A retrospective analysis of structure-activity relationships unveiled possibilities for utilizing DegraderTCM in the initial stages of designing novel PROTAC drugs.

High-throughput and reliable assessments of the ionization constant of monoprotic organic acids through an arginine based mixed mode HPLC.

The charge state of a molecule is the single most prominent attribute ruling out its interactions with the surrounding environment. In a previous study, the retention of acids on the new Celeris™ Arginine (ARG) column was found to be predominantly driven by electrostatics and, specifically, their charge state. Therefore, we analysed 41 compounds in liquid chromatography with ultraviolet detection to study possible relationships between the analytical retention on this phase and the pKa of the acidic solutes. Highly significant relationships were observed indicating either a linear (r2 = 0.86) or a quadratic (r2= 0.89) trend. To improve the throughput of the method, this was transferred to LC mass spectrometry, allowing the analysis of a molecule every 3 mins. The developed method was found to be fast, reliable, accurate, easily automatable and simple to set up. Finally, the analytical column's being industrially manufactured and commercially available offers broad applicability.

The Quest for Oral PROTAC drugs: Evaluating the Weaknesses of the Screening Pipeline.

A targeted bibliographic search exposed the deficiencies within existing PROTAC preclinical pipelines, including missing, poor-quality data and technical limitations in the experimental assays. Several recommendations are proposed to improve the efficiency of preclinical platforms for PROTACs.

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.

In Silico Tools to Extract the Drug Design Information Content of Degradation Data: The Case of PROTACs Targeting the Androgen Receptor.

Proteolysis-Targeting Chimeras (PROTACs) have recently emerged as a promising technology in the drug discovery landscape. Large interest in the degradation of the androgen receptor (AR) as a new anti-prostatic cancer strategy has resulted in several papers focusing on PROTACs against AR. This study explores the potential of a few in silico tools to extract drug design information from AR degradation data in the format often reported in the literature. After setting up a dataset of 92 PROTACs with consistent AR degradation values, we employed the Bemis-Murcko method for their classification. The resulting clusters were not informative in terms of structure-degradation relationship. Subsequently, we performed Degradation Cliff analysis and identified some key aspects conferring a positive contribution to activity, as well as some methodological limits when applying this approach to PROTACs. Linker structure degradation relationships were also investigated. Then, we built and characterized ternary complexes to validate previous results. Finally, we implemented machine learning classification models and showed that AR degradation for VHL-based but not CRBN-based PROTACs can be predicted from simple permeability-related 2D molecular descriptors.

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 pKa 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.

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.

Deleterious, protein-altering variants in the transcriptional coregulator ZMYM3 in 27 individuals with a neurodevelopmental delay phenotype.

Neurodevelopmental disorders (NDDs) result from highly penetrant variation in hundreds of different genes, some of which have not yet been identified. Using the MatchMaker Exchange, we assembled a cohort of 27 individuals with rare, protein-altering variation in the transcriptional coregulator ZMYM3, located on the X chromosome. Most (n = 24) individuals were males, 17 of which have a maternally inherited variant; six individuals (4 male, 2 female) harbor de novo variants. Overlapping features included developmental delay, intellectual disability, behavioral abnormalities, and a specific facial gestalt in a subset of males. Variants in almost all individuals (n = 26) are missense, including six that recurrently affect two residues. Four unrelated probands were identified with inherited variation affecting Arg441, a site at which variation has been previously seen in NDD-affected siblings, and two individuals have de novo variation resulting in p.Arg1294Cys (c.3880C>T). All variants affect evolutionarily conserved sites, and most are predicted to damage protein structure or function. ZMYM3 is relatively intolerant to variation in the general population, is widely expressed across human tissues, and encodes a component of the KDM1A-RCOR1 chromatin-modifying complex. ChIP-seq experiments on one variant, p.Arg1274Trp, indicate dramatically reduced genomic occupancy, supporting a hypomorphic effect. While we are unable to perform statistical evaluations to definitively support a causative role for variation in ZMYM3, the totality of the evidence, including 27 affected individuals, recurrent variation at two codons, overlapping phenotypic features, protein-modeling data, evolutionary constraint, and experimentally confirmed functional effects strongly support ZMYM3 as an NDD-associated gene.

Personalized Treatment for Infantile Ascending Hereditary Spastic Paralysis Based on In Silico Strategies.

Infantile onset hereditary spastic paralysis (IAHSP) is a rare neurological disease diagnosed in less than 50 children worldwide. It is transmitted with a recessive pattern and originates from mutations of the ALS2 gene, encoding for the protein alsin and involved in differentiation and maintenance of the upper motoneuron. The exact pathogenic mechanisms of IAHSP and other neurodevelopmental diseases are still largely unknown. However, previous studies revealed that, in the cytosolic compartment, alsin is present as an active tetramer, first assembled from dimer pairs. The C-terminal VPS9 domain is a key interaction site for alsin dimerization. Here, we present an innovative drug discovery strategy, which identified a drug candidate to potentially treat a patient harboring two ALS2 mutations: one truncation at lysine 1457 (not considered) and the substitution of arginine 1611 with a tryptophan (R1611W) in the C-terminus VPS9. With a protein modeling approach, we obtained a R1611W mutant model and characterized the impact of the mutation on the stability and flexibility of VPS9. Furthermore, we showed how arginine 1611 is essential for alsin's homo-dimerization and how, when mutated to tryptophan, it leads to an abnormal dimerization pattern, disrupting the formation of active tetramers. Finally, we performed a virtual screening, individuating an already therapy-approved compound (MK4) able to mask the mutant residue and re-establishing the alsin tetramers in HeLa cells. MK4 has now been approved for compassionate use.

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.

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 kwIAM (lipophilicity), EPSA, and Δ log kwIAM (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 kwIAM 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.

Application of the anthraquinone drug rhein as an axial ligand in bifunctional Pt(IV) complexes to obtain antiproliferative agents against human glioblastoma cells.

Octahedral Pt(IV) prodrugs are an effective way to combine cisplatin-like moieties and a second drug to obtain selective and stimuli responsive bifunctional antiproliferative compounds. Recently, two bifunctional Pt(IV) complexes have shown interesting in vitro and in vivo effects in glioblastoma, the most aggressive primary brain tumor. An interesting observation indicates that 4,5-dihydroxy-9,10-dioxo-9,10-dihydroanthracene-2-carboxylic acid (rhein) can inhibit in vivo glioma tumor progression. Furthermore, a prodrug in which cisplatin was combined with two molecules of rhein showed a potency higher than that of cisplatin toward cisplatin-resistant lung carcinoma cells. However, the high lipophilicity of this type of complex affects their solubility and bioavailability. To overcome these limits, in the present work, three Pt(IV) derivatives were obtained by differently linking one molecule of rhein and one acetato ligand at the axial position to a cisplatin core. The complexes proved to be similar to or more potent than the parent cisplatin and rhein, and the reference drug temozolomide on two human glioblastoma cell lines (U87-MG and T98G). They retained their activity under hypoxia and caused a significant reduction in the motility of both cell lines, which can be related to their ability to inhibit MMP2 and MMP9 matrix metalloproteinases. Finally, physicochemical and computational studies indicated that these Pt(IV) derivatives are more prone than rhein to cross the blood-brain barrier.

AI-based protein structure databases have the potential to accelerate rare diseases research: AlphaFoldDB and the case of IAHSP/Alsin.

Artificial intelligence (AI)-based protein structure databases are expected to have an impact on drug discovery. Here, we show how AlphaFold could support rare diseases research programs. We focus on Alsin, a protein responsible for rare motor neuron diseases, such as infantile-onset ascending hereditary spastic paralysis (IAHSP) and juvenile primary lateral sclerosis (JPLS), and involved in some cases of amyotrophic lateral sclerosis (ALS). First, we compared the AlphaFoldDB human Alsin model with homology models of Alsin domains. We then evaluated the flexibility profile of Alsin and of experimentally characterized mutants present in patients with IAHSP. Next, we compared preliminary models of dimeric/tetrameric Alsin responsible for its physiological action with hypothetical models reported in the literature. Finally, we suggest the best animal model for drug candidates testing. Overall, we computationally show that drug discovery efforts toward Alsin-involving diseases should be pursued.

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.

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.

Characterization of the new CelerisTM Arginine column: Retentive behaviour through a combination of chemometric tools and potential in drug analysis.

CelerisTM Arginine (ARG) is a mixed-mode stationary phase recently released on the market. To characterize its analytical behavior, the retention factors of a pool (n=100, of which 36 neutrals, 26 acids and 38 bases) of pharmaceutically relevant compounds have been measured on this phase over eight percentages (from 10 to 90% v/v) of acetonitrile (MeCN) as organic modifier. The ARG phase exhibited enhanced affinity for the molecules that are in their anionic form at the experimental pH, whilst basic compounds, albeit over a wide range of lipophilicity and pKa values, were on average poorly retained. To dissect the separation mechanism of the ARG phase, the overall analytical retention has been deconvoluted into the individual contributions of intermolecular forces by a QSPR/ Partial Least Square (PLS)/Block Relevance (BR) analysis tool recently developed by us. For the neutrals, the most relevant blocks were found to be Size, describing the interaction due to the dimension of the molecule, and O, representing the solute's hydrogen bond donor properties. The change in sign from positive to negative of the Size block, which occurs between 10% and 20% MeCN, allowed to visually appreciate the switch in the separation mode from reversed phase to normal phase. Some good statistic models for rationalizing the analytical behaviour of neutrals were developed from VS+ descriptors. However, their performance in modelling the analytical retention of acids was substandard, probably due to the intrinsic inefficacy of VS+ descriptors in handling electric charges. This instance was addressed by a complimentary MLR strategy, which led to successfully model the retention of acids on the ARG column and to shed light into their retention mechanism, which seemed to be substantially driven by electrostatics.

Permeability prediction in the beyond-Rule-of 5 chemical space: Focus on cyclic hexapeptides.

Cyclic peptides (CPs) are gaining more and more relevance in drug discovery. Since one of their main drawbacks is poor permeability, the discovery of new orally available CP drugs requires computational tools that predict CP permeability in very early drug discovery. In this study we used a literature dataset of 62 cyclic hexapeptides to evaluate the performances of a number of in silico tools based on different computational theory to model and rationalize PAMPA and Caco-2 permeability values. In particular, we submitted the dataset to a) online calculators, b) QSPR strategies, c) a physics-based tool, d) a mixed approach and e) a kinetic method. This latter is an emergent strategy in which a few relevant conformations retrieved from a set of molecular dynamics (MD) simulations by the Markov State Model (MSM) are used to establish the compounds permeability. Both free and commercial software were used. Results were compared with a model based on experimental physicochemical descriptors. All the computational approaches but online calculators performed quite well and show that lipophilicity and not polarity is the main determinant of the investigated event. A second major outcome of the study is that the impact of flexibility on the permeability of the considered dataset cannot be unambiguously assessed. Finally, our comparative analysis, which also included not common applied strategies, allowed a sound evaluation of the pros and cons of the applied computational approaches.

IκBα targeting promotes oxidative stress-dependent cell death.

BACKGROUND: Oxidative stress is a hallmark of many cancers. The increment in reactive oxygen species (ROS), resulting from an increased mitochondrial respiration, is the major cause of oxidative stress. Cell fate is known to be intricately linked to the amount of ROS produced. The direct generation of ROS is also one of the mechanisms exploited by common anticancer therapies, such as chemotherapy. METHODS: We assessed the role of NFKBIA with various approaches, including in silico analyses, RNA-silencing and xenotransplantation. Western blot analyses, immunohistochemistry and RT-qPCR were used to detect the expression of specific proteins and genes. Immunoprecipitation and pull-down experiments were used to evaluate protein-protein interactions. RESULTS: Here, by using an in silico approach, following the identification of NFKBIA (the gene encoding IκBα) amplification in various cancers, we described an inverse correlation between IκBα, oxidative metabolism, and ROS production in lung cancer. Furthermore, we showed that novel IκBα targeting compounds combined with cisplatin treatment promote an increase in ROS beyond the tolerated threshold, thus causing death by oxytosis. CONCLUSIONS: NFKBIA amplification and IκBα overexpression identify a unique cancer subtype associated with specific expression profile and metabolic signatures. Through p65-NFKB regulation, IκBα overexpression favors metabolic rewiring of cancer cells and distinct susceptibility to cisplatin. Lastly, we have developed a novel approach to disrupt IκBα/p65 interaction, restoring p65-mediated apoptotic responses to cisplatin due to mitochondria deregulation and ROS-production.

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.