Symbiosis and the origin of life

The paper uses chemical kinetic arguments and illustrations by computer modelling to discuss the origin and evolution of life. Complex self-reproducing chemical systems cannot arise spontaneously, whereas simple auto-catalytic systems can, especially in an irradiated aqueous medium. Self-reproducing chemical particles of any complexity, in an appropriate environment, have a self-regulating property which permits long-term survival. However, loss of materials from the environment can lead to continuing decay which is circumvented by physical union between different kinds of self-reproducing particles. The increasing complexity produced by such unions (symbioses) is irreversible so that the chemical system evolves. It is suggested that evolution by successive symbioses brough about the change from simple, spontaneously arising, auto-catalytic particles to complex prokaryotic cells.     

A model for oscillating chemical reactions.

A simple theoretical model for a class of oscillating chemical reactions is investigated, which is inspired by Degn's analysis of the mechanism of Belousov-Zhabotinsky reactions. Although the model involves only one auto-catalytic path, as distinct from models presented hitherto, a limit-cycle behaviour is found for the concentrations, corresponding to a definite region in the space of controllable parameters. Characteristics of the coherent oscillations are investigated together with the behaviour in the induction period. The critical slowing-down around the onset of instability is shown to behave very much like that found in the equilibrium phase transition. The existence of an "anti-Curie" point is found as well as the Curie point, corresponding to the fact that the domain for sustained oscillation is closed in certain directions. Possible extensions of the model and ways of improving the theory are discussed.     

Evolution of model proteins on a foldability landscape

We model the evolution of simple lattice proteins as a random walk in a fitness landscape, where the fitness represents the ability of the protein to fold. At higher selective pressure, the evolutionary trajectories are confined to neutral networks where the native structure is conserved and the dynamics are non self-averaging and nonexponential. The optimizability of the corresponding native structure has a strong effect on the size of these neutral networks and thus on the nature of the evolutionary process. Proteins 29:461–466, 1997. © 1997 Wiley-Liss, Inc.     

The mechanisms of hydrothermal deconstruction of lignocellulose: New insights from thermal-analytical and complementary studies

Differential Scanning Calorimetry, Dynamic Mechanical Thermal Analysis, gravimetric and chemical techniques have been used to study hydrothermal reactions of straw biomass. Exothermic degradation initiates above 195 °C, due to breakdown of the xylose ring from hemicellulose, which may be similar to reactions occurring during the early stage pyrolysis of dry biomass, though activated at lower temperature through water mediation. The temperature and magnitude of the exotherm reduce with increasing acid concentration, suggesting a reduction in activation energy and a change in the balance of reaction pathways. The presence of xylan oligomers in auto-catalytic hydrolysates is believed to be due to a low rate constant rather than a specific reaction mechanism. The loss of the lignin glass transition indicates that the lignin phase is reorganised under high temperature auto-catalytic conditions, but remains partially intact under lower temperature acid-catalytic conditions. This shows that lignin degradation reactions are activated thermally but are not effectively catalysed by aqueous acid.     

In silico two-hybrid system for the selection of physically interacting protein pairs

Deciphering the interaction links between proteins has become one of the main tasks of experimental and bioinformatic methodologies. Reconstruction of complex networks of interactions in simple cellular systems by integrating predicted interaction networks with available experimental data is becoming one of the most demanding needs in the postgenomic era. On the basis of the study of correlated mutations in multiple sequence alignments, we propose a new method (in silico two-hybrid, i2h) that directly addresses the detection of physically interacting protein pairs and identifies the most likely sequence regions involved in the interactions. We have applied the system to several test sets, showing that it can discriminate between true and false interactions in a significant number of cases. We have also analyzed a large collection of E. coli protein pairs as a first step toward the virtual reconstruction of its complete interaction network.     

Robustness and complexity co-constructed in multimodal signalling networks

