Spectral signatures of the pentagonal water cluster in bacteriorhodopsin

Nuclear Magnetic Resonance Detection of Hydrogen Bond Network in a Proton Pump Rhodopsin RxR and Its Alteration during the Cyclic Photoreaction

Hydrogen bond formation and deformation are crucial for the structural construction and functional expression of biomolecules. However, direct observation of exchangeable hydrogens, especially for oxygen-bound hydrogens, relevant to hydrogen bonds is challenging for current structural analysis approaches. Using solution-state NMR spectroscopy, this study detected the functionally important exchangeable hydrogens (i.e., Y49-ηOH and Y178-ηOH) involved in the pentagonal hydrogen bond network in the active site of R. xylanophilus rhodopsin (RxR), which functions as a light-driven proton pump. Moreover, utilization of the original light-irradiation NMR approach allowed us to detect and characterize the late photointermediate state (i.e., O-state) of RxR and revealed that hydrogen bonds relevant to Y49 and Y178 are still maintained during the photointermediate state. In contrast, the hydrogen bond between W75-εNH and D205-γCOO^- is strengthened and stabilizes the O-state.

Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin

An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water

The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.

pH-sensitive vibrational probe reveals a cytoplasmic protonated cluster in bacteriorhodopsin

Infrared spectroscopy has been used in the past to probe the dynamics of internal proton transfer reactions taking place during the functional mechanism of proteins but has remained mostly silent to protonation changes in the aqueous medium. Here, by selectively monitoring vibrational changes of buffer molecules with a temporal resolution of 6 µs, we have traced proton release and uptake events in the light-driven proton-pump bacteriorhodopsin and correlate these to other molecular processes within the protein. We demonstrate that two distinct chemical entities contribute to the temporal evolution and spectral shape of the continuum band, an unusually broad band extending from 2,300 to well below 1,700 cm^-1 The first contribution corresponds to deprotonation of the proton release complex (PRC), a complex in the extracellular domain of bacteriorhodopsin where an excess proton is shared by a cluster of internal water molecules and/or ionic E194/E204 carboxylic groups. We assign the second component of the continuum band to the proton uptake complex, a cluster with an excess proton reminiscent to the PRC but located in the cytoplasmic domain and possibly stabilized by D38. Our findings refine the current interpretation of the continuum band and call for a reevaluation of the last proton transfer steps in bacteriorhodopsin.

A delocalized proton-binding site within a membrane protein

The role of protein-bound water molecules in protein function and catalysis is an emerging topic. Here, we studied the solvation of an excess proton by protein-bound water molecules and the contribution of the surrounding amino acid residues at the proton release site of the membrane protein bacteriorhodopsin. It hosts an excess proton within a protein-bound water cluster, which is hydrogen bonded to several surrounding amino acids. Indicative of delocalization is a broad continuum absorbance experimentally observed by time-resolved Fourier transform infrared spectroscopy. In combination with site-directed mutagenesis, the involvement of several amino acids (especially Glu-194 and Glu-204) in the delocalization was elaborated. Details regarding the contributions of the glutamates and water molecules to the delocalization mode in biomolecular simulations are controversial. We carried out quantum mechanics/molecular mechanics (QM/MM) self-consistent charge density functional tight-binding simulations for all amino acids that have been experimentally shown to be involved in solvation of the excess proton, and systematically investigated the influence of the quantum box size. We compared calculated theoretical infrared spectra with experimental ones as a measure for the correct description of excess proton delocalization. A continuum absorbance can only be observed for small quantum boxes containing few amino acids and/or water molecules. Larger quantum boxes, including all experimentally shown involved amino acids, resulted in narrow absorbance bands, indicating protonation of a single binding site in contradiction to experimental results. We conclude that small quantum boxes seem to reproduce representative extreme cases of proton delocalization modes: proton delocalization only on water molecules or only between Glu-194 and Glu-204. Extending the experimental spectral region to lower wave numbers, a water-delocalized proton reproduces the observed continuum absorbance better than a glutamate-shared delocalized proton. However, a full agreement between QM simulations and experimental results on the delocalized excess proton will require a larger quantum box as well as more sophisticated QM/MM methods.

