The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

Peptide design by artificial neural networks and computer-based evolutionary search

A technique for systematic peptide variation by a combination of rational and evolutionary approaches is presented. The design scheme consists of five consecutive steps: (i) identification of a “seed peptide” with a desired activity, (ii) generation of variants selected from a physicochemical space around the seed peptide, (iii) synthesis and testing of this biased library, (iv) modeling of a quantitative sequence-activity relationship by an artificial neural network, and (v) de novo design by a computer-based evolutionary search in sequence space using the trained neural network as the fitness function. This strategy was successfully applied to the identification of novel peptides that fully prevent the positive chronotropic effect of anti-β1-adrenoreceptor autoantibodies from the serum of patients with dilated cardiomyopathy. The seed peptide, comprising 10 residues, was derived by epitope mapping from an extracellular loop of human β1-adrenoreceptor. A set of 90 peptides was synthesized and tested to provide training data for neural network development. De novo design revealed peptides with desired activities that do not match the seed peptide sequence. These results demonstrate that computer-based evolutionary searches can generate novel peptides with substantial biological activity.

rRNA complementarity within mRNAs: A possible basis for mRNA-ribosome interactions and translational control

Our recent demonstration that many eukaryotic mRNAs contain sequences complementary to rRNA led to the hypothesis that these sequences might mediate specific interactions between mRNAs and ribosomes and thereby affect translation. In the present experiments, the ability of complementary sequences to bind to rRNA was investigated by using photochemical cross-linking. RNA probes with perfect complementarity to 18S or 28S rRNA were shown to cross-link specifically to the corresponding rRNA within intact ribosomal subunits. Similar results were obtained by using probes based on natural mRNA sequences with varying degrees of complementarity to the 18S rRNA. RNase H cleavage localized four such probes to complementary regions of the 18S rRNA. The effects of complementarity on translation were assessed by using the mRNA encoding ribosomal protein S15. This mRNA contains a sequence within its coding region that is complementary to the 18S rRNA at 20 of 22 nucleotides. RNA from an S15-luciferase fusion construct was translated in a cell-free lysate and compared with the translation of four related constructs that were mutated to decrease complementarity to the 18S rRNA. These mutations did not alter the amino acid sequence or the codon bias. A correlation between complementarity and translation was observed; constructs with less complementarity increased the amount of translation up to 54%. These findings raised the possibility that direct base-pairing of particular mRNAs to rRNAs within ribosomes may function as a mechanism of translational control.

Structure and function in rhodopsin: Packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain are coupled*

A previous study of the retinitis pigmentosa mutation L125R and two designed mutations at this site, L125A and L125F, showed that these mutations cause partial or total misfolding of the opsins expressed in COS cells from the corresponding mutant opsin genes. We now report on expression and characterization of the opsins from the following retinitis pigmentosa mutants in the transmembrane domain of rhodopsin that correspond to six of the seven helices: G51A and G51V (helix A), G89D (helix B), A164V (helix D), H211P (helix E), P267L and P267R (helix F), and T297R (helix G). All the mutations caused partial misfolding of the opsins as observed by the UV/visible absorption characteristics and by separation of the expressed opsins into fractions that bound 11-cis-retinal to form the corresponding mutant rhodopsins and those that did not bind 11-cis-retinal. Further, all the mutant rhodopsins prepared from the above mutants, except for G51A, showed strikingly abnormal bleaching behavior with abnormal metarhodopsin II photointermediates. The results show that retinitis pigmentosa mutations in every one of the transmembrane helices can cause misfolding of the opsin. Therefore, on the basis of these and previous results, we conclude that defects in the packing of the transmembrane helices resulting from these mutations are relayed to the intradiscal domain, where they cause misfolding of the opsin by inducing the formation of a disulfide bond other than the native Cys-110—Cys-187 disulfide bond. Thus, there is coupling between packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain.

Incomplete penetrance of familial retinoblastoma linked to germ-line mutations that result in partial loss of RB function

