Modeling transcription factor binding events to DNA using a random walker/jumper representation on a 1D/2D lattice with different affinity sites

Surviving in a diverse environment requires corresponding organism responses. At the cellular level, such adjustment relies on the transcription factors (TFs) which must rapidly find their target sequences amidst a vast amount of non-relevant sequences on DNA molecules. Whether these transcription factors locate their target sites through a 1D or 3D pathway is still a matter of speculation. It has been suggested that the optimum search time is when the protein equally shares its search time between 1D and 3D diffusions. In this paper, we study the above problem using a Monte Carlo simulation by considering a very simple physical model. A 1D strip, representing a DNA, with a number of low affinity sites, corresponding to non-target sites, and high affinity sites, corresponding to target sites, is considered and later extended to a 2D strip. We study the 1D and 3D exploration pathways, and combinations of the two modes by considering three different types of molecules: a walker that randomly walks along the strip with no dissociation; a jumper that represents dissociation and then re-association of a TF with the strip at later time at a distant site; and a hopper that is similar to the jumper but it dissociates and then re-associates at a faster rate than the jumper. We analyze the final probability distribution of molecules for each case and find that TFs can locate their targets fast enough even if they spend 15% of their search time diffusing freely in the solution. This indeed agrees with recent experimental results obtained by Elf et al. 2007 and is in contrast with theoretical expectation.Comment: 24 pages, 9 figure

Classes of fast and specific search mechanisms for proteins on DNA

Problems of search and recognition appear over different scales in biological systems. In this review we focus on the challenges posed by interactions between proteins, in particular transcription factors, and DNA and possible mechanisms which allow for a fast and selective target location. Initially we argue that DNA-binding proteins can be classified, broadly, into three distinct classes which we illustrate using experimental data. Each class calls for a different search process and we discuss the possible application of different search mechanisms proposed over the years to each class. The main thrust of this review is a new mechanism which is based on barrier discrimination. We introduce the model and analyze in detail its consequences. It is shown that this mechanism applies to all classes of transcription factors and can lead to a fast and specific search. Moreover, it is shown that the mechanism has interesting transient features which allow for stability at the target despite rapid binding and unbinding of the transcription factor from the target.Comment: 65 pages, 23 figure

Salerno's model of DNA reanalysed: could solitons have biological significance?

We investigate the sequence-dependent behaviour of localised excitations in a toy, nonlinear model of DNA base-pair opening originally proposed by Salerno. Specifically we ask whether ``breather'' solitons could play a role in the facilitated location of promoters by RNA polymerase. In an effective potential formalism, we find excellent correlation between potential minima and {\em Escherichia coli} promoter recognition sites in the T7 bacteriophage genome. Evidence for a similar relationship between phage promoters and downstream coding regions is found and alternative reasons for links between AT richness and transcriptionally-significant sites are discussed. Consideration of the soliton energy of translocation provides a novel dynamical picture of sliding: steep potential gradients correspond to deterministic motion, while ``flat'' regions, corresponding to homogeneous AT or GC content, are governed by random, thermal motion. Finally we demonstrate an interesting equivalence between planar, breather solitons and the helical motion of a sliding protein ``particle'' about a bent DNA axis.Comment: Latex file 20 pages, 5 figures. Manuscript of paper to appear in J. Biol. Phys., accepted 02/09/0

Cytosine-to-Uracil Deamination by SssI DNA Methyltransferase

The prokaryotic DNA(cytosine-5)methyltransferase M.SssI shares the specificity of eukaryotic DNA methyltransferases (CG) and is an important model and experimental tool in the study of eukaryotic DNA methylation. Previously, M.SssI was shown to be able to catalyze deamination of the target cytosine to uracil if the methyl donor S-adenosyl-methionine (SAM) was missing from the reaction. To test whether this side-activity of the enzyme can be used to distinguish between unmethylated and C5-methylated cytosines in CG dinucleotides, we re-investigated, using a sensitive genetic reversion assay, the cytosine deaminase activity of M.SssI. Confirming previous results we showed that M.SssI can deaminate cytosine to uracil in a slow reaction in the absence of SAM and that the rate of this reaction can be increased by the SAM analogue 5’-amino-5’-deoxyadenosine. We could not detect M.SssI-catalyzed deamination of C5-methylcytosine (m5C). We found conditions where the rate of M.SssI mediated C-to-U deamination was at least 100-fold higher than the rate of m5C-to-T conversion. Although this difference in reactivities suggests that the enzyme could be used to identify C5-methylated cytosines in the epigenetically important CG dinucleotides, the rate of M.SssI mediated cytosine deamination is too low to become an enzymatic alternative to the bisulfite reaction. Amino acid replacements in the presumed SAM binding pocket of M.SssI (F17S and G19D) resulted in greatly reduced methyltransferase activity. The G19D variant showed cytosine deaminase activity in E. coli, at physiological SAM concentrations. Interestingly, the C-to-U deaminase activity was also detectable in an E. coli ung+ host proficient in uracil excision repair