Changes in the hydrogen-bonding strength of internal water molecules and cysteine residues in the conductive state of channelrhodopsin-1

Water plays an essential role in the structure and function of proteins, particularly in the less understood class of membrane proteins. As the first of its kind, channelrhodopsin is a light-gated cation channel and paved the way for the new and vibrant field of optogenetics, where nerve cells are activated by light. Still, the molecular mechanism of channelrhodopsin is not understood. Here, we applied time-resolved FT-IR difference spectroscopy to channelrhodopsin-1 from Chlamydomonas augustae. It is shown that the (conductive) P2(380) intermediate decays with τ ≈ 40 ms and 200 ms after pulsed excitation. The vibrational changes between the closed and the conductive states were analyzed in the X-H stretching region (X = O, S, N), comprising vibrational changes of water molecules, sulfhydryl groups of cysteine side chains and changes of the amide A of the protein backbone. The O-H stretching vibrations of "dangling" water molecules were detected in two different states of the protein using H2 (18)O exchange. Uncoupling experiments with a 1:1 mixture of H2O:D2O provided the natural uncoupled frequencies of the four O-H (and O-D) stretches of these water molecules, each with a very weakly hydrogen-bonded O-H group (3639 and 3628 cm(-1)) and with the other O-H group medium (3440 cm(-1)) to moderately strongly (3300 cm(-1)) hydrogen-bonded. Changes in amide A and thiol vibrations report on global and local changes, respectively, associated with the formation of the conductive state. Future studies will aim at assigning the respective cysteine group(s) and at localizing the "dangling" water molecules within the protein, providing a better understanding of their functional relevance in CaChR1.

Density functional tight binding: values of semi-empirical methods in an ab initio era

Semi-empirical (SE) methods are derived from Hartree-Fock (HF) or Density Functional Theory (DFT) by neglect and approximation of electronic integrals. Thereby, parameters are introduced which have to be determined from reference calculations and/or by fitting to available experimental data. This leads to computational methods that are about 2-3 orders of magnitude faster than the standard HF/DFT methods using medium sized basis sets while being about 3 orders of magnitude slower than empirical force field methods (Molecular Mechanics: MM). Therefore, SE methods are most appropriate for a specific range of applications. These include the study of systems that contain a large number of atoms and therefore being too large for ab initio or DFT methods and also problems where dynamic or entropic effects are particularly important. In the latter case, the errors made by considering a very limited number of molecular structures or neglecting entropic contributions can be much larger than the accuracy lost due to the use of SE methods. Another area where SE methods are attractive concerns the analysis of systems for which reliable MM models are not readily available. Therefore, even in an era when rapid progress is being made in ab initio methods, there is considerable interest in further developing SE methods. We illustrate this point by focusing on the discussion of recent development and application of the Density Functional Tight Binding method.

The role of protein-bound water molecules in microbial rhodopsins

Protein-bound internal water molecules are essential features of the structure and function of microbial rhodopsins. Besides structural stabilization, they act as proton conductors and even proton storage sites. Currently, the most understood model system exhibiting such features is bacteriorhodopsin (bR). During the last 20 years, the importance of water molecules for proton transport has been revealed through this protein. It has been shown that water molecules are as essential as amino acids for proton transport and biological function. In this review, we present an overview of the historical development of this research on bR. We furthermore summarize the recently discovered protein-bound water features associated with proton transport. Specifically, we discuss a pentameric water/amino acid arrangement close to the protonated Schiff base as central proton-binding site, a protonated water cluster as proton storage site at the proton-release site, and a transient linear water chain at the proton uptake site. We highlight how protein conformational changes reposition or reorient internal water molecules, thereby guiding proton transport. Last, we compare the water positions in bR with those in other microbial rhodopsins to elucidate how protein-bound water molecules guide the function of microbial rhodopsins. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

Convergence of QM/MM and Cluster Models for the Spectroscopic Properties of the Oxygen-Evolving Complex in Photosystem II

The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.

Conformational changes in the archaerhodopsin-3 proton pump: detection of conserved strongly hydrogen bonded water networks

