Dynamics of granzyme B-induced apoptosis: mathematical modeling

Apoptosis is mediated by an intracellular biochemical system that mainly includes proteins (procaspases, caspases, inhibitors, Bcl-2 protein family as well as substances released from mitochondrial intermembrane space). The dynamics of caspase activation and target cleavage in apoptosis induced by granzyme B in a single K562 cell was studied using a mathematical model of the dynamics of granzyme B-induced apoptosis developed in this work. Also the first application of optimization approach to determination of unknown kinetic constants of biochemical apoptotic reactions was presented. The optimization approach involves solving of two problems: direct and inverse. Solving the direct optimization problem, we obtain the initial (baseline) concentrations of procaspases for known kinetic constants through conditional minimization of a cost function based on the principle of minimum protein consumption by the apoptosis system. The inverse optimization problem is aimed at determination of unknown kinetic constants of apoptotic biochemical reactions proceeding from the condition that the optimal concentrations of procaspases resulting from the solution of the direct optimization problem coincide with the observed ones, that is, those determined by biochemical methods. The Multidimensional Index Method was used to perform numerical solution of the inverse optimization problem.     

Crystal structure of human factor VIII: implications for the formation of the factor IXa-factor VIIIa complex

Factor VIII is a procofactor that plays a critical role in blood coagulation, and is missing or defective in hemophilia A. We determined the X-ray crystal structure of B domain-deleted human factor VIII. This protein is composed of five globular domains and contains one Ca(2+) and two Cu(2+) ions. The three homologous A domains form a triangular heterotrimer where the A1 and A3 domains serve as the base and interact with the C2 and C1 domains, respectively. The structurally homologous C1 and C2 domains reveal membrane binding features. Based on biochemical studies, a model of the factor IXa-factor VIIIa complex was constructed by in silico docking. Factor IXa wraps across the side of factor VIII, and an extended interface spans the factor VIII heavy and light chains. This model provides insight into the activation of factor VIII and the interaction of factor VIIIa with factor IXa on the membrane surface.     

Non-linear dynamics of the complement system activation

The complement system (CS) plays a prominent role in the immune defense. The goal of this work is to study the dynamics of activation of the classic and alternative CS pathways based on the method of mathematical modeling. The principal difficulty that hinders modeling effort is the absence of the measured values of kinetic constants of many biochemical reactions forming the CS. To surmount this difficulty, an optimization procedure consisting of constrained minimization of the total protein consumption by the CS was designed. The constraints made use of published data on the in vitro kinetics of elimination of the Borrelia burgdorferi bacteria by the CS. Special features of the problem at hand called for a significant modification of the general constrained optimization procedure to include a mathematical model of the bactericidal effect of the CS in the iterative setting. Determination of the unknown kinetic constants of biochemical reactions forming the CS led to a fully specified mathematical model of the dynamics of cell killing induced by the CS. On the basis of the model, effects of the initial concentrations of complements and their inhibitors on the bactericidal action of the CS were studied. Proteins playing a critical role in the regulation of the bactericidal action of the CS were identified. Results obtained in this work serve as an important stepping stone for the study of functioning of the CS as a whole as well as for developing methods for control of pathogenic processes.     

Physiological casein and gluten protein requirements of growing rats

Using mounting casein and wheat gluten protein values (0-40%) in the animals' diet, the optimum and minimum physiological daily doses were determined in 49-day-old growing rats from changes in their body water, body nitrogen and protein intake. The optimum physiological doses were identical with the peak of linearity of the given parameters, which coincided with a 15% casein protein and a 20% gluten protein concentration in the diet. This was also confirmed by the maximum body amino acid values, which were found in animals given a 15% casein or 20% gluten protein diet. It was further confirmed by the finding of significantly elevated alanine aminotransferase and aspartate aminotransferase activity in the liver of animals with a higher intake of the above protein sources. The minimum physiological dose of the given protein was determined from the equations of the regression curves in the presence of zero changes in the body nitrogen or body water content. The optimum physiological daily doses of casein and wheat gluten protein were 3.25 g and 4.05 g respectively. The minimum physiological daily doses of casein protein were 268 mg (from body nitrogen changes) and 371 mg (from body water changes) and the minimum physiological daily doses of gluten protein were 892 mg (from body nitrogen changes) and 1,000 mg (from body water changes). The above indicators demonstrate, in the presence of higher and high dietary concentrations, that an intake of the given proteins over and above the optimum physiological daily dose is at the very least uneconomical (gluten), if not harmful (casein), making this a highly topical problem for further study.     

