Proton transfer reactions and hydrogen-bond networks in protein environments

Characterizing Protein Protonation Microstates Using Monte Carlo Sampling

Proteins are polyelectrolytes with acidic and basic amino acids Asp, Glu, Arg, Lys, and His, making up ≈25% of the residues. The protonation state of residues, cofactors, and ligands defines a "protonation microstate". In an ensemble of proteins some residues will be ionized and others neutral, leading to a mixture of protonation microstates rather than in a single one as is often assumed. The microstate distribution changes with pH. The protein environment also modifies residue proton affinity so microstate distributions change in different reaction intermediates or as ligands are bound. Particular protonation microstates may be required for function, while others exist simply because there are many states with similar energy. Here, the protonation microstates generated in Monte Carlo sampling in MCCE are characterized in HEW lysozyme as a function of pH and bacterial photosynthetic reaction centers (RCs) in different reaction intermediates. The lowest energy and highest probability microstates are compared. The ΔG, ΔH, and ΔS between the four protonation states of Glu35 and Asp52 in lysozyme are shown to be calculated with reasonable precision. At pH 7 the lysozyme charge ranges from 6 to 10, with 24 accepted protonation microstates, while RCs have ≈50,000. A weighted Pearson correlation analysis shows coupling between residue protonation states in RCs and how they change when the quinone in the Q_B site is reduced. Protonation microstates can be used to define input MD parameters and provide insight into the motion of protons coupled to reactions.

The Solvation-Induced Onsager Reaction Field Rather than the Double-Layer Field Controls CO2 Reduction on Gold

The selectivity and activity of the carbon dioxide reduction (CO₂R) reaction are sensitive functions of the electrolyte cation. By measuring the vibrational Stark shift of in situ-generated CO on Au in the presence of alkali cations, we quantify the total electric field present at catalytic active sites and deconvolute this field into contributions from (1) the electrochemical Stern layer and (2) the Onsager (or solvation-induced) reaction field. Contrary to recent theoretical reports, the CO₂R kinetics does not depend on the Stern field but instead is closely correlated with the strength of the Onsager reaction field. These results show that in the presence of adsorbed (bent) CO₂, the Onsager field greatly exceeds the Stern field and is primarily responsible for CO₂ activation. Additional measurements of the cation-dependent water spectra using vibrational sum frequency generation spectroscopy show that interfacial solvation strongly influences the CO₂R activity. These combined results confirm that the cation-dependent interfacial water structure and its associated electric field must be explicitly considered for accurate understanding of CO₂R reaction kinetics.

The Phosphorus Bond, or the Phosphorus-Centered Pnictogen Bond: The Covalently Bound Phosphorus Atom in Molecular Entities and Crystals as a Pnictogen Bond Donor

The phosphorus bond in chemical systems, which is an inter- or intramolecular noncovalent interaction, occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a covalently or coordinately bonded phosphorus atom in a molecular entity and a nucleophile in another, or the same, molecular entity. It is the second member of the family of pnictogen bonds, formed by the second member of the pnictogen family of the periodic table. In this overview, we provide the reader with a snapshot of the nature, and possible occurrences, of phosphorus-centered pnictogen bonding in illustrative chemical crystal systems drawn from the ICSD (Inorganic Crystal Structure Database) and CSD (Cambridge Structural Database) databases, some of which date back to the latter part of the last century. The illustrative systems discussed are expected to assist as a guide to researchers in rationalizing phosphorus-centered pnictogen bonding in the rational design of molecular complexes, crystals, and materials and their subsequent characterization.

Exploring Intra- and Intermolecular Interactions in Selected N-Oxides-The Role of Hydrogen Bonds

Intra- and intermolecular interactions have been explored in selected N-oxide derivatives: 2-(N,N-dimethylamino-N-oxymethyl)-4,6-dimethylphenyl (1) and 5,5'-dibromo-3-diethylaminomethyl-2,2'-biphenol N-oxide (2). Both compounds possess intramolecular hydrogen bonding, which is classified as moderate in 1 and strong in 2, and resonance-assisted in both cases. Density Functional Theory (DFT) in its classical formulation as well as Time-Dependent extension (TD-DFT) were employed to study proton transfer phenomena. The simulations were performed in the gas phase and with implicit and explicit solvation models. The obtained structures of the studied N-oxides were compared with experimental data available. The proton reaction path was investigated using scan with an optimization method, and water molecule reorientation in the monohydrate of 1 was found upon the proton scan progress. It was found that spontaneous proton transfer phenomenon cannot occur in the electronic ground state of the compound 1. An opposite situation was noticed for the compound 2. The changes of nucleophilicity and electrophilicity upon the bridged proton migration were analyzed on the basis of Fukui functions in the case of 1. The interaction energy decomposition of dimers and microsolvation models was investigated using Symmetry-Adapted Perturbation Theory (SAPT). The simulations were performed in both phases to introduce polar environment influence on the interaction energies. The SAPT study showed rather minor role of induction in the formation of homodimers. However, it is worth noticing that the same induction term is responsible for the preference of water molecules' interaction with N-oxide hydrogen bond acceptor atoms in the microsolvation study. The Natural Bond Orbital (NBO) analysis was performed for the complexes with water to investigate the charge flow upon the polar environment introduction. Finally, the TD-DFT was applied for isolated molecules as well as for microsolvation models showing that the presence of solvent affects excited states, especially when the N-oxide acceptor atom is microsolvated.

