Moiety conservation

Last updated

Moiety conservation is the conservation of a subgroup in a chemical species, which is cyclically transferred from one molecule to another. In biochemistry, moiety conservation can have profound effects on the system's dynamics. [1]

Contents

Moiety-conserved cycles in biochemistry

A typical example of a conserved moiety [2] in biochemistry is the Adenosine diphosphate (ADP) subgroup that remains unchanged when it is phosphorylated to create adenosine triphosphate (ATP) and then dephosphorylated back to ADP forming a conserved cycle. Moiety-conserved cycles in nature exhibit unique network control features which can be elucidated using techniques such as metabolic control analysis. Other examples in metabolism include NAD/NADH, NADP/NADPH, CoA/Acetyl-CoA. Conserved cycles also exist in large numbers in protein signaling networks when proteins get phosphorylated and dephosphorylated.

Phosphorylation Cycle. Moiety A is conserved. PhosphorylationCycle.png
Phosphorylation Cycle. Moiety A is conserved.

Most, if not all, of these cycles, are time-scale-dependent. For example, although a protein in a phosphorylation cycle is conserved during the interconversion, over a longer time scale, there will be low levels of protein synthesis and degradation, which change the level of protein moiety. The same applies to cycles involving ATP, NAD, etc. Thus, although the concept of a moiety-conserved cycle in biochemistry is a useful approximation, [3] over time scales that include significant net synthesis and degradation of the moiety, the approximation is no longer valid. When invoking the conserved-moiety assumption on a particular moiety, we are, in effect, assuming the system is closed to that moiety.

Identifying conserved cycles

Conserved cycles in a biochemical network can be identified by examination of the stoichiometry matrix, [4] [5] [6] . The stoichiometry matrix for a simple cycle with species A and AP is given by:

The rates of change of A and AP can be written using the equation:

Expanding the expression leads to:

Note that . This means that , where is the total mass of moiety .

Given an arbitrary system:

elementary row operations can be applied to both sides such that the stoichiometric matrix is reduced to its echelon form, giving:

The elementary operations are captured in the matrix. We can partition to match the echelon matrix where the zero rows begin such that:

By multiplying out the lower partition, we obtain:

The matrix will contain entries corresponding to the conserved cycle participants.

Conserved cycles and computer models

The presence of conserved moieties can affect how computer simulation models are constructed. [7] [8] Moiety-conserved cycles will reduce the number of differential equations required to solve a system. For example, a simple cycle has only one independent variable. The other variable can be computed using the difference between the total mass and the independent variable. The set of differential equations for the two-cycle is given by:

These can be reduced to one differential equation and one linear algebraic equation:

Related Research Articles

<span class="mw-page-title-main">Spherical coordinate system</span> Coordinates comprising a distance and two angles

In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a given point in space is specified by three numbers, : the radial distance of the radial liner connecting the point to the fixed point of origin ; the polar angle θ of the radial line r; and the azimuthal angle φ of the radial line r.

Kinematics is a subfield of physics and mathematics, developed in classical mechanics, that describes the motion of points, bodies (objects), and systems of bodies without considering the forces that cause them to move. Kinematics, as a field of study, is often referred to as the "geometry of motion" and is occasionally seen as a branch of both applied and pure mathematics since it can be studied without considering the mass of a body or the forces acting upon it. A kinematics problem begins by describing the geometry of the system and declaring the initial conditions of any known values of position, velocity and/or acceleration of points within the system. Then, using arguments from geometry, the position, velocity and acceleration of any unknown parts of the system can be determined. The study of how forces act on bodies falls within kinetics, not kinematics. For further details, see analytical dynamics.

The vorticity equation of fluid dynamics describes the evolution of the vorticity ω of a particle of a fluid as it moves with its flow; that is, the local rotation of the fluid. The governing equation is:

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition.

<span class="mw-page-title-main">Hooke's law</span> Physical law: force needed to deform a spring scales linearly with distance

In physics, Hooke's law is an empirical law which states that the force needed to extend or compress a spring by some distance scales linearly with respect to that distance—that is, Fs = kx, where k is a constant factor characteristic of the spring, and x is small compared to the total possible deformation of the spring. The law is named after 17th-century British physicist Robert Hooke. He first stated the law in 1676 as a Latin anagram. He published the solution of his anagram in 1678 as: ut tensio, sic vis. Hooke states in the 1678 work that he was aware of the law since 1660.

In linear algebra, a square matrix  is called diagonalizable or non-defective if it is similar to a diagonal matrix. That is, if there exists an invertible matrix  and a diagonal matrix such that . This is equivalent to . This property exists for any linear map: for a finite-dimensional vector space , a linear map  is called diagonalizable if there exists an ordered basis of  consisting of eigenvectors of . These definitions are equivalent: if  has a matrix representation as above, then the column vectors of  form a basis consisting of eigenvectors of , and the diagonal entries of  are the corresponding eigenvalues of ; with respect to this eigenvector basis,  is represented by .

<span class="mw-page-title-main">Euler equations (fluid dynamics)</span> Set of quasilinear hyperbolic equations governing adiabatic and inviscid flow

In fluid dynamics, the Euler equations are a set of partial differential equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. In particular, they correspond to the Navier–Stokes equations with zero viscosity and zero thermal conductivity.

In mathematics, the matrix exponential is a matrix function on square matrices analogous to the ordinary exponential function. It is used to solve systems of linear differential equations. In the theory of Lie groups, the matrix exponential gives the exponential map between a matrix Lie algebra and the corresponding Lie group.

