A molecular view of DNA flexibility

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.

Shape and size tunability of sheets of interlocked ring copolymers

Mechanically bonded membranes of interlocked ring polymers are a significant generalization of conventional elastic sheets, where connectivity is provided by covalent bonding, and represent a promising class of topological meta-materials. In this context, two open questions regard the large-scale reverberations of the heterogeneous composition of the rings and the inequivalent modes of interlocking neighboring rings. We address these questions with Langevin dynamics simulations of chainmails with honeycomb-lattice connectivity, where the rings are block copolymers with two segments of different rigidity. We considered various combinations of the relative lengths of the two segments and the patterns of the over- and under-passes linking neighboring rings. We find that varying ring composition and linking patterns have independent and complementary effects. While the former sets the overall size of the chainmail, the latter defines the shape, enabling the selection of starkly different conformation types. Notably, one of the considered linking patterns favors saddle-shaped membranes, providing a first example of spontaneous negative Gaussian curvature in mechanically bonded sheets. The results help establish the extent to which mechanically bonded membranes can differ from conventional elastic ones, particularly for the achievable shape and size tunability.

Effect of ethanol on the elasticities of double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The elasticities of double-stranded (ds) DNA and RNA, which are critical to their biological functions and applications in materials science, can be significantly modulated by solution conditions such as ions and temperature. However, there is still a lack of a comprehensive understanding of the role of solvents in the elasticities of dsRNA and dsDNA in a comparative way. In this work, we explored the effect of ethanol solvent on the elasticities of dsRNA and dsDNA by magnetic tweezers and all-atom molecular dynamics simulations. We found that the bending persistence lengths and contour lengths of dsRNA and dsDNA decrease monotonically with the increase in ethanol concentration. Furthermore, the addition of ethanol weakens the positive twist-stretch coupling of dsRNA, while promotes the negative twist-stretch coupling of dsDNA. Counter-intuitively, the lower dielectric environment of ethanol causes a significant re-distribution of counterions and enhanced ion neutralization, which overwhelms the enhanced repulsion along dsRNA/dsDNA, ultimately leading to the softening in bending for dsRNA and dsDNA. Moreover, for dsRNA, ethanol causes slight ion-clamping across the major groove, which weakens the major groove-mediated twist-stretch coupling, while for dsDNA, ethanol promotes the stretch-radius correlation due to enhanced ion binding and consequently enhances the helical radius-mediated twist-stretch coupling.

Effect of temperature on anisotropic bending elasticity of dsRNA: an all-atom molecular dynamics simulation

Employing all-atom molecular dynamics simulations, we examined the temperature-dependent behavior of bending elasticity in double-stranded RNA (dsRNA). Specifically, we focused on the bending persistence length and its constituent components, namely, the tilt and roll stiffness. Our results revealed a near-linear decrease in these stiffness components as a function of temperature, thereby highlighting the increased flexibility of dsRNA at elevated temperatures. Furthermore, our data revealed a significant anisotropy in dsRNA bending elasticity, which diminished with increasing temperature, attributable to marked disparities in tilt and roll stiffness components. We delineated the underlying biophysical mechanisms and corroborated our findings with extant literature. These observations offer salient implications for advancing our understanding of nucleic acid elasticity, and are pertinent to potential medical applications.

Polycations as Aptamer-Binding Modulators for Sensitive Fluorescence Anisotropy Assay of Aflatoxin B1

Fluorescence induced by the excitation of a fluorophore with plane-polarized light has a different polarization depending on the size of the fluorophore-containing reagent and the rate of its rotation. Based on this effect, many analytical systems have been implemented in which an analyte contained in a sample and labeled with a fluorophore (usually fluorescein) competes to bind to antibodies. Replacing antibodies in such assays with aptamers, low-cost and stable oligonucleotide receptors, is complicated because binding a fluorophore to them causes a less significant change in the polarization of emissions. This work proposes and characterizes the compounds of the reaction medium that improve analyte binding and reduce the mobility of the aptamer-fluorophore complex, providing a higher analytical signal and a lower detection limit. This study was conducted on aflatoxin B1 (AFB1), a ubiquitous toxicant contaminating foods of plant origins. Eight aptamers specific to AFB1 with the same binding site and different regions stabilizing their structures were compared for affinity, based on which the aptamer with 38 nucleotides in length was selected. The polymers that interact reversibly with oligonucleotides, such as poly-L-lysine and polyethylene glycol, were tested. It was found that they provide the desired reduction in the depolarization of emitted light as well as high concentrations of magnesium cations. In the selected optimal medium, AFB1 detection reached a limit of 1 ng/mL, which was 12 times lower than in the tris buffer commonly used for anti-AFB1 aptamers. The assay time was 30 min. This method is suitable for controlling almond samples according to the maximum permissible levels of their contamination by AFB1. The proposed approach could be applied to improve other aptamer-based analytical systems.

Systematic assessment of the flexibility of uracil damaged DNA

Uracil is a common DNA lesion which is recognized and removed by uracil DNA-glycosylase (UDG) as a part of the base excision repair pathway. Excision proceeds by base flipping, and UDG efficiency is thought to depend on the ease of deformability of the bases neighboring the lesion. We used molecular dynamics simulations to assess the flexibility of a large library of dsDNA strands, containing all tetranucleotide motifs with U:A, U:G, T:A or C:G base pairs. Our study demonstrates that uracil damaged DNA largely follows trends in flexibility of undamaged DNA. Measured bending persistence lengths, groove widths, step parameters and base flipping propensities demonstrate that uracil increases the flexibility of DNA, and that U:G base paired strands are more flexible than U:A strands. Certain sequence contexts are more deformable than others, with a key role for the 3' base next to uracil. Flexibilities are large when this base is an A or G, and repressed for a C or T. A 5' T adjacent to the uracil strongly promotes flexibility, but other 5' bases are less influential. DNA bending is correlated to step deformations and base flipping, and bending aids flipping. Our study implies that the link between substrate flexibility and UDG efficiency is widely valid, helps explain why UDG prefers to bind U:G base paired strands, and suggests that the DNA bending angle of the UDG-substrate complex is optimal for base flipping.Communicated by Ramaswamy H. Sarma.

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The subtle differences in the chemical structures of double-stranded (ds) RNA and DNA lead to significant variations in their biological roles and medical implications, largely due to their distinct biophysical properties, such as bending stiffness. Although it is well known that A-form dsRNA is stiffer than B-form dsDNA under physiological salt conditions, the underlying cause of this difference remains unclear. In this study, we employ high-precision magnetic-tweezer experiments along with molecular dynamics simulations and reveal that the relative bending stiffness between dsRNA and dsDNA is primarily determined by the structure- and salt-concentration-dependent ion distribution around their helical structures. At near-physiological salt conditions, dsRNA shows a sparser ion distribution surrounding its phosphate groups compared to dsDNA, causing its greater stiffness. However, at very high monovalent salt concentrations, phosphate groups in both dsRNA and dsDNA become fully neutralized by excess ions, resulting in a similar intrinsic bending persistence length of approximately 39 nm. This similarity in intrinsic bending stiffness of dsRNA and dsDNA is coupled to the analogous fluctuations in their total groove widths and further coupled to the similar fluctuation of base-pair inclination, despite their distinct A-form and B-form helical structures.

Melting and Bubble Formation in a Double-Stranded DNA: Microscopic Aspects of Early Base-Pair Opening Events and the Role of Water

Despite its rigid structure, DNA is a remarkably flexible molecule. Flexibility is essential for biological functions (such as transcription and gene repair), which require large-amplitude structural changes such as bubble formation. The bubbles thus formed are required to have a certain stability of their own and survive long on the time scale of molecular motions. A molecular understanding of fluctuations leading to quasi-stable structures is not available. Through extensive atomistic molecular dynamics simulations, we identify a sequence of microscopic events that culminate in local bubble formation, which is initiated by base-pair (BP) opening, resulting from the cleavage of native BP hydrogen bonds (HBs). This is followed by the formation of mismatched BPs with non-native contacts. These metastable structures can either revert to their original forms or undergo a flipping transition to form a local bubble that can span across 3-4 BPs. A substantial distortion of the DNA backbone and a disruption of BP stacking are observed because of the structural changes induced by these local perturbations. We also explored how water helps in the entire process. A small number of water molecules undergo rearrangement to stabilize the intermediate states by forming HBs with DNA bases. Water thus acts as a lubricant that counteracts the enthalpic penalty suffered from the loss of native BP contacts. Although the process of bubble formation is reversible, the sequence of steps involved poses an entropic barrier, preventing it from easily retracing the path to the native state.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

Temperature-dependent elasticity of DNA, RNA, and hybrid double helices

Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.

