Protein Environment Research Articles

Environment-dependent chlorophyll-chlorophyll charge transfer states in Lhca4 pigment-protein complex.

Photosystem I (PSI) light-harvesting antenna complexes LHCI contain spectral forms that absorb and emit photons of lower energy than that of its primary electron donor, P700. The most red-shifted fluorescence is associated with the Lhca4 complex. It has been suggested that this red emission is related to the inter-chlorophyll charge transfer (CT) states. In this work we present a systematic quantum-chemical study of the CT states in Lhca4, accounting for the influence of the protein environment by estimating the electrostatic interactions. We show that significant energy shifts result from these interactions and propose that the emission of the Lhca4 complex is related not only to the previously proposed a603+-a608- state, but also to the a602+-a603- state. We also investigate how different protonation patterns of protein amino acids affect the energetics of the CT states.

A Comprehensive Study of Factors Affecting the Prediction of the pKa Shift of Asp26 in Thioredoxin Protein.

The stable protonation state of ionizable amino acids in a protein can be predicted by computing the pKa shift of that residue within the protein environment. Thermodynamic Integration (TI) is an ideal molecular dynamics-based approach for predicting the pKa shift of ionizable protein residues. Here, we probe TI-based simulation protocols for their ability to accurately predict the pKa shift of Asp26 in thioredoxin. While implicit solvent models can predict the pKa shift accurately, explicit solvent models result in substantial errors. To understand the underlying reason for this surprising discrepancy, we investigate the role of various factors such as solvent models, conformational sampling, background charges, and polarization.

A Robust Induced Fit Docking Approach with the Combination of the Hybrid All-Atom/United-Atom/Coarse-Grained Model and Simulated Annealing.

Molecular docking remains an indispensable tool in computational biology and structure-based drug discovery. However, the correct prediction of binding poses remains a major challenge for molecular docking, especially for target proteins where a substrate binding induces significant reorganization of the active site. Here, we introduce an Induced Fit Docking (IFD) approach named AA/UA/CG-SA-IFD, which combines a hybrid All-Atom/United-Atom/Coarse-Grained model with Simulated Annealing. In this approach, the core region is represented by the All-Atom(AA) model, while the protein environment beyond the core region and the solvent are treated with either the United-Atom (UA) or the Coarse-Grained (CG) model. By combining the Elastic Network Model (ENM) for the CG region, the hybrid model ensures a reasonable description of ligand binding and the environmental effects of the protein, facilitating highly efficient and reliable sampling of ligand binding through Simulated Annealing (SA) at a high temperature. Upon validation with two testing sets, the AA/UA/CG-SA-IFD approach demonstrates remarkable accuracy and efficiency in induced fit docking, even for challenging cases where the docked poses significantly deviate from crystal structures.

Capturing the blue-light activated state of the Phot-LOV1 domain from Chlamydomonas reinhardtii using time-resolved serial synchrotron crystallography.

Light-oxygen-voltage (LOV) domains are small photosensory flavoprotein modules that allow the conversion of external stimuli (sunlight) into intracellular signals responsible for various cell behaviors (e.g. phototropism and chloroplast relocation). This ability relies on the light-induced formation of a covalent thioether adduct between a flavin chromophore and a reactive cysteine from the protein environment, which triggers a cascade of structural changes that result in the activation of a serine/threonine (Ser/Thr) kinase. Recent developments in time-resolved crystallography may allow the activation cascade of the LOV domain to be observed in real time, which has been elusive. In this study, we report a robust protocol for the production and stable delivery of microcrystals of the LOV domain of phototropin Phot-1 from Chlamydomonas reinhardtii (CrPhotLOV1) with a high-viscosity injector for time-resolved serial synchrotron crystallography (TR-SSX). The detailed process covers all aspects, from sample optimization to data collection, which may serve as a guide for soluble protein preparation for TR-SSX. In addition, we show that the crystals obtained preserve the photoreactivity using infrared spectroscopy. Furthermore, the results of the TR-SSX experiment provide high-resolution insights into structural alterations of CrPhotLOV1 from Δt = 2.5 ms up to Δt = 95 ms post-photoactivation, including resolving the geometry of the thioether adduct and the C-terminal region implicated in the signal transduction process.

