Protein conformational prediction

Substrates of the chloroplast small heat shock proteins 22E/F point to thermolability as a regulative switch for heat acclimation in Chlamydomonas reinhardtii

We have identified 39 proteins that interact directly or indirectly with high confidence with chloroplast HSP22E/F under heat stress thus revealing chloroplast processes affected by heat. Under conditions promoting protein unfolding, small heat shock proteins (sHsps) prevent the irreversible aggregation of unfolding proteins by integrating into forming aggregates. Aggregates containing sHsps facilitate the access of Hsp70 and ClpB/Hsp104 chaperones, which in ATP-dependent reactions disentangle individual proteins from the aggregates and assist in their refolding to the native state. Chlamydomonas reinhardtii encodes eight different sHsps (HSP22A to H). The goal of this work was to identify chloroplast-targeted sHsps in Chlamydomonas and to obtain a comprehensive list of the substrates with which they interact during heat stress in order to understand which chloroplast processes are disturbed under heat stress. We show that HSP22E and HSP22F are major chloroplast-targeted sHsps that have emerged from a recent gene duplication event resulting from the ongoing diversification of sHsps in the Volvocales. HSP22E/F strongly accumulate during heat stress and form high molecular mass complexes. Using differential immunoprecipitation, mass spectrometry and a stringent filtering algorithm we identified 39 proteins that with high-confidence interact directly or indirectly with HSP22E/F under heat stress. We propose that the apparent thermolability of several of these proteins might be a desired trait as part of a mechanism enabling Chlamydomonas chloroplasts to rapidly react to thermal stress.

Structural and Functional Characterization of the VQ Protein Family and VQ Protein Variants from Soybean

Proteins containing the FxxxVQxhTG or VQ motif interact with WRKY transcription factors. Although VQ proteins have been reported in several plants, knowledge about their structures, functions and evolution is still very limited. Here, we report structural and functional analysis of the VQ protein family from soybean. Like Arabidopsis homologues, soybean VQ proteins bind only Group I and IIc WRKY proteins and a substantial number of their genes are responsive to stress-associated phytohormones. Overexpression of some soybean VQ genes in Arabidopsis had strong effects on plant growth, development, disease resistance and heat tolerance. Phylogenetic analysis, sequence alignment and site-directed mutagenesis revealed that the region immediately upstream of the FxxxVQxhTG motif also affects binding to WRKY proteins. Consistent with a larger WRKY-binding VQ domain, soybean VQ22 protein from cultivated soybean contains a 4-amino acid deletion in the region preceding its VQ motif that completely abolishes its binding to WRKY proteins. By contrast, the 4-amino acid deletion is absent in the VQ22 protein from wild soybean species (Glycine soja). Overexpression of wild soybean VQ22 in cultivated soybean inhibited growth, particularly after cold treatment. Thus, the mutation of soybean VQ22 is associated with advantageous phenotypes and may have been positively selected during evolution.

The N-terminal tropomyosin- and actin-binding sites are important for leiomodin 2's function

Leiomodin is a potent actin nucleator related to tropomodulin, a capping protein localized at the pointed end of the thin filaments. Mutations in leiomodin-3 are associated with lethal nemaline myopathy in humans, and leiomodin-2-knockout mice present with dilated cardiomyopathy. The arrangement of the N-terminal actin- and tropomyosin-binding sites in leiomodin is contradictory and functionally not well understood. Using one-dimensional nuclear magnetic resonance and the pointed-end actin polymerization assay, we find that leiomodin-2, a major cardiac isoform, has an N-terminal actin-binding site located within residues 43-90. Moreover, for the first time, we obtain evidence that there are additional interactions with actin within residues 124-201. Here we establish that leiomodin interacts with only one tropomyosin molecule, and this is the only site of interaction between leiomodin and tropomyosin. Introduction of mutations in both actin- and tropomyosin-binding sites of leiomodin affected its localization at the pointed ends of the thin filaments in cardiomyocytes. On the basis of our new findings, we propose a model in which leiomodin regulates actin poly-merization dynamics in myocytes by acting as a leaky cap at thin filament pointed ends.

Comparing SILAC- and stable isotope dimethyl-labeling approaches for quantitative proteomics

