Using three-dimensional microfluidic networks for solving computationally hard problems

This paper describes the design of a parallel algorithm that uses moving fluids in a three-dimensional microfluidic system to solve a nondeterministically polynomial complete problem (the maximal clique problem) in polynomial time. This algorithm relies on (i) parallel fabrication of the microfluidic system, (ii) parallel searching of all potential solutions by using fluid flow, and (iii) parallel optical readout of all solutions. This algorithm was implemented to solve the maximal clique problem for a simple graph with six vertices. The successful implementation of this algorithm to compute solutions for small-size graphs with fluids in microchannels is not useful, per se, but does suggest broader application for microfluidics in computation and control.

Quantifying The Diffusion And Erosion Of Degradable Hydrogels Using A Microfluidic Technique

Hydrogels are composed of cross-linked networks of hydrophilic polymers that are biocompatible due to their high water content. Mass transfer through hydrogels has been suggested as an effective method of drug delivery, specifically in degradable polymers to minimize lasting effects within the body. Diffusion of small molecules in poly (ethylene glycol) diacrylate (PEG-DA) and dextran methacrylate (dex-MA) hydrogels was characterized in a microfluidic device and by complementary techniques. Microfluidic devices were prepared by crosslinking a formulation of hydrogel and photo-initiator, with and without visible dye, using photolithography to define a central microchannel. Channel sizes within the devices were approximately 600 ¿m to simulate vessels within the body. The microfluidic technique allows for both image and effluent analyses. To visualize the diffusive behavior within the dextran hydrogel, methylene blue and sulforhodamine 101 dyes were used in both elution and uptake experiments. Three analysis techniques for measuring diffusion coefficients were used to quantify the diffusion of solute in the hydrogel, including optical microscopy, characterization of device effluent, and NMR analyses. The optical microscopy technique analyzes images of the dye diffusion captured by a stereomicroscope to generate dye concentration v. position profiles. The data was fit to a diffusion model to determine diffusion coefficients and the dye release profile. In a typical elution experiment, aqueous solution is pumped through the microchannel and dye diffuses out of the hydrogel and into the aqueous phase. During elution, images are taken at regular time intervals and the effluent was collected. Analysis of the device effluent was performed using ultraviolet-visible (UV/Vis) spectroscopy to determine the effluent dye concentration and thus a short-time diffusion coefficient. Nuclear magnetic resonance (NMR) was used to determine a free diffusion coefficient of molecules in hydrogel without the effect of a concentration gradient. Diffusion coefficients for methylene blue and sulforhodamine 101 dyes in dex-MA hydrogel calculated using the three analysis methods all agree well. It was determined that utilizing a combination of the three techniques offers greater insight into molecular diffusion in hydrogels than employing each technique individually. The use of the same microfluidic devices used to measure diffusion is explored in the use of studying the degradation of dex-MA hydrogels. By combining what is known about the degradation rate in regards to the effect of pH and crosslinking and the ability to use a dye solution in contrast to establish the hydrogel boundaries could be a novel approach to studying hydrogel degradation.

Suitability of artificial neural networks for designing LoC circuits

he simulation of complex LoC (Lab-on-a-Chip) devices is a process that requires solving computationally expensive partial differential equations. An interesting alternative uses artificial neural networks for creating computationally feasible models based on MOR techniques. This paper proposes an approach that uses artificial neural networks for designing LoC components considering the artificial neural network topology as an isomorphism of the LoC device topology. The parameters of the trained neural networks are based on equations for modeling microfluidic circuits, analogous to electronic circuits. The neural networks have been trained to behave like AND, OR, Inverter gates. The parameters of the trained neural networks represent the features of LoC devices that behave as the aforementioned gates. This would mean that LoC devices universally compute.

Microfluidic platforms for the characterization of in meso membrane protein crystallization

Membrane proteins, which reside in the membranes of cells, play a critical role in many important biological processes including cellular signaling, immune response, and material and energy transduction. Because of their key role in maintaining the environment within cells and facilitating intercellular interactions, understanding the function of these proteins is of tremendous medical and biochemical significance. Indeed, the malfunction of membrane proteins has been linked to numerous diseases including diabetes, cirrhosis of the liver, cystic fibrosis, cancer, Alzheimer's disease, hypertension, epilepsy, cataracts, tubulopathy, leukodystrophy, Leigh syndrome, anemia, sensorineural deafness, and hypertrophic cardiomyopathy.1-3 However, the structure of many of these proteins and the changes in their structure that lead to disease-related malfunctions are not well understood. Additionally, at least 60% of the pharmaceuticals currently available are thought to target membrane proteins, despite the fact that their exact mode of operation is not known.4-6 Developing a detailed understanding of the function of a protein is achieved by coupling biochemical experiments with knowledge of the structure of the protein. Currently the most common method for obtaining three-dimensional structure information is X-ray crystallography. However, no a priori methods are currently available to predict crystallization conditions for a given protein.7-14 This limitation is currently overcome by screening a large number of possible combinations of precipitants, buffer, salt, and pH conditions to identify conditions that are conducive to crystal nucleation and growth.7,9,11,15-24 Unfortunately, these screening efforts are often limited by difficulties associated with quantity and purity of available protein samples. While the two most significant bottlenecks for protein structure determination in general are the (i) obtaining sufficient quantities of high quality protein samples and (ii) growing high quality protein crystals that are suitable for X-ray structure determination,7,20,21,23,25-47 membrane proteins present additional challenges. For crystallization it is necessary to extract the membrane proteins from the cellular membrane. However, this process often leads to denaturation. In fact, membrane proteins have proven to be so difficult to crystallize that of the more than 66,000 structures deposited in the Protein Data Bank,48 less than 1% are for membrane proteins, with even fewer present at high resolution (< 2Å)4,6,49 and only a handful are human membrane proteins.49 A variety of strategies including detergent solubilization50-53 and the use of artificial membrane-like environments have been developed to circumvent this challenge.43,53-55 In recent years, the use of a lipidic mesophase as a medium for crystallizing membrane proteins has been demonstrated to increase success for a wide range of membrane proteins, including human receptor proteins.54,56-62 This in meso method for membrane protein crystallization, however, is still by no means routine due to challenges related to sample preparation at sub-microliter volumes and to crystal harvesting and X-ray data collection. This dissertation presents various aspects of the development of a microfluidic platform to enable high throughput in meso membrane protein crystallization at a level beyond the capabilities of current technologies. Microfluidic platforms for protein crystallization and other lab-on-a-chip applications have been well demonstrated.9,63-66 These integrated chips provide fine control over transport phenomena and the ability to perform high throughput analyses via highly integrated fluid networks. However, the development of microfluidic platforms for in meso protein crystallization required the development of strategies to cope with extremely viscous and non-Newtonian fluids. A theoretical treatment of highly viscous fluids in microfluidic devices is presented in Chapter 3, followed by the application of these strategies for the development of a microfluidic mixer capable of preparing a mesophase sample for in meso crystallization at a scale of less than 20 nL in Chapter 4. This approach was validated with the successful on chip in meso crystallization of the membrane protein bacteriorhodopsin. In summary, this is the first report of a microfluidic platform capable of performing in meso crystallization on-chip, representing a 1000x reduction in the scale at which mesophase trials can be prepared. Once protein crystals have formed, they are typically harvested from the droplet they were grown in and mounted for crystallographic analysis. Despite the high throughput automation present in nearly all other aspects of protein structure determination, the harvesting and mounting of crystals is still largely a manual process. Furthermore, during mounting the fragile protein crystals can potentially be damaged, both from physical and environmental shock. To circumvent these challenges an X-ray transparent microfluidic device architecture was developed to couple the benefits of scale, integration, and precise fluid control with the ability to perform in situ X-ray analysis (Chapter 5). This approach was validated successfully by crystallization and subsequent on-chip analysis of the soluble proteins lysozyme, thaumatin, and ribonuclease A and will be extended to microfluidic platforms for in meso membrane protein crystallization. The ability to perform in situ X-ray analysis was shown to provide extremely high quality diffraction data, in part as a result of not being affected by damage due to physical handling of the crystals. As part of the work described in this thesis, a variety of data collection strategies for in situ data analysis were also tested, including merging of small slices of data from a large number of crystals grown on a single chip, to allow for diffraction analysis at biologically relevant temperatures. While such strategies have been applied previously,57,59,61,67 they are potentially challenging when applied via traditional methods due to the need to grow and then mount a large number of crystals with minimal crystal-to-crystal variability. The integrated nature of microfluidic platforms easily enables the generation of a large number of reproducible crystallization trials. This, coupled with in situ analysis capabilities has the potential of being able to acquire high resolution structural data of proteins at biologically relevant conditions for which only small crystals, or crystals which are adversely affected by standard cryocooling techniques, could be obtained (Chapters 5 and 6). While the main focus of protein crystallography is to obtain three-dimensional protein structures, the results of typical experiments provide only a static picture of the protein. The use of polychromatic or Laue X-ray diffraction methods enables the collection of time resolved structural information. These experiments are very sensitive to crystal quality, however, and often suffer from severe radiation damage due to the intense polychromatic X-ray beams. Here, as before, the ability to perform in situ X-ray analysis on many small protein crystals within a microfluidic crystallization platform has the potential to overcome these challenges. An automated method for collecting a "single-shot" of data from a large number of crystals was developed in collaboration with the BioCARS team at the Advanced Photon Source at Argonne National Laboratory (Chapter 6). The work described in this thesis shows that, even more so than for traditional structure determination efforts, the ability to grow and analyze a large number of high quality crystals is critical to enable time resolved structural studies of novel proteins. In addition to enabling X-ray crystallography experiments, the development of X-ray transparent microfluidic platforms also has tremendous potential to answer other scientific questions, such as unraveling the mechanism of in meso crystallization. For instance, the lipidic mesophases utilized during in meso membrane protein crystallization can be characterized by small angle X-ray diffraction analysis. Coupling in situ analysis with microfluidic platforms capable of preparing these difficult mesophase samples at very small volumes has tremendous potential to enable the high throughput analysis of these systems on a scale that is not reasonably achievable using conventional sample preparation strategies (Chapter 7). In collaboration with the LS-CAT team at the Advanced Photon Source, an experimental station for small angle X-ray analysis coupled with the high quality visualization capabilities needed to target specific microfluidic samples on a highly integrated chip is under development. Characterizing the phase behavior of these mesophase systems and the effects of various additives present in crystallization trials is key for developing an understanding of how in meso crystallization occurs. A long term goal of these studies is to enable the rational design of in meso crystallization experiments so as to avoid or limit the need for high throughput screening efforts. In summary, this thesis describes the development of microfluidic platforms for protein crystallization with in situ analysis capabilities. Coupling the ability to perform in situ analysis with the small scale, fine control, and the high throughput nature of microfluidic platforms has tremendous potential to enable a new generation of crystallographic studies and facilitate the structure determination of important biological targets. The development of platforms for in meso membrane protein crystallization is particularly significant because they enable the preparation of highly viscous mixtures at a previously unachievable scale. Work in these areas is ongoing and has tremendous potential to improve not only current the methods of protein crystallization and crystallography, but also to enhance our knowledge of the structure and function of proteins which could have a significant scientific and medical impact on society as a whole. The microfluidic technology described in this thesis has the potential to significantly advance our understanding of the structure and function of membrane proteins, thereby aiding the elucidation of human biology, the development of pharmaceuticals with fewer side effects for a wide range of diseases. References (1) Quick, M.; Javitch, J. A. P Natl Acad Sci USA 2007, 104, 3603. (2) Trubetskoy, V. S.; Burke, T. J. Am Lab 2005, 37, 19. (3) Pecina, P.; Houstkova, H.; Hansikova, H.; Zeman, J.; Houstek, J. Physiol Res 2004, 53, S213. (4) Arinaminpathy, Y.; Khurana, E.; Engelman, D. M.; Gerstein, M. B. Drug Discovery Today 2009, 14, 1130. (5) Overington, J. P.; Al-Lazikani, B.; Hopkins, A. L. Nat Rev Drug Discov 2006, 5, 993. (6) Dauter, Z.; Lamzin, V. S.; Wilson, K. S. Current Opinion in Structural Biology 1997, 7, 681. (7) Hansen, C.; Quake, S. R. Current Opinion in Structural Biology 2003, 13, 538. (8) Govada, L.; Carpenter, L.; da Fonseca, P. C. A.; Helliwell, J. R.; Rizkallah, P.; Flashman, E.; Chayen, N. E.; Redwood, C.; Squire, J. M. J Mol Biol 2008, 378, 387. (9) Hansen, C. L.; Skordalakes, E.; Berger, J. M.; Quake, S. R. P Natl Acad Sci USA 2002, 99, 16531. (10) Leng, J.; Salmon, J.-B. Lab Chip 2009, 9, 24. (11) Zheng, B.; Gerdts, C. J.; Ismagilov, R. F. Current Opinion in Structural Biology 2005, 15, 548. (12) Lorber, B.; Delucas, L. J.; Bishop, J. B. J Cryst Growth 1991, 110, 103. (13) Talreja, S.; Perry, S. L.; Guha, S.; Bhamidi, V.; Zukoski, C. F.; Kenis, P. J. A. The Journal of Physical Chemistry B 2010, 114, 4432. (14) Chayen, N. E. Current Opinion in Structural Biology 2004, 14, 577. (15) He, G. W.; Bhamidi, V.; Tan, R. B. H.; Kenis, P. J. A.; Zukoski, C. F. Cryst Growth Des 2006, 6, 1175. (16) Zheng, B.; Tice, J. D.; Roach, L. S.; Ismagilov, R. F. Angew Chem Int Edit 2004, 43, 2508. (17) Li, L.; Mustafi, D.; Fu, Q.; Tereshko, V.; Chen, D. L. L.; Tice, J. D.; Ismagilov, R. F. P Natl Acad Sci USA 2006, 103, 19243. (18) Song, H.; Chen, D. L.; Ismagilov, R. F. Angew Chem Int Edit 2006, 45, 7336. (19) van der Woerd, M.; Ferree, D.; Pusey, M. Journal of Structural Biology 2003, 142, 180. (20) Ng, J. D.; Gavira, J. A.; Garcia-Ruiz, J. M. Journal of Structural Biology 2003, 142, 218. (21) Talreja, S.; Kenis, P. J. A.; Zukoski, C. F. Langmuir 2007, 23, 4516. (22) Hansen, C. L.; Quake, S. R.; Berger, J. M. US, 2007. (23) Newman, J.; Fazio, V. J.; Lawson, B.; Peat, T. S. Cryst Growth Des 2010, 10, 2785. (24) Newman, J.; Xu, J.; Willis, M. C. Acta Crystallographica Section D 2007, 63, 826. (25) Collingsworth, P. D.; Bray, T. L.; Christopher, G. K. J Cryst Growth 2000, 219, 283. (26) Durbin, S. D.; Feher, G. Annu Rev Phys Chem 1996, 47, 171. (27) Talreja, S.; Kim, D. Y.; Mirarefi, A. Y.; Zukoski, C. F.; Kenis, P. J. A. J Appl Crystallogr 2005, 38, 988. (28) Yoshizaki, I.; Nakamura, H.; Sato, T.; Igarashi, N.; Komatsu, H.; Yoda, S. J Cryst Growth 2002, 237, 295. (29) Anderson, M. J.; Hansen, C. L.; Quake, S. R. P Natl Acad Sci USA 2006, 103, 16746. (30) Hansen, C. L.; Sommer, M. O. A.; Quake, S. R. P Natl Acad Sci USA 2004, 101, 14431. (31) Lounaci, M.; Rigolet, P.; Abraham, C.; Le Berre, M.; Chen, Y. Microelectron Eng 2007, 84, 1758. (32) Zheng, B.; Roach, L. S.; Ismagilov, R. F. J Am Chem Soc 2003, 125, 11170. (33) Zhou, X.; Lau, L.; Lam, W. W. L.; Au, S. W. N.; Zheng, B. Anal. Chem. 2007. (34) Cherezov, V.; Caffrey, M. J Appl Crystallogr 2003, 36, 1372. (35) Qutub, Y.; Reviakine, I.; Maxwell, C.; Navarro, J.; Landau, E. M.; Vekilov, P. G. J Mol Biol 2004, 343, 1243. (36) Rummel, G.; Hardmeyer, A.; Widmer, C.; Chiu, M. L.; Nollert, P.; Locher, K. P.; Pedruzzi, I.; Landau, E. M.; Rosenbusch, J. P. Journal of Structural Biology 1998, 121, 82. (37) Gavira, J. A.; Toh, D.; Lopez-Jaramillo, J.; Garcia-Ruiz, J. M.; Ng, J. D. Acta Crystallogr D 2002, 58, 1147. (38) Stevens, R. C. Current Opinion in Structural Biology 2000, 10, 558. (39) Baker, M. Nat Methods 2010, 7, 429. (40) McPherson, A. In Current Topics in Membranes, Volume 63; Volume 63 ed.; DeLucas, L., Ed.; Academic Press: 2009, p 5. (41) Gabrielsen, M.; Gardiner, A. T.; Fromme, P.; Cogdell, R. J. In Current Topics in Membranes, Volume 63; Volume 63 ed.; DeLucas, L., Ed.; Academic Press: 2009, p 127. (42) Page, R. In Methods in Molecular Biology: Structural Proteomics - High Throughput Methods; Kobe, B., Guss, M., Huber, T., Eds.; Humana Press: Totowa, NJ, 2008; Vol. 426, p 345. (43) Caffrey, M. Ann Rev Biophys 2009, 38, 29. (44) Doerr, A. Nat Methods 2006, 3, 244. (45) Brostromer, E.; Nan, J.; Li, L.-F.; Su, X.-D. Biochemical and Biophysical Research Communications 2009, 386, 634. (46) Li, G.; Chen, Q.; Li, J.; Hu, X.; Zhao, J. Anal Chem 2010, 82, 4362. (47) Jia, Y.; Liu, X.-Y. The Journal of Physical Chemistry B 2006, 110, 6949. (48) RCSB Protein Data Bank. http://www.rcsb.org/ (July 11, 2010). (49) Membrane Proteins of Known 3D Structure. http://blanco.biomol.uci.edu/Membrane_Proteins_xtal.html (July 11, 2010). (50) Michel, H. Trends Biochem Sci 1983, 8, 56. (51) Rosenbusch, J. P. Journal of Structural Biology 1990, 104, 134. (52) Garavito, R. M.; Picot, D. Methods 1990, 1, 57. (53) Kulkarni, C. V. 2010; Vol. 12, p 237. (54) Landau, E. M.; Rosenbusch, J. P. P Natl Acad Sci USA 1996, 93, 14532. (55) Pebay-Peyroula, E.; Rummel, G.; Rosenbusch, J. P.; Landau, E. M. Science 1997, 277, 1676. (56) Cherezov, V.; Liu, W.; Derrick, J. P.; Luan, B.; Aksimentiev, A.; Katritch, V.; Caffrey, M. Proteins: Structure, Function, and Bioinformatics 2008, 71, 24. (57) Cherezov, V.; Rosenbaum, D. M.; Hanson, M. A.; Rasmussen, S. G. F.; Thian, F. S.; Kobilka, T. S.; Choi, H. J.; Kuhn, P.; Weis, W. I.; Kobilka, B. K.; Stevens, R. C. Science 2007, 318, 1258. (58) Cherezov, V.; Yamashita, E.; Liu, W.; Zhalnina, M.; Cramer, W. A.; Caffrey, M. J Mol Biol 2006, 364, 716. (59) Jaakola, V. P.; Griffith, M. T.; Hanson, M. A.; Cherezov, V.; Chien, E. Y. T.; Lane, J. R.; IJzerman, A. P.; Stevens, R. C. Science 2008, 322, 1211. (60) Rosenbaum, D. M.; Cherezov, V.; Hanson, M. A.; Rasmussen, S. G. F.; Thian, F. S.; Kobilka, T. S.; Choi, H. J.; Yao, X. J.; Weis, W. I.; Stevens, R. C.; Kobilka, B. K. Science 2007, 318, 1266. (61) Wacker, D.; Fenalti, G.; Brown, M. A.; Katritch, V.; Abagyan, R.; Cherezov, V.; Stevens, R. C. J Am Chem Soc 2010, 132, 11443. (62) Höfer, N.; Aragão, D.; Caffrey, M. Biophys J 2010, 99, L23. (63) Li, L.; Ismagilov, R. F. Ann Rev Biophys 2010. (64) Pal, R.; Yang, M.; Lin, R.; Johnson, B. N.; Srivastava, N.; Razzacki, S. Z.; Chomistek, K. J.; Heldsinger, D. C.; Haque, R. M.; Ugaz, V. M.; Thwar, P. K.; Chen, Z.; Alfano, K.; Yim, M. B.; Krishnan, M.; Fuller, A. O.; Larson, R. G.; Burke, D. T.; Burns, M. A. Lab Chip 2005, 5, 1024. (65) Jayashree, R. S.; Gancs, L.; Choban, E. R.; Primak, A.; Natarajan, D.; Markoski, L. J.; Kenis, P. J. A. J Am Chem Soc 2005, 127, 16758. (66) Wootton, R. C. R.; deMello, A. J. Chem Commun 2004, 266. (67) McPherson, A. J Appl Crystallogr 2000, 33, 397.

