Biomedical subjects

Weiliang Zhu

Publications and source records attributed to Weiliang Zhu.

At least 19 recordsLinked to original sources

PubMed2007

Discovering potassium channel blockers from synthetic compound database by using structure-based virtual screening in conjunction with electrophysiological assay.

Potassium ion (K+) channels are attractive targets for drug discovery because of the essential roles played in biological systems. However, high-throughput screening (HTS) cannot be used to screen K+ channel blockers. To overcome this disadvantage of HTS, we have developed a virtual screening approach for discovering novel blockers of K+ channels. On the basis of a three-dimensional model of the eukaryotic K+ channels, molecular docking-based virtual screening was employed to search the chemical database MDL Available Chemicals Directory (ACD). Compounds were ranked according to their relative binding energy, favorable shape complementarity, and potential to form hydrogen bonds with the outer mouth of the K+ channel model. Twenty candidate compounds selected from the virtual screening were examined using the whole-cell voltage-clamp recording in rat dissociated hippocampal neurons. Among them, six compounds (5, 6, 8, 18-20) potently blocked both the delayed rectifier (IK) and fast transient K+ currents (IA). When applied externally, these six compounds preferentially blocked IK with potencies 2- to 500-fold higher than that of tetraethylammonium chloride. Intracellular application of the six compounds had no effect on both K+ currents. In addition, the interaction models and binding free energy calculations demonstrated that hydrophobic interaction and solvent effects play important roles in the inhibitory activities of these compounds. The results demonstrated that structure-based computer screening strategy could be used to identify novel, structurally diverse compounds targeting the pore binding pocket of the outer mouth of voltage-gated K+ channels. This study provides an alternative way of finding new blockers of voltage-gated K+ channels, while the techniques for high-throughput screening of K+ channel drugs remain in development.

Animals↗

PubMed2006

QSAR analyses on avian influenza virus neuraminidase inhibitors using CoMFA, CoMSIA, and HQSAR.

The recent wide spreading of the H5N1 avian influenza virus (AIV) in Asia, Europe and Africa and its ability to cause fatal infections in human has raised serious concerns about a pending global flu pandemic. Neuraminidase (NA) inhibitors are currently the only option for treatment or prophylaxis in humans infected with this strain. However, drugs currently on the market often meet with rapidly emerging resistant mutants and only have limited application as inadequate supply of synthetic material. To dig out helpful information for designing potent inhibitors with novel structures against the NA, we used automated docking, CoMFA, CoMSIA, and HQSAR methods to investigate the quantitative structure-activity relationship for 126 NA inhibitors (NIs) with great structural diversities and wide range of bioactivities against influenza A virus. Based on the binding conformations discovered via molecular docking into the crystal structure of NA, CoMFA and CoMSIA models were successfully built with the cross-validated q (2) of 0.813 and 0.771, respectively. HQSAR was also carried out as a complementary study in that HQSAR technique does not require 3D information of these compounds and could provide a detailed molecular fragment contribution to the inhibitory activity. These models also show clearly how steric, electrostatic, hydrophobicity, and individual fragments affect the potency of NA inhibitors. In addition, CoMFA and CoMSIA field distributions are found to be in well agreement with the structural characteristics of the corresponding binding sites. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing novel potent NA inhibitors.

Amino Acid Sequence↗

PubMed2006

Detailed mechanism for AmtB conducting NH4+/NH3: molecular dynamics simulations.

The mechanism by which the ammonium transporter, AmtB, conducts NH4+/NH3 into the cytoplasm was investigated using conventional molecular dynamics (MD) simulations. These simulations revealed that the neutral molecule, NH3, passes automatically through the channel upon its arrival at the Am2 site and that the function of the channel as a one-way valve for passage of NH3 is not determined by the cytoplasmic exit gate but, rather, by the periplasmic entrance gate of the channel. The NH3, produced by deprotonation of NH4+ at the entrance gate, is spontaneously conveyed to the central region of the channel via a hydrogen-bond network comprising His-168, His-318, Tyr-32, and the NH3 molecule. Finally, the NH3 molecule exits the channel through the exit gate formed by Phe-31, Ile-266, Val-314, and His-318. In addition, Ser-263 is shown to play a critical role in the final stages, acting as a pivoting arm to shunt the NH3 molecule from the cytoplasmic exit gate of the channel out into the cytoplasm. This finding is further supported by another simulation which shows that NH3 fails to be translocated through the channel formed by the Ser-263-Ala mutation. Thus, this study casts new insights on the mechanism of AmtB-mediated passage of NH3 across cellular membranes.

Cation Transport Proteins↗

PubMed2006

