Evaluation of the Gibbs Free Energy Changes and Melting Temperatures of DNA/DNA Duplexes Using Hybridization Enthalpy Calculated by Molecular Dynamics Simulation

PD-L1 DNA aptamers isolated from agarose-bead SELEX

Increased expression and activity of the PD-L1/PD-1 pathway suppresses the activation of cytotoxic T cells, which is vital in anti-tumour defence, allowing tumours to rise, expand and progress. Current strategies using antibodies to target PD-1/PD-L1 have been very effective in cancer therapeutics and companion diagnostics. Aptamers are a new class of molecules that offer an alternative to antibodies. Herein, the systematic evolution of ligands by exponential enrichment (SELEX) using agarose slurry beads was conducted to isolate DNA aptamers specific to recombinant human PD-L1 (rhPD-L1). Isolated aptamers were sequenced and analysed using MEGA X and structural features were examined using mFold. Three aptamer candidates (P33, P32, and P12) were selected for evaluation of binding affinity (dissociation constant, Kd) using ELONA and specificity and competitive inhibition assessment using the potentiostat-electrochemical method. Among those three, P32 displayed the highest specificity (8 nM) against PD-L1. However, P32 competes for the same binding site with the control antibody, 28-8. This study warrants further assessment of P32 aptamer as a potential, cost-effective alternative tool for targeting PD-L1.

Construction of a T m-value prediction model and molecular dynamics study of AmNA-containing gapmer antisense oligonucleotide

RNase H-dependent antisense oligonucleotides (gapmer ASOs) represent a class of nucleic acid therapeutics that bind to target RNA to facilitate RNase H-mediated RNA cleavage, thereby regulating the expression of disease-associated proteins. Integrating artificial nucleic acids into gapmer ASOs enhances their therapeutic efficacy. Among these, amido-bridged nucleic acid (AmNA) stands out for its potential to confer high affinity and stability to ASOs. However, a significant challenge in the design of gapmer ASOs incorporating artificial nucleic acids, such as AmNA, is the accurate prediction of their melting temperature (T m ) values. The T m is a critical parameter for designing effective gapmer ASOs to ensure proper functioning. However, predicting accurate T m values for oligonucleotides containing artificial nucleic acids remains problematic. We developed a T m prediction model using a library of AmNA-containing ASOs to address this issue. We measured the T m values of 157 oligonucleotides through differential scanning calorimetry, enabling the construction of an accurate prediction model. Additionally, molecular dynamics simulations were used to elucidate the molecular mechanisms by which AmNA modifications elevate T m , thereby informing the design strategies of gapmer ASOs.

Repurposing of Drug Bank Compounds against Plasmodium falciparum Dihydroorotate Dehydrogenase as novel anti malarial drug candidates by Computational approaches

This study aimed to repurpose Drug Bank Compounds against P. falciparum Dihydroorotate dehydrogenase (Pf-DHODH)a potential molecular target for antimalarial drug development due to its vital role in P. falciparum survival. Initially, the MATGEN server was used to screen drugs against Pf-DHODH (PDB ID 6GJG), followed by revalidating the results through docking by Autodock Vina through PyRx. Based on the docking results, three drugs namely, Talnifumate, Sulfaphenazole, and (3S)-N-[(2S)-1-[2-(1H-indol-3-yl)ethylamino]-1-oxopropan-2-yl]-1-(4-methoxyphenyl)-5-oxopyrrolidine-3-carboxamide-were subjected to molecular dynamics simulation for 100 ns. Molecular dynamics simulation results indicate that (3S)-N-[(2S)-1-[2-(1H-indol-3-yl)ethylamino]-1-oxopropan-2-yl]-1-(4-methoxyphenyl)-5-oxopyrrolidine-3-carboxamide- and Sulfaphenazole may target Pf-DHODH by forming a stable protein-ligand complex as they showed better free binding energy -130.58 kJ/mol, and -79.84 kJ/mol, respectively as compared to the free binding energy 116.255 kJ/mol of the reference compound; 3,6-dimethyl- ~ {N}-[4-(trifluoromethyl)phenyl]-[1,2]oxazolo[5,4-d]pyrimidin-4-amine. Although the studied compounds are drugs, still we applied Lipinski's rules and ADMET analysis that reconfirmed that these drugs have favorable drug-like properties. In conclusion, the results of the study show that Talniflumate and Sulfaphenazole may be potential antimalarial drug candidates.The derivatives of these drugs could be designed and tested to develop better drugs against Plasmodium species.

