Selective Filtering of Particles by the Extracellular Matrix: An Electrostatic Bandpass

The transport of microscopic particles such as growth factors, proteins, or drugs through the extracellular matrix (ECM) is based on diffusion, a ubiquitous mechanism in nature. The ECM shapes the local distribution of the transported macromolecules and at the same time constitutes an important barrier toward infectious agents. To fulfill these competing tasks, the hydrogels have to employ highly selective filtering mechanisms. Yet, the underlying microscopic principles are still an enigma in cell biology and drug delivery. Here, we show that the extracellular matrix presents an effective electrostatic bandpass, suppressing the diffusive motion of both positively and negatively charged objects. This mechanism allows uncharged particles to easily diffuse through the matrix, while charged particles are effectively trapped. However, by tuning the strength of this physical interaction of the particles with the biopolymer matrix, the microscopic mobility of formerly trapped particles can be rescued on demand. Moreover, we identify heparan sulfate chains to be one important key factor for the barrier function of the extracellular matrix. We propose that localized charge patches in the ECM are responsible for its highly unspecific but strongly selective filtering effect. Such localized interactions could also account for the observed tunability and selectivity of many other important permeability barriers that are established by biopolymer-based hydrogels, e.g., the mucus layer of endothelial cells or the hydrogel in the nuclear core complex.     

Particle Transport through Hydrogels Is Charge Asymmetric

Transport processes within biological polymer networks, including mucus and the extracellular matrix, play an important role in the human body, where they serve as a filter for the exchange of molecules and nanoparticles. Such polymer networks are complex and heterogeneous hydrogel environments that regulate diffusive processes through finely tuned particle-network interactions. In this work, we present experimental and theoretical studies to examine the role of electrostatics on the basic mechanisms governing the diffusion of charged probe molecules inside model polymer networks. Translational diffusion coefficients are determined by fluorescence correlation spectroscopy measurements for probe molecules in uncharged as well as cationic and anionic polymer solutions. We show that particle transport in the charged hydrogels is highly asymmetric, with diffusion slowed down much more by electrostatic attraction than by repulsion, and that the filtering capability of the gel is sensitive to the solution ionic strength. Brownian dynamics simulations of a simple model are used to examine key parameters, including interaction strength and interaction range within the model networks. Simulations, which are in quantitative agreement with our experiments, reveal the charge asymmetry to be due to the sticking of particles at the vertices of the oppositely charged polymer networks.     

Solute size effects on the diffusion in biofilms of Streptococcus mutans

The diffusion of poly(ethylene glycol) (PEG) (MW varying between 200 and 10,000), and of three different types of micelles was examined in Streptococcus mutans biofilms using infrared spectroscopy. PEGs were used because they show limited interactions with biological materials and their weight can be selected in order to cover a wide range of size. The study showed that a considerable fraction at the base of the biofilm was not accessible to the diffusing solute molecules and this inaccessible fraction was very dependent on the size of the diffusing molecules. In parallel, it was found that the diffusion coefficients of these solutes in the biofilms were less than those in water and this reduction was less pronounced for large macromolecules, an effect proposed to be related to their limited penetration. Triton X-100, a neutral detergent, forms micelles that behave like PEG, suggesting that the behaviour observed for neutral macromolecules can be extrapolated to neutral macroassemblies. However, the diffusion, as well as the penetration of sodium dodecylsulphate micelles (a negatively charged surfactant) and cetylpyridinium chloride micelles (positively charged), in the biofilms appeared to be significantly influenced by electrostatic interactions with biofilm components. The present findings provide useful insights associated with the molecular parameters required to efficiently penetrate bacterial biofilms. The study suggests a rationale for the limited bactericidal power of some antibiotics (the large ones). The restricted accessibility of macromolecules and macroassemblies to biofilms must be examined carefully in order to offer guidelines in the development of novel antibacterial treatments.     

The kinetics of chemically induced nonequilibrium swelling of articular cartilage and corneal stroma

An electromechanical model for charged, hydrated tissues is developed to predict the kinetics of changes in swelling and isometric compressive stress induced by changes in bath salt concentration. The model focuses on ionic transport as the rate limiting step in chemically modulating electrical interactions between the charged macromolecules of the extracellular matrix. The swelling response to such changes in local interaction forces is determined by the relative rates of chemical diffusion and fluid redistribution in the tissue sample. We have tested the model by comparing the experimentally observed salt-induced stress relaxation response in bovine articular cartilage and corneal stroma to the response predicted by the model using constitutive relations for the concentration dependent material properties of the tissues reported in a related study. The qualitatively good agreement between our experimental measurements and the predictions of the model supports the physical basis of the model and demonstrates the model's ability to discriminate between the two soft connective tissues that were examined.     

