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.

Domain and nucleotide dependence of the interaction between Saccharomyces cerevisiae translation elongation factors 3 and 1A

Eukaryotic translation elongation factor 3 (eEF3) is a fungal-specific ATPase proposed to catalyze the release of deacylated-tRNA from the ribosomal E-site. In addition, it has been shown to interact with the aminoacyl-tRNA binding GTPase elongation factor 1A (eEF1A), perhaps linking the E and A sites. Domain mapping demonstrates that amino acids 775-980 contain the eEF1A binding sites. Domain III of eEF1A, which is also involved in actin-related functions, is the site of eEF3 binding. The binding of eEF3 to eEF1A is enhanced by ADP, indicating the interaction is favored post-ATP hydrolysis but is not dependent on the eEF1A-bound nucleotide. A temperature-sensitive P915L mutant in the eEF1A binding site of eEF3 has reduced ATPase activity and affinity for eEF1A. These results support the model that upon ATP hydrolysis, eEF3 interacts with eEF1A to help catalyze the delivery of aminoacyl-tRNA at the A-site of the ribosome. The dynamics of when eEF3 interacts with eEF1A may be part of the signal for transition of the post to pre-translocational ribosomal state in yeast.

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.

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.

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.

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.

Regulating the expression of eEF1A to decipher its non canonical functions in eukaryotic cells

The canonical function of eEF1A is delivery of the aminoacylated tRNA to the A site of the ribosome during protein translation, however, it is also known to be an actin binding protein. As well as this actin binding function, eEF1A has been shown to be involved in other cellular processes such as cell proliferation and apoptosis. It has long been thought that the actin cytoskeleton and protein synthesis are linked and eEF1A has been suggested to be a candidate protein to form this link, though very little is understood about the relationship between its two functions. Overexpression of eEF1A has also been shown to be implicated in many different types of cancers, especially cancers that are metastatic, therefore it is important to further understand how eEF1A can affect both translation and the organisation of the actin cytoskeleton. To this end, we aimed to determine the effects of reduced expression of eEF1A on both translation and its non canonical functions in CHO cells. We have shown that reduced expression of eEF1A in this cell system results in no change in protein synthesis, however results in an increased number of actin stress fibres and other proteins associated with these fibres such as myosin IIA, paxillin and vinculin. Cell motility and attachment are also affected by this reduction in eEF1A protein expression. The organisational and motility phenotypes were found to be specific to eEF1A by transforming the cells with plasmids containing either human eEF1A1 or eEF1A2. Though the mechanisms by which these effects are regulated have not yet been established, this data provides evidence to show that the translation and actin binding functions of eEF1A are independent of each other as well as being suggestive of a role for eEF1A in cell motility as supported by the observation that overexpression of eEF1A protein tends to be associated with the cancer cells that are metastatic.

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.

Multiparameter flow cytometric analysis of C -reactive protein binding to human-derived cell lines and peripheral blood monocytes

C-reactive protein (CRP), a normally occurring human plasma protein may become elevated as much as 1,000 fold during disease states involving acute inflammation or tissue damage. Through its binding to phosphorylcholine in the presence of calcium, CRP has been shown to potentiate the activation of complement, stimulate phagocytosis and opsonize certain microorganisms. Utilizing a flow cytometric functional ligand binding assay I have demonstrated that a monocyte population in human peripheral blood and specific human-derived myelomonocytic cell lines reproducibly bind an evolutionarily conserved conformational pentraxin epitope on human CRP through a mechanism that does not involve its ligand, phosphorylcholine. ^ A variety of cell lines at different stages of differentiation were examined. The monocytic cell line, THP-1, bound the most CRP followed by U937 and KG-1a cells. The HL-60 cell line was induced towards either the granulocyte or monocyte pathway with DMSO or PMA, respectively. Untreated HL-60 cells or DMSO-treated cells did not bind CRP while cells treated with PMA showed increased binding of CRP, similar to U-937 cells. T cell and B-cell derived lines were negative. ^ Inhibition studies with Limulin and human SAP demonstrated that the binding site is a conserved pentraxin epitope. The calcium requirement necessary for binding to occur indicated that the cells recognize a conformational form of CRP. Phosphorylcholine did not inhibit the reaction therefore the possibility that CRP had bound to damaged membranes with exposed PC sites was discounted. ^ A study of 81 normal donors using flow cytometry demonstrated that a majority of peripheral blood monocytes (67.9 ± 1.3, mean ± sem) bound CRP. The percentage of binding was normally distributed and not affected by gender, age or ethnicity. Whole blood obtained from donors representing a variety of disease states showed a significant reduction in the level of CRP bound by monocytes in those donors classified with infection, inflammation or cancer. This reduction in monocyte populations binding CRP did not correlate with the concentration of plasma CRP. ^ The ability of monocytes to specifically bind CRP combined with the binding reactivity of the protein itself to a variety of phosphorylcholine containing substances may represent an important bridge between innate and adaptive immunity. ^

