Jianhua Shen

TarFisDock: a web server for identifying drug targets with docking approach

TarFisDock is a web-based tool for automating the procedure of searching for small molecule–protein interactions over a large repertoire of protein structures. It offers PDTD (potential drug target database), a target database containing 698 protein structures covering 15 therapeutic areas and a reverse ligand–protein docking program. In contrast to conventional ligand–protein docking, reverse ligand–protein docking aims to seek potential protein targets by screening an appropriate protein database. The input file of this web server is the small molecule to be tested, in standard mol2 format; TarFisDock then searches for possible binding proteins for the given small molecule by use of a docking approach. The ligand–protein interaction energy terms of the program DOCK are adopted for ranking the proteins. To test the reliability of the TarFisDock server, we searched the PDTD for putative binding proteins for vitamin E and 4H-tamoxifen. The top 2 and 10% candidates of vitamin E binding proteins identified by TarFisDock respectively cover 30 and 50% of reported targets verified or implicated by experiments; and 30 and 50% of experimentally confirmed targets for 4H-tamoxifen appear amongst the top 2 and 5% of the TarFisDock predicted candidates, respectively. Therefore, TarFisDock may be a useful tool for target identification, mechanism study of old drugs and probes discovered from natural products. TarFisDock and PDTD are available at .

Virtual Screening on Natural Products for Discovering Active Compounds and Target Information

Natural products, containing inherently large-scale structural diversity than synthetic compounds, have been the major resources of bioactive agents and will continually play as protagonists for discovering new drugs. However, how to access this diverse chemical space efficiently and effectively is an exciting challenge for medicinal chemists and pharmacologists. While virtual screening, which has shown a great promise in drug discovery, will play an important role in digging out lead (active) compounds from natural products. This review focuses on the strategy of virtual screening based on molecular docking and, with successful examples from our laboratory, illustrates the efficiency of virtual screening in discovering active compounds from natural products. On the other hand, the sequencing of the human genome and numerous pathogen genomes has resulted in an unprecedented opportunity for discovering potential new drug targets. Chemogenomics has appeared as a new technology to initiate target discovery by using active compounds as probes to characterize proteome functions. Natural products are the ideal probes for such research. Binding affinity fingerprint is a powerful chemogenomic descriptor to characterize both small molecules and pharmacologically relevant proteins. Therefore, this review also discusses binding affinity fingerprint strategy for identifying target information from the genomic data by using natural products as the probes.

pH-dependent Conformational Flexibility of the SARS-CoV Main Proteinase (M pro ) Dimer: Molecular Dynamics Simulations and Multiple X-ray Structure Analyses

n n The SARS coronavirus main proteinase (Mpro) is a key enzyme in the processing of the viral polyproteins and thus an attractive target for the discovery of drugs directed against SARS. The enzyme has been shown by X-ray crystallography to undergo significant pH-dependent conformational changes. Here, we assess the conformational flexibility of the Mpro by analysis of multiple crystal structures (including two new crystal forms) and by molecular dynamics (MD) calculations. The MD simulations take into account the different protonation states of two histidine residues in the substrate-binding site and explain the pH-activity profile of the enzyme. The low enzymatic activity of the Mpro monomer and the need for dimerization are also discussed.n n

Molecular Dynamics Simulations on SDF-1α: Binding with CXCR4 Receptor

Insights into the interacting mode of CXCR4 with SDF-1α are crucial in understanding the structural and functional characteristics of CXCR4 receptor. In this paper a computational pipeline, integrating protein structure prediction, molecular dynamics simulations, automated molecular docking, and Brownian dynamics simulations were employed to investigate the dynamic and energetic aspects of CXCR4 associating with SDF-1α. The entire simulation revealed the surface distribution feature of electrostatic potentials and conformational “open-close” process of the receptor. The possible binding conformation of CXCR4 was identified, and the CXCR4—SDF-1α binding complex was generated. Arg188-Glu277 salt bridge plays an important role for both the extracellular domain conformational change and SDF-1α binding. Two binding sites were mapped at the extracellular domain (Site 1) and inside the transmembrane domain (Site 2), which are composed of conserved residues. Sites 1 and 2 contribute ∼60% and 40% to the binding affinity with SDF-1α, respectively. The binding model is in agreement with most of the experimental data. Transmembrane VI has more significant motion in the harmonious conformational transition of CXCR4 during SDF-1α binding, which may be possibly associated with signal transduction. Based on the modeling and simulation, a binding mechanism hypothesis between CXCR4 and SDF-1α and its relationship to the signal transduction has been proposed.

