Restricted exchange microenvironments for cell culture.

Metabolite diffusion in tissues produces gradients and heterogeneous microenvironments that are not captured in standard 2D cell culture models. Here we describe restricted exchange environment chambers (REECs) in which diffusive gradients are formed and manipulated on length scales approximating those found in vivo. In REECs, cells are grown in 2D in an asymmetric chamber (<50 µL) formed between a coverglass and a glass bottom cell culture dish separated by a thin (~100 µm) gasket. Diffusive metabolite exchange between the chamber and bulk media occurs through one or more openings micromachined into the coverglass. Cell-generated concentration gradients form radially in REECs with a single round opening (~200 µm diameter). At steady state only cells within several hundred micrometers of the opening experience metabolite concentrations that permit survival which is analogous to diffusive exchange near a capillary in tissue. The chamber dimensions, the openings' shape, size, and number, and the cellular density and metabolic activity define the gradient structure. For example, two parallel slots above confluent cells produce the 1D equivalent of a spheroid. Using REECs, we found that fibroblasts align along the axis of diffusion while MDCK cells do not. MDCK cells do, however, exhibit significant morphological variations along the diffusive gradient.

Spatial information dynamics during early zebrafish development.

During development inherited information directs growth and specifies the complex spatial organization of cells and molecules. Here we show that a new information metric, the k-space information (kSI), captures the growth and emergence of spatial organization in a developing embryo. Using zebrafish as a model, we quantify the rate of development over the first 24h and demonstrate that important developmental landmarks are associated with well-defined transitions in information dynamics. The rate of development during this time is highest immediately before and after gastrulation, as well early in the segmentation period. We also find that the majority of the information arises from spatial correlations on the length scale of 20-80 µm, but there are contributions from many length scales that change over time. A comparison of the information dynamics in the maternal-zygotic one-eyed pinhead mutant, which is defective in mesoderm induction, with the wild-type embryo shows that the information dynamics diverge near the onset of gastrulation. Subsequently the mutant lacks a peak in the information dynamics that appears to be associated with the formation of trunk somites in the wild-type embryo. These findings provide a common and objective basis by which to quantify spatial organization, compare mutants and quantify developmental dynamics. The kSI can also be applied to any form of developmental data of arbitrary dimensions, and it offers a broad conceptual framework with which to organize the large amounts of data emerging from various sources.

Spatial information analysis of chemotactic trajectories.

During bacterial chemotaxis, a cell acquires information about its environment by sampling changes in the local concentration of a chemoattractant, and then uses that information to bias its motion relative to the source of the chemoattractant. The trajectory of a chemotaxing bacteria is thus a spatial manifestation of the information gathered by the cell. Here we show that a recently developed approach for computing spatial information using Fourier coefficient probabilities, the k-space information (kSI), can be used to quantify the information in such trajectories. The kSI is shown to capture expected responses to gradients of a chemoattractant. We then extend the k-space approach by developing an experimental probability distribution (EPD) that is computed from chemotactic trajectories collected under a reference condition. The EPD accounts for connectivity and other constraints that the nature of the trajectories imposes on the k-space computation. The EPD is used to compute the spatial information from any trajectory of interest, relative to the reference condition. The EPD-based spatial information also captures the expected responses to gradients of a chemoattractant, although the results differ in significant ways from the original kSI computation. In addition, the entropy calculated from the EPD provides a useful measure of trajectory space. The methods developed are highly general, and can be applied to a wide range of other trajectory types as well as non-trajectory data.

Laser inactivation protein patterning of cell culture microenvironments.

