Molecular thermodynamic model to predict the alpha-helical secondary structure of polypeptide chains in solution.

The native state of a protein molecule in aqueous solutions represents one of the lowest states of Gibbs energy [Anfinsen, C.B. (1973) Science 181, 223-230]. Much progress has been made about the rules of protein folding [King, J. (1989) Chem. Eng. News 67, 32-54] and the dominant forces in protein folding [Dill, K.A. (1990) Biochemistry 29, 7133-7155]. However, the quantitative contributions of different Gibbs energy terms to protein stability remains a controversial issue [Moult, J., & Unger, R. (1991) Biochemistry 30, 3816-3824]. A molecular thermodynamic model has been proposed for the Gibbs energy of folding a residue in aqueous homopolypeptides from a random-coiled state to either the alpha-helix state or the beta-sheet state [Chen, C.-C., Zhu, Y., King, J.A., & Evans, L.B. (1992) Biopolymers 32, 1375-1392]. In this work, we present a generalization of the molecular thermodynamic model for the Gibbs energy of folding natural and synthetic heteropolypeptides from random-coiled conformations into alpha-helical conformations. The generalized model incorporates the intrinsic folding potential due to residue-solvent interactions, the cooperative folding effect due to residue-residue interactions, and the location and length of alpha-helices. The utility of the model was demonstrated by examining the stability of alpha-helical conformations of a number of natural polypeptides including C-peptide (residues 1-13) and S-peptide (residues 1-20) of RNase A (bovine pancreatic ribonuclease A), the P alpha fragment in BPTI (bovine pancreatic trypsin inhibitor), and synthetic polypeptides (the copolymers of different amino acid residues) including alanine-based peptides (16 or 17 residues long) in water. The computed Gibbs energies correspond well with the experimental data on helicity. The results also accounted for the effects of amino acid substitution and temperature on the stability of alpha-helical conformations of the test polypeptides.     

Optimized atomic radii for protein continuum electrostatics solvation forces

Recently, we presented a Green's function approach for the calculation of analytic continuum electrostatic solvation forces based on numerical solutions of the finite-difference Poisson-Botzmann (FDPB) equation [Im et al., Comp. Phys. Comm. 111 (1998) 59]. In this treatment the analytic forces were explicitly defined as the first derivative of the FDPB continuum electrostatic free energy with respect to the coordinates of the solute atoms. A smooth intermediate region for the solute-solvent dielectric boundary needed to be introduced to avoid abrupt discontinuous variations in the solvation free energy and forces as a function of the atomic positions. In the present paper we extend the set of optimized radii, which was previously parametrized from molecular dynamics free energy simulations of the 20 standard amino acids with explicit solvent molecules [Nina et al., J. Phys. Chem. 101 (1997) 5239], to yield accurate solvation free energy by taking the influence of the smoothed dielectric region into account.     

Atomic solvation parameters applied to molecular dynamics of proteins in solution.

A solvation energy function for use in the molecular simulation of proteins is proposed. It is based on the accessible surface areas of atoms in the protein and on atomic solvation parameters derived from empirical vapor-to-water free energies of transfer of amino acid side-chain analogs. The energy function and its derivatives were added to the CHARMM molecular simulation program (Brooks, B.R., Bruccoleri, R.E., Olafson, B.D., States, D.J., Swaminathan, S., & Karplus, M., 1983, J. Comput. Chem. 4(2), 187-217). The effect of the added energy term was evaluated by 110 ps of molecular dynamics on the 26-residue protein melittin. The melittin monomer and tetramer were studied both with and without the added term. With the added energy term the monomer partially unfolded, while the secondary structure of the tetramer was preserved, in agreement with reported experiments (Brown, L.R., Lauterwein, J., & Wuethrich, K., 1980, Biochim. Biophys. Acta 622(2), 231-244; Lauterwein, J., Brown, L.R., & Wuethrich, K., 1980, Biochim. Biophys. Acta 622(2), 219-230).     

