Covalent Modification of the JH2 Domain of Janus Kinase 2.

Probe molecules that covalently modify the JAK2 pseudokinase domain (JH2) are reported. Selective targeting of JH2 domains over the kinase (JH1) domains is a necessary feature for ligands intended to evaluate JH2 domains as therapeutic targets. The JH2 domains of three Janus kinases (JAK1, JAK2, and TYK2) possess a cysteine residue in the catalytic loop that does not occur in their JH1 domains. Starting from a non-selective kinase binding molecule, computer-aided design directed attachment of substituents terminating in acrylamide warheads to modify Cys675 of JAK2 JH2. Successful covalent attachment was demonstrated first through observation of enhanced binding with increasing incubation time in fluorescence polarization experiments. Covalent binding also increased selectivity to as much as ca. 30-fold for binding the JAK2 JH2 domain over the JH1 domain after a 20-h incubation. Covalency was confirmed through HPLC electrospray quadrupole time-of-flight HRMS experiments, which revealed the expected mass shifts.

Insights on JAK2 Modulation by Potent, Selective, and Cell-Permeable Pseudokinase-Domain Ligands.

JAK2 is a non-receptor tyrosine kinase that regulates hematopoiesis through the JAK-STAT pathway. The pseudokinase domain (JH2) is an important regulator of the activity of the kinase domain (JH1). V617F mutation in JH2 has been associated with the pathogenesis of various myeloproliferative neoplasms, but JAK2 JH2 has been poorly explored as a pharmacological target. In light of this, we aimed to develop JAK2 JH2 binders that could selectively target JH2 over JH1 and test their capacity to modulate JAK2 activity in cells. Toward this goal, we optimized a diaminotriazole lead compound into potent, selective, and cell-permeable JH2 binders leveraging computational design, synthesis, binding affinity measurements for the JH1, JH2 WT, and JH2 V617F domains, permeability measurements, crystallography, and cell assays. Optimized diaminotriazoles are capable of inhibiting STAT5 phosphorylation in both WT and V617F JAK2 in cells.

Conversion of a False Virtual Screen Hit into Selective JAK2 JH2 Domain Binders Using Convergent Design Strategies.

The Janus kinase 2 (JAK2) pseudokinase domain (JH2) is an ATP-binding domain that regulates the activity of the catalytic tyrosine kinase domain (JH1). Dysregulation of JAK2 JH1 signaling caused by the V617F mutation in JH2 is implicated in various myeloproliferative neoplasms. To explore if JAK2 activity can be modulated by a small molecule binding to the ATP site in JH2, we have developed several ligand series aimed at selectively targeting the JAK2 JH2 domain. We report here the evolution of a false virtual screen hit into a new JAK2 JH2 series. Optimization guided by computational modeling has yielded analogues with nanomolar affinity for the JAK2 JH2 domain and >100-fold selectivity for the JH2 domain over the JH1 domain. A crystal structure for one of the potent compounds bound to JAK2 JH2 clarifies the origins of the strong binding and selectivity. The compounds expand the platform for seeking molecules to regulate JAK2 signaling, including V617F JAK2 hyperactivation.

Structure-guided design of a perampanel-derived pharmacophore targeting the SARS-CoV-2 main protease.

There is a clinical need for direct-acting antivirals targeting SARS-CoV-2, the coronavirus responsible for the COVID-19 pandemic, to complement current therapeutic strategies. The main protease (Mpro) is an attractive target for antiviral therapy. However, the vast majority of protease inhibitors described thus far are peptidomimetic and bind to the active-site cysteine via a covalent adduct, which is generally pharmacokinetically unfavorable. We have reported the optimization of an existing FDA-approved chemical scaffold, perampanel, to bind to and inhibit Mpro noncovalently with IC50s in the low-nanomolar range and EC50s in the low-micromolar range. Here, we present nine crystal structures of Mpro bound to a series of perampanel analogs, providing detailed structural insights into their mechanism of action and structure-activity relationship. These insights further reveal strategies for pursuing rational inhibitor design efforts in the context of considerable active-site flexibility and potential resistance mechanisms.

Structural investigation of 2-naphthyl phenyl ether inhibitors bound to WT and Y181C reverse transcriptase highlights key features of the NNRTI binding site.

