Prediction of protein structure

LRC: A new algorithm for prediction of conformational B-cell epitopes using statistical approach and clustering method

Identifying of B-cell epitopes from antigen is a challenging task in bioinformatics and applied in vaccine design and drug development. Recently, several methods have been presented to predict epitopes. The physicochemical or structural properties are used by these methods. In this paper, we propose a more appropriate epitope prediction method, LRC, that is based on a combination of physicochemical and structural properties. First, we construct a graph from the surface of antigen, then by using the logistic regression, we model the physicochemical and structural properties and weight each vertex of the graph. Finally, we utilize a clustering method, MCL, to cluster the graph. The effectiveness of the proposed method is benchmarked using several antibody-antigen PDB complexes. The results of LRC algorithm are compared with other methods (DiscoTope, SEPPA and Ellipro) in terms of sensitivity, specificity and other well-known measures. Results indicate that applying the LRC algorithm improves the precision of prediction epitopes in comparison with the mentioned methods. Our LRC program and supplementary material are freely available from http://bs.ipm.ir/softwares/LRC/.

RECOGNITION OF TRANSMEMBRANE SEGMENTS IN PROTEINS: REVIEW AND CONSISTENCY-BASED BENCHMARKING OF INTERNET SERVERS

Membrane proteins perform a number of crucial functions as transporters, receptors, and components of enzyme complexes. Identification of membrane proteins and prediction of their topology is thus an important part of genome annotation. We present here an overview of transmembrane segments in protein sequences, summarize data from large-scale genome studies, and report results of benchmarking of several popular internet servers.

Genetic, immunological and biochemical evidence for a Rnf complex in the acetogen Acetobacterium woodii

Acetogenic bacteria grow by the oxidation of various substrates coupled to the reduction of carbon dioxide (acetogenesis) or other electron acceptors but the mechanisms of energy conservation are still enigmatic. Here, we report the presence of a rnf gene cluster rnfCDGEAB in Acetobacterium woodii that is speculated to encode a novel, energy-conserving ferredoxin:NAD(+)-oxidoreductase complex composed of at least six different subunits. Transcriptional analysis revealed that the genes constitute an operon. RnfC and RnfG were heterologously produced and antibodies were generated. Western blot analyses demonstrated that these subunits were produced and are associated with the cytoplasmic membrane. The subunits were present in cells respiring with either carbon dioxide or caffeate. A preparation with NADH dehydrogenase activity was obtained from detergent solubilized membranes that contained RnfC and RnfG.

Structure of the genetic code suggested by the hydropathy correlation between anticodons and amino acid residues

The correlation between hydropathies of anticodons and amino acids, detected by other authors utilizing scales of amino acid molecules in solution, was improved with the utilization of scales of amino acid residues in proteins. Three partitions were discerned in the correlation plot with the principal dinucleotides of anticodons (pDiN, excluding the wobble position). (a) The set of outliers of the correlation: Gly-CC, Pro-GG, Ser-GA and Ser-CU. The amino acids are consistently small, hydro-apathetic, stabilizers of protein N-ends, preferred in aperiodic protein conformations and belong to synthetases class II. The pDiN sequences are representative of the homogeneous sector (triplets NRR and NYY), distinguished from the mixed sector (triplets NRY and NYR), that depict a 70% correspondence to the synthetases class II and I, respectively. The triplet pairs proposed to be responsible for the coherence in the set of outliers are of the palindromic kind, where the lateral bases are the same, CCC: GGG and AGA: UCU. This suggests that UCU previously belonged to Ser, adding to other indications that the attribution of Arg to YCU was due to an expansion of the Arg-tRNA synthetase specificity. The other attributions produced two correlation sets. (b) One corresponds to the remaining pDiN of the homogeneous sector, containing both synthetase classes; its regression line overlapped the one formed by the remaining attributions to class II. (c) The other contains the pDiN of the mixed sector and produced steeper slopes, especially with the class I attributions. It is suggested that the correlation was established when the amino acid composition of the protein synthetases became progressively enriched and that the set of outliers were the earliest to have been fixed.

Isolation and characterization of a haploid germ cell specific sperm associated antigen 9 (SPAG9) from the baboon

