Brownian dynamics simulations of probe and self-diffusion in concentrated protein and DNA solutions.

We have developed a Brownian dynamics algorithm for simulating probe and self-diffusion in concentrated solutions of DNA and protein. In these simulations, proteins are represented as spheres with radii given by their hydrodynamic radii, while DNA is modeled as a wormlike chain of hydrodynamically equivalent spherical frictional elements. The molecular interaction potentials employed by the program allow for intramolecular stretching and bending motions of the DNA chains, short-range Lennard-Jones interactions, and long-range electrostatic interactions. To test the program, we have carried out simulations of bovine serum albumin (BSA) probe diffusion and DNA self-diffusion in solutions of short-chain DNA as a function of both DNA concentration and solution ionic strength. In addition, we report on simulations of BSA self-diffusion as a function of BSA concentration and ionic strength. Based on a comparison to available experimental data, we find that our simulations accurately predict these transport properties under conditions of physiological salt concentration and predict the stronger concentration dependence observed at lower salt concentrations. These results are discussed in light of the nature of the intermolecular interactions in such systems and the approximations and limitations of the simulation algorithm.     

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.     

Computer modeling demonstrates that electrostatic attraction of nucleosomal DNA is mediated by histone tails

We conducted molecular dynamics computer simulations of charged histone tail-DNA interactions in systems mimicking nucleosome core particles (NCP) . In a coarse-grained model, the NCP is modeled as a negatively charged spherical particle with flexible polycationic histone tails attached to it in a dielectric continuum with explicit mobile counterions and added salt. The size, charge, and distribution of the tails relative to the core were built to mimick real NCP. In this way, we incorporate attractive ion-ion correlation effects due to fluctuations in the ion cloud and the attractive entropic and energetic tail-bridging effects. In agreement with experimental data, increase of monovalent salt content from salt-free to physiological concentration leads to the formation of NCP aggregates; likewise, in the presence of MgCl₂, the NCPs form condensed systems via histone-tail bridging and accumulation of counterions. More detailed mechanisms of the histone tail-DNA interactions and dynamics have been obtained from all-atom molecular dynamics simulations (including water), comprising three DNA 22-mers and 14 short fragments of the H4 histone tail (amino acids 5–12) carrying three positive charges on lysine⁺ interacting with DNA. We found correlation of the DNA-DNA distance with the presence and association of the histone tail between the DNA molecules.     

Interactions and Diffusion in Fine-Stranded β-lactoglobulin Gels Determined via FRAP and Binding

The effects of electrostatic interactions and obstruction by the microstructure on probe diffusion were determined in positively charged hydrogels. Probe diffusion in fine-stranded gels and solutions of β-lactoglobulin at pH 3.5 was determined using fluorescence recovery after photobleaching (FRAP) and binding, which is widely used in biophysics. The microstructures of the β-lactoglobulin gels were characterized using transmission electron microscopy. The effects of probe size and charge (negatively charged Na₂-fluorescein (376Da) and weakly anionic 70kDa FITC-dextran), probe concentration (50 to 200 ppm), and β-lactoglobulin concentration (9% to 12% w/w) on the diffusion properties and the electrostatic interaction between the negatively charged probes and the positively charged gels or solutions were evaluated. The results show that the diffusion of negatively charged Na₂-fluorescein is strongly influenced by electrostatic interactions in the positively charged β-lactoglobulin systems. A linear relationship between the pseudo-on binding rate constant and the β-lactoglobulin concentration for three different probe concentrations was found. This validates an important assumption of existing biophysical FRAP and binding models, namely that the pseudo-on binding rate constant equals the product of the molecular binding rate constant and the concentration of the free binding sites. Indicators were established to clarify whether FRAP data should be analyzed using a binding-diffusion model or an obstruction-diffusion model.     

DNA complexing with polyamidoamine dendrimers: implications for transfection.