Requirement of Chloride for the Downhill Electron Transfer Pathway from the Water-Splitting Center in Natural Photosynthesis

In photosystem II (PSII), Cl^- is a prerequisite for the second flash-induced oxidation of the Mn₄CaO₅ cluster (the S₂ to S₃ transition). We report proton transfer from the substrate water molecule via D1-Asp61 and electron transfer via redox-active D1-Tyr161 (TyrZ) to the chlorophyll pair in Cl^--depleted PSII using a quantum mechanical/molecular mechanical approach. The low-barrier H-bond formation between the substrate water molecule and D1-Asp61 remained unaffected upon the depletion of Cl^-. However, the binding site, D2-Lys317, formed a salt bridge with D1-Asp61, leading to the inhibition of the subsequent proton transfer. Remarkably, the redox potential (E_m) of S₂/S₃ increased significantly, making electron transfer from S₂ to TyrZ energetically uphill, as observed in Ca²⁺-depleted PSII. The uphill electron transfer pathway was induced by the significant increase in E_m(S₂/S₃) caused by the loss of charge compensation for D2-Lys317 upon the depletion of Cl^-, whereas it was induced by the significant decrease in E_m(TyrZ) caused by the rearrangement of the water molecules at the Ca²⁺ binding moiety upon the depletion of Ca²⁺.

Elucidation of the electron transfer environment in the MMOR FAD-binding domain from Methylosinus sporium 5

By facilitating electron transfer to the hydroxylase diiron center, MMOR-a reductase-serves as an essential component of the catalytic cycle of soluble methane monooxygenase. Here, the X-ray structure analysis of the FAD-binding domain of MMOR identified crucial residues and its influence on the catalytic cycle.

Low-Field Functional Group Resolved Nuclear Spin Relaxation in Mesoporous Silica

Solid-fluid interactions underpin the efficacy of functional porous materials across a diverse array of chemical reaction and separation processes. However, detailed characterization of interfacial phenomena within such systems is hampered by their optically opaque nature. Motivated by the need to bridge this capability gap, we report low-magnetic-field two-dimensional (2D) ¹H nuclear spin relaxation measurements as a noninvasive probe of adsorbate identity and interfacial dynamics, exploring the relaxation characteristics exhibited by liquid hydrocarbon adsorbates confined to a model mesoporous silica. For the first time, we demonstrate the capacity of this approach in distinguishing functional group-specific relaxation phenomena across a diverse range of alcohols and carboxylic acids employed as solvents, reagents, and liquid hydrogen carriers, with distinct relaxation responses assigned to the alkyl and hydroxyl moieties of each confined liquid. Uniquely, this relaxation behavior is shown to correlate with adsorbate acidity, with the observed relationship rationalized on the basis of surface-adsorbate proton-exchange dynamics. Our results demonstrate that nuclear spin relaxation provides a molecular-level perspective on sorbent/sorbate interactions, motivating the exploration of such measurements as a unique probe of adsorbate identity within optically opaque porous media.

Proton in the ring: spectroscopy and dynamics of proton bonding in macrocycle cavities

The proton bond is a paradigmatic quantum molecular interaction and a major driving force of supramolecular chemistry. The ring cavities of crown ethers provide an intriguing environment, promoting competitive proton sharing with multiple coordination anchors. This study shows that protons confined in crown ether cavities form dynamic bonds that migrate to varying pairs of coordinating atoms when allowed by the flexibility of the macrocycle backbone. Prototypic native crown ethers (12-crown-4, 15-crown-5 and 18-crown-6) and aza-crown ethers (cyclen, 1-aza-18-crown-6 and hexacyclen) are investigated. For each system, Infrared action spectroscopy experiments and ab initio Molecular Dynamics computations are employed to elucidate the structural effects associated with proton diffusion and its entanglement with the conformational and vibrational dynamics of the protonated host.

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller-Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers' stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm-1-1700 cm-1 and 2300 cm-1-3400 cm-1 in the gas phase and 600 cm-1-1800 cm-1 and 2200 cm-1-3400 cm-1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm-1-1700 cm-1 and 2300 cm-1-3300 cm-1 for the gas phase and one broad absorption region in the solid state between 700 cm-1 and 3100 cm-1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.

Formation of the Metal-Binding Core of the ZRT/IRT-like Protein (ZIP) Family Zinc Transporter

Zinc homeostasis in mammals is constantly and precisely maintained by sophisticated regulatory proteins. Among them, the Zrt/Irt-like protein (ZIP) regulates the influx of zinc into the cytoplasm. In this work, we have employed all-atom molecular dynamics simulations to investigate the Zn²⁺ transport mechanism in prokaryotic ZIP obtained from Bordetella bronchiseptica (BbZIP) in a membrane bilayer. Additionally, the structural and dynamical transformations of BbZIP during this process have been analyzed. This study allowed us to develop a hypothesis for the zinc influx mechanism and formation of the metal-binding site. We have created a model for the outward-facing form of BbZIP (experimentally only the inward-facing form has been characterized) that has allowed us, for the first time, to observe the Zn²⁺ ion entering the channel and binding to the negatively charged M2 site. It is thought that the M2 site is less favored than the M1 site, which then leads to metal ion egress; however, we have not observed the M1 site being occupied in our simulations. Furthermore, removing both Zn²⁺ ions from this complex resulted in the collapse of the metal-binding site, illustrating the "structural role" of metal ions in maintaining the binding site and holding the proteins together. Finally, due to the long Cd²⁺-residue bond distances observed in the X-ray structures, we have proposed the existence of an H₃O⁺ ion at the M2 site that plays an important role in protein stability in the absence of the metal ion.