DNA Hyperstructure

This study presents a new procedure to condense DNA molecules and precipitate them onto a glass slide. The resulting DNA molecules undergo autonomous self-assembly, creating closed superstructures on the micrometer scale, which are called DNA hyperstructures. These structures can be observed using low-magnification (4×) light microscopy. Precisely controlling the alcohol/glacial acetic acid ratio and DNA concentration during precipitation enabled the regulation of structure compaction on the slide. The alcohol/glacial acetic acid ratio is inversely proportional to the DNA concentration to achieve optimal compaction on the slide. Confocal microscopy fluorescence analysis of DNA extracts stained with DAPI shows that nucleic acids self-assemble to form structures during precipitation on the slide. This methodology is relevant since it facilitates the precipitation and visualization of DNA, regardless of its origin or molecular weight. To confirm its versatility, results with DNA extracted from human peripheral blood, the Lambda virus, and plasmid pBR322 are presented. The study examined the morphological features of DNA hyperstructures in both healthy individuals and those diagnosed with different medical conditions or illnesses, revealing distinct patterns specific to each case. This innovative technology has potential for disease detection in peripheral blood samples, ranging from cancer and Alzheimer's disease to determining the gender of the gestational product at an early stage.

Predicting DNA structure using a deep learning method

Understanding the mechanisms of protein-DNA binding is critical in comprehending gene regulation. Three-dimensional DNA structure, also described as DNA shape, plays a key role in these mechanisms. In this study, we present a deep learning-based method, Deep DNAshape, that fundamentally changes the current k-mer based high-throughput prediction of DNA shape features by accurately accounting for the influence of extended flanking regions, without the need for extensive molecular simulations or structural biology experiments. By using the Deep DNAshape method, DNA structural features can be predicted for any length and number of DNA sequences in a high-throughput manner, providing an understanding of the effects of flanking regions on DNA structure in a target region of a sequence. The Deep DNAshape method provides access to the influence of distant flanking regions on a region of interest. Our findings reveal that DNA shape readout mechanisms of a core target are quantitatively affected by flanking regions, including extended flanking regions, providing valuable insights into the detailed structural readout mechanisms of protein-DNA binding. Furthermore, when incorporated in machine learning models, the features generated by Deep DNAshape improve the model prediction accuracy. Collectively, Deep DNAshape can serve as versatile and powerful tool for diverse DNA structure-related studies.

Topology- and size-dependent binding of DNA nanostructures to the DNase I

The susceptibility of DNA nanomaterials to enzymatic degradation in biological environments is a significant obstacle limiting their broad applications in biomedicine. While DNA nanostructures exhibit some resistance to nuclease degradation, the underlying mechanism of this resistance remains elusive. In this study, the interaction of tetrahedral DNA nanostructures (TDNs) and double-stranded DNA (dsDNA) with DNase I is investigated using all-atom molecular dynamics simulations. Our results indicate that DNase I can effectively bind to all dsDNA molecules, and certain key residues strongly interact with the nucleic bases of DNA. However, the binding of DNase I to TDNs exhibits a non-monotonic behavior based on size; TDN15 and TDN26 interact weakly with DNase I (∼ - 75 kcal/mol), whereas TDN21 forms a strong binding with DNase I (∼ - 110 kcal/mol). Furthermore, the topological properties of the DNA nanostructures are analyzed, and an under-twisting (∼32°) of the DNA helix is observed in TDN15 and TDN26. Importantly, this under-twisting results in an increased width of the minor groove in TDN15 and TDN26, which primarily explains their reduced binding affinity to DNase I comparing to the dsDNA. Overall, this study demonstrated a novel mechanism for local structural control of DNA at the nanoscale by adjusting the twisting induced by length.

A genome-wide comprehensive analysis of nucleosome positioning in yeast

