DNA Gene's Basic Structure as a Nonperturbative Circuit Quantum Electrodynamics: Is RNA Polymerase II the Quantum Bus of Transcription?

1. Introduction

Structure-function transcription studies revealed a highly evolutionary conserved process for all three life divisions. Transcription is divided into three steps: promoter DNA binding and RNA initiation, processive RNA chain elongation, and termination [1]. Key points of each step were elucidated, including the nucleotide (nt) addition cycle, and the RNA proofreading mechanism after recognizing errors [2]. Also, a similar architecture, active center, and transcription mechanism of cellular RNA polymerases (RNA pol) have been described [3]. Using nucleoside triphosphate (NTP) substrates, RNA pol synthesizes the RNA complementary to the DNA template strand. The transcription process occurs spatially, forming biomolecular condensates, and RNA pol II can be incorporated into them [2,4].

The qubit, the elementary block for quantum information processing [5], is equivalent to a two-level quantum system (TLS) [6]. Some systems have naturally qubit properties, such as nuclear, atomic, molecular spins, electric dipoles [7], and DNA base pairs (bp). Previously, we carried out a theoretical approach to DNA as a quantum computer [8]. In DNA qubits, Adenine (A) and Thymine (T), Cytosine (C), and Guanine (G) are paired superconductors with the central Hydrogen bonds (H-bond) acting as a Josephson Junction (JJ). The quantum tunneling of H in DNA and RNA has been described [9]. Qubits have to store information, be able to be manipulated, and the data must be transferred effectively (quantum bus). Maintaining the quantum coherence of qubits is a critical challenge. The system must be insensitive to external disturbances [10].

The central idea for multiple qubit architectures is based on coupling them to a transmission line and selecting from the frequency at which qubits couple together [11,12]. These behave as non-linear resonators in which the fundamental and excited levels can be isolated. The general scheme of a coplanar microwave transmission line consists of a one-dimensional cable evaporated on an insulating substrate with two superconducting planes parallel to it on the same surface [13,14]. A finite transmission line is obtained by interrupting the central wire at two points to get a resonator that supports a quantized electromagnetic mode, behaving as a harmonic [15].

Cavity Quantum Electrodynamics (cavity QED) studies an atom as a TLS and the quantized electromagnetic field within an optical cavity [16]. Inspired by the circuit QED, Circuit Quantum Electrodynamics (circuit QED) emerges, in which superconducting qubits and microwave transmission lines called resonators replace the atoms and the cavity, respectively [13,17]. A TLS and a coupled harmonic oscillator exist in the cavity QED and the circuit QED [16,17,18].

Superconducting resonators are crucial elements in a circuit QED [17]. A microwave resonator can be an excellent quantum bus for qubit manipulation and memory [11]. In the cavity QED, the atom couples to the electromagnetic field through its electric dipole moment, while superconducting qubits have several coupling alternatives [19]. Circuit QED is a well-suited platform for quantum data due to its flexibility, scalability, and tunability [20]. Experimental works demonstrated superconducting qubits and a microwave resonator strong coupling [21].

Developing a molecular circuit engineering that controls the haphazard states of individual molecular events is a big challenge. In this work, an electrical representation of DNA is carried out following the Quantum Electrodynamics Theory of Circuits and considering the results previously published by our team [8]. A gene is described as a problem of a qubit capacitively coupled to two transmission guides subject to a harmonic forcing in the magnetic flux that passes through it for a frequency ${\omega }_r$ and an amplitude $A$. A-T, T-A, C-G, and G-C bp are structurally similar to Radio Frequency Superconducting Quantum Interference Device (RF-SQUID), where each level is coupled to a different transmission line. In a resonant regime, the separation between the two energy levels of the qubit is in resonance with the frequency of the transmission lines. Entangled qubit-photon states are generated, and both system parts could exchange excitations. We consider genes as a combination of n single-mode microwave resonators and cavities coupling in parallel. The Frequency matching can selectively obtain coupling between the phosphorylation pulse emitted by RNA pol II and any particular SQUID bp during transcription.

2. DNA and Transcription: Structural and Biochemical Studies

Genomic DNA is double-stranded. The H-bonded bases are assembled in parallel, perpendicular to the helix backbone. The x-ray analysis revealed a regular space of 0.34 nm between the bases and 36° of twist angle about the helical axis. Ten pairs of bases define the helix anatomy of a complete turn every 3.4 nm [22,23]. A repeating sugar (deoxyribose)- phosphate units form the strand backbone [23]. Each bp triplet constitutes a codon that carries on the genetic code. However, most of the human genome does not code for proteins. These regions may provide RNA molecules and transcription factors' specific sequences of binding sites and exert regulatory functions [24]. The minor and major grooves are two helical spaces with different widths. The first one is delimited by the space when the two antiparallel DNA strands run closest together, while the major groove is defined when they are furthest apart [22]. The above information is the DNA B schematic representation, the predominant type in most cell DNA stretches [23].

The eukaryotic transcription cycle arises during three phases. The RNA chain initiation undergoes various steps, recruiting multiple cellular factors, the mediator complex, and the RNA pol II to form a pre-initiation complex (PIC) at core promoter sequences. PIC composition is identical at all yeast promoters and is the rate-limiting step for transcription [25]. After, the RNAP II accesses the DNA template strand to begin the RNA synthesis. All cellular RNA pol II forms a stable elongation complex (EC) with the DNA template and RNA transcript for transcription elongation [26].

Eukaryotic polymerases comprise an evolutionarily conserved catalytically competent ten-subunit core and up to seven peripheral subunits. The C-terminal domain (CTD) of the RNA pol II largest subunit has a repeated heptad motif (Tyr₁Ser₂Pro₃Thr₄Ser₅Pro₆Ser₇) [27]. The copy number depends on the organism's complexity (26 in yeast vs. 52 in humans). CTD's dynamic phosphorylation (P) and dephosphorylation coordinate the specific association and release of transcriptional factors at every stage of transcription (Figure 1A) [28,29]. For example, transcription initiation and promoter escape at the 5′-end are facilitated by Ser₅-P for Kin28/Cdk7, while for the elongation Ser₅-P is by Srb10/Cdk8 and collaborates Ser₂-P by Ctk1/Cdk9. The termination and the 3′-end processing are because of the Ser₂-Ser₂-P transition by Bur1/Cdk9 [30,31]. A similar P-process regulates the biological activity of many proteins.

In most promoters, the two main recognition consensus sequences TTGACA and TATAAT (−35 and −10 hexameric boxes) are separated by a region of 17 ±1 bp called a spacer [32,33]. The −35 element interacts with the σ4 domain and the DNA major groove by inserting a helix-turn-helix (HTH). The −43 to −38 sequences interact with the DNA minor grove and the α-CTD [34]. During elongation, the transcription bubble has a 12±1 nt constant size and 8±1 bp RNA–DNA hybrid until RNA pol II reaches a termination codon [35]. RNA pol II selects the correct nt, adds it to the growing RNA chain 3′-end, and releases pyrophosphate (PP_i). The same residues contacting the nt exist in all RNA pol, consistent with a general nt addition cycle (NAC) [36]. The NTP substrate attaches transiently to the open active center enzyme conformation. Structural information on an RNA pol II EC revealed a nine-base pair DNA-RNA hybrid duplex from the active center cleft at the floor [2,37,38].

In the transcription bubble, DNA is unwound before the RNA pol II catalytic site, allowing the template single strand to reach it and form the DNA-RNA hybrid duplex. It is then rewound to assemble the exiting DNA duplex. Within the unwound region, the growing RNA is attached to the active center cleft with its 3′-end and forms the nine bp hybrid duplex with the DNA template strand [26]. The high EC stability depends on the DNA-RNA tight binding and facilitates transcription processivity. Again, CTD-P plays an important role. Spt₅ spans the pol II cleft from the clamp, and contacts upstream DNA, and the non-template DNA strand [39]. Spt₅ positive charges interact with the negatively charged DNA non-template strand to prevent the bubble's collapse. The Spt₅ NGN domain links to all polymerases clamp in the coiled-coil [36], indicating an evolutionary conserved binding mode.

