Immunophilins and HIV-1 V3 Loop For Structure-Based Anti-AIDS Drug Design

Biological evaluation and interaction of a newly designed anti-cancer Pd(II) complex and human serum albumin

The pharmacokinetics and pharmacodynamics of any drug will depend, largely, on the interaction that has with human serum albumin (HSA), the most abundant plasma protein. The interaction between newly synthesized Pd(II) complexe, 2,2'-bipyridin Butylglycinato Pd(II) nitrate, an anti-tumor component, with HSA was studied at different temperatures by fluorescence, far UV circular dichroism (CD), UV-visible spectrophotometry and theoretical approaches. The Pd(II) complex has a strong ability to quench the intrinsic fluorescence of HSA through a dynamic quenching procedure. The binding parameters and thermodynamic parameters, including δH°, δS° and δG° were calculated by fluorescence quenching method, indicated that hydrophobic forces play a major role in the interaction of Pd(II) complex with HSA. Based on Autodock, FRET (fluorescence resonance energy transfer) and fluorescence quenching data, it may be concluded that one of the binding sites in the complex of HSA is near the only one Trp of HSA (Trp214) in sub domain IIA of the protein. Far-UV-CD results indicated that Pd(II)-complex induced increase in the α-helical content of the protein. The anti-tumor property of the synthesized Pd(II) complex was studied by testing it on human tumor cell line K562. The 50% cytotoxic concentration (Cc₅₀) of complex was determined using MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay. Also, fluorescence staining with DAPI (4,6-diamidino-2-phenylindole) revealed some typical nuclear changes that are characteristic of apoptosis which is induced at Cc₅₀ concentration of Pd(II) complex in K562 cell line after 24 h incubation. Our results suggest that Pd(II) complex is a promising anti-proliferative agent and should execute its biological effects by inducing apoptosis.

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca²⁺ dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 Å, 5.5 Å, 6 Å, 6.5 Å, 7.0 Å and 7.5 Å were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.

Biphasic association of T7 RNA polymerase and a nucleotide analogue, cibacron blue as a model to understand the role of initiating nucleotide in the mechanism of enzyme action

T7 RNA polymerase (T7 RNAP) is an enzyme that utilizes ribonucleotides to synthesize the nascent RNA chain in a template-dependent manner. Here we have studied the interaction of T7 RNAP with cibacron blue, an anthraquinone monochlorotriazine dye, its effect on the function of the enzyme and the probable mode of binding of the dye. We have used difference absorption spectroscopy and isothermal titration calorimetry to show that the dye binds T7 RNAP in a biphasic manner. The first phase of the binding is characterized by inactivation of the enzyme. The second binding site overlaps with the common substrate-binding site of the enzyme. We have carried out docking experiment to map the binding site of the dye in the promoter bound protein. Competitive displacement of the dye from the high affinity site by labeled GTP and isothermal titration calorimetry of high affinity GTP bound enzyme with the dye suggests a strong correlation between the high affinity dye binding and the high affinity GTP binding in T7 RNAP reported earlier from our laboratory.

In search of allosteric modulators of a7-nAChR by solvent density guided virtual screening

Nicotinic acetylcholine receptors (nAChR) are pentameric ligand gated ion channels whose activity can be modulated by endogenous neurotransmitters as well as by synthetic ligands that bind the same or distinct sites from the natural ligand. The subtype of α7 nAChR has been considered as a potenial therapeutic target for Alzheimer's disease, schizophrenia and other neurological and psychiatric disorders. Here we have developed a homology model of α7 nAChR based on two high resolution crystal structures with Brookhaven Protein Data Bank (PDB) codes 2QC1 and 2WN9 for threading on one monomer and then for building a pentamer, respectively. A number of small molecule binding sites are identified using Pocket Finder (J. An, M. Tortov, and R. Abagyan, Molecular & Cellular Proteomics, 4.6, 752-761 (2005)) of Internal Coordinate Mechanics (ICM). Remarkably, these computer-identified sites match perfectly with ordered solvent densities found in the high-resolution crystal structure of α1 nAChR, suggesting that the surface cavities in the α7 nAChR model are likely binding sites of small molecules. A high throughput virtual screening by flexible ligand docking of 5008 small molecule compounds was performed at three potential allosteric modulator (AM) binding sites of α7 nAChR using Molsoft ICM software (R. Abagyan, M. Tortov and D. Kuznetsov, J Comput Chem 15, 488-506, (1994)). Some experimentally verified allosteric modulators of α7 like CCMI comp-6, LY 7082101, 5-HI, TQS, PNU-120596, genistein, and NS-1738 ranked among top 100 compounds, while the rest of the compounds in the list could guide further search for new allosteric modulators.

Screening from the world's largest TCM database against H1N1 virus

