PubMed HealthSearch

PubMed · 8866051

Self-organization in computer simulated selective systems.

Abstract

In certain biologic systems, the signal selects functional or numerical expansion of the recognizing elements. Examples of these systems include the immune system, brain cortex, and evolution. One common feature of these Darwinian-type systems is degenerate recognition, in which one signal can recognize several different elements, with different affinities and consequences. For example, T cell antigen receptors and antibodies demonstrate relative but not absolute specificity of recognition. Thus, the variables of dose of the signal and affinity of the recognizing element modulate the outcome. Another feature of these systems is the ability to create self-organized patterns, which do not mirror the incoming signals. The hypothesis of this study is that degenerate recognition with subsequent selection of recognizing elements can explain self-organization of these systems. An entirely numerical model was explored, using the cellular automata approach. Three intrinsic features of a common selective system were incorporated into this model: a large number of recognizing elements; degenerative recognition of stimuli by these elements; and subsequent selection. Different numerical patterns of incoming stimuli were tested. The model showed self-organizing dynamics. Usually, the population of recognizing elements demonstrated an initial period of equilibrium, then a chaotic transitional state, and, finally, the bifurcational appearance of a stable self-organized pattern. The final resolution into a stable pattern can be either gradual or quasi-saltational. I conclude that systems with a large number of recognizing elements, degenerative recognition, and selection of recognizing elements can self organize based upon the pattern of the incoming stimuli.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

S P Atamas. 1996. Self-organization in computer simulated selective systems.. https://doi.org/10.1016/0303-2647(96)01612-7

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Generating correlated data for omics simulation.

Simulation of realistic omics data is a key input for benchmarking studies that help users obtain optimal computational pipelines. Omics data involves large numbers of measured features on each sample and these measures are generally correlated with each other. However, simulation too often ignores these correlations, perhaps due to computational and statistical hurdles of doing so. To alleviate this, we describe three approaches for generating omics-scale data with correlated measures which mimic real datasets. These approaches are all based on a Gaussian copula approach with a covariance matrix that decomposes into a diagonal part and a low-rank part. This decomposition allows for extremely efficient simulation, overcoming a hurdle for adoption of past methods. We use these approaches to demonstrate the importance of including correlation in two benchmarking applications. First, we show that variance of results from the popular DESeq2 method increases when dependence is included. Second, we demonstrate that CYCLOPS, a method for inferring circadian time of collection from transcriptomics, improves in performance when given gene-gene dependencies in some circumstances. We provide an R package, dependentsimr, that has efficient implementations of these methods and can generate dependent data with arbitrary marginal distributions, including discrete (binary, ordered categorical, Poisson, negative binomial), continuous (normal), or with an empirical distribution.

Computer Simulation

Addressing current challenges in cancer immunotherapy with mathematical and computational modelling.

The goal of cancer immunotherapy is to boost a patient's immune response to a tumour. Yet, the design of an effective immunotherapy is complicated by various factors, including a potentially immunosuppressive tumour microenvironment, immune-modulating effects of conventional treatments and therapy-related toxicities. These complexities can be incorporated into mathematical and computational models of cancer immunotherapy that can then be used to aid in rational therapy design. In this review, we survey modelling approaches under the umbrella of the major challenges facing immunotherapy development, which encompass tumour classification, optimal treatment scheduling and combination therapy design. Although overlapping, each challenge has presented unique opportunities for modellers to make contributions using analytical and numerical analysis of model outcomes, as well as optimization algorithms. We discuss several examples of models that have grown in complexity as more biological information has become available, showcasing how model development is a dynamic process interlinked with the rapid advances in tumour-immune biology. We conclude the review with recommendations for modellers both with respect to methodology and biological direction that might help keep modellers at the forefront of cancer immunotherapy development.

Computer Simulation

Effective energy function for proteins in solution.

A Gaussian solvent-exclusion model for the solvation free energy is developed. It is based on theoretical considerations and parametrized with experimental data. When combined with the CHARMM 19 polar hydrogen energy function, it provides an effective energy function (EEF1) for proteins in solution. The solvation model assumes that the solvation free energy of a protein molecule is a sum of group contributions, which are determined from values for small model compounds. For charged groups, the self-energy contribution is accounted for primarily by the exclusion model. Ionic side-chains are neutralized, and a distance-dependent dielectric constant is used to approximate the charge-charge interactions in solution. The resulting EEF1 is subjected to a number of tests. Molecular dynamics simulations at room temperature of several proteins in their native conformation are performed, and stable trajectories are obtained. The deviations from the experimental structures are similar to those observed in explicit water simulations. The calculated enthalpy of unfolding of a polyalanine helix is found to be in good agreement with experimental data. Results reported elsewhere show that EEF1 clearly distinguishes correctly from incorrectly folded proteins, both in static energy evaluations and in molecular dynamics simulations and that unfolding pathways obtained by high-temperature molecular dynamics simulations agree with those obtained by explicit water simulations. Thus, this energy function appears to provide a realistic first approximation to the effective energy hypersurface of proteins.

Computer Simulation