Microfluidic Paper-based Devices For Bioanalytical Applications.

Paper has become increasingly recognized as a very interesting substrate for the construction of microfluidic devices, with potential application in a variety of areas, including health diagnosis, environmental monitoring, immunoassays and food safety. The aim of this review is to present a short history of analytical systems constructed from paper, summarize the main advantages and disadvantages of fabrication techniques, exploit alternative methods of detection such as colorimetric, electrochemical, photoelectrochemical, chemiluminescence and electrochemiluminescence, as well as to take a closer look at the novel achievements in the field of bioanalysis published during the last 2 years. Finally, the future trends for production of such devices are discussed.

Contrasting Effects Of Land Use Intensity And Exotic Host Plants On The Specialization Of Interactions In Plant-herbivore Networks.

Human land use tends to decrease the diversity of native plant species and facilitate the invasion and establishment of exotic ones. Such changes in land use and plant community composition usually have negative impacts on the assemblages of native herbivorous insects. Highly specialized herbivores are expected to be especially sensitive to land use intensification and the presence of exotic plant species because they are neither capable of consuming alternative plant species of the native flora nor exotic plant species. Therefore, higher levels of land use intensity might reduce the proportion of highly specialized herbivores, which ultimately would lead to changes in the specialization of interactions in plant-herbivore networks. This study investigates the community-wide effects of land use intensity on the degree of specialization of 72 plant-herbivore networks, including effects mediated by the increase in the proportion of exotic plant species. Contrary to our expectation, the net effect of land use intensity on network specialization was positive. However, this positive effect of land use intensity was partially canceled by an opposite effect of the proportion of exotic plant species on network specialization. When we analyzed networks composed exclusively of endophagous herbivores separately from those composed exclusively of exophagous herbivores, we found that only endophages showed a consistent change in network specialization at higher land use levels. Altogether, these results indicate that land use intensity is an important ecological driver of network specialization, by way of reducing the local host range of herbivore guilds with highly specialized feeding habits. However, because the effect of land use intensity is offset by an opposite effect owing to the proportion of exotic host species, the net effect of land use in a given herbivore assemblage will likely depend on the extent of the replacement of native host species with exotic ones.