Previously, we cloned and sequenced a sperm specific antigen, designated as HSS (EMBL nomenclature human sperm associated antigen 9: hSPAG9) from human testis (Shankar et al.: Biochem Biophys Res Commun 243:561-565, 1998). The present study was conducted to isolate baboon proteomic homologue in order to find out whether the baboon can provide a suitable model for examining its immunocontraception effects. Baboon SPAG9 (bSPAG9) was cloned and sequenced from the baboon testis cDNA library. The baboon cDNA contained open reading frame encoding 760 amino acids. A 90.6 and 96.8% homology between baboon and human SPAG9 was found at protein and DNA levels. Analysis for tissue specificity by Northern blot procedure using various baboon tissues indicated that bSPAG9 was specifically expressed only in the baboon testis. Further, cell type expression analysis by in situ hybridization in baboon testis demonstrated the expression of bSPAG9 mRNA transcript only in the round spermatid suggesting haploid germ cell expression. Anti-human SPAG9 antibodies recognized the acrosomal compartment region of baboon sperm in indirect immunofluorescence (IIF). Flow cytometry analysis showed surface localization of bSPAG9 in live baboon sperm. The amino acid sequence data for nonhuman primate SPAG9 suggest that antibodies generated by vaccinating baboon with hSPAG9 will recognize nonhuman primate SPAG9, supporting the testing of SPAG9 contraceptive vaccine based on hSPAG9 in the nonhuman primate model.

Comparative analysis of amino acid distributions in integral membrane proteins from 107 genomes

We have performed a comparative analysis of amino acid distributions in predicted integral membrane proteins from a total of 107 genomes. A procedure for identification of membrane spanning helices was optimized on a homology-reduced data set of 170 multi-spanning membrane proteins with experimentally determined topologies. The optimized method was then used for extraction of highly reliable partial topologies from all predicted membrane proteins in each genome, and the average biases in amino acid distributions between loops on opposite sides of the membrane were calculated. The results strongly support the notion that a biased distribution of Lys and Arg residues between cytoplasmic and extra-cytoplasmic segments (the positive-inside rule) is present in most if not all organisms.

Molecular cloning and characterization of the macaque sperm associated antigen 9 (SPAG9): an orthologue of human SPAG9 gene

The present study was conducted to isolate macaque proteomic homologue of human SPAG9 (EMBL nomenclature human sperm associated antigen 9: hSPAG9; Shankar et al., 1998: Biochem Biophys Res Commun 243:561-565) in order to find out whether the macaque can provide a suitable model for examining its immunocontraception effects. Macaque SPAG9 was cloned and sequenced from the macaque testis cDNA library. The macaque cDNA contained open reading frame encoding 712 amino acids. A 84.9% and 94% homology between macaque and human SPAG9 was found at protein and DNA levels. Northern analysis and RNA in situ hybridization experiments revealed testis- and stage-specific expressions of macaque SPAG9 mRNA, mainly confined to round spermatid suggesting haploid germ cell expression. Anti-human SPAG9 antibodies recognized native SPAG9 in macaque sperm extract in Western blotting and the acrosomal compartment region of macaque sperm in indirect immunofluorescence. Flow cytometry analysis further revealed surface localization of macaque SPAG9 in live macaque sperm. The amino acid sequence data for nonhuman primate SPAG9 suggest that antibodies generated by vaccinating macaque with hSPAG9 will recognize nonhuman primate SPAG9, supporting the testing of SPAG9 contraceptive vaccine based on hSPAG9 in the nonhuman primate model.

Mutation in mitochondrial complex I ND6 subunit is associated with defective response to hypoxia in human glioma cells

BACKGROUND

Hypoxia-tolerant human glioma cells reduce oxygen consumption rate in response to oxygen deficit, a defense mechanism that contributes to survival under moderately hypoxic conditions. In contrast, hypoxia-sensitive cells lack this ability. As it has been previously shown that hypoxia-tolerant (M006x, M006xLo, M059K) and -sensitive (M010b) glioma cells express differences in mitochondrial function, we investigated whether mitochondrial DNA-encoded mutations are associated with differences in the initial response to oxygen deficit.

RESULTS

