Side chain electrostatic interactions and pH‐dependent expansion of the intrinsically disordered,highly acidic carboxyl‐terminus of γ‐tubulin

Intramolecular electrostatic attraction and repulsion strongly influence the conformational sampling of intrinsically disordered proteins and domains (IDPs). In order to better understand this complex relationship, we have used nuclear magnetic resonance to measure side chain pK_a values and pH‐dependent translational diffusion coefficients for the unstructured and highly acidic carboxyl‐terminus of γ‐tubulin (γ‐CT), providing insight into how the net charge of an IDP relates to overall expansion or collapse of the conformational ensemble. Many of the pK_a values in the γ‐CT are shifted upward by 0.3–0.4 units and exhibit negatively cooperative ionization pH profiles, likely due to the large net negative charge that accumulates on the molecule as the pH is raised. pK_a shifts of this magnitude correspond to electrostatic interaction energies between the affected residues and the rest of the charged molecule that are each on the order of 1 kcal mol^?1. Diffusion of the γ‐CT slowed with increasing net charge, indicative of an expanding hydrodynamic radius (r_H). The degree of expansion agreed quantitatively with what has been seen from comparisons of IDPs with different charge content, yielding the general trend that every 0.1 increase in relative charge (|Q|/res) produces a roughly 5% increase in r_H. While γ‐CT pH titration data followed this trend nearly perfectly, there were substantially larger deviations for the database of different IDP sequences. This suggests that other aspects of an IDP's primary amino acid sequence beyond net charge influence the sensitivity of r_H to electrostatic interactions.     

Protein adaptations to temperature and pressure: complementary roles of adaptive changes in amino acid sequence and internal milieu

Retention of required structural and functional properties of proteins in species adapted to different temperatures and pressures is achieved through variation in amino acid sequence and accumulation of small organic solutes that stabilize protein traits. Conservation of ligand binding and catalytic rate can be achieved by minor differences in sequence. For orthologs of lactate dehydrogenase-A (A₄-LDH) temperature adaptation may involve only a single amino acid substitution. Adaptation involves changes in conformational mobility of regions of A₄-LDH that undergo movement during ligand binding, movements that are rate-limiting to catalysis. A model that integrates adaptations in sequence and intracellular milieu is developed on the basis of conformational microstates. Although orthologs of different thermally adapted species vary in stability, at physiological temperatures it is hypothesized that a similar ensemble of conformational microstates exists for all orthologs. Organic solutes stabilize this ensemble of microstates. Differences among orthologs in responses to organic solutes at a common temperature lead to similar responses at normal body temperatures. Because protein stability increases at high protein concentrations, intrinsic stabilities of proteins may reflect the protein concentrations of the cellular compartments in which they occur. Protein–stabilizing solutes like trimethylamine-N-oxide (TMAO) conserve protein function and structure at elevated hydrostatic pressures.     

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK_a values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK_a values measured for individual ionizable residues with NMR spectroscopy are actually population‐weighted averages of the pK_a in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK_a values of surface Glu and Asp residues are affected by pH‐sensitive fluctuations of the backbone between folded and locally unfolded conformations. ¹⁵N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps‐ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub‐global unfolding near carboxylic groups. This study has implications for structure‐based pK_a calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH‐sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co‐operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK_a values of ionizable residues in proteins. Proteins 2014; 82:3132–3143. © 2014 Wiley Periodicals, Inc.     

Calculating pKa values in the cAMP-dependent protein kinase: the effect of conformational change and ligand binding

The conformational change observed upon ligand binding and phosphorylation for the cAMP-dependent protein kinase (protein kinase A-PKA) is of high importance for the regulation of its activity. We calculate pKa values and net charges for 18 3D structures of PKA in various conformations and liganded states to examine the role of electrostatics in ligand binding and activation. We find that the conformational change of PKA takes place without any significant net proton uptake/release at all pH values, thus indicating that PKA has evolved to reduce any pH-dependent barriers to the conformational motion. We furthermore find that the binding of ligands induces large changes in the net charge of PKA at most pH values, but significantly, we find that the net charge difference at physiological pH is close to zero, thus indicating that the active-site pKa values have been preorganized for substrate binding. We are unable to unequivocally resolve the identity of the groups responsible for determining the pH-activity profile of PKA but speculate that the titration of Lys 168 or the titration of ATP itself could be responsible for the loss of activity at high pH values. Finally, we examine the effect of point mutations on the pKa values of the PKA catalytic residues and find these to be relatively insensitive to both noncharge-altering and charge-altering mutations.     

