The Stabilizing Effects of O-glycosylation on the Secondary Structural Integrity of the Designed α-loop-α motif by Molecular Dynamics Simulations

Abstract

In this study, various 400 ps molecular dynamics simulations were conducted to determine the stabilizing effect of O-glycosylation on the secondary structural integrity of the design α-loop-α motif, which has the optimal loop length of 7 Gly residues (denoted as N-A₁₆G₇A₁₆-C). In general, O-glycosylation stabilizes the structural integrity of the model peptide regardless of the length and position of glycosylation sites because it decreases the opportunity for water molecules to compete for the intramolecular hydrogen bonds. The designed peptide exhibits the highest helicity when residues 11 and 31 are replaced with Ser residues followed by O-linked with 3 galactose residues, representing the “face-to-face” glycosylation near the loop. In this case, the loop exhibits an extended conformation and several new hydrogen bonds are observed between the main chain of the loop and the galactose residues, resulting in decreasing the fluctuation and increasing the stability of the entire peptide. When the glycosylation are made close to the loop, the secondary structural integrity of the α-loop-α motif increases with the number of galactose residues. In addition, “face- to-face” glycosylation increases the structural integrity of this motif to a greater extent than “back-to-back” glycosylation. However, when the glycosylation are created away from the loop and near the N- and C-termini, no general rule is found for the stabilizing effect.     

The identification of conserved interactions within the SH3 domain by alignment of sequences and structures

The SH3 domain, comprised of approximately 60 residues, is found within a wide variety of proteins, and is a mediator of protein-protein interactions. Due to the large number of SH3 domain sequences and structures in the databases, this domain provides one of the best available systems for the examination of sequence and structural conservation within a protein family. In this study, a large and diverse alignment of SH3 domain sequences was constructed, and the pattern of conservation within this alignment was compared to conserved structural features, as deduced from analysis of eighteen different SH3 domain structures. Seventeen SH3 domain structures solved in the presence of bound peptide were also examined to identify positions that are consistently most important in mediating the peptide-binding function of this domain. Although residues at the two most conserved positions in the alignment are directly involved in peptide binding, residues at most other conserved positions play structural roles, such as stabilizing turns or comprising the hydrophobic core. Surprisingly, several highly conserved side-chain to main-chain hydrogen bonds were observed in the functionally crucial RT-Src loop between residues with little direct involvement in peptide binding. These hydrogen bonds may be important for maintaining this region in the precise conformation necessary for specific peptide recognition. In addition, a previously unrecognized yet highly conserved beta-bulge was identified in the second beta-strand of the domain, which appears to provide a necessary kink in this strand, allowing it to hydrogen bond to both sheets comprising the fold.     

Binding patterns of DTR-specific antibodies reveal a glycosylation-conditioned tumor-specific epitope of the epithelial mucin (MUC1)

Glycosylation determines essential biological functions of epithelial mucins in health and disease. We report on the influence of glycosylation of the immunodominant DTR motif of MUC1 on its antigenicity. Sets of novel glycopeptides were synthesized that enabled us to examine sole and combined effects of peptide length (number of repeats) and O-glycosylation with GalNAc at the DTR motif on the binding patterns of 22 monoclonal antibodies recognizing this motif. In case of unglycosylated peptides almost all antibodies bound better to multiple MUC1 tandem repeats. Glycosylation at the DTR led to enhanced binding in 11 cases, whereas 10 antibodies were not influenced in binding, and one was inhibited. In nine of the former cases both length and DTR glycosylation were additive in their influence on antibody binding, suggesting that both effects are different. Improved binding to the glycosylated DTR motif was exclusively found with antibodies generated against tumor-derived MUC1. Based on these data a tumor-specific MUC1 epitope is defined comprising the ...PDTRP... sequence in a particular conformation essentially determined by O-glycosylation at its threonine with either GalNAcalpha1 or a related short glycan. The results can find application in the field of MUC1-based immunotherapy.     

Molecular dynamics simulations to determine the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme

In this study, various molecular dynamics simulations were conducted to investigate the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme. The analyses of backbone root-mean-square deviation, radius of gyration, and secondary structure stability all show that supercritical CO(2) exhibits the ability to increase the stability of this protein, probably as a result of the solvent with less polarity, where hydrophobic interactions stabilizing the native structure are weakened and simultaneously the local hydrogen bonds are strengthened, resulting in stabilization of the secondary structures. The hydrophobic cores in the alpha- and beta-domains also play an important role in preventing this protein from thermal unfolding. As supercritical CO(2) has been attractive for biomedical applications because of the advantages of mild critical condition, nonflammability, nontoxity, and the purity of the resulting products, the structural stabilizing effect found in this study strongly suggests that it is possible to increase the thermostability of hen egg white lysozyme by pretreatment with supercritical CO(2), leading to better industrial applications of this protein.     

Analysis of post-translational modification of mycobacterial proteins using a cassette expression system

A recombinant expression system was developed to analyse sequence determinants involved in O-glycosylation of proteins in mycobacteria. By expressing peptide sequences corresponding to known glycosylation sites within a chimeric lipoprotein construct, amino acids flanking modified threonine residues were found to have an important influence on glycosylation. The expression system was used to screen mycobacterial sequences selected using a neural network (NetOglyc) trained on eukaryotic O-glycoproteins. Evidence of glycosylation was obtained for eight of 11 proteins tested. The results suggest that sites involved in O-glycosylation of mycobacterial and eukaryotic proteins share similar structural features.     

Molecular dynamics simulations to determine the optimal loop length in the helix-loop-helix motif

In this study, 100 ps molecular dynamics simulations were conducted to determine the optimal length of the Gly loop in the helix-loop-helix motif. The results revealed that the length of the loop significantly affected the stability of the two alpha-helices. In addition, the optimal loop length to maintain the highest helicity of the target peptide was found to be the one corresponding to 7 Gly residues. Longer loop length resulted in structural destabilization due to the entropic effect, while shorter loop length was not suitable to achieve intrachain packing of the two alpha-helices.     

Role of peptide sequence and neighboring residue glycosylation on the substrate specificity of the uridine 5'-diphosphate-alpha-N-acetylgalactosamine:polypeptide N-acetylgalactosaminyl transferases T1 and T2: kinetic modeling of the porcine and canine submaxillary gland mucin tandem repeats

A large family of uridine 5'-diphosphate (UDP)-alpha-N-acetylgalactosamine (GalNAc):polypeptide N-acetylgalactosaminyl transferases (ppGalNAc Ts) initiates mucin-type O-glycan biosynthesis at serine and threonine. The peptide substrate specificities of individual family members are not well characterized or understood, leaving an inability to rationally predict or comprehend sites of O-glycosylation. Recently, a kinetic modeling approach demonstrated neighboring residue glycosylation as a major factor modulating the O-glycosylation of the porcine submaxillary gland mucin 81 residue tandem repeat by ppGalNAc T1 and T2 [Gerken et al. (2002) J. Biol. Chem. 277, 49850-49862]. To confirm the general applicability of this model and its parameters, the ppGalNAc T1 and T2 glycosylation kinetics of the 80+ residue tandem repeat from the canine submaxillary gland mucin was obtained and characterized. To reproduce the glycosylation patterns of both mucins (comprising 50+ serine/threonine residues), specific effects of neighboring peptide sequence, in addition to the previously described effects of neighboring residue glycosylation, were required of the model. Differences in specificity of the two transferases were defined by their sensitivities to neighboring proline and nonglycosylated hydroxyamino acid residues, from which a ppGalNAc T2 motif was identified. Importantly, the model can approximate the previously reported ppGalNAc T2 glycosylation kinetics of the IgA1 hinge domain peptide [Iwasaki, et al. (2003) J. Biol. Chem. 278, 5613-5621], further validating both the approach and the ppGalNAc T2 positional weighting parameters. The characterization of ppGalNAc transferase specificity by this approach may prove useful for the search for isoform-specific substrates, the creation of isoform-specific inhibitors, and the prediction of mucin-type O-glycosylation sites.     

Solution structure of the PDZ2 domain from human phosphatase hPTP1E and its interactions with C-terminal peptides from the Fas receptor

