Systems biology最新文献

Mitochondria in excitable cells are recurrently exposed to pulsatile calcium gradients that activate cell function. Rapid calcium uptake by the mitochondria has previously been shown to cause uncoupling of oxidative phosphorylation. To test (i) if periodic nerve firing may cause oscillation of the cytosolic thermodynamic potential of ATP hydrolysis and (ii) if cytosolic adenylate (AK) and creatine kinase (CK) ATP buffering reactions dampen such oscillations, a lumped kinetic model of an excitable cell capturing major aspects of the physiology has been developed. Activation of ATP metabolism by low-frequency calcium pulses caused large oscillation of the cytosolic, but not mitochondrial ATP/ADP, ratio. This outcome was independent of net ATP synthesis or hydrolysis during mitochondrial calcium uptake. The AK/CK ATP buffering reactions dampened the amplitude and rate of cytosolic ATP/ADP changes on a timescale of seconds, but not milliseconds. These model predictions suggest that alternative sources of capacitance in neurons and striated muscles should be considered to protect ATP-free energy-driven cell functions.

The modelling of molecule-molecule interactions has been widely accepted as a tool for drug discovery and development studies. However, this powerful technique is unappreciated in physiological and biochemical studies, where it could be extremely useful for understanding the mechanisms of action of various compounds in cases when experimental data are controversial due to complexity of the investigated systems. In this study, based on the biochemical data suggesting involvement of mitochondrial ADP/ATP carrier in K+ and H+ transport to mitochondrial matrix molecular modelling is applied to elucidate the possible interactions between the ADP/ATP carrier and its putative ligands--K(ATP) channel blockers glybenclamide, tolbutamide and 5-hydroxydecanoate. Results revealed that K(ATP) channel blockers could bind to the specific location proximal to H1, H4, H5 and H6 transmembrane helices within the cavity of the ADP/ ATP carrier. Analysis of the predicted binding site suggests that K(ATP) channel blockers could interfere with both the ADP/ATP translocation and possible cation flux through the ADP/ATP carrier, and supports the hypothesis that the ADP/ATP carrier is a target of K(ATP) channel modulators.

Metabolic control analysis (MCA) was developed to quantify how system variables are affected by parameter variations in a system. In addition, MCA can express the global properties of a system in terms of the individual catalytic steps, using connectivity and summation theorems to link the control coefficients to the elasticity coefficients. MCA was originally developed for steady-state analysis and not all summation theorems have been derived for dynamic systems. A method to determine time-dependent flux and concentration control coefficients for dynamic systems by expressing the time domain as a function of percentage progression through any arbitrary fixed interval of time is reported. Time-dependent flux and concentration control coefficients of dynamic systems, provided that they are evaluated in this novel way, obey the same summation theorems as steady-state flux and concentration control coefficients, respectively.

Cells adapt to changes in their environment by the concerted action of many different regulatory mechanisms. Examples of such mechanisms are feedback inhibition by intermediates of metabolism, covalent modification of enzymes and changes in the abundance of mRNAs and proteins. These mechanisms act in parallel at different levels in the cellular hierarchy while regulating a single process. Existing hierarchical regulation analysis determines the relative importance of these mechanisms when the cell regulates a transition from one steady-state to another. Here, the analysis is extended to the regulation of time-dependent phenomena, for which two methods are introduced and illustrated with a kinetic model incorporating transcription and translation of metabolic enzymes.

A core model is presented for protein production in Escherichia coli to address the question whether there is an optimal ribosomal concentration for non-ribosome protein production. Analysing the steady-state solution of the model over a range of mRNA concentrations, indicates that such an optimum ribosomal content exists, and that the optimum shifts to higher ribosomal contents at higher specific growth rates.

In the post-genomic era, the biochemical information for individual compounds, enzymes, reactions to be found within named organisms has become readily available. The well-known KEGG and BioCyc databases provide a comprehensive catalogue for this information and have thereby substantially aided the scientific community. Using these databases, the complement of enzymes present in a given organism can be determined and, in principle, used to reconstruct the metabolic network. However, such reconstructed networks contain numerous properties contradicting biological expectation. The metabolic networks for a number of organisms are reconstructed from KEGG and BioCyc databases, and features of these networks are related to properties of their originating database.

An approach for analysis of biological networks is proposed. In this approach, named the connectivity matrix (CM) method, all the connectivities of interest are expressed in a matrix. Then, a variety of analyses are performed on GNU Octave or Matlab. Each node in the network is expressed as a row vector or numeral that carries information defining or characterising the node itself. Information about connectivity itself is also expressed as a row vector or numeral. Thus, connection of node n1 to node n2 through edge e is expressed as [n1, n2, e], a row vector formed by the combination of three row vectors or numerals, where n1, n2 and e indicate two different nodes and one connectivity, respectively. All the connectivities in any given network are expressed as a matrix, CM, each row of which corresponds to one connectivity. Using this CM method, intermetabolite atom-level connectivity is investigated in a model metabolic network composed of the reactions for glycolysis, oxidative decarboxylation of pyruvate, citric acid cycle, pentose phosphate pathway and gluconeogenesis.

Glycine is the major inhibitory neurotransmitter in the brainstem and spinal cord, where it participates in a variety of motor and sensory functions. It activates a special type of ligand-gated membrane receptor, which provides for Cl- ion conductance of the neuronal membrane. Computer simulations of a single-channel current through this receptor have been carried out on the basis of Brownian (Langevin) dynamics. The dependence of the currents on pore diameter and the location of the charged amino acid residues have been obtained. It has been shown that the presence and the symmetry of the filter-forming residues determined not only the ion-selectivity of the channel but also increased transmembrane anion current.

The evaluation of a generic simplified bi-substrate enzyme kinetic equation, whose derivation is based on the assumption of equilibrium binding of substrates and products in random order, is described. This equation is much simpler than the mechanistic (ordered and ping-pong) models, in that it contains fewer parameters (that is, no K(i) values for the substrates and products). The generic equation fits data from both the ordered and the ping-pong models well over a wide range of substrate and product concentrations. In the cases where the fit is not perfect, an improved fit can be obtained by considering the rate equation for only a single set of product concentrations. Due to its relative simplicity in comparison to the mechanistic models, this equation will be useful for modelling bi-substrate reactions in computational systems biology.