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.

Resistance gene-guided genome mining reveals the roseopurpurins as inhibitors of cyclin-dependent kinases.

With the significant increase in the availability of microbial genome sequences in recent years, resistance gene-guided genome mining has emerged as a powerful approach for identifying natural products with specific bioactivities. Here, we present the use of this approach to reveal the roseopurpurins as potent inhibitors of cyclin-dependent kinases (CDKs), a class of cell cycle regulators implicated in multiple cancers. We identified a biosynthetic gene cluster (BGC) with a putative resistance gene with homology to human CDK2. Using targeted gene disruption and transcription factor overexpression in Aspergillus uvarum, and heterologous expression of the BGC in Aspergillus nidulans, we demonstrated that roseopurpurin C (1) is produced by this cluster and characterized its biosynthesis. We determined the potency, specificity, and mechanism of action of 1 as well as multiple intermediates and shunt products produced from the BGC. We show that 1 inhibits human CDK2 with a Kiapp of 44 nM, demonstrates selectivity for clinically relevant members of the CDK family, and induces G1 cell cycle arrest in HCT116 cells. Structural analysis of 1 complexed with CDK2 revealed the molecular basis of ATP-competitive inhibition.

An Open Drug Discovery Competition: Experimental Validation of Predictive Models in a Series of Novel Antimalarials.

The Open Source Malaria (OSM) consortium is developing compounds that kill the human malaria parasite, Plasmodium falciparum, by targeting PfATP4, an essential ion pump on the parasite surface. The structure of PfATP4 has not been determined. Here, we describe a public competition created to develop a predictive model for the identification of PfATP4 inhibitors, thereby reducing project costs associated with the synthesis of inactive compounds. Competition participants could see all entries as they were submitted. In the final round, featuring private sector entrants specializing in machine learning methods, the best-performing models were used to predict novel inhibitors, of which several were synthesized and evaluated against the parasite. Half possessed biological activity, with one featuring a motif that the human chemists familiar with this series would have dismissed as "ill-advised". Since all data and participant interactions remain in the public domain, this research project "lives" and may be improved by others.

Locating ligand binding sites in G-protein coupled receptors using combined information from docking and sequence conservation.

GPCRs (G-protein coupled receptors) are the largest family of drug targets and share a conserved structure. Binding sites are unknown for many important GPCR ligands due to the difficulties of GPCR recombinant expression, biochemistry, and crystallography. We describe our approach, ConDockSite, for predicting ligand binding sites in class A GPCRs using combined information from surface conservation and docking, starting from crystal structures or homology models. We demonstrate the effectiveness of ConDockSite on crystallized class A GPCRs such as the beta2 adrenergic and A2A adenosine receptors. We also demonstrate that ConDockSite successfully predicts ligand binding sites from high-quality homology models. Finally, we apply ConDockSite to predict the ligand binding sites on a structurally uncharacterized GPCR, GPER, the G-protein coupled estrogen receptor. Most of the sites predicted by ConDockSite match those found in other independent modeling studies. ConDockSite predicts that four ligands bind to a common location on GPER at a site deep in the receptor cleft. Incorporating sequence conservation information in ConDockSite overcomes errors introduced from physics-based scoring functions and homology modeling.

Mechanisms of allosteric and mixed mode aromatase inhibitors.

Aromatase (CYP19) catalyzes the last biosynthetic step of estrogens in mammals and is a primary drug target for hormone-related breast cancer. However, treatment with aromatase inhibitors is often associated with adverse effects and drug resistance. In this study, we used virtual screening targeting a predicted cytochrome P450 reductase binding site on aromatase to discover four novel non-steroidal aromatase inhibitors. The inhibitors have potencies comparable to the noncompetitive tamoxifen metabolite, endoxifen. Our two most potent inhibitors, AR11 and AR13, exhibit both mixed-type and competitive-type inhibition. The cytochrome P450 reductase-CYP19 coupling interface likely acts as a transient binding site. Our modeling shows that our inhibitors bind better at different sites near the catalytic site. Our results predict the location of multiple ligand binding sites on aromatase. The combination of modeling and experimental results supports the important role of the reductase binding interface as a low affinity, promiscuous ligand binding site. Our new inhibitors may be useful as alternative chemical scaffolds that may show different adverse effects profiles than current clinically used aromatase inhibitors.

Structural studies of full-length receptor tyrosine kinases and their implications for drug design.

Receptor tyrosine kinases (RTKs) are important drug targets for cancer and immunological disorders. Crystal structures of individual RTK domains have contributed greatly to the structure-based drug design of clinically used drugs. Low-resolution structures from electron microscopy are now available for the RTKs, EGFR, PDGFR, and Kit. However, there are still no high-resolution structures of full-length RTKs due to the technical challenges of working with these complex, membrane proteins. Here, we review what has been learned from structural studies of these three RTKs regarding their mechanisms of ligand binding, activation, oligomerization, and inhibition. We discuss the implications for drug design. More structural data on full-length RTKs may facilitate the discovery of druggable sites and drugs with improved specificity and effectiveness against resistant mutants.

