| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Microparticle Photophysics Laboratory, Polytechnic University, Brooklyn, New York 11201
Correspondence: Address reprint requests to Iwao Teraoka, E-mail: teraoka{at}poly.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
A transparent microsphere can accommodate WGM in the vicinity of the sphere surface. The light propagates near the curved surface by total internal reflection. Resonance is achieved when the light path closes upon itself in phase after one cycle. If the diameter of the sphere is sufficiently large compared with the wavelength, the resonance can have a narrow width. The Q value of a silica microsphere in water at a 1.3-µm wavelength can be as large as 2 x 106 (12
). A much greater Q value, exceeding 108, is reported for a toroidal resonator at 680 nm (13
).
In each reflection along the circular path of WGM, the light seeps into the surroundings as an evanescent wave. The wavelength of the sharp resonance is sensitive to small changes in dielectric property in the immediate neighborhood of the transparent microsphere (14
). The changes include adsorption of molecules onto the microsphere and a change of refractive index (RI) in the surrounding medium. The shift of the wavelength upon adsorption of biomolecules onto the microsphere has been heralded as the most sensitive detector ever made possible without the necessity for fluorescent labeling (12
,15
18
). Detection of a single protein molecule is considered within reach (15
). Recent prediction (19
) and demonstration (20
) of enhanced sensitivity by a high RI coating has paved the way for the difficult detection. The sensor's capability is not limited to estimating the surface density of adsorbed molecules. Independent detection of the resonance shifts for two polarization modestransverse electric (TE) and transverse magnetic (TM)is expected to allow us to estimate the orientation of adsorbed anisotropic molecules (21
).
Each WGM is specified by l, m,
, and polarization (22
). l represents the number of waves in a circular orbit, m (= l, l + 1, ..., l) is the azimuthal index, and
is equal to the number of peaks in the radial function of the electric field intensity, thus specifying the radial mode. The polarization is either TE or TM. The wavelength at resonance is determined by l,
, and polarization. In a perfect spherical resonator, modes of different m are degenerate. The shift of resonance wavelength in response to the environmental changes depends also on l,
, and polarization (23
). It was recently demonstrated that the observed shifts of TE modes due to RI changes in the surroundings were in agreement with the shifts calculated using the indices evaluated for the microsphere used (24
).
More than a decade ago, Folan distinguished TE and TM shifts of WGM in a small polystyrene sphere levitated electrodynamically in air (25
). Folan examined the change in the scattering spectrum as water condensed onto the polymer sphere for the two polarizations, but the difference between the two shifts was insignificant within experimental error.
In this report, we use side coupling of a core-exposed single-mode fiber to induce both TE and TM polarizations in a silica microsphere and measure the wavelength shifts when proteins are added to the surrounding fluid to adsorb onto the sphere surface. We find that the shifts are different for TE and TM and the ratio of the two shifts provides information on the state of adsorbed protein. We confirmed the polarization-sensitive shifts by adding NaCl to the surroundings to cause a uniform increase of RI. The latter situation simulates adsorption of large objects such as mitochondria (17
).
Recently developed dual polarization interferometry (DPI) (26
29
) can provide information on the state of surface-bound proteins. DPI uses two polarizations of light transmitted through a pair of planar waveguides to find the RI and thickness of the adsorption layer, which in turn provide information on the protein conformation. To achieve a high sensitivity comparable to that of the surface plasmon resonance (SPR) instrument, DPI uses a large sensor area,
150 mm2. Our WGM sensor has a much smaller sensor area, typically <0.005 mm2, yet easily surpasses the sensitivity of DPI and SPR in terms of adsorbed mass per unit area while retaining the capability to find the state of the adsorbate. More importantly, the WGM sensor allows easy interpretation of the resonance shift in terms of molecular parameters (21
), without the need to assume an adsorption layer of uniform RI and thickness (26
). Neutron reflectivity (7
) is another method that macroscopically characterizes the adsorbed molecules as a whole, but its sensitivity and usefulness are limited.
Theoretical background
A plain microsphere of radius a and uniform relative permittivity
is placed in a uniform medium of
(n1 > n2). When the wavelength,
, of WGM is much shorter than a, the electric field, E(r), of the WGM is mostly confined to the interior of the microsphere. However, the evanescent field extends into the surroundings to the depth of
(
/2
)(
)1/2, which polarizes the molecules in the immediate neighborhood of the microsphere surface. The resonance wavelength of WGM shifts from
0 to
0 + 
, when small molecules (much smaller than
/n2) adsorb onto the sphere surface to displace a part of the surrounding medium. The adsorbed molecules are polarized by E(|r| = a+), where a+ indicates the exterior side of the sphere surface.
The fractional shift 
/
0 is equal to the ratio of the polarization energy in the adsorbed molecules to the total mode energy (12
,15
,21
). In general, the TE and TM modes exhibit different shifts, as the TE mode has only a tangential component, Et(a) as E(a+), whereas the TM mode has also a normal component, En(a+). For uniform adsorption of Np molecules of excess polarizability,
, at a low surface density, the fractional shift of either mode is given as (21
)
![]() | (1) |
tt and
nn are the polarizability tensor components in directions of Et and En, respectively,
...
indicates the average on the sphere surface (for Et and En) or the average with respect to the configuration of adsorbed molecules (for
tt and
nn),
0 is the vacuum permittivity, and the volume integral in the denominator covers the entire space (relative permittivity
r is
r1 in the sphere;
r2 elsewhere).
We showed earlier that the denominator in Eq. 1 is equal to 4
a3
0(
r1
r2)
[Et(a)]2
for the TE mode (21
,23
). Then, the fractional shift of the TE mode is given as
![]() | (2) |
The shift is identical for all radial modes (
= 1, 2, ...). Since Np
a2 for a given surface density of adsorbates, 
/
0
a1. There is a weak dependence of 
TE/
0 on
0 through wavelength dispersions of
tt,
r1, and
r2. For the TM mode, the denominator in Eq. 1 is calculated as 4
a3
0(
r1
r2)(
[Et(a)]2
+ (
r2/
r1)
[En(a+)]2
) (21
,23
). The expression for the TM shift, 
TM, is then obtained. The ratio 
TM/
TE is
![]() | (3) |
. Here, A+
[Et(a)]2
/
[En(a+)]2
is the intensity anisotropy ratio of the evanescent field of the TM mode (right on the sphere surface). When l >> 1, the following approximation is useful (21
![]() | (4) |
1
r2/
r1. The second approximation is not good except for the first-order radial modes (
= 1). The shift ratio given by Eq. 3 is insensitive to the dimension of the adsorbed molecule, as long as it is sufficiently smaller than
0/n2.
We now evaluate Eq. 3 for a molecule of volume Vp and uniform, isotropic relative permittivity
We consider five geometries of the moleculea sphere, a rod (cylinder) standing vertically on the surface, a rod lying on the surface, a disk standing on the surface (edge-on), and a disk lying on the surface (face-on), as illustrated in Fig. 1. For now, we assume a low surface density of adsorbed molecules. The effect of interference from the dipoles induced at nearby adsorbed molecules will be discussed toward the end of this section.
|
nn is definite: 
nn
=
nn. For a sphere, a standing rod and a lying disk,
tt, is also definite. For a lying rod and a standing disk, the orientation of the molecule on the surface relative to Et varies from molecule to molecule. If we assume random orientation of molecules in the directions tangential to the surface, 
tt
is the isotropic mean of the orthogonal tensor components in the directions parallel to the surface.
Table 1 lists 
tt
/(
0Vp) and 
nn
/(
0Vp) for the five geometries. These expressions were obtained from the boundary conditions for the electric field across the surface of the adsorbate. The expressions for the sphere and the two orientations of the rod were obtained earlier (21
). Our early works also showed that the image dipole induced within the sphere does not affect the polarization of the adsorbed particle (21
,30
). For rods lying on the surface, 
tt
/(
0Vp) is the isotropic mean of 2
r2(
rp +
r2)1(
rp
f2) and
rp
r2, which gives the listed expression. For edge-on adsorbed disks, 
tt
/(
0Vp) is the isotropic mean of (
r2/
rp)(
rp
r2) and
rp
r2. A thin uniform layer has the same
tt and
nn as those for isolated disks lying on the surface. The two geometries are listed in the last row of the table.
|
r1 = 1.4522 (silica),
r2 = 1.322 (water), a = 171 µm,
0 = 1.312 µm, l = 1170, and
rp = 1.552, where Eq. 4 was used for A+. Here, the RI of protein at 1.32 µm was estimated as 1.55 from the value of 1.57 at 589 nm (31
r2/
r1. For the other geometries, the ratio is 2
r2/
r1 = 1.18 at np = n2 and deviates from that value with an increasing np. The shift ratio >1.18 indicates
nn >
tt, which is likely due to molecules of an anisotropic shape standing on the surface. The ratio smaller than 1.18 indicates a geometry of the adsorbate extending parallel to the surface. The capability of the WGM sensors to provide information on the molecular orientation, independent of the size of the molecule, will be useful in studies on protein adsorption in different solutions and surface environments as well as conformational changes.
The above discussion applies to low surface coverages. With an increasing surface density, dipolar fields by nearby particles decrease
[En(a+)]2
but increase
[Et(a)]2
(21
). Earlier (30
), we used dipolar approximation to consider the effect for spherical molecules sequentially and randomly adsorbed onto the surface. Since we do not have a formula for geometries other than spheres, we adopt a formula for the spheres. Calculations for spheres of np = 1.55 show that the effect increases Et(a) by 0.9% and decreases En(a+) by 1.8% at 15% coverage of the projection of the spheres onto the surface. As a result, the TE shift is 0.9% greater than the estimate given by Eq. 2, and the TM shift is 1.3% less. The increase in A+ with an increasing surface density causes the ratio to drop to 1.07 at the highest surface coverage in random sequential adsorption (32
34
). Concomitantly, the criterion of the shift ratio for the anisotopic polarizability moves to a smaller value.
| MATERIALS AND METHODS |
|---|
|
|
|---|
WGM resonance shift measurement system
The optical part in our measurement system is shown in Fig. 2. We used a pigtailed butterfly laser (distributed feedback laser; DFB) from NTT Electronics (NLK1B5E1AA; Saddle Brook, NJ) operating at
1.31 µm as a light source. The laser is linearly polarized with an extinction ratio
1000 as observed by a photodiode (PDA400; Thorlabs, Newton, NJ) at the end of the single-mode fiber. Exposing the core of the fiber by etching in hydrofluoric acid solution decreased the extinction ratio to
200. However, neither etching nor contact with the microsphere skewed the polarization. Rotation of the laser mount around the axis of the output fiber changed the polarization direction. The angle of rotation of the polarizer at the photodiode to maximize the intensity of transmitted light was measured as the laser mount was rotated in a step of 10° up to ±90° from the unstrained direction. The required polarizer rotation was nearly identical to the laser rotation. The standard deviation of the difference between the two angles of rotation was 4.8°. The extinction ratio barely changed during the rotation. Thus, we know what angles of the laser mount cause vertical and horizontal polarizations in the etched section of the fiber, which in turn will excite TE and TM modes, respectively, within the microsphere that touches the fiber at its horizontal equator.
|
, was scanned at 10 Hz by changing the laser drive current, i, linearly with time. A sawtooth function generator was used for that purpose. The scan range was
0.2 nm. The relationship between
and i was evaluated using an interferometer (Agilent, Santa Clara, CA; HP3325A). In each scan, the laser intensity and wavelength increase almost linearly with time, and the two increases are nearly parallel to each other. When a microsphere is placed in contact with the core-exposed section of the fiber, destructive interference by the WGM causes dips in the light intensity at the photodetector. Each dip represents a WGM of a unique set of indices (l, m,
, polarization). | RESULTS AND DISCUSSION |
|---|
|
|
|---|
0.042 nm, ascribed to splitting of the degenerate azimuthal modes (m) by a spheroidal shape of the microsphere; Lai et al. predicted polarization-independent splitting by distortion of the meridional cross section of the microspheroid (35
= 1. Our scan range of
0.2 nm sees just a part of a cluster of the dips having the same l but different values of m. The minor dips are ascribed to higher order radial modes (
= 2, 3; l may be different) and a tail portion of adjacent clusters with
= 1 (l is different by 1).
|
|

/
0, where k0 = 2
/
0 (k0a is called a size parameter). From Eq. 2, we find that k0a
/
0 is proportional to the product of the number density of BSA on the surface and the polarizability. A similar relationship exists for the TM modes. In our TE mode experiments, 103k0a
/
0 is 7.7 ± 0.1 (mean ± SD), regardless of whether the dip is deep or shallow. For the TM modes, the reduced shift is 10.3 ± 0.2. The shift's independence of the radial mode agrees with the theoretical prediction for adsorption of small particles (21
The volume of a BSA molecule, Vp, is estimated from the molecular mass (1.10 x 1019 g) and the specific volume of BSA, 0.734 g/cm3 (37
), as 81.0 nm3. To estimate the surface coverage of BSA from our experimental data, we use below a picture of a standing disk with radius R and height H. First, we note that the above Vp can be equated to a disk of R = 4.1 nm, H = 1.5 nm, where the aspect ratio is close to the one proposed (38
) in a phosphorescence study as consistent with x-ray crystallographic data (36
). Then, from Table 1, we obtain
tt/
0 as 46.1 nm3, where
r2 = 1.322 and
rp = 1.552 were used. The surface density of the BSA molecules can be estimated using the formula
![]() | (5) |
Since k0 = 2
/
0 = 4.796 µm1,
r1 = 1.4522, and k0a(
/
0)TE = 7.7 x 103 in our measurement,
is estimated as 1.3 x 104 µm3. Therefore, the area fraction
of the projection of the rectangular cross section 2RH onto the surface is estimated as
The latter value is <1/3 of the highest possible value of
by spheres, 0.55 (30
,32
34
).
We could assume another geometry for the BSA molecule, for instance, a standing rod. The molecular dimension of the rod that gives Vp = 81.0 nm3 is R = 1.9 nm and H = 7.1 nm, for example. Then,
tt/
0 = 44.9 nm3, virtually identical to the one we obtained for the disk model.
As discussed earlier,
= 0.15 is too low for the dipoles induced in nearby particles to affect the estimate of
or the TM/TE shift ratio. Therefore, we do not need to change our discussion for the surface density and orientation of the adsorbed BSA molecules.
Refractive index change of the surroundings
We tested our polarization-sensitive WGM sensor for a uniform change of relative permittivity,
in the surroundings. The change mimics adsorption of particles with a linear dimension greater than the penetration depth of the evanescent field. The shift will be greater for the mode with a greater
, since its evanescent field penetrates deeper into the surroundings. Numerical calculation of the resonance conditions (23
) gives the reduced response, k0a
/(
), of the TE mode in a microsphere with a = 174 µm at
0 = 1.312 µm as 2.507 and 2.755 for
= 1 and 2, respectively. The reduced response of the TM mode will be 2.950 and 3.252 for
= 1 and 2, respectively. The reduced response is insensitive to a: At a = 196 µm, the response for
= 1 is
0.9% less than it is at a = 174 µm.
In adding NaCl to the surroundings three times, the shifts exhibited a complicated pattern, as each radial mode had a different shift. Two parts of Fig. 5 show 103k0a
/
0 for TE and TM modes as a function of NaCl concentration c in PBS buffer surrounding a plain silica microsphere (a = 174196 µm). The data were compiled from the shifts of different dips in a few measurements. Attention was paid not to include broad dips that apparently consisted of two or more dips at any stage of NaCl addition; the shape of these dips changed as more NaCl was introduced.
|
= 1 modes. The line gives 103k0a
/(
0
c) = 0.990 L/g. With k0a
/
0 = 2.507 x 2n2
n2, we obtain dn/dc = 0.150 mL/g, slightly less than 0.171 mL/g, the value reported for
= 589 nm at 25°C (39
In both of TE and TM plots, two sets of data are away from those for
= 1. They are ascribed to the
= 2 modes. The ratio of k0a
/(
0
c) for these sets to that for
= 1 is 1.11 in TE and TM. The ratio compares favorably with the theoretical values.
| CONCLUSIONS |
|---|
|
|
|---|
45° from the TE direction, splitting the light by the polarizations right before the photodetector, and measuring the fiber transmission spectra using two photodiodes. The high sensitivity of the WGM sensor will allow such measurements at extremely low coverages of small molecules. | ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This work was supported by the National Science Foundation through BES0522668.
| FOOTNOTES |
|---|
Submitted on December 15, 2006; accepted for publication February 15, 2007.
| REFERENCES |
|---|
|
|
|---|
2. Gray, J. J. 2004. The interaction of proteins with solid surfaces. Cur. Op. Str. Bio. 14:110115.[CrossRef]
3. Robers, M., I. J. A. M. Rensink, C. E. Hack, L. A. Aarden, C. P. M. Reutelingsperger, J. F. C. Glatz, and W. Th. Hermens. 1999. A new principle for rapid immunoassay of proteins based on in situ precipitate-enhanced ellipsometry. Biophys. J. 76:27692776.
4. Larsericsdotter, H., S. Oscarsson, and J. Buijs. 2005. Structure, stability, and orientation of BSA adsorbed to silica. J. Coll. Interf. Sci. 289:2635.[CrossRef][Medline]
5. Malmsten, M. 1995. Ellipsometry studies of fibronectin adsorption. Coll. Surf. B. 3:371381.[CrossRef]
6. Sukhishvili, S. A., and S. Granick. 1999. Adsorption of human serum albumin: dependence on molecular architecture of the oppositely charged surface. J. Chem. Phys. 110:1015310161.[CrossRef]
7. Lu, J. R. 1999. Neutron reflection study of globular protein adsorption at planar interfaces. Annu. Rep. Progr. Chem. Sec. C. 95:345.
8. Giacomelli, C. E., M. G. E. G. Bremer, and W. Norde. 1999. ATR-FTIR study of IgG adsorbed on different silica surfaces. J. Coll. Interf. Sci. 220:1323.[CrossRef][Medline]
9. Giacomelli, C. E., and W. Norde. 2001. The adsorption-desorption cycle. Reversibility of the BSA-silica system. J. Coll. Interf. Sci. 233:234240.[CrossRef][Medline]
10. Noinville, S., M. Revault, M.-H. Baron, A. Tiss, S. Yapoudjian, M. Ivanova, and R. Verger. 2002. Conformational changes and orientation of Humicola lanuginosa lipase on a solid hydrophobic surface: an in situ interface Fourier transform infrared-attenuated total reflection study. Biophys. J. 82:27092719.
11. Su, T. J., J. R. Lu, R. K. Thomas, Z. F. Cui, and J. Penfold. 1998. The conformational structure of bovine serum albumin layers adsorbed at the silica-water interface. J. Phys. Chem. B. 102:81008108.
12. Vollmer, F., D. Braun, A. Libchaber, M. Khoshsima, I. Teraoka, and S. Arnold. 2002. Protein detection by optical shift of a resonant microcavity. Appl. Phys. Lett. 80:40574059.[CrossRef]
13. Armani, A. M., D. K. Armani, B. Min, K. J. Vahala, and S. M. Spillane. 2005. Ultra-high-Q microcavity operation in H2O and D2O. Appl. Phys. Lett. 87:151118.[CrossRef]
14. Serpengüzel, A., S. Arnold, and G. Griffel. 1995. Excitation of resonances of microspheres on an optical fiber. Opt. Lett. 20:654656.
15. Arnold, S., M. Khoshsima, I. Teraoka, S. Holler, and F. Vollmer. 2003. Shift of whispering-gallery modes in microspheres by protein adsorption. Opt. Lett. 28:272274.[Medline]
16. Vollmer, F., S. Arnold, D. Braun, I. Teraoka, and A. Libchaber. 2003. Multiplexed DNA detection by optical resonances in microspheres. Biophys. J. 85:19741979.
17. Noto, M., F. Vollmer, D. Keng, I. Teraoka, and S. Arnold. 2005. Nanolayer characterization through wavelength multiplexing of a microsphere resonator. Opt. Lett. 30:510512.[CrossRef][Medline]
18. Noto, M., M. Khoshsima, D. Keng, I. Teraoka, V. Kolchenko, and S. Arnold. 2005. Molecular weight dependence of a whispering gallery mode biosensor. Appl. Phys. Lett. 87:223901.[CrossRef]
19. Teraoka, I., and S. Arnold. 2006. Enhancing sensitivity of a whispering gallery mode microsphere sensor by a high-refractive index surface layer. J. Opt. Soc. Am. B. 23:14341441.[CrossRef]
20. Gaathon, O., J. Culic-Viskota, M. Mihnev, I. Teraoka, and S. Arnold. 2006. Enhancing sensitivity of a whispering gallery mode bio-sensor by sub-wavelength confinement. Appl. Phys. Lett. 89:223901.[CrossRef]
21. Teraoka, I., and S. Arnold. 2006. Theory on resonance shifts in TE and TM whispering gallery modes by non-radial perturbations for sensing applications. J. Opt. Soc. Am. B. 23:13811389.[CrossRef]
22. Johnson, B. R. 1993. Theory of morphology-dependent resonances: shape resonances and width formulas. J. Opt. Soc. Am. A. 10:343352.
23. Teraoka, I., S. Arnold, and F. Vollmer. 2003. Perturbation approach to resonance shifts of whispering galley modes in a dielectric microsphere as a probe of a surrounding medium. J. Opt. Soc. Am. B. 20:19371946.
24. Hanumegowda, N. M., C. J. Stica, B. C. Patel, I. White, and X. Fan. 2005. Refractometric sensors based on microsphere resonators. Appl. Phys. Lett. 87:201107.[CrossRef]
25. Folan, L. M. 1992. Characterization of the accretion of material by microparticles using resonant ellipsometry. Appl. Opt. 31:20662071.
26. Swann, M. J., L. L. Peel, S. Carrington, and N. J. Freeman. 2004. Dual-polarization interferometry: an analytical technique to measure changes in protein structure in real time, to determine the stoichiometry of binding events, and to differentiate between specific and nonspecific interactions. Anal. Biochem. 329:190198.[CrossRef][Medline]
27. Freeman, N. J., L. L. Peel, M. J. Swann, G. H. Cross, A. Reeves, S. Brand, and J. R. Lu. 2004. Real time, high resolution studies of protein adsorption and structure at the solid-liquid interface using dual polarization interferometry. J. Phys. Condens. Matter. 16:S2493S2496.[CrossRef]
28. Berney, H., and K. Oliver. 2005. Dual polarization interferometry size and density characterisation of DNA immobilisation and hybridisation. Biosens. Bioelectron. 21:618626.[CrossRef][Medline]
29. Ricard-Blum, S., L. L. Peel, F. Ruggiero, and N. J. Freeman. 2006. Dual polarization interferometry characterization of carbohydrate-protein interactions. Anal. Biochem. 352:252259.[CrossRef][Medline]
30. Teraoka, I., and S. Arnold. 2007. Dielectric property of particles at interface in random sequential adsorption and its application to whispering gallery mode resonance-shift sensors. J. Appl. Phys. 101:023505.[CrossRef]
31. Vörös, J. 2004. The density and refractive index of adsorbing protein layers. Biophys. J. 87:553561.
32. Onoda, G. Y., and E. G. Liniger. 1986. Experimental determination of the random-parking limit in two dimensions. Phys. Rev. A. 33:715716.[CrossRef][Medline]
33. Wang, J. S. 1994. A fast algorithm for random sequential adsorption of discs. Int. J. Mod. Phys. C. 5:707715.[CrossRef]
34. Torquato, S. 1995. Mean nearest-neighbor distance in random packings of hard d-dimensional spheres. Phys. Rev. Lett. 74:21562159.[CrossRef][Medline]
35. Lai, H. M., P. T. Leung, K. Young, P. W. Barber, and S. C. Hill. 1990. Time-independent perturbation for leaking electromagnetic modes in open systems with application to resonances in microdroplets. Phys. Rev. A. 41:51875198.[CrossRef][Medline]
36. Carter, D. C., and J. X. Ho. 1994. Structure of serum albumin. Adv. Protein Chem. 45:153196.[Medline]
37. Brill, A. S., J. S. Olin, and B. M. Siegel. 1962. The specific volume of dry protein. J. Cell Biol. 13:249251.
38. Ferrer, M. L., R. Duchowicz, B. Carrasco, J. García de la Torre, and A. U. Acuña. 2001. The conformation of serum albumin in solution: a combined phosphorescence depolarization-hydrodynamic modeling study. Biophys. J. 80:24222430.
39. Russo, P. Differential Index of Refraction, dn/dc. http://macro.lsu.edu/HowTo/.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |