Lability of Nanoparticulate Metal Complexes at a ... - ACS Publications

The nanoparticle (NP) size domain has been shown to exhibit unique reactivity features towards target ions and small ... (nano)particulate metal compl...
0 downloads 8 Views 2MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Article

Lability of Nanoparticulate Metal Complexes at a Macroscopic Metal Responsive (Bio)Interface: Expression and Asymptotic Scaling Laws Jerome F.L. Duval, Raewyn M. Town, and Herman P. Van Leeuwen J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b11982 • Publication Date (Web): 19 Feb 2018 Downloaded from http://pubs.acs.org on February 25, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Lability of Nanoparticulate Metal Complexes at a Macroscopic Metal Responsive (Bio)Interface: Expression and Asymptotic Scaling Laws ∗ Jérôme F.L. Duval ,1,2 Raewyn M. Town,3,4 Herman P. van Leeuwen4

1

CNRS, Laboratoire Interdisciplinaire des Environnements Continentaux (LIEC), UMR 7360, Vandoeuvre-lès-Nancy, F-54501 France. 2 Université de Lorraine, LIEC, UMR 7360, Vandoeuvre-lès-Nancy, F-54501 France. 3 Systemic Physiological and Ecotoxicological Research (SPHERE), Department of Biology, University of Antwerp, Groenenborgerlaan 171, 2020 Antwerp, Belgium. 4 Physical Chemistry and Soft Matter, Wageningen University & Research, Stippeneng 4, 6708 WE Wageningen, The Netherlands. ∗

Corresponding author. Email: [email protected] Laboratoire Interdisciplinaire des Environnements Continentaux- UMR 7360 Université de Lorraine, 15 avenue du Charmois, 54501 Vandoeuvre-les-Nancy, France.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract. The lability of metal complexes expresses the extent of the dissociative contribution of the complex species to the flux of metal ions toward a macroscopic metal-responsive (bio)interface, e.g. an electrodic sensor or an organism. While the case of molecular ligands is well established, it is only recently that a definition was elaborated for the lability of metal complexes with nanoparticles (NPs) in aqueous dispersions. The definition includes the thickness of the non-equilibrium reaction layer operational at the (bio)interface and the extent of geometrical exclusion of NPs therefrom. In this work, an explicit expression is derived for the lability of nanoparticulate metal complexes (M-NP) towards a macroscopic reactive (bio)interface. Interpretation accounts for the M-NP chemodynamic properties that depend on NP size, electrostatics, metal diffusion and dehydration rates, and density of metal binding sites for various NPs types, i.e. soft/core-shell and hard NPs having volume and surface sites distribution, respectively. Computational examples under practical conditions illustrate how these factors jointly determine the remarkable non-monotonous dependence of M-NP lability parameter on NP size. The analysis is supported by the formulation of asymptotic scaling laws clarifying how local M-NP dissociation dynamics affect the lability parameter for M-NP complexes at the scale of the macroscopic (bio)interface.

2 ACS Paragon Plus Environment

Page 2 of 35

