The structural and functional diversity of naturally occurring antimicrobial peptides

Antimicrobial peptides occur in a diverse range of organisms from microorganisms to insects, plants and animals. Although they all have the common function of inhibiting or killing invading microorganisms they achieve this function using an extremely diverse range of structural motifs. Their sizes range from approximately 10-90 amino acids. Most carry an overall positive charge, reflecting a preferred mode of electrostatic interaction with negatively charged microbial membranes. This article describes the structural diversity of a representative set of antimicrobial peptides divided into five structural classes: those with agr-helical structure, those with bgr-sheet structure, those with mixed helical / bgr- sheet structure, those with irregular structure, and those incorporating a macrocyclic structure. There is a significant diversity in both the size and charge of molecules within each of these classes and between the classes. The common feature of their three-dimensional structures is, however, that they have a degree of amphipathic character in which there is separate localisation of hydrophobic regions and positively charged regions. An emerging trend amongst antimicrobial proteins is the discovery of more macrocyclic analogues. Cyclisation appears to impart an additional degree of stability on these molecules and minimizes proteolytic cleavage. In conclusion, there appear to be a number of promising opportunities for the development of novel clinically useful antimicrobial peptides based on knowledge of the structures of naturally occurring antimicrobial molecules.

Common E protein determinants for attenuation of glycosaminoglycan-binding variants of Japanese encephalitis and West Nile viruses

Natural isolates and laboratory strains of West Nile virus (WNV) and Japanese encephalitis virus (JEV) were attenuated for neuroinvasiveness in mouse models for flavivirus encephalitis by serial passage in human adenocarcinoma (SW13) cells. The passage variants displayed a small-plaque phenotype, augmented affinity for heparin-Sepharose, and a marked increase in specific infectivity for SW13 cells relative to the respective parental viruses, while the specific infectivity for Vero cells was not altered. Therefore, host cell adaptation of passage variants was most likely a consequence of altered receptor usage for virus attachment-entry with the involvement of cell surface glycosaminoglycans (GAG) in this process. In vivo blood clearance kinetics of the passage variants was markedly faster and viremia was reduced relative to the parental viruses, suggesting that affinity for GAG (ubiquitously present on cell surfaces and extracellular matrices) is a key determinant for the neuroinvasiveness of encephalitic flaviviruses. A difference in pathogenesis between WNV and JEV, which was reflected in more efficient growth in the spleen and liver of the WNV parent and passage variants, accounted for a less pronounced loss of neuroinvasiveness of GAG binding variants of WNV than JEV. Single gain-of-net-positive-charge amino acid changes at E protein residue 49, 138, 306, or 389/390, putatively positioned in two clusters on the virion surface, define molecular determinants for GAG binding and concomitant virulence attenuation that are shared by the JEV serotype flaviviruses.

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.

Mimicking the action of GroEL in molecular dynamics simulations: Application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Biosynthesis and insecticidal properties of plant cyclotides: The cyclic knotted proteins from Oldenlandia affinis

Several members of the Rubiaceae and Violaceae families produce a series of cycloticles or macrocyclic peptides of 29-31 amino acids with an embedded cystine knot. We aim to understand the mechanism of synthesis of cyclic peptides in plants and have isolated a cDNA clone that encodes the cyclotide kalata Ell as well as three other clones for related cycloticles from the African plant Olden-landia affinis. The cDNA clones encode prepropeptides with a 20-aa signal sequence, an N-terminal prosequence of 46-68 amino acids and one, two, or three cyclotide domains separated by regions of about 25 aa. The corresponding cycloticles have been isolated from plant material, indicating that the cyclotide domains are excised and cyclized from all four predicted precursor proteins. The exact processing site is likely to lie on the N-terminal side of the strongly conserved GlyLeuPro or SerLeuPro sequence that flanks both sides of the cyclotide domain. Cyclotides have previously been assigned an antimicrobial function; here we describe a potent inhibitory effect on the growth and development of larvae from the Lepidopteran species Helicoverpa punctigera.

Integrating JERS-1 imaging radar and elevation models for mapping tropical rainforest communities in far North Queensland, Australia