Graphical Abstract

Evaluation of miR-155 silencing using a molecular beacon in human lung adenocarcinoma cell line

Lung cancer (LC) is a leading cause of global cancer-related deaths, highlighting the development of innovative methods for biomarker detection improving the early diagnostics. microRNAs (miRs) alterations are known to be involved in the initiation and progression of human cancers and can act as biomarkers for diagnostics and treatment. Herein, we develop the application of molecular beacon (MB) technology to monitor miR-155-3p expression in human lung adenocarcinoma A549 cells without complementary DNA synthesis, amplification, or expensive reagents. Furthermore, we produced gold nanoparticles (AuNPs) for delivering antisense oligonucleotides into A549 cells to reduce miR-155-3p expression, which was subsequently detectable using the MB. The MB was designed and structural characterized by Förster Resonance Energy Transfer (FRET)-melting, Circular Dichroism (CD), Nuclear magnetic resonance (NMR), and fluorometric experiments, and then the hybridization conditions were optimized for an in vitro approach involving the detection of miR-155-3p in total RNA extracted from A549 cell line. The expression profile of miR-155-3p was obtained by RT-qPCR. The results demonstrated that MB was properly designed and showed efficacy in targeting miR-155-3p. Furthermore, a limit of detection down to nanomolar concentration was achieved and the specificity of the biosensor was proved. Moreover, the self-assembly of ASOs with AuNPs exhibited exceptional target specificity, effectively silencing miR-155-3p. Notably, compared to lipid-based transfection agent, AuNPs displayed superior silencing efficiency. We highlighted the ability of MB to detect changes in the target gene expression after gene silencing. Overall, this innovative approach represents a promising tool for detecting various biomarkers at the same time, with potential applications in clinical settings.

Diversity of Self-Assembled RNA Complexes: From Nanoarchitecture to Nanomachines

New tool development for various nucleic acid applications is an essential task in RNA nanotechnology. Here, we determined the ability of self-limited complex formation by a pair of oligoribonucleotides carrying two pairwise complementary blocks connected by a linker of different lengths in each chain. The complexes were analyzed using UV melting, gel shift assay analysis, and molecular dynamics (MD) simulations. We have demonstrated the spontaneous formation of various self-limited and concatemer complexes. The linear concatemer complex is formed by a pair of oligomers without a linker in at least one of them. Longer linkers resulted in the formation of circular complexes. The self-limited complexes formation was confirmed using the toehold strand displacement. The MD simulations indicate the reliability of the complexes' structure and demonstrate their dynamics, which increase with the rise of complex size. The linearization of 2D circular complexes into 1D structures and a reverse cyclization process were demonstrated using a toehold-mediated approach. The approach proposed here for the construction and directed modification of the molecularity and shape of complexes will be a valuable tool in RNA nanotechnology, especially for the rational design of therapeutic nucleic acids with high target specificity and the programmable response of the immune system of organisms.

Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes

To date, many derivatives and analogs of nucleic acids (NAs) have been developed. Some of them have found uses in scientific research and biomedical applications. Their effective use is based on the data about their properties. Some of the most important physicochemical properties of oligonucleotides are thermodynamic parameters of the formation of their duplexes with DNA and RNA. These parameters can be calculated only for a few NA derivatives: locked NAs, bridged oligonucleotides, and peptide NAs. Existing predictive approaches are based on an analysis of experimental data and the consequent construction of predictive models. The ongoing pilot studies aimed at devising methods for predicting the properties of NAs by computational modeling techniques are based only on knowledge about the structure of oligonucleotides. In this work, we studied the applicability of the weighted histogram analysis method (WHAM) in combination with umbrella sampling to the calculation of thermodynamic parameters of DNA duplex formation (changes in enthalpy ΔH°, entropy ΔS°, and Gibbs free energy ΔG° 37). A procedure was designed involving WHAM for calculating the hybridization properties of oligodeoxyribonucleotides. Optimal parameters for modeling and calculation of thermodynamic parameters were determined. The feasibility of calculation of ΔH°, ΔS°, and ΔG° 37 was demonstrated using a representative sample of 21 oligonucleotides 4-16 nucleotides long with a GC content of 14-100 %. Error of the calculation of the thermodynamic parameters was 11.4, 12.9, and 11.8 % for ΔH°, ΔS°, and ΔG° 37, respectively, and the melting temperature was predicted with an average error of 5.5 °C. Such high accuracy of computations is comparable with the accuracy of the experimental approach and of other methods for calculating the energy of NA duplex formation. In this paper, the use of WHAM for computation of the energy of DNA duplex formation was systematically investigated for the first time. Our results show that a reliable calculation of the hybridization parameters of new NA derivatives is possible, including derivatives not yet synthesized. This work opens up new horizons for a rational design of constructs based on NAs for solving problems in biomedicine and biotechnology.

