| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


* Solid State Institute, Technion-Israel Institute of Technology, Haifa 32000, Israel;
Biological Structure and Function Section, Division of Biomedical Sciences, Faculty of Medicine, Imperial College, London SW7 2AZ, United Kingdom; and
Laboratory of Structural and Supramolecular Chemistry, Institute of Physical Chemistry, National Centre for Scientific Research Demokritos, Aghia Paraskevi, 15310 Athens, Greece
Correspondence: Address reprint requests to S. Stolyarova, Solid State Institute, Technion-Israel Institute of Technology, Haifa 32000, Israel. Tel.: 972-4-8293797; Fax: 972-4-8235107; E-mail: ssstolya{at}tx.technion.ac.il.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
|
|
|
, i.e., w(kL) = k
w(L) and
is not an integer,
1,2,3.... The solution of this functional equation is w(L) = AL
where A is a constant and its dimensionality is [A] = [w]/(length)
. Scaling exponents
are different for different properties. If the values of
are independent of the dilation direction i (i = x, y, or z), we have an isotropic, or self-similar, fractal. If the values of
depend on the dilation direction, the object possesses an anisotropic dilation symmetry; it is a self-affine fractal (13
The fractal properties differ for various regions of PS since it has an intermediate structure between bulk self-similar fractals like silica gels and samples with rough self-affine fractal surface and monolithic bulk, such as appear as a result of deposition or erosion processes (14
). The pore interior surface is formed by external surfaces of nanocrystallites constituting the porous layer. This surface is a self-similar fractal (6
8
,12
), i.e., it is scale-invariant along any direction up to the correlation length Lc,pore. This scale is of the order of the pore diameter (6
8
) or the size of nanocrystallites constituting the porous layer (12
). This statement can be related to any individual pore, i.e., there are many different values of Lc,pore in accordance with the pore size distribution function.
In contrast with the pore interior surface, the external surface of the PS sample is a self-affine fractal (9
,10
), i.e., it has different scale-invariant properties along different directions (14
). The external surface of the PS sample has one preferred direction perpendicular to the wafer surface. So, the self-affine fractal surface has two upper correlation lengths: the longitudinal Lcl in the (x,y) plane and the perpendicular Lcp in the z direction. Lcp characterizes the scale of the surface roughness. The surface can be described by the function h(x,y), which gives its height at position (x,y) from a referent flat surface. The surface width w(L) of a section of the surface having a size L x L in the referent plane is defined by the root mean square of the height fluctuation by w(L) = (
h2
h
2)1/2. The averaging
...
means integration over this section:
h
= L2
h(x,y)dxdy.
The morphology of the surface is defined by the definite scaling rule: if L' = kL then w(L') = k
w(L) if L' and L < Lcl. The solution of this functional equation is well known:
![]() | (1) |
Here, the roughness exponent
is limited by the condition 0 <
< 1.
The longitude correlation length Lcl appears as a result of definite correlations during electrochemical etching of monocrystalline Si wafer and depends on etching conditions. Lcl is roughly 100 nm and can be expanded by an increase in etching time (9
,10
). The surface width w(L) may be regarded as a measure of the correlation scale in the direction perpendicular to the referent surface (14
). According to the definition of the self-affine surface, this scale depends on the size of the averaging section. Its maximal value is the correlation length Lcp in the direction perpendicular to the referent plane:
Thus, according to Eq. 1,
![]() | (2) |
< 1 and the minimal scale of self-similarity l << Lcl. For PS, Lcp is of the order of the mean pore size
Lc,pore
because the surface roughness is due to the pores etched on a fragment L x L. Using the above definitions, we will show that due to the fractal properties, local high supersaturation conditions can be created near the PS surface in the metastable solution. Supersaturation is a necessary, but not sufficient, condition for spontaneous nucleation: a free energy barrier of entropic origin must be overcome for a critical nucleus to form. An obvious way to promote nucleation is by inducing a local supersaturation spike, thus lowering the barrier at ideally one place in the medium, which will be the point where the critical nucleus forms. Naturally, the spatial size of this local spike should not be less than the size of the critical nucleus. A substrate that would locally increase the concentration of the molecules to be crystallized, with respect to the bulk, may therefore be an effective promoter of heterogeneous nucleation. The substrate does not need to be a postcritical nucleus (i.e., a seed crystal) itself, nor to be able to simultaneously constrain the minimal number of molecules that are necessary for criticality. A limited constraint, leading to a small increase in the probability of critical nucleus formation at the substrate, may be sufficient to induce crystallization if the system is not too far from labile conditions.
A high local supersaturation is possible near any surface if attractive forces are large enough. These conditions are realized in the cases of capillary condensation of neutral particles (15
,16
) and adsorption of charged particles on an oppositely charged surface (17
). The attractive force between opposite charges of a molecule and a surface can produce a sufficiently large difference in surface free energy to produce sufficient local supersaturation of the charged particles, resulting in heterogeneous nucleation (17
). Charged surfaces, including functionalized self-assembled monolayers (18
), doped silicon microfluidic devices (19
), poly-L-lysin covered glass surfaces (20
), and charged beads like Sephadex and various zeolites, which combine charge with an ordinary regular porous structure (21
), have all been tested as nucleation promoters. Most of these have been particularly successful with lysozyme, an easily crystallizable, fairly compact protein with a high surface charge density when at crystallization conditions. Success with a wider range of proteins has not been forthcoming and therefore such surfaces and beads have not gained widespread use.
Hydrophobic surfaces have also been shown to promote nucleation in some cases, possibly due to preferential binding of the few hydrophobic residues on the protein surface with the nucleant surface (22
). This, however, is not true in general: hydrophobic surfaces are commonly used as inert substrates for protein crystallization by homogeneous nucleation.
A nucleation promoter can also be one that provides a template for epitaxial crystal growth. This is the obvious mechanism of one of the most common methods of improving macromolecular crystals, namely, seeding metastable solutions with a previously obtained small crystal of the sample (homologous seeding). Heteroepitaxial growth on minerals has also been attempted for protein crystallization (23
). The presence of a heteroepitaxial growth mechanism has been demonstrated, but no single mineral has been found that can be used for a variety of different proteins. Hence, this arduous route has been largely abandoned by the macromolecular crystallographic community, for the moment at least.
In sum, it has been argued that a series of interactions may drive sufficient local supersaturation to lead to heterogeneous nucleation. These range from attraction between net opposite charges to short-range specific and nonspecific interactions (24
), dispersion forces (17
), hydrophobicity, etc.
For the most obvious case of attraction between net opposite charges, the signs of net molecular charge and surface charge are crucial for crystallization (17
,19
). However, the results of Chayen et al. (2
) are independent of charge sign since the proteins under investigation had different net charge signs, yet crystallized on the same PS samples. The precise origin of the attraction potential is thus likely to vary between proteins, and it is not the purpose of this work to discuss it further.
We consider here the case where there is a definite attraction potential between proteins and silicon surface resulting in protein accumulation on the surface layer at a concentration higher than in the solution bulk (positive adsorption (25
)) but still not sufficient for heterogeneous nucleation on a flat surface. We have focused on the fractal properties of PS surfaces, whichprovided that some kind of sufficiently attractive potential is presentmake it a more effective and wide-ranging nucleant than (functionalized) flat surfaces or surfaces with an ordinary regular porous structure.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
![]() | (3) |
L here refers to the shortest distance between two adjacent corners of the fragment, not to the length of the cross sectional fractal curve; D > 2 is the fractal dimension of the surface. In the case of self-affine surface, D = 3
. The correlation length Lc = Lc,pore and l < R < min(L, Lc,pore) for a self-similar pore surface. In the case of a self-affine surface, Lc = Lcl and l < R < w(L) < min(L, Lcp, Lcl).
Equation 3 accounts for the well-known effect that NR(L)/L2 for particles deposited on fractal surfaces with the same fractal dimension D has different values due to different upper limits Lc of the fractal regime (26
)the greater the correlation scale Lc, the higher that quantity:
![]() | (4) |
Note that the quantity NR(L)/L2 is an intensive thermodynamic variable (i.e., one independent of the size of averaging region L) for representative volume only (L >> Lc). We call this quantity projected effective density, nef. It corresponds to the density of the projection onto a flat surface (of area L x L) of a set of particles lying on the corrugated surface. As is clear from Eq. 4, nef is much higher than that on a regular surface: since D > 2 and Lc >> R, it results that nef >> 1/R2, where R2 would be the surface density on a flat surface for a complete monolayer.
We are exploiting this effect for the case of molecular ion adsorption from solution. We certainly do not have a complete monolayer due to mutual repulsion of the ions, and the mean distance between adsorbed ions can be greater than their size. Nevertheless, an effect analogous to that described by Eq. 3 is possible. Let the density of adsorbed molecules on a flat surface be ns and the mean distance
between absorbed particles be less than Lc: (
). The value of ns depends on the concentration of particles in the metastable solution, interaction between the particles, and interaction of the particles with the surface. Certainly, we suppose that positive adsorption takes place (25
), i.e., the mean distance
between absorbed particles is less than the mean distance between solute particles in the solution volume. So,
is a natural spatial scale for our model, and all distances will be measured in these units.
We can then speak of monolayer covering of the surface by "particles" of size
. Thus, the effective projected density nef of these particles lying on the fractal surface is much higher than their density ns on a regular surface:
![]() | (5) |
. The effect is absent if the molecule concentration in solution is not high enough and the mean distance
between absorbed particles on the flat surface is much more than Lc (
). Since the correlation scale Lc is defined by the PS preparation conditions, it is possible to select the value of Lc for any concentration of metastable solution for the local density nef to be above metastability. From simple geometric considerations, a higher nef implies a shorter average interparticle distance in physical (3D) space for the same average interparticle distance along the surface (i.e., on a two-dimensional metric attached to the surface). This can be intuitively pictured as follows: imagine a corrugated surface with particles rigidly fixed to it at given distances. This surface is now stretched out until completely flat, but no more. If the particles have remained fixed with respect to the surface, they will find themselves at greater distances from each other in 3D space than when the surface was corrugated. So, even if the actual surface density of the molecules on the fractal surface remains the same as on a flat surface, a higher nef means that the molecules are still brought much closer together than they would be if they were lying on the flat surface. Thus, in the vicinity of the fractal surface fragment, the effective surface density of protein molecules will be greater than in the vicinity of the flat surface.
Since we are interested in heterogeneous nucleation, let us show that the real local volume concentration nv,loc in the vicinity of the fractal surface is greater than in the vicinity of a flat surface. Moreover, this concentration should be no less than a critical concentration for spontaneous nucleation at given pressure P and temperature T, ncr(P, T). Naturally, the thickness z of this near-surface layer is no less than the size of the critical nucleus rcr, z
rcr at the same time or z
Lcp because we would like to use scaling relation Eq. 3. Here Lcp is the correlation length in the direction perpendicular to the referent plane of the self-affine surface if we are considering supersaturation in the vicinity of this surface and Lcp = Lc,pore if we are considering supersaturation inside the pore.
So, the relation of scales for the layer is
![]() | (6) |
Let us estimate the volume of this boundary layer Vb(z) and the number N(z) of particles that can be placed inside this layer.
According to the fractal surface definition (13
) and Eq. 3 we can calculate Vb(z) as a complete volume of (L/z)D boxes, each of volume equal to z3:
![]() | (7) |
We would now like to arrange the particles inside this volume. Since z is the distance in physical (3D) space, there will be
layers of particles. The number of particles in each layer is defined by the surface fractality and is equal to
. So, the total number of particles is
![]() | (8) |
The local volume density will therefore be
![]() | (9) |
Here we suppose that the distance between particles in the direction normal to the flat surface coincides with distance
along the surface, i.e., the volume density near the flat surface
. In a more realistic case, the distance between particles in the direction normal to the surface may be more than
since adsorption force (for example, van der Waals dispersion force) decreases with z. However, this is not important for our estimates since we use the same conditions for the flat surface, too.
So, although the boundary layer volume increases as z3D (3 D < 1 since D > 2), the number of particles inside the layer increases as z1 (Eqs. 7 and 8). This is the reason for the increase of local volume density. Recall that this is the effect of being in the vicinity of a fractal surface, since it is necessary that the condition in Eq. 6 be fulfilled. These conditions are rather rigid since in reality, fractal behavior is typically based on a scaling range that spans 0.52 decades (27
).
Strictly speaking, LD should be replaced by the Hausdorff D-measure mD(
) for the fractal surface (13
). Then the boundary layer volume Vb(z) is expressed as
![]() | (10) |
Here prefactor B is proportional to the Hausdorff D-measure mD(
) for a fractal (self-affine or self-similar) surface
(28
):
![]() | (11) |
The number of particles N(z) inside the volume Vb(z) is equal to the product of the number of layers with thickness
and the number of particles inside one layer,
:
![]() | (12) |
So, the local particle density nv(z) = N(z)/Vb(z) inside a boundary layer of thickness z near fractal surface with fractal dimension D depends on D and z only and is independent of the Hausdorff D-measure
![]() | (13) |
We can see that the results of Eqs. 9 and 13 are identical. It is more convenient to operate with a fragment of the surface with lateral size L. It is known that the Hausdorff D-measure mD(
)
LD. Since the effective density nV(z) is independent of mD(
) we prefer the simplest form for N(z) and Vb(z):
.
Thus, any fractal substrate yields a sufficient local concentration for nucleation if the following conditions are fulfilled:
![]() | (14) |
Here ncr (P, T) is a critical concentration for spontaneous nucleation at given pressure P and temperature T. The conditions in Eq. 14 let us make an optimal choice of crystallization and substrate-related parameters for heterogeneous nucleation of any protein. Indeed, ns is a function of solution metastability and silicon-protein interaction, whereas fractal dimension D and correlation length Lc are controlled by fabrication conditions. It was shown that the fractal dimension D, reflecting surface morphology of individual pore and the whole structure of porous network, is strongly porosity dependent (12
). For example, its value can be changed from 2.1 to 2.4 for porosity range 4070%, indicating an increase of nanocrystallite roughness and, consequently, the specific surface area with porosity (12
). In turn, the correlation length Lc was shown to be a function of anodization time (9
,10
). The pore size which serves as a correlation length for an individual fractal pore can be varied from nanometer to micrometer scale. Practically, the pore size as well as the porosity can be changed by varying the anodization current density, HF solution concentration as well as doping type and level of the starting silicon substrate (29
,30
). So, we can select favorable fabrication parameters to tailor fractal dimension D and correlation length Lc for a given protein to satisfy nucleation conditions. Assuming that useful pores for our purposes range in size from the diameter of a protein molecule to that of a critical nucleus, we would like to have a pore size distribution that ranges over a few protein diameters.
It remains however to be judged whether a very wide, "nontailored" and therefore more "universal" pore size distribution is to be preferred over a more protein-adapted tailored surface. A broad distribution of fractal pore sizes in PS sample may indeed be a crucial advantage, in providing for each of various proteins the pore size most suited to its molecular diameter and to the shape of the initial aggregates that it forms. On the same PS surface, different pores provide different supersaturation conditions for a given protein molecule because of the variation of the fractal parameters. Hence, the same PS can serve as a nucleant for various protein molecules. To reach conclusions on this matter, more crystallization experience with this and possibly other similar materials will doubtless be required.
| CONCLUSIONS |
|---|
|
|
|---|
Local supersaturation "spikes" are not in themselves sufficient to guarantee successful crystal growth, especially where macromolecular crystals are concerned. For crystal nucleation, rather than the triggering of phase separation or precipitation/amorphous aggregation, the protein solution has to be metastable with respect to a crystalline phase. In other words, the protein solution must be close to conditions where the protein would crystallize spontaneously in the absence of a nucleant. Practically, this means that the specific interactions between protein molecules (the ones that lead to 3D crystal formation) must be of comparable magnitude to the protein-nucleant attraction. If protein-protein interactions are too weak at the given conditions, then deposition of molecules on the nucleant surface will at best be achieved. If the interactions are not specific and conducive to ordered packing, then only amorphous aggregation can result.
These caveats should not detract from the importance of heterogeneous nucleation, which has two advantages. First, nucleation and growth at metastable conditions is highly advantageous because the effects of spontaneous nucleation are prevented. The slower growth and the lack of excess and secondary nucleation often provide for growth of larger, better diffracting crystals (17
). Second, the probability of a hit in the usual trial-and-error initial crystallization screening of a protein is increased, since metastable conditions will also yield crystals in the presence of the heterogeneous nucleant whereas they would remain clear in the absence of nucleant.
Our model shows that a sufficient local concentration of molecules for nucleation is possible inside and in the close vicinity of the pores, even when the average conditions in the bulk of the solution correspond to metastability. Whatever the mechanism of molecule adsorption, with correct choice of crystallization and of surface-related parameters, heterogeneous nucleation may occur due to the local supersaturation at the fractal surface. The PS substrate can serve as a rather universal nucleant for various protein molecules due to the wide distribution of fractal pore sizes. In addition, the PS technology is very flexible, allowing tailoring the pore size and concentration as well as the fractal properties to specific proteins by changing the fabrication conditions. It should be noted that heterogeneous nucleation of protein crystals never succeeded on oxidized PS, which had lost its fractal properties (11
). This fact is probably an additional demonstration of the role of fractality in the promotion of heterogeneous nucleation.
According to our theory, this nucleation-promoting property of PS is more general and is not limited to amphoteric macromolecules only. Indeed, it was also observed in the case of enhanced crystallization and bonding of ceramics and oxides to PS surfaces (4
,5
). On the other hand, our consideration of the phenomenon of local supersaturation caused by fractality does not rely upon specific parameters of PS and may be applied to any fractal surface.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Submitted on February 2, 2006; accepted for publication July 6, 2006.
| REFERENCES |
|---|
|
|
|---|
2. Chayen, N. E., E. Saridakis, A. El-Bahar, and Y. Nemirovsky. 2001. Porous silicon: a nucleation inducing material for protein crystallization. J. Mol. Biol. 312:591596.[CrossRef][Medline]
3. Canham, L. T. 1997. Biomedical applications of porous silicon. In Properties of Porous Silicon. L. Canham, editor. EMIS Datareviews series No. 18, INSPEC. Institution of Engineering and Technology, London, UK. 371376.
4. Stolyarova, S., A. El-Bahar, and Y. Nemirovsky. 2002. Unexpected room temperature growth of silicon dioxide crystallites on passivated porous silicon. J. Cryst. Growth. 237239:19201925.[CrossRef]
5. Stolyarova, S., B. Malic, S. Javoric, A. El-Bahar, M. Kosec, and Y. Nemirovsky. 2003. Integration of porous silicon with sol-gel derived ceramic films. Proc. Mater. Res. Soc. Sym. USA. 768:5762.
6. Di Francia, G., V. La Ferrara, P. Maddalena, D. Ninno, L. P. Odierna, and V. Cataudella. 1996. AC conductivity of porous silicon: a fractal and surface transport mechanism. Nuovo-Cimento. D18:11871196.
7. Matthai, C. C., J. L. Gavartin, and A. A. Cafolla. 1995. Structural and elastic properties of porous silicon. Thin Solid Films. 255:174176.[CrossRef]
8. Aleksandrow, L. N., and P. L. Novikov. 1998. Mechanisms of formation and topological analysis of porous silicon - computational modeling. Comput. Mater. Sci. 10:406410.[CrossRef]
9. Happo, N., M. Imamatsu, and K. Horii. 2000. Self-affine fractal of porous silicon surface before and after natural oxidation. Phys. Status Solidi A. 182:233237.[CrossRef]
10. Happo, N., M. Fujiwara, M. Iwamatsu, and K. Horii. 1998. Atomic force microscopy study of self-affine fractal roughness of porous silicon surfaces. Jpn. J. Appl. Phys. 37:39513953.[CrossRef]
11. Young, T. F., I. W. Huang, Y. L. Yang, W. C. Kuo, I. M. Jing, T. C. Chang, and C. Y. Chang. 1996. Atomic force microscopy study of the surface structure of oxidized porous silicon. Appl. Surf. Sci. 102:404407.[CrossRef]
12. Nychyporuk, T., V. Lysenko, and D. Barbier. 2005. Fractal structure of porous silicon nanocrystallites. Phys. Rev. B71:11540211154025.
13. Gouyet, J.-F. 1996. Physics and Fractal Structures. Springer-Verlag, Berlin/Heidelberg, Germany.
14. Family, F. 1990. Dynamic scaling and phase transitions in interface growth. Physica A. 168:561580.[CrossRef]
15. Donley, J. P., and A. J. Liu. 1997. Phase behavior of near-critical fluids confined in periodic gels. Phys. Rev. E. 55:539543.[CrossRef]
16. Kunzner, N., E. Gross, J. Diener, D. Kovalev, V. Y. Timoshenko, and D. Wallacher. 2003. Capillary condensation monitored in birefringent porous silicon layers. J. Appl. Phys. 94:49134917.[CrossRef]
17. Sear, R. P. 2003. Protein crystals and charged surfaces: interactions and heterogeneous nucleation. Phys. Rev E. 67:06190710619107.
18. Zhu, P., Y. Masuda, and K. Koumoto. 2004. The effect of surface charge on hydroxyapatite nucleation. Biomaterials. 25:39153921.[CrossRef][Medline]
19. Sanjoh, A., and T. Tsukihara. 1999. Spatiotemporal protein crystal growth studies using microfluidic silicon devices. J. Cryst. Growth. 196:691702.[CrossRef]
20. Nanev, C. N., and D. Tsekova. 2000. Heterogeneous nucleation of hen-egg-white lysozymemolecular approach. Cryst. Res. Technol. 35:189195.[CrossRef]
21. Chayen, N. E., and E. Saridakis. 2001. Is lysozyme really the ideal model protein? J. Cryst. Growth. 232:262264.[CrossRef]
22. Tsekova, D., S. Dimitrova, and C. N. Nanev. 1999. Heterogeneous nucleation (and adhesion) of lysozyme crystals. J. Cryst. Growth. 196:229233.[CrossRef]
23. McPherson, A., and P. Shlichta. 1988. Heterogeneous and epitaxial nucleation of protein crystals on mineral surfaces. Science. 239:385387.
24. Falini, G., S. Fermani, G. Conforti, and A. Ripamonti. 2002. Protein crystallization on chemically modified mica surfaces. Acta Crystallogr. D58:16491652.
25. Landau, L. D., and E. M. Lifshitz. 1993. Statistical Physics. Pergamon Press, Oxford/New York/Seoul/Tokyo.
26. Pfeifer, P., Y. J. Wu, M. W. Cole, and J. Krim. 1989. Multilayer adsorption on a fractally rough surface. Phys. Rev. Lett. 62:19972000.[CrossRef][Medline]
27. Malcai, O., D. A. Lidar, O. Biham, and D. Avnir. 1997. Scaling range and cutoffs in empirical fractals. Phys. Rev. E. 56:28172827.[CrossRef]
28. Wong, Pz., and A. J. Bray. 1988. Porod scattering from fractal surfaces. Phys. Rev. Lett. 60:1344.[CrossRef][Medline]
29. Halimaoui, J. 1997. Porous silicon formation by anodisation. In Properties of Porous Silicon. L. Canham, editor. EMIS Datareviews series No. 18, INSPEC. Institution of Engineering and Technology, London, UK. 1222.
30. Herino, R. 1997. Pore size distribution in porous silicon. In Properties of Porous Silicon. L. Canham, editor. EMIS Datareviews series No. 18, INSPEC. Institution of Engineering and Technology, London, UK. 8996.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |