Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Investigating wetting characteristics on microstructured surfaces for superhydrophobicity and metal microcasting

The engineering of liquid behavior on surfaces is important for infrastructure, transportation, manufacturing, and sensing. Surfaces can be rendered superhydrophobic by microstructuring, and superhydrophobic devices could lead to practical corrosion inhibition, self-cleaning, fluid flow control, and surface drag reduction. To more fully understand how liquid interacts with microstructured surfaces, this dissertation introduces a direct method for determining droplet solid-liquid-vapor interfacial geometry on microstructured surfaces. The technique performs metrology on molten metal droplets deposited onto microstructured surfaces and then frozen. Unlike other techniques, this visualization technique can be used on large areas of curved and opaque microstructured surfaces to determine contact line. This dissertation also presents measurements and models for how curvature and flexing of microstructured polymers affects hydrophobicity. Increasing curvature of microstructured surfaces leads to decreased slide angle for liquid droplets suspended on the surface asperities. For a surface with regularly spaced asperities, as curvature becomes more positive, droplets suspended on the tops of asperities are suspended on fewer asperities. Curvature affects superhydrophobicity because microscopic curvature changes solid-liquid interaction, pitch is altered, and curvature changes the shape of the three phase contact line. This dissertation presents a model of droplet interactions with curved microstructured surfaces that can be used to design microstructure geometries that maintain the suspension of a droplet when curved surfaces are covered with microstructured polymers. Controlling droplet dynamics could improve microfluidic devices and the shedding of liquids from expensive equipment, preventing corrosion and detrimental performance. This dissertation demonstrates redirection of dynamic droplet spray with anisotropic microstructures. Superhydrophobic microstructured surfaces can be economically fabricated using metal embossing masters, so this dissertation describes casting-based microfabrication of metal microstructures and nanostructures. Low melting temperature metal was cast into flexible silicone molds which were themselves cast from microfabricated silicon templates. The flexibility of the silicone mold permits casting of curved surfaces, which this dissertation demonstrates by fabricating a cylindrical metal roller with microstructures. The metal microstructures can be in turn used as a reusable molding tool. This dissertation also describes an industrial investment casting process to produce aluminum molds having integrated microstructures. Unlike conventional micromolding tools, the aluminum mold was large and had complex curved surfaces. The aluminum was cast into curved microstructured ceramic molds which were themselves cast from curved microstructured rubber. Many structures were successfully cast into the aluminum with excellent replication fidelity, including circular, square, and triangular holes. This dissertation demonstrates molding of large, curved surfaces having surface microstructures using the aluminum mold. This work contributes a more full understanding of the phenomenon of superhydrophobicity and techniques for the economic fabrication of superhydrophobic microstructures.

Measurements of small ice crystals in arctic cirrus: observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC)

In-situ observations on the size and shape of particles in arctic cirrus are less common than those in mid-latitude and tropical cirrus with considerable uncertainty about the contributions of small ice crystals (maximum dimension D<50 µm) to the mass and radiative properties that impact radiative forcing. In situ measurements of small ice crystals in arctic cirrus were made during the Indirect and Semi-Direct Aerosol Campaign (ISDAC) in April 2008 during transits of the National Research Council of Canada Convair-580 between Fairbanks and Barrow, Alaska and during Mixed Phase Arctic Cloud Experiment (MPACE) in October 2004 with the University of North Dakota (UND) Citation over Barrow, Alaska. Concentrations of small ice crystals with D < 50 μm from a Cloud and Aerosol Spectrometer (CAS), a Cloud Droplet Probe (CDP), a Forward Scattering Spectrometer Probe (FSSP), and a two-dimensional stereo probe (2DS) were compared as functions of the concentrations of crystals with D > 100 μm measured by a Cloud Imaging Probe (CIP) and two-dimensional stereo probe (2DS) in order to assess whether the shattering of large ice crystals on protruding components of different probes artificially amplified measurements of small ice crystal concentrations. The dependence of the probe comparison on other variables as CIP N>100 (number concentrations greater than diameter D>100 μm),temperature, relative humidity respect to ice (RHice), dominant habit from the Cloud Particle Imager (CPI), aircraft roll, pitch, true air speed and angle of attack was examined to understand potential causes of discrepancies between probe concentrations. Data collected by these probes were also compared against the data collected by a CAS, CDP and CIP during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) and by a CAS and 2DS during the Tropical Composition, Cloud and Climate Coupling (TC4) missions. During ISDAC, the CAS and FSSP both overestimated measurements of small ice crystals compared to both the CDP and 2DS by 1-2 orders of magnitude. Further, the amount of overestimation increased with the concentrations from the CIP2 (N>100 > 0.1 L-1). There was an unexplained discrepancy in concentrations of small crystals between the CDP and 2DS during ISDAC. In addition, there was a strong dependence on RHice of the average ratios of the N3-50, CAS/N3-50,CDP, N3-50, FSSP096/N3-50,CDP, N3-50, CAS/N3-50,FSSP096, N10-50, CDP/N3-50,2DS, N10-50, FSSP096/N10-50,2DS. Continued studies are needed to understand the discrepancy of these probes.

