by Robert A. Freitas Jr.
© Copyright 1998, Robert A. Freitas Jr.
Please send comments to
The human body consists of ~7 x 1027 atoms arranged in a highly aperiodic physical structure. Although 41 chemical elements are commonly found in the body's construction (Table 3-1), CHON comprises 99% of its atoms. Fully 87% of human body atoms are either hydrogen or oxygen.
Somatic atoms are generally present in combined form as molecules or ions, not individual atoms. The molecules of greatest nanomedical interest are incorporated into cells or circulate freely in blood plasma or the interstitial fluid. Table 3-2 summarizes the gross molecular contents of the typical human cell, which is 99.5% water and salts, by molecule count, and contains ~5000 different types of molecules. Appendix B lists 261 of the most common molecular and cellular constituents of human blood, and their normal concentrations in whole blood and plasma. This listing is far from complete. The human body is comprised of ~105 different molecular species, mostly proteinsa large but nonetheless finite molecular parts list. By 1997, at least ~104 of these proteins had been sequenced, ~103 had been spatially mapped, and ~7,000 structures (including proteins, peptides, viruses, protein/nucleic acid complexes, nucleic acids, and carbohydrates) had been registered in the Protein Data Bank maintained at Brookhaven National Laboratory . It is likely that the sequences and 3-D or tertiary structures of all human proteins will have been determined by the second decade of the 21st century, given the current accelerating pace of improving technology .
Transporting and sorting such a broad range of essential molecular species will be an important basic capability of many nanomedical systems. The three principal methods for distinguishing and conveying molecules that are most useful in nanomedicine are diffusion transport (Section 3.2), membrane filtration (Section 3.3), and receptor-based transport (Section 3.4). The chapter ends with a brief discussion of binding site engineering (Section 3.5).
Fluidic transfer of material, known as convective-diffusive transport, can occur either by convection due to bulk flow or by diffusion due to Brownian motion. In convective transport, material is carried along fluid streamlines at the mean velocity of the fluid, with a velocity distribution such as that in Poiseuille flow (Section 9.4.1.X). Bulk flow is customarily regarded as the most important physiological transport mechanism in the human circulation. Only for the smallest molecules, such as water or glucose, does the time required to diffuse across the width of a capillary roughly equal the time taken by a fluid element to flow the same distance (~0.02 sec). Larger molecules such as fibrinogen take ~100 times longer (~2 sec) to diffuse across one capillary width.
However, bulk flow in the body is usually laminar. Transported materials travel parallel to (and thus cannot reach) fluid/solid interfaces such as the surfaces of blood vessels or membranes. Wall interactions are made possible by diffusion, a random process in which particles can move transversely to fluid streamlines in response to molecular-scale collisions.
Additionally, the movement of micron-scale devices within a bulk fluid flow is dominated by viscous, not inertial, forces (Section 9.4.1.X). Molecular transport to and from such nanodevices is governed by diffusion, not by bulk flow.
A particle suspended in a fluid is subjected to continuous collisions, from all directions, with the surrounding molecules. If the velocities of all molecules were the same all the time, the particle would experience no net movement. However, molecules do not have a single velocity at a given temperature, but rather have a distribution of velocities of varying degrees of probability. Thus from time to time, a suspended particle receives a finite momentum of unpredictable direction and magnitude. The velocity vector of the particle changes continuously, resulting in an observable random zigzag motion, called Brownian movement.
where k = 1.381 x 10-23 joule/kelvin (K) or 0.01381 zJ/K (Boltzmann constant).1 Particles under bombardment also experience a rotational Brownian motion around randomly oriented axes, with the RMS angle of rotation:
although for < min = M / 15R, where M is particle mass (see below), rotation is ballistic.
In human blood plasma, with = 1.1 centipoise (1.1 x 10-3 kg/m-sec) and T = 310 K, a spherical 1-micron diameter nanodevice (R = 0.5 micron) translates ~1 micron in 1 sec (vbrownian ~ 10-6 m/sec) or ~8 microns (~the width of a capillary) after 77 sec (vbrownian ~ 10-7 m/sec), and rotates once in ~16 sec (min = 2 x 10-8 sec). In the same environment, a rigid 10-nm particle (roughly the diameter of a globular protein) would translate ~8 microns in one second (vbrownian ~ 10-5 m/sec) while rotating ~250 times, due to Brownian motion (min = 2 x 10-12 sec).
The instantaneous thermal velocity over one mean free path (the average distance between collisions) is much higher than the net Brownian translational velocity would suggest. For a particle of mass M = 4/3R3 with mean (working) density , the mean thermal velocity is
At T = 310 K, a spherical 1-micron diameter nanodevice of normal density (e.g. taking ~ H2O = 994.9 kg/m3 to minimize ballasting requirements; Section 10.X.X) has vthermal ~ 5 x 10-3 m/sec; for a spherical 10-nm diameter protein with ~ 1500 kg/m3, vthermal ~ 4 m/sec.
Medical nanodevices will frequently be called upon to absorb some particular material from the external aqueous operating environment. Molecular diffusion presents a fundamental limit to the speed at which this absorption can occur. (Once a block of solution has passed into the interior of a nanodevice, it may be divergently subdivided and transported at ~0.01-1 m/sec along internal pathways of characteristic dimension ~1 micron far faster than the <1 mm/sec diffusion velocity across 1 micron distances; Section 220.127.116.11.)
where J is the number of molecules/sec presented to the entire surface of the device, assumed to be 100% absorbed (but see 4.2.5), D (m2/sec) is the translational Brownian diffusion coefficient for the molecule to be absorbed, and C (molecules/m3) is the steady-state concentration of the molecule far from the device . (Blood concentrations in gm/cm3 from Appendix B are converted to molecules/m3 by multiplying Appendix B figures by (106 x NA/MW), where NA = 6.023 x 1023 molecules/mole (Avogadro's number), MW = molecular weight in gm/mole or daltons.) For rigid spherical particles of radius r, where r >> rH2O, the Einstein-Stokes equation  gives
though this is only an approximation because D varies slightly with concentration, with departure from molecular sphericalness, and other factors.
Measured diffusion coefficients in water for various molecules of physiological interest, converted to 310 K, are in Table 3-3. (Diffusion coefficient data for ionic salts such as NaCl and KCl, which dissociate in water and diffuse as independent ions, are for solvated electrolytes.) A 1-micron (diameter) spherical nanodevice suspended in arterial blood plasma at 310 K, with C = 7.3 x 1022 molecules/m3 of oxygen and D = 2.0 x 10-9 m2/sec, encounters a flow rate of J = 9.2 x 108 molecules/sec of O2 impinging upon its surface. (The same calculation applied to serum glucose yields J = 1.3 x 1010 molecules/sec.) The characteristic time for change mediated by diffusion in a region of size L scales as ~L2/D (Eqn. 3.9, below). Across the diameter of an L = 1 micron nanodevice, small molecules such as glucose diffuse in ~0.001 sec, small proteins like hemoglobin in ~0.01 sec, and virus particles diffuse in ~0.1 sec. (Diffusion coefficients of the same molecules in air at room temperature are a factor of ~60 higher, because air ~ 183 micropoise at 20 °C.)
In blood, the diffusivity of larger particles is significantly elevated because local fluid motions generated by individual red cell rotation lead to greater random excursions of the particles . The effective diffusivity De = D + Dr, where the rotation-induced increase in diffusivity Dr ~ 0.25 Rrbc2, with red cell radius Rrbc ~ 2.8 microns (taken for convenience as a spherical volume equivalent) and a typical blood shear rate ~ 500 sec-1, giving Dr ~10-9 m2/sec in normal whole blood. The elevation of diffusivity caused by red cell stirring is just 50% for O2 molecules. However, for large proteins and viruses the effective diffusivity increases 10-100 times, and the effective diffusivity of particles the size of platelets is a factor of 10,000 higher than for Brownian molecular diffusion.
The diffusion current to the surface of a nanodevice can also be estimated for various nonspherical configurations . For instance, the diffusion current to both sides of an isolated thin disk of radius R is given by J = 8RDC. The two-sided current to a square thin plate of area L2 is J = (8/1/2) LDC. The steady-state diffusion current to an isolated cylinder of length Lc and radius R is approximated by J = 2Lc DC/(ln (2Lc/R) - 1), for Lc >> R. The diffusion current through a circular hole of radius R in an impermeable wall separating regions of concentration c1 and c2 is J = 4RD (c1-c2).
Translational Brownian Diffusion Coefficients for
Foraging nanodevices operating in aqueous environments may only modestly exceed the maximum rates of passive diffusive intake described in Section 3.2.2 by engaging in active physical movements designed to increase access to the desired molecules.
The first strategy for active diffusive intake is local stirring. For this, the nanodevice is equipped with suitable active appendages used to manipulate the fluid in its vicinity. Transport by stirring is characterized by a velocity va, the speed of the appendage, and by a length La, its distance of travel, which together define a characteristic stirring frequency stir ~ va/La sec-1. Movement of molecules over a distance La by diffusion alone is scaled by a characteristic time ~La2/D (Section 3.2.2), which defines a characteristic diffusion frequency diff ~ D/La2 sec-1. Stirring will be more effective than diffusion only if stir > diff, that is, if va > D/La. For local stirring, La cannot be much larger than the size of the nanodevice itself. Assuming La = 1 micron and D = 10-9 m2/sec for small molecules, then va > 1000 microns/sec, a faster motion than is exhibited by bacterial cells but quite modest for nanomechanical devices (Section 9.3.1). With D = 10-11 m2/sec for large proteins and virus particles, va > 10 microns/sec, well within the normal microbiological range.
provides a dimensionless measure of the effectiveness of stirring vs. diffusion. For bacteria absorbing small molecules, NSh ~ 10-2. Micron-scale nanodevices with 1-micron appendages capable of 0.01-1 m/sec movement can achieve NSh ~ 10-1000 for small to large molecules, hence could be considerably more effective stirrers.
In a classic paper, Berg and Purcell  analyzed the viscous frictional energy cost of moving the stirring appendages so that the fluid surrounding a spherical object (e.g. a nanodevice) of radius R, out to some maximum stirring radius Rs, is maintained approximately uniform in concentration. The objective is to transfer fluid from a distant region of relatively high concentration to a place much closer to the nanodevice, thereby increasing the concentration gradient near the absorbing surface. To double the passive diffusion current by stirring, the minimum required power density
If = 1.1 x 10-3 kg/m-sec, R = 0.5 micron, D = 10-9 m2/sec for small molecules, and using a modest La = 1 micron stirring apparatus giving Rs = 3R, then Pd ~ 3 x 107 watts/m3. This greatly exceeds the 102-106 watts/m3 power density commonly available to biological cells (Table 6-9) but lies well within the normal range for nanomechanical systems which typically operate at up to ~109 watts/m3. (Nanomedically safe in vivo power densities are discussed at length in Sections 6.5.2 and 6.5.3.) For D ~ 10-11 m2/sec for large molecules, Pd ~ 3 x 103 watts/m3, which is reasonable even by biological standards. The maximum possible gain from stirring is ~ Rs/R, because the current is ultimately limited to what can diffuse into the stirred region.
Local heating due to stirring is minor. Given device volume V ~ 1 micron3, Pd = 3 x 107 watts/m3, mixing distance Lmix ~ 5 microns, and thermal conductivity Kt = 0.623 watts/m-K for water, then T ~ (PdV/LmixKt) = 10 microkelvins; taking heat capacity CV = 4.19 x 106 J/m3-K for water, thermal equilibration time tEQ ~ Lmix2CV/Kt = 0.2 millisec.
The second strategy for active diffusive intake is by swimming. Again, the nanodevice is equipped with suitable active propulsion equipment (Section 9.4) which enable it to move so as to continuously encounter the highest possible concentration gradient near its surface. Consider a spherical motile nanorobot of radius R propelled at constant velocity vswim through a fluid containing a desired molecule for which the surface of the device is essentially a perfect sink (Section 4.2.5). Applying the Stokes velocity field flow around the sphere to the standard diffusion equation, a numerical solution by Berg and Purcell  found that the fractional increase in the diffusion current due to swimming is proportional to vswim2 for vswim << D/R, and to vswim1/3 for vswim >> D/R.
Diffusive intake is doubled at a swimming speed vswim = 2.5 D/R, which for 1-micron devices is ~5000 microns/sec when absorbing small molecules, ~50 microns/sec for large molecules. The viscous frictional energy cost to drive the nanodevice through the fluid, derived from Stokes' law, requires an onboard power density of
If = 1.1 x 10-3 kg/m-sec, vswim = 2.5D/R, R = 0.5 micron, D = 10-9 m2/sec for small molecules, then Pd ~ 5 x 10^5 watts/m3. For large molecules with D = 10-11 m2/sec, Pd ~ 50 watts/m3. Thus the energy cost of diffusive swimming appears modest for nanomechanical systems; gains in diffusion by swimming for nanodevices will be restricted primarily by the maximum safe velocity that may be employed in vivo (Section 9.4.X).
In general, outswimming diffusion requires movement over a characteristic distance Ls ~ D/vswim . For bacteria moving at ~30 micron/sec and absorbing small molecules, then Ls ~ 30 microns, roughly the sprint distance exhibited by flagellar microbes such as E. coli. For micron-scale nanodevices moving at ~1 cm/sec (Section 9.4.X), Ls ~ 1-100 nm for large to small molecules.
Nanodevices may also use diffusion to sort molecules. One of the remarkable features of diffusive sortation is that an input sample consisting of a complex mixture of many different molecular species can sometimes be completely resolved into pure fractions without having any direct knowledge of the precise shapes or electrochemical characteristics of the molecules being sorted. This can be a tremendous advantage for nanodevices operating in environments containing a large number of unknown substances. Another major advantage is the ability to readily distinguish isomeric (though not chiral) molecules. As one example of many possible, molecules suspended in water will diffuse into an adjacent region of pure water at different speeds, giving rise to dissimilar time-dependent concentration gradients which may be exploited for sortation by interrupting the process before complete diffusive equilibrium is reached.
For simplicity, assume we wish to separate two molecular species initially present in solution in equal concentrations (c1 = c2), but having unequal diffusion coefficients (D1 < D2). Consider a separation apparatus with two chambers. Chamber A contains input sample concentrate. Chamber B contains pure water. A dilating gate (Section 3.3.2) separates the two chambers. The gate is opened for a time t approximated by
which relates the diffusion coefficient to the mean displacement X, taken here as L, the length of Chamber B. Table 3-4 gives an estimate of the time required for diffusion to reach 90% completion for glycine, a typical small molecule, in aqueous solution.
After t has elapsed, the gate is closed. (A gate with 10-nm sliding segments moving at 10 cm/sec closes in 0.1 microsec.) The faster-diffusing component D2 approaches diffusive equilibrium in Chamber B, but the slower-diffusing component does not; it is present only in smaller amounts. This gives a separation factor c2/c1 ~ D2/D1 for each diffusion sortation unit. If n units are connected in series, with each unit receiving as input the output of the previous unit, the net concentration achieved by the entire cascade is ~(D2/D1)n. Such cascades are commonplace in gaseous diffusion isotope separation  and other applications.
Figure 3-1 shows a 2-dimensional representation of an efficient design for a simple diffusion unit that might be used in a sortation cascade. Each unit consists of five chambers of equal volume, 7 dilating gates, 3 flap valves, 3 pistons, and two sieves which pass only water (or smaller) molecules. Each chamber is roughly cubical with L ~ 35 nm along the inside edge; including full piston throws and drives, controls, interunit piping and other support structures, each unit measures ~125 nm x 100 nm x 8 0 nm or ~0.001 micron3 with a mass of ~10-18 kg.
The following is a precise description of one complete cycle of operation for each unit:
Increasingly purified sample passes through a multi-unit sortation cascade as described above. For small molecules, a cascade of n ~ 1000 units (total device volume ~1 micron3) completely resolves two mixed molecular species with D2/D1 = 1.01. As a crude approximation, D ~ 1/MW1/3 for small spherical particles , so this cascade separates small molecules differing by the mass of one hydrogen atom, which should be sufficient for most purposes. Structural isomeric forms of the same molecule, such as -alanine and -alanine, often have slightly different diffusion coefficients, thus are also easily separable using a diffusion cascade. However, stereoisomeric (chiral) forms cannot be sorted by diffusion through an optically inactive solvent like water.
For large molecules, a 1 million-unit cascade (total device volume ~1000 micron3) provides D2/D1 ~ 1.00001, sufficient to completely separate large molecules differing by the mass of a single carbon atom. The fidelity of such fine resolutions depends strongly upon the ability to hold constant the temperature of the chamber, since D varies directly with temperature (Eqn. 3.5). Device temperature stability will be determined by at least three factors: (1) the accuracy of onboard thermal sensors in measuring T (T/T < 10-6; Section 4.6.1), (2) the rapidity with which the temperature measurement can be taken (10-9 to 10-6 sec; Section 4.6.1), and (3) the time that elapses between the temperature measurement and the end of the diffusive sortation process (which may be of the same order as the gate closing time, ~10-6 sec).
Most of the waste heat is generated in this device by forced water sieving (Section 3.3.1). To remain within biocompatible thermogenic limits (~109 watts/m3), each unit may be cycled once every ~3 millisec, a 0.8% duty cycle of a ~23 microsec sieving stroke. Subject to this restriction, each device would consume ~1 picowatt in continuous operation. A unit presented with a ~0.1 M concentration of small molecules processes ~106 molecules/sec (e.g. ~1 gm/hour of glucose using 1 cm3 of n = 1000-unit cascades), or ~104 molecules/sec for a unit presented with large molecules at ~0.001 M, circulating ~109 molecules/sec of water as working fluid while running at 340 cycles/sec. Additional chamber segments on each unit, combined with more complex diffusion circuits among the many units in a cascade, should permit the simultaneous complete fractionation of the input feedstock even if hundreds of distinct molecular species are present.
Nanoscale centrifuges offer yet another method for rapid molecular sortation, by biasing diffusive forces with a strong external field. The well-known effect of gravitational acceleration on spherical particles suspended in a fluid is described by Stokes' Law for sedimentation:
where vt is terminal velocity, g is the acceleration of gravity (9.81 m/sec2), R is particle radius, particle and fluid are the particle and fluid densities (kg/m3), and is coefficient of viscosity of the fluid. Particles which are more dense than the suspending liquid tend to fall. Those which are less dense tend to rise (particle/fluid ~ 0.8 for lipids, up to ~1.5 for proteins, and ~1.6 for carbohydrates).
where c2 is the concentration at distance r2 from the axis of a spinning centrifuge (molecules/m3), c1 is the concentration at distance r1 (nearer the axis), MWkg is the molecular weight of the desired molecule in kg/mole, is the angular velocity of the vessel (rad/sec), T is temperature (K) and the universal gas constant Rg = 8.31 joule/mole-K. The approximate spinning time ts required to reach equilibrium is
where Sd is the sedimentation coefficient, usually given in units of 10-13 sec or svedbergs (Table 3-5). Research ultracentrifuges have reached accelerations of ~109 g's.
Consider a cylindrical diamondoid vessel of density vessel = 3510 kg/m3, radius rc = 200 nm, height h = 100 nm, and wall thickness xwall = 10 nm, securely attached to an axial drive shaft of radius ra = 50 nm (schematic in Figure 3-2). A fluid sample containing desired molecules enters the vessel through a hollow conduit in the drive shaft, and the device is rapidly spun. If rim speed vr = 1000 m/sec (max), then = vr/rc = 5 x 109 rad/sec (/2 = 8 x 108 rev/sec). The maximum bursting force Fb ~ 0.5 vessel vr2 = 2 x 109 N/m2, well below the 50 x 109 N/m2 diamondoid tensile strength conservatively assumed by Drexler . Since Sd ranges from 0.1-200 x 10-13 sec for most particles of nanomedical interest (Table 3-5), minimum separation time using acceleration ar / g = vr2/g rc = 5 x 1011 g's, when r2 = rc and r1 = ra, is ts = 0.003-6.0 x 10-6 sec. Fluid sample components migrate at ~0.1 m/sec.
Maximum centrifugation energy per particle Ec = (MWkg / NA) ar (rc-ra) ~ 10,000 zJ/molecule, or ~10 zJ/bond for proteins, well below the 180-1800 zJ/bond range for covalent chemical bonds (Section 3.5.1). However, operating the nanocentrifuge at peak speed may disrupt the weakest noncovalent bonds (including hydrophobic, hydrogen, and van der Waals) which range from 4-50 zJ/bond. The nanocentrifuge has mass ~10-17 kg, requires ~3 picojoules to spin up to speed (bearing drag consumes ~10 nanowatts of power, and fluid drag through the internal plumbing contributes another ~ 5 nanowatts), completes each separation cycle in ~104 revs (~10-5 sec), and processes ~300 micron3/sec which is ~1013 small molecules/sec (at 1% input concentration) or ~109 large molecules/sec (at 0.1% input concentration).
From Eqn. 3.11, the nanocentrifuge separates salt from seawater with c2/c1 ~ 300 across the width of the vessel (rc - ra = 150 nm); extracting glucose from water at 310 K, c2/c1 ~105 over 150 nm. For proteins with particle ~ 1500 kg/m3, separation product removal ports may be spaced, say, 10 nm apart along the vessel radius while maintaining c2/c1 ~ 103 between each port. Vacuum isolation of the unit in an isothermal environment and operation in continuous-flow mode could permit exchange of contents while the vessel is still moving, sharply reducing remixing, vibrations, and thermal convection currents between product layers. A complete design specification of product removal ports, batch and continuous flow protocols, compression profiles, etc. is beyond the scope of this book.
Variable gradient density centrifugation may be used to trap molecules of a specific density in a specific zone for subsequent harvesting, allowing recovery of each molecular species from complex mixtures of substances that are close in density. The traditional method is a series of stratified layers of sucrose or cesium chloride solutions that increase in density from the top to the bottom of the tube. A continuous density gradient may also be used, with the density of the suspension fluid calibrated by physical compression. For example, the coefficient of isothermal compressibility = - (Vl / Vl ) / Pl = ( fluid / fluid) / Pl = 4.492 x 10-5 atm-1 for water at 1 atm and 310 K (compressibility is pressure- and temperature-dependent). Applying Pl = 12,000 atm to the vessel raises fluid density to 1250 kg/m3 , sufficient to partially regulate protein zoning. A multistage cascade (Section 3.2.4) may be necessary for complete compositional separation. Protein denaturation between 5000-15,000 atm  due to hydrogen bond disruption may limit nanocentrifugation rotational velocity. Protein compressibility may further reduce separability. The balance between the differential densities and the differential compressibilities will determine the equilibrium radius of the protein in the centrifuge; in the limiting case of equal compressibilities for a given target protein and water, there is no stable equilibrium radius.
The nanocentrifuge may also be useful in isotopic separations. For a D2O/H2O mixture, c2/c1 = 1.415 per pass through the device; c2/c1 = 106 is achieved in a 40-unit cascade. Tracer glycine containing one atom of 14C is separated from natural glycine using a 113-unit cascade, achieving c2/c1 = 106.
Foresight materials on the Web are ©1986–2013 Foresight Institute. All rights reserved. Legal Notices.