The structural integrity of anion binding exosite I of thrombin is required and sufficient for timely cleavage and activation of factor V and factor VIII

Alpha-thrombin has two separate electropositive binding exosites (anion binding exosite I, ABE-I and anion binding exosite II, ABE-II) that are involved in substrate tethering necessary for efficient catalysis. Alpha-thrombin catalyzes the activation of factor V and factor VIII following discrete proteolytic cleavages. Requirement for both anion binding exosites of the enzyme has been suggested for the activation of both procofactors by alpha-thrombin. We have used plasma-derived alpha-thrombin, beta-thrombin (a thrombin molecule that has only ABE-II available), and a recombinant prothrombin molecule rMZ-II (R155A/R284A/R271A) that can only be cleaved at Arg(320) (resulting in an enzymatically active molecule that has only ABE-I exposed, rMZ-IIa) to ascertain the role of each exosite for procofactor activation. We have also employed a synthetic sulfated pentapeptide (DY(SO(3)(-))DY(SO(3)(-))Q, designated D5Q1,2) as an exosite-directed inhibitor of thrombin. The clotting time obtained with beta-thrombin was increased by approximately 8-fold, whereas rMZ-IIa was 4-fold less efficient in promoting clotting than alpha-thrombin under similar experimental conditions. Alpha-thrombin readily activated factor V following cleavages at Arg(709), Arg(1018), and Arg(1545) and factor VIII following proteolysis at Arg(372), Arg(740), and Arg(1689). Cleavage of both procofactors by alpha-thrombin was significantly inhibited by D5Q1,2. In contrast, beta-thrombin was unable to cleave factor V at Arg(1545) and factor VIII at both Arg(372) and Arg(1689). The former is required for light chain formation and expression of optimum factor Va cofactor activity, whereas the latter two cleavages are a prerequisite for expression of factor VIIIa cofactor activity. Beta-thrombin was found to cleave factor V at Arg(709) and factor VIII at Arg(740), albeit less efficiently than alpha-thrombin. The sulfated pentapeptide inhibited moderately both cleavages by beta-thrombin. Under similar experimental conditions, membrane-bound rMZ-IIa cleaved and activated both procofactor molecules. Activation of the two procofactors by membrane-bound rMZ-IIa was severely impaired by D5Q1,2. Overall the data demonstrate that ABE-I alone of alpha-thrombin can account for the interaction of both procofactors with alpha-thrombin resulting in their timely and efficient activation. Because formation of meizothrombin precedes that of alpha-thrombin, our findings also imply that meizothrombin may be the physiological activator of both procofactors in vivo in the presence of a procoagulant membrane surface during the early stages of coagulation.     

A bipartite autoinhibitory region within the B-domain suppresses function in factor V

Activation of blood coagulation factor V (FV) is a key reaction of hemostasis. FV circulates in plasma as an inactive procofactor, and proteolytic removal of a large central B-domain converts it to an active cofactor (FVa) for factor Xa (FXa). Here we show that two short evolutionary conserved segments of the B-domain, together termed the procofactor regulatory region, serve an essential autoinhibitory function. This newly identified motif consists of a basic (963-1008) and an acidic (1493-1537) region and defines the minimal sequence requirements to maintain FV as a procofactor. Our data suggest that dismantling this autoinhibitory region via deletion or proteolysis is the driving force to unveil a high affinity binding site(s) for FXa. These findings document an unexpected sequence-specific role for the B-domain by negatively regulating FV function and preventing activity of the procofactor. These new mechanistic insights point to new ways in which the FV procofactor to cofactor transition could be modulated to alter hemostasis.     