In animal communication, signals are frequently emitted using different channels (e.g. frequencies in a vocalization) and different modalities (e.g. gestures can accompany vocalizations). We explore two explanations that have been provided for multimodality: (i) selection for high information transfer through dedicated channels and (ii) increasing fault tolerance or robustness through multichannel signals. Robustness relates to an accurate decoding of a signal when parts of a signal are occluded. We show analytically in simple feed-forward neural networks that while a multichannel signal can solve the robustness problem, a multimodal signal does so more effectively because it can maximize the contribution made by each channel while minimizing the effects of exclusion. Multimodality refers to sets of channels where within each set information is highly correlated. We show that the robustness property ensures correlations among channels producing complex, associative networks as a by-product. We refer to this as the principle of robust overdesign. We discuss the biological implications of this for the evolution of combinatorial signalling systems; in particular, how robustness promotes enough redundancy to allow for a subsequent specialization of redundant components into novel signals.     

A model of autocatalytic replication

Summary The catalytic effects that existing polymer chains have on the formation of new chains are modeled using ideas from spin glasses and neural networks. Computer simulation shows that isolated groups of chains in this model are capable of accurately replicating a wide variety of complex structures without templating. Replication in the model arises spontaneously and rapidly, leading to an extremely simple realization of a system exhibiting Darwinian evolution.     

The foldability landscape of model proteins

Molecular evolution may be considered as a walk in a multidimensional fitness landscape, where the fitness at each point is associated with features such as the function, stability, and survivability of these molecules. We present a simple model for the evolution of protein sequences on a landscape with a precisely defined fitness function. We use simple lattice models to represent protein structures, with the ability of a protein sequence to fold into the structure with lowest energy, quantified as the foldability, representing the fitness of the sequence. The foldability of the sequence is characterized based on the spin glass model of protein folding. We consider evolution as a walk in this foldability landscape and study the nature of the landscape and the resulting dynamics. Selective pressure is explicitly included in this model in the form of a minimum foldability requirement. We find that different native structures are not evenly distributed in interaction space, with similar structures and structures with similar optimal foldabilities clustered together. Evolving proteins marginally fulfill the selective criteria of foldability. As the selective pressure is increased, evolutionary trajectories become increasingly confined to “neutral networks,” where the sequence and the interactions can be significantly changed while a constant structure is maintained. © 1997 John Wiley & Sons, Inc. Biopoly 42: 427–438, 1997     

Super folds,networks, and barriers

Exhaustive enumeration of sequences and folds is conducted for a simple lattice model of conformations, sequences, and energies. Examination of all foldable sequences and their nearest connected neighbors (sequences that differ by no more than a point mutation) illustrates the following: (i) There exist unusually large number of sequences that fold into a few structures (super‐folds). The same observation was made experimentally and computationally using stochastic sampling and exhaustive enumeration of related models. (ii) There exist only a few large networks of connected sequences that are not restricted to one fold. These networks cover a significant fraction of fold spaces (super‐networks). (iii) There exist barriers in sequence space that prevent foldable sequences of the same structure to “connect” through a series of single point mutations (super‐barrier), even in the presence of the sequence connection between folds. While there is ample experimental evidence for the existence of super‐folds, evidence for a super‐network is just starting to emerge. The prediction of a sequence barrier is an intriguing characteristic of sequence space, suggesting that the overall sequence space may be disconnected. The implications and limitations of these observations for evolution of protein structures are discussed. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Epistasis can lead to fragmented neutral spaces and contingency in evolution

In evolution, the effects of a single deleterious mutation can sometimes be compensated for by a second mutation which recovers the original phenotype. Such epistatic interactions have implications for the structure of genome space--namely, that networks of genomes encoding the same phenotype may not be connected by single mutational moves. We use the folding of RNA sequences into secondary structures as a model genotype-phenotype map and explore the neutral spaces corresponding to networks of genotypes with the same phenotype. In most of these networks, we find that it is not possible to connect all genotypes to one another by single point mutations. Instead, a network for a phenotypic structure with n bonds typically fragments into at least 2(n) neutral components, often of similar size. While components of the same network generate the same phenotype, they show important variations in their properties, most strikingly in their evolvability and mutational robustness. This heterogeneity implies contingency in the evolutionary process.     

The relative efficiency of modular and non-modular networks of different size