Insight into polyproline II helical bundle stability in an antifreeze protein denatured state

The use of polyproline II (PPII) helices in protein design is currently hindered by limitations in our understanding of their conformational stability and folding. Recent studies of the snow flea antifreeze protein (sfAFP), a useful model system composed of six PPII helices, suggested that a low denatured state entropy contributes to folding thermodynamics. Here, circular dichroism spectroscopy revealed minor populations of PPII like conformers at low temperature. To get atomic level information on the conformational ensemble and entropy of the reduced, denatured state of sfAFP, we have analyzed its chemical shifts and {¹H}-¹⁵N relaxation parameters by NMR spectroscopy at four experimental conditions. No significant populations of stable secondary structure were detected. The stiffening of certain N-terminal residues at neutral versus acidic pH and shifted pK_a values leads us to suggest that favorable charge-charge interactions could bias the conformational ensemble to favor the formation the C1-C28 disulfide bond during nascent folding, although no evidence for preferred contacts between these positions was detected by paramagnetic relaxation enhancement under denaturing conditions. Despite a high content of flexible glycine residues, the mobility of the sfAFP denatured ensemble is similar for denatured α/β proteins both on fast ps/ns as well as slower μs/ms timescales. These results are in line with a conformational entropy in the denatured ensemble resembling that of typical proteins and suggest that new structures based on PPII helical bundles should be amenable to protein design.     

Theoretical and computational analysis of the membrane potential generated by cytochrome c oxidase upon single electron injection into the enzyme

We have developed theory and the computational scheme for the analysis of the kinetics of the membrane potential generated by cytochrome c oxidase upon single electron injection into the enzyme. The theory allows one to connect the charge motions inside the enzyme to the membrane potential observed in the experiments by using data from the “dielectric topography” map of the enzyme that we have created. The developed theory is applied for the analysis of the potentiometric data recently reported by the Wikström group [I. Belevich, D.A. Bloch, N. Belevich, M. Wikström and M.I. Verkhovsky, Exploring the proton pump mechanism of cytochrome c oxidase in real time, Proc. Natl. Acad. Sci. U. S. A. 104 (2007) 2685-2690] on the O to E transition in Paracoccus denitrificans oxidase. Our analysis suggests, that the electron transfer to the binuclear center is coupled to a proton transfer (proton loading) to a group just “above” the binuclear center of the enzyme, from which the pumped proton is subsequently expelled by the chemical proton arriving to the binuclear center. The identity of the pump site could not be determined with certainty, but could be localized to the group of residues His326 (His291 in bovine), propionates of heme a₃, Arg 473/474, and Trp164. The analysis also suggests that the dielectric distance from the P-side to Fe a is 0.4 or larger. The difficulties and pitfalls of quantitative interpretation of potentiometric data are discussed.     

Conditionally lethal mutants of bacteriophage T4 defective in production of a transfer RNA

The two buried carboxyls (Asp-102 and Asp-194) in both chymotrypsin and chymotrypsinogen are ionized at pH values greater than 4.2 and may be ionized even as low as pH 3.This was demonstrated by coupling most of the surface carboxyis of the proteins by a carbodi-imide with glycinamide or semicarbazide to diminish the groups ionizing at low pH and then titrating the proton uptake on denaturation by sodium dodecyl sulphate between pH 3.0 and 4.6. At pH values greater than 4.2 all unblocked carboxyls are ionized. The proton uptake during the conformational change on denaturation was determined by a stopped-flow procedure and found to be about 2H⁺/mol between pH 3.0 and 3.6. The rate constant for the uptake of protons is the same as that for the exposure of tryptophan and lies in the tens of millisecond region.The buried negative charge at the active site appears to be mainly on Asp-102 rather than on His-57, the pK_a of which must be raised by the buried charge. This enhances its efficacy as a base catalyst in the “charge relay system”.The presence of an intact charge relay system in the inactive zymogen illustrates the importance of stereochemical fit between enzyme and substrate. Enzyme catalysis could hardly be mediated by a catalyst which is uniquely reactive in the absence of correct enzyme-substrate orientation as this would be inconsistent with its specificity.     