Human immunodeficiency virus (HIV)-1 remains as a global health issue that is primarily treated with highly active antiretroviral therapy, a combination of drugs that target the viral life cycle. One class of these drugs are non-nucleoside reverse transcriptase inhibitors (NNRTIs) that target the viral reverse transcriptase (RT). First generation NNRTIs were troubled with poor pharmacological properties and drug resistance, incentivizing the development of improved compounds. One class of developed compounds are the 2-naphthyl phenyl ethers, showing promising efficacy against the Y181C RT mutation. Further biochemical and structural work demonstrated differences in potency against the Y181C mutation and binding mode of the compounds. This work aims to understand the relationship between the binding mode and ability to overcome drug resistance using macromolecular x-ray crystallography. Comparison of 2-naphthyl phenyl ethers bound to Y181C RT reveal that compounds that interact with the invariant W229 are more capable of retaining efficacy against the resistance mutation. Additional modifications to these compounds at the 4-position, computationally designed to compensate for the Y181C mutation, do not demonstrate improved potency. Ultimately, we highlight important considerations for the development of future HIV-1 drugs that are able to combat drug resistance.

Structure-Guided Identification of DNMT3B Inhibitors.

Methyltransferase 3 beta (DNMT3B) inhibitors that interfere with cancer growth are emerging possibilities for treatment of melanoma. Herein we identify small molecule inhibitors of DNMT3B starting from a homology model based on a DNMT3A crystal structure. Virtual screening by docking led to purchase of 15 compounds, among which 5 were found to inhibit the activity of DNMT3B with IC50 values of 13-72 µM in a fluorogenic assay. Eight analogues of 7, 10, and 12 were purchased to provide 2 more active compounds. Compound 11 is particularly notable as it shows good selectivity with no inhibition of DNMT1 and 22 µM potency toward DNMT3B. Following additional de novo design, exploratory synthesis of 17 analogues of 11 delivered 5 additional inhibitors of DNMT3B with the most potent being 33h with an IC50 of 8.0 µM. This result was well confirmed in an ultrahigh-performance liquid chromatography (UHPLC)-based analytical assay, which yielded an IC50 of 4.8 µM. Structure-activity data are rationalized based on computed structures for DNMT3B complexes.

Selective Janus Kinase 2 (JAK2) Pseudokinase Ligands with a Diaminotriazole Core.

Janus kinases (JAKs) are non-receptor tyrosine kinases that are essential components of the JAK-STAT signaling pathway. Associated aberrant signaling is responsible for many forms of cancer and disorders of the immune system. The present focus is on the discovery of molecules that may regulate the activity of JAK2 by selective binding to the JAK2 pseudokinase domain, JH2. Specifically, the Val617Phe mutation in JH2 stimulates the activity of the adjacent kinase domain (JH1) resulting in myeloproliferative disorders. Starting from a non-selective screening hit, we have achieved the goal of discovering molecules that preferentially bind to the ATP binding site in JH2 instead of JH1. We report the design and synthesis of the compounds and binding results for the JH1, JH2, and JH2 V617F domains, as well as five crystal structures for JH2 complexes. Testing with a selective and non-selective JH2 binder on the autophosphorylation of wild-type and V617F JAK2 is also contrasted.

Structural and pharmacological evaluation of a novel non-nucleoside reverse transcriptase inhibitor as a promising long acting nanoformulation for treating HIV.

Combination antiretroviral therapy (cART) has been proven effective in inhibiting human immunodeficiency virus type 1 (HIV-1) infection and has significantly improved the health outcomes in acquired immune deficiency syndrome (AIDS) patients. The therapeutic benefits of cART have been challenged because of the toxicity and emergence of drug-resistant HIV-1 strains along with lifelong patient compliance resulting in non-adherence. These issues also hinder the clinical benefits of non-nucleoside reverse transcriptase inhibitors (NNRTIs), which are one of the vital components of cART for the treatment of HIV-1 infection. In this study, using a computational and structural based drug design approach, we have discovered an effective HIV -1 NNRTI, compound I (Cmpd I) that is very potent in biochemical assays and which targets key residues in the allosteric binding pocket of wild-type (WT)-RT as revealed by structural studies. Furthermore, Cmpd I exhibited very potent antiviral activity in HIV-1 infected T cells, lacked cytotoxicity (therapeutic index >100,000), and no significant off-target effects were noted in pharmacological assays. To address the issue of non-adherence, we developed a long-acting nanoformulation of Cmpd I (Cmpd I-NP) using poly (lactide-coglycolide) (PLGA) particles. The pharmacokinetic studies of free and nanoformulated Cmpd I were carried out in BALB/c mice. Intraperitoneal administration of Cmpd I and Cmpd I-NP in BALB/c mice revealed prolonged serum residence time of 48 h and 30 days, respectively. The observed serum concentrations of Cmpd I in both cases were sufficient to provide >97% inhibition in HIV-1 infected T-cells. The significant antiviral activity along with favorable pharmacological and pharmacokinetic profile of Cmpd I, provide compelling and critical support for its further development as an anti-HIV therapeutic agent.

