Change In Protonation State Research Articles

The Buffer Capacity of Polyelectrolyte Microcapsules Depends on the Type of Template.

One of the key physicochemical parameters of polyelectrolyte microcapsules (PMCs) is their buffer capacity (BC). The BC of the microcapsules allows for an assessment of the change in protonation state across the entire polyelectrolyte system, which directly impacts the buffer barrier of PMCs, as well as the stability and physical properties of their shell. However, the buffer capacity of PMCs and their behavior under changes in ionic strength and temperature can differ depending on the type of core used to form the microcapsules. As part of this study, we revealed the buffer capacity (BC) of polyelectrolyte microcapsules formed on polystyrene cores (PMCPs) and studied the influence of ionic strength and environmental temperature on the BC of these capsules. We found that the buffer capacity of PMCPs differs from the BC of water at a pH above 8; the addition of sodium chloride leads to an increase in buffer capacity in alkaline conditions, and conversely, thermal treatment leads to its decrease at a pH of 9. The results obtained are different from the BC of polyelectrolyte microcapsules formed on CaCO3 cores, which suggests a difference in the physicochemical properties of these types of capsules. The buffer capacity of polyelectrolyte microcapsules depends on the type of template used.

Probing bilin-protein interaction in the protochromic photocycle of cyanobacteriochrome RcaE by site-directed mutagenesis.

Cyanobacteriochromes (CBCRs) are members of the phytochrome superfamily of photosensor proteins that bind a bilin chromophore. CBCRs exhibit substantial diversity in their absorption wavelengths through a variety of bilin-protein interactions. RcaE is the first discovered cyanobacteriochrome as a regulator of chromatic acclimation, where cyanobacteria optimize the absorption wavelength of their photosynthetic antenna. RcaE undergoes a reversible photoconversion between a green-absorbing (Pg) and a red-absorbing (Pr) states, where the bilin chromophore adopts a deprotonated C15-Z,anti and a protonated C15-E,syn structures, respectively. This photocycle is designated as "protochromic photocycle" as the change of the bilin protonation state is responsible for the large absorption shift. With the guidance of recently determined Pg and Pr structures of RcaE, in this study, we investigated bilin-chromophore interaction by site-directed mutagenesis on three key residues referred to as a protochromic triad and also other conserved residues interacting with the bilin. Among the protochromic triad residues, Glu217 and Lys261 are critical for the formation of the Pr state, while Leu249 is critical for the formation of both Pg and Pr states. Substitution in other conserved residues, including Val218, Phe219, and Pro220 in the wind-up helix and Phe252, Phe214, and Leu209 in a part of the bilin-binding pocket, had less substantial effects on the spectral sensitivity in RcaE. These data provide insights into our understanding of the bilin-chromophore interaction in the protochromic photocycle and also its evolution in the CBCRs.

Elucidating the Activation Mechanism of the Proton-sensing GPR68 Receptor

GPR68 is a proton-sensing G-protein Coupled Receptor (GPCR) involved in a variety of physiological processes and disorders including neoplastic pathologies. While GPR68 and few other GPCRs have been shown to be activated by a decrease in the extracellular pH, the molecular mechanism of their activation remains largely unknown. In this work, we used a combined computational and in vitro approach to provide new insight into the activation mechanism of the receptor. Molecular Dynamics simulations of GPR68 were used to model the changes in residue interactions and motions triggered by pH. Global and local rearrangements consistent with partial activation were observed upon protonation of the inactive state. Selected extracellular histidine and transmembrane acidic residues were found to have significantly upshifted pKa values during the simulations, consistently with their previously hypothesised role in activation through changes in protonation state. Moreover, a novel pairing between histidine and acidic residues in the extracellular region was highlighted by both sequence analyses and simulation data and tested through site-directed mutagenesis. At last, we identified a previously unknown hydrophobic lock in the extracellular region that might stabilise the inactive conformation and regulate the transition to the active state.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures.

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Proposed modification to a muscle-like acid-base switchable [2](2)rotaxane for improved force delivery