Minimization of intermediate concentrations as a suggested optimality principle for biochemical networks

The criterion of minimum intermediate concentrations in steady states is suggested to be of essential relevance in the evolution of biochemical reaction networks. This extremum principle is phrased in two different ways, firstly in terms of total osmolarity of intermediates and, secondly, as a multiple criterion problem. The relationships between the two problems are elucidated and a solving method for the latter is then given. It turns out that in each optimal state, the network can be subdivided into a slow and a fast subsystem. The notion of convex conservation relations is introduced and the implications of such relations for the optimization problem are investigated.     

Crystal structure of enterotoxigenic Escherichia coli colonization factor CS6 reveals a novel type of functional assembly

Coli surface antigen 6 (CS6) is a widely expressed enterotoxigenic Escherichia coli (ETEC) colonization factor that mediates bacterial attachment to the small intestinal epithelium. CS6 is a polymer of two protein subunits CssA and CssB, which are secreted and assembled on the cell surface via the CssC/CssD chaperone usher (CU) pathway. Here, we present an atomic resolution model for the structure of CS6 based on the results of X‐ray crystallographic, spectroscopic and biochemical studies, and suggest a mechanism for CS6‐mediated adhesion. We show that the CssA and CssB subunits are assembled alternately in linear fibres by the principle of donor strand complementation. This type of fibre assembly is novel for CU assembled adhesins. We also show that both subunits in the fibre bind to receptors on epithelial cells, and that CssB, but not CssA, specifically recognizes the extracellular matrix protein fibronectin. Taken together, structural and functional results suggest that CS6 is an adhesive organelle of a novel type, a hetero‐polyadhesin that is capable of polyvalent attachment to different receptors.     

Restoring the Procofactor State of Factor Va-like Variants by Complementation with B-domain Peptides

Coagulation factor V (FV) circulates as an inactive procofactor and is activated to FVa by proteolytic removal of a large inhibitory B-domain. Conserved basic and acidic sequences within the B-domain appear to play an important role in keeping FV as an inactive procofactor. Here, we utilized recombinant B-domain fragments to elucidate the mechanism of this FV autoinhibition. We show that a fragment encoding the basic region (BR) of the B-domain binds with high affinity to cofactor-like FV(a) variants that harbor an intact acidic region. Furthermore, the BR inhibits procoagulant function of the variants, thereby restoring the procofactor state. The BR competes with FXa for binding to FV(a), and limited proteolysis of the B-domain, specifically at Arg¹⁵⁴⁵, ablates BR binding to promote high affinity association between FVa and FXa. These results provide new insight into the mechanism by which the B-domain stabilizes FV as an inactive procofactor and reveal how limited proteolysis of FV progressively destabilizes key regulatory regions of the B-domain to produce an active form of the molecule.     

旱地小麦根际细菌中产生1-氨基环丙烷-1-羧酸(ACC)脱氨酶菌株的分离和鉴定

采用富集定向筛选法,从旱地小麦的根际土壤中分离到2株产生1-氨基环丙烷-1-羧酸(ACC)脱氨酶的菌株AS和CS。经测定菌株AS和CS的ACC脱氨酶的比活力分别为0.018 6 U/mg和0.016 7 U/mg蛋白。根据培养特征观察和生理生化指标测定结果,结合16S rDNA碱基序列测定和系统发育同源性分析,确定菌株AS和菌株CS分别属于霍氏肠杆菌(Enterobacter hormaechei)和变形斑沙雷氏菌(Serratia proteamaculans)。     

Petri net modeling of high-order genetic systems using grammatical evolution

