DNA Molecules Research Articles

Endemic Radiation of African Moonfish, Selene dorsalis (Gill 1863), in the Eastern Atlantic: Mitogenomic Characterization and Phylogenetic Implications of Carangids (Teleostei: Carangiformes).

This study offers an in-depth analysis of the mitochondrial genome of Selene dorsalis (Gill 1863), a species native to the Eastern Atlantic Ocean. The circular mitochondrial DNA molecule measures 16,541 base pairs and comprises 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and a control region (CR). The nucleotide composition exhibits a notable adenine-thymine (AT) bias, accounting for 53.13%, which aligns with other species in the Carangidae family. Most PCGs initiate with the ATG codon, with the exception of Cytochrome C oxidase subunit I, which starts with GTG. Analysis of relative synonymous codon usage reveals that leucine and serine are the most prevalent amino acids in the mitochondrial genome of S. dorsalis and its congeners (S. vomer and S. setapinnis). All tRNAs display the typical cloverleaf structure, though tRNA Serine (S1) lacks a dihydrouracil arm. Pairwise comparisons of synonymous and nonsynonymous substitutions for all PCGs yielded values below '1', indicating strong purifying selection. The CR spans 847 bp, representing 5.12% of the mitochondrial genome, and is characterized by high AT content (62.81%). It is situated between tRNA-Pro (TGG) and tRNA-Phe (GAA). The CR contains conserved sequence blocks, with CSB-1 being the longest at 22 bp and CSB-D the shortest at 18 bp. Phylogenetic analysis, using Bayesian and Maximum-likelihood trees constructed from concatenated PCGs across 72 species, successfully differentiates S. dorsalis from other carangids. This study also explores how ocean currents and gyres might influence lineage diversification and parapatric speciation of Selene species between the Atlantic and Pacific Oceans. These results highlight the importance of the mitochondrial genome in elucidating the structural organization and evolutionary dynamics of S. dorsalis and its relatives within marine ecosystems.

Genome Instability Induced by Topoisomerase Misfunction.

Topoisomerases alter DNA topology by making transient DNA strand breaks (DSBs) in DNA. The DNA cleavage reaction mechanism includes the formation of a reversible protein/DNA complex that allows rapid resealing of the transient break. This mechanism allows changes in DNA topology with minimal risks of persistent DNA damage. Nonetheless, small molecules, alternate DNA structures, or mutations in topoisomerase proteins can impede the resealing of the transient breaks, leading to genome instability and potentially cell death. The consequences of high levels of enzyme/DNA adducts differ for type I and type II topoisomerases. Top1 action on DNA containing ribonucleotides leads to 2-5 nucleotide deletions in repeated sequences, while mutant Top1 enzymes can generate large deletions. By contrast, small molecules that target Top2, or mutant Top2 enzymes with elevated levels of cleavage lead to small de novo duplications. Both Top1 and Top2 have the potential to generate large rearrangements and translocations. Thus, genome instability due to topoisomerase mis-function is a potential pathogenic mechanism especially leading to oncogenic progression. Recent studies support the potential roles of topoisomerases in genetic changes in cancer cells, highlighting the need to understand how cells limit genome instability induced by topoisomerases. This review highlights recent studies that bear on these questions.

Disparate requirements for RAD54L in replication fork reversal.

RAD54L is a DNA motor protein with multiple roles in homologous recombination DNA repair. In vitro, RAD54L was shown to also catalyze the reversal and restoration of model replication forks. In cells, however, little is known about how RAD54L may regulate the dynamics of DNA replication. Here, we show that RAD54L restrains the progression of replication forks and functions as a fork remodeler in human cancer cell lines and non-transformed cells. Analogous to HLTF, SMARCAL1and FBH1, and consistent with a role in fork reversal, RAD54L decelerates fork progression in response to replication stress and suppresses the formation of replication-associated ssDNA gaps. Interestingly, loss of RAD54L prevents nascent strand DNA degradation in both BRCA1/2- and 53BP1-deficient cells, suggesting that RAD54L functions in both pathways of RAD51-mediated replication fork reversal. In the HLTF/SMARCAL1 pathway, RAD54L is critical, but its ability to catalyze branch migration is dispensable, indicative of its function downstream of HLTF/SMARCAL1. Conversely, in the FBH1 pathway, branch migration activity of RAD54L is essential, and FBH1 engagement is dependent on its concerted action with RAD54L. Collectively, our results reveal disparate requirements for RAD54L in two distinct RAD51-mediated fork reversal pathways, positing its potential as a future therapeutic target.