Nucleocapsid protein of SARS coronavirus tightly binds to human cyclophilin A

n Abstractn n Severe acute respiratory syndrome coronavirus (SARS-CoV) is responsible for SARS infection. Nucleocapsid protein (NP) of SARS-CoV (SARS_NP) functions in enveloping the entire genomic RNA and interacts with viron structural proteins, thus playing important roles in the process of virus particle assembly and release. Protein–protein interaction analysis using bioinformatics tools indicated that SARS_NP may bind to human cyclophilin A (hCypA), and surface plasmon resonance (SPR) technology revealed this binding with the equilibrium dissociation constant ranging from 6 to 160nM. The probable binding sites of these two proteins were detected by modeling the three-dimensional structure of the SARS_NP–hCypA complex, from which the important interaction residue pairs between the proteins were deduced. Mutagenesis experiments were carried out for validating the binding model, whose correctness was assessed by the observed effects on the binding affinities between the proteins. The reliability of the binding sites derived by the molecular modeling was confirmed by the fact that the computationally predicted values of the relative free energies of the binding for SARS_NP (or hCypA) mutants to the wild-type hCypA (or SARS_NP) are in good agreement with the data determined by SPR. Such presently observed SARS_NP–hCypA interaction model might provide a new hint for facilitating the understanding of another possible SARS-CoV infection pathway against human cell.n n

Molecular docking and 3-D-QSAR studies on the possible antimalarial mechanism of artemisinin analogues.

Artemisinin (Qinghaosu) is a natural constituent found in Artemisia annua L, which is an effective drug against chloroquine-resistant Plasmodium falciparum strains and cerebral malaria. The antimalarial activities of artemisinin and its analogues appear to be mediated by the interactions of the drugs with hemin. In order to understand the antimalarial mechanism and the relationship between the physicochemical properties and the antimalarial activities of artemisinin analogues, we performed molecular docking simulations to probe the interactions of these analogues with hemin, and then performed three-dimensional quantitative structure-activity relationship (3-D-QSAR) studies on the basis of the docking models employing comparative molecular force fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Molecular docking simulations generated probable bioactive conformations of artemisinin analogues and provided a new insight into the antimalarial mechanism. The subsequent partial least squares (PLS) analysis indicates that the calculate binding energies correlate well with the experimental activity values. The CoMFA and CoMSIA models based on the bioactive conformations proved to have good predictive ability and in turn match well with the docking result, which further testified the reliability of the docking model. Combining these results, that is molecular docking and 3-D-QSAR, together, the binding model and activity of new synthesized artemisinin derivatives were well explained.

Steered Molecular Dynamics Simulation on the Binding of NNRTI to HIV-1 RT