DNA and polyamidamine (PAMAM) dendrimers form complexes on the basis of the electrostatic interactions between negatively charged phosphate groups of the nucleic acid and protonated (positively charged) amino groups of the polymers. Charge neutralization of both components and subsequent increases of the net positive charge of the complex result in changes in the physicochemistry and biological properties of the complexes. The formation of soluble, low-density and insoluble, high-density complexes was analyzed using UV light absorption and measurements of radioactive labeled DNA. Formation of high molecular weight and high-density complexes depended mainly on the DNA concentration and was enhanced by increasing the dendrimer-DNA charge ratio. Electrostatic charge related effects (attraction or repulsion of charged particles) appeared to be modulated by the generation of dendrimer (size of the polymer). With the progressive increases in the dendrimer-DNA charge ratio (above 20), an increase in the amount of low-density, soluble complexes was observed. Functional analysis revealed that the great majority (>90%) of transfection is carried by low-density, soluble, complexes which only represent approximately 10-20% of total complexed DNA. The ability of the dendrimer to complex and form aggregates with DNA is crucial for efficient transfection and the function of the complexed DNA.     

Probe diffusion in aqueous poly(vinyl alcohol) solutions studied by fluorescence correlation spectroscopy

We report fluorescence correlation spectroscopy measurements of the translational diffusion coefficient of various probe particles in dilute and semidilute aqueous poly(vinyl alcohol) solutions. The range of sizes of the particles (fluorescent molecules, proteins, and polymers) was chosen to explore various length scales of the polymer solutions as defined by the polymer-polymer correlation length. For particles larger than the correlation length, we find that the diffusion coefficient, D, decreases exponentially with the polymer concentration. This can be explained by an exponential increase in the solution viscosity, consistent with the Stokes-Einstein equation. For probes on the order of the correlation length, the decrease of the diffusion coefficient cannot be accounted for by the Stokes-Einstein equation, but can be fit by a stretched exponential, D approximately exp(-alphacn), where we find n = 0.73-0.84 and alpha is related to the probe size. These results are in accord with a diffusion model of Langevin and Rondelez (Polymer 1978, 19, 1875), where these values of n indicate a good solvent quality.     

Protein Diffusion on Charged Biopolymers: DNA versus Microtubule

Protein diffusion in lower-dimensional spaces is used for various cellular functions. For example, sliding on DNA is essential for proteins searching for their target sites, and protein diffusion on microtubules is important for proper cell division and neuronal development. On the one hand, these linear diffusion processes are mediated by long-range electrostatic interactions between positively charged proteins and negatively charged biopolymers and have similar characteristic diffusion coefficients. On the other hand, DNA and microtubules have different structural properties. Here, using computational approaches, we studied the mechanism of protein diffusion along DNA and microtubules by exploring the diffusion of both protein types on both biopolymers. We found that DNA-binding and microtubule-binding proteins can diffuse on each other’s substrates; however, the adopted diffusion mechanism depends on the molecular properties of the diffusing proteins and the biopolymers. On the protein side, only DNA-binding proteins can perform rotation-coupled diffusion along DNA, with this being due to their higher net charge and its spatial organization at the DNA recognition helix. By contrast, the lower net charge on microtubule-binding proteins enables them to diffuse more quickly than DNA-binding proteins on both biopolymers. On the biopolymer side, microtubules possess intrinsically disordered, negatively charged C-terminal tails that interact with microtubule-binding proteins, thus supporting their diffusion. Thus, although both DNA-binding and microtubule-binding proteins can diffuse on the negatively charged biopolymers, the unique molecular features of the biopolymers and of their natural substrates are essential for function.     

Interactions and Diffusion in Fine-Stranded β-lactoglobulin Gels Determined via FRAP and Binding

The effects of electrostatic interactions and obstruction by the microstructure on probe diffusion were determined in positively charged hydrogels. Probe diffusion in fine-stranded gels and solutions of β-lactoglobulin at pH 3.5 was determined using fluorescence recovery after photobleaching (FRAP) and binding, which is widely used in biophysics. The microstructures of the β-lactoglobulin gels were characterized using transmission electron microscopy. The effects of probe size and charge (negatively charged Na₂-fluorescein (376Da) and weakly anionic 70kDa FITC-dextran), probe concentration (50 to 200 ppm), and β-lactoglobulin concentration (9% to 12% w/w) on the diffusion properties and the electrostatic interaction between the negatively charged probes and the positively charged gels or solutions were evaluated. The results show that the diffusion of negatively charged Na₂-fluorescein is strongly influenced by electrostatic interactions in the positively charged β-lactoglobulin systems. A linear relationship between the pseudo-on binding rate constant and the β-lactoglobulin concentration for three different probe concentrations was found. This validates an important assumption of existing biophysical FRAP and binding models, namely that the pseudo-on binding rate constant equals the product of the molecular binding rate constant and the concentration of the free binding sites. Indicators were established to clarify whether FRAP data should be analyzed using a binding-diffusion model or an obstruction-diffusion model.     

