Permalink
Cannot retrieve contributors at this time
Fetching contributors…
Cannot retrieve contributors at this time
| <?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN" "JATS-archivearticle1.dtd"><article article-type="research-article" dtd-version="1.1d1" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="hwp">eLife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">00668</article-id><article-id pub-id-type="doi">10.7554/eLife.00668</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Biochemistry</subject></subj-group><subj-group subj-group-type="heading"><subject>Genomics and evolutionary biology</subject></subj-group></article-categories><title-group><article-title>The autoregulation of a eukaryotic DNA transposon</article-title></title-group><contrib-group><contrib contrib-type="author" equal-contrib="yes" id="author-4326"><name><surname>Claeys Bouuaert</surname><given-names>Corentin</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="fn" rid="pa1">‡</xref><xref ref-type="other" rid="par-1"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-4327"><name><surname>Lipkow</surname><given-names>Karen</given-names></name><xref ref-type="aff" rid="aff2"/><xref ref-type="aff" rid="aff3"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="other" rid="par-2"/><xref ref-type="other" rid="par-5"/><xref ref-type="other" rid="par-6"/><xref ref-type="other" rid="par-7"/><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-4328"><name><surname>Andrews</surname><given-names>Steven S</given-names></name><xref ref-type="aff" rid="aff4"/><xref ref-type="other" rid="par-3"/><xref ref-type="other" rid="par-4"/><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-4329"><name><surname>Liu</surname><given-names>Danxu</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="other" rid="par-1"/><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-1312"><name><surname>Chalmers</surname><given-names>Ronald</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="other" rid="par-1"/><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/></contrib><aff id="aff1"><institution content-type="dept">School of Biomedical Sciences</institution>, <institution>University of Nottingham</institution>, <addr-line><named-content content-type="city">Nottingham</named-content></addr-line>, <country>United Kingdom</country></aff><aff id="aff2"><institution>The Babraham Institute</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United Kingdom</country></aff><aff id="aff3"><institution content-type="dept">Department of Biochemistry and Cambridge Systems Biology Centre</institution>, <institution>University of Cambridge</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United Kingdom</country></aff><aff id="aff4"><institution>Fred Hutchinson Cancer Research Center</institution>, <institution content-type="dept">Basic Sciences Division</institution>, <addr-line><named-content content-type="city">Seattle</named-content></addr-line>, <country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Ferrell</surname><given-names>James</given-names></name><role>Reviewing editor</role><aff><institution>Stanford University</institution>, <country>United States</country></aff></contrib></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>chalmers@nottingham.ac.uk</email></corresp><fn fn-type="con" id="equal-contrib"><label>†</label><p>These authors contributed equally to this work</p></fn><fn fn-type="present-address" id="pa1"><label>‡</label><p>Molecular Biology Program, Howard Hughes Medical Institute, Memorial Sloan-Kettering Cancer Center, New York, United States</p></fn></author-notes><pub-date date-type="pub" publication-format="electronic"><day>18</day><month>06</month><year>2013</year></pub-date><pub-date pub-type="collection"><year>2013</year></pub-date><volume>2</volume><elocation-id>e00668</elocation-id><history><date date-type="received"><day>27</day><month>02</month><year>2013</year></date><date date-type="accepted"><day>13</day><month>05</month><year>2013</year></date></history><permissions><copyright-statement>© 2013, Claeys Bouuaert et al</copyright-statement><copyright-year>2013</copyright-year><copyright-holder>Claeys Bouuaert et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/3.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife00668.pdf"/><abstract><object-id pub-id-type="doi">10.7554/eLife.00668.001</object-id><p>How do DNA transposons live in harmony with their hosts? Bacteria provide the only documented mechanisms for autoregulation, but these are incompatible with eukaryotic cell biology. Here we show that autoregulation of Hsmar1 operates during assembly of the transpososome and arises from the multimeric state of the transposase, mediated by a competition for binding sites. We explore the dynamics of a genomic invasion using a computer model, supported by in vitro and in vivo experiments, and show that amplification accelerates at first but then achieves a constant rate. The rate is proportional to the genome size and inversely proportional to transposase expression and its affinity for the transposon ends. Mariner transposons may therefore resist post-transcriptional silencing. Because regulation is an emergent property of the reaction it is resistant to selfish exploitation. The behavior of distantly related eukaryotic transposons is consistent with the same mechanism, which may therefore be widely applicable.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.001">http://dx.doi.org/10.7554/eLife.00668.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.00668.002</object-id><title>eLife digest</title><p>Transposons are regions of mobile DNA that can jump from one location in the genome to another. This represents a genetic burden to the host because there is always the risk that the transposon will inactivate a cellular gene. However, a greater problem is that transposition is accompanied by an increase in the number of copies of the transposon. Since each new copy will be a source of further new copies, amplification of transposons is necessarily exponential. The fact that eukaryotic cells are able to tolerate DNA transposons suggests the existence of regulatory mechanisms to defuse the inevitable genomic melt-down. Host-mediated epigenetic modifications and RNA interference will provide some level of protection. However, they are by no means completely effective and a well-adapted genomic parasite, such as a transposon, might be expected to have its own mechanism of regulation.</p><p>Now, Claeys Bouuaert, Lipkow and colleagues have used a computer model in combination with in vivo and in vitro experiments to search for this mechanism. Their experiments reveal how a DNA transposon is down-regulated by its own transposase. The transposase is the enzyme that catalyzes the ‘jump’ or transposition. It binds to specific sites at either end of the transposon and brings these together to make up a nucleoprotein complex called the transpososome. It is within this complex that the chemical steps of the reaction take place. When the number of transposons increases, so does the concentration of transposase. Claeys Bouuaert et al. show that the binding sites become saturated at a relatively low transposase concentration and that negative regulation arises from the resulting competition. Thus, the rate of transposition decreases as the number of transposons increases. They further use the computer model to explore how the amplification of the transposon is affected by transposon-specific and cellular-specific factors.</p><p>Claeys Bouuaert, Lipkow and colleagues based their study predominantly on a resurrected copy of the Hsmar1 transposon, which was active in the human genome 50 million years ago. However, they also tested two distantly related eukaryotic transposons and observed that their behavior was similar, which suggests that this could be a general mechanism that controls the activity of jumping genes. They also note that their competition mechanism is conceptually similar to the immunological ‘prozone effect’. This is a recurrent theme in protein chemistry and demonstrates once again that less is in fact sometimes more.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.002">http://dx.doi.org/10.7554/eLife.00668.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>DNA recombination</kwd><kwd>evolution</kwd><kwd>Host–parasite interactions</kwd><kwd>systems modeling</kwd><kwd>allostery</kwd><kwd>allosteric</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd>Human</kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution>Wellcome Trust</institution></institution-wrap></funding-source><award-id>093160/Z/10/Z</award-id><principal-award-recipient><name><surname>Claeys Bouuaert</surname><given-names>Corentin</given-names></name><name><surname>Liu</surname><given-names>Danxu</given-names></name><name><surname>Chalmers</surname><given-names>Ronald</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution>Royal Society</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Lipkow</surname><given-names>Karen</given-names></name></principal-award-recipient></award-group><award-group id="par-3"><funding-source><institution-wrap><institution>National Institute of General Medical Sciences</institution></institution-wrap></funding-source><award-id>91687</award-id><principal-award-recipient><name><surname>Andrews</surname><given-names>Steven S</given-names></name></principal-award-recipient></award-group><award-group id="par-4"><funding-source><institution-wrap><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>R01GM086615-01</award-id><principal-award-recipient><name><surname>Andrews</surname><given-names>Steven S</given-names></name></principal-award-recipient></award-group><award-group id="par-5"><funding-source><institution-wrap><institution>Microsoft Research Faculty Fellowship</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Lipkow</surname><given-names>Karen</given-names></name></principal-award-recipient></award-group><award-group id="par-6"><funding-source><institution-wrap><institution>Apple Research & Technology Support</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Lipkow</surname><given-names>Karen</given-names></name></principal-award-recipient></award-group><award-group id="par-7"><funding-source><institution-wrap><institution>EU-IndiaGrid2</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Lipkow</surname><given-names>Karen</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>A DNA transposon, or ‘jumping gene’, controls its amplification within a genome through a competition between the enzyme multimers that are responsible for its mobility.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>Transposons are discrete sections of mobile DNA that are amplified as they move from one location in the genome to another. They have had a profound influence on our evolutionary history by providing DNA duplications and rearrangements as substrates for natural selection (e.g. <xref ref-type="bibr" rid="bib22">Hua-Van et al., 2011</xref>). Although retrotransposons, which transpose via an RNA intermediate, can persist for long periods in a given genome, DNA transposons appear to rely on regular horizontal transfer (<xref ref-type="bibr" rid="bib32">Lohe et al., 1995</xref>; <xref ref-type="bibr" rid="bib19">Hartl et al., 1997</xref>). In the period following horizontal transfer the transpositional activity of the element is under positive selection. However, once there are multiple copies in the genome, selection is relaxed because a freely diffusing pool of transposase acts on all copies of the element. Gradually, with time, as various copies of the element acquire detrimental mutations, the pool of active transposase is poisoned by dominant-negative complementation leading to the extinction of the transposon (<xref ref-type="bibr" rid="bib32">Lohe et al., 1995</xref>; <xref ref-type="bibr" rid="bib19">Hartl et al., 1997</xref>; <xref ref-type="bibr" rid="bib28">Le Rouzic and Capy, 2005</xref>; <xref ref-type="bibr" rid="bib27">Le Rouzic et al., 2007</xref>).</p><p>Whilst this mechanism likely dictates the ultimate fate of a eukaryotic DNA-transposon, it remains unknown how control is exerted in the short- to medium-term after an active transposon first arrives in a genome. This is a serious problem because transposition is inherently exponential and each new copy of the element is a source of further new copies. The fact that long-lived multicellular organisms tolerate transposons suggests the existence of control mechanisms. Host-mediated epigenetic responses may provide some level of protection. However, these are probably never entirely effective, particularly for euchromatic copies of the transposon (e.g. <xref ref-type="bibr" rid="bib55">Kelleher and Barbash, 2013</xref>). Indeed, modeling suggests that unregulated transposition results either in the demise of the transposon or the demise of the host (<xref ref-type="bibr" rid="bib28">Le Rouzic and Capy, 2005</xref>). Experiments with a non-autonomous Mos1 element in <italic>Drosophila</italic> provided the first experimental evidence of autoregulation (<xref ref-type="bibr" rid="bib31">Lohe and Hartl, 1996</xref>). The phenomenon, termed overproduction inhibition (OPI), was revealed by a reduction in the frequency of excision when multiple copies of the transposase gene are present, or when transposase is over-expressed from a heat-shock promoter. OPI was later shown to affect other elements, both in vivo and in vitro (<xref ref-type="bibr" rid="bib26">Lampe et al., 1998</xref>; <xref ref-type="bibr" rid="bib10">Clark et al., 2009</xref>; <xref ref-type="bibr" rid="bib8">Claeys Bouuaert and Chalmers, 2010</xref>).</p><p>Our present work focuses on a resurrected copy of the Hsmar1 transposon, which was active in the human genome about 50 million years ago (<xref ref-type="bibr" rid="bib45">Robertson and Zumpano, 1997</xref>; <xref ref-type="bibr" rid="bib11">Cordaux et al., 2006</xref>; <xref ref-type="bibr" rid="bib30">Liu et al., 2007</xref>; <xref ref-type="bibr" rid="bib34">Miskey et al., 2007</xref>). It is closely related to Mos1 and is a member of the mariner family, which is probably the most successful group of transposons in nature, as judged by the depth and breadth of its phylogenetic distribution (e.g. <xref ref-type="bibr" rid="bib44">Robertson and Lampe, 1995</xref>; <xref ref-type="bibr" rid="bib15">Feschotte and Wessler, 2002</xref>). The mariner family is a good model system to address autoregulation because of the demands placed by horizontal transfer on the mechanism of transposition and its control. Ideally, the mechanism should be independent of specific host factors and should allow a rate of transposition sufficiently high to protect the founding element from genetic drift. High activity also provides a high probability of integration into a vector, such as a virus or endosymbiont, which may mediate horizontal transfer. However, this must be balanced against a detrimental effect on host fitness (e.g. <xref ref-type="bibr" rid="bib28">Le Rouzic and Capy, 2005</xref>).</p><p>Here we present a mechanism for autoregulation based on a biochemical analysis of Hsmar1 transposition. The model reveals how the kinetics of a genomic invasion are dictated by the multimeric state of the transposase and the order in which the various components are recruited into the developing transpososome. Once a certain number of copies of the transposon are established in the genome, and the transposase concentration is above a critical threshold, the mechanism provides a steady-state rate, perfectly damping the exponential amplification which is a natural consequence of unregulated transposition. The mechanism is an emergent property of the transposition reaction, based on the competition of active transposase multimers for their primary binding sites at the transposon ends. It is therefore robust and resistant to selfish exploitation. The key prediction of the model is that doubling the transposon copy number, which is tantamount to doubling the transposase concentration, halves the rate of transposition. In vivo transposase dose-response curves for Hsmar1 and the distantly related Sleeping Beauty (SB) and piggyBac (PB) transposons fit this condition, suggesting they are all regulated in the same way.</p></sec><sec id="s2" sec-type="results"><title>Results</title><sec id="s2-1"><title>The synapsis-by-protein-dimerization mechanism for transpososome assembly</title><p>The documented mechanisms for autoregulation in the bacterial transposons Tn10 and Tn5 provided the starting point for our investigation (<xref ref-type="bibr" rid="bib25">Kleckner, 1990</xref>; <xref ref-type="bibr" rid="bib38">Reznikoff, 2008</xref>). The respective transposases belong to the same family of nucleotidyl-transferases as the Hsmar1 transposase, and use the same cut-and-paste mechanism of transposition (e.g. <xref ref-type="bibr" rid="bib20">Hickman et al., 2010</xref>; see <xref ref-type="fig" rid="fig2">Figure 2A</xref> below for an illustration of cut-and-paste transposition). The cut-and-paste mechanism does not increase the number of transposons directly. Amplification requires that the empty donor site is repaired by homologous recombination from a sister chromosome, or that the transposon is excised behind a replication fork and integrates in front. However, unless stated otherwise we will assume a maximally efficient reaction in which each transposition event converts one copy of the element into two copies.</p><p>In Tn10 and Tn5, the first stage of the reaction is when transposase monomers (T) bind to the ends of the element (A and B) and bring them into a synapsis (<xref ref-type="fig" rid="fig1">Figure 1A</xref>). This complex is known as a paired ends complex (PEC) or transpososome. We will refer to this mechanism of transpososome assembly as ‘synapsis-by-protein-dimerization’ (S-PD). A kinetic diagram for the model is provided in <xref ref-type="fig" rid="fig1">Figure 1B</xref>. The first catalytic steps of the reaction are double strand breaks at the transposon ends, which are irreversible owing to the loss of enthalpy (subsumed into <italic>k</italic><sub>3</sub> in the diagram). Any regulation must occur before the formation of a productive synapsis because almost all transposons that achieve the first nick go on to complete cleavage and integration (<xref ref-type="bibr" rid="bib6">Chalmers and Kleckner, 1994</xref>, <xref ref-type="bibr" rid="bib7">1996</xref>).<fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.003</object-id><label>Figure 1.</label><caption><title>Dynamic models for genome invasions.</title><p>(<bold>A</bold>) The mechanism of transpososome assembly in Tn10/5. In the bacterial elements Tn10 and Tn5 synapsis is mediated by the dimerization of monomers bound to either end of the transposon. T, transposase; T<sub>2</sub>, transposase dimer; A and B, transposon ends; PEC, paired ends complex. (<bold>B</bold>) The kinetic model embedded in the computer simulation for the S-PD mechanism of synapsis. Abbreviations are as given in part (<bold>A</bold>); T<sub>free</sub> and T<sub>nsb</sub>, free and non-specifically bound transposase. (<bold>C</bold>) Simulation of the S-PD model. Kinetic parameters: <italic>k</italic><sub>0</sub>, 9.9 × 10<sup>6</sup> M<sup>−1</sup>s<sup>−1</sup>; <italic>k</italic><sub>−0</sub>, 139 s<sup>−1</sup>; <italic>k</italic><sub>1</sub>, 3.8 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup>; <italic>k</italic><sub>−1</sub>, 1.2 × 10<sup>−2</sup> s<sup>−1</sup>; <italic>k</italic><sub>2</sub>, 12.7 s<sup>−1</sup>; <italic>k</italic><sub>−2</sub>, 1 × 10<sup>−10</sup> s<sup>−1</sup>; <italic>k</italic><sub>3</sub>, 1.4 × 10<sup>−3</sup> s<sup>−1</sup>. Sources of the values are given in <xref ref-type="table" rid="tbl2">Table 2</xref>. (<bold>D</bold>) As in part (<bold>C</bold>) but the transposase is 99% cis-acting. 1% of the transposase leaks into the bulk phase and acts on all copies of the element. (<bold>E</bold>) As in part (<bold>C</bold>) but the transposon expresses a trans-acting inhibitor which reduces the transposase concentration. The inhibitor is 1000-fold more active than the transposase, which approximates the situation in Tn10 where the inhibitor is an antisense RNA. (<bold>F</bold>) The cis-acting transposase and the trans-acting inhibitor from parts (<bold>D</bold>) and (<bold>E</bold>) are combined. (<bold>G</bold>) The mechanism of transpososome assembly in mariner. A transposase dimer bound to one transposon end recruits a naked transposon end. Abbreviations as in part (<bold>A</bold>). (<bold>H</bold>) The assembly-site occlusion model. When the transposase concentration is low, most transposons are occupied by only one transposase dimer, leading to productive synapsis (bottom left). When the transposase concentration is high, a transposon may be occupied by two transposase dimers, causing overproduction inhibition (OPI, top right). SEC2 is a transposon with a transposase dimer bound to one of the ends. (<bold>I</bold>) Kinetic model embedded in the computer simulation for the ASO mechanism. Abbreviations as given in parts (<bold>A</bold>) and (<bold>B</bold>). (<bold>J</bold>) Simulation of the S-NEC mechanism and comparison with S-PD. Parameters as in part (<bold>C</bold>). (<bold>K</bold>) The time scale of the S-NEC simulation in part (<bold>J</bold>) is extended.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.003">http://dx.doi.org/10.7554/eLife.00668.003</ext-link></p></caption><graphic xlink:href="elife00668f001"/></fig></p><p>To explore the dynamics of a genomic invasion by a transposon that uses the S-PD mechanism for transpososome assembly we used the kinetic diagram in <xref ref-type="fig" rid="fig1">Figure 1B</xref> to develop a computer model (‘Materials and methods’). At the start of the simulation a single transposon appears in the genome and initiates an exponential amplification (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). Changing the values of the kinetic constants governing the model changes the scales on the axes but does not alter the exponential character of the reaction. Note that in this case the scale on the X axis is unrealistically short because we allow a productive synapsis to form products almost instantly (the value of <italic>k</italic><sub>3</sub> in <xref ref-type="fig" rid="fig1">Figure 1B</xref> is large). We are thus ignoring the time that would be required in vivo for the maturation of the integration products and repair of the donor site. This allows us to focus on the earlier stages of the reaction where regulation takes place.</p><p>Amplification of the Tn10 and Tn5 transposons, which use the S-PD mechanism, is controlled by a combination of two effects (<xref ref-type="bibr" rid="bib25">Kleckner, 1990</xref>; <xref ref-type="bibr" rid="bib38">Reznikoff, 2008</xref>). Firstly, the transposases do not diffuse freely and tend to act in cis to the encoding element. This slows the amplification considerably because it prevents the effective concentration in transposase from rising with transposon copy number. However, it does not alter the exponential nature of the curve (<xref ref-type="fig" rid="fig1">Figure 1D</xref>). Secondly, the transposons express a trans-acting inhibitor, which in the case of Tn10 is an antisense RNA directed against the transposase mRNA. This strategy is also ineffective in dampening the exponential amplification, even if the inhibitor is much more abundant than the transposase (<xref ref-type="fig" rid="fig1">Figure 1E</xref>). In contrast, the combination of a cis-acting transposase and a trans-acting inhibitor is very effective (<xref ref-type="fig" rid="fig1">Figure 1F</xref>). This combined strategy is not viable in a eukaryotic host because of the physical separation of transcription and translation. Transposons using the S-PD mechanism would experience exponential amplification in a eukaryote.</p></sec><sec id="s2-2"><title>The assembly-site occlusion model for eukaryotic transposons</title><p>Whilst the Tn10/5 transposases are monomers in solution, the Hsmar1 enzyme behaves as a dimer in gel filtration experiments (not shown). Furthermore, several mariner transposases readily form a complex known as single end complex 2 (SEC2), which is composed of a transposase dimer bound to a single transposon end (<xref ref-type="bibr" rid="bib29">Lipkow et al., 2004</xref>; <xref ref-type="bibr" rid="bib2">Auge-Gouillou et al., 2005</xref>; CCB & RC data not shown). We would like to propose that SEC2 gives rise to the transpososome by capturing a naked transposon end (<xref ref-type="fig" rid="fig1">Figure 1G</xref>). We will refer to this mechanism as synapsis-by-naked-end-capture (S-NEC). Recruitment of a naked transposon end immediately suggests a mechanism that can account for OPI (<xref ref-type="fig" rid="fig1">Figure 1H</xref>). As the concentration of transposase rises, transposition is inhibited by the progressive reduction in the number of unbound ends available for synapsis. We will refer to this as an assembly-site-occlusion model (ASO). A kinetic diagram for the model is provided in <xref ref-type="fig" rid="fig1">Figure 1I</xref>, following the same conventions as before, except that T<sub>2</sub> represents the transposase dimer.</p><p>To explore the dynamics of mariner transposition, we simulated the S-NEC mechanism of transpososome assembly in a second computer model (‘Materials and methods’). Using an identical set of parameters, S-NEC has an advantage over the S-PD mechanism in terms of the rate at early time points (<xref ref-type="fig" rid="fig1">Figure 1J</xref>). However, as the invasion progresses, the S-NEC model is dampened by ASO, which eventually dominates the reaction (<xref ref-type="fig" rid="fig1">Figure 1K</xref>). The peak rate of the reaction is when the free-transposase concentration is equal to its equilibrium binding constant.</p></sec><sec id="s2-3"><title>Verification of the ASO model</title><p>The ASO model provides testable predictions. The simplest is that the rate of transposition is governed by the ratio of transposase to transposon ends, and not by the absolute concentrations. We tested this using the in vitro excision assay (<xref ref-type="fig" rid="fig2">Figure 2A</xref>, but please also consult <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref> for a guide to interpretation of the in vitro reactions). Excision of the transposon leaves behind the vector backbone, which is an end product of the reaction and provides a direct measure of the efficiency. With 10 nM substrate, the reaction peaked at 25 nM transposase and then declined sharply, as judged by the amount of backbone produced. When the concentration of substrate was increased or decreased by 10-fold, the amount of transposase required for peak activity changed in direct proportion. This is consistent with the ASO model, and excludes models in which OPI arises from a concentration-dependent mechanism of transposase aggregation.<fig-group><fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.004</object-id><label>Figure 2.</label><caption><title>Experimental test of the ASO model for mariner transposition.</title><p>(<bold>A</bold>) The efficiency of transposition depends on the ratio of transposase to DNA. The cartoon to the left of the gels illustrates the mechanism of cut-and-paste transposition using a supercoiled (SC) substrate. First strand nicking generates an open circular product (Nick.). Second strand nicking at one end yields the linear single-end break product (Lin.). A similar set of nicks at the other transposon end yield the plasmid backbone (BB) plus the excised transposon. Transposase was titrated into reactions with three different substrate concentrations. Reactions contained 1.5 µg of supercoiled substrate plasmid in a volume of 500 µl, 50 µl and 5 µl, which provided a final concentration of 1 nM, 10 nM and 100 nM, respectively. Transposase was diluted so that the addition of one tenth of the respective reaction volumes achieved the indicated concentrations. Reactions were incubated for 4 hr at 37 °C and deproteinated. Photographs of ethidium bromide stained agarose gels are shown. Consumption of the supercoiled substrate and production of the plasmid backbone both provide a measure of the efficiency of the reaction. Inter., product of intermolecular integration of the transposon into an unreacted substrate. (<bold>B</bold>) Preassembly of the paired-end complex protects from OPI. Four sets of staged reactions were assembled as indicated. Set A was a standard transposition reaction with 6.7 nM plasmid substrate and 20 nM of transposase except that the components were pre-incubated for 3 hr before the addition of the catalytic Mg<sup>2+</sup> ion. Set A′ was identical except that Ca<sup>2+</sup> was present during the 3 hr pre-incubation period, which enhances PEC assembly. Sets B and B′ were identical to A and A′ except that the mixture was supplemented with 50 nM transposase 30 s before the addition of Mg<sup>2+</sup>. OPI inhibited the reaction in Set B. In contrast, the stable PEC assembly supported by the presence of Ca<sup>2+</sup> protected Set B′ from the inhibitory effects of the excess transposase added just before the catalytic Mg<sup>2+</sup> ion. Photographs of ethidium bromide stained agarose gels are shown. See also <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.004">http://dx.doi.org/10.7554/eLife.00668.004</ext-link></p></caption><graphic xlink:href="elife00668f002"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00668.005</object-id><label>Figure 2—figure supplement 1.</label><caption><title>Kinetic analyses of OPI and interpretation of in vitro transposition reactions.</title><p>The kinetics of the transposition reaction were analyzed in standard reactions containing 6.7 nM plasmid substrate and the indicated concentrations of the transposases. Photographs of ethidium bromide stained agarose gels are shown. (<bold>A</bold>) The effects of the ASO mechanism complicates the interpretation of the gels, particularly the kinetic analyses. In the reactions with 20 nM transposase about half of the substrate is converted to the nicked intermediate in the first 2 min. However, it takes an hour or more to consume the remainder of the substrate. This biphasic behavior arises from the ASO mechanism, which operates in any reaction that contains more than one dimer of transposase. Thus, even if the transposase was 100% active and the reaction was performed with the optimum ratio of one dimer to one transposon, only half of the substrate would react initially. At the start of the reaction half of the transposons would be occupied by a dimer and would react, a quarter would be occupied by two dimers and would suffer OPI and the other quarter would be completely unoccupied. The slow phase of consumption corresponds to the redistribution of dimers from transposons which initially suffered OPI to those that were completely unbound. Even under OPI conditions, when both ends of the transposon are occupied by transposase dimers, their occasional dissociation provides a window of opportunity for synapsis. OPI can therefore be overcome by extending the incubation period so that the accumulated windows of opportunity eventually suffice to complete the reaction. Thus, the reaction with 40 nM transposase goes on to reach completion during the 8 hr incubation, despite initially suffering from OPI. With 200 nM transposase, the extended incubation is unable to overcome OPI because the windows of opportunity provided by unoccupied transposon ends are shorter at this concentration. (<bold>B</bold>) The nicked substrate is more sensitive to OPI and is inhibited at only 40 nM transposase. This is because any free transposon ends, made available by the dissociation of a transposase dimer, are re-bound before they can achieve synapsis. In other words, the slow synapsis of the nicked substrate requires a longer window of opportunity than is available at 40 nM transposase. Time points and transposase concentrations as in part (<bold>A</bold>). (<bold>C</bold>) The RA104 transposase mutant is resistant to OPI because the higher OFF-rate provides more windows of opportunity for synapsis (see <xref ref-type="fig" rid="fig4">Figure 4C</xref> and text for details). Time points and transposase concentrations as in part (<bold>A</bold>).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.005">http://dx.doi.org/10.7554/eLife.00668.005</ext-link></p></caption><graphic xlink:href="elife00668fs001"/></fig></fig-group></p><p>Another prediction is that pre-assembly of the PEC will sequester the transposon ends in productive interactions, and protect the reaction from the subsequent addition of excess transposase. To test this hypothesis we measured the kinetics of four transposition reactions, which we assembled in different stages prior to the addition of the catalytic Mg<sup>2+</sup> ions at time zero (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). Set A is a standard transposition reaction in which the components were incubated for 3 hr prior to the addition of Mg<sup>2+</sup>. Set A′ is identical except that Ca<sup>2+</sup> was present during the pre-incubation period. Although Ca<sup>2+</sup> does not support DNA cleavage, it supports assembly of the PEC (<xref ref-type="bibr" rid="bib8">Claeys Bouuaert and Chalmers, 2010</xref>). In Set A′ this is evident from the faster consumption of the supercoiled substrate. Upon addition of Mg<sup>2+</sup>, 50% of the substrate was converted from the supercoiled to the nicked form within 1 min.</p><p>Set B is identical to Set A, except that extra transposase was added 30 s before the catalytic metal ion at time zero (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). As expected, this inhibited the reaction and very little substrate was converted into product. Set B′ is identical to Set B except that Ca<sup>2+</sup> was present during the pre-incubation period, before the addition of excess transposase. The reaction kinetics were then identical to those of Set A′. Pre-assembly of the transpososome is completely effective in the relief of OPI because the competition for free transposon ends is irrelevant at this stage of the reaction. Transposition is thus regulated during assembly of the transpososome, before productive synapsis of the ends.</p></sec><sec id="s2-4"><title>Synapsis is very sensitive to the rate of diffusion</title><p>We have assumed that each copy of the transposon contributes a certain amount of transposase to the nucleus and that this is where transposition takes place. We therefore estimated the rate of diffusion in this compartment and calculated its effects on the association rate constants for specific and non-specific DNA (‘Materials and methods’). Compared to the in vitro rates, the values were reduced by about twofold, but their effects cancel each other out and the rate of transposition is unchanged as a result. However, the effect of diffusion on the rate of synapsis is much greater because of the drag experienced by a chromosomally-bound protein. By treating the DNA with the worm-like chain model we obtained an estimate for ‘segmental-diffusion’, which describes the rate of motion on a short scale, such as the distance between the transposon ends (‘Materials and methods’). Our calculated value is consistent with the experimentally determined rate of looping derived from lox/Cre recombination in mammalian cells (<xref ref-type="bibr" rid="bib41">Ringrose et al., 1999</xref>). When the calculated diffusion rates are applied in the computer model, a pseudo-steady-state rate of transposition is established when there are about 200 copies of the transposon (<xref ref-type="fig" rid="fig3">Figure 3A</xref>, compare with <xref ref-type="fig" rid="fig1">Figure 1K</xref>). This steady-state rate persists well beyond the copy numbers typically achieved during a genomic invasion (<xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1</xref>).<fig-group><fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.006</object-id><label>Figure 3.</label><caption><title>A semi-quantitative description of mariner transposition.</title><p>(<bold>A</bold>) The ASO model for mariner transposition was simulated after taking account of the slow diffusion which prevails in vivo: <italic>k</italic><sub>0</sub>, 5 × 10<sup>6</sup> M<sup>−1</sup>s<sup>−1</sup>; <italic>k</italic><sub>1</sub>, 1.9 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup>; <italic>k</italic><sub>2</sub>, 4.5 × 10<sup>−4</sup> s<sup>−1</sup>. Other parameters were as in <xref ref-type="fig" rid="fig1">Figure 1K</xref>. (<bold>B</bold>) Transposase binds the transposon end rapidly and tightly. Binding reactions (20 µl) contained 40 fmol of <sup>32</sup>P-labeled transposon end and 120 fmol of transposase. The reactions were incubated at 37 °C and separated on a native polyacrylamide gel. The lane indicated as time zero contains no transposase. The 30 s time interval was the time required to mix the sample, load the gel and apply the voltage. SEC1 and SEC2 represent a single transposon end bound by a transposase monomer and dimer, respectively. An autoradiogram is shown. (<bold>C</bold>) Binding reactions were set up and analyzed as in part (<bold>B</bold>). After 5 min incubation, to allow the complexes to form, a 20-fold molar excess of unlabeled transposon end ('cold ITR') was added. The lane indicated as time zero contains transposase but no cold competitor. The rate of transposase dissociation can be estimated from the amount of free DNA released from the complexes. (<bold>D</bold>) The rate of transposon excision from a nicked substrate provides an estimate of the rate of synapsis. The kinetics of a transposition reaction was analyzed in standard reactions containing 6.7 nM of nicked plasmid substrate and 20 nM transposase. A photograph of an ethidium bromide stained agarose gels is shown. (<bold>E</bold>) The ASO model was simulated as in <xref ref-type="fig" rid="fig1">Figure 1K</xref> but taking account of the effects of allostery, as described in the main text and ‘Materials and methods’. Parameters as in <xref ref-type="fig" rid="fig1">Figure 1K</xref> except <italic>k</italic><sub>2</sub> = 9.6 × 10<sup>−5</sup> s<sup>−1</sup>; <italic>k</italic><sub>−1</sub> = 5.8 × 10<sup>4</sup> s<sup>−1</sup>. (<bold>F</bold>) As in part (<bold>E</bold>) but the rate of transposition is divided by the transposon copy number. (<bold>G</bold>) As in part (<bold>A</bold>), but taking account of allostery and slow diffusion in vivo: <italic>k</italic><sub>0</sub>, 5.0 × 10<sup>6</sup> M<sup>−1</sup>s<sup>−1</sup>; <italic>k</italic><sub>1</sub>, 1.9 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup>; <italic>k</italic><sub>2</sub>, 3.4 × 10<sup>−9</sup> s<sup>−1</sup>. The effects of changing transposition efficiency are shown (see text for definition of efficiency). (<bold>H</bold>) The relationship between transposition efficiency and the relative rate of transposition is plotted with the maximum rate scaled to 1: <italic>y</italic> = 1.4427 ln(<italic>x</italic>). See also <xref ref-type="fig" rid="fig3s1 fig3s2">Figure 3—figure supplements 1 and 2</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.006">http://dx.doi.org/10.7554/eLife.00668.006</ext-link></p></caption><graphic xlink:href="elife00668f003"/></fig><fig id="fig3s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00668.007</object-id><label>Figure 3—figure supplement 1.</label><caption><title>The time scale of the graph in <xref ref-type="fig" rid="fig3">Figure 3A</xref> is extended.</title><p>The parameters take account of the slow diffusion in vivo, but take no account of allosteric interactions between transposase subunits. This plot therefore represents the in vivo dynamics of a generic homomeric DNA looping reaction in the absence of allostery.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.007">http://dx.doi.org/10.7554/eLife.00668.007</ext-link></p></caption><graphic xlink:href="elife00668fs002"/></fig><fig id="fig3s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00668.008</object-id><label>Figure 3—figure supplement 2.</label><caption><title>EMSA analysis of transpososome assembly shows that SEC1 arises from dissociation of the PEC.</title><p>Transposon ends were encoded on radiolabeled linear DNA fragments. Binding reactions were incubated at 37 °C for 2 hr, separated on a 5% native polyacrylamide gel and recorded by phosphoimaging. (<bold>A</bold>) SEC1 and SEC2 represent a single transposon end bound by a transposase monomer and dimer, respectively (see main text for details). SEC2 comes to dominate the reaction as the transposase concentration rises. There is a significant transition between 4 and 8 nM transposase when SEC1 largely disappears. This corresponds to the point at which the transposon ends become sub-stoichiometric to the transposase dimers. At this point no free transposon ends remain as they are all bound by transposase. According to the S-NEC mechanism (see <xref ref-type="fig" rid="fig1">Figure 1</xref> for details), SEC2 is converted to the transpososome (=PEC) by recruitment of a naked transposon end. OPI occurs when the transposon ends are sub-stoichiometric to transposase dimers and there is a shortage of free transposon ends available for recruitment. Note that the various species observed in these binding reactions are identical to those observed in reactions with the related Mos1 and Himar1 transposons, which display the same behavior for example (<xref ref-type="bibr" rid="bib12">Dawson and Finnegan, 2003</xref>; <xref ref-type="bibr" rid="bib29">Lipkow et al., 2004</xref>). The present data suggests that the PEC in all three of these related systems is unstable in the EMSA and dissociates into two SEC1 complexes soon after the start of electrophoresis. Thus, in agreement with the S-NEC mechanism, SEC1 disappears at the point in the titration when the transposon ends become sub-stoichiometric to transposase dimers. (<bold>B</bold>) A fixed amount of transposase was titrated with an increasing amount of free transposon ends. The appearance of SEC1 coincides exactly with the appearance of the free transposon ends, which are required for PEC assembly in the S-NEC model. As the amount of transposon ends is increased further, the amount of SEC1 increases. This reflects mass action, which drives PEC assembly by favoring the capture of a free transposon end (see part D below for confirmation). This supports the data in part (<bold>A</bold>) in suggesting that SEC1 arises from the dissociation of the PEC. (<bold>C</bold>) Binding reactions were with a single-chain transposase dimer, in which two monomers are concatenated by a linker peptide joining the C-terminus of one to the N-terminus of the other. Concatenation of the subunits stabilizes the PEC, which is now detected in the gel. As the transposase concentration increases, the PEC disappears at the same point as the free DNA and gives way to SEC2. This behavior is identical to SEC1 in parts (<bold>A</bold>) and (<bold>B</bold>). The single-chain dimer of transposase is fully proficient for the transposition reaction (not shown), demonstrating that SEC1 is not an obligate intermediate of the reaction. This supports the data in parts (<bold>A</bold>) and (<bold>B</bold>) further confirming that SEC1 arises from the dissociation of the PEC. (<bold>D</bold>) In vitro transposition reactions were performed with a plasmid substrate encoding a single transposon end. Reactions were stopped at the indicated times and deproteinated before analysis on a 1.1% agarose TBE gel stained with ethidium bromide. All three sets of transposition reactions contained the same amounts of transposase and supercoiled plasmid substrate. However, the respective reaction volumes were adjusted to 500 µl, 50 µl and 5 µl to achieve the indicated concentrations. Transpososome assembly requires bimolecular synapsis between ends located on different molecules, as illustrated below the gels. This is inefficient owing to the relatively low concentration of one transposon end with respect to another when on separate molecules. When such a transpososome performs cleavage, followed by integration into an unreacted substrate molecule, the product is a linear molecule three times the size of the substrate (<xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>). There is very little reaction when the substrate concentration is low. This reflects the inefficiency of second end recruitment. At high substrate concentration, mass action drives the reaction by favoring second-end recruitment.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.008">http://dx.doi.org/10.7554/eLife.00668.008</ext-link></p></caption><graphic xlink:href="elife00668fs003"/></fig></fig-group></p></sec><sec id="s2-5"><title>Hsmar1-specific factors</title><p>Up until now we have assumed that synapsis of the transposon ends is the product of two sequential, chemically identical, collision events: a transposase dimer collides with and binds first one transposon end and then the other. This system is useful to illustrate the fundamental dynamics of a generic DNA-looping reaction, and how it is governed by the concentration of the looping protein. However, the Hsmar1 transpososome is more complicated because of communication between the transposase subunits. The allostery was first noted in competition experiments using different topological forms of the substrate (<xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>). The experiments revealed that transposase binds quickly and tightly to the first transposon end it encounters, but has a much lower affinity for the second transposon end (<xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>). To see how this affects the rate of transposition we estimated the association and dissociation rate constants, and the rate of synapsis, in in vitro experiments. To estimate the association rate constants (<italic>k</italic><sub>1</sub> in <xref ref-type="fig" rid="fig1">Figure 1I</xref>) we used an electrophoretic mobility shift assay (EMSA). Assembly of SEC2 was complete after 30 s, which was the minimum time required to load the gel and apply the voltage (<xref ref-type="fig" rid="fig3">Figure 3B</xref>). Although this did not provide a direct measure of <italic>k</italic><sub>1</sub>, it was consistent with the rapid site-specific binding of a helix-turn-helix protein and we therefore retained the previous best estimate from the literature (‘Materials and methods’).</p><p>To estimate the dissociation rate constant (<italic>k</italic><sub>−1</sub> in <xref ref-type="fig" rid="fig1">Figure 1I</xref>) we challenged SEC2 with a 10-fold molar excess of unlabeled transposon end (<xref ref-type="fig" rid="fig3">Figure 3C</xref>). The complex dissociated slowly and about 50% remained after 20 min. We therefore estimate that <italic>k</italic><sub>−1</sub> is 5.8 × 10<sup>−4</sup> s<sup>−1</sup>. This is some 20-fold lower than the value for most helix-turn-helix proteins (‘Materials and methods’). It is possible that the slow dissociation reflects a second aspect of the aforementioned allostery. In this experiment we also observed single-end-complex 1 (SEC1), which contains a monomer of transposase bound to a single transposon end. It arises from decay of the PEC, which is unstable in the EMSA (<xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2A–C</xref>). SEC1 is highly abundant in this experiment because the excess of cold competitor ends drives PEC assembly by mass action (<xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2D</xref>).</p><p>To measure the rate of synapsis (<italic>k</italic><sub>2</sub>) we used a transposon encoded on an open circular plasmid (<xref ref-type="fig" rid="fig3">Figure 3D</xref>), which is probably the most relevant substrate because eukaryotic DNA contains very little free supercoiling. The half-time for synapsis was estimated at 2 hr, which corresponds to a pseudo-first order rate constant of 9.6 × 10<sup>−5</sup> s<sup>−1</sup> (‘Materials and methods’).</p><p>When the experimental estimates of <italic>k</italic><sub>−1</sub> and <italic>k</italic><sub>2</sub> were applied in the computer model for the generic DNA looping reaction in vitro, the rate of transposition was greatly reduced and the accelerating phase at the beginning of the reaction was largely eliminated (<xref ref-type="fig" rid="fig3">Figure 3E</xref>, compare with <xref ref-type="fig" rid="fig1 fig3">Figures 1K and 3A</xref>). Although the overall rate of transposition remains constant, the rate for any particular copy of the element decays exponentially (<xref ref-type="fig" rid="fig3">Figure 3F</xref>). Finally, when we also take account of the slow diffusion in vivo, the rate of transposition is reduced even further (<xref ref-type="fig" rid="fig3">Figure 3G</xref>, solid line).</p></sec><sec id="s2-6"><title>The efficiency of transposition in vivo</title><p>Cut-and-paste transposition relies on the host homologous recombination machinery to reinstate a copy of the transposon at the donor site following excision. The highest rate of amplification, which we define as an efficiency of two, is achieved if, furthermore, the transposon excises after passage of a replication fork and inserts into an un-replicated region ahead of the fork. In this situation, one element on one chromosome becomes four elements on two sister chromosomes. The rate of transposition is quite sensitive to the efficiency of the reaction (<xref ref-type="fig" rid="fig3">Figure 3G,H</xref>). If the transposon inserts into a replicated region of the chromosome, or if the donor site is not reinstated, or if the excised transposon fails to reintegrate, the efficiency is reduced to 1.5 and the rate of transposition is halved. The rate of transposition may also be negative (<xref ref-type="fig" rid="fig3">Figure 3H</xref>).</p></sec><sec id="s2-7"><title>Association and dissociation rate constants</title><p>The computer model reveals an inverse linear relationship between the association rate constant (<italic>k</italic><sub>1</sub>) and the rate of transposition (<xref ref-type="fig" rid="fig4">Figure 4A</xref>). Increasing the affinity of the transposase for the transposon end reduces the rate of transposition, and vice versa. This is because the transposase acts primarily as an inhibitor of its own activity once the steady-state is established. In principle, this behavior could cause the preferential amplification of transposons which had acquired mutations in the transposase binding sites. However, when we reduced the association rate constant by 100-fold it was clear that the advantage was accompanied by a penalty in the form of a lag phase at the start of the invasion (<xref ref-type="fig" rid="fig4">Figure 4B</xref>). A low-affinity founding element would therefore be vulnerable to genetic drift during the period before sufficient copies had been produced to sustain the transposase concentration required for efficient amplification.<fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.009</object-id><label>Figure 4.</label><caption><title>Association and dissociation rate constants.</title><p>(<bold>A</bold>) Simulation as in <xref ref-type="fig" rid="fig3">Figure 3E</xref>, with <italic>k</italic><sub>1</sub> fivefold up or down. We have retained the effects of allostery but here and in subsequent simulations we have ignored the effects of the slow diffusion in vivo. This allows the algorithm to run more smoothly and shortens the scale on the x axis, but does not affect the conclusions, which are based on the differential responses to changing the various parameters. (<bold>B</bold>) Simulation as in part (<bold>A</bold>) (solid line), with <italic>k</italic><sub>1</sub> 100-fold down. (<bold>C</bold>) Binding reactions with the RA104 mutant transposase were set up as in <xref ref-type="fig" rid="fig3">Figure 3B</xref>. In the lane with 20 nM transposase the smear between SEC2 and the position of free DNA is probably due to complexes that have dissociated during electrophoresis. Autoradiogram of an EMSA is shown. (<bold>D</bold>) The kinetics of the transposition reaction were analyzed in standard reactions containing 6.7 nM of supercoiled plasmid substrate and the indicated concentrations of the transposases. Photographs of ethidium bromide stained agarose gels are shown. With 200 nM wild-type transposase, the windows of opportunity for synapsis, provided by periods when one transposon end is unoccupied, are too short to allow for synapsis. The RA104 transposase mutant is resistant to OPI because the higher dissociating rate provides more windows of opportunity for synapsis. (<bold>E</bold>) Mutant and wild type transposase were assayed in HeLa cell culture with 8 ng or 1000 ng of helper plasmid and 500 ng of neomycin resistant reporter plasmid.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.009">http://dx.doi.org/10.7554/eLife.00668.009</ext-link></p></caption><graphic xlink:href="elife00668f004"/></fig></p><p>To verify this property of the model we mutagenized the transposase to lower its affinity for the transposon end. The RA104 mutation is located within the second helix-turn-helix motif, and abolishes an ionic interaction with the phosphodiester backbone (<xref ref-type="bibr" rid="bib46">Roman et al., 2007</xref>; <xref ref-type="bibr" rid="bib40">Richardson et al., 2009</xref>). In the EMSA the RA104 transposase produced very little SEC2, even when the standard transposase concentration was increased threefold (<xref ref-type="fig" rid="fig4">Figure 4C</xref>, compare with <xref ref-type="fig" rid="fig3">Figure 3B</xref>). However, in the in vitro plasmid transposition assay, RA104 consumed the supercoiled substrate faster than wild type (<xref ref-type="fig" rid="fig4">Figure 4D</xref>). Presumably, the mutant protein redistributes more quickly from those plasmids initially suffering OPI due to the double occupancy of their transposon ends. When the transposase concentration was increased 10-fold the wild type activity was almost abolished, whereas the mutant was almost unaffected (<xref ref-type="fig" rid="fig4">Figure 4D</xref>).</p><p>We next turned to an established eukaryotic cell culture assay. It is based on the co-transfection of a helper plasmid, expressing transposase, and a reporter plasmid, encoding a transposon that confers neomycin resistance when it integrates into a host chromosome (e.g. [<xref ref-type="bibr" rid="bib18">Grabundzija et al., 2010</xref>]). The rate of transposition is given by the number of stable transfectants obtained after drug selection. With the optimum amount of helper plasmid RA104 was only 25% as active as wild type. However, with higher levels RA104 was less severely inhibited (<xref ref-type="fig" rid="fig4">Figure 4E</xref>). Thus, the RA104 transposase reflects the behavior of the low affinity element in <xref ref-type="fig" rid="fig4">Figure 4B</xref>, which is at a disadvantage to wild type when the transposase concentration is low.</p></sec><sec id="s2-8"><title>Transposase expression level</title><p>Once a transposon is established in the genome, ASO provides a constant rate of transposition. With a given set of kinetic parameters for transposase binding and synapsis, the actual rate achieved depends on the amount of transposase expressed by each copy of the transposon (<xref ref-type="fig" rid="fig5">Figure 5A</xref>). Reducing the expression level by fivefold increases the steady-state rate by the same factor. However, it takes slightly longer to establish the final rate because of the difficulty locating the transposon ends at the start of the invasion when the transposase concentration is low. The limits of this effect are reached when the expression of transposase dimers is fractionally sub-stoichiometric to transposon ends (<xref ref-type="fig" rid="fig5">Figure 5B</xref>). Since ASO requires at least two dimers per transposon, the accelerating phase is never dampened and amplification, although slow at first, is exponential.<fig-group><fig id="fig5" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.010</object-id><label>Figure 5.</label><caption><title>Transposase expression level.</title><p>(<bold>A</bold>) Simulation as in <xref ref-type="fig" rid="fig3">Figure 3E</xref> showing the effects of changing the transposase expression fivefold up or down. Unless noted otherwise, each transposon copy produces 500 transposase dimers, which are considered as being contained within a 500 fl nucleus. (<bold>B</bold>) As in part (<bold>A</bold>) (1× line) but with a transposase expression level of 1.9 dimers per transposon. (<bold>C</bold>) Parameters as part (<bold>A</bold>) (1× line), plotting the relationship between transposase expression and the rate of transposition at the point in the invasion when there are 1000 copies of the transposon present in the genome. (<bold>D</bold>) HeLa cell assay for Hsmar1 transposition. HeLa cells were transfected with 500 ng of neomycin reporter plasmid and the indicated amount of helper plasmid expressing transposase. The rate of transposition is given by the number of neomycin resistant colonies recovered after drug selection. Bars indicate standard error of the mean, n = 3. R<sup>2</sup> is least squares goodness of fit to the line <italic>y = ax</italic><sup><italic>−1</italic></sup>. (<bold>E</bold>) As in part (<bold>D</bold>) but with isogenic Sleeping Beauty (SB100X) transposon reporter and helper plasmids. Data points are a mean of three experiments and were extracted and re-plotted from <xref ref-type="fig" rid="fig2">Figure 2A</xref> of (<xref ref-type="bibr" rid="bib18">Grabundzija et al., 2010</xref>). (<bold>F</bold>) As in part (<bold>D</bold>) but with isogenic piggyBac transposon reporter and helper plasmids. n = 3. See also <xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.010">http://dx.doi.org/10.7554/eLife.00668.010</ext-link></p></caption><graphic xlink:href="elife00668f005"/></fig><fig id="fig5s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00668.011</object-id><label>Figure 5—figure supplement 1.</label><caption><title>Transposase expression in HeLa cells.</title><p>HeLa cells were transfected with an increasing quantity of the helper plasmid, which expresses Hsmar1 transposase from a strong viral promoter. The titration was identical to that shown in <xref ref-type="fig" rid="fig5">Figure 5D</xref>. Cells were harvested after two days and analyzed by SDS-PAGE and Western blotting with antibodies against β-tubulin and the Hsmar1 transposase. Transposase expression increased throughout the titration. This indicates that transposase expression is not toxic. If transposase expression had been toxic the amount of transposase recovered from the cells would have plateaued or decreased with higher amounts of helper plasmid transfected.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.011">http://dx.doi.org/10.7554/eLife.00668.011</ext-link></p></caption><graphic xlink:href="elife00668fs004"/></fig></fig-group></p><p>We next considered a genome with 1000 copies of the transposon and calculated the rate of transposition over a range of transposase expression levels (<xref ref-type="fig" rid="fig5">Figure 5C</xref>). This revealed an inverse-exponential relationship between transposase expression and the rate of transposition. This is noteworthy because it suggests that once the steady-state is established, a post-transcriptional silencing response may potentiate the genomic invasion. It also suggests that titration of the transposase by non-autonomous elements may not contribute to the demise of a transposon as has been previously suggested.</p><p>To verify the inverse relationship between the transposase concentration and the rate of transposition, we titrated the helper plasmid in the HeLa cell culture assay (<xref ref-type="fig" rid="fig5">Figure 5D</xref>). The inhibitory response followed the power law relationship predicted from the ASO mechanism (<italic>y</italic> = <italic>ax</italic><sup>−1</sup>). A Western blot was used to confirm that transposase expression increased with increasing amounts of helper plasmid (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1</xref>). We also examined the transposase dose-response relationships for two other eukaryotic DNA transposons (<xref ref-type="fig" rid="fig5">Figure 5E,F</xref>). Sleeping Beauty (SB) is distantly related to Hsmar1 within the mariner/Tc1 superfamily (<xref ref-type="bibr" rid="bib17">Goodier and Davidson, 1994</xref>; <xref ref-type="bibr" rid="bib23">Ivics et al., 1997</xref>). The lepidopteran transposon piggyBac is very distant from Hsmar1 (<xref ref-type="bibr" rid="bib16">Fraser et al., 1996</xref>). In both cases the transposase dose-response curve matched the predicted power law relationship.</p></sec><sec id="s2-9"><title>Genome size and nuclear volume</title><p>Mariner transposons are widely distributed in plants and animals where genome sizes and nuclear volumes vary greatly. In the S-NEC computer model, the rate of transposition is barely sensitive to volume changes less than 100-fold (<xref ref-type="fig" rid="fig6">Figure 6A</xref>). This is because changes in the amount of transposase bound to specific and non-specific sites almost exactly cancel each other. The S-PD mechanism of synapsis responds in the opposite direction (<xref ref-type="fig" rid="fig6">Figure 6B</xref>).<fig id="fig6" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.012</object-id><label>Figure 6.</label><caption><title>Genome size and nuclear volume.</title><p>(<bold>A</bold>) Simulation as in <xref ref-type="fig" rid="fig3">Figure 3E</xref>, changing the nuclear volume by the indicated amounts. (<bold>B</bold>) Simulation of regulated S-PD mechanism as in <xref ref-type="fig" rid="fig1">Figure 1F</xref>, changing the nuclear volume by the indicated amounts. (<bold>C</bold>) Simulation as in part (<bold>A</bold>), changing the genome size by the indicated amounts. (<bold>D</bold>) Simulation as in part (<bold>B</bold>), changing the genome size by the indicated amounts.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.012">http://dx.doi.org/10.7554/eLife.00668.012</ext-link></p></caption><graphic xlink:href="elife00668f006"/></fig></p><p>Increasing the size of the genome in the S-NEC/ASO model increases the rate of transposition in direct proportion (<xref ref-type="fig" rid="fig6">Figure 6C</xref>). This is a result of the inverse relationship between the transposase concentration and the rate of transposition, which was verified in <xref ref-type="fig" rid="fig5">Figure 5D</xref>. The effect arises from the facts that the number of non-specific binding sites scales with the genomes size, and that these relieve OPI by absorbing free transposase. Whether or not non-specific sites are composed of DNA or chromatin does not matter. The ASO mechanism thus tailors the rate of transposition according to the genome size. The S-PD model responds in the opposite direction, and smaller genomes receive progressively higher rates of transposition (<xref ref-type="fig" rid="fig6">Figure 6D</xref>).</p></sec><sec id="s2-10"><title>Variations and alternative models</title><p>The rice transposon Osmar14 has a secondary transposase binding site close to one of the two principal binding sites at the transposon ends (<xref ref-type="bibr" rid="bib54">Yang et al., 2009</xref>). It was suggested that transposase binding at the secondary site may down-regulate transposition. To investigate the behavior of such a system, we introduced an inhibitory binding site into the S-NEC reaction scheme (<xref ref-type="fig" rid="fig7">Figure 7A</xref>). The steady-state rate of transposition was reduced in proportion to the binding affinity of the transposase. Note that the behavior of the system would be identical if the transposon instead expressed a transposase variant or a completely different protein, provided that binding interfered with transposition.<fig id="fig7" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.013</object-id><label>Figure 7.</label><caption><title>Variations and alternative mechanism of regulation.</title><p>(<bold>A</bold>) The S-NEC mechanism is simulated as in <xref ref-type="fig" rid="fig3">Figure 3E</xref> but with a secondary transposase binding site that inhibits transposition when occupied. 0, the secondary site has no affinity for the transposase; 1×, the secondary site has the same affinity for transposase as the primary binding sites at the transposon ends; 2×, secondary site with twice the affinity for transposase. (<bold>B</bold>) The dimerization end occlusion model for OPI. (<bold>C</bold>) All possible states of the DEO model are illustrated. A and B are transposon ends, T is an active transposase monomer, T<sub>2</sub> is an inhibitory transposase dimer. Solid lines indicate reaction pathways allowed in the model. The active species is shown in red. (<bold>D</bold>) Simulation of the DEO model using in vitro, non-allosteric parameters as in <xref ref-type="fig" rid="fig1">Figure 1C</xref>. Inhibitory dimers bind transposon ends with the same affinity as the active monomers. Active monomers dimerize with the same affinity as they bind transposon ends. (<bold>E</bold>) DEO model as in part (<bold>D</bold>) but with in vivo slow diffusion rate as in <xref ref-type="fig" rid="fig3">Figure 3E</xref>. Affinity of the monomers in the inhibitory dimer is doubled as indicated. (<bold>F</bold>) DEO model as in part (<bold>E</bold>) showing the effect of doubling in the genome size and nuclear volume. (<bold>G</bold>) DEO model as in part (<bold>E</bold>) showing the effect of changing the association rate constant of the active monomers and the inhibitory dimers by the indicated amount.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.013">http://dx.doi.org/10.7554/eLife.00668.013</ext-link></p></caption><graphic xlink:href="elife00668f007"/></fig></p><p>In the context of the S-PD reaction scheme, a secondary transposase binding site is identical to the trans-acting inhibitor modeled in <xref ref-type="fig" rid="fig1">Figure 1E</xref>. Thus, even if the inhibitory binding site is stronger than the principal binding sites, it will not halt the exponential amplification, but merely slow it down. This highlights the weakness of ‘secondary site’ models for regulation in the absence of the dose-dependent inhibition provided by the ASO mechanism.</p><p>ASO is reminiscent of the dimerization end-occlusion model (DEO) first proposed for Tn5 and later in more detail for the Mos1 mariner transposon (<xref ref-type="bibr" rid="bib53">Weinreich et al., 1994</xref>; <xref ref-type="bibr" rid="bib51">Townsend and Hartl, 2000</xref>). The model postulates an S-PD mechanism of transpososome assembly, accompanied by the concentration-dependent formation of inactive dimers, which bind and occlude the transposon ends (<xref ref-type="fig" rid="fig7">Figure 7B</xref>). Tn5 was later shown to be regulated by other factors (above), and the current biochemical analysis excludes the DEO mechanism for mariner transposition. Nevertheless, DEO is a viable model and we wished to explore its properties: firstly to compare them with ASO, and, secondly, to facilitate its recognition should it be present in any natural element that may be characterized in the future.</p><p>All potential states of the DEO model are illustrated in <xref ref-type="fig" rid="fig7">Figure 7C</xref>. However, in constructing a computer model we used a slightly simplified version in which only free monomers produce inactive dimers. In this scheme there is a single potentially active species, shown in red. As previously demonstrated, the model provides the dose-dependent inhibition needed to prevent exponential amplification (<xref ref-type="bibr" rid="bib51">Townsend and Hartl, 2000</xref>; <xref ref-type="fig" rid="fig7">Figure 7D,E</xref>). The eventual steady-state rate is inversely proportional to the stability of the inhibitory dimer (<xref ref-type="fig" rid="fig7">Figure 7E</xref>). As in the ASO model, the steady-state rate is also proportional to the genome size (<xref ref-type="fig" rid="fig7">Figure 7F</xref>). However, larger genomes experience a longer lag phase at the start of the invasion. This is an inherent disadvantage of the S-PD mechanism of transpososome assembly and reflects the improbability of transposase monomers binding to both ends of the element simultaneously when their concentration is low (<xref ref-type="fig" rid="fig1">Figure 1J</xref>). Interestingly, once established, the steady state rate is independent of the association and dissociation rate constants, <italic>k</italic><sub>1</sub> and <italic>k</italic><sub>−1</sub>, provided that dimerization does not change either of these values (<xref ref-type="fig" rid="fig7">Figure 7G</xref> and not shown).</p></sec></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>DNA transposons require regulation to balance their fitness at different stages of a genomic invasion (<xref ref-type="bibr" rid="bib28">Le Rouzic and Capy, 2005</xref>). Initially, a high rate of transposition is desirable as it guards against genetic drift. Later, as the copy number rises, a progressively lower activity per element would prevent excessive damage to the host (e.g. <xref ref-type="bibr" rid="bib51">Townsend and Hartl, 2000</xref>; <xref ref-type="bibr" rid="bib28">Le Rouzic and Capy, 2005</xref>). In the ASO model for mariner transposition, the rate of the reaction increases faster than the copy number at the start of a genomic invasion (<xref ref-type="fig" rid="fig3">Figure 3A,G,E</xref>). However, later when a number of copies exist in the genome, and the transposase concentration is above a certain threshold, the increased transpositional activity provided by each new copy is almost exactly balanced by the increased competition for binding sites at the transposon ends (<xref ref-type="fig" rid="fig3">Figure 3E,F</xref>). We note that the competition is conceptually similar to the immunological prozone phenomenon, which has been predicted to apply widely in biological systems that involve bridging interactions (<xref ref-type="bibr" rid="bib4">Bray and Lay, 1997</xref>).</p><sec id="s3-1"><title>The cell cycle, and the soma vs the germ line</title><p>Transposition in somatic cells leaves no trace in the genomic fossil record and does not help the vertical spread of the transposon through the population. It is also detrimental to the host, particularly in organisms where there is a risk of cancer. Furthermore, since somatic cells may spend a long time in the G1 phase of the cell cycle, when double strand breaks are processed by non-homologous end joining (NHEJ), transposition may cause a net decrease in the number of transposon copies in the event of a failure to reintegrate (<xref ref-type="fig" rid="fig3">Figure 3H</xref>). The only selective advantage of transposition in soma that we can envisage is that it provides an opportunity for integration into a vector, which may mediate horizontal transfer. Indeed, mariner transposes readily in the soma and is among the most widespread of elements in nature. In contrast, the P element has a limited phylogenetic distribution and its activity is restricted to the germ line (<xref ref-type="bibr" rid="bib14">Engels, 1983</xref>).</p><p>Transposition in the germ line aids the vertical spread of the transposon in the population and leaves a genomic fossil record of some of the events, which we still observe today. Cut-and-paste transposition is favored in the germ line because homologous recombination, which is required to reinstate the donor site, is more prevalent than in the soma throughout all phases of the cell cycle (<xref ref-type="bibr" rid="bib42">Robert et al., 2008</xref>; <xref ref-type="bibr" rid="bib50">Tichy et al., 2010</xref>; <xref ref-type="bibr" rid="bib48">Serrano et al., 2011</xref>). Synchronization of transposition with S phase of the cell cycle would provide a further advantage due to the general up-regulation of homologous recombination and the chance of integrating ahead of a replication fork. This strategy was first observed in the bacterial IS<italic>10</italic> (<xref ref-type="bibr" rid="bib43">Roberts et al., 1985</xref>). P element transposition may also be linked to the cell cycle because it preferentially integrates near origins of replication, which will ensure its early replication (<xref ref-type="bibr" rid="bib49">Spradling et al., 2011</xref>). A counter example is that expression of SB transposase is reported to extend G1 and promote the repair of excision sites by NHEJ in the soma (<xref ref-type="bibr" rid="bib52">Walisko et al., 2006</xref>). It is unknown whether mariner transposition is synchronized with the cell cycle. However, this seems plausible because Hsmar1 transposition is greatly stimulated by negative supercoiling (<xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>), which probably exists transiently behind an advancing replication fork.</p></sec><sec id="s3-2"><title>The rate of synapsis is a dominant factor</title><p>The values of the kinetic parameters determine the magnitude of the steady-state rate of transposition as well as the copy number at which it is achieved. The rate of synapsis is a dominant factor, and is in turn influenced by the multimeric state of the transposase, the rate of diffusion and, in the case of Hsmar1, by the allostery between the subunits. At the start of an invasion, a preformed transposase multimer overcomes the improbability of two monomers simultaneously binding opposite ends of the same element when their concentration is low. The advantage arises from the continuity of the DNA between the transposon ends, which ensures that the first dimer-bound-site has a high concentration with respect to its partner. The molar concentration of one DNA site with respect to another is abbreviated as <italic>j</italic><sub><italic>M</italic></sub>, and is about 55 nM for the ends of a 1.3 kb transposon, such as Hsmar1. The magnitude of this value dictates the effectiveness of the ASO mechanism. In dilute solutions, where synapsis is fast, OPI is not significant until the free transposase concentration is a significant fraction of <italic>j</italic><sub><italic>M</italic></sub>. This is evident in <xref ref-type="fig" rid="fig1">Figure 1K</xref> where 50% inhibition is not achieved until there are 10<sup>7</sup> copies of the transposon, which corresponds to 40 haploid human-genomes of DNA. However, under the slow diffusion regime that prevails in vivo, OPI becomes significant much earlier, and the steady state is achieved when there are about 200 copies of the transposon (<xref ref-type="fig" rid="fig3">Figure 3A</xref>). Allostery further increases the effectiveness of ASO by reducing the affinity of the developing transpososome for the second transposon end. This increases the competition for binding sites to such an extent that the steady-state rate is established when there are only a few copies of the transposon (<xref ref-type="fig" rid="fig3">Figure 3G</xref>).</p></sec><sec id="s3-3"><title>Evolutionary considerations</title><p>The influence of <italic>j</italic><sub><italic>M</italic></sub> on the rate of transposition helps to explain the success of miniaturized transposons, which often greatly outnumber their cognate master elements (e.g. <xref ref-type="bibr" rid="bib54">Yang et al., 2009</xref>). However, these relationships have additional complexities. Osmar14 is an autonomous transposon in the rice genome and provides the cognate transposase for Ost35, a more numerous miniature element (<xref ref-type="bibr" rid="bib54">Yang et al., 2009</xref>). This relationship is not obvious from a comparison of the respective transposon end sequences, which appear quite divergent. Our finding that OPI is relieved by reducing the affinity between the transposase and the transposon end provides a further explanation for the unexpected efficiency of this cross-mobilization in addition to their short length (<xref ref-type="fig" rid="fig4">Figure 4A</xref>).</p><p>Another factor that may explain why Osmar14 is less numerous than Ost35 is that it has a secondary transposase binding site near its 3′-end, which inhibits transposition and may represent a mechanism of autoregulation (<xref ref-type="bibr" rid="bib54">Yang et al., 2009</xref>). However, as noted already, a repressive binding site such as this does not change the underlying dynamics of ASO (<xref ref-type="fig" rid="fig7">Figure 7A</xref>), nor does it dampen exponential amplification in the absence of ASO or a similarly effective mechanism (<xref ref-type="fig" rid="fig1">Figure 1E</xref>). Another type of secondary-site model is one in which transposase represses its own expression. This is conceptually similar to the cis-action observed in some bacterial systems, in that it prevents the concentration of transposase rising in proportion to the transposon copy number. The weakness of such models is that they only prevent the transposase concentration from rising and can never cause a reduction (<xref ref-type="bibr" rid="bib51">Townsend and Hartl, 2000</xref>). For example, expression of the P element transposase in the soma is repressed in part by binding of a truncated splice variant. However, co-expression of the truncated transposase with the full-length version in the germ line is not sufficient to establish the repressive P-cytotype, which probably requires piwi RNA arising from specific telomeric copies of the element (<xref ref-type="bibr" rid="bib35">Misra and Rio, 1990</xref>; <xref ref-type="bibr" rid="bib5">Brennecke et al., 2008</xref>; <xref ref-type="bibr" rid="bib24">Jensen et al., 2008</xref>).</p><p>A more general argument against the effectiveness of ‘secondary-site’ models is that they are vulnerable to selfish exploitation by any element that lacks the site. Thus, the burden of Ost35 transposition, which lacks the secondary site, would presumably negate any advantage accruing to Osmar14 as a result of its self restraint. ASO is less vulnerable than the secondary-site models because it relies on a competition between active multimers of the recombinase for their primary binding sites at the transposon ends. Although it can, in principle, be subverted by the resistance of low-affinity elements to OPI, this effect can only operate within limits (<xref ref-type="fig" rid="fig4">Figure 4B</xref>, note the slow amplification when the copy number is low). Although not completely immune, the ASO mechanism is robust to this type of selfish exploitation. ASO also provides a homeostatic mechanism in that a doubling of the transposase concentration always halves the rate of transposition per element, irrespective of the values of the kinetic constants governing the reaction. A larger genome will thus always receive a proportionately higher rate of transposition. In the case of mariner, the relative genetic burden is thus independent of the host genome size and is dictated by adaptive features of the transposon itself. These features include the transposase binding affinity, the strength of the promoter in a given host and the strength of the allosteric coupling between subunits. The latter is a powerful determinant able to change the rate of synapsis over many orders of magnitude. Indeed, allostery provides significant adaptive flexibility, which is unavailable in the S-PD reaction mechanism. Finally, by buffering the rate of transposition against changes in the transposase concentration, ASO may also provide a degree of protection from post-transcriptional silencing and heterochromatinization responses.</p></sec></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><sec id="s4-1"><title>Protein expression, purification and in vitro assays</title><p>The wild type and RA104 transposases were expressed from pRC880 and pRC1128, respectively (<xref ref-type="bibr" rid="bib8">Claeys Bouuaert and Chalmers, 2010</xref>). Unless stated otherwise, each 50 µl transposition reaction contained 6.7 nM of plasmid substrate and 20 nM transposase in 20 mM Tris-HCl pH8, 100 mM NaCl, 2 mM DTT, 2.5 mM MgCl<sub>2</sub> and 10% glycerol. The standard substrate was pRC650, which encodes a mini-transposon with 30 bp Hsmar1 transposon ends. Plasmid pRC919 has a single 30 bp Hsmar1 end cloned into the pBluescript polylinker. Open circular substrates were prepared using the Nb.BsrDI endonuclease, which nicks the plasmids at several sites some distance away from the transposon ends. Transposition reactions were incubated at 37 °C and analyzed by loading 400 ng of the DNA on each lane of a TBE-buffered 1.1% agarose gel. After electrophoresis, the gel was stained with ethidium bromide, destained in water, and photographed on a 310 nm transilluminator using a DC290 camera (Kodak, Rochester, NY) with a 590 DF bandpass filter. Digital photographs were quantified using the Image Gauge software package ScienceLab 2003 (Fuji Corporation, Tokyo, Japan).</p><p>In the EMSA 97 bp DNA fragments carrying Hsmar1 transposon ends were prepared by digesting pRC919 with XmaI and labeled at both ends using [α-<sup>32</sup>P]dCTP and the Klenow enzyme. Unless stated otherwise each 20 μl reaction contained 250 ng of non-specific plasmid DNA as a carrier, 2 nM labeled substrate and 6 nM transposase. Complexes were assembled for the indicated times in a buffer containing 20 mM HEPES pH7.5, 100 mM NaCl, 2 mM DTT, 10% glycerol, 250 μg/ml BSA. Products were separated on 5 or 7% polyacrylamide Tris-acetate-EDTA gels.</p></sec><sec id="s4-2"><title>In vivo assays</title><p>Transposition assays in HeLa cells were performed as described by (<xref ref-type="bibr" rid="bib34">Miskey et al., 2007</xref>; <xref ref-type="bibr" rid="bib18">Grabundzija et al., 2010</xref>) using isogenic plasmids, except that the SB transposase gene and transposon ends were replaced by the respective sequences from PB and Hsmar1.</p></sec><sec id="s4-3"><title>Development and implementation of the computer model</title><p>The kinetic models for interactions between transposase and DNA is presented in <xref ref-type="fig" rid="fig1">Figure 1B,I</xref>. A and B represent the identical inverted repeat sequences flanking the transposon, which are the specific binding sites of the transposase. The hyphens joining A and B represent the transposon DNA between the inverted repeats. T<sub>2</sub> represents transposase, which is always dimeric in this model, T<sub>2,nsb</sub> represents non-specifically bound transposase, T<sub>2,free</sub> represents freely diffusing transposase, AT and BT represent transposase monomers bound to A or B, respectively, and AT<sub>2</sub> and BT<sub>2</sub> represent transposase dimers bound to A or B, respectively. The flux through the reaction was simulated in the Macintosh version of MATLAB R2010b using the ordinary differential equation solver ode15s.</p><p>For the mariner model, we begin by considering an idealized situation in which we make two significant assumptions. Firstly, that since the A and B sites are identical, they have identical binding kinetics. Secondly, that the two DNA-binding sites on a single transposase dimer have identical behavior and do not interact with each other, meaning that there are no allosteric interactions between the transposase subunits. In this idealized situation the first and second transposon end to be bound by the transposase, which ever this may be, are bound with equal affinity and the rate depends only on the relative concentration of the reactants. Later we will also consider a model in which allosteric interactions between the transposase subunits causes a kinetic barrier to recruitment of the second transposon end.</p><p>Because of these assumptions mentioned above, our model only needs to include the reaction rate constants shown in <xref ref-type="fig" rid="fig1">Figure 1I</xref>, which are <italic>k</italic><sub>0</sub>, <italic>k</italic><sub>−0</sub>, <italic>k</italic><sub>1</sub>, <italic>k</italic><sub>−1</sub>, <italic>k</italic><sub>2</sub>, <italic>k</italic><sub>−2</sub>, and <italic>k</italic><sub>3</sub>. Of these, <italic>k</italic><sub>0</sub> and <italic>k</italic><sub>1</sub> are second order rate constants (units of M<sup>−1</sup>s<sup>−1</sup>) and the rest are first order rate constants (units of s<sup>−1</sup>). The experimental literature provides reasonable estimates for <italic>k</italic><sub>−0</sub>/<italic>k</italic><sub>0</sub>, <italic>k</italic><sub>1</sub> and <italic>k</italic><sub>−1</sub>, given below, but not for <italic>k</italic><sub>2</sub>, which we go on to calculate. <italic>k</italic><sub>3</sub> is also relatively unknown, which we discuss later.</p><p>For the synapsis-by-protein-dimerization model, as exemplified by Tn10 and Tn5, we make the same assumptions as detailed above, except that the free transposase is monomeric.</p><p>In our model, protein expression is instantaneous, producing a specified number of transposase dimers per transposon copy. Unless stated otherwise this is 500 dimers in the mariner model and 1000 monomers in the prokaryotic model. Instantaneous expression is necessary to simulate the reaction under the fast diffusion regime which prevails in vitro. Under the in vivo slow-diffusion regime, the reaction is so slow that the dynamics of protein expression and degradation are irrelevant.</p></sec><sec id="s4-4"><title>Choice of the kinetic parameters for DNA binding–<italic>k</italic><sub>0</sub>, <italic>k</italic><sub>−0</sub>, <italic>k</italic><sub>1</sub>, and <italic>k</italic><sub>−1</sub></title><p>The Hsmar1 and Mos1 transposases are 37% identical and align throughout their entire length with only a single residue indel. The respective proteins will therefore have very similar three dimensional structures. A crystal structure for the Mos1 transposase revealed that it has an N-terminal DNA binding domain with two helix-turn-helix (H-T-H) motifs that are inserted into the major groove of the DNA close to the end of the transposon (<xref ref-type="bibr" rid="bib40">Richardson et al., 2009</xref>). We can therefore expect that the dissociation constant for transposon end binding (<italic>K</italic><sub>d</sub> = <italic>k</italic><sub>−1</sub>/<italic>k</italic><sub>1</sub>) will be in the mid-picomolar concentration range. The lactose repressor (LacI) is one of the most thoroughly characterized H-T-H proteins, with many studies addressing the kinetic parameters of DNA binding. The LacI dimer is similar to a monomer of Hsmar1 transposase in that it binds to operator sequences by a pair of H-T-H motifs, which are located in adjacent major grooves. We therefore elected to use the kinetic parameters for LacI to model transposase binding to the transposon end. Estimated parameters for the behavior of the transposase are perhaps not ideal. However, our conclusions do not depend on the absolute values of the parameters because they have minimal effect beyond changing the scale on the axes of the graphs generated by the simulation. The relationships between the various parameters, which provide the key insights of the work, remain unchanged, as explained further below and in the main text.</p><p>The values for the lactose repressor’s dissociation constant quoted in the textbooks and the literature range from about 10<sup>−9</sup> to 10<sup>−14</sup> M. However, it is now generally accepted that the <italic>K</italic><sub>d</sub> for LacI, and other dimeric H-T-H proteins, binding to their primary sites is in the mid to high picomolar concentration range. Many of the very low values quoted in the literature date from the late 1970s and early 1980s, and were determined at unphysiologically low salt concentrations. Another factor that complicates the lactose repressor literature is the complexity of the multiple binding sites encoded by the operator, which promote DNA looping and the mutual stabilization of dimers bound in cis to each other. Here we have used the kinetic parameters determined by Wells and Matthews and colleagues for LacI binding to a single site encoded on a 40 bp fragment of the lac operator under a physiological salt concentration (Table 3 of <xref ref-type="bibr" rid="bib21">Hsieh et al. (1987)</xref>).</p><p>In the course of our experiments with Hsmar1 transposase we have never recorded significant inhibition of the reaction by the moderate amounts of non-specific DNA used as a target in many of our experiments. This suggests that the affinity of the transposase for non-specific DNA is relatively low. We have therefore elected to adopt the equilibrium constant (<italic>K</italic><sub>d</sub> = <italic>k</italic><sub>−0</sub>/<italic>k</italic><sub>0</sub>) of LacI for non-specific DNA because it is probably fairly typical for trans-acting DNA binding proteins. We adopted the value provided by Von Hippel and colleagues because it was determined at a similar physiological salt concentration to the parameters noted above for site-specific binding (<xref ref-type="bibr" rid="bib37">Revzin and Von Hippel, 1977</xref>). We used this equilibrium constant to assign arbitrary association and dissociation rate constants, <italic>k</italic><sub>0</sub> and <italic>k</italic><sub>−0</sub>. It is difficult to define these rate constants precisely because DNA binding proteins have two distinct modes of non-specific interaction. Firstly, binding from the bulk solution is probably rapid with a rate equal to a sizable fraction of a diffusion-limited reaction. This is followed by a one-dimensional diffusion phase, during which the protein visits non-specific sites that are considered to be 1 bp apart. The equilibrium constant, <italic>K</italic><sub>d</sub>, is thus the product of one- and three-dimensional events, each of which has its own <italic>k</italic><sub><italic>on</italic></sub> and <italic>k</italic><sub><italic>off</italic></sub> components. Fortunately, the absolute values of <italic>k</italic><sub>−0</sub> and <italic>k</italic><sub>0</sub> are not important in the simulation because they serve only to define the equilibrium constant for non-specific sites, which in turn specifies the concentration of free transposase available for transposon end binding in a genome of a given size. The validity of this assumption was confirmed by running the simulation using several different pairs of values for <italic>k</italic><sub>0</sub> and <italic>k</italic><sub>−0</sub>: the output was unchanged provided that <italic>k</italic><sub>−0</sub>/<italic>k</italic><sub>0</sub> was always equal to the experimentally determined <italic>K</italic><sub>d</sub> (not shown).</p><p>The values for the specific and non-specific parameters are tabulated in <xref ref-type="table" rid="tbl1">Table 1</xref>, where they are also converted into the absolute units used in the simulation (i.e., numbers of molecules per cubic micrometer). The dissociation constant for non-specific DNA provided by Von Hippel and colleagues is expressed in units of per base pair. We have assumed a haploid genome size of 3 × 10<sup>9</sup> bp, which provides 6 × 10<sup>9</sup> non-specific binding sites per diploid genome. The volume of the nucleus is assumed to be 500 µm<sup>3</sup>.<table-wrap id="tbl1" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.014</object-id><label>Table 1.</label><caption><p>Kinetic parameters for DNA binding</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.014">http://dx.doi.org/10.7554/eLife.00668.014</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Reaction</th><th><italic>K</italic><sub>d</sub> (M)</th><th><italic>k</italic><sub>0</sub>, <italic>k</italic><sub>1</sub> (M<sup>−1</sup>s<sup>−1</sup>)</th><th><italic>k</italic><sub>0</sub>, <italic>k</italic><sub>1</sub> (µm<sup>3</sup> s<sup>−1</sup>)</th><th><italic>k</italic><sub>−0</sub>, <italic>k</italic><sub>−1</sub> (s<sup>−1</sup>)</th><th><italic>t</italic><sub>1/2</sub> (s)</th></tr></thead><tbody><tr><td>LacI specific binding</td><td>3.1 × 10<sup>−11</sup></td><td>3.8 × 10<sup>8</sup> <xref ref-type="table-fn" rid="tblfn1">*</xref></td><td align="char" char=".">0.633</td><td>1.2 × 10<sup>−2</sup> <xref ref-type="table-fn" rid="tblfn1">*</xref></td><td>58 <xref ref-type="table-fn" rid="tblfn2">†</xref></td></tr><tr><td>LacI non-specific binding</td><td>1.4 × 10<sup>−5</sup> <xref ref-type="table-fn" rid="tblfn3">‡</xref></td><td>9.9 × 10<sup>6</sup> <xref ref-type="table-fn" rid="tblfn4">§</xref></td><td align="char" char=".">0.0166</td><td>139 <xref ref-type="table-fn" rid="tblfn4">§</xref></td><td>5 × 10<sup>−3</sup> <xref ref-type="table-fn" rid="tblfn2">†</xref></td></tr></tbody></table><table-wrap-foot><fn id="tblfn1"><label>*</label><p>From Table 3 of <xref ref-type="bibr" rid="bib21">Hsieh et al. (1987)</xref>.</p></fn><fn id="tblfn2"><label>†</label><p>This is the experimental estimate for LacI bound to non-specific DNA dissociating to the bulk solution (<xref ref-type="bibr" rid="bib13">Elf et al., 2007</xref>). During this time it will have visited about 85 non-specific sites by one-dimensional diffusion. However, for the purposes of the simulation it is convenient to subsume the two phases into a single bulk behavior of the system (see text for details).</p></fn><fn id="tblfn3"><label>‡</label><p>From <xref ref-type="bibr" rid="bib37">Revzin and Von Hippel (1977)</xref>.</p></fn><fn id="tblfn4"><label>§</label><p>Estimated from <italic>K</italic><sub>d</sub> using the relationship <italic>K</italic><sub>d</sub> = <italic>k</italic><sub>−1</sub>/<italic>k</italic><sub>1</sub> or <italic>K</italic><sub>d</sub> = <italic>k</italic><sub>−0</sub>/<italic>k</italic><sub>0</sub>.</p></fn></table-wrap-foot></table-wrap></p></sec><sec id="s4-5"><title>Rate of synapsis, <italic>k</italic><sub>2</sub></title><p>We calculated the first-order reaction rate constant <italic>k</italic><sub>2</sub> as follows. Throughout this discussion, T<sub>2</sub> represents freely diffusing transposase, which is shown as T<sub>2,free</sub> in the model diagram.</p><p>Consider the two elementary reactions from our model that are shown at the bottom-left and the top-right regions of <xref ref-type="fig" rid="fig1">Figure 1I</xref>, respectively. Ignoring the transposon’s internal DNA sequences for now, these reactions are (1) T<sub>2</sub> + B → BT<sub>2</sub> and (2) AT<sub>2</sub> + B → AT<sub>2</sub>B. Both reactions involve the binding of transposase to DNA site B, so they are chemically identical. This might suggest that they would have the same reaction rates. However, in actuality, their reaction rates differ because the rate of diffusion for the free transposase is different from that of the DNA-bound transposase. Expressing their reaction rate constants as <italic>k</italic><sub>1</sub> for reaction 1 (as in the model diagram) and <inline-formula><mml:math id="inf1"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> for reaction 2, the respective reaction rates are<disp-formula id="equ1"><label>(1)</label><mml:math id="m1"><mml:mrow><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mrow><mml:mtext>BT</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:msub><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="italic">dt</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:msub><mml:mi>k</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mtext>T</mml:mtext><mml:mn>2</mml:mn></mml:msub><mml:mo>]</mml:mo></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mtext>B</mml:mtext><mml:mo>]</mml:mo></mml:mrow><mml:mtext>,</mml:mtext></mml:mrow></mml:math></disp-formula><disp-formula id="equ2"><label>(2)</label><mml:math id="m2"><mml:mrow><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mrow><mml:mtext>AT</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:msub><mml:mtext>B</mml:mtext><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="italic">dt</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mrow><mml:mtext>AT</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:msub><mml:mo>]</mml:mo></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mtext>B</mml:mtext><mml:mo>]</mml:mo></mml:mrow><mml:mtext>.</mml:mtext></mml:mrow></mml:math></disp-formula></p><p>The Collins and Kimball reaction rate theory (<xref ref-type="bibr" rid="bib39">Rice, 1985</xref>) enables us to expand these reaction rate constants into one contribution that arises from diffusion and a second that arises from the binding activation energy. According to this theory, the expanded reaction rate constants are<disp-formula id="equ3"><label>(3)</label><mml:math id="m3"><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mrow><mml:mi mathvariant="italic">act</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula><disp-formula id="equ4"><label>(4)</label><mml:math id="m4"><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mrow><mml:mi mathvariant="italic">act</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>k</italic><sub>1,<italic>diff</italic></sub> and <inline-formula><mml:math id="inf2"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> are diffusion-limited reaction rate constants and <italic>k</italic><sub><italic>act</italic></sub> is an activation-limited reaction rate constant. These equations include the same <italic>k</italic><sub><italic>act</italic></sub> value because we assumed above that the two reactions are chemically identical; under this assumption, they should have the same binding activation energy and thus the same activation-limited reaction rate constants.</p><p>The necessary diffusion-limited reaction rate constants can be calculated from equations that Berg derived for the diffusion-limited association rates of proteins to DNA (<xref ref-type="bibr" rid="bib3">Berg, 1984</xref>). Berg’s equation 16 is for the association rate of a freely diffusing protein to a specific DNA site, while his equation 27 is for the association rate of a DNA-bound protein to another specific DNA site. In our notation, these equations are<disp-formula id="equ5"><label>(5)</label><mml:math id="m5"><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>4</mml:mn><mml:mi>π</mml:mi><mml:msub><mml:mi>D</mml:mi><mml:mi>p</mml:mi></mml:msub><mml:mi>R</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>D</mml:mi><mml:mi>s</mml:mi></mml:msub><mml:mi>a</mml:mi><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mi>R</mml:mi><mml:mi>a</mml:mi></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula><disp-formula id="equ6"><label>(6)</label><mml:math id="m6"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow><mml:mo>′</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:mn>1.4</mml:mn><mml:msub><mml:mi>D</mml:mi><mml:mi>s</mml:mi></mml:msub><mml:mi>a</mml:mi><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mi>R</mml:mi><mml:mi>a</mml:mi></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p><p>These equations were derived by treating DNA with the worm-like chain model, in which DNA is represented as a long thin semi-flexible filament that bends smoothly over the course of its length. This model was recently shown to be reasonably accurate for modeling DNA dynamics (<xref ref-type="bibr" rid="bib36">Petrov et al., 2006</xref>). The central model parameter is the filament persistence length, <italic>a</italic>, which is the characteristic length for filament bending. According to the model, DNA binding sites diffuse rapidly within their local regions, while also gradually diffusing away to more distant regions. The rate of this local ‘segmental diffusion’ is quantified with the diffusion coefficient <italic>D</italic><sub><italic>s</italic></sub>. <italic>D</italic><sub><italic>s</italic></sub> is the translational diffusion coefficient of a free DNA fragment with length <italic>a</italic> (<xref ref-type="bibr" rid="bib3">Berg, 1984</xref>). Continuing with the above equations, <italic>R</italic> represents the distance over which a specific interaction occurs between transposase and its DNA binding sites. It arises from the assumption that the two reactants (the transposase binding domain and a DNA binding site, in this case) can be treated as hard spheres that react upon collision. Finally, <italic>D</italic><sub><italic>p</italic></sub> is the diffusion coefficient of free transposase relative to the center of mass of the DNA.</p><p>Two final equations need to be introduced to enable us to calculate the rate of synapsis, <italic>k</italic><sub>2</sub>. The reaction rate equation for reaction 2, given in <xref ref-type="disp-formula" rid="equ2">equation 2</xref>, can be rearranged by grouping <inline-formula><mml:math id="inf3"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> and [B],<disp-formula id="equ7"><label>(7)</label><mml:math id="m7"><mml:mrow><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mrow><mml:mtext>AT</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:msub><mml:mtext>B</mml:mtext><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="italic">dt</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup><mml:mrow><mml:mo>[</mml:mo><mml:mtext>B</mml:mtext><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mtext>AT</mml:mtext><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>k</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mtext>AT</mml:mtext><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p><p>The latter equality defines <italic>k</italic><sub>2</sub>, a pseudo-first order reaction rate constant, as <inline-formula><mml:math id="inf4"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup><mml:mrow><mml:mo>[</mml:mo><mml:mtext>B</mml:mtext><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula>. Here, [B] is the concentration of the B DNA binding site in the vicinity of the A binding site. We used the following empirical equation from <xref ref-type="bibr" rid="bib41">Ringrose et al. (1999)</xref> to find [B]:<disp-formula id="equ8"><label>(8)</label><mml:math id="m8"><mml:mrow><mml:msub><mml:mi>j</mml:mi><mml:mi>M</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>4</mml:mn><mml:mi>a</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mn>4</mml:mn></mml:msup><mml:mi>b</mml:mi></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>exp</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mo>−</mml:mo><mml:mn>460</mml:mn><mml:msup><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mn>6.25</mml:mn><mml:msup><mml:mi>b</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mfrac><mml:mrow><mml:mn>1.25</mml:mn><mml:mo>⋅</mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mn>5</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:msup><mml:mi>a</mml:mi><mml:mn>3</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p><p><italic>j</italic><sub><italic>M</italic></sub> is the local concentration of one DNA site in the vicinity of another in M, <italic>a</italic> is the DNA persistence length in nanometers, and <italic>b</italic> is the separation between the sites in base pairs.</p><p><xref ref-type="table" rid="tbl2">Table 2</xref> lists the values that we used and derived to estimate <italic>k</italic><sub>2</sub>. In brief, we used <xref ref-type="disp-formula" rid="equ3">equation 3</xref>, including an experimental value for <italic>k</italic><sub>1</sub> and a calculated value for <italic>k</italic><sub>1,<italic>diff</italic></sub> from <xref ref-type="disp-formula" rid="equ5">equation 5</xref>, to calculate <italic>k</italic><sub><italic>act</italic></sub>. We then substituted <italic>k</italic><sub><italic>act</italic></sub> into <xref ref-type="disp-formula" rid="equ4">equation 4</xref>, along with a calculated value for <inline-formula><mml:math id="inf5"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> from <xref ref-type="disp-formula" rid="equ6">equation 6</xref>, to estimate <inline-formula><mml:math id="inf6"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>. Finally, we used <xref ref-type="disp-formula" rid="equ7">equations 7</xref> and <xref ref-type="disp-formula" rid="equ8">8</xref> to estimate <italic>k</italic><sub>2</sub>, the rate of synapsis, from <inline-formula><mml:math id="inf7"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>.<table-wrap id="tbl2" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.015</object-id><label>Table 2.</label><caption><p>Estimation of <italic>k</italic><sub>2</sub> for transposase-DNA interactions in vitro</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.015">http://dx.doi.org/10.7554/eLife.00668.015</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Quantity</th><th>Symbol</th><th>Value</th><th>Source</th></tr></thead><tbody><tr><td>Reaction 1 rate constant</td><td><italic>k</italic><sub>1</sub></td><td>3.8 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup></td><td>Table 3 of <xref ref-type="bibr" rid="bib21">Hsieh et al. (1987)</xref>.</td></tr><tr><td>DNA persistence length</td><td><italic>a</italic></td><td>50 nm</td><td>Well established (e.g. <xref ref-type="bibr" rid="bib36">Petrov et al., 2006</xref>).</td></tr><tr><td>Reaction radius</td><td><italic>R</italic></td><td>2.5 nm</td><td>The transposase dimeric binding domain is 111 amino acids and weighs 26 kDa. From Note 3 of <xref ref-type="bibr" rid="bib1">Andrews (2012)</xref>, its radius is about 1.9 nm. Add to this about 0.6 nm to account for the 0.34 nm length of 1 DNA basepair (because transposase specificity is accurate to 1 bp), combined with the 1 nm radius of DNA.</td></tr><tr><td>Segmental diffusion coefficient</td><td><italic>D</italic><sub><italic>s</italic></sub></td><td>27 µm<sup>2</sup>s<sup>−1</sup></td><td>Figure 3 of <xref ref-type="bibr" rid="bib36">Petrov et al. (2006)</xref> for 147 bp, which is 1 DNA persistence length.</td></tr><tr><td>Transposase diffusion coefficient</td><td><italic>D</italic><sub><italic>p</italic></sub></td><td>61 µm<sup>2</sup>s<sup>−1</sup></td><td>Note 3 of <xref ref-type="bibr" rid="bib1">Andrews (2012)</xref>, using a molecular weight of 80.7 kDa for transposase. Assumes DNA center of mass is effectively stationary.</td></tr><tr><td>Diffusion-limited part of <italic>k</italic><sub>1</sub></td><td><italic>k</italic><sub>1,<italic>diff</italic></sub></td><td>1.5 × 10<sup>9</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ5">Equation 5</xref>.</td></tr><tr><td>Activation-limited reaction rate</td><td><italic>k</italic><sub><italic>act</italic></sub></td><td>5.2 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ3">Equation 3</xref>.</td></tr><tr><td>Diffusion-limited part of <inline-formula><mml:math id="inf8"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula></td><td><inline-formula><mml:math id="inf9"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula></td><td>4.2 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ6">Equation 6</xref>.</td></tr><tr><td>Reaction 2 rate constant</td><td><inline-formula><mml:math id="inf10"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula></td><td>2.3 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ4">Equation 4</xref>.</td></tr><tr><td>Effective local DNA site concentration</td><td>[B]</td><td>5.5 × 10<sup>−8</sup> M</td><td><xref ref-type="disp-formula" rid="equ8">Equation 8</xref>, using a 1287 bp transposon length.</td></tr><tr><td>Transposase-transposon end dissociation rate constant</td><td><italic>k</italic><sub>−1</sub></td><td>5.8 × 10<sup>−4</sup> s<sup>−1</sup></td><td>Estimated from <xref ref-type="fig" rid="fig3">Figure 3C</xref>.</td></tr><tr><td>Reaction 2, pseudo-first order rate constant for synapsis</td><td><italic>k</italic><sub>2</sub></td><td>12.7 s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ7">Equation 7</xref>.</td></tr><tr><td>Reaction 2, pseudo-first order rate constant for synapsis with an open circular substrate in vitro.</td><td><italic>k</italic><sub>2</sub></td><td>9.6 × 10<sup>−5</sup> s<sup>−1</sup></td><td>Estimated using the equation t<sub>1/2</sub> = ln2/<italic>k</italic><sub>2</sub>. t<sub>1/2</sub> is the time taken to consume one quarter of the substrate in <xref ref-type="fig" rid="fig3">Figure 3D</xref>. See <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref> for an explanation of why the time to consume one quarter of the substrate is used, rather than the time to consume half of the substrate.</td></tr></tbody></table></table-wrap></p><p>These theoretical considerations allow us to develop a model of an idealized transposition reaction, in which the monomers within a transposase dimer bind the first and second transposon ends with equal affinity. In this idealized situation, synapsis of the transposon ends is simply a product of the sequential binding of the transposase dimer to sites at opposite ends of the element. We then go on to consider a more realistic model in which allosteric interactions between the subunits reduces the affinity of the developing transpososome for the second transposon end. The magnitude of this effect is provided by experimental estimates of the rate of synapsis (<xref ref-type="fig" rid="fig3">Figure 3D</xref>, <xref ref-type="table" rid="tbl2">Table 2</xref>; <xref ref-type="bibr" rid="bib8">Claeys Bouuaert and Chalmers, 2010</xref>; <xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>).</p><p>In <xref ref-type="table" rid="tbl2">Table 2</xref>, the reaction radius, <italic>R</italic>, is a fairly rough estimate. However, this is not a major concern because reaction rates are relatively insensitive to the reaction radius (<xref ref-type="bibr" rid="bib3">Berg, 1984</xref>) (they scale as <italic>R</italic><sup>1/3</sup>). Other values are likely to be more accurate. A result of the calculations shown there is that <inline-formula><mml:math id="inf11"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> is only about a factor of 1.6 slower than <italic>k</italic><sub>1</sub>. This indicates that the slower diffusion of DNA-bound transposase than of free transposase has only modest impact for in vitro experiments.</p><p>For Tn10/5 there is no experimental data regarding the affinity of the transposon-end-bound monomers when they collide with each other and achieve synapsis (see kinetic model in <xref ref-type="fig" rid="fig1">Figure 1B</xref>). For the purposes of the simulation we assigned an association rate constant equal to that for the mariner transposase binding to a transposon end. In the respective prokaryotic and eukaryotic models, the rate of synapsis is therefore determined by the same association rate constant and the rate of diffusion, which is identical in both systems. The outputs of the respective simulations are therefore directly comparable and reveal differences in the underlying kinetic models.</p></sec><sec id="s4-6"><title>The effects of macromolecular crowding in vivo</title><p>Thus far we have considered reactions in dilute solution. To provide a more realistic description of the reaction in vivo we next had to account for macromolecular crowding, which reduces the rate of diffusion. We therefore computed the same numbers that are shown in <xref ref-type="table" rid="tbl2">Table 2</xref>, but for the in vivo situation, which are shown in <xref ref-type="table" rid="tbl3">Table 3</xref>. The only parameters sensitive to the high viscosity in vivo are the association rate constants, which includes the rate of synapsis. Our calculations show that the association rate constant, <italic>k</italic><sub>1</sub>, suffers a modest twofold penalty in vivo (<xref ref-type="table" rid="tbl3">Table 3</xref>). This relieves OPI and increases the rate of transposition twofold (not shown). However, this is precisely offset by the penalty on non-specific binding, <italic>k</italic><sub>0</sub>, which achieves the opposite. The results are very different for the rate of synapsis, which suffers an in vivo diffusion penalty of 2.8 × 10<sup>4</sup>-fold (<xref ref-type="table" rid="tbl2 tbl3">Tables 2 and 3</xref>). However, many DNA interactions are known to occur over much longer-ranges than those between transposon ends. It is therefore possible that other factors in vivo mitigate the severity of the diffusion penalty. For example, there is evidence that the factor <italic>j</italic><sub><italic>M</italic></sub> is doubled by the chromatinization of DNA (<xref ref-type="bibr" rid="bib41">Ringrose et al., 1999</xref>). The in vivo values in <xref ref-type="table" rid="tbl3">Table 3</xref> should therefore be seen as initial estimates rather than an exhaustive treatment of the topic. Nevertheless, our calculated rate of synapsis is very close to the rates of lox/cre recombination measured in vivo with similarly separated recombination sites (<xref ref-type="bibr" rid="bib41">Ringrose et al., 1999</xref>).<table-wrap id="tbl3" position="float"><object-id pub-id-type="doi">10.7554/eLife.00668.016</object-id><label>Table 3.</label><caption><p>Estimation of in vivo transposase-DNA interactions</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00668.016">http://dx.doi.org/10.7554/eLife.00668.016</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Quantity</th><th>Symbol</th><th>Value</th><th>Source</th></tr></thead><tbody><tr><td>DNA persistence length</td><td><italic>a</italic></td><td>50 nm</td><td>See <xref ref-type="table" rid="tbl2">Table 2</xref>.</td></tr><tr><td>Reaction radius</td><td><italic>R</italic></td><td>2.5 nm</td><td>See <xref ref-type="table" rid="tbl2">Table 2</xref>.</td></tr><tr><td>Segmental diffusion coefficient</td><td><italic>D</italic><sub><italic>s</italic></sub></td><td>5 × 10<sup>−4</sup> µm<sup>2</sup>s<sup>−1</sup></td><td>From Figure 2 caption of <xref ref-type="bibr" rid="bib33">Marshall et al. (1997)</xref>.</td></tr><tr><td>Transposase diffusion coefficient</td><td><italic>D</italic><sub><italic>p</italic></sub></td><td>15 µm<sup>2</sup>s<sup>−1</sup></td><td>¼ of <xref ref-type="table" rid="tbl2">Table 2</xref> value, based on Note 3 of <xref ref-type="bibr" rid="bib1">Andrews, (2012)</xref>.</td></tr><tr><td>Diffusion-limited part of <italic>k</italic><sub>1</sub></td><td><italic>k</italic><sub>1,<italic>diff</italic></sub></td><td>2.9 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ5">Equation 5</xref>.</td></tr><tr><td>Activation-limited reaction rate</td><td><italic>k</italic><sub><italic>act</italic></sub></td><td>5.2 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup></td><td>From <xref ref-type="table" rid="tbl2">Table 2</xref>.</td></tr><tr><td>Diffusion-limited part of <inline-formula><mml:math id="inf12"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula></td><td><inline-formula><mml:math id="inf13"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="italic">diff</mml:mi></mml:mrow><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula></td><td>7.8 × 10<sup>3</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ6">Equation 6</xref>.</td></tr><tr><td>Reaction 1 rate constant</td><td><italic>k</italic><sub>1</sub></td><td>1.9 × 10<sup>8</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ3">Equation 3</xref>.</td></tr><tr><td>Reaction 2 rate constant</td><td><inline-formula><mml:math id="inf14"><mml:mrow><mml:msubsup><mml:mi>k</mml:mi><mml:mn>1</mml:mn><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula></td><td>7.8 × 10<sup>3</sup> M<sup>−1</sup>s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ4">Equation 4</xref>.</td></tr><tr><td>Effective local DNA site concentration</td><td>[B]</td><td>5.5 × 10<sup>−8</sup> M</td><td>See <xref ref-type="table" rid="tbl2">Table 2</xref>.</td></tr><tr><td>Reaction 2 pseudo-first order rate constant</td><td><italic>k</italic><sub>2</sub></td><td>4.3 × 10<sup>−4</sup> s<sup>−1</sup></td><td><xref ref-type="disp-formula" rid="equ7">Equation 7</xref>.</td></tr></tbody></table></table-wrap></p></sec><sec id="s4-7"><title>Synapsis unbinding, <italic>k</italic><sub>−2</sub></title><p>In the simplest model for synapsis the rate of the reverse reaction, <italic>k</italic><sub>−2</sub>, would be equal to the rate of transposon end unbinding, <italic>k</italic><sub>−1</sub> (see <xref ref-type="fig" rid="fig1">Figure 1I</xref>). However, in practice we consider synapsis to be essentially irreversible for the following reasons. Transposition reaction mixtures in which the catalytic Mg<sup>2+</sup> ions are replaced by Ca<sup>2+</sup> support transpososome assembly, but none of the chemical steps of the reaction. After the addition of the catalytic metal ion, the first nick at the transposon end is detected immediately (<xref ref-type="bibr" rid="bib8">Claeys Bouuaert and Chalmers, 2010</xref>). This shows that nicking is very much faster than synapsis or synapsis unbinding. Synapsis in the presence of the catalytic metal ion is therefore essentially irreversible because of the loss of enthalpy associated with hydrolysis of the phosphodiester bond. In principle, it is possible that there exists some form of unstable synaptic complex that matures into the stable form that immediately precedes the first chemical reaction. If such a complex existed it would be irrelevant to our experiments, in which the rate of synapsis is estimated from the rate of the first nick. When we consider the rate of synapsis we are therefore restricting ourselves to ‘productive synapses’, which yield the first nick. For the purposes of the simulation we have set the rate constant for synapsis unbinding (<italic>k</italic><sub>−2</sub>) at 10<sup>−10</sup> s<sup>−1</sup>. Synapsis in the model is therefore essentially irreversible.</p></sec><sec id="s4-8"><title>Allostery and experimental estimates of <italic>k</italic><sub>2</sub> and <italic>k</italic><sub>−1</sub></title><p>When transposase dimer binds to the first transposon end it undergoes a conformational change that lowers the affinity of the unoccupied DNA binding domain for the second transposon end compared to the first (for details see main text and <xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>). This means that the actual value of <italic>k</italic><sub>2</sub> will be much lower than that calculated above for the idealized reaction, in which synapsis is the product of sequential, chemically identical, collisions between DNA binding domains and transposon ends. To account for the allosteric interactions between the transposase subunits we therefore require experimental estimates of <italic>k</italic><sub>2</sub> and <italic>k</italic><sub>−1</sub>.</p><p>Estimates for <italic>k</italic><sub>−1</sub> and <italic>k</italic><sub>2</sub> are provided by the kinetics of the in vitro reaction. In <xref ref-type="fig" rid="fig3">Figure 3C</xref> we see that about 50% of the transposon end present in SEC2 is released after 20 min incubation in the presence of cold competitor DNA. This corresponds to a dissociation rate constant of 5.8 × 10<sup>−4</sup> s<sup>−1</sup>. We have previously shown that first strand nicking depends on, and is much faster than, synapsis (<xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref> and ‘Discussion’ above). Consumption of the supercoiled substrate therefore provides an estimate of the rate of synapsis. Under optimal reaction conditions, with about one dimer of active transposase per transposon, only one half of the substrate is initially occupied by a single dimer and is therefore able to react. As explained above, the random association of transposase dimers and transposon ends at the start of the reaction means that a quarter would be occupied by two dimers and would suffer OPI and the other quarter would be completely unoccupied. The time required to consume one quarter of the substrate therefore approximates the half-time for synapsis. Since eukaryotic DNA has very little free supercoiling, the most relevant rate of synapsis is perhaps provided by an open circular substrate. In <xref ref-type="fig" rid="fig3">Figure 3D</xref>, the time taken to consume one quarter of the open circular substrate is about 2 hr, which corresponds to a pseudo-first order rate constant of 9.6 × 10<sup>−5</sup> s<sup>−1</sup>. Note that the reason we are able to estimate the rate of synapsis is owing to the optimal reaction condition in vitro where <italic>k</italic><sub>2</sub> is largely independent of <italic>k</italic><sub>0</sub>, <italic>k</italic><sub>−0</sub>, <italic>k</italic><sub>1</sub> and <italic>k</italic><sub>−1</sub>. This is because there is very little non-specific DNA present and the transposase concentration is such that most dimers are engaged in productive interactions with transposons.</p></sec><sec id="s4-9"><title>Rate of cleavage, integration and maturation of the product, <italic>k</italic><sub>3</sub></title><p>The rate limiting steps of the reaction with various substrates has been determined previously (<xref ref-type="bibr" rid="bib8">Claeys Bouuaert and Chalmers, 2010</xref>; <xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>). With relaxed plasmid and short linear substrates, synapsis is the rate limiting step. Supercoiling in the substrate accelerates synapsis because the transposon ends have a high relative concentration in the plectosome and a favorable angular distribution. In an in vitro reaction, where the plasmid has more than twice its natural level of free supercoiling, the rate of synapsis is faster than cleavage of the second strand at the transposon end. Note that second strand cleavage at the second transposon end releases the transposon from the donor site (illustrated in <xref ref-type="fig" rid="fig2">Figure 2A</xref>). In a kinetic analysis of staged in vitro reactions with a supercoiled substrate, a small amount of excised transposon is observed at early time points [this is below the region of the gels shown in <xref ref-type="fig" rid="fig2">Figure 2B</xref>, but can be seen in Figure 5B of <xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al. (2011)</xref>]. The rate of integration must therefore be faster than the rate of second strand cleavage. The half life of the excised transposon is approximated by the time required to convert a quarter of the supercoiled substrate into product. From <xref ref-type="fig" rid="fig2">Figure 2B</xref> and other experiments this appears to be about 8 min which corresponds to a rate constant for integration of 1.4 × 10<sup>−3</sup> s<sup>−1</sup>.</p><p>Although this value for <italic>k</italic><sub>3</sub> is appropriate in the simulations using in vitro parameters, it would be unrealistically short in an in vivo situation where the integration complex must be disassembled and the empty donor site restored. However, our reference in vitro simulation (<xref ref-type="fig" rid="fig3">Figure 3E</xref>) is relatively insensitive to the value of <italic>k</italic><sub>3</sub>. This is because the extent to which the Hsmar1 specific factors, namely the slow dissociation rate of SEC2 and the slow recruitment of the second end, slow synapsis. Furthermore, once a steady-state rate is established the length of the maturation process does not matter because the rates of pathway entry and exit are equal. Thus, if <italic>k</italic><sub>3</sub> is reduced by three orders of magnitude to 1.4 × 10<sup>−6</sup> s<sup>−1</sup> (t<sub>1/2</sub> = 5.7 days), the final rate of transposon amplification does not change significantly. However, it does take slightly longer to be achieved because of the time required for pathway entry and exit rates to balance (not shown). Note that we assume that every transposon that initiates catalysis goes on to complete excision and integration. This is supported by the in vitro reaction systems for Hsmar1 and Tn10 (<xref ref-type="bibr" rid="bib6">Chalmers and Kleckner, 1994</xref>, <xref ref-type="bibr" rid="bib7">1996</xref>; <xref ref-type="bibr" rid="bib8">Claeys Bouuaert and Chalmers, 2010</xref>; <xref ref-type="bibr" rid="bib9">Claeys Bouuaert et al., 2011</xref>). The efficiency of the reaction in vivo is unlikely to be so high and this is dealt with in <xref ref-type="fig" rid="fig3">Figure 3G</xref> and the section of text entitled ‘The efficiency of transposition in vivo.’</p></sec><sec id="s4-10"><title>Transposase expression level</title><p>Many transposons, and mariner in particular, have a wide host range. The promoters driving expression of transposase genes are therefore assumed to rely on general transcription factors. However, there is little published data pertaining to the possible strength of the promoters. Expression of the Hsmar1 promoter has been detected using a luciferase reporter system, but quantification is lacking for the contribution of a single copy of the element (<xref ref-type="bibr" rid="bib34">Miskey et al., 2007</xref>). A recent proteomic study quantified the abundance of 5000 vertebrate proteins in mouse cells (<xref ref-type="bibr" rid="bib47">Schwanhausser et al., 2011</xref>). The copy number per cell for most proteins ranged from about 100 to 1 million: The median was 16,000 and the mode was about 5000. In the simulation we set the transposase expression at 500 dimers per transposon. Although this is at the low end of the range for vertebrate proteins, it is not unusually low and allows for up to 1000 copies of the transposon before protein levels would reach the high end of the range. However, we note, once again, that our conclusions do not depend on the absolute values of the parameters because they have minimal effect beyond changing the scale on axes of the graphs generated by the simulation. The relationships between the various parameters, which provide the key insights of the work, remain unchanged.</p></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>We would like to thank Yves Bigot, Nancy Kleckner, Boyan Bonev, Steve Halford, Steve Kowalczykowski, Csaba Miskey, Zoltan Ivics, Deepesh Agarwal, and Andrew Howlett for helpful advice; Neil Walker for technical assistance; and Mark Calleja and Jenny Barna for maintaining CamGrid. We would also like to thank the referees for valuable contributions to this work. This publication is dedicated to the memory of Peter and Barbara Blakey whose generous bequest constructed the laboratory infrastructure.</p></ack><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>The authors declare that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>CCB, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con2"><p>DL, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con3"><p>RC, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con4"><p>KL, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con5"><p>SSA, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn></fn-group></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Andrews</surname><given-names>SS</given-names></name></person-group><year>2012</year><article-title>Spatial and stochastic cellular modeling with the Smoldyn simulator</article-title><source>Methods Mol Biol</source><volume>804</volume><fpage>519</fpage><lpage>42</lpage><pub-id pub-id-type="doi">10.1007/978-1-61779-361-5_26</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Auge-Gouillou</surname><given-names>C</given-names></name><name><surname>Brillet</surname><given-names>B</given-names></name><name><surname>Germon</surname><given-names>S</given-names></name><name><surname>Hamelin</surname><given-names>MH</given-names></name><name><surname>Bigot</surname><given-names>Y</given-names></name></person-group><year>2005</year><article-title>Mariner Mos1 transposase dimerizes prior to ITR binding</article-title><source>J Mol Biol</source><volume>351</volume><fpage>117</fpage><lpage>30</lpage><pub-id pub-id-type="doi">10.1016/j.jmb.2005.05.019</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berg</surname><given-names>OG</given-names></name></person-group><year>1984</year><article-title>Diffusion-controlled protein-DNA association: influence of segmental diffusion of the DNA</article-title><source>Biopolymers</source><volume>23</volume><fpage>1869</fpage><lpage>89</lpage><pub-id pub-id-type="doi">10.1002/bip.360231005</pub-id></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bray</surname><given-names>D</given-names></name><name><surname>Lay</surname><given-names>S</given-names></name></person-group><year>1997</year><article-title>Computer-based analysis of the binding steps in protein complex formation</article-title><source>Proc Natl Acad Sci USA</source><volume>94</volume><fpage>13493</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1073/pnas.94.25.13493</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brennecke</surname><given-names>J</given-names></name><name><surname>Malone</surname><given-names>CD</given-names></name><name><surname>Aravin</surname><given-names>AA</given-names></name><name><surname>Sachidanandam</surname><given-names>R</given-names></name><name><surname>Stark</surname><given-names>A</given-names></name><name><surname>Hannon</surname><given-names>GJ</given-names></name></person-group><year>2008</year><article-title>An epigenetic role for maternally inherited piRNAs in transposon silencing</article-title><source>Science</source><volume>322</volume><fpage>1387</fpage><lpage>92</lpage><pub-id pub-id-type="doi">10.1126/science.1165171</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chalmers</surname><given-names>RM</given-names></name><name><surname>Kleckner</surname><given-names>N</given-names></name></person-group><year>1994</year><article-title>Tn10/IS10 transposase purification, activation, and in vitro reaction</article-title><source>J Biol Chem</source><volume>269</volume><fpage>8029</fpage><lpage>35</lpage></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chalmers</surname><given-names>RM</given-names></name><name><surname>Kleckner</surname><given-names>N</given-names></name></person-group><year>1996</year><article-title>IS10/Tn10 transposition efficiently accommodates diverse transposon end configurations</article-title><source>EMBO J</source><volume>15</volume><fpage>5112</fpage><lpage>22</lpage></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Claeys Bouuaert</surname><given-names>C</given-names></name><name><surname>Chalmers</surname><given-names>R</given-names></name></person-group><year>2010</year><article-title>Transposition of the human Hsmar1 transposon: rate-limiting steps and the importance of the flanking TA dinucleotide in second strand cleavage</article-title><source>Nucleic Acids Res</source><volume>38</volume><fpage>190</fpage><lpage>202</lpage><pub-id pub-id-type="doi">10.1093/nar/gkp891</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Claeys Bouuaert</surname><given-names>C</given-names></name><name><surname>Liu</surname><given-names>D</given-names></name><name><surname>Chalmers</surname><given-names>R</given-names></name></person-group><year>2011</year><article-title>A simple topological filter in a eukaryotic transposon as a mechanism to suppress genome instability</article-title><source>Mol Cell Biol</source><volume>31</volume><fpage>317</fpage><lpage>27</lpage><pub-id pub-id-type="doi">10.1128/MCB.01066-10</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Clark</surname><given-names>KJ</given-names></name><name><surname>Carlson</surname><given-names>DF</given-names></name><name><surname>Leaver</surname><given-names>MJ</given-names></name><name><surname>Foster</surname><given-names>LK</given-names></name><name><surname>Fahrenkrug</surname><given-names>SC</given-names></name></person-group><year>2009</year><article-title>Passport, a native Tc1 transposon from flatfish, is functionally active in vertebrate cells</article-title><source>Nucleic Acids Res</source><volume>37</volume><fpage>1239</fpage><lpage>47</lpage><pub-id pub-id-type="doi">10.1093/nar/gkn1025</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cordaux</surname><given-names>R</given-names></name><name><surname>Udit</surname><given-names>S</given-names></name><name><surname>Batzer</surname><given-names>MA</given-names></name><name><surname>Feschotte</surname><given-names>C</given-names></name></person-group><year>2006</year><article-title>Birth of a chimeric primate gene by capture of the transposase gene from a mobile element</article-title><source>Proc Natl Acad Sci USA</source><volume>103</volume><fpage>8101</fpage><lpage>6</lpage><pub-id pub-id-type="doi">10.1073/pnas.0601161103</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dawson</surname><given-names>A</given-names></name><name><surname>Finnegan</surname><given-names>DJ</given-names></name></person-group><year>2003</year><article-title>Excision of the Drosophila mariner transposon Mos1. Comparison with bacterial transposition and V(D)J recombination</article-title><source>Mol Cell</source><volume>11</volume><fpage>225</fpage><lpage>35</lpage><pub-id pub-id-type="doi">10.1016/S1097-2765(02)00798-0</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Elf</surname><given-names>J</given-names></name><name><surname>Li</surname><given-names>GW</given-names></name><name><surname>Xie</surname><given-names>XS</given-names></name></person-group><year>2007</year><article-title>Probing transcription factor dynamics at the single-molecule level in a living cell</article-title><source>Science</source><volume>316</volume><fpage>1191</fpage><lpage>4</lpage><pub-id pub-id-type="doi">10.1126/science.1141967</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Engels</surname><given-names>WR</given-names></name></person-group><year>1983</year><article-title>The P family of transposable elements in Drosophila</article-title><source>Annu Rev Genet</source><volume>17</volume><fpage>315</fpage><lpage>44</lpage><pub-id pub-id-type="doi">10.1146/annurev.ge.17.120183.001531</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Feschotte</surname><given-names>C</given-names></name><name><surname>Wessler</surname><given-names>SR</given-names></name></person-group><year>2002</year><article-title>Mariner-like transposases are widespread and diverse in flowering plants</article-title><source>Proc Natl Acad Sci USA</source><volume>99</volume><fpage>280</fpage><lpage>5</lpage><pub-id pub-id-type="doi">10.1073/pnas.022626699</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fraser</surname><given-names>MJ</given-names></name><name><surname>Ciszczon</surname><given-names>T</given-names></name><name><surname>Elick</surname><given-names>T</given-names></name><name><surname>Bauser</surname><given-names>C</given-names></name></person-group><year>1996</year><article-title>Precise excision of TTAA-specific lepidopteran transposons piggyBac (IFP2) and tagalong (TFP3) from the baculovirus genome in cell lines from two species of Lepidoptera</article-title><source>Insect Mol Biol</source><volume>5</volume><fpage>141</fpage><lpage>51</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2583.1996.tb00048.x</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goodier</surname><given-names>JL</given-names></name><name><surname>Davidson</surname><given-names>WS</given-names></name></person-group><year>1994</year><article-title>Tc1 transposon-like sequences are widely distributed in salmonids</article-title><source>J Mol Biol</source><volume>241</volume><fpage>26</fpage><lpage>34</lpage><pub-id pub-id-type="doi">10.1006/jmbi.1994.1470</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grabundzija</surname><given-names>I</given-names></name><name><surname>Irgang</surname><given-names>M</given-names></name><name><surname>Mates</surname><given-names>L</given-names></name><name><surname>Belay</surname><given-names>E</given-names></name><name><surname>Matrai</surname><given-names>J</given-names></name><name><surname>Gogol-Doring</surname><given-names>A</given-names></name><etal/></person-group><year>2010</year><article-title>Comparative analysis of transposable element vector systems in human cells</article-title><source>Mol Ther</source><volume>18</volume><fpage>1200</fpage><lpage>9</lpage><pub-id pub-id-type="doi">10.1038/mt.2010.47</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hartl</surname><given-names>DL</given-names></name><name><surname>Lozovskaya</surname><given-names>ER</given-names></name><name><surname>Nurminsky</surname><given-names>DI</given-names></name><name><surname>Lohe</surname><given-names>AR</given-names></name></person-group><year>1997</year><article-title>What restricts the activity of mariner-like transposable elements</article-title><source>Trends Genet</source><volume>13</volume><fpage>197</fpage><lpage>201</lpage><pub-id pub-id-type="doi">10.1016/S0168-9525(97)01087-1</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hickman</surname><given-names>AB</given-names></name><name><surname>Chandler</surname><given-names>M</given-names></name><name><surname>Dyda</surname><given-names>F</given-names></name></person-group><year>2010</year><article-title>Integrating prokaryotes and eukaryotes: DNA transposases in light of structure</article-title><source>Crit Rev Biochem Mol Biol</source><volume>45</volume><fpage>50</fpage><lpage>69</lpage><pub-id pub-id-type="doi">10.3109/10409230903505596</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hsieh</surname><given-names>WT</given-names></name><name><surname>Whitson</surname><given-names>PA</given-names></name><name><surname>Matthews</surname><given-names>KS</given-names></name><name><surname>Wells</surname><given-names>RD</given-names></name></person-group><year>1987</year><article-title>Influence of sequence and distance between two operators on interaction with the lac repressor</article-title><source>J Biol Chem</source><volume>262</volume><fpage>14583</fpage><lpage>91</lpage></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hua-Van</surname><given-names>A</given-names></name><name><surname>Le Rouzic</surname><given-names>A</given-names></name><name><surname>Boutin</surname><given-names>TS</given-names></name><name><surname>Filee</surname><given-names>J</given-names></name><name><surname>Capy</surname><given-names>P</given-names></name></person-group><year>2011</year><article-title>The struggle for life of the genome’s selfish architects</article-title><source>Biol Direct</source><volume>6</volume><fpage>19</fpage><pub-id pub-id-type="doi">10.1186/1745-6150-6-19</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ivics</surname><given-names>Z</given-names></name><name><surname>Hackett</surname><given-names>PB</given-names></name><name><surname>Plasterk</surname><given-names>RH</given-names></name><name><surname>Izsvak</surname><given-names>Z</given-names></name></person-group><year>1997</year><article-title>Molecular reconstruction of sleeping beauty, a Tc1-like transposon from fish, and its transposition in human cells</article-title><source>Cell</source><volume>91</volume><fpage>501</fpage><lpage>10</lpage><pub-id>10.1016/S0092-8674(00)80436-5</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jensen</surname><given-names>PA</given-names></name><name><surname>Stuart</surname><given-names>JR</given-names></name><name><surname>Goodpaster</surname><given-names>MP</given-names></name><name><surname>Goodman</surname><given-names>JW</given-names></name><name><surname>Simmons</surname><given-names>MJ</given-names></name></person-group><year>2008</year><article-title>Cytotype regulation of P transposable elements in Drosophila melanogaster: repressor polypeptides or piRNAs?</article-title><source>Genetics</source><volume>179</volume><fpage>1785</fpage><lpage>93</lpage><pub-id pub-id-type="doi">10.1534/genetics.108.087072</pub-id></element-citation></ref><ref id="bib55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kelleher</surname><given-names>ES</given-names></name><name><surname>Barbash</surname><given-names>DA</given-names></name></person-group><year>2013</year><article-title>Analysis of piRNA-mediated silencing of active TEs in <italic>Drosophila melanogaster</italic> suggests limits on the evolution of host genome defense</article-title><source>Mol Biol Evol</source><pub-id pub-id-type="doi">10.1093/molbev/mst081</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kleckner</surname><given-names>N</given-names></name></person-group><year>1990</year><article-title>Regulating Tn10 and Is10 transposition</article-title><source>Genetics</source><volume>124</volume><fpage>449</fpage><lpage>54</lpage></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lampe</surname><given-names>DJ</given-names></name><name><surname>Grant</surname><given-names>TE</given-names></name><name><surname>Robertson</surname><given-names>HM</given-names></name></person-group><year>1998</year><article-title>Factors affecting transposition of the Himar1 mariner transposon in vitro</article-title><source>Genetics</source><volume>149</volume><fpage>179</fpage><lpage>87</lpage></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Le Rouzic</surname><given-names>A</given-names></name><name><surname>Boutin</surname><given-names>TS</given-names></name><name><surname>Capy</surname><given-names>P</given-names></name></person-group><year>2007</year><article-title>Long-term evolution of transposable elements</article-title><source>Proc Natl Acad Sci USA</source><volume>104</volume><fpage>19375</fpage><lpage>80</lpage><pub-id pub-id-type="doi">10.1073/pnas.0705238104</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Le Rouzic</surname><given-names>A</given-names></name><name><surname>Capy</surname><given-names>P</given-names></name></person-group><year>2005</year><article-title>The first steps of transposable elements invasion: parasitic strategy vs. genetic drift</article-title><source>Genetics</source><volume>169</volume><fpage>1033</fpage><lpage>43</lpage><pub-id pub-id-type="doi">10.1534/genetics.104.031211</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lipkow</surname><given-names>K</given-names></name><name><surname>Buisine</surname><given-names>N</given-names></name><name><surname>Lampe</surname><given-names>DJ</given-names></name><name><surname>Chalmers</surname><given-names>R</given-names></name></person-group><year>2004</year><article-title>Early intermediates of mariner transposition: catalysis without synapsis of the transposon ends suggests a novel architecture of the synaptic complex</article-title><source>Mol Cell Biol</source><volume>24</volume><fpage>8301</fpage><lpage>11</lpage><pub-id pub-id-type="doi">10.1128/MCB.24.18.8301-8311.2004</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname><given-names>D</given-names></name><name><surname>Bischerour</surname><given-names>J</given-names></name><name><surname>Siddique</surname><given-names>A</given-names></name><name><surname>Buisine</surname><given-names>N</given-names></name><name><surname>Bigot</surname><given-names>Y</given-names></name><name><surname>Chalmers</surname><given-names>R</given-names></name></person-group><year>2007</year><article-title>The human SETMAR protein preserves most of the activities of the ancestral Hsmar1 transposase</article-title><source>Mol Cell Biol</source><volume>27</volume><fpage>1125</fpage><lpage>32</lpage><pub-id pub-id-type="doi">10.1128/MCB.01899-06</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lohe</surname><given-names>AR</given-names></name><name><surname>Hartl</surname><given-names>DL</given-names></name></person-group><year>1996</year><article-title>Autoregulation of mariner transposase activity by overproduction and dominant-negative complementation</article-title><source>Mol Biol Evol</source><volume>13</volume><fpage>549</fpage><lpage>55</lpage><pub-id pub-id-type="doi">10.1093/oxfordjournals.molbev.a025615</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lohe</surname><given-names>AR</given-names></name><name><surname>Moriyama</surname><given-names>EN</given-names></name><name><surname>Lidholm</surname><given-names>DA</given-names></name><name><surname>Hartl</surname><given-names>DL</given-names></name></person-group><year>1995</year><article-title>Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements</article-title><source>Mol Biol Evol</source><volume>12</volume><fpage>62</fpage><lpage>72</lpage><pub-id pub-id-type="doi">10.1093/oxfordjournals.molbev.a040191</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Marshall</surname><given-names>WF</given-names></name><name><surname>Straight</surname><given-names>A</given-names></name><name><surname>Marko</surname><given-names>JF</given-names></name><name><surname>Swedlow</surname><given-names>J</given-names></name><name><surname>Dernburg</surname><given-names>A</given-names></name><name><surname>Belmont</surname><given-names>A</given-names></name><etal/></person-group><year>1997</year><article-title>Interphase chromosomes undergo constrained diffusional motion in living cells</article-title><source>Curr Biol</source><volume>7</volume><fpage>930</fpage><lpage>9</lpage><pub-id pub-id-type="doi">10.1016/S0960-9822(06)00412-X</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Miskey</surname><given-names>C</given-names></name><name><surname>Papp</surname><given-names>B</given-names></name><name><surname>Mates</surname><given-names>L</given-names></name><name><surname>Sinzelle</surname><given-names>L</given-names></name><name><surname>Keller</surname><given-names>H</given-names></name><name><surname>Izsvak</surname><given-names>Z</given-names></name><etal/></person-group><year>2007</year><article-title>The ancient mariner sails again: transposition of the human Hsmar1 element by a reconstructed transposase and activities of the SETMAR protein on transposon ends</article-title><source>Mol Cell Biol</source><volume>27</volume><fpage>4589</fpage><lpage>600</lpage><pub-id pub-id-type="doi">10.1128/MCB.02027-06</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Misra</surname><given-names>S</given-names></name><name><surname>Rio</surname><given-names>DC</given-names></name></person-group><year>1990</year><article-title>Cytotype control of Drosophila P element transposition: the 66 kd protein is a repressor of transposase activity</article-title><source>Cell</source><volume>62</volume><fpage>269</fpage><lpage>84</lpage><pub-id pub-id-type="doi">10.1016/0092-8674(90)90365-L</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Petrov</surname><given-names>EP</given-names></name><name><surname>Ohrt</surname><given-names>T</given-names></name><name><surname>Winkler</surname><given-names>RG</given-names></name><name><surname>Schwille</surname><given-names>P</given-names></name></person-group><year>2006</year><article-title>Diffusion and segmental dynamics of double-stranded DNA</article-title><source>Phys Rev Lett</source><volume>97</volume><fpage>258101</fpage><pub-id pub-id-type="doi">10.1103/PhysRevLett.97.258101</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Revzin</surname><given-names>A</given-names></name><name><surname>Von Hippel</surname><given-names>PH</given-names></name></person-group><year>1977</year><article-title>Direct measurement of association constants for the binding of <italic>Escherichia coli</italic> lac repressor to non-operator DNA</article-title><source>Biochemistry</source><volume>16</volume><fpage>4769</fpage><lpage>76</lpage><pub-id pub-id-type="doi">10.1021/bi00641a002</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reznikoff</surname><given-names>WS</given-names></name></person-group><year>2008</year><article-title>Transposon Tn5</article-title><source>Annu Rev Genet</source><volume>42</volume><fpage>269</fpage><lpage>86</lpage><pub-id pub-id-type="doi">10.1146/annurev.genet.42.110807.091656</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="book"><person-group person-group-type="editor"><name><surname>Rice</surname><given-names>SA</given-names></name></person-group><year>1985</year><source>Diffusion limited reactions</source><publisher-loc>Amsterdam</publisher-loc><publisher-name>Elsevier</publisher-name></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Richardson</surname><given-names>JM</given-names></name><name><surname>Colloms</surname><given-names>SD</given-names></name><name><surname>Finnegan</surname><given-names>DJ</given-names></name><name><surname>Walkinshaw</surname><given-names>MD</given-names></name></person-group><year>2009</year><article-title>Molecular architecture of the Mos1 paired-end complex: the structural basis of DNA transposition in a eukaryote</article-title><source>Cell</source><volume>138</volume><fpage>1096</fpage><lpage>108</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2009.07.012</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ringrose</surname><given-names>L</given-names></name><name><surname>Chabanis</surname><given-names>S</given-names></name><name><surname>Angrand</surname><given-names>PO</given-names></name><name><surname>Woodroofe</surname><given-names>C</given-names></name><name><surname>Stewart</surname><given-names>AF</given-names></name></person-group><year>1999</year><article-title>Quantitative comparison of DNA looping in vitro and in vivo: chromatin increases effective DNA flexibility at short distances</article-title><source>EMBO J</source><volume>18</volume><fpage>6630</fpage><lpage>41</lpage><pub-id pub-id-type="doi">10.1093/emboj/18.23.6630</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robert</surname><given-names>VJ</given-names></name><name><surname>Davis</surname><given-names>MW</given-names></name><name><surname>Jorgensen</surname><given-names>EM</given-names></name><name><surname>Bessereau</surname><given-names>JL</given-names></name></person-group><year>2008</year><article-title>Gene conversion and end-joining-repair double-strand breaks in the <italic>Caenorhabditis elegans</italic> germline</article-title><source>Genetics</source><volume>180</volume><fpage>673</fpage><lpage>9</lpage><pub-id pub-id-type="doi">10.1534/genetics.108.089698</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Roberts</surname><given-names>D</given-names></name><name><surname>Hoopes</surname><given-names>BC</given-names></name><name><surname>Mcclure</surname><given-names>WR</given-names></name><name><surname>Kleckner</surname><given-names>N</given-names></name></person-group><year>1985</year><article-title>IS10 transposition is regulated by DNA adenine methylation</article-title><source>Cell</source><volume>43</volume><fpage>117</fpage><lpage>30</lpage><pub-id pub-id-type="doi">10.1016/0092-8674(85)90017-0</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robertson</surname><given-names>HM</given-names></name><name><surname>Lampe</surname><given-names>DJ</given-names></name></person-group><year>1995</year><article-title>Distribution of transposable elements in arthropods</article-title><source>Annu Rev Entomol</source><volume>40</volume><fpage>333</fpage><lpage>57</lpage><pub-id pub-id-type="doi">10.1146/annurev.en.40.010195.002001</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robertson</surname><given-names>HM</given-names></name><name><surname>Zumpano</surname><given-names>KL</given-names></name></person-group><year>1997</year><article-title>Molecular evolution of an ancient mariner transposon, Hsmar1, in the human genome</article-title><source>Gene</source><volume>205</volume><fpage>203</fpage><lpage>17</lpage><pub-id pub-id-type="doi">10.1016/S0378-1119(97)00472-1</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Roman</surname><given-names>Y</given-names></name><name><surname>Oshige</surname><given-names>M</given-names></name><name><surname>Lee</surname><given-names>YJ</given-names></name><name><surname>Goodwin</surname><given-names>K</given-names></name><name><surname>Georgiadis</surname><given-names>MM</given-names></name><name><surname>Hromas</surname><given-names>RA</given-names></name><etal/></person-group><year>2007</year><article-title>Biochemical characterization of a SET and transposase fusion protein, Metnase: its DNA binding and DNA cleavage activity</article-title><source>Biochemistry</source><volume>46</volume><fpage>11369</fpage><lpage>76</lpage><pub-id pub-id-type="doi">10.1021/bi7005477</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schwanhausser</surname><given-names>B</given-names></name><name><surname>Busse</surname><given-names>D</given-names></name><name><surname>Li</surname><given-names>N</given-names></name><name><surname>Dittmar</surname><given-names>G</given-names></name><name><surname>Schuchhardt</surname><given-names>J</given-names></name><name><surname>Wolf</surname><given-names>J</given-names></name><etal/></person-group><year>2011</year><article-title>Global quantification of mammalian gene expression control</article-title><source>Nature</source><volume>473</volume><fpage>337</fpage><lpage>42</lpage><pub-id pub-id-type="doi">10.1038/nature10098</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Serrano</surname><given-names>L</given-names></name><name><surname>Liang</surname><given-names>L</given-names></name><name><surname>Chang</surname><given-names>Y</given-names></name><name><surname>Deng</surname><given-names>L</given-names></name><name><surname>Maulion</surname><given-names>C</given-names></name><name><surname>Nguyen</surname><given-names>S</given-names></name><etal/></person-group><year>2011</year><article-title>Homologous recombination conserves DNA sequence integrity throughout the cell cycle in embryonic stem cells</article-title><source>Stem Cells Dev</source><volume>20</volume><fpage>363</fpage><lpage>74</lpage><pub-id pub-id-type="doi">10.1089/scd.2010.0159</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Spradling</surname><given-names>AC</given-names></name><name><surname>Bellen</surname><given-names>HJ</given-names></name><name><surname>Hoskins</surname><given-names>RA</given-names></name></person-group><year>2011</year><article-title>Drosophila P elements preferentially transpose to replication origins</article-title><source>Proc Natl Acad Sci USA</source><volume>108</volume><fpage>15948</fpage><lpage>53</lpage><pub-id pub-id-type="doi">10.1073/pnas.1112960108</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tichy</surname><given-names>ED</given-names></name><name><surname>Pillai</surname><given-names>R</given-names></name><name><surname>Deng</surname><given-names>L</given-names></name><name><surname>Liang</surname><given-names>L</given-names></name><name><surname>Tischfield</surname><given-names>J</given-names></name><name><surname>Schwemberger</surname><given-names>SJ</given-names></name><etal/></person-group><year>2010</year><article-title>Mouse embryonic stem cells, but not somatic cells, predominantly use homologous recombination to repair double-strand DNA breaks</article-title><source>Stem Cells Dev</source><volume>19</volume><fpage>1699</fpage><lpage>711</lpage><pub-id pub-id-type="doi">10.1089/scd.2010.0058</pub-id></element-citation></ref><ref id="bib51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Townsend</surname><given-names>JP</given-names></name><name><surname>Hartl</surname><given-names>DL</given-names></name></person-group><year>2000</year><article-title>The kinetics of transposable element autoregulation</article-title><source>Genetica</source><volume>108</volume><fpage>229</fpage><lpage>37</lpage><pub-id pub-id-type="doi">10.1023/A:1004172703790</pub-id></element-citation></ref><ref id="bib52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Walisko</surname><given-names>O</given-names></name><name><surname>Izsvak</surname><given-names>Z</given-names></name><name><surname>Szabo</surname><given-names>K</given-names></name><name><surname>Kaufman</surname><given-names>CD</given-names></name><name><surname>Herold</surname><given-names>S</given-names></name><name><surname>Ivics</surname><given-names>Z</given-names></name></person-group><year>2006</year><article-title>Sleeping beauty transposase modulates cell-cycle progression through interaction with Miz-1</article-title><source>Proc Natl Acad Sci USA</source><volume>103</volume><fpage>4062</fpage><lpage>7</lpage><pub-id pub-id-type="doi">10.1073/pnas.0507683103</pub-id></element-citation></ref><ref id="bib53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Weinreich</surname><given-names>MD</given-names></name><name><surname>Gasch</surname><given-names>A</given-names></name><name><surname>Reznikoff</surname><given-names>WS</given-names></name></person-group><year>1994</year><article-title>Evidence that the cis preference of the Tn5 transposase is caused by nonproductive multimerization</article-title><source>Genes Dev</source><volume>8</volume><fpage>2363</fpage><lpage>74</lpage><pub-id pub-id-type="doi">10.1101/gad.8.19.2363</pub-id></element-citation></ref><ref id="bib54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>G</given-names></name><name><surname>Nagel</surname><given-names>DH</given-names></name><name><surname>Feschotte</surname><given-names>C</given-names></name><name><surname>Hancock</surname><given-names>CN</given-names></name><name><surname>Wessler</surname><given-names>SR</given-names></name></person-group><year>2009</year><article-title>Tuned for transposition: molecular determinants underlying the hyperactivity of a Stowaway MITE</article-title><source>Science</source><volume>325</volume><fpage>1391</fpage><lpage>4</lpage><pub-id pub-id-type="doi">10.1126/science.1175688</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.00668.017</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Ferrell</surname><given-names>James</given-names></name><role>Reviewing editor</role><aff><institution>Stanford University</institution>, <country>United States</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://www.elifesciences.org/the-journal/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for sending your work entitled “The autoregulation of DNA transposons in higher eukaryotes: insights from the human mariner element Hsmar1” for consideration at <italic>eLife</italic>. Your article has been favorably evaluated by a Senior editor and 2 reviewers, one of whom is a member of our Board of Reviewing Editors.</p><p>The Reviewing editor and the other reviewer discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>The fundamental problem at hand is how a transposable element manages to co-exist with its host without its intrinsic tendency for proliferation eventually rendering the host (and thus also the transposon) extinct. As the authors point out, a solution to this problem is known for certain prokaryotic transposons, but the elucidated mechanism does not apply to eukaryotes (because it requires cis-acting transposase). In the case of eukaryotic transposons, there have been only sporadic proposed models, which have not been analyzed or assessed either experimentally or by quantitative modeling (which, in every case is dependent on information regarding the basic biochemistry of the transposition reaction).</p><p>The present paper suggests an entirely new model for regulation of eukaryotic transposons, motivated by biochemical findings, with direct implications for the model from theoretical simulations in combination with both in vitro and in vivo analysis of transposition. Numerous extensions of the model, including features relating to the presence/accumulation/roles of non-autonomous elements, are also considered.</p><p>The proposed model is very interesting and very sensible, not only a priori in a qualitative way but also quantitatively, as illustrated by modeling. Furthermore, taken together, the presented observations provide very strong support for the validity of the model. The experimental observations are important and compelling, particularly the in vitro findings relating the responses to variations in transposase and transposon end concentrations, and the in vivo verification of several key predictions. No such study has ever been performed for any transposon, even in the prokaryotic cases, and the results presented here should be of broad interest to both the evolutionary biology community and the transposition mechanism community.</p><p>Major comments:</p><p>1) <xref ref-type="fig" rid="fig1">Figure 1</xref> and related text: We think the reader might appreciate an explanation of why the kinetic schemes shown in <xref ref-type="fig" rid="fig1">Figures 1A, G, and H</xref> result in duplication of the transposon. If we understand correctly, the idea is that the excision of the transposon occurs after the transposon has replicated. The site from which it was excised is repaired by homologous recombination using the other chromosome as a template, regenerating the transposon at the excision site, and the excised transposon goes on to integrate somewhere else. Depending upon whether it integrates before or after the replication fork, or fails to integrate at all, the number of transposons in the genome can increase by different amounts or decrease.</p><p>In any case, it would certainly be helpful if the authors were to explain their thinking here, and to clarify the relationships between the kinetic schemes in <xref ref-type="fig" rid="fig1">Figure 1</xref> and the discussion in <xref ref-type="fig" rid="fig3">Figure 3GH</xref>.</p><p>2) <xref ref-type="fig" rid="fig7">Figure 7</xref> and related text: we stumbled over this section. The writing makes it seem as though a previously-proposed dimerization end-occlusion model also accounts for stable transposon levels.</p><p>Assuming we have understood correctly, we think the situation is as follows: the first paragraph discusses a case with a secondary binding site. It is likely that this site was proposed to have the effect of repressing transposition. The important point, which is implicit but should be stated explicitly, is that this feature does not in any way prevent the accumulation over time, which is the problem at hand here – it simply slows things down, regardless of whether that is in the current mechanism or the S-PD mechanism. That is, you still need the current mechanism (and the S-PD mechanism doesn’t work) if this feature is present.</p><p>The second paragraph addresses the DEO model for Tn<italic>5</italic> and mariner. Again, the writing does not make clear what the authors wish to say. The case of Tn<italic>5</italic> is not relevant here because it involves a cis-acting transposase. But the case of mariner obviously is relevant. Here there are two points. First, the DEO model is fine, but the biochemistry presented in this paper for mariner suggests that it is not what is going on. This point should be made explicitly. Second, in principle, an element with a <italic>different</italic> biochemistry might use this mechanism; but, if so, the modeling suggests that there is a longer lag time before the effect kicks in. The authors are perhaps trying to be sure to be generous to the previous considerations. But they need to say (a) DEO can’t be true for mariner because the biochemistry says it isn’t and (b) it could potentially be true for some other element but with a disadvantage. Are there further candidate elements that might work in such a way?</p><p>3) This general class of mechanism is formally analogous to the “prozone effect”, which is the reason why complexes formed between a bivalent antibody and a bivalent antigen tend to be dimeric when either antibody or antigen is in excess, and highly polymeric when they are balanced. This mechanism pops up again and again in protein regulation. See, for example Bray E, Lay S, PNAS 1997.</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.00668.018</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p><italic>1) <xref ref-type="fig" rid="fig1">Figure 1</xref> and related text: We think the reader might appreciate an explanation of why the kinetic schemes shown in <xref ref-type="fig" rid="fig1">Figures 1A, G, and H</xref> result in duplication of the transposon. If we understand correctly, the idea is that the excision of the transposon occurs after the transposon has replicated. The site from which it was excised is repaired by homologous recombination using the other chromosome as a template, regenerating the transposon at the excision site, and the excised transposon goes on to integrate somewhere else. Depending upon whether it integrates before or after the replication fork, or fails to integrate at all, the number of transposons in the genome can increase by different amounts or decrease</italic>.</p><p><italic>In any case, it would certainly be helpful if the authors were to explain their thinking here, and to clarify the relationships between the kinetic schemes in <xref ref-type="fig" rid="fig1">Figure 1</xref> and the discussion in <xref ref-type="fig" rid="fig3">Figure 3GH</xref></italic>.</p><p>Agreed: we have now restructured the first two paragraphs of the Results section. The issues of replication and repair are now covered at the end of the first paragraph.</p><p><italic>2) <xref ref-type="fig" rid="fig7">Figure 7</xref> and related text: we stumbled over this section. The writing makes it seem as though a previously-proposed dimerization end-occlusion model also accounts for stable transposon levels</italic>.</p><p><italic>Assuming we have understood correctly, we think the situation is as follows: the first paragraph discusses a case with a secondary binding site. It is likely that this site was proposed to have the effect of repressing transposition. The important point, which is implicit but should be stated explicitly, is that this feature does not in any way prevent the accumulation over time, which is the problem at hand here – it simply slows things down, regardless of whether that is in the current mechanism or the S-PD mechanism. That is, you still need the current mechanism (and the S-PD mechanism doesn’t work) if this feature is present</italic>.</p><p>Agreed: the referees have understood correctly. We have substantially rewritten the first paragraph of this section as two new paragraphs, which make the suggested changes.</p><p><italic>The second paragraph addresses the DEO model for Tn</italic>5 <italic>and mariner. Again, the writing does not make clear what the authors wish to say. The case of Tn</italic>5 <italic>is not relevant here because it involves a cis-acting transposase. But the case of mariner obviously is relevant. Here there are two points. First, the DEO model is fine, but the biochemistry presented in this paper for mariner suggests that it is not what is going on. This point should be made explicitly. Second, in principle, an element with a</italic> different <italic>biochemistry might use this mechanism; but, if so, the modeling suggests that there is a longer lag time before the effect kicks in. The authors are perhaps trying to be sure to be generous to the previous considerations. But they need to say (a) DEO can’t be true for mariner because the biochemistry says it isn’t and (b) it could potentially be true for some other element but with a disadvantage. Are there further candidate elements that might work in such a way?</italic></p><p>Yes, thank you. We always like to try to be as generous as possible. However, we have now rewritten and restructured the last two paragraphs of the results to make the explicit statements suggested by the referees.</p><p><italic>3) This general class of mechanism is formally analogous to the “prozone effect”, which is the reason why complexes formed between a bivalent antibody and a bivalent antigen tend to be dimeric when either antibody or antigen is in excess, and highly polymeric when they are balanced. This mechanism pops up again and again in protein regulation. See, for example Bray E, Lay S, PNAS 1997</italic>.</p><p>Thank you for highlighting this excellent example of the less-is-more effect. We were familiar with it from the Ouchterlony double diffusion experiments of our youth. However, we had quite overlooked its connection to our current work. In response we have deleted the last sentence of the Introduction and added a sentence to the end of the first paragraph of the Discussion, in which we cite the Bray and Lay paper.</p></body></sub-article></article> |