A Coupled Ionization-Conformational Equilibrium Is Required To Understand the Properties of Ionizable Residues in the Hydrophobic Interior of Staphylococcal Nuclease

Mechanism of Radical Initiation and Transfer in Class Id Ribonucleotide Reductase Based on Density Functional Theory

Class Id ribonucleotide reductase (RNR) is a newly discovered enzyme, which employs the dimanganese cofactor in the superoxidized state (Mn^III/Mn^IV) as the radical initiator. The dimanganese cofactor of class Id RNR in the reduced state (inactive) is clearly based on the crystal structure of the Fj-β subunit. However, the state of the dimanganese cofactor of class Id RNR in the oxidized state (active) is not known. The X-band EPR spectra have shown that the activated Fj-β subunit exists in two distinct complexes, 1 and 2. In this work, quantum mechanical/molecular mechanical calculations were carried out to study class Id RNR. First, we have determined that complex 2 contains a Mn^III-(μ-oxo)₂-Mn^IV cluster, and complex 1 contains a Mn^III-(μ-hydroxo/μ-oxo)-Mn^IV cluster. Then, based on the determined dimanganese cofactors, the mechanism of radical initiation and transfer in class Id RNR is revealed. The Mn^III-(μ-oxo)₂-Mn^IV cluster in complex 2 has not enough reduction potential to initiate radical transfer directly. Instead, it needs to be monoprotonated into Mn^III-(μ-hydroxo/μ-oxo)-Mn^IV (complex 1) before the radical transfer. The protonation state of μ-oxo can be regulated by changing the protein microenvironment, which is induced by the protein aggregation and separation of β subunits with α subunits. The radical transfer between the cluster of Mn^III-(μ-hydroxo/μ-oxo)-Mn^IV and Trp30 in the radical-transfer chain of the Fj-β subunit (Mn^III/Mn^IV ↔ His100 ↔ Asp194 ↔ Trp30 ↔ Arg99) is a water-mediated tri-proton-coupled electron transfer, which transfers proton from the ε-amino group of Lys71 to the carboxyl group of Glu97 via the water molecule Wat551 and the bridging μ-hydroxo ligand through a three-step reaction. This newly discovered proton-coupled electron-transfer mechanism in class Id RNR is different from those reported in the known Ia-Ic RNRs. The ε-amino group of Lys71, which serves as a proton donor, plays an important role in the radical transfer.

Different pKa Shifts of Internal GLU8 in Human β-Endorphin Amyloid Revealing a Coupling of Internal Ionization and Stepwise Fibril Disassembly

As a functional amyloid, human β-endorphin amyloid fibril features a β-solenoid conformation and store peptide hormones within acidic secretory granules, which would be released into the blood through fibril disassembly when the cellular milieu pH increases from acidic to neutral level on exocytosis. To gain detailed atomic mechanism of β-endorphin amyloid fibrils' pH-responsive disassembly, we conduct constant pH molecular dynamics simulations to investigate the structural and dynamical properties of β-endorphin amyloid fibrils in experiencing the environmental pH changes. Our results demonstrate a clear pK_a shift of the internal ionizable residue of GLU8, and this shift becomes even more pronounced when it is buried more deeply in the amyloid fibrils. The unusual pK_a of GLU8 reveals that its protonation state changes from the protonated state in the acidic secretory granule to the deprotonated state in the neutral pH conditions in the blood, where the deprotonation of GLU8 leads to unfavorable interactions within the hydrophobic core of the amyloid and subsequent fibril disassembly. The different pK_a shifts of GLU8 relative to its positions in the amyloid fibril indicate that the β-endorphin amyloid fibril disassembly is a stepwise process, accounting for the experimental observation that the disassembly always initiates from the outermost layer. This study reveals the critical role of the protonation state of GLU8 in amyloid fibrils' pH-responsive disassembly, and provides clear insights for understanding the structural transitions of amyloids in hormone secretion. This study also provides theoretical basis for designing pH-sensitive biological tools for specific use with precise positioning of ionizable residues into the hydrophobic interior of proteins.

Active site center redesign increases protein stability preserving catalysis in thioredoxin

