Discovering properties of new DNA-binding activity of proteins

Protein-DNA interactions are an essential feature in the genetic activities of life, and the ability to predict and manipulate such interactions has applications in a wide range of fields. This Thesis presents the methods of modelling the properties of protein-DNA interactions. In particular, it investigates the methods of visualising and predicting the specificity of DNA-binding Cys2His2 zinc finger interaction. The Cys2His2 zinc finger proteins interact via their individual fingers to base pair subsites on the target DNA. Four key residue positions on the a- helix of the zinc fingers make non-covalent interactions with the DNA with sequence specificity. Mutating these key residues generates combinatorial possibilities that could potentially bind to any DNA segment of interest. Many attempts have been made to predict the binding interaction using structural and chemical information, but with only limited success. The most important contribution of the thesis is that the developed model allows for the binding properties of a given protein-DNA binding to be visualised in relation to other protein-DNA combinations without having to explicitly physically model the specific protein molecule and specific DNA sequence. To prove this, various databases were generated, including a synthetic database which includes all possible combinations of the DNA-binding Cys2His2 zinc finger interactions. NeuroScale, a topographic visualisation technique, is exploited to represent the geometric structures of the protein-DNA interactions by measuring dissimilarity between the data points. In order to verify the effect of visualisation on understanding the binding properties of the DNA-binding Cys2His2 zinc finger interaction, various prediction models are constructed by using both the high dimensional original data and the represented data in low dimensional feature space. Finally, novel data sets are studied through the selected visualisation models based on the experimental DNA-zinc finger protein database. The result of the NeuroScale projection shows that different dissimilarity representations give distinctive structural groupings, but clustering in biologically-interesting ways. This method can be used to forecast the physiochemical properties of the novel proteins which may be beneficial for therapeutic purposes involving genome targeting in general.

In silico modelling of the interaction of flavonoids with human P-glycoprotein nucleotide-binding domain

A three-dimensional model of human ABCB1 nucleotide-binding domain (NBD) was developed by homology modelling using the high-resolution human TAP1 transporter structure as template. Interactions between NBD and flavonoids were investigated using in silico docking studies. Ring-A of unmodified flavonoid was located within the NBD P-loop with the 5-hydroxyl group involved in hydrogen bonding with Lys1076. Ring-B was stabilised by hydrophobic stacking interactions with Tyr1044. The 3-hydroxyl group and carbonyl oxygen were extensively involved in hydrogen bonding interactions with amino acids within the NBD. Addition of prenyl, benzyl or geranyl moieties to ring-A (position-6) and hydrocarbon substituents (O-n-butyl to O-n-decyl) to ring-B (position-4) resulted in a size-dependent decrease in predicted docking energy which reflected the increased binding affinities reported in vitro.

Interaction of calcitonin-gene-related peptide with its receptors

The receptor for calcitonin-gene-related peptide (CGRP) is a heterodimer formed by calcitonin-receptor-like receptor (CRLR), a type II (family B) G-protein-coupled receptor, and receptor-activity-modifying protein 1 (RAMP1), a single-membrane-pass protein. It is likely that the first seven or so amino acids of CGRP (which form a disulphide-bonded loop) interact with the transmembrane domain of CRLR to cause receptor activation. The rest of the CGRP molecule falls into three domains. Residues 28-37 and 8-18 are normally required for high-affinity binding, while residues 19-27 form a hinge region. The 28-37 region is almost certainly in direct contact with the receptor; 8-18 may make additional receptor contacts or may stabilize an appropriate conformation of 28-37. It is likely that these regions of CGRP interact both with CRLR and with the extracellular domain of RAMP1.

The C1q binding activity of IgG is modified in vitro by reactive oxygen species: implications for rheumatoid arthritis

