Effect of Binding Pose and Modeled Structures on SVMGen and GlideScore Enrichment of Chemical Libraries

Virtual screening consists of docking libraries of small molecules to a target protein followed by rank-ordering of the resulting structures using scoring functions. The ability of scoring methods to distinguish between actives and inactives depends on several factors that include the accuracy of the binding pose during the docking step and the quality of the three-dimensional structure of the target. Here, we build on our previous work to introduce a new scoring approach (SVMGen) that uses machine learning trained with features from statistical pair potentials obtained from three-dimensional crystal structures. We use SVMGen and GlideScore to explore how enrichment or rank-ordering is affected by binding pose accuracy. To that end, we create a validation set that consists strictly of proteins whose crystal structure was solved in complex with their inhibitors. For the rank-ordering studies, we use crystal structures from PDBbind along with corresponding binding affinity data provided in the database. In addition to binding pose, we investigate the effect of using modeled structures for the target on the enrichment performance of SVMGen and GlideScore. To accomplish this, we generated homology models for protein kinases in DUD-E for which crystal structures are available to enable comparison of enrichment between modeled and crystal structure. We also generate homology models for kinases in SARfari for which there are many known small-molecule inhibitors but no known crystal structure. These models are used to assess the ability of SVMGen and GlideScore to distinguish between actives and decoys. We focus our work on protein kinases considering the wealth of structural and binding affinity data that exists for this family of proteins

Molecular Recognition in a Diverse Set of Protein-Ligand Interactions Studied with Molecular Dynamics Simulations and End-Point Free Energy Calculations

End-point free energy calculations using MM-GBSA and MM-PBSA provide a detailed understanding of molecular recognition in protein-ligand interactions. The binding free energy can be used to rank-order protein-ligand structures in virtual screening for compound or target identification. Here, we carry out free energy calculations for a diverse set of 11 proteins bound to 14 small molecules using extensive explicit-solvent MD simulations. The structure of these complexes was previously solved by crystallography and their binding studied with isothermal titration calorimetry (ITC) data enabling direct comparison to the MM-GBSA and MM-PBSA calculations. Four MM-GBSA and three MM-PBSA calculations reproduced the ITC free energy within 1 kcal•mol−1 highlighting the challenges in reproducing the absolute free energy from end-point free energy calculations. MM-GBSA exhibited better rank-ordering with a Spearman ρ of 0.68 compared to 0.40 for MM-PBSA with dielectric constant (ε = 1). An increase in ε resulted in significantly better rank-ordering for MM-PBSA (ρ = 0.91 for ε = 10). But larger ε significantly reduced the contributions of electrostatics, suggesting that the improvement is due to the non-polar and entropy components, rather than a better representation of the electrostatics. SVRKB scoring function applied to MD snapshots resulted in excellent rank-ordering (ρ = 0.81). Calculations of the configurational entropy using normal mode analysis led to free energies that correlated significantly better to the ITC free energy than the MD-based quasi-harmonic approach, but the computed entropies showed no correlation with the ITC entropy. When the adaptation energy is taken into consideration by running separate simulations for complex, apo and ligand (MM-PBSAADAPT), there is less agreement with the ITC data for the individual free energies, but remarkably good rank-ordering is observed (ρ = 0.89). Interestingly, filtering MD snapshots by pre-scoring protein-ligand complexes with a machine learning-based approach (SVMSP) resulted in a significant improvement in the MM-PBSA results (ε = 1) from ρ = 0.40 to ρ = 0.81. Finally, the non-polar components of MM-GBSA and MM-PBSA, but not the electrostatic components, showed strong correlation to the ITC free energy; the computed entropies did not correlate with the ITC entropy

A new class of orthosteric uPAR·uPA small-molecule antagonists are allosteric inhibitors of the uPAR·vitronectin interaction

