Protonation State Research Articles

PH Dependence of HSF1 trimerization is shaped by intramolecular interactions

Heat shock factor 1 (HSF1) primarily regulates various cellular stress responses. Previous studies have shown that low pH within the physiological range directly activates HSF1 function in vitro. However, the detailed molecular mechanisms remain unclear. This study proposes a molecular mechanism based on the trimerization behavior of HSF1 at different pH values. Extensive mutagenesis of human and goldfish HSF1 revealed that the optimal pH for trimerization depended on the identity of residue 103. In particular, when residue 103 was occupied by tyrosine, a significant increase in the optimal pH was observed, regardless of the rest of the sequence. This behavior can be explained by the protonation state of the neighboring histidine residues, His101 and His110. Residue 103 plays a key role in trimerization by forming disulfide or non-covalent bonds with Cys36. If tyrosine resides at residue 103 in an acidic environment, its electrostatic interactions with positively charged histidine residues prevent effective trimerization. His101 and His110 are neutralized at a higher pH, which releases Tyr103 to interact with Cys36 and drives the effective trimerization of HSF1. This study showed that the protonation state of a histidine residue can regulate the intramolecular interactions, which consequently leads to a drastic change in the oligomerization behavior of the entire protein.

1H, 15N and13C resonance assignments of S2A mutant of human carbonic anhydrase II.

In preparation for a detailed exploration of the structural and functional aspects of the Ser2Ala mutant of human carbonic anhydrase II, we present here almost complete sequence-specific resonance assignments for 1H, 15N, and 13C. The mutation of serine to alanine at position 2, located in the N-terminal region of the enzyme, significantly alters the hydrophilic nature of the site, rendering it hydrophobic. Consequently, there is an underlying assumption that this mutation would repel water from the site. However, intriguingly, comparative analysis of the mutant structure with the wild type reveals minimal discernible differences. These assignments serve as the basis for in-depth studies on histidine dynamics, protonation states, and its intricate role in protein-water interactions and catalysis.

In situ potassium and hydrogen ion exchange into a cubic zirconium silicate microporous material.

The selective separation of ions from aqueous systems, and even in the human body, is a crucial to overall environmental management and health. Nanoporous materials are widely known for their selective removal of cations from aqueous media, and therefore have been targeted for use as a pharmaceutical to treat hyperkalemia. This study investigated the detailed crystallographic molecular mechanisms that control the potassium ion selectivity in the nanoporous cubic zirconium silicate (CZS) related materials. Using time-resolved in situ Raman spectroscopy and time-resolved in situ X-ray diffraction, the selectivity mechanisms were determined to involve a synchronous cation-cation repulsion process that serves to open a favorable coordination bonding environment for potassium, not unlike the ion selectivity filter process found in potassium ion channels in proteins. Enhancement of ion exchange was observed when the CZS material was in a partial protonated state (≈3:1 Na:H), causing an expansion of the unit-cell volume, enlargement of the 7 member-ring window, and distortion of framework polyhedra, which allowed increased accessibility to the cage structures and resulted in rapid irreversible potassium ion exchange.

Histidine Protonation States Regulate the State Transition from R State Hemoglobin.

The objective of our work is to investigate the impact of pH on the structural changes of hemoglobin that affect its O2 affinity, known as the Bohr effect. We conducted molecular dynamics (MD) simulations to explore the transition between various hemoglobin states based on the protonation states (PSs) of two histidine residues (βHis143 and βHis146). We conducted the MD simulations from the R and R2 states with three sets of PSs assuming pH values of 7.0, 6.5, and 5.5, aiming to investigate the influence of pH on hemoglobin behavior. Our results demonstrated that the protonated His residues promote the state transition from the R state to the R2 state and encourage elongation of the distance between the β1-β2 subunits by weakening the inter-subunit interactions in the R state. These observations, aligning with the experimental evidence, indicate that the R2 state typically crystallizes under low pH conditions. Our findings suggest that the relationship between the PSs and the structural stability of the R state plays a role in the acid and alkaline Bohr effect.

Elucidating the Structure of the Eu‐EDTA Complex in Solution at various Protonation States