Most biological networks are modular but previous work with small model networks has indicated that modularity does not necessarily lead to increased functional efficiency. Most biological networks are large, however, and here we examine the relative functional efficiency of modular and non-modular neural networks at a range of sizes. We conduct a detailed analysis of efficiency in networks of two size classes: ‘small’ and ‘large’, and a less detailed analysis across a range of network sizes. The former analysis reveals that while the modular network is less efficient than one of the two non-modular networks considered when networks are small, it is usually equally or more efficient than both non-modular networks when networks are large. The latter analysis shows that in networks of small to intermediate size, modular networks are much more efficient that non-modular networks of the same (low) connective density. If connective density must be kept low to reduce energy needs for example, this could promote modularity. We have shown how relative functionality/performance scales with network size, but the precise nature of evolutionary relationship between network size and prevalence of modularity will depend on the costs of connectivity.     

Evolution of robustness to damage in artificial 3-dimensional development

GReaNs is an Artificial Life platform we have built to investigate the general principles that guide evolution of multicellular development and evolution of artificial gene regulatory networks. The embryos develop in GReaNs in a continuous 3-dimensional (3D) space with simple physics. The developmental trajectories are indirectly encoded in linear genomes. The genomes are not limited in size and determine the topology of gene regulatory networks that are not limited in the number of nodes. The expression of the genes is continuous and can be modified by adding environmental noise. In this paper we evolved development of structures with a specific shape (an ellipsoid) and asymmetrical pattering (a 3D pattern inspired by the French flag problem), and investigated emergence of the robustness to damage in development and the emergence of the robustness to noise. Our results indicate that both types of robustness are related, and that including noise during evolution promotes higher robustness to damage. Interestingly, we have observed that some evolved gene regulatory networks rely on noise for proper behaviour.     

Analysis of metabolic network based on conservation of molecular structure

Recent work has revealed much about chemical reactions inside hundreds of organisms as well as universal characteristics of metabolic networks, which shed light on the evolution of the networks. However, characteristics of individual metabolites have been neglected. For example, some carbohydrates have structures that are decomposed into small molecules by metabolic reactions, but coenzymes such as ATP are mostly preserved. Such differences in metabolite characteristics are important for understanding the universal characteristics of metabolic networks. To quantify the structure conservation of metabolites, we defined the "structure conservation index" (SCI) for each metabolite as the fraction of metabolite atoms restored to their original positions through metabolic reactions. As expected, coenzymes and coenzyme-like metabolites that have reaction loops in the network show a higher SCI. Using the index, we found that the sum of metabolic fluxes is negatively correlated with the structure preservation of metabolite. Also, we found that each reaction path around high SCI metabolites changes independently, while changes in reaction paths involving low SCI metabolites coincide through evolution processes. These correlations may provide a clue to universal properties of metabolic networks.     

Stochastic P systems and the simulation of biochemical processes with dynamic compartments

We introduce a sequential rewriting strategy for P systems based on Gillespie's stochastic simulation algorithm, and show that the resulting formalism of stochastic P systems makes it possible to simulate biochemical processes in dynamically changing, nested compartments. Stochastic P systems have been implemented using the spatially explicit programming language MGS. Implementation examples include models of the Lotka-Volterra auto-catalytic system, and the life cycle of the Semliki Forest virus.     

Evolutionary bioscience as regulatory systems biology

At present several entirely different explanatory approaches compete to illuminate the mechanisms by which animal body plans have evolved. Their respective relevance is briefly considered here in the light of modern knowledge of genomes and the regulatory processes by which development is controlled. Just as development is a system property of the regulatory genome, causal explanation of evolutionary change in developmental process must be considered at a system level. Here I enumerate some mechanistic consequences that follow from the conclusion that evolution of the body plan has occurred by alteration of the structure of developmental gene regulatory networks. The hierarchy and multiple additional design features of these networks act to produce Boolean regulatory state specification functions at upstream phases of development of the body plan. These are created by the logic outputs of network subcircuits, and in modern animals these outputs are impervious to continuous adaptive variation unlike genes operating more peripherally in the network.     

Receptors rather than signals change in expression in four physiological regulatory networks during evolutionary divergence in threespine stickleback

