Impact of Dynamically Exposed Polarity on Permeability and Solubility of Chameleonic Drugs Beyond the Rule of 5

Structural diversity and biological relevance of benzenoid and atypical ansamycins and their congeners

Covering: 2011 to 2021The structural division of ansamycins, including those of atypical cores and different lengths of the ansa chains, is presented. Recently discovered benzenoid and atypical ansamycin scaffolds are presented in relation to their natural source and biosynthetic routes realized in bacteria as well as their muta and semisynthetic modifications influencing biological properties. To better understand the structure-activity relationships among benzenoid ansamycins structural aspects together with mechanisms of action regarding different targets in cells, are discussed. The most promising directions for structural optimizations of benzenoid ansamycins, characterized by predominant anticancer properties, were discussed in view of their potential medical and pharmaceutical applications. The bibliography of the review covers mainly years from 2011 to 2021.

Modifications, biological origin and antibacterial activity of naphthalenoid ansamycins

Covering: 2011 to 2021Structural division of natural naphthalenoid ansamycins, regarding the type of the core and length of the ansa chain, and their biosynthetic pathways in microorganisms are discussed. The great biosynthetic plasticity of natural naphthalenoid ansamycins is reflected in their structural variety due to the alterations within ansa bridge or naphthalenoid core portions. A comparison between the biological potency of natural and semisynthetic naphthalenoid ansamycins was performed and discussed in relation to the molecular targets in cells. The antibacterial potency of naphthalenoid ansamycins seems to be dependent on the ansa chain length and conformational flexibility - the higher flexibility of the ansa chain the better biological outcome is noted.

Strategies for designing proteolysis targeting chimaeras (PROTACs)

Proteolysis targeting chimaeras (PROTACs) is a cutting edge and rapidly growing technique for new drug discovery and development. Currently, the largest challenge in the molecular design and drug development of PROTACs is efficient identification of potent and drug-like degraders. This review aims to comprehensively summarize and analyse state-of-the-art methods and strategies in the design of PROTACs. We provide a detailed illustration of the general principles and tactics for designing potent PROTACs, highlight representative case studies, and discuss the advantages and limitations of these strategies. Particularly, structure-based rational PROTAC design and emerging new types of PROTACs (e.g., homo-PROTACs, multitargeting PROTACs, photo-control PROTACs and PROTAC-based conjugates) will be focused on.

Relating Conformational Equilibria to Conformer-Specific Lipophilicities: New Opportunities in Drug Discovery

Efficient drug discovery is based on a concerted effort in optimizing bioactivity and compound properties such as lipophilicity, and is guided by efficiency metrics that reflect both aspects. While conformation–activity relationships and ligand conformational control are known strategies to improve bioactivity, the use of conformer‐specific lipophilicities (logp) is much less explored. Here we show how conformer‐specific logp values can be obtained from knowledge of the macroscopic logP value, and of the equilibrium constants between the individual species in water and in octanol. This is illustrated with fluorinated amide rotamers, with integration of rotamer ¹⁹F NMR signals as a facile, direct method to obtain logp values. The difference between logp and logP optimization is highlighted, giving rise to a novel avenue for lipophilicity control in drug discovery.

Emerging investigator series: the role of chemical properties in human exposure to environmental chemicals

One of the ultimate goals of environmental exposure science is to mechanistically understand how chemical properties and human behavior interactively determine human exposure to the wide spectrum of chemicals present in the environment. This comprehensive review assembles state-of-the-art knowledge of the role of partitioning, dissociation, mass transfer, and reactive properties in human contact with and absorption of organic chemicals via oral, dermal, and respiratory routes. Existing studies have revealed that chemicals with different properties vary greatly in mass distribution and occurrence among multiple exposure media, resulting in distinct patterns of human intake from the environment. On the other hand, these chemicals encounter different levels of resistance in the passage of intestinal, dermal, and pulmonary absorption barriers and demonstrate different levels of bioavailability, due to the selectivity of biochemical, anatomical and physiological structures of these absorption barriers. Moving forward, the research community needs to gain more in-depth mechanistic insights into the complex processes in human exposure, advance the technique to better characterize and predict chemical properties, generate and leverage experimental data for a more diverse range of chemicals, and describe better the interactions between chemical properties and human behavior.