To study the molecular basis for the clinical phenotype of incomplete penetrance of familial retinoblastoma, we have examined the functional properties of three RB mutations identified in the germ line of five different families with low penetrance. RB mutants isolated from common adult cancers and from classic familial retinoblastoma (designated as classic RB mutations) are unstable and generally do not localize to the nucleus, do not undergo cyclin-dependent kinase (cdk)-mediated hyperphosphorylation, show absent protein “pocket” binding activity, and do not suppress colony growth of RB(−) cells. In contrast, two low-penetrant alleles (661W and “deletion of codon 480”) retained the ability to localize to the nucleus, showed normal cdk-mediated hyperphosphorylation in vivo, exhibited a binding pattern to simian virus 40 large T antigen using a quantitative yeast two-hybrid assay that was intermediate between classic mutants (null) and wild-type RB, and had absent E2F1 binding in vitro. A third, low-penetrant allele, “deletion of RB exon 4,” showed minimal hyperphosphorylation in vivo but demonstrated detectable E2F1 binding in vitro. In addition, each low-penetrant RB mutant retained the ability to suppress colony growth of RB(−) tumor cells. These findings suggest two categories of mutant, low-penetrant RB alleles. Class 1 alleles correspond to promoter mutations, which are believed to result in reduced or deregulated levels of wild-type RB protein, whereas class 2 alleles result in mutant proteins that retain partial activity. Characterization of the different subtypes of class 2 low-penetrant genes may help to define more precisely functional domains within the RB product required for tumor suppression.

Suppression of an absolute defect in Type IV pilus biogenesis by loss-of-function mutations in pilT, a twitching motility gene in Neisseria gonorrhoeae

Type IV pili of Neisseria gonorrhoeae, the Gram-negative etiologic agent of gonorrhea, facilitate colonization of the human host. Gonococcal PilT, a protein belonging to a large family of molecules sharing a highly conserved nucleotide binding domain motif, has been shown to be dispensable for organelle biogenesis but essential for twitching motility and competence for genetic transformation. Here, we show that the defect in pilus biogenesis resulting from mutations in the pilC gene, encoding a putative pilus-associated adhesin for human tissue, can be suppressed by the absence of functional PilT. These data conclusively demonstrate that PilT influences the Type IV pilus biogenesis pathway and strongly suggest that organelle expression is a dynamic process. In addition, these findings imply that PilT antagonizes the process of organelle biogenesis and provide the basis for a model for how the counteractive roles of PilT and PilC might relate mechanistically to the phenomenon of twitching motility.

Molecular Cloning and Characterization of a Radial Spoke Head Protein of Sea Urchin Sperm Axonemes: Involvement of the Protein in the Regulation of Sperm Motility

Monoclonal antibodies raised against axonemal proteins of sea urchin spermatozoa have been used to study regulatory mechanisms involved in flagellar motility. Here, we report that one of these antibodies, monoclonal antibody D-316, has an unusual perturbating effect on the motility of sea urchin sperm models; it does not affect the beat frequency, the amplitude of beating or the percentage of motile sperm models, but instead promotes a marked transformation of the flagellar beating pattern which changes from a two-dimensional to a three-dimensional type of movement. On immunoblots of axonemal proteins separated by SDS-PAGE, D-316 recognized a single polypeptide of 90 kDa. This protein was purified following its extraction by exposure of axonemes to a brief heat treatment at 40°C. The protein copurified and coimmunoprecipitated with proteins of 43 and 34 kDa, suggesting that it exists as a complex in its native form. Using D-316 as a probe, a full-length cDNA clone encoding the 90-kDa protein was obtained from a sea urchin cDNA library. The sequence predicts a highly acidic (pI = 4.0) protein of 552 amino acids with a mass of 62,720 Da (p63). Comparison with protein sequences in databases indicated that the protein is related to radial spoke proteins 4 and 6 (RSP4 and RSP6) of Chlamydomonas reinhardtii, which share 37% and 25% similarity, respectively, with p63. However, the sea urchin protein possesses structural features distinct from RSP4 and RSP6, such as the presence of three major acidic stretches which contains 25, 17, and 12 aspartate and glutamate residues of 34-, 22-, and 14-amino acid long stretches, respectively, that are predicted to form α-helical coiled-coil secondary structures. These results suggest a major role for p63 in the maintenance of a planar form of sperm flagellar beating and provide new tools to study the function of radial spoke heads in more evolved species.