3. A-T, T-A, G-C, AND C-G Base Pairs: RF Squids—The Hamiltonian Scheme

Superconducting qubits steam as a leading technology platform for quantum computing [11,40]. Flux qubits exploit magnetic flux through a superconducting loop. A superconducting ring is interrupted by one or more JJs and crossed by a magnetic flux based on the Josephson Effect [11]. The qubit's quantum state is encoded in the persistent loop current and manipulated by altering the magnetic flux. Suppose an external oscillating flow of fixed frequency and variable amplitude is applied. In that case, the qubit can be forced by generating oscillations in the populations of its two energy levels through the Landau-Zener-Stuckelberg Transitions.

Phase qubits store quantum information in the difference of phases between the JJ on both sides. The quantum states coherent rotations can be induced by applying external microwave pulses [41]. One of the most used systems in Circuit QED is the charge qubit (Cooper pair box) coupled through the electrical component of the field to a superconducting coplanar waveguide resonator [40]. Recently, many studies have highlighted the central role of transmon qubits (Figure 2A) in quantum computing [19]. In them, the superconducting island is controlled by an external magnetic field that induces a magnetic flux across the JJs, affecting the superconducting phase difference [11,13].

An anharmonicity with a high energy difference between consecutive levels ensures well-defined qubits and long-lasting quantum calculations. In a previous study, our working group established the role of DNA as a perfect quantum computer by characterizing Adenine (A), Thymine (T), Cytosine (C), and Guanine (G) nitrogenous bases as superconducting islands. In addition, we defined DNA qubits by combining the system quantum state with classical information:

$$|0⟩\left[AT\right]=|0⟩\left[TA\right]=|0⟩\left[X\right],\mathrm{and}|1⟩\left[GC\right]=|1⟩\left[CG\right]=|1⟩\left[Y\right]$$

A-T and C-G are superconducting qubits with intrinsically quantum behavior. They have quantized energy levels, so their separation is more significant than thermal fluctuations and can be in superposition states. The JJ introduces anharmonicity to the system [11], thus obtaining non-equispaced energy levels. Each transition has a distinct frequency because of the decreasing energy spacing as the quantum number progressively increases. If externally applied a magnetic flux, the double-well potential and the two eigenfunctions become symmetrical while antisymmetrical superpositions of the two basis states. In these superconducting rings, a single total wave function describes the state of all electron and hole pairs in the condensate: ${\psi }_{\overrightarrow{p}}=\psi e^{i\overrightarrow{p}.\overrightarrow{r}∕\hslash }$.

$\overrightarrow{p}$ single electron and hole pair momentum.

${\left|\psi \right|}^2$is the probability density of finding an electron and hole pair in a given volume

$\overrightarrow{p}.\overrightarrow{r}∕\hslash$ of the phase associated with each electron and hole pair in space.

Because this wave function must be unvalued, the accumulated phase change ($∆\varnothing$) when traveling a closed path around the year must be a multiple of $2\pi$:

$∆\varnothing =2n\pi$ where $n$ is an integer number

Because the RNA pol II contacts only one element of the pair (classical information) in the template strand, we considered the number of aromatic rings to get one of the four Bell states:

Purines $=|0⟩$, Pyrimidines $=|1⟩$, then

$|\psi ⟩A=|\psi ⟩G=|0⟩$, and $|\psi ⟩T=|\psi ⟩C=|1⟩$, then

$|0⟩\left[X\right]=|0⟩\left[01\right]$ or $|0⟩\left[10\right]$, and $|1⟩\left[Y\right]=|1⟩\left[01\right]$or $|1⟩\left[10\right]$

$|\psi ⟩AT=\left[00\right]$, $|\psi ⟩TA=\left[01\right]$, $|\psi ⟩CG=\left[11\right]$, and $|\psi ⟩GC=\left[10\right]$

A JJ forms the foundation of a SQUID. However, these junctions also present an internal capacitive effect because they are constituted by two superconducting electrodes separated by an insulator [11]. The JJ can be considered a linear harmonic LC circuit with nonlinear inductance. Two types of SQUIDs are commonly used in quantum computer applications: a DC SQUID, which contains two parallel junctions, and an RF SQUID, which includes just one junction.

In DNA qubits, the central H-bond, acting like a JJ, forms the foundation of the RF SQUIDs (Figure 2B). The H-bond can induce strong inter- or intra-molecular electronic coupling by enhancing resonance, electron delocalization, or planarizing the conjugated backbone.

Nitrogenous (N) number three (N3) of a T in T-A bp has a lone pair of electrons as part of the $\pi$-cloud. The N1 in A uses the lone pair to attract the H attached to the N3 of T. Similarly, in C-G bp, the N3 in C uses its lone pair of electrons to form an H-bond with the N1 of G. Then, in our previous work, we explained, using physics approximations, how electrons interchange a biological quantum of energy to form electrons and hole pairs in the JJ. The electron ($e$) three ($e_3$) with the momentum $P+K$ emits a Biology boson ($B_b=\hslash \omega$) to $e_2$ losing the momentum $K$ and occupies the position of the hole ($h$) one ($h_1$) with momentum $P$. Then, $e_1,$ through the Tunnel Effect, passes through the barrier (NH--N) to the position of $e_3$, where there is now a hole. Please check Figure 4 in the reference [8].

The first Josephson Equation $I_s=I_0.sinsin\left({\varnothing }_2-{\varnothing }_1\right)$describes how the tunneling current depends on the two superconductors' phase difference. The second Josephson equation ${\displaystyle \frac{d}{dt}}\left({\varnothing }_2-{\varnothing }_1\right)={\displaystyle \frac{qV}{\hslash }}$ represents the time evolution of the phase difference given an external voltage [11]. $V$ is the potential difference in JJ. These two equations describe the temporal phase evolution of the two superconductors. If a constant voltage is applied to a JJ, the phase difference will evolve linearly in time (AC current). As opposed, if we use no external voltage, the phase difference becomes constant (constant supercurrent despite no external voltage) due to the phase coherence of the pairs. It is known as the DC Josephson Effect.

If $\varnothing ={\varnothing }_2-{\varnothing }_1$, then

$I_s=I_0.sinsin\left(\mathrm{\varnothing }\right)$Y $V={\displaystyle \frac{h}{4\pi e}}{\displaystyle \frac{d\mathrm{\varnothing }}{dt}}$

$P={\displaystyle \frac{d{E}_{jj}}{dt}}$ with $E_{jj}=\int Pdt=\int V{I}_{s}dt=\int {\displaystyle \frac{h}{4\pi e}}{\displaystyle \frac{d\mathrm{\varnothing }}{dt}}{I}_0.sinsin\left(\mathrm{\varnothing }\right)dt$

Then, $E_{jj}={\displaystyle \frac{h}{4\pi e}}$ $I_{0}coscos\mathrm{\varnothing }$, where $E_j={\displaystyle \frac{h}{4\pi e}}$ $I_0$

Then, $coscos\mathrm{\varnothing }=1-{\displaystyle \frac{1}{2}}{\mathrm{\varnothing }}^2-{\displaystyle \frac{1}{4!}}{\mathrm{\varnothing }}^4+O\left({\mathrm{\varnothing }}^6\right)$The first term $E_j$ can be disregarded because it is a constant factor while the terms of order $O\left({\varnothing }^6\right)$ can be disregarded because of being $\varnothing$ small.

The Hamiltonian is an energy operator containing kinetic energy $K={\displaystyle \frac{P^2}{2m}}$, and a potential energy of the form $V\left(r\right)$.

Them, $H={\displaystyle \frac{P^2}{2m}}+V\left(r\right)$The Hamiltonian of a qubit with a single JJ is given by

