BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael Margaliot

Michael Margaliot appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.DSmath.OCGenomicsmath.CAQuantitative MethodsSystems and ControlMolecular NetworksSubcellular ProcessesComputational ComplexityData Structures and Algorithmsmath.CONetworking and Internet Architecture

Catalog footprint

What is connected

33works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2025arXiv

Instability of equilibrium and convergence to periodic orbits in strongly 2-cooperative systems

We consider time-invariant nonlinear $n$-dimensional strongly $2$-cooperative systems, that is, systems that map the set of vectors with up to weak sign variation to its interior. Strongly $2$-cooperative systems enjoy a strong Poincare-Bendixson property: bounded solutions that maintain a positive distance from the set of equilibria converge to a periodic solution. For strongly $2$-cooperative systems whose trajectories evolve in a bounded and invariant set that contains a single unstable equilibrium, we provide a simple criterion for the existence of periodic trajectories. Moreover, we explicitly characterize a positive-measure set of initial conditions which yield solutions that asymptotically converge to a periodic trajectory. We demonstrate our theoretical results using two models from systems biology, the $n$-dimensional Goodwin oscillator and a $4$-dimensional biomolecular oscillator with RNA-mediated regulation, and provide numerical simulations that verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Negative feedback and oscillations in a model for mRNA translation

The ribosome flow model (RFM) is a phenomenological model for the unidirectional flow of particles along a 1D chain of $n$ sites. The RFM has been extensively used to study the dynamics of ribosome flow along a single mRNA molecule during translation. In this case, the particles model ribosomes and each site corresponds to a consecutive group of codons. Networks of interconnected RFMs have been used to model and analyze large-scale translation in the cell and, in particular, the effects of competition for shared resources. Here, we analyze the RFM with a negative feedback connection from the protein production rate to the initiation rate. This models, for example, the production of proteins that inhibit the translation of their own mRNA. Using tools from the theory of 2-cooperative dynamical systems, we provide a simple condition guaranteeing that the closed-loop system admits at least one non-trivial periodic solution. When this condition holds, we also explicitly characterize a large set of initial conditions such that any solution emanating from this set converges to a non-trivial periodic solution. Such a solution corresponds to a periodic pattern of ribosome densities along the mRNA, and to a periodic pattern of protein production.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Multiplicative and additive compounds via Kronecker products and Kronecker sums

Compound matrices play an important role in many fields of mathematics and have recently found new applications in systems and control theory. However, the explicit formulas for these compounds are non-trivial and not always easy to use. Here, we derive new formulas for the multiplicative and additive compounds of a matrix using Kronecker products and sums. This provides a new approach to matrix compounds based on the well-known and powerful theory of Kronecker products and sums. We demonstrate several applications of these new formulas, including deriving a new expression for the additive compound of the product of two matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Compound matrices in systems and control theory: a tutorial

The multiplicative and additive compounds of a matrix play an important role in several fields of mathematics including geometry, multi-linear algebra, combinatorics, and the analysis of nonlinear time-varying dynamical systems. There is a growing interest in applications of these compounds, and their generalizations, in systems and control theory. The goal of this tutorial paper is to provide a gentle and self-contained introduction to these topics with an emphasis on the geometric interpretation of the compounds, and to describe some of their recent applications including several non-trivial generalizations of positive systems, cooperative systems, contracting systems, and more.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Minimum effort decentralized control design for contracting network systems

We consider the problem of making a networked system contracting by designing minimal effort local controllers. Our method combines a hierarchical contraction characterization and a matrix-balancing approach to stabilizing a Metzler matrix via minimal diagonal perturbations. We demonstrate our approach by designing local controllers that render contractive a network of FitzHugh-Nagumo neurons with a general topology of interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On matrices whose exponential is a P-matrix

A matrix is called a P-matrix if all its principal minors are positive. P-matrices have found important applications in functional analysis, mathematical programming, and dynamical systems theory. We introduce a new class of real matrices denoted~$\EP$. A matrix is in~$\EP$ if and only if its matrix exponential is a P-matrix for all positive times. In other words, $A\in \EP$ if and only if the transition matrix of the linear system~$\dot x=Ax$ is a P-matrix for any positive time~$t$. We analyze the properties of this new class of matrices and describe an application of our theoretical results to opinion dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Verifying $k$-Contraction without Computing $k$-Compounds

