The evolution of controlled multitasked gene networks: The role of introns and other noncoding RNAs in the development of complex organisms

Eukaryotic phenotypic diversity arises from multitasking of a core proteome of limited size. Multitasking is routine in computers, as well as in other sophisticated information systems, and requires multiple inputs and outputs to control and integrate network activity. Higher eukaryotes have a mosaic gene structure with a dual output, mRNA (protein-coding) sequences and introns, which are released from the pre-mRNA by posttranscriptional processing. Introns have been enormously successful as a class of sequences and comprise up to 95% of the primary transcripts of protein-coding genes in mammals. In addition, many other transcripts (perhaps more than half) do not encode proteins at all, but appear both to be developmentally regulated and to have genetic function. We suggest that these RNAs (eRNAs) have evolved to function as endogenous network control molecules which enable direct gene-gene communication and multitasking of eukaryotic genomes. Analysis of a range of complex genetic phenomena in which RNA is involved or implicated, including co-suppression, transgene silencing, RNA interference, imprinting, methylation, and transvection, suggests that a higher-order regulatory system based on RNA signals operates in the higher eukaryotes and involves chromatin remodeling as well as other RNA-DNA, RNA-RNA, and RNA-protein interactions. The evolution of densely connected gene networks would be expected to result in a relatively stable core proteome due to the multiple reuse of components, implying,that cellular differentiation and phenotypic variation in the higher eukaryotes results primarily from variation in the control architecture. Thus, network integration and multitasking using trans-acting RNA molecules produced in parallel with protein-coding sequences may underpin both the evolution of developmentally sophisticated multicellular organisms and the rapid expansion of phenotypic complexity into uncontested environments such as those initiated in the Cambrian radiation and those seen after major extinction events.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Complex natural resonances of conducting planar objects buried in a dielectric half-space

A scheme is presented to incorporate a mixed potential integral equation (MPIE) using Michalski's formulation C with the method of moments (MoM) for analyzing the scattering of a plane wave from conducting planar objects buried in a dielectric half-space. The robust complex image method with a two-level approximation is used for the calculation of the Green's functions for the half-space. To further speed up the computation, an interpolation technique for filling the matrix is employed. While the induced current distributions on the object's surface are obtained in the frequency domain, the corresponding time domain responses are calculated via the inverse fast Fourier transform (FFT), The complex natural resonances of targets are then extracted from the late time response using the generalized pencil-of-function (GPOF) method. We investigate the pole trajectories as we vary the distance between strips and the depth and orientation of single, buried strips, The variation from the pole position of a single strip in a homogeneous dielectric medium was only a few percent for most of these parameter variations.

Complementation analysis of the flavivirus Kunjin NS3 and NS5 proteins defines the minimal regions essential for formation of a replication complex and shows a requirement of NS3 in cis for virus assembly

We have previously reported successful trans-complementation of defective Kunjin virus genomic RNAs with a range of large lethal deletions in the nonstructural genes NSI, NS3, and NS5 (A. A. Khromykh et al., J. Virol. 74:3253-3263, 2000). In this study we have mapped further the minimal region in the NS5 gene essential for efficient trans-complementation of genome-length RNAs in repBHK cells to the first 316 of the 905 codons. To allow amplification and easy detection of complemented defective RNAs with deletions apparently affecting virus assembly, we have developed a dual replicon complementation system. In this system defective replicon RNAs with a deletion(s) in the nonstructural genes also encoded the puromycin resistance gene (PAC gene) and the reporter gene for beta-galactosidase (beta-Gal). Complementation of these defective replicon RNAs in repBHK cells resulted in expression of PAC and beta-Gal which allowed establishment of cell lines stably producing replicating defective RNAs by selection with puromycin and comparison of replication efficiencies of complemented defective RNAs by beta-Gal assay. Using this system we demonstrated that deletions in the C-terminal 434 codons of NS3 (codons 178 to 611) were complemented for RNA replication, while any deletions in the first 178 codons were not. None of the genome-length RNAs containing deletions in NS3 shown to be complementable for RNA replication produced secreted defective viruses during complementation in repBHK cells. In contrast, structural proteins produced from these complemented defective RNAs were able to package helper replicon RNA. The results define minimal regions in the NS3 and NS5 genes essential for the formation of complementable replication complex and show a requirement of NS3 in cis for virus assembly.

Gene-Environment Interaction Effects on Behavioral Variation and Risk of Complex Disorders: The Example of Alcoholism and Other Psychiatric Disorders

There have been few replicated examples of genotype x environment interaction effects on behavioral variation or risk of psychiatric disorder. We review some of the factors that have made detection of genotype x environment interaction effects difficult, and show how genotype x shared environment interaction (GxSE) effects are commonly confounded with genetic parameters in data from twin pairs reared together. Historic data on twin pairs reared apart can in principle be used to estimate such GxSE effects, but have rarely been used for this purpose. We illustrate this using previously published data from the Swedish Adoption Twin Study of Aging (SATSA), which suggest that GxSE effects could account for as much as 25% of the total variance in risk of becoming a regular smoker. Since few separated twin pairs will be available for study in the future, we also consider methods for modifying variance components linkage analysis to allow for environmental interactions with linked loci.

