Extent and nature of contacts between protein molecules in crystal lattices and between subunits of protein oligomers

A survey was compiled of several characteristics of the intersubunit contacts in 58 oligomeric proteins, and of the intermolecular contacts in the lattice for 223 protein crystal structures. The total number of atoms in contact and the secondary structure elements involved are similar in the two types of interfaces. Crystal contact patches are frequently smaller than patches involved in oligomer interfaces. Crystal contacts result from more numerous interactions by polar residues, compared with a tendency toward nonpolar amino acids at oligomer interfaces. Arginine is the only amino acid prominent in both types of interfaces. Potentials of mean force for residue–residue contacts at both crystal and oligomer interfaces were derived from comparison of the number of observed residue–residue interactions with the number expected by mass action. They show that hydrophobic interactions at oligomer interfaces favor aromatic amino acids and methionine over aliphatic amino acids; and that crystal contacts form in such a way as to avoid inclusion of hydrophobic interactions. They also suggest that complex salt bridges with certain amino acid compositions might be important in oligomer formation. For a protein that is recalcitrant to crystallization, substitution of lysine residues with arginine or glutamine is a recommended strategy. Proteins 28:494–514, 1997. © 1997 Wiley-Liss, Inc.     

Surface, subunit interfaces and interior of oligomeric proteins

The solvent-accessible surface area (As) of 23 oligomeric proteins is calculated using atomic co-ordinates from high-resolution and well-refined crystal structures. As is correlated with the protein molecular weight, and a power law predicts its value to within 5% on average. The accessible surface of the average oligomer is similar to that of monomeric proteins in its hydropathy and amino acid composition. The distribution of the 20 amino acid types between the protein surface and its interior is also the same as in monomers. Interfaces, i.e. surfaces involved in subunit contacts, differ from the rest of the subunit surface. They are enriched in hydrophobic side-chains, yet they contain a number of charged groups, especially from Arg residues, which are the most abundant residues at interfaces except for Leu. Buried Arg residues are involved in H-bonds between subunits. We counted H-bonds at interfaces and found that several have none, others have one H-bond per 200 A2 of interface area on average (1 A = 0.1 nm). A majority of interface H-bonds involve charged donor or acceptor groups, which should make their contribution to the free energy of dissociation significant, even when they are few. The smaller interfaces cover about 700 A2 of the subunit surface. The larger ones cover 3000 to 10,000 A2, up to 40% of the subunit surface area in catalase. The lower value corresponds to an estimate of the accessible surface area loss required for stabilizing subunit association through the hydrophobic effect alone. Oligomers with small interfaces have globular subunits with accessible surface areas similar to those of monomeric proteins. We suggest that these oligomers assemble from preformed monomers with little change in conformation. In oligomers with large interfaces, isolated subunits should be unstable given their excessively large accessible surface, and assembly is expected to require major structural changes.     

Role of the IXI/V motif in oligomer assembly and function of StHsp14.0, a small heat shock protein from the acidothermophilic archaeon, Sulfolobus tokodaii strain 7

Small heat shock proteins (sHsps) are one of the most ubiquitous molecular chaperones. They are grouped together based on a conserved domain, the alpha-crystallin domain. Generally, sHsps exist as oligomers of 9-40 subunits, and the oligomers undergo reversible temperature-dependent dissociation into smaller species as dimers, which interact with denaturing substrate proteins. Previous studies have shown that the C-terminal region, especially the consensus IXI/V motif, is responsible for oligomer assembly. In this study, we examined deletions or mutations in the C-terminal region on the oligomer assembly and function of StHsp14.0, an sHsp from an acidothermophilic archaeon, Sulfolobus tokodaii strain 7. Mutated StHsp14.0 with C-terminal deletion or replacement of IIe residues in the IXI/V motif to Ala, Ser, or Phe residues could not form large oligomers and lost chaperone activity. StHsp14.0WKW, whose Ile residues in the IXI/V motif are changed to Trp, existed as an oligomer like that of the wild type. However, it dissociates to small oligomers and exhibits chaperone activity at relatively lowered temperature. Replacement of two Ile residues in the motif to relatively small residues, Ala or Ser, also resulted in the change of beta-sheet rich secondary structure and decrease of hydrophobicity. Interestingly, StHsp14.0 mutant with amino acid replacements to Phe kept almost the same secondary structure and relatively high hydrophobicity despite that it could not form an oligomeric structure. The results show that hydrophobicity and size of the amino acids in the IXI/V motif in the C-terminal region are responsible not only for assembly of the oligomer but also for the maintenance of beta-sheet rich secondary structure and hydrophobicity, which are important for the function of sHsp.     