Structure-based design and synthesis of conformationally constrained derivatives of methyl-piperidinopyrazole (MPP) with estrogen receptor (ER) antagonist activity

Nuclear Estrogen receptors (ER) are cytoplasmic proteins; translocated to the nucleus to induce transcriptional signals after getting bound to the estrogen hormone. ER activation implicated in cancer cell proliferation of female reproductive organs. Thus, the discovery of ER antagonists is a reliable strategy to combat estrogen-dependent breast cancer. Endometrial carcinoma is one of the complications encountered upon long-term therapy by selective estrogen receptor modulators (SERMs) like Tamoxifen (TMX) and methyl piperidinopyrazole (MPP). Thus, the ER-full antagonist is a solution to improve the safety of this class of therapeutics during the treatment of breast cancer. We selected MPP as a lead structure to design conformationally constrained analogs. Structural rigidification is a proven strategy to transform the SERMs into full antagonists. Accordingly, we synthesized 7-methoxy-3-(4-methoxyphenyl)-4,5-dihydro-2H-benzo[g]indazoles (4), (6a-c),(8-12) along with the biphenolic counterparts(13-19)that are the anticipated active metabolites. The 4-nitrophenyl derivative(4)is with the most balanced profile regardingthe in vivoanti-uterotrophic potential (EC₅₀ = 4.160 μM); and the cytotoxicity assay of the corresponding active metabolite(13)against ER+ breast cancer cell lines (MCF-7 IC₅₀ = 7.200 μM, T-47D IC₅₀ = 11.710 μM). The inconsiderable uterotrophic activities of the elaborated ER-antagonists and weak antiproliferative activity of the compound(13)against ovarian cancer (SKOV-3 IC₅₀ = 29.800 μM) highlighted it as a good start point to elaborate potential ER-full antagonists devoid of endometrial carcinoma. Extending the pendant chain that protrudes from the 2-(4-(substituted)-phenyl) ring of the new benzo-indazoles is recommended for enhancing the potency based on the binding mode of compound(13)in the ligand-binding domain (LBD) of ER.

Crossing of the Cystic Barriers of Toxoplasma gondii by the Fluorescent Coumarin Tetra-Cyclopeptide

FR235222 is a natural tetra-cyclopeptide with a strong inhibition effect on histone deacetylases, effective on mammalian cells as well as on intracellular apicomplexan parasites, such as Toxoplasma gondii, in the tachyzoite and bradyzoite stages. This molecule is characterized by two parts: the zinc-binding group, responsible for the binding to the histone deacetylase, and the cyclic tetrapeptide moiety, which plays a crucial role in cell permeability. Recently, we have shown that the cyclic tetrapeptide coupled with a fluorescent diethyl-amino-coumarin was able to maintain properties of cellular penetration on human cells. Here, we show that this property can be extended to the crossing of the Toxoplasma gondii cystic cell wall and the cell membrane of the parasite in its bradyzoite form, while maintaining a high efficacy as a histone deacetylase inhibitor. The investigation by molecular modeling allows a better understanding of the penetration mechanism.

Discovery of Umibecestat (CNP520): A Potent, Selective, and Efficacious β-Secretase (BACE1) Inhibitor for the Prevention of Alzheimer's Disease

After identification of lead compound 6, 5-amino-1,4-oxazine BACE1 inhibitors were optimized in order to improve potency, brain penetration, and metabolic stability. Insertion of a methyl and a trifluoromethyl group at the 6-position of the 5-amino-1,4-oxazine led to 8 (NB-360), an inhibitor with a pK_a of 7.1, a very low P-glycoprotein efflux ratio, and excellent pharmacological profile, enabling high central nervous system penetration and exposure. Fur color changes observed with NB-360 in efficacy studies in preclinical animal models triggered further optimization of the series. Herein, we describe the steps leading to the discovery of 3-chloro-5-trifluoromethyl-pyridine-2-carboxylic acid [6-((3R,6R)-5-amino-3,6-dimethyl-6-trifluoromethyl-3,6-dihydro-2H-[1,4]oxazin-3-yl)-5-fluoro-pyridin-2-yl]amide 15 (CNP520, umibecestat), an inhibitor with superior BACE1/BACE2 selectivity and pharmacokinetics. CNP520 reduced significantly Aβ levels in mice and rats in acute and chronic treatment regimens without any side effects and thus qualified for Alzheimer's disease prevention studies in the clinic.