Novel allosteric covalent inhibitors of bifunctional Cryptosporidium hominis TS-DHFR from parasitic protozoa identified by virtual screening.

Protozoans of the genus Cryptosporidium are the causative agent of the gastrointestinal disease, cryptosporidiosis, which can be fatal in immunocompromised individuals. Cryptosporidium hominis (C. hominis) bifunctional thymidylate synthase-dihydrofolate reductase (TS-DHFR) is an essential enzyme in the folate biosynthesis pathway and a molecular target for inhibitor design. Previous studies have demonstrated the importance of the ChTS-DHFR linker region "crossover helix" to the enzymatic activity and stability of the ChDHFR domain. We conducted a virtual screen of a novel non-active site pocket located at the interface of the ChDHFR domain and crossover helix. From this screen we have identified and characterized a noncompetitive inhibitor, compound 15, a substituted diphenyl thiourea. Through subsequent structure activity relationship studies, we have identified a time-dependent inhibitor lead, compound 15D17, a thiol-substituted 2-hydroxy-N-phenylbenzamide, which is selective for ChTS-DHFR, and whose effects appear to be mediated by covalent bond formation with a non-catalytic cysteine residue adjacent to the non-active site pocket.

Computationally-guided optimization of small-molecule inhibitors of the Aurora A kinase-TPX2 protein-protein interaction.

Free energy perturbation theory, in combination with enhanced sampling of protein-ligand binding modes, is evaluated in the context of fragment-based drug design, and used to design two new small-molecule inhibitors of the Aurora A kinase-TPX2 protein-protein interaction.

Improved Description of Sulfur Charge Anisotropy in OPLS Force Fields: Model Development and Parameterization.

The atomic point-charge model used in most molecular mechanics force fields does not represent well the electronic anisotropy that is featured in many directional noncovalent interactions. Sulfur participates in several types of such interactions with its lone pairs and σ-holes. The current study develops a new model, via the addition of off-atom charged sites, for a variety of sulfur compounds in the OPLS-AA and OPLS/CM5 force fields to address the lack of charge anisotropy. Parameter optimization is carried out to reproduce liquid-state properties, torsional and noncovalent energetics from reliable quantum mechanical calculations, and electrostatic potentials. Significant improvements are obtained for computed free energies of hydration, reducing the mean unsigned errors from ca. 1.4 to 0.4-0.7 kcal/mol. Enhanced accuracy in directionality and energetics is also obtained for molecular complexes with sulfur-containing hydrogen and halogen bonds. Moreover, the new model reproduces the unusual conformational preferences of sulfur-containing compounds with 1,4-intramolecular chalcogen bonds. Transferability of the new force field parameters to cysteine and methionine is verified via molecular dynamic simulations of blocked dipeptides. The study demonstrates the effectiveness of using off-atom charge sites to address electronic anisotropy, and provides a parametrization methodology that can be applied to other systems.

Identification and Characterization of JAK2 Pseudokinase Domain Small Molecule Binders.

Janus kinases (JAKs) regulate hematopoiesis via the cytokine-mediated JAK-STAT signaling pathway. JAKs contain tandem C-terminal pseudokinase (JH2) and tyrosine kinase (JH1) domains. The JAK2 pseudokinase domain adopts a protein kinase fold and, despite its pseudokinase designation, binds ATP with micromolar affinity. Recent evidence shows that displacing ATP from the JAK2 JH2 domain alters the hyperactivation state of the oncogenic JAK2 V617F protein while sparing the wild type JAK2 protein. In this study, small molecule binders of JAK2 JH2 were identified via an in vitro screen. Top hits were characterized using biophysical and structural approaches. Development of pseudokinase-selective compounds may offer novel pharmacological opportunities for treating cancers driven by JAK2 V617F and other oncogenic JAK mutants.

A Fluorescence Polarization Assay for Binding to Macrophage Migration Inhibitory Factor and Crystal Structures for Complexes of Two Potent Inhibitors.