In eukaryotic cells, the one-dimensional DNA molecules need to be tightly packaged into the spatially constraining nucleus. Folding is achieved on its lowest level by wrapping the DNA around nucleosomes. Their arrangement regulates other nuclear processes, such as transcription and DNA repair. Despite strong efforts to study nucleosome positioning using Next Generation Sequencing (NGS) data, the mechanism of their collective arrangement along the gene body remains poorly understood. Here, we classify nucleosome distributions of protein-coding genes in Saccharomyces cerevisiae according to their profile similarity and analyse their differences using functional Principal Component Analysis. By decomposing the NGS signals into their main descriptive functions, we compared wild type and chromatin remodeler-deficient strains, keeping position-specific details preserved whilst considering the nucleosome arrangement as a whole. A correlation analysis with other genomic properties, such as gene size and length of the upstream Nucleosome Depleted Region (NDR), identified key factors that influence the nucleosome distribution. We reveal that the RSC chromatin remodeler-which is responsible for NDR maintenance-is indispensable for decoupling nucleosome arrangement within the gene from positioning outside, which interfere in rsc8-depleted conditions. Moreover, nucleosome profiles in chd1Δ strains displayed a clear correlation with RNA polymerase II presence, whereas wild type cells did not indicate a noticeable interdependence. We propose that RSC is pivotal for global nucleosome organisation, whilst Chd1 plays a key role for maintaining local arrangement.

Optical Limiting Properties of DNA Biopolymer Doped with Natural Dyes

The high-power lasers have important implications for present and future light-based technologies; therefore, the protection measures against their high-intensity radiation are extremely important. Currently, a great deal of interest is directed towards the development of new nonlinear optical materials for passive optical limiters, which are used to protect the human eye and sensitive optical and optoelectronic devices from laser-induced damage. Biopolymers doped with natural dyes are emerging as a new class of optical materials with interesting photosensitive properties. In this paper, the optical limiting capability of deoxyribonucleic acid bio-polymer functionalized with Turmeric natural dye has been demonstrated for the first time, to the best of our knowledge. The experimental investigation of the optical limit has been done by the Intensity-scan method in the NIR spectral domain at the important telecommunication wavelength of 1550 nm, using ultrashort laser pulses (~120 fs). Several optical properties of this natural dye are presented and discussed. The values of the optical transmittance in the linear regime, the saturation intensity of the nonlinear transmittance curves, and the coefficient of the nonlinear absorption have been determined. The influence of the DNA biopolymer and natural dye concentration on the optical limiting properties of the investigated biomaterials is reported and discussed. The photostability and thermal stability of the investigated solutions have also been evaluated by monitoring the temporal decay of the normalized absorption spectra under illumination with UVA light and heating, respectively. Our results evidence the positive influence of the DNA, which embeds Turmeric natural dye, on the optical limiting functionality itself and on the photostability and thermal stability of this novel material. The performed study reveals the potential of the investigated novel biomaterial for applications in nonlinear photonics, in particular in optical limiting.

Predictions of DNA mechanical properties at a genomic scale reveal potentially new functional roles of DNA flexibility

Mechanical properties of DNA have been implied to influence many of its biological functions. Recently, a new high-throughput method, called loop-seq, which allows measuring the intrinsic bendability of DNA fragments, has been developed. Using loop-seq data, we created a deep learning model to explore the biological significance of local DNA flexibility in a range of different species from different kingdoms. Consistently, we observed a characteristic and largely dinucleotide-composition-driven change of local flexibility near transcription start sites. In the presence of a TATA-box, a pronounced peak of high flexibility can be observed. Furthermore, depending on the transcription factor investigated, flanking-sequence-dependent DNA flexibility was identified as a potential factor influencing DNA binding. Compared to randomized genomic sequences, depending on species and taxa, actual genomic sequences were observed both with increased and lowered flexibility. Furthermore, in Arabidopsis thaliana, mutation rates, both de novo and fixed, were found to be associated with relatively rigid sequence regions. Our study presents a range of significant correlations between characteristic DNA mechanical properties and genomic features, the significance of which with regard to detailed molecular relevance awaits further theoretical and experimental exploration.

Deep DNAshape: Predicting DNA shape considering extended flanking regions using a deep learning method

Understanding the mechanisms of protein-DNA binding is critical in comprehending gene regulation. Three-dimensional DNA shape plays a key role in these mechanisms. In this study, we present a deep learning-based method, Deep DNAshape, that fundamentally changes the current k -mer based high-throughput prediction of DNA shape features by accurately accounting for the influence of extended flanking regions, without the need for extensive molecular simulations or structural biology experiments. By using the Deep DNAshape method, refined DNA shape features can be predicted for any length and number of DNA sequences in a high-throughput manner, providing a deeper understanding of the effects of flanking regions on DNA shape in a target region of a sequence. Deep DNAshape method provides access to the influence of distant flanking regions on a region of interest. Our findings reveal that DNA shape readout mechanisms of a core target are quantitatively affected by flanking regions, including extended flanking regions, providing valuable insights into the detailed structural readout mechanisms of protein-DNA binding. Furthermore, when incorporated in machine learning models, the features generated by Deep DNAshape improve the model prediction accuracy. Collectively, Deep DNAshape can serve as a versatile and powerful tool for diverse DNA structure-related studies.