<span class="mw-page-title-main">Total least squares</span> Statistical technique

In applied statistics, total least squares is a type of errors-in-variables regression, a least squares data modeling technique in which observational errors on both dependent and independent variables are taken into account. It is a generalization of Deming regression and also of orthogonal regression, and can be applied to both linear and non-linear models.

A Hamiltonian system is a dynamical system governed by Hamilton's equations. In physics, this dynamical system describes the evolution of a physical system such as a planetary system or an electron in an electromagnetic field. These systems can be studied in both Hamiltonian mechanics and dynamical systems theory.

In applied mathematics, in particular the context of nonlinear system analysis, a phase plane is a visual display of certain characteristics of certain kinds of differential equations; a coordinate plane with axes being the values of the two state variables, say (x, y), or (q, p) etc. (any pair of variables). It is a two-dimensional case of the general n-dimensional phase space.

In linear algebra, an eigenvector or characteristic vector is a vector that has its direction unchanged by a given linear transformation. More precisely, an eigenvector of a linear transformation is scaled by a constant factor when the linear transformation is applied to it: . It is often important to know these vectors in linear algebra. The corresponding eigenvalue, characteristic value, or characteristic root is the multiplying factor .

<span class="mw-page-title-main">Tissot's indicatrix</span> Characterization of distortion in map projections

In cartography, a Tissot's indicatrix is a mathematical contrivance presented by French mathematician Nicolas Auguste Tissot in 1859 and 1871 in order to characterize local distortions due to map projection. It is the geometry that results from projecting a circle of infinitesimal radius from a curved geometric model, such as a globe, onto a map. Tissot proved that the resulting diagram is an ellipse whose axes indicate the two principal directions along which scale is maximal and minimal at that point on the map.

In control theory, we may need to find out whether or not a system such as is controllable, where , , and are, respectively, , , and matrices for a system with inputs, state variables and outputs.

The angular velocity tensor is a skew-symmetric matrix defined by:

The Cauchy momentum equation is a vector partial differential equation put forth by Cauchy that describes the non-relativistic momentum transport in any continuum.

<span class="mw-page-title-main">Stokes' theorem</span> Theorem in vector calculus

Stokes' theorem, also known as the Kelvin–Stokes theorem after Lord Kelvin and George Stokes, the fundamental theorem for curls or simply the curl theorem, is a theorem in vector calculus on . Given a vector field, the theorem relates the integral of the curl of the vector field over some surface, to the line integral of the vector field around the boundary of the surface. The classical theorem of Stokes can be stated in one sentence: The line integral of a vector field over a loop is equal to the surface integral of its curl over the enclosed surface. It is illustrated in the figure, where the direction of positive circulation of the bounding contour ∂Σ, and the direction n of positive flux through the surface Σ, are related by a right-hand-rule. For the right hand the fingers circulate along ∂Σ and the thumb is directed along n.

In numerical linear algebra, the conjugate gradient method is an iterative method for numerically solving the linear system

In mathematical optimization, the revised simplex method is a variant of George Dantzig's simplex method for linear programming.

The biochemical systems equation is a compact equation of nonlinear differential equations for describing a kinetic model for any network of coupled biochemical reactions and transport processes.

References

  1. Markevich, Nick I.; Hoek, Jan B.; Kholodenko, Boris N. (2 February 2004). "Signaling switches and bistability arising from multisite phosphorylation in protein kinase cascades". Journal of Cell Biology. 164 (3): 353–359. doi: 10.1083/jcb.200308060 . PMC   2172246 . PMID   14744999.
  2. Reich, JG (1981). Energy metabolism of the cell: a theoretical treatise. London: Academic Press.
  3. Fell, D. A. (1 September 1992). "Metabolic control analysis: a survey of its theoretical and experimental development". The Biochemical Journal. 286 ( Pt 2) (Pt 2): 313–330. doi:10.1042/bj2860313. ISSN   0264-6021. PMC   1132899 . PMID   1530563.
  4. Sauro, Herbert M.; Fell, David A. (1991). "SCAMP: A metabolic simulator and control analysis program". Mathematical and Computer Modelling. 15 (12): 15–28. doi: 10.1016/0895-7177(91)90038-9 .
  5. Cornish-Bowden, Athel; Hofmeyr, Jan-Hendrik S. (May 2002). "The Role of Stoichiometric Analysis in Studies of Metabolism: An Example". Journal of Theoretical Biology. 216 (2): 179–191. Bibcode:2002JThBi.216..179C. doi:10.1006/jtbi.2002.2547. PMID   12079370.
  6. Reder, Christine (November 1988). "Metabolic control theory: A structural approach". Journal of Theoretical Biology. 135 (2): 175–201. Bibcode:1988JThBi.135..175R. doi:10.1016/s0022-5193(88)80073-0. PMID   3267767.
  7. Sauro, Herbert M.; Fell, David A. (1991). "SCAMP: A metabolic simulator and control analysis program". Mathematical and Computer Modelling. 15 (12): 15–28. doi: 10.1016/0895-7177(91)90038-9 .
  8. Mendes, Pedro (1993). "GEPASI: a software package for modelling the dynamics, steady states and control of biochemical and other systems". Bioinformatics. 9 (5): 563–571. doi:10.1093/bioinformatics/9.5.563. PMID   8293329.