Microsecond exchange of internal water molecules 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.

Hydroxide Ion Mechanism for Long-Range Proton Pumping in the Third Proton Transfer of Bacteriorhodopsin

In bacteriorhodopsin, representative light-driven proton pump, five proton transfers yield vectorial active proton translocation, resulting in a proton gradient in microbes. Third proton transfer occurs from Asp96 to the Schiff base on the photocycle, which is expected to be a long-range proton transfer via the Grotthuss mechanism through internal water molecules. Here, large-scale quantum molecular dynamics simulations are performed for the third proton transfer, where all the atoms (~50000 atoms) are treated quantum-mechanically. The simulations demonstrate that two reaction paths exist along the water wire, namely, via hydronium and via hydroxide ions. The free energy analysis confirms that the path via hydroxide ions is considerably favorable and consistent with the observed lifetime of the transient water wire. Therefore, the proposed hydroxide ion mechanism, as in the first proton transfer, is responsible for the third long-range proton transfer.

Graphs of Hydrogen-Bond Networks to Dissect Protein Conformational Dynamics

Dynamic hydrogen bonds and hydrogen-bond networks are ubiquitous in proteins and protein complexes. Functional roles that have been assigned to hydrogen-bond networks include structural plasticity for protein function, allosteric conformational coupling, long-distance proton transfers, and transient storage of protons. Advances in structural biology provide invaluable insights into architectures of large proteins and protein complexes of direct interest to human physiology and disease, including G Protein Coupled Receptors (GPCRs) and the SARS-Covid-19 spike protein S, and give rise to the challenge of how to identify those interactions that are more likely to govern protein dynamics. This Perspective discusses applications of graph-based algorithms to dissect dynamical hydrogen-bond networks of protein complexes, with illustrations for GPCRs and spike protein S. H-bond graphs provide an overview of sites in GPCR structures where hydrogen-bond dynamics would be required to assemble longer-distance networks between functionally important motifs. In the case of spike protein S, graphs identify regions of the protein where hydrogen bonds rearrange during the reaction cycle and where local hydrogen-bond networks likely change in a virus variant of concern.

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

Mechanisms of long-distance allosteric couplings in proton-binding membrane transporters

Membrane transporters that use proton binding and proton transfer for function couple local protonation change with changes in protein conformation and water dynamics. Changes of protein conformation might be required to allow transient formation of hydrogen-bond networks that bridge proton donor and acceptor pairs separated by long distances. Inter-helical hydrogen-bond networks adjust rapidly to protonation change, and ensure rapid response of the protein structure and dynamics. Membrane transporters with known three-dimensional structures and proton-binding groups inform on general principles of protonation-coupled protein conformational dynamics. Inter-helical hydrogen bond motifs between proton-binding carboxylate groups and a polar sidechain are observed in unrelated membrane transporters, suggesting common principles of coupling protonation change with protein conformational dynamics.

Protein Structure Prediction in CASP13 Using AWSEM-Suite

Recently several techniques have emerged that significantly enhance the quality of predictions of protein tertiary structures. In this study, we describe the performance of AWSEM-Suite, an algorithm that incorporates template-based modeling and coevolutionary restraints with a realistic coarse-grained force field, AWSEM. With its roots in neural networks, AWSEM contains both physical and bioinformatical energies that have been optimized using energy landscape theory. AWSEM-Suite participated in CASP13 as a server predictor and generated reliable predictions for most targets. AWSEM-Suite ranked eighth in both the free-modeling category and the hard-to-model category and in one case provided the best submitted prediction. Here we critically discuss the prediction performance of AWSEM-Suite using several examples from different categories in CASP13. Structure prediction tests on these selected targets, two of them being hard-to-model targets, show that AWSEM-Suite can achieve high-resolution structure prediction after incorporating both template guidances and coevolutionary restraints even when homology is weak. For targets with reliable templates (template-easy category), introducing coevolutionary restraints sometimes damages the overall quality of the predictions. Free energy profile analyses demonstrate, however, that the incorporations of both of these evolutionarily informed terms effectively increase the funneling of the landscape toward native-like structures while still allowing sufficient flexibility to correct for discrepancies between the correct target structure and the provided guidance. In contrast to other predictors that are exclusively oriented toward structure prediction, the connection of AWSEM-Suite to a statistical mechanical basis and affiliated molecular dynamics and importance sampling simulations makes it suitable for functional explorations.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Pf1 bacteriophage hydration by magic angle spinning solid-state NMR

High resolution two- and three-dimensional heteronuclear correlation spectroscopy ((1)H-(13)C, (1)H-(15)N, and (1)H-(13)C-(13)C HETCOR) has provided a detailed characterization of the internal and external hydration water of the Pf1 virion. This long and slender virion (2000 nm × 7 nm) contains highly stretched DNA within a capsid of small protein subunits, each only 46 amino acid residues. HETCOR cross-peaks have been unambiguously assigned to 25 amino acids, including most external residues 1-21 as well as residues 39-40 and 43-46 deep inside the virion. In addition, the deoxyribose rings of the DNA near the virion axis are in contact with water. The sets of cross-peaks to the DNA and to all 25 amino acid residues were from the same hydration water (1)H resonance; some of the assigned residues do not have exchangeable side-chain protons. A mapping of the contacts onto structural models indicates the presence of water "tunnels" through a highly hydrophobic region of the capsid. The present results significantly extend and modify results from a lower resolution study, and yield a comprehensive hydration surface map of Pf1. In addition, the internal water could be distinguished from external hydration water by means of paramagnetic relaxation enhancement. The internal water population may serve as a conveniently localized magnetization reservoir for structural studies.