In vitro Microfluidic Model For The Study Of Vaso-occlusive Processes.

Vaso-occlusion, responsible for much of the morbidity of sickle-cell disease, is a complex multicellular process, apparently triggered by leukocyte adhesion to the vessel wall. The microcirculation represents a major site of leukocyte-endothelial interactions and vaso-occlusive processes. We have developed a biochip with subdividing interconnecting microchannels that decrease in size (40 μm to 10 μm in width), for use in conjunction with a precise microfluidic device, to mimic cell flow and adhesion through channels of sizes that approach those of the microcirculation. The biochips were utilized to observe the dynamics of the passage of neutrophils and red blood cells, isolated from healthy and sickle-cell anemia (SCA) individuals, through laminin or endothelial adhesion molecule-coated microchannels at physiologically relevant rates of flow and shear stress. Obstruction of E-selectin/intercellular adhesion molecule 1-coated biochip microchannels by SCA neutrophils was significantly greater than that observed for healthy neutrophils, particularly in the microchannels of 40-15 μm in width. Whereas SCA red blood cells alone did not significantly adhere to, or obstruct, microchannels, mixed suspensions of SCA neutrophils and red blood cells significantly adhered to and obstructed laminin-coated channels. Results from this in vitro microfluidic model support a primary role for leukocytes in the initiation of SCA occlusive processes in the microcirculation. This assay represents an easy-to-use and reproducible in vitro technique for understanding molecular mechanisms and cellular interactions occurring in subdividing microchannels of widths approaching those observed in the microvasculature. The assay could hold potential for testing drugs developed to inhibit occlusive mechanisms such as those observed in SCA and thrombotic diseases.