$H={\displaystyle \frac{Q^2}{2C}}+U\left(\Phi \right)$ with $U\left(\Phi \right)={\displaystyle \frac{{\Phi }^2}{2L}}-E_{j}coscos\left[{\displaystyle \frac{2\pi }{{\Phi }_0}}\left(\Phi -{\Phi }_{ext}\right)\right]$

$$\mathrm{\varnothing }=2\pi {\displaystyle \frac{\Phi }{{\Phi }_0}}$$

$n={\displaystyle \frac{Q}{2e}}$ and $E_C={\displaystyle \frac{e^2}{2C}}$

$$H={\displaystyle \frac{Q^2}{2C}}+{\displaystyle \frac{{\Phi }^2}{2L}}-E_{j}coscos\left[{\displaystyle \frac{2\pi }{{\Phi }_0}}\left(\Phi -{\Phi }_{ext}\right)\right]$$

$$H=4E_C{n}^2+E_j\left({\displaystyle \frac{1}{2L}}+{\displaystyle \frac{4\pi }{{{\Phi }_0}^2}}\right){\Phi }^2+E_j{\displaystyle \frac{2}{3{{\Phi }_0}^4}}{\Phi }^4$$

(1)

$\Phi$: flow in the JJ

${\Phi }_{ext}$: external magnetic flux, ${\Phi }_{ext}=0$

$C$: capacitance

$L$: Inductance

$E_j$: Josephson energy

${\Phi }_0$: magnetic flux quantum

$E_{jj}$: JJ energy

Equation 1 is the Hamilton equation. Following the first quantization procedure, the two conjugate variables $n$ and $\Phi$ are promoted to non-commuting operators

$n=\widehat{n}$ y $\Phi =\widehat{\Phi }$

$$\widehat{H}=4E_C{\widehat{n}}^2+E_j\left({\displaystyle \frac{1}{2L}}+{\displaystyle \frac{4\pi }{{{\Phi }_0}^2}}\right){\widehat{\Phi }}^2+E_j{\displaystyle \frac{2}{3{{\Phi }_0}^4}}{\widehat{\Phi }}^4$$

(2)

The first two terms in Equation 2 represent harmonic oscillators, and the third is the anharmonicity.

Introducing the standard annihilation ${\widehat{a}}^{-}$ and creation ${\widehat{a}}^{+}$ operators is useful.

$${\widehat{a}}^{+}={\displaystyle \frac{1}{\sqrt{2}}}\left({\displaystyle \frac{1}{\sqrt{\beta }}}\widehat{\Phi }+i\sqrt{\beta }\widehat{n}\right)$$

$${\widehat{a}}^{-}={\displaystyle \frac{1}{\sqrt{2}}}\left({\displaystyle \frac{1}{\sqrt{\beta }}}\widehat{\Phi }-i\sqrt{\beta }\widehat{n}\right)$$

With $\beta =\sqrt{8{\displaystyle \frac{E_C}{E_{j^{*}}}}}$ where $E_{j^{*}}={\displaystyle \frac{1}{2}}\left({\displaystyle \frac{1}{L}}+{\displaystyle \frac{8\pi }{{{\Phi }_0}^2}}\right)E_j$The action of ${\widehat{a}}^{+}$ is to create a quantized excitation of the flux and charge degrees of freedom of the magnetic and electric fields. It creates a photon of frequency ${\omega }_r$ stored in the circuit.

In function of ${\widehat{a}}^{+}$and ${\widehat{a}}^{-}$:

$\widehat{\Phi }=\sqrt{{\displaystyle \frac{\beta }{2}}}\left({\widehat{a}}^{-}+{\widehat{a}}^{+}\right)$ and $\widehat{n}={\displaystyle \frac{i}{\sqrt{2\beta }}}\left({\widehat{a}}^{-}-{\widehat{a}}^{+}\right)$

Then, the Hamiltonian in Equation 2 combines in:

$$\widehat{H}={\omega }_r\left({\widehat{a}}^{+}{\widehat{a}}^{-}+{\displaystyle \frac{1}{2}}\right)+{\displaystyle \frac{\alpha }{12}}{\left({\widehat{a}}^{+}+{\widehat{a}}^{-}\right)}^4$$

(3)

With ${\omega }_r=\sqrt{8E_C{E}_{j^{*}}}$ the new frequency and the anharmonicity $\alpha =-E_C$

Please note that the Hamiltonian represented in Equations 2 and 3 corresponds to our qubit proposal based on the DNA characteristics and differs from those described in other works (Figure 2B).

4. DNA Backbone as a Finite Transmission Line with Distributed Parameters—The Hamiltonian Scheme

A transmission line is a structure of uniform geometry used to transport electrical signals or RF energy efficiently from one point to another [11]. It is also a medium that propagates information through electromagnetic waves at very high frequencies [42]. Because it behaves like a cable that spreads energy, it can be analyzed like a quantum oscillator electrical circuit based on its frequency response. The circuit is considered a line with distributed parameters at high frequencies since it has dimensions comparable to the wavelength [43]. Thus, the circulated current's amplitude and phase differ at every point. The wavelength is very short for high frequencies, so the transmission lines behave like resonant circuits for a specific frequency range.

Two parallel conductors form pair lines. Its various variants are used in telephony and data transmission [44,45]. Electromagnetic energy propagates around the conductors of the line. If a voltage $Vs$ is applied to a pair line, an electric ($E$) field is generated between the conductors as opposite charges accumulate. The $Vs$ makes an electric current $\left({\rm I}\right)$ flow through the line's conductors terminated at a load impedance ($Z_L$). This current also generates a magnetic $H$ field around the conductors. The direction of the current and fields is reversed in each half-cycle of the voltage [45]. By multiplying $H$ by $\mu$, we obtain the induction magnetic vector ($B$), and then, by integrating $B$ over a plane parallel to the wires, we get the magnetic flux ($\Phi$) that “joins” the circuit.

The mathematical definition of a circuit inductance is $L={\displaystyle \frac{\Phi }{{\rm I}}}$. If we have a transmission line length $x_0$ units long, and if that line has a distributed inductance $L$, then the inductance is just $L=L{x}_0$. Two conductors with charge $\pm Q$ separated by a distance generate a potential difference ($V$). $Q=\pm \left(\rho x_0\right)$ in every section of a transmission line length $x_0$. Where $\rho$ is the lineal distribute charge. The basic transmission line parameters are its capacitance (C) and inductance (L) per unit length, the impedance ($Z_0$), and the mode frequency ($\omega$) [21]. These parameters are related by:

$$Z_0=\sqrt{{\displaystyle \frac{L}{C}}}\mathrm{and}\omega ={\displaystyle \frac{1}{\sqrt{LC}}}$$

The electric and magnetic fields are transverse to the direction of propagation. It is called Transverse Electromagnetic Waves (TEM) [46]. We can distinguish between infinite, semi-finite, and finite transmission lines. The latter is obtained by interrupting the central cable at two points so that a resonator that supports a quantized electromagnetic mode behaving like a harmonic oscillator is obtained. A resonator is a structure that encloses oscillatory electromagnetic energy without spreading it [11]. Boundary conditions restrict the wave configuration inside to patterns with integral numbers of half or quarter wavelengths along a propagation axis [45]. Therefore, specific discrete resonant frequencies ${\omega }_n$ are allowed. There is no energy dissipation if the resistive element is in series at zero current or parallel at zero voltage. Typical TEM resonators contain lossless elements at their end. Resonators are widely used to manipulate signals and energies, acting as bandpass filters allowing only frequencies close to the desired resonant frequency (${\omega }_n$) or as bands-stop filters letting all frequencies pass [21].

Characterizing single-stranded DNA as a transmission line is also of great interest. Single-stranded DNA can exhibit rich electrical properties. In the DNA backbone (formed by phosphate groups and deoxyriboses), C-N bonds attached deoxyriboses to the bp [47]. The negative charges at the phosphate groups connect the base pairs and give structural support for the double helix. Small polarons assisted hopping between the next nearest neighbors of the DNA molecular wire was suggested as the transport mechanism responsible for the solid high-temperature dependence of the electrical conductivity [48]. Our DNA wire model has connections between successive base pairs, the only sites included. Porath and Cols demonstrated that the sugar-phosphate backbones nor the specific bp sequence mediate the long-distance conduction in DNA [49]. Besides the multiple efforts to study the electrical current through single DNA fragments, the conduction mechanism, especially in long DNA molecules, is far from being established.