Protein micropatterned substrates have emerged as important tools for studying how cells interact with their environment, as well as allowing useful experimental control over, for example, cell shape and cell position on a surface. Here we present a new approach for protein micropatterning in which a focused laser is used to locally inactivate proteins on a protein-coated substrate. By translating the laser relative to the substrate, protein patterns of essentially arbitrary shape can be produced. This approach has a number of useful features. To begin, it is a maskless writing approach. Thus new patterns can be designed and implemented quickly. Laser inactivation can also be performed on a number of different substrate materials, ranging from glass to polydimethylsiloxane. Further, the inactivation is dose dependent, thus complex gradients and other non-uniform distributions of proteins can be produced. Because the focus of the laser can be changed quickly, laser-based patterning can also be applied to substrates with complex topographies or enclosed surfaces--as long as an optical path is available. To demonstrate this capability, protein patterns were made on the inside of small quartz capillary tubes. Patterned substrates produced using laser inactivation constrain cell shape in predictable ways, and we show that these substrates are compatible with a number of different eukaryotic cell lines.

Interactions between planar grafted neurofilament side-arms.

The side-arms of neurofilaments (NFs) have been proposed to be highly disordered, leading to an entropically and electrostatically based repulsion that modulates interfilament spacing. To characterize the behavior of two interacting polymer brushes in a system of this type, we performed molecular dynamics simulations of neurofilament brushes using a four bead reduced amino acid set coarse-grained model. In these simulations, we examined components of the neurofilament brush, NF-L, NF-M, and phosphorylated NF-H (NF-HP), individually. Each protein type was grafted to planar surfaces and simulations were performed for a range of separations of two apposed grafted surfaces. The calculated force-separation curves show the force increases as the reciprocal separation as predicted for polyelectrolyte brushes at high salt. All three systems can be overlapped on a single force-separation curve, which is not expected given the variation in amino acid sequence and charges on the polymers. Examination of structural properties shows scaling behavior in the average brush height, end-to-end distance, and the density interpenetration. Some of this scaling can be understood in terms of treating the NF proteins as effective polyelectrolytes, but some cannot suggesting a distinct polyampholyte behavior. Correlations are found between oppositely charged residues in opposite brushes. However, these correlations are weak in comparison to the strong correlations within each brush. In comparison with recent experimental data that observes condensed and expanded gel states, our results suggest that the condensed state structure involves significant interdigitation of the side-arms.

Computing spatial information from Fourier coefficient distributions.

The spatial relationships between molecules can be quantified in terms of information. In the case of membranes, the spatial organization of molecules in a bilayer is closely related to biophysically and biologically important properties. Here, we present an approach to computing spatial information based on Fourier coefficient distributions. The Fourier transform (FT) of an image contains a complete description of the image, and the values of the FT coefficients are uniquely associated with that image. For an image where the distribution of pixels is uncorrelated, the FT coefficients are normally distributed and uncorrelated. Further, the probability distribution for the FT coefficients of such an image can readily be obtained by Parseval's theorem. We take advantage of these properties to compute the spatial information in an image by determining the probability of each coefficient (both real and imaginary parts) in the FT, then using the Shannon formalism to calculate information. By using the probability distribution obtained from Parseval's theorem, an effective distance from the uncorrelated or most uncertain case is obtained. The resulting quantity is an information computed in k-space (kSI). This approach provides a robust, facile and highly flexible framework for quantifying spatial information in images and other types of data (of arbitrary dimensions). The kSI metric is tested on a 2D Ising model, frequently used as a model for lipid bilayer; and the temperature-dependent phase transition is accurately determined from the spatial information in configurations of the system.

Microelastic properties of lung cell-derived extracellular matrix.

The mechanical properties of the extracellular microenvironment regulate cell behavior, including migration, proliferation and morphogenesis. Although the elastic moduli of synthetic materials have been studied, little is known about the properties of naturally produced extracellular matrix. Here we have utilized atomic force microscopy to characterize the microelastic properties of decellularized cell-derived matrix from human pulmonary fibroblasts. This heterogeneous three-dimensional matrix had an average thickness of 5 ± 0.4 µm and a Young's modulus of 105 ± 14 Pa. Ascorbate treatment of the lung fibroblasts prior to extraction produced a twofold increase in collagen I content, but did not affect the stiffness of the matrices compared with matrices produced in standard medium. However, fibroblast-derived matrices that were crosslinked with glutaraldehyde demonstrated a 67% increase in stiffness. This work provides a microscale characterization of fibroblast-derived matrix mechanical properties. An accurate understanding of native three-dimensional extracellular microenvironments will be essential for controlling cell responses in tissue engineering applications.

