Dynamics of cell–matrix mechanical interactions in three dimensions

The forces cells apply to their surroundings control biological processes such as growth, adhesion, development, and migration. In the past 20 years, a number of experimental techniques have been developed to measure such cell tractions. These approaches have primarily measured the tractions applied by cells to synthetic two-dimensional substrates, which do not mimic in vivo conditions for most cell types. Many cell types live in a fibrous three-dimensional (3D) matrix environment. While studying cell behavior in such 3D matrices will provide valuable insights for the mechanobiology and tissue engineering communities, no experimental approaches have yet measured cell tractions in a fibrous 3D matrix.

This thesis describes the development and application of an experimental technique for quantifying cellular forces in a natural 3D matrix. Cells and their surrounding matrix are imaged in three dimensions with high speed confocal microscopy. The cell-induced matrix displacements are computed from the 3D image volumes using digital volume correlation. The strain tensor in the 3D matrix is computed by differentiating the displacements, and the stress tensor is computed by applying a constitutive law. Finally, tractions applied by the cells to the matrix are computed directly from the stress tensor.

The 3D traction measurement approach is used to investigate how cells mechanically interact with the matrix in biologically relevant processes such as division and invasion. During division, a single mother cell undergoes a drastic morphological change to split into two daughter cells. In a 3D matrix, dividing cells apply tensile force to the matrix through thin, persistent extensions that in turn direct the orientation and location of the daughter cells. Cell invasion into a 3D matrix is the first step required for cell migration in three dimensions. During invasion, cells initially apply minimal tractions to the matrix as they extend thin protrusions into the matrix fiber network. The invading cells anchor themselves to the matrix using these protrusions, and subsequently pull on the matrix to propel themselves forward.

Lastly, this thesis describes a constitutive model for the 3D fibrous matrix that uses a finite element (FE) approach. The FE model simulates the fibrous microstructure of the matrix and matches the cell-induced matrix displacements observed experimentally using digital volume correlation. The model is applied to predict how cells mechanically sense one another in a 3D matrix. It is found that cell-induced matrix displacements localize along linear paths. These linear paths propagate over a long range through the fibrous matrix, and provide a mechanism for cell-cell signaling and mechanosensing. The FE model developed here has the potential to reveal the effects of matrix density, inhomogeneity, and anisotropy in signaling cell behavior through mechanotransduction.

Synthesis and characterization of 1,1-di-tert-butyldiazene

The synthesis and direct observation of 1,1-di-tert-butyldiazene (16) at -127°C is described. The absorption spectrum of a red solution of 1,1-diazene 16 reveals a structured absorption band with λ max at 506 run (Me_2O, -125°C). The vibrational spacing in S_1 is about 1200 cm^(-1). The excited state of 16 emits weakly with a single maximum at 715 run observed in the fluorescence spectrum (Me_2O:CD_2Cl_2, -196°C). The proton NMR spectrum of 16 occurs as a singlet at 1.41 ppm. Monitoring this NMR absorption at -94^0 ± 2°C shows that 1,1-diazene 16 decomposes with a first-order rate of 1.8 x 10^(-3) sec(-1) to form isobutane, isobutylene and hexarnethylethane. This rate is 10^8 and 10^(34) times faster than the thermal decomposition of the corresponding cis and trans 1,2-di-tert-butyldiazene isomers. The free energy of activation for decomposition of 1,1-diazene 16 is found to be 12.5 ± 0.2 kcal/mol at -94°C which is much lower than the values of 19.1 and 19.4 kcal/lmole calculated at -94°C for N-(2,2,6,6- tetramethylpiperidyl)nitrene (3) and N-(2,2,5,5- tetrarnethylpyrrolidyl)nitrene (4), respectively. This difference between 16 and the cyclic-1,1-diazenes 3 and 4 can be attributed to a large steric interaction between the tert-butyl groups in 1,1-diazene 16.

In order to investigate the nature of the singlet-triplet gap in 1,1-diazenes, 2,5-di-tert-butyl-N-pyrrolynitrene (22) was generated but was found to be too reactive towards dimerization to be persistent. In the presence of dimethylsulfoxide, however, N-pyrrolynitrene (22) can be trapped as N-(2,5-di-tert-butyl- N'-pyrrolyl)dimethylsulfoxirnine (38). N-(2,5-di-tert-butyl-N'-pyrrolyl)dimethylsulfoximine (38-d^6) exchanges with free dimethylsulfoxide at 50°C in solution, presumably by generation and retrapping of pyrrolynitrene 22.

