Waddington redux: models and explanation in stem cell and systems biology

Stem cell biology and systems biology are two prominent new approaches to studying cell development. In stem cell biology, the predominant method is experimental manipulation of concrete cells and tissues. Systems biology, in contrast, emphasizes mathematical modeling of cellular systems. For scientists and philosophers interested in development, an important question arises: how should the two approaches relate? This essay proposes an answer, using the model of Waddington’s landscape to triangulate between stem cell and systems approaches. This simple abstract model represents development as an undulating surface of hills and valleys. Originally constructed by C. H. Waddington to visually explicate an integrated theory of genetics, development and evolution, the landscape model can play an updated unificatory role. I examine this model’s structure, representational assumptions, and uses in all three contexts, and argue that explanations of cell development require both mathematical models and concrete experiments. On this view, the two approaches are interdependent, with mathematical models playing a crucial but circumscribed role in explanations of cell development.     

Thiophosphorylated proteins in chromaffin cells are chromaffin vesicle matrix proteins

Bovine adrenal chromaffin cells were incubated with inorganic thiophosphate, using a protocol similar to experiments with inorganic phosphate, in order to determine the source of previously observed thiophosphoproteins. Incubation of cultured cells with [³⁵S]thiophosphate resulted in its incorporation into cell constituents within 2 min. SDS PAGE of the treated cells showed incorporation of label into a broad 97–121 kDa band that was evident after 5 min of treatment and increased progressively to the 40 min exposure limit. Monolayers of chronically treated cells were fractionated into subcellular constituents. The only particulate fraction containing radiolabelled proteins was the chromaffin vesicle fraction. Two-dimensional electrophoresis of the treated cells and isolated chromaffin vesicles showed a majority of proteins in the acidic region of the first dimension gel. A fluorogram of the gel revealed two regions of radiolabelled proteins at acidic and neutral regions of the 2-D gel. These were within the boundaries of the 97–121 kDa band. The thiophosphorylated proteins were released as soluble proteins upon osmotic or freeze-thaw lysis of the vesicles. Chromaffin vesicles isolated from either cultured cells or adrenal medulla tissue were energized by 2 mM ATP but not by the analog adenosine 5′-O-(3-thiotriphosphate). The 97–121 kDa proteins in intact or lysed vesicles prepared from adrenal medulla tissue were not thiophosphorylated by either inorganic thiophosphate or adenosine 5′-O-(3-thiotriphosphate) in the presence or absence of energization by ATP. Nearly complete loss of radiolabel from matrix proteins treated with chondroitinase ABC suggests that it is a component of vesicle proteoglycans.

The results demonstrate that chromaffin vesicle matrix proteins are rapidly and intensely thiophosphorylated in cultured chromaffin cells but not in isolated vesicles. The data suggest that phosphorylation must play an important role in the normal function of these vesicle proteins.     

Anisotropy decay associated fluorescence spectra and analysis of rotational heterogeneity. 2. 1,6-Diphenyl-1,3,5-hexatriene in lipid bilayers

The application of a new spectroscopic tool [Knutson, J. R., Davenport, L., & Brand, L. (1986) Biochemistry (preceding paper in this issue)] for studying rotational microheterogeneity of probe location in lipid bilayer systems is described. Anisotropy decay associated spectra are derived from experimentally obtained polarized emission components. "Early" difference spectra (IV - IH) contain contributions from both fast and slow rotors, while "late" difference spectra predominantly reflect the emission from slowly rotating fluorophores. Anisotropy decay associated spectra have been used to resolve the emission spectra of 1,6-diphenyl-1,3,5-hexatriene (DPH) imbedded within a known rotationally heterogeneous mixture of two vesicle types (L-alpha-dimyristoyllecithin and L-alpha-dipalmitoyllecithin). At 29 degrees C, diphenylhexatriene within pure dimyristoyllecithin vesicles rotates rapidly, with a small r infinity, while diphenylhexatriene in dipalmitoyllecithin vesicles exhibits a large r infinity. Spectra for diphenylhexatriene imbedded in the two vesicle types show small but significant spectral differences. A spectrum of a mixture of the two vesicle types with DPH lies between these characteristic component spectra. The spectrum extracted for "immobilized" probes in the mixture correctly overlays the dipalmitoyllecithin spectrum. Further studies have shown that diphenylhexatriene exhibits more than one emission anisotropy decay associated spectrum in vesicles of a single lipid type, when that lipid is near its phase transition temperature. Diphenylhexatriene apparently inhabits more than one rotational environment even in these "homogeneous" vesicle preparations.     

Sperm aging in the male after sexual rest: Contribution to chromosome anomalies

