Peptide design by artificial neural networks and computer-based evolutionary search

A technique for systematic peptide variation by a combination of rational and evolutionary approaches is presented. The design scheme consists of five consecutive steps: (i) identification of a “seed peptide” with a desired activity, (ii) generation of variants selected from a physicochemical space around the seed peptide, (iii) synthesis and testing of this biased library, (iv) modeling of a quantitative sequence-activity relationship by an artificial neural network, and (v) de novo design by a computer-based evolutionary search in sequence space using the trained neural network as the fitness function. This strategy was successfully applied to the identification of novel peptides that fully prevent the positive chronotropic effect of anti-β1-adrenoreceptor autoantibodies from the serum of patients with dilated cardiomyopathy. The seed peptide, comprising 10 residues, was derived by epitope mapping from an extracellular loop of human β1-adrenoreceptor. A set of 90 peptides was synthesized and tested to provide training data for neural network development. De novo design revealed peptides with desired activities that do not match the seed peptide sequence. These results demonstrate that computer-based evolutionary searches can generate novel peptides with substantial biological activity.

A hierarchical neuronal network for planning behavior

Planning a goal-directed sequence of behavior is a higher function of the human brain that relies on the integrity of prefrontal cortical areas. In the Tower of London test, a puzzle in which beads sliding on pegs must be moved to match a designated goal configuration, patients with lesioned prefrontal cortex show deficits in planning a goal-directed sequence of moves. We propose a neuronal network model of sequence planning that passes this test and, when lesioned, fails in a way that mimics prefrontal patients’ behavior. Our model comprises a descending planning system with hierarchically organized plan, operation, and gesture levels, and an ascending evaluative system that analyzes the problem and computes internal reward signals that index the correct/erroneous status of the plan. Multiple parallel pathways connecting the evaluative and planning systems amend the plan and adapt it to the current problem. The model illustrates how specialized hierarchically organized neuronal assemblies may collectively emulate central executive or supervisory functions of the human brain.

A model of multiplicative neural responses in parietal cortex.

Visual responses of neurons in parietal area 7a are modulated by a combined eye and head position signal in a multiplicative manner. Neurons with multiplicative responses can act as powerful computational elements in neural networks. In the case of parietal cortex, multiplicative gain modulation appears to play a crucial role in the transformation of object locations from retinal to body-centered coordinates. It has proven difficult to uncover single-neuron mechanisms that account for neuronal multiplication. Here we show that multiplicative responses can arise in a network model through population effects. Specifically, neurons in a recurrently connected network with excitatory connections between similarly tuned neurons and inhibitory connections between differently tuned neurons can perform a product operation on additive synaptic inputs. The results suggest that parietal responses may be based on this architecture.

Neural-network-based classification of cognitively normal, demented, Alzheimer disease and vascular dementia from single photon emission with computed tomography image data from brain.

Single photon emission with computed tomography (SPECT) hexamethylphenylethyleneamineoxime technetium-99 images were analyzed by an optimal interpolative neural network (OINN) algorithm to determine whether the network could discriminate among clinically diagnosed groups of elderly normal, Alzheimer disease (AD), and vascular dementia (VD) subjects. After initial image preprocessing and registration, image features were obtained that were representative of the mean regional tissue uptake. These features were extracted from a given image by averaging the intensities over various regions defined by suitable masks. After training, the network classified independent trials of patients whose clinical diagnoses conformed to published criteria for probable AD or probable/possible VD. For the SPECT data used in the current tests, the OINN agreement was 80 and 86% for probable AD and probable/possible VD, respectively. These results suggest that artificial neural network methods offer potential in diagnoses from brain images and possibly in other areas of scientific research where complex patterns of data may have scientifically meaningful groupings that are not easily identifiable by the researcher.

Emergent spindle oscillations and intermittent burst firing in a thalamic model: specific neuronal mechanisms.

