DNA under Force: Mechanics, Electrostatics, and Hydration

Supramolecular assembly of phenanthrene-DNA conjugates into light-harvesting nanospheres

The self-assembly of highly functionalized phenanthrene-DNA conjugates into supramolecular nanostructures is presented. DNA oligomers modified with phenanthrene residues at the 3'-end and internal positions self-assemble into spherical nanostructures. The nanospheres exhibit light-harvesting properties. Upon irradiation of phenanthrene, the excitation energy is transferred along phenanthrene units, resulting in phenanthrene-pyrene exciplex formation.

Developing an Implicit Solvation Machine Learning Model for Molecular Simulations of Ionic Media

Molecular dynamics (MD) simulations of biophysical systems require accurate modeling of their native environment, i.e., aqueous ionic solution, as it critically impacts the structure and function of biomolecules. On the other hand, the models should be computationally efficient to enable simulations of large spatiotemporal scales. Here, we present the deep implicit solvation model for sodium chloride solutions that satisfies both requirements. Owing to the use of the neural network potential, the model can capture the many-body potential of mean force, while the implicit water treatment renders the model inexpensive. We demonstrate our approach first for pure ionic solutions with concentrations ranging from physiological to 2 M. We then extend the model to capture the effective ion interactions in the vicinity and far away from a DNA molecule. In both cases, the structural properties are in good agreement with all-atom MD, showcasing a general methodology for the efficient and accurate modeling of ionic media.

Supramolecular assembly of pyrene-DNA conjugates: influence of pyrene substitution pattern and implications for artificial LHCs

The supramolecular self-assembly of pyrene-DNA conjugates into nanostructures is presented. DNA functionalized with different types of pyrene isomers at the 3'-end self-assemble into nano-objects. The shape of the nanostructures is influenced by the type of pyrene isomer appended to the DNA. Multilamellar vesicles are observed with the 1,6- and 1,8-isomers, whereas conjugates of the 2,7-isomer exclusively assemble into spherical nanoparticles. Self-assembled nano-spheres obtained with the 2,7-dialkynyl pyrene isomer were used for the construction of an artificial light-harvesting complex (LHC) in combination with Cy3 as the energy acceptor.

Sequence-Dependent Orientational Coupling and Electrostatic Attraction in Cation-Mediated DNA-DNA Interactions

Condensation of DNA is vital for its biological functions and controlled nucleic acid assemblies. However, the mechanisms of DNA condensation are not fully understood due to the inability of experiments to access cation distributions and the complex interplay of energetic and entropic forces during assembly. By constructing free energy surfaces using exhaustive sampling and detailed analysis of cation distributions, we elucidate the mechanism of DNA condensation in different salt conditions and with different DNA sequences. We found that DNA condensation is facilitated by the correlated dynamics of the localized cations at the grooves of DNA helices. These dynamics are strongly dependent on the salt conditions and DNA sequences. In the presence of magnesium ions, major groove binding facilitates attraction. In contrast, in the presence of polyvalent cations, minor groove binding serves to create charge patterns, leading to condensation. Our findings present a novel advancement in the field and have broad implications for understanding and controlling nucleic acid complexes in vivo and in vitro.

Unfolding of crumpled thin sheets

Crumpled thin sheets are complex fractal structures whose physical properties are influenced by a hierarchy of ridges. In this paper, we report experiments that measure the stress-strain relation and show the coexistence of phases in the stretching of crumpled surfaces. The pull stress showed a change from a linear Hookean regime to a sublinear scaling with an exponent of 0.65±0.03, which is identified with the Hurst exponent of the crumpled sheets. The stress fluctuations are studied. The statistical distribution of force peaks is analyzed. It is shown that the unpacking of crumpled sheets is guided by a hierarchical order of ridges.

Complex DNA Architectonics─Self-Assembly of Amphiphilic Oligonucleotides into Ribbons, Vesicles, and Asterosomes

