BIND—The Biomolecular Interaction Network Database

The Biomolecular Interaction Network Database (BIND; http://binddb.org) is a database designed to store full descriptions of interactions, molecular complexes and pathways. Development of the BIND 2.0 data model has led to the incorporation of virtually all components of molecular mechanisms including interactions between any two molecules composed of proteins, nucleic acids and small molecules. Chemical reactions, photochemical activation and conformational changes can also be described. Everything from small molecule biochemistry to signal transduction is abstracted in such a way that graph theory methods may be applied for data mining. The database can be used to study networks of interactions, to map pathways across taxonomic branches and to generate information for kinetic simulations. BIND anticipates the coming large influx of interaction information from high-throughput proteomics efforts including detailed information about post-translational modifications from mass spectrometry. Version 2.0 of the BIND data model is discussed as well as implementation, content and the open nature of the BIND project. The BIND data specification is available as ASN.1 and XML DTD.

Polar residues drive association of polyleucine transmembrane helices

Although many polar residues are directly involved in transmembrane protein functions, the extent to which they contribute to more general structural features is still unclear. Previous studies have demonstrated that asparagine residues can drive transmembrane helix association through interhelical hydrogen bonding [Choma, C., Gratkowski, H., Lear, J. D. & DeGrado, W. F. (2000) Nat. Struct. Biol. 7, 161–166; and Zhou, F. X., Cocco, M. J., Russ, W. P., Brunger, A. T. & Engelman, D. M. (2000) Nat. Struct. Biol. 7, 154–160]. We have studied the ability of other polar residues to promote helix association in detergent micelles and in biological membranes. Our results show that polyleucine sequences with Asn, Asp, Gln, Glu, and His, residues capable of being simultaneously hydrogen bond donors and acceptors, form homo- or heterooligomers. In contrast, polyleucine sequences with Ser, Thr, and Tyr do not associate more than the polyleucine sequence alone. The results therefore provide experimental evidence that interactions between polar residues in the helices of transmembrane proteins may serve to provide structural stability and oligomerization specificity. Furthermore, such interactions can allow structural flexibility required for the function of some membrane proteins.

The immunoprotective MHC II epitope of a chemically induced tumor harbors a unique mutation in a ribosomal protein

CD4+ T lymphocyte clones, generated from mice immunized with the methylcholanthrene-induced fibrosarcoma Meth A (H-2d), are restricted by I-Ed and recognize a unique antigen on Meth A. The antigen has been purified and characterized as the ribosomal protein L11. The antigenic epitope is contained within the sequence EYELRKHNFSDTG and is generated by substitution of Asn by His (italic) caused by a single point mutation. The tumor contains the wild-type and the mutated alleles. Immunization of BALB/cJ mice with the mutated epitope but not with the wild-type epitope protects mice against a subsequent challenge with the Meth A sarcoma. Adoptive transfer of CD4+ clones into BALB/c mice renders the mice specifically resistant to Meth A sarcoma. The mutated L11 epitope is thus shown to be an immunoprotective epitope in vivo by several criteria.

Prediction of protein deamidation rates from primary and three-dimensional structure

A method for the quantitative estimation of instability with respect to deamidation of the asparaginyl (Asn) residues in proteins is described. The procedure involves the observation of several simple aspects of the three-dimensional environment of each Asn residue in the protein and a calculation that includes these observations, the primary amino acid residue sequence, and the previously reported complete set of sequence-dependent rates of deamidation for Asn pentapeptides. This method is demonstrated and evaluated for 23 proteins in which 31 unstable and 167 stable Asn residues have been reported and for 7 unstable and 63 stable Asn residues that have been reported in 61 human hemoglobin variants. The relative importance of primary structure and three-dimensional structure in Asn deamidation is estimated.

Ca2+-binding activity of a COOH-terminal fragment of the Drosophila BK channel involved in Ca2+-dependent activation

Mutational and biophysical analysis suggests that an intracellular COOH-terminal domain of the large conductance Ca2+-activated K+ channel (BK channel) contains Ca2+-binding site(s) that are allosterically coupled to channel opening. However the structural basis of Ca2+ binding to BK channels is unknown. To pursue this question, we overexpressed the COOH-terminal 280 residues of the Drosophila slowpoke BK channel (Dslo-C280) as a FLAG- and His6-tagged protein in Escherichia coli. We purified Dslo-C280 in soluble form and used a 45Ca2+-overlay protein blot assay to detect Ca2+ binding. Dslo-C280 exhibits specific binding of 45Ca2+ in comparison with various control proteins and known EF-hand Ca2+-binding proteins. A mutation (D5N5) of Dslo-C280, in which five consecutive Asp residues of the “Ca-bowl” motif are changed to Asn, reduces 45Ca2+-binding activity by 56%. By electrophysiological assay, the corresponding D5N5 mutant of the Drosophila BK channel expressed in HEK293 cells exhibits lower Ca2+ sensitivity for activation and a shift of ≈+80 mV in the midpoint voltage for activation. This effect is associated with a decrease in the Hill coefficient (N) for activation by Ca2+ and a reduction in apparent Ca2+ affinity, suggesting the loss of one Ca2+-binding site per monomer. These results demonstrate a functional correlation between Ca2+ binding to a specific region of the BK protein and Ca2+-dependent activation, thus providing a biochemical approach to study this process.

Mutations within a furin consensus sequence block proteolytic release of ectodysplasin-A and cause X-linked hypohidrotic ectodermal dysplasia

X-linked hypohidrotic ectodermal dysplasia (XLHED) is a heritable disorder of the ED-1 gene disrupting the morphogenesis of ectodermal structures. The ED-1 gene product, ectodysplasin-A (EDA), is a tumor necrosis factor (TNF) family member and is synthesized as a membrane-anchored precursor protein with the TNF core motif located in the C-terminal domain. The stalk region of EDA contains the sequence -Arg-Val-Arg-Arg156-Asn-Lys-Arg159-, representing overlapping consensus cleavage sites (Arg-X-Lys/Arg-Arg↓) for the proprotein convertase furin. Missense mutations in four of the five basic residues within this sequence account for ≈20% of all known XLHED cases, with mutations occurring most frequently at Arg156, which is shared by the two consensus furin sites. These analyses suggest that cleavage at the furin site(s) in the stalk region is required for the EDA-mediated cell-to-cell signaling that regulates the morphogenesis of ectodermal appendages. Here we show that the 50-kDa EDA parent molecule is cleaved at -Arg156Asn-Lys-Arg159↓- to release the soluble C-terminal fragment containing the TNF core domain. This cleavage appears to be catalyzed by furin, as release of the TNF domain was blocked either by expression of the furin inhibitor α1-PDX or by expression of EDA in furin-deficient LoVo cells. These results demonstrate that mutation of a functional furin cleavage site in a developmental signaling molecule is a basis for human disease (XLHED) and raise the possibility that furin cleavage may regulate the ability of EDA to act as a juxtacrine or paracrine factor.

Induction of a Carbon-Starvation-Related Proteolysis in Whole Maize Plants Submitted to Light/Dark Cycles and to Extended Darkness1

Three-week-old maize (Zea mays L.) plants were submitted to light/dark cycles and to prolonged darkness to investigate the occurrence of sugar-limitation effects in different parts of the whole plant. Soluble sugars fluctuated with light/dark cycles and dropped sharply during extended darkness. Significant decreases in protein level were observed after prolonged darkness in mature roots, root tips, and young leaves. Glutamine and asparagine (Asn) changed in opposite ways, with Asn increasing in the dark. After prolonged darkness the increase in Asn accounted for most of the nitrogen released by protein breakdown. Using polyclonal antibodies against a vacuolar root protease previously described (F. James, R. Brouquisse, C. Suire, A. Pradet, P. Raymond [1996] Biochem J 320: 283–292) or the 20S proteasome, we showed that the increase in proteolytic activities was related to an enrichment of roots in the vacuolar protease, with no change in the amount of 20S proteasome in either roots or leaves. Our results show that no significant net proteolysis is induced in any part of the plant during normal light/dark cycles, although changes in metabolism and growth appear soon after the beginning of the dark period, and starvation-related proteolysis probably appears in prolonged darkness earlier in sink than in mature tissues.

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Electric-field-induced Schiff-base deprotonation in D85N mutant bacteriorhodopsin.

The application of an external electric field to dry films of Asp-85-->Asn mutant bacteriorhodopsin causes deprotonation of the Schiff base, resulting in a shift of the optical absorption maximum from 600 nm to 400 nm. This is in marked contrast to the case of wild-type bacteriorhodopsin films, in which electric fields produce a red-shifted product whose optical properties are similar to those of the acid-blue form of the protein. This difference is due to the much weaker binding of the Schiff-base proton in the mutant protein, as indicated by its low pK of approximately 9, as compared with the value pK approximately 13 in the wild type. Other bacteriorhodopsins with lowered Schiff-base pK values should also exhibit a field-induced shift in the protonation equilibrium of the Schiff base. We propose mechanisms to account for these observations.

A carbonic anhydrase from the nacreous layer in oyster pearls.

It is believed that the polymorphism observed in calcium carbonate crystals, such as aragonite and calcite in mollusk shells, is controlled by organic matrix proteins secreted from the mantle epithelia. However, the fine structures of these proteins are still unknown, and to understand the molecular mechanisms of mineralization process, detailed structural analyses of the organic matrix proteins are essential. For this, we have carried out purification, characterization, and cDNA cloning of nacrein, which is a soluble organic matrix protein in the nacreous layer of oyster pearls. Northern blot analysis showed that the nacrein transcript was specifically expressed in mantle pallial. Analysis of the deduced amino acid sequence revealed that the protein contained two functional domains: one was a carbonic anhydrase and another was a Gly-Xaa-Asn (Xaa = Asp, Asn, or Glu) repeat domain; however, the carbonic anhydrase domain was split into two subdomains with insertion of the Gly-Xaa-Asn repeat domain between them. Our findings suggest that nacrein actually functions as a matrix protein whose repeated Gly-Xaa-Asn domain possibly binds calcium and as a carbonic anhydrase that catalyzes the HCO3- formation, thus participating in calcium carbonate crystal formation of the nacreous layer.

Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element.

NM23-H2, a presumed regulator of tumor metastasis in humans, is a hexameric protein with both enzymatic (NDP kinase) and regulatory (transcriptional activation) activity. While the structure and catalytic mechanisms have been well characterized, the mode of DNA binding is not known. We examined this latter function in a site-directed mutational study and identified residues and domains essential for the recognition of a c-myc regulatory sequence. Three amino acids, Arg-34, Asn-69, and Lys-135, were found among 30 possibilities to be critical for DNA binding. Two of these, Asn-69 and Lys-135, are not conserved between NM23 variants differing in DNA-binding potential, suggesting that DNA recognition resides partly in nonconserved amino acids. All three DNA-binding defective mutant proteins are active enzymatically and appear to be stable hexamers, suggesting that they perform at the level of DNA recognition and that separate functional domains exist for enzyme catalysis and DNA binding. In the context of the known crystal structure of NM23-H2, the DNA-binding residues are located within distinct structural motifs in the monomer, which are exposed to the surface near the 2-fold axis of adjacent subunits in the hexamer. These findings are explained by a model in which NM23-H2 binds DNA with a combinatorial surface consisting of the "outer" face of the dimer. Chemical crosslinking data support a dimeric DNA-binding mode by NM23-H2.

Selenocysteine, identified as the penultimate C-terminal residue in human T-cell thioredoxin reductase, corresponds to TGA in the human placental gene.

The possible relationship of selenium to immunological function which has been suggested for decades was investigated in studies on selenium metabolism in human T cells. One of the major 75Se-labeled selenoproteins detected was purified to homogeneity and shown to be a homodimer of 55-kDa subunits. Each subunit contained about 1 FAD and at least 0.74 Se. This protein proved to be thioredoxin reductase (TR) on the basis of its catalytic activities, cross-reactivity with anti-rat liver TR antibodies, and sequence identities of several tryptic peptides with the published deduced sequence of human placental TR. Physicochemical characteristics of T-cell TR were similar to those of a selenocysteine (Secys)-containing TR recently isolated from human lung adenocarcinoma cells. The sequence of a 12-residue 75Se-labeled tryptic peptide from T-cell TR was identical with a C-terminal-deduced sequence of human placental TR except that Secys was present in the position corresponding to TGA, previously thought to be the termination codon, and this was followed by Gly-499, the actual C-terminal amino acid. The presence of the unusual conserved Cys-Secys-Gly sequence at the C terminus of TR in addition to the redox active cysteines of the Cys-Val-Asn-Val-Gly-Cys motif in the FAD-binding region may account for the peroxidase activity and the relatively low substrate specificity of mammalian TRs. The finding that T-cell TR is a selenoenzyme that contains Se in a conserved C-terminal region provides another example of the role of selenium in a major antioxidant enzyme system (i.e., thioredoxin-thioredoxin reductase), in addition to the well-known glutathione peroxidase enzyme system.

Potentiation of proton transfer function by electrostatic interactions in photosynthetic reaction centers from Rhodobacter sphaeroides: First results from site-directed mutation of the H subunit.

The x-ray crystallographic structure of the photosynthetic reaction center (RC) has proven critical in understanding biological electron transfer processes. By contrast, understanding of intraprotein proton transfer is easily lost in the immense richness of the details. In the RC of Rhodobacter (Rb.) sphaeroides, the secondary quinone (QB) is surrounded by amino acid residues of the L subunit and some buried water molecules, with M- and H-subunit residues also close by. The effects of site-directed mutagenesis upon RC turnover and quinone function have implicated several L-subunit residues in proton delivery to QB, although some species differences exist. In wild-type Rb. sphaeroides, Glu L212 and Asp L213 represent an inner shell of residues of particular importance in proton transfer to QB. Asp L213 is crucial for delivery of the first proton, coupled to transfer of the second electron, while Glu L212, possibly together with Asp L213, is necessary for delivery of the second proton, after the second electron transfer. We report here the first study, by site-directed mutagenesis, of the role of the H subunit in QB function. Glu H173, one of a cluster of strongly interacting residues near QB, including Asp L213, was altered to Gln. In isolated mutant RCs, the kinetics of the first electron transfer, leading to formation of the semiquinone, QB-, and the proton-linked second electron transfer, leading to the formation of fully reduced quinol, were both greatly retarded, as observed previously in the Asp L213 --> Asn mutant. However, the first electron transfer equilibrium, QA-QB <==> QAQB-, was decreased, which is opposite to the effect of the Asp L213 --> Asn mutation. These major disruptions of events coupled to proton delivery to QB were largely reversed by the addition of azide (N3-). The results support a major role for electrostatic interactions between charged groups in determining the protonation state of certain entities, thereby controlling the rate of the second electron transfer. It is suggested that the essential electrostatic effect may be to "potentiate" proton transfer activity by raising the pK of functional entities that actually transfer protons in a coupled fashion with the second electron transfer. Candidates include buried water (H3O+) and Ser L223 (serine-OH2+), which is very close to the O5 carbonyl of the quinone.

In vivo selection of basic region-leucine zipper proteins with altered DNA-binding specificities.

A transcription interference assay was used to generate mutant basic region-leucine zipper proteins with altered DNA-binding specificities. A library of mutants of a CCAAT/enhancer binding protein was constructed by randomizing five DNA-contacting amino acids in the basic region Asn-18, Ala-15, Val-14, Ser-11, and Arg-10. These mutants were then selected for their ability to bind mutant recognition sequences containing substitutions at the 2 and 3 positions of the wild-type sequence 5'-A5T4T3G2C1G1'C2'A3A4'T5'-3'. Mutants containing the sequence Leu-18Tyr-15Xaa-14Tyr-11Arg-10, in which four of the five contact residues are altered, were identified that recognize the palindromic sequence 5'-ATCYCGY'GAT-3' (Xaa = asparagine when Y = G; Xaa = methionine when Y = A). Moreover, in a selection against the sequence 5'-ATTACGTAAT-3', mutants were obtained containing substitutions not only in the basic region but also in the hinge region between the basic and leucine zipper regions. The mutant proteins showed high specificity in a functional transcription interference assay. A model for the interaction of these mutants with the target DNA sequences is discussed.

An integrated map of the genome of the tubercle bacillus, Mycobacterium tuberculosis H37Rv, and comparison with Mycobacterium leprae.

An integrated map of the genome of the tubercle bacillus, Mycobacterium tuberculosis, was constructed by using a twin-pronged approach. Pulsed-field gel electrophoretic analysis enabled cleavage sites for Asn I and Dra I to be positioned on the 4.4-Mb circular chromosome, while, in parallel, clones from two cosmid libraries were ordered into contigs by means of fingerprinting and hybridization mapping. The resultant contig map was readily correlated with the physical map of the genome via the landmarked restriction sites. Over 165 genes and markers were localized on the integrated map, thus enabling comparisons with the leprosy bacillus, Mycobacterium leprae, to be undertaken. Mycobacterial genomes appear to have evolved as mosaic structures since extended segments with conserved gene order and organization are interspersed with different flanking regions. Repetitive sequences and insertion elements are highly abundant in M. tuberculosis, but the distribution of IS6110 is apparently nonrandom.