Characterization and sequence prediction of structural variations in α-helix

Background

Centralised resources such as GenBank and UniProt are perfect examples of the major international efforts that have been made to integrate and share biological information. However, additional data that adds value to these resources needs a simple and rapid route to public access. The Distributed Annotation System (DAS) provides an adequate environment to integrate genomic and proteomic information from multiple sources, making this information accessible to the community. DAS offers a way to distribute and access information but it does not provide domain experts with the mechanisms to participate in the curation process of the available biological entities and their annotations.

Results

We designed and developed a Collaborative Annotation System for proteins called DAS Writeback. DAS writeback is a protocol extension of DAS to provide the functionalities of adding, editing and deleting annotations. We implemented this new specification as extensions of both a DAS server and a DAS client. The architecture was designed with the involvement of the DAS community and it was improved after performing usability experiments emulating a real annotation task.

Conclusions

We demonstrate that DAS Writeback is effective, usable and will provide the appropriate environment for the creation and evolution of community protein annotation.     

Accurate prediction of the burial status of transmembrane residues of α-helix membrane protein by incorporating the structural and physicochemical features

Predicting the burial status (the residue exposure to the lipid bilayer or buried within the protein core) of transmembrane (TM) residues of α-helix membrane protein (αHMP) is of great importance for genome-wide annotation and for experimental researchers to elucidate diverse physiological processes. In this work, we developed a new computational model that can be used for predicting the burial status of TM residues of αHMP. By incorporating physicochemical scales and conservation index, an efficient prediction model using least squares support vector machine (LS-SVM) was developed. The model was developed from 43 protein chains and its prediction ability was evaluated by an independent test set of other non-redundant ten protein chains. The prediction accuracy of our method was much better than the results of the reported works. Our results demonstrate that the LS-SVM prediction model incorporating structural and physicochemical features derived from sequence information could greatly improve the prediction accuracy.     

Characterization and strain distribution of four alleles at the hemoglobin α-chain structural locus in the mouse

In a survey of mice from 40 inbred strains, largely previously untested, four alleles were distinguished at the Hba locus, determining structure of adult mouse hemoglobin chain. This finding supports and extends previous sequence studies by others. Methods are given by which each phenotype was characterized by its solubility profile at varying pH and by its chromatography pattern. Concordance was complete between histidine-positive T-4 (defining Hba ^c) and high-intermediate solubility profile. In three inbred strains, a distinct new low-solubility profile, not associated with his-positive T-4, was recognized, and mice with this phenotype were classified as Hba ^d. Implications of observed widely distributed four-allele polymorphism of mouse hemoglobin -chain structure are discussed.This investigation was supported in part by NIH research grant CA-01074 from the National Cancer Institute and in part by the Virginia and D. K. Ludwig Foundation. The Jackson Laboratory is fully accredited by the American Association for Accreditation of Laboratory Animal Care.     

Hydrophobic organization of α-helix membrane bundle in bacteriorhodopsin

The hydrophobic organization of the intramembraneα-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM)α-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices—neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TMα-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.     

Modelling of transmembrane α-helix bundles

A strategy for modelling transmembrane -helix bundles has been investigated. Results concerning the rotational orientations of the helices are described and perspectives for extensions of the method are discussed.     

Thermal denaturation and structural changes of α-helical proteins in keratins

To gain insight into the thermal stability of intermediate filaments and matrix in the biological composite structure of α-keratins, the thermal denaturation performance of human hair fibers was investigated by Modulated Differential Scanning Calorimetry (MDSC) in the dry and the wet state. Denaturation enthalpy ΔH(D) in water was found to be independent of heating rate (11.5J/g) and to be approximately double as high as in the dry state (5.2J/g). The lower enthalpy (dry) and its dependency on heating rate are attributed to effects of pyrolysis. The stepwise change of reversing heat capacity ΔC(p) marks the denaturation process as a classic two-stage transition. The increase of ΔC(p) with heating rate reflects a continuous shift of the nature of the denaturation of the α-helical material, first, into random coil and then towards random β-materials for lower heating rates. Denaturation temperatures follow Arrhenius relationships with heating rate, yielding activation energies of 416kJ/mol (dry) and 263kJ/mol (wet), respectively. A decrease of activation energy (wet) for high heating rates supports the hypothesis of systematic changes of the pathway of denaturation.     