Proton-Coupled Electron Transfer upon Oxidation of Tyrosine in a De Novo Protein: Analysis of Proton Acceptor Candidates.

Redox-active residues, such as tyrosine and tryptophan, play important roles in a wide range of biological processes. The α3Y de novo protein, which is composed of three α helices and a tyrosine residue Y32, provides a platform for investigating the redox properties of tyrosine in a well-defined protein environment. Herein, the proton-coupled electron transfer (PCET) reaction that occurs upon oxidation of tyrosine in this model protein by a ruthenium photosensitizer is studied by using a vibronically nonadiabatic PCET theory that includes hydrogen tunneling and excited vibronic states. The input quantities to the analytical nonadiabatic rate constant expression, such as the diabatic proton potential energy curves and associated proton vibrational wave functions, reorganization energy, and proton donor-acceptor distribution functions, are obtained from density functional theory calculations on model systems and molecular dynamics simulations of the solvated α3Y protein. Two possible proton acceptors, namely, water or a glutamate residue in the protein scaffold, are explored. The PCET rate constant is greater when glutamate is the proton acceptor, mainly due to the more favorable driving force and shorter equilibrium proton donor-acceptor distance, although contributions from excited vibronic states mitigate these effects. Nevertheless, water could be the dominant proton acceptor if its equilibrium constant associated with hydrogen bond formation is significantly greater than that for glutamate. Although these calculations do not definitively identify the proton acceptor for this PCET reaction, they elucidate the conditions under which each proton acceptor can be favored. These insights have implications for tyrosine-based PCET in a wide variety of biochemical processes.

Real-time observation of a metal complex-driven reaction intermediate using a porous protein crystal and serial femtosecond crystallography

Determining short-lived intermediate structures in chemical reactions is challenging. Although ultrafast spectroscopic methods can detect the formation of transient intermediates, real-space structures cannot be determined directly from such studies. Time-resolved serial femtosecond crystallography (TR-SFX) has recently proven to be a powerful method for capturing molecular changes in proteins on femtosecond timescales. However, the methodology has been mostly applied to natural proteins/enzymes and limited to reactions promoted by synthetic molecules due to structure determination challenges. This work demonstrates the applicability of TR-SFX for investigations of chemical reaction mechanisms of synthetic metal complexes. We fix a light-induced CO-releasing Mn(CO)3 reaction center in porous hen egg white lysozyme (HEWL) microcrystals. By controlling light exposure and time, we capture the real-time formation of Mn-carbonyl intermediates during the CO release reaction. The asymmetric protein environment is found to influence the order of CO release. The experimentally-observed reaction path agrees with quantum mechanical calculations. Therefore, our demonstration offers a new approach to visualize atomic-level reactions of small molecules using TR-SFX with real-space structure determination. This advance holds the potential to facilitate design of artificial metalloenzymes with precise mechanisms, empowering design, control and development of innovative reactions.

Presence of low-energy chlorophylls d in photosystem I trimer and monomer cores isolated from Acaryochloris sp. NBRC 102871.

Acaryochloris species belong to a special category of cyanobacteria possessing chlorophyll (Chl) d. One of the photosynthetic characteristics of Acaryochloris marina MBIC11017 is that the absorption spectra of photosystem I (PSI) showed almost no bands and shoulders of low-energy Chls d over 740nm. In contrast, the absorption spectra of other Acaryochloris species showed a shoulder around 740nm, suggesting that low-energy Chls d within PSI are diversified among Acaryochloris species. In this study, we purified PSI trimer and monomer cores from Acaryochloris sp. NBRC 102871 and examined their protein and pigment compositions and spectral properties. The protein bands and pigment compositions of the PSI trimer and monomer of NBRC102871 were virtually identical to those of MBIC11017. The absorption spectra of the NBRC102871 PSIs exhibited a shoulder around 740nm, whereas the fluorescence spectra of PSI trimer and monomer displayed maximum peaks at 754 and 767nm, respectively. These spectral properties were different from those of MBIC11017, indicating the presence of low-energy Chls d within the NBRC102871 PSIs. Moreover, we analyzed the NBRC102871 genome to identify amino acid sequences of PSI proteins and compared them with those of the A. marina MBIC11017 and MBIC10699 strains whose genomes are available. The results showed that some of the sequences in NBRC102871 were distinct from those in MBIC11017 and MBIC10699. These findings provide insights into the variety of low-energy Chls d with respect to the protein environments of PSI cores among the three Acaryochloris strains.