The mitochondrial genome was sequenced and 23 mtDNA alterations were identified, one of which was an unreported mutation (T-C transition in base pair 14634) in the hypoxia-sensitive cell line, M010b, that resulted in a single amino acid change in the gene encoding the ND6 subunit of NADH:ubiquinone oxidoreductase (Complex I). The T14634C mutation did not abrogate ND6 protein expression, however, M010b cells were more resistant to rotenone, an agent used to screen for Complex I mutations, and adriamycin, an agent activated by redox cycling. The specific function of mtDNA-encoded, membrane-embedded Complex I ND subunits is not known at present. Current models suggest that the transmembrane arm of Complex I may serve as a conformationally driven proton channel. As cellular respiration is regulated, in part, by proton flux, we used homology-based modeling and computational molecular biology to predict the 3D structure of the wild type and mutated ND6 proteins. These models predict that the T14634C mutation alters the structure and orientation of the trans-membrane helices of the ND6 protein.

CONCLUSION

Complex I ND subunits are mutational hot spots in tumor mtDNA. Genetic changes that alter Complex I structure and function may alter a cell's ability to respond to oxygen deficit and consolidate hypoxia rescue mechanisms, and may contribute to resistance to chemotherapeutic agents that require redox cycling for activation.

Role of the sequence surrounding predicted transmembrane helix M4 in membrane association and function of the Ca(2+) release channel of skeletal muscle sarcoplasmic reticulum (ryanodine receptor isoform 1)

The role of the sequence surrounding M4 in ryanodine receptors (RyR) in membrane association and function was investigated. This sequence contains a basic, 19-amino acid M3/M4 loop, a hydrophobic 44-49 amino acid sequence designated M4 (or M4a/M4b), and a hydrophilic M4/M5 loop. Enhanced green fluorescent protein (EGFP) was inserted into RyR1 and truncated just after the basic sequence, just after M4, within the M4/M5 loop, just before M5 and just after M5. The A52 epitope was inserted into RyR2 and truncated just after M4a. Analysis of these constructs ruled out a M3/M4 transmembrane hairpin and narrowed the region of membrane association to M4a/M4b. EGFP inserted between M4a and M4b in full-length RyR2 was altered conformationally, losing fluorescence and gaining trypsin sensitivity. Although it was accessible to an antibody from the cytosolic side, tryptic fragments were membrane-bound. The expressed protein containing EGFP retained caffeine-induced Ca(2+) release channel function. These results suggest that M4a/M4b either forms a transmembrane hairpin or associates in an unorthodox fashion with the cytosolic leaflet of the membrane, possibly involving the basic M3/M4 loop. The expression of a mutant RyR1, Delta4274-4535, deleted in the sequence surrounding both M3 and M4, restored robust, voltage-gated L-type Ca(2+) currents and Ca(2+) transients in dyspedic myotubes, demonstrating that this sequence is not required for either orthograde (DHPR activation of sarcoplasmic reticulum Ca(2+) release) or retrograde (RyR1 increase in DHPR Ca(2+) channel activity) signals of excitation-contraction coupling. Maximal amplitudes of L-currents and Ca(2+) transients with Delta4274-4535 were larger than with wild-type RyR1, and voltage-gated sarcoplasmic reticulum Ca(2+) release was more sensitive to activation by sarcolemmal voltage sensors. Thus, this region may act as a negative regulatory module that increases the energy barrier for Ca(2+) release channel opening.

Membrane topology of the ZntB efflux system of Salmonella enterica serovar Typhimurium

The membrane topology of the ZntB Zn(2+) transport protein of Salmonella enterica serovar Typhimurium was determined by constructing deletion derivatives of the protein and genetically fusing them to blaM or lacZ cassettes. The enzymatic activities of the hybrid proteins indicate that ZntB is a bitopic integral membrane protein consisting largely of two independent domains. The first 266 amino acids form a large, highly charged domain within the cytoplasm, while the remaining 61 residues form a small membrane domain containing two membrane-spanning segments. The overall orientation towards the cytoplasm is consistent with the ability of ZntB to facilitate zinc efflux.

Geometrical analysis of Cys-Cys bridges in proteins and their prediction from incomplete structural information

Analysis of C-alpha atom positions from cysteines involved in disulphide bridges in protein crystals shows that their geometric characteristics are unique with respect to other Cys-Cys, non-bridging pairs. They may be used for predicting disulphide connections in incompletely determined protein structures, such as low resolution crystallography or theoretical folding experiments. The basic unit for analysis and prediction is the 3 x 3 distance matrix for Cx positions of residues (i - 1), Cys(i), (i +1) with (j - 1), Cys(j), (j + 1). In each of its columns, row and diagonal vector--outer distances are larger than the central distance. This analysis is compared with some analytical models.