The charge-relay system of serine proteinases: proton magnetic resonance titration studies of the four histidines of porcine trypsin.

Peaks corresponding to the C₍₂₎-protons of all four histidine residues of porcine β-trypsin were resolved in 250 MHz nuclear magnetic resonance spectra after deuteration of the slowly exchangeable N-H groups (whose resonances obscure the histidine peaks) by reversible unfolding of the protein in D₂O. One of the four peaks was assigned to the charge-relay histidine in the active site of trypsin (His(57) in the bovine chymotrypsinogen numbering system). Whereas the three other histidine C₍₂₎-peaks exhibited normal titration curves with single pK′ values of 7.20, 6.71 and 6.67, the peak assigned to His(57) had an abnormal titration curve showing two protonation steps in the pH range from 1 to 9. The first protonation with a pH′_mid of 5.0 is rapid on the nuclear magnetic resonance time-scale; the second with a pH′_mid of 4.5 is slow and apparently involves conformational transitions between two states having lifetimes of approximately 18 ms.In the complex between porcine β-trypsin and bovine pancreatic trypsin inhibitor (Kunitz) His(57) was found to be insensitive to pH over the range from 4 to 9 and its chemical shift resembles that of His(57) in the singly protonated charge relay of free trypsin. This result provides direct evidence that the trypsin charge relay acts as a proton acceptor in the initial catalytic step which leads to the formation of a tetrahedral complex. In the presence of equimolar bovine pancreatic trypsin inhibitor (Kunitz) the pH'_mid of the conformational transition that affects the charge-relay histidine is lowered from 4.5 to approximately 3.5.     

Does secondary structure determine tertiary structure in proteins?

Is highly approximate knowledge of a protein's backbone structure sufficient to successfully identify its family, superfamily, and tertiary fold? To explore this question, backbone dihedral angles were extracted from the known three‐dimensional structure of 2,439 proteins and mapped into 36 labeled, 60° × 60° bins, called mesostates. Using this coarse‐grained mapping, protein conformation can be approximated by a linear sequence of mesostates. These linear strings can then be aligned and assessed by conventional sequence‐comparison methods. We report that the mesostate sequence is sufficient to recognize a protein's family, superfamily, and fold with good fidelity. Proteins 2005. © 2005 Wiley‐Liss, Inc.     

Reaction of aspirin with cysteinyl residues of lens γ-crystallins: a mechanism for the proposed anti-cataract effect of aspirin

Incubation of lens crystallins with aspirin inhibits the development of opacities caused by cyanate. Cyanate-induced opacities are thought to be due to carbamylation of the lysyl residues which causes a decrease in the protein charge and subsequent conformational changes that permit disulfide bonding. Because aspirin can also react with lysyl residues, it has been proposed that the aspiring inhibition of cataractogonesis is due to acetylation of the lysyl which would block their reaction with cyanata. However, acetylation oflysyl residues also lowers the protein charge and would be expected to effect changes in protein conformation similar to those caused by carbamylation. Therefore, acetylation of the lysyl residues is not a satisfactory explanation for the inhibitory effect of aspirin on lens opacification. Our investigations of the reactions of cyanate and aspirin with bovine γ_II-crystallins show that the cysteinyl residues are also carbamylated and acetylated at pH 7.4. At this pH, the carbamylation at the cysteinyl residues is reversible, leading to regeneration of the thiol group and disulfide bonding. In contrast, the acetylation at cysteinyl residues is stable at pH 7.4 and can prevent disulfide bonding. This difference in stability explains how cyanate promotes, and aspirin inhibits, cataractogenesis.     

Does Lack of Secondary Structure Imply Intrinsic Disorder in Proteins? A Sequence Analysis