Understanding how DNA sequence variations impact human health through a hierarchy of biochemical and physiological systems is expected to improve the diagnosis, prevention, and treatment of common, complex human diseases. We have previously developed a hierarchical dynamic systems approach based on Petri nets for generating biochemical network models that are consistent with genetic models of disease susceptibility. This modeling approach uses an evolutionary computation approach called grammatical evolution as a search strategy for optimal Petri net models. We have previously demonstrated that this approach routinely identifies biochemical network models that are consistent with a variety of genetic models in which disease susceptibility is determined by nonlinear interactions between two DNA sequence variations. In the present study, we evaluate whether the Petri net approach is capable of identifying biochemical networks that are consistent with disease susceptibility due to higher order nonlinear interactions between three DNA sequence variations. The results indicate that our model-building approach is capable of routinely identifying good, but not perfect, Petri net models. Ideas for improving the algorithm for this high-dimensional problem are presented.     

Effect of cigarette smoking on nitric oxide, structural, and mechanical properties of mouse arteries

Cigarette smoking (CS) is a major risk factor for vascular disease. The aim of this study was to quantitatively assess the influence of CS on mouse arteries. We studied the effect of short-term (6 wk) and long-term (16 wk) CS exposure on structural and mechanical properties of coronary arteries compared with that of control mice. We also examined the reversibility of the deleterious effects of CS on structural [e.g., wall thickness (WT)], mechanical (e.g., stiffness), and biochemical [e.g., nitric oxide (NO) by-products] properties with the cessation of CS. The left and right coronary arteries were cannulated in situ and mechanically distended. The stress, strain, elastic modulus, and WT of coronary arteries were determined. Western blot analysis was used to analyze endothelial NO synthase (eNOS) in the femoral and carotid arteries of the same mice, and NO by-products were determined by measuring the levels of nitrite. Our results show that the mean arterial pressure was increased by CS. Furthermore, CS significantly increased the elastic modulus, decreased stress and strain, and increased the WT and WT-to-radius ratio compared with those of control mice. The reduction of eNOS protein expression was found only after long-term CS exposure. Moreover, the NO metabolite was markedly decreased in CS mice after short- and long-term exposure of CS. These findings suggest that 16 wk of CS exposure can cause an irreversible deterioration of structural and elastic properties of mouse coronary arteries. The decrease in endothelium-derived NO in CS mice was seen to significantly correlate with the remodeling of arterial wall.     

Effects of calcium and phosphate on the corpuscles of Stannius of the teleost fish,Oreochromis mossambicus

Summary Removal of the corpuscles of Stannius (CS) in Oreochromis mossambicus leads to hypercalcemia and hypophosphatemia. The effects on CS size and ultrastructure of different calcium and phosphate concentrations of the ambient water and of the food were investigated. A six-fold increase of the calcium concentration of the water leads to a four-fold increase in CS volume; this is mainly caused by an increase in the size and number of the type-1 cells. The effect of external calcium is most probably mediated by the calcium concentration of the blood plasma. Plasma ionic calcium may be the relevant factor. Changes in the calcium concentration of the food had no effect on the CS. Similarly, hyperphosphatemia or hypophosphatemia induced by high phosphate concentrations of the water or the food, or by a phosphate-deficient diet, had no noticeable effect on the CS. The results support the hypothesis that the type-1 cells produce the hypocalcemic factor of the CS. There is no evidence for the production by the CS of an endocrine factor involved in the control of phosphate metabolism.     

Impact of limited solvent capacity on metabolic rate, enzyme activities, and metabolite concentrations of S. cerevisiae glycolysis

