Hydrogen bonding, base stacking, and steric effects in dna replication

Xeno Nucleic Acids as Functional Materials: From Biophysical Properties to Application

Xeno nucleic acid (XNA) are artificial nucleic acids, in which the chemical composition of the sugar moiety is changed. These modifications impart distinct physical and chemical properties to XNAs, leading to changes in their biological, chemical, and physical stability. Additionally, these alterations influence the binding dynamics of XNAs to their target molecules. Consequently, XNAs find expanded applications as functional materials in diverse fields. This review provides a comprehensive summary of the distinctive biophysical properties exhibited by various modified XNAs and explores their applications as innovative functional materials in expanded fields.

Hydration Waters Make Up for the Missing Third Hydrogen Bond in the A·T Base Pair

Base pairing complementarity is central to DNA function. G·C and A·T pair specificity is thought to originate from the different number of hydrogen bonds the pairs make. Quantifying how many hydrogen bonds exist can be difficult because water molecules in the surrounding can make up for or disrupt direct hydrogen bonds, and the hydration structures around A·T and G·C pairs on duplex DNA are distinct. Large-scale computer simulations have been used here to create a detailed map for the hydration structure on A·T and G·C base pairs in water. The contributions of specific hydration waters to the free energy of each of the hydrogen bonds in the A·T and G·C pairs were computed. Using the equilibrium fractions of hydrated versus unhydrated states from the hydration profiles, the impact of specific bound waters on each hydrogen bond can be uniquely quantified using a thermodynamic construction. The findings suggest that hydration water in the minor groove of an A·T pair can provide up to about 2 kcal/mol of free energy advantage, effectively making up for the missing third hydrogen bond in the A·T pair compared to G·C, rendering the intrinsic thermodynamic stability of the A·T pair almost synonymous with G·C.

Theoretical insights into the photophysics of an unnatural base Z: A MS-CASPT2 investigation

We have employed the highly accurate multistate complete active space second-order perturbation theory (MS-CASPT2) method to investigate the photoinduced excited state relaxation properties of one unnatural base, namely Z. Upon excitation to the S2 state of Z, the internal conversion to the S1 state would be dominant. From the S1 state, two intersystem crossing paths leading to the T2 and T1 states and one internal conversion path to the S0 state are possible. However, considering the large barrier to access the S1 /S0 conical intersection and the strong spin-orbit coupling between S1 and T2 states (>40 cm-1 ), the intersystem crossing to the triplet manifolds is predicted to be more preferred. Arriving at the T2 state, the internal conversion to the T1 state and the intersystem crossing back to the S1 state are both possible considering the S1 /T2 /T1 three-state intersection near the T2 minimum. Upon arrival at the T1 state, the deactivation to S0 can be efficient after overcoming a small barrier to access T1 /S0 crossing point, where the spin-orbit coupling (SOC) is as large as 39.7 cm-1 . Our present work not only provides in-depth insights into the photoinduced process of unnatural base Z, but can also help the future design of novel unnatural bases with better photostability.

Single-molecule analysis of DNA base-stacking energetics using patterned DNA nanostructures

The DNA double helix structure is stabilized by base-pairing and base-stacking interactions. However, a comprehensive understanding of dinucleotide base-stacking energetics is lacking. Here we combined multiplexed DNA-based point accumulation in nanoscale topography (DNA-PAINT) imaging with designer DNA nanostructures and measured the free energy of dinucleotide base stacking at the single-molecule level. Multiplexed imaging enabled us to extract the binding kinetics of an imager strand with and without additional dinucleotide stacking interactions. The DNA-PAINT data showed that a single additional dinucleotide base stacking results in up to 250-fold stabilization for the DNA duplex nanostructure. We found that the dinucleotide base-stacking energies vary from -0.95 ± 0.12 kcal mol-1 to -3.22 ± 0.04 kcal mol-1 for C|T and A|C base-stackings, respectively. We demonstrate the application of base-stacking energetics in designing DNA-PAINT probes for multiplexed super-resolution imaging, and efficient assembly of higher-order DNA nanostructures. Our results will aid in designing functional DNA nanostructures, and DNA and RNA aptamers, and facilitate better predictions of the local DNA structure.

Single-Molecule Conductance of Intramolecular Hydrogen Bonding in Histamine on Gold

We perform single-molecule conductance measurements and DFT calculations on histamine, a biogenic amine that contains a flexible aliphatic linker and several nitrogen moieties with a potential for hydrogen bonding. Our study determines that junctions containing the free-base form of histamine can bridge through a molecular structure containing an intramolecular hydrogen bond. Conductance of this structure is higher than that through the saturated aliphatic linker. Flicker noise analysis of junction conductance confirms that transport occurs through the hydrogen bond and establishes a benchmark for noise measurements in hydrogen-bonded junctions. Overall, our work provides insights into the formation and conduction of intramolecular hydrogen bonding in single-molecule conductance measurements and into the conformations of the neurotransmitter histamine on noble metal surfaces.

π-Stacking among the Anthracenyl Groups of a Copper Complex Resulted in Doubling of Unit Cell Volume To Provide New Polymorphs

Two polymorphs of the 9-N-(3-imidazolylpropylamino)methylanthracene (Hanthraimmida) containing hydrated copper(II)-2,6-pyridinedicarboxylate complex are reported. The two polymorphs have either lamellar or Herringbone arrangements of π-stacks among the anthracenyl groups of organocation. The difference between the two polymorphs originated from having face-to-face stacking arrangements between the two anthracenyl groups of the symmetry independent cations within the unit cell in one of the polymorphs. The π-stacked anthracenyl groups in consecutive layers of the polymorphs are oriented in one direction in the polymorph designated as P1, whereas the polymorph designated as P2 has such orientations in opposite directions. The unit cell volume of the polymorph P2 (Z = 4) has approximately twice the volume of the polymorph P1 (Z = 2); it happend due to coalescence of two unit cells of P1 in the ab-crystallographic plane. A mixed methanol/water solvate of the copper complex is also reported. It has a channel-like arrangement of the cations; has the anions and the solvents within the cation embraced channel-like enclosures. This complex is unstable, once taken out from the methanol solvent, it transforms in real time to P2 by replacements of the methanol molecules by water molecules.

Influence of a Single Deuterium Substitution for Protium on the Frequency Generation of Different-Size Bubbles in IFNA17

The influence of a single 2H/1H replacement on the frequency generation of different-size bubbles in the human interferon alpha-17 gene (IFNA17) under various energies was studied by a developed algorithm and mathematical modeling without simplifications or averaging. This new approach showed the efficacy of researching DNA bubbles and open states both when all hydrogen bonds in nitrogenous base pairs are protium and after an 2H-substitution. After a single deuterium substitution under specific energies, it was demonstrated that the non-coding region of IFNA17 had a more significant regulatory role in bubble generation in the whole gene than the promoter had. It was revealed that a single deuterium substitution for protium has an influence on the frequency generation of DNA bubbles, which also depends on their size and is always higher for the smaller bubbles under the largest number of the studied energies. Wherein, compared to the natural condition under the same critical value of energy, the bigger raises of the bubble frequency occurrence (maximums) were found for 11-30 base pair (bp) bubbles (higher by 319%), 2-4 bp bubbles (higher by 300%), and 31 bp and over ones (higher by 220%); whereas the most significant reductions of the indicators (minimums) were observed for 11-30 bp bubbles (lower by 43%) and bubbles size over 30 bp (lower by 82%). In this study, we also analyzed the impact of several circumstances on the AT/GC ratio in the formation of DNA bubbles, both under natural conditions and after a single hydrogen isotope exchange. Moreover, based on the obtained data, substantial positive and inverse correlations were revealed between the AT/GC ratio and some factors (energy values, size of DNA bubbles). So, this modeling and variant of the modified algorithm, adapted for researching DNA bubbles, can be useful to study the regulation of replication and transcription in the genes under different isotopic substitutions in the nucleobases.

