Sequence-dependent response of DNA to torsional stress: a potential biological regulation mechanism

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.

Reactions of deoxyribonucleotide bases with sulfooxymethyl or halomethyl polycyclic aromatic hydrocarbons induce unwinding of DNA supercoils

Torsional stress in double-stranded DNA enables and regulates facets of chromosomal metabolism, replication, and transcription and requires regulatory enzymatic systems including topoisomerases and histone methyltransferases. As such, this machinery may be subject to deleterious effects from reactive mutagens, including ones from carcinogenic polycyclic aromatic hydrocarbon (PAH) adduct formation with DNA. Supercoiled plasmid DNA was investigated for its torsional responses to adducts formed in vitro from PAH benzylic carbocation reactive intermediates created spontaneously by release of leaving groups. PAH sulfate esters were found to (1) unwind DNA in a concentration dependent manner, and (2) provide maximum unwinding in a pattern consistent with known carcinogenicities of the parent PAHs, that is, 6-methylbenzo[a]pyrene > 7,12-methylbenz[a]anthracene > 3-methylcholanthrene > 9-methylanthracene > 7-methylbenz[a]anthracene > 1-methylpyrene. Supercoil unwinding was demonstrated to be dependent on the presence of sulfate or chloride leaving groups such that reactive carbocations were generated in situ by hydrolysis. In silico modeling of intercalative complex topology showed PAH benzylic carbocation reactive functional groups in alignment with target nucleophiles on guanine bases in a 5'-dCdG-3' pocket in agreement with known formation of nucleotide adducts. Inhibitory or modulatory effects on PAH-induced supercoil unwinding were seen with ascorbic acid and an experimental antineoplastic agent Antineoplaston A10 in agreement with their known anticarcinogenic properties. In summary, the reactive PAH intermediates studied here undoubtedly participate in well-known mutational mechanisms such as frameshifts and apurinic site generation. However, they are also capable of random disruption of chromosomal supercoiling in a manner consistent with the known carcinogenicities of the parent compounds, and this mechanism may represent an additional detrimental motif worthy of further study for a more complete understanding of chemical carcinogenicity.

Structural and dynamical aspect of DNA motif sequence specific binding of AP-1 transcription factor

Activator protein-1 (AP-1) comprises one of the largest and most evolutionary conserved families of ubiquitous eukaryotic transcription factors that act as a pioneer factor. Diversity in DNA binding interaction of AP-1 through a conserved basic-zipper (bZIP) domain directs in-depth understanding of how AP-1 achieves its DNA binding selectivity and consequently gene regulation specificity. Here, we address the structural and dynamical aspects of the DNA target recognition process of AP-1 using microsecond-long atomistic simulations based on the structure of the human AP-1 FosB/JunD bZIP-DNA complex. Our results show the unique role of DNA shape features in selective base specific interactions, characteristic ion population, and solvation properties of DNA grooves to form the motif sequence specific AP-1-DNA complex. The TpG step at the two terminals of the AP-1 site plays an important role in the structural adjustment of DNA by modifying the helical twist in the AP-1 bound state. We addressed the role of intrinsic motion of the bZIP domain in terms of opening and closing gripper motions of DNA binding helices, in target site recognition and binding of AP-1 factors. Our observations suggest that binding to the cognate motif in DNA is mainly accompanied with the precise adjustment of closing gripper motion of DNA binding helices of the bZIP domain.

Spectral analysis of DNA superhelical dynamics from molecular minicircle simulations

Torsional and bending deformations of DNA molecules often occur in vivo and are important for biological functions. DNA "under stress" is a conformational state, which is by far the most frequent state during DNA-protein and gene regulation. In DNA minicircles of length <100 base pairs (bp), the combined effect of torsional and bending stresses can cause local unusual conformations, with certain base pair steps often absorbing most of the stress, leaving other steps close to their relaxed conformation. To better understand the superhelical dynamics of DNA under stress, molecular simulations of 94 bp minicircles with different torsional linking numbers were interpreted using Fourier analyses and principal component analyses. Sharp localized bends of nearly 90° in the helical axis were observed, which in turn decreased fluctuations of the rotational register and helped redistribute the torsional stress into writhe, i.e., superhelical turn up to 360°. In these kinked minicircles, only two-thirds of the DNA molecule bends and writhes and the remaining segment stays close to straight and preserves a conformational flexibility typical of canonical B-DNA (bending of 39° ± 17° distributed parsimoniously across 36 bp), which was confirmed and visualized by principal component analysis. These results confirm that stressed DNA molecules are highly heterogeneous along their sequence, with segments designed to locally store and release stress so that nearby segments can stay relaxed.

