Mechanisms of DNA-Mediated Allostery

Temperature-induced swelling and unwinding of double-stranded DNA

We utilized all-atom molecular dynamics simulations to investigate the temperature-induced swelling and unwinding of double-stranded DNA (dsDNA). We adopted three helical parameters, specifically helical twist, helical rise, and diameter, to quantitatively describe the deformations and elastic properties associated with swelling and unwinding processes within an orthogonal cylindrical coordinate system. The results indicate that as temperature increases, dsDNA experiences a weak swelling accompanied by unwinding. This is associated with a slight increase in helical rise, while the helical diameter almost remains unchanged and the helical twist decreases. We evaluated all potential pathways for unwinding and elucidated that twist-diameter coupling drives the unwinding from an entropy perspective. On the other hand, we employed the rigid base pair model to examine the swelling and unwinding elasticities, with a focus on the stiffnesses of twist and diameter. The results suggest that the temperature induces variations in the local twist and diameter elasticities, as well as their couplings of dsDNA, which are closely related to the distance between the base pairs, attributed to its thermal fluctuations and correlations. The global twist elasticity reduces as the temperature rises; nonetheless, the global diameter elasticity and the twist-diameter coupling can be considered as constants, which indicate independence from the increasing temperature.

Stacking correlation length in single-stranded DNA

Base stacking is crucial in nucleic acid stabilization, from DNA duplex hybridization to single-stranded DNA (ssDNA) protein binding. While stacking energies are tiny in ssDNA, they are inextricably mixed with hydrogen bonding in DNA base pairing, making their measurement challenging. We conduct unzipping experiments with optical tweezers of short poly-purine (dA and alternating dG and dA) sequences of 20-40 bases. We introduce a helix-coil model of the stacking-unstacking transition that includes finite length effects and reproduces the force-extension curves. Fitting the model to the experimental data, we derive the stacking energy per base, finding the salt-independent value $\Delta G_0^{ST}=0.14(3)$ kcal/mol for poly-dA and $\Delta G_0^{ST}=0.07(3)$ kcal/mol for poly-dGdA. Stacking in these polymeric sequences is predominantly cooperative with a correlation length of ∼4 bases at zero force . The correlation length reaches a maximum of ∼10 and 5 bases at the stacking-unstacking transition force of ∼10 and 20 pN for poly-dA and poly-dGdA, respectively. The salt dependencies of the cooperativity parameter in ssDNA and the energy of DNA hybridization are in agreement, suggesting that double-helix stability is primarily due to stacking. Analysis of poly-rA and poly-rC RNA sequences shows a larger stacking stability but a lower stacking correlation length of ∼2 bases.

Mechanical response of double-stranded DNA: Bend, twist, and overwind

We employed all-atom molecular dynamics simulations to explore the mechanical response of bending, twisting, and overwinding for double-stranded DNA (dsDNA). We analyzed the bending and twisting deformations, as well as their stiffnesses, using the tilt, roll, and twist modes under stretching force. Findings indicate that the roll and twist angles vary linearly with the stretching force but show opposite trends. The tilt, roll, and twist elastic moduli are considered constants, while the coupling between roll and twist modes slightly decreases under stretching force. The effect of the stretching force on the roll and twist modes, including both their deformations and elasticities, exhibits sequence-dependence, with symmetry around the base pair step. Furthermore, we examined the overwinding path and mechanism of dsDNA from the perspective of the stiffness matrix, based on the tilt, roll, and twist modes. The correlations among tilt, roll, and twist angles imply an alternative overwinding pathway via twist-roll coupling when dsDNA is stretched, wherein entropic contribution prevails.

Role of Shape Deformation of DNA-Binding Sites in Regulating the Efficiency and Specificity in Their Recognition by DNA-Binding Proteins

Accurate transcription of genetic information is crucial, involving precise recognition of the binding motifs by DNA-binding proteins. While some proteins rely on short-range hydrophobic and hydrogen bonding interactions at binding sites, others employ a DNA shape readout mechanism for specific recognition. In this mechanism, variations in DNA shape at the binding motif resulted from either inherent flexibility or binding of proteins at adjacent sites are sensed and capitalized by the searching proteins to locate them specifically. Through extensive computer simulations, we investigate both scenarios to uncover the underlying mechanism and origin of specificity in the DNA shape readout mechanism. Our findings reveal that deformation in shape at the binding motif creates an entropy funnel, allowing information about altered shapes to manifest as fluctuations in minor groove widths. This signal enhances the efficiency of nonspecific search of nearby proteins by directing their movement toward the binding site, primarily driven by a gain in entropy. We propose this as a generic mechanism for DNA shape readout, where specificity arises from the alignment between the molecular frustration of the searching protein and the ruggedness of the entropic funnel governed by molecular features of the protein and arrangement of the DNA bases at the binding site, respectively.

Insights into elastic properties of coarse-grained DNA models: q-stiffness of cgDNA vs cgDNA

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs all-atom models using stiffness matrices in Fourier space (q-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version, cgDNA+. Our analysis shows that while cgDNA+ closely follows the q-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables, which are rather strong in the q → 0 (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects that manifest themselves at longer length scales. In particular, cgDNA performs poorly on twist stiffness, with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist are strongly length scale dependent: DNA is torsionally soft at a few base pair distances but becomes more rigid at distances of a few dozen base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated.