ABSTRACT
Structure-based drug design, which relies on precise understanding of the target protein and its interaction with the drug candidate, is dramatically expedited by advances in computational methods for candidate prediction. Yet, the accuracy needs to be improved with more structural data from high throughput experiments, which are challenging to generate, especially for dynamic and weak associations. Herein, we applied native mass spectrometry (native MS) to rapidly characterize ligand binding of an allosteric heterodimeric complex of SARS-CoV-2 nonstructural proteins (nsp) nsp10 and nsp16 (nsp10/16), a complex essential for virus survival in the host and thus a desirable drug target. Native MS showed that the dimer is in equilibrium with monomeric states in solution. Consistent with the literature, well characterized small cosubstrate, RNA substrate, and product bind with high specificity and affinity to the dimer but not the free monomers. Unsuccessfully designed ligands bind indiscriminately to all forms. Using neutral gas collision, the nsp16 monomer with bound cosubstrate can be released from the holo dimer complex, confirming the binding to nsp16 as revealed by the crystal structure. However, we observed an unusual migration of the endogenous zinc ions bound to nsp10 to nsp16 after collisional dissociation. The metal migration can be suppressed by using surface collision with reduced precursor charge states, which presumably resulted in minimal gas-phase structural rearrangement and highlighted the importance of complementary techniques. With minimal sample input (â¼µg), native MS can rapidly detect ligand binding affinities and locations in dynamic multisubunit protein complexes, demonstrating the potential of an "all-in-one" native MS assay for rapid structural profiling of protein-to-AI-based compound systems to expedite drug discovery.
Subject(s)
Mass Spectrometry , Methyltransferases , Protein Multimerization , SARS-CoV-2 , Viral Nonstructural Proteins , Viral Regulatory and Accessory Proteins , Viral Nonstructural Proteins/chemistry , Viral Nonstructural Proteins/metabolism , SARS-CoV-2/chemistry , Mass Spectrometry/methods , Allosteric Regulation , Protein Binding , Humans , Ligands , Models, MolecularABSTRACT
Identification of potential therapeutic candidates can be expedited by integrating computational modeling with domain aware machine learning (ML) models followed by experimental validation in an iterative manner. Generative deep learning models can generate thousands of new candidates, however, their physiochemical and biochemical properties are typically not fully optimized. Using our recently developed deep learning models and a scaffold as a starting point, we generated tens of thousands of compounds for SARS-CoV-2 Mpro that preserve the core scaffold. We utilized and implemented several computational tools such as structural alert and toxicity analysis, high throughput virtual screening, ML-based 3D quantitative structure-activity relationships, multi-parameter optimization, and graph neural networks on generated candidates to predict biological activity and binding affinity in advance. As a result of these combined computational endeavors, eight promising candidates were singled out and put through experimental testing using Native Mass Spectrometry and FRET-based functional assays. Two of the tested compounds with quinazoline-2-thiol and acetylpiperidine core moieties showed IC[Formula: see text] values in the low micromolar range: [Formula: see text] [Formula: see text]M and 3.41±0.0015 [Formula: see text]M, respectively. Molecular dynamics simulations further highlight that binding of these compounds results in allosteric modulations within the chain B and the interface domains of the Mpro. Our integrated approach provides a platform for data driven lead optimization with rapid characterization and experimental validation in a closed loop that could be applied to other potential protein targets.
Subject(s)
COVID-19 , SARS-CoV-2 , Humans , Molecular Docking Simulation , Molecular Dynamics Simulation , Protease Inhibitors/pharmacology , Antiviral Agents/pharmacology , Antiviral Agents/chemistryABSTRACT
Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 µM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.