{'fraction_non_alphanumeric': 0.037456920007255576, 'fraction_numerical': 0.02808513211197775, 'mean_word_length': 5.003448275862069, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Structure and function of seed storage proteins in faba bean (Vicia faba L.).

The protein subunit is the most important basic unit of protein, and its study can unravel the structure and function of seed storage proteins in faba bean. In this study, we identified six specific protein subunits in Faba bean (cv. Qinghai 13) combining liquid chromatography (LC), liquid chromatography-electronic spray ionization mass (LC-ESI-MS/MS) and bio-information technology. The results suggested a diversity of seed storage proteins in faba bean, and a total of 16 proteins (four GroEL molecular chaperones and 12 plant-specific proteins) were identified from 97-, 96-, 64-, 47-, 42-, and 38-kD-specific protein subunits in faba bean based on the peptide sequence. We also analyzed the composition and abundance of the amino acids, the physicochemical characteristics, secondary structure, three-dimensional structure, transmembrane domain, and possible subcellular localization of these identified proteins in faba bean seed, and finally predicted function and structure. The three-dimensional structures were generated based on homologous modeling, and the protein function was analyzed based on the annotation from the non-redundant protein database (NR database, NCBI) and function analysis of optimal modeling. The objective of this study was to identify the seed storage proteins in faba bean and confirm the structure and function of these proteins. Our results can be useful for the study of protein nutrition and achieve breeding goals for optimal protein quality in faba bean.


Introduction
The seeds of flowering plants usually accumulate large quantities of seed protein including storage protein and lectins. The amount of protein found in seed varies from 10% (in cereals) to 40% (in certain legumes and oilseeds) of the seed dry weight, which forms a major source of dietary protein (Shimada et al. 2003;Shewry et al. 1995;Shewry and Halford 2002). Faba bean is consumed worldwide as a rich source of plant protein for humans and animals. In addition to its high protein content, faba bean is also a globally important nitrogen-fixing legume (Sosulski and McCurdy 1987;Cazzato and Tufarelli 2012;Webb et al. 2016).
Seed storage proteins account for almost 80% of the total seed protein in weight (Shewry and Halford 2002). In faba bean, seed storage proteins mainly include globulin and albumin (Liu et al. 2012b;Shi et al. 2012Shi et al. , 2013. Globulins are widely distribution and can be divided into two groups based on their sedimentation coefficient (S 20.w ): the 7S vicilin-type globulins and the 11S legumin-type globulins. The 7-11S globulins tend to be deficient in the sulfur-containing amino acids cysteine and methionine and also tryptophan. Albumin, a functional protein composed of trypsin inhibitor and phytolectin, contains more sulfurcontaining amino acids than globulin (Utsumi et al. 1980;Mandal and Mandal 2000).
The amino acid composition of seed storage proteins is important for protein quality, and thus impacts the various end uses (Shewry et al. 1995). However, seed proteins in general are deficient in some essential amino acids which compromises nutritional quality. For example, legumes are deficient in sulfur-containing amino acids (cysteine and methionine), and cereal proteins are usually deficient in lysine and tryptophan which may limit the nutrition quality of seed (Shewry et al. 1995;Mandal and Mandal 2000). Therefore, exploring the amino acid composition is also important in evaluating the nutritional quality of seed storage proteins in faba bean.
Seed storage proteins have been widely studied in soybean, rice, Arabidopsis thaliana, wheat, and pea (Scott et al. 1992;Meinke et al. 1981;Choi et al. 2000;Shimada et al. 2003;Wu et al. 2016;Gatehouse et al. 1983). Scientists also evaluated the protein nutrition and agricultural importance of seed storage proteins in faba bean Chandra-Hioe et al. 2016;Makri et al. 2005;Carbonaro et al. 2000). Bailey and Boulter (1970) analyzed the structure of legumin in faba bean seed storage proteins using peptide mapping techniques after a tryptic digest. Valizadeh (2001) analyzed the seed storage protein profile of grain legumes grown in Iran using sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). However, the seed storage protein profile of faba bean has not been explored comprehensively, and the structure of the proteins is also unknown. Our previous work only analyzed the protein composition and protein subunits for different genotypes of faba bean (Liu et al. 2012a, b;Hou et al. 2015). Given the importance of the nutritional evaluation and protein consumption of seed storage proteins in faba bean, exploring the variation of the proteins in faba bean would help breeders in their selection efforts to develop new varieties. Our study may supply a reliable data for breeding high-quality faba bean varieties. Meanwhile, the identification of proteins would provide a window into complex cellular regulatory networks (Domon and Aebersold 2006).
Biological mass spectrometry (MS) possesses an advantage when exploring protein structure and function because of its ability to precisely predict the relative molecular weight, amino acid sequence, and Post-transcriptional modification in a polypeptide or protein (Smith et al. 1997). In addition, a complex protein mixture can be distinguished by combining liquid chromatography (LC) and tandem mass spectrometry (MS-MS) (Chen 2009;Link et al. 1999). Therefore, we used LC, LC-ESI-MS/ MS, and bio-information technology to analyze the seed storage proteins of faba bean in this study. We specifically evaluated as follows: (1) the seed storage protein composition and their abundance in faba bean; (2) their physicochemical characteristics and amino acid composition; (3) and the structure and function prediction of these seed storage proteins.

