Bo Wang

In Silico Docking and Electrophysiological Characterization of Lacosamide Binding Sites on Collapsin Response Mediator Protein-2 Identifies a Pocket Important in Modulating Sodium Channel Slow Inactivation

The anti-epileptic drug (R)-lacosamide ((2R)-2-(acetylamino)-N-benzyl-3-methoxypropanamide (LCM)) modulates voltage-gated sodium channels (VGSCs) by preferentially interacting with slow inactivated sodium channels, but the observation that LCM binds to collapsin response mediator protein 2 (CRMP-2) suggests additional mechanisms of action for LCM. We postulated that CRMP-2 levels affects the actions of LCM on VGSCs. CRMP-2 labeling by LCM analogs was competitively displaced by excess LCM in rat brain lysates. Manipulation of CRMP-2 levels in the neuronal model system CAD cells affected slow inactivation of VGSCs without any effects on other voltage-dependent properties. In silico docking was performed to identify putative binding sites in CRMP-2 that may modulate the effects of LCM on VGSCs. These studies identified five cavities in CRMP-2 that can accommodate LCM. CRMP-2 alanine mutants of key residues within these cavities were functionally similar to wild-type CRMP-2 as assessed by similar levels of enhancement in dendritic complexity of cortical neurons. Next, we examined the effects of expression of wild-type and mutant CRMP-2 constructs on voltage-sensitive properties of VGSCs in CAD cells: 1) steady-state voltage-dependent activation and fast-inactivation properties were not affected by LCM, 2) CRMP-2 single alanine mutants reduced the LCM-mediated effects on the ability of endogenous Na+ channels to transition to a slow inactivated state, and 3) a quintuplicate CRMP-2 alanine mutant further decreased this slow inactivated fraction. Collectively, these results identify key CRMP-2 residues that can coordinate LCM binding thus making it more effective on its primary clinical target.

Support vector regression scoring of receptor-ligand complexes for rank-ordering and virtual screening of chemical libraries

The community structure-activity resource (CSAR) data sets are used to develop and test a support vector machine-based scoring function in regression mode (SVR). Two scoring functions (SVR-KB and SVR-EP) are derived with the objective of reproducing the trend of the experimental binding affinities provided within the two CSAR data sets. The features used to train SVR-KB are knowledge-based pairwise potentials, while SVR-EP is based on physicochemical properties. SVR-KB and SVR-EP were compared to seven other widely used scoring functions, including Glide, X-score, GoldScore, ChemScore, Vina, Dock, and PMF. Results showed that SVR-KB trained with features obtained from three-dimensional complexes of the PDBbind data set outperformed all other scoring functions, including best performing X-score, by nearly 0.1 using three correlation coefficients, namely Pearson, Spearman, and Kendall. It was interesting that higher performance in rank ordering did not translate into greater enrichment in virtual screening assessed using the 40 targets of the Directory of Useful Decoys (DUD). To remedy this situation, a variant of SVR-KB (SVR-KBD) was developed by following a target-specific tailoring strategy that we had previously employed to derive SVM-SP. SVR-KBD showed a much higher enrichment, outperforming all other scoring functions tested, and was comparable in performance to our previously derived scoring function SVM-SP.

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 nonpolar and entropy components, rather than a better representation of the electrostatics. The 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 prescoring 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 nonpolar 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.

Exploring a structural protein-drug interactome for new therapeutics in lung cancer.

The pharmacology of drugs is often defined by more than one protein target. This property can be exploited to use approved drugs to uncover new targets and signaling pathways in cancer. Towards enabling a rational approach to uncover new targets, we expand a structural protein-ligand interactome () by scoring the interaction among 1000 FDA-approved drugs docked to 2500 pockets on protein structures of the human genome. This afforded a drug-target network whose properties compared favorably with previous networks constructed using experimental data. Among drugs with the highest degree and betweenness two are cancer drugs and one is currently used for treatment of lung cancer. Comparison of predicted cancer and non-cancer targets reveals that the most cancer-specific compounds were also the most selective compounds. Analysis of compound flexibility, hydrophobicity, and size showed that the most selective compounds were low molecular weight fragment-like heterocycles. We use a previously-developed screening approach using the cancer drug erlotinib as a template to screen other approved drugs that mimic its properties. Among the top 12 ranking candidates, four are cancer drugs, two of them kinase inhibitors (like erlotinib). Cellular studies using non-small cell lung cancer (NSCLC) cells revealed that several drugs inhibited lung cancer cell proliferation. We mined patient records at the Regenstrief Medical Record System to explore the possible association of exposure to three of these drugs with occurrence of lung cancer. Preliminary in vivo studies using the non-small cell lung cancer (NCLSC) xenograft model showed that losartan- and astemizole-treated mice had tumors that weighed 50 (p < 0.01) and 15 (p < 0.01) percent less than the treated controls. These results set the stage for further exploration of these drugs and to uncover new drugs for lung cancer therapy.

Structure-based computational approaches for small-molecule modulation of protein-protein interactions.

Three-dimensional structures of proteins offer an opportunity for the rational design of small molecules to modulate protein-protein interactions. The presence of a well-defined binding pocket on the surface of protein complexes, particularly at their interface, can be used for docking-based virtual screening of chemical libraries. Several approaches have been developed to identify binding pockets that are implemented in programs such as SiteMap, fpocket, and FTSite. These programs enable the scoring of these pockets to determine whether they are suitable to accommodate high-affinity small molecules. Virtual screening of commercial or combinatorial libraries can be carried out to enrich these libraries and select compounds for further experimental validation. In virtual screening, a compound library is docked to the target protein. The resulting structures are scored and ranked for the selection and experimental validation of top candidates. Molecular docking has been implemented in a number of computer programs such as AutoDock Vina. We select a set of protein-protein interactions that have been successfully inhibited with small molecules in the past. Several computer programs are applied to identify pockets on the surface, and molecular docking is conducted in an attempt to reproduce the binding pose of the inhibitors. The results highlight the strengths and limitations of computational methods for the design of PPI inhibitors.

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.

Building a Chemical-Protein Interactome on the Open Science Grid

The Structural Protein-Ligand Interactome (SPLINTER) project predicts the interaction of thousands of small molecules with thousands of proteins. These interactions are predicted using the three-dimensional structure of the bound complex between each pair of protein and compound that is predicted by molecular docking. These docking runs consist of millions of individual short jobs each lasting only minutes. However, computing resources to execute these jobs (which cumulatively take tens of millions of CPU hours) are not readily or easily available in a cost effective manner. By looking to National Cyberinfrastructure resources, and specifically the Open Science Grid (OSG), we have been able to harness CPU power for researchers at the Indiana University School of Medicine to provide a quick and efficient solution to their unmet computing needs. Using the job submission infrastructure provided by the OSG, the docking data and simulation executable was sent to more than 100 universities and research centers worldwide. These opportunistic resources provided millions of CPU hours in a matter of days, greatly reducing time docking simulation time for the research group. The overall impact of this approach allows researchers to identify small molecule candidates for individual proteins, or new protein targets for existing FDA-approved drugs and biologically active compounds.