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.

Glial hyaluronate-binding protein expression in aggregating brain cell cultures.

We analyzed the expression of glial hyaluronate-binding protein (GHAP), an integral component of the extracellular matrix, in aggregating brain cell cultures of fetal rat telencephalon using immunofluorescence. GHAP immunoreactivity appeared after 1 week in culture, simultaneous with the first deposits of myelin basic protein, and showed a development-dependent increase. Comparison of glia-enriched and neuron-enriched cultures showed that only glial cells express GHAP. Three peptide growth factors, epidermal growth factor, fibroblast growth factor and platelet-derived growth factor, which are known to stimulate the differentiation of glial cells, modulated the deposit of GHAP immunoreactivity. The 3-dimensional structure of aggregate cultures promoted GHAP deposition, suggesting that cell-cell interactions are required for extracellular matrix formation. Furthermore GHAP production seemed to depend on the developmental stage of the glial cells.

Inhibition of the shade avoidance response by formation of non-DNA binding bHLH heterodimers.

In shade-intolerant plants such as Arabidopsis, a reduction in the red/far-red (R/FR) ratio, indicative of competition from other plants, triggers a suite of responses known as the shade avoidance syndrome (SAS). The phytochrome photoreceptors measure the R/FR ratio and control the SAS. The phytochrome-interacting factors 4 and 5 (PIF4 and PIF5) are stabilized in the shade and are required for a full SAS, whereas the related bHLH factor HFR1 (long hypocotyl in FR light) is transcriptionally induced by shade and inhibits this response. Here we show that HFR1 interacts with PIF4 and PIF5 and limits their capacity to induce the expression of shade marker genes and to promote elongation growth. HFR1 directly inhibits these PIFs by forming non-DNA-binding heterodimers with PIF4 and PIF5. Our data indicate that PIF4 and PIF5 promote SAS by directly binding to G-boxes present in the promoter of shade marker genes, but their action is limited later in the shade when HFR1 accumulates and forms non-DNA-binding heterodimers. This negative feedback loop is important to limit the response of plants to shade.

Ligands for pheromone-sensing neurons are not conformationally activated odorant binding proteins.

Pheromones form an essential chemical language of intraspecific communication in many animals. How olfactory systems recognize pheromonal signals with both sensitivity and specificity is not well understood. An important in vivo paradigm for this process is the detection mechanism of the sex pheromone (Z)-11-octadecenyl acetate (cis-vaccenyl acetate [cVA]) in Drosophila melanogaster. cVA-evoked neuronal activation requires a secreted odorant binding protein, LUSH, the CD36-related transmembrane protein SNMP, and the odorant receptor OR67d. Crystallographic analysis has revealed that cVA-bound LUSH is conformationally distinct from apo (unliganded) LUSH. Recombinantly expressed mutant versions of LUSH predicted to enhance or diminish these structural changes produce corresponding alterations in spontaneous and/or cVA-evoked activity when infused into olfactory sensilla, leading to a model in which the ligand for pheromone receptors is not free cVA, but LUSH that is "conformationally activated" upon cVA binding. Here we present evidence that contradicts this model. First, we demonstrate that the same LUSH mutants expressed transgenically affect neither basal nor pheromone-evoked activity. Second, we compare the structures of apo LUSH, cVA/LUSH, and complexes of LUSH with non-pheromonal ligands and find no conformational property of cVA/LUSH that can explain its proposed unique activated state. Finally, we show that high concentrations of cVA can induce neuronal activity in the absence of LUSH, but not SNMP or OR67d. Our findings are not consistent with the model that the cVA/LUSH complex acts as the pheromone ligand, and suggest that pheromone molecules alone directly activate neuronal receptors.

