Mechanics of Twisted DNA Molecule Adsorbed on a Biological Membrane ()
1. Introduction
The Deoxyribose-Nucleic-Acid molecule (DNA) is considered among the most important biological systems for living beings. Since Watson and Crick unveiled its topological structure [1]. This molecule has become a vast field of research. All the genetic information that controls the cell reproduction is coded there, constituting the most sensitive molecule of the cellular nucleic.
A chromosome with protein components is a complete DNA molecule. Amino acids are assembling in the right order to produce the protein. This assembly is triggered by a chemical message carried by the DNA of each gene, tells the cell how to organize this assembly. Monomers that constitute the molecule of DNA are called nucleotides contained information. The structure of the DNA molecule consists of a base and a skeleton alternating phosphate ions and sugar molecules. In DNA, there are four different nucleotides, namely guanine (G), adenine (A), cytosine (C) and thymine (T), whose repeated stacks are formed by AT or GC. The stability of the double helix, i.e. its denaturation resistance depends on its sequence: the G-C rich sequences have a greater resistance to denaturation than the A-T rich one. Because the opening of the G-C base requires the breaking of three hydrogen bonds, while the opening of a pair of A-T involves only two hydrogen bridges. In the case of thermal denaturation, the mean denaturation temperature of a given chain is effectively related to the percentage of G-C contained in the sequence [2].
Several works have been interesting by the study of DNA-interface interactions considered as a subject of general nature. DNA interactions with surfaces can be repulsive leading to confinement or depletion [3] [4]. In this situation, there exist a repulsive interaction between negatively charged phosphates of DNA and phospholipids of membranes and contributes to enclose the genetic material within the cell [5]. The depletion forces also confine DNA within viral capsids [6]. To indicate the repulsive interactions are a common requirement in micro and nano-fluidics [7]. On the one hand, the attractive interactions lead to adsorption on a surface. The adsorption studied is done by x-ray reflectivity [8]. In some recent work [9], a theoretical analysis of such scenarios is based on the self-consistent field theory approach. In this work, they treated the importance of the charge correlation effect.
The thermodynamic magnitudes of DNA are affected by several factors, for example the thermodynamic properties with the interaction of the solvent on the nonlinear dynamical structure of a DNA segment [10], by using a time-independent perturbation approach. Also, for a short fragment of heterogeneous DNA with a stabilizing solvent interaction term, the technique used in this case is the imaginary time path integral formalism that is applied to a nonlinear Hamiltonian [11]. These nonlinearities in DNA dynamics were first emphasized by Englander et al. [12] who interpreted the formation of temporary open segments of base pairs as moving defects propagating coherently along the backbone of the molecule. This phenomenon has been treated like a bubble dynamics in a random DNA sequence, in a study of the localization of denaturation [13].
The denaturation transition is done by the hydrogen bond rupture between two pair bases, so the set of the two strands is transformed into a single strand. This phenomenon is done by several factors, for example, when a sufficiently heated solution of DNA. The temperature at which the DNA strands are half denatured, meaning half double-stranded, half single-stranded, is called the melting temperature (Tm) or denaturation temperature (Td). The amount of strand separation, or melting, is measured by the absorbance of the DNA solution at 260 nm (absorb in the ultraviolet). Nucleic acids absorb light at this wavelength because of the electronic structure in their bases, but when two strands of DNA come together, the proximity of the bases in the two strands quenches some of this absorbance. As the DNA becomes denatured, its absorption of light ultraviolet increases, when the two strands separate, this quenching disappears, and the absorbance rises 30% - 40%. This is called Hyperchromicity. The Hypochromic effect is the effect of stacked bases in a double helix absorbing less ultra-violet light.
For a quantitative study, Thierry Dauxois and Michel Peyrard presented a model for the dynamical structure of DNA that can be considered as an extension of the usual Ising-like statistical approach to the melting curves. Also, Thierry Dauxois and Michel showed by numerical simulation method at constrained temperature show that it provides a good qualitative description of the collective motions of the base pairs, including their large-amplitude fluctuational openings and the emergence of the denaturation bubbles from the thermal fluctuation [14].
In this paper, we focus on the statistical of a DNA molecule on the fluid membrane, our aim is to determine the quantity physics. Our system physics constituted by a fluid membrane fluctuating around a horizontal, plane on this the latter adsorbed a DNA molecule. Firstly, we determine the partition function, from which we establish the free energy of the system. Based on the latter, we derive the magnitudes thermodynamics such as heat energy. As results, we have showed that the constant of harmonic stacking is renormalized, and the DNA becomes more elastic. This constant depends crucially on the membrane parameters. The thermal fluctuations of membrane modify the depth of interaction potential of staking and inverse length. On the other hand, we concluded that the fluctuations of the membrane increase the denaturation temperature. In the transition point we constant that the average separation between base-pair and specific heat are diverging.
The remaining of the presentation proceeds as follows. In Section 2, we presented the physical system with its parameters, and the Section 3 is reserved to formalism for partition function and free energy, the exact study of the denaturation transition is presented in the Section 4, Some concluding remarks are drawn in Section 5.
2. Description of the Physical System and the Basic Equation
The denaturation transition can be achieved in vitro, while submitting the DNA to all chemicals or physical agent capabilities to destabilize links hydrogen, such as pH, temperature, certain solvents, alkaline agents, high ionic concentrations [10], in this paper, the authors investigated the effect of solution concentrations on physical quantities such as specific heat, entropy, melting temperature, and the mean hydrogen bonding stretching.
In the living cells, the DNA molecule is compacted in a biological membrane. Quantitative understanding and optimization of its functions require precise experimental characterization and accurate modelling of DNA properties on the biological membrane. The legitimate question, how the thermal fluctuations of the membrane affect the denaturation temperature of the DNA molecule, and its influence on thermodynamics quantities. Due to the nature of electrically charged DNA, as well as superficially charged cell membrane leads to the phenomenon of adsorption [9]. Therefore, in this investigation, we assume a DNA-molecule adsorbed on a fluctuating membrane, such as each base pair of index n, formed an angle
, with the plan of the membrane (see Figure 1). In this model, the interaction between DNA-molecule and membrane per area has an attraction, as well as hard wall repulsion for a base pair residing at a point on the membrane surface.
On the other hand, we introduce the twist angle
, for the pairs bases adsorbed, we take the equilibrium twist angle value, in our situation
, we choose the model of Dauxois-Peyrard-Bishop (DPB) modified by Zoli [11]. This author introduced the rotation angle between adjacent bases n, and
along the DNA backbone. The twist angle between adjacent base n, and
, we assume constant, therefore
(1)
Figure 1. Biological membrane with a twisted DNA adsorbed.
We consider the adsorption of DNA on the membrane as the adsorption of copolymer on the substrate, which forms loops, so there are a few base-pair are adsorbing on this membrane for each loop. In the other hand, some work show that [15] the distance between adjacent bases in closed case of base pair (~3.3 Å) is smaller than the distance of separation between base pairs (~18 Å), but it's clear that length between base pair in the opening state is greater than the distance between two base, so with these considerations we assume that the two bases adjacent are coplanar. We take for example the first adjacent bases (base 1 and 2), we can establish the first relation between the angels
and
, such as
and the second adjacent bases (base 2 and 3), we find
, etc. (see Figure 2), finally, we get the relationship between
and
(2)
where
is the first angle between adsorbed base pair and membrane, from this relationship, we can determine the maximum number of base pairs (bp) adsorbed on the membrane and contribute to the denaturation transition phenomena by the following relation:
(2a)
As long as adsorption takes place with a few monomers (we take for example
) to form the brush, we can choose for example
, and the number of base pair which contribute to adsorption is
, we deduct that
. In the absence of the membrane, the number of base pairs that contribute to the denaturation transition is of the order of 20 bp for
total number of base pairs [13].
Figure 2. Cross-section of a membrane with an adsorbed DNA molecule.
The double-strand is denoted by the base sequence
(where
and
are the nucleotides with
, of one of the strands, ordered from the 5' to 3' end. The distance between the mass center of the base pairs and a point situated on the surface of the membrane is given by barycenter technique
, where
is the reduced mass,
is the abscissa of the nucleotide
, and
is the distance between two nucleotides
and
(i.e. is the length of the hydrogen bond), since we are interested in the bases of pairs, that is adsorbed on the membrane, we take
or
equal the fluctuation amplitude of the membrane
2.1. The Hamiltonian of Membrane
The membrane is a flexible, continuous surface that has a size L and area
. Membrane local position vectors are, in the Monge representation,
, where
is the height of the membrane. The position of the (almost flat) membrane is specified through the displacement field
. The surface fluctuates around the horizontal plane
. The equilibrium statistical mechanic of the membrane is based on the Canham-Helfrich Hamiltonian [16] [17]
(3)
where
, is the membrane bending rigidity,
is the microscopic membrane surface tension, and
is the projection of r in the reference plane,
is the cartesian operator in the basic
, and
is the projected membrane area. Although we mainly refer to membranes, our development below with this Hamiltonian applies to a variety of membranes and interface including fluid membranes (
) and elastic interface (
).
2.2. Hamiltonian Model of DNA
In our investigation, we choose the model (DPB) modified by Zoli [11] [14], thus generalizing the DPB Hamiltonian, is given by
(4)
where m is the base pair reduced mass,
is the nucleotide velocity, K is the harmonic stacking,
and a are the pair dissociation energy and the inverse length setting the hydrogen bond potential range for the n-th base pair. The first term is the kinetic energy, the second term is the potential energy of the longitudinal links and the last term represents the dissociation energy of the pair (for a very deep reading see the Ref. [14]). The settings in this Hamiltonian are adjusted to match our problem (
,
,
, [18]).
2.3. Interaction between DNA Molecule and Membrane
The DNA molecule and membrane are subject to mutual interaction which we denote by
(5)
the sum of all the base pairs, and continued sum covers the surface of the membrane. The potential that can be used to model interactions such as the interaction between an atom and a surface is the Morse potential. Because of its simplicity, it is not used in modern spectroscopy [19]. However, its mathematical form has inspired the MLR (Morse/Long Range) potential, which is the most used potential energy function for spectroscopic data fitting [19]. The potential chose to mimic the adsorption potential by a simpler one, the Morse potential
(5b)
we replace
by its expression in this potential we find
(5c)
with
and
is the membrane’s surface. In our model, the interaction
between base pair and membrane per area has an attraction, as well as hard wall repulsion (impenetrability) for a base pair residing at a position
on the surface. Such impenetrability by DNA molecule is conceivable even in fluid membrane unless the striking barrier formed by the lipid bilayer self-assembly is disturbed by a strong attraction, (i.e. between hydrophobic segment in the DNA-chain and inner hydrophobic part of the bilayer) leading to chain penetration into the membrane, the choice of the Morse potential to ensure the adsorption phenomena, in a similar work [20], the potential chosen is the hard square-well with depth
and width
. The Parameters for the Morse potential are adjusted to have a good match at the minimum with the true potential (
,
).
3. The Partition Function of System and Free Energy
For a chain containing N units (nucleotide), the classical partition function of the system may be factored as
(6)
where we define as usual
, with T the absolute temperature,
Boltzmann’s constant. The momentums parts are readily integrated to give the familiar kinetic factor for N particles
. The sum on the variable
in Equation (6) is independent of the other variables, so we can compute this sum separately, therefore we get the following relation (see Appendix) [21]
(7)
we are restricted to the first order, and, we have been interested in the phenomenon of breathing, which is affected near of point of equilibrium, then
is weak. Thus, we get the potential expressed by the following relationship
(7a)
The partition function (6) becomes
(8)
where
is the partition function of the free membrane for developing the part related to the twist angle in the last equation, we are introducing sum and difference coordinates
and
. In this work, we also assume the twist angle is very low, which affects the harmonic part, becomes
. Therefore
can be expressed in the form
(9)
where
denotes the total potential energy in
. This integral can be evaluated exactly in the thermodynamic limit of a large system (
) using the eigenfunctions and eigenvalues of the transfer-integral operator [22] [23] [24]
(10)
To solve this integral we change the variable of integration from z to z' (where
), in what follows, we assume
, and introducing
, and performing the Gaussian integrals over
. From this equation, with the new parameter, in the continuum limit approximation, the transfer integral eigenvalue problem can be reduced the following Schrödinger-like equation
(11)
where the harmonic stacking
depend on the
defined by
,
is the new potential and s is the modified entropy by the presence the angle twist
(12)
with
is the entropy determined by T. Dauxois et al. [14], from this relation we can deduce that twist angle minimizes the entropy of the DNA-molecule and the harmonic stacking. The first term is harmonic potential, we can be transformed as a Morse potential at high temperature, we find potential with six pertinent parameters
given by the following relation
(13)
where
is the depth also, and
, this parameter is homogeneous to the inverse of a length. For a temperature
, we can estimate the values of the depth, and the range respectively,
and
. This potential is plotted in Figure 3, this potential reflects standard expression for chemical bonds and, moreover, it has the appropriate qualitative shape: 1) it includes a strong repulsive part for
corresponding to the steric hindrance 2) it has a minimum at the equilibrium position
, 3) it becomes flat for large z, giving a force between the bases that tends to vanish, as expected when the bases are very far apart, this feature allows a complete dissociation of the base pair, which would be forbidden if we had chosen a simple harmonic potential.
The shape of this potential is more adapted to our problem because it ensures the phenomenon of breathing of the DNA-molecule near its minimum has a harmonic form. So, this potential can be well fitted to the Morse potential, we can be rewritten as follows
(14)
Using the relation (13), to find the zero and minimum of this potential (
), then we can find the expression of the depth and the minimum (
). All the parameters depend on the six pertinent parameters, we find
(15)
is the depth of the Morse potential, and
defines the range parameter given by
Figure 3. Reduced Morse potential, versus distance, with the ranges parameters, are
, and the depths
,
,
.
(15b)
according to the values of the parameters, we notice that this last parameter is lower than a, so the width of the potential decreases. On the other hand, the dissociation energy of the pair D (the depth of potential) increases. Note this potential describes not only the hydrogen bonds and the repulsive interactions of phosphate groups in DNA but all interactions in the system [14]. Let comment on this form of potential: the value
corresponds to a closed base pair as in the Ising model, but now z can increase continuously to infinity if the two bases separate completely as in DNA denaturation. The variable z can even take negative values, corresponding to a compression of the bond linking the bases with respect to its equilibrium length. Large negative values will be forbidden by steric hindrance, which is introduced in the model by the potential linking the bases in a pair [15].
Now we determine the free energy. The calculation is similar to the one performed by Krumhansl and Schrieffer [23] for the statistical mechanics of the
field. It yields
where
is the lowest eigenvalue of the operator. Therefore, we can then compute the free energy of our model as the sum of the different contribution in
,
(16)
where
, is the free energy of the free membrane, with,
.
The following step is to determine
. The presence of the membrane and thermal fluctuations affects the eigenfunctions and eigenvalues of DNA molecules, so,
and
. [20] The ground state energy
, as is well-known in quantum mechanics, is alternatively obtained by a variational method, by varying a variational function
,
is given the following relation
(17)
Since we have interested, by the distance of separation of the base pair of the DNA molecule, one stands in space united dimensional, to facilitate our development further, we use the relation [20]
(18)
where
is the ground state variational function of the 1D free energy operator,
(19)
where
. With necessary ingredients, the expression of the minimum free energy
can be written as follows
(20)
where the pure DNA-molecule contribution is
(21a)
With
(21b
(21c)
Bay a Fourier transformation for
one can reduce the partition function to another form
(22)
Here
(23)
is the renormalized surface tension at the state
corresponding to the minimum of
, given by the Equation (20). From the partition function, we obtain the total free energy base per pair induced by its interaction with membrane
(24)
where
(25a)
With
(25b)
is the partition function of the nude membrane. Thus, we have the following expression of the free energy
(26)
where the first term is the membrane contribution. The condition that f is the minimum for
leads to the new following Schrödinger-like equation
(27)
this equation for
minimizing
, Equation (21), where
is the renormalized spring constant given by
(28)
it is easy to find
(28a)
where
is given by the following equation
(28b)
here l and L are the microscopic and macroscopic lengths that define the upper and lower wave number cutoffs. From the relationship (28a) we conclude that there is a supplementary force added to its previous
such as
,
and
, (
unit vector), the sign (-) signifies that this force is the inverse of the longitudinal force. These later have two different effects, one pulling the nucleotide, and the other returning it to its balanced position (similar to the nucleotide in each strand a mass attached by two springs). One deducts that the undulation of the membrane, decreases the elasticity constant K, so every strand in the double helix becomes more elastic.
4. Exact Study of Denaturation Transition
4.1. Eigenvalue and Ground State
Equation (27) is formally identical to the Schrödinger equation for a particle in a Morse potential so that it can be solved exactly [25]. So, the melting temperature of DNA adsorbed on a fluctuating membrane is given by
(29)
where
is the denaturation temperature in the absence of membrane [14], and the quantity
is given by the following equation:
(29a)
According to the values of the parameters
, it is deduced that the presence of the membrane increases the denaturation temperature of the DNA molecule. Also, this temperature depends on all the parameters of the problem.
Let's comment this result: like a set of external factors, the presence of the cell membrane leads to a change in the melting temperature of DNA, the presence of the membrane fluctuations increasing the melting temperature denaturation temperature. According to the parameters, this temperature increases by ~30%. Therefore, when the temperature increasing, the fluctuations of membrane increasing also, which implies that the membrane undergoes a phase transition from a flat phase to a crumpled phase, i.e. its orientations in different points of the membrane are decorrelated, where the persistence length
, is above the linear size L, [26]. A very recent work [27], show a good agreement between experience and numeric study, that the temperature change affects the DNA twist, and suggest that the untwisting of DNA with temperature is predominantly due to changes in DNA structure for defined backbone substates.
The eigenfunction of this system is given by
(30)
where
, and
it is the normalization constant, this equation can be rewritten according to the wave function of the DNA in the absence of the membrane we obtain
(30a)
where
is the eigenfunction of DNA with
, and
is the function given by the following relation
(30b)
with
, is a constant for a free DNA, and
, is given by
(30c)
We notice the eigenfunction is amplified, by a factor
. The eigenvalue of the associated ground state
(31)
where
, it is the eigenvalue of free DNA given by
(31a)
and
, is the excess energy of ground state its expression is
(31b)
It is noted that the thermal fluctuations of the membrane disturb the eigenvalue of the DNA molecule, as well as the associated wave function. According to the values of the parameters included in Equation (31a), the difference
, is very low, so the energy specter of
, stay near to that
, we conclude that membrane undulations cause a small perturbation of molecule DNA
4.2. Thermodynamic Magnitudes at the Critical Point
4.2.1. Eigenvalue and Ground State
We are recalled that to study of the denaturation transition of the DNA molecule, the important quantity which gives a measure of the extent of denaturation of the molecule is the mean stretching of the hydrogen bonds, which given by
(32)
We take the results obtained in the case of unbinding transition from two interacting manifolds (strings or bilayer membranes), we find that near the transition temperature (is the similar case) [28], this average distance scale exactly as follows
(33)
This average distance is plotted in Figure 4. One notices that the pace of the curve coincides with the one to find by the continuum approximation and the exact numerical calculation [14]. All methods show a divergence of the mean distance at
, corresponding to the stretch of links hydrogens, from a given temperature. We concluded that, for the weak values of
, the DNA-molecules are closed and opened for the strong values.
4.2.2. Surface Tension and Specific Heat
A quantity that reflects the ability of the molecule of DNA adsorbed on a fluctuating membrane to accumulate energy in thermal form, it is the specific heat. Before calculating the specific heat, we attempt to determine the expression of
. We consider below critical behaviors near
. In terms of small reduced temperature
, we have all ingredients to find the expression of the surface tension depending on the temperature difference
Figure 4. Variation of the reduced mean distance
versus reduced temperature.
(34a)
where
(34b)
this expression is according to the difference in the temperature. A simple dimensional analysis shows that
. It is necessary to note that
at
.
From the relationship (26) we can derive the specific heat
(35)
This specific heat includes a regular part, which is continuous at a transition point, and a singular part, which behaves as a power law with an exponent.
(35a)
with regular part given by
(35b)
Our objective is to study the critical comportment of the specific heat at denaturation temperature (singular part), therefore we get the following expression
(36)
with a factor, where
, is a critical exponent. Our system is immersed in three-dimensional space (
), but the study of the spherical model has explicitly showed the role of the space dimension on critical phenomena. This model allows calculating exactly the critical properties for an arbitrary space dimension. In this case
, (in spherical model
), since
we have
, it is a universal critical exponent, the specific heat plotted in Figure 5.
Figure 5. The variation of the specific heat
, versus reduced temperature.
5. Conclusions and Remarks
In this work, we have studied the influence of the thermal fluctuations membrane on the properties thermodynamics of DNA molecule adsorbed on the fluid membrane. We formulated the problem in terms of a mesoscopic-level partition function and free energy, using the approach variational method, which renormalizes the elasticity constant of DNA. The undulations of the membrane, decrease the elasticity constant K. We deduce that nucleotides in each strand vibrate in the longitudinal direction. We can deduce that the modes of vibration also change. This is because the force induced by the membrane, tends to bring the nucleotide back to the equilibrium position.
Consequently, the Schrödinger-like equation which describes the state of the DNA molecule in the absence of the membrane becomes a function of the parameters of our problem. This equation is solved exactly. The resolution of this equation gives the bound state and the associated wave function. From the Schrödinger-like equation, we have determined the melting temperature. This latter depends also on all problem parameters. We have used the eigenfunction associated the ground state to compute the average separation
between base pairs. This quantity diverges when
. Another physical quantity has been changed by the adsorption of DNA on the membrane; it is the surface tension
, where
vanishes in the vicinity of the denaturation temperature
. In this investigation, we have also found that the heat specific diverges at melting temperature with universal exponent
.
For the implementation of the denaturation transition of the DNA molecule adsorbed on a membrane, the experimental results concerning the absorption of ultraviolet (UV) by a pure membrane poly (acrylic) acid (PAA) doped with graphite [29], these authors determine the absorbance as a function of the wavelength. They ascended that the absorbance increases with the increase of the concentration of graphite (in the absence of graphite there is a weak absorption [29]). For our case it is necessary to use a pure membrane (absence of graphite and DNA) and UV with a fixed wavelength at 260 nm, and then determine the maximum absorbance for different temperature values, and then use a molecule of DNA adsorbed on a membrane (experiments have been performed on artificial DNAs, which are homopolymers, i.e. have only one type of base pairs [15]), and then measure the absorbance for each temperature value, in the case where the absorbent is maximizing for a certain temperature, that is called the melting or denaturation temperature. It should be noted that as the DNA becomes denatured, its ultraviolet light absorption increases.
The relation (28a) raises interesting questions is the influence of the membrane on the property’s DNA, especially the force-extension curve of a single DNA molecule. The first experience obtained in 1992 [30] and repeated with taller precision by several groups a few years ago [31] [32], can be repeated once again for different molecules of DNA with different values of the stiffness constant K, and look at the melting hysteresis comportment.
Acknowledgements
We are grateful to Professor M. Benhamou for illuminating discussions.
Appendix
In this appendix, we will try to develop the expression of following integral
(A1)
To develop this integral we make the first variable change
, we find
(A2)
we use also the second variable change as follow
, we see that
, so the integral becomes [21]
(A3)
We are restricting to the first order
, we find a simple development of our integral
(A4)
The term
, is as similar as Morse potential.
Parameters and Nomenclature
: Twist angle between each base pair.
: Angle formed between each base pair and membrane.
: The maximal number of monomers adsorbed on the membrane.
: The abscissa of the mass center of point
.
: The abscissa of nucleotide (
).
: Mass of nucleotide (
).
m: The reduced mass.
: Fluctuation amplitude of the membrane.
r: Vector position.
: Projection of r in the reference plane.
: Membrane bending rigidity.
: Microscopic membrane surface tension.
K: Harmonic stacking.
D0: Dissociation energy.
a: Inverse length.
U0: Dissociation energy (DNA-membrane).
b0: Inverse length.
s: Entropy.
: Partition function of system.
: Free energy.
: Free energy of the free membrane.
: Denaturation temperature.
: Denaturation temperature in the absence of membrane.
: Eigenfunction.
: Eigenvalue of the associated ground state.
: Mean stretching of the hydrogen bonds.
: Specific heat.