RORγ Structural Plasticity and Druggability.

Retinoic acid receptor-related orphan receptor γ (RORγ) is a transcription factor regulating the expression of the pro-inflammatory cytokine IL-17 in human T helper 17 (Th17) cells. Activating RORγ can induce multiple IL-17-mediated autoimmune diseases but may also be useful for anticancer therapy. Its deep immunological functions make RORɣ an attractive drug target. Over 100 crystal structures have been published describing atomic interactions between RORɣ and agonists and inverse agonists. In this review, we focus on the role of dynamic properties and plasticity of the RORɣ orthosteric and allosteric binding sites by examining structural information from crystal structures and simulated models. We discuss the possible influences of allosteric ligands on the orthosteric binding site. We find that high structural plasticity favors the druggability of RORɣ, especially for allosteric ligands.

Decoding allosteric communication pathways in protein lysine acetyltransferase.

In bacteria, protein lysine acetylation circuits can control core processes such as carbon metabolism. In E. coli, cyclic adenosine monophosphate (cAMP) controls the transcription level and activity of protein lysine acetyltransferase (PAT). The M. tuberculosis PatA (Mt-PatA) resides in two different conformations; the activated state and autoinhibited state. However, the mechanism of cAMP allosteric regulation of Mt-PatA remains mysterious. Here, we performed extensive all-atom molecular dynamics (MD) simulations (three independent run for each system) and built a residue-residue dynamic correlation network to show how cAMP mediates allosteric activation. cAMP binds at the regulatory site in the regulatory domain, which is 32 Å away from the catalytic site. An extensive conformational restructuring relieves autoinhibition caused by a molecular Lid (residues 161-203) that shelters the substrate-binding surface. In the activated state, the regulatory domain rotates (~40°) around Ser144, which links both domains. Rotation removes the C-terminus from the cAMP site and relieves the autoinhibited state. Also, the molecular Lid refolds and creates an activator binding site. A conserved residue, His173, was mutated into Lys in the Lid, and during an MD trajectory of the activated state, positioned itself near an acetyl donor molecule in the catalytic domain, suggesting a direct mechanism for acetylation. This study describes the allosteric framework for Mt-PatA and prerequisite intermediate states that permit long-distance signal transmission.

In-silico design of peptide inhibitors of K-Ras target in cancer disease.

Cancer is a leading cause of death, over one million individuals analyzed, and around 500,000 deaths happen due to cancer every year alone in the United States. The Ras is a significant protein in the signaling transduction pathways and has a leading role in cell proliferation. Above 30% of all human tumors arises due to the mutations in genes that encode a Ras protein that operate signaling cascades necessary for malignant transformation, tumor angiogenesis, and metastasis. The Ras gene family comprised of 36 total genes in human. The N-Ras, K-Ras, and H-Ras are accounted for to assume noticeable function in human cancer. The mutation in K-Ras protein is most commonly found in tumors. K-Ras is the most crucial driver in lung and pancreatic cancers. Among the mutations of N-Ras, H-Ras, and K-Ras, the mutant K-Ras is the most prevalent target for the development of Lungs, colon, and pancreatic cancers. The study aimed to develop the peptide inhibitors of the K-Ras G12D. The crystal structure of the mutant K-Ras/R11.1.6 G12D complex was retrieved from the protein databank. The protein R11.1.6 directly blocks interaction with Raf and diminishes signaling through the Raf-MEK-ERK signaling pathway. Here, in this study, we designed novel peptides from the truncated reference peptide (R11.1.6) through residue scan methodology. The top ten designed peptides (based on binding free energies) were subjected to molecular dynamics simulations using AMBER to evaluate stability. Our results indicate that the top ten selected peptides have strong interactions with K-Ras than the reference peptide (R11.1.6) and have the potency to prevent the binding of Raf and K-Ras.Communicated by Ramaswamy H. Sarma.

G Protein-Coupled Estrogen Receptor Production Using an Escherichia coli Cell-Free Expression System.

Heterologous expression of the G protein-coupled estrogen receptor (GPER) comes with a suite of challenges intrinsic to membrane proteins. This receptor's low expression levels and tendency to form insoluble aggregates in Escherichia coli and yeast make it a difficult receptor-target to study. In this unit, we detail steps to produce monomeric GPER using a precipitation-based cell-free system. We provide information on the DNA construct for expression, the pipetting scheme for the reaction supplements to generate a master mix, and the cell-free reaction setup. In the last portion of this unit, we outline steps for solubilization and purification, and we provide a viable method for qualitatively observing functionality by liquid chromatography-mass spectrometry detection. © 2019 by John Wiley & Sons, Inc.