Implications of the Conformationally Flexible, Macrocyclic Structure of the First-Generation, Direct-Acting Anti-Viral Paritaprevir on Its Solid Form Complexity and Chameleonic Behavior

Direct-acting antiviral regimens have transformed therapeutic management of hepatitis C across all prevalent genotypes. Most of the chemical matter in these regimens comprises molecules well outside the traditional drug development chemical space and presents significant challenges. Herein, the implications of high conformational flexibility and the presence of a 15-membered macrocyclic ring in paritaprevir are studied through a combination of advanced computational and experimental methods with focus on molecular chameleonicity and crystal form complexity. The ability of the molecule to toggle between high and low 3D polar surface area (PSA) conformations is underpinned by intramolecular hydrogen bonding (IMHB) interactions and intramolecular steric effects. Computational studies consequently show a very significant difference of over 75 Ų in 3D PSA between polar and apolar environments and provide the structural basis for the perplexingly favorable passive permeability of the molecule. Crystal packing and protein binding resulting in strong intermolecular interactions disrupt these intramolecular interactions. Crystalline Form I benefits from strong intermolecular interactions, whereas the weaker intermolecular interactions in Form II are partially compensated by the energetic advantage of an IMHB. Like Form I, no IMHB is observed within the receptor-bound conformation; instead, an intermolecular H-bond contributes to the potency of the molecule. The choice of metastable Form II is derisked through strategies accounting for crystal surface and packing features to manage higher form specific solid-state chemical reactivity and specific processing requirements. Overall, the results show an unambiguous link between structural features and derived properties from crystallization to dissolution, permeation, and docking into the protein pocket.

JNJ-67569762, A 2-Aminotetrahydropyridine-Based Selective BACE1 Inhibitor Targeting the S3 Pocket: From Discovery to Clinical Candidate

The discovery of a novel 2-aminotetrahydropyridine class of BACE1 inhibitors is described. Their pK_a and lipophilicity were modulated by a pending sulfonyl group, while good permeability and brain penetration were achieved via intramolecular hydrogen bonding. BACE1 selectivity over BACE2 was achieved in the S3 pocket by a novel bicyclic ring system. An optimization addressing reactive metabolite formation, cardiovascular safety, and CNS toxicity is described, leading to the clinical candidate JNJ-67569762 (12), which gave robust dose-dependent BACE1-mediated amyloid β lowering without showing BACE2-dependent hair depigmentation in preclinical models. We show that 12 has a favorable projected human dose and PK and hence presented us with an opportunity to test a highly selective BACE1 inhibitor in humans. However, 12 was found to have a QT effect upon repeat dosing in dogs and its development was halted in favor of other selective leads, which will be reported in the future.

Effect of Flexibility, Lipophilicity, and the Location of Polar Residues on the Passive Membrane Permeability of a Series of Cyclic Decapeptides

Cyclic peptides have received increasing attention over the recent years as potential therapeutics for "undruggable" targets. One major obstacle is, however, their often relatively poor bioavailability. Here, we investigate the structure-permeability relationship of 24 cyclic decapeptides that share the same backbone N-methylation pattern but differ in their side chains. The peptides cover a large range of values for passive membrane permeability as well as lipophilicity and solubility. To rationalize the observed differences in permeability, we extracted for each peptide the population of the membrane-permeable conformation in water from extensive explicit-solvent molecular dynamics simulations and used this as a metric for conformational rigidity or "prefolding." The insights from the simulations together with lipophilicity measurements highlight the intricate interplay between polarity/lipophilicity and flexibility/rigidity and the possible compensating effects on permeability. The findings allow us to better understand the structure-permeability relationship of cyclic peptides and extract general guiding principles.

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.

Large-Scale Membrane Permeability Prediction of Cyclic Peptides Crossing a Lipid Bilayer Based on Enhanced Sampling Molecular Dynamics Simulations

