Effects of mutations of Lys41 and Asp102 of bacteriorhodopsin

Bacteriorhodopsin (BR) is a retinal protein that functions as a light-driven proton pump. In this study, six novel mutants including K41E and D102K, were obtained to verify or rule out the possibility that residues Lys41 and Asp102 are determinants of the time order of proton release and uptake, because we found that the order was reversed in another retinal protein archaerhodopsin 4 (AR4), which had different 41th and 102th residues. Our results rule out that possibility and confirm that the pK(a) of the proton release complex (PRC) determines the time order. Nevertheless, mutations, especially D102K, were found to affect the kinetics of proton uptake substantially and the pK(a) of Asp96. Compared to the wild-type BR (BR-WT), the decay of the M intermediate and proton uptake in the photocycle was slowed about 3-fold in D102K. Hence those residues might be involved in proton uptake and delivery to the internal proton donor.     

Conformational sampling and polarization of Asp26 in pKa calculations of thioredoxin

Thioredoxin is a protein that has been used as model system by various computational methods to predict the pK_a of aspartate residue Asp26 which is 3.5 units higher than a solvent exposed one (eg, Asp20). Here, we use extensive atomistic molecular dynamics simulations of two different protonation states of Asp26 in combination with conformational analysis based on RMSD clustering and principle component analysis to identify representative conformations of the protein in solution. For each conformation, the Gibbs free energy of proton transfer between Asp26 and Asp20, which is fully solvated in a loop region of the protein, is calculated with the Amber99sb force field in alchemical transformations. The varying polarization of the two residues in different molecular environments and protonation states is described by Hirshfeld-I (HI) atomic charges obtained from the averaged polarized electron density. Our results show that the Gibbs free energy of proton transfer is dependent on the protein conformation, the proper sampling of the neighboring Lys57 residue orientations and on water molecules entering the hydrophobic cavity upon deprotonating Asp26. The inclusion of the polarization of both aspartate residues in the free energy cycle by HI atomic charges corrects the results from the non-polarizable force field and reproduces the experimental ΔpK_a value of Asp26.     

Factors That Differentiate the H-bond Strengths of Water Near the Schiff Bases in Bacteriorhodopsin and Anabaena Sensory Rhodopsin

Bacteriorhodopsin (BR) functions as a light-driven proton pump, whereas Anabaena sensory rhodopsin (ASR) is believed to function as a photosensor despite the high similarity in their protein sequences. In Fourier transform infrared (FTIR) spectroscopic studies, the lowest O-D stretch for D₂O was observed at ∼2200 cm⁻¹ in BR but was significantly higher in ASR (>2500 cm⁻¹), which was previously attributed to a water molecule near the Schiff base (W402) that is H-bonded to Asp-85 in BR and Asp-75 in ASR. We investigated the factors that differentiate the lowest O-D stretches of W402 in BR and ASR. Quantum mechanical/molecular mechanical calculations reproduced the H-bond geometries of the crystal structures, and the calculated O-D stretching frequencies were corroborated by the FTIR band assignments. The potential energy profiles indicate that the smaller O-D stretching frequency in BR originates from the significantly higher pK_a(Asp-85) in BR relative to the pK_a(Asp-75) in ASR, which were calculated to be 1.5 and −5.1, respectively. The difference is mostly due to the influences of Ala-53, Arg-82, Glu-194–Glu-204, and Asp-212 on pK_a(Asp-85) in BR and the corresponding residues Ser-47, Arg-72, Ser-188-Asp-198, and Pro-206 on pK_a(Asp-75) in ASR. Because these residues participate in proton transfer pathways in BR but not in ASR, the presence of a strongly H-bonded water molecule near the Schiff base ultimately results from the proton-pumping activity in BR.     

Bacteriorhodopsin: a high-resolution structural view of vectorial proton transport