Page 3 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1. Introduction. The nanoparticle (NP) size domain has been shown to exhibit unique reactivity features towards target ions and small molecules.1-3 The modified reactivity of various target species upon their association with NPs has broad significance for e.g. robust risk assessments of NPs,4,5 and design of NPs with tailored uptake/release kinetics for application in e.g. drug delivery.6-10 The interpretation calls for coupling the physicochemical features of NPs, e.g. size, charge density, chemical functionality, with their chemodynamic behaviour under the conditions encountered in environmental and biological media. The conceptual framework we develop herein quantitatively describes the reactivity of NP-target associates at a fundamental mechanistic level. The concepts are applicable to any target species: without loss of generality, we take the example of metal ion (M) complexes with negatively charged NPs. Recent work has elaborated a conceptual framework for evaluation of the chemodynamics of (nano)particulate metal complexes, i.e. the dynamics of interconversion between free and complexed metal forms.1,2,11-17 The work extends the original Eigen mechanism18 for formation of complexes between M and molecular ligands to the more intricate situation of charged complexing agents with nanoparticulate and colloidal dimensions. The confinement of metal binding sites within the volume of the permeable particle body or shell (3D soft/core-shell NPs) or at the particle surface (2D hard NPs) leads to metal complex chemodynamics that generally involve a balance between chemical reaction kinetic and diffusive metal transport contributions.1,2,11-17 Several reports demonstrated that these contributions are strongly affected by the electrostatic conditions prevailing in the vicinity, at the surface and/or in the volume of the NP body.2 Depending on medium salinity, NP charge and size, the particulate electric field may lead to a significant increase of the concentration of metal ions at the surface or in the volume of the NPs, and subsequently to a kinetics of M-NP complex formation that may be several orders of magnitude faster than that for complexes with molecular ligands.1 Given this knowledge on the chemodynamic properties of M-NP complexes, it was recently possible to elaborate a definition of their lability at a metal-responsive macroscopic interface such as a sensing electrode or the surface of an organism.13 The lability parameter is a measure of the extent to which the dissociation of the complexes contributes to the flux of metal ions toward the reactive sensing interface.13 Such formulation of M-NP lability is mandatory for e.g. a proper evaluation of the bioavailability of metal ions to organisms in the presence of complexants,19 for interpreting measurements obtained by dynamic metal speciation techniques such as various voltammetries20-22 and Diffusive Gradients in Thin Films (DGT),23-25 as well as in predicting release kinetics in the context of nanodrug formulations.7 At the labile limit, complexation of metal ions by nanoparticulate ligands is at equilibrium for any position from the sensing surface so that the overall flux of metal to the sensor results from the coupled diffusion-controlled transport of M-NP and free M.13 In the other limit of non-labile complexes, the complexation equilibrium 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

is not maintained and the relevant metal flux to the sensor is governed by the rate at which metal ions are released from the particle body.13 For intermediate situations, the lability parameter is determined by the ratio between the kinetic dissociation flux of complexes and the coupled diffusion flux of M-NP and M (Figure 1). A convenient way to formulate the dissociative kinetic flux of M-NP complexes at the sensor surface is to invoke the concept of a reaction layer, originally developed by Heyrovský and Brdička26-28 for molecular complexing ligands and recently detailed for nanoparticulate ligands.29 The thickness of this reaction layer, λ , basically derives from the relative mobilities of free and complexed metal species in the medium and their respective lifetimes as determined by the rate of re-association of M with NPs and the rate of M-NP dissociation.2,29 In turn, the spatial zone lying within a distance λ from the sensor defines a non-labile region where metal release from the complexes to the macroscopic interface is operational, i.e. λ denotes an effective reaction layer thickness at the macroscopic reactive interface.2,29 By its very nature, a reaction layer can only exist in the presence of complexing sites.13,30 Accordingly, for the case of NP complexants, the reaction layer at the sensing surface is an operational one that is related to the time-averaged presence of particle body volume/surface area and the corresponding time-averaged complexing site concentration. Unlike molecular ligands,25-28 the dissociation kinetic flux of M-NP complexes, or for that matter, the reaction layer thickness λ , generally depends on the degree of geometrical exclusion of the complexing NP’s body from the reaction layer, i.e. on the fraction of metal binding sites effectively involved in the dissociation of the complexes at the metal-sensing interface.2,13,29 Recent successful confrontation between the above lability framework and electrochemical data collected for Cd(II) in a dispersion of polymer nanoparticles evidenced the paramount importance of NP exclusion effects for a proper assessment of the lability of the nanoparticulate metal complexes considered.2,29 In contrast, the application of the conventional reaction layer approach - which neglects such effects25-28 was shown to overestimate the lability parameter of the complexes by more than 3 orders of magnitude.29 In the present work, an expression is developed for the lability parameter of M-NP complexes at a macroscopic metal sensor with explicit integration of the local M-NP chemodynamic features in conjunction with the macroscopic scale kinetic and transport-mediated exchange of metal species at the sensor interface. Such elaboration generalizes previous work13 and is further complemented by the derivation of useful tractable asymptotic scaling laws for analysis of the cumbersome dependence of lability of nanoparticulate complexes on their size. In turn, the results allow for a clear identification of the multi-factorial processes that control the magnitude of the M-NP lability parameter over the practical range of particle size as a function of e.g. the particle electrostatic features, the density of metal binding sites located within the volume or at the surface of the particle, the dehydration kinetics of the complexing metal ion (Eigen type association) and the intrinsic stability of the nanoparticulate metal complex.