The Wet Tropics World Heritage Area in Far North Queens- land, Australia consists predominantly of tropical rainforest and wet sclerophyll forest in areas of variable relief. Previous maps of vegetation communities in the area were produced by a labor-intensive combination of field survey and air-photo interpretation. Thus,. the aim of this work was to develop a new vegetation mapping method based on imaging radar that incorporates topographical corrections, which could be repeated frequently, and which would reduce the need for detailed field assessments and associated costs. The method employed G topographic correction and mapping procedure that was developed to enable vegetation structural classes to be mapped from satellite imaging radar. Eight JERS-1 scenes covering the Wet Tropics area for 1996 were acquired from NASDA under the auspices of the Global Rainforest Mapping Project. JERS scenes were geometrically corrected for topographic distortion using an 80 m DEM and a combination of polynomial warping and radar viewing geometry modeling. An image mosaic was created to cover the Wet Tropics region, and a new technique for image smoothing was applied to the JERS texture bonds and DEM before a Maximum Likelihood classification was applied to identify major land-cover and vegetation communities. Despite these efforts, dominant vegetation community classes could only be classified to low levels of accuracy (57.5 percent) which were partly explained by the significantly larger pixel size of the DEM in comparison to the JERS image (12.5 m). In addition, the spatial and floristic detail contained in the classes of the original validation maps were much finer than the JERS classification product was able to distinguish. In comparison to field and aerial photo-based approaches for mapping the vegetation of the Wet Tropics, appropriately corrected SAR data provides a more regional scale, all-weather mapping technique for broader vegetation classes. Further work is required to establish an appropriate combination of imaging radar with elevation data and other environmental surrogates to accurately map vegetation communities across the entire Wet Tropics.

Cycloviolacin H4, a hydrophobic cyclotide from Viola hederaceae

Cycloviolacin H4, a new macrocyclic miniprotein comprising 30 amino acid residues, was isolated from the underground parts of the Australian native violet Viola hederaceae. Its sequence, cyclo-(CAESCVWIPCTVTALLGCSCSNNVCYNGIP), was determined by nanospray tandem mass spectrometry and quantitative amino acid analysis. A knotted disuffide arrangement, which was designated as a cyclic cystine knot motif and characteristic to all known cyclotides, is proposed for stabilizing the molecular structure and folding. The cyclotide is classified in the bracelet subfamily of cyclotides due to the absence of a cis-Pro peptide bond in the circular peptide backbone. A model of its three-dimensional structure was derived based on the template of the homologous cyclotide vhr1 (Trabi et al. Plant Cell 2004, 16, 2204-2216). Cycloviolacin H4 exhibits the most potent hemolytic activity in cyclotides reported so far, and this activity correlates with the size of a surface-exposed hydrophobic patch. This work has thus provided insight into the factors that modulate the cytotoxic properties of cyclotides.

Predicting residue-wise contact orders in proteins by support vector regression

Background: The residue-wise contact order (RWCO) describes the sequence separations between the residues of interest and its contacting residues in a protein sequence. It is a new kind of one-dimensional protein structure that represents the extent of long-range contacts and is considered as a generalization of contact order. Together with secondary structure, accessible surface area, the B factor, and contact number, RWCO provides comprehensive and indispensable important information to reconstructing the protein three-dimensional structure from a set of one-dimensional structural properties. Accurately predicting RWCO values could have many important applications in protein three-dimensional structure prediction and protein folding rate prediction, and give deep insights into protein sequence-structure relationships. Results: We developed a novel approach to predict residue-wise contact order values in proteins based on support vector regression (SVR), starting from primary amino acid sequences. We explored seven different sequence encoding schemes to examine their effects on the prediction performance, including local sequence in the form of PSI-BLAST profiles, local sequence plus amino acid composition, local sequence plus molecular weight, local sequence plus secondary structure predicted by PSIPRED, local sequence plus molecular weight and amino acid composition, local sequence plus molecular weight and predicted secondary structure, and local sequence plus molecular weight, amino acid composition and predicted secondary structure. When using local sequences with multiple sequence alignments in the form of PSI-BLAST profiles, we could predict the RWCO distribution with a Pearson correlation coefficient (CC) between the predicted and observed RWCO values of 0.55, and root mean square error (RMSE) of 0.82, based on a well-defined dataset with 680 protein sequences. Moreover, by incorporating global features such as molecular weight and amino acid composition we could further improve the prediction performance with the CC to 0.57 and an RMSE of 0.79. In addition, combining the predicted secondary structure by PSIPRED was found to significantly improve the prediction performance and could yield the best prediction accuracy with a CC of 0.60 and RMSE of 0.78, which provided at least comparable performance compared with the other existing methods. Conclusion: The SVR method shows a prediction performance competitive with or at least comparable to the previously developed linear regression-based methods for predicting RWCO values. In contrast to support vector classification (SVC), SVR is very good at estimating the raw value profiles of the samples. The successful application of the SVR approach in this study reinforces the fact that support vector regression is a powerful tool in extracting the protein sequence-structure relationship and in estimating the protein structural profiles from amino acid sequences.