MUP-4 is a novel transmembrane protein with functions in epithelial cell adhesion in Caenorhabditis elegans

Tissue functions and mechanical coupling of cells must be integrated throughout development. A striking example of this coupling is the interactions of body wall muscle and hypodermal cells in Caenorhabditis elegans. These tissues are intimately associated in development and their interactions generate structures that provide a continuous mechanical link to transmit muscle forces across the hypodermis to the cuticle. Previously, we established that mup-4 is essential in embryonic epithelial (hypodermal) morphogenesis and maintenance of muscle position. Here, we report that mup-4 encodes a novel transmembrane protein that is required for attachments between the apical epithelial surface and the cuticular matrix. Its extracellular domain includes epidermal growth factor-like repeats, a von Willebrand factor A domain, and two sea urchin enterokinase modules. Its intracellular domain is homologous to filaggrin, an intermediate filament (IF)-associated protein that regulates IF compaction and that has not previously been reported as part of a junctional complex. MUP-4 colocalizes with epithelial hemidesmosomes overlying body wall muscles, beginning at the time of embryonic cuticle maturation, as well as with other sites of mechanical coupling. These findings support that MUP-4 is a junctional protein that functions in IF tethering, cell-matrix adherence, and mechanical coupling of tissues.

Class II antimicrobial peptides from lactic acid bacteria

Strains of lactic acid bacteria (LAB) produce a wide variety of antibacterial peptides. More than fifty of these so-called peptide bacteriocins have been isolated in the last few years. They contain 20-60 amino acids, and are cationic and hydrophobic in nature. Several of these bacteriocins consist of two complementary peptides. The peptide bacteriocins of LAB are inhibitory at concentrations in the nanomolar range, and cause membrane permeabilization and leakage of intracellular components in sensitive cells. The inhibitory spectrum is limited to gram-positive bacteria, and in many cases to bacteria closely related to the producing strain. Among the target organisms are food spoilage bacteria and pathogens such as Listeria, so that many of these antimicrobial peptides could have a potential as food preservatives as well as in medical applications.

The complete genome sequence of PM2, the first lipid-containing bacterial virus To Be isolated

Bacteriophage PM2 was isolated from the Pacific Ocean off the coast of Chile in the late 1960s. It was a new virus type, later classified as Corticoviridae, and also the first bacterial virus for which it was demonstrated that lipids are part of the virion structure. Here we report the determination and analysis of the 10, 079-bp circular dsDNA genome sequence. Noteworthy discoveries are the replication initiation system, which is related to the rolling circle mechanism described for phages such as φX174 and P2, and a 1.2-kb sequence that is similar to the maintenance region of a plasmid found in a marine Pseudoalteromonas sp. strain A28.

SOR1, a gene required for photosensitizer and singlet oxygen resistance in Cercospora fungi, is highly conserved in divergent organisms

Filamentous Cercospora fungi are resistant to photosensitizing compounds that generate singlet oxygen. C. nicotianae photosensitizer-sensitive mutants were restored to full resistance by transformation with SOR1 (Singlet Oxygen Resistance 1), a gene recovered from a wild-type genomic library. SOR1 null mutants generated via targeted gene replacement confirmed the requirement for SOR1 in photosensitizer resistance. SOR1 RNA is present throughout the growth cycle. Although resistance to singlet oxygen is rare in biological systems, SOR1, a gene with demonstrated activity against singlet-oxygen-generating photosensitizers, is highly conserved in organisms from widely diverse taxa. The characterization of SOR1 provides an additional phenotype to this large group of evolutionarily conserved genes.

Protein SRP54 of human signal recognition particle: cloning, expression, and comparative analysis of functional sites

Signal recognition particle (SRP) plays a critical role in the targeting of secretory proteins to cellular membranes. An essential component of SRP is the protein SRP54, which interacts not only with the nascent signal peptide, but also with the SRP RNA. To understand better how protein targeting occurs in the human system, the human SRP54 gene was cloned, sequenced, and the protein was expressed in bacteria and insect cells. Recombinant SRP54 was purified from both sources. The protein bound to SRP RNA in the presence of protein SRP19, and associated with the signal peptide of in vitro translated pre-prolactin. Comparative sequence analysis of human SRP54 with homologs from all three phylogenetic domains was combined with high-stringency protein secondary structure prediction. A conserved RNA-binding loop was predicted in the largely helical M-domain of SRP54. Contrary to general belief, the unusually high number of methionine residues clustered outside the predicted helices, thus indicating a mechanism of signal peptide recognition that may involve methionine-rich loops.