Neuronal mechanism of state control in Drosophila melanogaster

The changes in internal states, such as fear, hunger and sleep affect behavioral responses in animals. In most of the cases, these state-dependent influences are “pleiotropic”: one state affects multiple sensory modalities and behaviors; “scalable”: the strengths and choices of such modulations differ depending on the imminence of demands; and “persistent”: once the state is switched on the effects last even after the internal demands are off. These prominent features of state-control enable animals to adjust their behavioral responses depending on their internal demands. Here, we studied the neuronal mechanisms of state-controls by investigating energy-deprived state (hunger state) and social-deprived state of fruit flies, Drosophila melanogaster, as prototypic models. To approach these questions, we developed two novel methods: a genetically based method to map sites of neuromodulation in the brain and optogenetic tools in Drosophila.

These methods, and genetic perturbations, reveal that the effect of hunger to alter behavioral sensitivity to gustatory cues is mediate by two distinct neuromodulatory pathways. The neuropeptide F (NPF) – dopamine (DA) pathway increases sugar sensitivity under mild starvation, while the adipokinetic hormone (AKH)- short neuropeptide F (sNPF) pathway decreases bitter sensitivity under severe starvation. These two pathways are recruited under different levels of energy demands without any cross interaction. Effects of both of the pathways are mediated by modulation of the gustatory sensory neurons, which reinforce the concept that sensory neurons constitute an important locus for state-dependent control of behaviors. Our data suggests that multiple independent neuromodulatory pathways are underlying pleiotropic and scalable effects of the hunger state.

In addition, using optogenetic tool, we show that the neural control of male courtship song can be separated into probabilistic/biasing, and deterministic/command-like components. The former, but not the latter, neurons are subject to functional modulation by social experience, supporting the idea that they constitute a locus of state-dependent influence. Interestingly, moreover, brief activation of the former, but not the latter, neurons trigger persistent behavioral response for more than 10 min. Altogether, these findings and new tools described in this dissertation offer new entry points for future researchers to understand the neuronal mechanism of state control.

Elemental composition of solar energetic particles

The Low Energy Telescopes on the Voyager spacecraft are used to measure the elemental composition (2 ≤ Z ≤ 28) and energy spectra (5 to 15 MeV /nucleon) of solar energetic particles (SEPs) in seven large flare events. Four flare events are selected which have SEP abundance ratios approximately independent of energy/nucleon. The abundances for these events are compared from flare to flare and are compared to solar abundances from other sources: spectroscopy of the photosphere and corona, and solar wind measurements.

The selected SEP composition results may be described by an average composition plus a systematic flare-to-flare deviation about the average. For each of the four events, the ratios of the SEP abundances to the four-flare average SEP abundances are approximately monotonic functions of nuclear charge Z in the range 6 ≤ Z ≤ 28. An exception to this Z-dependent trend occurs for He, whose abundance relative to Si is nearly the same in all four events.

The four-flare average SEP composition is significantly different from the solar composition determined by photospheric spectroscopy: The elements C, N and O are depleted in SEPs by a factor of about five relative to the elements Na, Mg, Al, Si, Ca, Cr, Fe and Ni. For some elemental abundance ratios (e.g. Mg/O), the difference between SEP and photospheric results is persistent from flare to flare and is apparently not due to a systematic difference in SEP energy/nucleon spectra between the elements, nor to propagation effects which would result in a time-dependent abundance ratio in individual flare events.

The four-flare average SEP composition is in agreement with solar wind abundance results and with a number of recent coronal abundance measurements. The evidence for a common depletion of oxygen in SEPs, the corona and the solar wind relative to the photosphere suggests that the SEPs originate in the corona and that both the SEPs and solar wind sample a coronal composition which is significantly and persistently different from that of the photosphere.

Magnetospheric access of solar particles and the configuration of the distant geomagnetic field

The access of 1.2-40 MeV protons and 0.4-1.0 MeV electrons from interplanetary space to the polar cap regions has been investigated with an experiment on board a low altitude, polar orbiting satellite (OG0-4).