Transport of nucleosome core particles in semidilute DNA solutions

We studied the diffusion of native and trypsinized nucleosome core particles (NCPs), in aqueous solution and in concentrated DNA solutions (0.25-100 mg/ml) using fluorescence correlation spectroscopy (FCS). The highest DNA concentrations studied mimic the DNA density inside the cell nucleus. The diffusion coefficient of freely diffusing NCPs depends on the presence or absence of histone tails and is affected by the salt concentration due to the relaxation effect of counterions. NCPs placed in a network of long DNA molecules (30-50 kbp) reveal anomalous diffusion. We demonstrate that NCPs diffusion is in agreement with known particle transport in entangled macromolecular solutions as long as the histone tails are folded onto the particles. In contrast, when these tails are unfolded, the reversible adsorption of NCPs onto the DNA network has to be taken into account. This is confirmed by the fact that removal of the tails leads to reduction of the interaction between NCPs and the DNA network. The findings suggest that histone tail bridging plays an important role in chromatin dynamics.     

Diffusion of Particles in the Extracellular Matrix: The Effect of Repulsive Electrostatic Interactions

Diffusive transport of macromolecules and nanoparticles in charged fibrous media is of interest in many biological applications, including drug delivery and separation processes. Experimental findings have shown that diffusion can be significantly hindered by electrostatic interactions between the diffusing particle and charged components of the extracellular matrix. The implications, however, have not been analyzed rigorously. Here, we present a mathematical framework to study the effect of charge on the diffusive transport of macromolecules and nanoparticles in the extracellular matrix of biological tissues. The model takes into account steric, hydrodynamic, and electrostatic interactions. We show that when the fiber size is comparable to the Debye length, electrostatic forces between the fibers and the particles result in slowed diffusion. However, as the fiber diameter increases the repulsive forces become less important. Our results explain the experimental observations that neutral particles diffuse faster than charged particles. Taken together, we conclude that optimal particles for delivery to tumors should be initially cationic to target the tumor vessels and then change to neutral charge after exiting the blood vessels.     

Distribution and diffusivity of a hydrophobic probe molecule in the interior of a membrane: theory and simulation.

We propose a simple model for the distribution of position and orientation and the diffusion of a hydrophobic probe molecule embedded in a membrane. The molecule experiences both a Maier-Saupe orienting potential as well as an enclosing potential of repulsion from the membrane walls. A statistical thermodynamics treatment of the model provides predictions of the location and orientation of the molecule within the membrane. In particular, we evaluate the order parameter of the molecule in terms of the model constants. The diffusivity of the probe is studied by Brownian dynamics simulation. For rotational diffusion, we check an available analytical approximate treatment that allows for the prediction of the dynamics in terms of equilibrium quantities. We also pay attention to quantities related to the initial and mean reorientational rate of the probe. For translational diffusion, we use the simulation results to analyze some general aspects of lateral and transversal diffusion.     

Mass Transfer in the Cornea: I. Interacting Ion Flows in an Arbitrarily Charged Membrane

The formalisms of irreversible thermodynamics are used to describe multi-ionic nonconvective flow through an arbitrarily charged membrane. Interactions between oppositely charged ions are included and are measured by a single phenomenological coefficient. The consequent generalized Nernst-Planck flux equations are integrated to yield a relation between the species fluxes and the composition of the solutions bounding the membrane. It is assumed in the derivation that activity coefficient gradients within the membrane and direct interactions between ions of like charge are negligible. Some special cases are examined. To illustrate the use of the final equations, a single membrane separating solutions of differing composition is modeled, and the effect of ion-ion interactions on the membrane potential and the ion fluxes is demonstrated for several values of diffusion current density and membrane charge density.     

Kinetics of oligonucleotide hybridization to DNA probe arrays on high-capacity porous silica substrates