Leveraging Steric Moieties for Kinetic Control of DNA Strand Displacement Reactions

DNA strand displacement networks are a critical part of dynamic DNA nanotechnology and are proven primitives for implementing chemical reaction networks. Precise kinetic control of these networks is important for their use in a range of applications. Among the better understood and widely leveraged kinetic properties of these networks are toehold sequence, length, composition, and location. While steric hindrance has been recognized as an important factor in such systems, a clear understanding of its impact and role is lacking. Here, a systematic investigation of steric hindrance within a DNA toehold-mediated strand displacement network was performed through tracking kinetic reactions of reporter complexes with incremental concatenation of steric moieties near the toehold. Two subsets of steric moieties were tested with systematic variation of structures and reaction conditions to isolate sterics from electrostatics. Thermodynamic and coarse-grained computational modeling was performed to gain further insight into the impacts of steric hindrance. Steric factors yielded up to 3 orders of magnitude decrease in the reaction rate constant. This pronounced effect demonstrates that steric moieties can be a powerful tool for kinetic control in strand displacement networks while also being more broadly informative of DNA structural assembly in both DNA-based therapeutic and diagnostic applications that possess elements of steric hindrance through DNA functionalization with an assortment of chemistries.

Measuring the Enthalpy of an Individual Hydrogen Bond in a DNA Duplex with Nucleobase Isotope Editing and Variable-Temperature Infrared Spectroscopy

The level of interest in probing the strength of noncovalent interactions in DNA duplexes is high, as these weak forces dictate the range of suprastructures the double helix adopts under different conditions, in turn directly impacting the biological functions and industrial applications of duplexes that require making and breaking them to access the genetic code. However, few experimental tools can measure these weak forces embedded within large biological suprastructures in the native solution environment. Here, we develop experimental methods for detecting the presence of a single noncovalent interaction [a hydrogen bond (H-bond)] within a large DNA duplex in solution and measure its formation enthalpy (ΔH_f). We report that introduction of a H-bond into the TC2═O group from the noncanonical nucleobase 2-aminopurine produces an expected decrease ∼10 ± 0.76 cm^-1 (from ∼1720 cm^-1 in Watson-Crick to ∼1710 cm^-1 in 2-aminopurine), which correlates with an enthalpy of ∼0.93 ± 0.066 kcal/mol for this interaction.

Endogenous piRNAs Can Interact with the Omicron Variant of the SARS-CoV-2 Genome

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the COVID-19 pandemic, can still infect populations in many countries around the globe. The Omicron strain is the most mutated variant of SARS-CoV-2. The high transmissibility of the strain and its ability to evade immunity necessitate a priority study of its properties in order to quickly create effective means of preventing its spread. The current research aimed to examine the in silico interaction between PIWI-interacting RNAs (piRNAs) and the SARS-CoV-2 genome (gRNA) to identify endogenous piRNAs and propose synthetic piRNAs with strong antiviral activity for drug development. This study used validated bioinformatic approaches regarding the interaction of more than eight million piRNAs with the SARS-CoV-2 genome. The piRNAs’ binding sites (BSs) in the 5′UTR were located with overlapping nucleotide sequences termed clusters of BSs. Several BSs clusters have been found in the nsp3, nsp7, RNA-dependent RNA polymerase, endoRNAse, S surface glycoprotein, ORF7a, and nucleocapsid. Sixteen synthetic piRNAs that interact with gRNA have been proposed with free binding energy ranging from −170 kJ/mol to −175 kJ/mol, which can be used to create drugs that suppress the reproduction of SARS-CoV-2.

Novel genosensor for probing DNA mismatches and UV-induced DNA damage: Sequence-specific recognition

Human genome is continuously susceptible to changes that may lead to undesirable mutations causing various diseases and cancer. Vast majority of techniques has investigated the discrimination between base-pair mismatched nucleic acid, but many of these techniques are time-consuming, complex, expensive, and limited to the detection of specific type of dsDNA mismatches. In this study, we introduce a simple mix-and-read assay for the sensitive and cost-effective analysis of DNA base mismatches and UV-induced DNA damage using Hoechst genosensor dye (H258). This dye is a minor groove binder that undergoes a drastic conformational change upon binding with mismatch DNA. The difference in binding affinity between perfectly matched and mismatched DNA was studied for sequences at different base mismatch locations and finally, extended for the detection of dsDNA damage by UVC radiation in calf thymus DNA. In addition, a comparative DNA damage kinetic study was performed using H258 (minor groove binder) and EvaGreen (intercalating) dye to get insight on assay selectivity and sensitivity with dye binding mechanism. The result shows good reproducibility making H258 genosensor a cheaper alternative for DNA mismatch and damage studies with possibility of extension for in-vitro detection of hot spots of DNA mutations.

Inter- vs. Intra-Molecular Hydrogen Bond in Complexes of Nitrophthalic Acids with Pyridine

This study covers the analysis of isomeric forms of nitrophthalic acids with pyridine. This work dwells on the complementary experimental (X-ray, IR and Raman) and theoretical (Car-Parrinello Molecular Dynamics (CPMD) and Density Functional Theory (DFT)) studies of the obtained complexes. The conducted studies showed that steric repulsion between the nitro group in ortho-position and the carboxyl group causes significant isomeric changes. Modeling of the nitrophthalic acid—pyridine complex yielded a short strong intramolecular hydrogen bond (SSHB). The transition energy from the isomeric form with an intermolecular hydrogen bond to the isomeric form with an intramolecular hydrogen bond was estimated.

High-throughput single-molecule quantification of individual base stacking energies in nucleic acids

Base stacking interactions between adjacent bases in DNA and RNA are important for many biological processes and in biotechnology applications. Previous work has estimated stacking energies between pairs of bases, but contributions of individual bases has remained unknown. Here, we use a Centrifuge Force Microscope for high-throughput single molecule experiments to measure stacking energies between adjacent bases. We found stacking energies strongest between purines (G|A at -2.3 ± 0.2 kcal/mol) and weakest between pyrimidines (C|T at -0.5 ± 0.1 kcal/mol). Hybrid stacking with phosphorylated, methylated, and RNA nucleotides had no measurable effect, but a fluorophore modification reduced stacking energy. We experimentally show that base stacking can influence stability of a DNA nanostructure, modulate kinetics of enzymatic ligation, and assess accuracy of force fields in molecular dynamics simulations. Our results provide insights into fundamental DNA interactions that are critical in biology and can inform design in biotechnology applications.

Resolving the Chemical Formula of Nesquehonite via NMR Crystallography, DFT Computation, and Complementary Neutron Diffraction

