Molecular interactions of DNA with transfectants: a study based on infrared spectroscopy and quantum chemistry as aids to ﬂ uorescence spectroscopy and dynamic light scattering analyses †

Cationic polymers are promising non-viral agents because of their ease of use, inexpensiveness and favorable safety pro ﬁ le. Being inherently cationic, they do spontaneously assemble in water with nucleic acids to give polyplexes. The more popular techniques and methods routinely applied for the characterization of such nanosupensions, unfortunately, do not a ﬀ ord direct observation of anion – cation interactions. We thus tackle the issue of probing at the molecular level the interplay occurring between the PO 4 (cid:1) groups of double stranded DNA and the NH 3+ groups of the simple and representative transfectant poly( L -lysine). By means of IR spectroscopy and Density Functional Theory (DFT) calculations we assigned the red shift of the anti-symmetric PO stretching wavenumber to the intermolecular water-mediated NH 3+ – PO 4 (cid:1) interaction. These changes observed in moving from ine ﬀ ective, low nitrogen (N) to DNA phosphate (P) ratios (N/Ps) to more e ﬀ ective conditions are in overall good agreement with changes in polyplex behavior unveiled by means of other analytical techniques such as dynamic light scattering, laser Doppler micro-electrophoresis, and ﬂ uorescence spectroscopy, but provide direct insight into the basic intermolecular DNA:transfectant interplay. The application of IR spectroscopic analysis to other non-viral gene delivery vectors may become an appealing approach to disclose molecular details on PO 4 (cid:1) – cationic polymer interaction, taking into account conformational constraints and steric hindrance e ﬀ ects that the other techniques actually overlook.


Introduction
The use of nucleic acids as drugs for the treatment of a plethora of genetic and acquired diseases that conventional medicine fails to cure is collectively known as gene therapy.Since nucleic acids by themselves cannot be efficiently internalized by cells, one of the major challenges that still stands in the way of the complete success of gene therapy is the development of effective and non-toxic delivery agents, namely viral and non-viral vectors. 1 Although genetically engineered viral vectors are at present the most efficient delivery systems, serious safety concerns have hampered their widespread use and prompted the scientic community to develop non-viral gene vectors. 2Since their rst description as gene vectors in 1987 by Wu and Wu, 3 cationic polymers are one of the most used classes of non-viral agents because of their easiness of use, inexpensiveness and favorable safety prole.As these gene delivery vectors are inherently cationic, they do spontaneously interact with polyanionic nucleic acids to form nano-to-micro-scaled complexes called polyplexes.Polymer/DNA complexes, in turn, have to be taken up by cells and unpacked intracellularly to release their nucleic acid cargo that undergoes transcription and ultimately translation (i.e.transfection). 4 Unfortunately, the in vivo application of polymeric gene vectors is rather limited because of their still low transfection efficiency caused by one or more bottlenecks among those mentioned here above.3][14] Unfortunately, mechanistic insights into polyplex formation and disassembly have been hampered by a lack of atomic resolution, and structural information on polymer/DNA complexes. 10On these premises, nding an effective gene delivery solution remains mostly a trial and error process. 15t is known that vibrational states depend on molecular structure and intermolecular interactions; hence, within this framework, IR spectroscopy might provide useful molecular information on DNA complexation that other techniques overlook.Furthermore, the assignment and interpretation of changes in vibrational spectra is nowadays fostered by the use of reliable quantum chemical tools, among which Density Functional Theory (DFT) framework provides the most effective balance between accuracy and computational cost. 16This allows considering sizable molecular models to tackle the theoretical interpretation of vibrational spectra of biologically relevant molecules, as we describe in this work.Instead of the overall view of polyplexes that many other techniques for characterizing colloidal dispersions allow for, the joint use of DFT and IR spectroscopy provides a molecular picture of the interactions occurring within polyplexes, which is experimentally probed by the spectral changes of a selected IR transition (anti-symmetric PO stretching) caused by the interaction of PO 4 À groups of the DNA with the NH 3 + groups of poly(L-lysine) (PLL).Interestingly we noticed that this interaction is modulated by the presence of bridging water molecules.PLL was chosen as model transfectant because it is the oldest and a structurally simple (i.e. with only primary amines, the 3-amino groups) cationic polymer for gene delivery and it has been largely investigated in the literature over the last almost three decades. 3,17n this study, we compared different dye-based and label-free techniques to evaluate the DNA-complexation behavior of PLL.Our results lead to the proposal of using IR spectroscopy for the investigation of DNA complexation by cationic polymers, possibly paving the way toward the application of this technique for the evaluation of polyplex stability.Furthermore, DFT calculations were shown very successful in describing the spectroscopic effects of the interaction between the DNA and cationic polymers within polyplexes.

