Drug Discovery in Low Data Regimes: Leveraging a Computational Pipeline for the Discovery of Novel SARS-CoV-2 Nsp14-MTase Inhibitors.

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has led to significant global morbidity and mortality. A crucial viral protein, the non-structural protein 14 (nsp14), catalyzes the methylation of viral RNA and plays a critical role in viral genome replication and transcription. Due to the low mutation rate in the nsp region among various SARS-CoV-2 variants, nsp14 has emerged as a promising therapeutic target. However, discovering potential inhibitors remains a challenge. In this work, we introduce a computational pipeline for the rapid and efficient identification of potential nsp14 inhibitors by leveraging virtual screening and the NCI open compound collection, which contains 250,000 freely available molecules for researchers worldwide. The introduced pipeline provides a cost-effective and efficient approach for early-stage drug discovery by allowing researchers to evaluate promising molecules without incurring synthesis expenses. Our pipeline successfully identified seven promising candidates after experimentally validating only 40 compounds. Notably, we discovered NSC620333, a compound that exhibits a strong binding affinity to nsp14 with a dissociation constant of 427 ± 84 nM. In addition, we gained new insights into the structure and function of this protein through molecular dynamics simulations. We identified new conformational states of the protein and determined that residues Phe367, Tyr368, and Gln354 within the binding pocket serve as stabilizing residues for novel ligand interactions. We also found that metal coordination complexes are crucial for the overall function of the binding pocket. Lastly, we present the solved crystal structure of the nsp14-MTase complexed with SS148 (PDB:8BWU), a potent inhibitor of methyltransferase activity at the nanomolar level (IC50 value of 70 ± 6 nM). Our computational pipeline accurately predicted the binding pose of SS148, demonstrating its effectiveness and potential in accelerating drug discovery efforts against SARS-CoV-2 and other emerging viruses.

Expanded ensemble predictions of absolute binding free energies in the SAMPL9 host-guest challenge.

As part of the SAMPL9 community-wide blind host-guest challenge, we implemented an expanded ensemble workflow to predict absolute binding free energies for 13 small molecules against pillar[6]arene. Notable features of our protocol include consideration of a variety of protonation and enantiomeric states for both host and guests, optimization of alchemical intermediates, and analysis of free energy estimates and their uncertainty using large numbers of simulation replicates performed using distributed computing. Our predictions of absolute binding free energies resulted in a mean absolute error of 2.29 kcal mol-1 and an R2 of 0.54. Overall, results show that expanded ensemble calculations using all-atom molecular dynamics simulations are a valuable and efficient computational tool in predicting absolute binding free energies.

Metal Cation-Binding Mechanisms of Q-Proline Peptoid Macrocycles in Solution.

The rational design of foldable and functionalizable peptidomimetic scaffolds requires the concerted application of both computational and experimental methods. Recently, a new class of designed peptoid macrocycle incorporating spiroligomer proline mimics (Q-prolines) has been found to preorganize when bound by monovalent metal cations. To determine the solution-state structure of these cation-bound macrocycles, we employ a Bayesian inference method (BICePs) to reconcile enhanced-sampling molecular simulations with sparse ROESY correlations from experimental NMR studies to predict and design conformational and binding properties of macrocycles as functional scaffolds for peptidomimetics. Conformations predicted to be most populated in solution were then simulated in the presence of explicit cations to yield trajectories with observed binding events, revealing a highly preorganized all-trans amide conformation, whose formation is likely limited by the slow rate of cis/trans isomerization. Interestingly, this conformation differs from a racemic crystal structure solved in the absence of cation. Free energies of cation binding computed from distance-dependent potentials of mean force suggest Na+ has a higher affinity to the macrocycle than K+, with both cations binding much more strongly in acetonitrile than water. The simulated affinities are able to correctly rank the extent to which different macrocycle sequences exhibit preorganization in the presence of different metal cations and solvents, suggesting our approach is suitable for solution-state computational design.

Assigning confidence to molecular property prediction.