Quantifying sequence and structural features of protein–RNA interactions

Increasing awareness of the importance of protein–RNA interactions has motivated many approaches to predict residue-level RNA binding sites in proteins based on sequence or structural characteristics. Sequence-based predictors are usually high in sensitivity but low in specificity; conversely structure-based predictors tend to have high specificity, but lower sensitivity. Here we quantified the contribution of both sequence- and structure-based features as indicators of RNA-binding propensity using a machine-learning approach. In order to capture structural information for proteins without a known structure, we used homology modeling to extract the relevant structural features. Several novel and modified features enhanced the accuracy of residue-level RNA-binding propensity beyond what has been reported previously, including by meta-prediction servers. These features include: hidden Markov model-based evolutionary conservation, surface deformations based on the Laplacian norm formalism, and relative solvent accessibility partitioned into backbone and side chain contributions. We constructed a web server called aaRNA that implements the proposed method and demonstrate its use in identifying putative RNA binding sites.     

Sequence preference of α-helix N-terminal tetrapeptide

The α-helix is the most abundant secondary structure in proteins. Due to the specific i, i+4 hydrogen bond pattern, the two termini have unsatisfied hydrogen bonds, and are less constrained; in order to compensate for this, specific residues are preferred for the terminal positions. However, a naive combination of the statistically-preferred residues for each position may not result in a stable N-terminal helical sequence. In order to provide a set of preferable N-terminal peptides for α-helix design, we have studied the N-terminal tetrapeptide sequence motifs that are favorable for helix formation using statistical analysis and atomistic simulations. A set of tetrapeptide sequences including TEEE and TPEE were found to be favorable motifs. In addition to forming more hydrogen bonds in the helical conformation, the favorable motifs also tended to form more capping boxes. To empirically test our predictions, we obtained 10 peptides with different N-terminal motifs and measured their α-helical content by circular dichroism spectroscopy. The experimental results agreed qualitatively with the statistical and simulation results. Furthermore, some of the suggested preferable tetrapeptide sequences have been successfully applied in de novo protein design.     

Influence of minor backbone structural variations in modulating the in vitro gene transfer efficacies of α-tocopherol based cationic transfection lipids

Toward probing the influence of backbone structural variation in cationic lipid mediated gene delivery of α-tocopherol based lipids, two novel α-tocopherol based lipids 1 and 2 have been designed and synthesized. The only structural difference between the cationic amphiphiles 1 and 2 is the backbone structure, where lipid 1 has a non-glycerol backbone and lipid 2 has a glycerol backbone. The lipids 1 and 2 showed contrasting transfection efficiencies: lipid 1 showed high gene transfer efficacy in multiple cultured animals cell lines, whereas lipid 2 is transfection incompetent. In summary, the present findings demonstrate that in the case of α-tocopherol based lipids even minor structural variations like backbone can profoundly influence size, DNA binding characteristics, cellular uptake, and consequently gene delivery efficacies.     

Isolation and partial sequence of goat spleen prothymosin α

Goat prothymosin , a highly acidic polypeptide of pl 3.5, 109 amino acid residues, has been isolated from lymphoid and non-lymphoid tissues of young female goats. Unlike rat, murine and porcine prothymosins , goat prothymosin appears at a higher concentration in the spleen compared with the thymus. The sequence of segments of the polypeptide involving known mutations has been determined, by automatic sequencing of its tryptic peptide fragments. The acidic amino acid-rich segment in the middle of the molecule, including residues 49–83, has not been sequenced. Goat prothymosin closely resembles bovine prothymosin , with only one substitution, proline for alanine at position 85. It also resembles human prothymosin , with only three substitutions. It differs more significantly from rat and murine prothymosins , by two deletions and three substitutions. The results show the highly conserved nature of the molecule, with substitutions at given positions only.Abbreviations ProT Prothymosin - T₁ Thymosin ₁ - MLR Mixed Lymphocyte Response - HPLC High Performance Liquid Chromatography - RIA Radioimmunoassay - B Aspartic acid or Asparagine - Z Glutamic acid or Glutamine     