Derivatives of 1,3,5-triamino-1,3,5-trideoxy-cis-inositas versatile pentadentate ligands for protein labeling with Re-186/188. Prelabeling, biodistribution, and X-ray structural studies

The pentadentate H(3)bhci [1,3,5-trideoxy-1,3-bis((2-hydroxybenzyl)amino)-cis-inistol] and its bifunctionalized analogue H(3)bhci-glu-H [1,3,5-trideoxy-1,3-bis((2-hydroxybenzyl)amino)-5-glutaramido-cis-inositol] were synthesized, and their coordination chemistry was investigated with inactive rhenium, with no carrier added Re-188 and with carrier added Re-186. The neutral Re(V) complexes [ReO-(bhci)] and [ReO(bhci-glu-H)] are formed in good yields starting from [ReOCl3(P(C6H5)(3))(2)] or in quantitative yield directly from [(ReO4)-Re-186/188](-) in aqueous solution by reduction with Sn(II) or Sn(0). The X-ray structures of [ReO(bhci)] and [ReO(bhci-glu-H)] were elucidated revealing pentadentate side on coordination of the ligands to the Re=O core. The basic cyclohexane frame adopts a chair form in the case of [ReO(bhci)] and a twisted boat form in the case of [ReO(bhci-glu-H)]. [ReO(bhci)] crystallizes in the monoclinic space group C2/c with a = 27.425(3), b = 14.185(1), c = 19.047(2) Angstrom, and beta = 103.64(2)degrees and [ReO(bhci-glu-H)] in the monoclinic space group P2(1)/c with a = 13.056(3), b = 10.180(1), c = 22.378(5) Angstrom and beta = 98.205(9)degrees Both Re-188 complexes are stable in human serum for at least 3 days without decomposition. After injection into mice, [ReO(bhci-glu)](-) is readily excreted through the intestines, while [ReO(bhci)] is excreted by intestines, liver, and the kidneys. TLC investigations of the urine showed exclusively the complexes [ReO(bhci-glu-H)] and [ReO(bhci)], respectively, and no decomposition products. For derivatization of antibodies, the carboxylic group of [ReO(bhci-glu-H)] was activated with N-hydroxysuccinimide, which required unusually vigorous reaction conditions (heating). The anti colon cancer antibody mAb-35 [IgG and F(ab')(2) fragment] was labeled with [(ReO)-Re-186/188(bhci-glu)] to a specific activity of up to 1.5 mCi/mg (55 MBq/mg) with full retention of immunoreactivity. Labeling yields followed pseudo-first-order kinetics in antibody concentration with the ratio of rates between aminolysis and hydrolysis being about 2. Biodistributions of Re-186-labeled intact mAb-35 as well as of its F(ab')(2) fragment in tumor-bearing nude mice revealed good uptake by the tumor with only low accumulation of radioactivity in normal tissue.

Crystal structure of human T cell leukemia virus type 1 gp21 ectodomain crystallized as a maltose-binding protein chimera reveals structural evolution of retroviral transmembrane proteins

Retroviral entry into cells depends on envelope glycoproteins, whereby receptor binding to the surface-exposed subunit triggers membrane fusion by the transmembrane protein (TM) subunit. We determined the crystal structure at 2.5-Angstrom resolution of the ectodomain of gp21, the TM from human T cell leukemia virus type 1. The gp21 fragment was crystallized as a maltose-binding protein chimera, and the maltose-binding protein domain was used to solve the initial phases by the method of molecular replacement. The structure of gp21 comprises an N-terminal trimeric coiled coil, an adjacent disulfide-bonded loop that stabilizes a chain reversal, and a C-terminal sequence structurally distinct from HIV type 1/simian immunodeficiency virus gp41 that packs against the coil in an extended antiparallel fashion. Comparison of the gp21 structure with the structures of other retroviral TMs contrasts the conserved nature of the coiled coil-forming region and adjacent disulfide-bonded loop with the variable nature of the C-terminal ectodomain segment. The structure points to these features having evolved to enable the dual roles of retroviral TMs: conserved fusion function and an ability to anchor diverse surface-exposed subunit structures to the virion envelope and infected cell surface. The structure of gp21 implies that the N-terminal fusion peptide is in close proximity to the C-terminal transmembrane domain and likely represents a postfusion conformation.

Structural basis of recognition of monopartite and bipartite nuclear localization sequences by mammalian importin-alpha

Importin-alpha is the nuclear import receptor that recognizes cargo proteins which contain classical monopartite and bipartite nuclear localization sequences (NLSs), and facilitates their transport into the nucleus. To determine the structural basis of the recognition of the two classes of NLSs by mammalian importin-alpha, we co-crystallized an N-terminally truncated mouse receptor protein with peptides corresponding to the monopartite NLS from the simian virus 40 (SV40) large T-antigen, and the bipartite NLS from nucleoplasmin. We show that the monopartite SV40 large T-antigen NLS binds to two binding sites on the receptor, similar to what was observed in yeast importin-alpha. The nucleoplasmin NLS-importin-alpha complex shows, for the first time, the mode of binding of bipartite NLSs to the receptor. The two basic clusters in the NLS occupy the two binding sites used by the monopartite NLS, while the sequence linking the two basic clusters is poorly ordered, consistent with its tolerance to mutations. The structures explain the structural basis for binding of diverse NLSs to the sole receptor protein. (C) 2000 Academic Press.

Structural interpretation of mutations in phenylalanine hydroxylase protein aids in identifying genotype-phenotype correlations in phenylketonuria

Phenylalanine hydroxylase (PAH) is the enzyme that converts phenylalanine to tyrosine as a rate-limiting step in phenylalanine catabolism and protein and neurotransmitter biosynthesis. Over 300 mutations have been identified in the gene encoding PAH that result in a deficient enzyme activity and lead to the disorders hyperphenylalaninaemia and phenylketonuria. The determination of the crystal structure of PAH now allows the determination of the structural basis of mutations resulting in PAH deficiency. We present an analysis of the structural basis of 120 mutations with a 'classified' biochemical phenotype and/or available in vitro expression data. We find that the mutations can be grouped into five structural categories, based on the distinct expected structural and functional effects of the mutations in each category. Missense mutations and small amino acid deletions are found in three categories:'active site mutations', 'dimer interface mutations', and 'domain structure mutations'. Nonsense mutations and splicing mutations form the category of 'proteins with truncations and large deletions'. The final category, 'fusion proteins', is caused by frameshift mutations. We show that the structural information helps formulate some rules that will help predict the likely effects of unclassified and newly discovered mutations: proteins with truncations and large deletions, fusion proteins and active site mutations generally cause severe phenotypes; domain structure mutations and dimer interface mutations spread over a range of phenotypes, but domain structure mutations in the catalytic domain are more likely to be severe than domain structure mutations in the regulatory domain or dimer interface mutations.

Giant protein kinases: Domain interactions and structural basis of autoregulation

The myosin-associated giant protein kinases twitchin and titin are composed predominantly of fibronectin- and immunoglobulin-like modules, We report the crystal structures of two autoinhibited twitchin kinase fragments, one from Aplysia and a larger fragment from Caenorhabditis elegans containing an additional C-terminal immunoglobulin-like domain, The structure of the longer fragment shoes that the immunoglobulin domain contacts the protein kinase domain on the opposite side from the catalytic cleft, laterally exposing potential myosin binding residues, Together, the structures reveal the cooperative interactions between the autoregulatory region and the residues from the catalytic domain involved in protein substrate binding, ATP binding, catalysis and the activation loop, and explain the differences between the observed autoinhibitory mechanism and the one found in the structure of calmodulin-dependent kinase I.

Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR) -containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor. using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general. cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, interior/exterior orientations of side chains. and backbone conformation of the LRR structures can be predicted correctly. On the other hand. the analysis revealed that, in some cases. it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.