Search

The problem of coordinating genome replication with cell growth in bacteria was posed over four decades ago. Unlike for eukaryotes, this problem has not been completely solved even for Escherichia coli, which has been comprehensively studied by molecular biologists, to say nothing of other bacteria. Current models of the bacterial life cycle solve the coupling problem by introducing a phenomenological hypothesis that considers the dynamic coordination of growth and replication but does not unveil the underlying molecular mechanisms. Here we review the mechanisms regulating genome replication initiation with regards to their coupling to growth processes in the three best investigated bacterial species: E. coli, Bacillus subtilis, and Caulobacter crescentus. A putative correlation between the type of cell growth laws and the actual mechanisms regulating the replication of DNA formed during the process of evolution in various classes of bacteria, is discussed, including those intracellular parasites in which degenerative evolution his discarded most of their genomes. We contemplate the concept of a universal growth law for bacterial cells and some features in the formation of a primitive negative replication regulating mechanism in the context of the coupling problem.

A mathematical model of the intracellular reproduction of the influenza virus has been constructed, which describes the regulation of transcription, translation, replication, and assembly of the virus in a sensitive noninterferonogenic cell. A structural and parametric verification of the model has been performed, which involves an analysis of alternative mechanisms of regulation of the key events of the reproduction of the influenza virus at the stages of the primary and secondary transcription. It was shown that the dynamics of the concentrations of the system components is in qualitative agreement with the experimental data. An analysis of the model parameters showed a high sensitivity of the model to changes in the parameters that control the efficiency of the binding of the viral polymerase to virus-specific RNAs at the stages of initiation and elongation of transcription and replication of the viral genome. The results obtained suggest that the stages of the virus reproduction sensitive to changes in the parameters are potential targets for drugs. Using the model, a minimal multiplicity of infection at which the formation of viral particles is observed was predicted.

Regulation of plant growth and development by auxin is mediated by the hormone controlled distribution and dose-dependent mechanisms of its action. A mathematical model is proposed, which described the distribution of auxin in the cells extending along the central axis of the Arabidopsis thaliana root. This model reproduces qualitatively both auxin distribution in cells of the root central axis under the normal conditions and under the conditions of decreased active transport and the recovery of auxin distribution and related meristem restoration during root regeneration after ablation of its tip. Different types of distribution of the auxin concentration over the vertical root axis are described, possible variants of root growth and lateral roots formation are proposed, and biological interpretation is given to different regimes of model behavior.

Unlabelled: Gene expression is known to correlate with degree of codon bias in many unicellular organisms. However, such correlation is absent in some organisms. Recently we demonstrated that inverted complementary repeats within coding DNA sequence must be considered for proper estimation of translation efficiency, since they may form secondary structures that obstruct ribosome movement. We have developed a program for estimation of potential coding DNA sequence expression in defined unicellular organism using its genome sequence. The program computes elongation efficiency index. Computation is based on estimation of coding DNA sequence elongation efficiency, taking into account three key factors: codon bias, average number of inverted complementary repeats, and free energy of potential stem-loop structures formed by the repeats. The influence of these factors on translation is numerically estimated. An optimal proportion of these factors is computed for each organism individually. Quantitative translational characteristics of 384 unicellular organisms (351 bacteria, 28 archaea, 5 eukaryota) have been computed using their annotated genomes from NCBI GenBank. Five potential evolutionary strategies of translational optimization have been determined among studied organisms. A considerable difference of preferred translational strategies between Bacteria and Archaea has been revealed. Significant correlations between elongation efficiency index and gene expression levels have been shown for two organisms (S. cerevisiae and H. pylori) using available microarray data. The proposed method allows to estimate numerically the coding DNA sequence translation efficiency and to optimize nucleotide composition of heterologous genes in unicellular organisms., Availability: http://www.mgs.bionet.nsc.ru/mgs/programs/eei-calculator/.

The review describes several modules of the GeneExpress integrated computer system concerning the regulation of gene expression in eukaryotes. Approaches to the presentation of experimental data in databases are considered. The employment of GeneExpress in computer analysis and modeling of the organization and function of genetic systems is illustrated with examples. GeneExpress is available at http://wwwmgs.bionet.nsc.ru/mgs/gnw/.

Development of methods for mathematical simulation of biological systems and building specific simulations is an important trend of bioinformatics development. Here we describe the method of generalized chemokinetic simulation generating flexible and adequate simulations of various biological systems. Adequate simulations of complex nonlinear gene networks--control system of cholesterol by synthesis in the cell and erythrocyte differentiation and maturation--are given as the examples. The simulations were expressed in terms of unit processes--biochemical reactions. Optimal sets of parameters were determined and the systems were numerically simulated under various conditions. The simulations allow us to study possible functional conditions of these gene networks, calculate consequences of mutations, and define optimal strategies for their correction including therapeutic ones. Graphical user interface for these simulations is available at http://wwwmgs.bionet.nsc.ru/systems/MGL/GeneNet/.