In vitro transfection and cytotoxicity of PLL-based polyplexes
It is a matter of fact that the performance of a given gene delivery vector is a function of the nitrogen (N) to plasmid DNA phosphate (P) ratio (N/P) of use. 1,7,13Thus, we rst evaluated the transfection activity of PLL as a function of N/P.Transfection experiments were carried out on HeLa and SKOV-3 cell lines in culture medium supplemented with 10% FBS, commonly used in in vitro experiments to take into account the effect of serum on gene vectors prior to animal studies. 18In fact, it is now renown that non-viral gene vectors in biological uids such as plasma or serum, interact with proteins forming a dynamic protein shell called protein corona that can impair polyplex stability and oen adversely impacts their behaviour. 19,20As expected, transfection efficiency increased progressively with raising the N/P (Fig. 1).At N/P 0.5 the luciferase signal was below the threshold of sensitivity of the assay, while signicant transgene expression was observed only at N/P $ 1.5.The maximum transfection levels were obtained at N/P 3 (p < 0.05 vs. all lower N/Ps), beyond which a further increase in N/P resulted in a mere increase of the cytotoxicity levels (Fig. 1).Noteworthy, in all the transfection conditions tested, the cytotoxicity was very low in both the cell lines (<10%), thus not impairing transfection results.

Size and z-potential of PLL-based polyplexes
For efficient gene delivery it is fundamental that the polymers bind to and condense the extended structure of the DNA to form relatively tight complexes a few hundred nanometers in size and slightly positively charged on the surface.This is because aer interacting with negatively charged cell membrane components, particles smaller than 1 mm in diameter can be internalized by non-phagocytic cells via endocytosis. 1e thus evaluated the size (average hydrodynamic diameter, D H ) and the overall surface charge (z-potential, z P ) of PLL-based polyplexes by DLS and laser Doppler micro-electrophoresis, respectively, as a function of the N/P.These techniques unfortunately do not provide any information about PLL-DNA interactions at the very molecular level.
As shown in Fig. 2, z P curves as a function of N/P displayed a typical sigmoidal shape with the charge-inversion point (i.e.0 mV) at N/P z 1.25 at which a spike in the D H curve occurred. 21his neutralization or charge reversal phenomenon (also known as "overcharging" effect) results in an isoelectric point (IEP) and rapid aggregation of the colloidal uncharged polyplexes.
The slight displacement of charge-inversion point from the theoretical N/P of 1 to the experimental 1 < N/P < 1.5, suggests that more than a NH 3 + group on the polyelectrolyte chain is necessary to neutralize a single PO 4 À on the double-stranded DNA (dsDNA), and this can be ascribed to steric hindrance and conformational constraints of the PLL chain.Beyond the IEP, further addition of polycations leads to continuous polymer adsorption onto the polyplex surface, charge reversal, and restabilization of the particle suspension. 22Accordingly, we found that beyond N/P $ 2, a single unimodal population of polyplexes (polydispersity index z 0.2) was observed with a D H of z 100 nm and z P > +20 mV.Inversely, a reliable estimation of the diameter of uncomplexed plasmid DNA (pDNA) could not be determined by DLS, as previously indicated also by others. 23mplexation ability of PLL Fluorescence spectroscopy is a powerful technique with utility for studying the dynamic interactions and visualization of macromolecular complexes. 24Due to the excellent sensitivity of uorescence compared to other optical techniques and to the large variety of uorescent intercalating dyes available in the market, spectrouorimetric steady-state analyses of DNA-dyes are particularly suited to investigate the DNA-condensation ability of gene delivery vectors. 7,25We evaluated by uorophore-exclusion bio-analytical assay using SYBR Green I as DNA dye the ability of PLL to complex and condense nucleic acids as a function of N/P.Although very useful in providing fundamental information about the affinity between a polycation and the DNA, this technique does not allow direct identication of the interacting atoms and/or functionalities.
A dramatic increase in brightness (>1000-fold) is observed when the unsymmetrical cyanine dye SYBR Green I does intercalate the dsDNA helix. 24When dsDNA is condensed by a polymer, the intercalating dyes are hindered in their access to the intercalation sites.Thus, the extent of formation of polyplexes can be inferred from the relative uorescence of unassociated DNA aer exposure of the polymer/DNA complexes to SYBR Green I.
As shown in Fig. 3, the amount of free DNA decreased as the N/P increased.Specically, we observed that PLL can effectively complex the DNA already at very low N/P, with residual uorescence levels lower than 10% with respect to naked DNA already for N/P $ 1, in roughly good agreement with what we obtained by DLS analysis and laser Doppler micro-electrophoresis.It is worth noting that the maximum quenching of the uorescence signal (Fig. 3) was achieved immediately beyond the charge-neutralization point (i.e. at z P z 0 mV for N/P 1.25 and higher, Fig. 2) when the NH 3 + charges of the PLL did counterbalance and neutralize the PO 4 À on the dsDNA.