Splaying of aliphatic tails plays a central role in barrier crossing during liposome fusion.

The fusion between two lipid bilayers involves crossing a complicated energy landscape. The limiting barrier in the process appears to be between two closely opposed bilayers and the intermediate state where the outer leaflets are fused. We have performed molecular dynamics simulations to characterize the free energy barrier for the fusion of two liposomes and to examine the molecular details of barrier crossing. To capture the slow dynamics of fusion, a model using coarse-grained representations of lipids was used. The fusion between pairs of liposomes was simulated for four systems: DPPC, DOPC, a 3:1 mixture of DPPC/DPPE, and an asymmetric lipid tail system in which one tail of DPPC was reduced to half the length (ASTail). The weighted histogram method was used to compute the free energy as a function of separation distance. The relative barrier heights for these systems was found to be ASTail >> DPPC > DPPC/DPPE > DOPC, in agreement with experimental observations. Further, the free energy curves for all four can be overlaid on a single curve by plotting the free energy versus the surface separation (differing only in the point of fusion). These simulations also confirm that the two main contributions to the free energy barrier are the removal of water between the vesicles and the deformation of the vesicle. The most prominent molecular detail of barrier crossing in all cases examined was the splaying of lipid tails, where initially a single splayed lipid formed a bridge between the two outer leaflets that promotes additional lipid mixing between the vesicles and eventually leads to fusion. The tail splay appears to be closely connected to the energetics of the process. For example, the high barrier for the ASTail is the result of a smaller distance between terminal methyl groups in the splayed molecule. The shortening of this distance requires the liposomes to be closer together, which significantly increases the cost of water removal and bilayer deformation. Before tail splay can initiate fusion, contact must occur between a tail end and the external water. In isolated vesicles, the contact fraction is correlated to the fusogenicity difference between DPPC and DOPC. Moreover, for planar bilayers, the contact fraction is much lower for DPPC, which is consistent with its lack of fusion in giant vesicles. The simulation results show the key roles of lipid tail dynamics in governing the fusion energy landscape.

Conformational dynamics of neurofilament side-arms.

The side-arms of neurofilaments (NFs) have been proposed to be highly disordered, leading to entropic repulsion that modulates interfilament spacing. To gain further insight into the dynamics and organization of the side-arms, we performed molecular dynamics simulations of neurofilament brushes using a coarse-grained model. The density profiles for three NF proteins, NF-L, NF-M, and phosphorylated NF-H (NF-HP), grafted to planar surfaces were calculated and examined as a function of component (salt, residues) and as a function of charge. Analysis of these profiles reveals that the NF with the shortest side arm, NF-L, is disproportionately long compared to the other NFs. The reason for difference is that NF-L is effectively a strong polyelectrolyte, while NF-M and NF-HP are effectively weaker polyelectrolytes. Further, we find cross-correlations between neurofilament side-arms within the brush, even for the NF-L polymers. These correlations occur because of strong attractions between the long sequence repeats of negative residues and the long postive residue repeats and impart a time average structure of the neurofilament brush that deviates from an ideal polymer in a theta solvent.

Nanometer-scale embossing of polydimethylsiloxane.

Microstructured polydimethylsiloxane (PDMS) is an important and widely used material in biology and chemistry. Here we report that micrometer- and nanometer-scale features can be introduced into the surface of PDMS in a process that is functionally equivalent to embossing. We show that surface features <50 nm can be replicated onto the surface of previously cured PDMS at room temperature and at low pressure. This type of embossing can be performed on samples in solution. It also allows one template to be used for many different types of microstructures by changing the embossing time or serial embossing at different alignments. The balance between elastic and plastic properties of the PDMS has the effect of high-pass filtering the features that are captured and produces a sample that is suitable for sensitive surface characterization technologies such as atomic force microscopy. These findings extend the applications of PDMS as well as open the possibility for new uses.

