All issues

In this paper we study the rotational oscillations of the nitrous bases forming a central pair in a short DNA fragment consisting of three base pairs. A simple mechanical analog of the fragment where the bases are imitated by pendulums and the interactions between pendulums — by springs, has been constructed. We derived Lagrangian of the model system and the nonlinear equations of motions. We found solutions in the homogeneous case when the fragment considered consists of identical base pairs: Adenine-Thymine (AT- pair) or Guanine-Cytosine (GC-pair). The trajectories of the model system in the configuration space were also constructed.

It is known that in the native state the DNA molecule always contains some amount of locally unwound regions, often called the open states of DNA. It is believed that these states play an important role in DNA-protein recognition and that the study of the open states dynamics may shed further light on the mechanisms of regulation of transcription and replication. In this paper we consider the effect of the thermostat on the movement of the open states in the artificial sequence consisting of four homogeneous regions. We construct the energetic profile of the sequence and investigate the trajectories of the movement of the open states under the action of a random force.

A vertically distributed three-component model of marine ecosystem is considered. State of the plankton community with nutrients is analyzed under the active movement of zooplankton in a vertical column of water. The necessary conditions of the Turing instability in the vicinity of the spatially homogeneous equilibrium are obtained. Stability of the spatially homogeneous equilibrium, the Turing instability and the oscillatory instability are examined depending on the biological characteristics of zooplankton and spatial movement of plankton. It is shown that at low values of zooplankton grazing rate and intratrophic interaction rate the system is Turing instable when the taxis rate is low. Stabilization occurs either through increased decline of zooplankton either by increasing the phytoplankton diffusion. With the increasing rate of consumption of phytoplankton range of parameters that determine the stability is reduced. A type of instability depends on the phytoplankton diffusion. For large values of diffusion oscillatory instability is observed, with a decrease in the phytoplankton diffusion zone of Turing instability is increases. In general, if zooplankton grazing rate is faster than phytoplankton growth rate the spatially homogeneous equilibrium is Turing instable or oscillatory instable. Stability is observed only at high speeds of zooplankton departure or its active movements. With the increase in zooplankton search activity spatial distribution of populations becomes more uniform, increasing the rate of diffusion leads to non-uniform spatial distribution. However, under diffusion the total number of the population is stabilized when the zooplankton grazing rate above the rate of phytoplankton growth. In general, at low rate of phytoplankton consumption the spatial structures formation is possible at low rates of zooplankton decline and diffusion of all the plankton community. With the increase in phytoplankton predation rate the phytoplankton diffusion and zooplankton spatial movement has essential effect on the spatial instability.

It is known that the internal mobility of DNA molecules plays an important role in the functioning of these molecules. This explains the great interest of researchers in studying the internal dynamics of DNA. Complexity, laboriousness and high cost of research in this field stimulate the search and creation of simpler physical analogues, convenient for simulating the various dynamic regimes possible in DNA. One of the directions of such a search is connected with the use of a mechanical analogue of DNA — a chain of coupled pendulums. In this model, pendulums imitate nitrous bases, horizontal thread on which pendulums are suspended, simulates a sugarphosphate chain, and gravitational field simulates a field induced by a second strand of DNA. Simplicity and visibility are the main advantages of the mechanical analogue. However, the model becomes too cumbersome in cases where it is necessary to simulate long (more than a thousand base pairs) DNA sequences. Another direction is associated with the use of an electronic analogue of the DNA molecule, which has no shortcomings of the mechanical model. In this paper, we investigate the possibility of using the Josephson line as an electronic analogue. We calculated the coefficients of the direct and indirect transformations for the simple case of a homogeneous, synthetic DNA, the sequence of which contains only adenines. The internal mobility of the DNA molecule was modeled by the sine-Gordon equation for angular vibrations of nitrous bases belonging to one of the two polynucleotide chains of DNA. The second polynucleotide chain was modeled as a certain average field in which these oscillations occur. We obtained the transformation, allowing the transition from DNA to an electronic analog in two ways. The first includes two stages: (1) the transition from DNA to the mechanical analogue (a chain of coupled pendulums) and (2) the transition from the mechanical analogue to the electronic one (the Josephson line). The second way is direct. It includes only one stage — a direct transition from DNA to the electronic analogue.