Analysis of phosphate stretching signals in the IR absorption of DNA
7][28] Upon effective complexation, which, depending on the investigation technique considered was found to happen for N/P $ 1, we observed a red shi of the anti-symmetric PO stretching of about 4 cm À1 .Once complexation has occurred no further PO peak red shi was observed.Based on the discussion above, we may deduce that the local environment nearby the PO 4 À of DNA has reached a stable structure with respect to the interacting species which is not affected any further by an increase of the N/ P. Consequently, we can infer that IR measurements on PLL/ DNA can prociently recognise DNA complexation states adequate for transfection (i.e.N/P $ 1) but cannot identify polyplexes with the highest efficiency (i.e.N/P $ 3).On the other hand, IR is a probe-free technique that yields atomistic information not obtainable by uorescent dye measurements.We have also carried out DFT calculations on selected model molecules to investigate in some details the vibrational dynamics of the PO 4 À groups and correlate the observed wavenumber shis to interactions with specic molecular entities.This is relevant for the investigation of the basic interactions that exist in polyplexes between the inner molecular layer of cationic polymers facing the nucleic acids.
As starting model we have considered a section of the double helix of DNA comprising four bases, whose structure optimized with DFT is given below (Fig. 5).This model has a remarkable size for current quantum chemical methods and common computational facilities (125 atoms, 1697 basis functions).The expensive part of the calculation is the evaluation of the second Fig. 4 IR spectra of a sample of salmon sperm DNA (0.7 mg mL À1 ) in water solution of PLL at different nitrogen (N) to plasmid DNA phosphate (P) ratios (N/Ps) (reported for each spectrum).The red-shift of about 4 cm À1 upon complexation is observed for N/P > 1.  derivatives of the energy with respect to nuclear displacements, which is necessary for obtaining the vibrational normal modes and the simulation of the IR spectrum.Hence, this DNA model was not practical for investigating the spectroscopic changes induced by PO 4 À -PLL interactions, but it could be protably used to dene smaller sub-models focused on the PO 4 À group.
In Fig. 6 we present the comparison between the experimental IR spectrum of the DNA in water and the simulated spectrum obtained from DFT calculations carried out on the model of Fig. 5.The overall agreement between simulation and experimental observation is remarkably good, aer taking into due consideration the fact that the DNA model is a signicant simplication of the real sample.In particular, all main features in the experimentally accessible range (1400-950 cm À1 ) nd a theoretical counterpart, aer accepting an inevitable wavenumber shi and relative intensity adjustment.
The model also provides a detailed description of the normal modes, which foster the understanding of the spectroscopic features associated to the PO 4 À groups.In particular two aims are reached: (i) the assignment of the PO symmetric/anti-symmetric stretching bands shows that the PO stretching content is actually shared among several normal modes in the expected spectral regions.This is the effect of vibrational coupling existing between PO 4 À and the nearby groups (see Fig. 6); (ii) based on the inspection of the coupling existing between phosphate and nearby groups (i) one can use the structure of the DNA backbone (taken from Fig. 5) to dene a smaller submodel aimed at describing the effects of interaction of the PO 4 À with cationic species of interest in transfection.
Following point (ii) we have considered a section of DNA backbone comprising two furanose rings connected by one PO 4 À group.This model is shown in Fig. 7 and labeled (a).To model the interaction effect of the cationic transfectant considered in our experiments (PLL) we have considered just the NH 3 + and a section of the connected alkyl chain (C 4 H 8 ), namely CH 3 -CH 2 -NH 3 + , and placed it in close vicinity of the PO 4 À group of the backbone (model (b) in Fig. 7).The effect of anion-cation interaction on the anti-symmetric PO stretching is dramatic, as the comparison provided in Fig. 7 clearly shows.Upon interaction a marked red shi of about 40 cm À1 is observed.
From the experimental point of view, the red shi of the antisymmetric phosphate stretching is more limited (see Fig. 7).This is explained by the fact that not all PO 4 À groups, even at a N/P 1.2, are actually able to strongly interact with the NH 3 + groups as in model (b).Hence the experimental spectrum is the result of the average signal produced by interacting and noninteracting phosphates, which reduces the observed red shi with respect to the theoretical expectations of model (b).The water environment typical of transfection scenarios requires considering also the explicit role of water molecules at the anion-cation interaction site.For this reason, we have modied model (b) including one or two explicit water molecules in mediating the anion-cation interaction (model (W 1 ) and (W 2 ), respectively).
The presence of water molecules screening the direct NH 3 + -PO 4 À interaction has a signicant effect on the red shi of the anti-symmetric phosphate stretching, which is progressively reduced along the sequence (a) > (W 1 ) > (W 2 ) (see Fig. 7).Finally, we have considered a further model (W 3 ) (not shown in Fig. 7, but reported in ESI †) that corresponds to model (W 2 ) with the inclusion of one additional water molecule hydrogen-bonded to the phosphodiester oxygen of the PO 4 À group.Consistently with the trend discussed above, its associated anti-symmetric stretching lies between (W 2 ) and (a) (Fig. SI1 †).