Structure-based analysis of VDAC1 protein: defining oligomer contact sites

The outer mitochondrial membrane protein, the voltage-dependent anion channel (VDAC), is increasingly implicated in the control of apoptosis. Oligomeric assembly of VDAC1 was shown to be coupled to apoptosis induction, with oligomerization increasing substantially upon apoptosis induction and inhibited by apoptosis blockers. In this study, structure- and computation-based selection of the predicated VDAC1 dimerization site, in combination with site-directed mutagenesis, cysteine replacement, and chemical cross-linking, were employed to identify contact sites between VDAC1 molecules in dimers and higher oligomers. The predicted weakly stable β-strands were experimentally found to represent the interfaces between VDAC1 monomers composing the oligomer. Replacing hydrophobic amino acids with charged residues in β-strands 1, 2, and 19 interfered with VDAC1 oligomerization. The proximity of β-strands 1, 2, and 19 within the VDAC1 dimer and the existence of other association sites involving β-strand 16 were confirmed when a cysteine was introduced at defined positions in cysteineless VDAC1 mutants, together with the use of cysteine-specific cross-linker bis(maleimido)ethane. Moreover, the results suggest that VDAC1 also exists as a dimer that upon apoptosis induction undergoes conformational changes and that its oligomerization proceeds through a series of interactions involving two distinct interfaces. Dissection of VDAC1 dimerization/oligomerization as presented here provides structural insight into the oligomeric status of cellular VDAC1 under physiological and apoptotic conditions.     

Oligomeric protein structure networks: insights into protein-protein interactions

Background  

Protein-protein association is essential for a variety of cellular processes and hence a large number of investigations are being carried out to understand the principles of protein-protein interactions. In this study, oligomeric protein structures are viewed from a network perspective to obtain new insights into protein association. Structure graphs of proteins have been constructed from a non-redundant set of protein oligomer crystal structures by considering amino acid residues as nodes and the edges are based on the strength of the non-covalent interactions between the residues. The analysis of such networks has been carried out in terms of amino acid clusters and hubs (highly connected residues) with special emphasis to protein interfaces.     

Identification of an additional interaction domain in transmembrane domains 11 and 12 that supports oligomer formation in the human serotonin transporter

Na+/Cl--dependent neurotransmitter transporters form constitutive oligomers. The topological arrangement is not known, but a leucine heptad repeat in transmembrane domain (TM) 2 and a glycophorin-like motif in TM6 have been proposed to stabilize the oligomer. To determine the topology, we generated versions of the human serotonin transporter (hSERT) that carried cyan or yellow fluorescent proteins at their amino and/or carboxyl terminus. Appropriate pairs were coexpressed to measure fluorescence resonance energy transfer (FRET). Donor photobleaching FRET microscopy was employed to deduce the following arrangement: within the monomer, the amino and carboxyl termini are in close vicinity. In addition, in the oligomer, the carboxyl termini are closer to each other than the amino termini. Hence, a separate interaction domain (i.e. distinct from TM2 and TM6) must reside in the carboxyl-terminal half of hSERT. This was confirmed by expressing the amino- and carboxyl-terminal halves of hSERT. These were retained intracellularly; they also retained the coexpressed full-length transporter by forming export-deficient oligomers and, when cotransfected in all possible combinations, supported FRET. Hence, both the carboxyl and amino termini contain elements that drive oligomerization. By employing fragments comprising two neighboring TM helices, we unequivocally identified TM11/12 as a new contact site by donor photobleaching FRET and beta-lactamase protein fragment complementation assay. TM1/2 was also found to self-associate. Thus, oligomerization of hSERT involves at least two discontinuous interfaces. The currently identified interaction sites drive homophilic interactions. This is consistent with assembly of SERT oligomers in an array-like structure containing multimers of dimers.     

Beta-strand interfaces of non-dimeric protein oligomers are characterized by scattered charged residue patterns

