ABSTRACT
Histidine is a key amino-acid residue in proteins with unique properties engendered by its imidazole side chain that can exist in three different states: two different neutral tautomeric forms and a protonated, positively charged one with a pKa value close to physiological pH. Commonly, two or all three states coexist and interchange rapidly, enabling histidine to act as both donor and acceptor of hydrogen bonds, coordinate metal ions, and engage in acid/base catalysis. Understanding the exchange dynamics among the three states is critical for assessing histidine's mechanistic role in catalysis, where the rate of proton exchange and interconversion among tautomers might be rate limiting for turnover. Here, we determine the exchange kinetics of histidine residues with pKa values representative of the accessible range from 5 to 9 by measuring pH-dependent 15N, 13C, and 1H transverse relaxation rate constants for 5 nuclei in each imidazole. Proton exchange between the imidazole and the solvent is mediated by hydronium ions at acidic and neutral pH, whereas hydroxide mediated exchange becomes the dominant mechanism at basic pH. Proton transfer is very fast and reaches the diffusion limit for pKa values near neutral pH. We identify a direct pathway between the two tautomeric forms, likely mediated by a bridging water molecule or, in the case of high pH, hydroxide ion. For histidines with pKa 7, we determine all rate constants (lifetimes) involving protonation over the entire pH range. Our approach should enable critical insights into enzymatic acid/base catalyzed reactions involving histidines in proteins.
Subject(s)
Histidine , Protons , Histidine/chemistry , Hydrogen-Ion Concentration , Kinetics , Nuclear Magnetic Resonance, BiomolecularABSTRACT
Human carbonic anhydrase II catalyzes the reversible reaction of carbon dioxide and water to form bicarbonate and a proton. His64-mediated proton shuttling between the active site and the bulk solvent is rate limiting. Here we investigate the protonation behavior of His64 as well as its structural and dynamic features in a pH dependent way. We derive two pKa values for His64, 6.25 and 7.60, that we were able to assign to its inward and outward conformation. Furthermore, we show that His64 exists in both conformations equally, independent of pH. Both conformations display an equal distribution of their two neutral tautomeric states. The life time of each conformation is short and both states display high flexibility within their orientation. Therefore, His64 is never static, but rather poised to change conformation. These findings support an energetic, dynamic and solution ensemble-based framework for the high enzymatic activity of human carbonic anhydrase II.
Subject(s)
Carbonic Anhydrase II , Protons , Humans , Bicarbonates , Carbon Dioxide , WaterABSTRACT
In the human eye lenses, the crystallin proteins facilitate transparency, light refraction, as well as UV light protection. A deregulated balanced interplay between α-, ß-, and γ-crystallin can cause cataract. γD-crystallin (hγD) is involved in the energy dissipation of absorbed UV light by energy transfer between aromatic side chains. Early UV-B induced damage of hγD with molecular resolution is studied by solution NMR and fluorescence spectroscopy. hγD modifications are restricted to Tyr 17 and Tyr 29 in the N-terminal domain, where a local unfolding of the hydrophobic core is observed. None of the tryptophan residues assisting fluorescence energy transfer is modified and hγD is remained soluble over month. Investigating isotope-labeled hγD surrounded by eye lens extracts from cataract patients reveals very week interactions of solvent-exposed side chains in the C-terminal hγD domain and some remaining photoprotective properties of the extracts. Hereditary E107A hγD found in the eye lens core of infants developing cataract shows under the here used conditions a thermodynamic stability comparable to the wild type but an increased sensitivity toward UV-B irradiation.
Subject(s)
Cataract , Lens, Crystalline , gamma-Crystallins , Humans , gamma-Crystallins/chemistry , gamma-Crystallins/metabolism , Ultraviolet Rays , Protein Folding , Lens, Crystalline/metabolism , Cataract/metabolismABSTRACT
Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.
Subject(s)
Phenylalanine , Proteins , Proteins/chemistry , Magnetic Resonance Spectroscopy/methods , Phenylalanine/chemistry , Tyrosine/chemistry , Magnetic Resonance ImagingABSTRACT
Aromatic ring flips are a hallmark of protein dynamics. They are experimentally studied by NMR spectroscopy, where recent advances have led to improved characterization across a wide range of time scales. Results on different proteins have been interpreted as continuous diffusive ring rotations or jumplike flips, leading to diverging views of the protein interior as being fluidlike or solidlike, respectively. It is challenging to distinguish between these mechanisms and other types of conformational exchange because chemical-shift-mediated line broadening provides only conclusive evidence for ring flips only if the system can be moved from the slow- to intermediate/fast-exchange regime. Moreover, whenever the chemical shift difference between the two symmetry-related sites is close to zero, it is not generally possible to determine the exchange time scale. Here we resolve these issues by measuring residual dipolar coupling (RDC)-mediated exchange contributions using NMR relaxation dispersion experiments on proteins dissolved in dilute liquid crystalline media. Excellent agreement is found between the experimental difference in RDC between the two symmetry-related sites and the value calculated from high-resolution X-ray structures, demonstrating that dynamics measured for F52 in the B1 domain of protein G reports on distinct, jumplike flips rather than other types of conformational exchange.
Subject(s)
Proteins , Magnetic Resonance Spectroscopy , Molecular Conformation , Nuclear Magnetic Resonance, Biomolecular/methods , Proteins/chemistryABSTRACT
Peptidyl-prolyl isomerases (PPIases) catalyze intrinsically slow and often rate-limiting isomerization of prolyl-peptide bonds in unfolded or partially folded proteins, thereby speeding up the folding process and preventing misfolding. They often possess binding and chaperone domains in addition to the domain carrying the isomerization activity. Although generally, their substrates display no identity in their amino acid sequence upstream and downstream of the proline with 20 possibilities for each residue, PPIases are efficient enzymes. SlyD is a highly efficient PPIase consisting of an isomerase domain and an additional chaperone domain. The binding of peptide substrates to SlyD and its enzymatic activity depend to some extend on the proline-proximal residues, however, the impact of proline-distant residues has not been investigated so far. Here, we introduce a label-free NMR-based method to measure SlyD activity on different peptide substrates and analysed the data in the context of obtained binding affinities and several co-crystal structures. We show that especially charged and aromatic residues up to eight positions downstream and three positions upstream of the proline and outside the canonical region of similar conformations affect the activity and binding, although they rarely display distinct conformations in our crystal structures. We hypothesize that these positions primarily influence the association reaction. In the absence of the chaperone domain the isomerase activity strongly correlates with substrate affinity, whereas additional factors play a role in its presence. The mutual orientation of isomerase and chaperone domains depends on the presence of substrates in both binding sites, implying allosteric regulation of enzymatic activity.
Subject(s)
Escherichia coli Proteins/chemistry , Escherichia coli Proteins/metabolism , Escherichia coli/metabolism , Peptide Fragments/metabolism , Peptidylprolyl Isomerase/chemistry , Peptidylprolyl Isomerase/metabolism , Protein Folding , Binding Sites , Catalysis , Escherichia coli/genetics , Peptide Fragments/chemistry , Protein Conformation , Substrate SpecificityABSTRACT
NMR-spectroscopy has certain unique advantages for recording unfolding transitions of proteins compared e.g. to optical methods. It enables per-residue monitoring and separate detection of the folded and unfolded state as well as possible equilibrium intermediates. This allows a detailed view on the state and cooperativity of folding of the protein of interest and the correct interpretation of subsequent experiments. Here we summarize in detail practical and theoretical aspects of such experiments. Certain pitfalls can be avoided, and meaningful simplification can be made during the analysis. Especially a good understanding of the NMR exchange regime and relaxation properties of the system of interest is beneficial. We show by a global analysis of signals of the folded and unfolded state of GB1 how accurate values of unfolding can be extracted and what limits different NMR detection and unfolding methods. E.g. commonly used exchangeable amides can lead to a systematic under determination of the thermodynamic protein stability. We give several perspectives of how to deal with more complex proteins and how the knowledge about protein stability at residue resolution helps to understand protein properties under crowding conditions, during phase separation and under high pressure.
Subject(s)
Protein Folding , Protein Unfolding , Magnetic Resonance Spectroscopy , Nuclear Magnetic Resonance, Biomolecular/methods , Protein Denaturation , Proteins/chemistry , ThermodynamicsABSTRACT
Intrinsically disordered proteins (IDPs) are an important class of proteins which lack tertiary structure elements. Their dynamic properties can depend on reversible post-translational modifications and the complex cellular milieu, which provides a crowded environment. Both influences the thermodynamic stability and folding of globular proteins as well as the conformational plasticity of IDPs. Here we investigate the intrinsically disordered C-terminal region (amino acids 613-694) of human Grb2-associated binding protein 1 (Gab1), which binds to the disease-relevant Src homolog region 2 (SH2) domain-containing protein tyrosine phosphatase SHP2 (PTPN11). This binding is mediated by phosphorylation at Tyr 627 and Tyr 659 in Gab1. We characterize induced structure in Gab1613-694 and binding to SHP2 by NMR, CD and ITC under non-crowding and crowding conditions, employing chemical and biological crowding agents and compare the results of the non-phosphorylated and tyrosine phosphorylated C-terminal Gab1 fragment. Our results show that under crowding conditions pre-structured motifs in two distinct regions of Gab1 are formed whereas phosphorylation has no impact on the dynamics and IDP character. These structured regions are identical to the binding regions towards SHP2. Therefore, biological crowders could induce some SHP2 binding capacity. Our results therefore indicate that high concentrations of macromolecules stabilize the preformed or excited binding state in the C-terminal Gab1 region and foster the binding to the SH2 tandem motif of SHP2, even in the absence of tyrosine phosphorylation.
Subject(s)
Adaptor Proteins, Signal Transducing , Intrinsically Disordered Proteins , Protein Tyrosine Phosphatase, Non-Receptor Type 11 , Adaptor Proteins, Signal Transducing/chemistry , Adaptor Proteins, Signal Transducing/genetics , Humans , Intrinsically Disordered Proteins/chemistry , Phosphorylation , Protein Binding , Protein Conformation , Protein Tyrosine Phosphatase, Non-Receptor Type 11/chemistry , Tyrosine/chemistryABSTRACT
Eps15 homology (EH) domains are universal interaction domains to establish networks of protein-protein interactions in the cell. These networks mainly coordinate cellular functions including endocytosis, actin remodeling, and other intracellular signaling pathways. They are well characterized in structural terms, except for the internal EH domain from human γ-synergin (EHγ). Here, we complete the family of EH domain structures by determining the solution structure of the EHγ domain. The structural ensemble follows the canonical EH domain fold and the identified binding site is similar to other known EH domains. But EHγ differs significantly in the N- and C-terminal regions. The N-terminal α-helix is shortened compared to known homologues, while the C-terminal one is fully formed. A significant proportion of the remaining N- and C-terminal regions are well structured, a feature not seen in other EH domains. Single mutations in both the N-terminal and the C-terminal structured extensions lead to the loss of the distinct three-dimensional fold and turn EHγ into a molten globule like state. Therefore, we propose that the structural extensions in EHγ function as a clamp and are undoubtedly required to maintain its tertiary fold.
Subject(s)
Calcium-Binding Proteins , Phosphoproteins , Amino Acid Sequence , Binding Sites , Calcium-Binding Proteins/chemistry , Humans , Intracellular Signaling Peptides and Proteins/metabolism , Phosphoproteins/chemistryABSTRACT
Aromatic side chains are attractive probes of protein dynamic, since they are often key residues in enzyme active sites and protein binding sites. Dynamic processes on microsecond to millisecond timescales can be studied by relaxation dispersion experiments that attenuate conformational exchange contributions to the transverse relaxation rate by varying the refocusing frequency of applied radio-frequency fields implemented as either CPMG pulse trains or continuous spin-lock periods. Here we present an aromatic 1H R1ρ relaxation dispersion experiment enabling studies of two to three times faster exchange processes than achievable by existing experiments for aromatic side chains. We show that site-specific isotope labeling schemes generating isolated 1H-13C spin pairs with vicinal 2H-12C moieties are necessary to avoid anomalous relaxation dispersion profiles caused by Hartmann-Hahn matching due to the 3JHH couplings and limited chemical shift differences among 1H spins in phenylalanine, tyrosine and the six-ring moiety of tryptophan. This labeling pattern is sufficient in that remote protons do not cause additional complications. We validated the approach by measuring ring-flip kinetics in the small protein GB1. The determined rate constants, kflip, agree well with previous results from 13C R1ρ relaxation dispersion experiments, and yield 1H chemical shift differences between the two sides of the ring in good agreement with values measured under slow-exchange conditions. The aromatic1H R1ρ relaxation dispersion experiment in combination with the site-selective 1H-13C/2H-12C labeling scheme enable measurement of exchange rates up to kex = 2kflip = 80,000 s-1, and serve as a useful complement to previously developed 13C-based methods.
Subject(s)
Proteins , Protons , Binding Sites , Kinetics , Nuclear Magnetic Resonance, BiomolecularABSTRACT
Proteins are dynamic entities that intermittently depart from their ground-state structures and undergo conformational transitions as a critical part of their functions. Central to understanding such transitions are the structural rearrangements along the connecting pathway, where the transition state plays a special role. Using NMR relaxation at variable temperature and pressure to measure aromatic ring flips inside a protein core, we obtain information on the structure and thermodynamics of the transition state. We show that the isothermal compressibility coefficient of the transition state is similar to that of short-chain hydrocarbon liquids, implying extensive local unfolding of the protein. Our results further indicate that the required local volume expansions of the protein can occur not only with a net positive activation volume of the protein, as expected from previous studies, but also with zero activation volume by compaction of remote void volume, when averaged over the ensemble of states.
ABSTRACT
Gliding, a type of motility based on an actin-myosin motor, is specific to apicomplexan parasites. Myosin A binds two light chains which further interact with glideosome associated proteins and assemble into the glideosome. The role of individual glideosome proteins is unclear due to the lack of structures of larger glideosome assemblies. Here, we investigate the role of essential light chains (ELCs) in Toxoplasma gondii and Plasmodium falciparum and present their crystal structures as part of trimeric sub-complexes. We show that although ELCs bind a conserved MyoA sequence, P. falciparum ELC adopts a distinct structure in the free and MyoA-bound state. We suggest that ELCs enhance MyoA performance by inducing secondary structure in MyoA and thus stiffen its lever arm. Structural and biophysical analysis reveals that calcium binding has no influence on the structure of ELCs. Our work represents a further step towards understanding the mechanism of gliding in Apicomplexa.
Subject(s)
Apicomplexa , Multiprotein Complexes/chemistry , Multiprotein Complexes/metabolism , Myosin Light Chains/chemistry , Myosin Light Chains/metabolism , Amino Acid Sequence , Apicomplexa/metabolism , Calcium/chemistry , Calcium/metabolism , Conserved Sequence , Magnetic Resonance Spectroscopy , Models, Molecular , Nonmuscle Myosin Type IIA/chemistry , Nonmuscle Myosin Type IIA/metabolism , Protein Binding , Protein Conformation , Protein Multimerization , Protein Stability , Protozoan Proteins/chemistry , Protozoan Proteins/metabolism , Structure-Activity Relationship , ThermodynamicsABSTRACT
The nucleation of Alzheimer-associated Aß peptide monomers can be catalyzed by preexisting Aß fibrils. This leads to autocatalytic amplification of aggregate mass and underlies self-replication and generation of toxic oligomers associated with several neurodegenerative diseases. However, the nature of the interactions between the monomeric species and the fibrils during this key process, and indeed the ultrastructural localization of the interaction sites have remained elusive. Here we used NMR and optical spectroscopy to identify conditions that enable the capture of transient species during the aggregation and secondary nucleation of the Aß42 peptide. Cryo-electron microscopy (cryo-EM) images show that new aggregates protrude from the entire length of the progenitor fibril. These protrusions are morphologically distinct from the well-ordered fibrils dominating at the end of the aggregation process. The data provide direct evidence that self-replication through secondary nucleation occurs along the sides of fibrils, which become heavily decorated under the current solution conditions (14 µM Aß42, 20 mM sodium phosphate, 200 µM EDTA, pH 6.8).
Subject(s)
Amyloid beta-Peptides/metabolism , Amyloid beta-Peptides/ultrastructure , Peptide Fragments/metabolism , Peptide Fragments/ultrastructure , Alzheimer Disease/pathology , Amyloid/metabolism , Amyloid/ultrastructure , Amyloid beta-Peptides/chemistry , Benzothiazoles/chemistry , Benzothiazoles/metabolism , Cryoelectron Microscopy , Humans , Hydrogen-Ion Concentration , Kinetics , Magnetic Resonance Spectroscopy , Peptide Fragments/chemistry , Time-Lapse ImagingABSTRACT
Ring flips of phenylalanine and tyrosine are a hallmark of protein dynamics. They report on transient breathing motions of proteins. In addition, flip rates also depend on stabilizing interactions in the ground state, like aromatic stacking or cation-π interaction. So far, experimental studies of ring flips have almost exclusively been performed on aromatic rings without stabilizing interactions. Here we investigate ring flip dynamics of Phe and Tyr in the aromatic cluster in GB1. We found that all four residues of the cluster, Y3, F30, Y45 and F52, display slow ring flips. Interestingly, F52, the central residue of the cluster, which makes aromatic contacts with all three others, is flipping significantly faster, while the other rings are flipping with the same rates within margin of error. Determined activation enthalpies and activation volumes of these processes are in the same range of other reported ring flips of single aromatic rings. There is no correlation of the number of aromatic stacking interactions to the activation enthalpy, and no correlation of the ring's extent of burying to the activation volume. Because of these findings, we speculate that F52 is undergoing concerted ring flips with each of the other rings.
Subject(s)
Bacterial Proteins/chemistry , Carbon Isotopes/chemistry , Nuclear Magnetic Resonance, Biomolecular , Streptococcus/chemistryABSTRACT
Intracellular aggregates of superoxide dismutase 1 (SOD1) are associated with amyotrophic lateral sclerosis. In vivo, aggregation occurs in a complex and dense molecular environment with chemically heterogeneous surfaces. To investigate how SOD1 fibril formation is affected by surfaces, we used an in vitro model system enabling us to vary the molecular features of both SOD1 and the surfaces, as well as the surface area. We compared fibril formation in hydrophilic and hydrophobic sample wells, as a function of denaturant concentration and extraneous hydrophobic surface area. In the presence of hydrophobic surfaces, SOD1 unfolding promotes fibril nucleation. By contrast, in the presence of hydrophilic surfaces, increasing denaturant concentration retards the onset of fibril formation. We conclude that the mechanism of fibril formation depends on the surrounding surfaces and that the nucleating species might correspond to different conformational states of SOD1 depending on the nature of these surfaces.
Subject(s)
Amyloid/chemistry , Biocatalysis , Hydrophobic and Hydrophilic Interactions , Protein Aggregates , Protein Unfolding , Superoxide Dismutase-1/chemistry , Superoxide Dismutase-1/metabolism , Adsorption , Apoenzymes/chemistry , Apoenzymes/metabolism , Disulfides/chemistry , Surface PropertiesABSTRACT
Post-translational methylation of lysine side chains is of great importance for protein regulation, including epigenetic control. Here, we present specific 13CHD2 labeling of dimethylated lysines as a sensitive probe of the structure, interactions, and dynamics of these groups, and outline a theoretical and experimental framework for analyzing their conformational dynamics using 1H and 13C CPMG relaxation dispersion experiments. Dimethylated lysine side chains in calcium-loaded calmodulin show a marked pH dependence of their Carr-Purcell-Meiboom-Gill (CPMG) dispersion profiles, indicating complex exchange behavior. Combined analysis of 1H and 13C CPMG relaxation dispersions requires consideration of 12-state correlated exchange of the two methyl groups due to circular three-state rotamer jumps around the Cε-Nζ axis combined with proton exchange and amine inversion. Taking into account a number of fundamental constraints, the exchange model can be reduced to include only three fitted parameters, namely, the geometric average of the rotamer-jump rate constants, the rate constant of deprotonation of Nζ, and the chemical shift difference between the trans and gauge positions of the 13C or 1H nuclei. The pH dependence indicates that protonation of the end group dramatically slows down rotamer exchange for some lysine residues, whereas deprotonation leads to rapid amine inversion coupled with rotamer scrambling. The observed variation among residues in their exchange behavior appears to depend on the structural environment of the side chain. Understanding this type of exchange process is critical to correctly interpreting NMR spectra of methylated lysine side chains. The exchange model presented here forms the basis for studying the structure and dynamics of epigenetically modified lysine side chains and perturbations caused by changes in pH or interactions with target proteins.
Subject(s)
Amines/chemistry , Carbon-13 Magnetic Resonance Spectroscopy , Lysine/chemistry , Proteins/chemistry , Proton Magnetic Resonance Spectroscopy , Hydrogen-Ion Concentration , Methylation , Protons , Solvents/chemistryABSTRACT
Aromatic side chains are often key residues in enzyme active sites and protein binding sites, making them attractive probes of protein dynamics on the millisecond timescale. Such dynamic processes can be studied by aromatic 13C or 1H CPMG relaxation dispersion experiments. Aromatic 1H CPMG relaxation dispersion experiments in phenylalanine, tyrosine and the six-ring moiety of tryptophan, however, are affected by 3J 1H-1H couplings which are causing anomalous relaxation dispersion profiles. Here we show that this problem can be addressed by site-selective 1H/2H labeling of the aromatic side chains and that artifact-free relaxation dispersion profiles can be acquired. The method has been further validated by measuring folding-unfolding kinetics of the small protein GB1. The determined rate constants and populations agree well with previous results from 13C CPMG relaxation dispersion experiments. Furthermore, the CPMG-derived chemical shift differences between the folded and unfolded states are in excellent agreement with those obtained directly from the spectra. In summary, site-selective 1H/2H labeling enables artifact-free aromatic 1H CPMG relaxation dispersion experiments in phenylalanine and the six-ring moiety of tryptophan, thereby extending the available methods for studying millisecond dynamics in aromatic protein side chains.
Subject(s)
Amino Acids, Aromatic/chemistry , Nuclear Magnetic Resonance, Biomolecular/methods , Protein Folding , Proteins/chemistry , Staining and Labeling/methods , Deuterium , Kinetics , Molecular Dynamics Simulation , Protein Conformation , Staining and Labeling/standardsABSTRACT
The deposition of α-synuclein fibrils is one hallmark of Parkinson's disease. Here, we investigate how ganglioside lipids, present in high amounts in neurons and exosomes, influence the aggregation kinetics of α-synuclein. Gangliosides, as well as, other anionic lipid species with small or large headgroups were found to induce conformational changes of α-synuclein monomers and catalyse their aggregation at mildly acidic conditions. Although the extent of this catalytic effect was slightly higher for gangliosides, the results imply that charge interactions are more important than headgroup chemistry in triggering aggregation. In support of this idea, uncharged lipids with large headgroups were not found to induce any conformational change and only weakly catalyse aggregation. Intriguingly, aggregation was also triggered by free ganglioside headgroups, while these caused no conformational change of α-synuclein monomers. Our data reveal that partially folded α-synuclein helical intermediates are not required species in triggering of α-synuclein aggregation.
ABSTRACT
Aromatic side chains in proteins are often directly evolved in stabilizing the hydrophobic core, protein binding, or enzymatic activity. They are also responsible for specific local dynamic processes, such as histidine tautomerization or ring flips. Despite their importance, they are often not targeted directly by NMR spectroscopy, because of spectroscopic complications and challenges. This chapter addresses state-of-the-art site-selective 13C-labeling methods for aromatic side chains, and describes how they solve several of the spectroscopic issues. A special emphasis is put on thereby enabled protein dynamics experiments of aromatic side chains.
Subject(s)
Amino Acids, Aromatic/chemistry , Carbon Isotopes/chemistry , GTP-Binding Protein alpha Subunits/chemistry , Glucose/chemistry , Isotope Labeling/methods , Magnetic Resonance Spectroscopy/methods , Staining and Labeling/methods , Deuterium/chemistry , Humans , Hydrophobic and Hydrophilic Interactions , Magnetic Resonance Spectroscopy/instrumentation , Molecular Dynamics Simulation , Recombinant Proteins/chemistry , ThermodynamicsABSTRACT
Electrostatic interactions significantly contribute to the stability and function of proteins. The stabilizing or destabilizing effect of local charge is reflected in the perturbation of the pKa value of an ionizable group from the intrinsic pKa value. Herein, the charge network of a hyperstable dimeric protein (ribbon-helix-helix (rhh) protein from plasmid pRN1 from Sulfolobus islandicus) is studied through experimental determination of the pKa values of all ionizable groups. Transitions were monitored by multiple NMR signals per ionizable group between pHâ 0 and 12.5, prior to a global analysis, which accounted for the effects of neighboring residues. It is found that for several residues involved in salt bridges (four Asp and one Lys) the pKa values are shifted in favor of the charged state. Furthermore, the pKa values of residues C40 and Y47, both located in the hydrophobic dimer interface, are shifted beyond 13.7. The necessary energy for such a shift is about two-thirds of the total stability of the protein, which confirms the importance of the hydrophobic core to the overall stability of the rhh protein.