Complex networks: the key to systems biology

Though introduced recently, complex networks research has grown steadily because of its potential to represent, characterize and model a wide range of intricate natural systems and phenomena. Because of the intrinsic complexity and systemic organization of life, complex networks provide a specially promising framework for systems biology investigation. The current article is an up-to-date review of the major developments related to the application of complex networks in biology, with special attention focused on the more recent literature. The main concepts and models of complex networks are presented and illustrated in an accessible fashion. Three main types of networks are covered: transcriptional regulatory networks, protein-protein interaction networks and metabolic networks. The key role of complex networks for systems biology is extensively illustrated by several of the papers reviewed.

Neural networks and statistical analysis for classification of corneal videokeratography maps based on Zernike coefficients: a quantitative comparison

PURPOSE: The main goal of this study was to develop and compare two different techniques for classification of specific types of corneal shapes when Zernike coefficients are used as inputs. A feed-forward artificial Neural Network (NN) and discriminant analysis (DA) techniques were used. METHODS: The inputs both for the NN and DA were the first 15 standard Zernike coefficients for 80 previously classified corneal elevation data files from an Eyesys System 2000 Videokeratograph (VK), installed at the Departamento de Oftalmologia of the Escola Paulista de Medicina, São Paulo. The NN had 5 output neurons which were associated with 5 typical corneal shapes: keratoconus, with-the-rule astigmatism, against-the-rule astigmatism, "regular" or "normal" shape and post-PRK. RESULTS: The NN and DA responses were statistically analyzed in terms of precision ([true positive+true negative]/total number of cases). Mean overall results for all cases for the NN and DA techniques were, respectively, 94% and 84.8%. CONCLUSION: Although we used a relatively small database, results obtained in the present study indicate that Zernike polynomials as descriptors of corneal shape may be a reliable parameter as input data for diagnostic automation of VK maps, using either NN or DA.

