Docking studies suggest ligand-specific δ-opioid receptor conformations

An automated docking procedure was used to study binding of a series of δ-selective ligands to three models of the δ-opioid receptor. These models are thought to represent the three ligand-specific receptor conformations. Docking results are in agreement with point mutation studies and suggest that different ligands—agonists and antagonists—may bind to the same binding site under different receptor conformations. Docking to different receptor models (conformations) also suggests that by changing to a receptor-specific conformation, the receptor may open or close different binding sites to other ligands. Figure  Ligands 5 (green) and 6 (orange) in bindingpocket BP1 of the R1 δ-opioid receptor model     

Exploring the binding pocket for pyridopyrimidine ligands at the CCK1 receptor by molecular docking

Pyridopyrimidine-based analogues are among the most highly potent and selective antagonists of cholecystokinin receptor subtype-1 (CCK1R) described to date. To better understand the structural and chemical features responsible for the recognition mechanism, and to explore the binding pocket of these compounds, we performed automated molecular docking using GOLD2.2 software on some derivatives with structural diversity, and propose a putative binding conformation for each compound. The docking protocol was guided by the key role of the Asn333 residue, as revealed by site directed mutagenesis studies. The results suggest two putative binding modes located in the same pocket. Both are characterized by interaction with the main residues revealed by experiment, Asn333 and Arg336, and differ in the spatial position of the Boc-Trp moiety of these compounds. Hydrophobic contacts with residues Thr117, Phe107, Ile352 and Ile329 are also in agreement with experimental data. Despite the poor correlation obtained between the estimated binding energies and the experimental activity, the proposed models allow us to suggest a plausible explanation of the observed binding data in accordance with chemical characteristics of the compounds, and also to explain the observed diastereoselectivity of this family of antagonists towards CCK1R. The most reasonable selected binding conformations could be the starting point for future studies. Figure Superimposition of the two putative binding conformations revealed by molecular docking for pyridopyrimidine-based CCK1 antagonists     

The effect of a tightly bound water molecule on scaffold diversity in the computer-aided de novo ligand design of CDK2 inhibitors

We have determined the effects that tightly bound water molecules have on the de novo design of cyclin-dependent kinase-2 (CDK2) ligands. In particular, we have analyzed the impact of a specific structural water molecule on the chemical diversity and binding mode of ligands generated through a de novo structure-based ligand generation method in the binding site of CDK2. The tightly bound water molecule modifies the size and shape of the binding site and we have found that it also imposed constraints on the observed binding modes of the generated ligands. This in turn had the indirect effect of reducing the chemical diversity of the underlying molecular scaffolds that were able to bind to the enzyme satisfactorily.      

The vitamin riboflavin and its derivative lumichrome activate the LasR bacterial quorum-sensing receptor

Many bacteria use quorum sensing (QS) as an intercellular signaling mechanism to regulate gene expression in local populations. Plant and algal hosts, in turn, secrete compounds that mimic bacterial QS signals, allowing these hosts to manipulate QS-regulated gene expression in bacteria. Lumichrome, a derivative of the vitamin riboflavin, was purified and chemically identified from culture filtrates of the alga Chlamydomonas as a QS signal-mimic compound capable of stimulating the Pseudomonas aeruginosa LasR QS receptor. LasR normally recognizes the N-acyl homoserine lactone (AHL) signal, N-3-oxo-dodecanoyl homoserine lactone. Authentic lumichrome and riboflavin stimulated the LasR receptor in bioassays and lumichrome activated LasR in gel shift experiments. Amino acid substitutions in LasR residues required for AHL binding altered responses to both AHLs and lumichrome or riboflavin. These results and docking studies indicate that the AHL binding pocket of LasR recognizes both AHLs and the structurally dissimilar lumichrome or riboflavin. Bacteria, plants, and algae commonly secrete riboflavin or lumichrome, raising the possibility that these compounds could serve as either QS signals or as interkingdom signal mimics capable of manipulating QS in bacteria with a LasR-like receptor.     

Cancer-relevant biochemical targets of cytotoxic Lonchocarpus flavonoids: A molecular docking analysis