Before and beyond the micellization of n-alkyl glycosides. A water-1H NMR relaxation study

The interactions between the headgroups of n-alkyl glycoside (AG) and water molecules were studied by nuclear magnetic transverse relaxation times (T2) of the water protons before and beyond the micellization. Despite the low concentration of the surfactants (mM), their micellization induce strong effect on the T2 values of bulk water when the AG molecules self-aggregate into micelles. This is associated with the decreasing of the fraction of OH headgroups of AG to exchange protons with water molecules due to the OH headgroups intermolecular interactions of AG at the micelle surface. These findings support the computational results described in the literature, which indicate that the water hydrogen bonding to OH headgroups is perturbed at AG micelle surfaces.

Characteristics of halorhodopsin-bacterioruberin complex from Natronomonas pharaonis membrane in the solubilized system

Halorhodopsin is a retinal protein with a seven-transmembrane helix and acts as an inward light-driven Cl(-) pump. In this study, structural state of the solubilized halorhodopsin (NpHR) from the biomembrane of mutant strain KM-1 of Natronomonas pharaonis in nonionic detergent was investigated. A gel filtration chromatography monitored absorbances at 280 and 504 nm corresponding to the protein and a lipid soluble pigment of bacterioruberin (BR), respectively, has clearly detected an oligomer formation of the NpHRs and a complex formation between the NpHR and BR in the solubilized system. A molar ratio of NpHR:BR in the solubilized complex was close to 1:1. Further SDS-PAGE analysis of the solubilized NpHR cross-linked by 1% glutaraldehyde has revealed that the NpHR forms homotrimer in detergent system. Although this trimeric structure was stable in the presence of NaCl, it was dissociated to the monomer by the heat treatment at 45 °C in the desalted condition. The same tendency has been reported in the case of trimeric NpHR expressed heterologously on the E. coli membrane, leading to a conclusion that the change of strength of the trimeric association dependent on the ion binding is a universal feature of the NpHR. Interestingly, the trimer dissociation on the NpHR was accompanied by the complete dissociation of the BR molecule from the protein, indicated that the cavity formed by the NpHR protomers in the trimeric conformation is important for tight binding of the BR. Because the binding affinity for Cl(-) and the resistance to hydroxylamine under light illumination showed only minor differences between the NpHR in the solubilized state and that on the biomembrane, the influences of solubilization to the tertiary structure and function of the protein are thought to be minor. This NpHR-BR complex in the solubilized system has a potential to be a good model system to investigate the intermolecular interaction between the membrane protein and lipid.

Water-assisted proton transfer in ferredoxin I

The role of water molecules in assisting proton transfer (PT) is investigated for the proton-pumping protein ferredoxin I (FdI) from Azotobacter vinelandii. It was shown previously that individual water molecules can stabilize between Asp(15) and the buried [3Fe-4S](0) cluster and thus can potentially act as a proton relay in transferring H(+) from the protein to the μ(2) sulfur atom. Here, we generalize molecular mechanics with proton transfer to studying proton transfer reactions in the condensed phase. Both umbrella sampling simulations and electronic structure calculations suggest that the PT Asp(15)-COOH + H(2)O + [3Fe-4S](0) → Asp(15)-COO(-) + H(2)O + [3Fe-4S](0) H(+) is concerted, and no stable intermediate hydronium ion (H(3)O(+)) is expected. The free energy difference of 11.7 kcal/mol for the forward reaction is in good agreement with the experimental value (13.3 kcal/mol). For the reverse reaction (Asp(15)-COO(-) + H(2)O + [3Fe-4S](0)H(+) → Asp(15)-COOH + H(2)O + [3Fe-4S](0)), a larger barrier than for the forward reaction is correctly predicted, but it is quantitatively overestimated (23.1 kcal/mol from simulations versus 14.1 from experiment). Possible reasons for this discrepancy are discussed. Compared with the water-assisted process (ΔE ≈ 10 kcal/mol), water-unassisted proton transfer yields a considerably higher barrier of ΔE ≈ 35 kcal/mol.

Protein fluctuations as the possible origin of the thermal activation of rod photoreceptors in the dark