IgG can be denatured in vitro by reactive oxygen species (ROS). Native IgG activates the complement cascade through C1q. Using a modified ELISA, C1q binding activity of rheumatoid IgG has been compared to IgG denatured by neutrophil-derived ROS. The C1q binding activity of rheumatoid synovial fluid IgG is greater than the corresponding serum IgG (P < 0.01). Denaturation of IgG by activated polymorphs or the Fenton reaction decreased its C1q binding activity (P < 0.01). In vitro exposure of IgG to OH. and ROO. increased its interaction with C1q (P < 0.01). Hypochlorous acid had no effect. ROS-induced alteration to IgG-C1q binding activity may promote the inflammatory response in rheumatoid arthritis.

LacO-LacI interaction in affinity adsorption of plasmid DNA

Current approaches for purifying plasmids from bacterial production systems exploit the physiochemical properties of nucleic acids in non-specific capture systems. In this study, an affinity system for plasmid DNA (pDNA) purification has been developed utilizing the interaction between the lac operon (lacO) sequence contained in the pDNA and a 64mer synthetic peptide representing the DNA-binding domain of the lac repressor protein, LacI. Two plasmids were evaluated, the native pUC19 and pUC19 with dual lacO3/lacOs operators (pUC19lacO3/lacOs), where the lacOs operator is perfectly symmetrical. The DNA-protein affinity interaction was evaluated by surface plasmon resonance using a Biacore system. The affinity capture of DNA in a chromatography system was evaluated using LacI peptide that had been immobilized to Streamline™ adsorbent. The KD-values for double stranded DNA (dsDNA) fragments containing lacO1 and lacO3 and lacOs and lacO3 were 5.7 ± 0.3 × 10 -11 M and 4.1 ± 0.2 × 10-11 M respectively, which compare favorably with literature reports of 5 × 10-10 - 1 × 10-9 M for native laCO1 and 1-1.2 × 10-10 M for lacO1 in a saline buffer. Densitometric analysis of the gel bands from the affinity chromatography run clearly showed a significant preference for capture of the supercoiled fraction from the feed pDNA sample. The results indicate the feasibility of the affinity approach for pDNA capture and purification using native protein-DNA interaction. © 2006 Wiley Periodicals, Inc.

Thermodynamics of binding of regulatory ligands to tissue transglutaminase

The transamidating activity of tissue transglutaminase is regulated by the ligands calcium and GTP, via conformational changes which facilitate or interfere with interaction with the peptidyl-glutamine substrate. We have analysed binding of these ligands by calorimetric and computational approaches. In the case of GTP we have detected a single high affinity site (K (D) approximately 1 muM), with moderate thermal effects suggestive that binding GTP involves replacement of GDP, normally bound to the protein. On line with this possibility no significant binding was observed during titration with GDP and computational studies support this view. Titration with calcium at a high cation molar excess yielded a complex binding isotherm with a number of "apparent binding sites" in large excess over those detectable by equilibrium dialysis (6 sites). This binding pattern is ascribed to occurrence of additional thermal contributions, beyond those of binding, due to the occurrence of conformational changes and to catalysis itself (with protein self-crosslinking). In contrast only one site for binding calcium with high affinity (K (D) approximately 0.15 muM) is observed with samples of enzyme inactivated by alkylation at the active site (to prevent enzyme crosslinkage and thermal effects of catalysis). These results indicate an intrinsic ability of tissue transglutaminase to bind calcium with high affinity and the necessity of careful reassessment of the enzyme regulatory pattern in relation to the concentrations of ligands in living cells, taking also in account effects of ligands on protein subcellular compartimentation.

A functional chimeric DNA primase: the Cys4 zinc-binding domain of bacteriophage T3 primase fused to the helicase of bacteriophage T7