Ethylenediaminetetraacetic acid (EDTA), which has two amine and four carboxylate protonation sites, forms stable complexes with lanthanide ions. This work analyzes the coordination structure, in atomic resolution, of the Eu3+ ion complexed with EDTA in all its protonation states in aqueous solution. Eu‐EDTA complexes were modeled using classical molecular dynamics (MD) simulations using force field parameters optimized with ab initio molecular dynamics (AIMD) simulations. Structures from the MD simulations were used to predict extended X‐ray absorption fine structure (EXAFS) spectra and compared with EXAFS measurements of the Eu3+ aqua ion and Eu‐EDTA complexes at pH 3 and 11. This work details how Eu‐EDTA complex coordination structures change with increasing protonation of the EDTA ligand in the complex, from the tightly bound unprotonated complex to the unbinding of the fully protonated EDTA ligand from the Eu3+ ion as both become solvated by water. Agreement between predicted and measured EXAFS spectra supports the findings from simulation.

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.

Molecular Dynamics Simulations of Native Protein Charging via Proton Transfer during Electrospray Ionization with Grotthuss Diffuse H3O.

Unraveling the mechanism by which native proteins are charged through electrospray ionization (ESI) has been the focus of considerable research because observable charge states can be correlated to biophysical characteristics, such as protein folding and, thus, solution conformation. Difficulties in characterizing electrosprayed droplets have catalyzed the use of molecular dynamics (MD) to provide insights into the mechanisms by which proteins are charged and transferred to the gas phase. However, prior MD studies have utilized metal ions, primarily Na+, as charge carriers, even though proteins are primarily detected as protonated ions in the mass spectra. Here, we propose a modified MD protocol for simulating discrete Grotthuss diffuse H3O+ that is capable of dynamically altering amino-acid protonation states to model electrospray charging and gaseous ion formation of model proteins, ubiquitin, and myoglobin. Application of the protocol to the evaporation of acidic droplets enables a molecular perspective of H3O+ coordination and proton transfer to/from proteins, which is unfeasible with the metal charge carriers used in previous MD studies of ESI. Our protocol recreates experimentally observed charge-state distributions and supports the charge residue model (CRM) as the dominant mechanism of native protein ionization during ESI. Additionally, our results suggest that protonation is highly specific to individual residues and is correlated to the formation of localized hydrated regions on the protein surface as droplets desolvate. Considering the use of discrete H3O+ instead of Na+, the developed protocol is a necessary step toward developing a more comprehensive model of protein ionization during ESI.

QM/MM Study of the Catalytic Mechanism and Substrate Specificity of the Aromatic Substrate C-Methyltransferase Fur6.

In the field of medical chemistry and other organic chemistry, introducing a methyl group into a designed position has been difficult to achieve. However, owing to the vigorous developments in the field of enzymology, methyltransferases are considered potential tools for addressing this problem. Within the methyltransferase family, Fur6 catalyzes the methylation of C3 of 1,2,4,5,7-pentahydroxynaphthalene (PHN) using S-adenosyl-l-methionine (SAM) as the methyl donor. Here, we report the catalytic mechanism and substrate specificity of Fur6 based on computational studies. Our molecular dynamics (MD) simulation studies reveal the reactive form of PHN and its interactions with the enzyme. Our hybrid quantum mechanics/molecular mechanics (QM/MM) calculations suggest the reaction pathway of the methyl transfer step in which the energy barrier is 8.6 kcal mol-1. Our free-energy calculations with a polarizable continuum model (PCM) indicate that the final deprotonation step of the methylated intermediate occurs after it is ejected into the water solvent from the active center pocket of Fur6. Additionally, our studies on the protonation states, the highest occupied molecular orbital (HOMOs), and the energy barriers of the methylation reaction for the analogs of PHN demonstrate the mechanism of the specificity to PHN. Our study provides valuable insights into Fur6 chemistry, contributing to a deeper understanding of molecular mechanisms and offering an opportunity to engineer the enzyme to achieve high yields of the desired product(s).

Metal Complexes of Redox Non-Innocent Ligand N,N'-Bis(3,5-di-tertbutyl-2-hydroxy-phenyl)-1,2-phenylenediamine.

Redox non-innocent ligands react with metal precursors to form complexes where the oxidation states of the ligand and thus the metal atom cannot be easily defined. A well-known example of such ligands is bis(o-aminophenol) N,N'-bis(3,5-di-tertbutyl-2-hydroxy-phenyl)-1,2-phenylenediamine, previously developed by the Wieghardt group, which has a potentially tetradentate coordination mode and four distinct protonation states, whereas its electrochemical behavior allows for five distinct oxidation states. This rich redox chemistry, as well as the ability to coordinate to various transition metals, has been utilized in the syntheses of metal complexes with M2L, ML and ML2 stoichiometries, sometimes supported with other ligands. Different oxidation states of the ligand can adopt different coordination modes. For example, in the fully oxidized form, two N donors are sp2-hybridized, which makes the ligand planar, whereas in the fully reduced form, the sp3-hybridized N donors allow the formation of more flexible chelate structures. In general, the metal can be reduced during complexation, but redox processes of the isolated complexes typically occur on the ligand. Combination of this non-innocent ligand with redox-active transition metals may lead to complexes with interesting magnetic, electrochemical, photonic and catalytic properties.