A molecular docking investigation has been carried out on cytotoxic prenylated flavonoids from Lonchocarpus haberi with cancer-relevant chemotherapeutic targets known to be inhibited by flavonoids. Two molecular docking programs, Molegro and ArgusDock, were used to compare the binding energies of Lonchocarpus flavonoids with other flavonoids, inhibitors, or known ligands, to aromatase (CYP 19), fatty acid synthase (FAS), xanthine oxidase (XO), cyclooxygenases (COX-1 and COX-2), lipoxygenase (LOX-3), ornithine decarboxylase (ODC), protein tyrosine kinase (PTK), phosphoinositide 3-kinase (PI3K), protein kinase C (PKC), topoisomerase II (ATP binding site), ATP binding cassette (ABC) transporter, and phospholipase A₂ (PLA). The Lonchocarpus flavonoids examined in this study exhibited docking energies comparable to or stronger than other flavonoids that had been previously shown to be effective inhibitors of these enzymes. Furthermore, prenylated flavonoids, such as the Lonchocarpus flavonoids and xanthohumol, generally showed greater binding energies than the non-prenylated flavonoids. We conclude, therefore, that the Lonchocarpus flavonoids possibly owe their cytotoxic activity by inhibition of one or more of these enzymes.      

Molecular modeling of the three-dimensional structure of GLP-1R and its interactions with several agonists

Glucagon-like peptide-1 receptor (GLP-1R) is a promising molecular target for developing drugs treating type 2 diabetes. We have predicted the complete three-dimensional structure of GLP-1R and the binding modes of several GLP-1R agonists, including GLP-1, Boc5, and Cpd1, through a combination of homology modeling, molecular docking, and long-time molecular dynamics simulation on a lipid bilayer. Our model can reasonably interpret the results of a number of mutation experiments regarding GLP-1R as well as the successful modification to GLP-1 by Liraglutide. Our model is also validated by a recently revealed crystal structure of the extracellular domain of GLP-1R. An activation mechanism of GLP-1R agonists is proposed based on the principal component analysis and normal mode analysis on our predicted GLP-1R structure. Before the complete structure of GLP-1R is determined through experimental means, our model may serve as a valuable reference for characterizing the interactions between GLP-1R and its agonists. Figure Comparison of our predicted model of rGLP-1R (left) with the recently revealed crystal structure of hGLP-1R (right)     

ONIOM DFT/PM3 calculation on the interaction between STI-571 and abelson tyrosine kinase

The ONIOM2 (B3LYP/6–31G (d, p): PM3) and B3LYP/6–31G (d, p) methods were applied to investigate the interaction between STI-571 and abelson tyrosine kinase binding site. The complex of N-[4-methyl-3-(4-pyridin-3-yl-pyrimidin-2-ylamino)- phenyl]-benzamide (part of STI-571) and related 16 amino acid residues were found at B3LYP/6–31G (d, p) level to have hydrogen bonds and π....π stacking interaction, their binding energy via HAF optimization was −20.4 kcal mol⁻¹. The results derived from this study agreed well with the reported observation. Figure Optimized structure of STI-571 and Thr315 in abelson tyrosine kinase based on ONIOM2 method     

Toward the virtual screening of Cdc25A phosphatase inhibitors with the homology modeled protein structure

Cdc25 phosphatases have been considered as attractive drug targets for anticancer therapy due to the correlation of their overexpression with a wide variety of cancers. As a method for the discovery of novel inhibitors of Cdc25 phosphatases, we have evaluated the computer-aided drug design protocol involving the homology modeling of Cdc25A and virtual screening with the two docking tools: FlexX and the modified AutoDock program implementing the effects of ligand solvation in the scoring function. The homology modeling with the X-ray crystal structure of Cdc25B as a template provides a high-quality structure of Cdc25A that enables the structure-based inhibitor design. Of the two docking programs under consideration, AutoDock is found to be more accurate than FlexX in terms of scoring putative ligands. A detailed binding mode analysis of the known inhibitors shows that they can be stabilized in the active site of Cdc25A through the simultaneous establishment of the multiple hydrogen bonds and the hydrophobic interactions. The present study demonstrates the usefulness of the modified AutoDock program as a docking tool for virtual screening of new Cdc25 phosphatase inhibitors as well as for binding mode analysis to elucidate the activities of known inhibitors. Figure Structures and available IC50 values (in μM) of the twenty Cdc25 phosphatase inhibitors seeded in docking library     