Intrinsically disordered proteins (IDPs)/protein regions (IDPRs) lack unique three-dimensional structure at the level of secondary and/or tertiary structure and are represented as an ensemble of interchanging conformations. To investigate the role of presence/absence of secondary structures in promoting intrinsic disorder in proteins, a comparative sequence analysis of IDPs, IDPRs and proteins with minimal secondary structures (less than 5%) is required. A sequence analysis reveals proteins with minimal secondary structure content have high mean net positive charge, low mean net hydrophobicity and low sequence complexity. Interestingly, analysis of the relative local electrostatic interactions reveal that an increase in the relative repulsive interactions between amino acids separated by three or four residues lead to either loss of secondary structure or intrinsic disorder. IDPRs show increase in both local negative-negative and positive-positive repulsive interactions. While IDPs show a marked increase in the local negative-negative interactions, proteins with minimal secondary structure depict an increase in the local positive-positive interactions. IDPs and IDPRs are enriched in D, E and Q residues, while proteins with minimal secondary structure are depleted of these residues. Proteins with minimal secondary structures have higher content of G and C, while IDPs and IDPRs are depleted of these residues. These results confirm that proteins with minimal secondary structure have a distinctly different propensity for charge, hydrophobicity, specific amino acids and local electrostatic interactions as compared to IDPs/IDPRs. Thus we conclude that lack of secondary structure may be a necessary but not a sufficient condition for intrinsic disorder in proteins.     

Interplay of charge distribution and conformation in peptides: comparison of theory and experiment

We assessed the correlation between charge distribution and conformation of flexible peptides by comparing the theoretically calculated potentiometric-titration curves of two model peptides, Ac-Lys5-NHMe (a model of poly-L-lysine) and Ac-Lys-Ala11-Lys-Gly2-Tyr-NH2 (P1) in water and methanol, with the experimental curves. The calculation procedure consisted of three steps: (i) global conformational search of the peptide under study using the electrostatically driven Monte Carlo (EDMC) method with the empirical conformational energy program for peptides (ECEPP)/3 force field plus the surface-hydration (SRFOPT) or the generalized Born surface area (GBSA) solvation model as well as a molecular dynamics method with the assisted model building and energy refinement (AMBER)99/GBSA force field; (ii) reevaluation of the energy in the pH range considered by using the modified Poisson-Boltzmann approach and taking into account all possible protonation microstates of each conformation, and (iii) calculation of the average degree of protonation of the peptide at a given pH value by Boltzmann averaging over conformations. For Ac-Lys5-NHMe, the computed titration curve agrees qualitatively with the experimental curve of poly-L-lysine in 95% methanol. The experimental titration curves of peptide P1 in water and methanol indicate a remarkable downshift of the first pK(a) value compared to the values for reference compounds (n-butylamine and phenol, respectively), suggesting the presence of a hydrogen bond between the tyrosine hydroxyl oxygen and the H(epsilon) proton of a protonated lysine side chain. The theoretical titration curves agree well with the experimental curves, if conformations with such hydrogen bonds constitute a significant part of the ensemble; otherwise, the theory predicts too small a downward pH shift.     

Trans to Cis proton concentration gradients accelerate ionophore A23187-mediated net fluxes of Ca2+ across the human red cell membrane

Ionophore A23187-mediated net influx of Ca²⁺ in ATP-depleted human red cells was studied as a function of the pH and the proton concentration gradient across the membranes. Utilizing the Ca²⁺-induced increase in K⁺ conductance of the cell membranes, various CCCP-mediated proton gradients were raised across the membranes of cells suspended in unbuffered salt solutions with different K⁺ concentrations. In ionophore-mediated equilibrium the concentration ratios of ionized Ca between ATP-depleted, DIDS-treated cells and their suspension medium were equal to the concentration ratios of protons raised to the second power. With no proton concentration gradient across the membranes the net influxes of Ca²⁺ as a function of pH resembled a titration curve of a weak acid, with half maximal net influx at pH 7.3, at 100 μM extracellular Ca²⁺. With cellular pH fixed at various values, the net influx of Ca²⁺ was determined as a function of the proton concentration gradient. A linear relationship between the logarithm of net influx and the difference between extracellular and cellular pH was found at all cellular pH values tested, but the proton concentration gradient acceleration was a function of the cellular pH. Accelerations between 10- and 40- times per unit ΔpH were found and net effluxes were correspondingly decreased. The results are discussed in relation to present models of the mechanism of ionophore A23187-mediated Ca²⁺ transport. The importance of the proton concentration gradient dependency is discussed in relation to the induced oscillations in K⁺-conductance of human red cell membranes previously reported (Vestergaard-Bogind and Bennekou (1982) Biochim. Biophys. Acta 688, 37ndash;44).     

