Modelling the structure of latexin-carboxypeptidase A complex based on chemical cross-linking and molecular docking

We have determined the three-dimensional structure of the protein complex between latexin and carboxypeptidase A using a combination of chemical cross-linking, mass spectrometry and molecular docking. The locations of three intermolecular cross-links were identified using mass spectrometry and these constraints were used in combination with a speed-optimised docking algorithm allowing us to evaluate more than 3 x 10(11) possible conformations. While cross-links represent only limited structural constraints, the combination of only three experimental cross-links with very basic molecular docking was sufficient to determine the complex structure. The crystal structure of the complex between latexin and carboxypeptidase A4 determined recently allowed us to assess the success of this structure determination approach. Our structure was shown to be within 4 angstrom r.m.s. deviation of C alpha atoms of the crystal structure. The study demonstrates that cross-linking in combination with mass spectrometry can lead to efficient and accurate structural modelling of protein complexes.

Peptide binding prediction for the human class II MHC allele HLA-DP2:a molecular docking approach

MHC class II proteins bind oligopeptide fragments derived from proteolysis of pathogen antigens, presenting them at the cell surface for recognition by CD4+ T cells. Human MHC class II alleles are grouped into three loci: HLA-DP, HLA-DQ and HLA-DR. In contrast to HLA-DR and HLA-DQ, HLA-DP proteins have not been studied extensively, as they have been viewed as less important in immune responses than DRs and DQs. However, it is now known that HLA-DP alleles are associated with many autoimmune diseases. Quite recently, the X-ray structure of the HLA-DP2 molecule (DPA*0103, DPB1*0201) in complex with a self-peptide derived from the HLA-DR a-chain has been determined. In the present study, we applied a validated molecular docking protocol to a library of 247 modelled peptide-DP2 complexes, seeking to assess the contribution made by each of the 20 naturally occurred amino acids at each of the nine binding core peptide positions and the four flanking residues (two on both sides).

Peptidoglycan derived from Staphylococcus epidermidis induces Connexin43 hemichannel activity with consequences on the innate immune response in endothelial cells

Gram-positive bacterial cell wall components including PGN (peptidoglycan) elicit a potent pro-inflammatory response in diverse cell types, including endothelial cells, by activating TLR2 (Toll-like receptor 2) signalling. The functional integrity of the endothelium is under the influence of a network of gap junction intercellular communication channels composed of Cxs (connexins) that also form hemichannels, signalling conduits that are implicated in ATP release and purinergic signalling. PGN modulates Cx expression in a variety of cell types, yet effects in endothelial cells remain unresolved. Using the endothelial cell line b.End5, a 6 h challenge with PGN induced IL-6 (interleukin 6), TLR2 and Cx43 mRNA expression that was associated with enhanced Cx43 protein expression and gap junction coupling. Cx43 hemichannel activity, measured by ATP release from the cells, was induced following 15 min of exposure to PGN. Inhibition of hemichannel activity with carbenoxolone or apyrase prevented induction of IL-6 and TLR2 mRNA expression by PGN, but had no effect on Cx43 mRNA expression levels. In contrast, knockdown of TLR2 expression had no effect on PGN-induced hemichannel activity, but reduced the level of TLR2 and Cx43 mRNA expression following 6 h of PGN challenge. PGN also acutely induced hemichannel activity in HeLa cells transfected to express Cx43, but had no effect in Cx43-deficient HeLa OHIO cells. All ATP responses were blocked with Cx-specific channel blockers. We conclude that acute Cx43 hemichannel signalling plays a role in the initiation of early innate immune responses in the endothelium.

EpiDOCK:a molecular docking-based tool for MHC class II binding prediction

Cellular peptide vaccines contain T-cell epitopes. The main prerequisite for a peptide to act as a T-cell epitope is that it binds to a major histocompatibility complex (MHC) protein. Peptide MHC binder identification is an extremely costly experimental challenge since human MHCs, named human leukocyte antigen, are highly polymorphic and polygenic. Here we present EpiDOCK, the first structure-based server for MHC class II binding prediction. EpiDOCK predicts binding to the 23 most frequent human, MHC class II proteins. It identifies 90% of true binders and 76% of true non-binders, with an overall accuracy of 83%. EpiDOCK is freely accessible at http://epidock.ddg-pharmfac. net. © The Author 2013. Published by Oxford University Press. All rights reserved.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Peptide binding to HLA-DP proteins at pH 5.0 and pH 7.0:a quantitative molecular docking study

