Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.

On the origin of the selectivity of plasmidic H-NS towards horizontally acquired DNA: Linking H-NS oligomerization and cooperative DNA binding

The nucleoid-associated protein H-NS is a global modulator of the expression of genes associated with adaptation to environmental changes. A variant of H-NS expressed in the R27 plasmid was previously shown to selectively modulate the expression of horizontally acquired genes, with minimal effects on core genes that are repressed by the chromosomal form of H-NS. Both H-NS proteins are formed by an oligomerization domain and a DNA-binding domain, which are connected by a linker that is highly flexible in the absence of DNA. We studied DNA binding by means of oligomer-forming chimeric proteins in which domains of the chromosomal and plasmidic variants are exchanged, as well as in monomeric truncated forms containing the DNA-binding domain and variable portions of the linker. Point mutations in the linker were also examined in full-length and truncated H-NS constructs. These experiments show that the linker region contributes to DNA binding affinity and that it is a main component of the distinct DNA binding properties of chromosomal and plasmidic H-NS. We propose that interactions between the linker and DNA limit the flexibility of the connection between H- NS oligomerization and DNA binding and provide an allosteric indirect readout mechanism to detect long- range distortions of DNA, thus enabling discrimination between core and horizontally acquired DNA.

Distinct OGT-Binding Sites Promote HCF-1 Cleavage.

Human HCF-1 (also referred to as HCFC-1) is a transcriptional co-regulator that undergoes a complex maturation process involving extensive O-GlcNAcylation and site-specific proteolysis. HCF-1 proteolysis results in two active, noncovalently associated HCF-1N and HCF-1C subunits that regulate distinct phases of the cell-division cycle. HCF-1 O-GlcNAcylation and site-specific proteolysis are both catalyzed by O-GlcNAc transferase (OGT), which thus displays an unusual dual enzymatic activity. OGT cleaves HCF-1 at six highly conserved 26 amino acid repeat sequences called HCF-1PRO repeats. Here we characterize the substrate requirements for OGT cleavage of HCF-1. We show that the HCF-1PRO-repeat cleavage signal possesses particular OGT-binding properties. The glutamate residue at the cleavage site that is intimately involved in the cleavage reaction specifically inhibits association with OGT and its bound cofactor UDP-GlcNAc. Further, we identify a novel OGT-binding sequence nearby the first HCF-1PRO-repeat cleavage signal that enhances cleavage. These results demonstrate that distinct OGT-binding sites in HCF-1 promote proteolysis, and provide novel insights into the mechanism of this unusual protease activity.

Highly functional binding surface for miniaturised solid-phase immunoassays. - A spot story -

Several bioaffinity assays are based on the detection of an analyte which is bound on a solid substrate via biochemical interaction. These so called solid phase assays are based on the adhesion of the primary binding partner on a solid surface, which then binds the analyte to be detected. In this thesis work a novel solid phase based assay technology, known as spot technology, was developed. The spot technology is based on combination of high-capacity solid phases, concentrated in a spot format, utilising modified streptavidin molecules and recombinant antibody fragments. The reduction of the solid phase binding surface to a size of a spot enabled denser binding of the target molecules, providing improved signal intensities and signal-to-background ratio when applied in different solid phase immunoassays. Streptavidin-biotin interactions are commonly utilised in numerous different bioaffinity assays and the ultimate nature of streptavidin to bind biotin is among the strongest non-covalent interaction reported between two biomolecules. In this study native core streptavidin was chemically modified to provide polymerised streptavidin molecules with altered adsorption properties. These streptavidin conjugates, when coated onto polystyrene surface, provided enhanced biotin binding capacity and surface stability when compared to a reference coating constructed with native streptavidin. Furthermore, the combination of chemically modified streptavidin, sitespecifically biotinylated antibody fragments and the spot coating technology provided highly dense solid phase coating with improved binding properties. The performance of the spot assay technology was further demonstrated in different immunoassay configurations. Human thyroid stimulating hormone (TSH) and human cardiac troponin I (cTnI) were used as model analytes to show the applicability of the highly sensitive spot-based solid-phase immunoassay for detection of very low levels of analytes. It was demonstrated that the spot technology provided an assay concept with enhanced sensitivity and short turn-around times, characteristics that are highly suitable for point-of-care applications.