Stable isotope labeling is widely used to encode and quantify proteins in mass-spectrometry-based proteomics. We compared metabolic labeling with stable isotope labeling by amino acids in cell culture (SILAC) and chemical labeling by stable isotope dimethyl labeling and find that they have comparable accuracy and quantitative dynamic range in unfractionated proteome analyses and affinity pull-down experiments. Analyzing SILAC- and dimethyl-labeled samples together in single liquid chromatography–mass spectrometric analyses minimizes differences under analytical conditions, allowing comparisons of quantitative errors introduced during sample processing. We find that SILAC is more reproducible than dimethyl labeling. Because proteins from metabolically labeled populations can be combined before proteolytic digestion, SILAC is particularly suited to studies with extensive sample processing, such as fractionation and enrichment of peptides with post-translational modifications. We compared both methods in pull-down experiments using a kinase inhibitor, dasatinib, and tagged GRB2-SH2 protein as affinity baits. We describe a StageTip dimethyl-labeling protocol that we applied to in-solution and in-gel protein digests. Comparing the impact of post-digest isotopic labeling on quantitative accuracy, we demonstrate how specific experimental designs can benefit most from metabolic labeling approaches like SILAC and situations where chemical labeling by stable isotope-dimethyl labeling can be a practical alternative.

Sequence-based predictor of ATP-binding residues using random forest and mRMR-IFS feature selection

We develop a computational and statistical approach (ATPBR) for predicting ATP-binding residues in proteins from amino acid sequences by using random forests with a novel hybrid feature. The hybrid feature incorporates a new feature called PSSMPP, the predicted secondary structure and orthogonal binary vectors. The mRMR-IFS feature selection method is utilized to construct the best prediction model. At last, ATPBR achieves significantly improved performance over existing methods, with 87.53% accuracy and a Matthew׳s correlation coefficient of 0.554. In addition, our further analysis demonstrates that PSSMPP distinguishes more effectively between ATP-binding and non-binding residues. Besides, the optimal features selected by the mRMR-IFS method improve the prediction performance and may provide useful insights for revealing the mechanisms of ATP and proteins interactions.

Phosphorylation of tropomodulin1 contributes to the regulation of actin filament architecture in cardiac muscle

Tropomodulin1 (Tmod1) is an actin-capping protein that plays an important role in actin filament pointed-end dynamics and length in striated muscle. No mechanisms have been identified to explain how Tmod1's functional properties are regulated. The purpose of this investigation was to explore the functional significance of the phosphorylation of Tmod1 at previously identified Thr54. Rat cardiomyocytes were assessed for phosphorylation of Tmod1 using Pro-Q Diamond staining and (32)P labeling. Green fluorescent protein-tagged phosphorylation-mimic (T54E) and phosphorylation-deficient (T54A) versions of Tmod1 were expressed in cultured cardiomyocytes, and the ability of these mutants to assemble and restrict actin lengths was observed. We report for the first time that Tmod1 is phosphorylated endogenously in cardiomyocytes, and phosphorylation at Thr54 causes a significant reduction in the ability of Tmod1 to assemble to the pointed end compared with that of the wild type (WT; 48 vs. 78%, respectively). In addition, overexpression of Tmod1-T54E restricts actin filament lengths by only ∼3%, whereas Tmod1-WT restricts the lengths significantly by ∼8%. Finally, Tmod1-T54E altered the actin filament-capping activity in polymerization assays. Taken together, our data suggest that pointed-end assembly and Tmod1's thin filament length regulatory function are regulated by its phosphorylation state.

Alteration of tropomyosin-binding properties of tropomodulin-1 affects its capping ability and localization in skeletal myocytes

Tropomodulin (Tmod) is an actin-capping protein that binds to the two tropomyosins (TM) at the pointed end of the actin filament to prevent further actin polymerization and depolymerization. Therefore, understanding the role of Tmod is very important when studying actin filament dependent processes such as muscle contraction and intracellular transport. The capping ability of Tmod is highly influenced by TM and is 1000-fold greater in the presence of TM. There are four Tmod isoforms (Tmod1-4), three of which, Tmod1, Tmod3, and Tmod4, are expressed in skeletal muscles. The affinity of Tmod1 to skeletal striated TM (stTM) is higher than that of Tmod3 and Tmod4 to stTM. In this study, we tested mutations in the TM-binding sites of Tmod1, using circular dichroism (CD) and prediction analysis (PONDR). The mutations R11K, D12N, and Q144K were chosen because they decreased the affinity of Tmod1 to stTM, making it similar to that of affinity of Tmod3 and Tmod4 to stTM. Significant reduction of inhibition of actin pointed-end polymerization in the presence of stTM was shown for Tmod1 (R11K/D12N/Q144K) as compared with WT Tmod1. When GFP-Tmod1 and mutants were expressed in primary chicken skeletal myocytes, decreased assembly of Tmod1 mutants was revealed. This indicates a direct correlation between TM-binding and the actin-capping abilities of Tmod. Our data confirmed the hypothesis that assembly of Tmod at the pointed-end of the actin filament depends on its TM-binding affinity.

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.

SNOSite: exploiting maximal dependence decomposition to identify cysteine S-nitrosylation with substrate site specificity

