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.

Self-assembly of fibronectin into fibrillar networks underneath dipalmitoyl phosphatidylcholine monolayers: Role of lipid matrix and tensile forces

The cell-mediated assembly of fibronectin (Fn) into fibrillar matrices is a complex multistep process that is incompletely understood because of the chemical complexity of the extracellular matrix and a lack of experimental control over molecular interactions and dynamic events. We have identified conditions under which Fn assembles into extended fibrillar networks after adsorption to a dipalmitoyl phosphatidylcholine (DPPC) monolayer in contact with physiological buffer. We propose a sequential model for the Fn assembly pathway, which involves the orientation of Fn underneath the lipid monolayer by insertion into the liquid expanded (LE) phase of DPPC. Attractive interactions between these surface-anchored proteins and the liquid condensed (LC) domains leads to Fn enrichment at domain edges. Spontaneous self-assembly into fibrillar networks, however, occurs only after expansion of the DPPC monolayer from the LC phase though the LC/LE phase coexistence. Upon monolayer expansion, the domain boundaries move apart while attractive interactions among Fn molecules and between Fn and domain edges produce a tensile force on the proteins that initiates fibril assembly. The resulting fibrils have been characterized in situ by using fluorescence and light-scattering microscopy. We have found striking similarities between fibrils produced under DPPC monolayers and those found on cellular surfaces, including their assembly pathways.

Continuous and lurching traveling pulses in neuronal networks with delay and spatially decaying connectivity

Propagation of discharges in cortical and thalamic systems, which is used as a probe for examining network circuitry, is studied by constructing a one-dimensional model of integrate-and-fire neurons that are coupled by excitatory synapses with delay. Each neuron fires only one spike. The velocity and stability of propagating continuous pulses are calculated analytically. Above a certain critical value of the constant delay, these pulses lose stability. Instead, lurching pulses propagate with discontinuous and periodic spatio-temporal characteristics. The parameter regime for which lurching occurs is strongly affected by the footprint (connectivity) shape; bistability may occur with a square footprint shape but not with an exponential footprint shape. For strong synaptic coupling, the velocity of both continuous and lurching pulses increases logarithmically with the synaptic coupling strength gsyn for an exponential footprint shape, and it is bounded for a step footprint shape. We conclude that the differences in velocity and shape between the front of thalamic spindle waves in vitro and cortical paroxysmal discharges stem from their different effective delay; in thalamic networks, large effective delay between inhibitory neurons arises from their effective interaction via the excitatory cells which display postinhibitory rebound.

Functional reorganization in thalamocortical networks: Transition between spindling and delta sleep rhythms

Thalamic reticularis, thalamocortical, and cortical cells participate in the 7–14-hz spindling rhythm of early sleep and the slower delta rhythms of deeper sleep, with different firing patterns. In this case study, showing the interactions of intrinsic and synaptic properties, a change in the conductance of one kind of cell effectively rewires the thalamocortical circuit, leading to the transition from the spindling to the delta rhythm. The two rhythms make different uses of the fast (GABAA) and slow (GABAB) inhibition generated by the thalamic reticularis cells.

The interaction of Arp2/3 complex with actin: Nucleation, high affinity pointed end capping, and formation of branching networks of filaments

