Minimal cut sets in biochemical reaction networks

MOTIVATION: Structural studies of metabolic networks yield deeper insight into topology, functionality and capabilities of the metabolisms of different organisms. Here, we address the analysis of potential failure modes in metabolic networks whose occurrence will render the network structurally incapable of performing certain functions. Such studies will help to identify crucial parts in the network structure and to find suitable targets for repressing undesired metabolic functions. RESULTS: We introduce the concept of minimal cut sets for biochemical networks. A minimal cut set (MCS) is a minimal (irreducible) set of reactions in the network whose inactivation will definitely lead to a failure in certain network functions. We present an algorithm which enables the computation of the MCSs in a given network related to user-defined objective reactions. This algorithm operates on elementary modes. A number of potential applications are outlined, including network verifications, phenotype predictions, assessing structural robustness and fragility, metabolic flux analysis and target identification in drug discovery. Applications are illustrated by the MCSs in the central metabolism of Escherichia coli for growth on different substrates. AVAILABILITY: Computation and analysis of MCSs is an additional feature of the FluxAnalyzer (freely available for academic users upon request, special contracts for industrial companies; see web page below). Supplementary information: http://www.mpi-magdeburg.mpg.de/projects/fluxanalyzer     

Selecting SNPs for association studies based on population frequencies: a novel interactive tool and its application to polygenic diseases

Common complex polygenic diseases as autoimmune diseases have not been completely understood on a molecular level. While many genes are known to be involved in the pathways responsible for the phenotype, explicit causes for the susceptibility of the disease remain to be elucidated. The susceptibility to disease is thought to be the result of genetic epistatic interactions between common polymorphic genes. This polymorphism is mostly caused by single nucleotide polymorphisms (SNPs). Human subpopulations are known to differ in the susceptibility to the diseases and generally in the distribution of single nucleotide polymorphisms. The here presented approach retrieves SNPs with the most divergent frequencies for selected human subpopulations to help defining properties for the experimental verification of SNPs within defined regions. A web-accessible program implementing this approach was evaluated for multiple sclerosis (MS), a common human polygenic disease. A link to a summary of data from "The SNP Consortium" (TSC) with sex-dependencies of SNPs is available. Associations of SNPs to genes, genetic markers and chromosomal loci are retrieved from the Ensembl project. This tool is recommended to be used in conjunction with microarray analyses or marker association studies that link genes or chromosomal loci to particular diseases.     

Sera of lung cancer patients affect the release of Th1, Th2 and monocyte-derived cytokines, and the expression of IL-2Ralpha by normal, stimulated mononuclear cells

We have shown that the sera of lung cancer patients affect the response of ConA-stimulated normal peripheral blood mononuclear cells by decreasing the expression of IL-2Ralpha and inhibiting the release of IL-1beta and IL-2. A tendency to enhance the release of IL-6 was also observed. We conclude that an imbalance in the Th1/Th2 cytokine response, typical for cancer patients, may at least partly be related to soluble factors circulating in the patients' blood. We discuss a putative role of serum IL-10, IL-1ra, and soluble IL-2Ralpha in the effects observed.     

Abi1 is essential for the formation and activation of a WAVE2 signalling complex

WAVE2 belongs to a family of proteins that mediates actin reorganization by relaying signals from Rac to the Arp2/3 complex, resulting in lamellipodia protrusion. WAVE2 displays Arp2/3-dependent actin nucleation activity in vitro, and does not bind directly to Rac. Instead, it forms macromolecular complexes that have been reported to exert both positive and negative modes of regulation. How these complexes are assembled, localized and activated in vivo remains to be established. Here we use tandem mass spectrometry to identify an Abi1-based complex containing WAVE2, Nap1 (Nck-associated protein) and PIR121. Abi1 interacts directly with the WHD domain of WAVE2, increases WAVE2 actin polymerization activity and mediates the assembly of a WAVE2-Abi1-Nap1-PIR121 complex. The WAVE2-Abi1-Nap1-PIR121 complex is as active as the WAVE2-Abi1 sub-complex in stimulating Arp2/3, and after Rac activation it is re-localized to the leading edge of ruffles in vivo. Consistently, inhibition of Abi1 by RNA interference (RNAi) abrogates Rac-dependent lamellipodia protrusion. Thus, Abi1 orchestrates the proper assembly of the WAVE2 complex and mediates its activation at the leading edge in vivo.     

Material changes in osteoporotic human cancellous bone following infiltration with acrylic bone cement for a vertebral cement augmentation

Bone cement infiltration can be effective at mechanically augmenting osteoporotic vertebrae. While most published literature describes the gain in mechanical strength of augmented vertebrae, we report the first measurements of viscoelastic material changes of cancellous bone due to cement infiltration. We infiltrated cancellous core specimen harvested from osteoporotic cadaveric spines with acrylic bone cement. Bone specimen before and after cement infiltration were subjected to identical quasi-static and relaxation loading in confined and free compression. Testing data were fitted to a linear viscoelastic model of compressible material and the model parameters for cement, native cancellous bone, and cancellous bone infiltrated (composite) with cement were identified. The fitting demonstrated that the linear viscoelastic model presented in this paper accurately describes the mechanical behaviour of cement and bone, before and after infiltration. Although the composite specimen did not completely adopt the properties of bulk bone cement, the stiffening of cancellous bone due to cement infiltration is considerable. The composite was, for example, 8.5 times stiffer than native bone. The local stiffening of cancellous bone in patients may alter the load transfer of the augmented motion segment and may be the cause of subsequent fractures in the vertebrae adjacent to the ones infiltrated with cement. The material model and parameters in this paper, together with an adequate finite-element model, can be helpful to investigate the load shift, the mechanism for subsequent fractures, and filling patterns for ideal cement infiltration.     