Role of the structural context in selection of hydrophobic side-chain rotamers in a- and d-positions of α-helices

It was demonstrated for the first time that the distribution of side-chain rotamers in the a-and d-positions of α-helices of coiled-coil (cc) proteins follows a certain trend, rather then being random. For instance, most side chains adopt t rotamers in the a-positions and g^? rotamers in the d-positions of helical dimers. Vice versa, most side chains adopt g^? rotamers in the a-positions and t rotamers in the d-positions of tetramers. It was concluded that selection of the side-chain rotamers depends on the packing of α-helices and, consequently, depends on the structural context.     

Studies on molecular evolution and structural features of double-headed inhibitors of α-amylase and trypsin in plants

Plant seeds usually have high concentrations of proteinase and amylase inhibitors. These inhibitors exhibit a wide range of specificity, stability and oligomeric structure. In this communication, we report analysis of sequences that show statistically significant similarity to the double-headed α-amylase/trypsin inhibitor of ragi (Eleusine coracana). Our aim is to understand their evolutionary and structural features. The 14 sequences of this family that are available in the SWISSPROT database form three evolutionarily distinct branches. The branches relate to enzyme specificities and also probably to the oligomeric state of the proteins and not to the botanical class of the plant from which the enzymes are derived. This suggests that the enzyme specificities of the inhibitors evolved before the divergence of commercially cultivated cereals. The inhibitor sequences have three regions that display periodicity in hydrophobicity. It is likely that this feature reflects extended secondary structure in these segments. One of the most variable regions of the polypeptide corresponds to a loop, which is most probably exposed in the native structure of the inhibitors and is responsible for the inhibitory property.     

Characterization of Inhibitor-Bound α-Synuclein Dimer: Role of α-Synuclein N-Terminal Region in Dimerization and Inhibitor Binding

α-Synuclein is a major component of filamentous inclusions that are histological hallmarks of Parkinson's disease and other α-synucleinopathies. Previous analyses have revealed that several polyphenols inhibit α-synuclein assembly with low micromolar IC₅₀ values, and that SDS-stable, noncytotoxic soluble α-synuclein oligomers are formed in their presence. Structural elucidation of inhibitor-bound α-synuclein oligomers is obviously required for the better understanding of the inhibitory mechanism. In order to characterize inhibitor-bound α-synucleins in detail, we have prepared α-synuclein dimers in the presence of polyphenol inhibitors, exifone, gossypetin, and dopamine, and purified the products. Peptide mapping and mass spectrometric analysis revealed that exifone-treated α-synuclein monomer and dimer were oxidized at all four methionine residues of α-synuclein. Immunoblot analysis and redox-cycling staining of endoproteinase Asp-N-digested products showed that the N-terminal region (1-60) is involved in the dimerization and exifone binding of α-synuclein. Ultra-high-field NMR analysis of inhibitor-bound α-synuclein dimers showed that the signals derived from the N-terminal region of α-synuclein exhibited line broadening, confirming that the N-terminal region is involved in inhibitor-induced dimerization. The C-terminal portion still predominantly exhibited the random-coil character observed in monomeric α-synuclein. We propose that the N-terminal region of α-synuclein plays a key role in the formation of α-synuclein assemblies.     

Fold-recognition and comparative modeling of human α2,3-sialyltransferases reveal their sequence and structural similarities to CstII from Campylobacter jejuni

Background  

The 3-D structure of none of the eukaryotic sialyltransferases (SiaTs) has been determined so far. Sequence alignment algorithms such as BLAST and PSI-BLAST could not detect a homolog of these enzymes from the protein databank. SiaTs, thus, belong to the hard/medium target category in the CASP experiments. The objective of the current work is to model the 3-D structures of human SiaTs which transfer the sialic acid in α2,3-linkage viz., ST3Gal I, II, III, IV, V, and VI, using fold-recognition and comparative modeling methods. The pair-wise sequence similarity among these six enzymes ranges from 41 to 63%.     

Functional characterization of amphipathic α-helix in the osmoregulatory ABC transporter OpuA