Background: HLA-DPs are class II MHC proteins mediating immune responses to many diseases. Peptides bind MHC class II proteins in the acidic environment within endosomes. Acidic pH markedly elevates association rate constants but dissociation rates are almost unchanged in the pH range 5.0 - 7.0. This pH-driven effect can be explained by the protonation/deprotonation states of Histidine, whose imidazole has a pKa of 6.0. At pH 5.0, imidazole ring is protonated, making Histidine positively charged and very hydrophilic, while at pH 7.0 imidazole is unprotonated, making Histidine less hydrophilic. We develop here a method to predict peptide binding to the four most frequent HLA-DP proteins: DP1, DP41, DP42 and DP5, using a molecular docking protocol. Dockings to virtual combinatorial peptide libraries were performed at pH 5.0 and pH 7.0. Results: The X-ray structure of the peptide - HLA-DP2 protein complex was used as a starting template to model by homology the structure of the four DP proteins. The resulting models were used to produce virtual combinatorial peptide libraries constructed using the single amino acid substitution (SAAS) principle. Peptides were docked into the DP binding site using AutoDock at pH 5.0 and pH 7.0. The resulting scores were normalized and used to generate Docking Score-based Quantitative Matrices (DS-QMs). The predictive ability of these QMs was tested using an external test set of 484 known DP binders. They were also compared to existing servers for DP binding prediction. The models derived at pH 5.0 predict better than those derived at pH 7.0 and showed significantly improved predictions for three of the four DP proteins, when compared to the existing servers. They are able to recognize 50% of the known binders in the top 5% of predicted peptides. Conclusions: The higher predictive ability of DS-QMs derived at pH 5.0 may be rationalised by the additional hydrogen bond formed between the backbone carbonyl oxygen belonging to the peptide position before p1 (p-1) and the protonated ε-nitrogen of His 79β. Additionally, protonated His residues are well accepted at most of the peptide binding core positions which is in a good agreement with the overall negatively charged peptide binding site of most MHC proteins. © 2012 Patronov et al.; licensee BioMed Central Ltd.

MHC class II binding prediction by molecular docking

Proteins of the Major Histocompatibility Complex (MHC) bind self and nonself peptide antigens or epitopes within the cell and present them at the cell surface for recognition by T cells. All T-cell epitopes are MHC binders but not all MCH binders are T-cell epitopes. The MHC class II proteins are extremely polymorphic. Polymorphic residues cluster in the peptide-binding region and largely determine the MHC's peptide selectivity. The peptide binding site on MHC class II proteins consist of five binding pockets. Using molecular docking, we have modelled the interactions between peptide and MHC class II proteins from locus DRB1. A combinatorial peptide library was generated by mutation of residues at peptide positions which correspond to binding pockets (so called anchor positions). The binding affinities were assessed using different scoring functions. The normalized scoring functions for each amino acid at each anchor position were used to construct quantitative matrices (QM) for MHC class II binding prediction. Models were validated by external test sets comprising 4540 known binders. Eighty percent of the known binders are identified in the best predicted 15% of all overlapping peptides, originating from one protein. © 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Docking and simulation study of the enantiospecificity of chiral epoxidation reactions catalyzed by chloroperoxidase

Chloroperoxidase (CPO), a 298-residue glycosylated protein from the fungus Caldariomyces fumago, is probably the most versatile heme enzyme yet discovered. Interest in CPO as a catalyst is based on its power to produce enantiomerically enriched products. Recent research has focused its attention on the ability of CPO to epoxidize alkenes in high regioselectivity and enantioselectivity as an efficient and environmentally benign alternative to traditional synthetic routes. There has been little work on the nature of ligand binding, which probably controls the regio- and enantiospecifity of CPO. Consequently it is here that we focus our work. We report docking calculations and computer simulations aimed at predicting the enantiospecificity of CPO-catalyzed epoxidation of three model substrates. On the basis of this work candidate mutations to improve the efficiency of CPO are predicted. In order to accomplish these aims, a simulated annealing and molecular dynamics protocol is developed to sample potentially reactive substrate/CPO complexes.

Molecular docking and geographical information systems as tools to assess the potential impact of veterinary medicines on non-target organisms and the environment