Protein oligomers are formed either permanently, transiently or even by default. The protein chains are associated through intermolecular interactions constituting the protein interface. The protein interfaces of 40 soluble protein oligomers of stœchiometries above two are investigated using a quantitative and qualitative methodology, which analyzes the x-ray structures of the protein oligomers and considers their interfaces as interaction networks. The protein oligomers of the dataset share the same geometry of interface, made by the association of two individual β-strands (β-interfaces), but are otherwise unrelated. The results show that the β-interfaces are made of two interdigitated interaction networks. One of them involves interactions between main chain atoms (backbone network) while the other involves interactions between side chain and backbone atoms or between only side chain atoms (side chain network). Each one has its own characteristics which can be associated to a distinct role. The secondary structure of the β-interfaces is implemented through the backbone networks which are enriched with the hydrophobic amino acids favored in intramolecular β-sheets (MCWIV). The intermolecular specificity is provided by the side chain networks via positioning different types of charged residues at the extremities (arginine) and in the middle (glutamic acid and histidine) of the interface. Such charge distribution helps discriminating between sequences of intermolecular β-strands, of intramolecular β-strands and of β-strands forming β-amyloid fibers. This might open new venues for drug designs and predictive tool developments. Moreover, the β-strands of the cholera toxin B subunit interface, when produced individually as synthetic peptides, are capable of inhibiting the assembly of the toxin into pentamers. Thus, their sequences contain the features necessary for a β-interface formation. Such β-strands could be considered as ‘assemblons’, independent associating units, by homology to the foldons (independent folding unit). Such property would be extremely valuable in term of assembly inhibitory drug development.     

Mariner Mos1 transposase dimerizes prior to ITR binding

The mariner Mos1 synaptic complex consists of a tetramer of transposase molecules that bring together the two ends of the element. Such an assembly requires at least two kinds of protein-protein interfaces. The first is involved in "cis" dimerization, and consists of transposase molecules bound side-by-side on a single DNA molecule. The second, which is involved in "trans" dimerization, consists of transposase molecules bound to two different DNA molecules. Here, we used biochemical and genetic methods to enhance the definition of the regions involved in cis and trans-dimerization in the mariner Mos1 transposase. The cis and trans-dimerization interfaces were both found within the first 143 amino acid residues of the protein. The cis-dimerization activity was mainly contained in amino acids 1-20. The region spanning from amino acid residues 116-143, and containing the WVPHEL motif, was involved in the cis- to trans-shift as well as in trans-dimerization stabilization. Although the transposase exists mainly as a monomer in solution, we provide evidence that the transposase cis-dimer is the active species in inverted terminal repeat (ITR) binding. We also observed that the catalytic domain of the mariner Mos1 transposase modulates efficient transposase-transposase interactions in the absence of the transposon ends.     

Concentration and Extinction Coefficient Determination for Oligonucleotides and Analogs Using a General Phosphate Analysis

A general phosphate analysis (GPA) is developed which assays the concentration of nucleic acid oligomers and their analogs based on stoichiometric phosphorus in the sequence. The method involves complete digestion of the oligomer sample to orthophosphate using acid at high temperature and subsequent colorimetric analysis by phosphomolybdate complex formation. GPA is applied to oligomers having phosphodiester, methylphosphonate, and phosphorothioate backbone linkages. Given the absorption spectra of oligomers having these backbones, extinction coefficients are obtained and compared to other quantitative and predictive methods. In addition to sequences having the usual nucleoside residues found in naturally occurring nucleic acids, oligomers having base analog residues can be readily quantified by GPA.     

Glycoprotein B of herpes simplex virus type 1 oligomerizes through the intermolecular interaction of a 28-amino-acid domain.

Herpes simplex virus type 1 glycoprotein B (gB) is an envelope component that plays an essential role in virus infection. The biologically active form of gB is an oligomer that contributes to the process of viral envelope fusion with the cell surface membrane, resulting in viral penetration and initiation of the replication cycle. In previous studies, two discontinuous sites for oligomer formation were identified: a nonessential upstream site located between residues 93 and 282 and an essential downstream site located between residues 596 and 711. In this study, in vitro-transcribed and -translated gB test molecules were used to characterize the more active essential membrane-proximal domain. A series of gB test polypeptides mutated in this downstream oligomerization domain were assayed for their abilities to form oligomers with a mutant gB capture polypeptide containing the analogous wild-type domain. Detection of oligomers was achieved by coimmunoprecipitation of two gB mutant molecules by using a monoclonal antibody specific for a hemagglutinin epitope tag introduced into the coding sequence of the capture polypeptide. Analysis of the immune-precipitated products by sodium dodecyl sulfate-polyacrylamide gel electrophoresis demonstrated that the downstream oligomerization domain resided within residues 626 to 676. This region was further resolved into two segments, residues 626 to 653 and 653 to 675, each of which was independently sufficient to form oligomers. However, residues 626 to 653 provided for a stronger interaction between gB monomers. Moreover, this stretch of 28 amino acids was shown to form oligomers when introduced into the carboxy-terminal region of gB monomers lacking this domain at the normal site, thus indicating that this domain was functionally independent of its natural location within the gB molecule. Further analysis of the sequence within residues 596 to 653 by using mutant test polypeptides altered in individual amino acids revealed that cysteines 9 and 10 located at positions 596 and 633, respectively, were not required for oligomer formation but contributed to dimer formation and/or stabilization. The results of this study suggest that oligomerization of gB monomers is induced by interactions between contiguous residues localized within the ectodomain near the site of molecule insertion into the viral envelope membrane.     