The stabilization of natural proteins is a long-standing desired goal in protein engineering. Optimizing the hydrophobicity of the protein core often results in extensive stability enhancements. However, the presence of totally or partially buried catalytic charged residues, essential for protein function, has limited the applicability of this strategy. Here, focusing on the thioredoxin, we aimed to augment protein stability by removing buried charged residues in the active site without loss of catalytic activity. To this end, we performed a charged-to-hydrophobic substitution of a buried and functional group, resulting in a significant stability increase yet abolishing catalytic activity. Then, to simulate the catalytic role of the buried ionizable group, we designed a combinatorial library of variants targeting a set of seven surface residues adjacent to the active site. Notably, more than 50% of the library variants restored, to some extent, the catalytic activity. The combination of experimental study of 2% of the library with the prediction of the whole mutational space by partial least squares regression revealed that a single point mutation at the protein surface is sufficient to fully restore the catalytic activity without thermostability cost. As a result, we engineered one of the highest thermal stabilities reported for a protein with a natural occurring fold (137°C). Further, our hyperstable variant preserves the catalytic activity both in vitro and in vivo.

Ionization of D571 Is Coupled with SARS-CoV-2 Spike Up/Down Equilibrium Revealing the pH-Dependent Allosteric Mechanism of Receptor-Binding Domains

As a type I viral fusion protein, SARS-CoV-2 spike undergoes a pH-dependent switch to mediate the endosomal positioning of the receptor-binding domain to facilitate viral entry into cells and immune evasion. Gaps in our knowledge concerning the conformational transitions and key intramolecular motivations have hampered the development of effective therapeutics against the virus. To clarify the pH-sensitive elements on spike-gating the receptor-binding domain (RBD) opening and understand the details of the RBD opening transition, we performed microsecond-time scale constant pH molecular dynamics simulations in this study. We identified the deeply buried D571 with a clear pKa shift, suggesting a potential pH sensor, and showed the coupling of ionization of D571 with spike RBD-up/down equilibrium. We also computed the free-energy landscape for RBD opening and identified the crucial interactions that influence RBD dynamics. The atomic-level characterization of the pH-dependent spike activation mechanism provided herein offers new insights for a better understanding of the fundamental mechanisms of SARS-CoV-2 viral entry and infection and hence supports the discovery of novel therapeutics for COVID-19.

Protein pKa Prediction by Tree-Based Machine Learning

Protonation states of ionizable protein residues modulate many essential biological processes. For correct modeling and understanding of these processes, it is crucial to accurately determine their pKa values. Here, we present four tree-based machine learning models for protein pKa prediction. The four models, Random Forest, Extra Trees, eXtreme Gradient Boosting (XGBoost), and Light Gradient Boosting Machine (LightGBM), were trained on three experimental PDB and pKa datasets, two of which included a notable portion of internal residues. We observed similar performance among the four machine learning algorithms. The best model trained on the largest dataset performs 37% better than the widely used empirical pKa prediction tool PROPKA and 15% better than the published result from the pKa prediction method DelPhiPKa. The overall root-mean-square error (RMSE) for this model is 0.69, with surface and buried RMSE values being 0.56 and 0.78, respectively, considering six residue types (Asp, Glu, His, Lys, Cys, and Tyr), and 0.63 when considering Asp, Glu, His, and Lys only. We provide pKa predictions for proteins in human proteome from the AlphaFold Protein Structure Database and observed that 1% of Asp/Glu/Lys residues have highly shifted pKa values close to the physiological pH.

Prediction of protein pK a with representation learning

The behavior of proteins is closely related to the protonation states of the residues. Therefore, prediction and measurement of pK a are essential to understand the basic functions of proteins. In this work, we develop a new empirical scheme for protein pK a prediction that is based on deep representation learning. It combines machine learning with atomic environment vector (AEV) and learned quantum mechanical representation from ANI-2x neural network potential (J. Chem. Theory Comput. 2020, 16, 4192). The scheme requires only the coordinate information of a protein as the input and separately estimates the pK a for all five titratable amino acid types. The accuracy of the approach was analyzed with both cross-validation and an external test set of proteins. Obtained results were compared with the widely used empirical approach PROPKA. The new empirical model provides accuracy with MAEs below 0.5 for all amino acid types. It surpasses the accuracy of PROPKA and performs significantly better than the null model. Our model is also sensitive to the local conformational changes and molecular interactions.

Mechanism of DOPA radical generation and transfer in metal-free class Ie ribonucleotide reductase based on density functional theory

•

The mechanism of DOPA radical generation, transfer and regeneration is revealed.

•

The superoxide O₂^•− should be protonated to HO₂^• prior to oxidizing Tyr126 to DOPA radical.

