Search PubMed⌕ Search

Biomedical subjects

S A Kauffman

Publications and source records attributed to S A Kauffman.

At least 19 recordsLinked to original sources

Network growth models and genetic regulatory networks.

We study a class of growth algorithms for directed graphs that are candidate models for the evolution of genetic regulatory networks. The algorithms involve partial duplication of nodes and their links, together with the innovation of new links, allowing for the possibility that input and output links from a newly created node may have different probabilities of survival. We find some counterintuitive trends as the parameters are varied, including the broadening of the in-degree distribution when the probability for retaining input links is decreased. We also find that both the scaling of transcription factors with genome size and the measured degree distributions for genes in yeast can be reproduced by the growth algorithm if and only if a special seed is used to initiate the process.

Animals↗

Scaling in ordered and critical random boolean networks.

Random Boolean networks, originally invented as models of genetic regulatory networks, are simple models for a broad class of complex systems that show rich dynamical structures. From a biological perspective, the most interesting networks lie at or near a critical point in parameter space that divides "ordered" from "chaotic" attractor dynamics. We study the scaling of the average number of dynamically relevant nodes and the median number of distinct attractors in such networks. Our calculations indicate that the correct asymptotic scalings emerge only for very large systems.

Gene Expression↗

Prolegomenon to a general biology.

This chapter is a slightly modified version of Chapter 1 of my book, Investigations. I discuss "autonomous agents," the origin of life, autocatalytic sets of polymers, work cycles, puzzles-about evolutionary theory, and the birth of a general biology dealing with biospheres anywhere in the universe.

Behavior↗

Selection dynamics in autocatalytic systems: templates replicating through binary ligation.

The theory of autocatalytic binary ligation is reviewed within the context of a consistently applied Michaelis-Menten quasi-steady-state approximation to obtain explicit analytical results describing time-course data from experiments. A detailed protocol for the step-wise elucidation of a minimal set of experimental parameters is outlined. The kinetic equations are then generalized to cases of self- and cross-catalysis among an arbitrary number of different templates and applied to experiments involving just two templates. Depending on the values of various kinetic parameters such systems can display exclusionary Darwinian selection corresponding to an exponential growth law, selective coexistence or coexistence of all species characteristic of a parabolic growth law; the intermediate behaviour arises as a property of the full mechanism analysed here. Our results are applicable to the classical case of self-replicating nucleic acids and their analogues as well as to newly discovered self-replicating peptides.

Catalysis↗

Use of hearing protection and perceptions of noise exposure and hearing loss among construction workers.

The purpose of this study was to describe construction workers' use of hearing protection devices (HPDs) and determine their perceptions of noise exposure and hearing loss. Operating engineers, carpenters, and plumbers/pipe fitters in the Midwest (n = 400) completed a written questionnaire regarding their use of HPDs and their perceptions of noise exposure and hearing loss. Subjects were recruited through their trade union groups. Mean reported use of HPDs and mean perceived noise exposure were compared across trade groups. Bivariate and multivariate analysis techniques were used to assess relationships between use of HPDs and trade category, education, age, years of employment, noise exposure, and hearing loss. Bivariate analyses identified significant differences in mean use of HPDs by age, years of employment, and trade group. Multivariate logistic regression assessing the independent effects of these variables found significant differences only by trade group. Results indicate a need for significant improvement in all three trade groups' use of HPDs, and suggest a need to consider use and exposure levels, demographics, and trade group membership in designing hearing conservation programs.

Adult↗

Criticality and parallelism in combinatorial optimization.

Local search methods constitute one of the most successful approaches to solving large-scale combinatorial optimization problems. As these methods are increasingly parallelized, optimization performance initially improves but then abruptly degrades to no better than that of random search beyond a certain point. The existence of this transition is demonstrated for a family of generalized spin-glass models and the traveling salesman problem. Finite-size scaling is used to characterize size-dependent effects near the transition, and analytical insight is obtained through a mean-field approximation.

Algorithms↗

Search strategies for applied molecular evolution.