3D-QSAR studies of Dipeptidyl peptidase IV inhibitors using a docking based alignment

Dipeptidyl peptidase IV (DPP-IV) deactivates the incretin hormones GLP-1 and GIP by cleaving the penultimate proline or alanine from the N-terminal (P1-position) of the peptide. Inhibition of this enzyme will prevent the degradation of the incretin hormones and maintain glucose homeostasis; this makes it an attractive target for the development of drugs for diabetes. This paper reports 3D-QSAR analysis of several DPP-IV inhibitors, which were aligned by the receptor-based technique. The conformation of the molecules in the active site was obtained through docking methods. The QSAR models were generated on two training sets composed of 74 and 25 molecules which included phenylalanine, thiazolidine, and fluorinated pyrrolidine analogs. The 3D-QSAR models are robust with statistically significant r², q², and values. The CoMFA and CoMSIA models were used to design some new inhibitors with several fold higher binding affinity. Figure The CoMFA contours around molecule D1T155 (a) steric contours - favored (green); disfavored (yellow) (b) electrostatic contours - electropositive (blue); electronegative (red)     

PpoR is a conserved unpaired LuxR solo of Pseudomonas putida which binds N-acyl homoserine lactones

Background  

Only a small number of Pseudomonas putida strains possess the typical N-acyl homoserine lactone quorum sensing system (AHL QS) that consists of a modular LuxR family protein and its cognate LuxI homolog that produces the AHL signal. Moreover, AHL QS systems in P. putida strains are diverse in the type of AHLs they produce and the phenotypes that they regulate.     

The 3D structure of the defense-related rice protein Pir7b predicted by homology modeling and ligand binding studies

To better understand the ligand-binding mechanism of protein Pir7b, important part in detoxification of a pathogen-derived compound against Pyricularia oryzae, a 3D structure model of protein Pir7b was constructed based on the structure of the template SABP2. Three substrates were docking to this protein, two of them were proved to be active, and some critical residues are identified, which had not been confirmed by the experiments. His87 and Leu17 considered as ‘oxyanion hole’ contribute to initiating the Ser86 nucleophilic attack. Gln187 and Asp139 can form hydrogen bonds with the anilid group to maintain the active binding orientation with the substrates. The docking model can well interpret the specificity of protein Pir7b towards the anilid moiety of the substrates and provide valuable structure information about the ligand binding to protein Pir7b. Figure Ligand binding analysis based on the refined Pir7b model. Magenta dash line, hydrogen bond; Red dash line, distance label. (a) Docking of 2-naphthol AS-acetate to Pir7b model. A 3D figure of 2-naphthol AS-acetate-Pir7b complex is also attached (b) Docking of 2-naphthol AS-2-chlor-propionate to Pir7b model. (c) Docking of 2-naphthol-acetate to Pir7b model.     

Theoretical study of the intermolecular H-bonding and intermolecular proton transfer between isocytosine tautomeric forms and <Emphasis Type="Italic">R</Emphasis>,<Emphasis Type="Italic">S</Emphasis>-lactic acid

Eight H-bonded complexes between isocytosine (isoC) tautomeric forms and R/S-lactic acid (LA) have been studied at the B3LYP and HF levels of theory using 6–31+G(d) basis set. The energy barriers of the intermolecular proton transfers were also estimated as the results showed that they are several times lower than those of the intramolecular proton transfers of isoC in the gas phase. Furthermore, the energy barriers of the tautomerizations in which the carboxylic H-atom takes part are several times lower than those in which the LA OH group assists the proton transfer. Figure     

Computer-based design of novel HIV-1 entry inhibitors: neomycin conjugated to arginine peptides at two specific sites