The cell's cytoplasm is crowded by its various molecular components, resulting in a limited solvent capacity for the allocation of new proteins, thus constraining various cellular processes such as metabolism. Here we study the impact of the limited solvent capacity constraint on the metabolic rate, enzyme activities, and metabolite concentrations using a computational model of Saccharomyces cerevisiae glycolysis as a case study. We show that given the limited solvent capacity constraint, the optimal enzyme activities and the metabolite concentrations necessary to achieve a maximum rate of glycolysis are in agreement with their experimentally measured values. Furthermore, the predicted maximum glycolytic rate determined by the solvent capacity constraint is close to that measured in vivo. These results indicate that the limited solvent capacity is a relevant constraint acting on S. cerevisiae at physiological growth conditions, and that a full kinetic model together with the limited solvent capacity constraint can be used to predict both metabolite concentrations and enzyme activities in vivo.     

Chick Embryo Extract, Muscle Trophic Factor and Chick and Horse Sera as Environments for Chick Myogenic Cell Growth

Chick myogenic cells grew in a medium composed of Eagle's minimum essential medium (MEM), horse serum (HS), and one of the essential factors needed for myogenic cell growth (EFMG), that is, chick embryo extract (EE), chick serum (CS), or the muscle trophic factor (MTF). But they did not grow in the absence of the EFMG. In the absence of HS, they scarcely grew in a medium composed of MEM, and EE or MTF. They grew in a medium composed of MEM and CS; they grew much better in a medium composed of MEM, CS, and HS.
In the presence of one of the EFMG, the optimal HS concentration for growth varied depending on its lot. At higher HS concentrations, growth was suppressed. Further, it was suggested that an inhibitory substance(s) for myogenic cell growth was present in HS. The inhibitory effects can usually be minimized by diluting the serum with an artificial medium.     

Control analysis of metabolic networks. 2. Total differentials and general formulation of the connectivity relations

The mathematical background of the connectivity relations of metabolic control theory is analysed. The connectivity relations are shown to reflect general properties of total differentials of reaction rate vi, flux J, and metabolite concentration Xj. Connectivity relations hold for any metabolic network in which all vi are homogeneous functions of enzyme concentration Ei. This notion allows established algebraic methods to be used for the formulation of connectivity relations for metabolic systems in which numerous constraints are imposed on metabolite concentrations. A general procedure to derive connectivity relations for such metabolic systems is given. To encourage a broader audience to apply control theory to physiological systems, an easy-to-use graphical procedure is derived for formulating connectivity relations for biochemical systems in which no metabolite is involved in more than one constraint.     

Cold shock and cold acclimation proteins in the psychrotrophic bacterium Arthrobacter globiformis SI55.

The psychrotrophic bacterium Arthrobacter globiformis SI55 was grown at 4 and 25 degrees C, and the cell protein contents were analyzed by two-dimensional electrophoresis. Cells subjected to cold shocks of increasing magnitude were also analyzed. Correspondence analysis of protein appearance distinguished four groups of physiological significance. Group I contained cold shock proteins (Csps) overexpressed only after a large temperature downshift. Group II contained Csps with optimal expression after mild shocks. Group III contained proteins overexpressed after all cold shocks. These last proteins were also overexpressed in cells growing at 4 degrees C and were considered to be early cold acclimation proteins (Caps). Group IV contained proteins which were present at high concentrations only in 4 degrees C steady-state cells and appeared to be late Caps. A portion of a gene very similar to the Escherichia coli cspA gene (encoding protein CS7.4) was identified. A synthetic peptide was used to produce an antibody which detected a CS7.4-like protein (A9) by immunoblotting two-dimensional electrophoresis gels of A. globiformis SI55 total proteins. Unlike mesophilic microorganisms, this CS7.4-like protein was still produced during prolonged growth at low temperature, and it might have a particular adaptive function needed for balanced growth under harsh conditions. However, A9 was induced at high temperature by chloramphenicol, suggesting that CS7.4-like proteins have a more general role than their sole implication in cold acclimation processes.     

Direct measurement of the interactions of glycosaminoglycans and a heparin decasaccharide with the malaria circumsporozoite protein