The development of nucleic acids force fields: From an unchallenged past to a competitive future

Molecular dynamics simulations have strongly matured as a method to study biomolecular processes. Their validity, however, is determined by the accuracy of the underlying force fields that describe the forces between all atoms. In this article, we review the development of nucleic acids force fields. We describe the early attempts in the 1990s and emphasize their strong influence on recent force fields. State-of-the-art force fields still use the same Lennard-Jones parameters derived 25 years ago in spite of the fact that these parameters were in general not fitted for nucleic acids. In addition, electrostatic parameters also are deprecated, which may explain some of the current force field deficiencies. We compare different force fields for various systems and discuss new tests of the recently developed Tumuc1 force field. The OL-force fields and Tumuc1 are arguably the best force fields to describe the DNA double helix. However, no force field is flawless. In particular, the description of sugar-puckering remains a problem for nucleic acids force fields. Future refinements are required, so we review methods for force field refinement and give an outlook to the future of force fields.

Prediction of hub genes and key pathways associated with the radiation response of human hematopoietic stem/progenitor cells using integrated bioinformatics methods

Hematopoietic stem cells (HSCs) are indispensable for the maintenance of the entire blood program through cytokine response. However, HSCs have high radiosensitivity, which is often a problem during radiation therapy and nuclear accidents. Although our previous study has reported that the combination cytokine treatment (interleukin-3, stem cell factor, and thrombopoietin) improves the survival of human hematopoietic stem/progenitor cells (HSPCs) after radiation, the mechanism by which cytokines contribute to the survival of HSPCs is largely unclear. To address this issue, the present study characterized the effect of cytokines on the radiation-induced gene expression profile of human CD34⁺ HSPCs and explored the hub genes that play key pathways associated with the radiation response using a cDNA microarray, a protein-protein interaction-MCODE module analysis and Cytohubba plugin tool in Cytoscape. This study identified 2,733 differentially expressed genes (DEGs) and five hub genes (TOP2A, EZH2, HSPA8, GART, HDAC1) in response to radiation in only the presence of cytokines. Furthermore, functional enrichment analysis found that hub genes and top DEGs based on fold change were enriched in the chromosome organization and organelle organization. The present findings may help predict the radiation response and improve our understanding of this response of human HSPCs.

Manipulation of topoisomerase expression inhibits cell division but not growth and reveals a distinctive promoter structure in Synechocystis

In cyanobacteria DNA supercoiling varies over the diurnal cycle and is integrated with temporal programs of transcription and replication. We manipulated DNA supercoiling in Synechocystis sp. PCC 6803 by CRISPRi-based knockdown of gyrase subunits and overexpression of topoisomerase I (TopoI). Cell division was blocked but cell growth continued in all strains. The small endogenous plasmids were only transiently relaxed, then became strongly supercoiled in the TopoI overexpression strain. Transcript abundances showed a pronounced 5'/3' gradient along transcription units, incl. the rRNA genes, in the gyrase knockdown strains. These observations are consistent with the basic tenets of the homeostasis and twin-domain models of supercoiling in bacteria. TopoI induction initially led to downregulation of G+C-rich and upregulation of A+T-rich genes. The transcriptional response quickly bifurcated into six groups which overlap with diurnally co-expressed gene groups. Each group shows distinct deviations from a common core promoter structure, where helically phased A-tracts are in phase with the transcription start site. Together, our data show that major co-expression groups (regulons) in Synechocystis all respond differentially to DNA supercoiling, and suggest to re-evaluate the long-standing question of the role of A-tracts in bacterial promoters.

An updated definition of V(D)J recombination signal sequences revealed by high-throughput recombination assays

