Circular dichroism analysis of the secondary structures of bovine blood coagulation factor IX, factor X, and prothrombin

Analysis of the far-ultraviolet circular dichroism spectrum of bovine blood coagulation factor IX reveals the presence of approximately 14% helical structures 26% -sheets, 20% -turns, and 40% coils. These values are essentially the same for the activation products of this zymogen, factor IXa and factor IXa. Similar analysis for bovine factor X permits calculation of these secondary structural as approximately 11% helices, 31% -structures, 22% -turns, and 36% random structures. Bovine prothrombin contains approximately 12% helical structures, 35% -structures, 24% -turns, and 29% coils. None of these values is substantially altered as a result of increase of thepH from 7.4 to 10.5, or upon addition of Ca²⁺ to a concentration of at least 20 mM. Analysis of the near-ultraviolet spectra of factor IX and prothrombin suggests that several aromatic amino acid residues and the disulfide bond present in their -carboxyglutamic acid-containing regions are exposed to solvent and are perturbed by the abovepH adjustment and Ca²⁺ addition. Similar effects are observed in the case of factor X; in addition, the Trp residue at the amino terminus of the heavy chain appears to be influenced by the abovepH alteration. The results reported in this paper show that these vitamin K-dependent blood coagulation proteins are similar in their ordered secondary structures, which are dominated by -sheets and -turns. Their overall secondary structures are not influenced by Ca²⁺ binding and are stable to alkalinepH changes. However, these same environmental alterations appear to be effective probes of aromatic residues in the -carboxyglutamic acid regions.     

Distance-constraint approach to higher-order structures of globular proteins with empirically determined distances between amino acid residues

An analysis of higher-order structures of globular proteins by means of a distance-constraint approach is presented. Conformations are generated for each of 21 test proteins of small and medium sizes by optimizing an objective functionf=w _ij(d _ij–d _ij)², whered _ij is a distance between residuesi andj in a calculated conformation, d _ij is an assigned distance to the (ij) pair of residues which is determined based on the statistics of known three-dimensional structures of 14 proteins in the earlier study, andw _ij is a weighting factor. d _ij involves information about hydrophobicity and hydrophilicity of each amino acid residue and about connectivity of a polypeptide chain. In these calculations, only the amino acid sequence is used as input data specific to a calculated protein. With respect to higher-order structures regenerated in the optimized conformations, the following properties are analyzed: (a) N14 of a residue, defined as the number of residues surrounding the residue located within a sphere of radius of 14 Å; (b) root-mean-square differences of the global and local conformations from the corresponding X-ray conformations; (c) distance profiles in the short and medium ranges; and (d) distance maps. The effects of supplementary information about locations of secondary structures and disulfide bonds are also examined to discuss the potential ability of this methodology to predict the three-dimensional structures of globular proteins.     

Secretion of Escherichia coli DsbA and DsbC proteins from Brevibacillus choshinensis: stimulation of human epidermal growth factor production

DsbA and DsbC, members of the thioredoxin super-family of redox proteins, which are expressed in the periplasmic space of Escherichia coli, were cloned into and successfully secreted from Brevibacillus choshinensis at 100 g ml^–1. Both proteins were active in exchanging disulfide bonds of bovine insulin in vitro. Furthermore, DsbA secreted by B. choshinensis promoted the conversion of non-native human epidermal growth factor to the native form.     

Circular dichroism studies on conformational changes in protein molecules upon adsorption on ultrafine polystyrene particles

The conformational changes in well-characterized model proteins [bovine ribonuclease A (RNase A), horseradish peroxidase, sperm-whole myoglobin, human hemoglobin, and bovine serum albumin (BSA)] upon adsorption on ultrafine polystyrene (PS) particles have been studied using circular dichroism (CD) spectroscopy. These proteins were chosen with special attention to molecular flexibility. The ultrafine PS particles were negatively charged and have average diameters of 20 or 30 nm. Utilization of these ultrafine PS particles makes it possible to apply the CD technique to determine the secondary structure of proteins adsorbed on the PS surface. Effects of protein properties and adsorption conditions on the extent of the changes in the secondary structure of protein molecules upon adsorption on ultrafine PS particles were studied. The CD spectrum changes upon adsorption were significant in the "soft" protein molecules (myoglobin, hemoglobin, and BSA), while they were insingnificant in the "rigid" proteins (RNase A and peroxidase). The soft proteins sustained a marked decrease in alpha-helix content upon adsorption. Moreover, the native alpha-helix content, which is given as the percentage of the alpha-helix content in the free proteins, of adsorbed BSA was found to decrease with decreasing pH and increase with increasing adsorbed amount. These observations confirm some well-known hypotheses for the confirmational chages in protein molecules upon adsorption. (c) 1992 John Wiley & Sons, Inc.     