The urokinase receptor (uPAR) is a GPI-anchored cell surface receptor that is at the center of an intricate network of protein-protein interactions. Its immediate binding partners are the serine proteinase urokinase (uPA), and vitronectin (VTN), a component of the extracellular matrix. uPA and VTN bind at distinct sites on uPAR to promote extracellular matrix degradation and integrin signaling, respectively. Here, we report the discovery of a new class of pyrrolone small-molecule inhibitors of the tight ∼1 nM uPAR·uPA protein-protein interaction. These compounds were designed to bind to the uPA pocket on uPAR. The highest affinity compound, namely 7, displaced a fluorescently labeled α-helical peptide (AE147-FAM) with an inhibition constant Ki of 0.7 μM and inhibited the tight uPAR·uPAATF interaction with an IC50 of 18 μM. Biophysical studies with surface plasmon resonance showed that VTN binding is highly dependent on uPA. This cooperative binding was confirmed as 7, which binds at the uPAR·uPA interface, also inhibited the distal VTN·uPAR interaction. In cell culture, 7 blocked the uPAR·uPA interaction in uPAR-expressing human embryonic kidney (HEK-293) cells and impaired cell adhesion to VTN, a process that is mediated by integrins. As a result, 7 inhibited integrin signaling in MDA-MB-231 cancer cells as evidenced by a decrease in focal adhesion kinase (FAK) phosphorylation and Rac1 GTPase activation. Consistent with these results, 7 blocked breast MDA-MB-231 cancer cell invasion with IC50 values similar to those observed in ELISA and surface plasmon resonance competition studies. Explicit-solvent molecular dynamics simulations show that the cooperativity between uPA and VTN is attributed to stabilization of uPAR motion by uPA. In addition, free energy calculations revealed that uPA stabilizes the VTNSMB·uPAR interaction through more favorable electrostatics and entropy. Disruption of the uPAR·VTNSMB interaction by 7 is consistent with the cooperative binding to uPAR by uPA and VTN. Interestingly, the VTNSMB·uPAR interaction was less favorable in the VTNSMB·uPAR·7 complex suggesting potential cooperativity between 7 and VTN. Compound 7 provides an excellent starting point for the development of more potent derivatives to explore uPAR biology

Enrichment of Chemical Libraries Docked to Protein Conformational Ensembles and Application to Aldehyde Dehydrogenase 2

Molecular recognition is a complex process that involves a large ensemble of structures of the receptor and ligand. Yet, most structure-based virtual screening is carried out on a single structure typically from X-ray crystallography. Explicit-solvent molecular dynamics (MD) simulations offer an opportunity to sample multiple conformational states of a protein. Here we evaluate our recently developed scoring method SVMSP in its ability to enrich chemical libraries docked to MD structures of seven proteins from the Directory of Useful Decoys (DUD). SVMSP is a target-specific rescoring method that combines machine learning with statistical potentials. We find that enrichment power as measured by the area under the ROC curve (ROC-AUC) is not affected by increasing the number of MD structures. Among individual MD snapshots, many exhibited enrichment that was significantly better than the crystal structure, but no correlation between enrichment and structural deviation from crystal structure was found. We followed an innovative approach by training SVMSP scoring models using MD structures (SVMSPMD). The resulting models were applied to two difficult cases (p38 and CDK2) for which enrichment was not better than random. We found remarkable increase in enrichment power, particularly for p38, where the ROC-AUC increased by 0.30 to 0.85. Finally, we explored approaches for a priori identification of MD snapshots with high enrichment power from an MD simulation in the absence of active compounds. We found that the use of randomly selected compounds docked to the target of interest using SVMSP led to notable enrichment for EGFR and Src MD snapshots. SVMSP rescoring of protein–compound MD structures was applied for the search of small-molecule inhibitors of the mitochondrial enzyme aldehyde dehydrogenase 2 (ALDH2). Rank-ordering of a commercial library of 50 000 compounds docked to MD structures of ALDH2 led to five small-molecule inhibitors. Four compounds had IC50s below 5 μM. These compounds serve as leads for the design and synthesis of more potent and selective ALDH2 inhibitors

Small-molecule CaVα1⋅CaVβ antagonist suppresses neuronal voltage-gated calcium-channel trafficking