Compound matrices have found applications in many fields of science including systems and control theory. In particular, a sufficient condition for $k$-contraction is that a logarithmic norm (also called matrix measure) of the $k$-additive compound of the Jacobian is uniformly negative. However, this may be difficult to check in practice because the $k$-additive compound of an $n\times n$ matrix has dimensions $\binom{n}{k}\times \binom{n}{k}$. For an $n\times n$ matrix $A$, we prove a duality relation between the $k$ and $(n-k)$ compounds of $A$. We use this duality relation to derive a sufficient condition for $k$-contraction that does not require the computation of any $k$-compounds. We demonstrate our results by deriving a sufficient condition for $k$-contraction of an $n$-dimensional Hopfield network that does not require to compute any compounds. In particular, for $k=2$ this sufficient condition implies that the network is $2$-contracting and this implies a strong asymptotic property: every bounded solution of the network converges to an equilibrium point, that may not be unique. This is relevant, for example, when using the Hopfield network as an associative memory that stores patterns as equilibrium points of the dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Behavior of Totally Positive Differential Systems Near a Periodic Solution

A time-varying nonlinear dynamical system is called a totally positive differential system (TPDS) if its Jacobian admits a special sign pattern: it is tri-diagonal with positive entries on the super- and sub-diagonals. If the vector field of a TPDS is T-periodic then every bounded trajectory converges to a T-periodic solution. In particular, when the vector field is time-invariant every bounded trajectory of a TPDS converges to an equlbrium. Here, we use the spectral theory of oscillatory matrices to analyze the behavior near a periodic solution of a TPDS. This yields information on the perturbation directions that lead to the fastest and slowest convergence to or divergence from the periodic solution. We demonstrate the theoretical results using a model from systems biology called the ribosome flow model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Diagonal Stability of Discrete-time $k$-Positive linear Systems with Applications to Nonlinear Systems

A linear dynamical system is called $k$-positive if its dynamics maps the set of vectors with up to $k-1$ sign variations to itself. For $k=1$, this reduces to the important class of positive linear systems. Since stable positive linear time-invariant (LTI) systems always admit a diagonal quadratic Lyapunov function, i.e. they are diagonally stable, we may expect that this holds also for stable $k$-positive systems. We show that, in general, this is not the case both in the continuous-time (CT) and discrete-time (DT) case. We then focus on DT $k$-positive linear systems and introduce the new notion of DT $k$-diagonal stability. It is shown that this is a necessary condition for standard DT diagonal stability. We demonstrate an application of this new notion to the analysis of a class of DT nonlinear systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the Diagonal Stability of $k$-Positive Linear Systems

We consider $k$-positive linear systems, that is, systems that map the set of vectors with up to $k-1$ sign variations to itself. For $k=1$, this reduces to positive linear systems. It is well-known that stable positive linear time invariant (LTI) systems admit a diagonal Lyapunov function. This property has many important implications. A natural question is whether stable $k$-positive systemsalso admit a diagonal Lyapunov function. This paper shows that, in general, the answer is no. However, for both continuous-time and discrete-time $n$-dimensional systems that are $(n-1)$-positive we provide a sufficient condition for diagonal stability.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Variability in mRNA Translation: A Random Matrix Theory Approach

The rate of mRNA translation depends on the initiation, elongation, and termination rates of ribosomes along the mRNA. These rates depend on many "local" factors like the abundance of free ribosomes and tRNA molecules in the vicinity of the mRNA molecule. All these factors are stochastic and their experimental measurements are also noisy. An important question is how protein production in the cell is affected by this considerable variability. We develop a new theoretical framework for addressing this question by modeling the rates as identically and independently distributed random variables and using tools from random matrix theory to analyze the steady-state production rate. The analysis reveals a principle of universality: the average protein production rate depends only on the of the set of possible values that the random variable may attain. This explains how total protein production can be stabilized despite the overwhelming stochasticticity underlying cellular processes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Is my system of ODEs $k$-cooperative?