We have investigated the kinetics of DNA hybridization to oligonucleotide arrays on high-capacity porous silica films that were deposited by two techniques. Films created by spin coating pure colloidal silica suspensions onto a substrate had pores of approximately 23 nm, relatively low porosity (35%), and a surface area of 17 times flat glass (for a 0.3-microm film). In the second method, latex particles were codeposited with the silica by spin coating and then pyrolyzed, which resulted in larger pores (36 nm), higher porosity (65%), and higher surface area (26 times flat glass for a 0.3-microm film). As a result of these favorable properties, the templated silica hybridized more quickly and reached a higher adsorbed target density (11 vs. 8 times flat glass at 22 degrees C) than the pure silica. Adsorption of DNA onto the high-capacity films is controlled by traditional adsorption and desorption coefficients, as well as by morphology factors and transient binding interactions between the target and the probes. To describe these effects, we have developed a model based on the analogy to diffusion of a reactant in a porous catalyst. Adsorption values (k(a), k(d), and K) measured on planar arrays for the same probe/target system provide the parameters for the model and also provide an internally consistent comparison for the stability of the transient complexes. The interpretation of the model takes into account factors not previously considered for hybridization in three-dimensional films, including the potential effects of heterogeneous probe populations, partial probe/target complexes during diffusion, and non-1:1 binding structures. The transient complexes are much less stable than full duplexes (binding constants for full duplexes higher by three orders of magnitude or more), which may be a result of the unique probe density and distribution that is characteristic of the photolithographically patterned arrays. The behavior at 22 degrees C is described well by the predictive equations for morphology, whereas the behavior at 45 degrees C deviates from expectations and suggests that more complex phenomena may be occurring in that temperature regime.     

Multiscale Coarse-Grained Approach to Investigate Self-Association of Antibodies

Self-association of therapeutic monoclonal antibodies (mabs) are thought to modulate the undesirably high viscosity observed in their concentrated solutions. Computational prediction of such a self-association behavior is advantageous early during mab drug candidate selection when material availability is limited. Here, we present a coarse-grained (CG) simulation method that enables microsecond molecular dynamics simulations of full-length antibodies at high concentrations. The proposed approach differs from others in two ways: first, charges are assigned to CG beads in an effort to reproduce molecular multipole moments and charge asymmetry of full-length antibodies instead of only localized charges. This leads to great improvements in the agreement between CG and all-atom electrostatic fields. Second, the distinctive hydrophobic character of each antibody is incorporated through empirical adjustments to the short-range van der Waals terms dictated by cosolvent all-atom molecular dynamics simulations of antibody variable regions. CG simulations performed on a set of 15 different mabs reveal that diffusion coefficients in crowded environments are markedly impacted by intermolecular interactions. Diffusion coefficients computed from the simulations are in correlation with experimentally measured observables, including viscosities at a high concentration. Further, we show that the evaluation of electrostatic and hydrophobic characters of the mabs is useful in predicting the nonuniform effect of salt on the viscosity of mab solutions. This CG modeling approach is particularly applicable as a material-free screening tool for selecting antibody candidates with desirable viscosity properties.     

Stability, specificity and fluorescence brightness of multiply-labeled fluorescent DNA probes.

In this work, we studied the fluorescence and hybridization of multiply-labeled DNA probes which have the hydrophilic fluorophore 1-(straightepsilon-carboxypentynyl)-1'-ethyl- 3,3,3', 3'-tetramethylindocarbocyanine-5,5'-disulfonate (Cy3) attached via either a short or long linker at the C-5 position of deoxyuridine. We describe the effects of labeling density, fluorophore charge and linker length upon five properties of the probe: fluorescence intensity, the change in fluorescence upon duplex formation, the quantum yield of fluorescence (Phif), probe-target stability and specificity. For the hydrophilic dye Cy3, we have demonstrated that the fluorescence intensity andPhifare maximized when labeling every 6th base using the long linker. With a less hydrophilic dye, a labeling density this high could not be achieved without serious quenching of the fluorescence. The target specificity of multiply-labeled DNA probes was just as high as compared to the unmodified control probe, however, a less stable probe-target duplex is formed that exhibits a lower melting temperature. A mechanism that accounts for this destabilization is proposed which is consistent with our data. It involves dye-dye and dye-nucleotide interactions which appear to stabilize a single-stranded conformation of the probe.     