In the adaptive immune system, V(D)J recombination initiates the production of a diverse antigen receptor repertoire in developing B and T cells. Recombination activating proteins, RAG1 and RAG2 (RAG1/2), catalyze V(D)J recombination by cleaving adjacent to recombination signal sequences (RSSs) that flank antigen receptor gene segments. Previous studies defined the consensus RSS as containing conserved heptamer and nonamer sequences separated by a less conserved 12 or 23 base-pair spacer sequence. However, many RSSs deviate from the consensus sequence. Here, we developed a cell-based, massively parallel assay to evaluate V(D)J recombination activity on thousands of RSSs where the 12-RSS heptamer and adjoining spacer region contained randomized sequences. While the consensus heptamer sequence (CACAGTG) was marginally preferred, V(D)J recombination was highly active on a wide range of non-consensus sequences. Select purine/pyrimidine motifs that may accommodate heptamer unwinding in the RAG1/2 active site were generally preferred. In addition, while different coding flanks and nonamer sequences affected recombination efficiency, the relative dependency on the purine/pyrimidine motifs in the RSS heptamer remained unchanged. Our results suggest RAG1/2 specificity for RSS heptamers is primarily dictated by DNA structural features dependent on purine/pyrimidine pattern, and to a lesser extent, RAG:RSS base-specific interactions.

TripDesign: A DNA Triplex Design Approach Based on Interaction Forces

A DNA triplex has the advantages of improved nanostructure stability and pH environment responsiveness compared with single-stranded and double-stranded nucleic acids. However, sequence stability and low design efficiency hinder the application of DNA triplexes. Therefore, a DNA triplex design approach (TripDesign) based on interaction forces is proposed. First, we present the stacking force constraint, torsional stress constraint, and G-quadruplex motif constraint and then use an improved memetic algorithm to design triplex sequences under combinatorial constraints. Finally, to quantify the process of triplex formation, we also explore the minimum length of the triplex-forming oligos (TFOs) required to form the triplex and the factors that produce depletion in cyclic pH-jump experiments. The experimental results show that the sequences produced by TripDesign have high stability and reversibility, and the proposed approach achieves efficient and automatic sequence design. In addition, this study characterizes multiple basic parameters of DNA triplex formation and promotes the wider application of DNA triplexes in nanotechnology.

RNA kink-turns are highly anisotropic with respect to lateral displacement of the flanking stems

Kink-turns are highly bent internal loop motifs commonly found in the ribosome and other RNA complexes. They frequently act as binding sites for proteins and mediate tertiary interactions in larger RNA structures. Kink-turns have been a topic of intense research, but their elastic properties in the folded state are still poorly understood. Here we use extensive all-atom molecular dynamics simulations to parameterize a model of kink-turn in which the two flanking helical stems are represented by effective rigid bodies. Time series of the full set of six interhelical coordinates enable us to extract minimum energy shapes and harmonic stiffness constants for kink-turns from different RNA functional classes. The analysis suggests that kink-turns exhibit isotropic bending stiffness but are highly anisotropic with respect to lateral displacement of the stems. The most flexible lateral displacement mode is perpendicular to the plane of the static bend. These results may help understand the structural adaptation and mechanical signal transmission by kink-turns in complex natural and artificial RNA structures.

Homologous basic helix–loop–helix transcription factors induce distinct deformations of torsionally-stressed DNA: a potential transcription regulation mechanism

Changing torsional restraints on DNA is essential for the regulation of transcription. Torsional stress, introduced by RNA polymerase, can propagate along chromatin facilitating topological transitions and modulating the specific binding of transcription factors (TFs) to DNA. Despite the importance, the mechanistic details on how torsional stress impacts the TFs-DNA complexation remain scarce. Herein, we address the impact of torsional stress on DNA complexation with homologous human basic helix–loop–helix (BHLH) hetero- and homodimers: MycMax, MadMax and MaxMax. The three TF dimers exhibit specificity towards the same DNA consensus sequence, the E-box response element, while regulating different transcriptional pathways. Using microseconds-long atomistic molecular dynamics simulations together with the torsional restraint that controls DNA total helical twist, we gradually over- and underwind naked and complexed DNA to a maximum of ± 5°/bp step. We observe that the binding of the BHLH dimers results in a similar increase in DNA torsional rigidity. However, under torsional stress the BHLH dimers induce distinct DNA deformations, characterised by changes in DNA grooves geometry and a significant asymmetric DNA bending. Supported by bioinformatics analyses, our data suggest that torsional stress may contribute to the execution of differential transcriptional programs of the homologous TFs by modulating their collaborative interactions.