Yeast-derived recombinant human insulin-like growth factor I: Production, purification, and structural characterization

Recombinant human insulin-like growth factor I (IGF-I) is efficiently expressed and secreted fromSaccharomyces cerevisiae using a yeast -factor leader to direct secretion. However, approximately 10–20% of the IGF-I was in a monomeric form, the remaining materials being disulfide-linked aggregates. When the purified material was subjected to reverse-phase high-performance liquid chromatography (rp-HPLC), it gave two doublet peaks, I and II. Upon reduction, doublet peaks I and II converged to one doublet peak. This suggests that peaks I and II result from different disulfide structures, and the doublet feature of each peak results from other causes. Different disulfide structures between peaks I and II were also suggested from the near UV circular dichroism of these proteins. Only the peak II was biologically active, indicating that peak II has the correct disulfide structure. Concanavalin A affinity chromatography of the purified peak II doublet showed binding of the subpeak with an earlier rp-HPLC retention time, indicating that it was glycosylated. Sequence analysis of tryptic peptides suggested that Thr²⁹ was the site of glycosylation. Site-directed mutagenesis was used to convert Thr²⁹ to Asn²⁹. This substitution reduced, but did not eliminate IGF-I glycosylation, suggesting additional glycosylation sites. The site of carbohydrate addition was consistent with the model that O-glycosylations occur on hydroxyl amino acids near proline residues in -turns.     

Predicting secondary structures of membrane proteins with neural networks

Back-propagation, feed-forward neural networks are used to predict the secondary structures of membrane proteins whose structures are known to atomic resolution. These networks are trained on globular proteins and can predict globular protein structures having no homology to those of the training set with correlation coefficients (C) of 0.45, 0.32 and 0.43 for a-helix, -strand and random coil structures, respectively. When tested on membrane proteins, neural networks trained on globular proteins do, on average, correctly predict (Q_i) 62%, 38% and 69% of the residues in the -helix, -strand and random coil structures. These scores rank higher than those obtained with the currently used statistical methods and are comparable to those obtained with the joint approaches tested so far on membrane proteins. The lower success score for -strand as compared to the other structures suggests that the sample of -strand patterns contained in the training set is less representative than those of a-helix and random coil. Our analysis, which includes the effects of the network parameters and of the structural composition of the training set on the prediction, shows that regular patterns of secondary structures can be successfully extrapolated from globular to membrane proteins. Correspondence to: R. Casadio     

Complete nucleotide sequences of seven eubacterial genes coding for the elongation factor Tu: functional,structural and phylogenetic evaluations

The nucleotide sequences of cloned genes coding for the elongation factor Tu of seven eubacteria have been determined. These genes were fiom Anacystis nidulans, Bacillus subtilis, Bacteroides fragilis, Deinonema spec., Pseudomonas cepacia, Shewanella putrefaciens and Streptococcus oralis. The primary structures of the genes were compared to the available sequences of prokaryotic elongation factors Tu and eukaryotic elongation factors 1 alpha. A conservation profile was determined for homologous amino acid residues. Sites of known or putative functions are usually located at highly conserved positions or within highly conserved sequence stretches. The aligned 24 amino acid sequences were used as basis for a phylogenetic analysis. The phylogenetic tree corroborates the kingdom as well as phylum concept deduced from 16S rRNA data.Abbreviations EF-Tu elongation factor Tu - GDP guanosine 5-diphosphate - GTP guanosine 5-triphosphate; tuf gene, gene coding for elongation factor Tu     

Effect of lysine modification on the secondary structure of ovalbumin