Continuous molecular evolution of protein-domain structures by single amino acid changes

Protein structures cluster into families of folds that can result from extremely different amino acid sequences [1]. Because the enormous amount of genetic information generates a limited number of protein folds [2], a particular domain structure often assumes numerous functions. How new protein structures and new functions evolve under these limitations remains elusive. Molecular evolution may be driven by the ability of biomacromolecules to adopt multiple conformations as a bridge between different folds [3-6]. This could allow proteins to explore new structures and new tasks while part of the structural ensemble retains the initial conformation and function as a safeguard [7]. Here we show that a global structural switch can arise from single amino acid changes in cysteine-rich domains (CRD) of cnidarian nematocyst proteins. The ability of these CRDs to form two structures with different disulfide patterns from an identical cysteine pattern is distinctive [8]. By applying a structure-based mutagenesis approach, we demonstrate that a cysteine-rich domain can interconvert between two natively occurring domain structures via a bridge state containing both structures. Comparing cnidarian CRD sequences leads us to believe that the mutations we introduced to stabilize each structure reflect the birth of new protein folds in evolution.     

Deletion analysis of the Drosophila Inscuteable protein reveals domains for cortical localization and asymmetric localization

The Drosophila Inscuteable protein acts as a key regulator of asymmetric cell division during the development of the nervous system [1] [2]. In neuroblasts, Inscuteable localizes into an apical cortical crescent during late interphase and most of mitosis. During mitosis, Inscuteable is required for the correct apical-basal orientation of the mitotic spindle and for the asymmetric segregation of the proteins Numb [3] [4] [5], Prospero [5] [6] [7] and Miranda [8] [9] into the basal daughter cell. When Inscuteable is ectopically expressed in epidermal cells, which normally orient their mitotic spindle parallel to the embryo surface, these cells reorient their mitotic spindle and divide perpendicularly to the surface [1]. Like the Inscuteable protein, the inscuteable RNA is asymmetrically localized [10]. We show here that inscuteable RNA localization is not required for Inscuteable protein localization. We found that a central 364 amino acid domain - the Inscuteable asymmetry domain - was necessary and sufficient for Inscuteable localization and function. Within this domain, a separate 100 amino acid region was required for asymmetric localization along the cortex, whereas a 158 amino acid region directed localization to the cell cortex. The same 158 amino acid fragment could localize asymmetrically when coexpressed with the full-length protein, however, and could bind to Inscuteable in vitro, suggesting that this domain may be involved in the self-association of Inscuteable in vivo.     

Exhaustive mutagenesis in silico: multicoordinate free energy calculations on proteins and peptides

We have extended and applied a multicoordinate free energy method, chemical Monte Carlo/Molecular Dynamics (CMC/MD), to calculate the relative free energies of different amino acid side-chains. CMC/MD allows the calculation of the relative free energies for many chemical species from a single free energy calculation. We have previously shown its utility in host:guest chemistry (Pitera and Kollman, J Am Chem Soc 1998;120:7557-7567)1 and ligand design (Eriksson et al., J Med Chem 1999;42:868-881)2, and here demonstrate its utility in calculations of amino acid properties and protein stability. We first study the relative solvation free energies of N-methylated and acetylated alanine, valine, and serine amino acids. With careful inclusion of rotameric states, internal energies, and both the solution and vacuum states of the calculation, we calculate relative solvation free energies in good agreement with thermodynamic integration (TI) calculations. Interestingly, we find that a significant amount of the unfavorable solvation of valine seen in prior work (Sun et al., J Am Chem Soc 1992;114:6798-6801)3 is caused by restraining the backbone in an extended conformation. In contrast, the solvation free energy of serine is calculated to be less favorable than expected from experiment, due to the formation of a favorable intramolecular hydrogen bond in the vacuum state. These monomer calculations emphasize the need to accurately consider all significant conformations of flexible molecules in free energy calculations. This development of the CMC/MD method paves the way for computations of protein stability analogous to the biochemical technique of "exhaustive mutagenesis." We have carried out just such a calculation at position 133 of T4 lysozyme, where we use CMC/MD to calculate the relative stability of eight different side-chain mutants in a single free energy calculation. Our T4 calculations show good agreement with the prior free energy calculations of Veenstra et al. (Prot Eng 1997;10:789-807)4 and excellent agreement with the experiments of Mendel et al. (Science 1992;256:1798-1802).     