Novel FOXF1 Mutations in Sporadic and Familial Cases of Alveolar Capillary Dysplasia with Misaligned Pulmonary Veins Imply a Role for its DNA Binding Domain.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACD/MPV) is a rare and lethal developmental disorder of the lung defined by a constellation of characteristic histopathological features. Nonpulmonary anomalies involving organs of gastrointestinal, cardiovascular, and genitourinary systems have been identified in approximately 80% of patients with ACD/MPV. We have collected DNA and pathological samples from more than 90 infants with ACD/MPV and their family members. Since the publication of our initial report of four point mutations and 10 deletions, we have identified an additional 38 novel nonsynonymous mutations of FOXF1 (nine nonsense, seven frameshift, one inframe deletion, 20 missense, and one no stop). This report represents an up to date list of all known FOXF1 mutations to the best of our knowledge. Majority of the cases are sporadic. We report four familial cases of which three show maternal inheritance, consistent with paternal imprinting of the gene. Twenty five mutations (60%) are located within the putative DNA-binding domain, indicating its plausible role in FOXF1 function. Five mutations map to the second exon. We identified two additional genic and eight genomic deletions upstream to FOXF1. These results corroborate and extend our previous observations and further establish involvement of FOXF1 in ACD/MPV and lung organogenesis.

Identification and mutational analyses of phosphorylation sites of the calcineurin-binding protein CbpA and the identification of domains required for calcineurin binding in Aspergillus fumigatus.

Calcineurin is a key protein phosphatase required for hyphal growth and virulence in Aspergillus fumigatus, making it an attractive antifungal target. However, currently available calcineurin inhibitors, FK506 and cyclosporine A, are immunosuppressive, limiting usage in the treatment of patients with invasive aspergillosis. Therefore, the identification of endogenous inhibitors of calcineurin belonging to the calcipressin family is an important parallel strategy. We previously identified the gene cbpA as the A. fumigatus calcipressin member and showed its involvement in hyphal growth and calcium homeostasis. However, the mechanism of its activation/inhibition through phosphorylation and its interaction with calcineurin remains unknown. Here we show that A. fumigatus CbpA is phosphorylated at three distinct domains, including the conserved SP repeat motif (phosphorylated domain-I; PD-I), a filamentous fungal-specific domain (PD-II), and the C-terminal CIC motif (Calcipressin Inhibitor of Calcineurin; PD-III). While mutation of three phosphorylated residues (Ser208, Ser217, Ser223) in the PD-II did not affect CbpA function in vivo, mutation of the two phosphorylated serines (Ser156, Ser160) in the SP repeat motif caused reduced hyphal growth and sensitivity to oxidative stress. Mutational analysis in the key domains in calcineurin A (CnaA) and proteomic interaction studies confirmed the requirement of PxIxIT motif-binding residues (352-NIR-354) and the calcineurin B (CnaB)-binding helix residue (V371) for the binding of CbpA to CnaA. Additionally, while the calmodulin-binding residues (442-RVF-444) did not affect CbpA binding to CnaA, three mutations (T359P, H361L, and L365S) clustered between the CnaA catalytic and the CnaB-binding helix were also required for CbpA binding. This is the first study to analyze the phosphorylation status of calcipressin in filamentous fungi and identify the domains required for binding to calcineurin.

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: Proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of pH and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups. © 2011 American Institute of Physics.

Lipid rafts act as specialised domains for tetanus toxin binding and internalisation in neurons

Tetanus (TeNT) is a zinc protease that blocks neurotransmission by cleaving the synaptic protein vesicle-associated membrane protein/synaptobrevin. Although its intracellular catalytic activity is well established, the mechanism by which this neurotoxin interacts with the neuronal surface is not known. In this study, we characterize p15s, the first plasma membrane TeNT binding proteins and we show that they are glycosylphosphatidylinositol-anchored glycoproteins in nerve growth factor (NGF)-differentiated PC12 cells, spinal cord cells, and purified motor neurons. We identify p15 as neuronal Thy-1 in NGF-differentiated PC12 cells. Fluorescence lifetime imaging microscopy measurements confirm the close association of the binding domain of TeNT and Thy-1 at the plasma membrane. We find that TeNT is recruited to detergent-insoluble lipid microdomains on the surface of neuronal cells. Finally, we show that cholesterol depletion affects a raft subpool and blocks the internalization and intracellular activity of the toxin. Our results indicate that TeNT interacts with target cells by binding to lipid rafts and that cholesterol is required for TeNT internalization and/or trafficking in neurons.

Photocontrolled DNA Binding of a Receptor-Targeted Organometallic Ruthenium(II) Complex