Extracellular calcium flow through neuronal voltage-gated CaV2.2 calcium channels converts action potential-encoded information to the release of pronociceptive neurotransmitters in the dorsal horn of the spinal cord, culminating in excitation of the postsynaptic central nociceptive neurons. The CaV2.2 channel is composed of a pore-forming α1 subunit (CaVα1) that is engaged in protein-protein interactions with auxiliary α2/δ and β subunits. The high-affinity CaV2.2α1⋅CaVβ3 protein-protein interaction is essential for proper trafficking of CaV2.2 channels to the plasma membrane. Here, structure-based computational screening led to small molecules that disrupt the CaV2.2α1⋅CaVβ3 protein-protein interaction. The binding mode of these compounds reveals that three substituents closely mimic the side chains of hot-spot residues located on the α-helix of CaV2.2α1 Site-directed mutagenesis confirmed the critical nature of a salt-bridge interaction between the compounds and CaVβ3 Arg-307. In cells, compounds decreased trafficking of CaV2.2 channels to the plasma membrane and modulated the functions of the channel. In a rodent neuropathic pain model, the compounds suppressed pain responses. Small-molecule α-helical mimetics targeting ion channel protein-protein interactions may represent a strategy for developing nonopioid analgesia and for treatment of other neurological disorders associated with calcium-channel trafficking

Probing binding and cellular activity of pyrrolidinone and piperidinone small molecules targeting the urokinase receptor

The urokinase receptor (uPAR) is a cell-surface protein that is part of an intricate web of transient and tight protein interactions that promote cancer cell invasion and metastasis. Here, we evaluate the binding and biological activity of a new class of pyrrolidinone and piperidinone compounds, along with derivatives of previously-identified pyrazole and propylamine compounds. Competition assays revealed that the compounds displace a fluorescently labeled peptide (AE147-FAM) with inhibition constant (Ki ) values ranging from 6 to 63 μM. Structure-based computational pharmacophore analysis followed by extensive explicit-solvent molecular dynamics (MD) simulations and free energy calculations suggested the pyrazole-based and piperidinone-based compounds adopt different binding modes, despite their similar two-dimensional structures. In cells, pyrazole-based compounds showed significant inhibition of breast adenocarcinoma (MDA-MB-231) and pancreatic ductal adenocarcinoma (PDAC) cell proliferation, but piperidinone-containing compounds exhibited no cytotoxicity even at concentrations of 100 μM. One pyrazole-based compound impaired MDA-MB-231 invasion, adhesion, and migration in a concentration-dependent manner, while the piperidinone inhibited only invasion. The pyrazole derivative inhibited matrix metalloprotease-9 (gelatinase) activity in a concentration-dependent manner, while the piperidinone showed no effect suggesting different mechanisms for inhibition of cell invasion. Signaling studies further highlighted these differences, showing that pyrazole compounds completely inhibited ERK phosphorylation and impaired HIF1α and NF-κB signaling, while pyrrolidinones and piperidinones had no effect. Annexin V staining suggested that the effect of the pyrazole-based compound on proliferation was due to cell killing through an apoptotic mechanism. The compounds identified represent valuable leads in the design of further derivatives with higher affinities and potential probes to unravel the protein-protein interactions of uPAR

Small molecules inhibit STAT3 activation, autophagy, and cancer cell anchorage-independent growth

Triple-negative breast cancers (TNBCs) lack the signature targets of other breast tumors, such as HER2, estrogen receptor, and progesterone receptor. These aggressive basal-like tumors are driven by a complex array of signaling pathways that are activated by multiple driver mutations. Here we report the discovery of 6 (KIN-281), a small molecule that inhibits multiple kinases including maternal leucine zipper kinase (MELK) and the non-receptor tyrosine kinase bone marrow X-linked (BMX) with single-digit micromolar IC50s. Several derivatives of 6 were synthesized to gain insight into the binding mode of the compound to the ATP binding pocket. Compound 6 was tested for its effect on anchorage-dependent and independent growth of MDA-MB-231 and MDA-MB-468 breast cancer cells. The effect of 6 on BMX prompted us to evaluate its effect on STAT3 phosphorylation and DNA binding. The compound’s inhibition of cell growth led to measurements of survivin, Bcl-XL, p21WAF1/CIP1, and cyclin A2 levels. Finally, LC3B-II levels were quantified following treatment of cells with 6 to determine whether the compound affected autophagy, a process that is known to be activated by STAT3. Compound 6 provides a starting point for the development of small molecules with polypharmacology that can suppress TNBC growth and metastasis

The Importance of a Critical Protonation State and the Fate of the Catalytic Steps in Class A β-Lactamases and Penicillin-binding Proteins