Analysis of Shigella flexneri Wzz (Rol) function by mutagenesis and cross-linking: Wzz is able to oligomerize

The modal length or degree of polymerization (dp) of the Shigella flexneri O-antigen is determined in an unknown manner by the Wzz/Rol protein. The Wzz protein is anchored into the cytoplasmic membrane by two transmembrane domains (TM1 amino acids 32-52; TM2 amino acids 295-315) with the central loop of the protein located in the periplasm. Plasmids were constructed encoding hybrid Wzz proteins consisting of regions of S. flexneri Wzz (WzzSF) and Salmonella typhimurium Wzz (WzzST). These imparted O-antigen modal chain lengths that implied that the carboxy-terminal region of Wzz was involved in chain length determination. Site-directed mutagenesis was undertaken to investigate the functional significance of highly conserved residues in amino-/carboxy-terminal domains of WzzSF. Some of the WzzSF variants resulted in O-antigen modal chain lengths much shorter than those of wild-type WzzSF, whereas other mutants inactivated WzzSF function entirely and a third class had a longer O-antigen chain length distribution. The data indicate that amino acids throughout the length of the WzzSF protein are important in determination of O-antigen modal chain length. In vivo cross-linking experiments were performed to investigate the interactions between Wzz proteins. The experiments indicated that the WzzSF protein is able to form dimers and oligomers of at least six WzzSF proteins. A carboxy-terminal-truncated WzzSF protein having the amino terminal 194 amino acids was able to oligomerize, indicating that the amino-terminal region is sufficient for the Wzz-Wzz interaction observed. Shortened WzzSF proteins having internal deletions in the amino-terminal region were also able to oligomerize, suggesting that residues 59-194 are not essential for oligomerization. Cross-linking of WzzSF proteins with mutationally altered residues showed that loss of WzzSF function may be correlated to a reduced/altered ability to form oligomers, and that mutational alteration of glycine residues in the TM2 segment affects WzzSF-WzzSF dimer mobility in SDS polyacrylamide gels. These results provide the first evidence of protein-protein interactions for proteins involved in O-antigen polysaccharide biosynthesis.     

Direct incorporation of microtubule oligomers at high GTP concentrations

Chick brain microtubule protein consists primarily of a mixture of MAP2:tubulin oligomers and dimeric tubulin. The assembly of this protein is described by a single pseudofirst-order reaction at 20 microM GTP, but by the summation of two pseudofirst-order reactions at 1 mM GTP. The protein contains two GTP-binding species, corresponding to the tubulin dimers and the oligomers, and conditions which alter the dimer: oligomer equilibrium, affect the kinetics of microtubule assembly. The results indicate that the oligomers are only direct assembly intermediates at high GTP concentrations.     

Study of the assembly of vesicular stomatitis virus N protein: role of the P protein

To derive structural information about the vesicular stomatitis virus (VSV) nucleocapsid (N) protein, the N protein and the VSV phosphoprotein (P protein) were expressed together in Escherichia coli. The N and P proteins formed soluble protein complexes of various molar ratios when coexpressed. The major N/P protein complex was composed of 10 molecules of the N protein, 5 molecules of the P protein, and an RNA. A soluble N protein-RNA oligomer free of the P protein was isolated from the N/P protein-RNA complex using conditions of lowered pH. The molecular weight of the N protein-RNA oligomer, 513,879, as determined by analytical ultracentrifugation, showed that it was composed of 10 molecules of the N protein and an RNA of approximately 90 nucleotides. The N protein-RNA oligomer had the appearance of a disk with outer diameter, inner diameter, and thickness of 148 +/- 10 A, 78 +/- 9 A, and 83 +/- 8 A, respectively, as determined by electron microscopy. RNA in the complexes was protected from RNase digestion and was stable at pH 11. This verified that N/P protein complexes expressed in E. coli were competent for encapsidation. In addition to coexpression with the full-length P protein, the N protein was expressed with the C-terminal 72 amino acids of the P protein. This portion of the P protein was sufficient for binding to the N protein, maintaining it in a soluble state, and for assembly of N protein-RNA oligomers. With the results provided in this report, we propose a model for the assembly of an N/P protein-RNA oligomer.     