Mechanism of the formation of proton transfer pathways in photosynthetic reaction centers

The crystal structures of photosynthetic reaction centers from purple bacteria (PbRCs) and photosystem II show large structural similarity. However, the proposed mechanisms of proton transfer toward the terminal electron acceptor quinone (Q_B) are not consistent. In particular, not His-L190, which is an H-bond partner of Q_B, but rather Glu-L212, which is ∼6 Å away from Q_B, was assumed to be the direct proton donor for Q_B. We demonstrate that the H-bond between His-L190 and Q_B is a low-barrier H-bond, which facilitates proton transfer from singly protonated His-L190 to Q_B. Furthermore, Glu-L212 is not a direct H-bond donor for Q_B. However, it facilitates proton transfer toward deprotonated His-L190 via water molecules after Q_BH₂ forms and leaves the PbRC.

In photosynthetic reaction centers from purple bacteria (PbRCs) from Rhodobacter sphaeroides, the secondary quinone Q_B accepts two electrons and two protons via electron-coupled proton transfer (PT). Here, we identify PT pathways that proceed toward the Q_B binding site, using a quantum mechanical/molecular mechanical approach. As the first electron is transferred to Q_B, the formation of the Grotthuss-like pre-PT H-bond network is observed along Asp-L213, Ser-L223, and the distal Q_B carbonyl O site. As the second electron is transferred, the formation of a low-barrier H-bond is observed between His-L190 at Fe and the proximal Q_B carbonyl O site, which facilitates the second PT. As Q_BH₂ leaves PbRC, a chain of water molecules connects protonated Glu-L212 and deprotonated His-L190 forms, which serves as a pathway for the His-L190 reprotonation. The findings of the second pathway, which does not involve Glu-L212, and the third pathway, which proceeds from Glu-L212 to His-L190, provide a mechanism for PT commonly used among PbRCs.

Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient

Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.

Two Distinct Oxygen-Radical Conformations in the X-ray Free Electron Laser Structures of Photosystem II

We report the existence of two distinct oxygen-radical-containing Mn₄CaO_5/6 conformations with short O···O bonds in the crystal structures of the oxygen-evolving enzyme photosystem II (PSII), obtained using an X-ray free electron laser (XFEL). A short O···O distance of <2.3 Å between the O4 site of the Mn₄CaO₅ complex and the adjacent water molecule (W539) in the proton-conducting O4-water chain was observed in the second flash-induced (2F) XFEL structure (2F-XFEL), which may correspond to S₃. By use of a quantum mechanical/molecular mechanical approach, the OH^• formation at W539 and the short O4···O_W539 distance (<2.3 Å) were reproduced in S₂ and S₃ with reduced Mn1(III), which lacks the additional sixth water molecule O6. As the O^•- formation at O6 and the short O5···O6 distance (1.9 Å) have been reported in another 2F-XFEL structure with reduced Mn4(III), two distinct oxygen-radical conformations exist in the 2F-XFEL crystals.

CO2 Hydrogenation Catalyzed by a Ruthenium Protic N-Heterocyclic Carbene Complex

We describe the hydrogenation of CO₂ to formate catalyzed by a Ru(II) bis(protic N-heterocyclic carbene, p-NHC) phosphine complex [Ru(bpy)(MeCN)(P^Ph(p-NHC)₂)](PF₆)₂ (1). Under catalytic conditions (20 μmol catalyst, 20 bar CO₂, 60 bar H₂, 5 mL THF, 140 °C, 16 h), the activity of 1 is limited only by the amount of K₃PO₄ present in the reaction, yielding a nearly 1:1 ratio of turnover number (TON) to equivalents of K₃PO₄ (relative to 1), with the highest TON = 8040. Additionally, analysis of the reaction solution post-run reveals the catalyst intact with no free ligand observed. Stoichiometric studies, including examination of unique carbamate and hydride complexes as relevant intermediates, were carried out to probe the operative mechanism and understand the importance of metal-ligand cooperativity in this system.

Mitophagy in Human Diseases

Mitophagy is a selective autophagic process, essential for cellular homeostasis, that eliminates dysfunctional mitochondria. Activated by inner membrane depolarization, it plays an important role during development and is fundamental in highly differentiated post-mitotic cells that are highly dependent on aerobic metabolism, such as neurons, muscle cells, and hepatocytes. Both defective and excessive mitophagy have been proposed to contribute to age-related neurodegenerative diseases, such as Parkinson's and Alzheimer's diseases, metabolic diseases, vascular complications of diabetes, myocardial injury, muscle dystrophy, and liver disease, among others. Pharmacological or dietary interventions that restore mitophagy homeostasis and facilitate the elimination of irreversibly damaged mitochondria, thus, could serve as potential therapies in several chronic diseases. However, despite extraordinary advances in this field, mainly derived from in vitro and preclinical animal models, human applications based on the regulation of mitochondrial quality in patients have not yet been approved. In this review, we summarize the key selective mitochondrial autophagy pathways and their role in prevalent chronic human diseases and highlight the potential use of specific interventions.