In this paper, the system of equations of magnetic hydrodynamics (MHD) is considered. The exact solutions found describe fluid flows in a porous medium and are related to the development of a core simulator and are aimed at creating a domestic technology «digital deposit» and the tasks of controlling the parameters of incompressible fluid. The central problem associated with the use of computer technology is large-dimensional grid approximations and high-performance supercomputers with a large number of parallel microprocessors. Kinetic methods for solving differential equations and methods for «gluing» exact solutions on coarse grids are being developed as possible alternatives to large-dimensional grid approximations. A comparative analysis of the efficiency of computing systems allows us to conclude that it is necessary to develop the organization of calculations based on integer arithmetic in combination with universal approximate methods. A class of exact solutions of the Navier – Stokes system is proposed, describing three-dimensional flows for an incompressible fluid, as well as exact solutions of nonstationary three-dimensional magnetic hydrodynamics. These solutions are important for practical problems of controlled dynamics of mineralized fluids, as well as for creating test libraries for verification of approximate methods. A number of phenomena associated with the formation of macroscopic structures due to the high intensity of interaction of elements of spatially homogeneous systems, as well as their occurrence due to linear spatial transfer in spatially inhomogeneous systems, are highlighted. It is fundamental that the emergence of structures is a consequence of the discontinuity of operators in the norms of conservation laws. The most developed and universal is the theory of computational methods for linear problems. Therefore, from this point of view, the procedures of «immersion» of nonlinear problems into general linear classes by changing the initial dimension of the description and expanding the functional spaces are important. Identification of functional solutions with functions makes it possible to calculate integral averages of an unknown, but at the same time its nonlinear superpositions, generally speaking, are not weak limits of nonlinear superpositions of approximations of the method, i.e. there are functional solutions that are not generalized in the sense of S. L. Sobolev.

Population dynamics is considered in a modified chemostat model including diffusion, chemotaxis, and nonlocal competitive losses. To account for influence of the external environment on the population of the ecosystem, a random parameter is included into the model equations. Computer simulations reveal three dynamic modes depending on system parameters: the transition from initial state to a spatially homogeneous steady state, to a spatially inhomogeneous distribution of population density, and elimination of population density.

A model, describing the influence of external factors on temporal evolution of phytoplankton distribution in a horizontally-homogenous water layer, is presented. This model is based upon the reactiondiffusion equation and takes into account the main factors of influence: mineral nutrients, insolation and temperature. The mineral nutrients and insolation act oppositely on spatial phytoplankton distribution. The results of numerical modeling are presented and the prospect of applying this model to reconstruction of phytoplankton distribution from sea-surface satellite data is discussed. The model was used to estimate the chlorophyll content of the Peter the Great Bay (Sea of Japan).

To study nonlinear effects of biological species interactions numerical-analytical approach is being developed. The approach is based on the cosymmetry theory accounting for the phenomenon of the emergence of a continuous family of solutions to differential equations where each solution can be obtained from the appropriate initial state. In problems of mathematical ecology the onset of cosymmetry is usually connected with a number of relationships between the parameters of the system. When the relationships collapse families vanish, we get a finite number of isolated solutions instead of a continuum of solutions and transient process can be long-term, dynamics taking place in a neighborhood of a family that has vanished due to cosymmetry collapse.

We consider a model for spatiotemporal competition of predators or prey with an account for directed migration, Holling type II functional response and nonlinear prey growth function permitting Alley effect. We found out the conditions on system parameters under which there is linear with respect to population densities cosymmetry. It is demonstated that cosymmetry exists for any resource function in case of heterogeneous habitat. Numerical experiment in MATLAB is applied to compute steady states and oscillatory regimes in case of spatial heterogeneity.