ABSTRACTSwitchable rotaxanes are good candidates for molecular machines because of their ability to undergo reversible and repeatable co-conformational isomerism in response to an external stimulus. An acid-base switchable [2](2)rotaxane reported by Li et al. [J Org Chem. 79(15): [6996–7004] has been reported to exhibit a 48% change in length upon a stimulus-induced change in protonation state, suggesting that it could serve as the fundamental element of a ‘molecular muscle’. In the present work, the mechanical behaviour of this topologically complex system is explored with calculations based on semiempirical electronic structure theory. Our calculations predict that the original system will not work as envisioned. Although the calculations confirm a 48% change in length, they also reveal an energy barrier to ring shuttling upon deprotonation, preventing the immediate delivery of force. Therefore, a modification to the system is proposed. Similar semiempirical calculations were carried out on the modified system and the force calculated. The average force generated from the modified system is predicted to be approximately −86 pN. The results demonstrate that the proposed modification reduces the force generated by ring shuttling by 11% but largely eliminates the energy barrier to ring shuttling so that the force is immediately available.

Targeting Ultrafast Spectroscopic Insights into Red Fluorescent Proteins.

Red fluorescent proteins (RFPs) represent an increasingly popular class of genetically encodable bioprobes and biomarkers that can advance next-generation breakthroughs across the imaging and life sciences. Since the rational design of RFPs with improved functions or enhanced versatility requires a mechanistic understanding of their working mechanisms, while fluorescence is intrinsically an ultrafast event, a suitable toolset involving steady-state and time-resolved spectroscopic techniques has become powerful in delineating key structural features and dynamic steps which govern irreversible photoconverting or reversible photoswitching RFPs, and large Stokes shift (LSS)RFPs. The pertinent cis-trans isomerization and protonation state change of RFP chromophores in their local environments, involving key residues in protein matrices, lead to rich and complicated spectral features across multiple timescales. In particular, ultrafast excited-state proton transfer in various LSSRFPs showcases the resolving power of wavelength-tunable femtosecond stimulated Raman spectroscopy (FSRS) in mapping a photocycle with crucial knowledge about the red-emitting species. Moreover, recent progress in noncanonical RFPs with a site-specifically modified chromophore provides an appealing route for efficient engineering of redder and brighter RFPs, highly desirable for bioimaging. Such an effective feedback loop involving physical chemists, protein engineers, and biomedical microscopists will enable future successes to expand fundamental knowledge and improve human health.

Dynamic Protonation States Underlie Carbene Formation in ThDP-Dependent Enzymes: A Theoretical Study.

The activation mechanism of thiamine diphosphate (ThDP) in enzymes has long been the subject of intense research and controversial discussion. Particularly contentious is the formation of a carbene intermediate, the first one observed in an enzyme. For the formation of the carbene to take place, both intramolecular and intermolecular proton transfer pathways have been proposed. However, the physiologically relevant pH of ThDP-dependent enzymes around neutrality does not seem to be suitable for the formation of such reactive chemical species. Herein, we investigate the general mechanism of activation of the ThDP cofactor in human transketolase (TKT), by means of electronic structure methods. We show that in the case of the human TKT, the carbene species is accessible through a pKa shift induced by the electrostatics of a neighboring histidine residue (H110), whose protonation state change modulates the pKa of ThDP and suppresses the latter by more than 6 pH units. Our findings highlight that ThDP enzymes activate the cofactor beyond simple geometric constraints and the canonical glutamate. Such observations in nature can pave the way for the design of biomimetic carbene catalysts and the engineering of tailored enzymatic carbenes.

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube.

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU-), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU-, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU- inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU- is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU- to GLU0 to reduce the energy barriers. Overall, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

Enhancing Hit Discovery in Virtual Screening through Absolute Protein-Ligand Binding Free-Energy Calculations.

In the hit identification stage of drug discovery, a diverse chemical space needs to be explored to identify initial hits. Contrary to empirical scoring functions, absolute protein-ligand binding free-energy perturbation (ABFEP) provides a theoretically more rigorous and accurate description of protein-ligand binding thermodynamics and could, in principle, greatly improve the hit rates in virtual screening. In this work, we describe an implementation of an accurate and reliable ABFEP method in FEP+. We validated the ABFEP method on eight congeneric compound series binding to eight protein receptors including both neutral and charged ligands. For ligands with net charges, the alchemical ion approach is adopted to avoid artifacts in electrostatic potential energy calculations. The calculated binding free energies correlate with experimental results with a weighted average of R2 = 0.55 for the entire dataset. We also observe an overall root-mean-square error (RMSE) of 1.1 kcal/mol after shifting the zero-point of the simulation data to match the average experimental values. Through ABFEP calculations using apo versus holo protein structures, we demonstrated that the protein conformational and protonation state changes between the apo and holo proteins are the main physical factors contributing to the protein reorganization free energy manifested by the overestimation of raw ABFEP calculated binding free energies using the holo structures of the proteins. Furthermore, we performed ABFEP calculations in three virtual screening applications for hit enrichment. ABFEP greatly improves the hit rates as compared to docking scores or other methods like metadynamics. The good performance of ABFEP in rank ordering compounds demonstrated in this work confirms it as a useful tool to improve the hit rates in virtual screening, thus facilitating hit discovery.