Materials
The Faba bean (Vicia faba L.), cv. Qinghai 13 cultivated by Qinghai Academy of Agricultural and Forestry Science was used in this study. Endoproteinase trypsin was obtained from Promega (Madison, WI). All other chemicals used in proteolytic digestion and high-performance liquid chromatography (HPLC) were obtained from Sigma. The following instruments were used in this study: highspeed centrifuge (Thermo, Palo Alto, CA), high-performance liquid chromatograph (Agilent), and Quadrupole ion trap mass spectrometer (Thermo). The non-redundant protein database was from NCBI's GenBank.

Preparation and separation of protein
Seed storage proteins of faba bean were extracted according to Kumamaru et al. (1988) with some modification. Mature dry seeds were ground to a fine powder after strippng coat and screened through a mesh size of 80. To isolate the seed storage proteins, 50 mg of finely ground seed powder was transferred to a 1.5 mL microcentrifuge tube containing 35 lL b-mercaptoethanol and infiltrated for 20 min. Next, 675 lL protein extract buffer (8 M urea, 4% SDS, 20% glycerinum, 50 mM Tris-HCl, pH = 6.8) was added which was preheated at 40°C. The tube was vigorously vortexed for 30 s until the contents were transparent. After setting for 48 h at room temperature, the contents were centrifuged at 80009g for 5 min at 4°C. The clear supernatant was resolved by SDS-PAGE.
The seed storage proteins were first separated by SDS-PAGE and then stained with coomassie brilliant Blue. Then the band was excised from the gel for protein identification with MS/MS. The sample treatment followed a commonly used protocol (Shevchenko et al. 1996). In brief, the sample gel was cut from a target region of the gel and each protein gel band was destained with 50% MeOH in water containing 2% acetic acid for overnight. Reduction of disulfide linkages was carried out with dithiothreitol (DTT) followed by alkylation with iodoacetamide for 15 min at room temperature in the dark. The sample was then washed with water, and digested in-gel with sequencing grade-modified trypsin in the digestion buffer (ammonium bicarbonate 100 mM, pH 8.5). The peptides from the digestion were extracted with acetonitrile, and completely dried down in a SpeedVac device (Thermo). The dried sample was then redissolved in sample solution (2% acetonitrile, 97.5% water, 0.5% formic acid). A dissolved peptide sample was then analyzed by a NanoLC-ESI-MS/MS system.

NanoLC-ESI-MS/MS analysis
NanoLC-ESI-MS/MS analysis of digested protein samples was performed by an HPLC system with a 75 lm inner diameter and 8 cm in length reverse-phase C18 column. The particle size of the C18 was 3 lM with a pore size of 300 Ä . The injection time was 20 min. The HPLC Solvent A was 97.5% water, 2% acetonitrile, and 0.5% formic acid. The HPLC Solvent B was 9.5% water, 90% acetonitrile, and 0.5% formic acid. The gradation time was 60 min from 2% Solvent B to 90% solvent B, plus 20 min for sample loading and 20 min for column washing. The column flow rate was around 800 nL per minute after splitting. The typical injection volume was 3 lL.
The HPLC system was coupled with an ion trap mass spectrometer (LCQ DECA XP PLUS, Thermo) so that a sample eluted from the HPLC column was directly ionized by an electrospray ionization (ESI) process and entered into the mass spectrometer. The ionization voltage was optimized each time and was normally in a range of 1.2-1.8 kv. The capillary temperature was set at 110°C. The mass spectrometer was set at the data-dependent mode to acquire MS/MS data via a low-energy collision-induced dissociation (CID) process. The default collision energy was 33%, and the default charge state was 3. One full scan with one microscan with a mass range of 550 amu to 1800 amu was acquired, followed by one MS/MS scan of the most intense ion with a full mass range and three microscans. The dynamic exclusion feature was set as follows: repeat count of 1 within 0.3 min and exclusion duration of 0.4 min. The exclusion width was 4 Da.