Directed evolution of human dihydrofolate reductase: towards a better understanding of binding at the active site

La dihydrofolate réductase humaine (DHFRh) est une enzyme essentielle à la prolifération cellulaire, ce qui en fait une cible de choix pour le traitement de différents cancers. À cet effet, plusieurs inhibiteurs spécifiques de la DHFRh, les antifolates, ont été mis au point : le méthotrexate (MTX) et le pemetrexed (PMTX) en sont de bons exemples. Malgré l’efficacité clinique certaine de ces antifolates, le développement de nouveaux traitements s’avère nécessaire afin de réduire les effets secondaires liés à leur utilisation. Enfin, dans l’optique d’orienter la synthèse de nouveaux composés inhibiteurs des DHFRh, une meilleure connaissance des interactions entre les antifolates et leur enzyme cible est primordiale. À l’aide de l’évolution dirigée, il a été possible d’identifier des mutants de la DHFRh pour lesquels l’affinité envers des antifolates cliniquement actifs se voyait modifiée. La mutagenèse dite ¬¬de saturation a été utilisée afin de générer des banques de mutants présentant une diversité génétique au niveau des résidus du site actif de l’enzyme d’intérêt. De plus, une nouvelle méthode de criblage a été mise au point, laquelle s’est avérée efficace pour départager les mutations ayant entrainé une résistance aux antifolates et/ou un maintient de l’activité enzymatique envers son substrat natif, soient les phénotypes d’activité. La méthode de criblage consiste dans un premier temps en une sélection bactérienne à haut débit, puis dans un second temps en un criblage sur plaques permettant d’identifier les meilleurs candidats. Plusieurs mutants actifs de la DHFRh, résistants aux antifolates, ont ainsi pu être identifiés et caractérisés lors d’études de cinétique enzymatique (kcat et IC50). Sur la base de ces résultats cinétiques, de la modélisation moléculaire et des données structurales de la littérature, une étude structure-activité a été effectuée. En regardant quelles mutations ont les effets les plus significatif sur la liaison, nous avons commencé à construire un carte moléculaire des contacts impliqués dans la liaison des ligands. Enfin, des connaissances supplémentaires sur les propriétés spécifiques de liaison ont put être acquises en variant l’inhibiteur testé, permettant ainsi une meilleure compréhension du phénomène de discrimination du ligand.

Puroindoline-b mutations control the lipid binding interactions in mixed puroindoline-a: puroindoline-b systems

The interactions have been investigated of puroindoline-a (Pin-a) and mixed protein systems of Pin-a and wild-type puroindoline-b (Pin-b+) or puroindoline-b mutants (G46S mutation (Pin bH) or W44R mutation (Pin-bS)) with condensed phase monolayers of an anionic phospholipid (L-α-dipalmitoylphosphatidyl-dl-glycerol (DPPG)) at the air/water interface. The interactions of the mixed systems were studied at three different concentration ratios of Pin-a:Pin-b, namely 3:1, 1:1 and 1:3 in order to establish any synergism in relation to lipid binding properties. Surface pressure measurements revealed that Pin-a interaction with DPPG monolayers led to an equilibrium surface pressure increase of 8.7 ± 0.6 mN m-1. This was less than was measured for Pin-a:Pin-b+ (9.6 to 13.4 mN m-1), but was significantly more than was measured for Pin-a:Pin-bH (4.0 to 6.2 mN m-1) or Pin-a:Pin-bS (3.8 to 6.3 mN m-1) over the complete range of concentration ratio. Consequently, surface pressure increases were shown to correlate to endosperm hardness phenotype, with puroindolines present in hard-textured wheat varieties yielding lower equilibrium surface pressure changes. Integrated amide I peak areas from corresponding external reflectance Fourier-transform infrared (ER-FTIR) spectra, used to indicate levels of protein adsorption to the lipid monolayers, showed that differences in adsorbed amount were less significant. The data therefore suggest that Pin-b mutants having single residue substitutions within their tryptophan-rich loop that are expressed in some hard-textured wheat varieties influence the degree of penetration of Pin-a and Pin-b into anionic phospholipid films. These findings highlight the key role of the tryptophan-rich loop in puroindoline-lipid interactions.