Nanoconfinement Raises the Energy Barrier to Hydrogen Atom Exchange between Water and Glucose

In bulk aqueous environments, the exchange of protons between labile hydroxyl groups typically occurs easily and quickly. Nanoconfinement can dramatically change this normally facile process. Through exchange spectroscopy (EXSY) NMR measurements, we observe that nanoconfinement of glucose and water within AOT (sodium bis(2-ethylhexyl) sulfosuccinate) reverse micelles raises the energy barrier to labile hydrogen exchange, which suggests a disruption of the hydrogen bond network. Near room temperature, we measure barriers high enough to slow the process by as much as 2 orders of magnitude. Although exchange rates slow with decreasing temperatures in these nanoconfined environments, the barrier we measure below ∼285 K is 3-5 times lower than the barrier measured at room temperature, indicating a change in mechanism for the process. These findings suggest the possibility of hydrogen tunneling at a surprisingly high-temperature threshold. Furthermore, differences in exchange rates depend on the hydroxyl group position on the glucose pyranose ring and suggest a net orientation of glucose at the reverse micelle interface.

Deuterium Oxide (D2O) Induces Early Stress Response Gene Expression and Impairs Growth and Metastasis of Experimental Malignant Melanoma

There are two stable isotopes of hydrogen, protium (1H) and deuterium (2H; D). Cellular stress response dysregulation in cancer represents both a major pathological driving force and a promising therapeutic target, but the molecular consequences and potential therapeutic impact of deuterium (2H)-stress on cancer cells remain largely unexplored. We have examined the anti-proliferative and apoptogenic effects of deuterium oxide (D2O; 'heavy water') together with stress response gene expression profiling in panels of malignant melanoma (A375V600E, A375NRAS, G361, LOX-IMVI), and pancreatic ductal adenocarcinoma (PANC-1, Capan-2, or MIA PaCa-2) cells with inclusion of human diploid Hs27 skin fibroblasts. Moreover, we have examined the efficacy of D2O-based pharmacological intervention in murine models of human melanoma tumor growth and metastasis. D2O-induction of apoptosis was substantiated by AV-PI flow cytometry, immunodetection of PARP-1, and pro-caspase 3 cleavage, and rescue by pan-caspase inhibition. Differential array analysis revealed early modulation of stress response gene expression in both A375 melanoma and PANC-1 adenocarcinoma cells elicited by D2O (90%; ≤6 h) (upregulated: CDKN1A, DDIT3, EGR1, GADD45A, HMOX1, NFKBIA, or SOD2 (up to 9-fold; p < 0.01)) confirmed by independent RT-qPCR analysis. Immunoblot analysis revealed rapid onset of D2O-induced stress response phospho-protein activation (p-ERK, p-JNK, p-eIF2α, or p-H2AX) or attenuation (p-AKT). Feasibility of D2O-based chemotherapeutic intervention (drinking water (30% w/w)) was demonstrated in a severe combined immunodeficiency (SCID) mouse melanoma metastasis model using luciferase-expressing A375-Luc2 cells. Lung tumor burden (visualized by bioluminescence imaging) was attenuated by D2O, and inhibition of invasiveness was also confirmed in an in vitro Matrigel transwell invasion assay. D2O supplementation also suppressed tumor growth in a murine xenograft model of human melanoma, and median survival was significantly increased without causing adverse effects. These data demonstrate for the first time that systemic D2O administration impairs growth and metastasis of malignant melanoma through the pharmacological induction of deuterium (2H)-stress.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Unravelling the microhydration frameworks of prototype PAH by infrared spectroscopy: naphthalene–(water)1–3

Microhydration structures of the prototypical PAH, naphthalene, are probed by IR spectroscopy in helium droplets. The sequential water addition produces an extended hydrogen-bonded hydration network bound via π hydrogen bond to the aromatic ring.

The Nature of the Short Oxygen-Oxygen Distance in the Mn4CaO6 Complex of Photosystem II Crystals

The O···O distance for a typical H-bond is ∼2.8 Å, whereas the radiation-damage-free structures of photosystem II (PSII), obtained using the X-ray free electron laser (XFEL), shows remarkably short O···O distances of ∼2 Å in the oxygen-evolving Mn₄CaO_5/6 complex. Herein, we report the protonation/oxidation states of the short O···O atoms in the XFEL structures using a quantum mechanical/molecular mechanical approach. The O5···O6 distance of 1.9 Å is reproduced only when O6 is an unprotonated O radical (O^•-) with Mn(IV)₃Mn(III), i.e., the S₃ state. The potential energy profile shows a barrier-less energy minimum region when O5···O6 = 1.90-2.05 Å (O^•- ↓) or 2.05-2.20 Å (O^•- ↑). Formation of such a short O5···O6 distance is not possible when O6 is OH^- with Mn(IV)₄. In the case in which the O5···O6 distance is 1.9 Å, it seems likely that the O radical species exists in the oxygen-evolving complex of the XFEL-S₃ crystals.

Proton transfer pathway from the oxygen-evolving complex in photosystem II substantiated by extensive mutagenesis

