Structure of Petunia hybrida Defensin 1, a Novel Plant Defensin with Five Disulfide Bonds

The structure of a novel plant defensin isolated from the flowers of Petunia hybrida has been determined by H-1 NMR spectroscopy. P. hybrida defensin 1 (PhD1) is a basic, cysteine-rich, antifungal protein of 47 residues and is the first example of a new subclass of plant defensins with five disulfide bonds whose structure has been determined. PhD1 has the fold of the cysteine-stabilized alphabeta motif, consisting of an alpha-helix and a triple-stranded antiparallel beta-sheet, except that it contains a fifth disulfide bond from the first loop to the alpha-helix. The additional disulfide bond is accommodated in PhD1 without any alteration of its tertiary structure with respect to other plant defensins. Comparison of its structure with those of classic, four-disulfide defensins has allowed us to identify a previously unrecognized hydrogen bond network that is integral to structure stabilization in the family.

Protein contact prediction using patterns of correlation

We describe a new method for using neural networks to predict residue contact pairs in a protein. The main inputs to the neural network are a set of 25 measures of correlated mutation between all pairs of residues in two windows of size 5 centered on the residues of interest. While the individual pair-wise correlations are a relatively weak predictor of contact, by training the network on windows of correlation the accuracy of prediction is significantly improved. The neural network is trained on a set of 100 proteins and then tested on a disjoint set of 1033 proteins of known structure. An average predictive accuracy of 21.7% is obtained taking the best L/2 predictions for each protein, where L is the sequence length. Taking the best L/10 predictions gives an average accuracy of 30.7%. The predictor is also tested on a set of 59 proteins from the CASP5 experiment. The accuracy is found to be relatively consistent across different sequence lengths, but to vary widely according to the secondary structure. Predictive accuracy is also found to improve by using multiple sequence alignments containing many sequences to calculate the correlations. (C) 2004 Wiley-Liss, Inc.

Cardiac expression of the cystic fibrosis transmembrane conductance regulator involves novel Exon 1 usage to produce a unique amino-terminal protein

Cystic fibrosis is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which encodes a chloride channel present in many cells. In cardiomyocytes, we report that multiple exon 1 usage and alternative splicing produces four CFTR transcripts, with different 5'-untranslated regions, CFTRTRAD-139, CFTR-1C/-1A, CFTR-1C, and CFTR-1B. CFTR transcripts containing the novel upstream exons (exons -1C, -1B, and -1A) represent more than 90% of cardiac expressed CFTR mRNA. Regulation of cardiac CFTR expression, in response to developmental and pathological stimuli, is exclusively due to the modulation of CFTR-1C and CFTR-1C/-1A expression. Upstream open reading frames have been identified in the 5'-untranslated regions of all CFTR transcripts that, in conjunction with adjacent stem-loop structures, modulate the efficiency of translation initiation at the AUG codon of the main CFTR coding region in CFTRTRAD-139 and CFTR-1C/-1A transcripts. Exon(-1A), only present in CFTR-1C/-1A transcripts, encodes an AUG codon that is in-frame with the main CFTR open reading frame, the efficient translation of which produces a novel CFTR protein isoform with a curtailed amino terminus. As the expression of this CFTR transcript parallels the spatial and temporal distribution of the cAMP-activated whole-cell current density in normal and diseased hearts, we suggest that CFTR-1C/-1A provides the molecular basis for the cardiac cAMP-activated chloride channel. Our findings provide further insight into the complex nature of in vivo CFTR expression, to which multiple mRNA transcripts, protein isoforms, and post-transcriptional regulatory mechanisms are now added.

Better prediction of protein contact number using a support vector regression analysis of amino acid sequence

Background: Protein tertiary structure can be partly characterized via each amino acid's contact number measuring how residues are spatially arranged. The contact number of a residue in a folded protein is a measure of its exposure to the local environment, and is defined as the number of C-beta atoms in other residues within a sphere around the C-beta atom of the residue of interest. Contact number is partly conserved between protein folds and thus is useful for protein fold and structure prediction. In turn, each residue's contact number can be partially predicted from primary amino acid sequence, assisting tertiary fold analysis from sequence data. In this study, we provide a more accurate contact number prediction method from protein primary sequence. Results: We predict contact number from protein sequence using a novel support vector regression algorithm. Using protein local sequences with multiple sequence alignments (PSI-BLAST profiles), we demonstrate a correlation coefficient between predicted and observed contact numbers of 0.70, which outperforms previously achieved accuracies. Including additional information about sequence weight and amino acid composition further improves prediction accuracies significantly with the correlation coefficient reaching 0.73. If residues are classified as being either contacted or non-contacted, the prediction accuracies are all greater than 77%, regardless of the choice of classification thresholds. Conclusion: The successful application of support vector regression to the prediction of protein contact number reported here, together with previous applications of this approach to the prediction of protein accessible surface area and B-factor profile, suggests that a support vector regression approach may be very useful for determining the structure-function relation between primary sequence and higher order consecutive protein structural and functional properties.