Membrane permeability is a significant obstacle facing the development of cyclic peptide drugs. However, membrane permeation mechanisms are poorly understood. To investigate common features of permeable (and nonpermeable) designs, it is necessary to reproduce the membrane permeation process of cyclic peptides through the lipid bilayer. We simulated the membrane permeation process of 100 six-residue cyclic peptides across the lipid bilayer based on steered molecular dynamics (MD) and replica-exchange umbrella sampling simulations and predicted membrane permeability using the inhomogeneous solubility-diffusion model and a modified version of it. Furthermore, we confirmed the effectiveness of this protocol by predicting the membrane permeability of 56 eight-residue cyclic peptides with diverse chemical structures, including some confidential designs from a pharmaceutical company. As a result, a reasonable correlation between experimentally assessed and calculated membrane permeability of cyclic peptides was observed for the peptide libraries, except for strongly hydrophobic peptides. Our analysis of the MD trajectory demonstrated that most peptides were stabilized in the boundary region between bulk water and membrane and that for most peptides, the process of crossing the center of the membrane is the main obstacle to membrane permeation. The height of this barrier is well correlated with the electrostatic interaction between the peptide and the surrounding media. The structural and energetic features of the representative peptide at each vertical position within the membrane were also analyzed, revealing that peptides permeate the membrane by changing their orientation and conformation according to the surrounding environment.

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.

Global Analysis of Models for Predicting Human Absorption: QSAR, In Vitro, and Preclinical Models

Models intended to predict intestinal absorption are an essential part of the drug development process. Although many models exist for capturing intestinal absorption, many questions still exist around the applicability of these models to drug types like "beyond rule of 5" (bRo5) and low absorption compounds. This presents a challenge as current models have not been rigorously tested to understand intestinal absorption. Here, we assembled a large, structurally diverse dataset of ∼1000 compounds with known in vitro, preclinical, and human permeability and/or absorption data. In silico (quantitative structure-activity relationship), in vitro (Caco-2), and in vivo (rat) models were statistically evaluated for predictive performance against this human intestinal absorption dataset. We expect this evaluation to serve as a resource for DMPK scientists and medicinal/computational chemists to increase their understanding of permeability and absorption model utility and applications for academia and industry.

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.

Cell Permeability of Isomeric Macrocycles: Predictions and NMR Studies

Conformation-dependent 3D descriptors have been shown to provide better predictions of the physicochemical properties of macrocycles than 2D descriptors. However, the computational identification of relevant conformations for macrocycles is nontrivial. Herein, we report that the Caco-2 cell permeability difference between a pair of diastereomeric macrocycles correlated with their solvent accessible 3D polar surface area and radius of gyration. The descriptors were calculated from the macrocycles’ solution-phase conformational ensembles and independently from ensembles obtained by conformational sampling. Calculation of the two descriptors for three other stereo- and regioisomeric macrocycles also allowed the correct ranking of their cell permeability. Methods for conformational sampling may thus allow ranking of passive permeability for moderately flexible macrocycles, thereby contributing to the prioritization of macrocycles for synthesis in lead optimization.

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.

Applications of Virtual Screening in Bioprospecting: Facts, Shifts, and Perspectives to Explore the Chemo-Structural Diversity of Natural Products

Natural products are continually explored in the development of new bioactive compounds with industrial applications, attracting the attention of scientific research efforts due to their pharmacophore-like structures, pharmacokinetic properties, and unique chemical space. The systematic search for natural sources to obtain valuable molecules to develop products with commercial value and industrial purposes remains the most challenging task in bioprospecting. Virtual screening strategies have innovated the discovery of novel bioactive molecules assessing in silico large compound libraries, favoring the analysis of their chemical space, pharmacodynamics, and their pharmacokinetic properties, thus leading to the reduction of financial efforts, infrastructure, and time involved in the process of discovering new chemical entities. Herein, we discuss the computational approaches and methods developed to explore the chemo-structural diversity of natural products, focusing on the main paradigms involved in the discovery and screening of bioactive compounds from natural sources, placing particular emphasis on artificial intelligence, cheminformatics methods, and big data analyses.

Pt(iv) complexes based on cyclohexanediamines and the histone deacetylase inhibitor 2-(2-propynyl)octanoic acid: synthesis, characterization, cell penetration properties and antitumor activity