Mechanistic Insights on the Inhibition of C5 DNA Methyltransferases by Zebularine

In mammals DNA methylation occurs at position 5 of cytosine in a CpG context and regulates gene expression. It plays an important role in diseases and inhibitors of DNA methyltransferases (DNMTs)—the enzymes responsible for DNA methylation—are used in clinics for cancer therapy. The most potent inhibitors are 5-azacytidine and 5-azadeoxycytidine. Zebularine (1-(β-D-ribofuranosyl)-2(1H)- pyrimidinone) is another cytidine analog described as a potent inhibitor that acts by forming a covalent complex with DNMT when incorporated into DNA. Here we bring additional experiments to explain its mechanism of action. First, we observe an increase in the DNA binding when zebularine is incorporated into the DNA, compared to deoxycytidine and 5-fluorodeoxycytidine, together with a strong decrease in the dissociation rate. Second, we show by denaturing gel analysis that the intermediate covalent complex between the enzyme and the DNA is reversible, differing thus from 5-fluorodeoxycytidine. Third, no methylation reaction occurs when zebularine is present in the DNA. We confirm that zebularine exerts its demethylation activity by stabilizing the binding of DNMTs to DNA, hindering the methylation and decreasing the dissociation, thereby trapping the enzyme and preventing turnover even at other sites

Abstraction in ecology : reductionism and holism as complementary heuristics

In addition to their core explanatory and predictive assumptions, scientific models include simplifying assumptions, which function as idealizations, approximations, and abstractions. There are methods to investigate whether simplifying assumptions bias the results of models, such as robustness analyses. However, the equally important issue - the focus of this paper - has received less attention, namely, what are the methodological and epistemic strengths and limitations associated with different simplifying assumptions. I concentrate on one type of simplifying assumption, the use of mega parameters as abstractions in ecological models. First, I argue that there are two kinds of mega parameters qua abstractions, sufficient parameters and aggregative parameters, which have gone unnoticed in the literature. The two are associated with different heuristics, holism and reductionism, which many view as incompatible. Second, I will provide a different analysis of abstractions and the associated heuristics than previous authors. Reductionism and holism and the accompanying abstractions have different methodological and epistemic functions, strengths, and limitations, and the heuristics should be viewed as providing complementary research perspectives of cognitively limited beings. This is then, third, used as a premise to argue for epistemic and methodological pluralism in theoretical ecology. Finally, the presented taxonomy of abstractions is used to comment on the current debate whether mechanistic accounts of explanation are compatible with the use of abstractions. This debate has suffered from an abstract discussion of abstractions. With a better taxonomy of abstractions the debate can be resolved.Peer reviewe

The structure of SgrAI bound to DNA; recognition of an 8 base pair target

The three-dimensional X-ray crystal structure of the ‘rare cutting’ type II restriction endonuclease SgrAI bound to cognate DNA is presented. SgrAI forms a dimer bound to one duplex of DNA. Two Ca2+ bind in the enzyme active site, with one ion at the interface between the protein and DNA, and the second bound distal from the DNA. These sites are differentially occupied by Mn2+, with strong binding at the protein–DNA interface, but only partial occupancy of the distal site. The DNA remains uncleaved in the structures from crystals grown in the presence of either divalent cation. The structure of the dimer of SgrAI is similar to those of Cfr10I, Bse634I and NgoMIV, however no tetrameric structure of SgrAI is observed. DNA contacts to the central CCGG base pairs of the SgrAI canonical target sequence (CR|CCGGYG, | marks the site of cleavage) are found to be very similar to those in the NgoMIV/DNA structure (target sequence G|CCGGC). Specificity at the degenerate YR base pairs of the SgrAI sequence may occur via indirect readout using DNA distortion. Recognition of the outer GC base pairs occurs through a single contact to the G from an arginine side chain located in a region unique to SgrAI