Hyperpolarized long-lived spin state of methylene protons of 2-bromoethanol obtained from ethylene with non-equilibrium nuclear spin order

In this work we achieve a significant overpopulation (PLLS≈1%) of the long-lived spin state (LLS) of methylene protons in 2-bromoethan(2H)ol (BrEtOD) obtained in a reaction between ethylene with non-equilibrium nuclear spin order and bromine water. Given all protons in ethylene are magnetically equivalent, its nuclear states are classified into nuclear spin isomers (NSIM) with total spin I = 2,1,0. Addition of parahydrogen to acetylene produces ethylene with a population of only those NSIMs with I = 1,0. As a result of the reaction with bromine water the non-equilibrium spin order of ethylene is partly transferred to the singlet LLS involving the two methylene groups of BrEtOD. The 1H NMR signal enhancement (SE≈200) obtained as a result of the LLS readout is approximately equal to the SE of the hyperpolarized BrEtOD obtained with a single π/4 pulse. The LLS relaxation time (TLLS) was shown to be approximately 40 s (≈8T1) in the argon-bubbled sample.

Mapping the transporter-substrate interactions of the Trypanosoma cruzi NB1 nucleobase transporter reveals the basis for its high affinity and selectivity for hypoxanthine and guanine and lack of nucleoside uptake

Trypanosoma cruzi is a protozoan parasite and the etiological agent of Chagas disease, a debilitating and sometimes fatal disease that continues to spread to new areas. Yet, Chagas disease is still only treated with two related nitro compounds that are insufficiently effective and cause severe side effects. Nucleotide metabolism is one of the known vulnerabilities of T. cruzi, as they are auxotrophic for purines, and nucleoside analogues have been shown to have genuine promise against this parasite in vitro and in vivo. Since purine antimetabolites require efficient uptake through transporters, we here report a detailed characterisation of the T. cruzi NB1 nucleobase transporter with the aim of elucidating the interactions between TcrNB1 and its substrates and finding the positions that can be altered in the design of novel antimetabolites without losing transportability. Systematically determining the inhibition constants (Ki) of purine analogues for TcrNB1 yielded their Gibbs free energy of interaction, ΔG0. Pairwise comparisons of substrate (hypoxanthine, guanine, adenine) and analogues allowed us to determine that optimal binding affinity by TcrNB1 requires interactions with all four nitrogen residues of the purine ring, with N1 and N9, in protonation state, functioning as presumed hydrogen bond donors and unprotonated N3 and N7 as hydrogen bond acceptors. This is the same interaction pattern as we previously described for the main nucleobase transporters of Trypanosoma brucei spp. and Leishmania major and makes it the first of the ENT-family genes that is functionally as well as genetically conserved between the three main kinetoplast pathogens.

On the Validity of Constant pH Simulations.

Constant pH (cpH) simulations are now a standard tool for investigating charge regulation in coarse-grained models of polyelectrolytes and colloidal systems. Originally developed for studying solutions with implicit ions, extending this method to systems with explicit ions or solvents presents several challenges. Ensuring proper charge neutrality within the simulation cell requires performing titration moves in sync with the insertion or deletion of ions, a crucial aspect often overlooked in the literature. Contrary to the prevailing views, cpH simulations are inherently grand-canonical, meaning that the controlled pH is that of the reservoir. The presence of the Donnan potential between the implicit reservoir and the simulation cell introduces significant differences between titration curves calculated for open and closed systems; the pH of an isolated (closed) system is different from the pH of the reservoir for the same protonation state of the polyelectrolyte. To underscore this point, in this paper, we will compare the titration curves calculated using the usual cpH algorithm with those from the exact canonical simulation algorithm. In the latter case, titration moves adhere to the correct detailed balance condition, and pH is calculated using the recently introduced surface Widom insertion algorithm. Our findings reveal a very significant difference between the titration isotherms obtained using the standard cpH algorithm and the canonical titration algorithm, emphasizing the importance of using the correct simulation approach when studying charge regulation of polyelectrolytes, proteins, and colloidal particles.