Veterinary medicines (VMs) from agricultural industry can enter the environment in a number of ways. This includes direct exposure through aquaculture, accidental spillage and disposal, and indirect entry by leaching from manure or runoff after treatment. Many compounds used in animal treatments have ecotoxic properties that may have chronic or sometimes lethal effects when they come into contact with non-target organisms. VMs enter the environment in mixtures, potentially having additive effects. Traditional ecotoxicology tests are used to determine the lethal and sometimes reproductive effects on freshwater and terrestrial organisms. However, organisms used in ecotoxicology tests can be unrepresentative of the populations that are likely to be exposed to the compound in the environment. Most often the tests are on single compound toxicity but mixture effects may be significant and should be included in ecotoxicology testing. This work investigates the use, measured environmental concentrations (MECs) and potential impact of sea lice treatments on salmon farms in Scotland. Alternative methods for ecotoxicology testing including mixture toxicity, and the use of in silico techniques to predict the chronic impact of VMs on different species of aquatic organisms were also investigated. The Scottish Environmental Protection Agency (SEPA) provided information on the use of five sea lice treatments from 2008-2011 on Scottish salmon farms. This information was combined with the recently available data on sediment MECs for the years 2009-2012 provided by SEPA using ArcGIS 10.1. In depth analysis of this data showed that from a total of 55 sites, 30 sites had a MEC higher than the maximum allowable concentration (MAC) as set out by SEPA for emamectin benzoate and 7 sites had a higher MEC than MAC for teflubenzuron. A number of sites that were up to 16 km away from the nearest salmon farm reported as using either emamectin benzoate or teflubenzuron measured positive for the two treatments. There was no relationship between current direction and the distribution of the sea lice treatments, nor was there any evidence for alternative sources of the compounds e.g. land treatments. The sites that had MECs higher than the MAC could pose a risk to non-target organisms and disrupt the species dynamics of the area. There was evidence that some marine protected sites might be at risk of exposure to these compounds. To complement this work, effects on acute mixture toxicity of the 5 sea lice treatments, plus one major metabolite 3-phenoxybenzoic acid (3PBA), were measured using an assay using the bioluminescent bacteria Aliivibrio fischeri. When exposed to the 5 sea lice treatments and 3PBA A. fischeri showed a response to 3PBA, emamectin benzoate and azamethiphos as well as combinations of the three. In order to establish any additive effect of the sea lice treatments, the efficacy of two mixture prediction equations, concentration addition (CA) and independent action ii(IA) were tested using the results from single compound dose response curves. In this instance IA was the more effective prediction method with a linear regression confidence interval of 82.6% compared with 22.6% of CA. In silico molecular docking was carried out to predict the chronic effects of 15 VMs (including the five used as sea lice control). Molecular docking has been proposed as an alternative screening method for the chronic effects of large animal treatments on non-target organisms. Oestrogen receptor alpha (ERα) of 7 non-target bony fish and the African clawed frog Xenopus laevis were modelled using SwissModel. These models were then ‘docked’ to oestradiol, the synthetic oestrogen ethinylestradiol, two known xenoestrogens dichlorodiphenyltrichloroethane (DDT) and bisphenol A (BPA), the antioestrogen breast cancer treatment tamoxifen and 15 VMs using Auto Dock 4. Based on the results of this work, four VMs were identified as being possible xenoestrogens or anti-oestrogens; these were cypermethrin, deltamethrin, fenbendazole and teflubenzuron. Further investigation, using in vitro assays, into these four VMs has been suggested as future work. A modified recombinant yeast oestrogen screen (YES) was attempted using the cDNA of the ERα of the zebrafish Danio rerio and the rainbow trout Oncorhynchus mykiss. Due to time and difficulties in cloning protocols this work was unable to be completed. Use of such in vitro assays would allow for further investigation of the highlighted VMs into their oestrogenic potential. In conclusion, VMs used as sea lice treatments, such as teflubenzuron and emamectin benzoate may be more persistent and have a wider range in the environment than previously thought. Mixtures of sea lice treatments have been found to persist together in the environment, and effects of these mixtures on the bacteria A. fischeri can be predicted using the IA equation. Finally, molecular docking may be a suitable tool to predict chronic endocrine disrupting effects and identify varying degrees of impact on the ERα of nine species of aquatic organisms.

Docking of HIV-1 Vpr to the nuclear envelope is mediated by the interaction with the nucleoporin hCG1

The HIV-1 genome contains several genes coding for auxiliary proteins, including the small Vpr protein. Vpr affects the integrity of the nuclear envelope and participates in the nuclear translocation of the preintegration complex containing the viral DNA. Here, we show by photobleaching experiments performed on living cells expressing a Vpr-green fluorescent protein fusion that the protein shuttles between the nucleus and the cytoplasm, but a significant fraction is concentrated at the nuclear envelope, supporting the hypothesis that Vpr interacts with components of the nuclear pore complex. An interaction between HIV-1 Vpr and the human nucleoporin CG1 (hCG1) was revealed in the yeast two-hybrid system, and then confirmed both in vitro and in transfected cells. This interaction does not involve the FG repeat domain of hCG1 but rather the N-terminal region of the protein. Using a nuclear import assay based on digitonin-permeabilized cells, we demonstrate that hCG1 participates in the docking of Vpr at the nuclear envelope. This association of Vpr with a component of the nuclear pore complex may contribute to the disruption of the nuclear envelope and to the nuclear import of the viral DNA.

Low molecular weight compounds from mushrooms as potential Bcl-2 inhibitors: docking and virtual screening studies