The Pt(iv) complexes based on (SP-4-2)-dichlorido(cyclohexane-1,4-diamine)platinum(ii) (kiteplatin) and the histone deacetylase inhibitor 2-(2-propynyl)octanoic acid (POA) were investigated. Since POA contains a chiral carbon, all the possible Pt(iv) isomers were prepared and characterized, and their antiproliferative activity on six cancer cell lines was compared with that of the corresponding Pt(iv) complexes containing the cyclohexane-1R,2R-diamine equatorial ligand. To justify the very good antiproliferative activity (nanomolar IC₅₀), the polarity, lipophilicity, permeability, and cell accumulation of the complexes were studied. Overall, the two series of Pt(iv) complexes showed similar cell penetration properties, being significantly better than that of the Pt(ii) reference compounds. Finally, a representative compound of the whole set of complexes (i.e., that based on cyclohexane-1R,2R-diamine and racemic POA) was tested in vivo on mice bearing Lewis lung carcinoma, showing good tumor growth inhibition with negligible body weight loss.

Predicting cell-penetrating peptides using machine learning algorithms and navigating in their chemical space

Cell-penetrating peptides (CPPs) are naturally able to cross the lipid bilayer membrane that protects cells. These peptides share common structural and physicochemical properties and show different pharmaceutical applications, among which drug delivery is the most important. Due to their ability to cross the membranes by pulling high-molecular-weight polar molecules, they are termed Trojan horses. In this study, we proposed a machine learning (ML)-based framework named BChemRF-CPPred (beyond chemical rules-based framework for CPP prediction) that uses an artificial neural network, a support vector machine, and a Gaussian process classifier to differentiate CPPs from non-CPPs, using structure- and sequence-based descriptors extracted from PDB and FASTA formats. The performance of our algorithm was evaluated by tenfold cross-validation and compared with those of previously reported prediction tools using an independent dataset. The BChemRF-CPPred satisfactorily identified CPP-like structures using natural and synthetic modified peptide libraries and also obtained better performance than those of previously reported ML-based algorithms, reaching the independent test accuracy of 90.66% (AUC = 0.9365) for PDB, and an accuracy of 86.5% (AUC = 0.9216) for FASTA input. Moreover, our analyses of the CPP chemical space demonstrated that these peptides break some molecular rules related to the prediction of permeability of therapeutic molecules in cell membranes. This is the first comprehensive analysis to predict synthetic and natural CPP structures and to evaluate their chemical space using an ML-based framework. Our algorithm is freely available for academic use at http://comptools.linc.ufpa.br/BChemRF-CPPred .

Modulation of the Passive Permeability of Semipeptidic Macrocycles: N- and C-Methylations Fine-Tune Conformation and Properties

Incorporating small modifications to peptidic macrocycles can have a major influence on their properties. For instance, N-methylation has been shown to impact permeability. A better understanding of the relationship between permeability and structure is of key importance as peptidic drugs are often associated with unfavorable pharmacokinetic profiles. Starting from a semipeptidic macrocycle backbone composed of a tripeptide tethered head-to-tail with an alkyl linker, we investigated two small changes: peptide-to-peptoid substitution and various methyl placements on the nonpeptidic linker. Implementing these changes in parallel, we created a collection of 36 compounds. Their permeability was then assessed in parallel artificial membrane permeability assay (PAMPA) and Caco-2 assays. Our results show a systematic improvement in permeability associated with one peptoid position in the cycle, while the influence of methyl substitution varies on a case-by-case basis. Using a combination of molecular dynamics simulations and NMR measurements, we offer hypotheses to explain such behavior.

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.

Connecting Hydrophobic Surfaces in Cyclic Peptides Increases Membrane Permeability

N- or C-methylation in natural and synthetic cyclic peptides can increase membrane permeability, but it remains unclear why this happens in some cases but not others. Here we compare three-dimensional structures for cyclic peptides from six families, including isomers differing only in the location of an N- or Cα-methyl substituent. We show that a single methyl group only increases membrane permeability when it connects or expands hydrophobic surface patches. Positional isomers, with the same molecular weight, hydrogen bond donors/acceptors, rotatable bonds, calculated LogP, topological polar surface area, and total hydrophobic surface area, can have different membrane permeabilities that correlate with the size of the largest continuous hydrophobic surface patch. These results illuminate a key local molecular determinant of membrane permeability.