Revealing the allosteric mechanism of pH-dependence in the proton-activated chloride channel.

The ability of the Proton-activated chloride (PAC) channel to perform its chloride transfer function across the membrane was shown to be dependent on the pH in the extracellular environment. Thus, an increase of extracellular pH from pH 4.5 to pH 8 results in a conformational change in transmembrane domain (TMD) of PAC and abolishes chloride conductivity. Mutagenesis of various extracellular domain (ECD) regions of PAC, including the Acidic Pocket, the beta 3 sheet of the “finger region” and the distal end of the ECD, altered the pH-sensitivity of PAC without revealing how such changes affect the pH sensitivity and its effect on the structure and dynamics of the distal end of the TMD. To understand this mechanism we used molecular pH-specific molecular dynamics (MD) simulations of PAC protein with the ECpH-MD method, capturing the pH-dependent conformational transitions between the active and inactive forms of PAC protein. The effect of changes in protonation states of a number of experimentally determined pH-sensor residues revealed their specific structural roles in the mechanism of conformational change underlying the changes in PAC activity. The allosteric pathways connecting these pH-sensors in the ECD to the TMD at different pH values revealed and quantified from analyses of the various MD trajectories with the rigorous N-body information theory (NbIT) method. Together with the structural data determined for various states of PAC (inactive, active, desensitized), the atomistic level of the dynamics driving the allosteric communication offers a comprehensive mechanistic understanding of the pH-dependence in PAC function.

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 pKa 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 pKa 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 pKa 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.

Interaction of Chondroitin and Hyaluronan Glycosaminoglycans with Surfaces of Carboxylated Carbon Nanotubes Studied Using Molecular Dynamics Simulations

Interaction of β-D-glucopyranuronic acid (GlcA), N-acetyl-β-D-glucosamine (GlcNAc), N-acetyl-β-D-galactosamine (GalNAc) and two natural decameric glycosaminoglycans, hyaluronic acid (HA) and Chondroitin (Ch) with carboxylated carbon nanotubes, were studied using molecular dynamics simulations in a condensed phase. The force field used for carbohydrates was the GLYCAM-06j version, while functionalized carbon nanotubes (fCNT) were described using version two of the general amber force field. We found a series of significant differences in carbohydrate-fCNT adsorption strength depending on the monosaccharide molecule and protonation state of surface carboxyl groups. GlcNAc and GalNAc reveal a strong adsorption on fCNT with deprotonated carboxyl groups, and a slightly weaker adsorption on the fCNT with protonated carboxyl groups. On the contrary, GlcA weakly adsorbs on fCNT. The change in protonation state of surface carboxyl groups leads to the reversal orientation of GlcNAc and GalNAc in reference to the fCNT surface, while GlcA is not sensitive to that factor. Adsorption of decameric oligomers on the surface of fCNT weakens with the increasing number of monosaccharide units. Chondroitin adsorbs weaker than hyaluronic acid and incorporation of four Ch molecules leads to partial detachment of them from the fCNT surface. The glycan-fCNT interactions are strong enough to alter the conformation of carbohydrate backbone; the corresponding conformational changes act toward a more intensive contact of glycan with the fCNT surface. Structural and energetic features of the adsorption process suggest the CH-π interaction-driven mechanism.

Changes in Protonation States of In-Pathway Residues can Alter Ligand Binding Pathways Obtained From Spontaneous Binding Molecular Dynamics Simulations.

Protein-ligand binding processes often involve changes in protonation states that can be key to recognize and orient the ligand in the binding site. The pathways through which (bio)molecules interplay to attain productively bound complexes are intricate and involve a series of interconnected intermediate and transition states. Molecular dynamics (MD) simulations and enhanced sampling techniques are commonly used to characterize the spontaneous binding of a ligand to its receptor. However, the effect of protonation state changes of in-pathway residues in spontaneous binding MD simulations remained mostly unexplored. Here, we used molecular dynamics simulations to reconstruct the trypsin-benzamidine binding pathway considering different protonation states of His57. This residue is part of the trypsin catalytic triad and is located more than 10 Å away from Asp189, which is responsible for benzamidine binding in the trypsin S1 pocket. Our MD simulations showed that the binding pathways that benzamidine follow to target the S1 binding site are critically dependent on the His57 protonation state. Binding of benzamidine frequently occurs when His57 is protonated in the delta nitrogen while the binding process is significantly less frequent when His57 is positively charged. Constant-pH MD simulations retrieved the equilibrium populations of His57 protonation states at trypsin active pH offering a clearer picture of benzamidine recognition and binding. These results indicate that properly accounting for protonation states of distal residues can be important in spontaneous binding MD simulations.