Efficient retinal photoisomerization, signal transduction, and amplification contribute to single-photon electrical responses in vertebrates visual cells. However, spontaneous discrete electrical signals arising in the dark, with identical intensity and time profiles as those generated by genuine single photons (dark events), limit the potential capability of the rod visual system to discern single photons from thermal noise. It is accepted that the light and the thermal activation of the rod photoreceptor rhodopsin (Rho) triggers the light and the dark events, respectively. However the activation barrier for the dark events (80-110 kJ/mol) appears to be only half of the barrier for light-dependent activation of Rho (> or =180 kJ/mol). On the basis of these observations, it has been postulated that both processes should follow different pathways, but the molecular mechanism for the thermal activation process still remains an open question and subject of debate. Here, performing infrared difference spectroscopy measurements, we found that the -OH group of Thr118 from bovine Rho exhibits a slow but measurable hydrogen/deuterium exchange (HDX) under native conditions. Given the location of Thr118 in the X-ray structures, isolated from the aqueous phase and in steric contact with the buried retinal chromophore, we assume that a protein structural fluctuation must drive the retinal binding pocket (RBP) transiently open. We characterized the kinetics (rate and activation enthalpy) and thermodynamics (equilibrium constant and enthalpy) of this fluctuation from the global analysis of the HDX of Thr118-OH as a function of the temperature and pH. In parallel, using HPLC chromatography, we determined the kinetics of the thermal isomerization of the protonated 11-cis retinal in solution, as a model for retinal thermal isomerization in an open RBP. Finally, we propose a quantitative two-step model in which the dark activation of Rho is triggered by thermal isomerization of the retinal in a transiently opened RBP, which accurately reproduced both the experimental activation barrier and the rate of the dark events. We conclude that the absolute sensitivity threshold of our visual system is limited by structural fluctuations of the chromophore binding pocket rather than in the chromophore itself.

From protons to OXPHOS supercomplexes and Alzheimer's disease: structure-dynamics-function relationships of energy-transducing membranes

By the elucidation of high-resolution structures the view of the bioenergetic processes has become more precise. But in the face of these fundamental advances, many problems are still unresolved. We have examined a variety of aspects of energy-transducing membranes from large protein complexes down to the level of protons and functional relevant picosecond protein dynamics. Based on the central role of the ATP synthase for supplying the biological fuel ATP, one main emphasis was put on this protein complex from both chloroplast and mitochondria. In particular the stoichiometry of protons required for the synthesis of one ATP molecule and the supramolecular organisation of ATP synthases were examined. Since formation of supercomplexes also concerns other complexes of the respiratory chain, our work was directed to unravel this kind of organisation, e.g. of the OXPHOS supercomplex I(1)III(2)IV(1), in terms of structure and function. Not only the large protein complexes or supercomplexes work as key players for biological energy conversion, but also small components as quinones which facilitate the transfer of electrons and protons. Therefore, their location in the membrane profile was determined by neutron diffraction. Physico-chemical features of the path of protons from the generators of the electrochemical gradient to the ATP synthase, as well as of their interaction with the membrane surface, could be elucidated by time-resolved absorption spectroscopy in combination with optical pH indicators. Diseases such as Alzheimer's dementia (AD) are triggered by perturbation of membranes and bioenergetics as demonstrated by our neutron scattering studies.

Investigation of water bound to photosystem I with multiquantum filtered 17 O nuclear magnetic resonance

A new analytical approach was developed to characterize the properties of water molecules bound to macromolecules in solution using 17 O nuclear magnetic resonance (NMR) relaxation. A combination of conventional (single-quantum) and triple-quantum filtered Hahn echo and inversion recovery measurements was employed. From measured relaxation rate constants, the fraction and the correlation time of bound H2 17 O molecules and the relaxation rate constant of bulk water in solution were calculated. This was done by solving analytically a set of nonlinear equations describing the overall relaxation rate constants in the presence of chemical exchange between bulk and bound water. The analytical approach shows the uniqueness of the solution for a given set of three relaxation rate constants. This important result sheds light on the data reduction problem from 17 O NMR experiments on biological systems. Water bound in photosystem I isolated from the wild type and rubA variant of the cyanobacterium Synechocystis species PCC 7002 was investigated for the first time. The analysis revealed that photosystem I isolated from the wild type binds 1720+/-110 water molecules, whereas photosystem I isolated from the rubA variant binds only 1310+/-170. The accuracy of the method proposed can be increased by further 17 O enrichment. The methodology, established for the first time in this work, allows the study of a diverse range of biological samples regardless of their size and molecular weight. Applied initially to photosystem I, this novel method has important consequences for the future investigation of the assembly of biological molecules.

Metastable water clusters in the nonpolar cavities of the thermostable protein tetrabrachion

Water expulsion from the protein core is a key step in protein folding. Nevertheless, unusually large water clusters confined into the nonpolar cavities have been observed in the X-ray crystal structures of tetrabrachion, a bacterial protein that is thermostable up to at least 403 K (130 degrees C). Here, we use molecular dynamics (MD) simulations to investigate the structure and thermodynamics of water filling the largest cavity of the right-handed coiled-coil stalk of tetrabrachion at 365 K (92 degrees C), the temperature of optimal bacterial growth, and at room temperature (298 K). Hydrogen-bonded water clusters of seven to nine water molecules are found to be thermodynamically stable in this cavity at both temperatures, confirming the X-ray studies. Stability, as measured by the transfer free energy of the optimal size cluster, decreases with increasing temperature. Water filling is thus driven by the energy of transfer and opposed by the transfer entropy, both depending only weakly on temperature. Our calculations suggest that cluster formation becomes unfavorable at approximately 384 K (110 degrees C), signaling the onset of drying just slightly above the temperature of optimal growth. "Drying" thus precedes protein denaturation. At room temperature, the second largest cavity in tetrabrachion accommodates a five water molecule cluster, as reported in the X-ray studies. However, the simulations show that at 365 K the cluster is unstable and breaks up. We suggest that the large hydrophobic cavities may act as binding sites for two proteases, possibly explaining the unusual thermostability of the resulting protease-stalk complexes (up to approximately 393 K, 120 degrees C).