We report a structure-based biological approach to identify the proton-transfer pathway in photosystem II. First, molecular dynamics (MD) simulations were conducted to analyze the H-bond network that may serve as a Grotthuss-like proton conduit. MD simulations show that D1-Asp61, the H-bond acceptor of H₂O at the Mn₄CaO₅ cluster (W1), forms an H-bond via one water molecule with D1-Glu65 but not with D2-Glu312. Then, D1-Asp61, D1-Glu65, D2-Glu312, and the adjacent residues, D1-Arg334, D2-Glu302, and D2-Glu323, were thoroughly mutated to the other 19 residues, i.e., 114 Chlamydomonas chloroplast mutant cells were generated. Mutation of D1-Asp61 was most crucial. Only the D61E and D61C cells grew photoautotrophically and exhibit O₂-evolving activity. Mutations of D2-Glu312 were less crucial to photosynthetic growth than mutations of D1-Glu65. Quantum mechanical/molecular mechanical calculations indicated that in the PSII crystal structure, the proton is predominantly localized at D1-Glu65 along the H-bond with D2-Glu312, i.e., pK_a(D1-Glu65) > pK_a(D2-Glu312). The potential-energy profile shows that the release of the proton from D1-Glu65 leads to the formation of the two short H-bonds between D1-Asp61 and D1-Glu65, which facilitates downhill proton transfer along the Grotthuss-like proton conduit in the S₂ to S₃ transition. It seems possible that D1-Glu65 is involved in the dominant pathway that proceeds from W1 via D1-Asp61 toward the thylakoid lumen, whereas D2-Glu312 and D1-Arg334 may be involved in alternative pathways in some mutants.

Three-Step Spin State Transition and Hysteretic Proton Transfer in the Crystal of an Iron(II) Hydrazone Complex

A proton-electron coupling system, exhibiting unique bistability or multistability of the protonated state, is an attractive target for developing new switchable materials based on proton dynamics. Herein, we present an iron(II) hydrazone crystalline compound, which displays the stepwise transition and bistability of proton transfer at the crystal level. These phenomena are realized through the coupling with spin transition. Although the multi-step transition with hysteresis has been observed in various systems, the corresponding behavior of proton transfer has not been reported in crystalline systems; thus, the described iron(II) complex is the first example. Furthermore, because proton transfer occurs only in one of the two ligands and π electrons redistribute in it, the dipole moment of the iron(II) complexes changes with the proton transfer, wherein the total dipole moment in the crystal was canceled out owing to the antiferroelectric-like arrangement.

Microhydration of protonated biomolecular building blocks: protonated pyrimidine

Protonation and hydration of biomolecules govern their structure, conformation, and function. Herein, we explore the microhydration structure in mass-selected protonated pyrimidine-water clusters (H⁺Pym-W_n, n = 1-4) by a combination of infrared photodissociation spectroscopy (IRPD) between 2450 and 3900 cm^-1 and density functional theory (DFT) calculations at the dispersion-corrected B3LYP-D3/aug-cc-pVTZ level. We further present the IR spectrum of H⁺Pym-N₂ to evaluate the effect of solvent polarity on the intrinsic molecular parameters of H⁺Pym. Our combined spectroscopic and computational approach unequivocally shows that protonation of Pym occurs at one of the two equivalent basic ring N atoms and that the ligands in H⁺Pym-L (L = N₂ or W) preferentially form linear H-bonds to the resulting acidic NH group. Successive addition of water ligands results in the formation of a H-bonded solvent network which increasingly weakens the NH group. Despite substantial activation of the N-H bond upon microhydration, no intracluster proton transfer occurs up to n = 4 because of the balance of relative proton affinities of Pym and W_n and the involved solvation energies. Comparison to neutral Pym-W_n clusters reveals the drastic effects of protonation on microhydration with respect to both structure and interaction strength.

Low Operating Voltage Organic Field-Effect Transistors with Gelatin as a Moisture-Induced Ionic Dielectric Layer: The Issues of High Carrier Mobility

We have developed low-voltage (<2 V) flexible organic field-effect transistors (OFETs) with high carrier mobility using gelatin as a moisture-induced ionic gate dielectric system. Ionic concentration in the gelatin layer depends on the relative humidity condition during the measurement. The capacitance of the dielectric layer used for the calculation of field-effect carrier mobility for the OFETs crucially depends on the frequency at which the capacitance was measured. The results of frequency-dependent gate capacitance together with the anomalous bias-stress effect have been used to determine the exact frequency at which the carrier mobility should be calculated. The observed carrier mobility of the devices is 0.33 cm²/Vs with the capacitance measured at frequency 20 mHz. It can be overestimated to 14 cm²/Vs with the capacitance measured at 100 kHz. The devices can be used as highly sensitive humidity sensors. About three orders of magnitude variation in device current have been observed on the changes in relative humidity (RH) levels from 10 to 80%. The devices show a fast response with a response and recovery times of ∼100 and ∼110 ms, respectively. The devices are flexible up to a 5 mm bending radius.

The study of conformational changes in photosystem II during a charge separation