Characterization of [Fe(AMN3S3sarH)]3+ - a rigorously low-spin iron(II) complex

The iron(II) complex [Fe(AMN(3)S(3)sarH)](ClO4)(3).3H(2)O (AMN(3)S(3)sarH = 8-ammonio-1-methyl-3,13,16-trithia-6,10,19-triazabicyclo[6.6.6]icosane) has been synthesized and characterized by single crystal structure and spectroscopic methods. The Fe(II)-S(thiaether) bond lengths are short, indicative of a large degree of metal-ligand orbital mixing (pi-acceptor character) of the thiaether ligand. The complex is stable to metal centred oxidation. (C) 2002 Elsevier Science Ltd. All rights reserved.

The influence of cis/trans isomerism on the physical properties of a cyano-bridged dinuclear mixed valence complex

The cyano-bridged complexes cis-[L14CoIIINCFeII(CN)5]– and cis-[L14CoIIINCFeIII(CN)5] (L14= 6-methyl-1,4,8,11-tetraazacyclotetradecan-6-amine) are prepared and characterised spectroscopically, electrochemically and structurally: Na{cis-[L14CoIIINCFeII(CN)5]}·9H2O, monoclinic space group P21/c, a= 14.758(3), b= 10.496(1), c= 19.359(3) , = 92.00(2)°, Z= 4; cis-[L14CoIIINCFeIII(CN)5]·4H2O, orthorhombic space group P212121, a= 9.492(1), b= 14.709(2), c= 18.760(3) , Z= 4. In both complexes, the pendant amine is cis to the bridging cyanide ligand. An analysis of the metal-to-metal charge transfer (MMCT) transition in these systems with Hush theory has been carried out. This has revealed that the change in the configuration of the macrocycle both decreases the redox isomer energy difference (E1/2) and increases the reorganisational energy () of the cis-[L14CoIIINCFeII(CN)5]– complex with respect to the trans-[L14CoIIINCFeII(CN)5]– complex, the result being that both isomers display an MMCT transition of similar energy. The variation in redox isomer energy differences of the configurational isomers has been related to strain energy differences by molecular mechanics analysis of the [CoL14Cl]2+/+ precursor complexes.

A system for the heterologous expression of complex redox proteins in Rhodobacter capsulatus: characterisation of recombinant sulphite : cytochrome c oxidoreductase from Starkeya novella

The phototrophic purple non-sulfur bacterium Rhodobacter capsulatus expresses a wide variety of complex redox proteins in response to changing environmental conditions. Here we report the construction and evaluation of an expression system for recombinant proteins in that organism which makes use of the dor promoter from the same organism. A generic expression vector, pDorEX, was constructed and used to express sulphite:cytochrome c oxidoreductase from Starkeya novella, a heterodimeric protein containing both molybdenum and haem c. The recombinant protein was secreted to the periplasm and its biochemical properties were very similar to those of the native enzyme. The pDorEX system therefore seems to be potentially useful for heterologous expression of multi-subunit proteins containing complex redox cofactors. (C) 2002 Federation of European Biochemical Societies. Published by Elsevier Science B.V. All rights reserved.

Role of peripheral afference during acquisition of a complex coordination task

It has long been supposed that the interference observed in certain patterns of coordination is mediated, at least in part, by peripheral afference from the moving limbs. We manipulated the level of afferent input, arising from movement of the opposite limb, during the acquisition of a complex coordination task. Participants learned to generate flexion and extension movements of the right wrist, of 75degrees amplitude, that were a quarter cycle out of phase with a 1-Hz sinusoidal visual reference signal. On separate trials, the left wrist either was at rest, or was moved passively by a torque motor through 50degrees, 75degrees or 100degrees, in synchrony with the reference signal. Five acquisition sessions were conducted on successive days. A retention session was conducted I week later. Performance was initially superior when the opposite limb was moved passively than when it was static. The amplitude and frequency of active movement were lower in the static condition than in the driven conditions and the variation in the relative phase relation across trials was greater than in the driven conditions. In addition, the variability of amplitude, frequency and the relative phase relation during each trial was greater when the opposite limb was static than when driven. Similar effects were expressed in electromyograms. The most marked and consistent differences in the accuracy and consistency of performance (defined in terms of relative phase) were between the static condition and the condition in which the left wrist was moved through 50degrees. These outcomes were exhibited most prominently during initial exposure to the task. Increases in task performance during the acquisition period, as assessed by a number of kinematic variables, were generally well described by power functions. In addition, the recruitment of extensor carpi radialis (ECR), and the degree of co-contraction of flexor carpi radialis and ECR, decreased during acquisition. Our results indicate that, in an appropriate task context, afferent feedback from the opposite limb, even when out of phase with the focal movement, may have a positive influence upon the stability of coordination.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.