On the truncation of long-range electrostatic interactions in DNA

Long-range interactions are known to play an important role in highly polar biomolecules like DNA. In molecular dynamics simulations of nucleic acids and proteins, an accurate treatment of the long-range interactions are crucial for achieving stable nanosecond trajectories. In this report, we evaluate the structural and dynamic effects on a highly charged oligonucleotide in aqueous solution from different long-range truncation methods. Two group-based truncation methods, one with a switching function and one with a force-switching function were found to fail to give accurate stable trajectories close to the crystal structure. For these group-based truncation methods, large root mean square (rms) deviations from the initial structure were obtained and severe distortions of the oligonucleotide were observed. Another group-based truncation scheme, which used an abrupt truncation at 8. 0 A or at 12.0 A was also investigated. For the short cutoff distance, the conformations deviated far away from the initial structure and were significantly distorted. However, for the longer cutoff, where all necessary electrostatic interactions were included, the trajectory was quite stable. For the particle mesh Ewald (PME) truncation method, a stable DNA simulation with a heavy atom rms deviation of 1.5 A was obtained. The atom-based truncation methods also resulted in stable trajectories, according to the rms deviation from the initial B-DNA structure, of between 1.5 and 1.7 A for the heavy atoms. In these stable simulations, the heavy atom rms deviations were approximately 0.6-1.0 A lower for the bases than for the backbone. An increase of the cutoff radius from 8 to 12 A decreased the rms deviation by approximately 0.2 A for the atom-based truncation method with a force-shifting function, but increased the computational time by a factor of 2. Increasing the cutoff from 12 to 18 A for the atom-based truncation method with a force-shifting function requires 2-3 times more computational time, but did not significantly change the rms deviation. Similar rms deviations from the initial structure were found for the atom-based method with a force-shifting function and for the PME method. The computational cost was longer for the PME method with a cutoff of 12. 0 A for the direct space nonbonded calculations than for the atom-based truncation method with a force-shifting function and a cutoff of 12.0 A. If a nonperiodic boundary, e.g., a spherical boundary, was used, a considerable speedup could be achieved. From the rms fluctuations, the terminal nucleotides and especially the cytidines were found to be more flexible than the nonterminal nucleotides. The B-DNA form of the oligonucleotide was maintained throughout the simulations and is judged to depend on the parameters of the energy function and not on the truncation method used to handle the long-range electrostatic interactions. To perform accurate and stable simulations of highly charged biological macromolecules, we recommend that the atom-based force-shift method or the PME method should be used for the long-range electrostatics interactions.     

Translational and rotational motions of proteins in a protein crowded environment

Fluorescence correlation spectroscopy (FCS) was used to measure the translational diffusion of labeled apomyoglobin (tracer) in concentrated solutions of ribonuclease A and human serum albumin (crowders), as a quantitative model system of protein diffusive motions in crowded physiological environments. The ratio of the diffusion coefficient of the tracer protein in the protein crowded solutions and its diffusion coefficient in aqueous solution has been interpreted in terms of local apparent viscosities, a molecular parameter characteristic for each tracer-crowder system. In all protein solutions studied in this work, local translational viscosity values were larger than the solution bulk viscosity, and larger than rotational viscosities estimated for apomyoglobin in the same crowding solutions. Here we propose a method to estimate local apparent viscosities for the tracer translational and rotational diffusion directly from the bulk viscosity of the concentrated protein solutions. As a result of this study, the identification of protein species and the study of hydrodynamic changes and interactions in model crowded protein solutions by means of FCS and time-resolved fluorescence depolarization techniques may be expected to be greatly simplified.     

Weak interactions govern the viscosity of concentrated antibody solutions: high-throughput analysis using the diffusion interaction parameter

