Ultralarge Virtual Screening Identifies SARS-CoV-2 Main Protease Inhibitors with Broad-Spectrum Activity against Coronaviruses

Drugs targeting SARS-CoV-2 could have saved millions of lives during the COVID-19 pandemic, and it is now crucial to develop inhibitors of coronavirus replication in preparation for future outbreaks. We explored two virtual screening strategies to find inhibitors of the SARS-CoV-2 main protease in ultralarge chemical libraries. First, structure-based docking was used to screen a diverse library of 235 million virtual compounds against the active site. One hundred top-ranked compounds were tested in binding and enzymatic assays. Second, a fragment discovered by crystallographic screening was optimized guided by docking of millions of elaborated molecules and experimental testing of 93 compounds. Three inhibitors were identified in the first library screen, and five of the selected fragment elaborations showed inhibitory effects. Crystal structures of target-inhibitor complexes confirmed docking predictions and guided hit-to-lead optimization, resulting in a noncovalent main protease inhibitor with nanomolar affinity, a promising in vitro pharmacokinetic profile, and broad-spectrum antiviral effect in infected cells.

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.

Linker-Dependent Folding Rationalizes PROTAC Cell Permeability

Proteolysis-targeting chimeras (PROTACs) must be cell permeable to reach their target proteins. This is challenging as the bivalent structure of PROTACs puts them in chemical space at, or beyond, the outer limits of oral druggable space. We used NMR spectroscopy and molecular dynamics (MD) simulations independently to gain insights into the origin of the differences in cell permeability displayed by three flexible cereblon PROTACs having closely related structures. Both methods revealed that the propensity of the PROTACs to adopt folded conformations with a low solvent-accessible 3D polar surface area in an apolar environment is correlated to high cell permeability. The chemical nature and the flexibility of the linker were essential for the PROTACs to populate folded conformations stabilized by intramolecular hydrogen bonds, π-π interactions, and van der Waals interactions. We conclude that MD simulations may be used for the prospective ranking of cell permeability in the design of cereblon PROTACs.

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.

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.

Going Viral: An Investigation into the Chameleonic Behaviour of Antiviral Compounds

The ability to adjust conformations in response to the polarity of the environment, i.e. molecular chameleonicity, is considered to be important for conferring both high aqueous solubility and high cell permeability to drugs in chemical space beyond Lipinski's rule of 5. We determined the conformational ensembles populated by the antiviral drugs asunaprevir, simeprevir, atazanavir and daclatasvir in polar (DMSO-d₆ ) and non-polar (chloroform) environments with NMR spectroscopy. Daclatasvir was fairly rigid, whereas the first three showed large flexibility in both environments, that translated into major differences in solvent accessible 3D polar surface area within each conformational ensemble. No significant differences in size and polar surface area were observed between the DMSO-d₆ and chloroform ensembles of these three drugs. We propose that such flexible compounds are characterized as "partial molecular chameleons" and hypothesize that their ability to adopt conformations with low polar surface area contributes to their membrane permeability and oral absorption.

Conformation of the Macrocyclic Drug Lorlatinib in Polar and Nonpolar Environments: A MD Simulation and NMR Study

The replica exchange molecular dynamics (REMD) simulation is demonstrated to readily predict the conformations of the macrocyclic drug lorlatinib, as validated by solution NMR studies. In aqueous solution, lorlatinib adopts a conformer identical to its target bound structure. This conformer is stabilized by an extensive hydrogen bond network to the solvents. In chloroform, lorlatinib populates two conformers with the second one being less polar, which may contribute to lorlatinib's ability to cross cell membranes.

Antibacterial and cytotoxic prenylated dihydrochalcones from Eriosema montanum