A new approach to drug discovery is based on the generation of high diversity libraries of DNA, RNA, peptides or small molecules. Search of such libraries for useful molecules is an optimization problem on high-dimensional molecular fitness landscapes. We utilize a spin-glass-like model, the NK model, to analyze search strategies based on pooling, mutation, recombination and selective hill-climbing. Our results suggest that pooling followed by recombination and/or hill-climbing finds better candidate molecules than pooling alone on most molecular landscapes. Our results point to new experiments to assess the structure of molecular fitness landscapes and improve current models.

Drug Design↗

Libraries of random-sequence polypeptides produced with high yield as carboxy-terminal fusions with ubiquitin.

Libraries of random-sequence polypeptides have been shown to be valuable sources of novel molecules possessing a variety of useful biologic-like activities, some of which may hold promise as potential vaccines and therapeutics. Previous random peptide expression systems were limited to low levels of peptide production and often to short sequences. Here we describe a series of libraries designed for increased polypeptide length. Cloned as carboxy-terminal extensions of ubiquitin, the fusions were produced in E. coli at high levels, and were purified to homogeneity. The majority of the extension proteins examined could be cleaved from ubiquitin by treatment with a ubiquitin-fusion hydrolase. The libraries described here are appropriate sources of novel polypeptides with desired binding or catalytic function, as well as tools with which to examine inherent properties of proteins as a whole. Toward the latter goal, we have examined structural properties of random-sequence proteins purified from these libraries. Quite surprisingly, fluorescence emission spectra of intrinsic tryptophan residues in several purified fusion proteins, under native-like and denaturing conditions, often resemble those expected for folded and unfolded states, respectively. The results presented here detail an important expansion in the range of potential uses for random-sequence polypeptide libraries.

Amino Acids↗

Design of synthetic gene libraries encoding random sequence proteins with desired ensemble characteristics.

Libraries of random sequence polypeptides are useful as sources of unevolved proteins, novel ligands, and potential lead compounds for the development of vaccines and therapeutics. The expression of small random peptides has been achieved previously using DNA synthesized with equimolar mixtures of nucleotides. For many potential uses of random polypeptide libraries, concerns such as avoiding termination codons and matching target amino acid compositions make more complex designs necessary. In this study, three mixtures of nucleotides, corresponding to the three positions in the codon, were designed such that semirandom DNA synthesized by repeated cycles of the three mixtures created an open reading frame encoding random sequence polypeptides with desired ensemble characteristics. Two methods were used to design the nucleotide mixtures: the manual use of a spreadsheet and a refining grid search algorithm. Using design targets of less than or equal to 1% stop codons and an amino acid composition based on the average ratios observed in natural, globular proteins, the search methods yielded similar nucleotide ratios, Semirandom DNA, synthesized with a designed, three-residue repeat pattern, can encode libraries of very high diversity and represents an important tool for the construction of random polypeptide libraries.

Amino Acid Sequence↗

Coevolution to the edge of chaos: coupled fitness landscapes, poised states, and coevolutionary avalanches.

We introduce a broadened framework to study aspects of coevolution based on the NK class of statistical models of rugged fitness landscapes. In these models the fitness contribution of each of N genes in a genotype depends epistatically on K other genes. Increasing epistatic interactions increases the rugged multipeaked character of the fitness landscape. Coevolution is thought of, at the lowest level, as a coupling of landscapes such that adaptive moves by one player deform the landscapes of its immediate partners. In these models we are able to tune the ruggedness of landscapes, how richly intercoupled any two landscapes are, and how many other players interact with each player. All these properties profoundly alter the character of the coevolutionary dynamics. In particular, these parameters govern how readily coevolving ecosystems achieve Nash equilibria, how stable to perturbations such equilibria are, and the sustained mean fitness of coevolving partners. In turn, this raises the possibility that an evolutionary metadynamics due to natural selection may sculpt landscapes and their couplings to achieve coevolutionary systems able to coadapt well. The results suggest that sustained fitness is optimized when landscape ruggedness relative to couplings between landscapes is tuned such that Nash equilibria just tenuously form across the ecosystem. In this poised state, coevolutionary avalanches appear to propagate on all length scales in a power law distribution. Such avalanches may be related to the distribution of small and large extinction events in the record.

Biological Evolution↗

Drosophila segmentation: supercomputer simulation of prepattern hierarchy.

