Characterization of heparin-binding site of tissue transglutaminase:its importance in cell surface targeting, matrix deposition, and cell signaling

Tissue transglutaminase (TG2) is a multifunctional Ca2+ activated protein crosslinking enzyme secreted into the extracellular matrix (ECM), where it is involved in wound healing and scarring, tissue fibrosis, celiac disease and metastatic cancer. Extracellular TG2 can also facilitate cell adhesion important in wound healing through a non-transamidating mechanism via its association with fibronectin (FN), heparan sulphates (HS) and integrins. Regulating the mechanism how TG2 is translocated into the ECM therefore provides a strategy for modulating these physiological and pathological functions of the enzyme. Here, through molecular modelling and mutagenesis we have identified the HS binding site of TG2 202KFLKNAGRDCSRRSSPVYVGR222. We demonstrate the requirement of this binding site for translocation of TG2 into the ECM through a mechanism involving cell surface shedding of HS. By synthesizing a peptide NPKFLKNAGRDCSRRSS corresponding to the HS binding site within TG2, we also demonstrate how this mimicking peptide can in isolation compensate the RGD-induced loss of cell adhesion on FN via binding to syndecan-4, leading to activation of PKCa, pFAK-397 and ERK1/2 and the subsequent formation of focal adhesions and actin cytoskeleton organization. A novel regulatory mechanism for TG2 translocation into the extracellular compartment that depends upon TG2 conformation and the binding of HS is proposed.

PDNAsite:identification of DNA-binding site from protein sequence by incorporating spatial and sequence context

