Experimental pKa Value Determination of All Ionizable Groups of a Hyperstable Protein

Proton Occupancies in Histidine Side Chains of Carbonic Anhydrase II by Neutron Crystallography and NMR - Differences, Similarities and Opportunities

Histidine is a key amino-acid residues in proteins that can exist in three different protonation states: two different neutral tautomeric forms and a protonated, positively charged one. It can act as both donor and acceptor of hydrogen bonds, coordinate metal ions, and engage in acid/base catalysis. Human Carbonic Anhydrase II (HCA II) is a pivotal enzyme catalyzing the reversible hydration of carbon dioxide. It contains 12 histidine residues: six are surface exposed, two buried, three are active site zinc ion ligands, and one is a proton shuttle. Comparing results from NMR spectroscopy with previously determined neutron protein crystal structures enabled a side-by-side investigation of the proton occupancies and preferred tautomeric states of the histidine residues in HCA II. Buried and zinc coordinating histidines remain in one neutral tautomeric state across the entire pH range studied, as validated by both methods. In contrast, solvent-exposed histidines display high variability in proton occupancies. While the data were overall remarkably consistent between methods, some discrepancies were observed, shedding light on the limitations of each technique. Therefore, combining these methods with full awareness of the advantages and drawbacks of each, provides insights into the dynamic protonation landscape of HCA II histidines, crucial for elucidating enzyme catalytic mechanisms.

Discovery of alkaline laccases from basidiomycete fungi through machine learning-based approach

BACKGROUND

Laccases can oxidize a broad spectrum of substrates, offering promising applications in various sectors, such as bioremediation, biomass fractionation in future biorefineries, and synthesis of biochemicals and biopolymers. However, laccase discovery and optimization with a desirable pH optimum remains a challenge due to the labor-intensive and time-consuming nature of the traditional laboratory methods.

RESULTS

This study presents a machine learning (ML)-integrated approach for predicting pH optima of basidiomycete fungal laccases, utilizing a small, curated dataset against a vast metagenomic data. Comparative computational analyses unveiled the structural and pH-dependent solubility differences between acidic and neutral-alkaline laccases, helping us understand the molecular bases of enzyme pH optimum. The pH profiling of the two ML-predicted alkaline laccase candidates from the basidiomycete fungus Lepista nuda further validated our computational approach, showing the accuracy of this comprehensive method.

CONCLUSIONS

This study uncovers the efficacy of ML in the prediction of enzyme pH optimum from minimal datasets, marking a significant step towards harnessing computational tools for systematic screening of enzymes for biotechnology applications.

Proton Transfer Kinetics in Histidine Side Chains Determined by pH-Dependent Multi-Nuclear NMR Relaxation

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.

The Prenucleation Equilibrium of the Parathyroid Hormone Determines the Critical Aggregation Concentration and Amyloid Fibril Nucleation

Nucleation and growth of amyloid fibrils were found to only occur in supersaturated solutions above a critical concentration (ccrit ). The biophysical meaning of ccrit remained mostly obscure, since typical low values of ccrit in the sub-μM range hamper investigations of potential oligomeric states and their structure. Here, we investigate the parathyroid hormone PTH84 as an example of a functional amyloid fibril forming peptide with a comparably high ccrit of 67±21 μM. We describe a complex concentration dependent prenucleation ensemble of oligomers of different sizes and secondary structure compositions and highlight the occurrence of a trimer and tetramer at ccrit as possible precursors for primary fibril nucleation. Furthermore, the soluble state found in equilibrium with fibrils adopts to the prenucleation state present at ccrit . Our study sheds light onto early events of amyloid formation directly related to the critical concentration and underlines oligomer formation as a key feature of fibril nucleation. Our results contribute to a deeper understanding of the determinants of supersaturated peptide solutions. In the current study we present a biophysical approach to investigate ccrit of amyloid fibril formation of PTH84 in terms of secondary structure, cluster size and residue resolved intermolecular interactions during oligomer formation. Throughout the investigated range of concentrations (1 μM to 500 μM) we found different states of oligomerization with varying ability to contribute to primary fibril nucleation and with a concentration dependent equilibrium. In this context, we identified the previously described ccrit of PTH84 to mark a minimum concentration for the formation of homo-trimers/tetramers. These investigations allowed us to characterize molecular interactions of various oligomeric states that are further converted into elongation competent fibril nuclei during the lag phase of a functional amyloid forming peptide.

Energetics and dynamics of the proton shuttle of carbonic anhydrase II

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.

PKAD-2: New entries and expansion of functionalities of the database of experimentally measured pKa's of proteins

Almost all biological reactions are pH dependent and understanding the origin of pH dependence requires knowledge of the pKa's of ionizable groups. Here we report a new edition of PKAD, the PKAD-2, which is a database of experimentally measured pKa's of proteins, both wild type and mutant proteins. The new additions include 117 wild type and 54 mutant pKa values, resulting in total 1742 experimentally measured pKa's. The new edition of PKAD-2 includes 8 new wild type and 12 new mutant proteins, resulting in total of 220 proteins. This new edition incorporates a visual 3D image of the highlighted residue of interest within the corresponding protein or protein complex. Hydrogen bonds were identified, counted, and implemented as a search feature. Other new search features include the number of neighboring residues <4A from the heaviest atom of the side chain of a given amino acid. Here, we present PKAD-2 with the intention to continuously incorporate novel features and current data with the goal to be used as benchmark for computational methods.

Elucidating the Protonation State of the γ-Secretase Catalytic Dyad

γ-Secretase (GS) is an intramembrane aspartyl protease that participates in the sequential cleavage of C99 to generate different isoforms of the amyloid-β (Aβ) peptides that are associated with the development of Alzheimer's disease. Due to its importance in the proteolytic processing of C99 by GS, we performed pH replica exchange molecular dynamics (pH-REMD) simulations of GS in its apo and substrate-bound forms to sample the protonation states of the catalytic dyad. We found that the catalytic dyad is deprotonated at physiological pH in our apo form, but the presence of the substrate at the active site displaces its monoprotonated state toward physiological pH. Our results show that Asp257 acts as the general base and Asp385 as the general acid during the cleavage mechanism. We identified different amino acids such as Lys265, Arg269, and the PAL motif interacting with the catalytic dyad and promoting changes in its acid-base behavior. Finally, we also found a significant pK_a shift of Glu280 related to the internalization of TM6-CT in the GS-apo form. Our study provides critical mechanistic insight into the GS mechanism and the basis for future research on the genesis of Aβ peptides and the development of Alzheimer's disease.

Monitoring protein unfolding transitions by NMR-spectroscopy

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.

Heme Protein Binding of Sulfonamide Compounds: A Correlation Study by Spectroscopic, Calorimetric, and Computational Methods

Protein-ligand interaction studies are useful to determine the molecular mechanism of the binding phenomenon, leading to the establishment of the structure-function relationship. Here, we report the binding of well-known antibiotic sulfonamide drugs (sulfamethazine, SMZ; and sulfadiazine, SDZ) with heme protein myoglobin (Mb) using spectroscopic, calorimetric, ζ potential, and computational methods. Formation of a 1:1 complex between the ligand and Mb through well-defined equilibrium was observed. The binding constants obtained between Mb and SMZ/SDZ drugs were on the order of 104 M-1. SMZ with two additional methyl (-CH3) substitutions has higher affinity than SDZ. Upon drug binding, a notable loss in the helicity (via circular dichroism) and perturbation of the three-dimensional (3D) protein structure (via infrared and synchronous fluorescence experiments) were observed. The binding also indicated the dominance of non-polyelectrolytic forces between the amino acid residues of the protein and the drugs. The ligand-protein binding distance signified high probability of energy transfer between them. Destabilization of the protein structure upon binding was evident from differential scanning calorimetry results and ζ potential analyses. Molecular docking presented the best probable binding sites of the drugs inside protein pockets. Thus, the present study explores the potential binding characteristics of two sulfonamide drugs (with different substitutions) with myoglobin, correlating the structural and energetic aspects.

Signature of Glycylglutamic Acid Structure

BACKGROUND

Glutamate (Glu) is of great interest in biomedical research. It is considered a biomarker in diabetes, which may potentially contribute to the development of autism in genetically vulnerable human populations, and it is found in relation to advanced glycation end products (AGEs) [1]. Additionally, Glu plays an active role in the function of ligand-gated ion channel glutamate receptors, chloride channels capable of filtering glutamate, as well as Potassium (K+)-channel [2]. Glu attains α [3] and β [4] crystal forms and Cβ-CH2 show asymmetric 1H signal pattern in NMR spectra.

OBJECTIVES

The current study was undertaken to understand the signal patterns of Cβ-CH2 in Glu of the smallest dipeptide, Glycylglutamic Acid (GlyGlu), as well as the order, and planarity of the amide bond in the molecule.

MATERIALS AND METHODS

NMR spectra of GlyGlu were measured in D2O to deduce 1H and 13C chemical shifts and coupling constants. GlyGlu was crystallized from MeOH and the structure was determined by single crystal X-ray diffraction techniques.

RESULTS

The sidechain of Glu in the dipeptide dissimilates the β form. The amino group of Gly (Glycine) is protonated and exhibits hydrogen bonding with the main chain carboxylate group of a symmetry-related Glu that is deprotonated in the crystal packing of GlyGlu. The deprotonated main chain carboxylate of Glu is also in hydrogen-bonding distance from the side chain carboxylic acid group that is in the protonated form of a symmetry-related Glu of the dipeptide. The Cβ-CH2 geminal protons on the side chain of Glu have different chemical shifts and splitting pattern in 1H NMR reflecting their dissymmetric environment.

CONCLUSION

The results reported will be useful for monitoring changes that Glu and/or molecules in connection to Glu may undergo in in vivo, in situ, and in vitro conditions. This provides a valuable metric which will enable the examination of the metabolites relevant to the detection and diagnosis of disease or developmental conditions, as well as scrutinizing the effectiveness of treatment options.

Modeling and Experimental Validation of Algorithms for Maximum Quantity of Protein to be Immobilized on Solid Supports by Electrostatic Adsorption in the Strategy of Rational Design of Immobilized Derivatives

Protein immobilization by electrostatic adsorption to a support could represent a good option. On the other hand, lysozyme (EC 3.2.1.17) is a little and basic protein. The objective of this work was to test the functionality of the strategy of Rational Design of Immobilized Derivatives for the immobilization by electrostatic adsorption of egg white lysozyme on SP-Sepharose FastFlow support. The RDID_1.0 software was used to predict the superficial lysozyme clusters, the electrostatic configuration probability for each cluster, and the theoretical and estimated maximum quantity of protein to be immobilized. In addition, immobilization was performed and the experimental parameter practical maximum quantity of protein to be immobilized and the enzymatic activity of the immobilized derivative were assessed. The estimated maximum quantity of protein to be immobilized (9.49 protein mg/support g) was close to the experimental practical maximum quantity of protein to be immobilized (14.73 ± 0.09 protein mg/support g). The enzymatic activity assay with the chitosan substrate showed the catalytic functionality of the lysozyme-SP-Sepharose immobilized derivative (35.85 ± 3.07 U/support g), which preserved 78% functional activity. The used algorithm to calculate the estimated maximum quantity of protein to be immobilized works for other proteins, porous solid supports and immobilization methods, and this parameter has a high predictive value, useful for obtaining optimum immobilized derivatives. The applied methodology is valid to predict the most probable protein-support configurations and their catalytic competences, which concur with the experimental results. The produced biocatalyst had a high retention of functional activity. This indicates its functionality in enzymatic bioconversion processes.

L. major apo-acyl carrier protein forms ordered aggregates due to an exposed phenylalanine, while phosphopantetheine inhibits aggregation in the holo-form

L. major acyl carrier protein (ACP) is a mitochondrial protein, involved in fatty acid biosynthesis. The protein is expressed as an apo-protein, and post-translationally modified at Ser 37 by a 4'-Phosphopantetheinyl transferase. Crystal structure of the apo-form of the protein at pH 5.5 suggests a four helix bundle fold, typical of ACP's. However, upon lowering the pH to 5.0, it undergoes a conformational transition from α-helix to β-sheet, and displays amyloid like properties. When left for a few days at room temperature at this pH, the protein forms fibrils, visible under Transmission electron microscopy (TEM). Using an approach combining NMR, biophysical techniques, and mutagenesis, we have identified a Phe residue present on helix II of ACP, liable for this change. Phosphopantetheinylation of LmACP, or mutation of Phe 45 to the corresponding residue in E. coli ACP (methionine), slows down the conformational change. Conversely, substitution of methionine 44 of E. coli ACP with a phenylalanine, causes enhanced ThT binding. Thus, we demonstrate the unique property of an exposed Phe in inducing, and phophopantetheine in inhibiting amyloidogenesis. Taken together, our study adds L. major acyl carrier protein to the list of ACPs that act as pH sensors.

Experimental pKa Value Determination of All Ionizable Groups of a Hyperstable Protein

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.