A Detailed Protocol for Large-scale Recombinant Expression and Validation of Human FGFR2 with Its Transmembrane and Extracellular Domains in Escherichia coli.

Receptor tyrosine kinases (RTKs) are an important class of transmembrane receptors that mediate some of the most crucial biochemical pathways essential to the growth, differentiation, and survival of a cell and thus, are highly involved in cancers. Due to the complexity of RTKs having biochemically different domains including a transmembrane domain, an intact crystal structure of any of these proteins remain elusive as it is difficult to produce milligram amounts of intact functional RTKs for crystallography studies. A heavily studied RTK is fibroblast growth factor receptor 2 (FGFR2), which plays a key role in fibroblast growth regulation, differentiation, and oncogenesis. Previous studies have focused on expressing FGFR2's extracellular, transmembrane, and intracellular domains individually. For this protocol, we have focused on the extracellular and transmembrane domains of the FGFR2 protein. The function of the expressed protein is validated by demonstrating its ability to bind heparin and fibroblast growth factor 1 (FGF1). The primary contribution of our protocol is expressing two RTK domains together, including the transmembrane domain, in milligram quantities. Being able to express RTKs to define its crystal structures would enable pharmacologists to design cancer drugs that selectively target active conformations.

Mutation of Phenylalanine-223 to Leucine Enhances Transformation of Benzo[a]pyrene by Ring-Hydroxylating Dioxygenase of Sphingobium sp. FB3 by increasing Accessibility of the Catalytic Site.

Burning of agricultural biomass generates polycyclic aromatic hydrocarbons (PAHs) including the carcinogen benzo[a]pyrene, of which the catabolism is primarily initiated by a ring-hydroxylating dioxygenase (RHD). This study explores catalytic site accessibility and its role in preferential catabolism of some PAHs over others. The genes flnA1f, flnA2f, flnA3, and flnA4, encoding the oxygenase α and ß subunits, ferredoxin, and ferredoxin reductase, respectively, of the RHD enzyme complex (FlnA) were cloned from Sphingobium sp. FB3 and coexpressed in E. coli BL21. The FlnA effectively transformed fluoranthene but not benzo[a]pyrene. Substitution of the bulky phenylalanine-223 by leucine reduces the steric constraint in the substrate entrance to make the catalytic site of FlnA more accessible to large substrates, as visualized by 3D modeling, and allows the FlnA mutant to efficiently transform benzo[a]pyrene. Accessibility of the catalytic site to PAHs is a mechanism of RHD substrate specificity. The results shed light on why some PAHs are more recalcitrant than others.

Protection against ß-amyloid neurotoxicity by a non-toxic endogenous N-terminal ß-amyloid fragment and its active hexapeptide core sequence.

High levels (µM) of beta amyloid (Aß) oligomers are known to trigger neurotoxic effects, leading to synaptic impairment, behavioral deficits, and apoptotic cell death. The hydrophobic C-terminal domain of Aß, together with sequences critical for oligomer formation, is essential for this neurotoxicity. However, Aß at low levels (pM-nM) has been shown to function as a positive neuromodulator and this activity resides in the hydrophilic N-terminal domain of Aß. An N-terminal Aß fragment (1-15/16), found in cerebrospinal fluid, was also shown to be a highly active neuromodulator and to reverse Aß-induced impairments of long-term potentiation. Here, we show the impact of this N-terminal Aß fragment and a shorter hexapeptide core sequence in the Aß fragment (Aßcore: 10-15) to protect or reverse Aß-induced neuronal toxicity, fear memory deficits and apoptotic death. The neuroprotective effects of the N-terminal Aß fragment and Aßcore on Aß-induced changes in mitochondrial function, oxidative stress, and apoptotic neuronal death were demonstrated via mitochondrial membrane potential, live reactive oxygen species, DNA fragmentation and cell survival assays using a model neuroblastoma cell line (differentiated NG108-15) and mouse hippocampal neuron cultures. The protective action of the N-terminal Aß fragment and Aßcore against spatial memory processing deficits in amyloid precursor protein/PSEN1 (5XFAD) mice was demonstrated in contextual fear conditioning. Stabilized derivatives of the N-terminal Aßcore were also shown to be fully protective against Aß-triggered oxidative stress. Together, these findings indicate an endogenous neuroprotective role for the N-terminal Aß fragment, while active stabilized N-terminal Aßcore derivatives offer the potential for therapeutic application.

Recombinant expression in E. coli of human FGFR2 with its transmembrane and extracellular domains.

