Search PubMedSearch

PubMed · 2453282

Stoichiometric network analysis.

Abstract

Stoichiometric network analysis is a systematic, general approach to the qualitative, nonlinear dynamics of chemical reaction mechanisms and other systems with stoichiometry. The advantage of a qualitative approach is that no rate constants are needed to determine qualitative features of the dynamics. If one is interested in stability, the approach yields inequalities among the steady-state concentrations and the rate of flow through sequences of important reactions. These parameters are often the ones most easily measured experimentally. By comparing such experiments with the inequalities derived from stoichiometric network analysis, one can often prove that certain mechanisms cannot account for oscillations or other types of observed dynamics. The approach covers far more than stability. The existence of steady states of zero concentration has an interesting mathematics and applies to chemical evolution. The folding of the manifold of steady states can be found by direct calculation and plays a role in switching enzymes on and off. The approach leads to theorems showing that some steady states are globally attracting or, possibly, that a region containing chaos or an oscillation is globally attracting. The subject of sensitivity analysis has been reformulated in this context. Algorithms that apply many of the theoretical results to chemical networks have been developed and combined into a computer program package.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

B L Clarke. Stoichiometric network analysis.. https://doi.org/10.1007/bf02918360

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Preferential exclusion of sucrose from recombinant interleukin-1 receptor antagonist: role in restricted conformational mobility and compaction of native state.

Understanding the mechanism for sucrose-induced protein stabilization is important in many diverse fields, ranging from biochemistry and environmental physiology to pharmaceutical science. Timasheff and Lee [Lee, J. C. & Timasheff, S. N. (1981) J. Biol. Chem. 256, 7193-7201] have established that thermodynamic stabilization of proteins by sucrose is due to preferential exclusion of the sugar from the protein's surface, which increases protein chemical potential. The current study measures the preferential exclusion of 1 M sucrose from a protein drug, recombinant interleukin 1 receptor antagonist (rhIL-1ra). It is proposed that the degree of preferential exclusion and increase in chemical potential are directly proportional to the protein surface area and that, hence, the system will favor the protein state with the smallest surface area. This mechanism explains the observed sucrose-induced restriction of rhIL-1ra conformational fluctuations, which were studied by hydrogen-deuterium exchange and cysteine reactivity measurements. Furthermore, infrared spectroscopy of rhlL-1ra suggested that a more ordered native conformation is induced by sucrose. Electron paramagnetic resonance spectroscopy demonstrated that in the presence of sucrose, spin-labeled cysteine 116 becomes more buried in the protein's interior and that the hydrodynamic diameter of the protein is reduced. The preferential exclusion of sucrose from the protein and the resulting shift in the equilibrium between protein states toward the most compact conformation account for sucrose-induced effects on rhIL-1ra.

Chemical Phenomena

Comparison of hard-cylinder and screened Coulomb interactions in the modeling of supercoiled DNAs.

A 1000 base pair (bp) model supercoiled DNA is simulated using spherical screened Coulomb interactions between subunits on one hand and equivalent hard-cylinder interactions on the other. The amplitudes, or effective charges, of the spherical screened Coulomb electrostatic potentials are chosen so that the electrostatic potential surrounding the middle of a linear array of 2001 subunits (31.8 A diameter) closely matches the solution of the nonlinear Poisson-Boltzmann equation for a cylinder with 12 A radius and the full linear charge density of DNA at all distances beyond the 24 A hard-core diameter. This superposition of spherical screened Coulomb potentials is practically identical to the particular solution of the cylindrical linearized Poisson-Boltzmann equation that matches the solution of the nonlinear Poisson-Boltzmann equation at large distances. The interaction energy between subunits is reckoned from the effective charges according to the standard DLVO expression. The equivalent hard-cylinder diameter is chosen following Stigter's protocol for matching second virial coefficients, but for the full linear charge density of DNA. The electrostatic persistence length of the model with screened Coulomb interactions is extremely sensitive to the (arbitrarily) chosen subunit length at the higher salt concentrations. The persistence length of the hard-cylinder model is adjusted to match that of the screened Coulomb model for each ionic condition. Simulations for a superhelix density sigma = -0.05 using a spherical screened Coulomb interaction plus a 24 A hard-cylinder core (SCPHC) potential indicate that the radius of gyration of this 1000 bp DNA actually undergoes a slight increase as the NaCl concentration is raised from 0.01 to 1.0M. Thus, merely softening the potential from hard-cylinder to screened Coulomb form does not produce a large decrease in radius of gyration with increasing NaCl concentration for DNAs of this size. Radii of gyration, static structure factors, and diffusion coefficients obtained using the equivalent hard-cylinder (EHC) potential agree well with those obtained using the SCPHC potential in 1.0M NaCl, but in 0.1M NaCl the agreement is not as good, and in 0.01M NaCl the agreement is definitely unsatisfactory. These conclusions differ in significant respects from those obtained in previous studies.

Chemical Phenomena