Theoretical investigation of properties of assumed gene networks constructed from elementary units of two types, genetic elements and control links, was carried out. A test was formulated for a subclass of such networks with cyclic structure called S(n,k)-networks allowing calculation-free prediction of the network limiting properties (the presence/absence and number of stationery and/or cyclic functioning modes) from a graph of the network structure. The obtained data can be useful for constructing gene networks with predefined properties.

Four series of plasmids (pNSI, pNSII, pNLI, and pNLII) with artificial polycistrons containing the lacZ test gene were constructed. These plasmids coded for polycistronic mRNAs with two different types of cistron (orfZ and lacZ) coupling: in pNSI and pNLI, the orfZ termination codon and the lacZ initiation codon overlapped (type I); in pNSII and pNLII, the orfZ termination codon, was located upstream of the lacZ SD sequence. The length of the orfZ cistron was 60 bp in pNSI and pNSII or 300 bp in pNLI and pNLII. Plasmids with the same type of cistron coupling contained the same lacZ translation initiation region, whereas the structure of the orfZ translation initiation region varied, thereby providing varying efficiency of the orfZ gene translation. The effect of these variations on the efficiency of the lacZ gene translation was evaluated by direct measurement of the beta-galactosidase activity in Escherichia coli cells transformed with the corresponding plasmids. We found that the level of translation of the distal lacZ gene depended on the ribosome stream from the proximal gene and was maximal at the optimal ribosome stream level, which, in turn, depended on the type of cistron coupling.

The role of "stream" of ribosomes upon translation of polycistronic mRNAs has been studied using an artificial polycistron. It has been found that the levels of activation of cistron Ci + 1 out of two adjacent cistrons (Ci and Ci + 1) depends, in addition to earlier described effects of mutual arrangement of initiation and termination signals, also on efficiency of translation of the foregoing cistron Ci. The results obtained lead to the conclusion that in polycistronic systems the levels of translation of cistron Ci + 1 can be regulated by "stream" of ribosomes resulted from translation of the proximal cistron Ci.

Models of processive and distributive DNA synthesis and degradation catalysed by matrix enzymes were investigated. Distribution of polymer products dependent on the reaction model chosen, on the type of the matrix and on the enzyme-matrix initial concentration ratio were determined by methods of numerical modeling. Conditions were found where the scattering of the reaction polymer products was minimal. The homopolymer matrix choice of experimental condition may generate a distribution of product that obeys the Poisson distribution. Numerical investigations of polymerization (hydrolysis) processes showed that for a number of heteropolymer matrixes a distribution may exist with scattering much less than that for homopolymer matrixes. Conditions are found when processive and distributive models give different distributions of the reaction product.

The role of the translational terminator and initiator signals arrangement for two adjacent genes in polycistronic mRNA has been studied. Semisynthetic beta-galactosidase gene (lacZ) of E. coli and fragment of phage M13 DNA (with promoter PVIII, gene IX, and part of gene VIII) were used for constructing of the IX-VIII-lacZ artificial polycistronic operon. Cloning of the constructs into pBR322 vector resulted in a number of pLZ381N plasmids differing by the mutual arrangement of gene VIII translation terminator codon and SD site and initiator codon (SD-ATG-region) of lacZ gene. The mutual arrangement of gene VIII terminator codon and SDlacZ-ATG region has been altered by means of deletions and insertions that have not affected lacZ translation initiation signals. The beta-galactosidase (beta-Gal) synthesis in E. coli harbouring different types of pLZ381N plasmids has been found to depend on type of cistron coupling (gene VIII and lacZ). The overlapping of terminator and initiator codons (ATGA) for genes VIII and lacZ (type I of polycistrons) provide approximately equal translational level for both cistrons. On the other side, levels of beta-Gal synthesis in case of polycistrons type II (gene VIII stop-codon position at the beginning of SDlacZ or 10 nucleotides upstream) were 20-30 times as high as for type I. Differences in beta-Gal levels have also been found for variants of VIII-lacZ coupling in types IV and III polycistrons (the SDlacZ-ATG region in 27-50 nucleotides downstream from the proximal cistron VIII stop-codon, which, in turn, is 41 nucleotides upstream this terminator). These data cannot be explained on the basis of possible secondary structure including the SDlacZ-ATG region and other parts of polycistronic mRNA. In all these cases similarly stable stem-loop structures have been found. Therefore, the arrangement of the translation termination and initiation signals for two adjacent genes in essential for distal gene translation efficiency. One can imagine that ribosome or its 30S subpartical, stalling on the proximal gene terminator codon, affects the distal gene translation initiation.