b-Lactamases and penicillin-binding proteins are bacterial enzymes involved in antibiotic resistance to b-lactam antibiotics and biosynthetic assembly of cell wall, respectively. Members of these large families of enzymes all experience acylation by their respective substrates at an active-site serine as the first step in their catalytic activities. A Ser-X-X-Lys sequence motif is seen in all these proteins and crystal structures demonstrate that the side chain functions of the serine and lysine are in contact with one another. Three independent methods were used in this report to address the question of the protonation state of this important lysine (Lys73) in the TEM-1 b-lactamase from Escherichia coli. These techniques included perturbation of the pKa of Lys73 by the study of the g-thialysine-73 variant and the attendant kinetic analyses, investigation of the protonation state by titration of specifically labeled proteins by nuclear magnetic resonance and by computational treatment using the thermodynamic integration method. All three methods indicated that the pKa of Lys73 of this enzyme is attenuated to 8.0-8.5. It is argued herein that the unique ground-state ion pair of Glu166 and Lys73 of class A b-lactamases has actually raised the pKa of the active site lysine to 8.0-8.5 from that of the parental penicillin-binding protein. Whereas we cannot definitively rule out that Glu166 activates the active site water, which in turn promotes Ser70 for the acylation event, such as proposed earlier, we would like to propose as a plausible alternative for the acylation step the possibility that the ion pair would reconfigure to the protonated Glu166 and unprotonated Lys73. As such, unprotonated Lys73 could promote serine for acylation, a process that should be shared among all active-site-serine b-lactamases and penicillin-binding proteins

Small-Molecule Binding Sites to Explore New Targets in the Cancer Proteome

The Cancer Genome Atlas (TCGA) offers an unprecedented opportunity to identify small-molecule binding sites on proteins with overexpressed mRNA levels that correlate with poor survival. Here, we analyze RNA-seq and clinical data for 10 tumor types to identify genes that are both overexpressed and correlate with patient survival. Protein products of these genes were scanned for binding sites that possess shape and physicochemical properties that can accommodate small-molecule probes or therapeutic agents (druggable). These binding sites were classified as enzyme active sites (ENZ), protein-protein interaction sites (PPI), or other sites whose function is unknown (OTH). Interestingly, the overwhelming majority of binding sites were classified as OTH. We find that ENZ, PPI, and OTH binding sites often occurred on the same structure suggesting that many of these OTH cavities can be used for allosteric modulation of enzyme activity or protein-protein interactions with small molecules. We discovered several ENZ (PYCR1, QPRT, and HSPA6) and PPI (CASC5, ZBTB32, and CSAD) binding sites on proteins that have been seldom explored in cancer. We also found proteins that have been extensively studied in cancer that have not been previously explored with small molecules that harbor ENZ (PKMYT1, STEAP3, and NNMT) and PPI (HNF4A, MEF2B, and CBX2) binding sites. All binding sites were classified by the signaling pathways to which the protein that harbors them belongs using KEGG. In addition, binding sites were mapped onto structural protein-protein interaction networks to identify promising sites for drug discovery. Finally, we identify pockets that harbor missense mutations previously identified from analysis of the TCGA data. The occurrence of mutations in these binding sites provides new opportunities to develop small-molecule probes to explore their function in cancer.

Discovery and characterization of small molecules that target the Ral GTPase

The Ras-like GTPases RalA and B are important drivers of tumor growth and metastasis. Chemicals that block Ral function would be valuable as research tools and for cancer therapeutics. Here, we used protein structure analysis and virtual screening to identify drug-like molecules that bind a site on the GDP-form of Ral. Compounds RBC6, RBC8 and RBC10 inhibited Ral binding to its effector RalBP1, Ral-mediated cell spreading in murine fibroblasts and anchorage-independent growth of human cancer cell lines. Binding of RBC8 derivative BQU57 to RalB was confirmed by isothermal titration calorimetry, surface plasma resonance and 15N-HSQC NMR. RBC8 and BQU57 show selectivity for Ral relative to Ras or Rho and inhibit xenograft tumor growth similar to depletion of Ral by siRNA. Our results show the utility of structure-based discovery for development of therapeutics for Ral-dependent cancers