Two new prenylated dihydrochalcones (1,2) and eighteen known secondary metabolites (3-20) were isolated from the CH₂Cl₂-MeOH (1:1) extracts of the roots, the stem bark and the leaves of Eriosema montanum Baker f. (Leguminosae). The structures of the isolated compounds were characterized by NMR, IR, and UV spectroscopic and mass spectrometric analyses. The structures of compounds 5, 10, 11 and 13 were confirmed by single crystal X-ray diffraction. The antibacterial activity of the crude extracts and the isolated constituents were established against Gram-positive and Gram-negative bacteria. Among the tested compounds, 1-4 and 10 showed strong activity against the Gram-positive bacterium Bacillus subtilis with minimum inhibitory concentration (MIC) ranging from 3.1 to 8.9 μM, as did the leaf crude extract with an MIC of 3.0 μg/mL. None of the crude extracts nor the isolated compounds were active against Escherichia coli. Compounds 1, 3 and 4 showed higher cytotoxicity, evaluated against the human breast cancer cell line MCF-7, with EC₅₀ of 7.0, 18.0 and 18.0 μM, respectively. These findings contribute to the phytochemical understanding of the genus Eriosema, and highlight the pharmaceutical potential of prenylated dihydrochalcones.

Prenylated Flavonoids from the Roots of Tephrosia rhodesica

Five new compounds-rhodimer (1), rhodiflavan A (2), rhodiflavan B (3), rhodiflavan C (4), and rhodacarpin (5)-along with 16 known secondary metabolites, were isolated from the CH₂Cl₂-CH₃OH (1:1) extract of the roots of Tephrosia rhodesica. They were identified by NMR spectroscopic, mass spectrometric, X-ray crystallographic, and ECD spectroscopic analyses. The crude extract and the isolated compounds 2-5, 9, 15, and 21 showed activity (100% at 10 μg and IC₅₀ = 5-15 μM) against the chloroquine-sensitive (3D7) strain of Plasmodium falciparum.

Biflavanones, Chalconoids, and Flavonoid Analogues from the Stem Bark of Ochna holstii

Two new biflavanones (1 and 2), three new bichalconoids (3-5), and 11 known flavonoid analogues (6-16) were isolated from the stem bark extract (CH3OH-CH2Cl2, 7:3, v/v) of Ochna holstii. The structures of the isolated metabolites were elucidated by NMR spectroscopic and mass spectrometric analyses. The crude extract and the isolated metabolites were evaluated for antibacterial activity against Bacillus subtilis (Gram-positive) and Escherichia coli (Gram-negative) as well as for cytotoxicity against the MCF-7 human breast cancer cell line. The crude extract and holstiinone A (1) exhibited moderate antibacterial activity against B. subtilis with MIC values of 9.1 μg/mL and 14 μM, respectively. The crude extract and lophirone F (14) showed cytotoxicity against MCF-7 with EC50 values of 11 μg/mL and 24 μM, respectively. The other isolated metabolites showed no significant antibacterial activities (MIC > 250 μM) and cytotoxicities (EC50 ≥ 350 μM).

Antiplasmodial prenylated flavanonols from Tephrosia subtriflora

The CH₂Cl₂/MeOH (1:1) extract of the aerial parts of Tephrosia subtriflora afforded a new flavanonol, named subtriflavanonol (1), along with the known flavanone spinoflavanone B, and the known flavanonols MS-II (2) and mundulinol. The structures were elucidated by the use of NMR spectroscopy and mass spectrometry. The absolute configuration of the flavanonols was determined based on quantum chemical ECD calculations. In the antiplasmodial assay, compound 2 showed the highest activity against chloroquine-sensitive Plasmodium falciparum reference clones (D6 and 3D7), artemisinin-sensitive isolate (F32-TEM) as well as field isolate (KSM 009) with IC₅₀ values 1.4-4.6 μM without significant cytotoxicity against Vero and HEp2 cell lines (IC₅₀ > 100 μM). The new compound (1) showed weak antiplasmodial activity, IC₅₀ 12.5-24.2 μM, but also showed selective anticancer activity against HEp2 cell line (CC₅₀ 16.9 μM).

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.

Modified ent-Abietane Diterpenoids from the Leaves of Suregada zanzibariensis

The leaf extract of Suregada zanzibariensis gave two new modified ent-abietane diterpenoids, zanzibariolides A (1) and B (2), and two known triterpenoids, simiarenol (3) and β-amyrin (4). The structures of the isolated compounds were elucidated based on NMR and MS data analysis. Single-crystal X-ray diffraction was used to establish the absolute configurations of compounds 1 and 2. The crude leaf extract inhibited the infectivity of herpes simplex virus 2 (HSV-2, IC50 11.5 μg/mL) and showed toxicity on African green monkey kidney (GMK AH1) cells at CC50 52 μg/mL. The isolated compounds 1-3 showed no anti-HSV-2 activity and exhibited insignificant toxicity against GMK AH1 cells at ≥100 μM.

