Submit manuscript...
Journal of
eISSN: 2572-8466

Applied Biotechnology & Bioengineering

Research Article Volume 2 Issue 4

Structural hierarchy in DNA hydrogels

Bohidar HB,1,2 Pandey PK,1 Rawat K,2,3 Aswal VK,4 Kohlbrecher J5

1Department of Physical Sciences, Jawaharlal Nehru University, India
2Department of Nanosciences, Jawaharlal Nehru University, India
3Inter University Accelerator Centre, India
4Department of physics, Bhabha Atomic Research Centre, India
5Laboratory for Neutron Scattering, Switzerland

Correspondence: Bohidar HB, School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110067, India, Tel 91 11 26704637, Fax 91 11 2674 1837

Received: December 28, 2016 | Published: March 10, 2017

Citation: Pandey PK, Rawat K, Aswal VK, et al. Structural Hierarchy in DNA Hydrogels. J Appl Biotechnol Bioeng. 2017;2(4):144-150. DOI: 10.15406/jabb.2017.02.00038

Download PDF

Abstract

Sol-gel transition in aqueous DNA solutions was probed to determine the gelation time tgel and temperature Tgel. It was remarkable to observe that Tgel increased from 36 to 57°C and tgel linearly decreased from 700 to 350 s with increase in DNA concentration. The hydrogels were studied by using small angle neutron scattering (SANS) ([DNA]≤3%(w/v)) to determine the underlying hierarchy of length scales. From structure factor profile analysis, three distinct signatures were obtained:

  1. Radius of gyration Rg≈40±3nm, which assigned a length of 140nm to DNA strand
  2. Geometry of scattering moiety defined by the value of α≈2±0.1 indicating Gaussian chain behavior
  3. Correlation length, ξ increased from 0.50 to 3.30nm and the inter-chain spacing d0 decreased from 15.6 to 9.8nm, with increase in DNA concentration.

Physical networks were formed for concentration, c*≥2%(w/v) and the system involved at least five identifiable length scales. A revisit has shown that hierarchical structure of DNA hydrogel owes its origin to considerable self-organization at the molecular level dependent on biopolymer concentration.

Keywords: DNA hydrogels, neutron scattering, guinier-regime, power-law, broad peak region, correlation length, self-assembly

Abbreviations

SANS, small angle neutron scattering; EGDE, ethylene glycol diglycidyl ether

Introduction

DNA is found to be an extraordinary biopolymer because it can be used as molecular tool to construct many other objects such as nano devices, periodic arrays and genetic engineering templates. DNA plays a vital role in many biological processes like gene therapy. But due to its large extended chains it cannot be delivered directly into the cells.1 DNA hydrogels are eco-friendly because these are biocompatible and biodegradable. Unlike other hydrogels, DNA hydrogel are different in terms of efficiency and strong branched cross-linking.2 The benefit of it is one can achieve gelling under physiological conditions that allows for the encapsulation of drugs, proteins and mammalian cells in the sol phase instead of drug-loading step and denaturing conditions.3 A strongly interactive PVA and DNA cryogel was prepared and it showed a good mechanical resistance due to its heterogeneous porous structure.4 For natural DNA, formaldehyde and metal compounds such as arsenic, chromate and nickel are widely used as cross-linkers which are known carcinogens.5,6 DNA-cross-linked PAA m gel could be used to release nano particles for potential drug-delivery vehicles in response to right stimulus .7 DNA hydrogels reversibly respond to thermal stimulus, by switching between the gel and sol state (transition temperature) and enzymes when restriction sites are inserted into one of the building blocks.8