Predicting protein structural classes with pseudo amino acid composition: an approach using geometric moments of cellular automaton image

A novel approach was developed for predicting the structural classes of proteins based on their sequences. It was assumed that proteins belonging to the same structural class must bear some sort of similar texture on the images generated by the cellular automaton evolving rule [Wolfram, S., 1984. Cellular automation as models of complexity. Nature 311, 419-424]. Based on this, two geometric invariant moment factors derived from the image functions were used as the pseudo amino acid components [Chou, K.C., 2001. Prediction of protein cellular attributes using pseudo amino acid composition. Proteins: Struct., Funct., Genet. (Erratum: ibid., 2001, vol. 44, 60) 43, 246-255] to formulate the protein samples for statistical prediction. The success rates thus obtained on a previously constructed benchmark dataset are quite promising, implying that the cellular automaton image can help to reveal some inherent and subtle features deeply hidden in a pile of long and complicated amino acid sequences.     

Hydrophobic potential of mean force as a solvation function for protein structure prediction

We have developed a solvation function that combines a Generalized Born model for polarization of protein charge by the high dielectric solvent, with a hydrophobic potential of mean force (HPMF) as a model for hydrophobic interaction, to aid in the discrimination of native structures from other misfolded states in protein structure prediction. We find that our energy function outperforms other reported scoring functions in terms of correct native ranking for 91% of proteins and low Z scores for a variety of decoy sets, including the challenging Rosetta decoys. This work shows that the stabilizing effect of hydrophobic exposure to aqueous solvent that defines the HPMF hydration physics is an apparent improvement over solvent-accessible surface area models that penalize hydrophobic exposure. Decoys generated by thermal sampling around the native-state basin reveal a potentially important role for side-chain entropy in the future development of even more accurate free energy surfaces.     

Cetaceans on a molecular fast track to ultrasonic hearing

The early radiation of cetaceans coincides with the origin of?their defining ecological and sensory differences [1, 2]. Toothed whales (Odontoceti) evolved echolocation for hunting 36-34 million years ago, whereas baleen whales (Mysticeti) evolved filter feeding and do not echolocate [2]. Echolocation in toothed whales demands exceptional high-frequency hearing [3], and both echolocation and ultrasonic hearing have also evolved independently in bats [4, 5]. The motor protein Prestin that drives the electromotility of the outer hair cells (OHCs) is likely to be especially important in ultrasonic hearing, because it is the vibratory response of OHC to incoming sound waves that confers the enhanced sensitivity and selectivity of the mammalian auditory system [6, 7]. Prestin underwent adaptive change early in mammal?evolution [8] and also shows sequence convergence between bats and dolphins [9, 10], as well as within bats [11]. Focusing on whales, we show for the first time that the extent of protein evolution in Prestin can be linked directly to the evolution of high-frequency hearing. Moreover, we find that independent cases of sequence convergence in mammals have involved numerous identical amino acid site replacements. Our findings shed new light on the?importance of Prestin in the evolution of mammalian hearing.     

Aminoacyl-tRNA synthetases database

Aminoacyl-tRNA synthetases (AARSs) are at the center of the question of the origin of life. They constitute a family of enzymes integrating the two levels of cellular organization: nucleic acids and proteins. AARSs arose early in evolution and are believed to be a group of ancient proteins. They are responsible for attaching amino acid residues to their cognate tRNA molecules, which is the first step in the protein synthesis. The role they play in a living cell is essential for the precise deciphering of the genetic code. The analysis of AARSs evolutionary history was not possible for a long time due to a lack of a sufficiently large number of their amino acid sequences. The emerging picture of synthetases' evolution is a result of recent achievements in genomics [Woese,C., Olsen,G.J., Ibba,M. and S?ll,D. (2000) Microbiol. Mol. Biol. Rev., 64, 202-236]. In this paper we present a short introduction to the AARSs database. The updated database contains 1047 AARS primary structures from archaebacteria, eubacteria, mitochondria, chloroplasts and eukaryotic cells. It is the compilation of amino acid sequences of all AARSs known to date, which are available as separate entries via the WWW at http://biobases.ibch.poznan.pl/aars/.     