Aminoglycoside–arginine conjugates (AAC and APAC) are multi-target inhibitors of human immunodeficiency virus type-1 (HIV-1). Here, we predict new conjugates of neomycin with two arginine peptide chains binding at specific sites on neomycin [poly-arginine-neomycin-poly-arginine (PA-Neo-PA)]. The rationale for the design of such compounds is to separate two short arginine peptides with neomycin, which may extend the binding region of the CXC chemokine receptor type 4 (CXCR4). We used homology models of CXCR4 and unliganded envelope glycoprotein 120 (HIV-1_IIIB gp120) and docked PA-Neo-PAs and APACs to these using a multistep docking procedure. The results indicate that PA-Neo-PAs spread over two negatively charged patches of CXCR4. PA-Neo-PA–CXCR4 complexes are energetically more favorable than AACs/APAC–CXCR4 complexes. Notably, our CXCR4 model and docking procedure can be applied to predict new compounds that are either inhibitors of gp120–CXCR4 binding without affecting stromal cell-derived factor 1α (SDF-1α) chemotaxis activity, or inhibitors of SDF-1α–CXCR4 binding resulting in an anti-metastasis effect. We also predict that PA-Neo-PAs and APACs can interfere with CD4–gp120 binding in unliganded conformation. Figure The r5-Neo-r5-CXCR4 complex. CXCR4 is shown in CPK representation. The negatively charged residues are shown in red and positively charged residues in blue. The r5-Neo-r5 is shown in stick representation, neomycin core is colored yellow and arginine moieties are colored magenta. Two negatively charged patches separated by neutral and positively charged residues are visible.     

Structure-based 3D-QSAR studies on thiazoles as 5-HT<Subscript>3</Subscript> receptor antagonists

Structure-based 3D-QSAR studies were performed on 20 thiazoles against their binding affinities to the 5-HT₃ receptor with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The thiazoles were initially docked into the binding pocket of a human 5-HT_3A receptor homology model, constructed on the basis of the crystal structure of the snail acetylcholine binding protein (AChBP), using the GOLD program. The docked conformations were then extracted and used to build the 3D-QSAR models, with cross-validated values 0.785 and 0.744 for CoMFA and CoMSIA, respectively. An additional five molecules were used to validate the models further, giving satisfactory predictive values of 0.582 and 0.804 for CoMFA and CoMSIA, respectively. The results would be helpful for the discovery of new potent and selective 5-HT₃ receptor antagonists.      

N‐ACYL HOMOSERINE LACTONe LACTONASE,AiiA, INACTIVATION OF QUORUM‐SENSING AGONISTS PRODUCED BY CHLAMYDOMONAS REINHARDTII (CHLOROPHYTA) AND CHARACTERIZATION OF aiiA TRANSGENIC ALGAE1

Eukaryotes such as plants and the unicellular green alga Chlamydomonas reinhardtii P. A. Dang. produce and secrete compounds that mimic N‐acyl homoserine lactone (AHL) bacterial quorum‐sensing (QS) signals and alter QS‐regulated gene expression in the associated bacteria. Here, we show that the set of C. reinhardtii signal‐mimic compounds that activate the CepR AHL receptor of Burkholderia cepacia are susceptible to inactivation by AiiA, an AHL lactonase enzyme of Bacillus. Inactivation of these algal mimics by AiiA suggests that the CepR‐stimulatory class of mimics produced by C. reinhardtii may have a conserved lactone ring structure in common with AHL QS signals. To examine the role of AHL mimic compounds in the interactions of C. reinhardtii with bacteria, the aiiA gene codon optimized for Chlamydomonas was generated for the expression of AiiA as a chimeric fusion with cyan fluorescent protein (AimC). Culture filtrates of transgenic strains expressing the fusion protein AimC had significantly reduced levels of CepR signal‐mimic activities. When parental and transgenic algae were cultured with a natural pond water bacterial community, a morphologically distinct, AHL‐producing isolate of Aeromonas veronii was observed to colonize the transgenic algal cultures and form biofilms more readily than the parental algal cultures, indicating that secretion of the CepR signal mimics by the alga can significantly affect its interactions with bacteria it encounters in natural environments. The parental alga was also able to sequester and/or destroy AHLs in its growth media to further disrupt or manipulate bacterial QS.