Conclusions
By the joint use of several experimental and theoretical methods, this proof-of-concept study demonstrates that IR spectroscopy can be protably used to investigate the DNA complexation by cationic polymers, allowing identifying transfection-active DNA complexation states by a red shi in the anti-symmetric (z1220 cm À1 ) PO stretching IR band.Noteworthy, IR spectra allow denite discrimination between ineffective DNA binding (see N/P < 1 in Fig. 4) and transfectionactive polyplexes (see N/P $ 1 in Fig. 4) as a function of the polycation concentration, as demonstrated by comparison with DNA complexation assay, size and z P measurements.Hence we can conclude that IR spectroscopy can be used to probe DNAcation interaction for applications in gene delivery development and in the large segment of the research involving processes where DNA interactions with biologically active molecules are very relevant events, such as in gene regulation studies.Furthermore we were able to predict by a simple DFT model the effect of anion-cation interaction on the IR spectrum, thus paving the way toward the use of this quantum mechanical tool for the investigation of the interactions occurring within polyplexes and for the prediction of their behaviour and stability in different environments.
Based on the results presented here above, application of IR to other transfection scenarios becomes an appealing approach to disclose molecular details on the basic DNA-cation interaction.This is expected to depend on the selected cationic species and on the capability of the cationic polymer to effectively bind the DNA, which in turn depends on the molecular exibility of the transfectant and its steric hindrance.Through the modulation of the anti-symmetric PO stretching wavenumber as a function of the N/P ratio one can monitor the combined effect of all these molecular properties, which are at the basis of transfection.

Materials
PLL hydrobromide (M w : 4-15 kDa, as determined by viscosity) was purchased from Sigma-Aldrich (Milan, Italy).BCA Protein Assay Kit was from Thermo Fisher Scientic (Rockford, IL, USA).The Luciferase Assay System and the pDNA encoding for the modied rey luciferase pGL3-Control Vector (5.2 kbp) were all from Promega (Milan, Italy).Salmon sperm DNA was purchased from Life Technologies Italia (Monza, Italy) at a concentration of 10 mg mL À1 .All the other chemicals, solvents, and culture reagents used were purchased from Sigma-Aldrich if not differently specied.

Preparation of polyplexes
Polyplexes were prepared at room temperature (r.t.) by adding pGL3 or salmon sperm DNA to a solution of 4-15 kDa PLL in deionized water (dH 2 O) at the desired polymer concentration to yield different N/Ps, and incubated for 20 min at r.t.before use.

SYBR Green I-exclusion assay
The ability of 4-15 kDa PLL to complex DNA was evaluated by uorophore-exclusion titration assay, as previously described. 13olyplexes were prepared in different volumes, depending on the N/P, using 400 ng of pGL3 and then diluted to 100 mL with 10 mM HEPES pH 7.0 containing 2Â SYBR Green I (l ex ¼ 497 nm; l em ¼ 520 nm).Fluorescence was measured using a microplate reader GENios Plus (Tecan, Segrate, Italy) and normalized over the uorescence of naked pDNA.

