Hannah K.
Murnen‡
a,
Adrianne M.
Rosales‡
a,
Andrey V.
Dobrynin
b,
Ronald N.
Zuckermann
c and
Rachel A.
Segalman
*a
aDepartment of Chemical and Biomolecular Engineering, University of California, Berkeley, CA 94720, USA. E-mail: segalman@berkeley.edu
bInstitute of Materials Science, University of Connecticut, Storrs, Connecticut 06269, USA
cMolecular Foundry, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
First published on 16th October 2012
The conformation of polyelectrolytes in aqueous salt solutions is closely related to their self-assembly properties. In particular, the persistence length has a large impact on how the chain can arrange itself. In this work, biomimetic poly N-substituted glycines (polypeptoids) have been designed to position charged side chains at precise distances from each other to elucidate the relationship between the spacing of the charges along the backbone, the ionic strength, and the persistence length. Using small angle neutron scattering (SANS), it is shown that at low ionic strength, polypeptoids with charged groups located closer to each other along the polymer backbone are stiffer than those with the charged groups spaced further apart. At high ionic strength, the total persistence length decreases for both macromolecules because the electrostatic repulsions between ionized groups are screened. The measured persistence lengths were compared to those calculated using a discrete chain model with bending rigidity, and it is shown that the electrostatic persistence length scales quadratically with the Debye screening length. It is also shown that the bare persistence length of a molecule with alternating ionizable and hydrophilic groups is larger than that of a molecule containing 100% ionizable groups. This difference can be attributed to the longer hydrophilic side chains that may induce local chain stiffening.
One of the most common measures of polyelectrolyte chain conformation is the persistence length. Persistence length impacts a wide range of polymer properties, such as intrinsic viscosity,15 electrical conductivity,16 and in concentrated solutions or in the bulk, ion clustering.17 Thus, understanding chain conformation of these charged polymers lends insight into their functional properties. For example, Wang et al. measured the persistence length of a water-soluble conjugated polymer in semi-dilute solution to determine the effect of chain conformation on photoluminescence and found the presence of large aggregates.18,19 Upon the addition of a cationic surfactant, the electrostatic interactions are screened and the aggregates decrease in size, leading to a remarkable enhancement in photoluminescent efficiency. Further control over chain shape and charge spacing is therefore desirable to design highly functional materials. Recent research on charged polymers in the solid state has examined the influence of polymer architecture and monomer sequence on the spacing of ionic clusters.20,21 However, these parameters are not easily probed in synthetic polyelectrolytes, and to this end, biologically relevant polymers have proved useful, as they contain sequence specificity and defined secondary structures.
Experimental studies of biologically relevant polyelectrolytes in water have mostly focused on DNA. The effect of ionic strength on the persistence length of DNA has been studied using a variety of techniques, such as light scattering,22,23 force-extension24 experiments, dielectric spectroscopy,25,26 and fluorescence microscopy.27 However, the measurement of the persistence length can be convoluted by these techniques. For example, at high and moderate salt concentrations, a persistence length obtained by light scattering measurements for sufficiently long macromolecules can include an excluded volume contribution. In addition, force-extension experiments can underestimate values of the persistence length due to enthalpic effects associated with stretching the molecules at forces above 5 pN. Because of these issues, the relationship between ionic strength and DNA persistence length is difficult to quantify, and experimental observations often contradict each other. Various techniques have also been used to measure the ionic strength dependence of the persistence length for other biopolymers.28–31 In particular, light scattering has been used to determine a salt concentration dependence of the persistence length of hyaluronan32 and a polysaccharide,33 but aggregation can complicate measurements at low ionic strength, as seen in a neutron scattering study by Bonnet et al.34 Similar issues lead to very few studies for polypeptides.35,36
Theoretical studies on the chain conformation of polyelectrolytes also lead to much debate. The total persistence length for a polyelectrolyte contains contributions from both the bare persistence length of the chain and the electrostatic persistence length. Much of the debate on developing theoretical models for polyelectrolyte systems has centered on the flexibility of the underlying chain and how it impacts the scaling of the electrostatic persistence length with regards to ionic strength.37–39 On the one hand, theory originating with Kuhn, Kunzle, and Katchalsky40 suggests that persistence length scales with the inverse of the square root of the salt concentration (i.e., a linear dependence on the Debye screening length); this calculation corresponds to a flexible or semi-flexible chain with electrostatic corrections. On the other hand, OSF1,2,41–43 theory shows the persistence length scales with the inverse of the salt concentration (i.e., a quadratic dependence on the Debye screening length); OSF theories assume a locally stiff chain for which electrostatic interactions can modify the bending properties. In both cases, theory has been developed over the past several decades, and neither fully accounts for experimental observations. Given the difficulties in obtaining a definitive result for the verification of scaling laws, it is useful to investigate a simplified system in which factors such as chain length and charge distribution along the polymer backbone can be controlled with high accuracy.
The simplified system we will therefore focus on is polypeptoids, which are sequence-specific poly(N-substituted glycines).44 Biological polyelectrolytes have a sequence-specific arrangement of the charged groups along their backbone, and the sequence of charged groups is essential to the structure and function of the macromolecule. However, natural biopolymers also have many other interactions that impact their persistence length. For instance, polypeptides are inherently chiral and form intramolecular and intermolecular hydrogen bonds. They can also have strong hydrophobic interactions that drive a chain to collapse. Polypeptoids, however, do not have these competing interactions due in large part to their backbone chemistry. Polypeptoids have the same backbone as polypeptides, but the side chain is attached to the nitrogen rather than the backbone α-carbon. This N-substitution eliminates chirality and intra and intermolecular hydrogen bonding in the backbone, allowing the interactions in the system to be tuned by the introduction of different side chains. Furthermore, these materials are synthesized using a solid phase submonomer process that leads to sequence specific, monodisperse chains. In addition to serving as model biomimetic polymers, polypeptoids are an interesting self-assembly material in their own right.45 In solution, they can assemble into supramolecular helices46 and nanometer-thick sheets.47–49 To model and understand these assemblies, it is important to have a basic understanding of polypeptoid behavior with regard to ionic interactions.
In this paper, we explore the impact of charge density on chain stiffness using polypeptoids with either 100% carboxyl groups or 50% carboxyl groups (in an alternating sequence). The persistence lengths of these molecules have been measured using small angle neutron scattering (SANS). The results show that at low ionic strength, a polypeptoid in which every monomer contains an ionizable group has a higher persistence length than a polypeptoid in which ionizable and hydrophilic monomers alternate. At high ionic strength, this difference is far smaller, and both polypeptoids show a decrease in persistence length. The data has been analyzed using a semi-flexible chain model, and it is shown that the persistence length has a quadratic dependence on the Debye screening length. It is also shown that the bare bending rigidity is higher for the macromolecule with alternating charged group distribution. This result is attributed to the effect of the longer hydrophilic side chains that can induce local chain stiffening.
Polymer | Abbreviation | Purity | Obs/actual MW | Structure |
---|---|---|---|---|
Acetyl-poly(N-carboxyethyl glycine-N-methoxyethyl glycine)18 | Ac-p(NmeNce)18 | 98% | 4487.1/4488.2 |
![]() |
Acetyl-poly(N-carboxyethyl glycine)36 | Ac-p(Nce)36 | 97% | 4737.0/4737.6 |
![]() |
The experiments at NIST were conducted using a neutron wavelength of 6 Å and two sample-to-detector distances: 1.3 m and 4 m, which corresponded to a q-range of approximately 0.003–0.5 Å−1. Raw intensity measurements were corrected for the detector efficiency and the dark current, and absolute intensities (units of cm−1) were calculated using a direct beam measurement from an empty cell. Quartz banjo cells (Hellma USA, Plainview, NY) with a path length of 1 mm were used in a temperature controlled multiple position sample holder. The data were reduced using the NCNR SANS reduction macros52 in Igor Pro. For each sample, data from each q range was merged, and the scattering contribution from the solvent (deuterated water) was subtracted. Samples containing 1× and 3× charging levels were measured at NIST.
The experiments at ORNL were also conducted using a neutron wavelength of 6 Å and two sample-to-detector distances: 1.7 m and 14.5 m, which corresponded to a q-range of approximately 0.003–0.5 Å−1. Raw intensity measurements were corrected for the detector efficiency and the dark current, and absolute intensities (units of cm−1) were calculated using a direct beam measurement from an empty cell. Quartz banjo cells (Hellma USA, Plainview, NY) with a path length of 2 mm were used in a temperature controlled multiple position sample holder. The data were reduced using the Spice SANS reduction program in Igor Pro. For each sample, data from each q range was merged, and the scattering contribution from the solvent (deuterated water) was subtracted. Samples containing 10× and 30× charging levels were measured at ORNL.
The scattering intensity of the reduced sample data was fit to a semiflexible cylinder model with excluded volume using the Igor data analysis package provided by NCNR at NIST. In its simplest form, the scattering intensity can be represented as
I(q) = KP(q) + Iincoherent | (1) |
Sample | Polymer | NaOH equiv. | [NaOH] (M) | Overall fraction charged monomers, α | Measured Lp (nm) |
---|---|---|---|---|---|
1 | Ac-p(Nce)36 | 1 | 0.04 | 0.48 | 0.7 ± 0.2 |
2 | Ac-p(Nce)36 | 3 | 0.15 | 1 | 0.9 ± 0.3 |
3 | Ac-p(Nce)36 | 10 | 0.68 | 1 | 0.41 ± 0.06 |
4 | Ac-p(Nce)36 | 30 | 2.19 | 1 | 0.4 ± 0.2 |
5 | Ac-p(NmeNce)18 | 1 | 0.03 | 0.18 | 0.54 ± 0.02 |
6 | Ac-p(NmeNce)18 | 3 | 0.08 | 0.50 | 0.62 ± 0.2 |
7 | Ac-p(NmeNce)18 | 10 | 0.36 | 0.50 | 0.57 ± 0.03 |
8 | Ac-p(NmeNce)18 | 30 | 1.15 | 0.50 | 0.45 ± 0.03 |
The ionization of the carboxylic acid side chains is pH-dependent, and a titration with sodium hydroxide shows the relative charge states of each polymer. Interestingly, it is easier to deprotonate Ac-p(Nce)36 than Ac-p(NmeNce)18 as seen in Fig. 1. Fewer equivalents of NaOH were required to reach the equivalence point and eventually full deprotonation for Ac-p(Nce)36. This result may be because when the ionizable groups are directly next to each other, they can more easily share sodium ions, allowing the charges to be more evenly distributed. However, this may also be a kinetic effect, as both molecules tend to form aggregates over long times. To further examine this trend, additional molecules with carboxyl side chains spaced every third and every fourth monomer apart (Ac-p(Nme2Nce)12 and Ac-p(Nme3Nce)9) were studied and found to require the same amount of NaOH as Ac-p(NmeNce)18 if sufficient time was allowed for the system to stabilize (data shown in ESI, Fig. S1†). All solutions studied were made and filtered immediately prior to SANS experiments to avoid aggregation.
![]() | ||
Fig. 1 The titration curves of Ac-p(Nce)36 (open green squares) and Ac-p(NmeNce)18 (solid blue squares) as pH versus equivalents of NaOH. |
Solutions were made with each polymer and various concentrations of NaOH (Table 2). When the solutions contained 1 equivalent of NaOH, the polymers were only partially charged. From the titration curve equivalence point, it is calculated that the acid groups on Ac-p(Nce)36 (sample 1) are 48% charged while those on Ac-p(NmeNce)18 (sample 5) are 35% charged. Because only 50% of the monomers in Ac-p(NmeNce)18 are ionizable, the overall fraction of monomers that are charged, α, is actually 18%. The subsequent solutions were all made using higher concentrations of sodium hydroxide (3×, 10×, 30×), where all of the carboxylic acids are negatively charged. In these solutions, adding more sodium hydroxide simply increases the ionic strength and provides screening of the electrostatic repulsion between ionizable groups.
In order to obtain information about the conformation of these molecules at the various ionic strengths, small angle neutron scattering (SANS) experiments were performed. As a first assessment of the polymer conformation, the SANS data was plotted on a Kratky plot (Q2I versus Q), as shown in Fig. 2. The shape of this plot provides information about the polymer and its arrangement in solution. Fig. 2a shows the Kratky plot for Ac-p(Nce)36 at various equivalents of NaOH. The data for both the 1× and 3× solutions increase monotonically, indicating that the molecule is dispersed in solution. However, for both the 10× and 30× solutions, there is a distinct difference in the shape of the data curves. The data have a hump with a peak around 0.1–0.15 Å−1. This hump indicates that the macromolecule has a more coiled configuration in solution. This is likely due to the fact that when the amount of sodium in the solution is relatively low (in the 1× and 3× case), the ionized groups cause the macromolecule to remain extended. However, when more NaOH is added (in the case of 10× and 30× solutions), the macromolecule can also form complexes between the positive sodium and the negative carboxylic acids. This is known as the counterion condensation effect,4 and it pulls the molecule into a more coiled state. A similar trend is seen for Ac-p(NmeNce)18. The decrease in chain size due to counterion condensation has been observed experimentally54,55 and in computer simulations.56,57 While radius of gyration (Rg) measurements would be helpful in confirming this analysis, low signal to noise ratio in the small q range (Guinier region) of the scattering data makes it difficult to obtain reliable Rg values. In addition, the Rg measurement may be convoluted by any complexation of the polypeptoid with sodium ions.
![]() | ||
Fig. 2 Kratky plots for each polymer solution. (a) Includes the traces for Ac-p(Nce)36 while (b) includes the traces for Ac-p(NmeNce)18. For both polymers the solutions with 10× and 30× NaOH equivalents show peaks while the solutions with 1× and 3× NaOH equivalents monotonically increase. |
In order to understand the effect of electrostatic interactions on chain configurations, a semi-flexible chain model53 was used to fit the SANS data and to obtain a value of persistence length, Lp, for each ionic strength. These data are summarized in Table 2. In agreement with previous measurements of polypeptoid persistence lengths,58 these macromolecules are quite flexible with persistence lengths ranging from 0.4 nm up to approximately 1 nm. Dependence of the persistence length on NaOH concentration is shown in Fig. 3. The chain persistence length first increases with increasing NaOH concentration, then it begins to decrease as more NaOH is added. The initial increase is due to the NaOH behaving as a titrant, increasing the fraction of ionized groups along the polymer backbone. Electrostatic repulsion between these ionized groups stiffens the polymer chain, resulting in an increase in persistence length. As more sodium hydroxide is added, the persistence length decreases as expected because the added salt ions screen electrostatic interactions between ionized groups along the polymer backbone, making the macromolecules more flexible. This is in line with the qualitative results from the Kratky plots, where it was seen that the molecules adopted more coiled conformations (and therefore were more flexible) at higher ionic strengths. There is also a difference between the two polymers, particularly at low NaOH concentrations. The Ac-p(Nce)36 has a larger Lp at these concentrations, which is presumably due to the close proximity of the charged groups along the backbone and the larger overall amount of charge as compared to Ac-p(NmeNce)18. These two factors both contribute to increase the persistence length for Ac-p(Nce)36, demonstrating that the amount of charge has the ability to influence the polymer conformation in solution.
![]() | ||
Fig. 3 The persistence lengths for Ac-p(Nce)36 (open green squares) and Ac-p(NmeNce)18 (closed blue squares). The red dashed line is the best fit line using the expression Lp = f + g[NaOH]−0.5. |
As discussed in the introduction, the most common treatment for the theoretical scaling of persistence length with salt for polyelectrolytes is the work of Odijk1,41–43 and Skolnick and Fixman,59 who showed that for locally stiff polyelectrolytes, the persistence length should scale directly with the inverse of the salt concentration. This relationship holds well for low salt concentrations and for relatively stiff polyelectrolytes. However, at high salt concentrations and for more flexible polyelectrolytes, there is significant deviation from the OSF model. Several experimental60 and theoretical61,62 studies have shown that the persistence length of flexible polyelectrolytes is actually inversely proportional to the square root of the salt concentration, or chain persistence length scales linearly with the Debye screening length. These results hold especially well for polyelectrolytes that do not assume locally stiff conformations. Given the very low persistence length values measured here, it is likely that the polyelectrolytes described in this paper do not adopt locally stiff conformations. Therefore, it is possible to examine the scaling dependence of the persistence length on the ionic strength by fitting experimental data to the following equation
Lp = f + g[NaOH]−0.5 | (2) |
Below we will account for the finite chain length and map our polymers into a discrete chain model with internal bending rigidity.63 In the framework of this model,63 chain properties are described by two independent parameters: effective chain bond length b and chain bending constant K. For effective bond length, we will use a projection length of a bond in zig–zag conformation, which gives a value of b = 0.37 nm. In the case of charged polymers, a chain bending constant K has contributions from the bare chain bending rigidity K0 and from the electrostatic interactions between ionized groups along the polymer backbone (details shown in the ESI†):64
![]() | (3) |
![]() | (4) |
This equation was used to fit the measured persistence length as a function of NaOH concentration as shown in Fig. 4 by considering K0 in eqn (3) as fitting parameter. For both peptoids, the agreement between experimental and calculated values is very good. The biggest difference is observed for values of the persistence length at the highest ionic strength. As shown in Table 3, the value of the bare bending constant K0 is larger for the Ac-p(NmeNce)18 peptoid, which could be due to the longer side chain causing local chain stiffening. In Fig. 4, we plot the electrostatic contribution to the chain bending constant. All of the data collapse onto a universal line, confirming the quadratic dependence of the chain persistence length on the Debye screening length and illustrating the high quality of the fitting procedure.
![]() | ||
Fig. 4 Electrostatic contribution to the chain persistence length as a function of the reduced Debye screening length. Experimental data for Ac-p(NmeNce)18 are shown as open circles and those for Ac-p(Nce)36) shown as filled circles. Dashed red line corresponds to eqn (3). The solid black line illustrating a quadratic scaling serves as a guide to the eye. |
Sample | Polymer, NaOH equivalent | K 0 | K | L p,calc (nm) | L p,fit (nm) | Radius, Rfit (nm) |
---|---|---|---|---|---|---|
1 | Ac-p(Nce)36, 1× | 0.975 | 1.92 | 0.71 | 0.7 ± 0.2 | 0.65 ± 0.03 |
2 | Ac-p(Nce)36, 3× | 0.975 | 2.35 | 0.87 | 0.9 ± 0.3 | 0.45 ± 0.01 |
3 | Ac-p(Nce)36, 10× | 0.975 | 1.22 | 0.45 | 0.41 ± 0.06 | 0.98 ± 0.04 |
4 | Ac-p(Nce)36, 30× | 0.975 | 1.03 | 0.38 | 0.40 ± 0.2 | 0.46 ± 0.02 |
5 | Ac-p(NmeNce)18, 1× | 1.43 | 1.41 | 0.52 | 0.54 ± 0.02 | 0.54 ± 0.01 |
6 | Ac-p(NmeNce)18, 3× | 1.43 | 1.86 | 0.69 | 0.62 ± 0.2 | 0.60 ± 0.02 |
7 | Ac-p(NmeNce)18, 10× | 1.43 | 1.38 | 0.51 | 0.57 ± 0.03 | 0.69 ± 0.01 |
8 | Ac-p(NmeNce)18, 30× | 1.43 | 1.27 | 0.47 | 0.45 ± 0.03 | 0.69 ± 0.01 |
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c2sm26849c |
‡ Authors contributed equally. |
This journal is © The Royal Society of Chemistry 2013 |