Chemical synthesis and structure elucidation of bovine kappa-casein (1-44)

The caseins (alpha(s1), alpha(s2), beta, and kappa) are phosphoproteins present in bovine milk that have been studied for over a century and whose structures remain obscure. Here we describe the chemical synthesis and structure elucidation of the N-terminal segment (1-44) of bovine K-casein, the protein which maintains the micellar structure of the caseins. K-Casein (1-44) was synthesised by highly optimised Boc solid-phase peptide chemistry and characterised by mass spectrometry. Structure elucidation was carried out by circular dichroism and nuclear magnetic resonance spectroscopy. CD analysis demonstrated that the segment was ill defined in aqueous medium but in 30% trifluoroethanol it exhibited considerable helical structure. Further, NMR analysis showed the presence of a helical segment containing 26 residues which extends from Pro(8) to Arg(34). This is the first report which demonstrates extensive secondary structure within the casein class of proteins. (c) 2006 Elsevier Inc. All rights reserved.

Predicting the solvent accessibility of transmembrane residues from protein sequence

In this study, we propose a novel method to predict the solvent accessible surface areas of transmembrane residues. For both transmembrane alpha-helix and beta-barrel residues, the correlation coefficients between the predicted and observed accessible surface areas are around 0.65. On the basis of predicted accessible surface areas, residues exposed to the lipid environment or buried inside a protein can be identified by using certain cutoff thresholds. We have extensively examined our approach based on different definitions of accessible surface areas and a variety of sets of control parameters. Given that experimentally determining the structures of membrane proteins is very difficult and membrane proteins are actually abundant in nature, our approach is useful for theoretically modeling membrane protein tertiary structures, particularly for modeling the assembly of transmembrane domains. This approach can be used to annotate the membrane proteins in proteomes to provide extra structural and functional information.

The 2.1 angstrom crystal structure of copGFP, a representative member of the copepod clade within the green fluorescent protein superfamily

The green fluorescent protein (avGFP), its variants, and the closely related GFP-like proteins are characterized structurally by a cyclic tri-peptide chromophore located centrally within a conserved beta-can fold. Traditionally, these GFP family members have been isolated from the Cnidaria although recently, distantly related GFP-like proteins from the Bilateria, a sister group of the Cnidaria have been described, although no representative structure from this phylum has been reported to date. We have determined to 2.1 angstrom resolution the crystal structure of copGFP, a representative GFP-like protein from a copepod, a member of the Bilateria. The structure of copGFP revealed that, despite sharing only 19% sequence identity with GFP, the tri-peptide chromophore (Gly57-Tyr58-Gly59) of copGFP adopted a cis coplanar conformation within the conserved beta-can fold. However, the immediate environment surrounding the chromophore of copGFP was markedly atypical when compared to other members of the GFP-superfamily, with a large network of bulky residues observed to surround the chromophore. Arg87 and Glu222 (GFP numbering 96 and 222), the only two residues conserved between copGFP, GFP and GFP-like proteins are involved in autocatalytic genesis of the chromophore. Accordingly, the copGFP structure provides an alternative platform for the development of a new suite of fluorescent protein tools. Moreover, the structure suggests that the autocatalytic genesis of the chromophore is remarkably tolerant to a high degree of sequence and structural variation within the beta-can fold of the GFP superfamily. (c) 2006 Elsevier Ltd . All rights reserved.

Direct protein interaction underlies gene-for-gene specificity and coevolution of the flax resistance genes and flax rust avirulence genes

Plant resistance proteins (R proteins) recognize corresponding pathogen avirulence (Avr) proteins either indirectly through detection of changes in their host protein targets or through direct R-Avr protein interaction. Although indirect recognition imposes selection against Avr effector function, pathogen effector molecules recognized through direct interaction may overcome resistance through sequence diversification rather than loss of function. Here we show that the flax rust fungus AvrLS67 genes, whose products are recognized by the L5, L6, and L7 R proteins of flax, are highly diverse, with 12 sequence variants identified from six rust strains. Seven AvrL567 variants derived from Avr alleles induce necrotic responses when expressed in flax plants containing corresponding resistance genes (R genes), whereas five variants from avr alleles do not. Differences in recognition specificity between AvA567 variants and evidence for diversifying selection acting on these genes suggest they have been involved in a gene-specific arms race with the corresponding flax R genes. Yeast two-hybrid assays indicate that recognition is based on direct R-Avr protein interaction and recapitulate the interaction specificity observed in planta. Biochemical analysis of Escherichia coli-produced AvrL567 proteins shows that variants that escape recognition nevertheless maintain a conserved structure and stability, suggesting that the amino acid sequence differences directly affect the R-Avr protein interaction. We suggest that direct recognition associated with high genetic diversity at corresponding R and Avr gene loci represents an alternative outcome of plant-pathogen coevolution to indirect recognition associated with simple balanced polymorphisms for functional and nonfunctional R and Avr genes.