Micropatterns of an extracellular matrix protein with defined information content.

One powerful approach to understanding how cells process spatially variant signals is based on using micropatterned substrates to control the distribution of signaling molecules. However, quantifying spatially complex signals requires an appropriate metric. Here we propose that the Shannon information theory formalism provides a robust and useful way to quantify the organization of proteins in micropatterned systems. To demonstrate the use of informational entropy as a metric, we produced patterns of lines of fibronectin with varying information content. Fibroblasts grown on these patterns were sensitive to very small changes in informational entropy (6.6 bits), and the responses depended on the scale of the pattern.

High fidelity functional patterns of an extracellular matrix protein by electron beam-based inactivation.

Controlling the organization of proteins on surfaces provides a powerful biochemical tool for determining how cells interpret the spatial distribution of local signaling molecules. Here, we describe a general high fidelity approach based on electron beam writing to pattern the functional properties of protein-coated surfaces at length scales ranging from tens of nanometers to millimeters. A silicon substrate is first coated with the extracellular matrix protein fibronectin, which is then locally inactivated by exposure to a highly focused electron beam. Biochemical inactivation of the protein is established by the loss of antibody binding to the fibronectin. Functional inactivation is determined by the inability of cells to spread or form focal adhesions on the inactivated substrate, resulting in cell shapes constrained to the pattern, while they do both (and are unconstrained) on the remaining fibronectin. These protein patterns have very high fidelity, and typical patterns agree with the input dimensions of the pattern to within 2%. Further, the feature edges are well defined and approach molecular dimensions in roughness. Inactivation is shown to be dose dependent with observable suppression of the specific binding at 2 microC cm(-2) and complete removal of biochemical activity at approximately 50 microC cm(-2) for 5 keV electrons. The critical dose for inactivation also depends on accelerating voltage, and complete loss of antibody binding was achieved at approximately 4-7 microC cm(-2) for 1 keV electrons, which corresponds to approximately 50-90 electrons per cross-sectional area of a whole fibronectin dimer and ~2-4 electrons per type III fibronectin domain. AFM analysis of the pattern surfaces revealed that electron beam exposure does not remove appreciable amounts of material from the surface, suggesting that the patterning mechanism involves local inactivation rather than the ablation that has been observed in several organic thin film systems.

Insights into protein structure and function from disorder-complexity space.

Intrinsically disordered proteins have a wide variety of important functional roles. However, the relationship between sequence and function in these proteins is significantly different than that for well-folded proteins. In a previous work, we showed that the propensity to be disordered can be recognized based on sequence composition alone. Here that analysis is furthered by examining the relationship of disorder propensity to sequence complexity, where the metrics for these two properties depend only on composition. The distributions of 40 amino acid peptides from both ordered and disordered proteins are graphed in this disorder-complexity space. An analysis of Swiss-Prot shows that most peptides have high complexity and relatively low disorder. However, there are also an appreciable number of low complexity-high disorder peptides in the database. In contrast, there are no low complexity-low disorder peptides. A similar analysis for peptides in the PDB reveals a much narrower distribution, with few peptides of low complexity and high disorder. In this case, the bounds of the disorder-complexity distribution are well defined and might be used to evaluate the likelihood that a peptide can be crystallized with current methods. The disorder-complexity distributions of individual proteins and sets of proteins grouped by function are also examined. Among individual proteins, there is an enormous variety of distributions that in some cases can be rationalized with regard to function. Groups of functionally related proteins are found to have distributions that are similar within each group but show notable differences between groups. Finally, a pattern matching algorithm is used to search for proteins with particular disorder-complexity distributions. The results suggest that this approach might be used to identify relationships between otherwise dissimilar proteins.

Electron beam patterning of fibronectin nanodots that support focal adhesion formation.

Nanodots of fibronectin which have radii as small as 100 nm and are biofunctional at the cellular level, can be rapidly fabricated in arbitrary spatial patterns using a technique based on electron beam exposure of a protein monolayer with subsequent backfilling of a second protein species.