Recent 3-D structures of several intermediates in the photocycle of bacteriorhodopsin (bR) provide a detailed structural picture of this molecular proton pump in action. In this review, we describe the sequence of conformational changes of bR following the photoisomerization of its all-trans retinal chromophore, which is covalently bound via a protonated Schiff base to Lys216 in helix G, to a 13-cis configuration. The initial changes are localized near the protein's active site and a key water molecule is disordered. This water molecule serves as a keystone for the ground state of bR since, within the framework of the complex counter ion, it is important both for stabilizing the structure of the extracellular half of the protein, and for maintaining the high pK_a of the Schiff base (the primary proton donor) and the low pK_a of Asp85 (the primary proton acceptor). Subsequent structural rearrangements propagate out from the active site towards the extracellular half of the protein, with a local flex of helix C exaggerating an early movement of Asp85 towards the Schiff base, thereby facilitating proton transfer between these two groups. Other coupled rearrangements indicate the mechanism of proton release to the extracellular medium. On the cytoplasmic half of the protein, a local unwinding of helix G near the backbone of Lys216 provides sites for water molecules to order and define a pathway for the reprotonation of the Schiff base from Asp96 later in the photocycle. A steric clash of the photoisomerized retinal with Trp182 in helix F drives an outward tilt of the cytoplasmic half of this helix, opening the proton transport channel and enabling a proton to be taken up from the cytoplasm. Although bR is the first integral membrane protein to have its catalytic mechanism structurally characterized in detail, several key results were anticipated in advance of the structural model and the general framework for vectorial proton transport has, by and large, been preserved.     

Interaction between lysine 102 and aspartate 338 in the insect amino acid cotransporter KAAT1

KAAT1 is a lepidopteran neutral amino acid transporter belonging to the NSS super family (SLC6), which has an unusual cation selectivity, being activated by K⁺ and Li⁺ in addition to Na⁺. We have previously demonstrated that Asp338 is essential for KAAT1 activation by K⁺ and for the coupling of amino acid and driver ion fluxes. By comparing sequences of NSS family members, site-directed mutagenesis, and expression in Xenopus laevis oocytes, we identified Lys102 as a residue likely to interact with Asp338. Compared with wild type, the single mutants K102V and D338E each showed altered leucine uptake and transport-associated currents in the presence of both Na⁺ and K⁺. However, in K102V/D338E double mutant, the K102V mutation reversed both the inhibition of Na⁺-dependent transport and the block in K⁺-dependent transport that characterize the D338E mutant. K⁺-dependent leucine currents were not observed in any mutants with D338E. In the presence of the oxidant Cu(II) (1,10-phenanthroline)₃, we observed specific and reversible inhibition of K102C/D338C mutant, but not of the corresponding single cysteine mutants, suggesting that these residues are sufficiently close to form a disulfide bond. Thus both structural and functional evidence suggests that these two residues interact. Similar results have been obtained mutating the bacterial transporter homolog TnaT. Asp338 corresponds to Asn286, a residue located in the Na⁺ binding site in the recently solved crystal structure of the NSS transporter LeuT_Aa (41). Our results suggest that Lys102, interacting with Asp338, could contribute to the spatial organization of KAAT1 cation binding site and permeation pathway. residue interaction; oxidants; tertiary structure     

Molecular Dynamics Simulations Elucidate the Mechanism of Proton Transport in the Glutamate Transporter EAAT3

The uptake of glutamate in nerve synapses is carried out by the excitatory amino acid transporters (EAATs), involving the cotransport of a proton and three Na⁺ ions and the countertransport of a K⁺ ion. In this study, we use an EAAT3 homology model to calculate the pK_a of several titratable residues around the glutamate binding site to locate the proton carrier site involved in the translocation of the substrate. After identifying E374 as the main candidate for carrying the proton, we calculate the protonation state of this residue in different conformations of EAAT3 and with different ligands bound. We find that E374 is protonated in the fully bound state, but removing the Na2 ion and the substrate reduces the pK_a of this residue and favors the release of the proton to solution. Removing the remaining Na⁺ ions again favors the protonation of E374 in both the outward- and inward-facing states, hence the proton is not released in the empty transporter. By calculating the pK_a of E374 with a K⁺ ion bound in three possible sites, we show that binding of the K⁺ ion is necessary for the release of the proton in the inward-facing state. This suggests a mechanism in which a K⁺ ion replaces one of the ligands bound to the transporter, which may explain the faster transport rates of the EAATs compared to its archaeal homologs.     

Poisson-Boltzmann model studies of molecular electrostatic properties of the cAMP-dependent protein kinase