Archaerhodopsin-3 (AR3) is a light-driven proton pump from Halorubrum sodomense, but little is known about its photocycle. Recent interest has focused on AR3 because of its ability to serve both as a high-performance, genetically-targetable optical silencer of neuronal activity and as a membrane voltage sensor. We examined light-activated structural changes of the protein, retinal chromophore, and internal water molecules during the photocycle of AR3. Low-temperature and rapid-scan time-resolved FTIR-difference spectroscopy revealed that conformational changes during formation of the K, M, and N photocycle intermediates are similar, although not identical, to bacteriorhodopsin (BR). Positive/negative bands in the region above 3,600 cm( - 1), which have previously been assigned to structural changes of weakly hydrogen bonded internal water molecules, were substantially different between AR3 and BR. This included the absence of positive bands recently associated with a chain of proton transporting water molecules in the cytoplasmic channel and a weakly hydrogen bonded water (W401), which is part of a hydrogen-bonded pentagonal cluster located near the retinal Schiff base. However, many of the broad IR continuum absorption changes below 3,000 cm( - 1) assigned to networks of water molecules involved in proton transport through cytoplasmic and extracellular portions in BR were very similar in AR3. This work and subsequent studies comparing BR and AR3 structural changes will help identify conserved elements in BR-like proton pumps as well as bioengineer AR3 to optimize neural silencing and voltage sensing.

Proton storage site in bacteriorhodopsin: new insights from quantum mechanics/molecular mechanics simulations of microscopic pK(a) and infrared spectra

Identifying the group that acts as the proton storage/loading site is a challenging but important problem for understanding the mechanism of proton pumping in biomolecular proton pumps, such as bacteriorhodopsin (bR) and cytochrome c oxidase. Recent experimental studies of bR propelled the idea that the proton storage/release group (PRG) in bR is not an amino acid but a water cluster embedded in the protein. We argue that this idea is at odds with our knowledge of protein electrostatics, since invoking the water cluster as the PRG would require the protein to raise the pK(a) of a hydronium by almost 11 pK(a) units, which is difficult considering known cases of pK(a) shifts in proteins. Our recent quantum mechanics/molecular mechanics (QM/MM) simulations suggested an alternative "intermolecular proton bond" model in which the stored proton is shared between two conserved Glu residues (194 and 204). Here we show that this model leads to microscopic pK(a) values consistent with available experimental data and the functional requirement of a PRG. Extensive QM/MM simulations also show that, independent of a number of technical issues, such as the influence of QM region size, starting X-ray structure, and nuclear quantum effects, the "intermolecular proton bond" model is qualitatively consistent with available spectroscopic data. Potential of mean force calculations show explicitly that the stored proton strongly prefers the pair of Glu residues over the water cluster. The results and analyses help highlight the importance of considering protein electrostatics and provide arguments for why the "intermolecular proton bond" model is likely applicable to the PRG in biomolecular proton pumps in general.

Parameterization of Org27569: an allosteric modulator of the cannabinoid CB1 G protein-coupled receptor

The cannabinoid CB1 receptor is a class A G protein-coupled receptor (GPCR) that is the most widely expressed GPCR in the brain. Many GPCRs contain allosteric binding sites for endogenous and/or synthetic ligands, which are topographically distinct from the agonist-binding site that is known as the orthosteric site. While both endogenous and synthetic ligands that act at the CB1 orthosteric site have been known for some time, compounds that act at a CB1 allosteric site have only recently been discovered. The most studied of these is 5-chloro-3-ethyl-1H-indole-2-carboxylic acid [2-(4-piperidin-1-ylphenyl)ethyl]amide (Org27569). Because allosteric ligands are thought to act through conformational changes in the receptor that are transmitted from the allosteric to the orthosteric site, computational studies of the structural and dynamic interactions of Org27569 with the CB1 receptor are crucial to achieve a molecular level understanding of the basis of action of this important new class of compounds. To date, such computational studies have not been possible due to the lack of a complete set of molecular mechanics force field parameters for Org27569. Here, we present the development of missing CHARMM force field parameters for Org27569 using previously published methods and the validation and application of these new parameters using normal mode analysis and molecular dynamics simulations combined with experimental infrared measurements.

Proton transfer via a transient linear water-molecule chain in a membrane protein