Induction of an aromatic six-membered nitrogen ring via cation-pi interaction.

Nitrogen clusters have been intensively studied for their potential application as high-energy density materials, but a six-membered nitrogen ring (N6) was not found to be stable and aromatic. To explore the possibility of inducing an aromatic N6 ring via cation-pi interaction, quantum chemistry calculations were performed on the systems of Ca2N6, CaN6, CaN6(2-), N6, and N6(4-) at the B3LYP/6-311+G level. The optimized geometries reveal that the planar structure of the N6 ring is stable only in the Ca2N6 complex. The computed NBO and CHelpG charges demonstrate that the planar N6 moiety in the Ca2N6 complex is almost a 10pi-electron system. The predicted nucleus-independent chemical shift (NICS) values demonstrate that the N6 moiety is aromatic in comparison with the NICS values of benzene. The estimated enthalpy of formation for the Ca2N6 complex is 100.4 kcal/mol for the reaction of 2Ca and 3N2. The binding energy between the Ca2+ cation and the N6(4-) moiety is -1928.8 kcal/mol, with electrostatic interaction serving as the predominant component. When all the calculated results are taken into account, including the planar structure, 10pi-electron system, identical bond length, and negative NICS value of the N6(4-) moiety in the Ca2N6 complex, it is deduced that the alkaline earth metal Ca is capable of inducing an aromatic N6 ring through the cation-pi interaction formed by electron transfer from the Ca atom to the N6 ring.

Benzene↗

PubMed2006

Blocking of the nicotinic acetylcholine receptor ion channel by chlorpromazine, a noncompetitive inhibitor: A molecular dynamics simulation study.

A large series of pharmacological agents, distinct from the typical competitive antagonists, block in a noncompetitive manner the permeability response of the nicotinic acetylcholine receptor (nAChR) to the neurotransmitter acetylcholine. Taking the neuroleptic chlorpromazine (CPZ) as an example of such agents, the blocking mechanism of noncompetitive inhibitors to the ion channel pore of the nAChR has been explored at the atomic level using both conventional and steered molecular dynamics (MD) simulations. Repeated steered MD simulations have permitted calculation of the free energy (approximately 36 kJ/mol) of CPZ binding and identification of the optimal site in the region of the serine and leucine rings, at approximately 4 A from the pore entrance. Coulomb and the Lennard-Jones interactions between CPZ and the ion channel as well as the conformational fluctuations of CPZ were examined to assess the contribution of each to the binding of CPZ to the nAChR. The MD simulations disclose a dynamic interaction of CPZ binding to the nAChR ionic channel. The cationic ammonium head of CPZ forms strong hydrogen bonds with Glu262 (alpha), Asp268 (beta), Glu272 (beta), Ser276 (beta), Glu280 (delta), Gln271 (gamma), Glu275 (gamma), and Asn279 (gamma) nAChR residues. Finally, the conventional MD simulation of CPZ at its identified binding site demonstrates that the binding of CPZ not only blocks ion transport through the channel but also markedly inhibits the conformational transitions of the channel, necessary for nAChR to carry out its biological function.

Acetylcholine↗

PubMed2006

Binding interaction of quercetin-3-beta-galactoside and its synthetic derivatives with SARS-CoV 3CL(pro): structure-activity relationship studies reveal salient pharmacophore features.

The 3C-like protease (3CL(pro)) of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) is one of the most promising targets for discovery of drugs against SARS, because of its critical role in the viral life cycle. In this study, a natural compound called quercetin-3-beta-galactoside was identified as an inhibitor of the protease by molecular docking, SPR/FRET-based bioassays, and mutagenesis studies. Both molecular modeling and Q189A mutation revealed that Gln189 plays a key role in the binding. Furthermore, experimental evidence showed that the secondary structure and enzymatic activity of SARS-CoV 3CL(pro) were not affected by the Q189A mutation. With the help of molecular modeling, eight new derivatives of the natural product were designed and synthesized. Bioassay results reveal salient features of the structure-activity relationship of the new compounds: (1) removal of the 7-hydroxy group of the quercetin moiety decreases the bioactivity of the derivatives; (2) acetoxylation of the sugar moiety abolishes inhibitor action; (3) introduction of a large sugar substituent on 7-hydroxy of quercetin can be tolerated; (4) replacement of the galactose moiety with other sugars does not affect inhibitor potency. This study not only reveals a new class of compounds as potential drug leads against the SARS virus, but also provides a solid understanding of the mechanism of inhibition against the target enzyme.

Binding Sites↗

PubMed2006

Conformational transition pathway in the allosteric process of human glucokinase.