The swine influenza virus (H1N1) 2009 pandemic highlights the importance of having effective anti-viral strategies. Recently, oseltamivir (Tamiflu) resistant influenza viruses are identified; which further emphasizes the urgency in developing new antiviral agents. In influenza virus replication cycle, viral surface glycoprotein, hemagglutinin, is responsible for viral entry into host cells. Hence, a potentially effective antiviral strategy is to inhibit viral entry mechanism. To develop novel antiviral agent that inhibits viral entry, we analyzed 20,000 traditional Chinese medicine (TCM) ingredients in hemagglutinin subtype H1 sialic acid binding site found on H1N1 virus. We then performed molecular dynamics simulations to investigate receptor-ligand interaction of the candidates obtained from docking. Here, we report three TCM derivatives that have high binding affinities to H1 sialic acid binding site residues based on structure-based calculations. The top three derivatives, xylopine_2, rosmaricine_14 and rosmaricine_15, all have an amine group that interact with Glu83 and a pyridinium group that interact with Asp103. Molecular dynamics simulations show that these derivatives form strong hydrogen bonding with Glu83 but interact transiently with Asp103. We therefore suggest that an enhanced hemagglutinin inhibitor, based on our scaffold, should be designed to bind both Glu83 and Asp103 with high affinity.

[Computer modeling of the promising inhibitors of the HIV-1 subtype A replication as a framework for the rational anti-aids drug design]

The model of the structural complex of cyclophilin B belonging to the immunophilins family with the HIV-1 subtype A V3 loop presenting the principal neutralizing determinant of the virus gp120 envelope protein as well as determinants of cell tropism and syncutium formation was generated by molecular docking methods. Basing on the conformational and energy characteristics of the built complex, computer-aided design of the polypeptide able to block effectively the functionally crucial V3 segments was implemented. From the joint analysis of the results derived with the data of literature, the generated molecule was suggested to offer a promising pharmacological substance for making a reality of the protein engineering projects aimed at developing the anti-AIDS drugs able to stop the HIV's spread.

Melting of model HIV-1 stem-loop 1 RNA dimers monitored by 2-aminopurine fluorescence

Viral maturation of HIV-1 involves refolding of its genomic RNA, which is believed to include a rearrangement of the SL1 stem-loop from a metastable conformation called kissing loop dimer (KD) to a stable one termed extended dimer (ED). To investigate this rearrangement in vitro we have studied the thermal melting of the RNA dimers formed by slightly modified 23-nucleotide SL1 RNA of HIV-1 Mal. Local structural changes in the RNA dimers during the melting were monitored by changes in the fluorescence of 2-aminopurine (2AP) incorporated in predetermined positions of RNA. We have shown that the stem regions of both preformed KD and ED melt in the temperature interval from 75 ° C to 90 ° C. Kissing loop interface of the KD RNA is found to be disrupted at lower temperatures from 20 ° C to 55 ° C, at which the stem regions remain intact. Conversion of the preformed KD to ED overcoming the kinetic barrier occurs between 55 ° C and 65 ° C. The melting of "loop-loop" regions in both preformed and newly formed EDs takes place around 70 ° C. Our finding that thermoinduced KD-to-ED conversion is preceded by transient dissociation of loop-loop interface disagrees with a common idea of strand exchange without disruption of loop-loop-contact.

The key residues of active sites on the catalytic fragment for paclitaxel interacting with poly (ADP-ribose) polymerase

Poly(ADP-ribose) polymerase (PARP) is regarded as a target protein for paclitaxel (PTX) to bind. An important issue is to identify the key residues as active sites for PTX interacting with PARP, which will help to understand the potential drug activity of PTX against cancer cells. Using docking method and MD simulation, we have constructed a refined structure of PTX docked on the catalytic function domain of PARP (PDB code: 1A26). The residues Glu327(988), Tyr246(907), Lys242(903), His165(826), Asp105(766), Gln102(763) and Gln98(759) in PARP are identified as potential sites involved in interaction with PTX according to binding energy (E(b)) between PTX and single residue calculated with B3LYP/6-31G(d,p). These residues form an active binding pocket located on the surface of the catalytic fragment, possibly interacting with the required groups of PTX leading to its activity against cancer cells. It is noted that most of the active sites make conatct with the "southern hemisphere" of PTX except for one residue, Tyr246(907), which interacts with the "northern hemisphere" of PTX. The conformation of PTX in complex with the catalytic fragment is observed as being T-shaped, similar to that complexed with β-tubulin. The total Eb of -269.9 kJ/mol represents the potent interaction between PTX and the catalytic fragment, implying that PTX can readily bind to the active pocket. The tight association of PTX with the catalytic fragment would inhibit PARP activation, suggesting a potential application of PTX as an effective antineoplastic agent.

Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database

The rapid spread of influenza virus subtype H1N1 poses a great threat to million lives worldwide. To search for new anti-influenza compounds, we performed molecular docking and molecular dynamics simulation to identify potential traditional Chinese medicine (TCM) constituents that could block influenza M2 channel activity. Quinic acid, genipin, syringic acid, cucurbitine, fagarine, and methyl isoferulate all have extremely well docking results as compared to control amantadine. Further de novo drug design suggests that derivatives of genipin and methyl isoferulate could have enhanced binding affinity towards M2 channel. Selected molecular dynamics simulations of M2-derivative complexes show stable hydrogen bond interactions between the derivatives and M2 residues, Ser10 and Ala9. To our best knowledge, this is the first study on the anti-viral activity of the above listed TCM compounds.