Measurements of interbilayer forces and protein adsorption on uncharged lipid bilayers displaying poly(ethylene glycol) chains

Poly(ethylene glycol) (PEG)-stabilized liposomes were recently shown to exhibit differences in cell uptake that were linked to the liposome charge. To determine the differences and similarities between charged and uncharged PEG-decorated liposomes, we directly measured the forces between two supported, neutral bilayers with terminally grafted PEG chains. The measurements were performed with the surface force apparatus. The force profiles were similar to those measured with negatively charged PEG conjugates of 1, 2-distearoyl-sn-glycero-3-phosphatidyl ethanolamine (DSPE), except that they lacked the longer ranged electrostatic repulsion observed with the charged compound. Theories for simple polymers describe the forces between end-grafted polymer chains on neutral bilayers. The force measurements were complemented by surface plasmon resonance studies of protein adsorption onto these layers. The lack of electrostatic forces reduced the adsorption of positively charged proteins and enhanced the adsorption of negatively charged ones. The absence of charge also allowed us to determine how membrane charge and the polymer grafting density independently affect protein adsorption on the coated membranes. Such studies suggest the physical basis of the different interactions of charged and uncharged liposomes with proteins and cells.     

Nanoparticle diffusion measures bulk clot permeability

A clot's function is to achieve hemostasis by resisting fluid flow. Permeability is the measurement of a clot's hemostatic potential. It is sensitive to a wide range of biochemical parameters and pathologies. In this work, we consider the hydrodynamic phenomenon that reduces the mobility of fluid near the fiber surfaces. This no-slip boundary condition both defines the gel's permeability and suppresses nanoparticle diffusion in gel interstices. Here we report that, unlike previous work where steric effects also hindered diffusion, our system—nanoparticles in fibrin gel—was subject exclusively to hydrodynamic diffusion suppression. This result enabled an automated, high-throughput permeability assay that used small clot volumes. Permeability was derived from nanoparticle diffusion using the effective medium theory, and showed one-to-one correlation with measured permeability. This technique measured permeability without quantifying gel structure, and may therefore prove useful for characterizing similar materials (e.g., extracellular matrix) where structure is uncontrolled during polymerization and difficult to measure subsequently. We also report that PEGylation reduced, but did not eliminate, the population of immobile particles. We studied the forces required to pull stuck PEG particles free to confirm that the attachment is a result of neither covalent nor strong electrostatic binding, and discuss the relevance of this force scale to particle transport through physiological clots.     

Solute Size Effects on the Diffusion in Biofilms of Streptococcus mutans

The diffusion of poly(ethylene glycol) (PEG) (MW varying between 200 and 10,000), and of three different types of micelles was examined in Streptococcus mutans biofilms using infrared spectroscopy. PEGs were used because they show limited interactions with biological materials and their weight can be selected in order to cover a wide range of size. The study showed that a considerable fraction at the base of the biofilm was not accessible to the diffusing solute molecules and this inaccessible fraction was very dependent on the size of the diffusing molecules. In parallel, it was found that the diffusion coefficients of these solutes in the biofilms were less than those in water and this reduction was less pronounced for large macromolecules, an effect proposed to be related to their limited penetration. Triton X-100, a neutral detergent, forms micelles that behave like PEG, suggesting that the behaviour observed for neutral macromolecules can be extrapolated to neutral macroassemblies. However, the diffusion, as well as the penetration of sodium dodecylsulphate micelles (a negatively charged surfactant) and cetylpyridinium chloride micelles (positively charged), in the biofilms appeared to be significantly influenced by electrostatic interactions with biofilm components. The present findings provide useful insights associated with the molecular parameters required to efficiently penetrate bacterial biofilms. The study suggests a rationale for the limited bactericidal power of some antibiotics (the large ones). The restricted accessibility of macromolecules and macroassemblies to biofilms must be examined carefully in order to offer guidelines in the development of novel antibacterial treatments.     

Brownian dynamics simulation of probe diffusion in DNA: effects of probe size, charge and DNA concentration