The rhythmogenesis of 10-Hz sleep spindles is studied in a large-scale thalamic network model with two cell populations: the excitatory thalamocortical (TC) relay neurons and the inhibitory nucleus reticularis thalami (RE) neurons. Spindle-like bursting oscillations emerge naturally from reciprocal interactions between TC and RE neurons. We find that the network oscillations can be synchronized coherently, even though the RE-TC connections are random and sparse, and even though individual neurons fire rebound bursts intermittently in time. When the fast gamma-aminobutyrate type A synaptic inhibition is blocked, synchronous slow oscillations resembling absence seizures are observed. Near-maximal network synchrony is established with even modest convergence in the RE-to-TC projection (as few as 5-10 RE inputs per TC cell suffice). The hyperpolarization-activated cation current (Ih) is found to provide a cellular basis for the intermittency of rebound bursting that is commonly observed in TC neurons during spindles. Such synchronous oscillations with intermittency can be maintained only with a significant degree of convergence for the TC-to-RE projection.

Forming classes by stimulus frequency: Behavior and theory

Visual classification is the way we relate to different images in our environment as if they were the same, while relating differently to other collections of stimuli (e.g., human vs. animal faces). It is still not clear, however, how the brain forms such classes, especially when introduced with new or changing environments. To isolate a perception-based mechanism underlying class representation, we studied unsupervised classification of an incoming stream of simple images. Classification patterns were clearly affected by stimulus frequency distribution, although subjects were unaware of this distribution. There was a common bias to locate class centers near the most frequent stimuli and their boundaries near the least frequent stimuli. Responses were also faster for more frequent stimuli. Using a minimal, biologically based neural-network model, we demonstrate that a simple, self-organizing representation mechanism based on overlapping tuning curves and slow Hebbian learning suffices to ensure classification. Combined behavioral and theoretical results predict large tuning overlap, implicating posterior infero-temporal cortex as a possible site of classification.

Theory of orientation tuning in visual cortex.

The role of intrinsic cortical connections in processing sensory input and in generating behavioral output is poorly understood. We have examined this issue in the context of the tuning of neuronal responses in cortex to the orientation of a visual stimulus. We analytically study a simple network model that incorporates both orientation-selective input from the lateral geniculate nucleus and orientation-specific cortical interactions. Depending on the model parameters, the network exhibits orientation selectivity that originates from within the cortex, by a symmetry-breaking mechanism. In this case, the width of the orientation tuning can be sharp even if the lateral geniculate nucleus inputs are only weakly anisotropic. By using our model, several experimental consequences of this cortical mechanism of orientation tuning are derived. The tuning width is relatively independent of the contrast and angular anisotropy of the visual stimulus. The transient population response to changing of the stimulus orientation exhibits a slow "virtual rotation." Neuronal cross-correlations exhibit long time tails, the sign of which depends on the preferred orientations of the cells and the stimulus orientation.

Transgenic knockout mice exclusively expressing human hemoglobin S after transfer of a 240-kb βs-globin yeast artificial chromosome: A mouse model of sickle cell anemia

Sickle cell anemia (SCA) and thalassemia are among the most common genetic diseases worldwide. Current approaches to the development of murine models of SCA involve the elimination of functional murine α- and β-globin genes and substitution with human α and βs transgenes. Recently, two groups have produced mice that exclusively express human HbS. The transgenic lines used in these studies were produced by coinjection of human α-, γ-, and β-globin constructs. Thus, all of the transgenes are integrated at a single chromosomal site. Studies in transgenic mice have demonstrated that the normal gene order and spatial organization of the members of the human β-globin gene family are required for appropriate developmental and stage-restricted expression of the genes. As the cis-acting sequences that participate in activation and silencing of the γ- and β-globin genes are not fully defined, murine models that preserve the normal structure of the locus are likely to have significant advantages for validating future therapies for SCA. To produce a model of SCA that recapitulates not only the phenotype, but also the genotype of patients with SCA, we have generated mice that exclusively express HbS after transfer of a 240-kb βs yeast artificial chromosome. These mice have hemolytic anemia, 10% irreversibly sickled cells in their peripheral blood, reticulocytosis, and other phenotypic features of SCA.