Assembly of a Novel Cartilage Matrix Protein Filamentous Network: Molecular Basis of Differential Requirement of von Willebrand Factor A Domains

Cartilage matrix protein (CMP) is the prototype of the newly discovered matrilin family, all of which contain von Willebrand factor A domains. Although the function of matrilins remain unclear, we have shown that, in primary chondrocyte cultures, CMP (matrilin-1) forms a filamentous network, which is made up of two types of filaments, a collagen-dependent one and a collagen-independent one. In this study, we demonstrate that the collagen-independent CMP filaments are enriched in pericellular compartments, extending directly from chondrocyte membranes. Their morphology can be distinguished from that of collagen filaments by immunogold electron microscopy, and mimicked by that of self-assembled purified CMP. The assembly of CMP filaments can occur from transfection of a wild-type CMP transgene alone in skin fibroblasts, which do not produce endogenous CMP. Conversely, assembly of endogenous CMP filaments by chondrocytes can be inhibited specifically by dominant negative CMP transgenes. The two A domains within CMP serve essential but different functions during network formation. Deletion of the A2 domain converts the trimeric CMP into a mixture of monomers, dimers, and trimers, whereas deletion of the A1 domain does not affect the trimeric configuration. This suggests that the A2 domain modulates multimerization of CMP. Absence of either A domain from CMP abolishes its ability to form collagen-independent filaments. In particular, Asp22 in A1 and Asp255 in A2 are essential; double point mutation of these residues disrupts CMP network formation. These residues are part of the metal ion–dependent adhesion sites, thus a metal ion–dependent adhesion site–mediated adhesion mechanism may be applicable to matrilin assembly. Taken together, our data suggest that CMP is a bridging molecule that connects matrix components in cartilage to form an integrated matrix network.

Molecular basis for a link between complement and the vascular complications of diabetes

Activated terminal complement proteins C5b to C9 form the membrane attack complex (MAC) pore. Insertion of the MAC into endothelial cell membranes causes the release of growth factors that stimulate tissue growth and proliferation. The complement regulatory membrane protein CD59 restricts MAC formation. Because increased cell proliferation characterizes the major chronic vascular complications of human diabetes and because increased glucose levels in diabetes cause protein glycation and impairment of protein function, we investigated whether glycation could inhibit CD59. Glycation-inactivation of CD59 would cause increased MAC deposition and MAC-stimulated cell proliferation. Here, we report that (i) human CD59 is glycated in vivo, (ii) glycated human CD59 loses its MAC-inhibitory function, and (iii) inactivation of CD59 increases MAC-induced growth factor release from endothelial cells. We demonstrate by site-directed mutagenesis that residues K41 and H44 form a preferential glycation motif in human CD59. The presence of this glycation motif in human CD59, but not in CD59 of other species, may help explain the distinct propensity of humans to develop vascular proliferative complications of diabetes.

A structural basis for integrin activation by the cytoplasmic tail of the αIIb-subunit

A key step in the activation of heterodimeric integrin adhesion receptors is the transmission of an agonist-induced cellular signal from the short α- and/or β-cytoplasmic tails to the extracellular domains of the receptor. The structural details of how the cytoplasmic tails mediate such an inside-out signaling process remain unclear. We report herein the NMR structures of a membrane-anchored cytoplasmic tail of the αIIb-subunit and of a mutant αIIb-cytoplasmic tail that renders platelet integrin αIIbβ3 constitutively active. The structure of the wild-type αIIb-cytoplasmic tail reveals a “closed” conformation where the highly conserved N-terminal membrane-proximal region forms an α-helix followed by a turn, and the acidic C-terminal loop interacts with the N-terminal helix. The structure of the active mutant is significantly different, having an “open” conformation where the interactions between the N-terminal helix and C-terminal region are abolished. Consistent with these structural differences, the two peptides differ in function: the wild-type peptide suppressed αIIbβ3 activation, whereas the mutant peptide did not. These results provide an atomic explanation for extensive biochemical/mutational data and support a conformation-based “on/off switch” model for integrin activation.