Nesquehonite is a magnesium carbonate mineral relevant to carbon sequestration envisioned for carbon capture and storage of CO₂ . Its chemical formula remains controversial today, assigned as either a hydrated magnesium carbonate [MgCO₃  ⋅ 3H₂ O], or a hydroxy- hydrated- magnesium bicarbonate [Mg(HCO₃ )OH ⋅ 2H₂ O]. The resolution of this controversy is central to understanding this material's thermodynamic, phase, and chemical behavior. In an NMR crystallography study, using rotational-echo double-resonance ¹³ C{¹ H} (REDOR), ¹³ C-¹ H distances are determined with precision, and the combination of ¹³ C static NMR lineshapes and density functional theory (DFT) calculations are used to model different H atomic coordinates. [MgCO₃  ⋅ 3H₂ O] is found to be accurate, and evidence from neutron powder diffraction bolsters these assignments. Refined H positions can help understand how H-bonding stabilizes this structure against dehydration to MgCO₃ . More broadly, these results illustrate the power of NMR crystallography as a technique for resolving questions where X-ray diffraction is inconclusive.

π-Stacking Isomerism in Polycyclic Aromatic Hydrocarbons: The 2-Naphthalenethiol Dimer

π-Stacking is a common descriptor for face-to-face attractive forces between aromatic hydrocarbons. However, the physical origin of this interaction remains debatable. Here we examined π-stacking in a model homodimer formed by two thiol-substituted naphthalene rings. Two isomers of the 2-naphthalenethiol dimer were discovered using rotational spectroscopy, sharing a parallel-displaced crossed orientation and absence of thiol-thiol hydrogen bonds. One of the isomers presents C₂ symmetry, structurally analogous to the global minimum of the naphthalene dimer. The experimental data were rationalized with molecular orbital calculations, revealing a shallow potential energy surface. Noncovalent interactions are dominated by dispersion forces according to SAPT energy decomposition. In addition, the reduced electronic density shows a diffuse and extended region of inter-ring interactions, compatible with the description of π-stacking as a competition between dispersion and Pauli repulsion forces.

A two-residue nascent-strand steric gate controls synthesis of 2'-O-methyl- and 2'-O-(2-methoxyethyl)-RNA