A linear dynamical system is called positive if its flow maps the non-negative orthant to itself. More precisely, it maps the set of vectors with zero sign variations to itself. A linear dynamical system is called $k$-positive if its flow maps the set of vectors with up to $k-1$ sign variations to itself. A nonlinear dynamical system is called $k$-cooperative if its variational system, which is a time-varying linear dynamical system, is $k$-positive. These systems have special asymptotic properties. For example, it was recently shown that strongly $2$-cooperative systems satisfy a strong Poincaré-Bendixson property. Positivity and~$k$-positivity are easy to verify in terms of the sign-pattern of the matrix in the dynamics. However, these sign conditions are not invariant under a coordinate transformation. A natural question is to determine if a given~$n$-dimensional system is $k$-positive up to a coordinate transformation. We study this problem for two special kinds of transformations: permutations and scaling by a signature matrix. For any $n\geq 4$ and~$k\in\{2,\dots, n-2\}$, we provide a graph-theoretical necessary and sufficient condition for $k$-positivity up to such coordinate transformations. We describe an application of our results to a specific class of Lotka-Volterra systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximizing average throughput in oscillatory biological synthesis systems: an optimal control approach

A dynamical system entrains to a periodic input if its state converges globally to an attractor with the same period. In particular, for a constant input the state converges to a unique equilibrium point for any initial condition. We consider the problem of maximizing a weighted average of the system's output along the periodic attractor. The gain of entrainment is the benefit achieved by using a non-constant periodic input relative to a constant input with the same time average. Such a problem amounts to optimal allocation of resources in a periodic manner. We formulate this problem as a periodic optimal control problem which can be analyzed by means of the Pontryagin maximum principle or solved numerically via powerful software packages. We then apply our framework to a class of occupancy models that appear frequently in biological synthesis systems and other applications. We show that, perhaps surprisingly, constant inputs are optimal for various architectures. This suggests that the presence of non-constant periodic signals, which frequently appear in biological occupancy systems, is a signature of an underlying time-varying objective functional being optimized.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Exponent of Several Classes of Oscillatory Matrices

Oscillatory matrices were introduced in the seminal work of Gantmacher and Krein. An $n\times n$ matrix $A$ is called oscillatory if all its minors are nonnegative and there exists a positive integer $k$ such that all minors of $A^k$ are positive. The smallest $k$ for which this holds is called the exponent of the oscillatory matrix $A$. Gantmacher and Krein showed that the exponent is always smaller than or equal to $n-1$. An important and nontrivial problem is to determine the exact value of the exponent. Here we use the successive elementary bidiagonal factorization of oscillatory matrices, and its graph-theoretic representation, to derive an explicit expression for the exponent of several classes of oscillatory matrices, and a nontrivial upper-bound on the exponent for several other classes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Random attraction in the TASEP model

The totally asymmetric simple exclusion process (TASEP) is a basic model of statistical mechanics that has found numerous applications. We consider the case of TASEP with a finite chain where particles may enter from the left and leave to the right at prescribed rates. This model can be formulated as a Markov process with a finite number of states. Due to the irreducibility of the process it is well-known that the probability distribution on the states is globally attracted to a unique equilibrium distribution. We extend this result to the more detailed level of individual trajectories. To do so we formulate TASEP as a random dynamical system. Our main result is that the trajectories from all possible initial conditions contract to each other yielding the existence of a random attractor that consists of a single trajectory almost surely. This implies that in the long run TASEP ``filters out'' any perturbation that changes the state of the particles along the chain. In order to prove our main result we first establish that any random dynamical system on a finite state space possesses both a global random pullback attractor and a global random forward attractor. This observation appears to be missing in the literature. We then provide sufficient and necessary conditions for these attractors to be singletons. Finally, we show that TASEP satisfies one of these conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A Generalization of Smillie's Theorem on Strongly Cooperative Tridiagonal Systems