Conformational dynamics associated with ligand binding to vertebrate hexa-coordinate hemoglobins

Neuroglobin (Ngb) and cytoglobin (Cygb) are two new additions to the globin family, exhibiting heme iron hexa-coordination, a disulfide bond and large internal cavities. These proteins are implicated in cytoprotection under hypoxic-ischemic conditions, but the molecular basis of their cytoprotective function is unclear. Herein, a photothermal and spectroscopic study of the interactions of diatomic ligands with Ngb, Cygb, myoglobin and hemoglobin is presented. The impact of the disulfide bond in Ngb and Cygb and role of conserved residues in Ngb His64, Val68, Cys55, Cys120 and Tyr44 on conformational dynamics associated with ligand binding/dissociation were investigated. Transient absorption and photoacoustic calorimetry studies indicate that CO photo-dissociation from Ngb leads to a volume expansion (13.4±0.9 mL mol-1), whereas a smaller volume change was determined for Ngb with reduced Cys (ΔV=4.6±0.3 mL mol-1). Furthermore, Val68 side chain regulates ligand migration between the distal pocket and internal hydrophobic cavities since Val68Phe geminate quantum yield is ∼2.7 times larger than that of WT Ngb. His64Gln and Tyr44Phe mutations alter the thermodynamic parameters associated with CO photo-release indicating that electrostatic/hydrogen binding network that includes heme propionate groups, Lys 67, His64, and Tyr 44 in Ngb modulates the energetics of CO photo-dissociation. In Cygb, CO escape from the protein matrix is fast (< 40 ns) with a ΔH of 18±2 kcal mol-1 in Cygbred, whereas disulfide bridge formation promotes a biphasic ligand escape associated with an overall enthalpy change of 9±4 kcal mol-1. Therefore, the disulfide bond modulates conformational dynamics in Ngb and Cygb. I propose that in Cygb with reduced Cys the photo-dissociated ligand escapes through the hydrophobic tunnel as occurs in Ngb, whereas the CO preferentially migrates through the His64 gate in Cygbox. To characterize Cygb surface 1,8-ANS interactions with Cygb were investigated employing fluorescence spectroscopy, ITC and docking simulations. Two 1,8-ANS binding sites were identified. One binding site is located close to the extended N-terminus of Cygb and was also identified as a binding site for oleate. Furthermore, guanidinium hydrochloride-induced unfolding studies of Cygb reveal that the disulfide bond does not impact Cygb stability, whereas binding of cyanide slightly increases the protein stability.

Kinetic Analysis of Gill (Na+,K+)-ATPase Activity in Selected Ontogenetic Stages of the Amazon River Shrimp, (Decapoda, Palaemonidae): Interactions at ATP- and Cation-Binding Sites

