|Home | About | Journals | Submit | Contact Us | Français|
Microelectromechanical systems (MEMS), a technology that resulted from significant innovation in semiconductor fabrication, have recently been applied to the development of silicon nanopore membranes (SNM). In contrast to membranes fabricated from polymeric materials, SNM exhibit slit-shaped pores, monodisperse pore size, constant surface porosity, zero pore overlap, and sub-micron thickness. This development in membrane fabrication is applied herein for the validation of the XDLVO (extended Derjaguin, Landau, Verwey, and Overbeek) theory of membrane transport within the context of hemofiltration. In this work, the XDLVO model has been derived for the unique slit pore structure of SNM. Beta-2-microglobulin (B2M), a clinically relevant “middle molecular weight” solute in kidney disease, is highlighted in this study as the solute of interest. In order to determine interaction parameters within the XDLVO model for B2M and SNM, goniometric measurements were conducted, yielding a Hamaker constant of 4.61× 10−21 J and an acid-base Gibbs free energy at contact of 41 mJ/m2. The XDLVO model was combined with existing models for membrane sieving, with predictions of the refined model in good agreement with experimental data. Furthermore, the results show a significant difference between the XDLVO model and the simpler steric predictions typically applied in membrane transport. The refined model can be used as a tool to tailor membrane chemistry and maximize sieving or rejection of different biomolecules.
Theoretical models for solvent and solute transport have been at the core of membrane science since the development of highly selective filtration membranes in the 1960s. These models, which have been summarized extensively in the literature , have been used to describe solute transport in a variety of membrane types and processes. Most models investigating solute sieving through membranes have focused on size-based (steric) and electrostatic (EL) interactions, and even in cases where models have been exhaustively fit to experimental data, numerous assumptions are required to describe data obtained with polymeric membranes due to the broad pore size distribution and irregular pore shape in most polymeric membranes . Furthermore, the few studies investigating other interactions, such as Lifshitz van der Waals (LW) and acid-base (AB) interactions, have been only cursorily validated . There is, therefore, a desire to fabricate a membrane with highly tunable, uniform pore geometry for investigating transport model variations more rigorously than is possible with polymeric membranes.
Most polymeric membranes are made by a phase-inversion process, resulting in non-uniform lateral and axial pore dimensions, with significant uncertainties (and spatial variability) in the active layer thickness, surface porosity, and pore tortuosity. This makes it impossible to quantitatively validate models describing the nanoscale interactions between solute and membrane during sieving. In order to limit the variability in these membrane parameters, investigations have been performed with “track etched” membranes introduced in the 1970s, with the results showing rough agreement with membrane transport models . However, track etched membranes exhibit considerable limitations including variable pore number, inflexible pore design (with pore size greater than 15 nm), large thickness, and pore overlap .
Recently, silicon nanopore membranes (SNM) have been developed for use in numerous biomedical applications, including hemofiltration, the implantable renal assist device, the bioartificial pancreas, and extracorporeal membrane oxygenation . SNM possess nearly monodisperse pore size and tunable pore geometry. Unlike other silicon membranes with sub-micron pores, which are made using a long and laborious focused ion-bean (FIB) drilling fabrication process , SNM can be fabricated relatively rapidly using microelectromechanical systems (MEMS) fabrication techniques. Importantly, SNM can be fabricated with a very thin membrane thickness (~200–500 nm) and zero pore overlap. Figure 1 depicts a plan-view scanning electron microscopy (SEM) image of the membrane pores. Details on the manufacture of SNM membranes are provided elsewhere .
The purpose of this work is to apply the precisely known membrane porosity and pore dimensions of SNM to rigorously investigate, within the context of hemofiltration, the transport of macromolecules through membrane pores. This includes the full spectrum of interactions between membrane and solute, including LW, AB, and EL contributions (collectively, the extended Derjaguin, Landau, Verwey, and Overbeek, or XDLVO, theory). While the SNM technology itself has been well-characterized previously, the application of SNM to membrane pore model validation represents a new and significant opportunity. As currently fabricated SNM have a slit pore design, we first derive the equations for the interaction energy for a spherical solute in a slit pore and the osmotic reflection coefficient for slit pores in the following section. Since we are evaluating the validity of the XDLVO transport model from the perspective of kidney disease treatment, we select a clinically relevant middle molecular weight protein, Beta-2-microglobulin (B2M), as the solute of interest .
Solvent and mass transport through membranes can be determined by solving the relevant mass transfer equations. The rate of water transport, Jw, across the membrane can be written as
where Lp is the hydraulic permeability of the membrane, ΔP is the hydraulic pressure difference across the membrane, Δπ is the osmotic pressure difference, and σ is the osmotic reflection coefficient. Mass transport through membrane pores is due to both convective and diffusive contributions. By solving the relevant differential mass transport equation across the membrane, the following is obtained for the actual sieving coefficient, Sa,
where cf is the filtrate solute concentration, cb, m is the solute concentration in the feed solution at the membrane surface, S∞ is the asymptotic sieving coefficient corresponding to convective transport alone, and Pe is the Peclet number, i.e. the ratio of convective to diffusive transport in the membrane . The Peclet number itself can be calculated as
where Jw is the water flux, δm is the membrane thickness, ε is membrane porosity, ϕ is the solute-membrane partition coefficient, Kd is the diffusive hindrance factor, and D∞ is the bulk diffusion coefficient of the solute . The asymptotic sieving coefficient, S∞, is the sieving coefficient at very high Peclet number and is equal to the product of ϕ and the convective hindrance factor, Kc, i.e., S∞.= ϕKc.
Since it is difficult to experimentally determine the concentration of solute at the membrane surface, we define an observed sieving coefficient as
where cb,b is the bulk concentration of solute in the feed solution. Through applying film theory and solving the differential equations for mass transfer across the feed-side boundary layer, we obtain
where kb is the feed side external mass transfer coefficient, defined as kb =ShD∞/dh where Sh is the Sherwood number and dh is the hydraulic diameter of the flow channel . For a long, thin rectangular channel, the Sherwood number can be approximated by ,
where Re = ρsubdh/μ is the Reynolds number, i.e. the ratio of momentum forces to viscous forces, and Sc=μ/ρsD∞ is the Schmidt number, i.e. the ratio of momentum diffusivity to mass diffusivity. In both cases, ρs is the solvent density and μ is the solvent viscosity.
The osmotic reflection coefficient for slit pores can be derived in an analogous method to cylindrical pores . The Gibbs-Duhem equation describes thermodynamic equilibrium within the pore as
which states that the lateral (x direction) pressure gradient must be balanced by a solute chemical potential gradient at a given axial (z direction) position. Here, E represents the energy at a given position. We further make the assumption that the Boltzmann equation is valid such that c(x,z) = co (z)exp[−(E(x)−E(0))/kT], where co is the concentration at the pore centerline, k is the Boltzmann constant, and T is temperature. If we assume that the x and z components of the solute potential can be separated, incorporating this expression into Eq 7 yields
where R is the universal gas constant. The Navier-Stokes equations for fluid transport can be used to solve for the velocity profile within a single membrane pore
where μ is the solvent viscosity and vz is the velocity in the z direction. Note that Eq 9 neglects the effects of counterelectroosmosis arising from the electrical stresses in the fluid. Combining Eq 8 and Eq 9, we obtain
Application of the appropriate boundary conditions yields an expression for the velocity at a given position relative to the pore centerline:
where h is the pore half width and πo = coRT is the osmotic pressure. In order to determine the solvent flux through the pore, the velocity is averaged over the entire pore width, resulting in
Since Eq 12 is written in terms of intrapore concentrations and pressures, these values must be transformed into values just outside the pore entrance through relations as described in Anderson et al 
If we non-dimensionalize the pore position and solute radius (rs) by the pore half width we get ρ = x/h and λ=rs/h, respectively. Finally, substituting Eq 13 into Eq 12 and rearranging we obtain an expression for the water flux through each slit pore
The impact of the osmotic reflection coefficient may be particularly important in systems where the solute concentration is relatively high. In systems where solute concentrations are low, the osmotic reflection coefficient has no effect on the filtrate flux because the contribution from solute osmotic pressure is negligible. This simplifies the model considerably and may have relevance in a number of filtration modalities.
The partition coefficient for a slit pore is the ratio of the bulk membrane surface concentration just outside the pore to the average concentration just inside the pore co ,
where g(ρ) is the potential function at some non-dimensional position ρ is the Born in the pore and do is the Born repulsion limit of approximately 0.187 nm . Note that the upper limit of the integration in Eq 15 is less than unity due to steric hindrance and Born repulsion effects. In other words, there is a restriction that the solute center cannot sample all positions within the pore since the solute itself cannot pass through the pore wall. For the case where LW, AB, and EL interactions are neglected, the partition coefficient depends only on steric considerations such that ϕ = 1−λ−(do/h). The potential function can be represented by a Boltzmann distribution through g(ρ)=exp(−E(ρ)/kT), where E(ρ) is the energy at a given non-dimensional pore position .
The diffusive and convective hindrance factors can be calculated through existing correlations for the transport of solutes through slit pores . These correlations exist for both centerline positions and cross-sectionally averaged off-axis positions. Since repulsive interactions between solute and membrane are likely to position solutes closer to the pore centerline, we apply the centerline approximation in this study. The centerline hindrance factors are calculated through 
which are valid up to approximately λ = 0.5.
In order to determine the impact of LW, AB, and EL interactions on transport, the interaction energy between the solute and membrane surface (or pore wall) must be calculated. This energy has typically been determined through application of the Derjaguin approximation. However, as Bhattarcharjee at al  note, this approximation has been carefully validated only for sphere-sphere interactions. Therefore, an alternative approach, the surface element integration (SEI) method, is applied for the sphere-plane scenario discussed in this paper .
This method has been applied previously for circular pores and solutes approaching a plane from infinite distance, but not for a slit pore membrane geometry . Therefore, we develop a mathematical framework for the slit pore geometry shown in Figure 2. Applying an approach similar to Bhattacharjee , we use cylindrical coordinates to describe the solute as a group of cylindrical disks, each with infinitesimal thickness, dr. The distance of a cylindrical disk from the solute origin is defined
where r is the radius of a differential area element (dA = 2πrdr) at some position d from the solute origin, and rs is the solute radius. Using Eq 18, we can determine the distance of a differential area from the pore wall (L) by accounting for the pore width and the position of the solute laterally within the pore, such that
Eq 19 provides the distance of a differential area with the pore wall closest to the half sphere ABC. For the half sphere closest to the opposite pore wall, the distance is
Integrating over the entire sphere at a given position results in the following expression
which is the energy for the solute at a given position within the pore. However, each force is a function of the distance between the pore surface and solute surface element; thus, the interaction energy will also vary with distance as follows
for LW, AB, and EL components, respectively. Here, L is the distance between two differential area elements A, is the Hamaker constant for Lifshitz van der Waals interactions, is the AB interaction at contact, D is the AB decay length, εo is the vacuum permittivity, εr is the relative permittivity, ψ1 and ψ2 are the surface potentials of solute and membrane, and κ is the inverse Debye length. The Hamaker constant is itself a function of the LW free energy of interaction at contact such that . The AB (hydration repulsion) force can exhibit oscillatory or monotonic decrease with distance depending on surface composition and roughness . For the relatively rough polymer coated SNM surface used in our membranes, we assume monotonic decrease of the acid-base repulsion as shown in Eq 23. The value of the decay length can change due to the amount of surface bound water. Surfaces with either positive or negative charge can attract counterions to the interface, creating a hydration shell of water molecules, resulting in an increase in the decay length.
The contribution of LW forces to the total interaction energy can be represented by the Hamaker constant, by the surface potential for electrostatics, and by the acid-base free energy of interaction for the acid-base force. While values for the surface potential of protein and membrane can be easily determined using literature data for surface charge and zeta potential, LW and AB forces are more difficult to determine experimentally. The van Oss-Good-Chaudhury (VCG) method can be applied in order to isolate the Hamaker constant and the Gibbs free energy for acid-base interaction at contact using goniometric measurements. Briefly, this method assumes that the total surface tension (γT) of a solid or liquid is comprised of individual van der Waals (γLW), acid (γ+), and base (γ−) components. As the individual components are not known explicitly for the great majority of surfaces, the VCG method provides a framework for determining them experimentally using three liquids of known surface tension. Then, using thermodynamic combining relations, the surface tension values of the unknown surface can be calculated by solving
where Eq 25 applies for the membrane and Eq 26 applies for the solute of interest . Here, θ is the contact angle between the liquid and membrane or solute, and the subscripts l, m, and s represent the liquid, membrane, and solute, respectively. Applying Eq 25 and Eq 26 using three probe liquids (each with different γl values) results in a system of equations which, when solved, yield the individual LW and AB surface tension parameters. Finally, these surface tension parameters can be used to determine the interaction energy at contact between the solute and membrane as follows 
B2M is of particular interest in this study due to its molecular dimensions (MW 11800), which allow for partial sieving at SNM pore sizes (rs = 1.59 nm) . Partial sieving will allow us to elucidate the effects of XLDVO interactions between solute and pore. Furthermore, since B2M has been heavily researched as a kidney function middle molecular weight marker, there is published data available on its electrostatic properties .
The fabrication of SNM is detailed extensively in other published work . Briefly, microelectromechanical systems (MEMS) fabrication techniques are applied to grow a thin silicon oxide sacrificial layer on a silicon substrate. Upon etching of the oxide, nanometer scale pores are formed in the membrane. Virgin SNM have a roughly 11 nm pore size, and to facilitate membrane handling, membranes are fabricated with 0.95% surface porosity and a thickness of 400 nm. The porosity is precisely known since the photolithographic mask results in exactly 4.343 × 105 pores/mm2.
In order to prevent protein adhesion and pore blocking of the membrane during filtration, membranes are further functionalized through silane chemistry . SNM were immersed in 2 % w/v polyethylene glycol (PEG) silane (Gelest) dissolved in 75 mL of toluene for 2 hr at 70 °C. Following PEG treatment, the membrane pore size decreased from roughly 11 nm to 6 nm and water constant angles stabilized at 35 – 40 degrees, suggesting significant surface modification. The evaluation of the membrane pore size is discussed in more detail subsequently. Pure water permeability experiments on PEG-coated membranes showed linear increases in flux with applied pressure, suggesting a consistent pore size and uniform coverage of PEG within the membrane pores.
Filtration experiments were conducted using a custom experimental setup. A syringe pump (KD Scientific Inc, #KDS200) was used to pressurize a reservoir containing the feed solution. A peristaltic pump (Cole Palmer Masterflex, #7551-00) provided fluid flow (2 ml/min) through a pressure transducer (Omega Engineering Inc, #PX429) and specially designed aluminum flow cell. Experiments were automated using a connected computer, with hydraulic pressures taken over a range from 3 to 55 kPa (0.5 to 8 psi), which includes the typical physiological arterial pressure (~13.8 kPa). Aluminum tubing (McMaster-Carr, #5177K61) was used to attach the different system components and to prevent significant protein adsorption. The experimental system was assembled and sterilized by autoclave before each set of experiments. B2M (Lee Biosolutions Inc, #126-11) was dissolved in divalent cation-free phosphate buffered solution (UCSF Cell Culture Facility), and B2M immunoassays (performed on the SPA Plus, The Binding Site) were conducted at the UCSF Clinical Chemistry Laboratory in order to determine the solute concentrations in feed and filtrate.
Membrane pore sizes were determined before filtration experiments through hydraulic permeability measurements using a feed solution of pure water. SNM were placed in the testing flow cell and subjected to feed hydraulic pressures of 1, 2, and 4 psi (6.9, 13.8, and 27.6 kPa). The filtrate volume collected over a time interval was weighed in order to determine the water flux across the membrane at a given pressure. Using this water flux (Q) data, the membrane surface porosity, membrane area (Am) and known membrane thickness, the SNM slit half-width was then back-calculated using the follow equation:
which is derived from the Navier-Stokes equation for flow through a slit pore. As stated earlier, ε represents the membrane porosity, which is itself a function of the known pore number, pore length, and slit width h. The porosity can be represented through ε = 2NhL/Am, where L is pore length and N is the number of pores. The “effective” pore size is then determined through Eq 29 by incorporating the line of best fit for the water flux versus pressure data. For a membrane of a constant, known geometry (in this case, a slit pore), the pore size should not change across all tested pressures.
Using Eq 29, we determined the post-PEGylation membrane pore sizes as 6.1 ± 0.2 nm and 6.3 ± 0.2 nm for the two membranes used in this study. The reduced pore size necessarily results in a decrease in the membrane surface porosity. For pore sizes in this range, there is complete rejection of albumin as confirmed by previous studies . This method assumes that the PEG chains grafted on the membrane pore surface form a uniform, dense structure and the large majority of the fluid flows in the open pore region. With this assumption, polymer chain compression is likely to be limited due to the steric constraints arising from adjacent chains. Thus, the PEG polymer chains form an effective “hard wall” when modeled with Eq 29.
Contact angle measurements were conducted on both PEG-coated SNM and B2M samples using a commercially available goniometer (Biolin Scientific, Theta Lite). To prepare the B2M for goniometric analysis, a sample of B2M (10 mg) in 30 mL of solution was filtered at 140 kPa (20 psi) onto 5 kDa ultrafiltration membranes (Sterlitech Corp, Trisep UF5) placed in a stirred ultrafiltration cell (Amicon, 8050). Before filtration, each membrane was compacted at the operating pressure in order to prevent passage of solute through the membrane. Once the B2M solution was filtered onto the Trisep membrane, excess surface moisture was allowed to evaporate (approximately 30 minutes). Contact angles for water, diiodomethane (Sigma Aldrich, #158429), and ethylene glycol (Sigma Aldrich, #85978) were then measured on the hydrated protein layer. Further details about the procedure for performing these measurements is available in the literature  . Electrostatic parameters for the B2M  and PEG-coated SNM  are sourced from the literature.
As discussed above, B2M and PEG-SNM surface tension parameters were determined through contact angle analysis. Table 1 summarizes the surface tension parameters, interaction energy constants, and feed solution conditions. Comparing the B2M surface tension parameters to those for hydrated human serum albumin (HSA, from literature) , B2M exhibits a greater overall contribution from LW forces and a slightly lower polar AB component. B2M possesses similar LW and polar AB components to Immunoglobulin G (IgG, from literature) , but exhibits less apolar AB activity. The similarity to IgG is not surprising, as the amino acid sequence homology of B2M is related to the Ig constant domain . PEG-SNM surface tension values reported here are similar to values for PEG reported in the literature . Since the PEG-SNM is also highly polar, this results in large Gibbs free energy at contact with B2M. Similarly, as the PEG-SNM also possesses a large LW component, the Hamaker constant is high for the PEG-SNM/B2M pair. Note that the surface potential for B2M is significantly lower than that for BSA (−20 mV, from literature) . Therefore, there is less EL repulsion between B2M and PEG-SNM than there would be for a hypothetical BSA/PEG-SNM pair . The high ionic strength of the feed solution also limits the contribution of EL forces due to double layer shielding.
It is important to first consider the effect that each individual interaction energy component has on the total interaction between solute and membrane pore. The interaction energy (normalized by kT) for different positions relative to the pore centerline is shown in Figure 4 for LW, AB, and EL interactions. Here, the position relative to the centerline is a function of the pore half width, solute radius, and Born repulsion limit, such that xr = 1−ρ−do/h. Three pore half-widths are modeled, and in each case B2M is taken as the solute of interest. It is clear from the results that the biggest contribution to the interaction energy is the AB interaction. This is expected given the large value for AB interaction at contact between B2M and PEG-SNM as determined by contact angle measurements. As the solute moves farther from the closest point of contact, it rapidly approaches the opposing pore wall. Note that we use a decay length D = 0.1 nm for the acid-base interaction. Decay length values from literature are typically in the range of 0.1 to 0.6 nm, and we selected from the lower end of this range due to the non-ionic nature of the PEG-based coating on the SNM . Hydrophilic coatings with ionic constituents, such as zwitterions, can demonstrate much larger decay lengths due to the larger hydration shell around these moieties.
There is a considerably lower impact from LW and EL forces for the interaction of B2M and the PEG-SNM near the pore wall. However, since the AB interaction decays relatively quickly with D = 0.1 nm, the contribution of LW forces at the pore centerline remains much closer to the value at contact than for AB forces. For EL interactions, there is an almost negligible contribution to the interaction energy for all three pore sizes. This is a consequence of the small net charge on B2M at physiological pH in buffered solution. Proteins with a much greater net charge, such as albumin, will exhibit greater electrostatic repulsion, especially in cases where the pore wall itself is adsorbed with albumin. In non-physiological medium, such as very dilute salt solutions, the increased Debye length will result in a greater EL contribution.
It should be noted that the EL interaction does not appear to be monotonic close to the pore wall. This is the result of the behavior of the Debye-Huckel equation (Eq 24), which as presented is solved for the case of constant surface potential (CP). It has been shown previously that the behavior of the electrostatic force is markedly different when solved at CP compared to that at constant surface charge (CC) . This discrepancy is particularly pronounced when the difference between the potentials of interacting surfaces is relatively large, with the more strongly charged surface inducing an opposite charge on the other surface at small separations. Note that the “actual” interaction between the solute and pore is likely somewhere between the CP and CC solutions. Due the limited impact of EL interactions in our modeled scenario, we apply the CP case and assume that any deviations are minor relative to the AB and LW contributions.
As discussed earlier, the partition coefficient is a function of the integrated interaction energy over the pore width. Figure 5 shows the partition coefficient for a range of λ values (at constant rs= 1.59 nm for B2M) for different values of the AB decay constant. For the steric only condition, it is clear that ϕ decreases as the solute radius approaches the pore half width. As the AB decay constant increases, ϕ is shifted lower for a given λ value. At these larger decay constant values, the partition coefficient is essentially zero for even moderate values of λ. In practical terms, this means that shifting the AB decay constant through attachment of different modifiers to the silicon surface could be beneficial or detrimental to filtration performance. For example, in order to limit solute transmission, zwitterionic polymers with large AB decay constants could be covalently attached to the membrane. On the other hand, lower decay constant polymers such as PEG will only mildly prohibit sieving compared to that predicted based only on the steric exclusion. Note that, although certain surface modifiers such as zwitterionic polymers may have increased AB decay constants, there may be resultant increases in LW and EL contributions as well.
Like the partition coefficient, the osmotic reflection coefficient is a function of the interaction energy between the membrane and solute. As the solute radius approaches the pore half width, σ increases monotonically (Figure 6). It is clear from the model output that increasing the AB decay constant results in a much faster rise in σ at low values of λ. In this study, the osmotic reflection coefficient is unlikely to have a significant impact on sieving since the feed B2M concentration is relatively low. However, filtration of blood or analogous fluids with particularly high concentrations of certain proteins might exhibit decreased solvent flux compared to the pure water case. While the partition coefficient approaches zero at large D for relatively small values of λ, the osmotic reflection coefficient only approaches unity at large and λ. This behavior results from the differences in the way each parameter is calculated, with ϕ and σ requiring single and triple integration of the interaction energy, respectively.
Figure 7 depicts the modeled change in the observed and actual sieving coefficients for a variety of pore sizes. These simulations incorporate the hindrance factors described earlier (Eq 16 and Eq 17), which account only for steric effects. Note that the observed and actual sieving coefficients both initially decrease with increasing pressure, but exhibit different behavior at higher pressure. As convection begins to dominate diffusion at higher pressures, the actual sieving coefficient approaches its asymptotic value (at infinite convection). However, as convection dominates, concentration polarization causes a higher concentration of solute in solution at the feed-side membrane interface. Therefore, since the bulk solute concentration is relatively constant and lower than the increasingly higher interface concentration, the observed (or “apparent”) sieving coefficient increases to unity at higher pressures. Experimental validation for the minimum in the observed sieving coefficient at intermediate flux is available in the literature .
Experiments were conducted for B2M sieving using two membranes with different pore sizes after PEGylation (6.1 nm and 6.3 nm). With the exception of a data point at 1 psi for the membrane with 6.1 nm pore size, the experimental data match the XDLVO model well (Figure 8). The experimental data are also noticeably lower than predictions of the steric only model, which neglects the interaction energy between the pore and solute. Given the small net charge of B2M, the reduction in transmission compared to that predicted by the steric model is predominantly due to LW and AB interactions and their impact on the partition coefficient. Therefore, DLVO models, which inherently neglect these interactions, would have significantly overestimated the experimental sieving data. In our experiments the pure water and protein solution water permeability were approximately the same over all tested pressures, indicating that membrane fouling is unlikely to be responsible for the difference between the steric only and the XDLVO model predictions. This suggests that experimental and model differences are not the result of membrane changes resulting from pore occlusion by proteins. Note that the model predictions are calculated for D = 0.11 nm, which here acts a fitting parameter to model data for both membrane pore sizes tested. We believe this to be a reasonable value for the AB decay constant and consistent with the range of values predicted in the literature (see earlier discussion, Section 4.2.1). A shift in the value for D will indeed alter the model fit and agreement with experimental data.
The purpose of this work was to evaluate the accuracy of the XDLVO model for membrane transport using slit pore SNM with uniform pore size and surface porosity. Equations for the partition coefficient and osmotic reflection coefficient for slit pore membranes were developed accounting for LW, AB, and EL interactions using surface element integration. The results of the work demonstrate that the models can be used to accurately predict the sieving of B2M, a particularly relevant middle molecular weight molecule in applications of the artificial kidney.
A key conclusion of this work is that AB repulsion between solute and membrane can cause a marked decrease in the sieving coefficient of B2M. While the reduction in sieving coefficient due to XLDVO interactions is fairly small for the AB decay constant applied here (approximately 10–20%), other surface modifiers such as zwitterionic polymers (with larger decay constant) could cause a much more pronounced reduction in solute transmission.
The implications of this work for the design of biomedical devices and other processes are significant and should be investigated in future studies. If solutes are intended to be excluded from pores, large decay constant hydrophilic polymers could be grafted to the pore walls, allowing one to use membranes with larger pore size to maximize water permeability while still retaining the solutes of interest. On the other hand, if the aim is to increase solute transmission through the pores, polymers with small decay constants would facilitate penetration of solutes into the pores. Ultimately, this paper validates a model that can serve as a valuable tool in demonstrating how SNM can be optimized to achieve the desired removal of clinically relevant middle molecular weight solutes.
We would like to thank the National Institutes of Health (NIH) for funding this research work through the following grants: R01EB014315 and U01EB021214, both of which were provided through the National Institute of Biomedical Imaging and Bioengineering (NIBIB).
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.