Identification and molecular analysis of watermelon chlorotic stunt virus infecting snake gourd in Saudi Arabia

Snake gourd (Trichosanthes cucumerina) plants exhibiting typical begomovirus-like symptoms of stunted growth, leaf yellowing and mottling were observed at an open field in the Eastern region of Saudi Arabia. Sequencing analysis of the amplified complete DNA molecules revealed that the plants were infected with watermelon chlorotic stunt virus (WmCSV), which is a bipartite begomovirus prevalent mostly in the Old World and a serious threat to cucurbit production in the Arabian Peninsula, Middle East and Africa. The two WmCSV DNA-A isolates (SG31A and SG52A) were 98.9% identical and showed their highest nucleotide (nt) sequence identities (98.7%) with the isolates from Iran and Saudi Arabia. The DNA-B isolates (SG31B and SG52B), on the other hand, were 97.4% identical and exhibited their highest nt sequence identities (99.5 and 97%) with isolates reported from Iran and Oman. In the phylogenetic dendrograms the identified isolates clustered closely with previously reported WmCSV isolates from Iran and Saudi Arabia. Infectivity assays revealed that the DNA-A components alone could not induce infection in Nicotiana benthamiana plants however, together with DNA-B these isolates successfully caused typical begomovirus symptoms and both components were detected successfully using Southern blot hybridization. This study highlights the importance of conducting extensive future begomovirus surveillance to detect spillover events that could threaten native vegetable production in Saudi Arabia. This is crucial as begomoviruses pose a serious threat to vegetable cultivation throughout the Middle East.

The Nucleolus and Its Interactions with Viral Proteins Required for Successful Infection.

Nuclear bodies are structures in eukaryotic cells that lack a plasma membrane and are considered protein condensates, DNA, or RNA molecules. Known nuclear bodies include the nucleolus, Cajal bodies, and promyelocytic leukemia nuclear bodies. These bodies are involved in the concentration, exclusion, sequestration, assembly, modification, and recycling of specific components involved in the regulation of ribosome biogenesis, RNA transcription, and RNA processing. Additionally, nuclear bodies have been shown to participate in cellular processes such as the regulation of transcription of the cell cycle, mitosis, apoptosis, and the cellular stress response. The dynamics and functions of these bodies depend on the state of the cell. It is now known that both DNA and RNA viruses can direct their proteins to nuclear bodies, causing alterations in their composition, dynamics, and functions. Although many of these mechanisms are still under investigation, it is well known that the interaction between viral and nuclear body proteins is necessary for the success of the viral infection cycle. In this review, we concisely describe the interaction between viral and nuclear body proteins. Furthermore, we focus on the role of the nucleolus in RNA virus infections. Finally, we discuss the possible implications of the interaction of viral proteins on cellular transcription and the formation/degradation of non-coding RNAs.

Sequence-Unconstrained DNA Computing: DSN cycling and PER circuitry for dynamic miRNAs analysis and multifunctional logic operations

DNA strand displacement and enzyme-driven catalysis are pivotal in regulating molecular interactions and constructing molecular logic circuits. However, their application in intricate hierarchical cascading networks has been hindered by sequence crosstalk restrictions on the multiple inputs and outputs. In this study, we present a sequence-unconstrained strategy to implement versatile molecular logic circuits, relying on duplex-specific nuclease (DSN)-assisted cycling and primer exchange reaction (PER) amplification circuitry (DSN-PER). Employing linear capture probe (CP) conjugated with magnetic beads (MBs) significantly minimizes crosstalk reactions compared to the previous molecular beacon probe (MBP) for DNA barcoding. Importantly, the sequence design of “barcodes” is fully independent of target binding sites on MBs, allowing for the sensitive detection of a broad range of microRNAs (miRNAs) using the DSN-PER method. Specifically, this platform supports the programming of common logic computations (YES, NOT, OR, INHIBIT, AND, and NOR) and hierarchical circuit operation (OR-AND) in response to a specific pattern of varying miRNAs. These versatile logic circuits offer unique features with flexible sequence design, potentially broadening the applications of dynamic nucleic acid nanotechnologies in multiplexed bio-computation, intelligent diagnostics, and sophisticated DNA nanodevices.