Circumsporozoite (CS) protein is a predominant surface antigen of malaria sporozoites, the infective form of the parasite, and has been used for making anti-malaria vaccines. For the first time we have examined the interaction of CS protein with various glycosaminoglycans in real time using surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC). Heparin was the best binder among the glycosaminoglycans tested and bound to CS protein with nanomolar affinity. Using purified and structurally defined small heparin oligosaccharides, we identified a decasaccharide to be the minimum sized CS protein-binding sequence. In an indirect competition assay, this decasaccharide blocked the CS protein interaction with HepG2 cells with an ID(50) of less than 60 nM. The decasaccharide has a structure commonly found in hepatic heparan sulfate, and the same sequence has recently been shown to bind specifically to apolipoprotein E. Examination of porcine liver heparan sulfate in this indirect competition assay showed that it and heparin were the only glycosaminoglycans that could effectively block CS protein interaction with HepG2 cells in culture. These data support the hypothesis that the invasion of liver cells by the parasite shares a common mechanism with the hepatic uptake of lipoprotein remnants from the blood.     

Genomic Evidence for a Simpler Clotting Scheme in Jawless Vertebrates

Mammalian blood clotting involves numerous components, most of which are the result of gene duplications that occurred early in vertebrate evolution and after the divergence of protochordates. As such, the genomes of the jawless fish (hagfish and lamprey) offer the best possibility for finding systems that might have a reduced set of the many clotting factors observed in higher vertebrates. The most straightforward way of inventorying these factors may be through whole genome sequencing. In this regard, the NCBI Trace database ( http://www.ncbi.nlm.nih.gov/Traces/trace.cgi ) for the lamprey (Petromyzon marinus) contains more than 18 million raw DNA sequences determined by whole-genome shotgun methodology. The data are estimated to be about sixfold redundant, indicating that coverage is sufficiently complete to permit judgments about the presence or absence of particular genes. A search for 20 proteins whose sequences were determined prior to the trace database study found all 20. A subsequent search for specified coagulation factors revealed a lamprey system with a smaller number of components than is found in other vertebrates in that factors V and VIII seem to be represented by a single gene, and factor IX, which is ordinarily a cofactor of factor VIII, is not present. Fortuitously, after the completion of the survey of the Trace database, a draft assembly based on the same database was posted. The draft assembly allowed many of the identified Trace fragments to be linked into longer sequences that fully support the conclusion that lampreys have a simpler clotting scheme compared with other vertebrates. The data are also consistent with the hypothesis that a whole-genome duplication or other large scale block duplication occurred after the divergence of jawless fish from other vertebrates and allowed the simultaneous appearance of a second set of two functionally paired proteins in the vertebrate clotting scheme.     

Inference of biochemical network models in S-system using multiobjective optimization approach

MOTIVATION: The inference of biochemical networks, such as gene regulatory networks, protein-protein interaction networks, and metabolic pathway networks, from time-course data is one of the main challenges in systems biology. The ultimate goal of inferred modeling is to obtain expressions that quantitatively understand every detail and principle of biological systems. To infer a realizable S-system structure, most articles have applied sums of magnitude of kinetic orders as a penalty term in the fitness evaluation. How to tune a penalty weight to yield a realizable model structure is the main issue for the inverse problem. No guideline has been published for tuning a suitable penalty weight to infer a suitable model structure of biochemical networks. RESULTS: We introduce an interactive inference algorithm to infer a realizable S-system structure for biochemical networks. The inference problem is formulated as a multiobjective optimization problem to minimize simultaneously the concentration error, slope error and interaction measure in order to find a suitable S-system model structure and its corresponding model parameters. The multiobjective optimization problem is solved by the epsilon-constraint method to minimize the interaction measure subject to the expectation constraints for the concentration and slope error criteria. The theorems serve to guarantee the minimum solution for the epsilon-constrained problem to achieve the minimum interaction network for the inference problem. The approach could avoid assigning a penalty weight for sums of magnitude of kinetic orders.