We investigated modulation by ATP, Mg2+, Na+, K+ and NH4 (+) and inhibition by ouabain of (Na+,K+)-ATPase activity in microsomal homogenates of whole zoeae I and decapodid III (formerly zoea IX) and whole-body and gill homogenates of juvenile and adult Amazon River shrimps, . (Na+,K+)-ATPase-specific activity was increased twofold in decapodid III compared to zoea I, juveniles and adults, suggesting an important role in this ontogenetic stage. The apparent affinity for ATP ( (M) = 0.09 +/- A 0.01 mmol L-1) of the decapodid III (Na+,K+)-ATPase, about twofold greater than the other stages, further highlights this relevance. Modulation of (Na+,K+)-ATPase activity by K+ also revealed a threefold greater affinity for K+ ( (0.5) = 0.91 +/- A 0.04 mmol L-1) in decapodid III than in other stages; NH4 (+) had no modulatory effect. The affinity for Na+ ( (0.5) = 13.2 +/- A 0.6 mmol L-1) of zoea I (Na+,K+)-ATPase was fourfold less than other stages. Modulation by Na+, Mg2+ and NH4 (+) obeyed cooperative kinetics, while K+ modulation exhibited Michaelis-Menten behavior. Rates of maximal Mg2+ stimulation of ouabain-insensitive ATPase activity differed in each ontogenetic stage, suggesting that Mg2+-stimulated ATPases other than (Na+,K+)-ATPase are present. Ouabain inhibition suggests that, among the various ATPase activities present in the different stages, Na+-ATPase may be involved in the ontogeny of osmoregulation in larval The NH4 (+)-stimulated, ouabain-insensitive ATPase activity seen in zoea I and decapodid III may reflect a stage-specific means of ammonia excretion since functional gills are absent in the early larval stages.

Transport in the blood of an anti-tumor water-soluble rutheniun cyclopentadienyl complex: a fluorescence study on albumin binding

Dissertação de mestrado, Qualidade em Análises, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2014

The actinin family of actin crosslinking proteins: natural functions and potential applications in synthetic biology

Actinin and spectrin proteins are members of the Spectrin Family of Actin Crosslinking Proteins. The importance of these proteins in the cytoskeleton is demonstrated by the fact that they are common targets for disease causing mutations. In their most prominent roles, actinin and spectrin are responsible for stabilising and maintaining the muscle architecture during contraction, and providing shape and elasticity to the red blood cell in circulation, respectively. To carry out such roles, actinin and spectrin must possess important mechanical and physical properties. These attributes are desirable when choosing a building block for protein-based nanoconstruction. In this study, I assess the contribution of several disease-associated mutations in the actinin-1 actin binding domain that have recently been linked to a rare platelet disorder, congenital macrothrombocytopenia. I investigate the suitability of both actinin and spectrin proteins as potential building blocks for nanoscale structures, and I evaluate a fusion-based assembly strategy to bring about self-assembly of protein nanostructures. I report that the actinin-1 mutant proteins display increased actin binding compared to WT actinin-1 proteins. I find that both actinin and spectrin proteins exhibit enormous potential as nano-building blocks in terms of their stability and ability to self-assemble, and I successfully design and create homodimeric and heterodimeric bivalent building blocks using the fusion-based assembly strategy. Overall, this study has gathered helpful information that will contribute to furthering the advancement of actinin and spectrin knowledge in terms of their natural functions, and potential unnatural functions in protein nanotechnology.

Kinase Inhibitor Profile For Human Nek1, Nek6, And Nek7 And Analysis Of The Structural Basis For Inhibitor Specificity.

Human Neks are a conserved protein kinase family related to cell cycle progression and cell division and are considered potential drug targets for the treatment of cancer and other pathologies. We screened the activation loop mutant kinases hNek1 and hNek2, wild-type hNek7, and five hNek6 variants in different activation/phosphorylation statesand compared them against 85 compounds using thermal shift denaturation. We identified three compounds with significant Tm shifts: JNK Inhibitor II for hNek1(Δ262-1258)-(T162A), Isogranulatimide for hNek6(S206A), andGSK-3 Inhibitor XIII for hNek7wt. Each one of these compounds was also validated by reducing the kinases activity by at least 25%. The binding sites for these compounds were identified by in silico docking at the ATP-binding site of the respective hNeks. Potential inhibitors were first screened by thermal shift assays, had their efficiency tested by a kinase assay, and were finally analyzed by molecular docking. Our findings corroborate the idea of ATP-competitive inhibition for hNek1 and hNek6 and suggest a novel non-competitive inhibition for hNek7 in regard to GSK-3 Inhibitor XIII. Our results demonstrate that our approach is useful for finding promising general and specific hNekscandidate inhibitors, which may also function as scaffolds to design more potent and selective inhibitors.