Water as a cofactor in the unidirectional light-driven proton transfer steps in bacteriorhodopsin

Recent evidence for involvement of internal water molecules in the mechanism of bacteriorhodopsin is reviewed. Water O-H stretching vibration bands in the Fourier transform IR difference spectra of the L, M and N intermediates of bacteriorhodopsin were analyzed by photoreactions at cryogenic temperatures. A broad vibrational band in L was shown to be due to formation of a structure of water molecules connecting the Schiff base to the Thr46-Asp96 region. This structure disappears in the M intermediate, suggesting that it is involved in transient stabilization of the L intermediate prior to proton transfer from the Schiff base to Asp85. The interaction of the Schiff base with a water molecule is restored in the N intermediate. We propose that water is a critical mobile component of bacteriorhodopsin, forming organized structures in the transient intermediates during the photocycle and, to a large extent, determining the chemical behavior of these transient states.

Protein Dynamics Tightly Connected to the Dynamics of Surrounding and Internal Water Molecules

Proteins are key components of biological cells. For example, enzymes catalyze biochemical reactions, membrane transporters are responsible for uptake and release of critical and superfluous components from the cell environment, and structural proteins are responsible for the stability of the cell wall and cytoskeleton. Many of the diverse protein functions involve dynamic transitions ranging from small local atomic displacements up to large allosteric conformational changes. In any conformation, proteins are in contact with the universal solvent medium of cells, water. Water not only surrounds proteins but is often an integral part of proteins and also is involved in key mechanistic steps. This Minireview discusses recent experimental and theoretical results on the role of water for protein dynamics and function.

Effects of additives on surfactant phase behavior relevant to bacteriorhodopsin crystallization

The interactions leading to crystallization of the integral membrane protein bacteriorhodopsin solubilized in n-octyl-beta-D-glucoside were investigated. Osmotic second virial coefficients (B(22)) were measured by self-interaction chromatography using a wide range of additives and precipitants, including polyethylene glycol (PEG) and heptane-1,2,3-triol (HT). In all cases, attractive protein-detergent complex (PDC) interactions were observed near the surfactant cloud point temperature, and there is a correlation between the surfactant cloud point temperatures and PDC B(22) values. Light scattering, isothermal titration calorimetry, and tensiometry reveal that although the underlying reasons for the patterns of interaction may be different for various combinations of precipitants and additives, surfactant phase behavior plays an important role in promoting crystallization. In most cases, solution conditions that led to crystallization fell within a similar range of slightly negative B(22) values, suggesting that weakly attractive interactions are important as they are for soluble proteins. However, the sensitivity of the cloud point temperatures and resultant coexistence curves varied significantly as a function of precipitant type, which suggests that different types of forces are involved in driving phase separation depending on the precipitant used.

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.

Intermolecular electrostatic interactions and Brownian tumbling in protein solutions

It is often implicitly assumed that the long-range intermolecular electrostatic interactions in homogeneous protein solutions either are negligible for affecting protein Brownian tumbling or cause its deceleration without changing the shape of rotational auto-correlation function. This review presents a wide set of experimental data (NMR relaxation, dielectric spectroscopy and Brownian dynamics simulations) demonstrating that the interprotein electrostatic steering leads to a complication of the rotational correlation function. The key point of this effect is the rotational anisotropy caused by the interaction of the electric dipole moment of a protein with the external electric field produced by charges of neighboring proteins. Taking this effect into account in some cases might be of critical importance for the correct interpretation of various experimental data.

Detection of dynamic water molecules in a microcrystalline sample of the SH3 domain of alpha-spectrin by MAS solid-state NMR

Water molecules are a major determinant of protein stability and are important for understanding protein-protein interactions. We present two experiments which allow to measure first the effective T(2) decay rate of individual amide proton, and second the magnetization build-up rates for a selective transfer from H(2)O to H(N) using spin diffusion as a mixing element. The experiments are demonstrated for a uniformly (2)H, (15)N labeled sample of a microcrystalline SH3 domain in which exchangeable deuterons were back-substituted with protons. In order to evaluate the NMR experimental data, as X-ray structure of the protein was determined using the same crystallization protocol as for the solid-state NMR sample. The NMR experimental data are correlated with the dipolar couplings calculated from H(2)O-H(N) distances which were extracted from the X-ray structure of the protein. We find that the H(N) T(2) decay rates and H(2)O-H(N) build-up rates are sensitive to distance and dynamics of the detected water molecules with respect to the protein. We show that qualitative information about localization and dynamics of internal water molecules can be obtained in the solid-state by interpretation of the spin dynamics of a reporter amide proton.