The effect of water on the Fe(3+)/Fe(2+) reduction potential of heme

Hemeproteins can act as catalysts, oxygen carriers or electron conductors. The ferric/ferrous reduction potential E(m7) of iron in the center of the prosthetic group ranges from negative values for peroxidases to an extreme positive value for cytochrome a(3) with Hb and Mb in the middle [1]. Proteins exercise their influence on E(m7) in several ways: via substituents at the periphery of the chelate structure, via the proximal ligand, and via interaction with the surrounding medium, amino acid side chains, or polar solvents. Work on recombined proteins and 2,4-substituted free hemes documented that the first two effects are additive [2]. For the third effect, models of the dielectric media on a molecular level have been successfully applied [3-5]. E(m7) has also been empirically correlated to the degree of heme exposure to water [6-8]. The apoprotein/porphyrin and water/porphyrin interfaces are complementary since water molecules fill any empty space in the crevice and surround any pertinent part of heme outside the protein boundary. The present work links to this idea by a combination of statistical mechanics simulations and quantum mechanical calculations comparing heme in water with heme in an apolar environment. Our results show that polarization of the porphyrin pi-electron cloud by the field from water dipoles influences E(m7). The dominant effect of this and other determinates of iron electron availability is perturbations of delocalized electron density in the porphyrin chelate, reproduced by a model where the prosthetic group is treated as a disc of uniform electron density. The present work is also of interest since the interfacial energy constitutes the main barrier for heme-protein separation [9-11].     

An amino acid has two sides: a new 2D measure provides a different view of solvent exposure

The concept of amino acid solvent exposure is crucial for understanding and predicting various aspects of protein structure and function. The traditional measures of solvent exposure however suffer from various shortcomings, like for example the inability to distinguish exposed, partly exposed, buried, and deeply buried residues. This article introduces a new measure of solvent exposure called Half-Sphere Exposure that addresses many of the shortcomings of other methods. The new measure outperforms other measures with respect to correlation with protein stability, conservation among fold homologs, amino acid-type dependency and interpretation. The measure consists of the number of Calpha atoms in two half spheres around a residue's Calpha atom. Conceptually, one of the half spheres corresponds to the side chain's neighborhood, the other half sphere being in the opposite direction. We show here that the two half spheres correspond to two regions around an amino acid that are surprisingly distinct in terms of geometry and energy. This aspect of protein structure introduced here forms the basis of the Half-Sphere Exposure measure. The results strongly suggest that in many respects, a 2D measure is inherently much better suited to describe solvent exposure than the traditional 1D measures. Importantly, Half-Sphere Exposure can be calculated from the Calpha atom coordinates only, which abolishes the need for a full-atom model to calculate solvent exposure. Hence, the measure can be used in protein structure prediction methods that are based on various simplified models. Half-Sphere Exposure has great potential for use in protein structure prediction and analysis.     

Macromolecular solvation energies derived from small molecule crystal morphology.

The morphology of small molecule crystals provides a model for evaluating surface solvation energies in a system with similar packing density to that observed for amino acid residues in proteins. The solvation energies associated with the transfer of methylene and carboxyl groups between vacuum and aqueous phases are estimated to be approx. $40 and -260 cal/A2, respectively, from an analysis of the morphology of succinic acid crystals. These solvation energies predict values for contact angles in reasonable agreement with measurements determined from macroscopic monolayer surfaces. Transfer free energies between vapor and water phases for a series of carboxylic acids are also predicted reasonably well by these solvation energies, provided the surface exposure of different groups is quantitated with the molecular surface area rather than the more traditional accessible surface area. In general, molecular surfaces and molecular surface areas are seen to have important advantages for characterizing the structure and energetics of macromolecular surfaces. Crystal faces of succinic acid with the lowest surface energies in aqueous solution are characteristically smooth. Increasing surface roughness and apolarity are associated with higher surface energies, which suggests an approach for modifying the surface properties of proteins and other macromolecules.     