Photosystem II (PSII) is a multi-subunit pigment-protein complex and is one of several protein assemblies that function cooperatively in photosynthesis in plants and cyanobacteria. As more structural data on PSII become available, new questions arise concerning the nature of the charge separation in PSII reaction center (RC). The crystal structure of PSII RC from cyanobacteria Thermosynechococcus vulcanus was selected for the computational study of conformational changes in photosystem II associated to the charge separation process. The parameterization of cofactors and lipids for classical MD simulation with Amber force field was performed. The parametrized complex of PSII was embedded in the lipid membrane for MD simulation with Amber in Gromacs. The conformational behavior of protein and the cofactors directly involved in the charge separation were studied by MD simulations and QM/MM calculations. This study identified the most likely mechanism of the proton-coupled reduction of plastoquinone Q_B. After the charge separation and the first electron transfer to Q_B, the system undergoes conformational change allowing the first proton transfer to Q_B^- mediated via Ser264. After the second electron transfer to Q_BH, the system again adopts conformation allowing the second proton transfer to Q_BH^-. The reduced Q_BH₂ would then leave the binding pocket.

Thiol based mechanism internalises interacting partners to outer dense fibers in sperm

The sperm tail outer dense fibres (ODFs) contribute passive structural role in sperm motility. The level of disulphide cross-linking of ODFs and their structural thickness determines flagellar bending curvature and motility. During epididymal maturation, proteins are internalized to modify ODF disulphide cross-linking and enable motility. Sperm thiol status is further altered during capacitation in female tract. This suggests that components in female reproductive tract acting on thiol/disulphides could be capable of modulating the tail stiffness to facilitate modulation of the sperm tail rigidity and waveform en route to fertilization. Understanding the biochemical properties and client proteins of ODFs in reproductive tract fluids will help bridge this gap. Using recombinant ODF2 (aka Testis Specific Antigen of 70 kDa) as bait, we identified client proteins in male and female reproductive fluids. A thiol-based interaction and internalization indicates sperm can harness reproductive tract fluids for proteins that interact with ODFs and likely modulate the tail stiffness en route to fertilization.

How the Protein Environment Can Tune the Energy, the Coupling, and the Ultrafast Dynamics of Interacting Chlorophylls: The Example of the Water-Soluble Chlorophyll Protein

The interplay between active molecules and the protein environment in light-harvesting complexes tunes the photophysics and the dynamical properties of pigment-protein complexes in a subtle way, which is not fully understood. Here we characterized the photophysics and the ultrafast dynamics of four variants of the water-soluble chlorophyll protein (WSCP) as an ideal model system to study the behavior of strongly interacting chlorophylls. We found that when coordinated by the WSCP protein, the presence of the formyl group in chlorophyll b replacing the methyl group in chlorophyll a strongly affects the exciton energy and the dynamics of the system, opening up the possibility of tuning the photophysics and the transport properties of multichromophores by engineering specific interactions with the surroundings.

Dissociative photodetachment dynamics of the oxalate monoanion

The dissociative photodetachment (DPD) dynamics of the oxalate monoanion are studied using photoelectron-photofragment coincidence (PPC) spectroscopy. Following photodetachment of C₂O₄H^- at 4.66 eV HOCO + CO₂ products are observed, indicating the facile decarboxylation of the radical driven by the thermodynamic stability of CO₂. No evidence is seen for photodetachment to stable C₂O₄H or ionic photodissociation to produce HOCO^-. Calculations indicate the stabilizing presence of an intramolecular hydrogen bond in the anion via the formation of a strained five-membered ring. No intramolecular hydrogen bond is predicted in the radical due to the lower charge density on the oxygen atom. The PPC spectrum is consistent with a single direct two-body DPD channel that results in fragments of similar mass and is characterized by a large kinetic energy release (KER) and a broad photoelectron spectrum. The large KER is indicative of substantial repulsion in the radical following photodetachment. The form of the photoelectron spectrum is dominated by the bound to continuum Franck-Condon factors (BCFCF) and is suggestive of photodetachment to a repulsive potential energy surface. A lower bound for the electron affinity of C₂O₄H is reported as 4 eV. BCFCF calculations allow an approximate functional form of the repulsive surface along the C-C stretch coordinate to be extracted from the experimental photoelectron spectrum. PPC spectroscopy of the deuterated analogue (C₂O₄D^-), at higher anion beam energies is used to increase the detectability of any possible D atom products, but none are observed.

Catalytic mechanism and evolutionary characteristics of thioredoxin from Halobacterium salinarum NRC-1

Thioredoxin (TRX) is an important antioxidant against oxidative stress. TRX from the extremely halophilic archaeon Halobacterium salinarum NRC-1 (HsTRX-A), which has the highest acidic residue content [(Asp + Glu)/(Arg + Lys + His) = 9.0] among known TRXs, was chosen to elucidate the catalytic mechanism and evolutionary characteristics associated with haloadaptation. X-ray crystallographic analysis revealed that the main-chain structure of HsTRX-A is similar to those of homologous TRXs; for example, the root-mean-square deviations on C^α atoms were <2.3 Å for extant archaeal TRXs and <1.5 Å for resurrected Precambrian TRXs. A unique water network was located near the active-site residues (Cys45 and Cys48) in HsTRX-A, which may enhance the proton transfer required for the reduction of substrates under a high-salt environment. The high density of negative charges on the molecular surface (3.6 × 10^-3 e Å^-2) should improve the solubility and haloadaptivity. Moreover, circular-dichroism measurements and enzymatic assays using a mutant HsTRX-A with deletion of the long flexible N-terminal region (Ala2-Pro17) revealed that Ala2-Pro17 improves the structural stability and the enzymatic activity of HsTRX-A under high-salt environments (>2 M NaCl). The elongation of the N-terminal region in HsTRX-A accompanies the increased hydrophilicity and acidic residue content but does not affect the structure of the active site. These observations offer insights into molecular evolution for haloadaptation and potential applications in halophilic protein-related biotechnology.

