Kinetic pathways of topology simplification by Type-II topoisomerases in knotted supercoiled DNA

The topological state of covalently closed, double-stranded DNA is defined by the knot type $K$ and the linking-number difference $\\Delta Lk$ relative to unknotted relaxed DNA. DNA topoisomerases are essential enzymes that control the topology of DNA in all cells. In particular, type-II topoisomerases change both $K$ and $\\Delta Lk$ by a duplex-strand-passage mechanism and have been shown to simplify the topology of DNA to levels below thermal equilibrium at the expense of ATP hydrolysis. It remains a key question how small enzymes are able to preferentially select strand passages that result in topology simplification in much larger DNA molecules. Using numerical simulations, we consider the non-equilibrium dynamics of transitions between topological states $(K,\\Delta Lk)$ in DNA induced by type-II topoisomerases. For a biological process that delivers DNA molecules in a given topological state $(K,\\Delta Lk)$ at a constant rate we fully characterize the pathways of topology simplification by type-II topoisomerases in terms of stationary probability distributions and probability currents on the network of topological states $(K,\\Delta Lk)$. In particular, we observe that type-II topoisomerase activity is significantly enhanced in DNA molecules that maintain a supercoiled state with constant torsional tension. This is relevant for bacterial cells in which torsional tension is maintained by enzyme-dependent homeostatic mechanisms such as DNA-gyrase activity

DNA-Topology Simplification by Topoisomerases

The topological properties of DNA molecules, supercoiling, knotting, and catenation, are intimately connected with essential biological processes, such as gene expression, replication, recombination, and chromosome segregation. Non-trivial DNA topologies present challenges to the molecular machines that process and maintain genomic information, for example, by creating unwanted DNA entanglements. At the same time, topological distortion can facilitate DNA-sequence recognition through localized duplex unwinding and longer-range loop-mediated interactions between the DNA sequences. Topoisomerases are a special class of essential enzymes that homeostatically manage DNA topology through the passage of DNA strands. The activities of these enzymes are generally investigated using circular DNA as a model system, in which case it is possible to directly assay the formation and relaxation of DNA supercoils and the formation/resolution of knots and catenanes. Some topoisomerases use ATP as an energy cofactor, whereas others act in an ATP-independent manner. The free energy of ATP hydrolysis can be used to drive negative and positive supercoiling or to specifically relax DNA topologies to levels below those that are expected at thermodynamic equilibrium. The latter activity, which is known as topology simplification, is thus far exclusively associated with type-II topoisomerases and it can be understood through insight into the detailed non-equilibrium behavior of type-II enzymes. We use a non-equilibrium topologicalnetwork approach, which stands in contrast to the equilibrium models that are conventionally used in the DNA-topology field, to gain insights into the rates that govern individual transitions between topological states. We anticipate that our quantitative approach will stimulate experimental work and the theoretical/computational modeling of topoisomerases and similar enzyme systems

DNA cyclization and looping in the wormlike limit: normal modes and the validity of the harmonic approximation

For much of the last three decades, Monte Carlo-simulation methods have been the standard approach for accurately calculating the cyclization probability, J, or J factor, for DNA models having sequence-dependent bends or inhomogeneous bending flexibility. Within the last 10 years approaches based on harmonic analysis of semi-flexible polymer models have been introduced, which offer much greater computational efficiency than Monte Carlo techniques. These methods consider the ensemble of molecular conformations in terms of harmonic fluctuations about a well-defined elastic-energy minimum. However, the harmonic approximation is only applicable for small systems, because the accessible conformation space of larger systems is increasingly dominated by anharmonic contributions. In the case of computed values of the J factor, deviations of the harmonic approximation from the exact value of J as a function of DNA length have not been characterized. Using a recent, numerically exact method that accounts for both anharmonic and harmonic contributions to J for wormlike chains of arbitrary size, we report here the apparent error that results from neglecting anharmonic behavior. For wormlike chains having contour lengths less than four times the persistence length, the error in J arising from the harmonic approximation is generally small, amounting to free energies less than the thermal energy, kB T. For larger systems, however, the deviations between harmonic and exact J values increase approximately linearly with size

The thermodynamics of DNA loop formation, from J to Z