Equipping Uridine with Three Dipeptide Motifs for the Inhibition of Radical-Induced DNA Oxidation.

We report the discovery and characterization of antioxidative effects of uridine linked with three dipeptide motifs against DNA oxidation induced by peroxyl radicals. First, the dipeptide motifs are constructed by using the Ugi four-component reaction (Ugi 4CR), in which caffeic, ferulic, sinapic, and syringic acids are used as the carboxylic acid resources, vanillin, benzaldehyde, and p-hydroxybenzaldehyde are used as the aldehyde resources, tyramine- and dopamine-related isocyanides as well as ethyl isocyanoacetate are used as the isocyanide resources, and 2-(p-aminophenyl)ethanol is used as the amine component. We found that the antioxidative effects of the Ugi 4CR products are 1.3-2.8 times higher than those of caffeic, ferulic, sinapic, and syringic acids in the protection of DNA against peroxyl radical-induced oxidation. Moreover, when three Ugi 4CR products are linked with three hydroxyl groups of uridine by using three succinic anhydrides as the linkage, the inhibitory effects of the afforded uridine-dipeptide hybrids against the DNA oxidation increase 4.4-8.9 times (>3 times) compared to that of the Ugi 4CR product. This is due to the hybrid structure consisting of uridine and three motifs of the Ugi 4CR product enabling binding with the DNA strand more efficiently and quenching free radicals more rapidly. Therefore, the hybrid structure constructed by the nucleoside with antioxidative dipeptides offers an additional advantage for protecting DNA against radical-induced oxidation.

Predicting DNA Reactions with a Quantum Chemistry‐Based Deep Learning Model

AbstractIn this study, a deep learning model based on quantum chemistry is introduced to enhance the accuracy and efficiency of predicting DNA reaction parameters. By integrating quantum chemical calculations with self‐designed descriptor matrices, the model offers a comprehensive description of energy variations and considers a broad range of relevant factors. To overcome the challenge of limited labeled data, an active learning method is employed. The results demonstrate that this model outperforms existing methods in predicting DNA hybridization free energies and strand displacement rate constants, thus advancing the understanding of DNA molecular interactions, and aiding in the precise design and optimization of DNA‐based systems.

Strand dependent bypass of DNA lesions during fork reversal by ATP-dependent translocases SMARCAL1, ZRANB3, and HLTF.

During DNA replication, the replisome encounters obstacles including DNA lesions, transcription-replication conflicts, and other sources of replication stress. These obstacles must be efficiently overcome to complete DNA synthesis and minimize genome instability. One pathway to tolerate replication stress is replication fork reversal, in which parental template DNA strands are reannealed and a nascent-nascent DNA duplex is formed. Several enzymes promote replication fork reversal, including the ATP-dependent translocases SMARCAL1, ZRANB3, and HLTF. How these enzymes translocate on DNA that contains fork-stalling lesions is unknown. Here, we examined the abilities of SMARCAL1, ZRANB3, and HLTF to tolerate various lesions on leading or lagging template strands. We demonstrate that SMARCAL1 and ZRANB3 are selectively inhibited by lesions on the leading template strand, whereas HLTF is insensitive to bulky lesions on either strand. These results suggest that SMARCAL1 and ZRANB3 contact the leading strand during fork reversal and therefore are more sensitive to inhibition by bulky lesions on this strand. In contrast, HLTF DNA translocation is inherently insensitive to DNA lesions. These biochemical differences between the fork reversal enzymes provide insights into their mechanism of DNA remodeling and suggest they may act in lesion-specific contexts.

DciA secures bidirectional replication initiation in Vibrio cholerae.

Replication is initiated bidirectionally in the three domains of life by the assembly of two replication forks at an origin of replication. This is made possible by the recruitment of two replicative helicases to a nucleoprotein platform built at the origin of replication with the initiator protein. The reason why replication is initiated bidirectionally has never been experimentally addressed due to the lack of a suitable biological system. Using genetic and genomic approaches, we show that upon depletion of DciA, replication is no longer initiated bidirectionally at the origin of replication of Vibrio cholerae chromosome 1. We show that following unidirectional replication on the left replichore, nascent DNA strands at ori1 anneal to each other to form a double-stranded DNA end. While this DNA end can be efficiently resected in recB+ cells, only a few cells use it to trigger replication on the right replichore. In most DciA-depleted cells, chromosome 1 is degraded leading to cell death. Our results suggest that DciA is essential to ensuring bidirectional initiation of replication in bacteria, preventing a cascade of deleterious events following unidirectional replication initiation.