S-nitrosylation, the covalent attachment of a nitric oxide to (NO) the sulfur atom of cysteine, is a selective and reversible protein post-translational modification (PTM) that regulates protein activity, localization, and stability. Despite its implication in the regulation of protein functions and cell signaling, the substrate specificity of cysteine S-nitrosylation remains unknown. Based on a total of 586 experimentally identified S-nitrosylation sites from SNAP/L-cysteine-stimulated mouse endothelial cells, this work presents an informatics investigation on S-nitrosylation sites including structural factors such as the flanking amino acids composition, the accessible surface area (ASA) and physicochemical properties, i.e. positive charge and side chain interaction parameter. Due to the difficulty to obtain the conserved motifs by conventional motif analysis, maximal dependence decomposition (MDD) has been applied to obtain statistically significant conserved motifs. Support vector machine (SVM) is applied to generate predictive model for each MDD-clustered motif. According to five-fold cross-validation, the MDD-clustered SVMs could achieve an accuracy of 0.902, and provides a promising performance in an independent test set. The effectiveness of the model was demonstrated on the correct identification of previously reported S-nitrosylation sites of Bos taurus dimethylarginine dimethylaminohydrolase 1 (DDAH1) and human hemoglobin subunit beta (HBB). Finally, the MDD-clustered model was adopted to construct an effective web-based tool, named SNOSite (http://csb.cse.yzu.edu.tw/SNOSite/), for identifying S-nitrosylation sites on the uncharacterized protein sequences.

Graphical representation and mathematical characterization of protein sequences and applications to viral proteins

Graphical representation and numerical characterization (GRANCH) of nucleotide and protein sequences is a new field that is showing a lot of promise in analysis of such sequences. While formulation and applications of GRANCH techniques for DNA/RNA sequences started just over a decade ago, analyses of protein sequences by these techniques are of more recent origin. The emphasis is still on developing the underlying technique, but significant results have been achieved in using these methods for protein phylogeny, mass spectral data of proteins and protein serum profiles in parasites, toxicoproteomics, determination of different indices for use in QSAR studies, among others. We briefly mention these in this chapter, with some details on protein phylogeny and viral diseases. In particular, we cover a systematic method developed in GRANCH to determine conserved surface exposed peptide segments in selected viral proteins that can be used for drug and vaccine targeting. The new GRANCH techniques and applications for DNAs and proteins are covered briefly to provide an overview to this nascent field.

Transforming omics data into context: Bioinformatics on genomics and proteomics raw data

Differential gene expression analysis and proteomics have exerted significant impact on the elucidation of concerted cellular processes, as simultaneous measurement of hundreds to thousands of individual objects on the level of RNA and protein ensembles became technically feasible. The availability of such data sets has promised a profound understanding of phenomena on an aggregate level, expressed as the phenotypic response (observables) of cells, e.g., in the presence of drugs, or characterization of cells and tissue displaying distinct patho-physiological states. However, the step of transforming these data into context, i.e., linking distinct expression or abundance patterns with phenotypic observables - and furthermore enabling a sound biological interpretation on the level of reaction networks and concerted pathways, is still a major shortcoming. This finding is certainly based on the enormous complexity embedded in cellular reaction networks, but a variety of computational approaches have been developed over the last few years to overcome these issues. This review provides an overview on computational procedures for analysis of genomic and proteomic data introducing a sequential analysis workflow: Explorative statistics for deriving a first, from the purely statistical viewpoint, relevant candidate gene/protein list, followed by co-regulation and network analysis to biologically expand this core list toward functional networks and pathways. The review on these procedures is complemented by example applications tailored at identification of disease-associated proteins. Optimization of computational procedures involved, in conjunction with the continuous increase in additional biological data, clearly has the potential of boosting our understanding of processes on a cell-wide level.

The presence of valine at residue 129 in human prion protein accelerates amyloid formation

The polymorphism at residue 129 of the human PRNP gene modulates disease susceptibility and the clinico-pathological phenotypes in human transmissible spongiform encephalopathies. The molecular mechanisms by which the effect of this polymorphism are mediated remain unclear. It has been shown that the folding, dynamics and stability of the physiological, alpha-helix-rich form of recombinant PrP are not affected by codon 129 polymorphism. Consistent with this, we have recently shown that the kinetics of amyloid formation do not differ between protein containing methionine at codon 129 and valine at codon 129 when the reaction is initiated from the alpha-monomeric PrP(C)-like state. In contrast, we have shown that the misfolding pathway leading to the formation of beta-sheet-rich, soluble oligomer was favoured by the presence of methionine, compared with valine, at position 129. In the present work, we examine the effect of this polymorphism on the kinetics of an alternative misfolding pathway, that of amyloid formation using partially folded PrP allelomorphs. We show that the valine 129 allelomorph forms amyloids with a considerably shorter lag phase than the methionine 129 allelomorph both under spontaneous conditions and when seeded with pre-formed amyloid fibres. Taken together, our studies demonstrate that the effect of the codon 129 polymorphism depends on the specific misfolding pathway and on the initial conformation of the protein. The inverse propensities of the two allelomorphs to misfold in vitro through the alternative oligomeric and amyloidogenic pathways could explain some aspects of prion diseases linked to this polymorphism such as age at onset and disease incubation time.

Homology model of the structure of influenza B virus HA1

Influenza B virus is one of two types of influenza virus that cause substantial morbidity and mortality in humans, the other being influenza A virus. The inability to provide lasting protection to humans against influenza B virus infection is due, in part, to antigenic drift of the viral surface glycoprotein, haemagglutinin (HA). Studies of the antigenicity of the HA of influenza B virus have been hampered by lack of knowledge of its structure. To address this gap, two possible models have been inferred for this structure, based on two known structures of the homologous HA of the influenza A virus (subtypes H3 and H9). Statistical, structural and functional analyses of these models suggested that they matched important details of experimental observations and did not differ from each other in any substantive way. These models were used to investigate two HA sites at which viral variants appeared to carry a selective advantage. It was found that each of these sites coevolved with nearby sites to compensate for either size or charge changes.

Methionine 129 Variant of Human Prion Protein Oligomerizes More Rapidly than the Valine 129 Variant

The human PrP gene (PRNP) has two common alleles that encode either methionine or valine at codon 129. This polymorphism modulates disease susceptibility and phenotype of human transmissible spongiform encyphalopathies, but the molecular mechanism by which these effects are mediated remains unclear. Here, we compared the misfolding pathway that leads to the formation of beta-sheet-rich oligomeric isoforms of the methionine 129 variant of PrP to that of the valine 129 variant. We provide evidence for differences in the folding behavior between the two variants at the early stages of oligomer formation. We show that Met(129) has a higher propensity to form beta-sheet-rich oligomers, whereas Val(129) has a higher tendency to fold into alpha-helical-rich monomers. An equimolar mixture of both variants displayed an intermidate folding behavior. We show that the oligomers of both variants are initially a mixture of alpha- and beta-rich conformers that evolve with time to an increasingly homogeneous beta-rich form. This maturation process, which involves no further change in proteinase K resistance, occurs more rapidly in the Met(129) form than the Val(129) form. Although the involvement of such beta-rich oligomers in prion pathogenesis is speculative, the misfolding behavior could, in part, explain the higher susceptibility of individuals that are methionine homozygote to both sporadic and variant Creutzfeldt-Jakob disease.

Prediction of exposed domains of envelope glycoprotein in Indian HIV-1 isolates and experimental confirmation of their immunogenicity in humans

We describe the impact of subtype differences on the seroreactivity of linear antigenic epitopes in envelope glycoprotein of HIV-1 isolates from different geographical locations. By computer analysis, we predicted potential antigenic sites of envelope glycoprotein (gp120 and gp4l) of this virus. For this purpose, after fetching sequences of proteins of interest from data banks, values of hydrophilicity, flexibility, accessibility, inverted hydrophobicity, and secondary structure were considered. We identified several potential antigenic epitopes in a B subtype strain of envelope glycoprotein of HIV-1 (IIIB). Solid- phase peptide synthesis methods of Merrifield and Fmoc chemistry were used for synthesizing peptides. These synthetic peptides corresponded mainly to the C2, V3 and CD4 binding sites of gp120 and some parts of the ectodomain of gp41. The reactivity of these peptides was tested by ELISA against different HIV-1-positive sera from different locations in India. For two of these predicted epitopes, the corresponding Indian consensus sequences (LAIERYLKQQLLGWG and DIIGDIRQAHCNISEDKWNET) (subtype C) were also synthesized and their reactivity was tested by ELISA. These peptides also distinguished HIV-1-positive sera of Indians with C subtype infections from sera from HIV-negative subjects.

Origami in the outer membrane: the transmembrane arrangement of mitochondrial porins

Voltage-dependent anion-selective channels (VDAC), also known as mitochondrial porins, are key regulators of metabolite flow across the mitochondrial outer membrane. Porins from a wide variety of organisms share remarkably similar electrophysiological properties, in spite of considerable sequence dissimilarity, indicating that they share a common structure. Based on primary sequence considerations, analogy with bacterial porins, and circular dichroism analysis, it is agreed that VDAC spans the outer membrane as a beta-barrel. However, the residues that form the antiparallel beta-strands comprising this barrel remain unknown. Various predictive methods, largely based on the known structures of bacterial beta-barrels, have been applied to the primary sequences of VDAC. Refinement and confirmation of these predictions have developed through numerous investigations of wild-type and variant porins, both in mitochondria and in artificial membranes. These experiments have involved VDAC from several sources, precluding the generation of a unified model. Herein, using the Neurospora VDAC sequence as a template, the published structural information and predictions have been reassessed to delineate a model that satisfies most of the available data.

Identification of a "glycine-loop"-like coiled structure in the 34 AA Pro,Gly,Met repeat domain of the biomineral-associated protein, PM27

Fracture resistance in biomineralized structures has been linked to the presence of proteins, some of which possess sequences that are associated with elastic behavior. One such protein superfamily, the Pro,Gly-rich sea urchin intracrystalline spicule matrix proteins, form protein-protein supramolecular assemblies that modify the microstructure and fracture-resistant properties of the calcium carbonate mineral phase within embryonic sea urchin spicules and adult sea urchin spines. In this report, we detail the identification of a repetitive keratin-like "glycine-loop"- or coil-like structure within the 34-AA (AA: amino acid) N-terminal domain, (PGMG)(8)PG, of the spicule matrix protein, PM27. The identification of this repetitive structural motif was accomplished using two capped model peptides: a 9-AA sequence, GPGMGPGMG, and a 34-AA peptide representing the entire motif. Using CD, NMR spectrometry, and molecular dynamics simulated annealing/minimization simulations, we have determined that the 9-AA model peptide adopts a loop-like structure at pH 7.4. The structure of the 34-AA polypeptide resembles a coil structure consisting of repeating loop motifs that do not exhibit long-range ordering. Given that loop structures have been associated with protein elastic behavior and protein motion, it is plausible that the 34-AA Pro,Gly,Met repeat sequence motif in PM27 represents a putative elastic or mobile domain.

Role of SA-Le(a) and E-selectin in metastasis assessed with peptide antagonist

E-selectin ligand Sialyl-Lewis a (SA-Le(a)) carbohydrate is expressed on many carcinomas. Peptide mimicking SA-Le(a) (DLWDWVVGKPAG) was previously selected from a recombinant library by screening with monoclonal antibody (MAb) NS19-9. In this study, the residues critical for interaction with the NS19-9 were mapped using peptide array generated by substitution of various amino acid residues. The replacement of Trp 5 with Phe resulted in a change of peptide's secondary structure and increased binding with MAb and E-selectin, suggesting improved carbohydrate mimicry. Colonization of tumor cells expressing SA-Le(a) was blocked by the peptide and was completely abolished in E-selectin knock out mice. The data suggest the critical role of carbohydrate antigens and E-selectin in metastasis and that peptides mimicking carbohydrate antigens can function as antagonists of this process.

Prediction of beta-turns in proteins using the first-order Markov models

We present a method based on the first-order Markov models for predicting simple beta-turns and loops containing multiple turns in proteins. Sequences of 338 proteins in a database are divided using the published turn criteria into the following three regions, namely, the turn, the boundary, and the nonturn ones. A transition probability matrix is constructed for either the turn or the nonturn region using the weighted transition probabilities computed for dipeptides identified from each region. There are two such matrices constructed for the boundary region since the transition probabilities for dipeptides immediately preceding or following a turn are different. The window used for scanning a protein sequence from amino (N-) to carboxyl (C-) terminal is a hexapeptide since the transition probability computed for a turn tetrapeptide is capped at both the N- and C- termini with a boundary transition probability indexed respectively from the two boundary transition matrices. A sum of the averaged product of the transition probabilities of all the hexapeptides involving each residue is computed. This is then weighted with a probability computed from assuming that all the hexapeptides are from the nonturn region to give the final prediction quantity. Both simple beta-turns and loops containing multiple turns in a protein are then identified by the rising of the prediction quantity computed. The performance of the prediction scheme or the percentage (%) of correct prediction is evaluated through computation of Matthews correlation coefficients for each protein predicted. It is found that the prediction method is capable of giving prediction results with better correlation between the percent of correct prediction and the Matthews correlation coefficients for a group of test proteins as compared with those predicted using some secondary structural prediction methods. The prediction accuracy for about 40% of proteins in the database or 50% of proteins in the test set is better than 70%. Such a percentage for the test set is reduced to 30 if the structures of all the proteins in the set are treated as unknown.

Reconstitution of a chloroplast protein import channel

The chloroplastic outer envelope protein OEP75 with a molecular weight of 75 kDa probably forms the central pore of the protein import machinery of the outer chloroplastic membrane. Patch-clamp analysis shows that heterologously expressed, purified and reconstituted OEP75 constitutes a voltage-gated ion channel with a unit conductance of Lambda = 145pS. Activation of the OEP75 channel in vitro is completely dependent on the magnitude and direction of the voltage gradient. Therefore, movements of protein charges of parts of OEP75 in the membrane electric field are required either for pore formation or its opening. In the presence of precursor protein from only one side of the bilayer, strong flickering and partial closing of the channel was observed, indicating a specific interaction of the precursor with OEP75. The comparatively low ionic conductance of OEP75 is compatible with a rather narrow aqueous pore (dporeapproximately equal to 8-9 A). Provided that protein and ion translocation occur through the same pore, this implies that the environment of the polypeptide during the transit is mainly hydrophilic and that protein translocation requires almost complete unfolding of the precursor.

The use and misuse of FTIR spectroscopy in the determination of protein structure

Fourier transform infrared (FTIR) spectroscopy is an established tool for the structural characterization of proteins. However, many potential pitfalls exist for the unwary investigator. In this review we critically assess the application of FTIR spectroscopy to the determination of protein structure by (1) outlining the principles underlying protein secondary structure determination by FTIR spectroscopy, (2) highlighting the situations in which FTIR spectroscopy should be considered the technique of choice, (3) discussing the manner in which experiments should be conducted to derive as much physiologically relevant information as possible, and (4) outlining current methods for the determination of secondary structure from infrared spectra of proteins.

A quantum-mechanical study of the chain-length dependent stability of the extended and 3(10)-helix conformations in dehydroalanine oligopeptides

A quantum-chemical study of the chain-length dependent stability of the extended, 2-ribbon and 3(10)-helix conformations in dehydroalanine (delta Ala) oligopeptides has been performed by using both semiempirical AM1 and ab initio 4-31G methodologies. The validity of both methods in the study of the conformational properties of delta Ala oligopeptides was tested first on the dipeptide. The results of this test showed that 4-31G and AM1 calculations are in good agreement with 6-31G* calculations and experimental data. In order to monitor the conformational conversions, delta Ala oligopeptides comprising two to six residues were constructed. Molecular geometries were fully optimized using AM1, and the final conformations were verified to be minima by analysis of the corresponding second-derivative matrices. Conformational studies revealed that the 3(10)-helix is stabilized with respect to the 2(7)-ribbon when the number of residues is three or four, at the AM1 and ab initio 4-31G level respectively, while the extended form is the most stable in all the calculations performed. On the other hand, if a linear behaviour is assumed for longer chains, our calculations show a trend that would predict a conversion from extended form to 3(10)-helix in oligopeptides with around six (ab initio 4-31G) or eight (AM1) delta Ala residues. In order to explain these conformational changes, the cooperative effects for the different conformers were investigated. Large cooperative energy effects were found for the 3(10)-helix conformation.

Molecular cloning and evolutionary analysis of the calcium-modulated contractile protein, centrin, in green algae and land plants

Centrin (= caltractin) is a ubiquitous, cytoskeletal protein which is a member of the EF-hand superfamily of calcium-binding proteins. A centrin-coding cDNA was isolated and characterized from the prasinophyte green alga Scherffelia dubia. Centrin PCR amplification primers were used to isolate partial, homologous cDNA sequences from the green algae Tetraselmis striata and Spermatozopsis similis. Annealing analyses suggested that centrin is a single-copy-coding region in T. striata and S. similis and other green algae studied. Centrin-coding regions from S. dubia, S. similis and T. striata encode four colinear EF-hand domains which putatively bind calcium. Phylogenetic analyses, including homologous sequences from Chlamydomonas reinhardtii and the land plant Atriplex nummularia, demonstrate that the domains of centrins are congruent and arose from the two-fold duplication of an ancestral EF hand with Domains 1+3 and Domains 2+4 clustering. The domains of centrins are also congruent with those of calmodulins demonstrating that, like calmodulin, centrin is an ancient protein which arose within the ancestor of all eukaryotes via gene duplication. Phylogenetic relationships inferred from centrin-coding region comparisons mirror results of small subunit ribosomal RNA sequence analyses suggesting that centrin-coding regions are useful evolutionary markers within the green algae.

Cloning, sequencing and immunological characterization of the corrinoid-containing subunit of the N5-methyltetrahydromethanopterin: coenzyme-M methyltransferase from Methanobacterium thermoautotrophicum

A 3.5-kb EcoRI fragment of the Methanobacterium thermoautotrophicum chromosome contains five open reading frames, mtrA to mtrE. The deduced N-terminal amino acid sequence of mtrA is identical with 26 N-terminal amino acids of a corrinoid-containing membrane protein from Methanobacterium. Computer-aided analyses of mtrA predicts 237 amino acids with a molecular mass of 25,603 Da for its gene product. A hydropathy plot of this amino acid sequence indicates one hydrophobic helical conformation near the N-terminus of the peptide which represents a tentative membrane-spanning region. The main part of the protein, however, shows hydrophilic domains, suggesting a location outside the cytoplasmic membrane. These domains are probably accessible by monospecific polyclonal antibodies raised previously against the corrinoid-containing membrane protein. The immunogold-labeling technique revealed that the corrinoid-dependent membrane protein was detectable at the cytoplasmic face of the membranes and of vesicle preparations. No significant identity of the deduced amino acid sequence was found with sequences of several corrinoid-containing enzymes. In contrast to the hydrophilic gene product of mtrA, four other gene products from the gene cluster encode extremely hydrophobic proteins. The N-terminal sequences of mtrC and mtrD are identical with two peptides of the N5-methyltetrahydromethanopterin:coenzyme-M methyltransferase complex from Methanobacterium, indicating that the mtr genes encode this membrane protein.

The secondary structure of milk proteins and their biological function

A brief overview is given of the methods of determining and predicting secondary structure in proteins. The secondary structures of the milk serum proteins, lactoferrin, alpha-lactalbumin, lysozyme, and beta-lactoglobulin, as determined by x-ray crystallography, are compared with the results of a joint prediction method. This comparison evaluates critically the degree of success achieved and helps define what can reasonably be expected from a prediction in the absence of a known structure. The value of supplementary information from spectroscopic methods and the use of templates and sequence information from related proteins in improving the confidence of predictions are illustrated. One point that emerges is the general overprediction of helix content by the joint prediction method such that, for an all-beta protein such as beta-lactoglobulin, the method of Garnier, Osguthorpe, and Robson, applied with the correctly selected decision constants, provides a somewhat better approach. Secondary structure of the caseins can be predicted with less confidence than for globular proteins, and the results should be interpreted as evidence of a propensity to form transient secondary structures of the indicated type in view of the generally open and flexible conformation of caseins in solution.

Molecular cloning of a 47 kDa tissue-specific and differentiation-dependent urothelial cell surface glycoprotein

Despite the fact that bladder epithelium has many interesting biological features and is a frequent site of carcinoma formation, relatively little is known about its biochemical differentiation. We have shown recently that a 47 kDa glycoprotein, uroplakin III (UPIII), in conjunction with uroplakins I (27 kDa) and II (15 kDa), forms the asymmetric unit membrane (AUM)--a highly specialized biomembrane characteristic of the apical surface of bladder epithelium. Deglycosylation and cDNA sequencing revealed that UPIII contains up to 20 kDa of N-linked sugars attached to a core protein of 28.9 kDa. The presence of an N-terminal signal peptide sequence and a single transmembrane domain located near the C terminus, plus the N-terminal location of all the potential N-glycosylation sites, points to a type I (N-exo/C-cyto) configuration. Thus the mass of the extracellular domain (20 kDa plus up to 20 kDa of sugar) of UPIII greatly exceeds that of its intracellular domain (5 kDa). Such an asymmetrical mass distribution, a feature shared by the other two major uroplakins, provides a molecular explanation as to why the luminal leaflet of AUM is almost twice as thick as the cytoplasmic one. The fact that of the three major proteins of AUM only UPIII has a significant cytoplasmic domain suggests that this molecule may play an important role in AUM-cytoskeleton interaction in terminally differentiated urothelial cells.

An algorithm for protein engineering: simulations of recursive ensemble mutagenesis

An algorithm for protein engineering, termed recursive ensemble mutagenesis, has been developed to produce diverse populations of phenotypically related mutants whose members differ in amino acid sequence. This method uses a feedback mechanism to control successive rounds of combinatorial cassette mutagenesis. Starting from partially randomized "wild-type" DNA sequences, a highly parallel search of sequence space for peptides fitting an experimenter's criteria is performed. Each iteration uses information gained from the previous rounds to search the space more efficiently. Simulations of the technique indicate that, under a variety of conditions, the algorithm can rapidly produce a diverse population of proteins fitting specific criteria. In the experimental analog, genetic selection or screening applied during recursive ensemble mutagenesis should force the evolution of an ensemble of mutants to a targeted cluster of related phenotypes.

Evolution of EF-hand calcium-modulated proteins. II. Domains of several subfamilies have diverse evolutionary histories

In the first report in this series we described the relationships and evolution of 152 individual proteins of the EF-hand subfamilies. Here we add 66 additional proteins and define eight (CDC, TPNV, CLNB, LPS, DGK, 1F8, VIS, TCBP) new subfamilies and seven (CAL, SQUD, CDPK, EFH5, TPP, LAV, CRGP) new unique proteins, which we assume represent new subfamilies. The main focus of this study is the classification of individual EF-hand domains. Five subfamilies--calmodulin, troponin C, essential light chain, regulatory light chain, CDC31/caltractin--and three uniques--call, squidulin, and calcium-dependent protein kinase--are congruent in that all evolved from a common four-domain precursor. In contrast calpain and sarcoplasmic calcium-binding protein (SARC) each evolved from its own one-domain precursor. The remaining 19 subfamilies and uniques appear to have evolved by translocation and splicing of genes encoding the EF-hand domains that were precursors to the congruent eight and to calpain and to SARC. The rates of evolution of the EF-hand domains are slower following formation of the subfamilies and establishment of their functions. Subfamilies are not readily classified by patterns of calcium coordination, interdomain linker stability, and glycine and proline distribution. There are many homoplasies indicating that similar variants of the EF-hand evolved by independent pathways.

Dissection of an enzyme by protein engineering. The N and C-terminal fragments of barnase form a native-like complex with restored enzymic activity

A method is described for producing fragments of a protein suitable for studies of protein folding. The codon for a single methionine residue is introduced into the cloned gene of barnase, and the gene product cleaved with cyanogen bromide. The site of mutation was chosen to be at the surface of the protein in a region connecting segments of secondary structure in the native enzyme. The alpha + beta protein was mutated from Val36----Met, and split into two fragments, B(1-36) containing the alpha-helical regions and B(37-110), the beta-sheet. The fragments were purified by ion exchange chromatography. Neither retains catalytic activity. Fluorescence, circular dichroism, and 1H nuclear magnetic resonance data indicate that their structures are each close to that of random-coil peptides. The two fragments associate to form a tight complex (Kd = 0.2 to 0.6 microM), which displays spectroscopic properties similar to those of the uncleaved protein. The catalytic activity is restored in the complex with a value for Km similar to that for native enzyme but with kcat reduced about three- to fourfold. The second-order rate constant for association on mixing fragments in the concentration range 2.5 to 7.5 microM is 1 x 10(5) s-1 M-1.

Molecular self-assembly and nanochemistry: a chemical strategy for the synthesis of nanostructures

Molecular self-assembly is the spontaneous association of molecules under equilibrium conditions into stable, structurally well-defined aggregates joined by noncovalent bonds. Molecular self-assembly is ubiquitous in biological systems and underlies the formation of a wide variety of complex biological structures. Understanding self-assembly and the associated noncovalent interactions that connect complementary interacting molecular surfaces in biological aggregates is a central concern in structural biochemistry. Self-assembly is also emerging as a new strategy in chemical synthesis, with the potential of generating nonbiological structures with dimensions of 1 to 10(2) nanometers (with molecular weights of 10(4) to 10(10) daltons). Structures in the upper part of this range of sizes are presently inaccessible through chemical synthesis, and the ability to prepare them would open a route to structures comparable in size (and perhaps complementary in function) to those that can be prepared by microlithography and other techniques of microfabrication.

The combination of symbolic and numerical computation for three-dimensional modeling of RNA

Three-dimensional (3-D) structural models of RNA are essential for understanding of the cellular roles played by RNA. Such models have been obtained by a technique based on a constraint satisfaction algorithm that allows for the facile incorporation of secondary and other structural information. The program generates 3-D structures of RNA with atomic-level resolution that can be refined by numerical techniques such as energy minimization. The precision of this technique was evaluated by comparing predicted transfer RNA loop and RNA pseudoknot structures with known or consensus structures. The root-mean-square deviation (2.0 to 3.0 angstroms before minimization) between predicted and control structures reveal this system to be an effective method in modeling RNA.

Vaccines

Much progress has been made towards reaching an understanding of immune responses at the molecular level. This has provided much needed information for identifying the antigens which will afford protection against diseases such as rabies, malaria, whooping cough, hepatitis and acquired immune deficiency syndrome, and for presenting them to the immune system.

A class of amphipathic proteins associated with lipid storage bodies in plants. Possible similarities with animal serum apolipoproteins

The lipid-storing tissues of plants contain many small (0.2-1 microns) lipid (normally triacylglycerol) droplets which are surrounded and stabilized by a mixed phospholipid and protein annulus. The proteinaceous components of the lipid storage bodies are termed oleosins and are not associated with any other cellular structures. The major oleosins of rapeseed and radish have been isolated by preparative SDS-PAGE and are respectively classes of 19 kDa and 20 kDa proteins. Both protein classes were N-terminally blocked for direct sequencing, but were partially sequenced following limited proteolytic digestion. The major rapeseed oleosin was made up of at least two 19 kDa polypeptides, termed nap-I and nap-II, which have closely related but different amino acid sequences. A single 20 kDa oleosin, termed rad-I, was found in radish. A near full length cDNA clone for a major rapeseed oleosin was sequenced and found to correspond almost exactly to the sequence of nap-II. The sequences of nap-I and rad-I show very close similarity to one another, as do the sequences of nap-II and the previously determined sequence for the major oleosin from maize. All four oleosins have a large central hydrophobic domain flanked by polar N- and C-terminal domains. Secondary structure predictions for the four oleosins are similar and a novel model is proposed based on a central hydrophobic beta-strand region flanked by an N-terminal polar alpha-helix and a C-terminal amphipathic alpha-helix. The possibility that oleosins exhibit structural and functional similarities with some animal apolipoproteins is discussed.