Figure 2C shows a schematic representation of DNA based on the model of parameters distributed in a transmission line. The line length is proportional to the total capacitance and inductance. In this model, the backbone phosphate groups carry the inductive effect, while the bases have the capacitive effect. Since it is a superconducting circuit, we will ignore the electrical resistance and magnetic fields within the strands, which are expelled due to the Meissner effect [8,50].

Using the Telegrapher equations, we can obtain the transmission line equation and assume the DNA strand acts as a lossless line with $R=0$ and $G=0$, where R and G are the resistance and shunt admittance per unit length [51]. The first equation

$${\displaystyle \frac{\partial V\left(x,t\right)}{\partial x}}=-\left(L{\displaystyle \frac{\partial {\rm I}\left(x,t\right)}{\partial t}}\right)$$

explains the voltage dependence the distributed inductance $L$ multiplied by the time derived from the current flowing in the line at one point. The second equation

$${\displaystyle \frac{\partial {\rm I}\left(x,t\right)}{\partial x}}=-\left(C{\displaystyle \frac{\partial V\left(x,t\right)}{\partial t}}\right)$$

explains that the current loss as we go down the line is proportional to the distributed capacitance $C$ multiplied by the time rate of the voltage change in the line.

An LC oscillator is characterized by its inductance L, capacitance C, angular frequency ${\omega }_r={\displaystyle \frac{1}{\sqrt{LC}}}$, and impedance $Z_r=\sqrt{{\displaystyle \frac{L}{C}}}$. The oscillator energy is given by:

$$H_{LC}={\displaystyle \frac{Q^2}{2C}}+{\displaystyle \frac{{\Phi }^2}{2L}}$$

(4)

As $I={\displaystyle \frac{dQ}{dt}}$, and $Q\left(t\right)={\int }_0^{t}I\left(t^{\star }\right)d\left(t^{\star }\right)$. Then,

$$H_{LC}={\displaystyle \frac{Q^2}{2C}}+{\displaystyle \frac{1}{2}}C{W_r}^2{\Phi }^2$$

(5)

Where by Faraday's law:

$$\Phi \left(t\right)={\int }_0^{t}V\left(t\right)d\left(t\right)$$

Equation 5 is an LC harmonic oscillator equivalent to the coordinate mechanical harmonic oscillator $\Phi$

$$P=Q\mathrm{and}M=C$$

Then, $\left[\Phi ,Q\right]=i\mathrm{\hslash }$

Introducing the creation and annihilation operators (${\widehat{a}}^{+}$and ${\widehat{a}}^{-}$)

$\widehat{\Phi }={\Phi }_{zpf}\left({\widehat{a}}^{+}+{\widehat{a}}^{-}\right)$and $\widehat{Q}={iQ}_{zpf}\left({\widehat{a}}^{+}-{\widehat{a}}^{-}\right)$ with

${\Phi }_{zpf}=\sqrt{{\displaystyle \frac{\hslash }{2W_{r}C}}}=\sqrt{{\displaystyle \frac{\hslash Z_r}{2}}}$ and $Q_{zpf}=\sqrt{{\displaystyle \frac{\hslash W_{r}C}{2}}}=\sqrt{{\displaystyle \frac{\hslash }{2Z_r}}}$. Then,

$${\widehat{H}}_{LC}=\hslash W_r\left({\widehat{a}}^{+}{\widehat{a}}^{-}+{\displaystyle \frac{1}{2}}\right)$$

Where ${\widehat{a}}^{+}{\widehat{a}}^{-}|n⟩=n|n⟩$ with $n=\mathrm{0,1},2,\dots$and ${\displaystyle \frac{1}{2}}$ is the zero-order energy or null energy.

${\widehat{a}}^{+}=\sqrt{{\displaystyle \frac{\hslash Z_r}{2}}}\left(\widehat{\Phi }-i{Z}_r\widehat{Q}\right)$ creates a quantum excitation of flow and charge equivalent to the electric and magnetic Fields (a frequency photon ${\omega }_r$ inside the circuit).

We can generalize the result for a transmission line in the x direction and length d. The energy associated with each capacitor is ${\displaystyle \frac{{Q_n}^2}{2C_0}}$, and to each inductor ${\displaystyle \frac{{\left({\Phi }_{n+1}-{\Phi }_n\right)}^2}{2L_0}}$. Where ${\Phi }_n$ is the variable flow associated with node n and $Q_n$ the conjugate variable is the charge at node n.

Then, the Hamiltonian for the discrete linear transmission resonator associated with Figure 2C is

$$H=\sum _{n=0}^{N-1}\left[{\displaystyle \frac{1}{2C_0}}{Q_n}^2+{\displaystyle \frac{1}{2L_0}}{\left({\Phi }_{n+1}-{\Phi }_n\right)}^2\right]$$

$L_0$and $C_0$: the variable inductance and capacitance associated with each node n and flow ${\Phi }_n$.

In the continuous limit, considering ${\delta }_x$ the size of a unit cell that tends to zero:

$C_0={\delta }_x{c}_0$ and $L_0={\delta }_x{l}_0$ with $c_0$ y $l_0$ the linear density of capacitance and inductance.

Defining a continuous flow field as $\Phi \left(x_n\right)={\Phi }_n$ and the charge $Q\left(x_n\right)={\displaystyle \frac{Q_n}{{\delta }_x}}$. Taking the ${\delta }_x$and $d=N∆X$ constant, the Hamiltonian for a continuous linear transmission resonator:

$$H={\int }_0^{d}dx\left[{\displaystyle \frac{1}{2c_0}}{Q}^2\left(x\right)+{\displaystyle \frac{1}{2l_0}}{\left({\partial }_x{\Phi }_x\right)}^2\right]$$

(6)

Where ${\partial }_x{\Phi }_x={\displaystyle \frac{\left({\Phi }_{n+1}-{\Phi }_n\right)}{{\delta }_x}}$

The charge $Q\left(x,t\right)=c_0{\partial }_t\Phi \left(x,t\right)$ is the generalized canonical moment of the flow $\Phi \left(x,t\right)={\int }_0^{t}V\left(x,t\right)dt$ with $V\left(x,t\right)$ the voltage.

Using the information above, we obtain the wave propagation equation along the linear transmission resonator

$${v_0}^2{\displaystyle \frac{{\partial }^2\Phi \left(x,t\right)}{\partial x^2}}-{\displaystyle \frac{{\partial }^2\Phi \left(x,t\right)}{\partial t^2}}=0$$

(7)

${v_0}^2={\displaystyle \frac{1}{l_0{c}_0}}$ It is the speed of light in the biological medium of DNA.

The Equation 7 solution is expressed as a linear combination of normal nodes of a size ${\displaystyle \frac{\lambda }{2}}$.

$$\Phi \left(x,t\right)=\sum _{m=0}^{\mathrm{\infty }}{U}_m\left(x\right){\Phi }_m\left(t\right)$$

with ${\ddot{\Phi }}_m=-{{\omega }_m}^2{\Phi }_m$ a differential equation of an oscillating function in time at the frequency node ${\omega }_m$

$$U_m\left(x\right)=A_{m}coscos\left(K_{m}x+{\phi }_m\right)$$

with amplitude $A_m$, wave vector $K_m={\displaystyle \frac{{\omega }_m}{v_0}}$, and phase ${\phi }_m$.

The Hamiltonian given by Equation 6 can be expressed as a discrete Hamiltonian:

$$H=\sum _{m=0}^{\mathrm{\infty }}\left[{\displaystyle \frac{{Q_m}^2}{2C_r}}+{\displaystyle \frac{1}{2}}{C}_r{{\omega }_m}^2{{\Phi }_m}^2\right]$$