Solid-phase synthesis of functionalized bis-peptides

We demonstrate the first solid-phase synthesis of highly functionalized bis-peptides. Bis-peptides are ladder oligomers composed of stereochemically pure, cyclic bis-amino acids joined by substituted diketopiperazine linkages. They have a shape-programmable backbone that is controlled by controlling the stereochemistry and sequence of the monomers within each oligomer. Functionalized bis-peptides are assembled using a new amide bond forming reaction (acyl-transfer coupling) that we have previously developed and a novel activation strategy that allows the sequential formation of penta- and hexa-substituted diketopiperazines from extremely hindered N-alkyl-alpha,alpha-disubstituted amino acids. We present mechanistic evidence that acyl-transfer coupling is competitive with direct acylation in the formation of hindered amide bonds. We also detail the synthesis of four functionalized bis-peptides, and that by combining bis-peptides with amino acids through diketopiperazine linkages, bis-peptides can mimic the display of residues i, i+4, i+7 of an alpha-helical peptide.     

Intra-membrane oligomerization and extra-membrane oligomerization of amyloid-β peptide are competing processes as a result of distinct patterns of motif interplay

Soluble oligomers of amyloid-β peptide (Aβ) are emerging as the primary neurotoxic species in Alzheimer disease, however, whether the membrane is among their direct targets that mediate the downstream adverse effects remains elusive. Herein, we show that multiple soluble oligomeric Aβ preparations, including Aβ-derived diffusible ligand, protofibril, and zinc-induced Aβ oligomer, exhibit much weaker capability to insert into the membrane than Aβ monomer. Aβ monomers prefer incorporating into membrane rather than oligomerizing in solution, and such preference can be reversed by the aggregation-boosting factor, zinc ion. Further analyses indicate that the membrane-embedded oligomers of Aβ are derived from rapid assembly of inserted monomers but not due to the insertion of soluble Aβ oligomers. By comparing the behavior of a panel of Aβ truncation variants, we demonstrate that the intra- and extra-membrane oligomerization are mutually exclusive processes that proceed through distinct motif interplay, both of which require the action of amino acids 37-40/42 to overcome the auto-inhibitory interaction between amino acids 29-36 and the N-terminal portion albeit via different mechanisms. These results indicate that intra- and extra-membrane oligomerization of Aβ are competing processes and emphasize a critical regulation of membrane on the behavior of Aβ monomer and soluble oligomers, which may determine distinct neurotoxic mechanisms.     

Intermolecular β-Strand Networks Avoid Hub Residues and Favor Low Interconnectedness: A Potential Protection Mechanism against Chain Dissociation upon Mutation

Altogether few protein oligomers undergo a conformational transition to a state that impairs their function and leads to diseases. But when it happens, the consequences are not harmless and the so-called conformational diseases pose serious public health problems. Notorious examples are the Alzheimer''s disease and some cancers associated with a conformational change of the amyloid precursor protein (APP) and of the p53 tumor suppressor, respectively. The transition is linked with the propensity of β-strands to aggregate into amyloid fibers. Nevertheless, a huge number of protein oligomers associate chains via β-strand interactions (intermolecular β-strand interface) without ever evolving into fibers. We analyzed the layout of 1048 intermolecular β-strand interfaces looking for features that could provide the β-strands resistance to conformational transitions. The interfaces were reconstructed as networks with the residues as the nodes and the interactions between residues as the links. The networks followed an exponential decay degree distribution, implying an absence of hubs and nodes with few links. Such layout provides robustness to changes. Few links per nodes do not restrict the choices of amino acids capable of making an interface and maintain high sequence plasticity. Few links reduce the “bonding” cost of making an interface. Finally, few links moderate the vulnerability to amino acid mutation because it entails limited communication between the nodes. This confines the effects of a mutation to few residues instead of propagating them to many residues via hubs. We propose that intermolecular β-strand interfaces are organized in networks that tolerate amino acid mutation to avoid chain dissociation, the first step towards fiber formation. This is tested by looking at the intermolecular β-strand network of the p53 tetramer.     