Protein-DNA interactions are involved in many fundamental biological processes essential for cellular function. Most of the existing computational approaches employed only the sequence context of the target residue for its prediction. In the present study, for each target residue, we applied both the spatial context and the sequence context to construct the feature space. Subsequently, Latent Semantic Analysis (LSA) was applied to remove the redundancies in the feature space. Finally, a predictor (PDNAsite) was developed through the integration of the support vector machines (SVM) classifier and ensemble learning. Results on the PDNA-62 and the PDNA-224 datasets demonstrate that features extracted from spatial context provide more information than those from sequence context and the combination of them gives more performance gain. An analysis of the number of binding sites in the spatial context of the target site indicates that the interactions between binding sites next to each other are important for protein-DNA recognition and their binding ability. The comparison between our proposed PDNAsite method and the existing methods indicate that PDNAsite outperforms most of the existing methods and is a useful tool for DNA-binding site identification. A web-server of our predictor (http://hlt.hitsz.edu.cn:8080/PDNAsite/) is made available for free public accessible to the biological research community.

Interaction of peptides, peptidomimetics and other drug candidates with the DI/tripeptide transport system in intestinal epithelial cells, using the in vitro caco-2 cell culture system

An uptake system was developed using Caco-2 cell monolayers and the dipeptide, glycyl-[3H]L-proline, as a probe compound. Glycyl-[3H]L-proline uptake was via the di-/tripeptide transport system (DTS) and, exhibited concentration-, pH- and temperature-dependency. Dipeptides inhibited uptake of the probe, and the design of the system allowed competitors to be ranked against one another with respect to affinity for the transporter. The structural features required to ensure or increase interaction with the DTS were defined by studying the effect of a series of glycyl-L-proline and angiotensin-converting enzyme (ACE)-inhibitor (SQ-29852) analogues on the uptake of the probe. The SQ-29852 structure was divided into six domains (A-F) and competitors were grouped into series depending on structural variations within specific regions. Domain A was found to prefer a hydrophobic function, such as a phenyl group, and was intolerant to positive charges and H+ -acceptors and donors. SQ-29852 analogues were more tolerant of substitutions in the C domain, compared to glycyl-L-proline analogues, suggesting that interactions along the length of the SQ-29852 molecule may override the effects of substitutions in the C domain. SQ-29852 analogues showed a preference for a positive function, such as an amine group in this region, but dipeptide structures favoured an uncharged substitution. Lipophilic substituents in domain D increased affinity of SQ-29852 analogues with the DTS. A similar effect was observed for ACE-NEP inhibitor analogues. Domain E, corresponding to the carboxyl group was found to be tolerant of esterification for SQ-29852 analogues but not for dipeptides. Structural features which may increase interaction for one series of compounds, may not have the same effect for another series, indicating that the presence of multiple recognition sites on a molecule may override the deleterious effect of anyone change. Modifying current, poorly absorbed peptidomimetic structures to fit the proposed hypothetical model may improve oral bioavailability by increasing affinity for the DTS. The stereochemical preference of the transporter was explored using four series of compounds (SQ-29852, lysylproline, alanylproline and alanylalanine enantiomers). The L, L stereochemistry was the preferred conformation for all four series, agreeing with previous studies. However, D, D enantiomers were shown in some cases to be substrates for the DTS, although exhibiting a lower affinity than their L, L counterparts. All the ACE-inhibitors and β-lactam antibiotics investigated, produced a degree of inhibition of the probe, and thus show some affinity for the DTS. This contrasts with previous reports that found several ACE inhibitors to be absorbed via a passive process, thus suggesting that compounds are capable of binding to the transporter site and inhibiting the probe without being translocated into the cell. This was also shown to be the case for oligodeoxynucleotide conjugated to a lipophilic group (vitamin E), and highlights the possibility that other orally administered drug candidates may exert non-specific effects on the DTS and possibly have a nutritional impact. Molecular modelling of selected ACE-NEP inhibitors revealed that the three carbonyl functions can be oriented in a similar direction, and this conformation was found to exist in a local energy-minimised state, indicating that the carbonyls may possibly be involved in hydrogen-bond formation with the binding site of the DTS.

Behavioural and neurochemical correlates of drugs acting at imidazoline and a2-adrenoceptor binding sites

A number of agents with differing selectivity profiles for the non-a2 adrenoceptor binding site (NAIBS), imidazoline preferring receptor (IPR) and a2-adrenoceptor were employed in a series of behavioural and neurochemical experiments to determine a functional role for the former two sites. The highly selective NAIBS ligand RX801 077 produced an increase in rat brain extracellular noradrenaline (NA) levels, as determined by the technique of in vivo microdialysis, which may underlie its ability to produce a discriminable cue in the same species. This increase in NA may be due to a suggested link between the NAIBS and the monoamine oxidase inhibitor (MAOI) activity of RX801 077. For instance, the RX801 077 cue was substituted for by the MAOI drugs pargyline and moclobemide, which themselves down regulate NAIBS when administered chronically. RX811 059 substituted for the RX801 077 cue which may be due its ability to stimulate NA release via its activity as a highly selective a2-adrenoceptor antagonist. An effect upon NA output may also explain the ability of RX801 077 to 'mimic' the anti-immobility effect of the antidepressant drug desmethylimipramine (DMJ) in the forced swimming test. Further studies are therefore required to examine a possible role for the NAIBS in the treatment of depression. Discriminable cues were also produced by RX811 059 and the a2- adrenoceptor agonist clonidine, probably as a consequence of their respective ability to stimulate and inhibit NA output via their opposing activity at a2-adrenoceptors. The IPR has been suggested to play a role in mediating the hypotensive effect of clonidine, although a precise role was unable to be established for this site in the present studies due to the unavailability of highly selective IPA agents.

Characterization of metal ion-induced [3H]inositol hexakisphosphate binding to rat cerebellar membranes

The binding of [3H]inositol hexakisphosphate ([3H] InsP6) to rat cerebellar membranes has been characterized with the objective of establishing the role, if any, of a membrane protein receptor. In the presence of EDTA, we have previously identified an InsP6-binding site with a capacity of approximately 20 pmol/mg protein (Hawkins, P. T., Reynolds, D. J. M., Poyner, D. R., and Hanley, M. R. (1990) Biochem. Biophys. Res. Commun. 167, 819-827). However, in the presence of 1 mM Mg2+, the capacity of [3H]InsP6 binding to membranes was increased approximately 9-fold. This enhancing effect of Mg2+ was reversed by addition of 10 microM of several cation chelators, suggesting that the increased binding required trace quantities of other metal cations. This is supported by experiments where it was possible to saturate binding by addition of excess membranes, despite not significantly depleting radioligand, pointing to removal of some other factor. Removal of endogenous cations from the binding assay by pretreatment with chelex resin also prevents the Mg(2+)-induced potentiation. Consideration of the specificity of the chelators able to abolish this potentiation suggested involvement of Fe3+ or Al3+. Both these ions (but not several others) were able to increase [3H]InsP6 binding to chelex-pretreated membranes at concentrations of 1 microM. It is possible to demonstrate synergy between Fe3+ and Mg2+ under these conditions. We propose that [3H]InsP6 may interact with membranes through non-protein recognition possibly via phospholipids, in a manner dependent upon trace metals. The implications of this for InsP6 biology are considered.

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

In vivo and in vitro aspects of imidazoline-2(1/2)sites within the central nervous system

The in vivo and in vitro characteristics of the I2 binding site were probed using the technique of drug discrimination and receptor autoradiography. Data presented in this thesis indicates the I2 ligand 2-BFI generates a cue in drug discrimination. Further studies indicated agmatine, a proposed endogenous imidazoline ligand, and a number of imidazoline and imidazole analogues of 2-BFI substitute significantly for 2-BFI. In addition to specific I2 ligands the administration of NRl's (noradrenaline reuptake inhibitors), the sympathomimetic d-amphetamine, the α1-adrenoceptor agonist methoxamine, but not the β1 agonist dobutamine or the β2 agonist salbutamol, gave rise to significant levels of substitution for the 2-BFI cue. The administration of the α1-adrenoceptor antagonist WB4101, prior to 2- BFI itself significantly reduced levels of 2-BFI appropriate responding. Administration of the reversible MAO-A inhibitors moclobemide and Ro41-1049, but not the reversible MAO-B inhibitors lazabemide and Ro16-6491, gave rise to potent dose dependent levels of substitution for the 2-BFI cue. Further studies indicated the administration of a number of β-carbolines and the structurally related indole alkaloid ibogaine also gave rise to dose dependent significant levels of substitution. Due to the relationship of indole alkaloids to serotonin the 5-HT releaser fenfluramine and a number of SSRI's (selective serotonin reuptake inhibitor) were also administered and these compounds gave rise to significant partial (20-80% responses to the 2-BFI lever) levels of substitution. The autoradiographical studies reported here indicate [3H]2-BFI labels I2 sites within the rat arcuate nucleus, area postrema, pineal gland, interpeduncular nucleus and subfornical organ. Subsequent experiments confirmed that the drug discrimination dosing schedule significantly increases levels of [3H]2-BFI 12 binding within two of these nuclei. However, levels of [3H]2-BFI specific binding were significantly reduced within four of these nuclei after chronic treatment with the irreversible MAO inhibitors deprenyl and tranylcypromine but not pargyline, which only reduced levels significantly in two. Further autoradiographical studies indicated that the distribution of [3H]2-BFI within the C57/B mouse compares favourably to that within the rat. Comparison of these levels of binding to those from transgenic mice who over-express MAO-B indicates two possibly distinct populations of [3H]2-BFI 12 sites exist in mouse brain. The data presented here indicates the 2-BFI cue is associated with the selective activation of α1-adrenoceptors and possibly 5-HT receptors. 2-BFI trained rats recognise reversible MAO-A but not MAO-B inhibitors. However, data within this thesis indicates the autoradiographical distribution of I2 sites bears a closer resemblance to that of MAO-B not MAO-A and further studies using transgenic mice that over-express MAO-B suggests a non-MAO-B I2 site exists in mouse brain.

Design and synthesis of potential antimicrobial agents

Tuberculosis is one of the most devastating diseases in the world primarily due to several decades of neglect and an emergence of multidrug-resitance strains (MDR) of M. tuberculosis together with the increased incidence of disseminated infections produced by other mycobacterium in AIDS patients. This has prompted the search for new antimycobacterial drugs. A series of pyridine-2-, pyridine-3-, pyridine-4-, pyrazine and quinoline-2-carboxamidrazone derivatives and new classes of carboxamidrazone were prepared in an automated fashion and by traditional synthesis. Over nine hundred synthesized compounds were screened for their anti mycobacterial activity against M. fortutium (NGTG 10394) as a surrogate for M. tuberculosis. The new classes of amidrazones were also screened against tuberculosis H37 Rv and antimicrobial activities against various bacteria. Fifteen tested compounds were found to provide 90-100% inhibition of mycobacterium growth of M. tuberculosis H37 Rv in the primary screen at 6.25 μg mL-1. The most active compound in the carboxamidrazone amide series had an MIG value of 0.1-2 μg mL-1 against M. fortutium. The enzyme dihydrofolate reductase (DHFR) has been a drug-design target for decades. Blocking of the enzymatic activity of DHFR is a key element in the treatment of many diseases, including cancer, bacterial and protozoal infection. The x-ray structure of DHFR from M. tuberculosis and human DHFR were found to have differences in substrate binding site. The presence of glycerol molecule in the Xray structure from M. tuberculosis DHFR provided opportunity to design new antifolates. The new antifolates described herein were designed to retain the pharmcophore of pyrimethamine (2,4- diamino-5(4-chlorophenyl)-6-ethylpyrimidine), but encompassing a range of polar groups that might interact with the M. tuberculosis DHFR glycerol binding pockets. Finally, the research described in this thesis contributes to the preparation of molecularly imprinted polymers for the recognition of 2,4-diaminopyrimidine for the binding the target. The formation of hydrogen bonding between the model functional monomer 5-(4-tert-butyl-benzylidene)-pyrimidine-2,4,6-trione and 2,4-diaminopyrimidine in the pre-polymerisation stage was verified by 1H-NMR studies. Having proven that 2,4-diaminopyrimidine interacts strongly with the model 5-(4-tert-butylbenzylidene)- pyrimidine-2,4,6-trione, 2,4-diaminopyrimidine-imprinted polymers were prepared using a novel cyclobarbital derived functional monomer, acrylic acid 4-(2,4,6-trioxo-tetrahydro-pyrimidin-5- ylidenemethyl)phenyl ester, capable of multiple hydrogen bond formation with the 2,4- diaminopyrimidine. The recognition property of the respective polymers toward the template and other test compounds was evaluated by fluorescence. The results demonstrate that the polymers showed dose dependent enhancement of fluorescence emissions. In addition, the results also indicate that synthesized MIPs have higher 2,4-diaminopyrimidine binding ability as compared with corresponding non-imprinting polymers.

Extracellular loops 1 and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function

The first and third extracellular loops (ECL) of G protein-coupled receptors (GPCRs) have been implicated in ligand binding and receptor function. This study describes the results of an alanine/leucine scan of ECLs 1 and 3 and loop-associated transmembrane (TM) domains of the secretin-like GPCR calcitonin receptor-like receptor which associates with receptor activity modifying protein 1 to form the CGRP receptor. Leu195Ala, Val198Ala and Ala199Leu at the top of TM2 all reduced aCGRP-mediated cAMP production and internalization; Leu195Ala and Ala199Leu also reduced aCGRP binding. These residues form a hydrophobic cluster within an area defined as the "minor groove" of rhodopsin-like GPCRs. Within ECL1, Ala203Leu and Ala206Leu influenced the ability of aCGRP to stimulate adenylate cyclase. In TM3, His219Ala, Leu220Ala and Leu222Ala have influences on aCGRP binding and cAMP production; they are likely to indirectly influence the binding site for aCGRP as well as having an involvement in signal transduction. On the exofacial surfaces of TMs 6 and 7, a number of residues were identified that reduced cell surface receptor expression, most noticeably Leu351Ala and Glu357Ala in TM6. The residues may contribute to the RAMP1 binding interface. Ile360Ala impaired aCGRP-mediated cAMP production. Ile360 is predicted to be located close to ECL2 and may facilitate receptor activation. Identification of several crucial functional loci gives further insight into the activation mechanism of this complex receptor system and may aid rational drug design.

Role of GSH in estrone sulfate binding and translocation by the multidrug resistance protein 1 (MRP1/ABCC1)

Multidrug resistance protein 1 (MRP1/ABCC1) is an ATP-dependent efflux pump that can confer resistance to multiple anticancer drugs and transport conjugated organic anions. Unusually, transport of several MRP1 substrates requires glutathione (GSH). For example, estrone sulfate transport by MRP1 is stimulated by GSH, vincristine is co-transported with GSH, or GSH can be transported alone. In the present study, radioligand binding assays were developed to investigate the mechanistic details of GSH-stimulated transport of estrone sulfate by MRP1. We have established that estrone sulfate binding to MRP1 requires GSH, or its non-reducing analogue S-methyl GSH (S-mGSH), and further that the affinity (Kd) of MRP1 for estrone sulfate is 2.5-fold higher in the presence of S-mGSH than GSH itself. Association kinetics show that GSH binds to MRP1 first, and we propose that GSH binding induces a conformational change, which makes the estrone sulfate binding site accessible. Binding of non-hydrolyzable ATP analogues to MRP1 decreases the affinity for estrone sulfate. However, GSH (or S-mGSH) is still required for estrone sulfate binding, and the affinity for GSH is unchanged. Estrone sulfate affinity remains low following hydrolysis of ATP. The affinity for GSH also appears to decrease in the post-hydrolytic state. Our results indicate ATP binding is sufficient for reconfiguration of the estrone sulfate binding site to lower affinity and argue for the presence of a modulatory GSH binding site not associated with transport of this tripeptide. A model for the mechanism of GSH-stimulated estrone sulfate transport is proposed.

The coupling mechanism of P-glycoprotein involves residue L339 in the sixth membrane spanning segment

The transmembrane (TM) domains in P-glycoprotein (P-gp) contain the drug binding sites and undergo conformational changes driven by nucleotide catalysis to effect translocation. However, our understanding of exactly which regions are involved in such events remains unclear. A site-directed labelling approach was used to attach thiol-reactive probes to cysteines introduced into transmembrane segment 6 (TM6) in order to perturb function and infer involvement of specific residues in drug binding and/or interdomain communication. Covalent attachment of coumarin-maleimide at residue 339C within TM6 resulted in impaired ATP hydrolysis by P-gp. The nature of the effect was to reduce the characteristic modulation of basal activity caused by transported substrates, modulators and the potent inhibitor XR9576. Photoaffinity labelling of P-gp with [(3)H]-azidopine indicated that residue 339C does not alter drug binding per se. However, covalent modification of this residue appears to prevent conformational changes that lead to drug stimulation of ATP hydrolysis.

Toward prediction of class II mouse major histocompatibility complex peptide binding affinity:in silico bioinformatic evaluation using partial least squares, a robust multivariate statistical technique

The accurate identification of T-cell epitopes remains a principal goal of bioinformatics within immunology. As the immunogenicity of peptide epitopes is dependent on their binding to major histocompatibility complex (MHC) molecules, the prediction of binding affinity is a prerequisite to the reliable prediction of epitopes. The iterative self-consistent (ISC) partial-least-squares (PLS)-based additive method is a recently developed bioinformatic approach for predicting class II peptide−MHC binding affinity. The ISC−PLS method overcomes many of the conceptual difficulties inherent in the prediction of class II peptide−MHC affinity, such as the binding of a mixed population of peptide lengths due to the open-ended class II binding site. The method has applications in both the accurate prediction of class II epitopes and the manipulation of affinity for heteroclitic and competitor peptides. The method is applied here to six class II mouse alleles (I-Ab, I-Ad, I-Ak, I-As, I-Ed, and I-Ek) and included peptides up to 25 amino acids in length. A series of regression equations highlighting the quantitative contributions of individual amino acids at each peptide position was established. The initial model for each allele exhibited only moderate predictivity. Once the set of selected peptide subsequences had converged, the final models exhibited a satisfactory predictive power. Convergence was reached between the 4th and 17th iterations, and the leave-one-out cross-validation statistical terms - q2, SEP, and NC - ranged between 0.732 and 0.925, 0.418 and 0.816, and 1 and 6, respectively. The non-cross-validated statistical terms r2 and SEE ranged between 0.98 and 0.995 and 0.089 and 0.180, respectively. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made freely available online (http://www.jenner.ac.uk/MHCPred).

Cortical oscillatory dynamics and benzodiazepine-site modulation of tonic inhibition in fast spiking interneurons

Tonic conductance mediated by extrasynaptic GABAA receptors has been implicated in the modulation of network oscillatory activity. Using an in vitro brain slice to produce oscillatory activity and a kinetic model of GABAA receptor dynamics, we show that changes in tonic inhibitory input to fast spiking interneurons underlie benzodiazepine-site mediated modulation of neuronal network synchrony in rat primary motor cortex. We found that low concentrations (10 nM) of the benzodiazepine site agonist, zolpidem, reduced the power of pharmacologically-induced beta-frequency (15–30 Hz) oscillatory activity. By contrast, higher doses augmented beta power. Application of the antagonist, flumazenil, also increased beta power suggesting endogenous modulation of the benzodiazepine binding site. Voltage-clamp experiments revealed that pharmacologically-induced rhythmic inhibitory postsynaptic currents were reduced by 10 nM zolpidem, suggesting an action on inhibitory interneurons. Further voltage -clamp studies of fast spiking cells showed that 10 nM zolpidem augmented a tonic inhibitory GABAA receptor mediated current in fast spiking cells whilst higher concentrations of zolpidem reduced the tonic current. A kinetic model of zolpidem-sensitive GABAA receptors suggested that incubation with 10 nM zolpidem resulted in a high proportion of GABAA receptors locked in a kinetically slow desensitized state whilst 30 nM zolpidem favoured rapid transition into and out of desensitized states. This was confirmed experimentally using a challenge with saturating concentrations of GABA. Selective modulation of an interneuron-specific tonic current may underlie the reversal of cognitive and motor deficits afforded by low-dose zolpidem in neuropathological states.

Structural basis for receptor activity-modifying protein-dependent selective peptide recognition by a G protein-coupled receptor

Association of receptor activity-modifying proteins (RAMP1-3) with the G protein-coupled receptor (GPCR) calcitonin receptor-like receptor (CLR) enables selective recognition of the peptides calcitonin gene-related peptide (CGRP) and adrenomedullin (AM) that have diverse functions in the cardiovascular and lymphatic systems. How peptides selectively bind GPCR:RAMP complexes is unknown. We report crystal structures of CGRP analog-bound CLR:RAMP1 and AM-bound CLR:RAMP2 extracellular domain heterodimers at 2.5 and 1.8 Å resolutions, respectively. The peptides similarly occupy a shared binding site on CLR with conformations characterized by a β-turn structure near their C termini rather than the α-helical structure common to peptides that bind related GPCRs. The RAMPs augment the binding site with distinct contacts to the variable C-terminal peptide residues and elicit subtly different CLR conformations. The structures and accompanying pharmacology data reveal how a class of accessory membrane proteins modulate ligand binding of a GPCR and may inform drug development targeting CLR:RAMP complexes.

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.