The ATP-binding-cassette transporter OpuA from Lactococcus lactis is composed of two ATPase subunits (OpuAA) and two subunits (OpuABC) with the transmembrane domain fused to an extracellular substrate-binding protein. Of the almost 1900 homologues of OpuA known to date, a subset has an amino-terminal amphipathic helix (plus extra transmembrane segment) fused to the core of the transmembrane domain of the OpuABC subunit. FRET measurements indicate that the amphipathic α-helix is located close to the membrane surface, where its hydrophobic face interacts with the transport protein rather than the membrane lipids. Next, we determined the functional role of this accessory region by engineering the amphipathic α-helix. We analyzed the consequence of the mutations in intact cells by monitoring growth and transport of glycine betaine under normal and osmotic stress conditions. More detailed studies were performed in hybrid membrane vesicles, proteoliposomes, and bilayer nanodisks. We show that the amphipathic α-helix of OpuA is necessary for high activity of OpuA but is not critical for the biogenesis of the protein or the ionic regulation of transport.     

cDNA sequence and chromosome localization of pig α1,3 galactosyltransferase

Human serum contains natural antibodies (NAb), which can bind to endothelial cell surface antigens of other mammals. This is believed to be the major initiating event in the process of hyperacute rejection of pig to primate xenografts. Recent work has implicated galoctosyl 1,3 galactosyl 1,4 N-acetyl-glucosaminyl carbohydrate epitopes, on the surface of pig endothelial cells as a major target of human natural antibodies. This epitope is made by a specific galactosyltransferase (1,3 GT) present in pigs but not in higher primates. We have now cloned and sequenced a full-length pig 1,3 GT cDNA. The predicted 371 amino acid protein sequence shares 85% and 76% identity with previously characterized cattle and mouse 1,3 GT protein sequences, respectively. By using fluorescence and isotopic in situ hybridization, the GGTA1 gene was mapped to the region q2.10–q2.11 of pig chromosome 1, providing further evidence of homology between the subterminal region of pig chromosome 1q and human chromosome 9q, which harbors the locus encoding the AB0 blood group system, as well as a human pseudogene homologous to the pig GGTA1 gene.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number L36152     

Characterization of the mutant α-mannosidase in bovine mannosidosis

Residual acidic α-mannosidase, varying in amount up to approx. 15% of normal values, can be measured in various organs of a calf with mannosidosis. The highest specific activity and relative proportion of residual activity were found in the liver. Chromatography on DEAE-cellulose showed that the residual activity was associated with two components, which were eluted at comparable positions with those found in normal tissues. The residual activity had a lower thermal stability and a higher K_m value for a synthetic substrate than did the normal enzyme. No differences in molecular weight or electrophoretic mobility between normal acidic α-mannosidase and the residual activity were observed by gel filtration and electrophoresis on cellulose acetate respectively. The isoelectric focusing profiles for the α-mannosidase in the normal and pathological livers were very similar. It is suggested that a mutant enzyme, resulting from a mutation in a structural gene, accounts for the residual acidic α-mannosidase in mannosidosis. The mutant enzyme, which cross-reacts with antiserum raised against normal bovine acidic α-mannosidase, is present at a decreased concentration compared with the normal enzyme. There is a correlation between the concentrations of residual activity and cross-reacting material in mannosidosis. α-Mannosidase with a pH optimum of 5.75 and which is activated by Zn²⁺ was also detected in the liver of the calf with mannosidosis. However, it is probably not a product of the defective gene because addition of Zn²⁺ indicated that it was also present in normal tissues.     

Characterization of the α-cyanocinnamate binding site in rat heart mitochondria and in submitochondrial particles

The effect of pH and substrates on the binding of radiolabelled α-cyanocinnamate to mitochondria and submitochondrial particles has been investigated. It has been found that the binding is strongly influenced by the pH of the medium (it decreases on increasing the pH of the medium). The inhibition of pyruvate oxidation by this inhibitor follows the same pH dependence. The pH affects only the affinity of the α-cyanocinnamate binding site without changing their total number. A similar pH dependence has been found in inside-out submitochondrial particles where the binding sites are directly accessible. The quantitative parameters of the binding of α-cyanocinnamate in submitochondrial particles have been determined. The binding can be prevented or displaced by pyruvate and other substrates of the carrier. The turnover number for pyruvate transport in rat-heart mitochondria has been determined.