Subtle Change in the Charge Distribution of Surface Residues May Affect the Secondary Functions of Cytochrome c

Although the primary function of cytochrome c (cyt c) is electron transfer, the protein caries out an additional secondary function involving its interaction with membrane cardiolipin (CDL), its peroxidase activity, and the initiation of apoptosis. Whereas the primary function of cyt c is essentially conserved, its secondary function varies depending on the source of the protein. We report here a detailed experimental and computational study, which aims to understand, at the molecular level, the difference in the secondary functions of cyt c obtained from horse heart (mammalian) and Saccharomyces cerevisiae (yeast). The conformational landscape of cyt c has been found to be heterogeneous, consisting of an equilibrium between the compact and extended conformers as well as the oligomeric species. Because the determination of relative populations of these conformers is difficult to obtain by ensemble measurements, we used fluorescence correlation spectroscopy (FCS), a method that offers single-molecule resolution. The population of different species is found to depend on multiple factors, including the protein source, the presence of CDL and urea, and their concentrations. The complex interplay between the conformational distribution and oligomerization plays a crucial role in the variation of the pre-apoptotic regulation of cyt c observed from different sources. Finally, computational studies reveal that the variation in the charge distribution at the surface and the charge reversal sites may be the key determinant of the conformational stability of cyt c.     

Internal charge transfer in cytochrome c oxidase at a limited proton supply: Proton pumping ceases at high pH

Background

In the membrane-bound enzyme cytochrome c oxidase, electron transfer from cytochrome c to O₂ is linked to proton uptake from solution to form H₂O, resulting in a charge separation across the membrane. In addition, the reaction drives pumping of protons across the membrane.

Methods

In this study we have measured voltage changes as a function of pH during reaction of the four-electron reduced cytochrome c oxidase from Rhodobacter sphaeroides with O₂. These electrogenic events were measured across membranes containing purified enzyme reconstituted into lipid vesicles.

Results

The results show that the pH dependence of voltage changes (primarily associated with proton transfer) during O₂ reduction does not match that of the previously studied absorbance changes (primarily associated with electron transfer). Furthermore, the voltage changes decrease with increasing pH.

Conclusions

The data indicate that cytochrome c oxidase does not pump protons at high pH (10.5) (or protons are taken from the “wrong” side of the membrane) and that at this pH the net proton-uptake stoichiometry is ∼ 1/2 of that at pH 8. Furthermore, the results provide a basis for interpretation of results from studies of mutant forms of the enzyme.

General significance

These results provide new insights into the function of cytochrome c oxidase.     

Conformational response to charge clustering in synthetic intrinsically disordered proteins

Background

Recent theoretical and computational studies have shown that the charge content and, most importantly, the linear distribution of opposite charges are major determinants of conformational properties of intrinsically disordered proteins (IDPs). Charge segregation in a sequence can be measured through κ, which represents a normalized measure of charge asymmetry. A strong inverse correlation between κ and radius of gyration has been previously demonstrated for two independent sets of permutated IDP sequences.

Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways

The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.     

The Hexose-Proton Cotransport System of Chlorella : pH-Dependent Change in Km Values and Translocation Constants of the Uptake System

The proton concentration in the medium affects the maximal velocity of sugar uptake with a K_m of 0.3 mM (high affinity uptake). By decreasing the proton concentration a decrease in high affinity sugar uptake is observed, in parallel the activity of a low affinity uptake system (K_m of 50 mM) rises. Both systems add up to 100%. The existence of the carrier in two conformational states (protonated and unprotonated) has been proposed therefore, the protonated form with high affinity to 6-deoxyglucose, the unprotonated form with low affinity. A plot of extrapolated V_max values at low substrate concentration versus proton concentration results in a K_m for protons of 0.14 µM, i.e. half-maximal protonation of the carrier is achieved at pH 6.85. The stoichiometry of protons cotransported per 6-deoxyglucose is close to 1 at pH 6.0–6.5. At higher pH values the stoichiometry continuously decreases; at pH 8.0 only one proton is cotransported per four molecules of sugar. Whereas the translocation of the protonated carrier is strictly dependent on sugar this coupling is less strict for the unprotonated form. Therefore at alkaline pH a considerable net efflux of accumulated sugar can occur. The dependence of sugar accumulation on pH has been measured. The decrease in accumulation with higher pH values can quantitatively be explained by the decrease in the amount of protonated carrier. The properties of the unprotonated carrier resemble strikingly the properties of carrier at the inner side of the membrane. The inside pH of Chlorella was measured with the weak acid 5,5-dimethyl-2, 4-oxazolidinedion (DMO). At an outside pH of 6.5 the internal pH was found to be 7.2. To explain the extent of sugar accumulation it has to be assumed that the membrane potential also contributes to active sugar transport in this alga.     