Database search and validation
As a default, the mass spectrometric data were used to search against the most recent non-redundant protein database (NR database, NCBI) with ProtTech's ProtQuest software suite. The output from the database search was manually validated by a senior scientist before reporting. The parameters of the data acquisition and data processing were as following: Data acquisition parameters: Instrument: LCQ DECA XP Plus Software: Xcalibur 2.0 Centroid mode Minimum MS signal for Precursor-ion: 5 9 10 4 counts Processing parameters: Software: ProtQuest 2.0 Signal-to-noise: C5 De-isotoped: Yes Filtered: Yes Missed cleavages: B5 Modifications: Cys ?57 Database: NR database from NCBI's GenBank (update on 9/9/2011) Validation parameters: Peptide = 1, score C13 Peptide [ 1, score C15

Label-free protein quantitation
Relative protein abundance was calculated based on the method published by Griffin et al. (2010). The scoring function was based on MS abundance recorded in both MS and MS/MS datasets, spectral count (SC, number of MS/ MS spectra per peptide), unique peptide number (PN), and fragment ion (MS/MS) intensities. The normalized spectral index (SI N ) was calculated for each gel slice based on the following equation. The final relative protein abundance is the percentage of each SI N in the total SI N in each sample.
where pn is the unique peptide number, sc is the spectral count, i is the fragment ion intensity of peptide k, n is the: total number of proteins identified, and L is the number of amino acids in a protein.

Mass spectrometry analysis of different protein subunits
The protein composition, peptide mass, and peptide sequence for subunits of seed storage proteins in faba bean were obtained using LC-ESI-MS/MS. Six protein subunits of 97, 96, 64, 47, 42, and 38 kD were identified in faba bean. In total, 16 proteins were identified in these six protein subunits, and the molecule weight ranged from 20 to 57 kD (Table S1). All peptide sequences were listed and each peptide was sequenced by LC-MS independently, so even a single peptide could lead to a highly confident protein identification. Each protein was hydrolyzed into 1-10 peptide fragments through enzymatic hydrolysis at the arginine and lysine sites, and the calculated relative abundance of peptides ranged from 5.0 to 100%. In total, 42 peptide fragments were identified in the 16 seed storage proteins of faba bean (Table S1). Each of the peptide fragments was composed of 8-23 amino acids (Table S1).

Basic information of identified seed storage proteins
For further analysis, amino acid sequence alignment was performed by searching the NR database of NCBI (http:// www.ncbi.nlm.nih.gov/protein) based on identified peptide sequences. Table S2 lists 16 identified proteins from each sample. Since our results are based on MS/MS peptide sequencing of each individual peptide, all proteins listed here are actual proteins present in the sample instead of ''potential candidates'', which is often seen in a MALDI-TOF-based peptide mapping. For any identified protein here, there should be [98% certainty if the identification is based on LC-MS/MS sequence of one peptide, and it will have [99.9% certainty if it is based on sequencing of two or more peptides.
In total, 16 proteins, four GroEL molecular chaperones, and 12 plant-specific proteins (phosphate ABC transporter periplasmic substrate-binding protein, Tu elongation factor Tu, electron transfer flavoprotein subunit alpha, alkyl hydroperoxide reductase C22 subunit, cirtrate synthase, motA/TolQ/ExbB proton channel family protein, convicilin, legumin A precursor, putative sucrose-binding protein, LEGB7_VICFA, and Htlv-1 Gp21 Ectodomain MALTOSE-Binding Protein Chimera) were identified (Table S2). The 97-kD subunit was recognized as GroEL molecular chaperone, and its relative abundance reached up to 100%, which was much higher than in the 64-, 42-, and 38-kD subunits (Table S2). Among those 12 identified plant-specific proteins, convicilin in the 64-kD subunit, legumin (a precursor in the 47-and 42-kD subunits), putative sucrose-binding protein in the 47-kD subunit, and LEGB7_VICFC in the 38-kD subunit were recognized as specific proteins which had been identified in leguminous plants like Vicia faba, Vicia, and others (Table S2). Six proteins from the 96-kD subunit were all identified as specific proteins, and the sequence coverage ranged from 4.01 to 15.84% (Table S2).
The results from CELLO revealed various subcellular location patterns of those identified proteins. The results exhibited GroEL chaperonin, elongation factor Tu, electron transfer flavoprotein subunit alpha, alkyl hydroperoxide reductase C22 subunit, citrate synthase, convicilin, and putative sucrose-binding protein all located in the cytoplasm; phosphate ABC transporter periplasmic substratebinding protein and Htlv-1 Gp21 ectodomain MALTOSEbinding protein chimera all located periplasmic; motA/ TolQ/ExbB proton channel family protein located at the inner membrane; legumin A precursor located at the outer membrane; and LEGB7_VICFA which showed a wide distribution to the outer membrane, extracellular region, and cytoplasm (Table S2). The diverse subcellular localization of those seed proteins in faba bean imply a variety of functions.

Analysis of protein primary structure
The sequence of amino acids is the foundation of the native state of a protein, so we first analyzed the protein primary structure using Expert Protein Analysis (ExPASy, http:// expasy.org/tools) in order to precisely predict the potential function and basic structure of seed storage proteins in faba bean. As shown in Table 1, the number of amino acids in the 16 identified proteins ranged from 187 to 497 (Table 1). The maximal number and minimal number of amino acids was 497 in convicilin and 187 in alkyl hydroperoxide reductase C22 subunit, respectively, which resulted in a different molecular weight of 57501.1 and 20525.2 D, respectively. The theoretical isoelectric point of identified proteins ranged from 4.90 to 9.55, which is very important to protein purification by the isoelectric focusing method. The instability index ranged from 18.59 in electron transfer flavoprotein subunit alpha to 61.92 in the legumin A precursor. The aliphatic index ranged from 75.49 to 112.6. The grand average of hydropathicity (GRAVY) was -0.957 to 0.374. Convicilin and the alkyl hydroperoxide reductase C22 subunit showed the maximum and minimum value, respectively, in the total number of negatively/positively charged residues (Table 1).

Amino acid composition assay of different proteins
There was a vast difference among different proteins both in amino acid amount and amino acid composition. In legumes, the low content of cysteine and methionine in Table 1 The primary structure of seed storage proteins in faba bean  Htlv-1 Gp21 ectodomain MALTOSE-binding protein chimera major storage proteins may limit the nutritional quality of the seed for monogastric animals (Shewry et al. 1995), so the content of sulfur-containing amino acids (cysteine and methionine) is an important criteria when evaluating protein quality. The amino acid composition and percentage in specific proteins are shown in Table 2. As shown, the convilicin contained 62 glutamic and 46 leucine, which made up 12.5 and 9.3% of the total amino acid amount, respectively. In the legumin A precursor, the content of arginine was the highest and next was leucine, which accounted for 10 and 9.3% of the total amino acid amount, respectively. There were 51 lysine and 51 glutamic in putative sucrose-binding protein, accounting for the highest percentage at 10.6% each. Leucine was the second-most abundant amino acid and made up 10.2% of the total amino acid amount in putative sucrose-binding protein. In LEGB7_VICFA, there were 37 glutamic and 29 leucine, accounting for 11.0 and 8.7%, respectively. The sulfur-containing amino acid number in the four proteins above was 3, 8, 15, and 1, which corresponds with 0.6, 1.6, 3.2, and 0.3%, respectively. In addition, elongation factor Tu and citrate synthase showed a high sulfur-containing amino acid content with percentages of 3.9 and 4.4%, respectively. These results suggest that the proteins in faba bean are diverse in amino acid amount and composition.

Secondary structure prediction of different proteins
Protein secondary structure is the basic spatial structure of a protein. As shown in Table 3, alpha helix, extended strand, beta turn, and random coil are the most common  secondary structures in identified proteins. Abundance of alpha helices in the identified proteins ranged from 25.26% in elongation factor Tu to 58.65% in the motA/TolQ/ExbB proton channel family. Electron transfer flavoprotein subunit alpha contained the most beta turns with an abundance of 13.50%, and the content of beta turns in the motA/TolQ/ ExbB proton channel family was the lowest. Elongation factor Tu possessed the most extended strands with an abundance of 28.57%, whereas the percentage of extended strands in the motA/TolQ/ExbB proton channel family was only 11.28%. The amount of random coils ranged from 23.28% in GroEL chaperonin to 38.43% in convicilin.

Transmembrane domain prediction
To characterize a protein to be a membrane protein, we further detected the transmembrane domain of those identified proteins depending on the online severs TMHMM and HMMTOP (data from HMMTOP not shown), and acquired consistent results. The results suggest that the motA/TolQ/ExbB proton channel family protein contained three transmembrane spanning domains, consistent with the roles that the motA/TolQ/ExbB proton channel family protein undertakes in intracellular trafficking, secretion, and vesicular transport. The results also suggest that convicilin and putative sucrose-binding protein possess one transmembrane domain each (Fig. 1). There were no transmembrane domains found in GroEL chaperonin, phosphate ABC transporter periplasmic substrate-binding protein, elongation factor Tu, electron transfer flavoprotein subunit alpha, alkyl hydroperoxide reductase C22 subunit, citrate synthase, legumin A precursor, legumin A precursor, LEGB7_VICFA, or Htlv-1 Gp21 ectodomain MAL-TOSE-binding protein chimera (Fig. 1). The various transmembrane domains of proteins in faba bean seeds are closely related to the diverse functions.

Prediction of three-dimensional structure for identified proteins
Model quality estimation is an essential component of protein structure prediction, as the accuracy of a model determines its usefulness for practical applications (Biasini et al. 2014). A reliable three-dimensional structure based on homologous modeling was constructed using SWISS-MODEL and is shown in Fig. 2. The optimum templates information is listed in Table 4. We used homology modeling for structure prediction of GroEL chaperonin, elongation factor Tu, alkyl hydroperoxide reductase C22 subunit, citrate synthase, and legumin A precursor, where all the sequence identities were [60% between the target and template sequences (Table 4). Htlv-1 Gp21 ectodomain MALTOSE-binding protein chimera even shared a high sequence identity of 100% with the template (PDB ID 1mg1), which is consistent with the three-dimensional structure characteristic predicted by Kobe et al. (1999). LEGB7_VICFA shared a sequence identity of 45.90% with the template (Table 4). As shown in Table 4, we also used homology modeling for structure prediction of phosphate ABC transporter periplasmic substrate-binding protein, motA/TolQ/ExbB proton channel family protein, and putative sucrose-binding protein, where the sequence identity was 23.26, 27.44, and 34.61%, respectively. The GMEQ (global quality estimation score), in which the most likely structural similarity is the value at which the joint distribution is maximized, is an important indicator expressed as a number between zero and one (Biasini et al. 2014). GroEL chaperonin, elongation factor Tu, electron transfer flavoprotein subunit alpha, alkyl hydroperoxide reductase C22 subunit, citrate synthase, and Htlv-1 Gp21 ectodomain MALTOSE-binding protein chimera showed high GMEQ values of 0.8, 0.79, 0.8, 0.9, 0.82, and 0.99, respectively, meaning a high reliability of the hypothetical three-dimensional structure. Phosphate ABC transporter periplasmic substrate-binding protein, motA/TolQ/ExbB proton channel family protein, convicilin, legumin A precursor, putative sucrose-binding protein, and LEGB7_VICFA showed relatively low values of 0.51,0.58, 0.53, 0.73, 0.56, and 0.65, respectively. Coverage, another property of model quality, ranged from 0.71 to 1.
In addition, structural evaluation and stereochemical analyses were assessed using RAMPAGE Ramachandran plot analysis (http://mordred.bioc.cam.ac.uk/*rapper/rampage.php) (Lovell et al. 2003;Wang et al. 2016). The refined models showed good proportions of residues in favored, allowed, and outliner regions except LEGB7_VICFA (Online Resource 1 and Table 4). The models of GroEL chaperonin, elongation factor Tu, electron transfer flavoprotein subunit alpha, alkyl hydroperoxide reductase C22 subunit, citrate synthase, motA/TolQ/ExbB proton channel family protein, and Htlv-1 Gp21 ectodomain MALTOSE-binding protein chimera showed [95% of the residues in the allowed region of the Ramachandran plot (Online Resource 1 and Table 4).

Functional annotation and prediction of identified proteins
Structure similarity is more consistent than sequence similarity (Taylor and Orengo 1989), and the structures of homologous proteins are more conserved than their sequences (Chothia and Lesk 1986). The biological function cannot be predicted by comparison of sequence similarity alone (Illergard et al. 2009). In this study, protein function were explicated combined the protein structures and function annotation from the NR database. The objective was to assign a precise function to the proteins in faba bean. The major function of these proteins is listed in Online Resource 2 (Table S3).
The GroEL chaperonin was identified in the 97-, 64-, 42-, and 38-kD subunits with a relative abundance of 100, 8.4, 13.2, and 5.9%, respectively. Sequence similarity searches showed that GroEL chaperonin belongs to the HSP60 family. It was reported that GroEL chaperonin can l Htlv-1 Gp21 ectodomain MALTOSE-binding protein chimera capture the ploypeptide substrate and bind co-chaperone GroES and ATP to promote productive folding (Bartesaghi et al. 2012;Walter 2002). Rusanganwa and Gupta (1993) revealed a possible role of the GroEL or HSP60 chaperonin in the nodulation (symbiosis) and nitrogen fixation processes (Rusanganwa and Gupta 1993).
Phosphate ABC transporter periplasmic substrate-binding protein belongs to the periplasmic-binding protein type 2 family, which plays a role in the uptake of a variety of soluble substrates such as phosphate, sulfate, polysaccharides, lysine/ arginine/ornithine, and histidine. Besides transport proteins, the family includes ionotropic glutamate receptors and unorthodox sensor proteins involved in signal transduction. The above function was confirmed based on the homology modeling of 4omb, which was classified as a transport protein and was a component of the conserved bacterial phosphatespecific transport system (Neznansky et al. 2014).
Elongation factor Tu is a multifunctional protein and has been extensively studied in several species. This protein is a GTP-binding protein and plays an important role in protein biosynthesis as it promotes the binding of aminoacyl-tRNAs to ribosomes (Dale et al. 2004;Kim et al. 2000;Daviter et al. 2003;Carlin et al. 1992). The structure of Echerichia coli ribosome elongation factor Tu complex was available at \3 Å resolution by C s -corrected cryo-EM (Fischer et al. 2015).
The electron transfer flavoprotein serves as a specific electron acceptor for various mitochondrial dehydrogenases. Electron transfer flavoprotein transfers electrons to the main respiratory chain via electron transfer flavoprotein-ubiquinone oxidoreductase. Electron transfer flavoprotein is a heterodimer that consists of an alpha and a beta subunit which binds one molecule of FAD per dimer. The homologous pair of proteins (FixA/FixB) are essential for nitrogen fixation (Tsai and Saier 1995;Roberts et al. 1996). Based on the model of 2a1u, the crystal structure of this protein has been constructed (Toogood et al. 2005). Alkyl hydroperoxide reductase C22 subunit, a kind of oxidoreductase which belongs to the family of antioxidant proteins, directly reduces organic hydroperoxides in its reduced dithiol form. It is reported that alkyl hydroperoxide reductase acts as a specific alkyl hydroperoxide-scavenging enzyme for protection against radical oxygen damage (Jacobson et al. 1989), though elimination of reactive nitrogen intermediates has also been demonstrated to occur (Chen et al. 1998). These data argue for a significant role of alkyl hydroperoxide reductase in the resistance to peroxide damage.
Citrate synthase is a hexameric enzyme that catalyzes the entry of carbon into the citric acid cycle (Gerike et al. 1998;Textor et al. 1997). It is often regarded as the first enzyme in the citric acid cycle (Patton et al. 1993). It is also reported that the crystal structure of the closed form of citrate synthase from the hyperthermophilic Archaeon Pyrococcus furiosus is 1.9 Å (Russell et al. 1997).
The MotA/TolQ/ExbB proton channel family protein is mainly localized in the membrane and functions in intracellular trafficking, secretion, and vesicular transport using the motive force at the inner membrane (Jarosik and Hansen 1995;Sun and Webster 1987). Studies have also shown that the Ton subcomplex (including motA/TolQ/ ExbB) forms pH-sensitive cation-selective channels, providing insight into the mechanism by which it may harness the proton motive force to produce energy (Celia et al. 2016).
Convicilins, a set of seed storage proteins, differ from vicilins, a related group of seed storage proteins, mainly because of the presence of the N-terminal extension, an additional sequence of amino acids in the sequence that corresponds to the first exon. Convicilins have been described only in species of the legume tribe Vicieae. Legumin A precursor is a seed storage globulin found in Vicia narbonensis and provides the major nitrogen source for the developing plant (Tandang-Silvas et al. 2010).
Putative sucrose-binding protein has been described as a sucrose-binding and sucrose-transporting protein as well as an outgroup member of the large vicilin storage protein family (Hei et al. 2001).
The LEGB7_VICFA protein, found in the seeds of many leguminous and non-leguminous plants, is an important source of sulfur-containing amino acids in seed meals. The LEGB7_VICFA protein belongs to the 11S seed storage protein (globulin) family (Heim et al. 1989).
It is hypothesized that the Htlv-1 Gp21 ectodomain MALTOSE-binding protein chimera possesses dual roles of conserved fusion function and an ability to anchor diverse surface-exposed subunit structures to the virion envelope and infected cell surface (Kobe et al. 1999;Center et al. 1998).

Discussion and conclusion
Seed storage proteins are vital for plant growth and development as they (1) provide nutrients for seed respiratory metabolism, seed germination, and seedling establishment; (2) regulate physiological and biochemical reactions and metabolic processes. Studies show that the seed storage proteins of faba bean mainly consist of albumin and globulin (Utsumi et al. 1980;Bewley 1997;Li et al. 2003;Liu et al. 2012a). In this study, 16 proteins (Four GroEL molecular chaperones and 12 plant specific proteins) were identified. The convicilin, legumin A precursor, putative sucrose-binding protein, and LEGB7_-VICFC are all globulins and specifically found within leguminous plants. Meanwhile, the convicilin in the 64-kD subunit, legumin A precursor in the 42-and 47-kD subunits, putative sucrose-binding protein in the 47-kD subunit, and LEGB7_VICFA in the 38-kD subunit were found at a relatively high abundance of 91.6, 86.7, 69.2, 30.8, and 89.1% respectively, whereas the abundance of other proteins ranged from 4.6 to 12.3% except GroEL chaperonin and phosphate ABC transporter periplasmic substratebinding protein. This suggests that in faba bean seeds, globulin occupies a dominate part of seed storage proteins and provides the main nutrition and amino acid supply.
It is interesting that the relative abundance of GroEL chaperones in the 97-kD subunit reached up to 100%. This suggests that the GroEL chaperones existed as a much high amount in faba bean. It is reported that HSP60 has a 15-30% contribution in cellular protein folding and maintenance and also plays a vital role in mitochondrial protein transportation (Trivedi et al. 2016). In addition, the elevated expression level of HSP/molecular chaperones assist in the assembly of the highly abundant secretory storage proteins in the endoplasmic reticulum lumen during seed development (Trivedi et al. 2016). Taken together, that may be why there is a high abundance of GroEL chaperones in seed storage proteins.
The seed storage proteins in faba bean exhibited various functions depending on their diverse three-dimensional structures. In addition to those mentioned above, the proteins with a relatively low abundance possessed different functions such as protein transport, signal transduction, protein biosynthesis, electron transfer, peroxide resistance, catalyzing, and intracellular trafficking in order to maintain the physiological and biochemical activities in faba bean seeds. Together, protein composition and abundance is closely associated with protein function in faba bean seeds.
Amino acid composition, as the assessment standard for protein quality, often varies between different seed storage proteins. The 16 proteins identified in faba bean were mainly composed of glutamic, arginine, lysine, glycine, leucine, cysteine, and methionine, which are essential to human life. In this study, the number of sulfur-containing amino acids (cysteine and methionine), which may limit the nutritional quality of the legume seed, reached up to 15, 19, and 15 in elongation factor Tu, citrate synthase, and putative sucrose-binding protein, respectively. The percentages of sulfur-containing amino acid were 3.9, 4.4, and 3.2% out of the total number of amino acids in elongation factor Tu, citrate synthase, and putative sucrose-binding protein, respectively, which was higher than the percentages in the other proteins. This result suggests that generation of the ideal faba bean using traditional breeding methods or genetic engineering may need to increase the cysteine and methionine levels in seed storage proteins. In the future, we can purposefully promote the amount of the elongation factor Tu, citrate synthase, and putative sucrosebinding proteins which contained a relative high amount of sulfur-containing amino acids to optimize the protein quality of faba bean seeds.