(8)

This Hamiltonian (Equation 8) represents the sum of simple harmonic oscillators. With $C_r=d{c}_0$ as the total capacitance of the resonator and $Q_m=C_r{\dot{\Phi }}_m$ the conjugate charge ${\Phi }_m$.

Following the quantization procedure with the variables ${\Phi }_m$ and $Q_m$, we define the creation and annihilation operators:

$${\widehat{\Phi }}_m=\sqrt{{\displaystyle \frac{\hslash Z_m}{2}}}\left({{\widehat{a}}_m}^{+}+{{\widehat{a}}_m}^{-}\right)$$

$${\widehat{Q}}_m=i\sqrt{{\displaystyle \frac{\hslash }{2Z_m}}}\left({{\widehat{a}}_m}^{+}-{{\widehat{a}}_m}^{-}\right)$$

Where $Z_m=\sqrt{{\displaystyle \frac{L_m}{C_r}}}$ the characteristic impedance of the m node and ${\displaystyle \frac{1}{L_m}}=C_r{{\omega }_m}^2$.

Then, the Hamiltonian given in Equation 8 finally transforms into:

$$\widehat{H}=\sum _{m=0}^{\infty }\hslash {\omega }_m{{\widehat{a}}_m}^{+}{{\widehat{a}}_m}^{-}$$

With ${\omega }_m=\left(m+1\right){\omega }_0$, and ${\displaystyle \frac{{\omega }_0}{2\pi }}={\displaystyle \frac{v_0}{2d}}$ is the frequency.

5. DNA as a Quantum Electrodynamic Circuit—The Hamiltonian Scheme

A coplanar waveguide's (CPW) geometry allows for increased electromagnetic energy density and coupling [20]. The strong coupling permits operations to be carried out on the qubit and even entanglement of multiple qubits in the same cavity. Using two-dimensional CPWs (TD-CPWs) and circuit QED, researchers have realized high-efficiency entanglement up to 20 qubits and high-fidelity two-qubit gates for achieving up to 64 qubits. The electromagnetic field inside conventional superconducting resonators in circuit QED, (including TD-CPWs and 3D-cavities), are standing waves determined by boundary conditions, resulting in a topographic dependence of the fields’ amplitude [21]. In circuit QED, resonant microwave cavities are replaced by sections of microwave transmission lines (integrated circuits) [19]. If superconducting materials are used, the system's losses and dephasing would be limited only by the system properties [11]. Another way to mitigate losses is to lower the electric field energy stored at interfaces and surfaces to the energy stored in a vacuum. That is why using 3D-microwave cavities or 3D- resonators rather than planar circuits is preferred [17,52]. This is a new quantum information processing prototype with information stored in a cavity.

Here, the effective Hamiltonian for our DNA archetype can be represented as

$$H=H_{RF-S}+H_{R1}+H_{R2}+H_{I_{RFS-R_1-R_2}}+H_{\kappa }+H_{\gamma }$$

Where, $H_{RF-S}$ is the Hamiltonian of a flux qubit with a single JJ (A-T, T-A, G-C, or C-G). $H_{R1}$ and $H_{R2}$ represented the Hamiltonian of the two single-stranded DNA separately. To complete the effective representation of the Hamiltonian in DNA, $H_{I_{RFS-R_1-R_2}}$represents the interaction between the qubit and two-stranded DNA.

$H_{\kappa }$, the cavity losses, and the decay rate.

$H_{\gamma }$, decoherence, or decay of the two-level system.

As A-T and C-G are superconductors, their properties would limit the system's losses and dephasing ($H_{\kappa }=H_{\gamma }=0)$.

According to the Hamiltonian, a strong coupling qubit-cavity produces a special footprint when entry into resonance. The coherent oscillation between two quantum states is known as Rabi oscillation and the frequency as the Rabi frequency ${\Omega }_R$ (Figure 3). The ${\Omega }_R$ is proportional to the energy of the TLS-electromagnetic field interaction [53]. For a flux qubit coupled to a single mode of the electromagnetic field of a transmission line that acts as a resonator, the system's Hamiltonian consists of three parts: the qubit term ($H_1$), the electromagnetic field term confined in the resonator ($H_2$), and the Rabi-type interaction term ($H_3$).

$$\widehat{H}={\widehat{H}}_1+{\widehat{H}}_2+{\widehat{H}}_3={\displaystyle \frac{1}{2}}\left(ϵ\left(t\right){\sigma }_z+\Delta {\sigma }_x\right)+\hslash {\omega }_r\left({\widehat{b}}^{+}{\widehat{b}}^{-}+{\displaystyle \frac{1}{2}}\right)+\hslash g{\sigma }_y\left({\widehat{b}}^{+}+{\widehat{b}}^{-}\right)$$

(9)

${\omega }_r$ resonator resonance frequency.

$ϵ\left(t\right)={ϵ}_0+Acoscos\omega t$, but in DNA under the ground conditions $Acoscos\omega t=0$. At this point, we consider the oscillator at a fixed instant.

Figure 3. Model of DNA as a 3D nonperturbative cavity quantum electrodynamics defined by two resonators. In the zero-detuning case, a qubit is strongly coupled to a cavity and produces a unique signal when tuned into resonance. We hypothesize that the H-bond provides maximum coupling. RAHB: the resonance-assisted hydrogen bond. The Jaynes-Cummings Hamiltonian assumes the rotating wave approximation ($g\ll {\omega }_r,\Omega$, ${\omega }_r~\Omega )$. The speed of quantum operations is represented by $g$. To maximize and maintain the coupling $g∕\pi$ larger than the two main dephasing rates of the system $(g\gg k,\gamma )$ is the target.

Figure 3. Model of DNA as a 3D nonperturbative cavity quantum electrodynamics defined by two resonators. In the zero-detuning case, a qubit is strongly coupled to a cavity and produces a unique signal when tuned into resonance. We hypothesize that the H-bond provides maximum coupling. RAHB: the resonance-assisted hydrogen bond. The Jaynes-Cummings Hamiltonian assumes the rotating wave approximation ($g\ll {\omega }_r,\Omega$, ${\omega }_r~\Omega )$. The speed of quantum operations is represented by $g$. To maximize and maintain the coupling $g∕\pi$ larger than the two main dephasing rates of the system $(g\gg k,\gamma )$ is the target.

Where ${ϵ}_0$corresponds to the minimum energy or the fundamental state (${ϵ}_0={\displaystyle \frac{1}{2}}\hslash \omega$)

$\Delta$ corresponds to the two-level difference in energy for $ϵ\left(t\right)={ϵ}_0$, $\Delta =\hslash {\Omega }_R$

$$\widehat{H}={\displaystyle \frac{\hslash {\Omega }_R}{2}}{\sigma }_x+\hslash {\omega }_r\left({\widehat{b}}^{+}{\widehat{b}}^{-}+{\displaystyle \frac{1}{2}}\right)+\hslash g{\sigma }_y\left({\widehat{b}}^{+}+{\widehat{b}}^{-}\right)$$

In the absence of damping, the Hamiltonian can be readily diagonalized. The detuning parameter is $ϵ\left(t\right)={\Omega }_R-{\omega }_r$. In the zero-detuning case ${\Omega }_R={\omega }_r$.

${\sigma }_x$, ${\sigma }_y$ and ${\sigma }_z$ are the Pauli matrices.

${\widehat{b}}^{+}$ and ${\widehat{b}}^{-}$photon creation and annihilation operators.

$g$ Intensity of the qubit-electromagnetic field coupling rate. It depends on the strength of the electric or magnetic field at the system site, its dipole moment, and the system’s properties [52].

$$g={\displaystyle \frac{d{A}_e}{\sqrt{2}\hslash }}\mathrm{or}g={\displaystyle \frac{\mu A_m}{\sqrt{2}\hslash }}$$

(10)

where $A_e$ and $A_m$ are the electric and magnetic field amplitudes.