Smillie (1984) proved an interesting result on the stability of nonlinear, time-invariant, strongly cooperative, and tridiagonal dynamical systems. This result has found many applications in models from various fields including biology, ecology, and chemistry. Smith (1991) has extended Smillie's result and proved entrainment in the case where the vector field is time-varying and periodic. We use the theory of linear totally nonnegative differential systems developed by Schwarz (1970) to give a generalization of these two results. This is based on weakening the requirement for strong cooperativity to cooperativity, and adding an additional observability-type condition.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Output Selection and Observer Design for Boolean Control Networks: A Sub-Optimal Polynomial-Complexity Algorithm

Using a graph-theoretic approach, we derive a new sufficient condition for observability of a Boolean control network (BCN). Based on this condition, we describe two algorithms: the first selects a set of nodes so that observing this set makes the BCN observable. The second algorithm builds an observer for the observable BCN. Both algorithms are sub-optimal, as they are based on a sufficient but not necessary condition for observability. Yet their time-complexity is linear in the length of the description of the BCN, rendering them feasible for large-scale networks. We discuss how these results can be used to provide a sub-optimal yet polynomial-complexity algorithm for the minimal observability problem in BCNs. Some of the theoretical results are demonstrated using a BCN model of the core network regulating the mammalian cell cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

A Polynomial-Time Algorithm for Solving the Minimal Observability Problem in Conjunctive Boolean Networks

Many complex systems in biology, physics, and engineering include a large number of state-variables, and measuring the full state of the system is often impossible. Typically, a set of sensors is used to measure part of the state-variables. A system is called observable if these measurements allow to reconstruct the entire state of the system. When the system is not observable, an important and practical problem is how to add a \emph{minimal} number of sensors so that the system becomes observable. This minimal observability problem is practically useful and theoretically interesting, as it pinpoints the most informative nodes in the system. We consider the minimal observability problem for an important special class of Boolean networks, called conjunctive Boolean networks (CBNs). Using a graph-theoretic approach, we provide a necessary and sufficient condition for observability of a CBN with $n$ state-variables, and an efficient~$O(n^2)$-time algorithm for solving the minimal observability problem. We demonstrate the usefulness of these results by studying the properties of a class of random CBNs.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Minimal Controllability of Conjunctive Boolean Networks is NP-Complete

Given a conjunctive Boolean network (CBN) with $n$ state-variables, we consider the problem of finding a minimal set of state-variables to directly affect with an input so that the resulting conjunctive Boolean control network (CBCN) is controllable. We give a necessary and sufficient condition for controllability of a CBCN; an $O(n^2)$-time algorithm for testing controllability; and prove that nonetheless the minimal controllability problem for CBNs is NP-hard.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimal Down Regulation of mRNA Translation

Down regulation of mRNA translation is an important problem in various bio-medical domains ranging from developing effective medicines for tumors and for viral diseases to developing attenuated virus strains that can be used for vaccination. Here, we study the problem of down regulation of mRNA translation using a mathematical model called the ribosome flow model (RFM). In the RFM, the mRNA molecule is modeled as a chain of $n$ sites. The flow of ribosomes between consecutive sites is regulated by $n+1$ transition rates. Given a set of feasible transition rates, that models the outcome of all possible mutations, we consider the problem of maximally down regulating the translation rate by altering the rates within this set of feasible rates. Under certain conditions on the feasible set, we show that an optimal solution can be determined efficiently. We also rigorously analyze two special cases of the down regulation optimization problem. Our results suggest that one must focus on the position along the mRNA molecule where the transition rate has the strongest effect on the protein production rate. However, this rate is not necessarily the slowest transition rate along the mRNA molecule. We discuss some of the biological implications of these results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Model for Competition for Ribosomes in the Cell