Cryopreservation of avocado (Persea Americana Mill.) using somatic embryos

Avocado genetic resources are currently maintained in the form of field repositories at great cost and risk of natural disasters, pest and diseases. Cryopreservation offers a necessary, complimentary method that is safe, cost-effective and long-term. However, long-term maintenance and regeneration of plantlets from avocado somatic embryos has been a major barrier in the development of new avocado cultivars. In this study, two protocols for vitrification-based cryopreservation of avocado somatic embryos were investigated. Globular somatic embryos of two avocado cultivars were tested, revealing cultivar-dependent differences in desiccation tolerance and subsequent freezing resistance, possibly attributed to their size and culture age. A two-step regeneration system, involving an intermediate liquid phase step between subcultures in solid medium, significantly enhanced shoot development from somatic embryo tissue. This work will add considerable value towards cryopreservation of avocado somatic embryos for germplasm conservation and the generation of new and improved avocado cultivars.

Computational modelling of non-equilibrium condensing steam flows in low-pressure steam turbines

The steam turbines play a significant role in global power generation. Especially, research on low pressure (LP) steam turbine stages is of special importance for steam turbine man- ufactures, vendors, power plant owners and the scientific community due to their lower efficiency than the high pressure steam turbine stages. Because of condensation, the last stages of LP turbine experience irreversible thermodynamic losses, aerodynamic losses and erosion in turbine blades. Additionally, an LP steam turbine requires maintenance due to moisture generation, and therefore, it is also affecting on the turbine reliability. Therefore, the design of energy efficient LP steam turbines requires a comprehensive analysis of condensation phenomena and corresponding losses occurring in the steam tur- bine either by experiments or with numerical simulations. The aim of the present work is to apply computational fluid dynamics (CFD) to enhance the existing knowledge and understanding of condensing steam flows and loss mechanisms that occur due to the irre- versible heat and mass transfer during the condensation process in an LP steam turbine. Throughout this work, two commercial CFD codes were used to model non-equilibrium condensing steam flows. The Eulerian-Eulerian approach was utilised in which the mix- ture of vapour and liquid phases was solved by Reynolds-averaged Navier-Stokes equa- tions. The nucleation process was modelled with the classical nucleation theory, and two different droplet growth models were used to predict the droplet growth rate. The flow turbulence was solved by employing the standard k-ε and the shear stress transport k-ω turbulence models. Further, both models were modified and implemented in the CFD codes. The thermodynamic properties of vapour and liquid phases were evaluated with real gas models. In this thesis, various topics, namely the influence of real gas properties, turbulence mod- elling, unsteadiness and the blade trailing edge shape on wet-steam flows, are studied with different convergent-divergent nozzles, turbine stator cascade and 3D turbine stator-rotor stage. The simulated results of this study were evaluated and discussed together with the available experimental data in the literature. The grid independence study revealed that an adequate grid size is required to capture correct trends of condensation phenomena in LP turbine flows. The study shows that accurate real gas properties are important for the precise modelling of non-equilibrium condensing steam flows. The turbulence modelling revealed that the flow expansion and subsequently the rate of formation of liquid droplet nuclei and its growth process were affected by the turbulence modelling. The losses were rather sensitive to turbulence modelling as well. Based on the presented results, it could be observed that the correct computational prediction of wet-steam flows in the LP turbine requires the turbulence to be modelled accurately. The trailing edge shape of the LP turbine blades influenced the liquid droplet formulation, distribution and sizes, and loss generation. The study shows that the semicircular trailing edge shape predicted the smallest droplet sizes. The square trailing edge shape estimated greater losses. The analysis of steady and unsteady calculations of wet-steam flow exhibited that in unsteady simulations, the interaction of wakes in the rotor blade row affected the flow field. The flow unsteadiness influenced the nucleation and droplet growth processes due to the fluctuation in the Wilson point.

Shape and chemically anisotropic colloidal particles: A theoretical study of their structure, phase behavior, and slow dynamics

