Patricia
Montes
a,
Tania
Chopra
b,
Rafał
Konefał
b,
Pavla
Hájovská
c,
Igor
Lacík
c,
Vladimír
Raus
b,
Miroslav
Šlouf
b,
Mariusz
Uchman
a and
Miroslav
Štěpánek
*a
aDepartment of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, Hlavova 2030, 128 40, Prague 2, Czech Republic. E-mail: stepanek@natur.cuni.cz
bInstitute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovského nám. 2, 162 06, Prague 6, Czech Republic
cDepartment for Biomaterials Research, Polymer Institute of the Slovak Academy of Sciences, Dúbravská cesta 9, 845 41, Bratislava, Slovakia
First published on 18th November 2024
Biguanide-based cationic polyelectrolytes are used as key components of interpolyelectrolyte complexes bolstering alginate hydrogel microcapsules employed in cell therapies. Nevertheless, electrostatic complexation of these unique polycations has not been studied before. In this study, the interaction between biguanide condensates and anionic polyelectrolytes with carboxylate groups was studied on a model system of a metformin condensate (MFC) and an anionic diblock polyelectrolyte poly(methacrylic acid)-block-poly(ethylene oxide) (PMAA–PEO). The formation of MFC/PMAA–PEO core–corona nanoparticles was followed by static, dynamic and electrophoretic light scattering and by isothermal titration calorimetry and their internal structure was investigated by small angle neutron scattering and cryogenic transmission electron microscopy. It was found that the aggregation of PMAA–PEO chains induced by MFC occurred at much lower MFC/PMAA–PEO ratios that would correspond to the isoelectric point, thus yielding strongly negatively charged nanoparticles, suggesting the role of specific (non-electrostatic) interactions in the stabilization of the complex between PMAA and MFC.
It is generally assumed that properties of PEs in solution, including their electrostatic complexation, are dictated mostly by Coulombic interactions between charges on the PE backbones and on small ions (counterions, co-ions) in the solution. This means that two PEs differing only in charged groups should have the same properties as long as the charges of the groups and their mutual distances remain the same. In this approach, the solvent is described as a continuous medium with given dielectric permittivity, and all effects associated with the solvation of polymer chains and small ions in the solution are neglected.12 While it is known that the behaviour of IPECs can be affected by the PE chain length,13 by modifications of the PE backbone substitution (e.g., by hydrophobic/hydrophilic side groups),13,14 or, because of specific solvation effects, by the choice of counterions,15 less attention has been paid so far to the effect of non-Coulombic interactions of charged groups of the PE on the IPEC formation.
A few recent studies stressed the difference in the behavior of IPECs formed by a PE with ammonium groups and by a PE with guanidinium (Gnd+) groups.16,17 From the structural point of view, Gnd+ is an intriguing chemical species belonging to the class of planar Y-conjugated quasi-aromatic structures.18 The ion is partly hydrophobic as water molecules can form hydrogen bonds only from the sides of the planar ion. Therefore Gnd+ sticks to hydrophobic surfaces or proteins and forms like-charge contact ion pairs despite the obvious electrostatic repulsion. It was found that Gnd+-based cationic PE exhibited much stronger interaction with PEs than the ammonium-based ones, affording solid IPECs rather than a liquid coacervate phase as in the case of ammonium and showing much higher stability at increased ionic strength of the solution.17 A quasilelastic neutron scattering study19 showed that chain dynamics in IPECs with Gnd+-based cationic PE did not change with ionic strength of the solution so that the interactions between oppositely charged PE chains were less prone to electrostatic screening, indicating major contribution from specific (non-electrostatic) interactions to the stability of the complex.
Recently, the technological significance of guanidine-based PEs has been increasingly recognized, with these materials exploited in a number of different applications including potent antimicrobial polymers,20 high modulus supramolecular hydrogels,21 drug delivery vehicles,22 ultra-stable polyelectrolyte multilayers,23 or wastewater treatment agents.24 One of these applications is in the “next generation” treatment of type 1 diabetes via the transplantation of pancreatic cells.25–27 These cells are immunoprotected by their encapsulation in multicomponent alginate microcapsules that are reinforced by the IPEC formed using the polycationic poly(methylene-co-cyanoguanidine) (PMCG). The PMCG repeat unit is based on a biguanide moiety that can be seen as two condensed Gnd groups sharing one of their nitrogens. While the biguanide motif is formed in situ via a post-polymerization modification reaction during the PMCG synthesis, we proposed that similar biguanide-based polycations should be accessible also through an alternative route exploiting the direct condensation of N-substituted biguanides with formaldehyde. Conveniently, certain substituted biguanides are commercially available, such as 1,1-dimethylbiguanide (commonly known as metformin), which is widely used as a medication for the type 2 diabetes treatment.28
In this study, we synthesized and characterized a novel type of cationic biguanide condensate based on metformin (further referred to as metformin-formaldehyde condensate, MFC, Scheme 1) and investigated its interaction with a double hydrophilic diblock anionic polyelectrolyte poly(methacrylic acid)-b-poly(ethylene oxide) (PMAA–PEO) in basic aqueous solution, ensuring full ionization of the PMAA block. While previous studies focused on guanidine-based cationic PEs containing Gnd+ in polymer side groups,17,19 biguanide condensates represent a structurally different PE with guanidine motifs not only in the side groups but also in the backbone, whose ability to form IPECs has not yet been investigated. Using scattering techniques, ITC and Cryo-TEM, we followed the formation of the IPEC complex of PMAA–PEO and MFC resulting in their co-assembly within NPs composed of the PEO corona and the MFC/PMAA core.
Scheme 1 Representative chemical structure of the metformin-formaldehyde condensate (MFC). Note that the charge can be delocalized throughout the biguanide moiety. |
The g(2)(q,τ) functions were fitted using an inverse Laplace transform by means of a constrained regularization algorithm (CONTIN) and the obtained distribution functions of relaxation times A(τR) were recalculated to the distribution functions of apparent hydrodynamic radii using the relation RHapp = kTτRq2/6πη, where k is the Boltzmann constant, T temperature and η is solvent viscosity. The g(2)(q,τ) functions exhibiting monomodal RHapp distributions were further fitted using the 2nd cumulant expansion, g(2)(q,τ) = 1 + βexp[−2Γ1(q)τ + Γ2(q)τ2], where β is the coherence factor Γ1(q) and Γ2(q) are the 1st and the 2nd cumulant.
The front detector was positioned at 1.4 m from the sample at an angle of 20°, the rear detector was at 17.6 m or 5.6 m, the beam was collimated at 17.6 or 5.6 m with a source aperture of 40 × 55 mm2 and a sample aperture of 7 × 10 (h × v) mm2. A wavelength of λ = 0.6 nm (relative full width at half maximum 10%) was selected, allowing to cover a continuous q-range of 0.025–6.5 nm−1 for 17.6 m and 0.085–6.5 nm−1 for 5.6 m, where q is the magnitude of the scattering vector, q = 4π/λsin(θ/2), θ being the scattering angle. Samples were kept in quartz cuvettes (type 100-QS, Hellma GmbH, Müllheim, Germany) of 2 mm pathway, on a thermalized rack at 20 °C. Data were reduced with the program Grasp V.10.25f30 normalizing with monitor, subtracting the contribution from the empty cell, taking into account noise from the measurement with a sintered 10B4C piece at the sample position, and using for transmission the intensity transmitted by a semi-transparent beam stop. Parallax from the detector and from the sample attenuation were corrected for, and a flat field was used. The absolute scale was obtained from the measurement of the direct beam with a calibrated attenuator (flux of 1.64 MHz for 17.6 m and 13.8 MHz for 5.6 m). After reduction, the scattering from the buffer was subtracted.
The solution behaviour of MFC was studied by SANS (0.1–4 nm−1) and SLS (11–25 μm−1); Both SANS and SLS scattering curves (Fig. 1) were fitted with the Guinier equation, I(q) = I(0)exp(−Rg2q2/3) + Ib, where Ib is the background from the incoherent scattering in SANS; SLS data were fitted with the set value Ib = 0. The fits provided gyration radii, RgSANS = 0.6 nm and RgSLS = 24 nm.
The presence of two types of scatterers differing in the size was confirmed by DLS measurement which exhibited a bimodal distribution of relaxation times, τR (Fig. 2a). The mean relaxation rates of both modes, Γfast = 1/〈τR〉fast and Γslow = 1/〈τR〉slow are shown in Fig. 2b as functions of q2. The fast mode was purely diffusive, Γfast(q) = Dfastq2, with the corresponding apparent diffusion coefficient, Dfast = Γfast(q)/q2 = 4.1 × 10−10 m2 s−1, and the corresponding hydrodynamic radius of RH = kT/6πηDfast = 0.6 nm (k being the Boltzmann constant, T temperature and η solvent viscosity), equal to Rg obtained from the SANS data. This confirmed that the fast mode can be attributed to the diffusion of free MFC chains.
The slow mode (reflected also by the Guinier regime observed in static light scattering) was most likely caused by MFC dynamics slowed by attractive interactions between polyelectrolyte chains. A similar slow mode has frequently been observed in DLS from aqueous solutions of polyelectrolytes and the origin of this behaviour has been extensively debated,31–33 with its presence explained as diffusive motion of metastable domains of concentrated polyelectrolyte solution, stabilized by shared clouds of counterions. In the case of MFC, the behaviour in DLS was more complex as the dependence of Γslow on q2 was not linear, suggesting a contribution of non-diffusive fluctuations to the slow relaxation (most likely internal dynamics in the domains). We would like to note that Cryo-TEM measurement (Fig. S11a, ESI†) was not able to reveal any inhomogeneities in the MFC solution, indicating that MFC concentration in the domains was only slightly higher than its average concentration in the solution, so the formed aggregates were very loose.
A similar behaviour, that is, a bimodal distribution of relaxation times, was found for DLS of PMAA–PEO aqueous solution (Fig. 3a, curve 0; τR scale is recalculated to apparent hydrodynamic radii RHapp = kTq2τR/6πη). The mean RH for the fast mode, corresponding to the diffusion of free PMAA–PEO chains, was 10 nm.
Mixing PMAA–PEO with MFC led to the IPEC formation, which, when followed by DLS (Fig. 3a), was manifested by the increase in the amplitude and mean RHapp of the fast mode, indicating the aggregation of PMAA–PEO with MFC while the slow mode of PMAA–PEO gradually disappeared. For w > 0.6, the RHapp distribution became monomodal, dominated by the IPEC particles of PMAA–PEO and MFC. The transformation of loose aggregates to small compact NPs was apparent also from the static light scattering data (Fig. 3b) where the loose aggregates manifested themselves by the steeply decreasing scattering intensity with q in the low region up to w = 0.6. For w > 0.6, scattering became nearly isotropic and its intensity increased by an order of magnitude, indicating that the formed NPs were small but had high molar mass.
Gyration radius, Rg, and hydrodynamic radius, RH (Fig. 4) of the NPs formed at MFC-to-PMAA–PEO ratios above 0.6 were evaluated from SLS (from the Zimm equation, I(0)/I(q) = 1 + Rg2q2/3) and DLS (from the diffusion coefficient obtained as linear extrapolation of Γ1(q)/q2 to q = 0). The concomitant increase in Rg and decrease in RH with w could be explained by different contributions to Rg and to RH: The hydrodynamic radius reflected rather the overall radius of the NPs which decreased both with the collapse of PMAA chains and with the decrease in the aggregation number at higher w. On the contrary, dominant contributions to the gyration radius originated from the NP cores, thus growing as the amount of MFC (with higher dn/dc than that of both PMAA and PEO) in the core increased.
Fig. 4 Gyration radius, Rg (curve 1), and hydrodynamic radius, RH (curve 2) of MFC/PMAA–PEO particles, as functions of w. |
Scattering intensity at θ = 90°, I(90°), is shown as a function of w in Fig. 5, curve 1. It exhibited a maximum at w = 0.8, indicating that the stability of the formed particles reached its maximum at this composition; the further addition of MFC led to a decrease in aggregation number.
It is worth mentioning that the stoichiometric ratio between the MFC and the PMAA monomeric units would correspond to the mass ratio w = 0.75, assuming that 25% of MFC monomeric units bear CH2OH groups (see Scheme 1) as would correspond to the used stoichiometric ratio between metformin and CH2O monomeric units of paraformaldehyde, 1:1.25, used for the MFC synthesis. It was thus much lower than the mass ratio corresponding to the isoelectric point of the complex, which was reached at w = ∼2 (Fig. 5, curve 2). Moreover, the zeta potential of MFC/PMAA–PEO was decreasing with the increasing content of MFC up to w = 0.6. Similar, albeit less pronounced behaviour, was observed during DNA condensation by cationic surfactants,34 where the conformational changes of the DNA chain induced by the bound surfactant ions exposed charges which were initially screened, thus making the zeta potential more negative. In the case of the MFC/PMAA–PEO system, the aggregation of PMAA chains with MFC was triggered while the net charge of the complex was still strongly negative, and the increased electrophoretic mobility of MFC/PMAA–PEO particles in comparison with loose aggregates of PMAA–PEO chains caused the increase in the negative zeta potential.
To check whether the compact particles formed by the MFC/PMAA–PEO complex had a core/corona structure with the core formed by the MFC/PMAA complex and the corona of PEO chains, we conducted a SANS measurement for the mass ratio, w = 4. The scattering curve (Fig. 6) fitted well to the Pedersen–Gerstenberg (PG) model of the homogeneous sphere surrounded by Gaussian chains35 with an additional term for incoherent scattering. The fit, performed with scattering lengths of the core block (that is, a complex of a PMAA chain with MFC chains), bsph, and the corona block, bchain, provided the following results: sphere radius, rsph = 18 ± 3 nm, the chain gyration radius, rchain = 9 ± 2 nm, and bsph/bchain = 0.8 ± 0.1. The overall particle radius from SANS was thus, rsph + 2rchain = 36 ± 7 nm which was in accordance with the light scattering data.
Fig. 6 SANS curve of MFC/PMAA–PEO mixture (PMAA–PEO concentration 2 mg mL−1, w = 4). The red line is the fit of the data with the PG model. |
MFC/PMAA–PEO NP morphology was further investigated by cryogenic transmission electron microscopy for MFC-to-PMAA mass ratios, w = 1, 2 and 4. Cryo-TEM measurements revealed no substantial differences among NPs with different compositions. Fig. 7 shows the Cryo-TEM image for w = 1 together with the frequency histogram of NP sizes. The images for w = 2 and w = 4 are shown in ESI† (Fig. S11b and c). Assuming that the hydrated coronas did not have sufficient contrast to be imaged and only NP cores were visible, the results for w = 4 could be compared with the SANS analysis which provided the radius of the spherical core from the PG model, rsph = 18 ± 3 nm. The core radius from Cryo-TEM was twice lower, 9 ± 2 nm, which was the difference that could not explained only by the increased weight of larger NPs in the SANS data in contrast to number-averaged value from Cryo-TEM.
Fig. 7 (a) Cryo-TEM micrograph of MFC/PMAA–PEO aqueous solution (w = 4). (b) Frequency histogram of the radii of nanoparticles shown in Fig. 7a. |
In order to investigate thermodynamics of MFC/PMAA–PEO complex formation, we conducted an ITC experiment. The binding isotherm is shown in Fig. 8. The interpretation of the thermodynamic behaviour was complicated by the enthalpy change associated with the dilution of MFC (Fig. 8 inset, curve 2), which was not negligible as compared with the enthalpy of MFC mixing with PMAA–PEO.
The enthalpy curve showed that the process was exothermic up to w = 1, followed by the endothermic process for w > 1. A similar two-step behavior was reported for coacervation of oppositely charged polyelectrolytes where the exothermic process was associated with pairing of oppositely charged polyions while the following segregation of the complex into coacervate was endothermic.36 In the case of the MFC/PMAA–PEO system, the comparison with light scattering data showed that the exothermic-to-endothermic behaviour transition occurred after the formation of the NPs so that the endotherm corresponded to their rearrangement (swelling) and disruption.
Zeta potential measurements showed that aggregation of PMAA chains upon the interaction with MFC occurred at low MFC amounts, thus yielding particles with high charge density in the core and hence strongly negative zeta potential. This finding confirms the role of non-electrostatic interactions (such as the formation of guanidinium ion pairs, allowing for cooperative binding of MFC to PMAA chains) in the stability of interpolyelectrolyte complexes involving biguanide polycations.
Footnote |
† Electronic supplementary information (ESI) available: Molar mass distribution of MFC from size exclusion chromatography, 1H, 13C, 1H–13C HSBC and 1H–13C HSQC NMR spectra of metformin hydrochloride, 1H, 13C, 1H–1H COSY, 1H–13C HSBC and 1H–13C HSQC NMR spectra of MFC, discussion of NMR spectra and Cryo-TEM images of MFC and MFC/PMAA–PEO aqueous solutions. See DOI: https://doi.org/10.1039/d4sm00851k |
This journal is © The Royal Society of Chemistry 2024 |