We have used Brownian dynamics simulation to study probe diffusion in solutions of short chain DNA using our previously developed simulation algorithm. We have examined the effect of probe size, charge, and DNA concentration on the probe diffusion coefficient, with the aim of gaining insight into the diffusion of proteins in a concentrated DNA environment. In these simulations, DNA was modeled as a worm-like chain of hydrodynamically equivalent spherical frictional elements while probe particles were modeled as spheres of given charge and hydrodynamic radius. The simulations allowed for both short range Lennard-Jones interactions and long ranged electrostatic interactions between charged particles. For uncharged systems, we find that the effects of probe size and DNA concentration on the probe diffusion coefficient are consistent with excluded volume models and we interpret our results in terms of both empirical scaling laws and the predictions of scaled particle theory. For charged systems, we observe that the effects of probe size and charge are most pronounced for the smallest probes and interpret the results in terms of the probe charge density. For an ionic strength of 0.1 M we find that, below a critical probe surface charge density, the probe diffusion coefficient is largely independent of probe charge and only weakly dependent on the DNA charge. These effects are discussed in terms of the interactions between the probe and the DNA matrix and are interpreted in terms of both the underlying physics of transport in concentrated solutions and the assumptions of the simulation model.     

Diffusion Regulation in the Vitreous Humor

The efficient treatment of many ocular diseases depends on the rapid diffusive distribution of solutes such as drugs or drug delivery vehicles through the vitreous humor. However, this multicomponent hydrogel possesses selective permeability properties, which allow for the diffusion of certain molecules and particles, whereas others are immobilized. In this study, we perform an interspecies comparison showing that the selective permeability properties of the vitreous are conserved across several mammalian species. We identify the polyanionic glycosaminoglycans hyaluronic acid and heparan sulfate as two key macromolecules that establish this selective permeability. We show that electrostatic interactions between the polyanionic macromolecules and diffusing solutes can be weakened by charge screening or enzymatic glycosaminoglycan digestion. Furthermore, molecule penetration into the vitreous is also charge-dependent and only efficient as long as the net charge of the molecule does not exceed a certain threshold.     

Effects of Hydration on Steric and Electric Charge-Induced Interstitial Volume Exclusion—a Model

The presence of collagen and charged macromolecules like glycosaminoglycans (GAGs) in the interstitial space limits the space available for plasma proteins and other macromolecules. This phenomenon, known as interstitial exclusion, is of importance for interstitial fluid volume regulation. Physical/mathematical models are presented for calculating the exclusion of electrically charged and neutral macromolecules that equilibrate in the interstitium under various degrees of hydration. Here, a central hypothesis is that the swelling of highly electrically charged GAGs with increased hydration shields parts of the neutral collagen of the interstitial matrix from interacting with electrically charged macromolecules, such that exclusion of charged macromolecules exhibits change due to steric and charge effects. GAGs are also thought to allow relatively small neutral, but also charged macromolecules neutralized by a very high ionic strength, diffuse into the interior of GAGs, whereas larger macromolecules may not. Thus, in the model, relatively small electrically charged macromolecules, such as human serum albumin, and larger neutral macromolecules such as IgG, will have quite similar total volume exclusion properties in the interstitium. Our results are in agreement with ex vivo and in vivo experiments, and suggest that the charge of GAGs or macromolecular drugs may be targeted to increase the tissue uptake of macromolecular therapeutic agents.     

Diffusional anisotropy in collagenous tissues: fluorescence imaging of continuous point photobleaching

Molecular transport in avascular collagenous tissues such as articular cartilage occurs primarily via diffusion. The presence of ordered structures in the extracellular matrix may influence the local transport of macromolecules, leading to anisotropic diffusion depending on the relative size of the molecule and that of extracellular matrix structures. Here we present what we believe is a novel photobleaching technique for measuring the anisotropic diffusivity of macromolecules in collagenous tissues. We hypothesized that macromolecular diffusion is anisotropic in collagenous tissues, depending on molecular size and the local organization of the collagen structure. A theoretical model and experimental protocol for fluorescence imaging of continuous point photobleaching was developed to measure diffusional anisotropy. Significant anisotropy was observed in highly ordered collagenous tissues such as ligament, with diffusivity ratios >2 along the fiber direction compared to the perpendicular direction. In less-ordered tissues such as articular cartilage, diffusional anisotropy was dependent on site in the tissue and size of the diffusing molecule. Anisotropic diffusion was also dependent on the size of the diffusing molecule, with greatest anisotropy observed for larger molecules. These findings suggest that diffusional transport of macromolecules is anisotropic in collagenous tissues, with higher rates of diffusion along primary orientation of collagen fibers.     