Cytotoxicity of pilosulin 1, a peptide from the venom of the jumper ant Myrmecia pilosula

The synthetic peptide pilosulin 1, corresponding to the largest defined allergenic polypeptide found in the venom of the jumper ant Myrmecia pilosula, inhibited the incorporation of [methyl-H-3]thymidine into proliferating Epstein-Barr transformed (EBV) B-cells. The LD50 was four-fold lower in concentration than melittin, a cytotoxic peptide found in honey bee venom. Loss of cell viability was assessed by flow cytometry by measuring the proportion of cells that fluoresced in the presence of the fluorescent dye 7-aminoactinomycin D. Examination of proliferating EBV B-cells indicated that the cells lost viability within a few minutes exposure to pilosulin 1. Partial peptides of pilosulin 1 were less efficient in causing loss of cell viability and the results suggest that the 22 N-terminal residues are critical to the cytotoxic activity of pilosulin 1. Normal blood white cells were also labile to pilosulin I. T- and B-lymphocytes, monocytes and natural killer cells, however, were more labile than granulocytes. Analysis of pilosulin I using circular dichroism indicated that, in common with melittin and other Hymenoptera venom toxins, it had the potential to adopt an cc-helical secondary structure. (C) 1998 Elsevier Science B.V, All rights reserved.

The dimerization and topological specificity functions of MinE reside in a structurally autonomous C-terminal domain

Correct placement of the division septum in Escherichia coli requires the co-ordinated action of three proteins, MinC, MinD and MinE. MinC and MinD interact to form a non-specific division inhibitor that blocks septation at all potential division sites. MinE is able to antagonize MinCD in a topologically sensitive manner, as it restricts MinCD activity to the unwanted division sites at the cell poles, Here, we show that the topological specificity function of MinE residues in a structurally autonomous, trypsin-resistant domain comprising residues 31-88, Nuclear magnetic resonance (NMR) and circular dichroic spectroscopy indicate that this domain includes both alpha and beta secondary structure, while analytical ultracentrifugation reveals that it also contains a region responsible for MinE homodimerization. While trypsin digestion indicates that the anti-MinCD domain of MinE (residues 1-22) does not form a tightly folded structural domain, NMR analysis of a peptide corresponding to MinE(1-22) indicates that this region forms a nascent helix in which the peptide rapidly interconverts between disordered (random coil) and alpha-helical conformations, This suggests that the N-terminal region of MinE may be poised to adopt an alpha-helical conformation when it interacts with the target of its anti-MinCD activity, presumably MinD.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Sausage: protein threading with flexible force fields

Sausage is a protein sequence threading program, but with remarkable run-time flexibility. Using different scripts, it can calculate protein sequence-structure alignments, search structure libraries, swap force fields, create models form alignments, convert file formats and analyse results. There are several different force fields which might be classed as knowledge-based, although they do not rely on Boltzmann statistics. Different force fields are used for alignment calculations and subsequent ranking of calculated models.

Solution structure of χ-conopeptide MrIA, a modulator of the human norepinephrine transporter

The chi-conopeptides MrIA and MrIB are 13-residue peptides with two disulfide bonds that inhibit human and rat norepinephrine transporter systems and are of significant interest for the design of novel drugs involved in pain treatment. In the current study we have determined the solution structure of MrIA using NMR spectroscopy. The major element of secondary structure is a hairpin with the two strands connected by an inverse gamma-turn. The residues primarily involved in activity have previously been shown to be located in the turn region (Sharpe, I. A.; Palant, E.: Schroder, C. L; Kaye, D. M.; Adams, D. I.; Alewood, P. F.; Lewis, R. J. J Biol Client 2003, 278, 40317-40323), which appears to be more flexible than the beta-strands based on disorder in the ensemble of calculated structures. Analogues of MrIA with N-terminal truncations indicate that the N-terminal residues play a role in defining a stable conformation and the native disulfide connectivity. In particular, noncovalent interactions between Val3 and Hypl2 are likely to be involved in maintaining a stable conformation. The N-terminus also affects activity, as a single N-terminal deletion introduced additional pharmacology at rat vas deferens, while deleting the first two amino acids reduced chi-conopeptide potency. This article was originally published online as an accepted preprint. The Published Online date corresponds to the preprint version. You can request a copy of the preprint by entailing the Biopolymers editorial office at biopolymers@wiley.com (c) 2005 Wiley Periodicals, Inc.