A double-ring of human RAD52 remodels replication forks restricting fork reversal.

Using cryo-EM, biochemical and single-molecule approaches we show that the structure of stalled DNA replication fork promotes a unique two-ring organization of human RAD52 protein which remodels the fork via DNA strand exchange.

Breast Tumor Diagnosis Based on Molecular Learning Vector Quantization Neural Networks.

DNA nanotechnology plays a crucial role in precise cancer medicine. Currently, molecular logic circuits are applied to detect tumor-specific biomarkers and control the release of therapeutic drugs. However, these systems lack self-learning capabilities for intelligent diagnostics in biological samples, and their data processing capabilities are limited. Here, a molecular learning vector quantization neural network (LVQNN) model based on DNA strand displacement (DSD) technology for breast tumor diagnosis is developed. Compared to previous work, the molecular LVQNN boasts powerful computing abilities, handling high-dimensional data for intelligent cancer diagnosis. To verify the feasibility and versatility of the network, two distinct typical datasets are selected: one from a single source with cell morphology data from 569 cases, and a more extensive one spanning different populations and ages, with miRNA gene expression data from 1881 cases. By using the molecular LVQNN, diagnostic experiments are conducted on 50 and 120 public individuals from these two datasets, respectively, achieving accuracy rates of 94% and 97.5%. This study demonstrates that the LVQNN model exhibits significant advantages in breast cancer diagnosis and enhances diagnostic accuracy while introducing new approaches for intelligent cancer diagnosis, anticipated to bring significant breakthroughs and application prospects to precise cancermedicine.

Solvation free energy in governing equations for DNA hybridization, protein-ligand binding, and protein folding.

This work examines the thermodynamics of model biomolecular interactions using a governing equation that accounts for the participation of bulk water in the equilibria. In the first example, the binding affinities of two DNA duplexes, one of nine and one of 10 base pairs in length, are measured and characterized by isothermal titration calorimetry (ITC) as a function of concentration. The results indicate that the change in solvation free energy that accompanies duplex formation (ΔGS) is large and unfavorable. The duplex with the larger number of G:C pairings yields the largest change in solvation free energy, ΔGS = +460 kcal·mol-1per base pair at 25 °C. A van't Hoff analysis of the data is complicated by the varying degree of intramolecular base stacking within each DNA strand as a function of temperature. A modeling study reveals how the solvation free energy alters the output of a typical ITC experiment and leads to a good, though misleading, fit to the classical equilibrium equation. The same thermodynamic framework is applied to a model protein-ligand interaction, the binding of ribonuclease A with the nucleotide inhibitor 3'-UMP, and to a conformational equilibrium, the change in tertiary structure of α-lactalbumin in molar guanidinium chloride solutions. The ribonuclease study yields a value of ΔGS = +160 kcal·mol-1, whereas the folding equilibrium yields ΔGS ≈ 0, an apparent characteristic of hydrophobic interactions. These examples provide insight on the role of solvation energy in binding equilibria and suggest a pivot in the fundamental application of thermodynamics to solution chemistry.

An Improved Method and Device for Nucleic Acid Isolation Using a High-Salt Gel Electroelution Trap.

The success of DNA analytical methods, including long-read sequencing, depends on the availability of high-quality, purified DNA. Previously, we developed a method and device for isolating high-molecular-weight (HMW) DNA for long-read sequencing using a high-salt gel electroelution trap. Here, we present an improved version of this method for purifying nucleic acids with high yield and purity from even the most challenging biological samples. The proposed method is a significant improvement over the previously published procedure, offering a simple, fast, and efficient solution for isolating HMW DNA and smaller DNA and RNA molecules. The method utilizes vertical gel electrophoresis in two nested, partially overlapping electrophoretic columns. The upper, smaller-diameter column has a thin layer of agarose gel at the bottom, which separates nucleic acids from impurities, and an electrophoresis buffer on top. After the target nucleic acid has been gel-purified on the upper column, a larger-diameter column with a layer of high-salt gel overlaid with electrophoresis buffer is inserted from below. The purified nucleic acid is then electroeluted into the buffer-filled gap between the separating gel and the high-salt gel, where excess counterions from the high-salt gel slow its migration and cause it to accumulate. The proposed vertical purification system outperforms the previously described horizontal system in terms of ease of use, speed, scalability, and compatibility with high-throughput workflows. Furthermore, the vertical system allows for the sequential purification of several nucleic acid species from the same sample using interchangeable salt-gel columns.