Large-scale simultaneous mRNA translation and the resulting competition for the available ribosomes has important implications to the cell's functioning and evolution. Developing a better understanding of the intricate correlations between these simultaneous processes, rather than focusing on the translation of a single isolated transcript, should help in gaining a better understanding of mRNA translation regulation and the way elongation rates affect organismal fitness. A model of simultaneous translation is specifically important when dealing with highly expressed genes, as these consume more resources. In addition, such a model can lead to more accurate predictions that are needed in the interconnection of translational modules in synthetic biology. We develop and analyze a general model for large-scale simultaneous mRNA translation and competition for ribosomes. This is based on combining several ribosome flow models (RFMs) interconnected via a pool of free ribosomes. We prove that the compound system always converges to a steady-state and that it always entrains or phase locks to periodically time-varying transition rates in any of the mRNA molecules. We use this model to explore the interactions between the various mRNA molecules and ribosomes at steady-state. We show that increasing the length of an mRNA molecule decreases the production rate of all the mRNAs. Increasing any of the codon translation rates in a specific mRNA molecule yields a local effect: an increase in the translation rate of this mRNA, and also a global effect: the translation rates in the other mRNA molecules all increase or all decrease. These results suggest that the effect of codon decoding rates of endogenous and heterologous mRNAs on protein production is more complicated than previously thought.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analyzing Linear Communication Networks using the Ribosome Flow Model

The Ribosome Flow Model (RFM) describes the unidirectional movement of interacting particles along a one-dimensional chain of sites. As a site becomes fuller, the effective entry rate into this site decreases. The RFM has been used to model and analyze mRNA translation, a biological process in which ribosomes (the particles) move along the mRNA molecule (the chain), and decode the genetic information into proteins. Here we propose the RFM as an analytical framework for modeling and analyzing linear communication networks. In this context, the moving particles are data-packets, the chain of sites is a one dimensional set of ordered buffers, and the decreasing entry rate to a fuller buffer represents a kind of decentralized backpressure flow control. For an RFM with homogeneous link capacities, we provide closed-form expressions for important network metrics including the throughput and end-to-end delay. We use these results to analyze the hop length and the transmission probability (in a contention access mode) that minimize the end-to-end delay in a multihop linear network, and provide closed-form expressions for the optimal parameter values.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Contraction After Small Transients

Contraction theory is a powerful tool for proving asymptotic properties of nonlinear dynamical systems including convergence to an attractor and entrainment to a periodic excitation. We consider three generalizations of contraction with respect to a norm that allow contraction to take place after small transients in time and/or amplitude. These generalized contractive systems (GCSs) are useful for several reasons. First, we show that there exist simple and checkable conditions guaranteeing that a system is a GCS, and demonstrate their usefulness using several models from systems biology. Second, allowing small transients does not destroy the important asymptotic properties of contractive systems like convergence to a unique equilibrium point, if it exists, and entrainment to a periodic excitation. Third, in some cases as we change the parameters in a contractive system it becomes a GCS just before it looses contractivity with respect to a norm. In this respect, generalized contractivity is the analogue of marginal stability in Lyapunov stability theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ribosome Flow Model on a Ring

The asymmetric simple exclusion process (ASEP) is an important model from statistical physics describing particles that hop randomly from one site to the next along an ordered lattice of sites, but only if the next site is empty. ASEP has been used to model and analyze numerous multiagent systems with local interactions including the flow of ribosomes along the mRNA strand. In ASEP with periodic boundary conditions a particle that hops from the last site returns to the first one. The mean field approximation of this model is referred to as the ribosome flow model on a ring (RFMR). The RFMR may be used to model both synthetic and endogenous gene expression regimes. We analyze the RFMR using the theory of monotone dynamical systems. We show that it admits a continuum of equilibrium points and that every trajectory converges to an equilibrium point. Furthermore, we show that it entrains to periodic transition rates between the sites. We describe the implications of the analysis results to understanding and engineering cyclic mRNA translation in-vitro and in-vivo.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Nonlinear Consensus Algorithm Derived from Statistical Physics