Human macrophage migration inhibitory factor (MIF) is both a keto-enol tautomerase and a cytokine associated with numerous inflammatory diseases and cancer. Consistent with observed correlations between inhibition of the enzymatic and biological activities, discovery of MIF inhibitors has focused on monitoring the tautomerase activity using l-dopachrome methyl ester or 4-hydroxyphenyl pyruvic acid as substrates. The accuracy of these assays is compromised by several issues including substrate instability, spectral interference, and short linear periods for product formation. In this work, we report the syntheses of fluorescently labeled MIF inhibitors and their use in the first fluorescence polarization-based assay to measure the direct binding of inhibitors to the active site. The assay allows the accurate and efficient identification of competitive, noncompetitive, and covalent inhibitors of MIF in a manner that can be scaled for high-throughput screening. The results for 22 compounds show that the most potent MIF inhibitors bind with Kd values of ca. 50 nM; two are from our laboratory, and the other is a compound from the patent literature. X-ray crystal structures for two of the most potent compounds bound to MIF are also reported here. Striking combinations of protein-ligand hydrogen bonding, aryl-aryl, and cation-π interactions are responsible for the high affinities. A new chemical series was then designed using this knowledge to yield two more strong MIF inhibitors/binders.

OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins.

The parametrization and validation of the OPLS3 force field for small molecules and proteins are reported. Enhancements with respect to the previous version (OPLS2.1) include the addition of off-atom charge sites to represent halogen bonding and aryl nitrogen lone pairs as well as a complete refit of peptide dihedral parameters to better model the native structure of proteins. To adequately cover medicinal chemical space, OPLS3 employs over an order of magnitude more reference data and associated parameter types relative to other commonly used small molecule force fields (e.g., MMFF and OPLS_2005). As a consequence, OPLS3 achieves a high level of accuracy across performance benchmarks that assess small molecule conformational propensities and solvation. The newly fitted peptide dihedrals lead to significant improvements in the representation of secondary structure elements in simulated peptides and native structure stability over a number of proteins. Together, the improvements made to both the small molecule and protein force field lead to a high level of accuracy in predicting protein-ligand binding measured over a wide range of targets and ligands (less than 1 kcal/mol RMS error) representing a 30% improvement over earlier variants of the OPLS force field.

Improved Peptide and Protein Torsional Energetics with the OPLSAA Force Field.

The development and validation of new peptide dihedral parameters are reported for the OPLS-AA force field. High accuracy quantum chemical methods were used to scan φ, ψ, χ1, and χ2 potential energy surfaces for blocked dipeptides. New Fourier coefficients for the dihedral angle terms of the OPLS-AA force field were fit to these surfaces, utilizing a Boltzmann-weighted error function and systematically examining the effects of weighting temperature. To prevent overfitting to the available data, a minimal number of new residue-specific and peptide-specific torsion terms were developed. Extensive experimental solution-phase and quantum chemical gas-phase benchmarks were used to assess the quality of the new parameters, named OPLS-AA/M, demonstrating significant improvement over previous OPLS-AA force fields. A Boltzmann weighting temperature of 2000 K was determined to be optimal for fitting the new Fourier coefficients for dihedral angle parameters. Conclusions are drawn from the results for best practices for developing new torsion parameters for protein force fields.

Discovery and crystallography of bicyclic arylaminoazines as potent inhibitors of HIV-1 reverse transcriptase.

Non-nucleoside inhibitors of HIV-1 reverse transcriptase (HIV-RT) are reported that incorporate a 7-indolizinylamino or 2-naphthylamino substituent on a pyrimidine or 1,3,5-triazine core. The most potent compounds show below 10 nanomolar activity towards wild-type HIV-1 and variants bearing Tyr181Cys and Lys103Asn/Tyr181Cys resistance mutations. The compounds also feature good aqueous solubility. Crystal structures for two complexes enhance the analysis of the structure-activity data.

Design, synthesis, and protein crystallography of biaryltriazoles as potent tautomerase inhibitors of macrophage migration inhibitory factor.