Steric exclusion is a key element of enzyme substrate specificity, including in polymerases. Such substrate specificity restricts the enzymatic synthesis of 2'-modified nucleic acids, which are of interest in nucleic-acid-based drug development. Here we describe the discovery of a two-residue, nascent-strand, steric control 'gate' in an archaeal DNA polymerase. We show that engineering of the gate to reduce steric bulk in the context of a previously described RNA polymerase activity unlocks the synthesis of 2'-modified RNA oligomers, specifically the efficient synthesis of both defined and random-sequence 2'-O-methyl-RNA (2'OMe-RNA) and 2'-O-(2-methoxyethyl)-RNA (MOE-RNA) oligomers up to 750 nt. This enabled the discovery of RNA endonuclease catalysts entirely composed of 2'OMe-RNA (2'OMezymes) for the allele-specific cleavage of oncogenic KRAS (G12D) and β-catenin CTNNB1 (S33Y) mRNAs, and the elaboration of mixed 2'OMe-/MOE-RNA aptamers with high affinity for vascular endothelial growth factor. Our results open up these 2'-modified RNAs-used in several approved nucleic acid therapeutics-for enzymatic synthesis and a wider exploration in directed evolution and nanotechnology.

piRNA and miRNA Can Suppress the Expression of Multiple Sclerosis Candidate Genes

Multiple sclerosis (MS) is a common inflammatory demyelinating disease with a high mortality rate. MS is caused by many candidate genes whose specific involvement has yet to be established. The aim of our study was to identify endogenous miRNAs and piRNAs involved in the regulation of MS candidate gene expression using bioinformatic methods. A program was used to quantify the interaction of miRNA and piRNA nucleotides with mRNA of the target genes. We used 7310 miRNAs from three databases and 40,000 piRNAs. The mRNAs of the candidate genes revealed miRNA binding sites (BSs), which were located separately or formed clusters of BSs with overlapping nucleotide sequences. The miRNAs from the studied databases were generally bound to mRNAs in different combinations, but miRNAs from only one database were bound to the mRNAs of some genes. For the first time, a direct interaction between the complete sequence of piRNA nucleotides and the nucleotides of their mRNA BSs of target genes was shown. One to several clusters of BSs of miRNA and piRNA were identified in the mRNA of ADAM17, AHI1, CD226, EOMES, EVI5, IL12B, IL2RA, KIF21B, MGAT5, MLANA, SOX8, TNFRSF1A, and ZBTB46 MS candidate genes. These piRNAs form the expression regulation system of the MS candidate genes to coordinate the synthesis of their proteins. Based on these findings, associations of miRNAs, piRNAs, and candidate genes for MS diagnosis are recommended.

The role of filamentation in activation and DNA sequence specificity of the sequence-specific endonuclease SgrAI

Filament formation by metabolic, biosynthetic, and other enzymes has recently come into focus as a mechanism to fine-tune enzyme activity in the cell. Filamentation is key to the function of SgrAI, a sequence-specific DNA endonuclease that has served as a model system to provide some of the deepest insights into the biophysical characteristics of filamentation and its functional consequences. Structure-function analyses reveal that, in the filamentous state, SgrAI stabilizes an activated enzyme conformation that leads to accelerated DNA cleavage activity and expanded DNA sequence specificity. The latter is thought to be mediated by sequence-specific DNA structure, protein-DNA interactions, and a disorder-to-order transition in the protein, which collectively affect the relative stabilities of the inactive, non-filamentous conformation and the active, filamentous conformation of SgrAI bound to DNA. Full global kinetic modeling of the DNA cleavage pathway reveals a slow, rate-limiting, second-order association rate constant for filament assembly, and simulations of in vivo activity predict that filamentation is superior to non-filamenting mechanisms in ensuring rapid activation and sequestration of SgrAI's DNA cleavage activity on phage DNA and away from the host chromosome. In vivo studies demonstrate the critical requirement for accelerated DNA cleavage by SgrAI in its biological role to safeguard the bacterial host. Collectively, these data have advanced our understanding of how filamentation can regulate enzyme structure and function, while the experimental strategies used for SgrAI can be applied to other enzymatic systems to identify novel functional roles for filamentation.

Chemical Tools for the Study of DNA Repair

DNA repair enzymes continuously provide surveillance throughout our cells, protecting the enclosed DNA from the damage that is constantly arising from oxidation, alkylating species, and radiation. Members of this enzyme class are intimately linked to pathways controlling cancer and inflammation and are promising targets for diagnostics and future therapies. Their study is benefiting widely from the development of new tools and methods aimed at measuring their activities. Here, we provide an Account of our laboratory's work on developing chemical tools to study DNA repair processes in vitro, as well as in cells and tissues, and what we have learned by applying them.We first outline early work probing how DNA repair enzymes recognize specific forms of damage by use of chemical analogs of the damage with altered shapes and H-bonding abilities. One outcome of this was the development of an unnatural DNA base that is incorporated selectively by polymerase enzymes opposite sites of missing bases (abasic sites) in DNA, a very common form of damage.We then describe strategies for design of fluorescent probes targeted to base excision repair (BER) enzymes; these were built from small synthetic DNAs incorporating fluorescent moieties to engender light-up signals as the enzymatic reaction proceeds. Examples of targets for these DNA probes include UDG, SMUG1, Fpg, OGG1, MutYH, ALKBH2, ALKBH3, MTH1, and NTH1. Several such strategies were successful and were applied both in vitro and in cellular settings; moreover, some were used to discover small-molecule modulators of specific repair enzymes. One of these is the compound SU0268, a potent OGG1 inhibitor that is under investigation in animal models for inhibiting hyperinflammatory responses.To investigate cellular nucleotide sanitation pathways, we designed a series of "two-headed" nucleotides containing a damaged DNA nucleotide at one end and ATP at the other; these were applied to studying the three human sanitation enzymes MTH1, dUTPase, and dITPase, some of which are therapeutic targets. The MTH1 probe (ARGO) was used in collaboration with oncologists to measure the enzyme in tumors as a disease marker and also to develop the first small-molecule activators of the enzyme.We proceed to discuss the development of a "universal" probe of base excision repair processes (UBER), which reacts covalently with abasic site intermediates of base excision repair. UBER probes light up in real time as the reaction occurs, enabling the observation of base excision repair as it occurs in live cells and tissues. UBER probes can also be used in efficient and simple methods for fluorescent labeling of DNA. Finally, we suggest interesting directions for the future of this field in biomedicine and human health.

Nucleobase-modified nucleosides and nucleotides: Applications in biochemistry, synthetic biology, and drug discovery

. DNA is often referred to as the "molecule of life" since it contains the genetic blueprint for all forms of life on this planet. The core building blocks composing DNA are deoxynucleotides. While the deoxyribose sugar and phosphate group are ubiquitous, it is the composition and spatial arrangement of the four natural nucleobases, adenine (A), cytosine (C), guanine (G), and thymine (T), that provide diversity in the coding information present in DNA. The ability of DNA to function as the genetic blueprint has historically been attributed to the formation of proper hydrogen bonding interactions made between complementary nucleobases. However, recent chemical and biochemical studies using nucleobase-modified nucleotides that contain "non-hydrogen bonding" functional groups have challenged many of the dogmatic views for the necessity of hydrogen-bonding interactions for DNA stability and function. Based on years of exciting research, this area has expanded tremendously and is thus too expansive to provide a comprehensive review on the topic. As such, this review article provides an opinion highlighting how nucleobase-modified nucleotides are being applied in diverse biomedical fields, focusing on three exciting areas of research. The first section addresses how these analogs are used as mechanistic probes for DNA polymerase activity and fidelity during replication. This section outlines the synthetic logic and medicinal chemistry approaches used to replace hydrogen-bonding functional groups to examine the contributions of shape/size, nucleobase hydrophobicity, and pi-electron interactions. The second section extends these mechanistic studies to provide insight into how nucleobase-modified nucleosides are used in synthetic biology. One example is through expansion of the genetic code in which changing the composition of DNA makes it possible to site-specifically incorporate unnatural amino acids bearing unique functional groups into enzymes and receptors. The final section describes results of pre-clinical studies using nucleobase-modified nucleosides as potential therapeutic agents against diseases such as cancer.

Nucleotide Composition of Ultra-Conserved Elements Shows Excess of GpC and Depletion of GG and CC Dinucleotides

The public UCNEbase database, comprising 4273 human ultra-conserved noncoding elements (UCNEs), was thoroughly investigated with the aim to find any nucleotide signals or motifs that have made these DNA sequences practically unchanged over three hundred million years of evolution. Each UCNE comprises over 200 nucleotides and has at least 95% identity between humans and chickens. A total of 31,046 SNPs were found within the UCNE database. We demonstrated that every human has over 300 mutations within 4273 UCNEs. No association of UCNEs with non-coding RNAs, nor preference of a particular meiotic recombination rate within them were found. No sequence motifs associated with UCNEs nor their flanking regions have been found. However, we demonstrated that UCNEs have strong nucleotide and dinucleotide sequence abnormalities compared to genome averages. Specifically, UCNEs are depleted for CC and GG dinucleotides, while GC dinucleotides are in excess of 28%. Importantly, GC dinucleotides have extraordinarily strong stacking free-energy inside the DNA helix and unique resistance to dissociation. Based on the adjacent nucleotide stacking abnormalities within UCNEs, we conjecture that peculiarities in dinucleotide distribution within UCNEs may create unique 3D conformation and specificity to bind proteins. We also discuss the strange dynamics of multiple SNPs inside UCNEs and reasons why these sequences are extraordinarily conserved.

Formation Mechanism of Inter-Crosslink in DNA by Nitrogen Oxides Pollutants through A Diazonium Intermediate

Outdoor air pollution is a mixture of multiple atmospheric pollutants, among which nitrogen oxide (NOx) stands out due to its association with several diseases. NOx reactivity can conduct to DNA damage as severe as interstrand crosslinks (ICL) formation, that in turn is able to block DNA replication and transcription. Experimental studies have suggested that the ICL formation due to NOx is realized through a diazonium intermediate (DI). In this work, we have modeled the DI structure, including a DNA double-strand composed of two base pairs GC/CG, being diazotized as one of the guanine nucleotides. The structural stability of DNA with DI lesion was essayed through 500 ns molecular dynamics simulations. It was found that the DNA structure of the oligonucleotide is stable when the DI is present since the loss of a Guanine–Cytosine hydrogen bond is replaced by the presence of two cation-π interactions. Additionally, we have studied the mechanism of formation of a crosslink between the two guanine nucleobases from the modeled DI by carrying out DFT calculations at the M06-L/DNP+ level of theory. Our results show that the mechanism is thermodynamically favored by a strong stabilization of the ICL product, and the process is kinetically viable since its limiting stage is accessible.

In Silico Study of piRNA Interactions with the SARS-CoV-2 Genome

A prolonged pandemic with numerous human casualties requires a rapid search for means to control the various strains of SARS-CoV-2. Since only part of the human population is affected by coronaviruses, there are probably endogenous compounds preventing the spread of these viral pathogens. It has been shown that piRNA (PIWI-interacting RNAs) interact with the mRNA of human genes and can block protein synthesis at the stage of translation. Estimated the effects of piRNA on SARS-CoV-2 genomic RNA (gRNA) in silico. A cluster of 13 piRNA binding sites (BS) in the SARS-CoV-2 gRNA region encoding the oligopeptide was identified. The second cluster of BSs 39 piRNAs also encodes the oligopeptide. The third cluster of 24 piRNA BS encodes the oligopeptide. Twelve piRNAs were identified that strongly interact with the gRNA. Based on the identified functionally important endogenous piRNAs, synthetic piRNAs (spiRNAs) are proposed that will suppress the multiplication of the coronavirus even more strongly. These spiRNAs and selected endogenous piRNAs have little effect on human 17494 protein-coding genes, indicating a low probability of side effects. The piRNA and spiRNA selection methodology created for the control of SARS-CoV-2 (NC_045512.2) can be used to control all strains of SARS-CoV-2.

Assessing the Experimental Hydrogen Bonding Energy of the Cyclic Water Dimer Transition State with a Cyclooctatetraene-Based Molecular Balance

We have conducted an experimental and computational study of cyclooctatetraene-1,4/1,6-dimethanol (1,4 and 1,6) as a molecular balance with the goal in mind to determine the otherwise inaccessible hydrogen bonding energy (HBE) of the cyclic water dimer, which constitutes a transition state. The 1,4/1,6 folding equilibrium is governed by an intramolecular hydrogen bond in the folded 1,6-isomer, in which the OH groups adopt a cyclic planar geometry, akin to the structure of the cyclic water dimer transition state. We characterized hydrogen bonding in 1,6 and reference complexes utilizing SAPT2 + (3)δMP2/aug-cc-pVTZ and selected quantum theory of atoms in molecule descriptors at M06-2XD3(0)/ma-def2-TZVPP. Additionally, we computed HBEs at the DLPNO-CCSD(T)/aug-cc-pVQZ level of theory. We find that hydrogen bonding in 1,6 is very similar to the interaction in the C_i symmetric cyclic water dimer TS, both in magnitude and character. We experimentally determined the Gibbs free energy of the folding process (ΔG_eq) in a variety of organic solvents via nuclear magnetic resonance spectroscopy measurements at room temperature. By combining experimentally obtained ΔG_eq values with corrections derived from accurate computational methods, we provide estimates for the HBE of cyclic water dimers and the cyclic water dimer TS, as the most stable cyclic water dimer.

Introducing a New Bond-Forming Activity in an Archaeal DNA Polymerase by Structure-Guided Enzyme Redesign

DNA polymerases have evolved to feature a highly conserved activity across the tree of life: formation of, without exception, internucleotidyl O-P linkages. Can this linkage selectivity be overcome by design to produce xenonucleic acids? Here, we report that the structure-guided redesign of an archaeal DNA polymerase, 9°N, exhibits a new activity undetectable in the wild-type enzyme: catalyzing the formation of internucleotidyl N-P linkages using 3'-NH₂-ddNTPs. Replacing a metal-binding aspartate in the 9°N active site with asparagine was key to the emergence of this unnatural enzyme activity. MD simulations provided insights into how a single substitution enhances the productive positioning of a 3'-amino nucleophile in the active site. Further remodeling of the protein-nucleic acid interface in the finger subdomain yielded a quadruple-mutant variant (9°N-NRQS) displaying DNA-dependent NP-DNA polymerase activity. In addition, the engineered promiscuity of 9°N-NRQS was leveraged for one-pot synthesis of DNA─NP-DNA copolymers. This work sheds light on the molecular basis of substrate fidelity and latent promiscuity in enzymes.

Nearest-Neighbor dsDNA Stability Analysis Using Alchemical Free-Energy Simulations

The thermodynamic stability of double-stranded (ds)DNA depends on its sequence. It is influenced by the base pairing and stacking with neighboring bases along DNA molecules. Semiempirical schemes are available that allow us to predict the thermodynamic stability of DNA sequences based on empirically derived nearest-neighbor contributions of base pairs formed in the context of all possible nearest-neighbor base pairs. Current molecular dynamics (MD) simulations allow one to simulate the dynamics of DNA molecules in good agreement with experimentally obtained structures and available data on conformational flexibility. However, the suitability of current force field methods to reproduce dsDNA stability and its sequence dependence has been much less well tested. We have employed alchemical free-energy simulations of whole base pair transversions in dsDNA and in unbound single-stranded partner molecules. Such transversions change the sequence context but not the nucleotide content or base pairing in dsDNA and allow a direct comparison with the empirical nearest-neighbor dsDNA stability model. For the alchemical free-energy changes in the unbound single-stranded (ss)DNA partner molecules, we tested different setups assuming either complete unstacking or unrestrained simulations with partial stacking in the unbound ssDNA. The free-energy simulations predicted nearest-neighbor effects of similar magnitude, as observed experimentally but showed overall limited correlation with experimental data. An inaccurate description of stacking interactions and other possible reasons such as the neglect of electronic polarization effects are discussed. The results indicate the need to improve the realistic description of stacking interactions in current molecular mechanic force fields.

Aromatic interactions directing peptide nano-assembly

Self-assembly is a process of spontaneous organization of molecules as a result of non-covalent interactions. Organized self-assembly at the nano level is emerging as a powerful tool in the bottom-up fabrication of functional nanostructures for targeted applications. Aromatic π-π stacking plays a significant role by facilitating the persistent supramolecular association of individual subunits to the self-assembled structures of high stability. Understanding, the supramolecular chemistry of the materials interacting through aromatic interactions, is of tremendous interest in not only constructing functional materials but also in revealing the mechanism of molecular assembly in living organisms. This chapter aims to focus on understanding the potential role of π-π interactions in directing and regulating the self-assembly of peptide nanostructures. The scope of the chapter starts with an outline of the history and mechanism of the aromatic π-π interactions. It progresses through the design strategy for the assembly of peptides containing aromatic rings, the conditions affecting the aromatic stacking interactions, their resulting nanoassemblies, properties, and applications. The properties and applications of the supramolecular materials formed through the aromatic stacking interactions are highlighted to provide an increased understanding of the role of weak interactions in the design and construction of novel functional materials.

Pre-Transition State and Apo Structures of the Filament-Forming Enzyme SgrAI Elucidate Mechanisms of Activation and Substrate Specificity

Enzyme filamentation is a widespread phenomenon that mediates enzyme regulation and function. For the filament-forming sequence-specific DNA endonuclease SgrAI, the process of filamentation both accelerates its DNA cleavage activity and expands its DNA sequence specificity, thus allowing for many additional DNA sequences to be rapidly cleaved. Both outcomes—the acceleration of DNA cleavage and the expansion of sequence specificity—are proposed to regulate critical processes in bacterial innate immunity. However, the mechanistic bases underlying these events remain unclear. Herein, we describe two new structures of the SgrAI enzyme that shed light on its catalytic function. First, we present the cryo-EM structure of filamentous SgrAI bound to intact primary site DNA and Ca²⁺ resolved to ∼2.5 Å within the catalytic center, which represents the trapped enzyme–DNA complex prior to the DNA cleavage reaction. This structure reveals important conformational changes that contribute to the catalytic mechanism and the binding of a second divalent cation in the enzyme active site, which is expected to contribute to increased DNA cleavage activity of SgrAI in the filamentous state. Second, we present an X-ray crystal structure of DNA-free (apo) SgrAI resolved to 2.0 Å resolution, which reveals a disordered loop involved in DNA recognition. Collectively, these multiple new observations clarify the mechanism of expansion of DNA sequence specificity of SgrAI, including the indirect readout of sequence-dependent DNA structure, changes in protein–DNA interactions, and the disorder-to-order transition of a crucial DNA recognition element.

Self-Assembly of All-DNA Rods with Controlled Patchiness

Double-stranded DNA (dsDNA) fragments exhibit noncovalent attractive interactions between their tips. It is still unclear how DNA liquid crystal self-assembly is affected by such blunt-end attractions. It is demonstrated that stiff dsDNA fragments with moderate aspect ratio can specifically self-assemble in concentrated aqueous solutions into different types of smectic mesophases on the basis of selectively screening of blunt-end DNA stacking interactions. To this end, this type of attractions are engineered at the molecular level by constructing DNA duplexes where the attractions between one or both ends are screened by short hairpin caps. All-DNA bilayer and monolayer smectic-A type of phases, as well as a columnar phase, can be stabilized by controlling attractions strength. The results imply that the so far elusive smectic-A in DNA rod-like liquid crystals is a thermodynamically stable phase. The existence of the bilayer smectic phase is confirmed by Monte-Carlo simulations of hard cylinders decorated with one attractive terminal site. This work demonstrates that DNA blunt-ends behave as well-defined monovalent attractive patches whose strength and position can be potentially precisely tuned, highlighting unique opportunities concerning the stabilization of nonconventional DNA-based lyotropic liquid crystal phases assembled by all-DNA patchy particles with arbitrary geometry and composition.

Identification of Bovine miRNAs with the Potential to Affect Human Gene Expression

Milk and other products from large mammals have emerged during human evolution as an important source of nutrition. Recently, it has been recognized that exogenous miRNAs (mRNA inhibited RNA) contained in milk and other tissues of the mammalian body can enter the human body, which in turn have the ability to potentially regulate human metabolism by affecting gene expression. We studied for exogenous miRNAs from Bos taurus that are potentially contain miRNAs from milk and that could act postprandially as regulators of human gene expression. The interaction of 17,508 human genes with 1025 bta-miRNAs, including 245 raw milk miRNAs was studied. The milk bta-miR-151-5p, bta-miR-151-3p, bta-miRNA-320 each have 11 BSs (binding sites), and bta-miRNA-345-5p, bta-miRNA-614, bta-miRNA-1296b and bta-miRNA-149 has 12, 14, 15 and 26 BSs, respectively. The bta-miR-574-5p from cow’s milk had 209 human genes in mRNAs from one to 25 repeating BSs. We found 15 bta-miRNAs that have 100% complementarity to the mRNA of 13 human target genes. Another 12 miRNAs have BSs in the mRNA of 19 human genes with 98% complementarity. The bta-miR-11975, bta-miR-11976, and bta-miR-2885 BSs are located with the overlap of nucleotide sequences in the mRNA of human genes. Nucleotide sequences of BSs of these miRNAs in 5′UTR mRNA of human genes consisted of GCC repeats with a total length of 18 nucleotides (nt) in 18 genes, 21 nt in 11 genes, 24 nt in 14 genes, and 27–48 nt in nine genes. Nucleotide sequences of BSs of bta-miR-11975, bta-miR-11976, and bta-miR-2885 in CDS mRNA of human genes consisted of GCC repeats with a total length of 18 nt in 33 genes, 21 nt in 13 genes, 24 nt in nine genes, and 27–36 nt in 11 genes. These BSs encoded polyA or polyP peptides. In only one case, the polyR (SLC24A3 gene) was encoded. The possibility of regulating the expression of human genes by exogenous bovine miRNAs is discussed.

Dielectricity of a molecularly crowded solution accelerates NTP misincorporation during RNA-dependent RNA polymerization by T7 RNA polymerase

In biological systems, the synthesis of nucleic acids, such as DNA and RNA, is catalyzed by enzymes in various aqueous solutions. However, substrate specificity is derived from the chemical properties of the residues, which implies that perturbations of the solution environment may cause changes in the fidelity of the reaction. Here, we investigated non-promoter-based synthesis of RNA using T7 RNA polymerase (T7 RNAP) directed by an RNA template in the presence of polyethylene glycol (PEG) of various molecular weights, which can affect polymerization fidelity by altering the solution properties. We found that the mismatch extensions of RNA propagated downstream polymerization. Furthermore, PEG promoted the polymerization of non-complementary ribonucleoside triphosphates, mainly due to the decrease in the dielectric constant of the solution. These results indicate that the mismatch extension of RNA-dependent RNA polymerization by T7 RNAP is driven by the stacking interaction of bases of the primer end and the incorporated nucleotide triphosphates (NTP) rather than base pairing between them. Thus, proteinaceous RNA polymerase may display different substrate specificity with changes in dielectricity caused by molecular crowding conditions, which can result in increased genetic diversity without proteinaceous modification.

A new glimpse on the active site of SARS-CoV-2 3CLpro, coupled with drug repurposing study

Coronavirus disease 2019 (COVID-19) is caused by novel severe acute respiratory syndrome coronavirus (SARS-CoV-2). Its main protease, 3C-like protease (3CLpro), is an attractive target for drug design, due to its importance in virus replication. The analysis of the radial distribution function of 159 3CLpro structures reveals a high similarity index. A study of the catalytic pocket of 3CLpro with bound inhibitors reveals that the influence of the inhibitors is local, perturbing dominantly only residues in the active pocket. A machine learning based model with high predictive ability against SARS-CoV-2 3CLpro is designed and validated. The model is used to perform a drug-repurposing study, with the main aim to identify existing drugs with the highest 3CLpro inhibition power. Among antiviral agents, lopinavir, idoxuridine, paritaprevir, and favipiravir showed the highest inhibition potential. Enzyme - ligand interactions as a key ingredient for successful drug design.

DNA origami nano-mechanics

Invention of DNA origami has transformed the fabrication and application of biological nanomaterials. In this review, we discuss DNA origami nanoassemblies according to their four fundamental mechanical properties in response to external forces: elasticity, pliability, plasticity and stability. While elasticity and pliability refer to reversible changes in structures and associated properties, plasticity shows irreversible variation in topologies. The irreversible property is also inherent in the disintegration of DNA nanoassemblies, which is manifested by its mechanical stability. Disparate DNA origami devices in the past decade have exploited the mechanical regimes of pliability, elasticity, and plasticity, among which plasticity has shown its dominating potential in biomechanical and physiochemical applications. On the other hand, the mechanical stability of the DNA origami has been used to understand the mechanics of the assembly and disassembly of DNA nano-devices. At the end of this review, we discuss the challenges and future development of DNA origami nanoassemblies, again, from these fundamental mechanical perspectives.

Stacked but not Stuck: Unveiling the Role of π→π* Interactions with the Help of the Benzofuran-Formaldehyde Complex

The 1:1 benzofuran–formaldehyde complex has been chosen as model system for analyzing π→π* interactions in supramolecular organizations involving heteroaromatic rings and carbonyl groups. A joint “rotational spectroscopy–quantum chemistry” strategy unveiled the dominant role of π→π* interactions in tuning the intermolecular interactions of such adduct. The exploration of the intermolecular potential energy surface led to the identification of 14 low‐energy minima, with 4 stacked isomers being more stable than those linked by hydrogen bond or lone‐pair→π interactions. All energy minima are separated by loose transition states, thus suggesting an effective relaxation to the global minimum under the experimental conditions. This expectation has been confirmed by the experimental detection of only one species, which was unambiguously assigned owing to the computation of accurate spectroscopic parameters and the characterization of 11 isotopologues. The large number of isotopic species opened the way to the determination of the first semi‐experimental equilibrium structure for a molecular complex of such a dimension.

Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

Herein, we provide new insights into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of a large portion of human B-DNA by using advanced quantum mechanical methods. Our results indicate that (i) the effect of non-neighboring bases on the inter-strand interaction is negligibly small, (ii) London dispersion effects are essential for the stability of the duplex structure, (iii) the largest contribution to the stability of the duplex structure is the Watson-Crick base pairing - consistent with previous computational investigations, (iv) the effect of stacking between adjacent bases is relatively small but still essential for the duplex structure stability and (v) there are no cooperativity effects between intra-strand stacking and inter-strand base pairing interactions. These results are consistent with atomic force microscope measurements and provide the first theoretical validation of nearest neighbor approaches for predicting thermodynamic data of arbitrary DNA sequences.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Bifacial PNAs Destabilize MALAT1 by 3' A-Tail Displacement from the U-Rich Internal Loop

We report herein a new class of synthetic reagents for targeting the element for nuclear expression (ENE) in MALAT1, a long noncoding RNA upregulated in many cancers. The cis-acting ENE contains a U-rich internal loop (URIL) that forms an 11 base UAU-rich triplex stem with the truncated 3' oligo-A tail of MALAT1, protecting the terminus from exonuclease digestion and greatly extending transcript lifetime. Bifacial peptide nucleic acids (bPNAs) similarly bind URILs via base triple formation between two uracil bases and a synthetic base, melamine. We synthesized a set of low molecular weight bPNAs composed of α-linked peptide, isodipeptide, and diketopiperazine backbones and evaluated their ENE binding efficacy in vitro via oligo-A strand displacement and consequent exonuclease sensitivity. Degradation was greatly enhanced by bPNA treatment in the presence of exonucleases, with ENE half-life plunging to 6 min from >24 h. RNA digestion kinetics could clearly distinguish between bPNAs with similar URIL affinities, highlighting the utility of functional assays for evaluating synthetic RNA binders. In vitro activity was mirrored by a 50% knockdown of MALAT1 expression in pancreatic cancer (PANC-1) cells upon treatment with bPNAs, consistent with intracellular digestion triggered by a similar ENE A-tail displacement mechanism. Pulldown from PANC-1 total RNA with biotinylated bPNA enriched MALAT1 > 4000× , supportive of bPNA-URIL selectivity. Together, these experiments establish the feasibility of native transcript targeting by bPNA in both in vitro and intracellular contexts. Reagents such as bPNAs may be useful tools for the investigation of transcripts stabilized by cis-acting poly(A) binding RNA elements.

SAMPDI-3D: predicting the effects of protein and DNA mutations on protein-DNA interactions

MOTIVATION

Mutations that alter protein-DNA interactions may be pathogenic and cause diseases. Therefore, it is extremely important to quantify the effect of mutations on protein-DNA binding free energy to reveal the molecular origin of diseases and to assist the development of treatments. Although several methods that predict the change of protein-DNA binding affinity upon mutations in the binding protein were developed, the effect of DNA mutations was not considered yet.

RESULTS

Here, we report a new version of SAMPDI, the SAMPDI-3D, which is a gradient boosting decision tree machine learning method to predict the change of the protein-DNA binding free energy caused by mutations in both the binding protein and the bases of the corresponding DNA. The method is shown to achieve Pearson correlation coefficient of 0.76 and 0.80 in a benchmarking test against experimentally determined change of the binding free energy caused by mutations in the binding protein or DNA, respectively. Furthermore, three datasets collected from literature were used to do blind benchmark for SAMPDI-3D and it is shown that it outperforms all existing state-of-the-art methods. The method is very fast allowing for genome-scale investigations.

AVAILABILITY

It is available as a web server and a stand-code at http://compbio.clemson.edu/SAMPDI-3D/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pyrophosphate release acts as a kinetic checkpoint during high-fidelity DNA replication by the Staphylococcus aureus replicative polymerase PolC

Bacterial replication is a fast and accurate process, with the bulk of genome duplication being catalyzed by the α subunit of DNA polymerase III within the bacterial replisome. Structural and biochemical studies have elucidated the overall properties of these polymerases, including how they interact with other components of the replisome, but have only begun to define the enzymatic mechanism of nucleotide incorporation. Using transient-state methods, we have determined the kinetic mechanism of accurate replication by PolC, the replicative polymerase from the Gram-positive pathogen Staphylococcus aureus. Remarkably, PolC can recognize the presence of the next correct nucleotide prior to completing the addition of the current nucleotide. By modulating the rate of pyrophosphate byproduct release, PolC can tune the speed of DNA synthesis in response to the concentration of the next incoming nucleotide. The kinetic mechanism described here would allow PolC to perform high fidelity replication in response to diverse cellular environments.

Mechanical characterization of base analogue modified nucleic acids by force spectroscopy

We use mechanical unfolding of single DNA hairpins with modified bases to accurately assess intra- and intermolecular forces in nucleic acids. As expected, the modification stabilizes the hybridized hairpin, but we also observe intriguing stacking interactions in the unfolded hairpin. Our study highlights the benefit of using base-modified nucleic acids in force-spectroscopy.

In acid-aminopyrimidine continuum: experimental and computational studies of furan tetracarboxylate-2-aminopyrimidinium salt

Salts and cocrystals are the two important solid forms when a carboxylic acid crystallizes with an aminopyrimidine base such that the extent of proton transfer distinguishes between them. The ΔpK_a value (pK_a(base) − pK_a(acid)) predicts whether the proton transfer will occur or not. However, the ΔpK_a range, 0 < ΔpK_a < 3, is elusive where the formation of cocrystal or salt cannot be predicted. The current study has been done to obtain a generalization in this elusive range with the Cambridge Structural Database (CSD). Based on the generalization, a novel salt (FTCA)⁻(2-AP)⁺ of furantetracarboxylic acid (FTCA) with 2-aminopyrimidine (2-AP) is obtained. The structural confirmation was done by single-crystal X-ray diffraction (SCXRD). Density functional theory (DFT) calculations were performed at the IEF-PCM-B3LYP-D3/6-311G(d,p) level to optimize the geometrical coordinates of salt for frontier molecular orbitals (FMOs) and molecular electrostatic potential (MESP). The geometrical parameters of most of the atoms of the optimized salt structure were comparable with SCXRD data. Additionally, results of other computational methods such as ab initio (Hartree–Fock; HF and second-order-Møller–Plesset perturbation; MP2) and semi-empirical were also compared with experimental results of the salt. Quantum theory of atoms in molecules (QTAIM), reduced density gradient (RDG), and natural bond orbital (NBO) analyses were done to calculate the strength and nature of non-covalent interactions present in the salt. Furthermore, Hirshfeld surface analysis, interaction energy calculations, and total energy frameworks were performed for qualitative and quantitative estimations of strong and weak intermolecular interactions.

Generalization in the elusive ΔpK_a range, experimental and computational studies of furan tetracarboxylate-2-aminopyrimidinium salt.

17O NMR Spectroscopy: A Novel Probe for Characterizing Protein Structure and Folding

Oxygen is a key atom that maintains biomolecular structures, regulates various physiological processes, and mediates various biomolecular interactions. Oxygen-17 (¹⁷O), therefore, has been proposed as a useful probe that can provide detailed information about various physicochemical features of proteins. This is attributed to the facts that (1) ¹⁷O is an active isotope for nuclear magnetic resonance (NMR) spectroscopic approaches; (2) NMR spectroscopy is one of the most suitable tools for characterizing the structural and dynamical features of biomolecules under native-like conditions; and (3) oxygen atoms are frequently involved in essential hydrogen bonds for the structural and functional integrity of proteins or related biomolecules. Although ¹⁷O NMR spectroscopic investigations of biomolecules have been considerably hampered due to low natural abundance and the quadruple characteristics of the ¹⁷O nucleus, recent theoretical and technical developments have revolutionized this methodology to be optimally poised as a unique and widely applicable tool for determining protein structure and dynamics. In this review, we recapitulate recent developments in ¹⁷O NMR spectroscopy to characterize protein structure and folding. In addition, we discuss the highly promising advantages of this methodology over other techniques and explain why further technical and experimental advancements are highly desired.

Dipole moment enhanced π-π stacking in fluorophenylacetylenes is carried over from gas-phase dimers to crystal structures propagated through liquid like clusters

The aggregates of monofluorinated phenylacetylenes in the gas-phase, investigated using the IR-UV double resonance spectroscopic method in combination with extensive structural search and electronic structure calculations, reveal the formation of liquid-like clusters with a π-stacked dimeric core. The structural assignment based on the IR spectra in the acetylenic and aromatic C-H stretching regions suggests that, unlike the parent non-fluorinated phenylacetylene, the substitution of a F atom on the phenyl ring increases the dipole moment, leading to robustness in the formation of a ππ stacked dimer, which propagates incorporating C-Hπ_{Ar/Ac} and C-HF interactions involving both acetylenic and aromatic C-H groups. The structural evolution of fluorophenylacetylene aggregates in the gas phase shows marginal effects due to fluorine atom position on the phenyl ring, with substitution in the para-position tending towards phenylacetylene. The present study signifies that the ππ stacked dimers act as a nucleus for the growth of higher clusters to which other molecular units are added predominantly via the {Ar}_C-Hπ_{Ar} type of interaction and the dominant interactions present in the crystal structures gradually emerge with increasing cluster size. Based on these features, gas-phase clusters of fluorophenylacetylene are hypothesized as "liquid-like clusters" acting as intermediates in the generation of various polymorphic forms starting from a ππ stacked dimer as the core molecular unit.

DNA nanotechnology provides an avenue for the construction of programmable dynamic molecular systems

Self-assembled supramolecular structures in living cells and their dynamics underlie various cellular events, such as endocytosis, cell migration, intracellular transport, cell metabolism, and gene expression. Spatiotemporally regulated association/dissociation and generation/degradation of assembly components is one of the remarkable features of biological systems. The significant advancement in DNA nanotechnology over the last few decades has enabled the construction of various-shaped nanostructures via programmed self-assembly of sequence-designed oligonucleotides. These nanostructures can further be assembled into micrometer-sized structures, including ordered lattices, tubular structures, macromolecular droplets, and hydrogels. In addition to being a structural material, DNA is adopted to construct artificial molecular circuits capable of activating/inactivating or producing/decomposing target DNA molecules based on strand displacement or enzymatic reactions. In this review, we provide an overview of recent studies on artificially designed DNA-based self-assembled systems that exhibit dynamic features, such as association/dis-sociation of components, phase separation, stimulus responsivity, and DNA circuit-regulated structural formation. These biomacromolecule-based, bottom-up approaches for the construction of artificial molecular systems will not only throw light on bio-inspired nano/micro engineering, but also enable us to gain insights into how autonomy and adaptability of living systems can be realized.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Enzymatic construction of metal-mediated nucleic acid base pairs

Artificial metal base pairs have become increasingly important in nucleic acids chemistry due to their high thermal stability, water solubility, orthogonality to natural base pairs, and low cost of production. These interesting properties combined with ease of chemical and enzymatic synthesis have prompted their use in several practical applications, including the construction of nanomolecular devices, ions sensors, and metal nanowires. Chemical synthesis of metal base pairs is highly efficient and enables the rapid screening of novel metal base pair candidates. However, chemical synthesis is limited to rather short oligonucleotides and requires rather important synthetic efforts. Herein, we discuss recent progress made for the enzymatic construction of metal base pairs that can alleviate some of these limitations. First, we highlight the possibility of generating metal base pairs using canonical nucleotides and then describe how modified nucleotides can be used in this context. We also provide a description of the main analytical techniques used for the analysis of the nature and the formation of metal base pairs together with relevant examples of their applications.

Evolutionary Changes in the Interaction of miRNA With mRNA of Candidate Genes for Parkinson's Disease

Parkinson’s disease (PD) exhibits the second-highest rate of mortality among neurodegenerative diseases. PD is difficult to diagnose and treat due to its polygenic nature. In recent years, numerous studies have established a correlation between this disease and miRNA expression; however, it remains necessary to determine the quantitative characteristics of the interactions between miRNAs and their target genes. In this study, using novel bioinformatics approaches, the quantitative characteristics of the interactions between miRNAs and the mRNAs of candidate PD genes were established. Of the 6,756 miRNAs studied, more than one hundred efficiently bound to mRNA of 61 candidate PD genes. The miRNA binding sites (BS) were located in the 5′-untranslated region (5′UTR), coding sequence (CDS) and 3′-untranslated region (3′UTR) of the mRNAs. In the mRNAs of many genes, the locations of miRNA BS with overlapping nucleotide sequences (clusters) were identified. Such clusters substantially reduced the proportion of nucleotide sequences of miRNA BS in the 5′UTRs, CDSs, and 3′UTRs. The organization of miRNA BS into clusters leads to competition among miRNAs to bind mRNAs. Differences in the binding characteristics of miRNAs to the mRNAs of genes expressed at different rates were identified. Single miRNA BS, polysites for the binding for one miRNA, and multiple BS for two or more miRNAs in one mRNA were identified. Evolutionary changes in the BS of miRNAs and their clusters in 5′UTRs, CDSs and 3′UTRs of mRNA of orthologous candidate PD genes were established. Based on the quantitative characteristics of the interactions between miRNAs and mRNAs candidate PD genes, several associations recommended as markers for the diagnosis of PD.

Mismatch detection in homologous strand exchange amplified by hydrophobic effects

In contrast to DNA replication and transcription where nucleotides are added and matched one by one, homologous recombination by DNA strand exchange tests whole sequences for complementarity, which requires elimination of mismatched yet thermodynamically stable intermediates. To understand the remarkable sequence specificity of homologous recombination, we have studied strand exchange between a 20-mer duplex containing one single mismatch (placed at varied positions) with the matching single strand in presence of poly(ethylene glycol) representing a semi-hydrophobic environment. A FRET-based assay shows that rates and yields of strand exchange from mismatched to matched strands rapidly increase with semi-hydrophobic co-solute concentration, contrasting previously observed general strand exchange accelerating effect of ethyl glycol ethers. We argue that this effect is not caused simply by DNA melting or solvent-induced changes of DNA conformation but is more complex involving several mechanisms. The catalytic effects, we propose, involve strand invasion facilitated by reduced duplex stability due to weakened base stacking ("longitudinal breathing"). Secondly, decreased water activity makes base-pair hydrogen bonds stronger, increasing the relative energy penalty per mismatch. Finally, unstacked mismatched bases (gaps) are stabilized through partly intercalated hydrophobic co-solvent molecules, assisting nucleation of strand invasion at the point of mismatch. We speculate that nature long ago discovered, and now exploits in various enzymes, that sequence recognition power of nucleic acids may be modulated in a hydrophobic environment.