Controlled destabilization of caged circularized DNA oligonucleotides predicted by replica exchange molecular dynamics simulations

Spatiotemporal control is a critical issue in the design of strategies for the photoregulation of oligonucleotide activity. Efficient uncaging, i.e., activation by removal of photolabile protecting groups (PPGs), often necessitates multiple PPGs. An alternative approach is based on circularization strategies, exemplified by intrasequential circularization, also denoted photo-tethering, as introduced in [Seyfried et al., Angew. Chem., Int. Ed., 2017, 56, 359]. Here, we develop a computational protocol, relying on replica exchange molecular dynamics (REMD), in order to characterize the destabilization of a series of circularized, caged DNA oligonucleotides addressed in the aforementioned study. For these medium-sized (32 nt) oligonucleotides, melting temperatures are computed, whose trend is in good agreement with experiment, exhibiting a large destabilization and, hence, reduction of the melting temperature of the order of ΔTm ∼ 30 K as compared with the native species. The analysis of free energy landscapes confirms the destabilization pattern experienced by the circularized oligonucleotides. The present study underscores that computational protocols that capture controlled destabilization and uncaging of oligonucleotides are promising as predictive tools in the tailored photocontrol of nucleic acids.

Dissociative electron attachment to gold(I)-based compounds: 4,5-dichloro-1,3-diethyl-imidazolylidene trifluoromethyl gold(I)

With the use of proton-NMR and powder XRD (XRPD) studies, the suitability of specific Au-focused electron beam induced deposition (FEBID) precursors has been investigated with low electron energy, structure, excited states and resonances, structural crystal modifications, flexibility, and vaporization level. 4,5-Dichloro-1,3-diethyl-imidazolylidene trifluoromethyl gold(I) is a compound that is a uniquely designed precursor to meet the needs of focused electron beam-induced deposition at the nanostructure level, which proves its capability in creating high purity structures, and its growing importance in other AuIm_x and AuCl_nB (where x and n are the number of radicals, B = CH, CH₃, or Br) compounds in the radiation cancer therapy increases the efforts to design more suitable bonds in processes of SEM (scanning electron microscopy) deposition and in gas-phase studies. The investigation performed of its powder shape using the XRPD XPERT³ panalytical diffractometer based on CoK_α lines shows changes to its structure with change in temperature, level of vacuum, and light; the sensitivity of this compound makes it highly interesting in particular to the radiation research. Used in FEBID, though its smaller number of C, H, and O atoms has lower levels of C contamination in the structures and on the surface, it replaces these bonds with C-Cl and C-N bonds that have lower bond-breaking energy. However, it still needs an extra purification step in the deposition process, either H₂O, O₂, or H jets.

Controlled Release of DNA Binding Anticancer Drugs from Gold Nanoparticles with Near-Infrared Radiation

Traditional chemotherapies target rapidly developing cells in the human body resulting in harsh side effects including fatigue, immune system suppression, and nausea, among others. Delivery systems to focus the active pharmaceutical ingredients (APIs) to the diseased tissue can diminish the negative side effects while improving treatment outcomes. Gold nanoparticles (AuNP) offer many unique advantages as drug delivery vehicles, including being biologically inert, easily adaptable to various shapes and sizes, able to create a strong Au-thiol bond, and able to generate heat upon the absorption of near-infrared light. To this end, a AuNP delivery vehicle was engineered to load and release two DNA binding anti-cancer drugs, mithramycin and doxorubicin, in a controlled fashion. The drugs were loaded onto the surface of the AuNP with temperature sensitive linkages. The amount of heat generated, and subsequent release of the drugs was controlled by the irradiation time with a near-infrared laser. By modulating the linkage used to load the drugs three different release profiles were able to be achieved, indicating the feasibility of such a system for combinational therapy requiring sequential release of APIs.