Do Drug-likeness Rules Apply to Oral Prodrugs?

This paper describes a comparative analysis of the physicochemical and structural properties of prodrugs and their corresponding drugs with regard to drug-likeness rules. The dataset used in this work was obtained from the DrugBank. Sixty-five pairs of prodrugs/drugs were retrieved and divided into the following categories: carrier-linked to increase hydrophilic character, carrier-linked to increase absorption, and bioprecursors. We compared the physicochemical properties related to drug-likeness between prodrugs and drugs. Our results show that prodrugs do not always follow Lipinski's Rule of 5, especially as we observed 15 prodrugs with more than 10 hydrogen bond acceptors and 18 with a molecular weight greater than 500 Da. This fact highlights the importance of extending Lipinski's rules to encompass other parameters as both strategies (filtering of drug-like chemical libraries and prodrug design) aim to improve the bioavailability of compounds. Therefore, critical reasoning is fundamental to determine whether a structure has drug-like properties or could be considered a potential orally active compound in the drug-design pipeline.

Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations

Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.

Development of a Hierarchical Support Vector Regression-Based In Silico Model for Caco-2 Permeability

Drug absorption is one of the critical factors that should be taken into account in the process of drug discovery and development. The human colon carcinoma cell layer (Caco-2) model has been frequently used as a surrogate to preliminarily investigate the intestinal absorption. In this study, a quantitative structure–activity relationship (QSAR) model was generated using the innovative machine learning-based hierarchical support vector regression (HSVR) scheme to depict the exceedingly confounding passive diffusion and transporter-mediated active transport. The HSVR model displayed good agreement with the experimental values of the training samples, test samples, and outlier samples. The predictivity of HSVR was further validated by a mock test and verified by various stringent statistical criteria. Consequently, this HSVR model can be employed to forecast the Caco-2 permeability to assist drug discovery and development.

Employing complementary spectroscopies to study the conformations of an epimeric pair of side-chain stapled peptides in aqueous solution

Understanding the conformational preferences of free ligands in solution is often necessary to rationalize structure-activity relationships in drug discovery. Herein, we examine the conformational behavior of an epimeric pair of side-chain stapled peptides that inhibit the FAD dependent amine oxidase lysine specific demethylase 1 (LSD1). The peptides differ only at a single stereocenter, but display a major difference in binding affinity. Their Raman optical activity (ROA) spectra are most likely dominated by the C-terminus, obscuring the analysis of the epimeric macrocycle. By employing NMR spectroscopy, we show a difference in conformational behavior between the two compounds and that the LSD1 bound conformation of the most potent compound is present to a measurable extent in aqueous solution. In addition, we illustrate that Molecular Dynamics (MD) simulations produce ensembles that include the most important solution conformations, but that it remains problematic to identify relevant conformations with no a priori knowledge from the large conformational pool. Furthermore, this work highlights the importance of understanding the scope and limitations of the available techniques for conducting conformational analyses. It also emphasizes the importance of conformational selection of a flexible ligand in molecular recognition.

Solution Conformations Shed Light on PROTAC Cell Permeability

Proteolysis targeting chimeras (PROTACs) induce intracellular degradation of target proteins. Their bifunctional structure puts degraders in a chemical space where ADME properties often complicate drug discovery. Herein we provide the first structural insight into PROTAC cell permeability obtained by NMR studies of a VHL-based PROTAC (1), which is cell permeable despite having a high molecular weight and polarity and a large number of rotatable bonds. We found that 1 populates elongated and polar conformations in solutions that mimic extra- and intracellular compartments. Conformations were folded and had a smaller polar surface area in chloroform, mimicking a cell membrane interior. Formation of intramolecular and nonclassical hydrogen bonds, π–π interactions, and shielding of amide groups from solvent all facilitate cell permeability by minimization of size and polarity. We conclude that molecular chameleonicity appears to be of major importance for 1 to enter into target cells.

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.

Membrane Permeability in Cyclic Peptides is Modulated by Core Conformations