The asymmetric simple exclusion process (ASEP) is an important model from statistical physics describing particles that hop randomly from one site to the next along an ordered lattice of sites, but only if the next site is empty. ASEP has been used to model and analyze numerous multiagent systems with local interactions ranging from ribosome flow along the mRNA to pedestrian traffic. In ASEP with periodic boundary conditions a particle that hops from the last site returns to the first one. The mean field approximation of this model is referred to as the ribosome flow model on a ring (RFMR). We analyze the RFMR using the theory of monotone dynamical systems. We show that it admits a continuum of equilibrium points and that every trajectory converges to an equilibrium point. Furthermore, we show that it entrains to periodic transition rates between the sites. When all the transition rates are equal all the state variables converge to the same value. Thus, the RFMR with homogeneous transition rates is a nonlinear consensus algorithm. We describe an application of this to a simple formation control problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-order maximum principles for the stability analysis of positive bilinear control systems

We consider a continuous-time positive bilinear control system (PBCS), i.e. a bilinear control system with Metzler matrices. The positive orthant is an invariant set of such a system, and the corresponding transition matrix C(t) is entrywise nonnegative for all time t>0. Motivated by the stability analysis of positive linear switched systems (PLSSs) under arbitrary switching laws, we fix a final time T>0 and define a control as optimal if it maximizes the spectral radius of C(T). A recent paper developed a first-order necessary condition for optimality in the form of a maximum principle (MP). In this paper, we derive higher-order necessary conditions for optimality for both singular and bang-bang controls. Our approach is based on combining results on the second-order derivative of the spectral radius of a nonnegative matrix with the generalized Legendre-Clebsch condition and the Agrachev-Gamkrelidze second-order optimality condition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Maximizing Protein Translation Rate in the Nonhomogeneous Ribosome Flow Model: A Convex Optimization Approach

Translation is an important stage in gene expression. During this stage, macro-molecules called ribosomes travel along the mRNA strand linking amino-acids together in a specific order to create a functioning protein. An important question is how to maximize protein production. Indeed, translation is known to consume most of the cell's energy and it is natural to assume that evolution shaped this process so that it maximizes the protein production rate. If this is indeed so then one can estimate various parameters of the translation machinery by solving an appropriate mathematical optimization problem. The same problem also arises in the context of synthetic biology, namely, re-engineer heterologous genes in order to maximize their translation rate in a host organism. We consider the problem of maximizing the protein production rate using a computational model for translation-elongation called the ribosome flow model (RFM). This model describes the flow of the ribosomes along an mRNA chain of length n using a set of n first-order nonlinear ODEs. It also includes n+1 positive parameters: the ribosomal initiation rate into the mRNA chain, and n elongation rates along the chain sites. We show that the steady-state translation rate in the RFM is a strictly concave function of its parameters. This means that the problem of maximizing the translation rate under a suitable constraint always admits a unique solution, and that this solution can be determined using highly-efficient algorithms for solving convex optimization problems even for large values of n. Furthermore, our analysis shows that the optimal translation rate can be computed based only on the optimal initiation rate and the elongation rate of the codons near the beginning of the ORF. We discuss some applications of the theoretical results to synthetic biology, molecular evolution, and functional genomics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Maximizing Protein Translation Rate in the Ribosome Flow Model: the Homogeneous Case

Gene translation is the process in which intracellular macro-molecules, called ribosomes, decode genetic information in the mRNA chain into the corresponding proteins. Gene translation includes several steps. During the elongation step, ribosomes move along the mRNA in a sequential manner and link amino-acids together in the corresponding order to produce the proteins. The homogeneous ribosome flow model(HRFM) is a deterministic computational model for translation-elongation under the assumption of constant elongation rates along the mRNA chain. The HRFM is described by a set of n first-order nonlinear ordinary differential equations, where n represents the number of sites along the mRNA chain. The HRFM also includes two positive parameters: ribosomal initiation rate and the (constant) elongation rate. In this paper, we show that the steady-state translation rate in the HRFM is a concave function of its parameters. This means that the problem of determining the parameter values that maximize the translation rate is relatively simple. Our results may contribute to a better understanding of the mechanisms and evolution of translation-elongation. We demonstrate this by using the theoretical results to estimate the initiation rate in M. musculus embryonic stem cell. The underlying assumption is that evolution optimized the translation mechanism. For the infinite-dimensional HRFM, we derive a closed-form solution to the problem of determining the initiation and transition rates that maximize the protein translation rate. We show that these expressions provide good approximations for the optimal values in the n-dimensional HRFM already for relatively small values of n. These results may have applications for synthetic biology where an important problem is to re-engineer genomic systems in order to maximize the protein production rate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Boolean Control Networks with Maximal Topological Entropy