Slow Proton Transfer Coupled to Unfolding Explains the Puzzling Results of Single-Molecule Experiments on BBL,a Paradigmatic Downhill Folding Protein

A battery of thermodynamic, kinetic, and structural approaches has indicated that the small α-helical protein BBL folds-unfolds via the one-state downhill scenario. Yet, single-molecule fluorescence spectroscopy offers a more conflicting view. Single-molecule experiments at pH 6 show a unique half-unfolded conformational ensemble at mid denaturation, whereas other experiments performed at higher pH show a bimodal distribution, as expected for two-state folding. Here we use thermodynamic and laser T-jump kinetic experiments combined with theoretical modeling to investigate the pH dependence of BBL stability, folding kinetics and mechanism within the pH 6–11 range. We find that BBL unfolding is tightly coupled to the protonation of one of its residues with an apparent pK_a of ∼7. Therefore, in chemical denaturation experiments around neutral pH BBL unfolds gradually, and also converts in binary fashion to the protonated species. Moreover, under the single-molecule experimental conditions (denaturant midpoint and 279 K), we observe that proton transfer is much slower than the ∼15 microseconds folding-unfolding kinetics of BBL. The relaxation kinetics is distinctly biphasic, and the overall relaxation time (i.e. 0.2–0.5 ms) becomes controlled by the proton transfer step. We then show that a simple theoretical model of protein folding coupled to proton transfer explains quantitatively all these results as well as the two sets of single-molecule experiments, including their more puzzling features. Interestingly, this analysis suggests that BBL unfolds following a one-state downhill folding mechanism at all conditions. Accordingly, the source of the bimodal distributions observed during denaturation at pH 7–8 is the splitting of the unique conformational ensemble of BBL onto two slowly inter-converting protonation species. Both, the unprotonated and protonated species unfold gradually (one-state downhill), but they exhibit different degree of unfolding at any given condition because the native structure is less stable for the protonated form.     

Protein charge ladders reveal that the net charge of ALS-linked superoxide dismutase can be different in sign and magnitude from predicted values

This article utilized “protein charge ladders”—chemical derivatives of proteins with similar structure, but systematically altered net charge—to quantify how missense mutations that cause amyotrophic lateral sclerosis (ALS) affect the net negative charge (Z) of superoxide dismutase-1 (SOD1) as a function of subcellular pH and Zn²⁺ stoichiometry. Capillary electrophoresis revealed that the net charge of ALS-variant SOD1 can be different in sign and in magnitude—by up to 7.4 units per dimer at lysosomal pH—than values predicted from standard pK_a values of amino acids and formal oxidation states of metal ions. At pH 7.4, the G85R, D90A, and G93R substitutions diminished the net negative charge of dimeric SOD1 by up to +2.29 units more than predicted; E100K lowered net charge by less than predicted. The binding of a single Zn²⁺ to mutant SOD1 lowered its net charge by an additional +2.33 ± 0.01 to +3.18 ± 0.02 units, however, each protein regulated net charge when binding a second, third, or fourth Zn²⁺ (ΔZ < 0.44 ± 0.07 per additional Zn²⁺). Both metalated and apo-SOD1 regulated net charge across subcellular pH, without inverting from negative to positive at the theoretical pI. Differential scanning calorimetry, hydrogen-deuterium exchange, and inductively coupled plasma mass spectrometry confirmed that the structure, stability, and metal content of mutant proteins were not significantly affected by lysine acetylation. Measured values of net charge should be used when correlating the biophysical properties of a specific ALS-variant SOD1 protein with its observed aggregation propensity or clinical phenotype.