The chromosome complement was studied in first-cleavage metaphases of mouse zygotes resulting from sperm aged in the male physiologically, after sexual rest. Females were inseminated by control males mating at 3-day intervals while experimentals mated to males that had had a sexual rest of 14 or more days. A total of 1954 eggs were collected 33–35 h post-HCG from 101 superovulated females mated to 42 controls and 43 experimental males. The fertilization rate was similar in both groups, being 84% and 85%, respectively. G-banded or Q-banded chromosomes were analyzed in 301 (68.3%) controls and 392 (49%) experimental first-cleavage metaphases. The overall rate of chromosome anomalies in controls was 4.45% as compared to 10.94% in experimentals, a highly significant difference. In the experimental group compared to controls, the frequency of trisomy, triploidy, structural rearrangements, and tetraploidy increased from 3.9% to 6.9%, 0% to 1.6%, 0.8% to 2.8%, and 0% to 1.3%, respectively. The genomic source of origin of the abnormalities was determined on the basis of differential condensation of the genomes. In the experimentals, grossly unbalanced sperm (diploids, disomics, double disomics, and those with large fragments) fertilized significantly more oocytes compared to controls. Our results implicate an advantage either in numbers or fertilizing capability for chromosomally abnormal sperm in a physiologically aged population.     

Decreased Fc receptor avidity and degradative function of monocytes from patients with systemic lupus erythematosus

We studied the binding and degradation of stable, soluble heat aggregates of 125I-IgG (A-IgG) by monocytes from 30 patients with systemic lupus erythematosus (SLE) and 30 normals. Relative avidities (KE) for Fc receptor (FcR) binding of A-IgG and maximal binding of A-IgG by monocytes were determined from Scatchard plots of binding data obtained at 4 degrees C. Rates of degradation (Vmax) of A-IgG at 37 degrees C were calculated from Lineweaver-Burke plots of the Michaelis-Menton equation. KE were decreased in SLE monocytes (15.5 X 10(-9) L/M) as compared with normals (20.1 X 10(-9) L/M, p less than 0.005) and Vmax were decreased for SLE (0.89 ng/hr) as compared with normals (1.11 ng/hr, p less than 0.005). The maximal FcR binding by SLE monocytes was not statistically different in SLE patients and normals, but monocytes from SLE patients with active disease showed a lower maximal binding capacity for A-IgG (4.9 ng/10(5) cells) than normals (5.4 ng/10(5) cells, p less than 0.05). KE and Vmax in SLE were also lower for patients with active disease than for normal subjects. KE in patients whose anti-ssDNA binding was greater than 20% were lower than for those with DNA binding of less than 20% (p less than 0.005). These data suggest that patients with active SLE have diminished numbers of available FcR on their circulating monocytes, possibly due to interiorization of FcR during endocytosis of endogenous circulating immune complexes.     

Molecular evolution of enzyme structure: Construction of a hybrid hamster/Escherichia coli aspartate transcarbamoylase

Summary Aspartate transcarbamoylase (ATCase, EC 2.1.3.2) is the first unique enzyme common to de novo pyrimidine biosynthesis and is involved in a variety of structural patterns in different organisms. InEscherichia coli, ATCase is a functionally independent, oligomeric enzyme; in hamster, it is part of a trifunctional protein complex, designated CAD, that includes the preceding and subsequent enzymes of the biosynthetic pathway (carbamoyl phosphate synthetase and dihydroorotase). The complete complementary DNA (cDNA) nucleotide sequence of the ATCase-encoding portion of the hamster CAD gene is reported here. A comparison of the deduced amino acid sequences of the hamster andE. coli catalytic peptides revealed an overall 44% amino acid similarity, substantial conservation of predicted secondary structure, and complete conservation of all the amino acids implicated in the active site of theE. coli enzyme. These observations led to the construction of a functional hybrid ATCase formed by intragenic fusion based on the known tertiary structure of the bacterial enzyme. In this fusion, the amino terminal half (the “polar domain”) of the fusion protein was provided by a hamster ATCase cDNA subclone, and the carboxyl terminal portion (the “equatorial domain”) was derived from a clonedpyrBI operon ofE. coli K-12. The recombinant plasmid bearing the hybrid ATCase was shown to satisfy growth requirements of transformedE. coli pyrB ⁻ cells. The functionality of thisE. coli-hamster hybrid enzyme confirms conservation of essential structure-function relationships between evolutionarily distant and structurally divergent ATCases.     

The stratum corneum lipid thermotropic phase behavior.

The stratum corneum, the outermost layer of mammalian skin, is considered the least permeable skin layer to the diffusion of water and other solutes. It is generally accepted that the intercellular lipid multilayer domain is the diffusional pathway for most lipophilic solutes. Fluidization of the lipid multilayers is believed to result in the loss of barrier properties of the stratum corneum. Current investigations address the lipid thermotropic phase behavior in terms of lipid alkyl chain packing, mobility and conformational order as measured by Fourier transform infrared (FTIR) spectroscopy. A solid-solid phase transition is observed with increased alkyl chain mobility followed by a gel to liquid-crystalline phase transition near 65 degrees C. These results further elucidate the role of lipid fluidity that may contribute to the transport properties of the stratum corneum.