In order to determine the effect of chemical modification of the -amino groups on the secondary structure of ovalbumin, we prepared six acetylated (17, 36, 54, 70, 82, and 98%) and four succinylated derivatives (25, 50, 72, and 97%) of the protein. Native ovalbumin and the acylated derivatives were homogeneous as revealed by the electrophoretic pattern. The UV-absorption and fluorescence spectra changed progressively with the extent of modification. However, circular dichroic (CD) studies indicated that acylation of 15 of the 20 lysine residues had little effect on the secondary structure of ovalbumin. Acylation of the remaining five lysine residues resulted in a fairly severe change in the secondary structure. The -helical content decreased from about 31% in the native state to 16.5% in the 97% succinylated ovalbumin and to 21.5% in the 98% acetylated derivative. A comparison of these data with the spectral and hydrodynamic data of Qasim and Salahuddin (1978) suggested that the secondary structure of ovalbumin is more resistant to acylation than is the tertiary structure and, thus, the tertiary and the secondary structures are, to some extent, mutually independent. Raising thepH to 11.2 did not alter the secondary structure of ovalbumin and increasing the ionic strength by more than 20-fold did not reverse the loss of helical structure in 97% succinylated protein. These two observations suggest that the change in secondary structure upon maximal acylation may not only involve electrostatic effects, but also certain other factors, such as steric hindrance due to the entering bulky groups.     

Prediction of probable pathways of folding in globular proteins

A method is described for the prediction of probable folding pathways of globular proteins, based on the analysis of distance maps. It is applicable to proteins of unknown spatial structure but known amino acid sequence as well as to proteins of known structure. It is based on an objective procedure for the determination of the boundary of compact regions that contain high densities of interresidue contacts on the distance map of a globular protein. The procedure can be used both with contact maps derived from a known three-dimensional protein structure and with predicted contact maps computed by means of a statistical procedure from the amino acid sequence alone. The computed contact map can also be used to predict the location of compact short-range structures, viz. -helices and -turns, thereby complementing other statistical predictive procedures. The method provides an objective basis for the derivation of a theoretically predicted pathway of protein folding, proposed by us earlier [Tanaka and Scheraga (1977) Macromolecules10, 291–304; Némethy and Scheraga (1979) Proc. Natl. Acad. Sci., U.S.A.76, 6050–6054].     

Resonance assignments, secondary structure and topology of leukaemia inhibitory factor in solution

The chemical shift assignments and secondary structure of a murine–human chimera,MH35, of leukaemia inhibitory factor (LIF), a 180-residue protein of molecular mass 20 kDa,have been determined from multidimensional heteronuclear NMR spectra acquired on auniformly 13C,15N-labelled sample. Secondary structure elements were defined on the basisof chemical shifts, NH-CH coupling constants, medium-range NOEs and the location ofslowly exchanging amide protons. The protein contains four -helices, the relativeorientations of which were determined on the basis of long-range, interhelical NOEs. The fourhelices are arranged in an up-up-down-down orientation, as found in other four-helical bundlecytokines. The overall topology of MH35-LIF is similar to that of the X-ray crystallographicstructure for murine LIF [Robinson et al. (1994) Cell, 77, 1101–1116]. Differencesbetween the X-ray structure and the solution structure are evident in the N-terminal tail, wherethe solution structure has a trans-Pro17 compared with the cis-Pro17 found in the crystalstructure and the small antiparallel -sheet encompassing residues in the N-terminus andCD loop in the crystal structure is less stable.     

Structure-function studies of human ciliary neurotrophic factor

Ciliary neurotrophic factor (CNTF) is a polypeptide that promotes the survival and/or differentiation of a number of neural cell types. Here we present a structural and functional analysis of the human CNTF molecule. Variant proteins were synthesized byEscherichia coli trnsformmed with mutant cDNA constructs, and purified by SDS-polyacrylamide gel electrophoresis and reverse phase high pressure liquid chromatography. Most variant CNTF proteins lacked neurotrophic activity, but two N-and C-terminal deletions (2–14 and 173–200, respectively) actually displayed a several-fold increase in specific activity. Loss of biological activity was accompanied by changes in the alphahelical nature of CNTF as measured by circular dichroism. These data strengthen the proposed similarity between CNTF and the family of hematopoietic cytokines.     

Structural domains of phytochrome deduced from homologies in amino acid sequences