HIV-1 reverse transcriptase (RT) is the primary target for anti-AIDS chemotherapy. Nonnucleoside RT inhibitors (NNRTIs) are very potent and most promising anti-AIDS drugs that specifically inhibit HIV-1 RT. The binding and unbinding processes of alpha-APA, an NNRTI, have been studied using nanosecond conventional molecular dynamics and steered molecular dynamics simulations. The simulation results show that the unbinding process of alpha-APA consists of three phases based on the position of alpha-APA in relation to the entrance of the binding pocket. When alpha-APA is bound in the binding pocket, the hydrophobic interactions between HIV-1 RT and alpha-APA dominate the binding; however, the hydrophilic interactions (both direct and water-bridged hydrogen bonds) also contribute to the stabilizing forces. Whereas Tyr-181 makes significant hydrophobic interactions with alpha-APA, Tyr-188 forms a strong hydrogen bond with the acylamino group (N14) of alpha-APA. These two residues have very flexible side chains and appear to act as two flexible clamps discouraging alpha-APA to dissociate from the binding pocket. At the pocket entrance, two relatively inflexible residues, Val-179 and Leu-100, gauge the openness of the entrance and form the bottleneck of the inhibitor-unbinding pathway. Two special water molecules at the pocket entrance appear to play important roles in inhibitor recognition of binding and unbinding. These water molecules form water bridges between the polar groups of the inhibitor and the residues around the entrance, and between the polar groups of the inhibitor themselves. The water-bridged interactions not only induce the inhibitor to adopt an energetically favorable conformation so the inhibitor can pass through the pocket entrance, but also stabilize the binding of the inhibitor in the pocket to prevent the inhibitors dissociation. The complementary steered molecular dynamics and conventional molecular dynamics simulation results strongly support the hypothesis that NNRTIs inhibit HIV-1 RT polymerization activity by enlarging the DNA-binding cleft and restricting the flexibility and mobility of the p66 thumb subdomain that are believed to be essential during DNA translocation and polymerization.

Ligand-binding regulation of LXR/RXR and LXR/PPAR heterodimerizations: SPR technology-based kinetic analysis correlated with molecular dynamics simulation.

Liver X receptor (LXR) and peroxisome proliferator‐activated receptor (PPAR) are two members of nuclear receptors involved in the nutrient metabolisms of dietary fatty acid and cholesterol. They are found to be of cross‐talk function in that LXR regulates fatty acid synthesis and PPAR controls fatty acid degradation. LXRs (LXRα and LXRβ) function by forming obligate heterodimers with the retinoid X receptor (RXR), and subsequently binding to specific DNA response elements within the regulatory regions of their target genes. In this work, the kinetic features concerning LXR/RXR and LXR/PPAR interactions have been fully investigated using surface plasmon resonance (SPR) technology. It is found that LXRs could bind to all the three PPAR subtypes, PPARα, PPARγ and PPARδ with different binding affinities, and such receptor/receptor interactions could be regulated by ligand binding. Moreover, molecular dynamics (MD) simulations were performed on six typical complex models. The results revealed that ligands may increase the interaction energies between the receptor interfaces of the simulated receptor/receptor complexes. The MD results are in agreement with the SPR data. Further analyses on the MD results indicated that the ligand binding might increase the hydrogen bonds between the interfaces of the receptor/receptor complex.

Brownian Dynamics Simulations of the Recognition of the Scorpion Toxin Maurotoxin with the Voltage-Gated Potassium Ion Channels