Energization of amino acid uptake by system A in cultured human fibroblasts

The energization of System A in cultured human fibroblasts has been studied by measuring the energy transfer from the electrochemical gradient of Na+ to the chemical gradient of the site A-specific substrate amino acid 2-methylaminoisobutyric acid. The co-transport Na+/amino acid, studied by kinetic analysis and radiochemical measurements, showed a coupling ratio of 1:1. The assessment of the Na+ electrochemical gradient in cultured adherent cells relied on the development of noninvasive procedures as follows: the membrane electrical potential was estimated from the accumulation of L-arginine at equilibrium (Bussolati, O., Laris, P. C., Nucci, F. A., Dall'Asta, V., Longo, N., Guidotti, G. G., and Gazzola, G. C. (1987) Am. J. Physiol. 253, C391-C397); the chemical gradient of Na+ was determined from spectrometric measurements of Na+. The accumulation of 2-methylaminoisobutyric acid was strongly sensitive to changes of Na+ gradient and of membrane electrical potential, indicating that the electrochemical gradient of Na+ contributed energy for the uphill transport of the amino acid through System A. Changes in the Na+ electrochemical gradient were obtained by: (i) alterations of extracellular concentration of Na+; (ii) changes of membrane electrical potential obtained by variation of extracellular [K+]; and (iii) changes of [Na+]in and membrane electrical potential upon incubation of the cells in serum-free saline solutions (Dall'Asta, V., Gazzola, G. C., Longo, N., Bussolati, O., Franchi-Gazzola, R., and Guidotti, G. G. (1986) Biochim. Biophys. Acta 860, 1-8). The correlation between the chemical gradient of 2-methylaminoisobutyric acid and the Na+ electrochemical potential followed a straight line with a yield close to the thermodynamic equilibrium, thus suggesting that the energy stored in the gradient of Na+ electrochemical potential is fully adequate to energize the intracellular accumulation of site A-reactive amino acids in human fibroblasts.     

Energetics of protein folding

The energetics of protein folding determine the 3D structure of a folded protein. Knowledge of the energetics is needed to predict the 3D structure from the amino acid sequence or to modify the structure by protein engineering. Recent developments are discussed: major factors are reviewed and auxiliary factors are discussed briefly. Major factors include the hydrophobic factor (burial of non-polar surface area) and van der Waals interactions together with peptide hydrogen bonds and peptide solvation. The long-standing model for the hydrophobic factor (free energy change proportional to buried non-polar surface area) is contrasted with the packing-desolvation model and the approximate nature of the proportionality between free energy and apolar surface area is discussed. Recent energetic studies of forming peptide hydrogen bonds (gas phase) are reviewed together with studies of peptide solvation in solution. Closer agreement is achieved between the 1995 values for protein unfolding enthalpies in vacuum given by Lazaridis-Archontis-Karplus and Makhatadze-Privalov when the solvation enthalpy of the peptide group is taken from electrostatic calculations. Auxiliary factors in folding energetics include salt bridges and side-chain hydrogen bonds, disulfide bridges, and propensities to form alpha-helices and beta-structure. Backbone conformational entropy is a major energetic factor which is discussed only briefly for lack of knowledge.     

Complete amino acid sequence of the heavy-chain variable region from an A/J mouse antigen-nonbinding monoclonal antibody bearing the predominant arsonate idiotype