Cyclic peptides have the potential to bind to challenging targets, which are undruggable with small molecules, but their application is limited by low membrane permeability. Here, using a series of cyclic pentapeptides, we showed that established physicochemical criteria of permeable peptides are heavily violated. We revealed that a dominant core conformation, stabilized by amides' shielding pattern, could guide the design of novel compounds. As a result, counter-intuitive strategies, such as incorporation of polar residues, can be beneficial for permeability. We further find that core globularity is a promising descriptor, which can extend the capability of standard predictive models.

Understanding the Conformational Properties of Fluorinated Polypeptides: Molecular Modelling of Unguisin A

In this work, we investigate the conformational properties of unguisin A, a natural macrocyclic heptapeptide that incorporates a γ-aminobutyric acid (Gaba), and four of its difluorinated stereoisomers at the Gaba residue. According to nuclear magnetic resonance (NMR) experiments, their secondary structure depends dramatically on the stereochemistry of the fluorinated carbon atoms. However, many molecular details of the structure and flexibility of these systems remain unknown, so that a rationale of the conformational changes induced by the fluorine atoms in the macrocycle is still missing. To fill this gap, we apply enhanced molecular dynamics (MD) techniques to explore the peptide conformational space in dimethyl sulfoxide solution followed by 4-8 μs of conventional MD simulations that provide extensive equilibrium sampling. The simulations, which compare reasonably well with the NMR-based observations, show that the secondary structure of the macrocycle is altered substantially upon fluorination, except for the (S,S) diastereomer. It also turns out that the conformations of the fluorinated peptides are visited during the enhanced MD simulation of natural unguisin A, suggesting thus that conformations accessible to the unsubstituted macrocyclic peptide may be selected by fluorination. Therefore, computational characterization of the macrocyclic peptides could be helpful in the rational design of stereoselective fluorinated peptides with fine-tuned conformation and activity.

Benchmark of Generic Shapes for Macrocycles

Macrocycles target proteins that are otherwise considered undruggable because of a lack of hydrophobic cavities and the presence of extended featureless surfaces. Increasing efforts by computational chemists have developed effective software to overcome the restrictions of torsional and conformational freedom that arise as a consequence of macrocyclization. Moloc is an efficient algorithm, with an emphasis on high interactivity, and has been constantly updated since 1986 by drug designers and crystallographers of the Roche biostructural community. In this work, we have benchmarked the shape-guided algorithm using a dataset of 208 macrocycles, carefully selected on the basis of structural complexity. We have quantified the accuracy, diversity, speed, exhaustiveness, and sampling efficiency in an automated fashion and we compared them with four commercial (Prime, MacroModel, molecular operating environment, and molecular dynamics) and four open-access (experimental-torsion distance geometry with additional “basic knowledge” alone and with Merck molecular force field minimization or universal force field minimization, Cambridge Crystallographic Data Centre conformer generator, and conformator) packages. With three-quarters of the database processed below the threshold of high ring accuracy, Moloc was identified as having the highest sampling efficiency and exhaustiveness without producing thousands of conformations, random ring splitting into two half-loops, and possibility to interactively produce globular or flat conformations with diversity similar to Prime, MacroModel, and molecular dynamics. The algorithm and the Python scripts for full automatization of these parameters are freely available for academic use.

First total synthesis of versicotide A, B and C

The syntheses of versicotides A-C, natural products containing anthranilic acid and NMe-Ala, were achieved by solid phase peptide synthesis on 2-chlorotrityl resin followed by solution phase macrocyclization. Using an oxyma additive, the difficult coupling reactions to the deactivated aromatic amine of o-aminobenzoic acid, were performed in high yield, avoiding anthranilic rearrangements or side reactions.

Drug-Like Properties in Macrocycles above MW 1000: Backbone Rigidity versus Side-Chain Lipophilicity

Large macrocyclic peptides can achieve surprisingly high membrane permeability, although the properties that govern permeability in this chemical space are only beginning to come into focus. We generated two libraries of cyclic decapeptides with stable cross-β conformations, and found that peptoid substitutions within the β-turns of the macrocycle preserved the rigidity of the parent scaffold, whereas peptoid substitutions in the opposing β-strands led to "chameleonic" species that were rigid in nonpolar media but highly flexible in water. Both rigid and chameleonic compounds showed high permeability over a wide lipophilicity range, with peak permeabilities differing significantly depending on scaffold rigidity. Our findings indicate that modulating lipophilicity can be used to engineer favorable ADME properties into both rigid and flexible macrocyclic peptides, and that scaffold rigidity can be used to tune optimal lipophilicity.