Molecular Mechanisms of DNA Sequence Selectivity in V(D)J Recombination

Antigen receptor (AgR) diversity is central to the ability of adaptive immunity in jawed vertebrates to protect against pathogenic agents. The production of highly diverse AgR repertoires is initiated during B and T cell lymphopoiesis by V(D)J recombination, which assembles the receptor genes from component gene segments in a cut-and-paste recombination reaction. Recombination activating proteins, RAG1 and RAG2 (RAG1/2), catalyze V(D)J recombination by cleaving adjacent to recombination signal sequences (RSSs) that flank AgR 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, and the molecular mechanism for semiselective V(D)J recombination specificity is unknown. The modulation of chromatin structure during V(D)J recombination is essential in the formation of diverse AgRs in adaptive immunity while also reducing the likelihood for off-target recombination events that can result in chromosomal aberrations and genomic instability. Here we review what is presently known regarding mechanisms that facilitate assembly of RAG1/2 with RSSs, the ensuing conformational changes required for DNA cleavage activity, and how the readout of the RSS sequence affects reaction efficiency.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Force-dependent elasticity of nucleic acids

The functioning of double-stranded (ds) nucleic acids (NAs) in cellular processes is strongly mediated by their elastic response. These processes involve proteins that interact with dsDNA or dsRNA and distort their structures. The perturbation of the elasticity of NAs arising from these deformations is not properly considered by most theoretical frameworks. In this work, we introduce a novel method to assess the impact of mechanical stress on the elastic response of dsDNA and dsRNA through the analysis of the fluctuations of the double helix. Application of this approach to atomistic simulations reveals qualitative differences in the force dependence of the mechanical properties of dsDNA with respect to those of dsRNA, which we relate to structural features of these molecules by means of physically-sound minimalistic models.

Rigid Base Biasing in Molecular Dynamics Enables Enhanced Sampling of DNA Conformations

All-atom simulations have become increasingly popular to study conformational and dynamical properties of nucleic acids as they are accurate and provide high spatial and time resolutions. This high resolution, however, comes at a heavy computational cost, and, within the time scales of simulations, nucleic acids weakly fluctuate around their ideal structure exploring a limited set of conformations. We introduce the RBB-NA algorithm (available as a package in the Open Source Library PLUMED), which is capable of controlling rigid base parameters in all-atom simulations of nucleic acids. With suitable biasing potentials, this algorithm can "force" a DNA or RNA molecule to assume specific values of the six rotational (tilt, roll, twist, buckle, propeller, opening) and/or the six translational parameters (shift, slide, rise, shear, stretch, stagger). The algorithm enables the use of advanced sampling techniques to probe the structure and dynamics of locally strongly deformed nucleic acids. We illustrate its performance showing some examples in which DNA is strongly twisted, bent, or locally buckled. In these examples, RBB-NA reproduces well the unconstrained simulations data and other known features of DNA mechanics, but it also allows one to explore the anharmonic behavior characterizing the mechanics of nucleic acids in the high deformation regime.

Efficient modification and preparation of circular DNA for expression in cell culture

DNA plasmids are an essential tool for delivery and expression of RNAs and proteins in cell culture experiments. The preparation of plasmids typically involves a laborious process of bacterial cloning, validation, and purification. While the expression plasmids can be designed and ordered from the contract manufacturers, the cost may be prohibitive when a large number of plasmids is required. We have developed an efficient fully synthetic method and protocol that enables the production of circularized DNA containing expression elements ready for transfection in as little as 3 hours, thereby eliminating the bacterial cloning steps. The protocol describes how to take a linear double-stranded DNA fragment and efficiently circularize and purify this DNA fragment with minimal hands-on time. As proof of the principle, we applied Circular Vector expressing engineered prime editing guide RNA (epegRNA) in cell culture, and demonstrated matching and even exceeding performance of this method as compared to guides expressed by plasmids. The method's speed of preparation, low cost, and ease of use will make it a useful tool in applications requiring the expression of short RNAs and proteins.

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.

