Logo: to the web site of Uppsala University

Many proteins that bind specific DNA sequences search the genome by combining three-dimensional diffusion with one-dimensional sliding on nonspecific DNA(1-5). Here we combine resonance energy transfer and fluorescence correlation measurements to characterize how individual lac repressor (LacI) molecules explore the DNA surface during the one-dimensional phase of target search. To track the rotation of sliding LacI molecules on the microsecond timescale, we use real-time single-molecule confocal laser tracking combined with fluorescence correlation spectroscopy (SMCT-FCS). The fluctuations in fluorescence signal are accurately described by rotation-coupled sliding, in which LacI traverses about 40 base pairs (bp) per revolution. This distance substantially exceeds the 10.5-bp helical pitch of DNA; this suggests that the sliding protein frequently hops out of the DNA groove, which would result in the frequent bypassing of target sequences. We directly observe such bypassing using single-molecule fluorescence resonance energy transfer (smFRET). A combined analysis of the smFRET and SMCT-FCS data shows that LacI hops one or two grooves (10-20 bp) every 200-700 mu s. Our data suggest a trade-off between speed and accuracy during sliding: the weak nature of nonspecific protein-DNA interactions underlies operator bypassing, but also speeds up sliding. We anticipate that SMCT-FCS, which monitors rotational diffusion on the microsecond timescale while tracking individual molecules with millisecond resolution, will be applicable to the real-time investigation of many other biological interactions and will effectively extend the accessible time regime for observing these interactions by two orders of magnitude. Single-molecule fluorescence resonance energy transfer and real-time confocal laser tracking with fluorescence correlation spectroscopy together characterize how individual lac repressor molecules bypass operator sites while exploring the DNA surface at microsecond timescales.

It has become increasingly clear that low levels of antibiotics present in many environments can select for resistant bacteria, yet the evolutionary pathways for resistance development during exposure to low amounts of antibiotics remain poorly defined. Here we show that Salmonella enterica exposed to sub-MIC levels of streptomycin evolved high-level resistance via novel mechanisms that are different from those observed during lethal selections. During lethal selection only rpsL mutations are found, whereas at sub-MIC selection resistance is generated by several small-effect resistance mutations that combined confer high-level resistance via three different mechanisms: (i) alteration of the ribosomal RNA target (gidB mutations), (ii) reduction in aminoglycoside uptake (cyoB, nuoG, and trkH mutations), and (iii) induction of the aminoglycoside-modifying enzyme AadA (znuA mutations). These results demonstrate how the strength of the selective pressure influences evolutionary trajectories and that even weak selective pressures can cause evolution of high-level resistance.

Random mutations in genes from disparate protein classes may have different distributions of fitness effects (DFEs) depending on different structural, functional, and evolutionary constraints. We measured the fitness effects of 156 single mutations in the genes encoding AraC (transcription factor), AraD (enzyme), and AraE (transporter) used for bacterial growth on L-arabinose. Despite their different molecular functions these genes all had bimodal DFEs with most mutations either being neutral or strongly deleterious, providing a general expectation for the DFE. This contrasts with the unimodal DFEs previously obtained for ribosomal protein genes where most mutations were slightly deleterious. Based on theoretical considerations, we suggest that the 33-fold higher average mutational robustness of ribosomal proteins is due to stronger selection for reduced costs of translational and transcriptional errors. Whereas the large majority of synonymous mutations were deleterious for ribosomal proteins genes, no fitness effects could be detected for the AraCDE genes. Four mutations in AraC and AraE increased fitness, suggesting that slightly advantageous mutations make up a significant fraction of the DFE, but that they often escape detection due to the limited sensitivity of commonly used fitness assays. We show that the fitness effects of amino acid substitutions can be predicted based on evolutionary conservation, but those weakly deleterious mutations are less reliably detected. This suggests that large-effect mutations and the fraction of highly deleterious mutations can be computationally predicted, but that experiments are required to characterize the DFE close to neutrality, where many mutations ultimately fixed in a population will occur.

The apparent dissociation constant (K-d) for specific binding of glucocorticoid receptor (GR) and androgen receptor (AR) to DNA was determined in vivo in Xenopus oocytes. The total nuclear receptor concentration was quantified as specifically retained [H-3]-hormone in manually isolated oocyte nuclei. DNA was introduced by nuclear microinjection of single stranded phagemid DNA, chromatin is then formed during second strand synthesis. The fraction of DNA sites occupied by the expressed receptor was determined by dimethylsulphate in vivo footprinting and used for calculation of the receptor-DNA binding affinity. The forkhead transcription factor FoxA1 enhanced the DNA binding by GR with an apparent K-d of similar to 1 mu M and dramatically stimulated DNA binding by AR with an apparent K-d of similar to 0.13 mu M at a composite androgen responsive DNA element containing one FoxA1 binding site and one palindromic hormone receptor binding site known to bind one receptor homodimer. FoxA1 exerted a weak constitutive- and strongly cooperative DNA binding together with AR but had a less prominent effect with GR, the difference reflecting the licensing function of FoxA1 at this androgen responsive DNA element.