4 ACS Paragon Plus Environment

Page 4 of 35

Page 5 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

2. Theory: setting the stage. In the following, we consider charged spherical nanoparticles of radius rp dispersed in a medium ∗ containing aqueous metal ions M (valence zM ) with bulk concentration cM and ions from a monovalent

background electrolyte in large excess over M. In line with usual practice, the electric potential profile at the NP/medium interface is assumed to be pre-established by the indifferent electrolyte ions.2 The binding of metal ions to NPs occurs via reactive sites S, that are present in excess over the metal ions in solution, and distributed either at the NP surface or within the particle body. The former case refers to 2D hard complexing NPs with surface local site density denoted as ΓS and the latter to soft/permeable and coreshell particulate systems with intraparticulate volume site density cs .17,29 We introduce the ratio z = a / rp with a being the NP core radius so that soft, core-shell and hard complexing NPs refer to z = 0 ,

0 < z < 1 and z = 1 , respectively. The volume fraction of particles in dispersion is denoted as ϕ . For many practical cases with ϕ > J diff (ℒ >> 1 ), i.e. the metal flux at the sensor surface is determined by coupled diffusion of M

and M-NP, whereas the inequality J kin DM / DML ,

(

∗ DM / ka' the metal flux at the macroscopic surface in the non-labile limit is given by J kin = kd cML

(

with ka' = ka cL∗ . The latter expression is generally written J kin = kd c∗ML µ where µ = DM / ka'

)

1/2

)

1/2

is the

thickness of the conventional reaction layer involving the kinetic constant ka' for M association with L, and kd identifies with kdis according to Eigen association scheme between M and L.18 In the case where the transport of the released M from inside to outside of the NP is rate-limiting the overall M-NP

8 ACS Paragon Plus Environment

Page 9 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

dissociation process, the system is labile, as argued in previous work.30 This expression of J kin can be ∗ λ used for the situation of NP complexants provided that cML is replaced by cMS , the concentration of

inner-sphere complexes MS averaged over the entire dispersion volume with account of the only fraction of the reactive shell volume located within the reaction layer of thickness λ (Figure 1, see §3). This thickness is the pendant of µ defined above for molecular L (see §3).2,13,29 It is important to realize that the condition Kc∗L >> DM / DML argued to derive J kin (and subsequent µ ) is generally not satisfied for NP complexants because NP diffusion coefficient Dp (the analogue of DML ) is lower than DM to an extent that depends on the NP radius rp . This implies that the dissociative term DM / ka' alone is not sufficient for defining the reaction layer thickness λ relevant for NP complexants at a reactive sensor/solution interface. In turn, an MS formation term should be included in the definition of λ to account for the association between free M and S. This additional (associative) term corresponds to

kd / Dp in eq. (2). In that equation, kd (the reciprocal life time of MS) is the composite kinetic constant

(

kd = 1 / kdis + 1 / kd,p

)

−1

that should be considered in order to cover the entire range of M-NP dissociation

dynamics, i.e. from slowly (e.g. Ni) to rapidly (e.g. Pb) dehydrating metal ions.12,17 Finally, the

dissociative term DM / ka' involved in the defining expression of µ should be replaced - for the NP case -

(

)

by DM / ka,pcsλ . This term involves the concentration csλ of binding sites effectively present within the reaction layer (see §3) and smeared out over the entire dispersion volume, and the diffusion-component

ka,p for M-NP association. It is stressed here that it is ka,p and not the composite constant ka that must be considered because the relevant timescale is the one required to get M inside the volume of the NP complexant (Figure 1). All these elements motivate the form of eq. (2) that defines λ . It is further emphasized that the original paper51 that elaborated an extended definition of the reaction layer thickness used a notation based on the involved rate constants, which has carried into subsequent work.2,17,29 Herein we adopt a process-based nomenclature and define the association reaction layer

(

thickness as kd / Dp

)

−1/2

(

(

and dissociation layer thickness as DM / ka,p csλ

))

1/2

. The latter enters the