Distinct electrostatic frequency tuning rates for amide I and amide I′ vibrations

Amide I spectroscopy probes the backbone C=O stretch vibrations of peptides and proteins. Amide I spectra are often collected in deuterated water (D2O) since this provides a cleaner background in the amide I frequency range; such data are often referred to as amide I' spectra since deuteration induces changes in the mode structure, including a roughly ∼10 cm-1 redshift. For biological samples, however, deuteration is often not possible. As amide I frequency maps are increasingly applied to quantitative protein structural analysis, this raises the interesting challenge of drawing direct connections between amide I and amide I' data. We here analyze amide I and amide I' peak frequencies for a series of dipeptides and related compounds. Changes in protonation state induce large electrostatic shifts in the peak frequencies, allowing us to amass a sizable library of data points for direct amide I/amide I' comparison. While we find an excellent linear correlation between amide I and amide I' peak frequencies, the deuteration-induced shift is smaller for more red-shifted vibrations, indicating different electrostatic tuning rates in the two solvents. H2O/D2O shifts were negligible for proline-containing dipeptides that lack exchangeable amide hydrogens, indicating that the intrinsic properties of the solvent do not strongly influence the H/D shift. These findings indicate that the distinct tuning rates observed for the two vibrations arise from modifications to the intrinsic properties of the amide bond and provide (at least for solvated dipeptides) a simple, linear "map" for translating between amide I and amide I' frequencies.

A Computational Protocol Combining DFT and Cheminformatics for Prediction of pH-Dependent Redox Potentials.

Discovering new materials for energy storage requires reliable and efficient protocols for predicting key properties of unknown compounds. In the context of the search for new organic electrolytes for redox flow batteries, we present and validate a robust procedure to calculate the redox potentials of organic molecules at any pH value, using widely available quantum chemistry and cheminformatics methods. Using a consistent experimental data set for validation, we explore and compare a few different methods for calculating reaction free energies, the treatment of solvation, and the effect of pH on redox potentials. We find that the B3LYP hybrid functional with the COSMO solvation method, in conjunction with thermal contributions evaluated from BLYP gas-phase harmonic frequencies, yields a good prediction of pH = 0 redox potentials at a moderate computational cost. To predict how the potentials are affected by pH, we propose an improved version of the Alberty-Legendre transform that allows the construction of a more realistic Pourbaix diagram by taking into account how the protonation state changes with pH.

Structural basis of the protochromic green/red photocycle of the chromatic acclimation sensor RcaE

Cyanobacteriochromes (CBCRs) are bilin-binding photosensors of the phytochrome superfamily that show remarkable spectral diversity. The green/red CBCR subfamily is important for regulating chromatic acclimation of photosynthetic antenna in cyanobacteria and is applied for optogenetic control of gene expression in synthetic biology. It is suggested that the absorption change of this subfamily is caused by the bilin C15-Z/C15-E photoisomerization and a subsequent change in the bilin protonation state. However, structural information and direct evidence of the bilin protonation state are lacking. Here, we report a high-resolution (1.63Å) crystal structure of the bilin-binding domain of the chromatic acclimation sensor RcaE in the red-absorbing photoproduct state. The bilin is buried within a "bucket" consisting of hydrophobic residues, in which the bilin configuration/conformation is C5-Z,syn/C10-Z,syn/C15-E,syn with the A- through C-rings coplanar and the D-ring tilted. Three pyrrole nitrogens of the A- through C-rings are covered in the α-face with a hydrophobic lid of Leu249 influencing the bilin pKa, whereas they are directly hydrogen bonded in the β-face with the carboxyl group of Glu217. Glu217 is further connected to a cluster of waters forming a hole in the bucket, which are in exchange with solvent waters in molecular dynamics simulation. We propose that the "leaky bucket" structure functions as a proton exit/influx pathway upon photoconversion. NMR analysis demonstrated that the four pyrrole nitrogen atoms are indeed fully protonated in the red-absorbing state, but one of them, most likely the B-ring nitrogen, is deprotonated in the green-absorbing state. These findings deepen our understanding of the diverse spectral tuning mechanisms present in CBCRs.

Direct Observation of Protonation State Modulation in SARS-CoV-2 Main Protease upon Inhibitor Binding with Neutron Crystallography.