Protonation equilibria of residues important in the catalytic mechanism of a protein kinase were analyzed on the basis of the Poisson-Boltzmann electrostatic model along with a cluster-based treatment of the multiple titration state problem. Calculations were based upon crystallographic structures of the mammalian cAMP-dependent protein kinase, one representing the so called closed form of the enzyme and the other representing an open conformation. It was predicted that at pH 7 the preferred form of the phosphate group at the catalytically essential threonine 197 (P-Thr197) in the closed form is dianionic, whereas in the open form a monoanionic ionization state is preferred. This dianionic state of P-Thr197, in the closed form, is stabilized by interactions with ionizable residues His87, Arg165, and Lys189. Our calculations predict that the hydroxyl of the Ser residue in the peptide substrate is very difficult to ionize, both in the closed and open structures of the complex. Also, the supposed catalytic base, Asp166, does not seem to have a pK _a appropriate to remove the hydroxyl group proton of the peptide substrate. However, when Ser of the peptide substrate is forced to remain ionized, the predicted pK _a of Asp166 increases strongly, which suggests that the Asp residue is a likely candidate to attract the proton if the Ser residue becomes deprotonated, possibly during some structural change preceding formation of the transition state. Finally, in accord with suggestions made on the basis of the pH-dependence of kinase kinetics, our calculations predict that Glu230 and His87 are the residues responsible for the molecular pK _a values of 6.2 and 8.5, observed in the experiment. Received: 19 October 1998 / Revised version: 12 April 1999 / Accepted: 15 April 1999     

Breaking the Carboxyl Rule: LYSINE 96 FACILITATES REPROTONATION OF THE SCHIFF BASE IN THE PHOTOCYCLE OF A RETINAL PROTEIN FROM EXIGUOBACTERIUM SIBIRICUM*?

A lysine instead of the usual carboxyl group is in place of the internal proton donor to the retinal Schiff base in the light-driven proton pump of Exiguobacterium sibiricum (ESR). The involvement of this lysine in proton transfer is indicated by the finding that its substitution with alanine or other residues slows reprotonation of the Schiff base (decay of the M intermediate) by more than 2 orders of magnitude. In these mutants, the rate constant of the M decay linearly decreases with a decrease in proton concentration, as expected if reprotonation is limited by the uptake of a proton from the bulk. In wild type ESR, M decay is biphasic, and the rate constants are nearly pH-independent between pH 6 and 9. Proton uptake occurs after M formation but before M decay, which is especially evident in D₂O and at high pH. Proton uptake is biphasic; the amplitude of the fast phase decreases with a pK_a of 8.5 ± 0.3, which reflects the pK_a of the donor during proton uptake. Similarly, the fraction of the faster component of M decay decreases and the slower one increases, with a pK_a of 8.1 ± 0.2. The data therefore suggest that the reprotonation of the Schiff base in ESR is preceded by transient protonation of an initially unprotonated donor, which is probably the ϵ-amino group of Lys-96 or a water molecule in its vicinity, and it facilitates proton delivery from the bulk to the reaction center of the protein.     

Molecular Dynamics Simulations Elucidate the Mechanism of Proton Transport in the Glutamate Transporter EAAT3

The uptake of glutamate in nerve synapses is carried out by the excitatory amino acid transporters (EAATs), involving the cotransport of a proton and three Na⁺ ions and the countertransport of a K⁺ ion. In this study, we use an EAAT3 homology model to calculate the pK_a of several titratable residues around the glutamate binding site to locate the proton carrier site involved in the translocation of the substrate. After identifying E374 as the main candidate for carrying the proton, we calculate the protonation state of this residue in different conformations of EAAT3 and with different ligands bound. We find that E374 is protonated in the fully bound state, but removing the Na2 ion and the substrate reduces the pK_a of this residue and favors the release of the proton to solution. Removing the remaining Na⁺ ions again favors the protonation of E374 in both the outward- and inward-facing states, hence the proton is not released in the empty transporter. By calculating the pK_a of E374 with a K⁺ ion bound in three possible sites, we show that binding of the K⁺ ion is necessary for the release of the proton in the inward-facing state. This suggests a mechanism in which a K⁺ ion replaces one of the ligands bound to the transporter, which may explain the faster transport rates of the EAATs compared to its archaeal homologs.     