Evaluation of follicular lymphoid depletion in the Bursa of Fabricius: an alternative methodology using digital image analysis and artificial neural networks

Fifty Bursa of Fabricius (BF) were examined by conventional optical microscopy and digital images were acquired and processed using Matlab® 6.5 software. The Artificial Neuronal Network (ANN) was generated using Neuroshell® Classifier software and the optical and digital data were compared. The ANN was able to make a comparable classification of digital and optical scores. The use of ANN was able to classify correctly the majority of the follicles, reaching sensibility and specificity of 89% and 96%, respectively. When the follicles were scored and grouped in a binary fashion the sensibility increased to 90% and obtained the maximum value for the specificity of 92%. These results demonstrate that the use of digital image analysis and ANN is a useful tool for the pathological classification of the BF lymphoid depletion. In addition it provides objective results that allow measuring the dimension of the error in the diagnosis and classification therefore making comparison between databases feasible.

Classification of breast masses using a committee machine of artificial neural networks

This work proposes a new approach using a committee machine of artificial neural networks to classify masses found in mammograms as benign or malignant. Three shape factors, three edge-sharpness measures, and 14 texture measures are used for the classification of 20 regions of interest (ROIs) related to malignant tumors and 37 ROIs related to benign masses. A group of multilayer perceptrons (MLPs) is employed as a committee machine of neural network classifiers. The classification results are reached by combining the responses of the individual classifiers. Experiments involving changes in the learning algorithm of the committee machine are conducted. The classification accuracy is evaluated using the area A. under the receiver operating characteristics (ROC) curve. The A, result for the committee machine is compared with the A, results obtained using MLPs and single-layer perceptrons (SLPs), as well as a linear discriminant analysis (LDA) classifier Tests are carried out using the student's t-distribution. The committee machine classifier outperforms the MLP SLP, and LDA classifiers in the following cases: with the shape measure of spiculation index, the A, values of the four methods are, in order 0.93, 0.84, 0.75, and 0.76; and with the edge-sharpness measure of acutance, the values are 0.79, 0.70, 0.69, and 0.74. Although the features with which improvement is obtained with the committee machines are not the same as those that provided the maximal value of A(z) (A(z) = 0.99 with some shape features, with or without the committee machine), they correspond to features that are not critically dependent on the accuracy of the boundaries of the masses, which is an important result. (c) 2008 SPIE and IS&T.