Accurate Quantum-Mechanically Derived Force-Fields through a Fragment-Based Approach: Balancing Specificity and Transferability in the Prediction of Self-Assembly in Soft Matter

The wide range of time/length scales covered by self-assembly in soft matter makes molecular dynamics (MD) the ideal candidate for simulating such a supramolecular phenomenon at an atomistic level. However, the reliability of MD outcomes heavily relies on the accuracy of the adopted force-field (FF). The spontaneous re-ordering in liquid crystalline materials stands as a clear example of such collective self-assembling processes, driven by a subtle and delicate balance between supramolecular interactions and single-molecule flexibility. General-purpose transferable FFs often dramatically fail to reproduce such complex phenomena, for example, the error on the transition temperatures being larger than 100 K. Conversely, quantum-mechanically derived force-fields (QMD-FFs), specifically tailored for the target system, were recently shown (J. Phys. Chem. Lett.2022,13, 243) to allow for the required accuracy as they not only well reproduced transition temperatures but also yielded a quantitative agreement with the experiment on a wealth of structural, dynamic, and thermodynamic properties. The main drawback of this strategy stands in the computational burden connected to the numerous quantum mechanical (QM) calculations usually required for a target-specific parameterization, which has undoubtedly hampered the routine application of QMD-FFs. In this work, we propose a fragment-based strategy to extend the applicability of QMD-FFs, in which the amount of QM calculations is significantly reduced, being a single-molecule-optimized geometry and its Hessian matrix the only QM information required. To validate this route, a new FF is assembled for a large mesogen, exploiting the parameters obtained for two smaller liquid crystalline molecules, in this and previous work. Lengthy MD simulations are carried out with the new transferred QMD-FF, observing again a spontaneous re-orientation in the correct range of temperatures, with good agreement with the available experimental measures. The present results strongly suggest that a partial transfer of QMD-FF parameters can be invoked without a significant loss of accuracy, thus paving the way to exploit the method's intrinsic predictive capabilities in the simulation of novel soft materials.

The convergent xenogeneic silencer MucR predisposes α-proteobacteria to integrate AT-rich symbiosis genes

Bacterial adaptation is largely shaped by horizontal gene transfer, xenogeneic silencing mediated by lineage-specific DNA bridgers (H-NS, Lsr2, MvaT and Rok), and various anti-silencing mechanisms. No xenogeneic silencing DNA bridger is known for α-proteobacteria, from which mitochondria evolved. By investigating α-proteobacterium Sinorhizobium fredii, a facultative legume microsymbiont, here we report the conserved zinc-finger bearing MucR as a novel xenogeneic silencing DNA bridger. Self-association mediated by its N-terminal domain (NTD) is required for DNA–MucR–DNA bridging complex formation, maximizing MucR stability, transcriptional silencing, and efficient symbiosis in legume nodules. Essential roles of NTD, CTD (C-terminal DNA-binding domain), or full-length MucR in symbiosis can be replaced by non-homologous NTD, CTD, or full-length protein of H-NS from γ-proteobacterium Escherichia coli, while NTD rather than CTD of Lsr2 from Gram-positive Mycobacterium tuberculosis can replace the corresponding domain of MucR in symbiosis. Chromatin immunoprecipitation sequencing reveals similar recruitment profiles of H-NS, MucR and various functional chimeric xenogeneic silencers across the multipartite genome of S. fredii, i.e. preferring AT-rich genomic islands and symbiosis plasmid with key symbiosis genes as shared targets. Collectively, the convergently evolved DNA bridger MucR predisposed α-proteobacteria to integrate AT-rich foreign DNA including symbiosis genes, horizontal transfer of which is strongly selected in nature.

Sequence-Specific Features of Short Double-Strand, Blunt-End RNAs Have RIG-I- and Type 1 Interferon-Dependent or -Independent Anti-Viral Effects