A photoactivated ruthenium(II) arene complex has been conjugated to two receptor-binding peptides, a dicarba analogue of octreotide and the Arg-Gly-Asp (RGD) tripeptide. These peptides can act as"tumor-targeting devices" since their receptors are overexpressed on the membranes of tumor cells. Both ruthenium-peptide conjugates are stable in aqueous solution in the dark, but upon irradiation with visible light, the pyridyl-derivatized peptides were selectively photodissociated from the ruthenium complex, as inferred by UV-vis and NMR spectroscopy. Importantly, the reactive aqua species generated from the conjugates, [(η6-p-cym)Ru(bpm)(H2O)]2+, reacted with the model DNA nucleobase 9-ethylguanine as well as with guanines of two DNA sequences, 5′dCATGGCT and 5′dAGCCATG. Interestingly, when irradiation was performed in the presence of the oligonucleotides, a new ruthenium adduct involving both guanines was formed as a consequence of the photodriven loss of p-cymene from the two monofunctional adducts. The release of the arene ligand and the formation of a ruthenated product with a multidentate binding mode might have important implications for the biological activity of such photoactivated ruthenium(II) arene complexes. Finally, photoreactions with the peptide-oligonucleotide hybrid, Phac-His-Gly-Met-linker-p5′dCATGGCT, also led to arene release and to guanine adducts, including a GG chelate. The lack of interaction with the peptide fragment confirms the preference of such organometallic ruthenium(II) complexes for guanine over other potential biological ligands, such as histidine or methionine amino acids.

Differential Dimerization of Variants Linked to Enhanced S-Cone Sensitivity Syndrome (ESCS) Located in the NR2E3 Ligand-Binding Domain.

NR2E3 encodes the photoreceptor-specific nuclear hormone receptor that acts as a repressor of cone-specific gene expression in rod photoreceptors, and as an activator of several rod-specific genes. Recessive variants located in the ligand-binding domain (LBD) of NR2E3 cause enhanced short wavelength sensitive- (S-) cone syndrome (ESCS), a retinal degeneration characterized by an excess of S-cones and non-functional rods. We analyzed the dimerization properties of NR2E3 and the effect of disease-causing LBD missense variants by bioluminescence resonance energy transfer (BRET(2) ) protein interaction assays. Homodimerization was not affected in presence of p.A256V, p.R039G, p.R311Q, and p.R334G variants, but abolished in presence of p.L263P, p.L336P, p.L353V, p.R385P, and p.M407K variants. Homology modeling predicted structural changes induced by NR2E3 LBD variants. NR2E3 LBD variants did not affect interaction with CRX, but with NRL and rev-erbα/NR1D1. CRX and NRL heterodimerized more efficiently together, than did either with NR2E3. NR2E3 did not heterodimerize with TLX/NR2E1 and RXRα/NR2C1. The identification of a new compound heterozygous patient with detectable rod function, who expressed solely the p.A256V variant protein, suggests a correlation between LBD variants able to form functional NR2E3 dimers and atypical mild forms of ESCS with residual rod function.

Microscopic-macroscopic approach for binding energies with the Wigner-Kirkwood method. II. Deformed nuclei

The binding energies of deformed even-even nuclei have been analyzed within the framework of a recently proposed microscopic-macroscopic model. We have used the semiclassical Wigner-Kirkwood ̄h expansion up to fourth order, instead of the usual Strutinsky averaging scheme, to compute the shell corrections in a deformed Woods-Saxon potential including the spin-orbit contribution. For a large set of 561 even-even nuclei with Z 8 and N 8, we find an rms deviation from the experiment of 610 keV in binding energies, comparable to the one found for the same set of nuclei using the finite range droplet model of Moller and Nix (656 keV). As applications of our model, we explore its predictive power near the proton and neutron drip lines as well as in the superheavy mass region. Next, we systematically explore the fourth-order Wigner-Kirkwood corrections to the smooth part of the energy. It is found that the ratio of the fourth-order to the second-order corrections behaves in a very regular manner as a function of the asymmetry parameter I=(N−Z)/A. This allows us to absorb the fourth-order corrections into the second-order contributions to the binding energy, which enables us us to simplify and speed up the calculation of deformed nuclei.