Role of Putative Anion-Binding Sites in Cytoplasmic and Extracellular Channels ofNatronomonas pharaonisHalorhodopsin†

Natronomonas (Natronobacterium) pharaonis halorhodopsin (NpHR) is an inward light-driven Cl(-) ion pump. For efficient Cl(-) transport, the existence of Cl(-)-binding or -interacting sites in both extracellular (EC) and cytoplasmic (CP) channels is postulated. Candidates include Arg123 and Thr126 in EC channels and Lys215 and Thr218 in CP channels. The roles played by these amino acid residues in anion binding and in the photocycle have been investigated by mutation of the amino acid residues at these positions. Anion binding was assayed by changes in circular dichroism and the shift in the absorption maximum upon addition of Cl(-) to anion-free NpHR. The binding affinity was affected in mutants in which certain EC residues had been replaced; this finding revealed the importance of Arg123. On the other hand, mutants in which certain residues in the CP channel were replaced (CP mutants) did not show changes in their dissociation constants. The photocycles of these mutants were also examined, and in the case of the EC mutants, the transition to the last step was greatly delayed; on the other hand, in the CP mutants, L2-photointermediate decay was significantly prolonged, except in the case of K215Q, which lacked the O-photointermediate. The importance of Thr218 for binding of Cl(-) to the CP channel was indicated by these results. On the basis of these observations, the possible anion transport mechanism of NpHR was discussed.

Water molecules and hydrogen-bonded networks in bacteriorhodopsin--molecular dynamics simulations of the ground state and the M-intermediate

Protein crystallography provides the structure of a protein, averaged over all elementary cells during data collection time. Thus, it has only a limited access to diffusive processes. This article demonstrates how molecular dynamics simulations can elucidate structure-function relationships in bacteriorhodopsin (bR) involving water molecules. The spatial distribution of water molecules and their corresponding hydrogen-bonded networks inside bR in its ground state (G) and late M intermediate conformations were investigated by molecular dynamics simulations. The simulations reveal a much higher average number of internal water molecules per monomer (28 in the G and 36 in the M) than observed in crystal structures (18 and 22, respectively). We found nine water molecules trapped and 19 diffusive inside the G-monomer, and 13 trapped and 23 diffusive inside the M-monomer. The exchange of a set of diffusive internal water molecules follows an exponential decay with a 1/e time in the order of 340 ps for the G state and 460 ps for the M state. The average residence time of a diffusive water molecule inside the protein is approximately 95 ps for the G state and 110 ps for the M state. We have used the Grotthuss model to describe the possible proton transport through the hydrogen-bonded networks inside the protein, which is built up in the picosecond-to-nanosecond time domains. Comparing the water distribution and hydrogen-bonded networks of the two different states, we suggest possible pathways for proton hopping and water movement inside bR.

Photobleaching of Bacteriorhodopsin Solubilized with Triton X-100

In the current studies, we examined the effects of hexagonal lattice formation with lipid membranes on the structural stability of native bacteriorhodopsin (bR). Denaturation kinetic measurements for bR solubilized with the mild nonionic detergent Triton X-100 (TX100) were performed in the dark and under illumination by visible light. The solubilized bR was stable in the dark over a wide concentration range of TX100 (1 to 200 mM). In purple membranes, a bilobed band was observed in visible circular dichroism spectra due to interactions between neighboring chromophores. At all concentrations of TX100, this was replaced by a single positive band. Upon illumination with visible light, TX100-solubilized bR clearly showed photobleaching to bacterioopsin. These experimental results suggest that photobleaching is due to a lack of intermolecular interactions inside the purple membrane lattice. Extensive kinetic measurements further revealed that the rate constant of photobleaching is strongly dependent on the detergent concentration, although the activation energy for photobleaching does not significantly change with the TX100 concentration. The mechanism of photobleaching for the solubilized bR is discussed with respect to detergent micelle properties.

Role of water mediated interactions in protein-protein recognition landscapes

The energy landscape picture of protein folding and binding is employed to optimize a number of pair potentials for direct and water-mediated interactions in protein complex interfaces. We find that water-mediated interactions greatly complement direct interactions in discriminating against various types of trap interactions that model those present in the cell. We highlight the context dependent nature of knowledge-based binding potentials, as contrasted with the situation for autonomous folding. By performing a Principal Component Analysis (PCA) of the corresponding interaction matrixes, we rationalize the strength of the recognition signal for each combination of the contact type and reference trap states using the differential in the idealized "canonical" amino acid compositions of native and trap layers. The comparison of direct and water-mediated contact potential matrixes emphasizes the importance of partial solvation in stabilizing charged groups in the protein interfaces. Specific water-mediated interresidue interactions are expected to influence significantly the kinetics as well as thermodynamics of protein association.

Zinc-porphyrin Solvation in Folded and Unfolded States of Zn-cytochrome c