Spontaneous prepattern formation in a two level hierarchy of reaction-diffusion systems is simulated in three space co-ordinates and time, mimicking gap gene and primary pair-rule gene expression. The model rests on the idea of Turing systems of the second kind, in which one prepattern generates position dependent rate constants for a subsequent reaction-diffusion system. Maternal genes are assumed responsible for setting up gradients from the anterior and posterior ends, one of which is needed to stabilize a double period prepattern suggested to underly the read out of the gap genes. The resulting double period pattern in turn stabilizes the next prepattern in the hierarchy, which has a short wavelength with many characteristics of the stripes seen in actual primary pair-rule gene expression. Without such hierarchical stabilization, reaction-diffusion mechanisms yield highly patchy short wave length patterns, and thus unreliable stripes. The model yields seven stable stripes located in the middle of the embryo, with the potential for additional expression near the poles, as observed experimentally. The model does not rely on specific chemical reaction kinetics, rather the effect is general to many such kinetic schemes. This makes it robust to parameter changes, and it has good potential for adapting to size and shape changes as well. The study thus suggests that the crucial organizing principle in early Drosophila embryogenesis is based on global field mechanisms, not on particular local interactions.

Animals↗

Spatial harmonics and pattern specification in early Drosophila development. Part I. Bifurcation sequences and gene expression.

Molecular probes have now provided an unprecedented wealth of detail revealing the changing spatial patterns of gene products in early Drosophila development. This is examined for dynamic properties which might provide insights into the underlying behaviour of the patterning process. What emerges is that transcripts and protein products of members of the major categories of zygotically active genes involved in segmentation pass through transient spatial patterns that are suggestive of harmonic sequences arising from spatial frequency-doubling bifurcations. That is to say, these patterns are typically periodic in space and show a doubling in the number of domains of spatial expression as development proceeds. One of these patterns reflects the primary functional role of the gene in the establishment of the spatial pattern. The different categories of segmentation gene pass through these transients at different rates, those with the longest functional wavelength progressing most slowly. Each gene in a category has its own unique phase relationship to other members, as well as particular variations on the harmonic sequence theme. The result is that the developing embryo experiences a spatial hierarchy of phase-shifted patterning influences that span the range from the whole embryo to single segments, providing progressively more spatial resolution in the patterning process. The characteristic transients and the dynamic relationships between genes of the different categories suggest that gene products expressed in longer-wavelength patterns act as bifurcation parameters on the dynamic system generating the next shorter wavelength category. Such parametric influences are known to result in frequency-doubling bifurcations in Turing reaction-diffusion systems. A general model is proposed of a hierarchically-nested set of quasi-autonomous dynamic systems involving gene activities that can generate the progressively finer spatial order that emerges during embryogenesis. This model has implications for the general stability properties of evolving epigenetic systems.

Animals↗

Spatial harmonics and pattern specification in early Drosophila development. Part II. The four colour wheels model.

We review the evidence presented in Part I showing that transcripts and protein products of maternal, gap, pair-rule, and segment polarity genes exhibit increasingly complex, multipeaked longitudinal waveforms in the early Drosophila embryo. The central problem we address in Part II is the use the embryo makes of these wave forms to specify longitudinal pattern. Based on the fact that mutants of many of these genes generate deletions and mirror symmetrical duplications of pattern elements on length scales ranging from about half the egg to within segments, we propose that position is specified by measuring a "phase angle" by use of the ratios of two or more variables. Pictorially, such a phase angle can be thought of as a colour on a colour wheel. Any such model contains a phaseless singularity where all or many phases, or colours, come together. We suppose as well that positional values sufficiently close to the singularity are meaningless, hence a "dead zone". Duplications and deletions are accounted for by deformation of the cycle of morphogen values occurring along the antero-posterior axis. If the cycle of values surrounds the singularity and lies outside the dead zone, pattern is normal. If the curve transects the dead zone, pattern elements are deleted. If the curve lies entirely on one side of the singularity, pattern elements are deleted and others are duplicated with mirror symmetry. The existence of different wavelength transcript patterns in maternal, gap, pair-rule, and segment polarity genes and the roles of those same genes in generating deletions and mirror symmetrical duplications on a variety of length scales lead us to propose that position is measured simultaneously on at least four colour wheels, which cycle different numbers of times along the anterior-posterior axis. These yield progressively finer grained positional information. Normal pattern specification requires a unique angle, outside of the dead zone, from each of the four wheels. Deformations of the cycle of gene product concentrations yield the deletions and mirror symmetric duplications observed in the mutants discussed. The alternative familiar hypothesis that longitudinal position is specified in an "on" "off" combinatorial code does not readily account for the duplication deletion phenomena.