Antibacterial and cytotoxic biflavonoids from the root bark of Ochna kirkii

The new isoflavonoid kirkinone A (1) and biflavonoid kirkinone B (2) along with six known compounds (3-8) were isolated from the methanolic extract of the root bark of Ochna kirkii. The compounds were identified by NMR spectroscopic and mass spectrometric analyses. Out of the eight isolated natural products, calodenin B (4) and lophirone A (6) showed significant antibacterial activity against the Gram-positive bacterium Bacillus subtilis with MIC values of 2.2 and 28 μM, and cytotoxicity against the MCF-7 human breast cancer cell line with EC₅₀ values of 219.3 and 19.2 μM, respectively. The methanolic crude extract of the root bark exhibited cytotoxicity at EC₅₀ 8.4 μg/mL. The isolated secondary metabolites and the crude extract were generally inactive against the Gram-negative Escherichia coli (MIC ≥400 μg/mL). Isolation of biflavonoids and related secondary metabolites from O. kirkii demonstrates their chemotaxonomic significance to the genus Ochna and to other members of the family Ochnaceae.

Antiviral Rotenoids and Isoflavones Isolated from Millettia oblata ssp. teitensis

Three new (1-3) and six known rotenoids (5-10), along with three known isoflavones (11-13), were isolated from the leaves of Millettia oblata ssp. teitensis. A new glycosylated isoflavone (4), four known isoflavones (14-18), and one known chalcone (19) were isolated from the root wood extract of the same plant. The structures were elucidated by NMR and mass spectrometric analyses. The absolute configuration of the chiral compounds was established by a comparison of experimental ECD and VCD data with those calculated for the possible stereoisomers. This is the first report on the use of VCD to assign the absolute configuration of rotenoids. The crude leaves and root wood extracts displayed anti-RSV (human respiratory syncytial virus) activity with IC50 values of 0.7 and 3.4 μg/mL, respectively. Compounds 6, 8, 10, 11, and 14 showed anti-RSV activity with IC50 values of 0.4-10 μM, while compound 3 exhibited anti-HRV-2 (human rhinovirus 2) activity with an IC50 of 4.2 μM. Most of the compounds showed low cytotoxicity for laryngeal carcinoma (HEp-2) cells; however compounds 3, 11, and 14 exhibited low cytotoxicity also in primary lung fibroblasts. This is the first report on rotenoids showing antiviral activity against RSV and HRV viruses.

Impact of Linker Composition on VHL PROTAC Cell Permeability

The discovery of cell permeable and orally bioavailable von Hippel-Lindau (VHL) proteolysis targeting chimeras (PROTACs) is challenging as their structures locates them at, or beyond, the outer limits of oral druggable space. We have designed a set of nine VHL PROTACs and found that the linker had a profound impact on passive cell permeability. Determination of the solution ensembles in a nonpolar solvent revealed that high permeability was correlated to the ability of the PROTACs to adopt folded conformations that have a low solvent accessible 3D polar surface area. Our results suggest that the design of cell permeable VHL PROTACs could focus on linkers that facilitate shielding of polar surface area in the VHL ligand in a nonpolar but not in a polar environment. In addition, we found that not only intramolecular hydrogen bonds, but also NH-π and π-π interactions contribute to the stabilization of low-polarity conformations, and thereby to high cell permeability.

Antiviral iridoid glycosides from Clerodendrum myricoides

The methanol root extract of Clerodendrum myricoides (Hochst.) Vatke afforded two new (1, 2) and two known (3, 4) iridoid glycosides. The structures of the isolated compounds were established based on NMR, IR, UV and MS data analyses. The crude extract and the isolated constituents were assayed for antiviral activity against the human respiratory syncytial virus (RSV) in human laryngeal epidermoid carcinoma (HEp-2) cells. The crude extract inhibited RSV infectivity at EC₅₀ = 0.21 μg/ml, while it showed cytotoxicity against HEp-2 cells with CC₅₀ = 9 μg/ml. Compound 2 showed 43.2% virus inhibition at 100 μM, while compounds 1 as well as 3 and 4 had only weak antiviral and cytotoxic activities.

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.