Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using Only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.

Binding of an RNA trafficking response element to heterogeneous nuclear ribonucleoproteins A1 and A2

Heterogeneous nuclear ribonucleoprotein (hnRNP) A2 binds a 21-nucleotide myelin basic protein mRNA response element, the A2RE, and A2RE-like sequences in other localized mRNAs, and is a trans-acting factor in oligodendrocyte cytoplasmic RNA trafficking. Recombinant human hnRNPs A1 and A2 were used in a biosensor to explore interactions with A2RE and the cognate oligodeoxyribonucleotide. Both proteins have a single site that bound oligonucleotides with markedly different sequences but did not bind in the presence of heparin. Both also possess a second, specific site that bound only A2RE and was unaffected by heparin, hnRNP A2 bound A2RE in the latter site with a K-d near 50 nM, whereas the K-d for hnRNP A1 was above 10 muM. UV cross-linking assays led to a similar conclusion. Mutant A2RE sequences, that in earlier qualitative studies appeared not to bind hnRNP A2 or support RNA trafficking in oligodendrocytes, had dissociation constants above 5 muM for this protein. The two concatenated RNA recognition motifs (RRMs), but not the individual RRMs, mimicked the binding behavior of hnRNP A2. These data highlight the specificity of the interaction of A2RE with these hnRNPs and suggest that the sequence-specific A2RE-binding site on hnRNP A2 is formed by both RRMs acting in cis.

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.

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.

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

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.

Digestive Peptidase Evolution In Holometabolous Insects Led To A Divergent Group Of Enzymes In Lepidoptera.

Trypsins and chymotrypsins are well-studied serine peptidases that cleave peptide bonds at the carboxyl side of basic and hydrophobic l-amino acids, respectively. These enzymes are largely responsible for the digestion of proteins. Three primary processes regulate the activity of these peptidases: secretion, precursor (zymogen) activation and substrate-binding site recognition. Here, we present a detailed phylogenetic analysis of trypsins and chymotrypsins in three orders of holometabolous insects and reveal divergent characteristics of Lepidoptera enzymes in comparison with those of Coleoptera and Diptera. In particular, trypsin subsite S1 was more hydrophilic in Lepidoptera than in Coleoptera and Diptera, whereas subsites S2-S4 were more hydrophobic, suggesting different substrate preferences. Furthermore, Lepidoptera displayed a lineage-specific trypsin group belonging only to the Noctuidae family. Evidence for facilitated trypsin auto-activation events were also observed in all the insect orders studied, with the characteristic zymogen activation motif complementary to the trypsin active site. In contrast, insect chymotrypsins did not seem to have a peculiar evolutionary history with respect to their mammal counterparts. Overall, our findings suggest that the need for fast digestion allowed holometabolous insects to evolve divergent groups of peptidases with high auto-activation rates, and highlight that the evolution of trypsins led to a most diverse group of enzymes in Lepidoptera.