•

The protonation of Asp88 is the prerequisite for the DOPA radical generation and radical transfer.

•

Lys213 is a key residue for the transfer of the DOPA radical.

Quantum mechanical/molecular mechanical (QM/MM) calculations were carried out to investigate the mechanisms of the generation, transfer, and regeneration of the DOPA radical for metal-free class Ie ribonucleotide reductase. The crystal structure of MfR2 (Nature, 2018, 563, 416–420) was adopted for the calculations. The QM/MM calculations have revealed several key points that are vital for understanding the mechanisms. The superoxide O₂^•− provided by the flavoprotein NrdI cannot directly oxidize the residue Tyr126 to the DOPA radical. It should be protonated to HO₂^•. The calculation results suggest that the covalent modification of Tyr126 and the DOPA radical generation can be carried out with no involvement of metal cofactors. This addresses the concerns of the articles (Nature, 2018, 563, 416–420; PNAS, 2018, 115, 10022–10027). Another concern from the articles is that how the DOPA radical is transferred from the radical trap. The DFT calculations have demonstrated that Lys213 is a key residue for the radical transfer from the DOPA radical. The ε-amino group of Lys213 is used not only as a bridge for the electron transfer but also as a proton donor. It can provide a proton to DOPA126 via a water molecule, and thus the radical transfer from DOPA126 to Trp52 is facilitated. It has also been revealed that the protonation of Asp88 is the prerequisite for the DOPA radical generation and the radical transfer in class Ie. Once the radical is quenched, it can be regenerated via the oxidations by superoxide O₂^•− and hydroperoxyl radical HO₂^•.

pH-Dependent Conformational Changes Lead to a Highly Shifted pKa for a Buried Glutamic Acid Mutant of SNase

Ionizable residues are rarely present in the hydrophobic interior of proteins, but when they are, they play important roles in biological processes such as energy transduction and enzyme catalysis. Internal ionizable residues have anomalous experimental pK_a values with respect to their pK_a in bulk water. This work investigates the atomistic cause of the highly shifted pK_a of the internal Glu23 in the artificially mutated variant V23E of Staphylococcal Nuclease (SNase) using pH replica exchange molecular dynamics (pH-REMD) simulations. The pK_a of Glu23 obtained from our calculations is 6.55, which is elevated with respect to the glutamate pK_a of 4.40 in bulk water. The calculated value is close to the experimental pK_a of 7.10. Our simulations show that the highly shifted pK_a of Glu23 is the product of a pH-dependent conformational change, which has been observed experimentally and also seen in our simulations. We carry out an analysis of this pH-dependent conformational change in response to the protonation state change of Glu23. Using a four-state thermodynamic model, we estimate the two conformation-specific pK_a values of Glu23 and describe the coupling between the conformational and ionization equilibria.

Protonation-Induced Dynamic Allostery in PDZ Domain: Evidence of Perturbation-Independent Universal Response Network

Dynamic allostery is a relatively new paradigm where certain external perturbations may lead to modulation of conformational dynamics at a distant part of a protein without significant changes in the overall structure. While most well-characterized examples of dynamic allostery involve binding with other entities like small molecules, peptides, or nucleic acids, in this work we demonstrate that chemical modifications like protonation may lead to significant dynamical allosteric response in a PDZ domain protein. Tuning the protonation states of two histidine residues (H317 and H372), we identify the allosteric pathways responsible for the dynamic response. Interestingly, the same set of residues that constitute the allosteric response network upon ligand binding seem to be responsible for protonation-induced dynamic allostery. Thus, we propose the existence of an inherent universal response network in signaling proteins, where the same set of residues can respond to varying types of external perturbations in terms of rearrangement of hydrogen-bonded network and redistribution of electrostatic interaction energies.

Predicting Reactive Cysteines with Implicit-Solvent-Based Continuous Constant pH Molecular Dynamics in Amber