Requirement of phosphatidylglycerol for photosynthetic function in thylakoid membranes

To investigate the role of phosphatidylglycerol (PG) in photosynthesis, we constructed a mutant defective in the CDP-diacylglycerol synthase gene from a cyanobacterium, Synechocystis sp. PCC6803. The mutant, designated as SNC1, required PG supplementation for growth. Growth was repressed in PG-free medium concomitantly with the decrease in cellular content of PG. These results indicate that PG is essential, and that SNC1 is defective in PG synthesis. Decrease in PG content was accompanied by a reduction in the cellular content of chlorophyll, but with little effect on the contents of phycobilisome pigments, which showed that levels of chlorophyll–protein complexes decreased without alteration of those of phycobilisomes. Regardless of the decrease in the PG content, CO2-dependent photosynthesis by SNC1 was similar to that by the wild type on a chlorophyll basis, but consequently became lower on a cell basis. Simultaneously, the ratio of oxygen evolution of photosystem II (PSII) measured with p-benzoquinone to that of CO2-dependent photosynthesis, which ranged between 1.3 and 1.7 in the wild type. However, it was decreased in SNC1 from 1.3 to 0.4 during the early growth phase where chlorophyll content and CO2-dependent photosynthesis were little affected, and then finally to 0.1, suggesting that PSII first lost its ability to reduce p-benzoquinone and then decreased in its level and actual activity. These results indicate that PG contributes to the accumulation of chlorophyll–protein complexes in thylakoid membranes, and also to normal functioning of PSII.

Molecular basis of lutropin recognition by the mannose/GalNAc-4-SO4 receptor

The circulatory half-life of the glycoprotein hormone lutropin (LH) is precisely regulated by the mannose (Man)/GalNAc-4-SO4 receptor expressed in hepatic endothelial cells. Rapid clearance from the circulation contributes to the episodic rise and fall of LH levels that is essential for maximal stimulation of the G protein-coupled LH receptor. We have defined two molecular forms of the Man/GalNAc-4-SO4 receptor that differ in ligand specificity, cell and tissue expression, and function. The form expressed by hepatic endothelial cells binds GalNAc-4-SO4-bearing ligands and regulates hormone circulatory half-life, whereas the form expressed by macrophages binds Man-bearing ligands and may play a role in innate immunity. We demonstrate that the GalNAc-4-SO4-specific form in hepatic endothelial cells is dimeric whereas the Man-specific form in lung macrophages is monomeric, accounting for the different ligand specificities of the receptor expressed in these tissues. Two cysteine-rich domains, each of which binds a single GalNAc-4-SO4, are required to form stable complexes with LH. The kinetics of LH binding by the GalNAc-4-SO4-specific form of the receptor in conjunction with its rate of internalization from the cell surface make it likely that only two of the four terminal GalNAc-4-SO4 moieties present on native LH are engaged before receptor internalization. As a result, the rate of hormone clearance will remain constant over a wide range of LH concentrations and will not be sensitive to variations in the number of terminal GalNAc-4-SO4 moieties as long as two or more are present on multiple oligosaccharides.

Viral mediated expression of insulin-like growth factor I blocks the aging-related loss of skeletal muscle function

During the aging process, mammals lose up to a third of their skeletal muscle mass and strength. Although the mechanisms underlying this loss are not entirely understood, we attempted to moderate the loss by increasing the regenerative capacity of muscle. This involved the injection of a recombinant adeno-associated virus directing overexpression of insulin-like growth factor I (IGF-I) in differentiated muscle fibers. We demonstrate that the IGF-I expression promotes an average increase of 15% in muscle mass and a 14% increase in strength in young adult mice, and remarkably, prevents aging-related muscle changes in old adult mice, resulting in a 27% increase in strength as compared with uninjected old muscles. Muscle mass and fiber type distributions were maintained at levels similar to those in young adults. We propose that these effects are primarily due to stimulation of muscle regeneration via the activation of satellite cells by IGF-I. This supports the hypothesis that the primary cause of aging-related impairment of muscle function is a cumulative failure to repair damage sustained during muscle utilization. Our results suggest that gene transfer of IGF-I into muscle could form the basis of a human gene therapy for preventing the loss of muscle function associated with aging and may be of benefit in diseases where the rate of damage to skeletal muscle is accelerated.