Electrostatic analysis of the hepatitis C virus NS3 helicase reveals both active and allosteric site locations

Multi-conformation continuum electrostatics (MCCE) was used to analyze various structures of the NS3 RNA helicase from the hepatitis C virus in order to determine the ionization state of amino acid side chains and their pK_as. In MCCE analyses of HCV helicase structures that lacked ligands, several active site residues were identified to have perturbed pK_as in both the nucleic acid binding site and in the distant ATP-binding site, which regulates helicase movement. In all HCV helicase structures, Glu493 was unusually basic and His369 was abnormally acidic. Both these residues are part of the HCV helicase nucleic acid binding site, and their roles were analyzed by examining the pH profiles of site-directed mutants. Data support the accuracy of MCCE predicted pK_a values, and reveal that Glu493 is critical for low pH enzyme activation. Several key residues, which were previously shown to be involved in helicase-catalyzed ATP hydrolysis, were also identified to have perturbed pK_as including Lys210 in the Walker-A motif and the DExD/H-box motif residues Asp290 and His293. When DNA was present in the structure, the calculated pK_as shifted for both Lys210 and Asp290, demonstrating how DNA binding might lead to electrostatic changes that stimulate ATP hydrolysis.     

Long-range interaction between the enzyme active site and a distant allosteric site in the human mitochondrial NAD(P)-dependent malic enzyme

Our previous study has suggested that mutation of the amino acid residue Asp102 has a significant effect on the fumarate-mediated activation of human mitochondrial NAD(P)⁺-dependent malic enzyme (m-NAD(P)-ME). In this paper, we examine the cationic amino acid residue Arg98, which is adjacent to Asp102 and is highly conserved in most m-NAD(P)-MEs. A series of R98/D102 mutants were created to examine the possible interactions between Arg98 and Asp102 using the double-mutant cycle analysis. Kinetic analysis revealed that the catalytic efficiency of the enzyme was severely affected by mutating both Arg98 and Asp102 residues. However, the binding energy of these mutant enzymes to fumarate as determined by analysis of the K_A,Fum values, show insignificant differences, indicating that the mutation of Arg98 and Asp102 did not cause a significant decrease in the binding affinity of fumarate. The overall coupling energies for R98K/D102N as determined by analysis of the k_cat/K_m and K_A,Fum values were −2.95 and −0.32 kcal/mol, respectively. According to these results, we conclude that substitution of both Arg98 and Asp102 residues has a synergistic effect on the catalytic ability of the enzyme.     

Crystal structure of the O intermediate of the Leu93→Ala mutant of bacteriorhodopsin

The lifetime of the O intermediate of bacteriorhodopsin (BR) is extended by a factor of ～250 in the Leu93‐to‐Ala mutant (BR_L93A). To clarify the structural changes occurring in the last stage of the proton pumping cycle of BR, we crystallized BR_L93A into a hexagonal P622 crystal. Diffraction data from the unphotolyzed state showed that the deletion of three carbon atoms from Leu93 is compensated by the insertion of four water molecules in the cytoplasmic vicinity of retinal. This insertion of water is suggested to be responsible for the blue‐shifted λ_max (540 nm) of the mutant. A long‐lived substate of O with a red‐shifted λ_max (～565 nm) was trapped when the crystal of BR_L93A was flash‐cooled after illumination with green light. This substate (O_slow) bears considerable similarity to the M intermediate of native BR; that is, it commonly shows deformation of helix C and the FG loop, downward orientation of the side chain of Arg82, and disruption of the Glu194/Glu204 pair. In O_slow, however, the main chain of Lys216 is less distorted and retinal takes on the 13‐cis/15‐syn configuration. Another significant difference is seen in the pH dependence of the structure of the proton release group, the pK_a value of which is suggested to be much lower in O_slow than in M. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK_a values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK_a values measured for individual ionizable residues with NMR spectroscopy are actually population‐weighted averages of the pK_a in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK_a values of surface Glu and Asp residues are affected by pH‐sensitive fluctuations of the backbone between folded and locally unfolded conformations. ¹⁵N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps‐ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub‐global unfolding near carboxylic groups. This study has implications for structure‐based pK_a calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH‐sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co‐operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK_a values of ionizable residues in proteins. Proteins 2014; 82:3132–3143. © 2014 Wiley Periodicals, Inc.     

