Computational investigation of the impact of core sequence on immobile DNA four-way junction structure and dynamics

3dDNAscoreA: A scoring function for evaluation of DNA 3D structures

DNA molecules are vital macromolecules that play a fundamental role in many cellular processes and have broad applications in medicine. For example, DNA aptamers have been rapidly developed for diagnosis, biosensors, and clinical therapy. Recently, we proposed a computational method of predicting DNA 3D structures, called 3dDNA. However, it lacks a scoring function to evaluate the predicted DNA 3D structures, and so they are not ranked for users. Here, we report a scoring function, 3dDNAscoreA, for evaluation of DNA 3D structures based on a deep learning model ARES for RNA 3D structure evaluation but using a new strategy for training. 3dDNAscoreA is benchmarked on two test sets to show its ability to rank DNA 3D structures and select the native and near-native structures.

Sculpting photoproducts with DNA origami

Natural light-harvesting systems spatially organize densely packed dyes in different configurations to either transport excitons or convert them into charge photoproducts, with high efficiency. In contrast, artificial photosystems like organic solar cells and light-emitting diodes lack this fine structural control, limiting their efficiency. Thus, biomimetic multi-dye systems are needed to organize dyes with the sub-nanometer spatial control required to sculpt resulting photoproducts. Here, we synthesize 11 distinct perylene diimide (PDI) dimers integrated into DNA origami nanostructures and identify dimer architectures that offer discrete control over exciton transport versus charge separation. The large structural-space and site-tunability of origami uniquely provides controlled PDI dimer packing to form distinct excimer photoproducts, which are sensitive to interdye configurations. In the future, this platform enables large-scale programmed assembly of dyes mimicking natural systems to sculpt distinct photophysical products needed for a broad range of optoelectronic devices, including solar energy converters and quantum information processors.

Site-Specific Investigation of DNA Holliday Junction Dynamics and Structure with 6-Methylisoxanthopterin, a Fluorescent Guanine Analog

DNA Holliday Junction (HJ) formation and resolution is requisite for maintaining genomic stability in processes such as replication fork reversal and double-strand break repair. If HJs are not resolved, chromosome disjunction and aneuploidy result, hallmarks of tumor cells. To understand the structural features that lead to processing of these four-stranded joint molecule structures, we seek to identify structural and dynamic features unique to the central junction core. We incorporate the fluorescent guanine analog 6-methylisoxanthopterin (6-MI) at ten different locations throughout a model HJ structure to obtain site-specific information regarding the structure and dynamics of bases relative to those in a comparable sequence context in duplex DNA. These comparisons were accomplished through measuring fluorescence lifetime, relative brightness, fluorescence anisotropy, and thermodynamic stability, along with fluorescence quenching assays. These time-resolved and steady-state fluorescence measurements demonstrate that the structural distortions imposed by strand crossing result in increased solvent exposure, less stacking of bases and greater extrahelical nature of bases within the junction core. The 6-MI base analogs in the junction reflect these structural changes through an increase in intensity relative to those in the duplex. Molecular dynamics simulations performed using a model HJ indicate the primary sources of deformation are in the shift and twist parameters of the bases at the central junction step. These results suggest that junction-binding proteins may use the unique structure and dynamics of the bases at the core for recognition.

Prediction and Control in DNA Nanotechnology

DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.

The Rule of Thirds: Controlling Junction Chirality and Polarity in 3D DNA Tiles

The successful self-assembly of tensegrity triangle DNA crystals heralded the ability to programmably construct macroscopic crystalline nanomaterials from rationally-designed, nanoscale components. This 3D DNA tile owes its "tensegrity" nature to its three rotationally stacked double helices locked together by the tensile winding of a center strand segmented into 7 base pair (bp) inter-junction regions, corresponding to two-thirds of a helical turn of DNA. All reported tensegrity triangles to date have employed ( Z + 2 / 3 ) \[\left( {Z{\bm{ + }}2{\bf /}3} \right)\] turn inter-junction segments, yielding right-handed, antiparallel, "J1" junctions. Here a minimal DNA triangle motif consisting of 3-bp inter-junction segments, or one-third of a helical turn is reported. It is found that the minimal motif exhibits a reversed morphology with a left-handed tertiary structure mediated by a locally-parallel Holliday junction-the "L1" junction. This parallel junction yields a predicted helical groove matching pattern that breaks the pseudosymmetry between tile faces, and the junction morphology further suggests a folding mechanism. A Rule of Thirds by which supramolecular chirality can be programmed through inter-junction DNA segment length is identified. These results underscore the role that global topological forces play in determining local DNA architecture and ultimately point to an under-explored class of self-assembling, chiral nanomaterials for topological processes in biological systems.

Molecular Dynamic Studies of Dye-Dye and Dye-DNA Interactions Governing Excitonic Coupling in Squaraine Aggregates Templated by DNA Holliday Junctions

Dye molecules, arranged in an aggregate, can display excitonic delocalization. The use of DNA scaffolding to control aggregate configurations and delocalization is of research interest. Here, we applied Molecular Dynamics (MD) to gain an insight on how dye-DNA interactions affect excitonic coupling between two squaraine (SQ) dyes covalently attached to a DNA Holliday junction (HJ). We studied two types of dimer configurations, i.e., adjacent and transverse, which differed in points of dye covalent attachments to DNA. Three structurally different SQ dyes with similar hydrophobicity were chosen to investigate the sensitivity of excitonic coupling to dye placement. Each dimer configuration was initialized in parallel and antiparallel arrangements in the DNA HJ. The MD results, validated by experimental measurements, suggested that the adjacent dimer promotes stronger excitonic coupling and less dye-DNA interaction than the transverse dimer. Additionally, we found that SQ dyes with specific functional groups (i.e., substituents) facilitate a closer degree of aggregate packing via hydrophobic effects, leading to a stronger excitonic coupling. This work advances a fundamental understanding of the impacts of dye-DNA interactions on aggregate orientation and excitonic coupling.