High-resolution protein ground-state structures of proton pumps and channels have revealed internal protein-bound water molecules. Their possible active involvement in protein function has recently come into focus. An illustration of the formation of a protonated protein-bound water cluster that is actively involved in proton transfer was described for the membrane protein bacteriorhodopsin (bR) [Garczarek F, Gerwert K (2006) Nature 439:109-112]. Here we show through a combination of time-resolved FTIR spectroscopy and molecular dynamics simulations that three protein-bound water molecules are rearranged by a protein conformational change that resulted in a transient Grotthuss-type proton-transfer chain extending through a hydrophobic protein region of bR. This transient linear water chain facilitates proton transfer at an intermediate conformation only, thereby directing proton transfer within the protein. The rearrangement of protein-bound water molecules that we describe, from inactive positions in the ground state to an active chain in an intermediate state, appears to be energetically favored relative to transient incorporation of water molecules from the bulk. Our discovery provides insight into proton-transfer mechanisms through hydrophobic core regions of ubiquitous membrane spanning proteins such as G-protein coupled receptors or cytochrome C oxidases.

DFTB3: Extension of the self-consistent-charge density-functional tight-binding method (SCC-DFTB)

The self-consistent-charge density-functional tight-binding method (SCC-DFTB) is an approximate quantum chemical method derived from density functional theory (DFT) based on a second-order expansion of the DFT total energy around a reference density. In the present study we combine earlier extensions and improve them consistently with, first, an improved Coulomb interaction between atomic partial charges, and second, the complete third-order expansion of the DFT total energy. These modifications lead us to the next generation of the DFTB methodology called DFTB3, which substantially improves the description of charged systems containing elements C, H, N, O, and P, especially regarding hydrogen binding energies and proton affinities. As a result, DFTB3 is particularly applicable to biomolecular systems. Remaining challenges and possible solutions are also briefly discussed.

Real-time structural transitions are coupled to chemical steps in ATP hydrolysis by Eg5 kinesin

At the biochemical level, motor proteins are enzymatic molecules that function by converting chemical energy into mechanical motion. The key element for energy transduction and a major unresolved question common for all motor proteins is the coordination between the chemical and conformational steps in ATP hydrolysis. Here we show time-lapse monitoring of an in vitro ATP hydrolysis reaction by the motor domain of a human Kinesin-5 protein (Eg5) using difference Fourier transform infrared spectroscopy and UV photolysis of caged ATP. In this first continuous observation of a biological reaction coordinate from substrate to product, direct spectral markers for two catalytic events are measured: proton abstraction from nucleophilic water by the catalytic base and formation of the inorganic phosphate leaving group. Simultaneous examination of conformational switching in Eg5, in parallel with catalytic steps, shows structural transitions in solution consistent with published crystal structures of the prehydrolytic and ADP-bound states. In addition, we detect structural modifications in the Eg5 motor domain during bond cleavage between the beta- and gamma-phosphates of ATP. This conclusion challenges mechanochemical models for motor proteins that utilize only two stages of the catalytic cycle to drive force and motion.

Potential proton-release channels in bacteriorhodopsin

The protein bacteriorhodopsin pumps protons across a bacterial membrane; its pumping cycle is triggered by the photoisomerization of a retinal cofactor and involves multiple proton-transfer reactions between intermittent protonation sites. These transfers are either direct or mediated by hydrogen-bonded networks, which may include internal water molecules. The terminal step of the proton-transfer sequence is the proton release from a pocket near Glu194 and Glu204 to the extracellular bulk during the transition from the L to the M photointermediate states. The polar and charged side chains connecting these two regions in the crystal structures show no structural changes between the initial bR state and the L/M states, and no intermittent protonation changes have been detected so far in this region. Based on biomolecular simulations, we propose two potential proton-release channels, which connect the release pocket to the extracellular medium. In simulations of the L photointermediate we observe bulk water entering these channels and forming transient hydrogen-bonded networks, which could serve as fast deprotonation pathways from the release pocket to the bulk via a Grotthuss mechanism. For the first channel, we find that the triple Arg7, Glu9, and Tyr79 acts as a valve, thereby gating water uptake and release. The second channel has two release paths, which split at the position Asn76/Pro77 underneath the release group. Here, water molecules either exchange directly with the bulk or diffuse within the protein towards Arg 134/Lys129, where the exchange with the bulk occurs.