A detailed non-equilibrium state diagram of shape-anisotropic particle fluids is constructed. The effects of particle shape are explored using Naive Mode Coupling Theory (NMCT), and a single particle Non-linear Langevin Equation (NLE) theory. The dynamical behavior of non-ergodic fluids are discussed. We employ a rotationally frozen approach to NMCT in order to determine a transition to center of mass (translational) localization. Both ideal and kinetic glass transitions are found to be highly shape dependent, and uniformly increase with particle dimensionality. The glass transition volume fraction of quasi 1- and 2- dimensional particles fall monotonically with the number of sites (aspect ratio), while 3-dimensional particles display a non-monotonic dependence of glassy vitrification on the number of sites. Introducing interparticle attractions results in a far more complex state diagram. The ideal non-ergodic boundary shows a glass-fluid-gel re-entrance previously predicted for spherical particle fluids. The non-ergodic region of the state diagram presents qualitatively different dynamics in different regimes. They are qualified by the different behaviors of the NLE dynamic free energy. The caging dominated, repulsive glass regime is characterized by long localization lengths and barrier locations, dictated by repulsive hard core interactions, while the bonding dominated gel region has short localization lengths (commensurate with the attraction range), and barrier locations. There exists a small region of the state diagram which is qualified by both glassy and gel localization lengths in the dynamic free energy. A much larger (high volume fraction, and high attraction strength) region of phase space is characterized by short gel-like localization lengths, and long barrier locations. The region is called the attractive glass and represents a 2-step relaxation process whereby a particle first breaks attractive physical bonds, and then escapes its topological cage. The dynamic fragility of fluids are highly particle shape dependent. It increases with particle dimensionality and falls with aspect ratio for quasi 1- and 2- dimentional particles. An ultralocal limit analysis of the NLE theory predicts universalities in the behavior of relaxation times, and elastic moduli. The equlibrium phase diagram of chemically anisotropic Janus spheres and Janus rods are calculated employing a mean field Random Phase Approximation. The calculations for Janus rods are corroborated by the full liquid state Reference Interaction Site Model theory. The Janus particles consist of attractive and repulsive regions. Both rods and spheres display rich phase behavior. The phase diagrams of these systems display fluid, macrophase separated, attraction driven microphase separated, repulsion driven microphase separated and crystalline regimes. Macrophase separation is predicted in highly attractive low volume fraction systems. Attraction driven microphase separation is charaterized by long length scale divergences, where the ordering length scale determines the microphase ordered structures. The ordering length scale of repulsion driven microphase separation is determined by the repulsive range. At the high volume fractions, particles forgo the enthalpic considerations of attractions and repulsions to satisfy hard core constraints and maximize vibrational entropy. This results in site length scale ordering in rods, and the sphere length scale ordering in Janus spheres, i.e., crystallization. A change in the Janus balance of both rods and spheres results in quantitative changes in spinodal temperatures and the position of phase boundaries. However, a change in the block sequence of Janus rods causes qualitative changes in the type of microphase ordered state, and induces prominent features (such as the Lifshitz point) in the phase diagrams of these systems. A detailed study of the number of nearest neighbors in Janus rod systems reflect a deep connection between this local measure of structure, and the structure factor which represents the most global measure of order.

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.

Adipocyte size fluctuation, mechano-active lipid droplets and caveolae

International audience

Studies on bovine oocytes and ovarian tissues cryopreservation : fundamental and applied aspects

Tese de Doutoramento, Ciências Agrárias (Reprodução Animal), 26 de Junho de 2013, Universidade dos Açores.

Chemical and physical effects of ultrasound: sonoluminescence and materials