Fluorescent aptasensor based on conformational switch-induced hybridization for facile detection of β-amyloid oligomers

Aβ oligomers (AβO) are a dominant biomarker for early Alzheimer's disease diagnosis. A fluorescent aptasensor coupled with conformational switch-induced hybridization was established to detect AβO. The fluorescent aptasensor is based on the interaction of fluorophore-labeled AβO-specific aptamer (FAM-Apt) against its partly complementary DNA sequence on the surface of magnetic beads (cDNA-MBs). Once the FAM-Apt binds to AβO, the conformational switch of FAM-Apt increases the tendency to be captured by cDNA-MBs. This causes a descending fluorescence of supernatant, which can be utilized to determine the levels of AβO. Thus, the base-pair matching above 12 between FAM-Apt and cDNA-MBs with increasing hybridizing free energies reached the ascending fluorescent signal equilibrium. The optimized aptasensor showed linearity from 1.7 ng mL^-1 to 85.1 (R = 0.9977) with good recoveries (79.27-109.17%) in plasma. Furthermore, the established aptasensor possesses rational selectivity in the presence of monomeric Aβ, fibrotic Aβ, and interferences. Therefore, the developed aptasensor is capable of quantifying AβO in human plasma and possesses the potential to apply in clinical cases.

Protein DEK and DTA Aptamers: Insight Into the Interaction Mechanisms and the Computational Aptamer Design

By blocking the DEK protein, DEK-targeted aptamers (DTAs) can reduce the formation of neutrophil extracellular traps (NETs) to reveal a strong anti-inflammatory efficacy in rheumatoid arthritis. However, the poor stability of DTA has greatly limited its clinical application. Thus, in order to design an aptamer with better stability, DTA was modified by methoxy groups (DTA_OMe) and then the exact DEK–DTA interaction mechanisms were explored through theoretical calculations. The corresponding 2′-OCH3-modified nucleotide force field was established and the molecular dynamics (MD) simulations were performed. It was proved that the 2′-OCH3-modification could definitely enhance the stability of DTA on the premise of comparative affinity. Furthermore, the electrostatic interaction contributed the most to the binding of DEK–DTA, which was the primary interaction to maintain stability, in addition to the non-specific interactions between positively-charged residues (e.g., Lys and Arg) of DEK and the negatively-charged phosphate backbone of aptamers. The H-bond network analysis reminded that eight bases could be mutated to probably enhance the affinity of DTA_OMe. Therein, replacing the 29th base from cytosine to thymine of DTA_OMe was theoretically confirmed to be with the best affinity and even better stability. These research studies imply to be a promising new aptamer design strategy for the treatment of inflammatory arthritis.

Pairing nanoarchitectonics of oligodeoxyribonucleotides with complex diversity: concatemers and self-limited complexes

The development of approaches to the design of two- and three-dimensional self-assembled DNA-based nanostructures with a controlled shape and size is an essential task for applied nanotechnology, therapy, biosensing, and bioimaging. We conducted a comprehensive study on the formation of various complexes from a pair of oligonucleotides with two transposed complementary blocks that can be linked through a nucleotide or non-nucleotide linker. A methodology is proposed to prove the formation of a self-limited complex and to determine its molecularity. It is based on the "opening" of a self-limited complex with an oligonucleotide that effectively binds to a duplex-forming block. The complexes assembled from a pair of oligonucleotides with different block length and different linker sizes and types were investigated by theoretical analysis, several experimental methods (a gel shift assay, atomic force microscopy, and ultraviolet melting analysis), and molecular dynamics simulations. The results showed a variety of complexes formed by only a pair of oligonucleotides. Self-limited associates, concatemer complexes, or mixtures thereof can arise if we change the length of a duplex and loop-forming blocks in oligonucleotides or via introduction of overhangs and chemical modifications. We postulated basic principles of rational design of native self-limited DNA complexes of desired structure, shape, and molecularity. Our foundation makes self-limited complexes useful tools for nanotechnology, biological studies, and therapeutics.