Introduction: Computational modeling has rapidly advanced over the last decades. Recently, machine learning has emerged as a powerful and cost-effective strategy to learn from existing datasets and perform predictions on unseen molecules. Accordingly, the explosive rise of data-driven techniques raises an important question: What confidence can be assigned to molecular property predictions and what techniques can be used?Areas covered: The authors discuss popular strategies for predicting molecular properties, their corresponding uncertainty sources and methods to quantify uncertainty. First, the authors' considerations for assessing confidence begin with dataset bias and size, data-driven property prediction and feature design. Next, the authors discuss property simulation via computations of binding affinity in detail. Lastly, they investigate how these uncertainties propagate to generative models, as they are usually coupled with property predictors.Expert opinion: Computational techniques are paramount to reduce the prohibitive cost of brute-force experimentation during exploration. The authors believe that assessing uncertainty in property prediction models is essential whenever closed-loop drug design campaigns relying on high-throughput virtual screening are deployed. Accordingly, considering sources of uncertainty leads to better-informed validations, more reliable predictions and more realistic expectations of the entire workflow. Overall, this increases confidence in the predictions and, ultimately, accelerates drug design.

SARS-CoV-2 simulations go exascale to predict dramatic spike opening and cryptic pockets across the proteome.

SARS-CoV-2 has intricate mechanisms for initiating infection, immune evasion/suppression and replication that depend on the structure and dynamics of its constituent proteins. Many protein structures have been solved, but far less is known about their relevant conformational changes. To address this challenge, over a million citizen scientists banded together through the Folding@home distributed computing project to create the first exascale computer and simulate 0.1 seconds of the viral proteome. Our adaptive sampling simulations predict dramatic opening of the apo spike complex, far beyond that seen experimentally, explaining and predicting the existence of 'cryptic' epitopes. Different spike variants modulate the probabilities of open versus closed structures, balancing receptor binding and immune evasion. We also discover dramatic conformational changes across the proteome, which reveal over 50 'cryptic' pockets that expand targeting options for the design of antivirals. All data and models are freely available online, providing a quantitative structural atlas.

Metal-Binding Q-Proline Macrocycles.

We introduce the efficient Fmoc-SPPS and peptoid synthesis of Q-proline-based, metal-binding macrocycles (QPMs), which bind metal cations and display nine functional groups. Metal-free QPMs are disordered, evidenced by NMR and a crystal structure of QPM-3 obtained through racemic crystallization. Upon addition of metal cations, QPMs adopt ordered structures. Notably, the addition of a second functional group at the hydantoin amide position (R2) converts the proline ring from Cγ-endo to Cγ-exo, due to steric interactions.

SARS-CoV-2 Simulations Go Exascale to Capture Spike Opening and Reveal Cryptic Pockets Across the Proteome.

SARS-CoV-2 has intricate mechanisms for initiating infection, immune evasion/suppression, and replication, which depend on the structure and dynamics of its constituent proteins. Many protein structures have been solved, but far less is known about their relevant conformational changes. To address this challenge, over a million citizen scientists banded together through the Folding@home distributed computing project to create the first exascale computer and simulate an unprecedented 0.1 seconds of the viral proteome. Our simulations capture dramatic opening of the apo Spike complex, far beyond that seen experimentally, which explains and successfully predicts the existence of 'cryptic' epitopes. Different Spike homologues modulate the probabilities of open versus closed structures, balancing receptor binding and immune evasion. We also observe dramatic conformational changes across the proteome, which reveal over 50 'cryptic' pockets that expand targeting options for the design of antivirals. All data and models are freely available online, providing a quantitative structural atlas.

Fluorinated Aromatic Monomers as Building Blocks To Control α-Peptoid Conformation and Structure.

Peptoids are peptidomimetics of interest in the fields of drug development and biomaterials. However, obtaining stable secondary structures is challenging, and designing these requires effective control of the peptoid tertiary amide cis/trans equilibrium. Herein, we report new fluorine-containing aromatic monomers that can control peptoid conformation. Specifically, we demonstrate that a fluoro-pyridine group can be used to circumvent the need for monomer chirality to control the cis/trans equilibrium. We also show that incorporation of a trifluoro-methyl group ( NCF3Rpe) rather than a methyl group ( NRpe) at the α-carbon of a monomer gives rise to a 5-fold increase in cis-isomer preference.