Substrate effects in poly(ethylene glycol) self-assembled monolayers on granular and flame-annealed gold.

Poly(ethylene glycol) (PEG) self-assembled monolayers (SAMs) are surface coatings that efficiently prevent nonspecific adhesion of biomolecules to surfaces. Here, we report on SAM formation of the PEG thiol CH3O(CH2CH2O)17NHCO(CH2)2SH (PEG(17)) on three types of Au films: thermally evaporated granular Au and two types of Au films from hydrogen flame annealing of granular Au, Au(111), and Au silicide. The different Au surfaces clearly affects the morphology and mechanical properties of the PEG(17) SAM, which is shown by AFM topographs and force distance curves. The two types of SAMs found on flame-annealed Au were denoted "soft" and "hard" due to their difference in stiffness and resistance to scratching by the AFM probe. With the aim of nanometer scale patterning of the PEG(17), the SAMs were exposed by low energy (1 kV) electron beam lithography (EBL). Two distinctly different types of behaviour were observed on the different types of SAM; the soft PEG(17) SAM was destroyed in a self-developing process while material deposition was dominant for the hard PEG(17) SAM.

Directed immobilization of protein-coated nanospheres to nanometer-scale patterns fabricated by electron beam lithography of poly(ethylene glycol) self-assembled monolayers.

Controlling the spatial organization of biomolecules on solid supports with high resolution is important for a wide range of scientific and technological problems. Here we report a study of electron beam lithography (EBL) patterning of a self-assembled monolayer (SAM) of the amide-containing poly(ethylene glycol) (PEG) thiol CH(3)O(CH(2)CH(2)O)(17)NHCO(CH(2))(2)SH on Au and demonstrate the patterning of biomolecular features with dimensions approaching 40 nm. The electron dose dependence of feature size and pattern resolution is studied in detail by atomic force microscopy (AFM), which reveals two distinct patterning mechanisms. At low doses, the pattern formation occurs by SAM ablation in a self-developing process where the feature size is directly dose-dependent. At higher doses, electron beam-induced deposition of material, so-called contamination writing, is seen in the ablated areas of the SAM. The balance between these two mechanisms is shown to depend on the geometry of the pattern. The patterned SAMs were backfilled with fluorescent 40-nm spheres coated with NeutrAvidin. These protein-coated spheres adhered to exposed areas in the SAM with high selectivity. This direct writing approach for patterning bioactive surfaces is a fast and efficient way to produce patterns with a resolution approaching that of single proteins.

Amyloid-beta aggregates formed at polar-nonpolar interfaces differ from amyloid-beta protofibrils produced in aqueous buffers.

The deposition of aggregated amyloid-beta (Abeta) peptides in the brain as senile plaques is a pathological hallmark of Alzheimer's disease (AD). Several lines of evidence indicate that fibrillar and, in particular, soluble aggregates of these 40- and 42-residue peptides are important in the etiology of AD. Recent studies also stress that amyloid aggregates are polymorphic and that a single polypeptide can fold into multiple amyloid conformations. Here we review our recent reports that Abeta(1-40) in vitro can form soluble aggregates with predominant beta-structures that differ in stability and morphology. One class of aggregates involved soluble Abeta protofibrils, prepared by vigorous overnight agitation of monomeric Abeta(1-40) in low ionic strength buffers. These aggregates were quite stable and disaggregated to only a limited extent on dilution. A second class of soluble Abeta aggregates was generated at polar-nonpolar interfaces. Aggregation in a two-phase system of buffer over chloroform occurred more rapidly than in buffer alone. In buffered 2% hexafluoroisopropanol (HFIP), microdroplets of HFIP were formed and the half-time for aggregation was less than 10 minutes. Like Abeta protofibrils, these interfacial aggregates showed increased thioflavin T fluorescence and were rich in beta-structure by circular dichroism. However, electron microscopy and atomic force microscopy revealed very different morphologies. The HFIP aggregates formed initial globular clusters that progressed over several days to soluble fibrous aggregates. When diluted out of HFIP these aggregates initially were very unstable and disaggregated completely within 2 minutes. However, their stability increased as they progressed to fibers. It is important to determine whether similar interfacial Abeta aggregates are produced in vivo.