Periplakin, a novel component of cornified envelopes and desmosomes that belongs to the plakin family and forms complexes with envoplakin

The cornified envelope is a layer of transglutaminase cross-linked protein that is assembled under the plasma membrane of keratinocytes in the outermost layers of the epidermis. We have determined the cDNA sequence of one of the proteins that becomes incorporated into the cornified envelope of cultured epidermal keratinocytes, a protein with an apparent molecular mass of 195 kD that is encoded by a mRNA with an estimated size of 6.3 kb. The protein is expressed in keratinizing and nonkeratinizing stratified squamous epithelia and in a number of other epithelia. Expression of the protein is upregulated during the terminal differentiation of epidermal keratinocytes in vivo and in culture. Immunogold electron microscopy was used to demonstrate an association of the 195-kD protein with the desmosomal plaque and with keratin filaments in the differentiated layers of the epidermis. Sequence analysis showed that the 195-kD protein is a member of the plakin family of proteins, to which envoplakin, desmoplakin, bullous pemphigoid antigen 1, and plectin belong. Envoplakin and the 195-kD protein coimmunoprecipitate. Analysis of their rod domain sequences suggests that the formation of both homodimers and heterodimers would be energetically favorable. Confocal immunofluorescent microscopy of cultured epidermal keratinocytes revealed that envoplakin and the 195-kD protein form a network radiating from desmosomes, and we speculate that the two proteins may provide a scaffolding onto which the cornified envelope is assembled. We propose to name the 195-kD protein periplakin.

Identification of an RNA-binding-loop in the N-terminal region of signal-recognition-particle protein SRP19

Protein SRP19 is a 144-amino-acid polypeptide that associates intimately with the signal-recognition particle RNA (SRP RNA) and serves as an important structural and functional component of the SRP. We investigated the structure and RNA-binding activity of the human SRP19 protein by the use of comparative sequence analysis, high-stringency structure prediction, proteolytic susceptibility, and site-directed mutagenesis. SRP19 was found to consist of two distinct regions (called N-terminal and C-terminal regions) that are separated by a boundary of approximately 12-15 amino acid residues. Both regions contain an alpha-helix and several beta-strands that are connected by loops or turns. In agreement with the hypothetical model, proteolytic susceptibility demonstrated the predominant accessibility of two sites: one in a surface loop of the N-terminal region (YLNNKKTIAEGR33), and another site in the C-terminal tail at residues L129 and E133. The RNA-binding activities of mutant polypeptides with changes of conserved lysines and arginines (mutants K27Q, R33Q and R34Q) demonstrated that the proteolytically accessible loop of the N-terminal region is in direct contact with the SRP RNA. In contrast, alteration of a certain basic amino acid residues in the C-terminal region (R83, K116 and R118), as well as a deletion of four amino acid residues located at the boundary between the two regions, had no effect on the RNA-binding ability. The structural model that emerges from our data is thematically similar to that of ribosomal protein S5, the N-domain of which contains a loop motif believed to interact with double-stranded RNA. The presence of a similar structural feature in protein SRP19 has significant implications for the structure and function of the SRP19-RNA complex.

The rkpGHI and -J genes are involved in capsular polysaccharide production by Rhizobium meliloti

The first complementation unit of the fix-23 region of Rhizobium meliloti, which comprises six genes (rkpAB-CDEF) exhibiting similarity to fatty acid synthase genes, is required for the production of a novel type of capsular polysaccharide that is involved in root nodule development and structurally analogous to group II K antigens found in Escherichia coli (G. Petrovics, P. Putnoky, R. Reuhs, J. Kim, T. A. Thorp, K. D. Noel, R. W. Carlson, and A. Kondorosi, Mol. Microbiol. 8:1083-1094, 1993; B. L. Reuhs, R. W. Carlson, and J. S. Kim, J. Bacteriol. 175:3570-3580, 1993). Here we present the nucleotide sequence for the other three complementation units of the fix-23 locus, revealing the presence of four additional open reading frames assigned to genes rkpGHI and -J. The putative RkpG protein shares similarity with acyltransferases, RkpH is homologous to short-chain alcohol dehydrogenases, and RkpJ shows significant sequence identity with bacterial polysaccharide transport proteins, such as KpsS of E. coli. No significant homology was found for RkpI. Biochemical and immunological analysis of Tn5 derivatives for each gene demonstrated partial or complete loss of capsular polysaccharides from the cell surface; on this basis, we suggest that all genes in the fix-23 region are required for K-antigen synthesis or transport.