Impact of kV-cone beam computed tomography dose on DNA repair mechanisms: A pilot study

PurposeIn head and neck squamous cell carcinoma (HNSCC), early complications of the radiotherapy (RT) are observed from the beginning of the treatment to a few months after its end. During external radiotherapy treatment, several patient-dependent parameters can cause a modification of the dose distribution compared to the planned distribution due to variation in patient positioning, anatomy, or intra-fractional movements for example. To verify these parameters during treatment sessions, one of the most commonly used solutions is the cone-beam computed tomography (CBCT). Nowadays, the use of CBCT may constitutes a significant part of the total dose at the end of treatment (up to 10 cGy per session) and more often the volume irradiated by imaging is larger than the one irradiated by the treatment, leading to unintentional irradiation of nearby organs.In this study, we asked whether the imaging low dose added to a following fraction dose (2Gy) may affect the biological response in terms of DNA repair. Material and methodsUsing an IVInomad dosimeter and scintillating fiber probes specially designed for this exploratory study, we exposed fibroblasts cells from head and neck cancer (HNC) patients to a CBCT dose followed by a radiotherapy fraction dose. DNA double strand breaks and DNA repair were assessed by immunofluorescence using the biomarkers gamma H2AX (γH2AX) and pATM. ResultsThe median dose of CBCT was measured between 17 to 21 mGy per session. The kinetics of both biomarkers were found to be strongly dependent on the individual factor in radiosensitive patients. For HNC patients, a prior CBCT dose applied few minutes before the 2Gy dose may have a sublinear effect on the DNA repair mechanisms and potentially on observed health tissue toxicity. ConclusionThe preliminary results obtained highlight the importance of individual and tissue factors for recognizing and repairing DSB during a treatment by radiotherapy using CBCT.

Monitoring SARS-CoV-2 IgA, IgM and IgG antibodies in dried blood and saliva samples using antibody proximity extension assays (AbPEA)

Using a modified proximity extension assay, total and immunoglobulin (Ig) class-specific anti-SARS-CoV-2 antibodies were sensitively and conveniently detected directly from ø1.2 mm discs cut from dried blood and saliva spots (DBS and DSS) without the need for elution. For total Ig detection, antigen probes were prepared by conjugating recombinant spike protein subunit 1 (S1-RBD) to a pair of oligonucleotides. To detect isotype-specific antibody reactivity, one antigen probe was replaced with oligonucleotide-conjugated antibodies specific for antibody isotypes. Binding of pairs of oligonucleotide-conjugated probes to antibodies in patient samples brings oligonucleotides in proximity. An added DNA polymerase uses a transient hybridization between the oligonucleotides to prime synthesis of a DNA strand, which serves as a DNA amplicon that is quantified by real-time PCR. The S1-RBD-specific IgG, IgM, and IgA antibodies in DBS samples collected over the course of a first and second vaccination exhibited kinetics consistent with previous reports. Both DBS and DSS collected from 42 individuals in the autumn of 2023 showed significant level of total S1-RBD antibodies with a correlation of R = 0.70. However, levels in DSS were generally 10 to 100-fold lower than in DBS. Anti-S1-RBD IgG and IgA in DSS demonstrated a correlation of R = 0.6.

Theoretical predictions to facilitate the method development in slalom chromatography for the separation of large DNA molecules