Two colinear bacteriophage T7 gene 4 proteins provide helicase and primase functions in vivo. T7 primase differs from T7 helicase by an additional 63 residues at the amino terminus. This terminal domain contains a zinc-binding motif which mediates an interaction with the basic primase recognition sequence 3'-CTG-5'. We have generated a chimeric primase in which the 81 amino-terminal residues are derived from the primase of phage T3 and the 484 carboxyl-terminal residues are those of phage T7 helicase. The amino-terminal domain of T3 primase is 50% homologous with that of T7 primase. The resulting T3/T7 chimeric protein is a functional primase in vivo. While the primase activity of the purified protein is about one-third that of T7 primase, the recognition sites used and the oligoribonucleotides synthesized from these sites are identical. We conclude that the residues responsible for the interaction with the sequence 3'-CTG-5' are conserved between the chimeric and T7 proteins.

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.

Acanthamoeba polyphaga trophozoite binding of representative fungal single cell forms

Acanthamoeba polyphaga trophozoites bind yeast cells of Candida albicans isolates within a few hours, leaving few cells in suspension or still attached to trophozoite surfaces. The nature of yeast cell recognition, mediated by an acanthamoebal trophozoite mannose binding protein is confirmed by experiments utilizing concentration dependent mannose hapten blocking. Similarly, acapsulate cells of Cryptococcus neoformans are also bound within a relatively short timescale. However, even after protracted incubation many capsulate cells of Cryptococcus remain in suspension, suggesting that the capsulate cell form of this species is not predated by acanthamoebal trophozoites. Further aspects of the association of Acanthamoeba and fungi are apparent when studying their interaction with conidia of the biocontrol agent Coniothyrium minitans. Conidia which readily bind with increasing maturity of up to 42 days, were little endocytosed and even released. Cell and conidial surface mannose as determined by FITC-lectin binding, flow cytometry with associated ligand binding analysis and hapten blocking studies demonstrates the following phenomena. Candida isolates and acapsulate Cryptococcus expose most mannose, while capsulate Cryptococcus cells exhibit least exposure commensurate with yeast cellular binding or lack of trophozoites. Conidia of Coniothyrium, albeit in a localized fashion, also manifest surface mannose exposure but as shown by Bmax values, in decreasing amounts with increasing maturity. Contrastingly such conidia experience greater trophozoite binding with maturation, thereby questioning the primacy of a trophozoite mannose-binding-protein recognition model.

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.

Static energy analysis of MHC class I and class II peptide-binding affinity

Antigenic peptide is presented to a T-cell receptor (TCR) through the formation of a stable complex with a major histocompatibility complex (MHC) molecule. Various predictive algorithms have been developed to estimate a peptide's capacity to form a stable complex with a given MHC class II allele, a technique integral to the strategy of vaccine design. These have previously incorporated such computational techniques as quantitative matrices and neural networks. A novel predictive technique is described, which uses molecular modeling of predetermined crystal structures to estimate the stability of an MHC class II-peptide complex. The structures are remodeled, energy minimized, and annealed before the energetic interaction is calculated.

Evidence that interaction between conserved residues in transmembrane helices 2, 3, and 7 are crucial for human VPAC1 receptor activation

The VPAC(1) receptor belongs to family B of G protein-coupled receptors (GPCR-B) and is activated upon binding of the vasoactive intestinal peptide (VIP). Despite the recent determination of the structure of the N terminus of several members of this receptor family, little is known about the structure of the transmembrane (TM) region and about the molecular mechanisms leading to activation. In the present study, we designed a new structural model of the TM domain and combined it with experimental mutagenesis experiments to investigate the interaction network that governs ligand binding and receptor activation. Our results suggest that this network involves the cluster of residues Arg(188) in TM2, Gln(380) in TM7, and Asn(229) in TM3. This cluster is expected to be altered upon VIP binding, because Arg(188) has been shown previously to interact with Asp(3) of VIP. Several point mutations at positions 188, 229, and 380 were experimentally characterized and were shown to severely affect VIP binding and/or VIP-mediated cAMP production. Double mutants built from reciprocal residue exchanges exhibit strong cooperative or anticooperative effects, thereby indicating the spatial proximity of residues Arg(188), Gln(380), and Asn(229). Because these residues are highly conserved in the GPCR-B family, they can moreover be expected to have a general role in mediating function.