Sequence-dependent twist-bend coupling in DNA minicircles

Looping of double-stranded DNA molecules with 100-200 base pairs into minicircles, catenanes, and rotaxanes has been suggested as a potential tool for DNA nanotechnologies. However, sharp DNA bending into a minicircle with a diameter of several to ten nanometers occurs with alterations in the DNA helical structure and may lead to defective kink formation that hampers the use of DNA minicircles, catenanes, and rotaxanes in nanoscale DNA applications. Here, we investigated local variations of a helical twist in sharply bent DNA using microsecond-long all-atom molecular dynamics simulations of six different DNA minicircles, focusing on the sequence dependence of the coupling between DNA bending and its helical twist. Twist angles between consecutive base pairs were analyzed at different locations relative to the direction of DNA bending and, among 10 unique dinucleotide steps, we identified four dinucleotide steps with strong twist-bend coupling, the pyrimidine-purine dinucleotide steps of TA/TA, CG/CG, and CA/TG and the purine-purine dinucleotide step of GA/TC. This work suggests the sequence-dependent structural responses of DNA to strong mechanical deformation, providing new molecular-level insights into the structure and stability of sharply bent DNA minicircles for nanoscale applications.

Membrane Activity of a DNA-Based Ion Channel Depends on the Stability of Its Double-Stranded Structure

DNA nanotechnology has emerged as a promising method for designing spontaneously inserting and fully controllable synthetic ion channels. However, both insertion efficiency and stability of existing DNA-based membrane channels leave much room for improvement. Here, we demonstrate an approach to overcoming the unfavorable DNA-lipid interactions that hinder the formation of a stable transmembrane pore. Our all-atom MD simulations and experiments show that the insertion-driving cholesterol modifications can cause fraying of terminal base pairs of nicked DNA constructs, distorting them when embedded in a lipid bilayer. Importantly, we show that DNA nanostructures with no backbone discontinuities form more stable conductive pores and insert into membranes with a higher efficiency than the equivalent nicked constructs. Moreover, lack of nicks allows design and maintenance of membrane-spanning helices in a tilted orientation within the lipid bilayer. Thus, reducing the conformational degrees of freedom of the DNA nanostructures enables better control over their function as synthetic ion channels.

Loop-seq: A high-throughput technique to measure the mesoscale mechanical properties of DNA

The local mechanical properties of the DNA polymer influence molecular processes in biology that require mechanical deformations of DNA. Lack of suitable high-throughput experimental techniques had precluded measuring how these properties might vary with sequence along the vast lengths of genomes. Here, we present a detailed protocol for a recently developed experimental technique called loop-seq, which measures at least one local mechanical property of DNA-its propensity to cyclize-in genome-scale throughput. Loop-seq has been used to obtain experimentally derived genome-wide maps of a physical property of DNA. Such measurements have revealed that diverse DNA-deforming processes involved in chromatin organization at various genomic loci are regulated by the genetically encoded, sequence-dependent variations in the mechanical properties of DNA.

Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors

Activator proteins 1 (AP-1) comprise one of the largest families of eukaryotic basic leucine zipper transcription factors. Despite advances in the characterization of AP-1 DNA-binding sites, our ability to predict new binding sites and explain how the proteins achieve different gene expression levels remains limited. Here we address the role of sequence-specific DNA flexibility for stability and specific binding of AP-1 factors, using microsecond-long molecular dynamics simulations. As a model system, we employ yeast AP-1 factor Yap1 binding to three different response elements from two genetic environments. Our data show that Yap1 actively exploits the sequence-specific flexibility of DNA within the response element to form stable protein–DNA complexes. The stability also depends on the four to six flanking nucleotides, adjacent to the response elements. The flanking sequences modulate the conformational adaptability of the response element, making it more shape-efficient to form specific contacts with the protein. Bioinformatics analysis of differential expression of the studied genes supports our conclusions: the stability of Yap1–DNA complexes, modulated by the flanking environment, influences the gene expression levels. Our results provide new insights into mechanisms of protein–DNA recognition and the biological regulation of gene expression levels in eukaryotes.