Oligomers of D-2 dopamine receptors - Evidence from ligand binding

There is increasing evidence that G protein-coupled receptors form oligomers and that this might be important for their function. We have studied this phenomenon for the D-2 dopamine receptor and have shown-using a variety of biochemical and biophysical techniques-that this receptor forms dimers or higher-order oligomers. Using ligand-binding studies, we have also found evidence that this oligomer formation has functional relevance. Thus, for the receptor expressed in either CHO cells or Sf 9 insect cells, the binding properties of several radioligands (in saturation, competition, and dissociation assays) do not conform to those expected for a monomeric receptor with a single binding site. We propose that the receptors exist in oligomers with homotropic and heterotropic negatively cooperative interactions between ligands.

Oligomers of D-2 dopamine receptors: evidence from ligand binding

There is increasing evidence that G protein-coupled receptors form oligomers and that this might be important for their function. We have studied this phenomenon for the D-2 dopamine receptor and have shown-using a variety of biochemical and biophysical techniques-that this receptor forms dimers or higher-order oligomers. Using ligand-binding studies, we have also found evidence that this oligomer formation has functional relevance. Thus, for the receptor expressed in either CHO cells or Sf 9 insect cells, the binding properties of several radioligands (in saturation, competition, and dissociation assays) do not conform to those expected for a monomeric receptor with a single binding site. We propose that the receptors exist in oligomers with homotropic and heterotropic negatively cooperative interactions between ligands

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Structure and Calcium-Binding Activity of LipL32, the Major Surface Antigen of Pathogenic Leptospira sp.

Leptospixosis, a spirochaetal zoonotic disease caused by Leptospira, has been recognized as an important emerging infectious disease. LipL32 is the major exposed outer membrane protein found exclusively in pathogenic leptospires, where it accounts for up to 75% of the total outer membrane proteins. It is highly immunogenic, and recent studies have implicated LipL32 as an extracellular matrix binding protein, interacting with collagens, fibronectin, and laminin. In order to better understand the biological role and the structural requirements for the function of this important lipoprotein, we have determined the 2.25-angstrom-resolution structure of recombinant LipL32 protein corresponding to residues 21-272 of the wild-type protein (LipL32(21-272)). The LipL32(21-272) monomer is made of a jelly-roll fold core from which several peripheral secondary structures protrude. LipL32(21-272) is structurally similar to several other jelly-roll proteins, some of which bind calcium ions and extracellular matrix proteins. Indeed, spectroscopic data (circular dichroism, intrinsic tryptophan fluorescence, and extrinsic 1-amino-2-naphthol-4-sulfonic acid fluorescence) confirmed the calcium-binding properties of LipL32(21-272). Ca(2+) binding resulted in a significant increase in the thermal stability of the protein, and binding was specific for Ca(2+) as no structural or stability perturbations were observed for Mg(2+), Zn(2+), or Cu(2+). Careful examination of the crystal lographic structure suggests the locations of putative regions that could mediate Ca(2+) binding as well as binding to other interacting host proteins, such as collagens, fibronectin, and lamixidn. (C) 2009 Elsevier Ltd. All rights reserved.

The Role of Nanometer-Scaled Ligand Patterns in Polyvalent Binding by Large Mannan-Binding Lectin Oligomers