Weak protein-protein interactions are thought to modulate the viscoelastic properties of concentrated antibody solutions. Predicting the viscoelastic behavior of concentrated antibodies from their dilute solution behavior is of significant interest and remains a challenge. Here, we show that the diffusion interaction parameter (k(D)), a component of the osmotic second virial coefficient (B(2)) that is amenable to high-throughput measurement in dilute solutions, correlates well with the viscosity of concentrated monoclonal antibody (mAb) solutions. We measured the k(D) of 29 different mAbs (IgG(1) and IgG(4)) in four different solvent conditions (low and high ion normality) and found a linear dependence between k(D) and the exponential coefficient that describes the viscosity concentration profiles (|R| ≥ 0.9). Through experimentally measured effective charge measurements, under low ion normality where the electroviscous effect can dominate, we show that the mAb solution viscosity is poorly correlated with the mAb net charge (|R| ≤ 0.6). With this large data set, our results provide compelling evidence in support of weak intermolecular interactions, in contrast to the notion that the electroviscous effect is important in governing the viscoelastic behavior of concentrated mAb solutions. Our approach is particularly applicable as a screening tool for selecting mAbs with desirable viscosity properties early during lead candidate selection.     

Monoazido analog of ethidium as a chromatin probe: Binding to DNA

Two photoaffinity analogs of ethidium, 8-azido-3-amino, and 3-azido-8-amino-5-ethyl-6-phenylphenanthridinium chloride, have been used to probe the structure of mammalian chromatin and its interactions with the ethidium moiety. The monoazido analogs were established as suitable probes by comparing their interactions with chromatin and pure DNA prepared from chromatin to those of the parent ethidium bromide. Scatchard analysis of the binding data determined from spectrophotometric titrations showed that the analogs interacted with both nucleic acids in a manner similar to the parent compound. The effect of chromatin proteins on the interaction of the ethidium moiety with intact chromatin was investigated directly. By exposing the noncovalent complex to visible light, the monoazido analog was attached covalently in its interaction sites within chromatin, and the amount of drug bound covalently to DNA was determined for both protein-free DNA and chromatin. Using saturating concentrations of drug, DNA within intact chromatin was found to be associated with only half as much drug as DNA extracted from its protein prior to drug exposure. The distribution of drug bound within chromatin was determined following the attachment of the monoazido analog (by photoactivation) to chromatin that had undergone limited nuclease digestion. Several distinct populations isolated by size fractionation and quantitative measurements revealed that (1) both the core particles and the spacer-containing particles contained bound drug, reflecting high-affinity binding sites; and (2) chromatin particles containing 150 DNA base pairs (putatively nucleosome core structures) contained less total bound drug at high drug concentrations than those particles having intact spacer DNA.     

The biophysics of DNA hybridization with immobilized oligonucleotide probes.

A mathematical model based on receptor-ligand interactions at a cell surface has been modified and further developed to represent heterogeneous DNA-DNA hybridization on a solid surface. The immobilized DNA molecules with known sequences are called probes, and the DNA molecules in solution with unknown sequences are called targets in this model. Capture of the perfectly complementary target is modeled as a combined reaction-diffusion limited irreversible reaction. In the model, there are two different mechanisms by which targets can hybridize with the complementary probes: direct hybridization from the solution and hybridization by molecules that adsorb nonspecifically and then surface diffuse to the probe. The results indicate that nonspecific adsorption of single-stranded DNA on the surface and subsequent two-dimensional diffusion can significantly enhance the overall reaction rate. Heterogeneous hybridization depends strongly on the rate constants for DNA adsorption/desorption in the non-probe-covered regions of the surface, the two-dimensional (2D) diffusion coefficient, and the size of probes and targets. The model shows that the overall kinetics of DNA hybridization to DNA on a solid support may be an extremely efficient process for physically realistic 2D diffusion coefficients, target concentrations, and surface probe densities. The implication for design and operation of a DNA hybridization surface is that there is an optimal surface probe density when 2D diffusion occurs; values above that optimum do not increase the capture rate. Our model predicts capture rates in agreement with those from recent experimental literature. The results of our analysis predict that several things can be done to improve heterogeneous hybridization: 1) the solution phase target molecules should be about 100 bases or less in size to speed solution-phase and surface diffusion; 2) conditions should be created such that reversible adsorption and two-dimensional diffusion occur in the surface regions between DNA probe molecules; 3) provided that 2) is satisfied, one can achieve results with a sparse probe coverage that are equal to or better than those obtained with a surface totally covered with DNA probes.