Animals↗

The NK model of rugged fitness landscapes and its application to maturation of the immune response.

Adaptive evolution is, to a large extent, a complex combinatorial optimization process. Such processes can be characterized as "uphill walks on rugged fitness landscapes". Concrete examples of fitness landscapes include the distribution of any specific functional property such as the capacity to catalyze a specific reaction, or bind a specific ligand, in "protein space". In particular, the property might be the affinity of all possible antibody molecules for a specific antigenic determinant. That affinity landscape presumably plays a critical role in maturation of the immune response. In this process, hypermutation and clonal selection act to select antibody V region mutant variants with successively higher affinity for the immunizing antigen. The actual statistical structure of affinity landscapes, although knowable, is currently unknown. Here, we analyze a class of mathematical models we call NK models. We show that these models capture significant features of the maturation of the immune response, which is currently thought to share features with general protein evolution. The NK models have the important property that, as the parameter K increases, the "ruggedness" of the NK landscape varies from a single peaked "Fujiyama" landscape to a multi-peaked "badlands" landscape. Walks to local optima on such landscapes become shorter as K increases. This fact allows us to choose a value of K that corresponds to the experimentally observed number of mutational "steps", 6-8, taken as an antibody sequence matures. If the mature antibody is taken to correspond to a local optimum in the model, tuning the model requires that K be about 40, implying that the functional contribution of each amino acid in the V region is affected by about 40 others. Given this value of K, the model then predicts several features of "antibody space" that are in qualitative agreement with experiment: (1) The fraction of fitter variants of an initial "roughed in" germ line antibody amplified by clonal selection is about 1-2%. (2) Mutations at some sites of the mature antibody hardly affect antibody function at all, but mutations at other sites dramatically decrease function. (3) The same "roughed in" antibody sequence can "walk" to many mature antibody sequences. (4) Many adaptive walks can end on the same local optimum. (5) Comparison of different mature sequences derived from the same initial V region shows evolutionary hot spots and parallel mutations. All these predictions are open to detailed testing by obtaining monoclonal antibodies early in the immune response and carrying out in vitro mutagenesis and adaptive hill climbing with respect to affinity for the immunizing antigen.

Amino Acid Sequence↗

Modeling adaptive biological systems.

During the evolution of many systems found in nature, both the system composition and the interactions between components will vary. Equating the dimension with the number of different components, a system which adds or deletes components belongs to a class of dynamical systems with a finite dimensional phase space of variable dimension. We present two models of biochemical systems with a variable phase space, a model of autocatalytic reaction networks in the prebiotic soup and a model of the idiotypic network of the immune system. Each model contains characteristic meta-dynamical rules for constructing equations of motion from component properties. The simulation of each model occurs on two levels. On one level, the equations of motion are integrated to determine the state of each component. On a second level, algorithms which approximate physical processes in the real system are employed to change the equations of motion. Models with meta-dynamical rules possess several advantages for the study of evolving systems. First, there are no explicit fitness functions to determine how the components of the model rank in terms of survivability. The success of any component is a function of its relationship to the rest of the system. A second advantage is that since the phase space representation of the system is always finite but continually changing, we can explore a potentially infinite phase space which would otherwise be inaccessible with finite computer resources. Third, the enlarged capacity of systems with meta-dynamics for variation allows us to conduct true evolution experiments. The modeling methods presented here can be applied to many real biological systems. In the two studies we present, we are investigating two apparent properties of adaptive networks. With the simulation of the prebiotic soup, we are most interested in how a chemical reaction network might emerge from an initial state of relative disorder. With the study of the immune system, we study the self-regulation of the network including its ability to distinguish between species which are part of the network and those which are not.

Adaptation, Biological↗