The dynamics of three population interactions (two predators and a prey, two prey and a predator) are considered. The onset of families of stationary distributions and limit cycle branching out of equlibria of a family that lose stability are investigated in case of homogeneous habitat. The study of the system for two prey and a predator gave a wonderful result of species coexistence. We have found out parameter regions where three families of stable solutions can be realized: coexistence of two prey in absence of a predator, stationary and oscillatory distributions of three coexisting species. Cosymmetry collapse is analyzed and long-term transient dynamics leading to solutions with the exclusion of one of prey or extinction of a predator is established in the numerical experiment.

Our purpose is to study the spatio-temporal population wave behavior observed in the predator-prey system. It is assumed that predators move both directionally and randomly, and prey spread only diffusely. The model does not take into account demographic processes in the predator population; it’s total number is constant and is a parameter. The variables of the model are the prey and predator densities and the predator speed, which are connected by a system of three reaction – diffusion – advection equations. The system is considered on an annular range, that is the periodic conditions are set at the boundaries of the interval. We have studied the bifurcations of wave modes arising in the system when two parameters are changed — the total number of predators and their taxis acceleration coefficient.

The main research method is a numerical analysis. The spatial approximation of the problem in partial derivatives is performed by the finite difference method. Integration of the obtained system of ordinary differential equations in time is carried out by the Runge –Kutta method. The construction of the Poincare map, calculation of Lyapunov exponents, and Fourier analysis are used for a qualitative analysis of dynamic regimes.

It is shown that, population waves can arise as a result of existence of directional movement of predators. The population dynamics in the system changes qualitatively as the total predator number increases. А stationary homogeneous regime is stable at low value of parameter, then it is replaced by self-oscillations in the form of traveling waves. The waveform becomes more complicated as the bifurcation parameter increases; its complexity occurs due to an increase in the number of temporal vibrational modes. A large taxis acceleration coefficient leads to the possibility of a transition from multi-frequency to chaotic and hyperchaotic population waves. A stationary regime without preys becomes stable with a large number of predators.

In this work, the problem of constructing an electronic analogue of heterogeneous DNA is solved with the help of the methods of mathematical modeling. Electronic analogs of that type, along with other physical models of living systems, are widely used as a tool for studying the dynamic and functional properties of these systems. The solution to the problem is based on an algorithm previously developed for homogeneous (synthetic) DNA and modified in such a way that it can be used for the case of inhomogeneous (native) DNA. The algorithm includes the following steps: selection of a model that simulates the internal mobility of DNA; construction of a transformation that allows you to move from the DNA model to its electronic analogue; search for conditions that provide an analogy of DNA equations and electronic analogue equations; calculation of the parameters of the equivalent electrical circuit. To describe inhomogeneous DNA, the model was chosen that is a system of discrete nonlinear differential equations simulating the angular deviations of nitrogenous bases, and Hamiltonian corresponding to these equations. The values of the coefficients in the model equations are completely determined by the dynamic parameters of the DNA molecule, including the moments of inertia of nitrous bases, the rigidity of the sugar-phosphate chain, and the constants characterizing the interactions between complementary bases in pairs. The inhomogeneous Josephson line was used as a basis for constructing an electronic model, the equivalent circuit of which contains four types of cells: A-, T-, G-, and C-cells. Each cell, in turn, consists of three elements: capacitance, inductance, and Josephson junction. It is important that the A-, T-, G- and C-cells of the Josephson line are arranged in a specific order, which is similar to the order of the nitrogenous bases (A, T, G and C) in the DNA sequence. The transition from DNA to an electronic analog was carried out with the help of the A-transformation which made it possible to calculate the values of the capacitance, inductance, and Josephson junction in the A-cells. The parameter values for the T-, G-, and C-cells of the equivalent electrical circuit were obtained from the conditions imposed on the coefficients of the model equations and providing an analogy between DNA and the electronic model.