Torsional stress can regulate the unwrapping of two outer half superhelical turns of nucleosomal DNA

Torsional stress has a significant impact on the structure and stability of the nucleosome. RNA polymerase imposes torsional stress on the DNA in chromatin and unwraps the DNA from the nucleosome to access the genetic information encoded in the DNA. To understand how the torsional stress affects the stability of the nucleosome, we examined the unwrapping of two half superhelical turns of nucleosomal DNA from either end of the DNA under torsional stress with all-atom molecular dynamics simulations. The free energies for unwrapping the DNA indicate that positive stress that overtwists DNA facilitates a large-scale asymmetric unwrapping of the DNA without a large extension of the DNA. During the unwrapping, one end of the DNA was dissociated from H3 and H2A-H2B, while the other end of the DNA stably remained wrapped. The detailed analysis indicates that this asymmetric dissociation is facilitated by the geometry and bendability of the DNA under positive stress. The geometry stabilized the interaction between the major groove of the twisted DNA and the H3 αN-helix, and the straightened DNA destabilized the interaction with H2A-H2B. Under negative stress, the DNA became more bendable and flexible, which facilitated the binding of the unwrapped DNA to the octamer in a stable state. Consequently, we conclude that the torsional stress has a significant impact on the affinity of the DNA and the octamer through the inherent nature of the DNA and can change the accessibility of regulatory proteins.

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

Energetics of twisted DNA topologies

Our goal is to review the main theoretical models used to calculate free energy changes associated with common, torsion-induced conformational changes in DNA and provide the resulting equations hoping to facilitate quantitative analysis of both in vitro and in vivo studies. This review begins with a summary of work regarding the energy change of the negative supercoiling-induced B- to L-DNA transition, followed by a discussion of the energetics associated with the transition to Z-form DNA. Finally, it describes the energy changes associated with the formation of DNA curls and plectonemes, which can regulate DNA-protein interactions and promote cross talk between distant DNA elements, respectively. The salient formulas and parameters for each scenario are summarized in table format to facilitate comparison and provide a concise, user-friendly resource.

DNA mechanics and its biological impact

Almost all nucleoprotein interactions and DNA manipulation events involve mechanical deformations of DNA. Extraordinary progresses in single-molecule, structural, and computational methods have characterized the average mechanical properties of DNA, such as bendability and torsional rigidity, in high resolution. Further, the advent of sequencing technology has permitted measuring, in high-throughput, how such mechanical properties vary with sequence and epigenetic modifications along genomes. We review these recent technological advancements, and discuss how they have contributed to the emerging idea that variations in the mechanical properties of DNA play a fundamental role in regulating, genome-wide, diverse processes involved in chromatin organization.

Specifically bound BZIP transcription factors modulate DNA supercoiling transitions

Torsional stress on DNA, introduced by molecular motors, constitutes an important regulatory mechanism of transcriptional control. Torsional stress can modulate specific binding of transcription factors to DNA and introduce local conformational changes that facilitate the opening of promoters and nucleosome remodelling. Using all-atom microsecond scale molecular dynamics simulations together with a torsional restraint that controls the total twist of a DNA fragment, we address the impact of torsional stress on DNA complexation with a human BZIP transcription factor, MafB. We gradually over- and underwind DNA alone and in complex with MafB by 0.5° per dinucleotide step, starting from the relaxed state to a maximum of 5° per dinucleotide step, monitoring the evolution of the protein-DNA contacts at different degrees of torsional strain. Our computations show that MafB changes the DNA sequence-specific response to torsional stress. The dinucleotide steps that are susceptible to absorbing most of the torsional stress become more torsionally rigid, as they are involved in protein-DNA contacts. Also, the protein undergoes substantial conformational changes to follow the stress-induced DNA deformation, but mostly maintains the specific contacts with DNA. This results in a significant asymmetric increase of free energy of DNA twisting transitions, relative to free DNA, where overtwisting is more energetically unfavourable. Our data suggest that specifically bound BZIP factors could act as torsional stress insulators, modulating the propagation of torsional stress along the chromatin fibre, which might promote cooperative binding of collaborative DNA-binding factors.