Boolean control networks (BCNs) are discrete-time dynamical systems with Boolean state-variables and inputs that are interconnected via Boolean functions. BCNs are recently attracting considerable interest as computational models for genetic and cellular networks with exogenous inputs. The topological entropy of a BCN with m inputs is a nonnegative real number in the interval [0,m*log 2]. Roughly speaking, a larger topological entropy means that asymptotically the control is "more powerful". We derive a necessary and sufficient condition for a BCN to have the maximal possible topological entropy. Our condition is stated in the framework of Cheng's algebraic state-space representation of BCNs. This means that verifying this condition incurs an exponential time-complexity. We also show that the problem of determining whether a BCN with n state variables and m=n inputs has a maximum topological entropy is NP-hard, suggesting that this problem cannot be solved in general using a polynomial-time algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Structure of the Initiation and Elongation Rates that Maximize Protein Production in the Ribosome Flow Model

Translation is a crucial step in gene expression. During translation, macromolecules called ribosomes "read" the mRNA strand in a sequential manner and produce a corresponding protein. Translation is known to consume most of the cell's energy. Maximizing the protein production rate in mRNA translation, subject to the bounded biomolecullar budget, is thus an important problem in both biology and biotechnology. We consider this problem using a mathematical model for mRNA translation called the ribosome flow model (RFM). For an mRNA strand with $n$ sites the RFM includes $n$ state-variables that encode the normalized ribosomal density at each site, and $n+1$ positive parameters: the initiation rate and elongation rates along the chain. An affine constraint on these rates is used to model the bounded cellular budget. We show that for a homogeneous constraint the rates that maximize the steady-state protein production rate have a special structure. They are symmetric with respect to the middle of the chain, and monotonically increase as we move towards the center of the chain. The ribosomal densities corresponding to the optimal rates monotonically decrease along the chain. We discuss some of the biological implications of these results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Three Generalizations of Contraction

We introduce three forms of generalized contraction (GC). Roughly speaking, these are motivated by allowing contraction to take place after small transients in time and/or amplitude. Indeed, contraction is usually used to prove asymptotic properties, like convergence to an attractor or entrainment to a periodic excitation, and allowing initial transients does not affect this asymptotic behavior. We provide sufficient conditions for GC, and demonstrate their usefulness using examples of systems that are not contractive, with respect to any norm, yet are GC.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sensitivity of mRNA Translation

Using the dynamic mean-field approximation of the totally asymmetric simple exclusion process (TASEP), we investigate the effect of small changes in the initiation, exit, and elongation rates along the mRNA strand on the steady state protein translation rate. We focus on two special cases where exact closed-form expressions for the translation rate sensitivity can be derived. We discuss the ramifications of our results in the context of functional genomics, molecular evolution, and synthetic biology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Switching Between Linear Consensus~Protocols: A Variational Approach

We consider a linear consensus system with n agents that can switch between r different connectivity patterns. A natural question is which switching law yields the best (or worst) possible speed of convergence to consensus? We formulate this question in a rigorous manner by relaxing the switched system into a bilinear consensus control system, with the control playing the role of the switching law. A best (or worst) possible switching law then corresponds to an optimal control. We derive a necessary condition for optimality, stated in the form of a maximum principle (MP). Our approach, combined with suitable algorithms for numerically solving optimal control problems, may be used to obtain explicit lower and upper bounds on the achievable rate of convergence to consensus. We also show that the system will converge to consensus for any switching law if and only if a certain (n-1) dimensional linear switched system converges to the origin for any switching law. For the case n=3 and r=2, this yields a necessary and sufficient condition for convergence to consensus that admits a simple graph-theoretic interpretation.

0 citations0 reviewsNo code linkedOpen access