Measurement of size and surface charge of polyplexes
Size (D H ) analyses and surface charge (z P ) measurements of polyplexes were performed by means of DLS and laser Doppler micro-electrophoresis using a Malvern Zetasizer Nano ZS instrument (Malvern, UK) tted with a 633 nm laser at a xed scattering angle of 173 .Polyplexes (100 mL, 2 mg of pGL3) were prepared as described above and were diluted 1 : 4 in dH 2 O and le equilibrating at 25 C prior to measurements.

Cell culture and transfection
HeLa and SKOV-3 cells were cultured in humidied atmosphere of 5% CO 2 at 37 C. Cell culture medium was Dulbecco's Modied Eagle Medium (DMEM) containing 10 mM HEPES buffer, 1 mM sodium pyruvate, 100 U mL À1 penicillin, 0.1 mg mL À1 streptomycin, and 2 mM glutamine and supplemented with 10% (v/v) fetal bovine serum (FBS).The day before transfections, cells were plated in 96-well cell culture plates at a density of 2 Â 10 4 cells per cm 2 .102.4 ng per well of pGL3 were complexed as described above with 4-15 kDa PLL at the desired N/P and then added to cells in a nal volume of 100 mL per well of medium.Twenty-four hours later cell viability was assessed using AlamarBlue cell viability assay (Life Technologies Italia, Monza, Italy) according to manufacturer's guidelines.Cytotoxicity was dened as 100%viability [%].Results were expressed as % relative to the untreated controls.Aerwards, cells were lysed with Cell Culture Lysis Reagent (Promega, Milan, Italy) and luciferase expression was quantied by Luciferase Assay System and normalized to the protein content in each sample, as determined by bicinchoninic acid (BCA) assay.Results were expressed as relative light units (RLU) per mg of proteins.
Vibrational spectroscopy and quantum chemical calculations IR spectra have been collected with Nicolet Nexus FT-IR interferometer.Water solutions of salmon sperm DNA and PLL have been measured using sealed cells with CaF 2 windows and path length of 12 mm.Starting from the purchased salmon sperm DNA solution at 10 mg mL À1 , we have prepared more diluted solutions at 0.7 mg mL À1 using deionized water.This choice has been dictated by the trade-off between two contrasting aims: the will of more closely approach the low concentration regime adopted in transfection, while keeping the IR signal of salmon sperm DNA strong enough to allow reliable measurements.Omnic soware has been used for baseline subtraction and spectral averaging.
Quantum chemical calculations have been carried out with the Gaussian09 quantum chemistry computer code 29 on a small computer cluster operated at the Dept. of Chemistry, Materials and Chemical Engineering (Politecnico di Milano).We have considered a DFT method including dispersion forces, B97D/6-31+G(d,p), to account for p-p stacking interactions of the selected model DNA structure (Fig. 4).For models denoted in Fig. 6 as (a, b, W 1 , W 2 ) we have included IEFPCM mean eld correction for water solvation. 30This turns out to be essential for correctly modelling the interaction of the formal proton that is shared between the phosphate oxygen and the quaternary nitrogen.
By following the methodology described in the work by Radice et al., 31 we have analysed the Cartesian nuclear displacements associated to the normal modes of the selected molecular models in terms of valence internal coordinates, which directly relate the computed IR signals to vibrations of specic moieties in the chemical structure of the model molecule (such as, for instance, the stretching of a given PO bond).The output of DFT calculations typically reports normal modes information in terms of Cartesian nuclear displacements.This requires the application of a linear transformation (B) that converts Cartesian normal modes (L x ) into their internal coordinates counterpart (L R ¼ BL x ). 31 In this way, for each given k-th normal mode, it is possible to determine the contribution to the displacement of the j-th valence internal coordinate as the matrix element (L R ) jk .Further details on this procedure are given in ESI.† In order to unequivocally assign symmetric PO stretching and anti-symmetric PO stretching we have considered the free PO bonds (R 1 , R 2 ) of the two PO 4 À groups (A, B) of the DFT model reported in Fig. 4.This allows to dene symmetric (+) and anti-symmetric (À) combinations, i.e. (R 1 AE R 2 ), locally to A and B. For a given k-th normal mode, the associated content of say, symmetric/anti-symmetric PO stretching on PO 4 À A, will be given by |(L R ) ka AE (L R ) kb | where a and b denote the indices of the PO stretching coordinates R 1 and R 2 of PO 4 À A (respectively).A completely similar approach is taken for the coordinates of PO 4 À B.
The nal quantity, which is used in Fig. 5 to display the assignment of the computed IR features is given by the sum of the contributions from both PO 4 À groups, namely: |(L R ) ka AE (L R ) kb | A + |(L R ) ka AE (L R ) kb | B (where the + sign stands for the symmetric PO stretching and thesign stands for the antisymmetric PO stretching).Furthermore we also consider the other two PO stretching coordinates of both PO 4 À (we name them R 3 , R 4 ) that are directly connected to DNA backbone.Similarly to (R 1 , R 2 ) we dene in this case as well the symmetric and anti-symmetric combinations.It is possible to see that also the "backbone" PO bonds are participating to some extent to the symmetric PO stretching signal.Geometry optimizations of models (a) and (b) have been carefully carried out to ensure that the structure of the backbone was kept the same as in the DNA model, since our aim is focused on the normal modes of the phosphate group.Hence just the atoms of the phosphate group and counter cation were le free to move during geometry optimizations.It is worth to mention that the models (a and b) do not allow to investigate with the same accuracy the other vibrational features of DNA other than the anti-symmetric PO stretching modes.For the former a greater involvement of furanose vibrations is expected and actually found.This is due to the fact that models (a) and (b), for the sake of simplicity, do not consider any DNA base connected to the sugar rings.Actually, in both (a) and (b) model we replace one nitrogen atom of the two DNA bases (CA) with a hydrogen atom.This is also fully relaxed together with the connected carbon during geometry optimization.Geometry optimization of models (W 1 ), (W 2 ), (W 3 ) has been carried out similarly as models (a) and (b) with the additional optimization of the position of the explicit water molecules.