On the antibacterial photodynamic inactivation mechanism of Emodin and Dermocybin natural photosensitizers: A theoretical investigation.

A DFT and TDDFT study has been carried out on monomeric anthraquinones Emodin and Dermocybin (Em, Derm) recently proposed as natural antibacterial photosensitizers able to act also against gram-negative microbes. The computational study has been performed considering the relative amount of neutral and ionic forms of each compound in water, with the variation of pH. The occurrence of both Type I and Type II photoreactions has been explored computing the absorption properties of each species, the spin-orbit coupling constants (SOC), the vertical ionization potentials and the vertical electron affinities. The most plausible deactivation channels leading to the population of excited triplet states have been proposed. Our data indicate Emodin as more active than Dermocybin in antimicrobial photodynamic therapy throughout the Type II mechanism. Our data support a dual TypeI/II activity of the monomeric anthraquinones Emodin and Dermccybin in water, in all the considered protonation states.

Determination of Amino Acids' pKa: Importance of Cavity Scaling within Implicit Solvation Models and Choice of DFT Functionals.

Protonation states of molecules significantly influence the thermodynamics and kinetics of chemical reactions. This is especially important in biochemical processes, where appropriate protonation states of amino acids control the exo/endoergicity of practically all biochemical cycles. This paper is focused on appraisal of the impact of DFT functionals and PCM solvation models on the accuracy of pKa evaluations for all proteinogenic amino acids. Eight functionals (B3LYP, PBE0, revPBE0, M06-2X, M11, M11-L, TPSSh, and ωB97X-D) and four basis sets are considered, together with four kinds of implicit solvation models when additional attention is paid to a cavity construction. An influence of nonelectrostatic contributions and Wertz's corrections on Gibbs free energy is investigated together with accuracy of provided proton solvation energy. The best model is based on the M06-2X/6-311++G**/D-PCM/UAKS computational level. The fitting procedure is utilized to improve the accuracy of the evaluated models. All of these results are also compared with values obtained from the COSMOtherm program and CCSD(T) calculations. Results for cysteine and histidine are discussed individually, as they can be found in different protonation states at neutral pH.

The Interaction Mechanism of Picolinamide Fungicide Targeting on the Cytochrome bc1 Complex and Its Structural Modification.

Picolinamide fungicides, structurally related to UK-2A and antimycin-A, bind into the Qi-site in the bc1 complex. However, the detailed binding mode of picolinamide fungicides remains unknown. In the present study, antimycin-A and UK-2A were selected to study the binding mode of picolinamide inhibitors with four protonation states in the Qi-site by integrating molecular dynamics simulation, molecular docking, and molecular mechanics Generalized Born surface area (MM/GBSA) calculations. Subsequently, a series of new picolinamide derivatives were designed and synthesized to further understand the effects of substituents on the tail phenyl ring. The computational results indicated that the substituted aromatic rings in antimycin-A and UK-2A were the pharmacophore fragments and made the primary contribution when bound to a protein. Compound 9g-hydrolysis formed H-bonds with Hie201 and Ash228 and showed an IC50 value of 6.05 ± 0.24 μM against the porcine bc1 complex. Compound 9c, with a simpler chemical structure, showed higher control effects than florylpicoxamid against cucumber downy mildew and expanded the fungicidal spectrum of picolinamide fungicides. The structural and mechanistic insights obtained from the present study will provide a valuable clue for the future designing of new promising Qi-site inhibitors.

Scalar Relativistic All-Electron and Pseudopotential Ab Initio Study of a Minimal Nitrogenase [Fe(SH)4H]- Model Employing Coupled-Cluster and Auxiliary-Field Quantum Monte Carlo Many-Body Methods.

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to organisms. The detailed mechanism of this enzyme is currently not known, and computational studies are complicated by the fact that different density functional theory (DFT) methods give very different energetic results for calculations involving nitrogenase models. Recently, we designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S to mimic different possible protonation states of the nitrogenase active site. We showed that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical methods. Based on nonrelativistic single-reference coupled-cluster (CC) calculations, we estimated that the ΔE is 101 kJ/mol. In this study, we demonstrate that scalar relativistic effects play an important role and significantly affect ΔE. Our best revised single-reference CC estimates for ΔE are 85-91 kJ/mol, including energy corrections to account for contributions beyond triples, core-valence correlation, and basis-set incompleteness error. Among coupled-cluster approaches with approximate triples, the canonical CCSD(T) exhibits the largest error for this problem. Complementary to CC, we also used phaseless auxiliary-field quantum Monte Carlo calculations (ph-AFQMC). We show that with a Hartree-Fock (HF) trial wave function, ph-AFQMC reproduces the CC results within 5 ± 1 kJ/mol. With multi-Slater-determinant (MSD) trials, the results are 82-84 ± 2 kJ/mol, indicating that multireference effects may be rather modest. Among the DFT methods tested, τ-HCTH, r2SCAN with 10-13% HF exchange with and without dispersion, and O3LYP/O3LYP-D4, and B3LYP*/B3LYP*-D4 generally perform the best. The r2SCAN12 (with 12% HF exchange) functional mimics both the best reference MSD ph-AFQMC and CC ΔE results within 2 kJ/mol.