Cysteines existing in the deprotonated thiolate form or having a tendency to become deprotonated are important players in enzymatic and cellular redox functions and frequently exploited in covalent drug design; however, most computational studies assume cysteines as protonated. Thus, developing an efficient tool that can make accurate and reliable predictions of cysteine protonation states is timely needed. We recently implemented a generalized Born (GB) based continuous constant pH molecular dynamics (CpHMD) method in Amber for protein pKa calculations on CPUs and GPUs. Here we benchmark the performance of GB-CpHMD for predictions of cysteine pKa's and reactivities using a data set of 24 proteins with both down- and upshifted cysteine pKa's. We found that 10 ns single-pH or 4 ns replica-exchange CpHMD titrations gave root-mean-square errors of 1.2-1.3 and correlation coefficients of 0.8-0.9 with respect to experiment. The accuracy of predicting thiolates or reactive cysteines at physiological pH with single-pH titrations is 86 or 81% with a precision of 100 or 90%, respectively. This performance well surpasses the traditional structure-based methods, particularly a widely used empirical pKa tool that gives an accuracy less than 50%. We discuss simulation convergence, dependence on starting structures, common determinants of the pKa downshifts and upshifts, and the origin of the discrepancies from the structure-based calculations. Our work suggests that CpHMD titrations can be performed on a desktop computer equipped with a single GPU card to predict cysteine protonation states for a variety of applications, from understanding biological functions to covalent drug design.

Fully Atomistic Multiscale Approach for pKa Prediction

The ionization state of titratable amino acids strongly affects proteins structure and functioning in a large number of biological processes. It is therefore essential to be able to characterize the pK_a of ionizable groups inside proteins and to understand its microscopic determinants in order to gain insights into many functional properties of proteins. A big effort has been devoted to the development of theoretical approaches for the prediction of deprotonation free energies, yet the accurate theoretical/computational calculation of pK_a values is recognized as a current challenge. A methodology based on a hybrid quantum/classical approach is here proposed for the computation of deprotonation free energies. The method is applied to calculate the pK_a of formic acid, methylammonium, and methanethiol, providing results in good agreement with the corresponding experimental estimates. The pK_a is also calculated for aspartic acid and lysine as single residues in solution and for three aspartic/glutamic acids inside a well-characterized protein: hen egg white lysozyme. While for small molecules the method is able to deal with multiple protonation states of all titratable groups, this becomes computationally very expensive for proteins. The calculated pK_a values for the single amino acids (except for the zwitterionic aspartic acid) and inside the protein display a systematic shift with respect to the experimental values that suggests that the fine balance between hydrophobic and polar interactions might be not accurately reproduced by the usual classical force-fields, thus affecting the computation of deprotonation free energies. The calculated pK_a shifts inside the protein are in good agreement with the corresponding experimental ones (within 1 pK_a unit), well reproducing the pK_a changes due to the protein environment even in the case of large pK_a shifts.

Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor

The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.

Chemical rescue of H+ delivery in proton transfer mutants of reaction center of photosynthetic bacteria

In the native and most mutant reaction centers of bacterial photosynthesis, the electron transfer is coupled to proton transfer and is rate limiting for the second reduction of Q_B^- → Q_BH₂. In the presence of divalent metal ions (e.g. Cd²⁺) or in some ("proton transfer") mutants (L210DN/M17DN or L213DN), the proton delivery to Q_B^- is made rate limiting and the properties of the proton pathway can be directly examined. We found that small weak acids and buffers in large concentrations (up to 1 M) were able to rescue the severely impaired proton transfer capability differently depending on the location of the defects: lesions at the protein surface (proton gate H126H/H128H + Cd²⁺), beneath the surface (M17DN + Cd²⁺, L210DN/M17DN) or deep inside the protein (L213DN) could be completely, partially or to very small extent recovered, respectively. Small zwitterionic acids (azide/hydrazoic acid) and buffers (tricine) proved to be highly effective rescuers consistent with their enhanced binding affinity and access to any of the proton acceptors (including Q_B^- itself) in the pathway. As a consequence, back titration of the protons at L212Glu could be observed as a pH-dependence of the rate constant of the charge recombination in the presence of azide or formate. Model calculations support the collective influence of the acid cluster on the change of the protonation states upon extension of the cluster with the bound small acid. In proton transfer mutants, the rescuing agents decreased the free energy of activation together with their enthalpic and entropic components. This is in agreement with the hypothesis that they function as protein-penetrating protonophores delivering protons into the chain and select dominating paths out of many alternate routes. We estimate that the proton delivery will be accelerated in one pathway out of 100-200 alternate pathways. The implications for design of the chemical recovery of impaired intra-protein proton transfer pathways in proton transfer mutants are discussed.

Conformational sampling and polarization of Asp26 in pKa calculations of thioredoxin