A total of 333 quiet time observations of the electron polar cap boundary give a mapping of the boundary between open and closed geomagnetic field lines which is an order of magnitude more comprehensive than previously available.

Persistent features (north/south asymmetries) in the polar cap proton flux, which are established as normal during solar proton events, are shown to be associated with different flux levels on open geomagnetic field lines than on closed field lines. The pole in which these persistent features are observed is strongly correlated to the sector structure of the interplanetary magnetic field and uncorrelated to the north/south component of this field. The features were observed in the north (south) pole during a negative (positive) sector 91% of the time, while the solar field had a southward component only 54% of the time. In addition, changes in the north/south component have no observable effect on the persistent features.

Observations of events associated with co-rotating regions of enhanced proton flux in interplanetary space are used to establish the characteristics of the 1.2 - 40 MeV proton access windows: the access window for low polar latitudes is near the earth, that for one high polar latitude region is ~250 R_⊕ behind the earth, while that for the other high polar latitude region is ~1750 R_⊕ behind the earth. All of the access windows are of approximately the same extent (~120 R_⊕). The following phenomena contribute to persistent polar cap features: limited interplanetary regions of enhanced flux propagating past the earth, radial gradients in the interplanetary flux, and anisotropies in the interplanetary flux.

These results are compared to the particle access predictions of the distant geomagnetic tail configurations proposed by Michel and Dessler, Dungey, and Frank. The data are consistent with neither the model of Michel and Dessler nor that of Dungey. The model of Frank can yield a consistent access window configuration provided the following constraints are satisfied: the merging rate for open field lines at one polar neutral point must be ~5 times that at the other polar neutral point, related to the solar magnetic field configuration in a consistent fashion, the migration time for open field lines to move across the polar cap region must be the same in both poles, and the open field line merging rate at one of the polar neutral points must be at least as large as that required for almost all the open field lines to have merged in 0 (one hour). The possibility of satisfying these constraints is investigated in some detail.

The role played by interplanetary anisotropies in the observation of persistent polar cap features is discussed. Special emphasis is given to the problem of non-adiabatic particle entry through regions where the magnetic field is changing direction. The degree to which such particle entry can be assumed to be nearly adiabatic is related to the particle rigidity, the angle through which the field turns, and the rate at which the field changes direction; this relationship is established for the case of polar cap observations.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

Mechanics of sediment transport and bedrock erosion in steep landscapes

Erosion is concentrated in steep landscapes such that, despite accounting for only a small fraction of Earth’s total surface area, these areas regulate the flux of sediment to downstream basins, and their rugged morphology records transient changes (or lack thereof) in geologic and climatic forcing. Steep landscapes are geomorphically active; large sediment fluxes and rapid landscape evolution rates can create or destroy habitat for humans and wildlife alike, and landslides, debris flows, and floods common in mountainous areas represent a persistent natural and structural hazard. Despite the central role that steep landscapes play in the geosciences and in landscape management, the processes controlling their evolution have been poorly studied compared to lower-gradient areas. This thesis focuses on the basic mechanics of sediment transport and bedrock incision in steep landscapes, as these are the fundamental processes which set the pace and style of landscape evolution. Chapter 1 examines the spatial distribution of slow-moving landslides; these landslides can dominate sediment fluxes to river networks, but the controls on their occurrence are poorly understood. Using a case-study along the San Andreas Fault, California, I show that slow-moving landslides preferentially occur near the fault, suggesting a rock-strength control on landslide distribution. Chapter 2 provides the first field-measurements of incipient sediment motion in streams steeper than 14% and shows a large influence of slope-dependent flow hydraulics and grain-scale roughness on particle motion. Chapter 3 presents experimental evidence for bedrock erosion by suspended sediment, suggesting that, in contrast to prevailing theoretical predictions, suspension-regime transport in steep streams can be the dominant erosion agent. Steep streams are often characterized by the presence of waterfalls and bedrock steps which can have locally high rates of erosion; Chapters 4 and 5 present newly developed, experimentally validated theory on sediment transport through and bedrock erosion in waterfall plunge pools. Finally, Chapter 6 explores the formation of a bedrock slot canyon where interactions between sediment transport and bedrock incision lead to the formation of upstream-propagating bedrock step-pools and waterfalls.