Envoplakin, a novel precursor of the cornified envelope that has homology to desmoplakin

The cornified envelope is a layer of transglutaminase cross-linked protein that is deposited under the plasma membrane of keratinocytes in the outermost layers of the epidermis. We present the sequence of one of the cornified envelope precursors, a protein with an apparent molecular mass of 210 kD. The 210-kD protein is translated from a 6.5-kb mRNA that is transcribed from a single copy gene. The mRNA was upregulated during suspension-induced terminal differentiation of cultured human keratinocytes. Like other envelope precursors, the 210-kD protein became insoluble in SDS and beta-mercaptoethanol on activation of transglutaminases in cultured keratinocytes. The protein was expressed in keratinizing and nonkeratinizing stratified squamous epithelia, but not in simple epithelia or nonepithelial cells. Immunofluorescence staining showed that in epidermal keratinocytes, both in vivo and in culture, the protein was upregulated during terminal differentiation and partially colocalized with desmosomal proteins. Immunogold EM confirmed the colocalization of the 210-kD protein and desmoplakin at desmosomes and on keratin filaments throughout the differentiated layers of the epidermis. Sequence analysis showed that the 210-kD protein is homologous to the keratin-binding proteins desmoplakin, bullous pemphigoid antigen 1, and plectin. These data suggest that the 210-kD protein may link the cornified envelope to desmosomes and keratin filaments. We propose that the 210-kD protein be named "envoplakin."

Cl- channels in basolateral renal medullary vesicles. X. Cloning of a Cl- channel from rabbit outer medulla

These experiments were intended to identify candidate cDNAs which might encode basolateral membrane Cl- channels of the mTAL using a homology-based cloning strategy. We prepared a cDNA library using a 1.8 to 3.2 kb mRNA fraction from rabbit outer medulla that induces a Cl- conductance in cellular membranes of Xenopus laevis oocytes. The cDNA library was screened with two 32P-oligonucleotide probes corresponding to highly conserved sequences in other Cl- channels. We isolated two cDNAs: rbClC-Ka and rbClC-Kb. The protein sequences deduced from these two cDNAs had 99% homology. Using RT-PCR technology, cultured mouse mTAL cells were found to contain mRNA corresponding to those two cDNAs. Expression of the mRNAs corresponding to these two cDNAs was kidney-specific and was greater in rabbit renal medulla than rabbit renal cortex. Finally, by using RT-PCR technology in combination with microdissected glomeruli or tubule segments, we found mRNA for rbClC-Ka in glomeruli, proximal convoluted tubules, mTAL and cortical collecting tubules.

GenStar: a method for de novo drug design

A novel method, which we call GenStar, has been developed to suggest chemically reasonable structures which fill the active sites of enzymes. The proposed molecules provide good steric contact with the enzyme and exist in low-energy conformations. These structures are composed entirely of sp3 carbons which are grown sequentially, but which can also branch or form rings. User-selected enzyme seed atoms may be used to determine the area in which structure generation begins. Alternatively, GenStar may begin with a predocked 'inhibitor core' from which atoms are grown. For each new atom generated by the program, several hundred candidate positions representing a range of reasonable bond lengths, bond angles, and torsion angles are considered. Each of these candidates is scored, based on a simple enzyme contact model. The selected position is chosen at random from among the highest scoring cases. Duplicate structures may be removed using a variety of criteria. The compounds may be energy minimized and displayed using standard modeling programs. Also, it is possible to analyze the collection of all structures created by GenStar and locate binding motifs for common fragments such as benzene and naphthylene. Tests of the method using HIV protease, FK506 binding protein (FKBP-12) and human carbonic anhydrase (HCA-II) demonstrated that structures similar to known potent inhibitors may be generated with GenStar.

Effects of calcium on recombinant bovine chromogranin A

