What matters for lac repressor search in vivo--sliding, hopping, intersegment transfer, crowding on DNA or recognition?

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 Osym 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.

Transcription-factor binding and sliding on DNA studied using micro- and macroscopic models.

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 k(B)T 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.

Lost in presumption: stochastic reactions in spatial models.

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.

MesoRD 1.0: Stochastic reaction-diffusion simulations in the microscopic limit.

SUMMARY: MesoRD is a tool for simulating stochastic reaction-diffusion systems as modeled by the reaction diffusion master equation. The simulated systems are defined in the Systems Biology Markup Language with additions to define compartment geometries. MesoRD 1.0 supports scale-dependent reaction rate constants and reactions between reactants in neighbouring subvolumes. These new features make it possible to construct physically consistent models of diffusion-controlled reactions also at fine spatial discretization. AVAILABILITY: MesoRD is written in C++ and licensed under the GNU general public license (GPL). MesoRD can be downloaded at http://mesord.sourceforge.net. The MesoRD homepage, http://mesord.sourceforge.net, contains detailed documentation and news about recently implemented features. CONTACT: johan.elf@icm.uu.se.

The Infectious Dose Shapes Vibrio cholerae Within-Host Dynamics.

During infection, the rates of pathogen replication, death, and migration affect disease progression, dissemination, transmission, and resistance evolution. Here, we follow the population dynamics of Vibrio cholerae in a mouse model by labeling individual bacteria with one of >500 unique, fitness-neutral genomic tags. Using the changes in tag frequencies and CFU numbers, we inform a mathematical model that describes the within-host spatiotemporal bacterial dynamics. This allows us to disentangle growth, death, forward, and retrograde migration rates continuously during infection. Our model has robust predictive power across various experimental setups. The population dynamics of V. cholerae shows substantial spatiotemporal heterogeneity in replication, death, and migration. Importantly, we find that the niche available to V. cholerae in the host increases with inoculum size, suggesting cooperative effects during infection. Therefore, it is not enough to consider just the likelihood of exposure (50% infectious dose) but rather the magnitude of exposure to predict outbreaks. IMPORTANCE Determining the rates of bacterial migration, replication, and death during infection is important for understanding how infections progress. Separately measuring these rates is often difficult in systems where multiple processes happen simultaneously. Here, we use next-generation sequencing to measure V. cholerae migration, replication, death, and niche size along the mouse gastrointestinal tract. We show that the small intestine of the mouse is a heterogeneous environment, and the population dynamic characteristics change substantially between adjacent gut sections. Our approach also allows us to characterize the effect of inoculum size on these processes. We find that the niche size in mice increases with the infectious dose, hinting at cooperative effects in larger inocula. The dose-response relationship between inoculum size and final pathogen burden is important for the infected individual and is thought to influence the progression of V. cholerae epidemics.

RESTAMP - Rate estimates by sequence-tag analysis of microbial populations.

Microbial division rates determine the speed of mutation accumulation and thus the emergence of antimicrobial resistance. Microbial death rates are affected by antibiotic action and the immune system. Therefore, measuring these rates has advanced our understanding of host-pathogen interactions and antibiotic action. Several methods based on marker-loss or few inheritable neutral markers exist that allow estimating microbial division and death rates, each of which has advantages and limitations. Technical bottlenecks, i.e., experimental sampling events, during the experiment can distort the rate estimates and are typically unaccounted for or require additional calibration experiments. In this work, we introduce RESTAMP (Rate Estimates by Sequence Tag Analysis of Microbial Populations) as a method for determining bacterial division and death rates. This method uses hundreds of fitness neutral sequence barcodes to measure the rates and account for experimental bottlenecks at the same time. We experimentally validate RESTAMP and compare it to established plasmid loss methods. We find that RESTAMP has a number of advantages over plasmid loss or previous marker based techniques. (i) It enables to correct the distortion of rate estimates by technical bottlenecks. (ii) Rate estimates are independent of the sequence tag distribution in the starting culture allowing the use of an arbitrary number of tags. (iii) It introduces a bottleneck sensitivity measure that can be used to maximize the accuracy of the experiment. RESTAMP allows studying microbial population dynamics with great resolution over a wide dynamic range and can thus advance our understanding of host-pathogen interactions or the mechanisms of antibiotic action.