The main protease (3CL Mpro) from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is an essential enzyme for viral replication with no human counterpart, making it an attractive drug target. To date, no small-molecule clinical drugs are available that specifically inhibit SARS-CoV-2 Mpro. To aid rational drug design, we determined a neutron structure of Mpro in complex with the α-ketoamide inhibitor telaprevir at near-physiological (22 °C) temperature. We directly observed protonation states in the inhibitor complex and compared them with those in the ligand-free Mpro, revealing modulation of the active-site protonation states upon telaprevir binding. We suggest that binding of other α-ketoamide covalent inhibitors can lead to the same protonation state changes in the Mpro active site. Thus, by studying the protonation state changes induced by inhibitors, we provide crucial insights to help guide rational drug design, allowing precise tailoring of inhibitors to manipulate the electrostatic environment of SARS-CoV-2 Mpro.

Solution structure of a europium-nicotianamine complex supports that phytosiderophores bind lanthanides.

We report the solution structure of a europium-nicotianamine complex predicted from ab initio molecular dynamics simulations with density functional theory. Emission and excitation spectroscopy show that the Eu3+ coordination environment changes in the presence of nicotianamine, suggesting complex formation, such as what is seen for the Eu3+-nicotianamine complex structure predicted from computation. We modeled Eu3+-ligand complexes with explicit water molecules in periodic boxes, effectively simulating the solution phase. Our simulations consider possible chemical events (e.g. coordination bond formation, protonation state changes, charge transfers), as well as ligand flexibility and solvent rearrangements. Our computational approach correctly predicts the solution structure of a Eu3+-ethylenediaminetetraacetic acid complex within 0.05 Å of experimentally measured values, backing the fidelity of the predicted solution structure of the Eu3+-nicotianamine complex. Emission and excitation spectroscopy measurements were also performed on the well-known Eu3+-ethylenediaminetetraacetic acid complex to validate our experimental methods. The electronic structure of the Eu3+-nicotianamine complex is analyzed to describe the complexes in greater detail. Nicotianamine is a metabolic precursor of, and structurally very similar to, phytosiderophores, which are responsible for the uptake of metals in plants. Although knowledge that nicotianamine binds europium does not determine how plants uptake rare earths from the environment, it strongly supports that phytosiderophores bind lanthanides.

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 pKa values with respect to their pKa in bulk water. This work investigates the atomistic cause of the highly shifted pKa of the internal Glu23 in the artificially mutated variant V23E of Staphylococcal Nuclease (SNase) using pH replica exchange molecular dynamics (pH-REMD) simulations. The pKa of Glu23 obtained from our calculations is 6.55, which is elevated with respect to the glutamate pKa of 4.40 in bulk water. The calculated value is close to the experimental pKa of 7.10. Our simulations show that the highly shifted pKa 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 pKa values of Glu23 and describe the coupling between the conformational and ionization equilibria.

An acidic residue buried in the dimer interface of isocitrate dehydrogenase 1 (IDH1) helps regulate catalysis and pH sensitivity.

Isocitrate dehydrogenase 1 (IDH1) catalyzes the reversible NADP+-dependent conversion of isocitrate to α-ketoglutarate (αKG) to provide critical cytosolic substrates and drive NADPH-dependent reactions like lipid biosynthesis and glutathione regeneration. In biochemical studies, the forward reaction is studied at neutral pH, while the reverse reaction is typically characterized in more acidic buffers. This led us to question whether IDH1 catalysis is pH-regulated, which would have functional implications under conditions that alter cellular pH, like apoptosis, hypoxia, cancer, and neurodegenerative diseases. Here, we show evidence of catalytic regulation of IDH1 by pH, identifying a trend of increasing kcat values for αKG production upon increasing pH in the buffers we tested. To understand the molecular determinants of IDH1 pH sensitivity, we used the pHinder algorithm to identify buried ionizable residues predicted to have shifted pKa values. Such residues can serve as pH sensors, with changes in protonation states leading to conformational changes that regulate catalysis. We identified an acidic residue buried at the IDH1 dimer interface, D273, with a predicted pKa value upshifted into the physiological range. D273 point mutations had decreased catalytic efficiency and, importantly, loss of pH-regulated catalysis. Based on these findings, we conclude that IDH1 activity is regulated, at least in part, by pH. We show this regulation is mediated by at least one buried acidic residue ∼12 Å from the IDH1 active site. By establishing mechanisms of regulation of this well-conserved enzyme, we highlight catalytic features that may be susceptible to pH changes caused by cell stress and disease.