How global DNA unwinding causes non-uniform stress distribution and melting of DNA

DNA unwinding is an important process that controls binding of proteins, gene expression and melting of double-stranded DNA. In a series of all-atom MD simulations on two DNA molecules containing a transcription start TATA-box sequence we demonstrate that application of a global restraint on the DNA twisting dramatically changes the coupling between helical parameters and the distribution of deformation energy along the sequence. Whereas only short range nearest-neighbor coupling is observed in the relaxed case, long-range coupling is induced in the globally restrained case. With increased overall unwinding the elastic deformation energy is strongly non-uniformly distributed resulting ultimately in a local melting transition of only the TATA box segment during the simulations. The deformation energy tends to be stored more in cytidine/guanine rich regions associated with a change in conformational substate distribution. Upon TATA box melting the deformation energy is largely absorbed by the melting bubble with the rest of the sequences relaxing back to near B-form. The simulations allow us to characterize the structural changes and the propagation of the elastic energy but also to calculate the associated free energy change upon DNA unwinding up to DNA melting. Finally, we design an Ising model for predicting the local melting transition based on empirical parameters. The direct comparison with the atomistic MD simulations indicates a remarkably good agreement for the predicted necessary torsional stress to induce a melting transition, for the position and length of the melted region and for the calculated associated free energy change between both approaches.

Structure-mechanics statistical learning unravels the linkage between local rigidity and global flexibility in nucleic acids

The mechanical properties of nucleic acids underlie biological processes ranging from genome packaging to gene expression, but tracing their molecular origin has been difficult due to the structural and chemical complexity. We posit that concepts from machine learning can help to tackle this long-standing challenge. Here, we demonstrate the feasibility and advantage of this strategy through developing a structure-mechanics statistical learning scheme to elucidate how local rigidity in double-stranded (ds)DNA and dsRNA may lead to their global flexibility in bend, stretch, and twist. Specifically, the mechanical parameters in a heavy-atom elastic network model are computed from the trajectory data of all-atom molecular dynamics simulation. The results show that the inter-atomic springs for backbone and ribose puckering in dsRNA are stronger than those in dsDNA, but are similar in strengths for base-stacking and base-pairing. Our analysis shows that the experimental observation of dsDNA being easier to bend but harder to stretch than dsRNA comes mostly from the respective B- and A-form topologies. The computationally resolved composition of local rigidity indicates that the flexibility of both nucleic acids is mostly due to base-stacking. But for properties like twist-stretch coupling, backbone springs are shown to play a major role instead. The quantitative connection between local rigidity and global flexibility sets foundation for understanding how local binding and chemical modification of genetic materials effectuate longer-ranged regulatory signals.

The mechanical properties of nucleic acids underlie biological processes ranging from genome packaging to gene expression. We devise structural mechanics statistical learning method to reveal their molecular origin in terms of chemical interactions.

The emergence of sequence-dependent structural motifs in stretched, torsionally constrained DNA

The double-helical structure of DNA results from canonical base pairing and stacking interactions. However, variations from steady-state conformations resulting from mechanical perturbations in cells have physiological relevance but their dependence on sequence remains unclear. Here, we use molecular dynamics simulations showing sequence differences result in markedly different structural motifs upon physiological twisting and stretching. We simulate overextension on different sequences of DNA ((AA)12, (AT)12, (CC)12 and (CG)12) with supercoiling densities at 200 and 50 mM salt concentrations. We find that DNA denatures in the majority of stretching simulations, surprisingly including those with over-twisted DNA. GC-rich sequences are observed to be more stable than AT-rich ones, with the specific response dependent on the base pair order. Furthermore, we find that (AT)12 forms stable periodic structures with non-canonical hydrogen bonds in some regions and non-canonical stacking in others, whereas (CG)12 forms a stacking motif of four base pairs independent of supercoiling density. Our results demonstrate that 20-30% DNA extension is sufficient for breaking B-DNA around and significantly above cellular supercoiling, and that the DNA sequence is crucial for understanding structural changes under mechanical stress. Our findings have important implications for the activities of protein machinery interacting with DNA in all cells.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