Slalom chromatography (SC) re-emerged in 2024 due to the availability of low adsorption ultra-high pressure liquid chromatography (UHPLC) packed columns/instruments and large modalities being investigated in the context of cell and gene therapies. The physico-chemical principles of SC retention combined with hydrodynamic chromatography (HDC) exclusion have been recently reported. In SC, DNA macromolecules are retarded because: (1) they can be stretched to lengths comparable to the particle diameter, and (2) their elastic relaxation time is long enough to maintain them in non-equilibrium extended conformations under regular UHPLC shear flow conditions. Here, a quantitative HDC-SC retention model is consolidated. A general plate height model accounting for the band broadening of long DNA biopolymers along packed beds is also derived for supporting method development and predicting speed-resolution performance in SC.For illustration, the chromatographic speed-resolution properties in SC are predicted for the separation of specific critical pairs (4.0/4.5, 10/11, and 25/27 kbp) of linear dsDNA polymers. The calculations are performed for two available custom-made particle sizes, dp= 1.7 and 2.5μm, at a constant pressure of 10,000 psi. The predictions are directly validated from experimental data acquired using low adsorption MaxPeakTM 4.6 mm i.d. Columns packed with 1.7μm BEHTM 45 Å (15 cm long column) and 2.5μm BEH 125 Å (30 cm long column) Particles, and by injecting six linear dsDNAs (λ DNA-Hind III Digest). The LC system is very low dispersion ACQUITYTM UPLCTM I-class PLUS System, and the mobile phase is a 100 mM phosphate buffer at pH 8. Maximum resolution is always achieved when the average extended lengths of linear dsDNAs are equal to a critical length, which is proportional to the particle diameter and to the square root of the applied shear rate. Most advantageously, the experimental results reveal that the relaxation times of linear dsDNAs observed under shear flow conditions are two orders of magnitude shorter than those expected in the absence of flow: this enables the detection of the longest linear dsDNAs up to 25 kbp without irremediable loss in column performance. Finally, the retention-efficiency model elaborated in this work can be used to rapidly anticipate and develop methods (selection of particle size, column length, and operating pressure) for any targeted DNA and time-resolution constraints.

Processing and Visualization of Protec-Seq Data for the Generation of Calibrated, Genome-Wide Double Double-Strand Break(dDSB) Maps.

This chapter describes the computational pipeline for the processing and visualization of Protec-Seq data, a method for purification and genome-wide mapping of double-stranded DNA protected by a specific protein at both ends. In the published case, the protein of choice was Saccharomyces cerevisiae Spo11, a conserved topoisomerase-like enzyme that makes meiotic double-strand breaks (DSBs) to initiate homologous recombination, ensuring proper segregation of homologous chromosomes and fertility. The isolated DNA molecules were thus termed double DSB (dDSB) fragments and were found to represent 34 to several hundred base-pair long segments that are generated by Spo11 and are enriched at DSB hotspots, which are sites of topological stress. In order to allow quantitative comparisons between dDSB profiles across experiments, we implemented calibrated chromatin immunoprecipitation sequencing (ChIP-Seq) using the meiosis-competent yeast species Saccharomyces kudriavzevii as calibration strain. Here, we provide a detailed description of the computational methods for processing, analyzing, and visualizing Protec-Seq data, comprising the download of the raw data, the calibrated genome-wide alignments, and the scripted creation of either arc plots or Hi-C-style heatmaps for the illustration of chromosomal regions of interest. The workflow is based on Linux shell scripts (including wrappers for publicly available, open-source software) as well as R scripts and is highly customizable through its modular structure.

DNA and its Mosaic Structure

We consider that the DNA is an information processing system, receiving, registering and transferring information. In DNA we encounter four different nucleotides, which may be considered as the “letters” of an alphabet. Each distinct DNA base may be represented equally well by a specific number and a DNA strand by a unique Gödel number G, where G is a product of prime numbers raised to appropriate powers. We studied the statistical distribution of g, the logarithm of G, in the case of equal probability for the four DNA bases. In the present work we extend our analysis in two different directions We consider unequal probabilities for the DNA bases We raise the issue of scaling. For relatively small DNA strands we use the prime number theorem and we obtain a Poisson distribution for the average g. The Poisson distribution depends upon a parameter λ, which connects together the DNA pattern and the prime number density. Different values for λ lead to a mosaic structure of DNA. Our work indicates a strong link between DNA structure, its representation by Gödel’s numbers and the distribution of the prime numbers.

Majorana zero modes in a G-quadruplex DNA with an s-wave superconductor

We theoretically investigate topological properties of guanine-quadruplex (G4) DNA molecules with an s-wave superconductor. We show the emergence of Majorana zero modes (MZMs) at the ends of G4-DNA with topological nontrivial phase where a Zeeman field along x direction is considered. The topological superconducting G4-DNA can host one pair of MZMs or two pairs of MZMs, which can be regulated by chemical potential. The topological superconducting G4-DNAs with one pair of MZMs are more stable. The differential conductance of topological superconducting G4-DNA shows zero bias peak structure. MZMs can be identified by measuring the crossed Andreev reflection. It can be used to probe and control the MZMs in G4-DNA.