Photobiomodulation of Na,K‐ATPase in native membrane fraction and reconstituted in DPPC:DPPE‐liposome

AbstractStudies focusing on how photobiomodulation (PBM) can affect the structure and function of proteins are scarce in the literature. Few previous studies have shown that the enzymatic activity of Na,K‐ATPAse (NKA) can be photo‐modulated. However, the variability of sample preparation and light irradiation wavelengths have not allowed for an unequivocal conclusion about the PBM of NKA. Here, we investigate minimal membrane models containing NKA, namely, native membrane fraction and DPPC:DPPE proteoliposome upon laser irradiation at wavelengths 532, 650, and 780 nm. Interestingly, we show that the PBM on the NKA enzymatic activity has a bell‐shaped profile with a stimulation peak (~15% increase) at around 20 J.cm−2 and 6 J.cm−2 for the membrane‐bound and the proteoliposome samples, respectively, and are practically wavelength independent. Further, by normalizing the enzymatic activity by the NKA enzyme concentration, we show that the PBM response is related to the protein amount with small influence due to protein's environment. The stimulation decays over time reaching the basal level around 6 h after the irradiation for the three lasers and both NKA samples. Our results demonstrate the potential of using low‐level laser therapy to modulate NKA activity, which may have therapeutic implications and benefits.

The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications.

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

Fermi Resonance of the N-D Stretching Mode Probing the Local Hydrogen-Bonding Environment in Proteins.

Local H-bonding interactions are crucial for proteins to undergo various structural transitions and form different secondary structures. However, identifying slight distinctions in the local H-bonding of proteins is rather challenging. Here, we demonstrate that the Fermi resonance of the N-D stretching mode can provide an effective probe for the localized H-bonding environment of proteins both at the surface/interface and in the bulk. Using sum frequency generation vibrational spectroscopy and infrared spectroscopy, we established a correlation between the Fermi resonance of the N-D mode and protein secondary structures. The H-bond of N-D···C═O splits the N-D modes into two peaks (∼2410 and ∼2470 cm-1). The relative strength ratio (R) between the ∼2410 cm-1 peak and the ∼2470 cm-1 peak is very sensitive to H-bond strength and protein secondary structure. R is less than 1 for α-helical peptides, while R is greater than 1 for β-sheet peptides. For R < 2.5, both α-helical/loop structures and β-sheet structures exhibit almost identical Fermi coupling strengths (W = 28 cm-1). For R > 2.5, W decreases from 28 to 14 cm-1 and depends on the aggregation degree of the β-sheet oligomers or fibrils. The initial local H-bonding status impacts the misfolding dynamics of proteins at the lipid bilayer interface.

THEORETICAL DESCRIPTION OF PUCKER STATES OF CYCLIC GROUPS OF CAROTENOIDS

Carotenoids are one of the largest classes among biological pigments. They are found both in free form (dissolved in hydrophobic environments), and bound in stoichiometric complexes with protein (carotenoproteins). The interaction between the carotenoid and the protein environment is determined by the molecular conformation and dynamic behavior of the cofactor. In this work, using molecular modeling methods, we analyzed the entire set of alternative conformational states of cyclic carotenoids. The result of the theoretical study is an energysorted catalog of alternative conformational states for biogenic cyclic carotenoids and their optical isomers. The conformational states characterized in this study can be used to correctly select initial conditions when solving problems of molecular modeling of carotenoproteins

THEORETICAL DESCRIPTION OF PUCKER STATES OF CYCLIC GROUPS OF CAROTENOIDS