Characterizing the interplay between gene nucleotide composition bias and splicing

BACKGROUND

Nucleotide composition bias plays an important role in the 1D and 3D organization of the human genome. Here, we investigate the potential interplay between nucleotide composition bias and the regulation of exon recognition during splicing.

RESULTS

By analyzing dozens of RNA-seq datasets, we identify two groups of splicing factors that activate either about 3200 GC-rich exons or about 4000 AT-rich exons. We show that splicing factor-dependent GC-rich exons have predicted RNA secondary structures at 5' ss and are dependent on U1 snRNP-associated proteins. In contrast, splicing factor-dependent AT-rich exons have a large number of decoy branch points, SF1- or U2AF2-binding sites and are dependent on U2 snRNP-associated proteins. Nucleotide composition bias also influences local chromatin organization, with consequences for exon recognition during splicing. Interestingly, the GC content of exons correlates with that of their hosting genes, isochores, and topologically associated domains.

CONCLUSIONS

We propose that regional nucleotide composition bias over several dozens of kilobase pairs leaves a local footprint at the exon level and induces constraints during splicing that can be alleviated by local chromatin organization at the DNA level and recruitment of specific splicing factors at the RNA level. Therefore, nucleotide composition bias establishes a direct link between genome organization and local regulatory processes, like alternative splicing.

Hydrophobic catalysis and a potential biological role of DNA unstacking induced by environment effects

The main stabilizer of the DNA double helix is not the base-pair hydrogen bonds but coin-pile stacking of base pairs, whose hydrophobic cohesion, requiring abundant water, indirectly makes the DNA interior dry so that hydrogen bonds can exert full recognition power. We report that certain semihydrophobic agents depress the stacking energy (measurable in single-molecule experiments), leading to transiently occurring holes in the base-pair stack (monitorable via binding of threading intercalators). Similar structures observed in DNA complexes with RecA and Rad51, and previous observations of spontaneous strand exchange catalyzed in semihydrophobic model systems, make us propose that some hydrophobic protein residues may have roles in catalyzing homologous recombination. We speculate that hydrophobic catalysis is a general phenomenon in DNA enzymes.

Hydrophobic base stacking is a major contributor to DNA double-helix stability. We report the discovery of specific unstacking effects in certain semihydrophobic environments. Water-miscible ethylene glycol ethers are found to modify structure, dynamics, and reactivity of DNA by mechanisms possibly related to a biologically relevant hydrophobic catalysis. Spectroscopic data and optical tweezers experiments show that base-stacking energies are reduced while base-pair hydrogen bonds are strengthened. We propose that a modulated chemical potential of water can promote “longitudinal breathing” and the formation of unstacked holes while base unpairing is suppressed. Flow linear dichroism in 20% diglyme indicates a 20 to 30% decrease in persistence length of DNA, supported by an increased flexibility in single-molecule nanochannel experiments in poly(ethylene glycol). A limited (3 to 6%) hyperchromicity but unaffected circular dichroism is consistent with transient unstacking events while maintaining an overall average B-DNA conformation. Further information about unstacking dynamics is obtained from the binding kinetics of large thread-intercalating ruthenium complexes, indicating that the hydrophobic effect provides a 10 to 100 times increased DNA unstacking frequency and an “open hole” population on the order of 10⁻² compared to 10⁻⁴ in normal aqueous solution. Spontaneous DNA strand exchange catalyzed by poly(ethylene glycol) makes us propose that hydrophobic residues in the L2 loop of recombination enzymes RecA and Rad51 may assist gene recombination via modulation of water activity near the DNA helix by hydrophobic interactions, in the manner described here. We speculate that such hydrophobic interactions may have catalytic roles also in other biological contexts, such as in polymerases.

A sequence environment modulates the impact of methylation on the torsional rigidity of DNA

We explain, with atomic-level detail, how DNA methylation may contribute to the regulation of biological processes. We show that DNA sequence environment modulates the impact of methylation on the physical properties of the double helix.