Thioredoxin is a protein that has been used as model system by various computational methods to predict the pK_a of aspartate residue Asp26 which is 3.5 units higher than a solvent exposed one (eg, Asp20). Here, we use extensive atomistic molecular dynamics simulations of two different protonation states of Asp26 in combination with conformational analysis based on RMSD clustering and principle component analysis to identify representative conformations of the protein in solution. For each conformation, the Gibbs free energy of proton transfer between Asp26 and Asp20, which is fully solvated in a loop region of the protein, is calculated with the Amber99sb force field in alchemical transformations. The varying polarization of the two residues in different molecular environments and protonation states is described by Hirshfeld-I (HI) atomic charges obtained from the averaged polarized electron density. Our results show that the Gibbs free energy of proton transfer is dependent on the protein conformation, the proper sampling of the neighboring Lys57 residue orientations and on water molecules entering the hydrophobic cavity upon deprotonating Asp26. The inclusion of the polarization of both aspartate residues in the free energy cycle by HI atomic charges corrects the results from the non-polarizable force field and reproduces the experimental ΔpK_a value of Asp26.

Engineering Allostery into Proteins

Our ability to engineer protein structure and function has grown dramatically over recent years. Perhaps the next level in protein design is to develop proteins whose function can be regulated in response to various stimuli, including ligand binding, pH changes, and light. Endeavors toward these goals have tested and expanded on our understanding of protein function and allosteric regulation. In this chapter, we provide examples from different methods for developing new allosterically regulated proteins. These methods range from whole insertion of regulatory domains into new host proteins, to covalent attachment of photoswitches to generate light-responsive proteins, and to targeted changes to specific amino acid residues, especially to residues identified to be important for relaying allosteric information across the protein framework. Many of the examples we discuss have already found practical use in medical and biotechnology applications.

Switching Protein Conformational Substates by Protonation and Mutation

Protein modules that regulate the availability and conformational status of transcription factors determine the rapidity, duration, and magnitude of cellular response to changing conditions. One such system is the single-gene product Cnu, a four-helix bundle transcription co-repressor, which acts as a molecular thermosensor regulating the expression of virulence genes in enterobacteriaceae through modulation of its native conformational ensemble. Cnu and related genes have also been implicated in pH-dependent expression of virulence genes. We hypothesize that protonation of a conserved buried histidine (H45) in Cnu promotes large electrostatic frustration, thus disturbing the H-NS, a transcription factor, binding face. Spectroscopic and calorimetric methods reveal that H45 exhibits a suppressed p K_a of ∼5.1, the protonation of which switches the conformation to an alternate native ensemble in which the fourth helix is disordered. The population redistribution can also be achieved through a mutation H45V, which does not display any switching behavior at pH values greater than 4. The Wako-Saitô-Muñoz-Eaton (WSME) statistical mechanical model predicts specific differences in the conformations and fluctuations of the fourth and first helices of Cnu determining the observed pH response. We validate these predictions through fluorescence lifetime measurements of a sole tryptophan, highlighting the presence of both native and non-native interactions in the regions adjoining the binding face of Cnu. Our combined experimental-computational study thus shows that Cnu acts both as a thermo- and pH-sensor orchestrated via a subtle but quantifiable balance between the weak packing of a structural element and protonation of a buried histidine that promotes electrostatic frustration.

Hydronium Ions Accompanying Buried Acidic Residues Lead to High Apparent Dielectric Constants in the Interior of Proteins

Internal ionizable groups are known to play important roles in protein functions. A mystery that has attracted decades of extensive experimental and theoretical studies is the apparent dielectric constants experienced by buried ionizable groups, which are much higher than values expected for protein interiors. Many interpretations have been proposed, such as water penetration, conformational relaxation, local unfolding, protein intrinsic backbone fluctuations, etc. However, these interpretations conflict with many experimental observations. The virtual mixture of multiple states (VMMS) simulation method developed in our lab provides a direct approach for studying the equilibrium of multiple chemical states and can monitor p K_a values along simulation trajectories. Through VMMS simulations of staphylococcal nuclease (SNase) variants with internal Asp or Glu residues, we discovered that cations were attracted to buried deprotonated acidic groups and the presence of the nearby cations were essential to reproduce experimentally measured p K_a values. This finding, combined with structural analysis and validation simulations, suggests that the proton released from a deprotonation process stays near the deprotonated group inside proteins, possibly in the form of a hydronium ion. The existence of a proton near a buried charge has many implications in our understanding of protein functions.