After a brief review of the use of photochemical triggers and heme metal substitution to probe the folding dynamics of cytochrome c, we present new results on the photophysics and photochemistry of folded and unfolded states of the zinc-substituted protein (Zn-cyt c). Our measurements of Zn-cyt c triplet state decay kinetics reveal a systematic isotope effect on lifetimes: the decay in the folded protein (tau(H)2(O) approximately 10 ms) is only modestly affected by isotopically substituted buffers (k(H)2(O)/k(D)2(O) = 1.2), whereas a reduced triplet lifetime (approximately 1.3 ms) and greater isotope effect (1.4) were found for the chemically denatured, fully unfolded protein. The shortest lifetime (0.1-0.4 ms) and greatest isotope effect (1.5) were found for a fully exposed model compound, zinc-substituted N-acetyl-microperoxidase-8 (ZnAcMP8), implying that the unfolded protein provides some protection to the Zn-porphyrin group even under fully denaturing conditions. Further evidence for partial structure in unfolded Zn-cyt c comes from bimolecular quenching experiments using Ru(NH(3))(6)(3+) as an external Zn-porphyrin triplet state quencher. In the presence of quencher, partially unfolded protein at midpoint guanidinium chloride (GdmCl) and urea concentrations exhibits biphasic triplet decay kinetics, a fast component corresponding to an extended, solvent-exposed state (6.6 x 10(8) M(-1) s(-1) in GdmCl, 6.3 x 10(8) M(-1) s(-1) in urea) and a slow component attributable to a compact, relatively solvent-inaccessible, state (5.9 x 10(7) M(-1) s(-1) in GdmCl, 8.6 x 10(6) M(-1) s(-1) in urea). The variation in Zn-porphyrin solvation for the compact states in the two denaturants reveals that the cofactor in the partially unfolded protein is better protected in urea solutions.

Competitive Na(+) and Rb(+) binding in the minor groove of DNA

Sequence-dependent coordination of alkali ions to the nucleotide bases in the minor groove of AT-tract B-DNA has recently been inferred from X-ray crystallography, solution NMR and computer simulations. Here, we present new (23)Na and (87)Rb magnetic relaxation dispersion (MRD) data that demonstrate competitive and long-lived binding of Na(+) and Rb(+) ions in the minor groove of the B-DNA duplex [d(CGCGAATTCGCG)](2). The Na(+)/Rb(+) selectivity of the minor groove is found to be weak, consistent with local structural flexibility. The ion occupancies derived from the MRD data are substantially higher than previously reported, suggesting that groove-bound ions significantly influence the energetics and structural polymorphism of DNA in vivo. For example, in the presence of 0.20 M Na(+) and 0.56 M Rb(+) at 4 degrees C, the ApT site in the minor groove is occupied by a Rb(+) ion, a Na(+) ion, or a water molecule 40, 10, and 50% of the time, respectively. In the absence of Rb(+), the Na(+) occupancy increases to 50%. At 4 degrees C, the mean residence time of groove-bound ions is 0.2 +/- 0.1 micros for Rb(+) and 10 ns to 100 micros for Na(+). A shorter correlation time of 2 ns is attributed to counterions bridging cross-strand phosphate groups.

Protein hydration dynamics in solution: a critical survey

The properties of water in biological systems have been studied for well over a century by a wide range of physical techniques, but progress has been slow and erratic. Protein hydration--the perturbation of water structure and dynamics by the protein surface--has been a particularly rich source of controversy and confusion. Our aim here is to critically examine central concepts in the description of protein hydration, and to assess the experimental basis for the current view of protein hydration, with the focus on dynamic aspects. Recent oxygen-17 magnetic relaxation dispersion (MRD) experiments have shown that the vast majority of water molecules in the protein hydration layer suffer a mere twofold dynamic retardation compared with bulk water. The high mobility of hydration water ensures that all thermally activated processes at the protein-water interface, such as binding, recognition and catalysis, can proceed at high rates. The MRD-derived picture of a highly mobile hydration layer is consistent with recent molecular dynamics simulations, but is incompatible with results deduced from intermolecular nuclear Overhauser effect spectroscopy, dielectric relaxation and fluorescence spectroscopy. It is also inconsistent with the common view of hydration effects on protein hydrodynamics. Here, we show how these discrepancies can be resolved.

Biomolecular cryocrystallography: structural changes during flash-cooling

To minimize radiation damage, crystal structures of biological macromolecules are usually determined after rapid cooling to cryogenic temperatures, some 150-200 K below the normal physiological range. The biological relevance of such structures relies on the assumption that flash-cooling is sufficiently fast to kinetically trap the macromolecule and associated solvent in a room-temperature equilibrium state. To test this assumption, we use a two-state model to calculate the structural changes expected during rapid cooling of a typical protein crystal. The analysis indicates that many degrees of freedom in a flash-cooled protein crystal are quenched at temperatures near 200 K, where local conformational and association equilibria may be strongly shifted toward low-enthalpy states. Such cryoartifacts should be most important for strongly solvent-coupled processes, such as hydration of nonpolar cavities and surface regions, conformational switching of solvent-exposed side chains, and weak ligand binding. The dynamic quenching that emerges from the model considered here can also rationalize the glass transition associated with the atomic fluctuations in the protein.

The interaction of water molecules with purple membrane suspension using2H double-quantum filter,1H and2H diffusion nuclear magnetic resonance