Expanding the interaction term in terms of the rise and fall operators ${\sigma }_y={\sigma }_{+}+{\sigma }_{-}$

$$\hslash g\left({\sigma }_{+}{\widehat{b}}^{+}+{\sigma }_{+}{\widehat{b}}^{-}+{\sigma }_{-}{\widehat{b}}^{+}+{\sigma }_{-}+{\widehat{b}}^{-}\right)$$

A TLS and a harmonic oscillator are coupled in a circuit QED [21]. The circuit can operate in two limits: on-resonant and off-resonant dispersive regimes. On resonance, the separation between the two energy levels of the qubit resonates with the transmission line frequency. Entangled qubit-photon states are generated, and both system parts can exchange excitations without energy loss [11]. In an ultra-strong coupling regime, the coupling intensity reaches values comparable to the resonator frequency, rendering the Jaynes-Cummings Model invalid. In the dispersive regime, the transmission line frequency differs from the qubit energy levels separation. In this situation, the energy levels of the qubit will depend on the photons and vice versa. Consequently, the resonator could obtain information about the qubit's state and couple multiple qubits together [54].

Molecular circuits have been developed to control complex functions in biological and biochemical systems. In the DNA ground state, each nitrogenous base is coupled to a transmission line in resonance with the H-bond that joins it to its complementary base connected to a second transmission line. When RNA pol II emits pulses at the end of the transmission line, the system enters a dispersive regime. It is a sub-regime of the Jaynes-Cummings Model that allows measurements of the qubit's state acting on the resonator. The energy levels of the qubit are normalized depending on the state of the resonator and vice versa [54]. The qubit state is determined by varying the waveguide frequency until it matches the new frequency of the resonator. It also allows manipulation of the qubit by forcing the resonator with a signal sent through the waveguide similar to the natural qubit frequency.

This work approaches the DNA gene structure as the lineal combination of circuits-QED. Each qubit level of a flux qubit is coupled to an independent resonator. The frequency of the qubit is similar to that of the resonators it has been connected to. For a flux qubit coupled to two modes of the electromagnetic field of a transmission line that acts as a resonator, the system's Hamiltonian consists of three parts: the qubit term ($H_1$), the electromagnetic field term confined in the two resonators, and the Rabi-type interaction term ($H_2$). (modified from Equation 9):

$$\widehat{H}={\widehat{H}}_1+{\widehat{H}}_2={\displaystyle \frac{1}{2}}\left(ϵ\left(t\right){\sigma }_z+\Delta {\sigma }_x\right)+\sum _{i=1}^2\left(\hslash {\omega }_{r_i}\left({{\widehat{b}}_i}^{+}{{\widehat{b}}_i}^{-}+{\displaystyle \frac{1}{2}}\right)+\hslash g_i{\sigma }_y\left({{\widehat{b}}_i}^{+}+{{\widehat{b}}_i}^{-}\right)\right)$$

In the DNA ground state, $ϵ\left(t\right)={\Omega }_R-{\omega }_r=o$

$$\widehat{H}={\displaystyle \frac{\hslash {\Omega }_R}{2}}{\sigma }_x+\sum _{i=1}^2\left(\hslash {\omega }_{r_i}\left({{\widehat{b}}_i}^{+}{{\widehat{b}}_i}^{-}+{\displaystyle \frac{1}{2}}\right)+\hslash g_i{\sigma }_y\left({{\widehat{b}}_i}^{+}+{{\widehat{b}}_i}^{-}\right)\right)$$

In DNA, A-T, T-A, C-G, and G-C qubits are coupled to a virtual cavity defined by two resonators. This constitutes the DNA basic unit. The difference between the system’s energy levels in the DNA ground state is $\hslash {\Omega }_R$ due to the movement of electron and hole pairs in the aromatic ring without an external force [8]. A-T and C-G are strongly coupled to the resonators, which allow operations to be performed during the lifetime of the quantum state. The coherent oscillation between the two bp quantum states is equivalent to the phenomenon known as Rabi oscillations. The oscillation frequency is the Rabi frequency. It is related to the energy of the qubit-electromagnetic field interaction. A perturbation can tune the bp levels such that a certain level splitting coincides with the cavity photon energy Figure 4.

Analyzing Equation (10), the average field strength is estimated using the integral of the field amplitude over the space volume delimited by the two resonators. DNA spiral diameter is around 2 nm. Thus, we assume that the space dimensions give us $V={\displaystyle \frac{4}{3}}\pi r^3$, where $V$ is the sphere´s volume, and $r$ is the radio. Taking a space average and setting that $ϵ={ϵ}_0$ and $\mu ={\mu }_0$ [52]:

$${\displaystyle \frac{1}{2}}\int ϵ{\left|A_e\right|}^{2}dV={\displaystyle \frac{\hslash {\omega }_r}{2}}\mathrm{or}{\displaystyle \frac{1}{2}}\int {\displaystyle \frac{{\left|A_m\right|}^2}{\mu }}dV={\displaystyle \frac{\hslash {\omega }_r}{2}}$$

$${\displaystyle \frac{1}{2}}{ϵ}_0{A}_e^{2}V={\displaystyle \frac{\mathrm{\hslash }{\omega }_r}{2}}\mathrm{or}{\displaystyle \frac{1}{2{\mu }_0}}{A}_m^{2}V={\displaystyle \frac{\mathrm{\hslash }{\omega }_r}{2}}$$

$${\displaystyle \frac{2}{3}}{ϵ}_0{A}_e^2\pi r^3={\displaystyle \frac{\mathrm{\hslash }{\omega }_r}{2}}\mathrm{or}{\displaystyle \frac{2}{3{\mu }_0}}{A}_m^2\pi r^3={\displaystyle \frac{\mathrm{\hslash }{\omega }_r}{2}}$$

$$A_e^2={\displaystyle \frac{3\mathrm{\hslash }{\omega }_r}{4\pi {ϵ}_0{r}^3}}\mathrm{or}{A}_m^2={\displaystyle \frac{3{\mu }_0\mathrm{\hslash }{\omega }_r}{4\pi r^3}}$$

Then, $A_e={\displaystyle \frac{1}{2r}}\sqrt{{\displaystyle \frac{3\hslash {\omega }_r}{\pi {ϵ}_{0}r}}}$ for the electric field and $A_m={\displaystyle \frac{1}{2r}}\sqrt{{\displaystyle \frac{3{\mu }_0\hslash {\omega }_r}{\pi r}}}$ for the magnetic fiel.

In the following section, we will describe the theoretical basis for DNA Transcription and compare this process to the dispersive detuning regime.

6. Quantum Description of Transcription

Despite the progress in elucidating the critical events during transcription, what needs to be addressed is that the process takes place strictly in a spatially defined manner. How the transcription spatial organization is achieved is still an unsolved question. A physiochemical phenomenon called biomolecular condensates is partly formed through phase separation to achieve spatiotemporal regulation of some biological processes [55,56,57,58]. The phase separation is crucial in regulating the transcriptional events in every step of a transcription cycle (Figure 1A) [59,60].

Multivalent interactions constitute driving forces underlying phase separation in biomolecular condensates. They include disordered regions and structured, repeated modular oligomerization domains [61]. Some molecules function as scaffolds to assemble phase-separated structures, whereas client components are recruited by directly binding to scaffolds [55]. In general, modular interaction domains in evolutionarily conserved protein segments create multivalency that facilitates the engagement of multiple binding complexes. The dynamic interacting networks initiate and promote phase separation [60]. In many proteins, the succession of structural components such as tyrosine, glycine, and serine have provided multivalency to undergo phase separation and enable the formation of biomolecular condensates [29]. Unlike the classical transcription describing the RNA pol II tracks along the relatively static DNA template strand to synthesize the nascent RNA, the phase separation model allows for the description of a quantum model of transcription. RNA pol II can be incorporated into transcriptional condensates [58,61]. Transcription bubble translocation along the DNA requires rotation of the RNA pol around the DNA axis. The linking number density has been calculated using the equation:

