Protein parameters

What are protein parameters?

Protein parameters are quantitative biochemical properties calculated from amino acid sequences to characterize protein behavior. These computational metrics include molecular weight, isoelectric point, stability indices, and hydrophobicity measures that provide insights into protein structure and function.

Parameters are derived from primary sequence data using established algorithms based on individual amino acid contributions. The calculations encompass several categories:

• Physical properties: molecular weight and sequence length
• Electrochemical properties: isoelectric point
• Stability measures: instability index, aliphatic index for thermostability
• Hydrophobicity measures: GRAVY score and residue-class percentages
• Spectroscopic properties: extinction coefficients (reduced and oxidized)
• Compositional properties: aromaticity, secondary structure fractions (helix, turn, sheet)

Prediction reliability depends on experimental data quality used in computational models. Most algorithms utilize experimental datasets from the 1960s-1990s, with periodic refinements from newer structural data.

Molecular weight

Molecular weight is the total mass of a protein molecule, expressed in Daltons (Da) or kilodaltons (kDa). It equals the sum of all amino acid residue masses plus 18.015 Da for the water molecule accounting for free amino and carboxyl termini.

The molecular weight follows this formula:

$$MW = \sum_{i=1}^{n} m_{residue,i} + 18.015\ \mathrm{Da}$$

where $n$ is the total number of amino acids and $m_{residue,i}$ represents the mass of the $i$-th amino acid residue.

Standard atomic masses follow IUPAC atomic weights. Non-standard amino acids or post-translational modifications require manual mass corrections.

Molecular weight applications include SDS-PAGE migration prediction, mass spectrometry interpretation, enzyme assay calculations, and size-exclusion chromatography profiling.

Isoelectric point

The isoelectric point (pI) is the pH at which a protein has zero net charge. Below pI, proteins are positively charged through protonation of basic residues (lysine, arginine, histidine). Above pI, proteins become negatively charged as acidic residues (aspartic acid, glutamic acid) lose protons.

The pI calculation applies the Henderson-Hasselbalch equation to all ionizable groups:

$$pH = pK_a + \log\left( ([A^-]) / ([HA]) \right)$$

Amino acids contribute via established pKa values: aspartic acid (3.9), glutamic acid (4.3), histidine (6.0), cysteine (8.3), tyrosine (10.9), lysine (10.5), and arginine (12.5). N-terminal amino groups (pKa ~9.6) and C-terminal carboxyl groups (pKa ~2.3) also affect charge.

Isoelectric point knowledge is essential for purification strategies. Ion-exchange chromatography, isoelectric focusing, and crystallization depend on accurate pI predictions. Proteins show minimum solubility at pI due to reduced electrostatic repulsion.

Instability index

The instability index predicts cellular protein stability through dipeptide composition analysis. It assigns Dipeptide Instability Weight Values (DIWV) to all 400 amino acid pairs based on experimental in vivo half-life data.

The calculation follows:

$$II = \frac{10}{L}\sum_{i=1}^{L-1} DIWV(x_i, x_{i+1})$$

where $L$ represents protein length, $x_i$ indicates the amino acid at position $i$, and $DIWV$ represents the instability contribution of each adjacent pair.

Values below 40 indicate stability; above 40 suggests instability. This threshold distinguishes proteins with half-lives under 5 hours (unstable) from those exceeding 16 hours (stable).

Interpretation requires considering protein localization, post-translational modifications, and cellular context, which significantly influence biological stability.

Aliphatic index

The aliphatic index quantifies relative volume of aliphatic amino acids (alanine, valine, isoleucine, leucine) as a thermostability indicator. Higher values correlate with thermal stability through enhanced hydrophobic interactions at elevated temperatures.

The calculation assigns differential weights to each aliphatic residue:

$$AI = X(Ala) + a \cdot X(Val) + b \cdot [X(Ile) + X(Leu)]$$

where $X$ represents the mole percent of each amino acid, with empirical coefficients $a = 2.9$ and $b = 3.9$ reflecting relative thermostability contributions.

Applications include thermostable enzyme engineering and studying high-temperature adaptations. Thermophilic proteins consistently show higher aliphatic indices than mesophilic counterparts, enabling optimal temperature prediction.

GRAVY score

The Grand Average of Hydropathicity (GRAVY) quantifies protein hydrophobic character by averaging Kyte-Doolittle hydropathy values across all amino acids.

$$GRAVY = \frac{1}{L} \sum_{i=1}^{L} h_i$$

where $L$ is the protein length and $h_i$ represents the hydropathy value of amino acid $i$.

GRAVY ranges from -2.0 (hydrophilic) to +2.0 (hydrophobic):

• > +1.0: Integral membrane proteins with transmembrane domains
• 0 to +1.0: Mixed regions, often peripheral membrane proteins
• < 0: Soluble hydrophilic proteins, typically cytoplasmic enzymes
• < -1.0: Highly soluble proteins in transport or signaling

GRAVY scores predict localization, membrane association, and purification behavior.

Aromaticity

Aromaticity measures the fraction of aromatic amino acids (phenylalanine, tryptophan, tyrosine) in a protein sequence. It is calculated as the sum of the mole percentages of F, W, and Y residues, following the method described by Lobry (1994).

$$\text{Aromaticity} = \frac{n(Phe) + n(Trp) + n(Tyr)}{L}$$

where $L$ is the total sequence length. Values typically range from 0.05 to 0.15 for globular proteins. High aromaticity correlates with UV absorbance intensity and contributes to hydrophobic core formation and protein–protein interaction interfaces.

Secondary structure fractions

Secondary structure fractions estimate the proportion of amino acids associated with helix, turn, and sheet conformations based on residue propensities. Unlike Chou-Fasman prediction, which applies sliding-window algorithms to predict per-residue structure, this method simply groups amino acids by their general structural tendency:

Structure	Amino acids
Helix	E, M, A, L, K
Turn	N, P, G, S, D
Sheet	V, I, Y, F, W, L, T

Each fraction equals the sum of mole percentages for that group's amino acids. Note that leucine (L) appears in both helix and sheet groups, so the three fractions do not sum to 1.0. Amino acids not listed (R, C, Q, H) do not contribute to any fraction.

These fractions provide a quick compositional indicator of structural tendency rather than a positional prediction.

Extinction coefficients

Extinction coefficients quantify 280 nm light absorption for spectrophotometric concentration determination, depending primarily on aromatic amino acids (tryptophan, tyrosine) and disulfide-bonded cysteine.

Two coefficients account for cysteine oxidation states:

• Reduced: All cysteines as free thiols
• Oxidized: Complete disulfide bond formation

The calculation employs established molar absorptivity constants:

$$\varepsilon_{280} = n(Trp) \times 5500 + n(Tyr) \times 1490 + n(Cys-Cys) \times 125$$

where $n$ represents the number of each residue type and coefficients are expressed in M⁻¹cm⁻¹.

Extinction coefficients enable Beer-Lambert law ($A = \varepsilon c l$) application for concentration determination in enzyme kinetics, interaction assays, and biochemical analysis. Coefficient choice depends on disulfide status, confirmed by DTNB assay or mass spectrometry.