Optimization is reported for biaryltriazoles as inhibitors of the tautomerase activity of human macrophage migration inhibitory factor (MIF), a proinflammatory cytokine associated with numerous inflammatory diseases and cancer. A combined approach was taken featuring organic synthesis, enzymatic assaying, crystallography, and modeling including free-energy perturbation (FEP) calculations. X-ray crystal structures for 3a and 3b bound to MIF are reported and provided a basis for the modeling efforts. The accommodation of the inhibitors in the binding site is striking with multiple hydrogen bonds and aryl-aryl interactions. Additional modeling encouraged pursuit of 5-phenoxyquinolinyl analogues, which led to the very potent compound 3s. Activity was further enhanced by addition of a fluorine atom adjacent to the phenolic hydroxyl group as in 3w, 3z, 3aa, and 3bb to strengthen a key hydrogen bond. It is also shown that physical properties of the compounds can be modulated by variation of solvent-exposed substituents. Several of the compounds are likely the most potent known MIF tautomerase inhibitors; the most active ones are more than 1000-fold more active than the well-studied (R)-ISO-1 and more than 200-fold more active than the chromen-4-one Orita-13.

Structural studies provide clues for analog design of specific inhibitors of Cryptosporidium hominis thymidylate synthase-dihydrofolate reductase.

Cryptosporidium is the causative agent of a gastrointestinal disease, cryptosporidiosis, which is often fatal in immunocompromised individuals and children. Thymidylate synthase (TS) and dihydrofolate reductase (DHFR) are essential enzymes in the folate biosynthesis pathway and are well established as drug targets in cancer, bacterial infections, and malaria. Cryptosporidium hominis has a bifunctional thymidylate synthase and dihydrofolate reductase enzyme, compared to separate enzymes in the host. We evaluated lead compound 1 from a novel series of antifolates, 2-amino-4-oxo-5-substituted pyrrolo[2,3-d]pyrimidines as an inhibitor of Cryptosporidium hominis thymidylate synthase with selectivity over the human enzyme. Complementing the enzyme inhibition compound 1 also has anti-cryptosporidial activity in cell culture. A crystal structure with compound 1 bound to the TS active site is discussed in terms of several van der Waals, hydrophobic and hydrogen bond interactions with the protein residues and the substrate analog 5-fluorodeoxyuridine monophosphate (TS), cofactor NADPH and inhibitor methotrexate (DHFR). Another crystal structure in complex with compound 1 bound in both the TS and DHFR active sites is also reported here. The crystal structures provide clues for analog design and for the design of ChTS-DHFR specific inhibitors.

Biochemical assays for the discovery of TDP1 inhibitors.

Drug screening against novel targets is warranted to generate biochemical probes and new therapeutic drug leads. TDP1 and TDP2 are two DNA repair enzymes that have yet to be successfully targeted. TDP1 repairs topoisomerase I-, alkylation-, and chain terminator-induced DNA damage, whereas TDP2 repairs topoisomerase II-induced DNA damage. Here, we report the quantitative high-throughput screening (qHTS) of the NIH Molecular Libraries Small Molecule Repository using recombinant human TDP1. We also developed a secondary screening method using a multiple loading gel-based assay where recombinant TDP1 is replaced by whole cell extract (WCE) from genetically engineered DT40 cells. While developing this assay, we determined the importance of buffer conditions for testing TDP1, and most notably the possible interference of phosphate-based buffers. The high specificity of endogenous TDP1 in WCE allowed the evaluation of a large number of hits with up to 600 samples analyzed per gel via multiple loadings. The increased stringency of the WCE assay eliminated a large fraction of the initial hits collected from the qHTS. Finally, inclusion of a TDP2 counter-screening assay allowed the identification of two novel series of selective TDP1 inhibitors.

Bifunctional inhibition of human immunodeficiency virus type 1 reverse transcriptase: mechanism and proof-of-concept as a novel therapeutic design strategy.

Human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) is a major target for currently approved anti-HIV drugs. These drugs are divided into two classes: nucleoside and non-nucleoside reverse transcriptase inhibitors (NRTIs and NNRTIs). This study illustrates the synthesis and biochemical evaluation of a novel bifunctional RT inhibitor utilizing d4T (NRTI) and a TMC-derivative (a diarylpyrimidine NNRTI) linked via a poly(ethylene glycol) (PEG) linker. HIV-1 RT successfully incorporates the triphosphate of d4T-4PEG-TMC bifunctional inhibitor in a base-specific manner. Moreover, this inhibitor demonstrates low nanomolar potency that has 4.3-fold and 4300-fold enhancement of polymerization inhibition in vitro relative to the parent TMC-derivative and d4T, respectively. This study serves as a proof-of-concept for the development and optimization of bifunctional RT inhibitors as potent inhibitors of HIV-1 viral replication.