A conformationally constrained peptidomimetic binds to the extracellular region of HER2 protein

Human epidermal growth factor receptor 2 (HER2) is a member of the human epidermal growth factor receptor kinases (other members include EGFR or HER1, HER3, and HER4) that are involved in signaling cascades for cell growth and differentiation. It is well established that HER2-mediated heterodimerization has important implications in cancer. Deregulation of signaling pathways and overexpression of HER2 is known to occur in cancer cells, indicating a role of HER2 in tumorigenesis. Therefore, blocking HER2-mediated signaling has potential therapeutic value. We have designed several peptidomimetics to inhibit HER2-mediated signaling for cell growth. One of the compounds (HERP5, Arg-beta Naph-Phe) exhibited antiproliferative activity with IC(50) values in the micromolar-to-nanomolar range in breast cancer cell lines. Binding of fluorescently labeled HERP5 to HER2 protein was evaluated by fluorescence assay, microscopy, and circular dichroism spectroscopy. Results indicated that HERP5 binds to the extracellular region of the HER2 protein. Structure of the peptidomimetic HERP5 was studied by NMR and molecular dynamics simulations. Based on these results a model was proposed for HER2-EGFR dimerization and possible blocking by HERP5 peptidomimetic using a protein-protein docking method.

Anti-infectious bronchitis virus (IBV) activity of 1,8-cineole: effect on nucleocapsid (N) protein

In the present study, anti-IBV (infectious bronchitis virus) activity of 1,8-cineole was studied by MTT assay, as well as docking and molecular dynamic (MD) simulations. The CC50 of 1,8-cineole was above 10 mM. And the maximum noncytotoxic concentration (TD0) of 1,8-cineole was determined to be 3.90 ± 0.22 mM, which was much higher than that of ribavirin (0.78 ± 0.15 mM). 1,8-cineole could inhibit IBV with an IC(50) of 0.61 mM. MTT assay showed that the inhibition of IBV by 1, 8-cineole appears to occur moderately before entering the cell but much strongly after penetration of the virus into the cell. In silico simulations indicated that the binding site of 1,8-cineole was located at the N terminus of phosphorylated nucleocapsid (N) protein, with interaction energy equaling -40.33 kcal mol(-1). The residues TyrA92, ProA134, PheA137, AspA138 and TyrA140 had important roles during the binding process and are fully or partially conserved in various IBV strains. Based on spatial and energetic criteria, 1,8-cineole interfered with the binding between RNA and IBV N-protein. Results presented here may suggest that 1,8-cineole possesses anti-IBV properties, and therefore is a potential source of anti-IBV ingredients for the pharmaceutical industry.

Structure-based and ligand-based drug design for HER 2 receptor

Human epidermal growth factor receptor 2, HER2, is a commonly over-expressed tyrosine kinase receptor found in many types of carcinoma. Despite that there are several HER2 inhibitors, namely Iressa, Tarceva and Tykerb, currently in clinical trials, all can cause several side effects. In this study, both structure-based and ligand-based drug design were employed to design novel HER2 inhibitors from traditional Chinese medicine (TCM). The HER2 structure model was built in homology modeling based on known receptors of the same family. Docking and de novo evolution experiments were performed to identify candidates and to build derivatives. A training set of 32 compounds with inhibitory activities to HER2 was used to formulate the pharmacophore hypotheses that were subsequently used to examine candidates obtained from the docking study. Hydrogen bond interactions, salt-bridge formations and pi-stacking were observed between the ligands and Phe731, Lys753, Asp863 and Asp808 of HER2 protein. Combining results from both docking and pharmacophore mapping analysis, CLC015-5, CLC604-11 and CLC604-18 were well accepted and consistent in both approaches and were considered as the most potential HER2 inhibitors.

Virtual screening and drug design for PDE-5 receptor from traditional Chinese medicine database

Erectile dysfunction (ED) is a sexual disorder mainly caused by decrease in cellular concentration of cyclic guanosine monophosphate (cGMP), which is degraded by phosphodiesterase type-5 (PDE-5). As a potent therapeutic target, inhibitors such as Viagra , Cialis, and Levitra have already been developed to target PDE-5 for treating ED; traditional Chinese medicine, Epimedium sagittatum, also has shown prominent results as well. To developed new PDE-5 inhibitors, we performed a virtual screening of traditional Chinese medicine (TCM) database and docking analyses to identify candidates. Known PDE-5 inhibitors were used to construct a three dimensional quantitative structure-activity relationship (3D QSAR) model by HypoGen program. From docking analyses, isochlorogenic acid b was identified as the most potential inhibitory compound. De novo evolution designed 47 derivatives. Of the 47 derivatives, seven were able to map into the pharmacophore model, and these seven compounds were suggested to be the most promising leads for inhibiting PDE-5. An analysis of the hydrogen bond interactions formed between the docked ligands and PDE-5 identified ASN662, SER663 and GLN817 as the most frequently interacting residues. A total of eight novel leading compounds were identified to have favorable interaction with PDE-5. These compounds all had hydrogen bond interactions with three key residues that could be further investigated for understanding of PDE-5 and ligands interaction.