DNA hydrogels have a wide range of biomedical application in tissue engineering and drug-delivery system.9 Lee et al inspired by silkworm and spiders gave the evidence for the formation of DNA hydrogel fiber with self-entanglement prepared in the presence of ionic liquid.10 Moran MC et al.9 prepared DNA gel particles and studied swelling, dissolution behavior and DNA release.9 They made DNA gel particles in the presence surfactant CTAB and protein lyzozyme, but without using any cross-linker or organic solvents. Orakdogen N et al.11 showed evidence of strain hardening in both physical and chemical DNA gels.11 Physical gels were prepared by physical cross-linking of DNA with heating-cooling cycle, while chemical gels were prepared by adding a cross linker like EGDE (ethylene glycol diglycidyl ether). Visco elastic properties of DNA gels and the crossover frequency for DNA in solutions that vary with concentration have been reported.12,13 Liu Y et al.14 described the relaxation time of double and single stranded DNA on the basis of viscosity measurement in an oscillatory flow and by stress relaxation measurements.14 DNA hydrogel relaxation dynamics has been found to be hierarchical in concentration.15 The differential shrinking of hydrogels in presence of NaBr and cationic surfactant has been reported.16 Rhelogical behavior of DNA hydrogels indicated that storage modulus G’ and loss modulus G’’ are frequency independent which is due to dominance of viscoelastic relaxation of the networks at lower frequencies.17

Hierarchical self-assembly in DNA gels is the primary factor which distinguishes and makes it versatile. At appropriate concentration and controlled flexibility, one-pot self-assembly of DNA gives rise to very small nanometer range tetrahedral, dedecahedra and bucky-ball size complex structures. Therefore, self-assembly of DNA allows formation of complex 3-D structures.18 DNA gels exhibit universal scaled stiffening behavior that can be reproduced by a wormlike chain model keeping in mind the entropic elasticity of DNA strands.11 For tissue engineering applications, due to biocompatibility and permeability of hydrogel, one can design the required three-dimensional constructions.19 In spite of the availability of a large volume of literature on DNA, the micro-structural probing of their hydrogels by SANS remains poorly explored. In this report, we have comprehensively studied DNA hydrogels using SANS to develop a general understanding of hierarchical self-organization of these networks. Considering the importance of the nucleic acids in genetic engineering and pharmaceutics, these results are of importance.

Materials and methods

The salmon testis dsDNA sodium salt (2000 bp, molecular weight 1.3×106 g/mol) and D2O (base solvent) were purchased from Sigma-Aldrich, USA. In the medium of (0.15M) sodium chloride and 0.015M sodium citrate, it has 41.2% G-C content and 87°C melting temperature. DNA was dissolved in D2O in different amounts (1-3% (w/v)) at 60°C, followed by heating to 90°C under stirring to denature the biopolymer. Finally, samples were gradually cooled to room temperature (25°C) while their viscosity was measured in real time (viscometer model SV-50, A & D Company, Japan) and the gels were subjected to SANS studies. Rigid gels were formed for DNA concentration, c*≥2% (w/v). The SANS experiments were carried out using the SINQ neutron source at PSI, Villegen, Switzerland. These measurements were carried out with 1.6 nm neutron wavelength where the scattering wave vector (q) was scanned in the range of 0.03-3.4 nm-1.20 The scattering profile, I(q) vs q, exhibited characteristic signature pertaining to the internal micro-structure of the samples. All the measurements were performed using a 2D position sensitive detector of size (96cm x 96cm), source wavelength of 0.8nm, and sample to detector distance of 28 and 20m and path length of quartz cuvette of 2mm. All concentrations are in (w/v) unless otherwise stated.

Results and discussion

Sol-gel Transition

DNA solutions prepared with different DNA concentrations were analyzed for their characteristic sol-gel transition behavior. Time dependent evolution of viscosity η is shown by Figure 1 for different samples. The time dependent growth of viscosity could be described by the scaling relation given by 21

η r = ε k ; t  t gel MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacqaH3oaApaWaaSbaaeaapeGaamOCaaWdaeqaa8qacqGH9aqp cqaH1oqzpaWaaWbaaeqabaWdbiabgkHiTiaadUgaaaGaai4oaiaabc cacaWG0bGaeyOKH4QaaeiiaiaadshapaWaaSbaaeaapeGaam4zaiaa dwgacaWGSbaapaqabaaaaa@497C@ (1)