A method of semiempirical identification of structural domains is proposed. The procedure is based on the comparison of amino acid sequences in groups of homologous proteins. This approach was tested using 32 known protein sequences from different cytochromeb ₅, cytochromec, lysozyme, hemoglobin, and myoglobin proteins. The method presented was able to identify all structural domains of these reference proteins. A consensus secondary structure provided information on structural content of these domains predicting correctly 21 of 23 (91%) of -helices. We applied this method to six homologous phytochrome sequences fromAvena, Arabadopsis, Cucurbita, Maize, Oryza, andPisum. Some of the identified domains can be assigned to the known tertiary structure categories. For example, an / domain is localized in the region known to stabilize the phytochrome chromophore in the red light absorbing form (Pr). One -helical and one / domains are localized in regions important for the chromophore stabilization in the far-red absorbing form (Pfr). From an analysis of noncovalent interaction patterns in another domain it is proposed that a phytochrome dimer contact involves two segments localized between residues 730 and 821 (using numbering of aligned sequences). Also, a possible antiparallel -sheet structure of this region has been suggested. According to this model, the long axis of the interacting structures is perpendicular to a twofold symmetry axis of the phytochrome dimer.     

A use of Ramachandran potentials in protein solution structure determinations

A strategy is developed to use database-derived - constraints during simulated annealing procedures for protein solution structure determination in order to improve the Ramachandran plot statistics, while maintaining the agreement with the experimental constraints as the sole criterion for the selection of the family. The procedure, fully automated, consists of two consecutive simulated annealing runs. In the first run, the database-derived - constraints are enforced for all aminoacids (but prolines and glycines). A family of structures is then selected on the ground of the lowest violations of the experimental constraints only, and the - values for each residue are examined. In the second and final run, the database-derived - constraints are enforced only for those residues which in the first run have ended in one and the same favored - region. For residues which are either spread over different favored regions or concentrated in disallowed regions, the constraints are not enforced. The final family is then selected, after the second run, again only based on the agreement with the experimental constraints. This automated approach was implemented in DYANA and was tested on as many as 12 proteins, including some containing paramagnetic metals, whose structures had been previously solved in our laboratory. The quality of the structures, and of Ramachandran plot statistics in particular, was notably improved while preserving the agreement with the experimental constraints.     

Molecular dynamics simulations of epidermal growth factor and transforming growth factor-α structures in water

AMBER v. 4.1 force field in 1.5 ns NPT molecular dynamics simulations of murine epidermal growth factor (mEGF), human epidermal growth factor (hEGF), and human transforming growth factor-α (hTGF-α) structures with explicit TIP3P solvation were used to investigate differences in backbone stability, changes in secondary structure, interdomain flexibility, and weakly polar interactions. Backbone root mean square deviations of sections of each peptide show that the most stable regions in mEGF and hEGF are the A-, B-, and C-loops, whereas the most stable regions in hTGF-α are the A- and B-loops. The secondary structure in the B-loops of mEGF and hEGF differ significantly from the nuclear magnetic resonance (NMR) structures of mEGF and hEGF. The position and type of turns in the B-loop of mEGF and hEGF increase the interstrand distance of the antiparallel β-sheets thereby disrupting their structure. The interdomain flexibility of simulated hTGF-α structure is greater than in either mEGF or hEGF. The ϕ, ψ dihedrals of hTGF-α occupy two distinct populations of phase space corresponding to either a C or an α-helical conformation. This change in dihedral angle is stabilized by Phe¹⁵ with Arg⁴² and Phe¹⁷ with Arg⁴² N-π weakly polar interactions that are present only in hTGF-α but not in mEGF or hEGF. Proteins 33:396–407, 1998. © 1998 Wiley-Liss, Inc.     

Secondary structural studies of bovine caseins: structure and temperature dependence of beta-casein phosphopeptide (1-25) as analyzed by circular dichroism,FTIR spectroscopy,and analytical ultracentrifugation

The defining structural feature of all of the caseins is their common phosphorylation sequence. In milk, these phosphoserine residues combine with inorganic calcium and phosphate to form colloidal complexes. In addition, nutritional benefits have been ascribed to the phosphopeptides from casein. To obtain a molecular basis for the functional, chemical, and biochemical properties of these casein peptides, the secondary structure of the phosphopeptide of bovine -casein (1–25) was reexamined using Fourier transform infrared (FTIR) and circular dichroism (CD) spectroscopies. Both methods predict secondary structures for the peptide which include polyproline II elements as well as -extended sheet and turn-like elements. These structural elements were highly stable from 5° to 70°C. Reexamination of previously published ¹H NMR data using chemical shift indices suggests structures in accord with the CD and FTIR data. Dephosphorylation showed little or no secondary structural changes, as monitored by CD and FTIR, but the modified peptide demonstrated pronounced self-association. The polymers formed were not highly temperature sensitive, but were pressure sensitive as judged by analytical ultracentrifugation at selected rotor speeds. Molecular dynamics (MD) simulations demonstrated relatively large volume changes for the dephosphorylated peptide, in accord with the pressure dependent aggregation observed in the analytical ultracentrifuge data. In contrast the native peptide in MD remained relatively rigid. The physical properties of the peptide suggest how phosphorylation can alter its biochemical and physiological properties.     