The 1F6 hybridoma protein, exhibiting the predominant cross-reactive idiotype (CRI) associated with the immune response to p-azophenylarsonate in A/J mice but failing to bind the hapten arsonate, was elicited following immunization with rat anti-CRI [Wysocki, L.J., & Sato, V. (1981) Eur. J. Immunol. 11, 832-839]. The dissociation of idiotype and antigen binding in this hybridoma provides an opportunity to determine structural features involved in antigen binding and idiotypic sites. The complete heavy-chain variable region (VH) amino acid sequence was obtained by automated Edman degradation of the intact chain and fragments due to CNBr cleavage, trypsin digestion, mild acid hydrolysis, and carboxypeptidase A digestion of a CNBr fragment. Comparison of the CRI+ arsonate-nonbinding 1F6 sequence with the CRI+ germ-line VH gene sequence reveals that the 1F6 heavy chain differs from the germ-line-encoded amino acid sequence at seven positions within VH [Siekevitz, M., Gefter, M. L., Brodeur, P., Riblet, R., & Marshak-Rothstein, A. (1982) Eur. J. Immunol. 12, 1023-1032]. The 1F6 VH appears to arise from the CRI+ germ-line VH by somatic mutation at at least seven amino acid residues, each of which could be due to a single nucleotide base change. The diversity (D) gene-encoded segment of 1F6 is similar to that of the CRI+ antigen-binding hybridoma 36-65 except for two amino acid substitutions. Further, the idiotype (CRI) is preserved despite use of a JH4 gene segment in 1F6 as compared to JH2 in all CRI+ arsonate-binding hybridomas examined to date.(ABSTRACT TRUNCATED AT 250 WORDS)     

Prediction of membrane protein types from sequences and position-specific scoring matrices

Membrane protein plays an important role in some biochemical process such as signal transduction, transmembrane transport, etc. Membrane proteins are usually classified into five types [Chou, K.C., Elrod, D.W., 1999. Prediction of membrane protein types and subcellular locations. Proteins: Struct. Funct. Genet. 34, 137-153] or six types [Chou, K.C., Cai, Y.D., 2005. J. Chem. Inf. Modelling 45, 407-413]. Designing in silico methods to identify and classify membrane protein can help us understand the structure and function of unknown proteins. This paper introduces an integrative approach, IAMPC, to classify membrane proteins based on protein sequences and protein profiles. These modules extract the amino acid composition of the whole profiles, the amino acid composition of N-terminal and C-terminal profiles, the amino acid composition of profile segments and the dipeptide composition of the whole profiles. In the computational experiment, the overall accuracy of the proposed approach is comparable with the functional-domain-based method. In addition, the performance of the proposed approach is complementary to the functional-domain-based method for different membrane protein types.     

Complete amino acid sequence of human liver cytosolic alanine aminotransferase (GPT) determined by a combination of conventional and mass spectral methods

The complete amino acid sequence of human liver cytosolic alanine aminotransferase (GPT) (EC 2.6.1.2) is presented. Two primary sets of overlapping fragments were obtained by cleavage of the pyridylethylated protein at methionyl and lysyl bonds with cyanogen bromide and Achromobacter protease I, respectively. Isolated peptides were analyzed with a protein sequencer or with a plasma desorption time of flight mass spectrometer and placed in the sequence on the basis of their molecular mass and homology to the sequence of rat GPT. The protein was found to be acetylated at the amino terminus and contained 495 amino acid residues. The Mr of the subunit was calculated to be 54,479, which was in good agreement with a Mr of 55,000 estimated by SDS-PAGE, and also indicated that the active enzyme with a Mr of 114,000 was a homodimer composed of two identical subunits. The amino acid sequence is highly homologous to that of rat GPT (87.9% identity) recently determined [Ishiguro, M., Suzuki, M., Takio, K., Matsuzawa, T., & Titani, K. (1991) Biochemistry 30, 6048-6053]. All of the crucial amino acid residues are conserved in human GPT, which seem to be hydrogen bonding to pyridoxal 5'-phosphate in rat GPT by the sequence homology to other alpha-aminotransferases with known tertiary structures.     

