Random fluctuations and noise are ubiquitous within nature and will generally occur on some length scale. Consider long linear DNA strands suspended within a solution at some temperature and salt concentration Or in-vivo, compactified within a living cell, within a biological environment with average ambient temperature T. Thermal fluctuations will cause the DNA ribbon to fluctuate weakly about its average geometry. In addition, the thermal fluctuations will also result in the spontaneous opening and closing of the bonds or bridges spanning the edges of the DNA ribbon. The DNA ribbon is therefore not a static structure but fluctuates randomly, essentially becoming a stochastic or random geometry.} The aim therefore, will be to describe the fluctuating DNA ribbon mathematically, in terms of random geometrical and topological terms. First, some statistical mechanical issues are considered. The following proposition follows the arguments in Huang. The following assumptions are made:

A DNA suspension within a liquid at temperature and some given salt concentration, or in=vivo within a cell. It can be considered within the canonical ensemble as a subsystem within a heat bath at fixed temperature. The DNA strand of length and sugar-phosphate backbone edges has N bonds or weak hydrogen bridges across the ribbon binding the A-T and G-C nucleic acids. At temperature T the ribbon weakly fluctuates and the bridges open and close. A open bond can be assigned an energy and an open bond has energy . A link can be open only if all the other links along the ribbon are already open. The possible states are labeled by the number of open bridges or links so that the energy with open links is The at low temperatures , and so there are few open links in that . At high temperatures , and so most of the links are open in that

The partition function is

where is the inverse temperature. The expected number of open links at temperature is then given by the thermodynamic relation

so that

Then at low temperatures , and so there are few open links in that . At high temperatures , and so most of the links are open in that

This suggests that DNA will exist and function only within certain temperature ranges. The molecule dissociates at high temperatures. (Which also implies that DNA-based life cannot arise and evolve in high-temperature conditions in the universe.) Superhelicity also plays a role in that the free energy associated with supercoiling is of the quadratic form , where K is a proporitonality constant. This quadratic relation has been experimentally shown (refs.) DNA molecules that differ only in are topoisomers. The concentration of a given topoisomer with linking number at temperature T, is the Boltzmann distribution

which is a normal or Gaussian distribution centred on , the relaxed linking number, and with standard deviation . Proteins which regulate the supercoiling and maintain this distribution are known as topoisomerases. For high temperatures .

Taking into account thermodynamic, mechanical and biophysical properties of fluctuating DNA molecules within a thermal bath is complex. However, the aim of this section is to focus on a tractable stochastic mathematical description of a fluctuating DNA ribbon, in terms of a ‘stochastic geometry and topology’ which arises from contact with a thermal bath, or more precisely, a thermal noise bath.

\item This introduces analogies with Brownian motion and turbulence, which are studied and described on length scales as stochastic processes and where the deeper underlying statistical mechanics can be ignored. Indeed, the stochastic fluctuations of the DNA double-helix structure when immersed in a thermal bath whether in-vitro or in vivo, is a form of Brownian motion, arising from bombardment of the double helix by a ‘warm sea’ of smaller surrounding molecules.

For proteins, which consist of mobile helices, strands and sheets, Brownian bombardment from water molecules will also play an important role in the folding of the protein. Simple Brownian motion of a particle in a thermal bath at temperature T is described by a linear SDE of the for

where is the position of the Brownian particle. The random field is a thermal white noise which described the effect of the thermal Brownian bombardment of the particle by the smaller water molecules. Generally

and by the

equipartition theory , where is the statistical expectation, k is Boltzmann constant. The solution is the well-defined stochastic Ito integral

Suppose by analogy, we consider the tangent vectors and along the edges of the DNA ribbon. If the DNA ribbon undergoes a Brownian bombardment from smaller molecules, these vectors will be randomly perturbed as the ribbon fluctuates randomly about its average geometry/topology.

where is a Gaussian random vector field defined at all points in a domain containing the ribbon. GRVFs are discussed in the Appendix. The linking number will not be affected by thermal fluctuations of the DNA ribbon but the linking number $\mathcal{LK}=\mathcal{TW}+\mathcal{WR}$ will be randomly partitioned into twist and writhe which will randomly change. The total remains constant however. This can be made more precise as follows:

Let be a DNA ribbon with linking number and superhelicity and let

and be tangent vectors at and . The Wilson lines along the helices or ribbon edges are

If the domain is a thermal bath at some temperature then the tangent vectors can be interpreted as random Gaussian vector fields and such that

where the GRVF obeys the statistics

and where is a 2-point function or covariance regulated at and is a correlation length such that for .

(Appendix A.)

The stochastic Wilson lines along the helices are then

The stochastic expectations of the Wilson lines are

\

The stochastic integral

can be evaluated by a cluster integral decomposition, exploiting the Gaussian property of the random field . This is briefly stated and proved in the following theorem. Let be a random Gaussian vector field defined for all with respect to a probability space and having a 2-point function . Then the expectations of the stochastic Wilson loops along h and are then

latex \mathbf{\mathcal{W}}[\bar{h}]=\mathbf{E}\bigg(exp\bigg(i\int_{\bar{h}}\widehat{\mathcal{F}}_{i}(x)dx^{i}\bigg)\bigg)= exp\bigg(-\int_{h}\int_{\bar{h}}C_{ij}(x,y;\zeta)dx^{i}dy^{j}\bigg)$

where

. If for then as

## Leave a comment

Comments feed for this article