$${\omega }_{0}x=\theta +\mathrm{\varnothing }\mathrm{with}{\omega }_0=1.85{nm}^{-1}$$

Where ${\omega }_{0}x$ is the accumulative rotational angle, $x$ is the distance in $nm$, $\theta$ is the RNA pol rotation angle, and $\varnothing$ the DNA rotation at the RNA pol site. Based on DNA supercoiling-transcription interplay, experiments have shown the maximum RNA pol II translocation rate or velocity is ${\upsilon }_0=20nm.s^{-1}\approx \left(60bp\right).s^{-1}$. [26].

The ability of pol II as the scaffold for the PIC arises from its CTD, which is flexible and highly disordered. CTD and phase separation close relationship is vital in initiating efficient transcription. The role of CTD-P in activating co-transcriptional capping is well-established [31]. The favorable regulation of CTD in condensates is reversibly controlled depending on CTD-P [29]. It determines pol II positioning into different condensates as transcription progresses. For example, the transition between initiation and elongation is associated with Kin₂₈-mediated Ser₅-P within the heptapeptide repeats [27]. Hence, it raises an essential question of whether P is the transcription's driving pulse.

In P-reactions, the side chain hydroxyl groups on specific proteins' serine, threonine, and tyrosine residues are modified with the gamma phosphate, producing a theatrical change in the protein's architecture and function [62]. ATP is almost always the phosphate donor [63]. It can be represented as A-P~P~P, where °~°°~° are the high-energy acid anhydride bonds. One of those is hydrolyzed. Breaking these bonds is an endergonic process that consumes energy rather than releasing it [64,65]. The negative change in the Gibbs free energy is due to the bonds formed after the hydrolysis process having a lower energy than before. The new bonds in hydrolysis allow the release of enormous energy, exactly 7.7 kcal/mol: ΔG= -7,7 kcal/mol o -31 KJ/mol [66].

Unlike most biological processes in which only one phosphate group ($P_i$) is separated, the last two phosphate groups are released in the form of the pyrophosphate group (${PP}_i$) in replication. This process can be summarized in a chemical equation:

${\left(DNA\right)}_n+dNTP\leftrightarrow {\left(DNA\right)}_{n+1}+{PP}_i$ [67]

Similarly, during transcription elongation, the overall reaction equation:

${\left(NMP\right)}_n+NTP\leftrightarrow {\left(NMP\right)}_{n+1}+{PP}_i$ [68]

The H-bond is a robust fixed dipole-dipole electrostatic force when many molecules join. It provides excellent stability but is weaker than the covalent or ionic bonds. The energy of a H-bond is typically 5 to 30 kJ/mol, but it can vary from very weak to very strong. For N—H...N, it is approximately 13 kJ/mol. Intermolecular H-bond is responsible for the high boiling point of water (100°C). ATP to ADP hydrolysis releases 30.5 kJ/mol, with a free energy change of 3.4 kJ/mol [66].

A transition between energy levels caused by forcing one of the system parameters' Hamiltonian is called a Landau-Zener transition. A series of successive Landau-Zener transitions occur when a TLS is subject to periodic forcing on one of its parameters with a sufficiently intense amplitude [69]. Detuning can be modified externally

$$ϵ\left(t\right)={ϵ}_0+Acoscos{\omega }_{ext}t$$

where $A$ and ${\omega }_{ext}$ are the forcing amplitude and frequency, respectively. The qubit levels are anti-symmetrically coupled

$$g_1={-g}_2,\mathrm{and}{ϵ}_1\left(t\right)=-{ϵ}_2\left(t\right).$$

The qubit eigenstates can be defined as $|\uparrow ⟩$, $|\downarrow ⟩$.

The dispersive detuning regime, where $ϵ\left(t\right)\gg g$ is helpful in quantum computing because it allows quantum measurements of the qubit state and its manipulation [11]. Irradiating at the qubit frequency $\Omega$ handle the qubit coherently. Irradiating at the resonator frequency ${\omega }_r$ introduce photons in the cavity that become entangled with the qubit states. The law of conservation of linear momentum requires that at least two photons be created so that the resulting linear momentum can be equal to zero [46,70]. The photon is the gauge boson of the electromagnetic interaction and, therefore, is a boson that acts as a carrier of a fundamental interaction of nature. Conservation laws can determine the energies of the two photons or, equivalently, their frequencies. A driving radiation field of frequency ${\omega }_{ext}$ acting on the qubit-resonator system during a specific time can be described by the Hamiltonian [71]:

$$\widehat{H}\left(t\right)=\mathrm{\hslash }∆\left(t\right)\left({\widehat{b}}^{+}{q}^{-i{\omega }_{ext}t}+{\widehat{b}}^{-}{q}^{i{\omega }_{ext}t}\right)$$

(11)

where $q$is the electron charge.

Every qubit is represented as a point on the Bloch sphere, where the poles represent the eigenstates [72]. The colatitude angle, $\theta$, runs through the angles from 0 to $\pi$, while the azimuthal angle, $\varphi$, runs from 0 to 2$\pi$ [73]. The effect of a quantum gate on a single qubit is a rotation in the Bloch sphere, giving rise to a new qubit. In the Bloch sphere, we visualize the action of different logic gates or the temporal evolution of the state of a TLS described by a Hamiltonian. An axis and a rotation angle define a rotation operator [74]. Its action produces a rotation of the point on the sphere concerning the rotation axis at the rotation angle [73,75]. The relationship between the angle of rotation in the x-axis and the pulse parameters in the Bloch sphere is represented in the equation:

$$\theta \left(t\right)=\Omega V_0{\int }_0^{t}f\left(t^{´}\right)d{t}^{´}$$

Since the Rabi oscillations follow an expression of sin[2], we can obtain their period and determine the Rabi frequency and, even more importantly, the value of the variable parameter of the pulse makes it go from the state $|0⟩$ to the state $|1⟩$, a logic gate $\pi$ [75]. To perform rotations around the y-axis, multiply the pulse amplitude by the imaginary unit (a phase of $\varnothing ={\displaystyle \frac{\pi }{2}}$ that is added to the pulse).

$$U\left(t\right)=expexp\left(-i{\int }_0^t{H}_d^I\left(t^{´}\right)d{t}^{´}\right)=exp\left({\displaystyle \frac{i}{2}}\Omega V_0\left(I{\sigma }_x-Q{\sigma }_y\right){\int }_0^{t}f\left(t^{´}\right)d{t}^{´}\right)$$

By combining Equation 11 with the Hamiltonian of the system in Equation 9 and applying the unitary transformation [52]:

$$\widehat{U}=exp\left({\displaystyle \frac{g}{ϵ\left(t\right)}}\left({\widehat{b}}^{-}{\sigma }_{+}-{\widehat{b}}^{+}{\sigma }_{-}\right)\right)$$

we get the single qubit Hamiltonian in the frame rotating at the driving radiation field of frequency ${\omega }_{ext}$ acting on the qubit-resonator system during a specific time:

$$\widehat{H}={\displaystyle \frac{\mathrm{\hslash }}{2}}\left(\Omega +2{\displaystyle \frac{g^2}{ϵ\left(t\right)}}\left({\widehat{b}}^{+}{\widehat{b}}^{-}+{\displaystyle \frac{1}{2}}\right)-{\omega }_{ext}\right){\sigma }_z+\mathrm{\hslash }{\displaystyle \frac{g∆}{ϵ\left(t\right)}}{\sigma }_x+\mathrm{\hslash }\left({\omega }_r-{\omega }_{ext}\right){\widehat{b}}^{+}{\widehat{b}}^{-}+\mathrm{\hslash }∆\left({\widehat{b}}^{+}+{\widehat{b}}^{-}\right)$$

The current in a transmission line is proportional to the relationship between the voltage's magnitude and the line's impedance (I=V/Z). This simple Ohm's Law will hold for a transient period [76].

Three effects can generate a superconducting phase difference around the circuit [11].