Carotenoids are one of the largest classes among biological pigments. They are found both in free form (dissolved in hydrophobic environments), and bound in stoichiometric complexes with protein (carotenoproteins). The interaction between the carotenoid and the protein environment is determined by the molecular conformation and dynamic behavior of the cofactor. In this work, using molecular modeling methods, we analyzed the entire set of alternative conformational states of cyclic carotenoids. The result of the theoretical study is an energysorted catalog of alternative conformational states for biogenic cyclic carotenoids and their optical isomers. The conformational states characterized in this study can be used to correctly select initial conditions when solving problems of molecular modeling of carotenoproteins

Design, Synthesis and Catalytic Activity of Protein Containing Thiotyrosine as an Active Site Residue.

Native chemical ligation is a key reaction in the toolbox of chemical methods for the synthesis of native and modified proteins. The catalysis of ligation is commonly performed by using small aryl-thiol molecules added at high concentrations. In this work, we incorporated thiotyrosine, a non-canonical amino acid containing an aryl-thiol moiety, into a designed cyclic protein « sans queue ni tête ». Importantly, the protein environment reduced the pKa of the thiol group to 5.8-5.9, which is significantly lower than the previously reported value for thiotyrosine in a short peptide (pKa 6.4). Furthermore, we demonstrated the catalytic activity of this protein both as hydrolase and in native chemical ligation of peptides. These results will be useful for the development of efficient protein catalysts (enzymes) for protein synthesis and modification.

Protein Effects on the Excitation Energies and Exciton Dynamics of the CP24 Antenna Complex.

In this study, the site energy fluctuations, energy transfer dynamics, and some spectroscopic properties of the minor light-harvesting complex CP24 in a membrane environment were determined. For this purpose, a 3 μs-long classical molecular dynamics simulation was performed for the CP24 complex. Furthermore, using the density functional tight binding/molecular mechanics molecular dynamics (DFTB/MM MD) approach, we performed excited state calculations for the chlorophyll a and chlorophyll b molecules in the complex starting from five different positions of the MD trajectory. During the extended simulations, we observed variations in the site energies of the different sets as a result of the fluctuating protein environment. In particular, a water coordination to Chl-b 608 occurred only after about 1 μs in the simulations, demonstrating dynamic changes in the environment of this pigment. From the classical and the DFTB/MM MD simulations, spectral densities and the (time-dependent) Hamiltonian of the complex were determined. Based on these results, three independent strongly coupled chlorophyll clusters were revealed within the complex. In addition, absorption and fluorescence spectra were determined together with the exciton relaxation dynamics, which reasonably well agrees with experimental time scales.

Unraveling Electronic and Vibrational Coherences Following a Charge Transfer Process in a Photosystem II Reaction Center

A reaction center is a unique biological system that performs the initial charge separation within a Photosystem II (PSII) multiunit enzyme, which eventually drives the catalytic water-splitting in plants and algae. The possible role of quantum coherences coinciding with the energy and charge transfer processes in PSII reaction center is one of the active areas of research. Here, we study these quantum coherences by using a numerically exact method on an excitonic dimer model, including linear vibronic coupling and employing optimal parameters from experimental two-dimensional coherent spectroscopic measurements. This enables us to precisely capture the excitonic interaction between pigments and the dissipation of the energy from electronic and charge-transfer (CT) states to the protein environment. We employ the time nonlocal (TNL) quantum master equation to calculate the population dynamics, which yields numerically reliable results. The calculated results show that, due to the strong dissipation, the lifetime of electronic coherence is too short to have direct participation in the charge transfer processes. However, there are long-lived vibrational coherences present in the system at frequencies close to the excitionic energy gap. These are strongly coupled with the electronic coherences, which makes the detection of the electronic coherences with conventional techniques very challenging. Additionally, we unravel the strong excitonic interaction of radical pair (PD1 and PD2) in the reaction center, which results in a long-lived electronic coherence of &gt;100 fs, even at room temperature. Our work provide important physical insight to the charge separation process in PSII reaction center, which may be helpful for better understanding of photophysical processes in other natural and artificial light-harvesting systems.

Raman Vibrational Signatures of Excited States of Echinenone in the Orange Carotenoid Protein (OCP) and Implications for its Photoactivation Mechanism