Conformational differences between two amyloid β oligomers of similar size and dissimilar toxicity

Several protein conformational disorders (Parkinson and prion diseases) are linked to aberrant folding of proteins into prefibrillar oligomers and amyloid fibrils. Although prefibrillar oligomers are more toxic than their fibrillar counterparts, it is difficult to decouple the origin of their dissimilar toxicity because oligomers and fibrils differ both in terms of structure and size. Here we report the characterization of two oligomers of the 42-residue amyloid β (Aβ42) peptide associated with Alzheimer disease that possess similar size and dissimilar toxicity. We find that Aβ42 spontaneously forms prefibrillar oligomers at Aβ concentrations below 30 μm in the absence of agitation, whereas higher Aβ concentrations lead to rapid formation of fibrils. Interestingly, Aβ prefibrillar oligomers do not convert into fibrils under quiescent assembly conditions but instead convert into a second type of oligomer with size and morphology similar to those of Aβ prefibrillar oligomers. Strikingly, this alternative Aβ oligomer is non-toxic to mammalian cells relative to Aβ monomer. We find that two hydrophobic peptide segments within Aβ (residues 16-22 and 30-42) are more solvent-exposed in the more toxic Aβ oligomer. The less toxic oligomer is devoid of β-sheet structure, insoluble, and non-immunoreactive with oligomer- and fibril-specific antibodies. Moreover, the less toxic oligomer is incapable of disrupting lipid bilayers, in contrast to its more toxic oligomeric counterpart. Our results suggest that the ability of non-fibrillar Aβ oligomers to interact with and disrupt cellular membranes is linked to the degree of solvent exposure of their central and C-terminal hydrophobic peptide segments.     

Evolvability of yeast protein-protein interaction interfaces

The functional importance of protein-protein interactions indicates that there should be strong evolutionary constraint on their interaction interfaces. However, binding interfaces are frequently affected by amino acid replacements. Change due to coevolution within interfaces can contribute to variability but is not ubiquitous. An alternative explanation for the ability of surfaces to accept replacements may be that many residues can be changed without affecting the interaction. Candidates for these types of residues are those that make interchain interaction only through the protein main chain, β-carbon, or associated hydrogen atoms. Since almost all residues have these atoms, we hypothesize that this subset of interface residues may be more easily substituted than those that make interactions through other atoms. We term such interactions "residue type independent." Investigating this hypothesis, we find that nearly a quarter of residues in protein interaction interfaces make exclusively interchain residue-type-independent contacts. These residues are less structurally constrained and less conserved than residues making residue-type-specific interactions. We propose that residue-type-independent interactions allow substitutions in binding interfaces while the specificity of binding is maintained.     

Prediction of protein homo-oligomer types by pseudo amino acid composition: Approached with an improved feature extraction and Naive Bayes Feature Fusion

Summary. The interaction of non-covalently bound monomeric protein subunits forms oligomers. The oligomeric proteins are superior to the monomers within the scope of functional evolution of biomacromolecules. Such complexes are involved in various biological processes, and play an important role. It is highly desirable to predict oligomer types automatically from their sequence. Here, based on the concept of pseudo amino acid composition, an improved feature extraction method of weighted auto-correlation function of amino acid residue index and Naive Bayes multi-feature fusion algorithm is proposed and applied to predict protein homo-oligomer types. We used the support vector machine (SVM) as base classifiers, in order to obtain better results. For example, the total accuracies of A, B, C, D and E sets based on this improved feature extraction method are 77.63, 77.16, 76.46, 76.70 and 75.06% respectively in the jackknife test, which are 6.39, 5.92, 5.22, 5.46 and 3.82% higher than that of G set based on conventional amino acid composition method with the same SVM. Comparing with Chou’s feature extraction method of incorporating quasi-sequence-order effect, our method can increase the total accuracy at a level of 3.51 to 1.01%. The total accuracy improves from 79.66 to 80.83% by using the Naive Bayes Feature Fusion algorithm. These results show: 1) The improved feature extraction method is effective and feasible, and the feature vectors based on this method may contain more protein quaternary structure information and appear to capture essential information about the composition and hydrophobicity of residues in the surface patches that buried in the interfaces of associated subunits; 2) Naive Bayes Feature Fusion algorithm and SVM can be referred as a powerful computational tool for predicting protein homo-oligomer types.