Poly(ethylene glycol) self-assembled monolayer island growth.

Here, we report a study of the morphology and growth dynamics of a self-assembled monolayer (SAM) of the amide containing poly(ethylene glycol) (PEG) thiol (CH3O(CH2CH2O)17NHCO(CH2)2SH) on atomically flat Au(111) surfaces. SAM growth from a 20 muM ethanolic solution reveals island growth through three distinct steps: island nucleation, island growth, and coalescence. The coalescence-step, filling voids in the SAM, is by far slowest. The fine structure study reveals dendritic island formation, an observation which can be explained by attractive intermolecular interactions and surface diffusion-limited aggregation. We have also observed a change in the island height, which peaks during the island growth phase. This height change can be associated with a molecular conformational transition.

Rapid assembly of amyloid-beta peptide at a liquid/liquid interface produces unstable beta-sheet fibers.

Accumulation of aggregated amyloid-beta peptide (Abeta) in the brain is a pathological hallmark of Alzheimer's disease (AD). In vitro studies indicate that the 40- to 42-residue Abeta peptide in solution will undergo self-assembly leading to the transient appearance of soluble protofibrils and ultimately to insoluble fibrils. The Abeta peptide is amphiphilic and accumulates preferentially at a hydrophilic/hydrophobic interface. Solid surfaces and air-water interfaces have been shown previously to promote Abeta aggregation, but detailed characterization of these aggregates has not been presented. In this study Abeta(1-40) introduced to aqueous buffer in a two-phase system with chloroform aggregated 1-2 orders of magnitude more rapidly than Abeta in the buffer alone. The interface-induced aggregates were released into the aqueous phase and persisted for 24-72 h before settling as a visible precipitate at the interface. Thioflavin T fluorescence and circular dichroism analyses confirmed that the Abeta aggregates had a beta-sheet secondary structure. However, these aggregates were far less stable than Abeta(1-40) protofibrils prepared in buffer alone and disaggregated completely within 3 min on dilution. Atomic force microscopy revealed that the aggregates consisted of small globules 4-5 nm in height and long flexible fibers composed of these globules aligned roughly along a longitudinal axis, a morphology distinct from that of Abeta protofibrils prepared in buffer alone. The relative instability of the fibers was supported by fiber interruptions apparently introduced by brief washing of the AFM grids. To our knowledge, unstable aggregates of Abeta with beta-sheet structure and fibrous morphology have not been reported previously. Our results provide the clearest evidence yet that the intrinsic beta-sheet structure of an in vitro Abeta aggregate depends on the aggregation conditions and is reflected in the stability of the aggregate and the morphology observed by atomic force microscopy. Resolution of these structural differences at the molecular level may provide important clues to the further understanding of amyloid formation in vivo.

Evidence for a highly elastic shell-core organization of cochlear outer hair cells by local membrane indentation.

Cochlear outer hair cells (OHCs) are thought to play an essential role in the high sensitivity and sharp frequency selectivity of the hearing organ by generating forces that amplify the vibrations of this organ at frequencies up to several tens of kHz. This tuning process depends on the mechanical properties of the cochlear partition, which OHC activity has been proposed to modulate on a cycle-by-cycle basis. OHCs have a specialized shell-core ultrastructure believed to be important for the mechanics of these cells and for their unique electromotility properties. Here we use atomic force microscopy to investigate the mechanical properties of isolated living OHCs and to show that indentation mechanics of their membrane is consistent with a shell-core organization. Indentations of OHCs are also found to be highly nonhysteretic at deformation rates of more than 40 microm/s, which suggests the OHC lateral wall is a highly elastic structure, with little viscous dissipation, as would appear to be required in view of the very rapid changes in shape and mechanics OHCs are believed to undergo in vivo.