Glucokinase (GK) is an important enzyme for regulating blood glucose levels and a potentially attractive target for diabetes of the young type 2 and persistent hyperinsulinemic hypoglycemia of infancy. To characterize the conformational transition of GK from the closed state to the superopen state, a series of conventional molecular dynamics (MD) and target MD (TMD) simulations were performed on both the wild-type enzyme and its mutants. Two 10-ns conventional MD simulations showed that, although the allosteric site of GK is approximately 20 A away from the active site, the activator is able to enhance the activity of the enzyme through conformational restriction. Fourteen TMD simulations on GK and five of its mutants revealed a reliably conformational transition pathway. The overall conformational transition includes three stages, and three likely stable intermediate states were identified by free energy scanning for the snapshots throughout the pathway. The conformational transition feature revealed by our TMD simulations rationalized several important mutagenesis and kinetic data. Remarkably, the TMD simulations predicted that Y61S, I159A, A201R, V203E, and V452S mutations, which have not been investigated so far, may facilitate the opening process of GK. These predictions also have been verified by mutagenesis and kinetic analyses in this study. These observations are beneficial to understanding the mechanism of GK regulation and designing the compounds for treating metabolic diseases.

Allosteric Regulation↗

PubMed2006

Towards discovering dual functional inhibitors against both wild type and K103N mutant HIV-1 reverse transcriptases: molecular docking and QSAR studies on 4,1-benzoxazepinone analogues.

To find useful information for discovering dual functional inhibitors against both wild type (WT) and K103N mutant reverse transcriptases (RTs) of HIV-1, molecular docking and 3D-QSAR approaches were applied to a set of twenty-five 4,1-benzoxazepinone analogues of efavirenz (SUSTIVA), some of them are active against the two RTs. 3D-QSAR models were constructed, based on their binding conformations determined by molecular docking, with r(2)(cv) values ranging from 0.656 to 0.834 for CoMFA and CoMSIA, respectively. The models were then validated to be highly predictive and extrapolative by inhibitors in two test sets with different molecular skeletons. Furthermore, CoMFA models were found to be well matched with the binding sites of both WT and K103N RTs. Finally, a reasonable pharmacophore model of 4,1-benzoxazepinones were established. The application of the model not only successfully differentiated the experimentally determined inhibitors from non-inhibitors, but also discovered two potent inhibitors from the compound database SPECS. On the basis of both the 3D-QSAR and pharmacophore models, new clues for discovering and designing potent dual functional drug leads against HIV-1 were proposed: (i) adopting positively charged aliphatic group at the cis-substituent of C3; (ii) reducing the electronic density at the position of O4; (iii) positioning a small branched aliphatic group at position of C5; (iv) using the negatively charged bulky substituents at position of C7.

Drug Design↗

PubMed2006

Syntheses of triazole-modified zanamivir analogues via click chemistry and anti-AIV activities.

Sixteen novel 4-triazole-modified zanamivir (1) analogues were synthesized using the click reactions, and their inhibitory activities against avian influenza virus (AIV, H5N1) were determined. Compound 3b exerts promising inhibitory activity with EC(50) of 6.4 microM, which is very close to that of zanamivir (EC(50) = 2.8 microM). Molecular modeling provided the information about the binding model between inhibitors and neuraminidase, which are in good agreement with inhibitory activities.

Antiviral Agents↗

PubMed2006

Mutagenic probability estimation of chemical compounds by a novel molecular electrophilicity vector and support vector machine.