Understanding the self-assembly dynamics of A/T absent 'four-way DNA junctions with sticky ends' at altered physiological conditions through molecular dynamics simulations

Elucidation of structure and dynamics of alternative higher-order structures of DNA such as in branched form could be targeted for therapeutics designing. Herein, we are reporting the intrinsically dynamic and folds transitions of an unusual DNA junction with sequence d(CGGCGGCCGC)4 which self-assembles into a four-way DNA junction form with sticky ends using long interval molecular simulations under various artificial physiological conditions. The original crystal structure coordinates (PDB ID: 3Q5C) for the selected DNA junction was considered for a total of 1.1 μs molecular dynamics simulation interval, including different temperature and pH, under OPLS-2005 force field using DESMOND suite. Following, post-dynamics structure parameters for the DNA junction were calculated and analyzed by comparison to the crystal structure. We show here that the self-assembly dynamics of DNA junction is mitigated by the temperature and pH sensitivities, and discloses peculiar structural properties as function of time. From this study it can be concluded on account of temperature sensitive and pH dependent behaviours, DNA junction periodic arrangements can willingly be synthesized and redeveloped for multiple uses like genetic biomarkers, DNA biosensor, DNA nanotechnology, DNA Zipper, etc. Furthermore, the pH dis-regulation behaviour may be used to trigger the functionality of DNA made drug-releasing nanomachines.

Exciton Chirality Inversion in Dye Dimers Templated by DNA Holliday Junction

While only one enantiomer of chiral biomolecules performs a biological function, access to both enantiomers (or enantiomorphs) proved to be advantageous for technology. Using dye covalent attachment to a DNA Holliday junction (HJ), we created two pairs of dimers of bis(chloroindolenine)squaraine dye that enabled strongly coupled molecular excitons of opposite chirality in solution. The exciton chirality inversion was achieved by interchanging single covalent linkers of unequal length tethering the dyes of each dimer to the HJ core. Dimers in each pair exhibited profound exciton-coupled circular dichroism (CD) couplets of opposite signs. Dimer geometries, modeled by simultaneous fitting absorption and CD spectra, were related in each pair as nonsuperimposable and nearly exact mirror images. The origin of observed exciton chirality inversion was explained in the view of isomerization of the stacked Holliday junction. This study will open new opportunities for creating excitonic DNA-based materials that rely on programmable system chirality.

Spectroscopic details on the molecular structure of pyrimidine‑2‑thiones heterocyclic compounds: computational and antiviral activity against the main protease enzyme of SARS-CoV-2

Computational tools in investigating of spectral heterocyclic compounds ranges based on pyrimidine‑2‑thiones, take some importance in identifying their molecular and electronic behavior. Some charcoal heterocyclic compounds were previously synthesized in our laboratory and their experimental results were compared with the computational evaluation. Computational spectroscopic analytical items (IR, NMR and UV–Vis) were calculated using the more popular DFT methods and the predicted results were compared with the reported experimental ones. Quantum and chemical parameters were calculated and molecular electrostatic surface potential (MEP) was studied which predicted the highly electronic sites around the compounds. Some molecular properties (ionization energy, electron affinity, energy gap, hardness, electronegativity, electrophilicity index, static dipole moment and average linear polarizability) of these Schiff bases which were computed at B3LYP/6-31G(d,p) level in aqueous phase. Benchmark analysis was performed for three ab initio functionals such B3LYP, BPV86 and B3PW91 methods to explain the data resulted from NMR spectra. The docking study of some selected previously synthesized compounds was performed using the viral Mpro enzyme protein in compared to a k36 reference ligand inhibitor. The study indicated the ability of the synthesized compounds to form H-bond and hydrophobic (VDW, π-alkyl and π-sulfur) interactions with Mpro enzyme receptor with high inhibition effect of compound L2.

Predicting accurate ab initio DNA electron densities with equivariant neural networks

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

Molecular dynamics simulations of cyanine dimers attached to DNA Holliday junctions

Dye aggregates and their excitonic properties are of interest for their applications to organic photovoltaics, non-linear optics, and quantum information systems. DNA scaffolding has been shown to be effective at promoting the aggregation of dyes in a controllable manner. Specifically, isolated DNA Holliday junctions have been used to achieve strongly coupled cyanine dye dimers. However, the structural properties of the dimers and the DNA, as well as the role of Holliday junction isomerization are not fully understood. To study the dynamics of cyanine dimers in DNA, molecular dynamics simulations were carried out for adjacent and transverse dimers attached to Holliday junctions in two different isomers. It was found that dyes attached to adjacent strands in the junction exhibit stronger dye-DNA interactions and larger inter-dye separations compared to transversely attached dimers, as well as end-to-end arrangements. Transverse dimers exhibit lower inter-dye separations and more stacked configurations. Furthermore, differences in Holliday junction isomer are analyzed and compared to dye orientations. For transverse dyes exhibiting the smaller inter-dye separations, excitonic couplings were calculated and shown to be in agreement with experiment. Our results suggested that dye attachment locations on DNA Holliday junctions affect dye-DNA interactions, dye dynamics, and resultant dye orientations which can guide the design of DNA-templated cyanine dimers with desired properties.

Molecular dynamics simulations reveal dye attachment and DNA Holliday junction isomer effects on dye dimer orientations and excitonic couplings. These simulations can guide synthesis and experiments of dye-DNA structures for excitonic applications.