Fully Connected PLL Networks: How Filter Determines the Number of Nodes

Synchronization plays an important role in telecommunication systems, integrated circuits, and automation systems. Formerly, the masterslave synchronization strategy was used in the great majority of cases due to its reliability and simplicity. Recently, with the wireless networks development, and with the increase of the operation frequency of integrated circuits, the decentralized clock distribution strategies are gaining importance. Consequently, fully connected clock distribution systems with nodes composed of phase-locked loops (PLLs) appear as a convenient engineering solution. In this work, the stability of the synchronous state of these networks is studied in two relevant situations: when the node filters are first-order lag-lead low-pass or when the node filters are second-order low-pass. For first-order filters, the synchronous state of the network shows to be stable for any number of nodes. For second-order filter, there is a superior limit for the number of nodes, depending on the PLL parameters. Copyright (C) 2009 Atila Madureira Bueno et al.

A reliable measure of similarity based on dependency for short time series: an application to gene expression networks

Background: Microarray techniques have become an important tool to the investigation of genetic relationships and the assignment of different phenotypes. Since microarrays are still very expensive, most of the experiments are performed with small samples. This paper introduces a method to quantify dependency between data series composed of few sample points. The method is used to construct gene co-expression subnetworks of highly significant edges. Results: The results shown here are for an adapted subset of a Saccharomyces cerevisiae gene expression data set with low temporal resolution and poor statistics. The method reveals common transcription factors with a high confidence level and allows the construction of subnetworks with high biological relevance that reveals characteristic features of the processes driving the organism adaptations to specific environmental conditions. Conclusion: Our method allows a reliable and sophisticated analysis of microarray data even under severe constraints. The utilization of systems biology improves the biologists ability to elucidate the mechanisms underlying celular processes and to formulate new hypotheses.