Structural mechanism of specific ligand recognition by a lipocalin tailored for the complexation of digoxigenin

DigA16 is an artificial digoxigenin-binding protein, which was derived from the bilin-binding protein, a lipocalin of Pieris brassicae, via reshaping of its natural ligand pocket. Here we report the crystal structures of DigA16 in the presence of either digoxigenin or digitoxigenin and for the apo-protein at resolutions below 1.9A. As a consequence of the altogether 17 amino acid substitutions within the binding site significant structural changes have occurred in the four loops that form the entrance to the ligand pocket on top of the structurally conserved beta-barrel framework. For example, one loop adopts a new alpha-helical backbone structure, which seems to be induced by few critical side-chain contacts. Digoxigenin becomes almost fully buried (by 95%) upon complexation, whereby specificity for the hydrophilic steroid is maintained through hydrogen-bonding networks and shape complementarity. The differential binding of the related steroid digitoxigenin is mainly governed by an internal histidine residue, whose side-chain undergoes significant induced fit. Among those amino acids that line the ligand pocket two tyrosine and one tryptophan residue provide the largest contacts. Interestingly, corresponding three side-chains are found with the same mutual orientation in the anti-digoxigenin antibody 26-10, even though the hapten orientation is quite different there and only 66% of the steroid surface is buried in the combining site. Hence, in the case of the engineered lipocalin DigA16 an example of convergent in vitro evolution is observed. Generally, the remarkable structural plasticity of the loop region and the role of polar residues in the binding site illustrate the potential of the lipocalin scaffold for the generation of specific receptor proteins towards a variety of ligands.     

Oligodendrocyte-specific expression of human immunodeficiency virus type 1 Nef in transgenic mice leads to vacuolar myelopathy and alters oligodendrocyte phenotype in vitro

Vacuolar myelopathy (VM) is a frequent central nervous system complication of human immunodeficiency virus type 1 (HIV-1) infection. We report here that transgenic (Tg) mice expressing even low levels of Nef in oligodendrocytes under the regulation of the myelin basic protein (MBP) promoter (MBP/HIV(Nef)) developed VM similar to the human disease in its appearance and topography. The spinal cords of these Tg mice showed lower levels of the myelin proteins MAG and CNPase and of the 21-kDa isoform of MBP prior to the development of vacuoles. In addition, Tg oligodendrocytes in primary in vitro cultures appeared morphologically more mature but, paradoxically, exhibited a less mature phenotype based on O4, O1, CNPase, and MBP staining. In particular, mature CNPase(+) MBP(+) Tg oligodendrocytes were less numerous than non-Tg oligodendrocytes. Therefore, Nef appears to affect the proper differentiation of oligodendrocytes. These data suggest that even low levels of Nef expression in human oligodendrocytes may be responsible for the development of VM in HIV-1-infected individuals.     

Three distinct epitopes on the extracellular face of the glucagon receptor determine specificity for the glucagon amino terminus

The glucagon and glucagon-like peptide-1 (GLP-1) receptors are homologous family B seven-transmembrane (7TM) G protein-coupled receptors, and they selectively recognize the homologous peptide hormones glucagon (29 amino acids) and GLP-1 (30-31 amino acids), respectively. The amino-terminal extracellular domain of the glucagon and GLP-1 receptors (140-150 amino acids) determines specificity for the carboxyl terminus of glucagon and GLP-1, respectively. In addition, the glucagon receptor core domain (7TM helices and connecting loops) strongly determines specificity for the glucagon amino terminus. Only 4 of 15 residues are divergent in the glucagon and GLP-1 amino termini; Ser2, Gln3, Tyr10, and Lys12 in glucagon and the corresponding Ala8, Glu9, Val16, and Ser18 in GLP-1. In this study, individual substitution of these four residues of glucagon with the corresponding residues of GLP-1 decreased the affinity and potency at the glucagon receptor relative to glucagon. Substitution of distinct segments of the glucagon receptor core domain with the corresponding segments of the GLP-1 receptor rescued the affinity and potency of specific glucagon analogs. Site-directed mutagenesis identified the Asp385 --> Glu glucagon receptor mutant that specifically rescued Ala2-glucagon. The results show that three distinct epitopes of the glucagon receptor core domain determine specificity for the N terminus of glucagon. We suggest a glucagon receptor binding model in which the extracellular ends of TM2 and TM7 are close to and determine specificity for Gln3 and Ser2 of glucagon, respectively. Furthermore, the second extracellular loop and/or proximal segments of TM4 and/or TM5 are close to and determine specificity for Lys12 of glucagon.