Statistical analysis
GraphPad version 5.04 (GraphPad soware, La Jolla, CA, USA) was used for statistical analysis.Comparisons among groups were carried out by one-way analysis of variance (ANOVA).Signicance was retained when p < 0.05.Data are expressed as mean AE standard error of the mean (SEM).

Fig. 1
Fig. 1 Transfection efficiency (white bars) and cytotoxicity (striped bars) of PLL-based polyplexes over a wide range of nitrogen (N) to plasmid DNA phosphate (P) ratios (N/Ps) and assayed on (a) HeLa and (b) SKOV-3 cell lines.Results are expressed as mean AE standard error of the mean (n ¼ 4).

Fig. 2
Fig. 2 Average hydrodynamic diameter (D H , grey asterisks and dotted lines) and z-potential (z P , full black squares and lines) of PLL-based polyplexes determined over a wide range of nitrogen (N) to plasmid DNA phosphate (P) ratios (N/Ps).Results are expressed as mean AE standard deviation (n ¼ 4).

Fig. 3
Fig. 3 Fluorophore-exclusion titration of PLL as a function of nitrogen (N) to plasmid DNA phosphate (P) ratio (N/P).DNA complexation ability of PLL was evaluated by monitoring the exclusion of SYBR Green I from polyplexes.Results are expressed as mean AE standard deviation (n ¼ 3).

Fig. 5
Fig.5The fully optimized equilibrium structure of the model DNA section considered for investigating the phosphate stretching vibrations.The model represents two short sections of interacting DNA strands.One strand is a cytosine-adenine (CA) sequence and the other is the complementary sequence (guanine-thymine -GT).F denotes furanose rings.‡

Fig. 7
Fig. 7 Simulated IR spectrum of (a) the non-interacting and (b) interacting DNA backbone model extracted from the molecular model represented in Fig.4(from DFT calculations, simulated FWHM ¼ 30 cm À1 ).To ease comparison with experiment a scaling factor of 1.04 has been applied to the computed wavenumber values.The experimental IR absorption spectra of a sample of salmon sperm DNA in water solution of PLL at N/P 0 (grey line) and N/P 1.2 (magenta line) are also reported for comparison.For simplicity of representation just a section of models W 1 and W 2 is represented, but the models considered in the calculations are similar to (a) and (b) on the backbone side.‡ The starting structure used as input for subsequent DFT calculations as been obtained using tools available from: http://casegroup.rutgers.edu.