The Missing Part of Seed Dispersal Networks: Structure and Robustness of Bat-Fruit Interactions

Mutualistic networks are crucial to the maintenance of ecosystem services. Unfortunately, what we know about seed dispersal networks is based only on bird-fruit interactions. Therefore, we aimed at filling part of this gap by investigating bat-fruit networks. It is known from population studies that: (i) some bat species depend more on fruits than others, and (ii) that some specialized frugivorous bats prefer particular plant genera. We tested whether those preferences affected the structure and robustness of the whole network and the functional roles of species. Nine bat-fruit datasets from the literature were analyzed and all networks showed lower complementary specialization (H(2)' = 0.3760.10, mean 6 SD) and similar nestedness (NODF = 0.5660.12) than pollination networks. All networks were modular (M=0.32 +/- 0.07), and had on average four cohesive subgroups (modules) of tightly connected bats and plants. The composition of those modules followed the genus-genus associations observed at population level (Artibeus-Ficus, Carollia-Piper, and Sturnira-Solanum), although a few of those plant genera were dispersed also by other bats. Bat-fruit networks showed high robustness to simulated cumulative removals of both bats (R = 0.55 +/- 0.10) and plants (R = 0.68 +/- 0.09). Primary frugivores interacted with a larger proportion of the plants available and also occupied more central positions; furthermore, their extinction caused larger changes in network structure. We conclude that bat-fruit networks are highly cohesive and robust mutualistic systems, in which redundancy is high within modules, although modules are complementary to each other. Dietary specialization seems to be an important structuring factor that affects the topology, the guild structure and functional roles in bat-fruit networks.

Social Networks Shape the Transmission Dynamics of Hepatitis C Virus

Hepatitis C virus (HCV) infects 170 million people worldwide, and is a major public health problem in Brazil, where over 1% of the population may be infected and where multiple viral genotypes co-circulate. Chronically infected individuals are both the source of transmission to others and are at risk for HCV-related diseases, such as liver cancer and cirrhosis. Before the adoption of anti-HCV control measures in blood banks, this virus was mainly transmitted via blood transfusion. Today, needle sharing among injecting drug users is the most common form of HCV transmission. Of particular importance is that HCV prevalence is growing in non-risk groups. Since there is no vaccine against HCV, it is important to determine the factors that control viral transmission in order to develop more efficient control measures. However, despite the health costs associated with HCV, the factors that determine the spread of virus at the epidemiological scale are often poorly understood. Here, we sequenced partial NS5b gene sequences sampled from blood samples collected from 591 patients in Sao Paulo state, Brazil. We show that different viral genotypes entered Sao Paulo at different times, grew at different rates, and are associated with different age groups and risk behaviors. In particular, subtype 1b is older and grew more slowly than subtypes 1a and 3a, and is associated with multiple age classes. In contrast, subtypes 1a and 3b are associated with younger people infected more recently, possibly with higher rates of sexual transmission. The transmission dynamics of HCV in Sao Paulo therefore vary by subtype and are determined by a combination of age, risk exposure and underlying social network. We conclude that social factors may play a key role in determining the rate and pattern of HCV spread, and should influence future intervention policies.