Zeptomole detection of DNA based on microparticle dissociation from a glass plate in a combined acoustic-gravitational field

In this study, we propose a novel detection principle based on the dissociation of microparticles immobilized on a glass plate through weak hybridization involving 4-6 base pairs (bps) in a combined acoustic-gravitational field. Particle dissociation from the glass plate occurs when the resultant of the acoustic radiation force (F_ac) and the sedimentation force (F_sed) exerted on the particle exceeds the binding force owing to the weak hybridization (F_bind). Because F_ac and F_sed can be controlled by the microparticle density, and F_ac is a function of the applied voltage to the transducer (V), an increase in V induces particle dissociation. The binding of gold nanoparticles (AuNPs) onto silica microparticles (SPs) resulting from the strong hybridization of 20 bps induces an increase in the density of SPs, leading to an increase in F_ac and F_sed; consequently, the voltage V required for dissociation becomes lower than that required without AuNP binding. We demonstrate that the dependence of the binding number of AuNPs per SP on V follows the theoretical prediction. The binding of 7500 AuNPs per SP can be detected as a 10 V change in V. The present approach allows the detection of 2000 DNA molecules involved in the strong hybridization between AuNPs and SP.

Calculation of Energy for RNA/RNA and DNA/RNA Duplex Formation by Molecular Dynamics Simulation

Abstract

The development of approaches for predictive calculation of hybridization properties of various nucleic acid (NA) derivatives is the basis for the rational design of the NA-based constructs. Modern advances in computer modeling methods provide the feasibility of these calculations. We have analyzed the possibility of calculating the energy of DNA/RNA and RNA/RNA duplex formation using representative sets of complexes (65 and 75 complexes, respectively). We used the classical molecular dynamics (MD) method, the MMPBSA or MMGBSA approaches to calculate the enthalpy (ΔH°) component, and the quasi-harmonic approximation (Q-Harm) or the normal mode analysis (NMA) methods to calculate the entropy (ΔS°) contribution to the Gibbs energy ($$\Delta G_{{37}}^{^\circ }$$ ) of the NA complex formation. We have found that the MMGBSA method in the analysis of the MD trajectory of only the NA duplex and the empirical linear approximation allow calculation of the enthalpy of formation of the DNA, RNA, and hybrid duplexes of various lengths and GC content with an accuracy of 8.6%. Within each type of complex, the combination of rather efficient MMGBSA and Q-Harm approaches being applied to the trajectory of only the bimolecular complex makes it possible to calculate the $$\Delta G_{{37}}^{^\circ }$$ of the duplex formation with an error value of 10%. The high accuracy of predictive calculation for different types of natural complexes (DNA/RNA, DNA/RNA, and RNA/RNA) indicates the possibility of extending the considered approach to analogs and derivatives of nucleic acids, which gives a fundamental opportunity in the future to perform rational design of new types of NA-targeted sequence-specific compounds.

Recognition and removal of clustered DNA lesions via nucleotide excision repair

Clustered damage of DNA consists of two or more lesions located within one or two turns of the DNA helix. Clusters consisting of lesions of various structures can arise under the influence of strong damaging factors, especially if the cells have a compromised repair status. In this work, we analyzed how the presence of an analog of the apurinic/apyrimidinic site - a non-nucleoside residue consisting of diethylene glycol phosphodiester (DEG) - affects the recognition and removal of a bulky lesion (a non-nucleoside site of the modified DNA strand containing a fluorescein residue, nFlu) from DNA by a mammalian nucleotide excision repair system. Here we demonstrated that the efficiency of nFlu removal decreases in the presence of DEG in the complementary strand and is completely suppressed when the DEG is located opposite the nFlu. By contrast, protein factor XPC-RAD23B, which initiates global genomic nucleotide excision repair, has higher affinity for DNA containing clustered damage as compared to DNA containing a single bulky lesion; the affinity of XPC strengthens as the positions of DEG and nFlu become closer. The changes in the double-stranded DNA's geometry caused by the presence of clustered damage were also assessed. The obtained experimental data together with the results of molecular dynamics simulations make it possible to get insight into the structural features of DNA containing clustered lesions that determine the efficiency of repair. Speaking more broadly, this study should help to understand the probable fate of bulky adduct-containing clusters of various topologies in the mammalian cell.