Mannan-binding lectin (MBL) is an important protein of the innate immune system and protects the body against infection through opsonization and activation of the complement system on surfaces with an appropriate presentation of carbohydrate ligands. The quaternary structure of human MBL is built from oligomerization of structural units into polydisperse complexes typically with three to eight structural units, each containing three lectin domains. Insight into the connection between the structure and ligand-binding properties of these oligomers has been lacking. In this article, we present an analysis of the binding to neoglycoprotein-coated surfaces by size-fractionated human MBL oligomers studied with small-angle x-ray scattering and surface plasmon resonance spectroscopy. The MBL oligomers bound to these surfaces mainly in two modes, with dissociation constants in the micro to nanomolar order. The binding kinetics were markedly influenced by both the density of ligands and the number of ligand-binding domains in the oligomers. These findings demonstrated that the MBL-binding kinetics are critically dependent on structural characteristics on the nanometer scale, both with regard to the dimensions of the oligomer, as well as the ligand presentation on surfaces. Therefore, our work suggested that the surface binding of MBL involves recognition of patterns with dimensions on the order of 10-20 nm. The recent understanding that the surfaces of many microbes are organized with structural features on the nanometer scale suggests that these properties of MBL ligand recognition potentially constitute an important part of the pattern-recognition ability of these polyvalent oligomers. The Journal of Immunology, 2012, 188: 1292-1306.

Three-dimensional electron microscopy of myriapod hemocyanin, planorbid snail acetylcholine-binding protein and cyanobacterial Vipp1

Three-dimensional electron microscopy (3-D EM) provides a framework for the analysis of large protein quaternary structures. The advantage over the generally higher resolving meth- od of X-ray crystallography is the embedding of the proteins in their physiological environ- ment. However, results of the two methods can be combined to obtain superior structural information. In this work, three different protein types – (i) Myriapod hemocyanin, (ii) vesi- cle-inducing protein in plastids 1 (Vipp1) and (iii) acetylcholine-binding protein (AChBP) – were structurally analyzed by 2-D and 3-D EM and, where possible, functionally interpreted.rnMyriapod hemocyanins have been previously shown to be 6x6-meric assemblies that, in case of Scutigera coleoptrata hemocyanin (ScoHc), show two 3x6-mer planes whith a stag- gering angle of approximately 60°. Here, previously observed structural differences between oxy- and deoxy-ScoHc could be substantiated. A 4° rotation between hexamers of two dif- ferent 3x6-mer planes was measured, which originates at the most central inter-hexamer in- terface. Further information about allosteric behaviour in myriapod hemocyanin was gained by analyzing Polydesmus angustus hemocyanin (PanHc), which shows a stable 3x6-mer and divergent histidine patterns in the inter-hexamer interfaces when compared to ScoHc. Both findings would conclusively explain the very different oxygen binding properties of chilopod and diplopod hemocyanin.rnVipp1 is a protein found in cyanobacteria and higher plants which is essential for thyla- koid membrane function and forms highly variable ring-shaped structures. In the course of this study, the first 3-D analysis of Vipp1 was conducted and yielded reconstructions of six differently sized Vipp1 rings from negatively stained images at resolutions between 20 to 30 Å. Furthermore, mutational analyses identified specific N-terminal amino acids that are essential for ring formation. On the basis of these analyses and previously published results, a hypothetical model of the Vipp1 tertiary and quaternary structure was generated.rnAChBP is a water-soluble protein in the hemolymph of mollusks. It is a structural and functional homologue of the ligand-binding domain of nicotinic acetylcholine receptors. For the freshwater snail Biomphalaria glabrata, we previously described two types of AChBP (BgAChBP1 and BgAChBP2). In this work, a 6 Å 3-D reconstruction of native BgAChBP is presented, which shows a dodecahedral assembly that is unprecedented for an AChBP. Single particle analysis of recombinantely expressed BgAChBP types led to preliminary results show- ing a dodecahedral assembly of BgAChBP1 and a dipentameric assembly of BgAChBP2. This indicates divergent biological functions of the two types.