When a liquid is irradiated with ultrasound, acoustic cavitation (the formation, growth, and implosive collapse of bubbles in liquids irradiated with ultrasound) generally occurs. This is the phenomenon responsible for the driving of chemical reactions (sonochemistry) and the emission of light (sonoluminescence). The implosive collapse of bubbles in liquids results in an enormous concentration of sound energy into compressional heating of the bubble contents. Therefore, extreme chemical and physical conditions are generated during cavitation. The study of multibubble sonoluminescence (MBSL) and single-bubble sonoluminescence (SBSL) in exotic liquids such as sulfuric acid (H2SO4) and phosphoric acid (H3PO4) leads to useful information regarding the intracavity conditions during bubble collapse. Distinct sonoluminescing bubble populations were observed from the intense orange and blue-white emissions by doping H2SO4 and H3PO4 with sodium salts, which provides the first experimental evidence for the injected droplet model over the heated-shell model for cavitation. Effective emission temperatures measured based on excited OH• and PO• emission indicate that there is a temperature inhomogeneity during MBSL in 85% H3PO4. The formation of a temperature inhomogeneity is due to the existence of different cavitating bubble populations: asymmetric collapsing bubbles contain liquid droplets and spherical collapsing bubbles do not contain liquid droplets. Strong molecular emission from SBSL in 65% H3PO4 have been obtained and used as a spectroscopic probe to determine the cavitation temperatures. It is found that the intracavity temperatures are dependent on the applied acoustic pressures and the thermal conductivities of the dissolved noble gases. The chemical and physical effects of ultrasound can be used for materials synthesis. Highly reactive species, including HO2•, H•, and OH• (or R• after additives react with OH•), are formed during aqueous sonolysis as a consequence of the chemical effects of ultrasound. Reductive species can be applied to synthesis of water-soluble fluorescent silver nanoclusters in the presence of a suitable stabilizer or capping agent. The optical and fluorescent properties of the Ag nanoclusters can be easily controlled by the synthetic conditions such as the sonication time, the stoichiometry of the carboxylate groups to Ag+, and the polymer molecular weight. The chemical and physical effects of ultrasound can be combined to prepare polymer functionalized graphenes from graphites and a reactive solvent, styrene. The physical effects of ultrasound are used to exfoliate graphites to graphenes while the chemical effects of ultrasound are used to induce the polymerization of styrene which can then functionalize graphene sheets via radical coupling. The prepared polymer functionalized graphenes are highly stable in common organic solvents like THF, CHCl3, and DMF. Ultrasonic spray pyrolysis (USP) is used to prepare porous carbon spheres using energetic alkali propiolates as the carbon precursors. In this synthesis, metal salts are generated in situ, introducing porous structures into the carbon spheres. When different alkali salts or their mixtures are used as the precursor, carbon spheres with different morphologies and structures are obtained. The different precursor decomposition pathways are responsible for the observed structural difference. Such prepared carbon materials have high surface area and are thermally stable, making them potentially useful for catalytic supports, adsorbents, or for other applications by integrating other functional materials into their pores.

Solid-liquid interactions in microscale structures and devices

Liquid-solid interactions become important as dimensions approach mciro/nano-scale. This dissertation focuses on liquid-solid interactions in two distinct applications: capillary driven self-assembly of thin foils into 3D structures, and droplet wetting of hydrophobic micropatterned surfaces. The phenomenon of self-assembly of complex structures is common in biological systems. Examples include self-assembly of proteins into macromolecular structures and self-assembly of lipid bilayer membranes. The principles governing this phenomenon have been applied to induce self-assembly of millimeter scale Si thin films into spherical and other 3D structures, which are then integrated into light-trapping photovoltaic (PV) devices. Motivated by this application, we present a generalized analytical study of the self-folding of thin plates into deterministic 3D shapes, through fluid-solid interactions, to be used as PV devices. This study consists of developing a model using beam theory, which incorporates the two competing components — a capillary force that promotes folding and the bending rigidity of the foil that resists folding into a 3D structure. Through an equivalence argument of thin foils of different geometry, an effective folding parameter, which uniquely characterizes the driving force for folding, has been identified. A criterion for spontaneous folding of an arbitrarily shaped 2D foil, based on the effective folding parameter, is thus established. Measurements from experiments using different materials and predictions from the model match well, validating the assumptions used in the analysis. As an alternative to the mechanics model approach, the minimization of the total free energy is employed to investigate the interactions between a fluid droplet and a flexible thin film. A 2D energy functional is proposed, comprising the surface energy of the fluid, bending energy of the thin film and gravitational energy of the fluid. Through simulations with Surface Evolver, the shapes of the droplet and the thin film at equilibrium are obtained. A critical thin film length necessary for complete enclosure of the fluid droplet, and hence successful self-assembly into a PV device, is determined and compared with the experimental results and mechanics model predictions. The results from the modeling and energy approaches and the experiments are all consistent. Superhydrophobic surfaces, which have unique properties including self-cleaning and water repelling are desired in many applications. One excellent example in nature is the lotus leaf. To fabricate these surfaces, well designed micro/nano- surface structures are often employed. In this research, we fabricate superhydrophobic micropatterned Polydimethylsiloxane (PDMS) surfaces composed of micropillars of various sizes and arrangements by means of soft lithography. Both anisotropic surfaces, consisting of parallel grooves and cylindrical pillars in rectangular lattices, and isotropic surfaces, consisting of cylindrical pillars in square and hexagonal lattices, are considered. A novel technique is proposed to image the contact line (CL) of the droplet on the hydrophobic surface. This technique provides a new approach to distinguish between partial and complete wetting. The contact area between droplet and microtextured surface is then measured for a droplet in the Cassie state, which is a state of partial wetting. The results show that although the droplet is in the Cassie state, the contact area does not necessarily follow Cassie model predictions. Moreover, the CL is not circular, and is affected by the micropatterns, in both isotropic and anisotropic cases. Thus, it is suggested that along with the contact angle — the typical parameter reported in literature quantifying wetting, the size and shape of the contact area should also be presented. This technique is employed to investigate the evolution of the CL on a hydrophobic micropatterned surface in the cases of: a single droplet impacting the micropatterned surface, two droplets coalescing on micropillars, and a receding droplet resting on the micropatterned surface. Another parameter which quantifies hydrophobicity is the contact angle hysteresis (CAH), which indicates the resistance of the surface to the sliding of a droplet with a given volume. The conventional methods of using advancing and receding angles or tilting stage to measure the resistance of the micropatterned surface are indirect, without mentioning the inaccuracy due to the discrete and stepwise motion of the CL on micropillars. A micronewton force sensor is utilized to directly measure the resisting force by dragging a droplet on a microtextured surface. Together with the proposed imaging technique, the evolution of the CL during sliding is also explored. It is found that, at the onset of sliding, the CL behaves as a linear elastic solid with a constant stiffness. Afterwards, the force first increases and then decreases and reaches a steady state, accompanied with periodic oscillations due to regular pinning and depinning of the CL. Both the maximum and steady state forces are primarily dependent on area fractions of the micropatterned surfaces in our experiment. The resisting force is found to be proportional to the number of pillars which pin the CL at the trailing edge, validating the assumption that the resistance mainly arises from the CL pinning at the trailing edge. In each pinning-and-depinning cycle during the steady state, the CL also shows linear elastic behavior but with a lower stiffness. The force variation and energy dissipation involved can also be determined. This novel method of measuring the resistance of the micropatterned surface elucidates the dependence on CL pinning and provides more insight into the mechanisms of CAH.