Prediction of PCR amplification from primer and template sequences using recurrent neural network

We have developed a novel method to predict the success of PCR amplification for a specific primer set and DNA template based on the relationship between the primer sequence and the template. To perform the prediction using a recurrent neural network, the usual double-stranded formation between the primer and template nucleotide sequences was herein expressed as a five-lettered word. The set of words (pseudo-sentences) was placed to indicate the success or failure of PCR targeted to learn recurrent neural network (RNN). After learning pseudo-sentences, RNN predicted PCR results from pseudo-sentences which were created by primer and template sequences with 70% accuracy. These results suggest that PCR results could be predicted using learned RNN and the trained RNN could be used as a replacement for preliminary PCR experimentation. This is the first report which utilized the application of neural network for primer design and prediction of PCR results.

Effects of Phosphoryl Guanidine Modification of Phosphate Residues on the Structure and Hybridization of Oligodeoxyribonucleotides

Phosphoryl guanidine oligonucleotides (PGOs) are promising tools for biological research and development of biosensors and therapeutics. We performed structural and hybridization analyses of octa-, deca-, and dodecamers with all phosphate residues modified by 1,3-dimethylimidazolidine-2-imine moieties. Similarity of the B-form double helix between native and modified duplexes was noted. In PGO duplexes, we detected a decrease in the proportion of C2'-endo and an increased proportion of C1'-exo sugar conformations of the modified chain. Applicability of the two-state model to denaturation transition of all studied duplexes was proved for the first time. Sequence-dependent effects of this modification on hybridization properties were observed. The thermal stability of PGO complexes is almost native at 100 mM NaCl and slightly increases with decreasing ionic strength. An increase in water activity and dramatic changes in interaction with cations and in solvation of PGOs and their duplexes were noted, resulting in slight elevation of the melting temperature after an ionic-strength decrease from 1 M NaCl down to deionized water. Decreased binding of sodium ions and decreased water solvation were documented for PGOs and their duplexes. In contrast to DNA, the PGO duplex formation leads to a release of several cations. The water shell is significantly more disordered near PGOs and their complexes. Nevertheless, changes in solvation during the formation of native and PGO complexes are similar and indicate that it is possible to develop models for predictive calculations of the thermodynamic properties of phosphoryl guanidine oligomers. Our results may help devise an approach for the rational design of PGOs as novel improved molecular probes and tools for many modern methods involving oligonucleotides.

Thermodynamics of DNA Hybridization from Atomistic Simulations

Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.

Dynamic and scalable DNA-based information storage

The physical architectures of information storage systems often dictate how information is encoded, databases are organized, and files are accessed. Here we show that a simple architecture comprised of a T7 promoter and a single-stranded overhang domain (ss-dsDNA), can unlock dynamic DNA-based information storage with powerful capabilities and advantages. The overhang provides a physical address for accessing specific DNA strands as well as implementing a range of in-storage file operations. It increases theoretical storage densities and capacities by expanding the encodable sequence space and simplifies the computational burden in designing sets of orthogonal file addresses. Meanwhile, the T7 promoter enables repeatable information access by transcribing information from DNA without destroying it. Furthermore, saturation mutagenesis around the T7 promoter and systematic analyses of environmental conditions reveal design criteria that can be used to optimize information access. This simple but powerful ss-dsDNA architecture lays the foundation for information storage with versatile capabilities.

The physical architectures of information storage dictate how data is encoded, organised and accessed. Here the authors use DNA with a single-strand overhang as a physical address to access specific data and do in-storage file operations in a scalable and reusuable manner.

Spectroscopy, Thermodynamics and Molecular Docking of Fraxinellone with DNA

Fraxinellone is an important botanical lactone compound and has been demonstrated to have insecticidal activity. To provide theoretical support to the assessment on the safety of utilizing fraxinellone as a natural insecticidal agent, the interactions between fraxinellone and armyworm DNA, salmon sperm DNA and calf thymus DNA were investigated using UV-Vis absorption spectroscopy, isothermal titration calorimetry, and molecular docking. Results showed that there were two types of combinations between fraxinellone and three kinds of DNA. Type I combination had an equilibrium constant of combination (Ka₁) of about 10⁵ and binding sites (n₁) of 0.40-0.70, while type II combination had an equilibrium constant of combination (Ka₂) of 10³ and binding sites (n₂) of 1.35-3.15. Results of molecular docking showed that there were non-classical embedding type interactions between fraxinellone and three kinds of DNA, with the reaction taking place in small groove areas of the DNA structure, resulting in relatively weak interactive forces.