The precise arrangement of structural subunits is a key factor for the proper shape and function of natural and artificial supramolecular assemblies. In DNA nanotechnology, the geometrically well-defined double-stranded DNA scaffold serves as an element of spatial control for the precise arrangement of functional groups. Here, we describe the supramolecular assembly of chemically modified DNA hybrids into diverse types of architectures. An amphiphilic DNA duplex serves as the sole structural building element of the nanosized supramolecular structures. The morphology of the assemblies is governed by a single subunit of the building block. The chemical nature of this subunit, i.e., polyethylene glycols of different chain length or a carbohydrate moiety, exerts a dramatic influence on the architecture of the assemblies. Cryo-electron microscopy revealed the arrangement of the individual DNA duplexes within the different constructs. Thus, the morphology changes from vesicles to ribbons with increasing length of a linear polyethylene glycol. Astoundingly, attachment of a N-acetylgalactosamine carbohydrate to the DNA duplex moiety produces an unprecedented type of star-shaped architecture. The novel DNA architectures presented herein imply an extension of the current concept of DNA materials and shed new light on the fast-growing field of DNA nanotechnology.

High-sensitivity small-molecule detection of microcystin-LR cyano-toxin using a terahertz-aptamer biosensor

We demonstrate the rapid and highly sensitive detection of a small molecule, microcystin-LR (MC-LR) toxin using an aptasensor based on a terahertz (THz) emission technique named the terahertz chemical microscope (TCM). The main component of the TCM is the sensing plate, which consists of a thin silicon layer deposited on a sapphire substrate, with a natural SiO₂ layer formed on the top of the Si layer. The DNA aptamer is linked to the oxidized top surface of the silicon layer by a one-step reaction (click chemistry) between the DBCO-labeled aptamer and an azido group that binds to the surface. Using density functional theory (DFT) calculations, the number of active sites on the surface has been estimated to be 3.8 × 10¹³ cm^-2. Aptamer immobilization and MC-LR binding have been optimized by adjusting the aptamer concentration and the binding buffer composition. When MC-LR binds with the DNA aptamer, it causes a change in the chemical potential at the surface of the sensing plate, which leads to a change in the amplitude of the THz signal. Compared with other bio-sensing methods such as surface plasmon resonance (SPR), TCM is a rapid assay that can be completed in 15 min (10 min incubation and 5 min data acquisition). Moreover, our results show that the aptamer-based TCM can detect MC-LR with an excellent detection limit of 50 ng L^-1, which is 20 times more sensitive compared with SPR measurements of MC-LR.

Parallel DNA Extraction From Whole Blood for Rapid Sample Generation in Genetic Epidemiological Studies

Large-scale genetic epidemiological studies require high-quality analysis of samples such as blood or saliva from multiple patients, which is challenging at the point of care. To expand these studies’ impact, minimal sample storage time and less complex extraction of a substantial quantity and good purity of DNA or RNA for downstream applications are necessary. Here, a simple microfluidics-based system that performs genomic DNA (gDNA) extraction from whole blood was developed. In this system, a mixture of blood lysate, paramagnetic beads, and binding buffer are first placed into the input well. Then, the gDNA-bound paramagnetic beads are pulled using a magnet through a central channel containing a wash buffer to the output well, which contains elution buffer. The gDNA is eluted at 55°C off the chip. The 40-minute microfluidic protocol extracts gDNA from six samples simultaneously and requires an input of 4 μL of diluted blood and a total reagent volume of 75 μL per reaction. Techniques including quantitative PCR (qPCR) and spectrofluorimetry were used to test the purity and quantity of gDNA eluted from the chip following extraction. Bead transport and molecular diffusional analysis showed that an input of less than 4 ng of gDNA (∼667 white blood cells) is optimal for on-chip extraction. There was no observable transport of inhibitors into the eluate that would greatly affect qPCR, and a sample was successfully prepared for next-generation sequencing (NGS). The microfluidics-based extraction of DNA from whole blood described here is paramount for future work in DNA-based point-of-care diagnostics and NGS library workflows.

Molecular chains under tension: Thermal and mechanical activation of statistically interacting extension and contraction particles

This work introduces a methodology for the statistical mechanical analysis of polymeric chains under tension controlled by optical or magnetic tweezers at thermal equilibrium with an embedding fluid medium. The response of single bonds between monomers or of entire groups of monomers to tension is governed by the activation of statistically interacting particles representing quanta of extension or contraction. This method of analysis is capable of describing thermal unbending of the freely jointed or wormlike chain kind, linear or nonlinear contour elasticity, and structural transformations including effects of cooperativity. The versatility of this approach is demonstrated in an application to double-stranded DNA undergoing torsionally unconstrained stretching across three regimes of mechanical response including an overstretching transition. The three-regime force-extension characteristic, derived from a single free-energy expression, accurately matches empirical evidence.