Investigation of differences in the tertiary structures of food proteins by small-angle X-ray scattering

Summary With current emphasis in bioengineering on developing new and better structure-function relationships for proteins (e.g., the need for predictability of expected properties prior to cloning), practical and reliable methodology for providing characterization of appropriate features has become of increasing importance. The most potent and detailed technique, X-ray crystallography, has severe limitations: it is so demanding and time-consuming that X-ray coordinates are frequently unavailable for materials of interest; its data relate to static and essentially unhydrated structures, whereas proteins exhibit a variety of dynamic features and function in an aqueous environment; and many proteins of technological importance may never be crystallized. Small-angle X-ray scattering, however, is particularly suitable as a methodology that can provide a substantial number of significant geometric parameters consistent with crystallographic results, that can readily show tertiary structural changes occurring under varying conditions, and that can deal with solutions and gels. Results are presented here from small-angle X-ray scattering investigations of the apo and holo forms of chicken egg-white riboflavin-binding protein, chicken egg-white lysozyme, bovine milk-whey -lactalbumin and -lactoglobulin, and bovine ribonuclease. We utilize these observations to compare tertiary structures of these proteins as well as conformational changes in these structures, and to provide a basis for discussion of their physical and biological significance.Agricultural Research Service, U.S. Department of Agriculture. Reference to brand or firm name does not constitute endorsement by the U.S. Department of Agriculture over others of a similar nature not mentioned.     

Conserved and nonconserved structures in the secretory proteins encoded in the Balbiani ring genes ofChironomus tentans

Summary The large, repetitive Balbiani ring (BR) genes, BR 1, 2, and 6, inChironomus tentans originated from a short ancestral sequence and have all evolved according to analogous amplification schemes. We analyzed the structures of the BR-encoded secretory proteins and defined the parts that have been conserved during the evolutionary process. The BR products show striking similarities, with the BR 1 and BR 2 products being more similar to each other than to the BR 6 product. In the constant (C) region of the repeat units, 7 of the 30 amino acid residues are strictly conserved; 4 of these are the cysteine residues. The subrepeat (SR) regions of all the BR products are dominated by repeated tripeptide elements rich in proline and charged amino acid residues. Most of the amino acid replacements in both regions are conservative. Secondary structure predictions suggested that the C regions of the BR 1 and BR2 products have several elements of secondary structure: an -helix, a -strand, and one or two reverse turns, as in globular structures. The prediction for the C region of the BR 6 product is similar but lacks a -strand. The predictions for the intervening SR regions appear less conclusive, but are clearly different from those for the C regions, and suggest regular structures not differing in their conformational elements. The SR regions evolved from an ancestor sequence similar to the C region; thus, the BR products seem to represent an example of evolution from one structure to two differently folded products. It is proposed that the alignment and polymerization of the long BR proteins could be promoted by the repetitive structure of the molecules, due to the possibility of forming disulfide bridges between half-cystine residues and electrostatic interactions between the charged residues of the SR regions. The divergence among the BR products is discussed in relation to possible functional differences among the members of the BR gene family.     

Aminosalicylic acid reduces the antiproliferative effect of hyperglycaemia,advanced glycation endproducts and glycated basic fibroblast growth factor in cultured bovine aortic endothelial cells: Comparison with aminoguanidine