We have investigated which aspects of transcription factor DNA interactions are most important to account for the recent in vivo search time measurements for the dimeric lac repressor. We find the best agreement for a sliding model where non-specific binding to DNA is improbable at first contact and the sliding LacI protein binds at high probability when reaching the specific O-sym operator. We also find that the contribution of hopping to the overall search speed is negligible although physically unavoidable. The parameters that give the best fit reveal sliding distances, including hopping, close to what has been proposed in the past, i.e. similar to 40 bp, but with an unexpectedly high 1D diffusion constant on non-specific DNA sequences. Including a mechanism of inter-segment transfer between distant DNA segments does not bring down the 1D diffusion to the expected fraction of the in vitro value. This suggests a mechanism where transcription factors can slide less hindered in vivo than what is given by a simple viscosity scaling argument or that a modification of the model is needed. For example, the estimated diffusion rate constant would be consistent with the expectation if parts of the chromosome, away from the operator site, were inaccessible for searching.

We have investigated which aspects of transcription factor DNA interactions are most important to account for the recent in vivo search time measurements for the dimeric lac repressor. We find the best agreement for a sliding model where non-specific binding to DNA is improbable at first contact and the sliding LacI protein binds at high probability when reaching the specific O_sym operator. We also find that the contribution of hopping to the overall search speed is negligible although physically unavoidable. The parameters that give the best fit reveal sliding distances, including hopping, close to what has been proposed in the past, i.e. ∼40 bp, but with an unexpectedly high 1D diffusion constant on non-specific DNA sequences. Including a mechanism of inter-segment transfer between distant DNA segments does not bring down the 1D diffusion to the expected fraction of the in vitro value. This suggests a mechanism where transcription factors can slide less hindered in vivo than what is given by a simple viscosity scaling argument or that a modification of the model is needed. For example, the estimated diffusion rate constant would be consistent with the expectation if parts of the chromosome, away from the operator site, were inaccessible for searching.

An important mechanism for generation of new genes is by duplication-divergence of existing genes. Duplication-divergence includes several different sub-models, such as subfunctionalization where after accumulation of neutral mutations the original function is distributed between two partially functional and complementary genes, and neofunctionalization where a new function evolves in one of the duplicated copies while the old function is maintained in another copy. The likelihood of these mechanisms depends on the longevity of the duplicated state, which in turn depends on the fitness cost and genetic stability of the duplications. Here, we determined the fitness cost and stability of defined gene duplications/amplifications on a low copy number plasmid. Our experimental results show that the costs of carrying extra gene copies are substantial and that each additional kbp of DNA reduces fitness by approximately 0.15%. Furthermore, gene amplifications are highly unstable and rapidly segregate to lower copy numbers in absence of selection. Mathematical modelling shows that the fitness costs and instability strongly reduces the likelihood of both sub- and neofunctionalization, but that these effects can be off-set by positive selection for novel beneficial functions.

Transcription factors search for specific operator sequences by alternating rounds of 3D diffusion with rounds of 1D diffusion (sliding) along the DNA. The details of such sliding have largely been beyond direct experimental observation. For this purpose we devised an analytical formulation of umbrella sampling along a helical coordinate, and from extensive and fully atomistic simulations we quantified the free-energy landscapes that underlie the sliding dynamics and dissociation kinetics for the LacI dimer. The resulting potential of mean force distributions show a fine structure with an amplitude of 1 k(B)T for sliding and 12 kBT for dissociation. Based on the free-energy calculations the repressor slides in close contact with DNA for 8 bp on average before making a microscopic dissociation. By combining the microscopic molecular-dynamics calculations with Brownian simulation including rotational diffusion from the microscopically dissociated state we estimate a macroscopic residence time of 48 ms at the same DNA segment and an in vitro sliding distance of 240 bp. The sliding distance is in agreement with previous in vitro sliding-length estimates. The in vitro prediction for the macroscopic residence time also compares favorably to what we measure by single-molecule imaging of nonspecifically bound fluorescently labeled LacI in living cells. The investigation adds to our understanding of transcription-factor search kinetics and connects the macro-/mesoscopic rate constants to the microscopic dynamics.

The classical Meselson-Stahl density-shift method was used to study replication of pOU71, a runaway-replication derivative of plasmid R1 in Escherichia coli. The miniplasmid maintained the normal low copy number of R1 during steady growth at 30 degrees C, but as growth temperatures were raised above 34 degrees C, the copy number of the plasmid increased to higher levels, and at 42 degrees C, it replicated without control in a runaway replication mode with lethal consequences for the host. The eclipse periods (minimum time between successive replication of the same DNA) of the plasmid shortened with rising copy numbers at increasing growth temperatures (Olsson et al., 2003). In this work, eclipse periods were measured during downshifts in copy number of pOU71 after it had replicated at 39 and 42 degrees C, resulting in 7- and 50-fold higher than normal plasmid copy number per cell, respectively. Eclipse periods for plasmid replication, measured during copy number downshift, suggested that plasmid R1, normally selected randomly for replication, showed a bias such that a newly replicated DNA had a higher probability of replication compared to the bulk of the RI population. However, even the unexpected nonrandom replication followed the copy number kinetics such that every generation, the plasmids underwent the normal inherited number of replication, n, independent of the actual number of plasmid copies in a newborn cell.

Physical modeling is increasingly important for generating insights into intracellular processes. We describe situations in which combined spatial and stochastic aspects of chemical reactions are needed to capture the relevant dynamics of biochemical systems.