The recognition of the scorpion toxin maurotoxin (MTX) by the voltage-gated potassium (Kv1) channels, Kv1.1, Kv1.2, and Kv1.3, has been studied by means of Brownian dynamics (BD) simulations. All of the 35 available structures of MTX in the Protein Data Bank (http://www.rcsb.org/pdb) determined by nuclear magnetic resonance were considered during the simulations, which indicated that the conformation of MTX significantly affected both the recognition and the binding between MTX and the Kv1 channels. Comparing the top five highest-frequency structures of MTX binding to the Kv1 channels, we found that the Kv1.2 channel, with the highest docking frequencies and the lowest electrostatic interaction energies, was the most favorable for MTX binding, whereas Kv1.1 was intermediate, and Kv1.3 was the least favorable one. Among the 35 structures of MTX, the 10th structure docked into the binding site of the Kv1.2 channel with the highest probability and the most favorable electrostatic interactions. From the MTX-Kv1.2 binding model, we identified the critical residues for the recognition of these two proteins through triplet contact analyses. MTX locates around the extracellular mouth of the Kv1 channels, making contacts with its beta-sheets. Lys23, a conserved amino acid in the scorpion toxins, protrudes into the pore of the Kv1.2 channel and forms two hydrogen bonds with the conserved residues Gly401(D) and Tyr400(C) and one hydrophobic contact with Gly401(C) of the Kv1.2 channel. The critical triplet contacts for recognition between MTX and the Kv1.2 channel are Lys23(MTX)-Asp402(C)(Kv1), Lys27(MTX)-Asp378(D)(Kv1), and Lys30(MTX)-Asp402(A)(Kv1). In addition, six hydrogen-bonding interactions are formed between residues Lys23, Lys27, Lys30, and Tyr32 of MTX and residues Gly401, Tyr400, Asp402, Asp378, and Thr406 of Kv1.2. Many of them are formed by side chains of residues of MTX and backbone atoms of the Kv1.2 channel. Five hydrophobic contacts exist between residues Pro20, Lys23, Lys30 and Tyr32 of MTX and residues Asp402, Val404, Gly401, and Arg377 of the Kv1.2 channel. The simulation results are in agreement with the previous molecular biology experiments and explain the binding phenomena between MTX and Kv1 channels at the molecular level. The consistency between the results of the BD simulations and the experimental data indicated that our three-dimensional model of the MTX-Kv1.2 channel complex is reasonable and can be used in additional biological studies, such as rational design of novel therapeutic agents blocking the voltage-gated channels and in mutagenesis studies in both the toxins and the Kv1 channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp378 of the Kv1.2 channel is critical for its recognition and binding functionality toward MTX. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating this might be a new clue for additional functional study of Kv1 channels.

Brownian Dynamics Simulations of the Recognition of the Scorpion Toxin P05 with the Small-conductance Calcium-activated Potassium Channels

The recognition of the scorpion toxin P05 and the small-conductance, calcium-activated potassium (SK) channels, rsk1, rsk2, and rsk3, has been studied by means of the Brownian dynamics (BD) method. All of the 25 available structures of P05 in the RCSB Protein Data Bank determined by NMR were considered during the simulation, which indicated that the conformation of P05 affects both the recognition and binding between the two proteins significantly. Comparing the top four high-frequency structures of P05 binding to the SK channels, we found that the rsk2 channel, with high frequencies and lowest electrostatic interaction energies (E (int)(ES)), is the most favorable for P05 binding, while rsk3 is intermediate, and rsk1 is the least favorable. Among the 25 structures of P05, the 13th structure docks into the binding site of the rsk2 channel with the highest probability and most favorable electrostatic interactions. From the P05-rsk2 channel binding model, we identified the residues critical for the recognition of these two proteins through triplet contact analyses. P05 locates around the extracellular mouth of the SK channels and contacts the SK channels using its alpha-helix rather than beta-sheets. The critical triplet contacts for recognition between P05 and the rsk2 channel are Arg6 (P05)-Asp364 (SK), Arg7 (P05)-Asn368 (SK), and Arg13 (P05)-Asp341 (SK). The structure of the P05-rsk2 complex with the most favorable electrostatic interaction energy was further refined by molecular mechanics, showing that six hydrogen bonding interactions exist between P05 and the rsk2 channel: one hydrogen bond is formed between Arg6 (P05) and Asp364(D) (rsk2); Arg7 (P05) forms three hydrogen bonds with Asp341(B) (rsk2)) and Asp364(C) (rsk2); two hydrogen bonds are formed by Arg13 (P05) with Asp341(A) (rsk2) and Asp364(B) (rsk2). The simulation results are in good agreement with the previous molecular biological experiments and can explain the binding phenomena between P05 and SK channels at the level of molecular structure. The consistency between the results of the BD simulations and the experimental data indicated that our 3D model of the P05-rsk2 channel complex is reasonable and can be employed in further biological studies, such as rational design of the novel therapeutic agents blocking the small-conductance, calcium-activated and apamin-sensitive potassium channels, and for mutagenesis studies in both toxins and SK channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp364 of the rsk2 channel is critical for its recognition and binding functionality towards P05. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating that this might be a new clue for further functional study of SK channels.