original Heyrovsky formulation26 of the kinetic flux due to dissociation of MS. The former involves kd that defines the lifetime of the complex MS. The corresponding flux is an associative one because it

(

counts the total production of MS in the layer of thickness kd / Dp

)

9 ACS Paragon Plus Environment

−1/2

.

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Finally, the basis of eqs. (1)-(3) is supported by the excellent agreement between the lability parameter measured by electrochemistry for Cd(II)-latex particles (Table 3 in Ref. [29], and Figure 10 in Ref. [2]) and the treatment of the lability parameter outlined above for NP complexants at a macroscopic metal sensor surface. It is stressed that this treatment did not require any adjustable parameters given the known physicochemical properties of the latex NPs considered, both in terms of size and density of sites S derived from independent electrokinetic analysis.2,29 In particular, the comparison between experiments and theory highlighted the dramatic impact of the geometric exclusion of NP from the reaction layer (thickness λ ) on the lability parameter ℒ. Adopting the standard reaction layer formulation

(

µ = DM / ka'

)

1/2

with ignoring NP exclusion process, ℒ was found to exceed by more than 3 orders of

magnitude the experimental value. In addition, evaluation of ℒ for Cd-latex NP from eqs. (1)-(3) with ∗ arbitrarily ignoring the DM cM term in eq. (3) results in a lability parameter that is ca. 4 times larger than

the experimental value. These elements further comfort our defining expressions here for J kin , J diff and

λ . Within the context of lability evaluation for molecular complexes ML, some authors52 define J diff as ∗ the maximum diffusive supply of the complexes without any contribution from free M (i.e. DM cM = 0 in

eq. (3)). As previously argued, such definition would lead to a poorer comparison to experimental data on Cd-latex NP system, and to ignoring the necessarily coupled diffusions of free M and complexed metal species to the sensor surface as evidenced by van Leeuwen.50 Last, other operational definitions of the reaction layer thickness are proposed in the literature for the case of molecular ligands.25 At this stage, it is difficult to appreciate whether or not the proposed definitions can be easily transposed to the intricate situation of NP complexants with proper integration of their defining physicochemical features and of e.g. the paramount contribution of NP exclusion from the reaction layer.2,29 The solution proposed in this report and elsewhere2,13,29 achieves such a transition between molecular to NP complexant situations.

In the sections below (§3-4), starting from the eqs. (1)-(3) an expression for the lability parameter ℒ = J kin / J diff of M-NP complexes at a reactive (bio)interface is explicitly derived as a function of (i) the

relevant physicochemical features of the nanoparticulate ligands, which includes their size, their electrostatic characteristics and the density of metal binding sites S they carry, (ii) the metal diffusion and dehydration kinetic properties, and (iii) the stability of the intraparticulate or surface metal complexes MS. For that purpose, differentiation is made between the case of 3D soft/core-shell complexing nanoparticles (§3) and that of 2D hard nanoparticulate ligands (§4). It is emphasized that the developments below hold for cases in line with the applicability of the generalized Eigen

10 ACS Paragon Plus Environment

Page 10 of 35

Page 11 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

scheme1,2,12,14,16,17 (Figure 1) in the high NP charge density regime ( κ l c 1 . Previous work17,29

elaborated diagnosis for chemodynamic behaviour for M-NP association and dissociation from criteria involving the ratios l a / rp and µd / rp , respectively. Unlike the choice originally made in Ref. [17], the term

(1 + Kint cs ) Kint cs

could have been left outside the defining expression of µd , then resulting in

µd / l a = 1 . The dimensionless scalars fel,a involved in Φ el,a and f el,d = f el,a / f B appearing in µd correct the convergent (subscript ‘a’) or divergent (subscript ‘d’) diffusion of M to or from the NPs for particle electrostatic effects, respectively. Quantitatively, these scalars are obtained from solution of the non-linear Poisson-Boltzmann equation as detailed elsewhere,17,53 and their deviation from unity indicates the extent to which diffusion of M is accelerated or impeded by the particle electric double layer field. For practical situations where NPs are negatively charged, in contrast to the behaviour of f B , the factor fel,a does not exceed ca. one order of magnitude,12,17 so that the inequality Φ el,a