The solution structure of the second PDZ domain (PDZ2) from human phosphatase hPTP1E has been determined using 2D and 3D heteronuclear NMR experiments. The binding of peptides derived from the C-terminus of the Fas receptor to PDZ2 was studied via changes in backbone peptide and protein resonances. The structure is based on a total of 1387 nonredundant experimental NMR restraints including 1261 interproton distance restraints, 45 backbone hydrogen bonds, and 81 torsion angle restraints. Analysis of 30 lowest-energy structures resulted in rmsd values of 0.41 +/- 0.09 A for backbone atoms (N, Calpha, C') and 1.08 +/- 0.10 A for all heavy atoms, excluding the disordered N- and C-termini. The hPTP1E PDZ2 structure is similar to known PDZ domain structures but contains two unique structural features. In the peptide binding domain, the first glycine of the GLGF motif is replaced by a serine. This serine appears to replace a bound water observed in PDZ crystal structures that hydrogen bonds to the bound peptide's C-terminus. The hPTP1E PDZ2 structure also contains an unusually large loop following strand beta2 and proximal to the peptide binding site. This well-ordered loop folds back against the PDZ domain and contains several residues that undergo large amide chemical shift changes upon peptide binding. Direct observation of peptide resonances demonstrates that as many as six Fas peptide residues interact with the PDZ2 domain.     

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.     

Frequent occurrence of the T-loop RNA folding motif in ribosomal RNAs

Analysis of atomic resolution structures of the rRNAs within the context of the 50S and the 30S ribosomal subunits have revealed the presence of nine examples of a recurrent structural motif, first observed in the TpsiC loop of tRNAs. The key component of this T-loop motif is a UA trans Watson-Crick/Hoogsteen base pair stacked on a Watson-Crick pair on one side. This motif is stabilized by several noncanonical hydrogen bonds, facilitating RNA-RNA as well as RNA-protein interactions. In particular, the sugar edge of the purine on the 3' side of the pivotal uridine in the UA pair frequently forms a noncanonical base pair with a distant residue. The bulged-out bases, usually seen as part of the motif, also use their Watson-Crick edges to interact with nearby residues via base-specific hydrogen bonds. In certain occurrences, a backbone reversal is stabilized by specific hydrogen bonds as is observed in the U-turn motifs and the adenosine residue of the key UA pair interacts with a third base via its Watson-Crick edge, essentially generating a base triple.     

'Wave-type' structure of a synthetic hexaglycosylated decapeptide: A part of the extracellular domain of human glycophorin A

The three-dimensional structure of a glycopeptide, His-Thr*-Ser*-Thr*-Ser*-Ser*-Ser*-Val-Thr-Lys, with 2-acetamido-2-deoxy--D-galactose (GalNAc) residues linked to six adjacent amino acids from Thr-10 to Ser-15, was studied by NMR spectroscopy and molecular dynamics (MD) simulations. The hexaglycosylated decapeptide is part of the extracellular domain of human glycophorin A and shows an extended structure of the peptide backbone due to O-glycosylation. Furthermore, each GalNAc residue exhibits one and only one NOE contact from the NHAc proton to the backbone amide proton of the amino acid that the sugar is directly bound to. This indicates a strong preference for the orientation of all GalNAc residues towards the N-terminus. NOE build-up curves were used to determine 42 inter-proton distances that, in connection with angles of the peptide backbone obtained from 3J-coupling constants, resulted in constraints for a MD simulation in water. The NMR data and the MD simulations show a preference for an extended backbone structure. The GalNAc residues are located alternatingly on opposite sides of the backbone and reduce the flexibility of the peptide backbone. The conformation of the molecule is relatively rigid and shows a 'wave-type' 3D structure of the peptide backbone within the glycosylation cluster. This new structural element is also supported by the unusual CD spectrum of the glycopeptide.     

Disruption of the coenzyme binding site and dimer interface revealed in the crystal structure of mitochondrial aldehyde dehydrogenase "Asian" variant

Mitochondrial aldehyde dehydrogenase (ALDH2) is the major enzyme that oxidizes ethanol-derived acetaldehyde. A nearly inactive form of the enzyme, ALDH2*2, is found in about 40% of the East Asian population. This variant enzyme is defined by a glutamate to lysine substitution at residue 487 located within the oligomerization domain. ALDH2*2 has an increased Km for its coenzyme, NAD+, and a decreased kcat, which lead to low activity in vivo. Here we report the 2.1 A crystal structure of ALDH2*2. The structure shows a large disordered region located at the dimer interface that includes much of the coenzyme binding cleft and a loop of residues that form the base of the active site. As a consequence of these structural changes, the variant enzyme exhibits rigid body rotations of its catalytic and coenzyme-binding domains relative to the oligomerization domain. These structural perturbations are the direct result of the inability of lysine 487 to form important stabilizing hydrogen bonds with arginines 264 and 475. Thus, the elevated Km for coenzyme exhibited by this variant probably reflects the energetic penalty for reestablishing this site for productive coenzyme binding, whereas the structural alterations near the active site are consistent with the lowered Vmax.     

Conformational studies on the MUC1 tandem repeat glycopeptides: implication for the enzymatic O-glycosylation of the mucin protein core

The tandem repeat of the MUC1 protein core is a major site of O-glycosylation that is catalyzed by several polypeptide GalNAc-transferases. To define structural features of the peptide substrates that contribute to acceptor substrate efficiency, solution structures of the 21-residue peptide AHGVTSAPDTRPAPGSTAPPA (AHG21) from the MUC1 protein core and four isoforms, glycosylated with alpha-N-acetylgalactosamine on corresponding Thr residues, AHG21 (T5), AHG21 (T10), AHG21 (T17), and AHG21 (T5,T17), were investigated by NMR spectroscopy and computational methods. NMR studies revealed that sugar attachment affected the conformational equilibrium of the peptide backbone near the glycosylated Thr residues. The clustering of the low-energy conformations for nonglycosylated and glycosylated counterparts within the VTSA, DTR, and GSTA fragments (including all sites of potential glycosylation catalyzed by GalNAc-T1, -T2, and -T4 transferases) showed that the glycosylated peptides display distinct structural propensities that may explain, in part, the differences in substrate specificities exhibited by these polypeptide GalNAc-transferases.     

Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

N for AsN – O for StrOcture? A strand–loop–strand motif for prokaryotic O‐glycosylation

So far, it has not been possible to identify a general sequence motif for O-glycosylation in bacteria. In this issue, Charbonneau et al. (2012) demonstrate why O-glycosylation is mediated by a 13-residue strand-loop-strand motif which is part of a 19-residue imperfect repeat in the passenger domain of bacterial autotransporters. This motif provides a convenient 'glycosylation cassette' and raises intriguing questions about the structural regulation of the glycosylation pathway.     

Energetics and cooperativity of the hydrogen bonding and anchor interactions that bind peptides to MHC class II protein

The complexity of the interaction between major histocompatibility complex class II (MHC II) proteins and peptide ligands has been revealed through structural studies and crystallographic characterization. Peptides bind through side-chain "anchor" interactions with MHC II pockets and an extensive array of genetically conserved hydrogen bonds to the peptide backbone. Here we quantitatively investigate the kinetic hierarchy of these interactions. We present results detailing the impact of single side-chain mutations of peptide anchor residues on dissociation rates, utilizing two I-A(d)-restricted peptides, one of which has a known crystal structure, and 24 natural and non-natural amino acid mutant variants of these peptides. We find that the N-terminal P1, P4 and P6 anchor-pocket interactions can make significant contributions to binding stability. We also investigate the interactions of these peptides with four I-A(d) MHC II proteins, each mutated to disrupt conserved hydrogen bonds to the peptide backbone. These complexes exhibit kinetic behavior suggesting that binding energy is disproportionately invested near the peptide N terminus for backbone hydrogen bonds. We then evaluate the effects of simultaneously modifying both anchor and hydrogen bonding interactions. A quantitative analysis of 71 double mutant cycles reveals that there is little apparent cooperativity between anchor residue interactions and hydrogen bonds, even when they are directly adjacent (<5A).     

Mucin core O-glycosylation is modulated by neighboring residue glycosylation status. Kinetic modeling of the site-specific glycosylation of the apo-porcine submaxillary mucin tandem repeat by UDP-GalNAc:polypeptide N-acetylgalactosaminyltransferases T1 and T2

The influence of peptide sequence and environment on the initiation and elongation of mucin O-glycosylation is not well understood. The in vivo glycosylation pattern of the porcine submaxillary gland mucin (PSM) tandem repeat containing 31 O-glycosylation sites (Gerken, T. A., Gilmore, M., and Zhang, J. (2002) J. Biol. Chem. 277, 7736-7751) reveals a weak inverse correlation with hydroxyamino acid density (and by inference the density of glycosylation) with the extent of GalNAc glycosylation and core-1 substitution. We now report the time course of the in vitro glycosylation of the apoPSM tandem repeat by recombinant UDP-GalNAc:polypeptide alpha-GalNAc transferases (ppGalNAc transferase) T1 and T2 that confirm these findings. A wide range of glycosylation rates are found, with several residues showing apparent plateaus in glycosylation. An adjustable kinetic model that reduces the first-order rate constants proportional to neighboring glycosylation status, plus or minus three residues of the site of glycosylation, was found to reasonably reproduce the experimental rate data for both transferases, including apparent plateaus in glycosylation. The unique, transferase-specific, positional weighting constants reveal information on the peptide/glycopeptide recognition site for each transferase. Both transferases displayed high sensitivities to neighboring Ser/Thr glycosylation, whereas ppGalNAc T2 displayed additional high sensitivities to the presence of nonglycosylated Ser/Thr residues. This is the first demonstration of the ability to model mucin O-glycosylation kinetics, confirming that under the appropriate conditions neighboring glycosylation status can be a significant factor modulating the first step of mucin O-glycan biosynthesis.     

Plant aquaporins with non-aqua functions: deciphering the signature sequences

Research in recent years on plant Major Intrinsic Proteins (MIPs), commonly referred to as 'aquaporins', has seen a vast expansion in the substrates found to be transported via these membrane channels. The diversity in sizes, chemical nature and physiological significance of these substrates has meant a need to critically analyse the possible structural and biochemical properties of MIPs that transport these, in order to understand their roles. In this work we have undertaken a comprehensive analysis of all plant MIPs, coming from different families, that have been proven to transport ammonia, boron, carbon dioxide, hydrogen peroxide, silicon and urea. The sequences were analysed for all primary selectivity-related motifs (NPA motifs, ar/R filter, P1-P5 residues). In addition, the putative regulatory phosphorylation and glycosylation sites and mechanistic regulators such as loop lengths have been analysed. Further, nine specificity-determining positions (SDPs) were predicted for each group. The results show the ar/R filter residues, P2-P4 positions and some of the SDPs are characteristic for certain groups, and O-glycosylation sites are unique to a subgroup while N-glycosylation was characteristic of the other MIPs. Certain residues, especially in loop C, and structural parameters such as loop lengths also contribute to the uniqueness of groups. The comprehensive analysis makes significant inroads into appraising the intriguing diversity of plant MIPs and their roles in fundamental life processes, and provides tools for plant selections, protein engineering and transgenics.     

Primary and secondary structure of the pore-forming peptide of pathogenic Entamoeba histolytica.

A pore-forming peptide is implicated in the potent cytolytic activity of pathogenic Entamoeba histolytica. Using NH2-terminal sequence information of this peptide, the corresponding cDNA was isolated. The cDNA-deduced amino acid sequence revealed a putative signal peptide and a mature peptide of 77 amino acids including six cysteine residues. Computer-aided secondary structure analysis predicted that the peptide would be composed of four adjacent alpha-helices, and CD spectroscopy indicated an all alpha-helical conformation. The tertiary structure appears to be stabilized by three disulfide bonds; the pore-forming activity was not sensitive to heat but was lost in the presence of reducing agents. Sequence homology was found to the saposins and to surfactant-associated protein B, both mammalian polypeptides of similar size and secondary structure but of non-lytic function. In particular, the six cysteine residues were found to be conserved, suggesting a common motif for stabilizing a favourable tertiary structure. Compared with previously characterized toxic peptides also containing three disulfide bonds, the amoeba peptide may represent a distinct class of biologically active peptides.