The Arp2/3 complex is a stable assembly of seven protein subunits including two actin-related proteins (Arp2 and Arp3) and five novel proteins. Previous work showed that this complex binds to the sides of actin filaments and is concentrated at the leading edges of motile cells. Here, we show that Arp2/3 complex purified from Acanthamoeba caps the pointed ends of actin filaments with high affinity. Arp2/3 complex inhibits both monomer addition and dissociation at the pointed ends of actin filaments with apparent nanomolar affinity and increases the critical concentration for polymerization at the pointed end from 0.6 to 1.0 μM. The high affinity of Arp2/3 complex for pointed ends and its abundance in amoebae suggest that in vivo all actin filament pointed ends are capped by Arp2/3 complex. Arp2/3 complex also nucleates formation of actin filaments that elongate only from their barbed ends. From kinetic analysis, the nucleation mechanism appears to involve stabilization of polymerization intermediates (probably actin dimers). In electron micrographs of quick-frozen, deep-etched samples, we see Arp2/3 bound to sides and pointed ends of actin filaments and examples of Arp2/3 complex attaching pointed ends of filaments to sides of other filaments. In these cases, the angle of attachment is a remarkably constant 70 ± 7°. From these in vitro biochemical properties, we propose a model for how Arp2/3 complex controls the assembly of a branching network of actin filaments at the leading edge of motile cells.

Simple neural networks for the amplification and utilization of small changes in neuron firing rates

I describe physiologically plausible “voter-coincidence” neural networks such that secondary “coincidence” neurons fire on the simultaneous receipt of sufficiently large sets of input pulses from primary sets of neurons. The networks operate such that the firing rate of the secondary, output neurons increases (or decreases) sharply when the mean firing rate of primary neurons increases (or decreases) to a much smaller degree. In certain sensory systems, signals that are generally smaller than the noise levels of individual primary detectors, are manifest in very small increases in the firing rates of sets of afferent neurons. For such systems, this kind of network can act to generate relatively large changes in the firing rate of secondary “coincidence” neurons. These differential amplification systems can be cascaded to generate sharp, “yes–no” spike signals that can direct behavioral responses.

Visualization of biochemical networks in living cells

Functional annotation of novel genes can be achieved by detection of interactions of their encoded proteins with known proteins followed by assays to validate that the gene participates in a specific cellular function. We report an experimental strategy that allows for detection of protein interactions and functional assays with a single reporter system. Interactions among biochemical network component proteins are detected and probed with stimulators and inhibitors of the network. In addition, the cellular location of the interacting proteins is determined. We used this strategy to map a signal transduction network that controls initiation of translation in eukaryotes. We analyzed 35 different pairs of full-length proteins and identified 14 interactions, of which five have not been observed previously, suggesting that the organization of the pathway is more ramified and integrated than previously shown. Our results demonstrate the feasibility of using this strategy in efforts of genomewide functional annotation.

The generation of oscillations in networks of electrically coupled cells

In several biological systems, the electrical coupling of nonoscillating cells generates synchronized membrane potential oscillations. Because the isolated cell is nonoscillating and electrical coupling tends to equalize the membrane potentials of the coupled cells, the mechanism underlying these oscillations is unclear. Here we present a dynamic mechanism by which the electrical coupling of identical nonoscillating cells can generate synchronous membrane potential oscillations. We demonstrate this mechanism by constructing a biologically feasible model of electrically coupled cells, characterized by an excitable membrane and calcium dynamics. We show that strong electrical coupling in this network generates multiple oscillatory states with different spatio-temporal patterns and discuss their possible role in the cooperative computations performed by the system.

Intrinsic noise in gene regulatory networks

Cells are intrinsically noisy biochemical reactors: low reactant numbers can lead to significant statistical fluctuations in molecule numbers and reaction rates. Here we use an analytic model to investigate the emergent noise properties of genetic systems. We find for a single gene that noise is essentially determined at the translational level, and that the mean and variance of protein concentration can be independently controlled. The noise strength immediately following single gene induction is almost twice the final steady-state value. We find that fluctuations in the concentrations of a regulatory protein can propagate through a genetic cascade; translational noise control could explain the inefficient translation rates observed for genes encoding such regulatory proteins. For an autoregulatory protein, we demonstrate that negative feedback efficiently decreases system noise. The model can be used to predict the noise characteristics of networks of arbitrary connectivity. The general procedure is further illustrated for an autocatalytic protein and a bistable genetic switch. The analysis of intrinsic noise reveals biological roles of gene network structures and can lead to a deeper understanding of their evolutionary origin.