Multi-responsiveness of single anterior pituitary cells to hypothalamic-releasing hormones: A cellular basis for paradoxical secretion

The classic view for hypothalamic regulation of anterior pituitary (AP) hormone secretion holds that release of each AP hormone is controlled specifically by a corresponding hypothalamic-releasing hormone (HRH). In this scenario, binding of a given HRH (thyrotropin-, growth hormone-, corticotropin-, and luteinizing hormone-releasing hormones) to specific receptors in its target cell increases the concentration of cytosolic Ca2+ ([Ca2+]i), thereby selectively stimulating the release of the appropriate hormone. However, “paradoxical” responses of AP cells to the four well-established HRHs have been observed repeatedly with both in vivo and in vitro systems, raising the possibility of functional overlap between the different AP cell types. To explore this possibility, we evaluated the effects of HRHs on [Ca2+]i in single AP cells identified immunocytochemically by the hormone they stored. We found that each of the five major AP cell types contained discrete subpopulations that were able to respond to several HRHs. The relative abundance of these multi-responsive cells was 59% for lactotropes, 33% for thyrotropes, and in the range of 47–55% for gonadotropes, corticotropes, and somatotropes. Analysis of prolactin release from single living cells revealed that each of the four HRHs tested were able to induce hormone release from a discrete lactotrope subpopulation, the size of which corresponded closely to that in which [Ca2+]i changes were induced by the same secretagogues. When viewed as a whole, our diverse functional measurements of multi-responsiveness suggest that hypothalamic control of pituitary function is more complicated than previously envisioned. Moreover, they provide a cellular basis for the so-called “paradoxical” behavior of pituitary cells to hypothalamic hypophysiotropic agents.

Defective trafficking and function of KATP channels caused by a sulfonylurea receptor 1 mutation associated with persistent hyperinsulinemic hypoglycemia of infancy

The ATP-sensitive potassium channel (KATP) regulates insulin secretion in pancreatic β cells. Loss of functional KATP channels because of mutations in either the SUR1 or Kir6.2 channel subunit causes persistent hyperinsulinemic hypoglycemia of infancy (PHHI). We investigated the molecular mechanism by which a single phenylalanine deletion in SUR1 (ΔF1388) causes PHHI. Previous studies have shown that coexpression of ΔF1388 SUR1 with Kir6.2 results in no channel activity. We demonstrate here that the lack of functional expression is due to failure of the mutant channel to traffic to the cell surface. Trafficking of KATP channels requires that the endoplasmic reticulum-retention signal, RKR, present in both SUR1 and Kir6.2, be shielded during channel assembly. To ask whether ΔF1388 SUR1 forms functional channels with Kir6.2, we inactivated the RKR signal in ΔF1388 SUR1 by mutation to AAA (ΔF1388 SUR1AAA). Inactivation of similar endoplasmic reticulum-retention signals in the cystic fibrosis transmembrane conductance regulator has been shown to partially overcome the trafficking defect of a cystic fibrosis transmembrane conductance regulator mutation, ΔF508. We found that coexpression of ΔF1388 SUR1AAA with Kir6.2 led to partial surface expression of the mutant channel. Moreover, mutant channels were active. Compared with wild-type channels, the mutant channels have reduced ATP sensitivity and do not respond to stimulation by MgADP or diazoxide. The RKR → AAA mutation alone has no effect on channel properties. Our results establish defective trafficking of KATP channels as a molecular basis of PHHI and show that F1388 in SUR1 is critical for normal trafficking and function of KATP channels.