Structural investigation of bacteriorhodopsin and some of its photoproducts by polarized Fourier transform infrared spectroscopic methods-difference spectroscopy and photoselection

The direction of selected IR-transition moments of the retinal chromophore of bacteriorhodopsin (BR) and functional active amino acid residues are determined for light- and dark-adapted BR and for the intermediates K and L of the photocycle. Torsions around single bonds of the chromophore are found to be present in all the investigated BR states. The number of twisted single bonds and the magnitude of these torsions decreases in the order K, L, light-adapted BR, dark-adapted BR. In the last, only the C₁₄—C₁₅ single bond is twisted. The orientation of molecular planes and chemical bonds of such protein side chains, which are perturbed during the transition of light-adapted BR to the respective intermediates, are deduced and the results compared with the current three dimensional model of BR. Trp 86 and Trp 185 are found to form a rigid part of the protein, whereas Asp 96 and Asp 115 perform molecular rearrangements upon formation of the L-intermediate.     

Predicting protein pKa by environment similarity

A statistical method to predict protein pK_a has been developed by using the 3D structure of a protein and a database of 434 experimental protein pK_a values. Each pK_a in the database is associated with a fingerprint that describes the chemical environment around an ionizable residue. A computational tool, MoKaBio, has been developed to identify automatically ionizable residues in a protein, generate fingerprints that describe the chemical environment around such residues, and predict pK_a from the experimental pK_a values in the database by using a similarity metric. The method, which retrieved the pK_a of 429 of the 434 ionizable sites in the database correctly, was crossvalidated by leave‐one‐out and yielded root mean square error (RMSE) = 0.95, a result that is superior to that obtained by using the Null Model (RMSE 1.07) and other well‐established protein pK_a prediction tools. This novel approach is suitable to rationalize protein pK_a by comparing the region around the ionizable site with similar regions whose ionizable site pK_a is known. The pK_a of residues that have a unique environment not represented in the training set cannot be predicted accurately, however, the method offers the advantage of being trainable to increase its predictive power. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Energy cost for the proper ionization of active site residues in 6-phosphogluconate dehydrogenase from T. brucei

The catalytic mechanism of 6-phosphogluconate dehydrogenase requires the inversion of a Lys/Glu couple from its natural ionization state. The pK_a of these residues in free and substrate bound enzymes has been determined measuring by ITC the proton release/uptake induced by substrate binding at different pH values. Wt 6-phosphogluconate dehydrogenase from Trypanosoma brucei and two active site enzyme mutants, K185H and E192Q were investigated. Substrate binding was accompanied by proton release and was dependent on the ionization of a group with pK_a 7.07 which was absent in the E192Q mutant. Kinetic data highlighted two pK_a, 7.17 and 9.64, in the enzyme–substrate complex, the latter being absent in the E192Q mutant, suggesting that the substrate binding shifts Glu192 pK_a from 7.07 to 9.64. A comparison of wt and E192Q mutant appears to show that the substrate binding shifts Lys185 pK_a from 9.9 to 7.17. By comparing differences in proton release and the binding enthalpy of wt and mutant enzymes, the enthalpic cost of the change in the protonation state of Lys185 and Glu192 was estimated at ≈ 6.1 kcal/mol. The change in protonation state of Lys185 and Glu192 has little effect on Gibbs free energy, 240–325 cal/mol. However proton balance evidences the dissociation of other group(s) that can be collectively described by a single pK_a shift from 9.1 to 7.54. This further change in ionization state of the enzyme causes an increase of free energy with a total cost of 1.2–2.3 kcal/mol to set the enzyme into a catalytically competent form.     

pH-rate profiles support a general base mechanism for galactokinase (Lactococcus lactis)