Pathogen-associated molecular patterns, including cytoplasmic DNA and double-strand (ds)RNA trigger the induction of interferon (IFN) and antiviral states protecting cells and organisms from pathogens. Here we discovered that the transfection of human airway cell lines or non-transformed fibroblasts with 24mer dsRNA mimicking the cellular micro-RNA (miR)29b-1* gives strong anti-viral effects against human adenovirus type 5 (AdV-C5), influenza A virus X31 (H3N2), and SARS-CoV-2. These anti-viral effects required blunt-end complementary RNA strands and were not elicited by corresponding single-strand RNAs. dsRNA miR-29b-1* but not randomized miR-29b-1* mimics induced IFN-stimulated gene expression, and downregulated cell adhesion and cell cycle genes, as indicated by transcriptomics and IFN-I responsive Mx1-promoter activity assays. The inhibition of AdV-C5 infection with miR-29b-1* mimic depended on the IFN-alpha receptor 2 (IFNAR2) and the RNA-helicase retinoic acid-inducible gene I (RIG-I) but not cytoplasmic RNA sensors MDA5 and ZNFX1 or MyD88/TRIF adaptors. The antiviral effects of miR29b-1* were independent of a central AUAU-motif inducing dsRNA bending, as mimics with disrupted AUAU-motif were anti-viral in normal but not RIG-I knock-out (KO) or IFNAR2-KO cells. The screening of a library of scrambled short dsRNA sequences identified also anti-viral mimics functioning independently of RIG-I and IFNAR2, thus exemplifying the diverse anti-viral mechanisms of short blunt-end dsRNAs.

Mechanical properties of nucleic acids and the non-local twistable wormlike chain model

Mechanical properties of nucleic acids play an important role in many biological processes that often involve physical deformations of these molecules. At sufficiently long length scales (say, above ∼20-30 base pairs), the mechanics of DNA and RNA double helices is described by a homogeneous Twistable Wormlike Chain (TWLC), a semiflexible polymer model characterized by twist and bending stiffnesses. At shorter scales, this model breaks down for two reasons: the elastic properties become sequence-dependent and the mechanical deformations at distal sites get coupled. We discuss in this paper the origin of the latter effect using the framework of a non-local Twistable Wormlike Chain (nlTWLC). We show, by comparing all-atom simulations data for DNA and RNA double helices, that the non-local couplings are of very similar nature in these two molecules: couplings between distal sites are strong for tilt and twist degrees of freedom and weak for roll. We introduce and analyze a simple double-stranded polymer model that clarifies the origin of this universal distal couplings behavior. In this model, referred to as the ladder model, a nlTWLC description emerges from the coarsening of local (atomic) degrees of freedom into angular variables that describe the twist and bending of the molecule. Different from its local counterpart, the nlTWLC is characterized by a length-scale-dependent elasticity. Our analysis predicts that nucleic acids are mechanically softer at the scale of a few base pairs and are asymptotically stiffer at longer length scales, a behavior that matches experimental data.

Stretching of Long Double-Stranded DNA and RNA Described by the Same Approach

We propose an approach to help interpret polymer force-extension curves that exhibit plateau regimes. When coupled to a bead-spring dynamic model, the approach accurately reproduces a variety of experimental force-extension curves of long double-stranded DNA and RNA, including torsionally constrained and unconstrained DNA and negatively supercoiled DNA. A key feature of the model is a specific nonconvex energy function of the spring. We provide an algorithm to obtain the five required parameters of the model from experimental force-extension curves. The applicability of the approach to the force-extension curves of double-stranded (ds) DNA of variable GC content as well as to a DNA/RNA hybrid structure is explored and confirmed. We use the approach to explain counterintuitive sequence-dependent trends and make predictions. In the plateau region of the force-extension curves, our molecular dynamics simulations show that the polymer separates into a mix of weakly and strongly stretched states without forming macroscopically distinct phases. The distribution of these states is predicted to depend on the sequence.

Bending of Canonical and G/T Mismatched DNAs

Mismatched base pairs alter the flexibility and intrinsic curvature of DNA. The role of such DNA features is not fully understood in the mismatch repair pathway. MutS/DNA complexes exhibit DNA bending, PHE intercalation, and changes of base-pair parameters near the mismatch. Recently, we have shown that base-pair opening in the absence of MutS can discriminate mismatches from canonical base pairs better than DNA bending. However, DNA bending in the absence of MutS was found to be rather challenging to describe correctly. Here, we present a computational study on the DNA bending of canonical and G/T mismatched DNAs. Five types of geometric parameters covering template-based bending toward the experimental DNA structure, global, and local geometry parameters were employed in biased molecular dynamics in the absence of MutS. None of these parameters showed higher discrimination than the base-pair opening. Only roll could induce a sharply localized bending of DNA as observed in the experimental MutS/DNA structure. Further, we demonstrated that the intercalation of benzene mimicking PHE decreases the energetic cost of DNA bending without any effect on mismatch discrimination.