Bacteriorhodopsin is a membrane protein of the purple membrane (PM) of Halobacterium salinarum, which is isolated as sheets of highly organized two-dimensional hexagonal microcrystals and for which water molecules play a crucial role that affects its function as a proton pump. In this paper we used single- and double-quantum (2)H NMR as well as (1)H and (2)H diffusion NMR to characterize the interaction of water molecules with the PM in D(2)O suspensions. We found that, under the influence of a strong magnetic field on a concentrated PM sample (0.61 mM), the PM sheets affect the entire water population and a residual quadrupolar splitting (upsilon(q) approximately 5.5 Hz, 298 K, at 11.7 T) is observed for the D(2)O molecules. We found that the residual quadrupolar coupling, the creation time in which a maximal DQF signal was obtained (tau(max)), and the relative intensity of the (2)H DQF spectrum of the water molecules in the PM samples (referred to herein as NMR order parameters) are very sensitive to temperature, dilution, and chemical modifications of the PM. In concentrated PM samples in D(2)O, these NMR parameters seem to reflect the relative organization of the PM. Interestingly, we have observed that some of these parameters are sensitive to the efficiency of the trimer packing, as concluded from the apo-membrane behavior. The data for dionized blue membrane, partially delipidated sample, and detergent-treated PM show that these D(2)O NMR order parameters, which are magnetic field dependent, are sensitive to the structural integrity of the PM. In addition, we revealed that heating the PM sample inside or outside the NMR magnet has, after cooling, a different effect on the NMR characteristics of the water molecules in the concentrated PM suspensions. The difference in the D(2)O NMR order parameters for the PM samples, which were heated and cooled in the presence and in the absence of a strong magnetic field, corroborates the conclusions that the above D(2)O order parameters are indirect reflections of both microscopic and macroscopic order of the PM samples. In addition, (1)H NMR diffusion measurements showed that at least three distinct water populations could be identified, based on their diffusion coefficients. These water populations seem to correlate with different water populations previously reported for the PM system.

A discrete water exit pathway in the membrane protein cytochrome c oxidase

By using the non-redox-active Mg2+/Mn2+ site of cytochrome c oxidase as a probe, water access from the outside of the enzyme and water escape from the buried active site were studied. Water movement was time-resolved by monitoring the magnetic interaction of the oxygen isotope 17O with the Mn2+ by using a rapid freeze-quench-electron spin echo envelope modulation technique. Rapid (msec) access of water from the bulk phase to the Mn2+ was demonstrated by mixing cytochrome c oxidase with H217O. To determine whether a channel involving the Mn2+ was used for water exit from the active site, samples incubated in 17O2 were allowed to turn over approximately five times before freezing. The 17O, now in the form of H217O, was detected at the Mn2+. The significant broadening of the Mn2+ signal after the limited number of turnovers strongly suggests that the water exits the protein by means of one discrete pathway, not by random diffusion.

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.

Water molecule rearrangements around Leu93 and Trp182 in the formation of the L intermediate in bacteriorhodopsin's photocycle

After the chromophore's isomerization in the initial photochemical event in bacteriorhodopsin, the primary photoproduct K makes a thermal transition to the L intermediate, which prepares the pigment for Schiff base deprotonation in the following step (L --> M). Substantial changes in the hydrogen bonding of internal water molecules take place upon L formation. Some of these mobile waters are probably involved in changing the pK of the Schiff base and perhaps that of the proton acceptor Asp85 to allow proton movement [Maeda, A. (2001) Biochemistry (Moscow) 66, 1555-1569]. Here we show that mutations of Leu93 and Trp182, residues close to the 13-methyl group of the chromophore, allow the formation of L at much lower temperatures than in the wild type (80 K instead of 140 K). Moreover, an intense band due to weakly bound water that is peculiar for L was already present in the initial (unphotolyzed) state of each mutant at 2632 cm(-1) (in D2O) but not in the wild type. This unique, intense water band is shifted compared to the L band at 2589 cm(-1) but coincides with the band seen in L', the all-trans photoproduct of wild-type L formed at 80 K. We propose that the L93M and W182F mutations induce changes in the hydrogen bonding of one or more water molecules in the unphotolyzed states of these pigments, which are similar to those H-bonding changes that take place upon formation of L in the wild type, and thus facilitate the formation of L even at 80 K. We infer that L formation involves perturbation of a site which includes retinal, Trp182, and Leu93, and this structure is temporarily stabilized by rearranged hydrogen bonds with water molecules.

Structural proton diffusion along lipid bilayers

For H(+) transport between protein pumps, lateral diffusion along membrane surfaces represents the most efficient pathway. Along lipid bilayers, we measured a diffusion coefficient of 5.8 x 10(-5) cm(2) s(-1). It is too large to be accounted for by vehicle diffusion, considering proton transport by acid carriers. Such a speed of migration is accomplished only by the Grotthuss mechanism involving the chemical exchange of hydrogen nuclei between hydrogen-bonded water molecules on the membrane surface, and the subsequent reorganization of the hydrogen-bonded network. Reconstitution of H(+)-binding sites on the membrane surface decreased the velocity of H(+) diffusion. In the absence of immobile buffers, structural (Grotthuss) diffusion occurred over a distance of 100 micro m as shown by microelectrode aided measurements of the spatial proton distribution in the immediate membrane vicinity and spatially resolved fluorescence measurements of interfacial pH. The efficiency of the anomalously fast lateral diffusion decreased gradually with an increase in mobile buffer concentration suggesting that structural diffusion is physiologically important for distances of approximately 10 nm.