Galactokinase (GALK), a member the Leloir pathway for normal galactose metabolism, catalyzes the conversion of α-d-galactose to galactose-1-phosphate. For this investigation, we studied the kinetic mechanism and pH profiles of the enzyme from Lactococcus lactis. Our results show that the mechanism for its reaction is sequential in both directions. Mutant proteins D183A and D183N are inactive (<10 000 fold), supporting the role of Asp183 as a catalytic base that deprotonates the C-1 hydroxyl group of galactose. The pH-k_cat profile of the forward reaction has a pK_a of 6.9 ± 0.2 that likely is due to Asp183. The pH-k_cat/K_Gal profile of the reverse reaction further substantiates this role as it is lacking a key pK_a required for a direct proton transfer mechanism. The R36A and R36N mutant proteins show over 100-fold lower activity than that for the wild-type enzyme, thus suggesting that Arg36 lowers the pK_a of the C-1 hydroxyl to facilitate deprotonation.     

The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Theoretical studies on the redox-Bohr effect in cytochrome c 3 from Desulfovibrio vulgaris Hildenborough

 The pH dependence of the redox potentials in the tetrahemic cytochrome c ₃ from Desulfovibrio vulgaris Hildenborough (redox-Bohr effect) is here investigated using continuum electrostatics methods. The redox-Bohr effect seems to be associated with changes in the protonation state of charged residues in the protein, but the exact residues had not been identified. The global pK _a of this phenomenon is dependent on the redox state of the molecule, and the influence of the pH on the microscopic potential of each heme has been experimentally quantified. The availability of detailed experimental data provides us with important and unique guides to the performance of ab initio pK _a calculations aiming at the identification of the groups involved. These calculations were performed in several redox states along the reduction pathway, with the double objective of finding groups with redox-linked pK _a shifts, and absolute pK _as compatible with the redox-Bohr effect. The group with the largest pK _a shift along the reduction pathway is propionate D from heme I. Its effect on the redox potential of individual hemes, as calculated by electrostatic calculations, correlates very well with the experimental order of influence, making it a likely candidate. Abnormal titration of the same propionate has been experimentally observed on a homologous cytochrome c ₃ from a different strain, thus strengthening the theoretical result. However, its absolute calculated pK _a in the fully oxidised cytochrome is outside the zone where the phenomenon is known to occur, but the calculation shows a strong dependence on small conformational changes, suggesting large uncertainties in the calculated value. A group with a pK _a value within the experimentally observed range is propionate D from heme IV. Its influence on the redox potential of the hemes does not correlate with the experimental order, indicating that, although it may be one of the possible players on the phenomenon, it cannot be solely responsible for it. Mutation of the Lys45 residue is suggested as an indirect way of probing the importance of the propionate D from heme I in the mechanism. Non-heme groups may also be involved in this process; our calculations indicate His67 and the N-terminal as groups that may play a role. Accuracy and applicability of current continuum electrostatic methods are discussed in the context of this system. Received: 27 March 1997 / Accepted: 19 August 1997     

Crystal structure of the eukaryotic light-driven proton-pumping rhodopsin, Acetabularia rhodopsin II, from marine alga

Acetabularia rhodopsin (AR) is a rhodopsin from the marine plant Acetabularia acetabulum. The opsin-encoding gene from A. acetabulum, ARII, was cloned and found to be novel but homologous to that reported previously. ARII is a light-driven proton pump, as demonstrated by the existence of a photo-induced current through Xenopus oocytes expressing ARII. The photochemical reaction of ARII prepared by cell-free protein synthesis was similar to that of bacteriorhodopsin (BR), except for the lack of light-dark adaptation and the different proton release and uptake sequence. The crystal structure determined at 3.2 Å resolution is the first structure of a eukaryotic member of the microbial rhodopsin family. The structure of ARII is similar to that of BR. From the cytoplasmic side to the extracellular side of the proton transfer pathway in ARII, Asp92, a Schiff base, Asp207, Asp81, Arg78, Glu199, and Ser189 are arranged in positions similar to those of the corresponding residues directly involved in proton transfer by BR. The side-chain carboxyl group of Asp92 appears to interact with the sulfhydryl group of Cys218, which is unique to ARII and corresponds to Leu223 of BR and to Asp217 of Anabaena sensory rhodopsin. The orientation of the Arg78 side chain is opposite to the corresponding Arg82 of BR. The putative absence of water molecules around Glu199 and Arg78 may disrupt the formation of the low-barrier hydrogen bond at Glu199, resulting in the “late proton release”.