PH-dependent reaction triggering in PmHMGR crystals for time-resolved crystallography

Serial crystallography and time-resolved data collection can readily be employed to investigate the catalytic mechanism of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl (HMG)-coenzyme-A (CoA) reductase (PmHMGR) by changing the environmental conditions in the crystal and so manipulating the reaction rate. This enzyme uses a complex mechanism to convert mevalonate to HMG-CoA using the co-substrate CoA and cofactor NAD+. The multi-step reaction mechanism involves an exchange of bound NAD+ and large conformational changes by a 50-residue subdomain. The enzymatic reaction can be run in both forward and reverse directions in solution and is catalytically active in the crystal for multiple reaction steps. Initially, the enzyme was found to be inactive in the crystal starting with bound mevalonate, CoA, and NAD+. To observe the reaction from this direction, we examined the effects of crystallization buffer constituents and pH on enzyme turnover, discovering a strong inhibition in the crystallization buffer and a controllable increase in enzyme turnover as a function of pH. The inhibition is dependent on ionic concentration of the crystallization precipitant ammonium sulfate but independent of its ionic composition. Crystallographic studies show that the observed inhibition only affects the oxidation of mevalonate but not the subsequent reactions of the intermediate mevaldehyde. Calculations of the pKa values for the enzyme active site residues suggest that the effect of pH on turnover is due to the changing protonation state of His381. We have now exploited the changes in ionic inhibition in combination with the pH-dependent increase in turnover as a novel approach for triggering the PmHMGR reaction in crystals and capturing information about its intermediate states along the reaction pathway.

Controlling the Fluctuating Tip-Enhanced Raman Spectra of Chloramben on Silver Nanocubes.

Tip-enhanced Raman (TER) scattering from molecules residing at plasmonic junctions can be used to detect, identify, and image single molecules. This is most evident for flat molecules interrogated under conditions of extreme temperatures and pressure. It is also the case for (bio)molecular systems that feature preferred orientations/conformations under ambient laboratory conditions. More complex molecules that can adopt multiple conformations and/or feature different protonation or charge states give rise to complex TER spectra. We illustrate how the latter can be controlled in the case of chloramben molecules coated onto plasmonic silver nanocubes. We show that characteristic molecular Raman spectra cannot be obtained when tunneling plasmons are operative, i.e., when the tip is in direct contact with the chemically functionalized plasmonic nanoparticles. We rationalize these observations and propose an approach to less invasive and hence more analytical TER spectral imaging.

Inhomogenous broadening of antenna protein chlorophyll site energies due to side chain conformational and protonation states determined by Monte Carlo microstate analysis

Inhomogenous broadening of antenna protein chlorophyll site energies due to side chain conformational and protonation states determined by Monte Carlo microstate analysis

Surface Charge Can Modulate Phase Separation of Multidomain Proteins.

Phase separation has emerged as an important mechanism explaining the formation of certain biomolecular condensates. Biological phase separation is often driven by the multivalent interactions of modular protein domains. Beyond valency, the physical features of folded domains that promote phase separation are poorly understood. We used a model system─the small ubiquitin modifier (SUMO) and its peptide ligand, the SUMO interaction motif (SIM)─to examine how domain surface charge influences multivalency-driven phase separation. Phase separation of polySUMO and polySIM was altered by pH via a change in the protonation state of SUMO surface histidines. These effects were recapitulated by histidine mutations, which modulated SUMO solubility and polySUMO-polySIM phase separation in parallel and were quantitatively explained by atomistic modeling of weak interactions among proteins in the system. Thus, surface charge can tune the phase separation of multivalent proteins, suggesting a means of controlling phase separation biologically, evolutionarily, and therapeutically.