MOTIVATION: Mutagenicity is among the toxicological end points that pose the highest concern. The accelerated pace of drug discovery has heightened the need for efficient prediction methods. Currently, most available tools fall short of the desired degree of accuracy, and can only provide a binary classification. It is of significance to develop a discriminative and informative model for the mutagenicity prediction. RESULTS: Here we developed a mutagenic probability prediction model addressing the problem, based on datasets covering a large chemical space. A novel molecular electrophilicity vector (MEV) is first devised to represent the structure profile of chemical compounds. An extended support vector machine (SVM) method is then used to derive the posterior probabilistic estimation of mutagenicity from the MEVs of the training set. The results show that our model gives a better performance than TOPKAT (http://www.accelrys.com) and other previously published methods. In addition, a confidence level related to the prediction can be provided, which may help people make more flexible decisions on chemical ordering or synthesis. AVAILABILITY: The binary program (ZGTOX_1.1) based on our model and samples of input datasets on Windows PC are available at http://dddc.ac.cn/adme upon request from the authors.

Algorithms↗

PubMed2006

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 http://www.dddc.ac.cn/tarfisdock/.

Databases, Protein↗

PubMed2006

How does ammonium dynamically interact with benzene in aqueous media? A first principle study using the Car-Parrinello molecular dynamics method.

The Car-Parrinello molecular dynamics (CPMD) method was used to study the dynamic characteristics of the cation-pi interaction between ammonium and benzene in gaseous and aqueous media. The results obtained from the CPMD calculation on the cation-pi complex in the gaseous state were very similar to those calculated from the Gaussian98 program with DFT and MP2 algorithms, demonstrating that CPMD is a valid approach for studying this system. Unlike the interaction in the gaseous state, our 12-ps CPMD simulation showed that the geometry of the complex in aqueous solution changes frequently in terms of the interaction angles and distances. Furthermore, the simulation revealed that the ammonium is constantly oscillating above the benzene plane in an aqueous environment and interacts with benzene mostly through three of its hydrogen atoms. In contrast, the interaction of the cation with the aromatic molecule in the gaseous state involves two hydrogen atoms. In addition, the free energy profile in aqueous solution was studied using constrained CPMD simulations, resulting in a calculated binding free energy of -5.75 kcal/mol at an optimum interaction distance of approximately 3.25 A, indicating that the cation-pi interaction between ammonium and benzene is stable even in aqueous solution. Thus, this CPMD study suggested that the cation-pi interaction between an ammonium (group) and an aromatic structure could take place even on surfaces of protein or nucleic acids in solution.

Benzene↗

PubMed2006

Neuraminidase pharmacophore model derived from diverse classes of inhibitors.

A three-dimensional pharmacophore model was developed based on 22 currently available inhibitors, which were carefully selected with great diversity in both molecular structure and bioactivity, for discovering new potent neuraminidase (NA) inhibitors to fight against avian influenza virus. The best hypothesis (Hypo1), consisting of five features, namely, one positive ionizable group, one negative ionizable group, one hydrophobic point, and two hydrogen-bond donors, has a correlation coefficient of 0.902, a root mean square deviation of 1.392, and a cost difference of 72.88, suggesting that a highly predictive pharmacophore model was successfully obtained. The application of the model shows great success in predicting the activities of 88 known NA inhibitors in our test set with a correlation coefficient of 0.818 with a cross-validation of 98% confidence level. Accordingly, our model should be reliable in identifying structurally diverse compounds with desired biological activity.

Enzyme Inhibitors↗

PubMed2006

Essential structural profile of a dual functional inhibitor against cyclooxygenase-2 (COX-2) and 5-lipoxygenase (5-LOX): molecular docking and 3D-QSAR analyses on DHDMBF analogues.

It is recently proposed that compounds with equal capabilities of inhibiting COX and 5-LOX, both are key enzymes involved in the arachidonic acid (AA) cascade, are expected to be safer non-steroidal anti-inflammatory drugs (NSAIDs). To dig out helpful information in designing dual functional inhibitors against the two enzymes, homology modeling, molecular dynamics (MD) simulations, automated docking, and 3D-QSAR analyses were performed in this study on 21 COX-2/5-LOX dual inhibitors, namely, 7-tert-butyl-2,3-dihydro-3,3-dimethylbenzofuran (DHDMBF) analogues. A 3D-model of 5-LOX was built based on the high-resolution X-ray structure of rabbit reticulocyte 15-lipoxygenase. Molecular docking was then applied to locate the binding orientations and conformations of DHDMBF analogues with COX-2 and 5-LOX, respectively, leading to highly predictive CoMFA models constructed on the basis of the binding conformations with q2 values of 0.782 and 0.634 for COX-2 and 5-LOX, respectively. In addition, CoMFA field distributions were found in good agreement with the structural characteristics of the corresponding binding sites. Both the docking simulations and QSAR analyses suggest that new potent dual inhibitors should share a structural feature with a moderately bulky group at R2 position and a rather negatively charged group around the position of the carbonyl group of DHDMBFs. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing new NSAIDs as anti-inflammation drugs with favorable safety profile.

Amino Acid Sequence↗

PubMed2006

An overall picture of SARS coronavirus (SARS-CoV) genome-encoded major proteins: structures, functions and drug development.

A severe atypical pneumonia designated as severe acute respiratory syndrome (SARS) by The World Health Organization broke out in China and menaced to more than other 30 countries between the end of the year 2002 and June of the year 2003. A novel coronavirus called severe acute respiratory syndrome coronavirus (SARS-CoV) has been recently identified as the etiological agent responsible for the infectious SARS disease. Based on extensively scientific cooperation and almost two-year's studies, remarkable achievements have been made in the understanding of the phylogenetic property and the genome organization of SARS-CoV, as well as the detailed characters of the major proteins involved in SARS-CoV life cycle. In this review, we would like to summarize the substantial scientific progress that has been made towards the structural and functional aspects of SARS-CoV associated key proteins. The progress focused on the corresponding key proteins' structure-based drug and vaccine developments has been also highlighted. The concerted and cooperative response for the treatment of the SARS disease has been proved to be a triumph of global public health and provides a new paradigm for the detection and control of future emerging infectious disease threats.

Antiviral Agents↗

PubMed2006

Molecular dynamics of nicotinic acetylcholine receptor correlating biological functions.

The nicotinic acetylcholine receptor (nAChR) that mediates fast intercellular communication in response to neurotransmitters is a paradigm of ligand-gated ion channels. Molecular dynamics (MD) simulations are valuable in understanding membrane protein function at atomic level, providing useful clues for further experimental/theoretical studies. In this brief review, recent progress in MD simulations of the nAChR has been illustrated, mainly focusing on the latest simulation of the whole transmembrane domain of the receptor. On the basis of MD simulations, asymmetrical and asynchronous motions of five subunits were observed both in the ligand binding and transmembrane domains; a closed-to-open conformational shift of the gate was captured in different simulation systems; the contributions from the lipid molecules and other transmembrane segments rather than M2 to the gate switch as well as the conformational change of the whole channel were assessed; the dynamic behavior and related physical/chemical properties of the water molecules and cations within the ion channel were examined; and an experimentally comparable single-channel conductance and ion selectivity were obtained.

Kinetics↗

PubMed2005

Dynamic mechanism of E2020 binding to acetylcholinesterase: a steered molecular dynamics simulation.

The unbinding process of E2020 ((R,S)-1-benzyl-4-[(5,6-dimethoxy-1-indanon)-2-yl]-methylpiperidine) leaving from the long active site gorge of Torpedo californica acetylcholinesterase (TcAChE) was studied by using steered molecular dynamics (SMD) simulations on a nanosecond scale with different velocities, and unbinding force profiles were obtained. Different from the unbinding of other AChE inhibitors, such as Huperzine A that undergoes the greatest barrier located at the bottleneck of the gorge, the major resistance preventing E2020 from leaving the gorge is from the peripheral anionic site where E2020 interacts intensively with several aromatic residues (e.g., Tyr70, Tyr121, and Trp279) through its benzene ring and forms a strong direct hydrogen bond and a water bridge with Ser286 via its O24. These interactions cause the largest rupture force, approximately 550 pN. It was found that the rotatable bonds of the piperidine ring to the benzene ring and dimethoxyindanone facilitate E2020 to pass the bottleneck through continuous conformation change by rotating those bonds to avoid serious conflict with Tyr121 and Phe330. The aromatic residues lining the gorge wall are the major components contributing to hydrophobic interactions between E2020 and TcAChE. Remarkably, these aromatic residues, acting in three groups as "sender" and "receiver", compose a "conveyer belt" for E2020 entering and leaving the TcAChE gorge.

Acetylcholinesterase↗

PubMed2005

Dynamic mechanism for the autophosphorylation of CheA histidine kinase: molecular dynamics simulations.

The two-component system (TCS) is an important signal transduction component for most bacteria. This signaling pathway is mediated by histidine kinases via autophosphorylation between P1 and P4 domains. Taking chemotaxis protein CheA as a model of TCS, the autophosphorylation mechanism of the TCS histidine kinases has been investigated in this study by using a computational approach integrated homology modeling, ligand-protein docking, protein-protein docking, and molecular dynamics (MD) simulations. Four nanosecond-scale MD simulations were performed on the free P4 domain, P4-ATP, P4-TNPATP, and P1-P4-ATP complexes, respectively. Upon its binding to the binding pocket of P4 with a folded conformation, ATP gradually extends to an open state with help from a water molecule. Meanwhile, ATP forms two hydrogen bonds with His413 and Lys494 at this state. Because of the lower energy of the folded conformations, ATP shrinks back to its folded conformations, leading to the rupture of the hydrogen bond between ATP and Lys494. Consequently, Lys494 moves away from the pocket entrance, resulting in an open of the ATP lid of P4. It is the open state of P4 that can bind tightly to P1, where the His45 of P1 occupies a favorable position for its autophosphorylation from ATP. This indicates that ATP is not only a phosphoryl group donor but also an activator for CheA phosphorylation. Accordingly, a mechanism of the autophosphorylation of CheA is proposed as that the ATP conformational switch triggers the opening of the ATP lid of P4, leading to P1 binding tightly, and subsequently autophosphorylation from ATP to P1.

Adenosine Triphosphate↗