Synergistic Computational Modeling Approaches as Team Players in the Game of Solubility Predictions

Several approaches to predict and model drug solubility have been used in the drug discovery and development processes during the last decades. Each of these approaches have their own benefits and place, and are typically used as standalone approaches rather than in concert. The synergistic effects of these are often overlooked, partly due to the need of computational experts to perform the modeling and simulations as well as analyzing the data obtained. Here we provide our views on how these different approaches can be used to retrieve more information on drug solubility, ranging from multivariate data analysis over thermodynamic cycle modeling to molecular dynamics simulations. We are discussing aqueous solubility as well as solubility in more complex mixed solvents and media with colloidal structures present. We conclude that the field of computational pharmaceutics is in its early days but with a bright future ahead. However, education of computational formulators with broad knowledge of modeling and simulation approaches is imperative if computational pharmaceutics is to reach its full potential.

Cyclic Peptides as Drugs for Intracellular Targets: The Next Frontier in Peptide Therapeutic Development

Developing macrocyclic peptides that can reach intracellular targets is a significant challenge. This review discusses the most recent strategies used to develop cell permeable cyclic peptides that maintain binding to their biological target inside the cell. Macrocyclic peptides are unique from small molecules because traditional calculated physical properties are unsuccessful for predicting cell membrane permeability. Peptide synthesis and experimental membrane permeability is the only strategy that effectively differentiates between cell permeable and cell impermeable molecules. Discussed are chemical strategies, including backbone N-methylation and stereochemical changes, which have produced molecular scaffolds with improved cell permeability. However, these improvements often come at the expense of biological activity as chemical modifications alter the peptide conformation, frequently impacting the compound's ability to bind to the target. Highlighted is the most promising approach, which involves side-chain alterations that improve cell permeability without impact binding events.

Predicting the Permeability of Macrocycles from Conformational Sampling - Limitations of Molecular Flexibility

Macrocycles constitute superior ligands for targets that have flat binding sites but often require long synthetic routes, emphasizing the need for property prediction prior to synthesis. We have investigated the scope and limitations of machine learning classification models and of regression models for predicting the cell permeability of a set of denovo-designed, drug-like macrocycles. 2D-Based classification models, which are fast to calculate, discriminated between macrocycles that had low-medium and high permeability and may be used as virtual filters in early drug discovery projects. Importantly, stereo- and regioisomer were correctly classified. QSPR studies of two small sets of comparator drugs suggested that use of 3D descriptors, calculated from biologically relevant conformations, would allow development of more precise regression models for late phase drug projects. However, a 3D permeability model could only be developed for a rigid series of macrocycles. Comparison of NMR based conformational analysis with in silico conformational sampling indicated that this shortcoming originates from the inability of the molecular mechanics force field to identify the relevant conformations for flexible macrocycles. We speculate that a Kier flexibility index of ≤10 constitutes a current upper limit for reasonably accurate 3D prediction of macrocycle cell permeability.

The Critical Role of Passive Permeability in Designing Successful Drugs

Passive permeability is a key property in drug disposition and delivery. It is critical for gastrointestinal absorption, brain penetration, renal reabsorption, defining clearance mechanisms and drug-drug interactions. Passive diffusion rate is translatable across tissues and animal species, while the extent of absorption is dependent on drug properties, as well as in vivo physiology/pathophysiology. Design principles have been developed to guide medicinal chemistry to enhance absorption, which combine the balance of aqueous solubility, permeability and the sometimes unfavorable compound characteristic demanded by the target. Permeability assays have been implemented that enable rapid development of structure-permeability relationships for absorption improvement. Future advances in assay development to reduce nonspecific binding and improve mass balance will enable more accurately measurement of passive permeability. Design principles that integrate potency, selectivity, passive permeability and other ADMET properties facilitate rapid advancement of successful drug candidates to patients.

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.