Bovine chromogranin A, the acidic calcium-binding protein characteristic of endocrine secretory vesicles, has been expressed in Escherichia coli using the pET3a vector system under T7 polymerase control. The expressed protein is located in the bacterial cytosol and can be purified from bacterial proteins by a heat treatment step, followed by gel filtration, anion-exchange, and reversed-phase chromatography. The purified recombinant chromogranin A has an apparent M(r) of ca. 72,000 by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, in spite of its 432-amino acid polypeptide chain, consistent with observations on natural chromogranin A. The primary structure has been confirmed by mass spectral analysis of tryptic peptides, by Edman degradation of the intact protein, and by immunoreactivity with sequence-specific antibodies. Analysis by circular dichroism spectroscopy shows pH- and concentration-dependent spectra. The spectra are Ca2(+)-dependent from 5 to 40 microM.

Structural domains of phytochrome deduced from homologies in amino acid sequences

A method of semiempirical identification of structural domains is proposed. The procedure is based on the comparison of amino acid sequences in groups of homologous proteins. This approach was tested using 32 known protein sequences from different cytochrome b5, cytochrome c, lysozyme, hemoglobin, and myoglobin proteins. The method presented was able to identify all structural domains of these reference proteins. A consensus secondary structure provided information on structural content of these domains predicting correctly 21 of 23 (91%) of alpha-helices. We applied this method to six homologous phytochrome sequences from Avena, Arabadopsis, Cucurbita, Maize, Oryza, and Pisum. Some of the identified domains can be assigned to the known tertiary structure categories. For example, an alpha/beta domain is localized in the region known to stabilize the phytochrome chromophore in the red light absorbing form (Pr). One alpha-helical and one alpha/beta domains are localized in regions important for the chromophore stabilization in the far-red absorbing form (Pfr). From an analysis of noncovalent interaction patterns in another domain it is proposed that a phytochrome dimer contact involves two segments localized between residues 730 and 821 (using numbering of aligned sequences). Also, a possible antiparallel beta-sheet structure of this region has been suggested. According to this model, the long axis of the interacting structures is perpendicular to a twofold symmetry axis of the phytochrome dimer.

Structure of the reovirus cell-attachment protein: a model for the domain organization of sigma 1