The reactions of l nm particles of plutonium-238 dioxide and curium-244 dioxide with lung fluid

The reactions of 1 nm particles of plutonium-238 dioxide and curium-244 dioxide with rat lung fluid have been studied both in vitro and in vivo. The plutonium-238 particles are positively charged and combine by electrostatic attraction with the negative pulmonary surfactant which mediates the transfer of radioactivity to the blood. In contrast the curium-244 particles are negative and are assumed to diffuse passively through the alveolar walls. The results emphasise that electrostatic charge is an important factor governing the reaction of 1 nm actinide oxide particles with macromolecules.     

The height of biomolecules measured with the atomic force microscope depends on electrostatic interactions.

In biological applications of atomic force microscopy, the different surface properties of the biological sample and its support become apparent. Observed height differences between the biomolecule and its supporting surface are thus not only of structural origin, but also depend on the different sample-tip and support-tip interactions. This can result in negative or positive contributions to the measured height, effects that are described by the DLVO (Derjaguin, Landau, Verwey, Overbeek) theory. Experimental verification shows that the electrostatic interactions between tip and sample can strongly influence the result obtained. To overcome this problem, pH and electrolyte concentration of the buffer solution have to be adjusted to screen out electrostatic forces. Under these conditions, the tip comes into direct contact with the surface of support and biological system, even when low forces required to prevent sample deformation are applied. In this case, the measured height can be related to the thickness of the native biological structure. The observed height dependence of the macromolecules on electrolyte concentration makes it possible to estimate surface charge densities.     

Direct force measurements between cellulose surfaces and colloidal silica particles

The interaction of cellulose layers with colloidal silica particles was investigated by direct force measurements with the atomic force microscope (AFM). Upon approach, repulsive forces were found between the negatively charged silica particles and the cellulose surface. The forces were interpreted quantitatively in terms of electrostatic interactions due to overlap of diffuse layers originating from negatively charged carboxylic groups on the cellulose surface. The diffuse layer charge density of cellulose was estimated to be 0.80 mC/m2 at pH 9.5 and 0.21 mC/m2 at pH 4. The forces upon retraction are characterized by molecular adhesion events, whereby individual cellulose chains desorb from the probe surface. The retraction profiles are dominated by well-defined force plateaus, which correspond to single-chain desorption forces of 35-42 pN. We surmise that adsorption of cellulose to probe surfaces is dominated by nonelectrostatic forces, probably originating from hydrogen bonding. Electrostatic contributions to desorption force could be detected only at high pH, where the silica surface is highly charged.     

Enhanced viscoelasticity of human cystic fibrotic sputum correlates with increasing microheterogeneity in particle transport

Current biochemical characterizations of cystic fibrosis (CF) sputum do not address the high degree of microheterogeneity in the rheological properties of the mucosal matrix and only provide bulk-average particle diffusion coefficients. The viscoelasticity of CF sputum greatly reduces the diffusion rates of colloidal particles, limiting the effectiveness of gene delivery to underlying lung cells. We determine diffusion coefficients of hundreds of individual amine-modified and carboxylated polystyrene particles (diameter 100-500 nm) embedded in human CF sputum with 5 nm and 33 ms of spatiotemporal resolution. High resolution multiple particle tracking is used to calculate the effective viscoelastic properties of CF sputum at the micron scale, which we relate to its macroscopic viscoelasticity. CF sputum microviscosity, as probed by 100- and 200-nm particles, is an order of magnitude lower than its macroviscosity, suggesting that nanoparticles dispersed in CF sputum are transported primarily through lower viscosity pores within a highly elastic matrix. Multiple particle tracking provides a non-destructive, highly sensitive method to quantify the high heterogeneity of the mucus pore network. The mean diffusion coefficient becomes dominated by relatively few but fast-moving particles as particle size is reduced from 500 to 100 nm. Neutrally charged particles with a diameter <200 nm undergo more rapid transport in CF sputum than charged particles. Treatment with recombinant human DNase (Pulmozyme) reduces macroviscoelastic properties of CF sputum by up to 50% and dramatically narrows the distribution of individual particle diffusion rates but surprisingly does not significantly alter the ensemble-average particle diffusion rate.     

Intranasal M Cell Uptake of Nanoparticles Is Independently Influenced by Targeting Ligands and Buffer Ionic Strength