Fibroblast growth factor receptors (FGFRs) are a family of receptor tyrosine kinases containing three domains: an extracellular receptor domain, a single transmembrane helix, and an intracellular tyrosine kinase domain. FGFRs are activated by fibroblast growth factors (FGFs) as part of complex signal transduction cascades regulating angiogenesis, skeletal formation, cell differentiation, proliferation, cell survival, and cancer. We have developed the first recombinant expression system in E. coli to produce a construct of human FGFR2 containing its transmembrane and extracellular receptor domains. We demonstrate that the expressed construct is functional in binding heparin and dimerizing. Size exclusion chromatography demonstrates that the purified FGFR2 does not form a complex with FGF1 or adopts an inactive dimer conformation. Progress towards the successful recombinant production of intact FGFRs will facilitate further biochemical experiments and structure determination that will provide insight into how extracellular FGF binding activates intracellular kinase activity.

Simulations reveal increased fluctuations in estrogen receptor-alpha conformation upon antagonist binding.

Molecular dynamics (MD) simulations have been used to model dynamic fluctuations in the structure of estrogen receptor-alpha (ER-α) upon binding to the natural agonist 17ß-estradiol (E2) and to the active metabolite of the breast cancer drug and antagonist, 4-hydroxytamoxifen (OHT). We present the most extensive MD simulations to date of ER-α, with over 1µs of combined simulations for the monomer and dimer forms. Simulations reveal that the antagonist-bound complex includes significant fluctuations while the agonist-bound complex is tightly restrained. OHT increases dynamic disorder in the loops located to either side of the tail H12 helix; H12 has been associated with the activation status of ER-α. We also report that fluctuations near H12 lead to greater conformational variation in the binding mode of the ethylamine tail of OHT. Both the agonist and antagonist conformations are stable throughout the 240ns simulations, supporting the hypothesis that there are no transitions between these two states or into intermediate states. The stable position of H12 in the OHT-bound conformation suggests that OHT stabilizes a well-defined antagonist conformational ensemble rather than merely blocking the agonist-driven activation of ER-α. Simultaneously, the increased dynamic properties of the OHT-bound complex is a potential source of binding entropy.

Structure of a TCR-Mimic Antibody with Target Predicts Pharmacogenetics.

Antibody therapies currently target only extracellular antigens. A strategy to recognize intracellular antigens is to target peptides presented by immune HLA receptors. ESK1 is a human, T-cell receptor (TCR)-mimic antibody that binds with subnanomolar affinity to the RMF peptide from the intracellular Wilms tumor oncoprotein WT1 in complex with HLA-A*02:01. ESK1 is therapeutically effective in mouse models of WT1(+) human cancers. TCR-based therapies have been presumed to be restricted to one HLA subtype. The mechanism for the specificity and high affinity of ESK1 is unknown. We show in a crystal structure that ESK1 Fab binds to RMF/HLA-A*02:01 in a mode different from that of TCRs. From the structure, we predict and then experimentally confirm high-affinity binding with multiple other HLA-A*02 subtypes, broadening the potential patient pool for ESK1 therapy. Using the crystal structure, we also predict potential off-target binding that we experimentally confirm. Our results demonstrate how protein structure information can contribute to personalized immunotherapy.

An unusual sugar conformation in the structure of an RNA/DNA decamer of the polypurine tract may affect recognition by RNase H.

Retroviral conversion of single-stranded RNA into double-stranded DNA requires priming for each strand. While host cellular t-RNA serves as primer for the first strand, the viral polypurine tract (PPT) is primer for the second. Therefore, polypurine tracts of retroviruses are essential for viral replication by reverse transcriptase (RT). These purine tracts are resistant to cleavage during first strand synthesis. In obtaining the primer for second strand synthesis, the RNase H function of RT must cleave the PPT exactly for in vivo transcription to proceed efficiently and proper integration to occur. At the RNase H active site the protein makes contacts primarily along the backbone, with hydrogen bonds to the sugar-phosphate oxygen atoms. A high-resolution structure (1.10A) of the first ten base-pairs of the RNA/DNA hybrid PPT, r-(c-a-a-a-g-a-a-a-a-g)/d-(C-T-T-T-T-C-T-T-T-G), contains the highly deformable r-(a-g-a) steps found in retroviral polypurine tracts. This r-(a-g-a) motif is utilized in the "unzipping" or unpairing of bases that occurs when RT binds a malleable PPT. Another unusual feature found in our high-resolution PPT structure is the sugar switch at RNA adenine 2. All the RNA sugars are the expected C3'-endo, except sugar 2, which is C2'-endo, characteristic of B-form sugars. This local A-to-B conversion adversely affects the pattern of hydrogen bonds from protein to sugar-phosphate backbone, disrupting the catalytic site. Disruption could cause the enzyme to pause at the 5'-end of the PPT, leaving it intact. Pyrimidine-purine (YR) steps are most deformable and the T-A step especially can undergo A-to-B transitions readily.