This report describes a model for the structure of the reovirus cell-attachment protein sigma 1. S1 gene nucleotide sequences were determined for prototype strains of the three serotypes of mammalian reoviruses. Deduced amino acid sequences of the S1-encoded sigma 1 proteins were then compared in order to identify conserved features of these sequences. Discrete regions in the amino-terminal two-thirds of sigma 1 sequence share characteristics with the fibrous domains of other cellular and viral proteins. Most of the amino-terminal one-third of sigma 1 sequence is predicted to form an alpha-helical coiled coil like that of myosin. The middle one-third of sigma 1 sequence appears more heterogeneous; it is predicted to form a large region of beta-sheet that is followed by a region which contains two short alpha-helical coiled coils separated by a smaller region of beta-sheet. The two beta-sheet regions are each proposed to form a cross-beta sandwich like that suggested for the rod domain of the adenovirus fiber protein (N. M. Green, N. G. Wrigley, W. C. Russell, S. R. Martin, and A. D. McLachlan, EMBO J. 2:1357-1365, 1983). The remaining carboxy-terminal one-third of sigma 1 sequence is predicted to form a structurally complex globular domain. A model is suggested in which the discrete regions of sigma 1 sequence are ascribed to morphologic regions seen in computer-processed electron micrographic images of the protein (R. D. B. Fraser, D. B. Furlong, B. L. Trus, M. L. Nibert, B. N. Fields, and A. C. Steven, J. Virol. 64:2990-3000, 1990.

Secondary structure of the entomocidal toxin from Bacillus thuringiensis subsp. kurstaki HD-73

The secondary structure of the toxin from Bacillus thuringiensis subsp. kurstaki (Btk) HD-73 was estimated by Raman, infrared, and circular dichroism spectroscopy, and by predictive methods. Circular dichroism and infrared spectroscopy gave an estimate of 33-40% alpha-helix, whereas Raman and predictive methods gave approximately 20%. Raman and circular dichroism spectra, as well as predictive methods, indicated that the toxin contains 32-40% beta-sheet structure, whereas infrared spectroscopy gave a slightly lower estimate. Thus, all of these approaches are in agreement that the native conformation of Btk HD-73 toxin is highly folded and contains considerable amounts of both alpha-helical and beta-sheet structures. No significant differences were detected in the secondary structure of the toxin either in solution or as a hydrated pellet.

Prediction of the secondary structure of myosin light chains from comparison of homologous sequences. Implications for the interaction between myosin heavy and light chains

The primary sequences of seventeen essential and seventeen regulatory myosin light chains were analyzed and compared, using algorithms based on the different structural properties of their amino acid residues. This process allowed estimation of the structural homology between the proteins studied, and improved the prediction of their mean secondary structure and functionally important segments or residues. On the basis of the crystal structure of troponin C, a model of the myosin essential light chain with a fairly compact form is proposed. The possible sites of interaction between myosin light and heavy chains from rabbit skeletal muscle were also investigated by a complementarity method adapted to helix-rich proteins. Segments 139-149 and 65-75 in the essential light chain and segments 27-37, 67-77 and 97-107 in the regulatory light chain are suggested to constitute some of these sites, as most of them were found to have the features of surface-seeking helices.

A structural and dynamic model for the nicotinic acetylcholine receptor

Folding of the five polypeptide subunits (alpha 2 beta gamma delta) of the nicotinic acetylcholine receptor (AChR) into a functional structural model is described. The principles used to arrange the sequences into a structure include: (1) hydrophobicity----membrane-crossing segments; (2) amphipathic character----ion-carrying segments (ion channel with single group rotations); (3) molecular shape (elongated, pentagonal cylinder)----folding dimensions of exobilayer portion; (4) choice of acetylcholine binding sites----specific folding of exobilayer segments; (5) location of reducible disulfides (near agonist binding site)----additional specification of exobilayer arrangement; (6) genetic homology----consistency of functional group choices; (7) noncompetitive antagonist labeling----arrangement of bilayer helices. The AChR model is divided into three parts: (a) exobilayer consisting of 11 antiparallel beta-strands from each subunit; (b) bilayer consisting of four hydrophobic and one amphiphilic alpha-helix from each subunit; (c) cytoplasmic consisting of one (folded) loop from each subunit. The exobilayer strands can form a closed 'flower' (the 'resting state') which is opened ('activated') by agonists bound perpendicular to the strands. Rearrangement of the agonists to a strand-parallel position and partial closing of the 'flower' leads to a desensitized receptor. The actions of acetylcholine and succinoyl and suberoyl bis-cholines are clarified by the model. The opening and closing of the exobilayer 'flower' controls access to the ion channel which is composed of the five amphiphilic bilayer helices. A molecular mechanism for ion flow in the channel is given. Openings interrupted by short duration closings (50 microseconds) depend upon channel group motions. The unusual photolabeling of intrabilayer serines in alpha, beta and delta subunits but not in gamma subunits near the binding site for non-competitive antagonists is explained along with a mechanism for the action of these antagonists such as phencyclidine. The unusual alpha 192Cys-193Cys disulfide may have a special peptide arrangement, such as a cis-peptide bond to a following proline (G.A. Petsko and E.M. Kosower, unpublished results). The position of phosphorylatable sites and proline-rich segments are noted for the cytoplasmic loops. The dynamic behavior of the AChR channel and many different experimental results can be interpreted in terms of the model. An example is the lowering of ionic conductivity on substitution of bovine for Torpedo delta M2 segment. The model represents a useful construct for the design of experiments on AChR.

An integrated prediction of secondary, tertiary and quaternary structure of glucose dehydrogenase

Based on homology of partial sequences, on physico-chemical evidence and on studies using chemical modification, we came to the tentative conclusion that tetrameric glucose dehydrogenases from Bacillus megaterium and B. subtilis should have a structure closely related to that of lactate dehydrogenase. The overall homology of primary structures was found to be very low, however, and independent predictions of secondary structure produced a clearly different pattern of beta-strands and alpha-helices. We nevertheless tried a manual prediction based on the hydrophobic nature of internal beta-sheet and on the amphiphilic character of external helices. This treatment led to the identification of analogues of all the beta-strands present in lactate dehydrogenase with the exception of beta C. Six amphiphilic helices were identified corresponding to alpha B, alpha C, alpha D, alpha 1F, alpha 2F and alpha 3G in lactate dehydrogenase. Conserved functional residues were found at analogous positions. The Q and R intersubunit contacts could be identified and partial proteolysis was found to occur on the outer surface of the tetramer. The structure was found to explain the better binding of NADP as compared to NAD+ and offered a rationalization of the role of the essential lysine at position 201.