Strong proton-shared hydrogen bonding in a methyl imidazole⋯HCl complex: evidence from matrix isolation infrared spectroscopy and ab initio computations

Evidence for proton-shared hydrogen bonding is provided in a methyl imidazole⋯HCl complex using matrix isolation infrared spectroscopy and ab initio computations.

A subtle structural change in the distal haem pocket has a remarkable effect on tuning hydrogen peroxide reactivity in dye decolourising peroxidases fromStreptomyces lividans

A subtle positional shift of the distal haem pocket aspartate in two dye decolourising peroxidase homologs has a remarkable effect on their reactivity with H₂O₂.

Revealing electronic features governing hydrolysis of cephalosporins in the active site of the L1 metallo-β-lactamase

The QM/MM simulations followed by electron density feature analysis are carried out to deepen the understanding of the reaction mechanism of cephalosporin hydrolysis in the active site of the L1 metallo-β-lactamase. The differences in reactivity of ten similar cephalosporin compounds are explained by using an extended set of bonding descriptors. The limiting step of the reaction is characterized by the proton transfer to the nitrogen atom of the cephalosporin thiazine ring accompanied with formation of the C₄ Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> C₃ double bond in its N–C₄–C₃ fragment. The temporary N⋯H–O_w hydrogen bond, which is formed in the transition state of the limiting step of the reaction was recognized as a key atomic interaction governing the reactivity of various cephalosporins. Non-local real-space bonding descriptors show that different extent of localization of electron lone pair at N atom in the transition state affect the reactivity of compounds: smaller electron localization is typical for the less reactive species. In particular, the Fermi hole analysis shows how exchange electron correlation in the N⋯H–O_w fragment control electron lone pair localization. Delocalization tensor, linear response kernel and source function indicate that features of electron delocalization in the N–C₄–C₃ fragment of cephalosporins in the transition state complexes determine the differences in C₄–C₃ bond for substrates with high and low rate constants. The C₄–C₃ bond of the N–C₄–C₃ fragment at the transition state is similar to that of the preceding intermediate for the less reactive species and resembles the features of the enzyme–product complex for more reactive compounds. The power and limitations of the descriptors applied for solving the problem are discussed and the generality of approach is stressed.

Combination of QM/MM and modern bonding descriptors explains different reactivity of cephalosporins in the active site of the L1 metallo-β-lactamase.

X-ray crystallographic studies on the hydrogen isotope effects of green fluorescent protein at sub-ångström resolutions

Hydrogen atoms are critical to the nature and properties of proteins, and thus deuteration has the potential to influence protein function. In fact, it has been reported that some deuterated proteins show different physical and chemical properties to their protiated counterparts. Consequently, it is important to investigate protonation states around the active site when using deuterated proteins. Here, hydrogen isotope effects on the S65T/F99S/M153T/V163A variant of green fluorescent protein (GFP), in which the deprotonated B form is dominant at pH 8.5, were investigated. The pH/pD dependence of the absorption and fluorescence spectra indicates that the protonation state of the chromophore is the same in protiated GFP in H₂O and protiated GFP in D₂O at pH/pD 8.5, while the pK_a of the chromophore became higher in D₂O. Indeed, X-ray crystallographic analyses at sub-ångström resolution revealed no apparent changes in the protonation state of the chromophore between the two samples. However, detailed comparisons of the hydrogen OMIT maps revealed that the protonation state of His148 in the vicinity of the chromophore differed between the two samples. This indicates that protonation states around the active site should be carefully adjusted to be the same as those of the protiated protein when neutron crystallographic analyses of proteins are performed.

Possible Mechanisms of Biological Effects Observed in Living Systems during 2H/1H Isotope Fractionation and Deuterium Interactions with Other Biogenic Isotopes

This article presents the original descriptions of some recent physics mechanisms (based on the thermodynamic, kinetic, and quantum tunnel effects) providing stable 2H/1H isotope fractionation, leading to the accumulation of particular isotopic forms in intra- or intercellular space, including the molecular effects of deuterium interaction with 18O/17O/16O, 15N/14N, 13C/12C, and other stable biogenic isotopes. These effects were observed mainly at the organelle (mitochondria) and cell levels. A new hypothesis for heavy nonradioactive isotope fractionation in living systems via neutron effect realization is discussed. The comparative analysis of some experimental studies results revealed the following observation: "Isotopic shock" is highly probable and is observed mostly when chemical bonds form between atoms with a summary odd number of neutrons (i.e., bonds with a non-compensated neutron, which correspond to the following equation: Nn - Np = 2k + 1, where k ϵ Z, k is the integer, Z is the set of non-negative integers, Nn is number of neutrons, and Np is number of protons of each individual atom, or in pair of isotopes with a chemical bond). Data on the efficacy and metabolic pathways of the therapy also considered 2H-modified drinking and diet for some diseases, such as Alzheimer's disease, Friedreich's ataxia, mitochondrial disorders, diabetes, cerebral hypoxia, Parkinson's disease, and brain cancer.

How [FeFe]-Hydrogenase Facilitates Bidirectional Proton Transfer