The molecular mechanisms underlying behavioural evolution following colonization of novel environments are largely unknown. Molecules that interact to control equilibrium within an organism form physiological regulatory networks. It is essential to determine whether particular components of physiological regulatory networks evolve or if the network as a whole is affected in populations diverging in behavioural responses, as this may affect the nature, amplitude and number of impacted traits. We studied the regulation of four physiological regulatory networks in freshwater and marine populations of threespine stickleback raised in a common environment, which were previously characterized as showing evolutionary divergence in behaviour and stress reactivity. We measured nineteen components of these networks (ligands and receptors) using mRNA and monoamine levels in the brain, pituitary and interrenal gland, as well as hormone levels. Freshwater fish showed higher expression in the brain of adrenergic (adrb2a), serotonergic (htr2a) and dopaminergic (DRD2) receptors, but lower expression of the htr2b receptor. Freshwater fish also showed higher expression of the mc2r receptor of the glucocorticoid axis in the interrenals. Collectively, our results suggest that the inheritance of the regulation of these networks may be implicated in the evolution of behaviour and stress reactivity in association with population divergence. Our results also suggest that evolutionary change in freshwater threespine stickleback may be more associated with the expression of specific receptors rather than with global changes of all the measured constituents of the physiological regulatory networks.     

Degree distribution in plant–animal mutualistic networks: forbidden links or random interactions?

Recent studies show that ecological interaction networks depart from the "scale-free" topologies observed in many other real world networks. Such a departure has been hypothesized to result from non-matching biological attributes of species, such as phenology or morphology, that prevent the occurrence of certain interactions ("forbidden links"). Here I compare the topology of 17 plant–animal mutualistic networks with that predicted by a simple null model that assumes that a species' degree (number of interspecific interactions) is a function of its frequency of interaction. The topology predicted by this null model is strikingly close to that observed in the real networks. Thus, this null model provides a simple alternative interpretation of patterns observed in ecological interaction networks that does not require the existence of non-matching species traits.     

Emergence of complex social networks from spatial structure and rules of thumb: a modelling approach

Individual-based computer models show that simple heuristic governing individuals’ behavior may suffice to generate complex patterns of social behavior at the group level such as those observed in animal societies. ‘GrooFiWorld’ is an example of such kind of computer models. In this model, self-organization and simple behavioral rules generate complex patterns of social behavior like those described in tolerant and intolerant societies of macaques. Social complexity results from the socio-spatial structure of the group, the nature of which is, in turn, a side-effect of intensity of aggression. The model suggests that a similar mechanism may give rise to complex social structures in macaques. It is, however, unknown if the spatial structure of the model and that of macaques are indeed similar. Here we used social networks analysis as a proxy for spatial structure of the group. Our findings show that the social networks of the model share similar qualitative features with those of macaques. As group size increases, the density and the average individual eigenvector centrality decrease and the modularity and centralization of the network increase. In social networks emerging from simulations resembling intolerant societies the density is lower, the modularity and centralization are higher, and the individuals ranking higher in the dominance hierarchy are more central than in the social networks emerging from simulations resembling egalitarian societies. Given the qualitative similarity between the social networks of the model and that of empirical data, our results suggest that the spatial structure of macaques is similar to that of the model. It seems thus plausible that, as in the model, the spatial structure combined with simple behavioral rules plays a role in the emergence of complex social networks and complex social behavior in macaques.     

Highly covarying residues have a functional role in antibody constant domains

The ability to generate and design antibodies recognizing specific targets has revolutionized the pharmaceutical industry and medical imaging. Engineering antibody therapeutics in some cases requires modifying their constant domains to enable new and altered interactions. Engineering novel specificities into antibody constant domains has proved challenging due to the complexity of inter‐domain interactions. Covarying networks of residues that tend to cluster on the protein surface and near binding sites have been identified in some proteins. However, the underlying role these networks play in the protein resulting in their conservation remains unclear in most cases. Resolving their role is crucial, because residues in these networks are not viable design targets if their role is to maintain the fold of the protein. Conversely, these networks of residues are ideal candidates for manipulating specificity if they are primarily involved in binding, such as the myriad interdomain interactions maintained within antibodies. Here, we identify networks of evolutionarily‐related residues in C‐class antibody domains by evaluating covariation, a measure of propensity with which residue pairs vary dependently during evolution. We computationally test whether mutation of residues in these networks affects stability of the folded antibody domain, determining their viability as design candidates. We find that members of covarying networks cluster at domain‐domain interfaces, and that mutations to these residues are diverse and frequent during evolution, precluding their importance to domain stability. These results indicate that networks of covarying residues exist in antibody domains for functional reasons unrelated to thermodynamic stability, making them ideal targets for antibody design. Proteins 2013. © 2012 Wiley Periodicals, Inc.