The formation of DNA loops is a ubiquitous theme in biological processes, including DNA replication, recombination and repair, and gene regulation. These loops are mediated by proteins bound at specific sites along the contour of a single DNA molecule, in some cases many thousands of base pairs apart. Loop formation incurs a thermodynamic cost that is a sensitive function of the length of looped DNA as well as the geometry and elastic properties of the DNA-bound protein. The free energy of DNA looping is logarithmically related to a generalization of the Jacobson-Stockmayer factor for DNA cyclization, termed the J factor. In the present article, we review the thermodynamic origins of this quantity, discuss how it is measured experimentally and connect the macroscopic interpretation of the J factor with a statistical-mechanical description of DNA looping and cyclization

Free-energy calculations for semi-flexible macromolecules: Applications to DNA knotting and looping

We present a method to obtain numerically accurate values of configurational free energies of semiflexible macromolecular systems, based on the technique of thermodynamic integration combined with normal-mode analysis of a reference system subject to harmonic constraints. Compared with previous free-energy calculations that depend on a reference state, our approach introduces two innovations, namely the use of internal coordinates to constrain the reference states and the ability to freely select these reference states. As a consequence, it is possible to explore systems that undergo substantially larger fluctuations than those considered in previous calculations, including semiflexible biopolymers having arbitrary ratios of contour length L to persistence length P. To validate the method, high accuracy is demonstrated for free energies of prime DNA knots with L/P=20 and L/P=40, corresponding to DNA lengths of 3000 and 6000 base pairs, respectively. We then apply the method to study the free-energy landscape for a model of a synaptic nucleoprotein complex containing a pair of looped domains, revealing a bifurcation in the location of optimal synapse (crossover) sites. This transition is relevant to target-site selection by DNA-binding proteins that occupy multiple DNA sites separated by large linear distances along the genome, a problem that arises naturally in gene regulation, DNA recombination, and the action of type-II topoisomerases.Comment: Main article: 24 pages, 8 figures; Supplemental Material: 21 pages, 6 figures. Typos correcte

Analysis of In-Vivo LacR-Mediated Gene Repression Based on the Mechanics of DNA Looping

Interactions of E. coli lac repressor (LacR) with a pair of operator sites on the same DNA molecule can lead to the formation of looped nucleoprotein complexes both in vitro and in vivo. As a major paradigm for loop-mediated gene regulation, parameters such as operator affinity and spacing, repressor concentration, and DNA bending induced by specific or non-specific DNA-binding proteins (e.g., HU), have been examined extensively. However, a complete and rigorous model that integrates all of these aspects in a systematic and quantitative treatment of experimental data has not been available. Applying our recent statistical-mechanical theory for DNA looping, we calculated repression as a function of operator spacing (58–156 bp) from first principles and obtained excellent agreement with independent sets of in-vivo data. The results suggest that a linear extended, as opposed to a closed v-shaped, LacR conformation is the dominant form of the tetramer in vivo. Moreover, loop-mediated repression in wild-type E. coli strains is facilitated by decreased DNA rigidity and high levels of flexibility in the LacR tetramer. In contrast, repression data for strains lacking HU gave a near-normal value of the DNA persistence length. These findings underscore the importance of both protein conformation and elasticity in the formation of small DNA loops widely observed in vivo, and demonstrate the utility of quantitatively analyzing gene regulation based on the mechanics of nucleoprotein complexes

Fractal Dimension Analysis of Transient Visual Evoked Potentials: Optimisation and Applications

Purpose The visual evoked potential (VEP) provides a time series signal response to an external visual stimulus at the location of the visual cortex. The major VEP signal components, peak latency and amplitude, may be affected by disease processes. Additionally, the VEP contains fine detailed and non-periodic structure, of presently unclear relevance to normal function, which may be quantified using the fractal dimension. The purpose of this study is to provide a systematic investigation of the key parameters in the measurement of the fractal dimension of VEPs, to develop an optimal analysis protocol for application. Methods VEP time series were mathematically transformed using delay time, τ, and embedding dimension, m, parameters. The fractal dimension of the transformed data was obtained from a scaling analysis based on straight line fits to the numbers of pairs of points with separation less than r versus log(r) in the transformed space. Optimal τ, m, and scaling analysis were obtained by comparing the consistency of results using different sampling frequencies. The optimised method was then piloted on samples of normal and abnormal VEPs. Results Consistent fractal dimension estimates were obtained using τ = 4 ms, designating the fractal dimension = D2 of the time series based on embedding dimension m = 7 (for 3606 Hz and 5000 Hz), m = 6 (for 1803 Hz) and m = 5 (for 1000Hz), and estimating D2 for each embedding dimension as the steepest slope of the linear scaling region in the plot of log(C(r)) vs log(r) provided the scaling region occurred within the middle third of the plot. Piloting revealed that fractal dimensions were higher from the sampled abnormal than normal achromatic VEPs in adults (p = 0.02). Variances of fractal dimension were higher from the abnormal than normal chromatic VEPs in children (p = 0.01). Conclusions A useful analysis protocol to assess the fractal dimension of transformed VEPs has been developed