Hydrogenases are metalloenzymes that catalyze the conversion of protons and molecular hydrogen, H₂. [FeFe]-hydrogenases show particularly high rates of hydrogen turnover and have inspired numerous compounds for biomimetic H₂ production. Two decades of research on the active site cofactor of [FeFe]-hydrogenases have put forward multiple models of the catalytic proceedings. In comparison, our understanding of proton transfer is poor. Previously, residues were identified forming a hydrogen-bonding network between active site cofactor and bulk solvent; however, the exact mechanism of catalytic proton transfer remained inconclusive. Here, we employ in situ infrared difference spectroscopy on the [FeFe]-hydrogenase from Chlamydomonas reinhardtii evaluating dynamic changes in the hydrogen-bonding network upon photoreduction. While proton transfer appears to be impaired in the oxidized state (Hox), the presented data support continuous proton transfer in the reduced state (Hred). Our analysis allows for a direct, molecular unique assignment to individual amino acid residues. We found that transient protonation changes of glutamic acid residue E141 and, most notably, arginine R148 facilitate bidirectional proton transfer in [FeFe]-hydrogenases.

Does Oxygen Feature Chalcogen Bonding?

Using the second-order Møller-Plesset perturbation theory (MP2), together with Dunning's all-electron correlation consistent basis set aug-cc-pVTZ, we show that the covalently bound oxygen atom present in a series of 21 prototypical monomer molecules examined does conceive a positive (or a negative) σ-hole. A σ-hole, in general, is an electron density-deficient region on a bound atom M along the outer extension of the R-M covalent bond, where R is the reminder part of the molecule, and M is the main group atom covalently bonded to R. We have also examined some exemplar 1:1 binary complexes that are formed between five randomly chosen monomers of the above series and the nitrogen- and oxygen-containing Lewis bases in N2, PN, NH3, and OH2. We show that the O-centered positive σ-hole in the selected monomers has the ability to form the chalcogen bonding interaction, and this is when the σ-hole on O is placed in the close proximity of the negative site in the partner molecule. Although the interaction energy and the various other 12 characteristics revealed from this study indicate the presence of any weakly bound interaction between the monomers in the six complexes, our result is strongly inconsistent with the general view that oxygen does not form a chalcogen-bonded interaction.

Proton Transfer Pathways between Active Sites and Proximal Clusters in the Membrane-Bound [NiFe] Hydrogenase

[NiFe] hydrogenases are enzymes that catalyze the splitting of molecular hydrogen according to the reaction H₂ → 2H⁺ + 2e^-. Most of these enzymes are inhibited even by low traces of O₂. However, a special group of O₂-tolerant hydrogenases exists. A member of this group is the membrane-bound [NiFe] hydrogenase from Ralstonia eutropha ( ReMBH). The ReMBH harbors an unusual iron sulfur cluster with composition 4Fe3S(6Cys) that is able to undergo structural changes triggering the flow of two electrons to the [NiFe] active site. These electrons promote oxygen reduction at the active site, preventing, in this way, aerobic inactivation of the enzyme. In the superoxidized state, the [4Fe3S] cluster binds to a hydroxyl group that originates from either molecular oxygen or water reaching the site. Both reactions, oxygen reduction to water at the [NiFe]- or [4Fe3S]-centers and oxygen evolution from water at the proximal cluster, require the delivery of protons regulated by a subtle communication mechanism between these metal centers. In this work, we sequentially apply multiscale modeling techniques as quantum mechanical/molecular mechanics methods and classical molecular dynamics simulations to investigate the role of two distinct proton transfer pathways connecting the [NiFe] active site and the [4Fe3S] proximal cluster of ReMBH in the protection mechanism against an oxygen attack. Although the "glutamate" pathway is preferred by protons migrating toward the active site to avoid inactivation by O₂, the "histidine" pathway plays an essential role in delivering protons for O₂ reduction at the proximal cluster. The results obtained in this work not only provide new pieces to the puzzling catalytic mechanisms governing O₂-tolerant hydrogenases but also highlight the relevance of dynamics in the proper description of biochemical reactions in general.

Catalysis by Networks of Cooperative Hydrogen Bonds

The main paradigm of today's chemistry is sustainability. In pursuing sustainability, we need to learn from chemical processes carried out by Nature and realize that Nature does not use either strong acids, or strong bases or fancy reagents to achieve outstanding chemical processes. Instead, enzyme activity leans on the cooperation of several chemical entities to avoid strong acids or bases or to achieve such an apparently simple goal as transferring a proton from an NuH unit to an E unit (NuH + E → Nu–EH). Hydrogen bond catalysis emerged strongly two decades ago in trying to imitate Nature and avoid metal catalysis. Now to mount another step in pursuing the goal of sustainability, the focus is upon cooperativity between the different players involved in catalysis. This chapter looks at the concept of cooperativity and, more specifically, (a) examines the role of cooperative hydrogen bonded arrays of the general type NuH⋯(NuH)n⋯NuH (i.e. intermolecular cooperativity) to facilitate general acid–base catalysis, not only in the solution phase but also under solvent-free and catalyst-free conditions, and, most important, (b) analyzes the capacity of designer chiral organocatalysts displaying intramolecular networks of cooperative hydrogen bonds (NCHBs) to facilitate enantioselective synthesis by bringing conformational rigidity to the catalyst in addition to simultaneously increasing the acidity of key hydrogen atoms so to achieve better complementarity in the highly polarized transition states.