In silico exploration of the fructose-6-phosphate phosphorylation step in glycolysis: genomic evidence of the coexistence of an atypical ATP-dependent along with a PPi-dependent phosphofructokinase in Propionibacterium freudenreichii subsp. shermanii

We performed a detailed bioinformatic study of the catalytic step of fructose-6-phosphate phosphorylation in glycolysis based on the raw genomic draft of Propionibacterium freudenreichii subsp. shermanii (P. shermanii) ATCC9614 [Meurice et al., 2004]. Our results provide the first in silico evidence of the coexistence of genes coding for an ATP-dependent phosphofructokinase (ATP-PFK) and a PPi-dependent phosphofructokinase (PPi-PFK), whereas the fructose-1,6-bisphosphatase (FBP) and ADP-dependent phosphofructokinase (ADP-PFK) are absent. The deduced amino acid sequence corresponding to the PPi-PFK (AJ508922) shares 100% similarity with the already characterised propionibacterial protein (P29495; Ladror et al., 1991]. The unexpected ATP-PFK gene (AJ509827) encodes a protein of 373 aa which is highly similar (50% positive residues) along at least 95% of its sequence length to different well-characterised ATP-PFKs. The characteristic PROSITE pattern important for the enzyme function of ATP-PFKs (PS00433) was conserved in the putative ATP-PFK sequence: 8 out of 9 amino acid residues. According to the recent evolutionary study of PFK proteins with different phosphate donors [Bapteste et al., 2003], the propionibacterial ATP-PFK harbours a G104-K124 residue combination, which strongly suggested that this enzyme belongs to the group of atypical ATP-PFKs. According to our phylogenetic analyses the amino acid sequence of the ATP-PFK is clustered with the atypical ATP-PFKs from group III of the Siebers classification [Siebers et al., 1998], whereas the expected PPi-PFK protein is closer to the PPi-PFKs from clade P [Müller et al., 2001]. The possible significance of the co-existence of these two PFKs and their importance for the regulation of glycolytic pathway flux in P. shermanii is discussed.     

Streptococcal phosphoenolpyruvate-sugar phosphotransferase system: amino acid sequence and site of ATP-dependent phosphorylation of HPr

The amino acid sequence of histidine-containing protein (HPr) from Streptococcus faecalis has been determined by direct Edman degradation of intact HPr and by amino acid sequence analysis of tryptic peptides, V8 proteolytic peptides, thermolytic peptides, and cyanogen bromide cleavage products. HPr from S. faecalis was found to contain 89 amino acid residues, corresponding to a molecular weight of 9438. The amino acid sequence of HPr from S. faecalis shows extended homology to the primary structure of HPr proteins from other bacteria. Besides the phosphoenolpyruvate-dependent phosphorylation of a histidyl residue in HPr, catalyzed by enzyme I of the bacterial phosphotransferase system, HPr was also found to be phosphorylated at a seryl residue in an ATP-dependent protein kinase catalyzed reaction [Deutscher, J., & Saier, M. H., Jr. (1983) Proc. Natl. Acad. Sci. U.S.A. 80, 6790-6794]. The site of ATP-dependent phosphorylation in HPr of S. faecalis has now been determined. [32P]P-Ser-HPr was digested with three different proteases, and in each case, a single labeled peptide was isolated. Following digestion with subtilisin, we obtained a peptide with the sequence -(P)Ser-Ile-Met-. Using chymotrypsin, we isolated a peptide with the sequence -Ser-Val-Asn-Leu-Lys-(P)Ser-Ile-Met-Gly-Val-Met-. The longest labeled peptide was obtained with V8 staphylococcal protease. According to amino acid analysis, this peptide contained 36 out of the 89 amino acid residues of HPr. The following sequence of 12 amino acid residues of the V8 peptide was determined: -Tyr-Lys-Gly-Lys-Ser-Val-Asn-Leu-Lys-(P)Ser-Ile-Met-.(ABSTRACT TRUNCATED AT 250 WORDS)