Statistical practices of educational researchers: An analysis of their ANOVA, MANOVA, and ANCOVA analyses

Articles published in several prominent educational journals were examined to investigate the use of data-analysis tools by researchers in four research paradigms: between-subjects univariate designs, between-subjects multivariate designs, repeated measures designs, and covariance designs. In addition to examining specific details pertaining to the research design (e.g., sample size, group size equality/inequality) and methods employed for data analysis, we also catalogued whether: (a) validity assumptions were examined, (b) effect size indices were reported, (c) sample sizes were selected based on power considerations, and (d) appropriate textbooks and/or articles were cited to communicate the nature of the analyses that were performed. Our analyses imply that researchers rarely verify that validity assumptions are satisfied and accordingly typically use analyses that are nonrobust to assumption violations. In addition, researchers rarely report effect size statistics, nor do they routinely perform power analyses to determine sample size requirements. We offer many recommendations to rectify these shortcomings.Social Sciences and Humanities Research Counci

Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity.

Global dispersal and increasing frequency of the SARS-CoV-2 spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large dataset, well represented by both spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant

Changes in symptomatology, reinfection, and transmissibility associated with the SARS-CoV-2 variant B.1.1.7: an ecological study

Background The SARS-CoV-2 variant B.1.1.7 was first identified in December, 2020, in England. We aimed to investigate whether increases in the proportion of infections with this variant are associated with differences in symptoms or disease course, reinfection rates, or transmissibility. Methods We did an ecological study to examine the association between the regional proportion of infections with the SARS-CoV-2 B.1.1.7 variant and reported symptoms, disease course, rates of reinfection, and transmissibility. Data on types and duration of symptoms were obtained from longitudinal reports from users of the COVID Symptom Study app who reported a positive test for COVID-19 between Sept 28 and Dec 27, 2020 (during which the prevalence of B.1.1.7 increased most notably in parts of the UK). From this dataset, we also estimated the frequency of possible reinfection, defined as the presence of two reported positive tests separated by more than 90 days with a period of reporting no symptoms for more than 7 days before the second positive test. The proportion of SARS-CoV-2 infections with the B.1.1.7 variant across the UK was estimated with use of genomic data from the COVID-19 Genomics UK Consortium and data from Public Health England on spike-gene target failure (a non-specific indicator of the B.1.1.7 variant) in community cases in England. We used linear regression to examine the association between reported symptoms and proportion of B.1.1.7. We assessed the Spearman correlation between the proportion of B.1.1.7 cases and number of reinfections over time, and between the number of positive tests and reinfections. We estimated incidence for B.1.1.7 and previous variants, and compared the effective reproduction number, Rt, for the two incidence estimates. Findings From Sept 28 to Dec 27, 2020, positive COVID-19 tests were reported by 36 920 COVID Symptom Study app users whose region was known and who reported as healthy on app sign-up. We found no changes in reported symptoms or disease duration associated with B.1.1.7. For the same period, possible reinfections were identified in 249 (0·7% [95% CI 0·6–0·8]) of 36 509 app users who reported a positive swab test before Oct 1, 2020, but there was no evidence that the frequency of reinfections was higher for the B.1.1.7 variant than for pre-existing variants. Reinfection occurrences were more positively correlated with the overall regional rise in cases (Spearman correlation 0·56–0·69 for South East, London, and East of England) than with the regional increase in the proportion of infections with the B.1.1.7 variant (Spearman correlation 0·38–0·56 in the same regions), suggesting B.1.1.7 does not substantially alter the risk of reinfection. We found a multiplicative increase in the Rt of B.1.1.7 by a factor of 1·35 (95% CI 1·02–1·69) relative to pre-existing variants. However, Rt fell below 1 during regional and national lockdowns, even in regions with high proportions of infections with the B.1.1.7 variant. Interpretation The lack of change in symptoms identified in this study indicates that existing testing and surveillance infrastructure do not need to change specifically for the B.1.1.7 variant. In addition, given that there was no apparent increase in the reinfection rate, vaccines are likely to remain effective against the B.1.1.7 variant. Funding Zoe Global, Department of Health (UK), Wellcome Trust, Engineering and Physical Sciences Research Council (UK), National Institute for Health Research (UK), Medical Research Council (UK), Alzheimer's Society