In this study, the vibrational characteristics of optically excited echinenone in various solvents and the Orange Carotenoid Protein (OCP) in red and orange states are systematically investigated through steady-state and time-resolved spectroscopy techniques. Time-resolved experiments, employing both Transient Absorption (TA) and Femtosecond Stimulated Raman Spectroscopy (FSRS), reveal different states in the OCP photoactivation process. The time-resolved studies indicate vibrational signatures of exited states positioned above the S1 state during the initial 140 fs of carotenoid evolution in OCP, an absence of a vibrational signature for the relaxed S1 state of echinenone in OCP, and more robust signatures of a highly excited ground state (GS) in OCP. Differences in S1 state vibration population signatures between OCP and solvents are attributed to distinct conformations of echinenone in OCP and hydrogen bonds at the keto group forming a short-lived intramolecular charge transfer (ICT) state. The vibrational dynamics of the hot GS in OCP show a more pronounced red shift of ground state CC vibration compared to echinenone in solvents, thus suggesting an unusually hot form of GS. The study proposes a hypothesis for the photoactivation mechanism of OCP, emphasizing the high level of vibrational excitation in longitudinal stretching modes as a driving force. In conclusion, the comparison of vibrational signatures reveals unique dynamics of energy dissipation in OCP, providing insights into the photoactivation mechanism and highlighting the impact of the protein environment on carotenoid behavior. The study underscores the importance of vibrational analysis in understanding the intricate processes involved in early phase OCP photoactivation.

Brain Delivery of Protein Therapeutics by Cell Matrix-Inspired Biomimetic Nanocarrier.

Protein therapeutics are anticipated to offer significant treatment options for central nervous system (CNS) diseases. However, the majority of proteins are unable to traverse the blood-brain barrier (BBB) and reach their CNS target sites. Inspired by the natural environment of active proteins, the cell matrix components hyaluronic acid (HA) and protamine (PRTM) are used to self-assemble with proteins to form a protein-loaded biomimetic core and then incorporated into ApoE3-reconstituted high-density lipoprotein (rHDL) to form a protein-loaded biomimetic nanocarrier (Protein-HA-PRTM-rHDL). This cell matrix-inspired biomimetic nanocarrier facilitates the penetration of protein therapeutics across the BBB and enables their access to intracellular target sites. Specifically, CAT-HA-PRTM-rHDL facilitates rapid intracellular delivery and release of catalase (CAT) via macropinocytosis-activated membrane fusion, resulting in improved spatial learning and memory in traumatic brain injury (TBI) model mice (significantly reduces the latency of TBI mice and doubles the number of crossing platforms), and enhances motor function and prolongs survival in amyotrophic lateral sclerosis (ALS) model mice (extended the median survival of ALS mice by more than 10 days). Collectively, this cell matrix-inspired nanoplatform enables the efficient CNS delivery of protein therapeutics and provides a novel approach for the treatment of CNS diseases.

Delineating three-dimensional behavior of uveal melanoma cells under anchorage independent or dependent conditions