Mushrooms have the ability to promote apoptosis in tumor cell lines, but the mechanism of action is not quite well understood. Inhibition of the interaction between Bcl-2 and pro-apoptotic proteins could be an important step that leads to apoptosis. Therefore, the discovery of compounds with the ability to inhibit Bcl-2 is an ongoing research topic in drug discovery. In this study, we started by analyzing Bcl-2 experimental structures that are currently available in Protein Data Bank database. After analysis of the more relevant Bcl-2 structures, 4 were finally selected. An analysis of the best docking methodology was then performed using a cross-docking and re-docking approach while testing 2 docking softwares: AutoDock 4 and AutoDock Vina. Autodock4 provided the best docking results and was selected to perform a virtual screening study applied to a dataset of 40 Low Molecular Weight (LMW) compounds present in mushrooms, using the selected Bcl-2 structures as target. Results suggest that steroid are the more promising family, among the analyzed compounds, and may have the ability to interact with Bcl-2 and this way promoting tumor apoptosis. The steroids that presented lowest estimated binding energy (ΔG) were: Ganodermanondiol, Cerevisterol, Ganoderic Acid X and Lucidenic Lactone; with estimated ΔG values between -8,45 and -8,23 Kcal/mol. A detailed analysis of the docked conformation of these 4 top ranked LMW compounds was also performed and illustrates a plausible interaction between the 4 top raked steroids and Bcl-2, thus substantiating the accuracy of the predicted docked poses. Therefore, tumoral apoptosis promoted by mushroom might be related to Bcl-2 inhibition mediated by steroid family of compounds.

Investigating mushroom LMW compounds as potential Bcl-2 inhibitors: docking studies using AutoDock4

The B cell CLL/lymphoma-2 (Bcl-2) family is functionally classified as either anti-apoptotic or pro-apoptotic, and the regulation of its interactions dictates survival or commitment to apoptosis. Bcl-2 family is also implicated in a wide range of diseases. In some types of cancers, including lymphomas and epithelial cancers, protein overexpression of anti-apoptotic Bcl-2 family, such as the Bcl-2 protein is indicative of cancer in an advanced stage, with a poor prognosis and resistant to chemotherapy [1]. Several reports indicate that mushrooms have the ability to promote apoptosis in tumour cell lines, but the mechanism of action is not fully understood. Inhibition of the interaction between Bcl-2 (anti-apoptotic protein) and proapoptotic proteins could be an important step in the mechanism of mushroom induced apoptosis. Therefore, the discovery of compounds with the capacity to inhibit Bcl-2 is an ongoing research topic on cancer therapy. In this work, docking studies were performed using a dataset of 40 low molecular weight (LMW) compounds present in mushrooms. The docking software AutoDock 4 was used and docking studies were performed using 5 selected Bcl-2 crystal structures as targets. Compounds with the lowest predicted binding energy (predΔG) are expected to be the more potent inhibitors. Among the tested compounds, steroids presented the lowest predΔG with several exhibiting values below -9 kcal/mol. The results are corroborated by several reports that state that steroids induce apoptosis in several tumor cells. It is thus feasible that they might act by preventing Bcl-2 from forming complexes with the respective proapoptotic protein interaction partners, namely Bak, Bax, and Bim. Moreover, previous studies on our research group demonstrated that 48 h treatment of MCF-7 cells (breast carcinoma) with Suillus collinitus methanolic extract caused a decrease in Bcl-2, highlighting the antitumor potential of this mushroom species [2]. In conclusion, the process of apoptosis promoted by mushroom extracts may be related to the inhibition of Bcl-2 by the steroid derivatives herein studied. However, further studies are needed to confirm this hypothesis.

A New Multi-Objective Approach for Molecular Docking Based on RMSD and Binding Energy

Ligand-protein docking is an optimization problem based on predicting the position of a ligand with the lowest binding energy in the active site of the receptor. Molecular docking problems are traditionally tackled with single-objective, as well as with multi-objective approaches, to minimize the binding energy. In this paper, we propose a novel multi-objective formulation that considers: the Root Mean Square Deviation (RMSD) difference in the coordinates of ligands and the binding (intermolecular) energy, as two objectives to evaluate the quality of the ligand-protein interactions. To determine the kind of Pareto front approximations that can be obtained, we have selected a set of representative multi-objective algorithms such as NSGA-II, SMPSO, GDE3, and MOEA/D. Their performances have been assessed by applying two main quality indicators intended to measure convergence and diversity of the fronts. In addition, a comparison with LGA, a reference single-objective evolutionary algorithm for molecular docking (AutoDock) is carried out. In general, SMPSO shows the best overall results in terms of energy and RMSD (value lower than 2A for successful docking results). This new multi-objective approach shows an improvement over the ligand-protein docking predictions that could be promising in in silico docking studies to select new anticancer compounds for therapeutic targets that are multidrug resistant.