Measuring substrate binding and affinity of purified membrane transport proteins using the scintillation proximity assay

The scintillation proximity assay (SPA) is a rapid radioligand binding assay. Upon binding of radioactively labeled ligands (here L-[(3)H]arginine or D-[(3)H]glucose) to acceptor proteins immobilized on fluoromicrospheres (containing the scintillant), a light signal is stimulated and measured. The application of SPA to purified, detergent-solubilized membrane transport proteins allows substrate-binding properties to be assessed (e.g., substrate specificity and affinity), usually within 1 d. Notably, the SPA makes it possible to study specific transporters without interference from other cellular components, such as endogenous transporters. Reconstitution of the target transporter into proteoliposomes is not required. The SPA procedure allows high sample throughput and simple sample handling without the need for washing or separation steps: components are mixed in one well and the signal is measured directly after incubation. Therefore, the SPA is an excellent tool for high-throughput screening experiments, e.g., to search for substrates and inhibitors, and it has also recently become an attractive tool for drug discovery.

Antisense properties of tricyclo-DNA.

Tricyclo (tc)-DNA belongs to the class of conformationally constrained DNA analogs that show enhanced binding properties to DNA and RNA. We prepared tc-oligonucleotides up to 17 nt in length, and evaluated their binding efficiency and selectivity towards complementary RNA, their biological stability in serum, their RNase H inducing potential and their antisense activity in a cellular assay. Relative to RNA or 2'-O-Me-phosphorothioate (PS)-RNA, fully modified tc-oligodeoxynucleotides, 10-17 nt in length, show enhanced selectivity and enhanced thermal stability by approximately 1 degrees C/modification in binding to RNA targets. Tricyclodeoxyoligonucleotides are completely stable in heat-deactivated fetal calf serum at 37 degree C. Moreover, tc-DNA-RNA duplexes are not substrates for RNase H. To test for antisense effects in vivo, we used HeLa cell lines stably expressing the human beta-globin gene with two different point mutations in the second intron. These mutations lead to the inclusion of an aberrant exon in beta-globin mRNA. Lipofectamine-mediated delivery of a 17mer tc-oligodeoxynucleotide complementary to the 3'-cryptic splice site results in correction of aberrant splicing already at nanomolar concentrations with up to 100-fold enhanced efficiency relative to a 2'-O-Me-PS-RNA oligonucleotide of the same length and sequence. In contrast to 2'-O-Me-PS-RNA, tc-DNA shows antisense activity even in the absence of lipofectamine, albeit only at much higher oligonucleotide concentrations.

Positive modulation of synaptic and extrasynaptic GABAA receptors by an antagonist of the high affinity benzodiazepine binding site.

GABAA receptors are the major inhibitory neurotransmitter receptors in the brain and are the target for many clinically important drugs such as the benzodiazepines. Benzodiazepines act at the high-affinity binding site at the α+/γ- subunit interface. Previously, an additional low affinity binding site for diazepam located in the transmembrane (TM) domain has been described. The compound SJM-3 was recently identified in a prospective screening of ligands for the benzodiazepine binding site and investigated for its site of action. We determined the binding properties of SJM-3 at GABAA receptors recombinantly expressed in HEK-cells using radioactive ligand binding assays. Impact on function was assessed in Xenopus laevis oocytes with electrophysiological experiments using the two-electrode voltage clamp method. SJM-3 was shown to act as an antagonist at the α+/γ- site. At the same time it strongly potentiated GABA currents via the binding site for diazepam in the transmembrane domain. Mutation of a residue in M2 of the α subunit strongly reduced receptor modulation by SJM-3 and a homologous mutation in the β subunit abolished potentiation. SJM-3 acts as a more efficient modulator than diazepam at the site in the trans-membrane domain. In contrast to low concentrations of benzodiazepines, SJM-3 modulates both synaptic and extrasynaptic receptors. A detailed exploration of the membrane site may provide the basis for the design and identification of subtype-selective modulatory drugs.