چكيده لاتين
Following the discovery of ribozymes as nucleic acid-based enzymes, researchers were directed toward the identification and development of DNA-based catalytic molecules, known as deoxyribozymes. Since the recognition of the first DNAzyme with cleavage activity, various deoxyribozymes with diverse catalytic functions against different substrates have been identified and developed under in vitro conditions, and have subsequently been employed as powerful molecular tools in a wide range of applications. Similar to other enzymes, DNAzymes must exhibit a specific degree of structural stability in order to achieve optimal catalytic activity under different reaction conditions. Due to the intrinsic relationship between structure and catalytic activity, the double-stranded structure formed between the binding arms of the DNAzyme—designed as the reverse complement of the substrate sequence—must be arranged both spatially and thermodynamically in such a way that proper folding of the active site facilitates the transition state of the reaction with lower activation energy, thereby enhancing the catalytic efficiency of the DNAzyme.
The stability of the double-stranded structures formed between DNAzyme binding arms and the substrate depends on multiple parameters. One of the most widely examined thermodynamic parameters, predicted and measured in numerous studies as an indicator of structural stability, is the total binding free energy of the duplex (ΔGbinding). Accordingly, the stability of Watson–Crick base-paired duplexes is typically described and calculated in terms of the binding free energy (ΔG). Factors influencing ΔG include the nearest-neighbor nucleotide composition within the duplex sequence, the overall length of both the substrate and the binding arms, the relative positioning of nucleotides with respect to the catalytic domain, the presence of one or more unpaired nucleotides at the enzyme–substrate interface, and their positional context relative to the catalytic core.
In this study, to simultaneously investigate the influence of binding-arm length and the presence of single-nucleotide substitutions within the duplex structure on binding free energy—considered as two critical thermodynamic parameters affecting catalytic activity—two variants of the F8 DNAzyme, differing in arm length, were designed against eight single-nucleotide mutant substrates (previously employed in the Nucleic Acid Laboratory at the University of Isfahan under the supervision of Jowadi). Furthermore, the NUPACK Python API was employed to calculate thermodynamic parameters, including the total binding free energy (ΔG) of mutant-containing duplexes with different lengths, the change in free energy (ΔΔG), and the normalized index |ΔΔG / ΔG(wt)|, which was used to quantitatively assess the sensitivity (tolerance or intolerance) of the duplex structures to single-nucleotide substitutions. This approach enabled a precise evaluation of the extent to which different mutations, as well as sequence length, influence duplex stability and binding strength, and consequently their impact on catalytic efficiency.