Ligand binding promotes the entropy-driven oligomerization of integrin alpha IIb beta 3

Integrin alpha(IIb)beta(3) clusters on the platelet surface after binding adhesive proteins in a process that regulates signal transduction. However, the intermolecular forces driving integrin self-association are poorly understood. This work provides new insights into integrin clustering mechanisms by demonstrating how temperature and ligand binding interact to affect the oligomeric state of alpha(IIb)beta(3). The ligand-free receptor, solubilized in thermostable octyl glucoside micelles, exhibited a cooperative transition at approximately 43 degrees C, monitored by changes in intrinsic fluorescence and circular dichroism. Both signals changed in a direction opposite to that for global unfolding, and both were diminished upon binding the fibrinogen gamma-chain ligand-mimetic peptide cHArGD. Free and bound receptors also exhibited differential sensitivity to temperature-enhanced oligomerization, as measured by dynamic light scattering, sedimentation velocity, and sedimentation equilibrium. Van't Hoff analyses of dimerization constants for alpha(IIb)beta(3) complexed with cHArGD, cRGD, or eptifibatide yielded large, favorable entropy changes partly offset by unfavorable enthalpy changes. Transmission electron microscopy showed that ligand binding and 37 degrees C incubation enhanced assembly of integrin dimers and larger oligomers linked by tail-to-tail contacts. Interpretation of these images was aided by threading models for alpha(IIb)beta(3) protomers and dimers based on the ectodomain structure of alpha(v)beta(3). We propose that entropy-favorable nonpolar interactions drive ligand-induced integrin clustering and outside-in signaling.

Pressure dependence of the photocycle kinetics of bacteriorhodopsin

The pressure dependence of the photocycle kinetics of bacteriorhodopsin from Halobacterium salinarium was investigated at pressures up to 4 kbar at 25 degrees C and 40 degrees C. The kinetics can be adequately modeled by nine apparent rate constants, which are assigned to irreversible transitions of a single relaxation chain of nine kinetically distinguishable states P(1) to P(9). All states except P(1) and P(9) consist of two or more spectral components. The kinetic states P(2) to P(6) comprise only the two fast equilibrating spectral states L and M. From the pressure dependence, the volume differences DeltaV(o)(LM) between these two spectral states could be determined that range from DeltaV(o)(LM) = -11.4 +/- 0.7 ml/mol (P(2)) to DeltaV(o)(LM) = 14.6 +/- 2.8 mL/mol (P(6)). A model is developed that explains the dependence of DeltaV(o)(LM) on the kinetic state by the electrostriction effect of charges, which are formed and neutralized during the L/M transition.

Structural evidence that propofol stabilizes different GABA(A) receptor states at potentiating and activating concentrations

The GABA(A) receptor is a target of many general anesthetics, such as propofol. General anesthetic binding sites are distinct from the GABA binding sites. At low concentrations, the anesthetics potentiate the currents induced by submaximal GABA concentrations. At higher concentrations the anesthetics directly activate GABA(A) receptors. In contrast, benzodiazepines, such as diazepam, only potentiate currents induced by submaximal GABA concentrations. Channel kinetic studies suggest that these drugs stabilize different receptor states. We previously showed that the accessibility of the anionic sulfhydryl reagent p-chloromercuribenzenesulfonate (pCMBS(-)) applied extracellularly to cysteines substituted for residues in the GABA(A) alpha1 subunit M3 membrane-spanning segment was state-dependent. The subset of pCMBS(-)-accessible, M3 segment cysteine mutants acts as a reporter for receptor conformation. Here we show that pCMBS(-), applied in the presence of a potentiating concentration of propofol, reacts with a subset of alpha1 subunit, M3 segment, cysteine-substitution mutants (Y294C, V297C, I302C, F304C). In the presence of a directly activating concentration of propofol pCMBS(-) reacts with a different subset of the M3 cysteine-substitution mutants (Y294C, S299C, I302C, E303C, A305C). These subsets are distinct from the subsets of M3 cysteine-substitution mutants that are reactive with pCMBS(-) in the absence and presence of GABA and in the presence of diazepam. We hypothesize that distinct subsets of reactive residues represent distinct conformations or ensembles of conformations of the receptor. These results provide structural evidence for at least five distinct receptor states, three nonconducting states, resting, diazepam-bound and potentiating propofol-bound, and two conducting-desensitized states, the activating propofol-bound and GABA-bound states.

Protein-water interactions in a dynamic world

Protein-water interactions are key to biological function. They have an underlying dynamic component that pervades the functional roles associated both with particular systems and with the properties of proteins in general. This article focuses on the specific ways in which the dynamics of water are important to protein structure, motion and adaptability to changes in the protein environment.