Oxidized phospholipid inhibition of toll-like receptor (TLR) signaling is restricted to TLR2 and TLR4 roles for cd14, lps-binding protein, and md2 as targets for specificity of inhibition

The generation of reactive oxygen species is a central feature of inflammation that results in the oxidation of host phospholipids. Oxidized phospholipids, such as 1-palmitoyl-2-arachidonyl-sn-glycero-3-phosphorylcholine (OxPAPC), have been shown to inhibit signaling induced by bacterial lipopeptide or lipopolysac-charide (LPS), yet the mechanisms responsible for the inhibition of Toll-like receptor (TLR) signaling by OxPAPC remain incompletely understood. Here, we examined the mechanisms by which OxPAPC inhibits TLR signaling induced by diverse ligands in macrophages, smooth muscle cells, and epithelial cells. OxPAPC inhibited tumor necrosis factor- production, IB degradation, p38 MAPK phosphorylation, and NF-B-dependent reporter activation induced by stimulants of TLR2 and TLR4 (Pam3CSK4 and LPS) but not by stimulants of other TLRs (poly(I·C), flagellin, loxoribine, single-stranded RNA, or CpG DNA) in macrophages and HEK-293 cells transfected with respective TLRs and significantly reduced inflammatory responses in mice injected subcutaneously or intraperitoneally with Pam3CSK4. Serum proteins, including CD14 and LPS-binding protein, were identified as key targets for the specificity of TLR inhibition as supplementation with excess serum or recombinant CD14 or LBP reversed TLR2 inhibition by OxPAPC, whereas serum accessory proteins or expression of membrane CD14 potentiated signaling via TLR2 and TLR4 but not other TLRs. Binding experiments and functional assays identified MD2 as a novel additional target of OxPAPC inhibition of LPS signaling. Synthetic phospholipid oxidation products 1-palmitoyl-2-(5-oxovaleryl)-sn-glycero-3-phosphocholine and 1-palmitoyl-2-glutaryl-sn-glycero-3-phosphocholine inhibited TLR2 signaling from 30 µM. Taken together, these results suggest that oxidized phospholipid-mediated inhibition of TLR signaling occurs mainly by competitive interaction with accessory proteins that interact directly with bacterial lipids to promote signaling via TLR2 or TLR4.

Mapping interaction sites within the N-terminus of the calcitonin gene-related peptide receptor; the role of residues 23-60 of the calcitonin receptor-like receptor

The calcitonin receptor-like receptor (CLR) acts as a receptor for the calcitonin gene-related peptide (CGRP) but in order to recognize CGRP, it must form a complex with an accessory protein, receptor activity modifying protein 1 (RAMP1). Identifying the protein/protein and protein/ligand interfaces in this unusual complex would aid drug design. The role of the extreme N-terminus of CLR (Glu23-Ala60) was examined by an alanine scan and the results were interpreted with the help of a molecular model. The potency of CGRP at stimulating cAMP production was reduced at Leu41Ala, Gln45Ala, Cys48Ala and Tyr49Ala; furthermore, CGRP-induced receptor internalization at all of these receptors was also impaired. Ile32Ala, Gly35Ala and Thr37Ala all increased CGRP potency. CGRP specific binding was abolished at Leu41Ala, Ala44Leu, Cys48Ala and Tyr49Ala. There was significant impairment of cell surface expression of Gln45Ala, Cys48Ala and Tyr49Ala. Cys48 takes part in a highly conserved disulfide bond and is probably needed for correct folding of CLR. The model suggests that Gln45 and Tyr49 mediate their effects by interacting with RAMP1 whereas Leu41 and Ala44 are likely to be involved in binding CGRP. Ile32, Gly35 and Thr37 form a separate cluster of residues which modulate CGRP binding. The results from this study may be applicable to other family B GPCRs which can associate with RAMPs.