Engineering the transition from non-living to living matter

Re-creating and understanding the origin of life represents one of the major challenges facing the scientific community. We will never know exactly how life started on planet Earth, however, we can reconstruct the most likely chemical pathways that could have contributed to the formation of the first living systems. Traditionally, prebiotic chemistry has investigated the formation of modern life’s precursors and their self-organisation under very specific conditions thought to be ‘plausible’. So far, this approach has failed to produce a living system from the bottom-up. In the work presented herein, two different approaches are employed to explore the transition from inanimate to living matter. The development of microfluidic technology during the last decades has changed the way traditional chemical and biological experiments are performed. Microfluidics allows the handling of low volumes of reagents with very precise control. The use of micro-droplets generated within microfluidic devices is of particular interest to the field of Origins of Life and Artificial Life. Whilst many efforts have been made aiming to construct cell-like compartments from modern biological constituents, these are usually very difficult to handle. However, microdroplets can be easily generated and manipulated at kHz rates, making it suitable for high-throughput experimentation and analysis of compartmentalised chemical reactions. Therefore, we decided to develop a microfluidic device capable of manipulating microdroplets in such a way that they could be efficiently mixed, split and sorted within iterative cycles. Since no microfluidic technology had been developed before in the Cronin Group, the first chapter of this thesis describes the soft lithographic methods and techniques developed to fabricate microfluidic devices. Also, special attention is placed on the generation of water-in-oil microdroplets, and the subsequent modules required for the manipulation of the droplets such as: droplet fusers, splitters, sorters and single/multi-layer micromechanical valves. Whilst the first part of this thesis describes the development of a microfluidic platform to assist chemical evolution, finding a compatible set of chemical building blocks capable of reacting to form complex molecules with endowed replicating or catalytic activity was challenging. Abstract 10 Hence, the second part of this thesis focuses on potential chemistry that will ultimately possess the properties mentioned above. A special focus is placed on the formation of peptide bonds from unactivated amino acids, despite being one of the greatest challenges in prebiotic chemistry. As opposed to classic prebiotic experiments, in which a specific set of conditions is studied to fit a particular hypothesis, we took a different approach: we explored the effects of several parameters at once on a model polymerisation reaction, without constraints on hypotheses on the nature of optimum conditions or plausibility. This was facilitated by development of a new high-throughput automated platform, allowing the exploration of a much larger number of parameters. This led us to discover that peptide bond formation is less challenging than previously imagined. Having established the right set of conditions under which peptide bond formation was enhanced, we then explored the co-oligomerisation between different amino acids, aiming for the formation of heteropeptides with different structure or function. Finally, we studied the effect of various environmental conditions (rate of evaporation, presence of salts or minerals) in the final product distribution of our oligomeric products.