Hyperglycaemia reduces proliferation of bovine aortic endothelial cells in vitro. A similar effect in vivo may contribute to long-term complications of diabetes such as impaired wound-healing and retinopathy.We report the effect of increased glucose concentrations, glycated basic fibroblast growth factor (FGF-2) and bovine serum albumin-derived advanced glycation endproducts (BSA-AGE) on the proliferation of bovine aortic endothelial cells.Glucose (30 and 50 mmol/l) had an antiproliferative effect on endothelial cells. This effect may be mediated through reduced mitogenic activity of FGF-2. The glycation of FGF-2 with 250 mmol/l glucose-6-phosphate led to reduced mitogenic activity compared to native FGF-2. BSA-AGE at concentrations of 10, 50 and 250 g/ml had an antiproliferative effect on cultured endothelial cells.Aminosalicylic acid at a concentration of 200 mol/l proved to be more effective than equimolar concentrations of aminoguanidine in protecting endothelial cells against the antiproliferative effects of both high (30 mmol/l) glucose and 50 g/ml BSA-AGE. FGF-2 glycated in the presence of 4 mmol/l aminosalicylic acid or aminoguanidine retained mitogenic activity compared to that glycated in their absence.Compounds like aminoguanidine and, in particular, aminosalicylic acid protect endothelial cells against glucose-mediated toxicity and may therefore have therapeutic potential.     

Quenching by acrylamide and temperature of a fluorescent probe attached to the active site of Ribonuclease

The fluorescence properties of ribonuclease labelled at its active site with N-(iodoacetylamino)-ethyl-5-naphthylamine-1-sulfonic acid have been studied at different temperatures and in the presence of acrylamide. The rate constant for the quenching of the fluorescence of labelled ribonuclease by acrylamide is apparently not limited by the accessibility of the probe: similar values are obtained for the native and denatured states of the protein. Instead, acrylamide seems to be a rather inefficient quencher of this fluorescent group ((acetamidoamino) ethyl-5-naphtylamine-1-sulfonic acid), as shown by non-linear Stern-Volmer representations, biphasic decay kinetics, and a low value of the rate constant.The fluorescence intensity of the native state of the labelled protein is highly sensitive to temperature and exhibits a 20% decrease for an increase of temperature of from 10°C to 30°C, independent of solvent viscosity. This thermal quenching is specific for the native conformation and disappears when the protein is unfolded. When the fluorescence life-time of the label is shortened by addition of acrylamide, the effect of temperature becomes identical for native and unfolded structures. This suggests that the cause of the thermal quenching is the presence of conformational fluctuations within the native protein which apparently take place in the time range from 35 to 200 ns.Abbreviations used 1,5-IAEDANS N-(iodoacetylamino)ethyl-5-naphthylamine-1-sulfonic acid - AEDANS (acetamidoamine)-ethyl-5-naphthylamine-1-sulfonic acid - RNase bovine pancreatic ribonuclease - AEDANS-RNase RNase labelled with AEDANS - ME-AEDANS (hydroxyethylthioacetamido)ethyl-5-naphthylamine-1-sulfonic acid: the product of the reaction between 1,5-IAEDANS and -mercaptoethanol (Hudson and Weber 1973) - Gu-HCl guanidine hydrochloride     

Comparison of the three-dimensional molecular models of bovine submicellar caseins with small-angle X-ray scattering. Influence of protein hydration

To test the applicability of two energy-minimized, three-dimensional structures of the bovine casein submicelle, theoretical small-angle X-ray scattering curves in the presence and absence of water were compared to experimental data. The published method simulates molecular dynamics of proteins in solution by employing adjustable Debye-Waller temperature factors (B factors) for the protein, for the solvent, and for protein-bound water. The programs were first tested upon bovine pancreatic trypsin inhibitor beginning with its known X-ray crystal structure. To approximate the degree of protein hydration previously determined by NMR relaxation experiments (0.01 g water/g protein), 120 water molecules were docked into the large void of the-casein portion of the structure for both the symmetric and asymmetric casein submicelle models. To approximate hydrodynamic hydration (0.244 g water/g protein), 2703 water molecules were added to each of the above structures using the droplet algorithm in the Sybyl molecular modeling package. All structures were then energy-minimized and their solvation energies calculated. Theoretical small-angle X-ray scattering curves were calculated for all unhydrated and hydrated structures and compared with experimentally determined scattering profiles for submicellar casein. Best results were achieved with the 120-bound-water structure for both the symmetric and asymmetric submicelle models. Comparison of results for the protein submicelle models with those for the theoretical and literature values of bovine pancreatic trypsin inhibitor demonstrates the applicability of the methodology.Reference to a brand or firm name does not constitute endorsement by the U.S. Department of Agriculture over others of a similar nature not mentioned.