Genomic DNA is constantly subjected to various mechanical stresses arising from its biological functions and cell packaging. If the local mechanical properties of DNA change under torsional and tensional stress, the activity of DNA-modifying proteins and transcription factors can be affected and regulated allosterically. To check this possibility, appropriate steady forces and torques were applied in the course of all-atom molecular dynamics simulations of DNA with AT-and GC-alternating sequences. It is found that the stretching rigidity grows with tension as well as twisting. The torsional rigidity is not affected by stretching, but it varies with twisting very strongly, and differently for the two sequences. Surprisingly, for AT-alternating DNA it passes through a minimum with the average twist close to the experimental value in solution. For this fragment, but not for the GC-alternating sequence, the bending rigidity noticeably changes with both twisting and stretching. The results have important biological implications and shed light upon earlier experimental observations. PACS numbers: 87.14.gk 87.15.H-87.15.ap 87.15.ak
IntroductionInternal mechanical stress is ubiquitous in the biologically active state of double helical DNA. In eucaryotic cells, DNA is densely packed in chromosomes and forced to bend, twist and stretch by numerous protein factors involved in genome regulation [1,2]. In procaryotes, DNA is subjected to a constitutive unwinding torque maintained by special enzymes, which leads to supercoiling, as in a long rope with bending and twisting elasticity [3,4]. The supercoiling and, more generally, stressinduced DNA forms are key factors in a variety of cellular processes [5]. For instance, the degree of supercoiling in bacteria changes systematically during the cell cycle and in response to environmental conditions, which is accompanied by activation or suppression of certain genes [6]. The promoter sensitivity to supercoiling stems from the recognition of short promoter elements by RNA polymerase [7]. Detailed studies indicate that it probably does not require DNA melting nor transitions to alternative forms [6]. In E. coli, relaxation of the superhelical stress simultaneously alters the activity of 306 genes (7% of the genome), with 106 genes activated and others deactivated [8]. The genes concerned are functionally diverse and widely dispersed throughout the chromosome, and the effect is dose-dependent.The physical mechanisms of such effects are understood only partially. Long DNA is well described by the coarse-grained worm-like chain (WLC) model [9,10] supplemented with harmonic twisting and stretching elasticity [11][12][13][14][15][16]. This model nicely explains the stressmodulated probability of looping, wrapping around proteins, and juxtaposition of distant protein binding sites [17]. However, it cannot account for the promoter sen- * Electronic address: alexey@ibpc.fr sitivity to supercoiling, for instance, because in this and many other cases the gene regulation has a strong local...