Modelling DNA extension and fragmentation in contractive microfluidic devices: a Brownian dynamics and computational fluid dynamics approach

Fragmenting DNA into short pieces is an essential manipulation in many biological studies, ranging from genome sequencing to molecular diagnosis. Among various DNA fragmentation methods, microfluidic hydrodynamic DNA fragmentation has huge advantages especially in terms of handling small-volume samples and being integrated into automatic and all-in-one DNA analysis equipment. Despite the fast progress in experimental studies and applications, a systematic understanding of how DNA molecules are distributed, stretched and fragmented in a confined microfluidic field is still lacking. In this work, we investigate the extension and fragmentation of DNA in a typical contractive microfluidic field, which consists of a shear flow-dominated area and an elongational flow-dominated area, using the Brownian dynamics-computational fluid dynamics method. Our results show that the shear flow at the straight part of the microfluidic channel and the elongational flow at the contractive bottleneck together determine the performance of DNA fragmentation. The average fragment size of DNA decreases with the increase of the strain rate of the elongational flow, and the upstream shear flow can significantly precondition the conformation of DNA to produce shorter and more uniform fragments. A systematic study of the dynamics of DNA fragmentation shows that DNA tends to break at the mid-point when the strain rate of elongational flow is small, and the breakage point largely deviates from the midpoint as the strain rate increases. Our simulation of the thorough DNA fragmentation process in a realistic microfluidic field agrees well with experimental results. We expect that our study can shed new light on the development of future microfluidic devices for DNA fragmentation and integrated DNA analysis devices.

Charged nanorods at heterogeneously charged surfaces

We study the spatial and orientational distribution of charged nanorods (rodlike counterions) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The nanorods are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude.

Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Double-stranded (ds) RNAs play essential roles in many processes of cell metabolism. The knowledge of three-dimensional (3D) structure, stability, and flexibility of dsRNAs in salt solutions is important for understanding their biological functions. In this work, we further developed our previously proposed coarse-grained model to predict 3D structure, stability, and flexibility for dsRNAs in monovalent and divalent ion solutions through involving an implicit structure-based electrostatic potential. The model can make reliable predictions for 3D structures of extensive dsRNAs with/without bulge/internal loops from their sequences, and the involvement of the structure-based electrostatic potential and corresponding ion condition can improve the predictions for 3D structures of dsRNAs in ion solutions. Furthermore, the model can make good predictions for thermal stability for extensive dsRNAs over the wide range of monovalent/divalent ion concentrations, and our analyses show that the thermally unfolding pathway of dsRNA is generally dependent on its length as well as its sequence. In addition, the model was employed to examine the salt-dependent flexibility of a dsRNA helix, and the calculated salt-dependent persistence lengths are in good accordance with experiments.

Divalent Ion-Mediated DNA-DNA Interactions: A Comparative Study of Triplex and Duplex

Ion-mediated interaction between DNAs is essential for DNA condensation, and it is generally believed that monovalent and nonspecifically binding divalent cations cannot induce the aggregation of double-stranded (ds) DNAs. Interestingly, recent experiments found that alkaline earth metal ions such as Mg²⁺ can induce the aggregation of triple-stranded (ts) DNAs, although there is still a lack of deep understanding of the surprising findings at the microscopic level. In this work, we employed all-atom dynamic simulations to directly calculate the potentials of mean force (PMFs) between tsDNAs, between dsDNAs, and between tsDNA and dsDNA in Mg²⁺ solutions. Our calculations show that the PMF between tsDNAs is apparently attractive and becomes more strongly attractive at higher [Mg²⁺], although the PMF between dsDNAs cannot become apparently attractive even at high [Mg²⁺]. Our analyses show that Mg²⁺ internally binds into grooves and externally binds to phosphate groups for both tsDNA and dsDNA, whereas the external binding of Mg²⁺ is much stronger for tsDNA. Such stronger external binding of Mg²⁺ for tsDNA favors more apparent ion-bridging between helices than for dsDNA. Furthermore, our analyses illustrate that bridging ions, as a special part of external binding ions, are tightly and positively coupled to ion-mediated attraction between DNAs.