BackgroundAlthough rare, uveal melanoma (UM) is a life-threatening malignancy. Understanding its biology is necessary to improve disease outcome. Three-dimensional (3D) in vitro culture methods have emerged as tools that incorporate physical and spatial cues that better mimic tumor biology and in turn deliver more predictive preclinical data. Herein, we comprehensively characterize UM cells under different 3D culture settings as a suitable model to study tumor cell behavior and therapeutic intervention.MethodsSix UM cell lines were tested in two-dimensional (2D) and 3D-culture conditions. For 3D cultures, we used anchorage-dependent (AD) methods where cells were embedded or seeded on top of basement membrane extracts and anchorage-free (AF) methods where cells were seeded on agarose pre-coated plates, ultra-low attachment plates, and on hanging drops, with or without methylcellulose. Cultures were analyzed for multicellular tumor structures (MCTs) development by phase contrast and confocal imaging, and cell wellbeing was assessed based on viability, membrane integrity, vitality, apoptotic features, and DNA synthesis. Vascular endothelial growth factor (VEGF) production was evaluated under hypoxic conditions for cell function analysis.ResultsUM cells cultured following anchorage-free methods developed MCTs shaped as spheres. Regardless of their sizes and degree of compaction, these spheres displayed an outer ring of viable and proliferating cells, and a core with less proliferating and apoptotic cells. In contrast, UM cells maintained under anchorage-dependent conditions established several morphological adaptations. Some remained isolated and rounded, formed multi-size irregular aggregates, or adopted a 2D-like flat appearance. These cells invariably conserved their metabolic activity and conserved melanocytic markers (i.e., expression of Melan A/Mart-1 and HMB45). Notably, under hypoxia, cells maintained under 3D conditions secrete more VEGF compared to cells cultured under 2D conditions.ConclusionsUnder an anchorage-free environment, UM cells form sphere-like MCTs that acquire attributes reminiscent of abnormal vascularized solid tumors. UM cells behavior in anchorage-dependent manner exposed diverse cells populations in response to cues from an enriched extracellular matrix proteins (ECM) environment, highlighting the plasticity of UM cells. This study provides a 3D cell culture platform that is more predictive of the biology of UM. The integration of such platforms to explore mechanisms of ECM-mediated tumor resistance, metastatic abilities, and to test novel therapeutics (i.e., anti-angiogenics and immunomodulators) would benefit UM care.

(2S,4S)-5-Fluoroleucine, (2S,4R)-5-Fluoroleucine, and 5,5'-Difluoroleucine in Escherichia coli PpiB: Protein Production, 19F NMR, and Ligand Sensing Enhanced by the γ-Gauche Effect.

Global substitution of leucine for analogues containing CH2F instead of methyl groups delivers proteins with multiple sites for monitoring by 19F nuclear magnetic resonance (NMR) spectroscopy. The 19 kDa Escherichia coli peptidyl-prolyl cis-trans isomerase B (PpiB) was prepared with uniform high-level substitution of leucine by (2S,4S)-5-fluoroleucine, (2S,4R)-5-fluoroleucine, or 5,5'-difluoroleucine. The stability of the samples toward thermal denaturation was little altered compared to the wild-type protein. 19F nuclear magnetic resonance (NMR) spectra showed large chemical shift dispersions between 6 and 17 ppm. The 19F chemical shifts correlate with the three-bond 1H-19F couplings (3JHF), providing the first experimental verification of the γ-gauche effect predicted by [Feeney, J. J. Am. Chem. Soc. 1996, 118, 8700-8706] and establishing the effect as the predominant determinant of the 19F chemical shifts of CH2F groups. Individual CH2F groups can be confined to single rotameric states by the protein environment, but most CH2F groups exchange between different rotamers at a rate that is fast on the NMR chemical shift scale. Interactions between fluorine atoms in 5,5'-difluoroleucine bias the CH2F rotamers in agreement with results obtained previously for 1,3-difluoropropane. The sensitivity of the 19F chemical shift to the rotameric state of the CH2F groups potentially renders them particularly sensitive for detecting allosteric effects.

Energetics of the H-Bond Network in Exiguobacterium sibiricum Rhodopsin.

Exiguobacterium sibiricum rhodopsin (ESR) functions as a light-driven proton pump utilizing Lys96 for proton uptake and maintaining its activity over a wide pH range. Using a combination of methodologies including the linear Poisson-Boltzmann equation and a quantum mechanical/molecular mechanical approach with a polarizable continuum model, we explore the microscopic mechanisms underlying its pumping activity. Lys96, the primary proton uptake site, remains deprotonated owing to the loss of solvation in the ESR protein environment. Asp85, serving as a proton acceptor group for Lys96, does not form a low-barrier H-bond with His57. Instead, deprotonated Asp85 forms a salt-bridge with protonated His57, and the proton is predominantly located at the His57 moiety. Glu214, the only acidic residue at the end of the H-bond network exhibits a pKa value of ∼6, slightly elevated due to solvation loss. It seems likely that the H-bond network [Asp85···His57···H2O···Glu214] serves as a proton-conducting pathway toward the protein bulk surface.