Here, relative viscosity η r = η η 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacqaH3oaApaWaaSbaaeaapeGaamOCaaWdaeqaa8qacqGH9aqp daWcaaqaaiabeE7aObqaaiabeE7aO9aadaWgaaqaa8qacaaIWaaapa qabaaaaaaa@412E@ , η0 is solvent viscosity, ε =( 1 t t gel ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacqaH1oqzcaqGGaGaeyypa0ZaaeWaa8aabaWdbiaaigdacqGH sisldaWcaaWdaeaapeGaaeiDaaWdaeaapeGaaeiDa8aadaWgaaqaa8 qacaqGNbGaaeyzaiaabYgaa8aabeaaaaaapeGaayjkaiaawMcaaaaa @44B1@ and tgel is gelation time.

The plot between ηr vs. t shown in Figure 1A, provided the characteristic time tgel obtained from the change in slope which was found to decrease with DNA concentration. It was found that as the DNA concentration was increased from 1 to 3%, the viscosity of DNA sol increased gradually due to the formation of inter connected physical networks. To fit relative viscosity data to time dependent power-law function was used and corresponding power-law exponent k was determined. Theoretical value for Rouse dynamics corresponds to k= 0.7 for conducting and k= 1.3 for percolating gel network.22,23 In our case, k= 0.13 ±0.02, implying invariance of the network growth process with DNA content. Although it was an indicative description of gelation, but on the basis of experimental data, it was possible to give a reasonable description to gelation kinetics. The gelation time was tgel ≈ 700 s for 2% gel, which reduced to 350 s when the DNA concentration was 3%. This clearly indicated that at higher concentration due to the propensity of DNA strands the physical entanglement was highly favorable that caused earlier gelation (Figure 2).

No gelation transition was noticed when the biopolymer concentration was less than 2%. The temperature of gelation was 36°C for 2% sample which enhanced to 57°C for 3% DNA gels. It is remarkable to note that just by raising the concentration it was possible to increase the gelation temperature by 63 %. At the same time the gelation time reduced by almost 50%. Further, the data shown in Figure 3, could be fitted to a straight line (chi-squared= 0.99) with a slope of »23 (°C (v/w)). The linear dependence of gelation temperature with DNA concentration was not reported hither to. This result will have significant bearing on the applications of DNA hydrogels.


Figure 1

  1. Temporal dependence of viscosity. The arrows indicate gelation time tgel.
  2. Relative viscosity ηr, as function of reduced time ε =( 1 t t gel ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacqaH1oqzcaqGGaGaeyypa0ZdamaabmaabaWdbiaaigdacqGH sisldaWcaaqaaiaadshaaeaacaWG0bWdamaaBaaabaWdbiaadEgaca WGLbGaamiBaaWdaeqaaaaaaiaawIcacaGLPaaaaaa@446D@ . Solid lines are fitting of data to Equation (1).

Figure 2 Variation of time of gelation tgel as function of DNA concentration.



Figure 3

  1. Variation of viscosity with temperature shown for different DNA concentrations.
  2. Variation of gelation temperature Tgel with DNA concentrations and
  3. Temporal temperature dependence. Note the Newtonian cooling rate 0.069°C/s.

Static Structure Factor

The static structure factor profiles obtained from different DNA hydrogel samples are shown in Figure 4. The raw data was examined by checking their 1/I(q) vs q2 and 1 q MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaSaaae aacaaIXaaabaWaaOaaaeaacaWGXbaabeaaaaaaaa@3A66@ vs q2 behavior which correspond to Ornstein-Zernike 24 and Debye-Bueche functions.25 The statistical accuracy of this fitting was poor. The best fit to the data was found when it was fitted to an empirical function given by equation (2). In this formalism, we have identified three ranges of scattering vectors q:

  1. Low q-range (0.05nm-1≤ q ≤ 0.11nm-1) was fitted to Guinier function given by I( q ) ~   I G  exp( R g 2 q 2 3 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaqGjbWaaeWaa8aabaWdbiaabghaaiaawIcacaGLPaaacaqG GcGaaiOFaiaabckacaqGGcGaaeysa8aadaWgaaqaa8qacaqGhbaapa qabaWdbiaabckacaqGLbGaaeiEaiaabchadaqadaqaamaalaaabaGa eyOeI0IaaeOua8aadaqhaaqaa8qacaqGNbaapaqaa8qacaaIYaaaai aabghapaWaaWbaaeqabaWdbiaaikdaaaaabaGaaG4maaaaaiaawIca caGLPaaaaaa@4E86@

where Rg is known as radius of gyration of scattering moiety.

  1. In the intermediate q-range (0.12 nm-1≤ q ≤0.21 nm-1), power-law function was used for data fitting IPL (q) ~ q where α is related to the geometry of scattering moiety and
  1. In the high q-range (q ≥ 0.22 nm-1) the broad peak was fitted to a modified Lorentzian function given by, I( q )= I BP ( 1+( q q 0   )ξ ) m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaqGjbWaaeWaa8aabaWdbiaabghaaiaawIcacaGLPaaacqGH 9aqpdaWcaaqaaiaadMeadaWgaaqaaiaadkeacaWGqbaabeaaaeaada qadaWdaeaapeGaaGymaiabgUcaRmaabmaapaqaa8qacaqGXbGaeyOe I0IaaeyCa8aadaWgaaqaa8qacaaIWaaapaqabaWdbiaabckaaiaawI cacaGLPaaacaqG+oaacaGLOaGaayzkaaWdamaaCaaabeqaa8qacaqG Tbaaaaaaaaa@4C32@  where q0 is related to peak position of scattered peak and ξ is correlation length for polymer chain and m is a fitting parameter (Figure 4). Thus, the structure factor profile could be described by

I( q )= I G  exp( R g 2 q 2 3 )+ I PL ( q ) ~  q α + I BP ( 1+( q q 0   )ξ ) m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGjbWdamaabmaabaWdbiaadghaa8aacaGLOaGaayzkaaWd biabg2da9iaabMeapaWaaSbaaeaapeGaae4raaWdaeqaa8qacaqGGc GaaeyzaiaabIhacaqGWbWaaeWaaeaadaWcaaqaaiabgkHiTiaabkfa paWaa0baaeaapeGaae4zaaWdaeaapeGaaGOmaaaacaqGXbWdamaaCa aabeqaa8qacaaIYaaaaaqaaiaaiodaaaaacaGLOaGaayzkaaGaey4k aSIaamysa8aadaWgaaqaa8qacaWGqbGaamitaaWdaeqaamaabmaaba Wdbiaadghaa8aacaGLOaGaayzkaaWdbiaabccacaGG+bGaaeiiaiaa dghapaWaaWbaaeqabaWdbiabgkHiTiabeg7aHbaacqGHRaWkdaWcaa qaaiaabMeapaWaaSbaaeaapeGaaeOqaiaabcfaa8aabeaaa8qabaWa aeWaa8aabaWdbiaaigdacqGHRaWkdaqadaWdaeaapeGaaeyCaiabgk HiTiaabghapaWaaSbaaeaapeGaaGimaaWdaeqaa8qacaqGGcaacaGL OaGaayzkaaGaaeOVdaGaayjkaiaawMcaa8aadaahaaqabeaapeGaae yBaaaaaaaaaa@67C9@ (2)

Figure 4 SANS correlation between static structure factor I(q) and scattering vector q for different DNA concentrations. These correlations are best fitted with Guinier, Power-Law and Broad peak (Lorentzian) functions in different q-ranges.

This mode of analysis yielded excellent chi-squared values (>0.97) in all the cases. The cross-over between low and intermediate and intermediate and large-q values are identified as q* and q**, respectively. These defined the signature length scales given by L * ( = 2π q * ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGmbWdamaaCaaabeqaa8qacaGGQaaaamaabmaapaqaa8qa cqGH9aqpdaWcaaWdaeaapeGaaGOmaiaabc8aa8aabaWdbiaabghapa WaaWbaaeqabaWdbiaabQcaaaaaaaGaayjkaiaawMcaaaaa@4155@  and L ** ( = 2π q ** ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGmbWdamaaCaaabeqaa8qacaGGQaGaaiOkaaaadaqadaWd aeaapeGaeyypa0ZaaSaaa8aabaWdbiaaikdacaqGapaapaqaa8qaca qGXbWdamaaCaaabeqaa8qacaqGQaGaaeOkaaaaaaaacaGLOaGaayzk aaaaaa@42B0@

which plotted for different hydrogel samples. This biopolymer had an overlap concentration of c*(≥ 2%) above which the solution was in gel state (Figure 5). Remarkably, the value of L*(=50nm) coincides with the persistence length of DNA.26,27 The L** value was 26 nm in the sol and about 22 nm in gel state (Figure 5). Thus, to a very good approximation, L* ≈ 2L**. In the low q-region ( q <3 R g (1) ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qadaqadaqaaiaadghacaqGGaGaeyipaWJaaG4maiaadkfapaWa a0baaeaapeGaam4zaaWdaeaapeGaaiikaiabgkHiTiaaigdacaGGPa aaaaGaayjkaiaawMcaaaaa@42BC@ the scattering factor was approximated by Guinier-regime.

I( q ) = I G  exp( R g 2 q 2 3 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGjbWdamaabmaabaWdbiaadghaa8aacaGLOaGaayzkaaWd biaabccacqGH9aqpcaqGjbWdamaaBaaabaWdbiaabEeaa8aabeaape GaaeiOaiaabwgacaqG4bGaaeiCamaabmaabaWaaSaaaeaacqGHsisl caqGsbWdamaaDaaabaWdbiaabEgaa8aabaWdbiaaikdaaaGaaeyCa8 aadaahaaqabeaapeGaaGOmaaaaaeaacaaIZaaaaaGaayjkaiaawMca aaaa@4BE6@  (3)


Figure 5

  1. Variation of q* and q** with different DNA concentrations.
  2. Variation of L* and L** with different DNA concentrations. Overlap concentration c* is clearly identified here (arrow).

Where IG is the structure factor at very low q and the fitting yielded value of radius of gyration Rg ≈ 40±3nm in the sol state independent of DNA concentration (Figure 6). According to scaling laws correlation length dependence of ξ with concentration is given as.28

ξ=  R g ( c c * ) v ( 13v ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaqG+oGaeyypa0JaaeiOaiaabkfapaWaaSbaaeaapeGaae4z aaWdaeqaamaabmaabaWaaSaaaeaapeGaae4yaaWdaeaapeGaae4ya8 aadaahaaqabeaapeGaaeOkaaaaaaaapaGaayjkaiaawMcaamaaCaaa beqaamaalaaabaGaamODaaqaamaabmaabaGaaGymaiabgkHiTiaaio dacaWG2baacaGLOaGaayzkaaaaaaaaaaa@48D7@ (4)

Where c*(»2 %) is the typical overlap concentration and ν is excluded volume exponent of the system. From this scaling relation, using ξ = 0.5- 3.3 nm (deduced from equation 6), Rg= 40nm and c= 2% and scaling exponent obtained was v~0.40-0.45. Thus, in the sol state DNA strands were in the poor solvent environment, which changed to marginal solvent situation as the gelation temperature was reached. Therefore, the sol-gel transition was concomitant with change in the solvent quality for this biopolymer. Since, this change was in the favorable direction of solubility, there was no phase separation, but a clear transition to gel phase occurred due to a propensity in the physical entanglement of the DNA strands. We use equation (4) with caution because it strictly applies to chemically cross-linked networks and gels. In the intermediate q-regime (0.12 nm-1≤ q ≤0.21 nm-1), power- law dependence of structure factor was defined by

I PL ( q )~ q α MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGjbWdamaaBaaabaWdbiaadcfacaWGmbaapaqabaWaaeWa aeaapeGaamyCaaWdaiaawIcacaGLPaaapeGaaiOFaiaadghapaWaaW baaeqabaWdbiabgkHiTiabeg7aHbaaaaa@42EB@ (5)


Figure 6

  1. Correlation for I(q) versus q fitted by Guinier and Power-Law functions for different DNA concentrations at room temperature.
  2. Variation of radius of gyration and geometrical scattering factor α with different DNA concentration. Arrow demarcates the two regions.

The least-squares fitting of data yielded the value for α ≈ 2.2±0.1 independent of DNA concentration (Figure 6B), which attributed Gaussian chain behavior (theoretical value=2) to these strands. SANS data in the high q-region was described by a signature broad-peak, a profile normally noticed in scattering from amorphous soft materials. In this formalism, the typical distance between scattering heterogeneities is explained by d-spacing, d0which is the peak position of the broad correlation peak located at q0.Theq0value signifies the peak position of different morphologies such as in lamellar, cylindrical, spherical or for bi-continuous structures.29 Here in the case of hydrogels one can see the significance of such a model. We identify d0 as average distance between the polymer-poor and polymer-rich regions. The following functional form (modified Lorentzian function) described the broad peak feature seen in our system Peak position q0 is related to d0 through the relation

d 0 = 2Π q 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamizam aaBaaabaGaaGimaaqabaGaeyypa0ZaaSaaaeaacaaIYaGaeuiOdafa baGaamyCamaaCaaabeqaaiaaicdaaaaaaaaa@3F7B@ (7)

Figure 7 data clearly shows that for the sample with 1% DNA concentration there was no correlation peak. As we increased the DNA concentration, the existence of broad peak became more prominent and its position shifted to lower-q0 value. Thus the value for d0 decreased from 16.0±1.0 to 10.0 ±0.8 nm which signified that the hydrogel system had become compact and dense. Incidentally, the value of the exponent m≈ 2±0.2. Secondly, we found that with increase in DNA concentration (1 to 3%), the correlation length, ξ increased from 0.50 to 3.3nm. It can be inferred that as DNA concentration increased there was formation of network structures due to physical entanglement of DNA strands, which was predominant at higher DNA concentrations. This corresponded to dominant scattering arising from these structures. This behavior is depicted in Figure 8, which, in addition, defines the overlap concentration c*. Based on the SANS data, a schematic sol-gel phase diagram can be constructed which is shown in Figure 9. In the sol state, the mean inter strand distance was d0 = 16±1 nm, Rg =40±2 nm and correlation length (distance of minimum approach) ξ= 0.5±0.04 nm. These values changed to d0 = 10.0±1.0 nm and ξ=3.3±0.4 nm at the highest concentration used in this study. Thus, it was clearly seen that the self-assembly leading to the gelation involved at least five length scales. The other three parameters showed strong concentration dependence.

Figure 7 Correlation between I (q) versus q obtained from SANS measurement for high q regime performed on different DNA concentrations at room temperature. Note the red-shift of q0 (arrows).


Figure 8 Variation of

  1. Correlation length with different DNA concentration.
  2. D- Spacing with different DNA concentrations.

Figure 9 Schematic of sol- gel transition in DNA solution.

Conclusion

It was found that self-assembly of DNA strands precedes gelation transition. It was possible to extract micro-structural information about DNA strands with SANS measurement. The structure factor profile, I(q) vs q data, was divided into three different q-regions and were treated with specific fitting functions:

  1. low-q region (0.05nm-1≤ q≤ 0.11 nm-1), I( q )=  I G  exp( R g 2 q 2 3 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaqGjbWaaeWaa8aabaWdbiaabghaaiaawIcacaGLPaaacqGH 9aqpcaqGGcGaaeysa8aadaWgaaqaa8qacaqGhbaapaqabaWdbiaabc kacaqGLbGaaeiEaiaabchadaqadaqaamaalaaabaGaeyOeI0IaaeOu a8aadaqhaaqaa8qacaqGNbaapaqaa8qacaaIYaaaaiaabghapaWaaW baaeqabaWdbiaaikdaaaaabaGaaG4maaaaaiaawIcacaGLPaaaaaa@4C44@  (Guinier- regime)
  2. The intermediate q-region (0.12 nm-1≤ q≤ 0.21 nm-1), I(q)~ q (Power-law) and
  3. The asymptotic region (q ≥ 0.22 nm-1), I( q )=  I L ( 1+( q q 0 )ξ ) m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaqGjbWaaeWaa8aabaWdbiaabghaaiaawIcacaGLPaaacqGH 9aqpcaqGGcWaaSaaaeaacaWGjbWaaSbaaeaacaWGmbaabeaaaeaada qadaqaaiaaigdacqGHRaWkdaqadaqaaiaadghacqGHsislcaWGXbWa aSbaaeaacaaIWaaabeaaaiaawIcacaGLPaaacaqG+oaacaGLOaGaay zkaaWaaWbaaeqabaGaamyBaaaaaaaaaa@4AD2@  (Broad Peak region where q0is related to peak position of scattered peak and ξ is correlation length).

The crossover between regions between (i) and (ii) and between (ii) and (iii) designated as q* and q** defined two other length scales, L * = 2Π q* MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGmbWdamaaCaaabeqaa8qacaGGQaaaa8aacqGH9aqpdaWc aaqaaiaaikdacqqHGoauaeaacaWGXbGaaiOkaaaaaaa@3F78@  (=50 nm independent of concentration) and L ** = 2Π q** MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGmbWdamaaCaaabeqaa8qacaGGQaGaaiOkaaaacqGH9aqp daWcaaqaaiaaikdacqqHGoauaeaacaWGXbGaaiOkaiaacQcaaaaaaa@40C5@  (= 26 to22 nm decreasing with increase in DNA concentration). From statistical analysis of the structure factor data, we obtained the radius of gyration value, Rg≈ 40±3nm. Assuming, the strands to be rod-like, the length of the strands L can be determined from, R g 2 = L 2 12 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaieaaaaaa aaa8qacaWGsbWdamaaDaaabaWdbiaadEgaa8aabaWdbiaaikdaaaWd aiabg2da9maalaaabaWdbiaadYeapaWaaWbaaeqabaWdbiaaikdaaa aapaqaaiaaigdacaaIYaaaaaaa@401E@ which gives L»140 nm which is quite realistic. In the high-q region a broad correlation peak was detected and the broad peak fitting analysis yielded correlation length, ξ increasing from 0.50nm to 3.30nm and the inter strand spacingd0 decreasing from 15.6 to 9.8 nm as DNA content was increased from 1% to 3%. A clear overlap concentration c* = 2% (w/v) was found. DNA hydrogel was found to be consisting of hierarchy of structures with well-defined length scales that contribute uniquely to the scattering process. In summary; it has been shown that hierarchical structure of DNA hydrogels is replete with several self-organized assemblies, whose spatial extension is dependent on the biopolymer concentration.

Acknowledgements

PKP acknowledges University Grants Commission, Government of India for a Junior Research Fellowship. KR is very thankful to Department of Science and Technology, Government of India-Inspire Faculty Award. This work was supported by the Department of Science and Technology (DST), India.

Conflict of interest

The author declares no conflict of interest.

References

  1. Horkay F, Basser PJ. Osmotic observations on chemically cross-linked DNA gels in physiological salt solutions. Biomacromolecules. 2004;5(1):232–237.
  2. Feiyue Xu, Lei Xu, Qi Wang, et al. A three-dimensional dynamic coculture system enabling facile cell separation for chondrogenesis of mesenchymal stem cells. Biochemical Engineering Journal. 2015;103:68–76.
  3. Jiang H, Pan V, Vivek S, et al. Programmable DNA hydrogels assembled from multidomain DNA strands. Chembiocheml. 2016;17(12):1156–1162.
  4. Papancea A, Valente AJ, Patachia S, et al. PVA-DNA cryogel membranes: characterization, swelling, and transport studies. Langmuir. 2008;24(1):273–279.
  5. Cho YJ, Kim HY, Huang H, et al. Spectroscopic characterization of interstrand carbinolamine cross-links formed in the 5‘-CpG-3‘ sequence by the acrolein-derived γ-OH-1,N 2-propano-2‘-deoxyguanosine DNA Adduct. J Am Chem Soc. 2005;127(50):17686–17696.
  6. Dutta S, Chowdhury G, Gates KS. Interstrand cross-links generated by abasic sites in duplex DNA. J Am Chem Soc. 2007;129(7):1852–1853.
  7. Liedl T, Dietz H, Yurke B, et al. A pH-triggered, controlled trapping and release of quantum dots in a DNA‐Switchable Hydrogel. Small. 2007;3(10):1688–1693.
  8. Xing Y, Cheng E, Yang Y, et al. Self‐assembled DNA hydrogels with designable thermal and enzymatic responsiveness. Adv Mater. 2011;23(9):1117–1121.
  9. Moran MC, Miguel MG, Lindman B. DNA gel particles: particle preparation and release characteristics. Langmuir. 2007;23(12):6478–6481.
  10. Lee CK, Shin SR, Lee SH, et al. DNA hydrogel fiber with self‐entanglement prepared by using an ionic liquid. Angew Chem Int Ed Engl. 2008;47(13):2470–2474.
  11. Orakdogen N, ErmanB, Okay O. Evidence of strain hardening in DNA gels. Macromolecules. 2010;43 (3):1530–1538.
  12. Mason TG, Dhople A, Wirtz D. Linear viscoelastic moduli of concentrated DNA solutions. Macromolecules. 1998;31(11):3600–3603.
  13. Murakami Y, Maeda M. DNA-responsive hydrogels that can shrink or swell. Biomacromolecules. 2005;6(6):2927–2929.
  14. Liu Y, Jun Y, Steinberg V. Longest relaxation times of double-stranded and single-stranded DNA. Macromolecules. 2007;40(6):2172–2176.
  15. Arfin N, Aswal VK, Kohlbrecher J, et al. Relaxation dynamics and structural changes in DNA soft gels. Polymer. 2015;65:175–182.
  16. Costa D, Hansson P, Schneider S, et al. Interaction between covalent DNA gels and a cationic surfactant. Biomacromolecules. 2006;7(4):1090–1095.
  17. Okay O. DNA Hydrogels: New functional soft materials. Journal Of Polymer Science Part B: Polymer Physics. 2011;49:551–555.
  18. He Y, Ye T, Su M, et al. Hierarchical self-assembly of DNA into symmetric supramolecular polyhedra. Nature. 2008;452(7184):198–201.
  19. Li C, Jones AF, Dun AR, et al. Rapid formation of a supramolecular polypeptide-dna hydrogel for in situ three‐dimensional multilayer bioprinting. Angew Chem Int Ed Engl. 2015;54(13):3957–3961.
  20. Kohlbrecher J, Wagner W. The new SANS instrument at the swiss spallation source SINQ. J Appl Cryst. 2000;33:804–806.
  21. Stauffer D, Aharony A. Introduction to Percolation Theory. 2nd ed. London, UK: Taylor and Francis; 1994.
  22. Stauffer D, Coniglio A, Adam M. Gelation and critical phenomena. In: Dušek K editor. polymer networks. Advances in Polymer Science. 1982;44:103–158.
  23. De Arcangelis L, Del Gado E, Coniglio A. Complex dynamics in gelling systems. Eur Phys J E. 2002;9(3):277–282.
  24. Debye P, Bueche AM. Scattering by an inhomogeneous solid. Journal of Applied Physics. 1949;20:518–525.
  25. Mohanty B, Bohidar HB. Internal pressure fluctuations in coacervates and syneresis. EPL (Europhysics Letters). 2006;76(5):965–971.
  26. Rawat K, Pathak J, Bohidar HB. Effect of persistence length on binding of DNA to polyions and overcharging of their intermolecular complexes in aqueous and in 1-methyl-3-octyl imidazolium chloride ionic liquid solutions. Phys Chem Chem Phys. 2013;15(29):12262–12273.
  27. Mrevlishvili GM, Svintradze DV. Complex between triple helix of collagen and double helix of DNA in aqueous solution. Int J Biol Macromol. 2005;35(5):243–245.
  28. Mohanty B, Aswal VK, Kohlbrecher J, et al. Length scale hierarchy in sol, gel, and coacervate phases of gelatin. J Polym Sci B Polym Phys. 2006;44(12):1653–1667.
  29. Simmons BA, Irvin GC, Agarwal V, et al. Small angle neutron scattering study of microstructural transitions in a surfactant-based gel mesophase. Langmuir. 2002;18(3):624–632.
Creative Commons Attribution License

©2017 Pandey, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.