Chest x-ray resolution after SARS-CoV-2 infection.

Aim To analyse the resolution of chest X-ray findings in relation to laboratory parameters in patients infected with acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in a two- month followup. Analysis of chest X-ray findings in the first few months after the disease is the main goal of our work. Methods Out of the total of 343 patients chest X-ray findings were followed in 269 patients. Patients were divided into groups according to the severity of findings. D-dimer, inflammatory markers, blood cell count, neutrophil lymphocyte ratio (NLR) were analysed. Chest X-ray was analysed during the hospitalization on the day of admission, on the third, the seventh and the fourteenth day (scoring method was used). After discharge chest X-ray was performed in a two-week follow-up, then after one and two months, and after three months if necessary. Results Incomplete chest X-ray resolution was identified in 24 (39.34%) patients with severe, 27 (22.31 %) patients with moderate and in three (3.91%) patients with mild findings. Statistical significance was established in overall score by comparison between all groups (p<0.001), and in the moderate compared to the mild group (p=0.0051). The difference of NLR in the severe compared to the moderate group was observed (p=0.0021) and in the severe group compared to the mild group (p=0.00013). Conclusion Chest X-ray findings persisted mostly in the severe group followed by the moderate and mild ones. Long-term followup is necessary for the appropriate treatment and prevention of fibrosis, and reduction of symptoms.

Selection or drift: The population biology underlying transposon insertion sequencing experiments.

Transposon insertion sequencing methods such as Tn-seq revolutionized microbiology by allowing the identification of genomic loci that are critical for viability in a specific environment on a genome-wide scale. While powerful, transposon insertion sequencing suffers from limited reproducibility when different analysis methods are compared. From the perspective of population biology, this may be explained by changes in mutant frequency due to chance (drift) rather than differential fitness (selection). Here, we develop a mathematical model of the population biology of transposon insertion sequencing experiments, i.e. the changes in size and composition of the transposon-mutagenized population during the experiment. We use this model to investigate mutagenesis, the growth of the mutant library, and its passage through bottlenecks. Specifically, we study how these processes can lead to extinction of individual mutants depending on their fitness and the distribution of fitness effects (DFE) of the entire mutant population. We find that in typical in vitro experiments few mutants with high fitness go extinct. However, bottlenecks of a size that is common in animal infection models lead to so much random extinction that a large number of viable mutants would be misclassified. While mutants with low fitness are more likely to be lost during the experiment, mutants with intermediate fitness are expected to be much more abundant and can constitute a large proportion of detected hits, i.e. false positives. Thus, incorporating the DFEs of randomly generated mutations in the analysis may improve the reproducibility of transposon insertion experiments, especially when strong bottlenecks are encountered.

The lac repressor displays facilitated diffusion in living cells.

Transcription factors (TFs) are proteins that regulate the expression of genes by binding sequence-specific sites on the chromosome. It has been proposed that to find these sites fast and accurately, TFs combine one-dimensional (1D) sliding on DNA with 3D diffusion in the cytoplasm. This facilitated diffusion mechanism has been demonstrated in vitro, but it has not been shown experimentally to be exploited in living cells. We have developed a single-molecule assay that allows us to investigate the sliding process in living bacteria. Here we show that the lac repressor slides 45 ± 10 base pairs on chromosomal DNA and that sliding can be obstructed by other DNA-bound proteins near the operator. Furthermore, the repressor frequently (>90%) slides over its natural lacO(1) operator several times before binding. This suggests a trade-off between rapid search on nonspecific sequences and fast binding at the specific sequence.