In mucosal tissues, epithelial M cells capture and transport microbes across the barrier to underlying immune cells. Previous studies suggested that high affinity ligands targeting M cells may be used to deliver mucosal vaccines; here, we show that particle composition and dispersion buffer ionic strength can independently influence their uptake in vivo. First, addition of a poloxamer 188 to nanoparticle formulations increased uptake of intranasally administered nanoparticles in vivo, but the effect was dependent on the presence of the M cell-targeting ligand. Second, solvent ionic strength is known to effect electrostatic interactions; accordingly, reduced ionic strength increased the electrostatic potential between the epithelium and the particles. Interestingly, below a critical ionic strength, intranasal particle uptake in vivo significantly was increased even when controlled for osmolarity. Similar results were obtained for uptake of bacterial particles. Surprisingly, at low ionic strength, the specific enhancement effect by the targeting peptide was negligible. Modeling of the electrostatic forces predicted that the enhancing effects of the M cell-targeting ligand only are enabled at high ionic strength, as particle electrostatic forces are reduced through Debye screening. Thus, electrostatic forces can have a dramatic effect on the in vivo M cell particle uptake independent of the action of targeting ligands. Examination of these forces will be helpful to optimizing mucosal vaccine and drug delivery.     

Prediction of charge-induced molecular alignment of biomolecules dissolved in dilute liquid-crystalline phases

Alignment of macromolecules in nearly neutral aqueous lyotropic liquid-crystalline media such as bicelles, commonly used in macromolecular NMR studies, can be predicted accurately by a steric obstruction model (Zweckstetter and Bax, 2000). A simple extension of this model is described that results in improved predictions for both the alignment orientation and magnitude of protein and DNA solutes in charged nematic media, such as the widely used medium of filamentous phage Pf1. The extended model approximates the electrostatic interaction between a solute and an ordered phage particle as that between the solute's surface charges and the electric field of the phage. The model is evaluated for four different proteins and a DNA oligomer. Results indicate that alignment in charged nematic media is a function not only of the solute's shape, but also of its electric multipole moments of net charge, dipole, and quadrupole. The relative importance of these terms varies greatly from one macromolecule to another, and evaluation of the experimental data indicates that these terms scale differently with ionic strength. For several of the proteins, the calculated alignment is sensitive to the precise position of the charged groups on the protein surface. This suggests that NMR alignment measurements can potentially be used to probe protein electrostatics. Inclusion of electrostatic interactions in addition to steric effects makes the extended model applicable to all liquid crystals used in biological NMR to date.     

How does protein phosphorylation regulate photosynthesis?

Phosphorylation of light-harvesting antenna proteins redirects absorbed light energy between reaction centres of photosynthetic membranes. A generally accepted explanation for this is that electrostatic forces drive the more negatively charged, phosphorylated antenna proteins between membrane domains that differ in surface charge. However, structural studies on soluble phosphoproteins indicate that phosphorylated amino acid side chains have specific effects on molecular recognition, by ligand blocking or by intramolecular interactions which alter protein structure. These studies suggest alternative mechanisms for phosphorylation in control of pairwise protein-protein interactions in biological membranes. Thus, in photosynthesis, the surface charge model is only one possible interpretation.     

Charge grouping approaches to calculation of electrostatic forces in molecular dynamics of macromolecules.

Calculation of long-range electrostatic interactions is the most time-consuming step in theoretical simulation of the structure and dynamics of macromolecules. In practice very short cutoff distances are used, which may distort the behavior of the model system. We describe two accurate approaches to calculation of electrostatic forces based on hierarchical grouping of charges into cubes. The first is similar to the O(NlogN) algorithm developed by Barnes, J. and Hut, P., Nature (London) 324, 446-449 (1986), for simulation of a gravitational motion of N bodies. The second approach we formulate for a system with periodic boundary conditions in the nearest image approximation. The calculation of electrostatic interactions and a charge grouping procedure are faster than O(N2). The average inaccuracy in the force introduced by the grouping does not exceed 1%. We describe a small modification of the same approach which makes it suitable for long strongly charged polymers as well. This accurate approach to calculation of electrostatic interactions is illustrated with an example of the dynamics of ions near DNA. Quick equilibration of the ionic distribution is observed during molecular dynamics simulation if electrostatic forces are properly calculated, while the behavior and distribution of ions are less realistic when the conventional cutoff distances are used.