DNA complexes with human apurinic/apyrimidinic endonuclease 1: structural insights revealed by pulsed dipolar EPR with orthogonal spin labeling

A DNA molecule is under continuous influence of endogenous and exogenous damaging factors, which produce a variety of DNA lesions. Apurinic/apyrimidinic sites (abasic or AP sites) are among the most common DNA lesions. In this work, we applied pulse dipolar electron paramagnetic resonance (EPR) spectroscopy in combination with molecular dynamics (MD) simulations to investigate in-depth conformational changes in DNA containing an AP site and in a complex of this DNA with AP endonuclease 1 (APE1). For this purpose, triarylmethyl (TAM)-based spin labels were attached to the 5' ends of an oligonucleotide duplex, and nitroxide spin labels were introduced into APE1. In this way, we created a system that enabled monitoring the conformational changes of the main APE1 substrate by EPR. In addition, we were able to trace substrate-to-product transformation in this system. The use of different (orthogonal) spin labels in the enzyme and in the DNA substrate has a crucial advantage allowing for detailed investigation of local damage and conformational changes in AP-DNA alone and in its complex with APE1.

Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides

MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.

The Role of Structural Enthalpy in Spherical Nucleic Acid Hybridization

DNA hybridization onto DNA-functionalized nanoparticle surfaces (e.g., in the form of a spherical nucleic acid (SNA)) is known to be enhanced relative to hybridization free in solution. Surprisingly, via isothermal titration calorimetry, we reveal that this enhancement is enthalpically, as opposed to entropically, dominated by ∼20 kcal/mol. Coarse-grained molecular dynamics simulations suggest that the observed enthalpic enhancement results from structurally confining the DNA on the nanoparticle surface and preventing it from adopting enthalpically unfavorable conformations like those observed in the solution case. The idea that structural confinement leads to the formation of energetically more stable duplexes is evaluated by decreasing the degree of confinement a duplex experiences on the nanoparticle surface. Both experiment and simulation confirm that when the surface-bound duplex is less confined, i.e., at lower DNA surface density or at greater distance from the nanoparticle surface, its enthalpy of formation approaches the less favorable enthalpy of duplex formation for the linear strand in solution. This work provides insight into one of the most important and enabling properties of SNAs and will inform the design of materials that rely on the thermodynamics of hybridization onto DNA-functionalized surfaces, including diagnostic probes and therapeutic agents.

Direct Comparison of Amino Acid and Salt Interactions with Double-Stranded and Single-Stranded DNA from Explicit-Solvent Molecular Dynamics Simulations

Given the ubiquitous nature of protein-DNA interactions, it is important to understand the interaction thermodynamics of individual amino acid side chains for DNA. One way to assess these preferences is to perform molecular dynamics (MD) simulations. Here we report MD simulations of 20 amino acid side chain analogs interacting simultaneously with both a 70-base-pair double-stranded DNA and with a 70-nucleotide single-stranded DNA. The relative preferences of the amino acid side chains for dsDNA and ssDNA match well with values deduced from crystallographic analyses of protein-DNA complexes. The estimated apparent free energies of interaction for ssDNA, on the other hand, correlate well with previous simulation values reported for interactions with isolated nucleobases, and with experimental values reported for interactions with guanosine. Comparisons of the interactions with dsDNA and ssDNA indicate that, with the exception of the positively charged side chains, all types of amino acid side chain interact more favorably with ssDNA, with intercalation of aromatic and aliphatic side chains being especially notable. Analysis of the data on a base-by-base basis indicates that positively charged side chains, as well as sodium ions, preferentially bind to cytosine in ssDNA, and that negatively charged side chains, and chloride ions, preferentially bind to guanine in ssDNA. These latter observations provide a novel explanation for the lower salt dependence of DNA duplex stability in GC-rich sequences relative to AT-rich sequences.