• Due to the Josephson effect. When the pairs coherently pass through a JJ, the supercurrent ${\rm I}$ related to the phase difference between electrodes can be represented by ${\rm I}={{\rm I}}_{0}sinsin\varphi$. ${{\rm I}}_0$ is the maximum critical current the junction can maintain. In our model of DNA, the maximum crucial current of the H-bond is ${{\rm I}}_0$.
• Due to the supercurrent, an intrinsic phase shift appears when computing the line integral of momentum $\overrightarrow{p}$ pairs moving along the curve.
• Due to external magnetic flux. Classical electromagnetism and Maxwell's laws show that a magnetic field that changes with time induces an electric field. This, in turn, modifies the momentum of the charged particles by exerting a force on them. The moment of the electron pairs is given by equation $\overrightarrow{p}=2m_e\overrightarrow{v}+2e\overrightarrow{A}$. The first component corresponds to the kinetic part, and the second to the field contribution—the potential vector is $\overrightarrow{A}$.

In our DNA model, RNA pol II changes the external magnetic flux through a P-mechanism and generates a phase difference. RNA pol II sends the system pair of pulses with the same parameters many times (n). The total number of pulses required may explain why the enzyme's catalytic site has 52 repeats (Figure 1A). Then, the number of times (m) the system was in state $|0⟩$ or $|1⟩$ is extracted. With these data, the probability that the system is in a specific state is calculated by applying Laplace's rule. For example, to obtain the probability that it is in state $|0⟩$: $P\left(|0⟩\right)={\displaystyle \frac{m}{n}}$. We can obtain the Rabi oscillations by representing the probabilities as a function of times. At a computational level, a pulse is a finite time series of complex values, where each value represents a pulse amplitude at a given time. In DNA, a pulse is a finite time series of P-reactions.

7. Discussion

The spin 1/2 particles are TLS useful for quantum computation. In the DNA ground state, the $\pi$-ring controlled the hyperfine and electron-mediated nuclear spin interactions. Previously, close connections between classical and quantum information were described in DNA [8]. In this work, we characterized A-T and C-G bp as flow qubits, as they constitute a superconducting ring interrupted by one JJ and crossed by a magnetic flux. Recent studies treated the H-bond between DNA bp as adiabatic systems with spin-orbit coupling. Hubač and Cols. extensively defined qubits formed by Majorana fermions in the H-bond and discussed the entangled states in the bp [78]. The Hamiltonian obtained in this work for the DNA qubits differs from that of the transmon. However, in the limit when L tends to be infinite (omitting the term ${\displaystyle \frac{1}{2L}}$ in Equation 2), our Hamiltonian (Figure 2B) tends to be the transmon Hamiltonian (Figure 2A).

In DNA, each bp is connected to the adjacent by the sugar-phosphate backbone. This structure is similar to a transmission line with distributed parameters. Changes in the magnetic field by the RNA pol II during transcription could flip nuclear spins at resonance. An applied magnetic field breaks electrons’ twofold spin degeneracy. A pulse set detuned from the resonant frequency and enables the state of the qubit to be determined. One question emerged: What is the energy to get the lowest donor-excited state? This work considers high intra-complex P-activity manifested in processive P-reactions of alcohol residues in RNA pol II using ATP as the phosphate donor as the RNA pol II driving force. Multisite P-networks emerge as a new level of signal based on routes encoded into disordered regions of proteins [79]. PPi release during transcription is a signature step in each nt addition cycle [68]. A DNA circuit can be represented as a symmetric and antisymmetric sin function: $y=Asinsinx$and $y=-Asinsinx$. We proposed the graphic representation of the DNA strands’ deformation due to the H-bond break to form the Bloch sphere. The rotations of the qubit are carried out during transcription (Figure 1B).

The coupling of each spin to its neighbors and the magnetic field enables numerous operations on each spin. Measurements are performed by determining the electron spin state effect on the electrons' orbital wavefunction after transferring nuclear spin polarization to the electrons [11]. Also, quantum computation requires a two-qubit controlled rotation proceeding. The result is a rotation of a target qubit when the control qubit is oriented in a specified direction and still unchanged [20]. The measurable signal decreases with the number of qubits [80]. Previous studies reported that more than ten qubits will be challenging. Looking at the spins inside a DNA loop, they are arranged in parallel. In principle, logical operations and measurements can be performed independently and in parallel on each spin in the array.

Circuit QED points to multiple superconducting qubits and 3D cavities with high coherence [21]. Distributed quantum computation arises across numerous single-mode resonators [20]. Single-step multi-qubit phase gates on multiple single-mode resonators mediated by a superconducting bus in circuit QED have been described [81]. The superconducting bus may be possible single-mode resonator interactions, and quantum information is encoded in various single-mode resonators' vacuum and single-photon states [11]. In addition, the pulse engineering technique adjusts the coupling strength between resonators and the bus enhancing fidelity and robustness. The tunable coupling strength plays a broad role in the circuit QED system [20].

Based on the above information, we consider that a DNA gene's basic structure is the nonperturbative circuit QED, and a complete gene would be the linear combination of $n$ circuits in parallel. This parallel arrangement, together with its superconducting nature, contributes to the total resistance being zero, and there is no heat dissipation. DNA gene transcription process includes the combination of $n$ cavities hosting nine qubits and coupled to a shared machine. Many qubits in a single cavity may increase the unwanted qubit-qubit interaction and the cavity decay while decreasing the qubit-cavity coupling strength [20,54,80]. Hence, multiple cavities and quantum operations on qubits distributed in different cavities are mandatory for large-scale quantum information processing.

The RNA pol II bus holds a level structure formed by ground and excited levels, denoted by $|g_1,g_2\dots g_n⟩$ and $|e_1,e_2\dots e_n⟩$, respectively. Adjusting the qubit level spacings achieves the coupling and decoupling of each qubit from its cavity. The P-pulses applied to the bus drive the transitions resonantly between $|e_x⟩$ and $|g_{x+1}⟩$ with Rabi frequency ${\Omega }_x$. Realizing quantum information transfer (QIT) and entanglement with SQUIDs in a microwave cavity is feasible in artificial atom models [82].

8. Conclusions

Biomolecular condensates concentrate functionally related components through multivalent and dynamic interactions in biopolymers without a bounding membrane. Using the transcription condensate architecture and the Electrodinamyc Quantum Theory of Circuits, we further developed ideas about the electrical design of quantum information in DNA. This Theory is an excellent model for characterizing the transcription process and describing the Hamiltonians that define the function in secondary quantification using the methods of Quantum Physics. In this work, we address the graphic and physical-mathematical representation of the DNA qubit as a nonlinear quantum circuit that introduces anharmonity through the H-bond. There are similarities between DNA qubits and RF-SQUIDs coupled to a cavity mode defined by two resonators. However, DNA qubits differ from the transmon qubits.

To characterize DNA as a transmission line or a wire model, the backbone structure of single-stranded DNA must be analyzed. It has connections between successive bp. Hence, the line's superconductors in a given gene are the bp. Each single DNA strand backbone works like a TEM with discrete distributed parameters. In this work, we described the Hamiltonian of one qubit coupled to two resonators. We concluded that every single gene in DNA can be modeled as a combination of $n$ circuit. The RNA pol II, acting like a multifunctional biomolecular machine, produces the system's decoherence like an artificially engineered device. A processive phosphorylation circuit with multiple kinase inputs and in the 52 repeated heptad motif (Tyr₁Ser₂Pro₃Thr₄Ser₅Pro₆Ser₇) controls multi-target qubits. The PPi release step is crucial to the mechano-chemical coupling mechanism during transcription elongation because it is equivalent to the computation pulse. The transcription process involves coupling n circuits, each containing nine qubits.

Evolutionarily conserved biological processes are critical and contain essential secrets that allow a specific function to be carried out. Knowing these deep secrets helps us understand fundamental human and natural development processes. Understanding basic quantum interactions is essential to discovering these biological secrets. Quantum Molecular Biology is an emerging and fascinating field of study. The quantum information in DNA is related to life evolution and preservation.