Molecular dynamics simulations of bovine rhodopsin: influence of protonation states and different membrane-mimicking environments

G-protein coupled receptors (GPCRs) are a protein family of outstanding pharmaceutical interest. GPCR homology models, based on the crystal structure of bovine rhodopsin, have been shown to be valuable tools in the drug-design process. The initial model is often refined by molecular dynamics (MD) simulations, a procedure that has been recently discussed controversially. We therefore analyzed MD simulations of bovine rhodopsin in order to identify contacts that could serve as constraints in the simulation of homology models. Additionally, the effect of an N-terminal truncation, the nature of the membrane mimic, the influence of varying protonation states of buried residues and the importance of internal water molecules was analyzed. All simulations were carried out using the program-package GROMACS. While N-terminal truncation negatively influenced the overall protein stability, a stable simulation was possible in both solvent environments. As regards the protonation state of titratable sites, the experimental data could be reproduced by the program UHBD (University of Houston Brownian Dynamics), suggesting its application for studying homology models of GPCRs. A high flexibility was observed for internal water molecules at some sites. Finally, interhelical hydrogen-bonding interactions could be derived, which can now serve as constraints in the simulations of GPCR homology models.     

Asp83, Glu113 and Glu134 are not specifically involved in Schiff base protonation or wavelength regulation in bovine rhodopsin

Site-specific mutagenesis was employed to investigate the proposed contribution of proton-donating residues (Glu, Asp) in the membrane domains of bovine rhodopsin to protonation of the Schiff base-linking protein and chromophore or to wavelength modulation of this visual pigment. Three point-mutations were introduced to replace the highly conserved residues Asp83 by Asn (D83N), Glu113 by Gln (E113 Q) or Glu134 by Asp (E134D), respectively. All 3 substitutions had only marginal effects on the spectral properties of the final pigment (less than or equal to 3 nm blue-shift relative to native rhodopsin). Hence, none of these residues by itself is specifically involved in Schiff base protonation or wavelength modulation of bovine rhodopsin.     

Simulating proteins at constant pH: An approach combining molecular dynamics and Monte Carlo simulation

For the structure and function of proteins, the pH of the solution is one of the determining parameters. Current molecular dynamics (MD) simulations account for the solution pH only in a limited way by keeping each titratable site in a chosen protonation state. We present an algorithm that generates trajectories at a Boltzmann distributed ensemble of protonation states by a combination of MD and Monte Carlo (MC) simulation. The algorithm is useful for pH-dependent structural studies and to investigate in detail the titration behavior of proteins. The method is tested on the acidic residues of the protein hen egg white lysozyme. It is shown that small structural changes may have a big effect on the pK(A) values of titratable residues.     

Protonation changes upon ligand binding to trypsin and thrombin: structural interpretation based on pK(a) calculations and ITC experiments

The protonation states of a protein and a ligand can be altered upon complex formation. Such changes can be detected experimentally by isothermal titration calorimetry (ITC). For a series of ligands binding to the serine proteases trypsin and thrombin, we previously performed an extensive ITC and crystallographic study and were able to identify protonation changes for four complexes. However, since ITC measures only the overall proton exchange, it does not provide structural insights into the functional groups involved in the proton transfer. Using Poisson-Boltzmann calculations based on our recently developed PEOE_PB charges, we compute pK(a) values for all complexes of our former study in order to reveal the residues with altered protonation states. The results indicate that His57, a member of the catalytic triad, is responsible for the most relevant pK(a) shifts leading to the experimentally detected protonation changes. This finding is in contrast to our previous assumption that the observed protonation changes occur at the carboxylic group of the ligands. The newly detected proton acceptor is used for a revised factorization of the ITC data, which is necessary whenever the protonation inventory changes upon complexation. The pK(a) values of complexes showing no protonation change in the ITC experiment are reliably predicted in most cases, whereas predictions of strongly coupled systems remain problematic.     

Long-range nature of the interactions between titratable groups in Bacillus agaradhaerens family 11 xylanase: pH titration of B. agaradhaerens xylanase

Xylanase from Bacillus agaradhaerens belongs to a large group of glycosyl hydrolases which catalyze the degradation of xylan. The protonation behavior of titratable groups of the uniformly (15)N- and (13)C-labeled xylanase was investigated by multinuclear NMR spectroscopy. A total of 224 chemical shift titration curves corresponding to (1)H, (13)C, and (15)N resonances revealed pK(a) values for all aspartic and glutamic acid residues, as well as for the C-terminal carboxylate and histidine residues. Most of the titratable groups exhibit a complex titration behavior, which is most likely due to the mutual interactions with other neighboring groups or due to an unusual local microenvironment. Subsite -1 containing the catalytic dyad shows a long-range interaction over 9 A with Asp21 via two hydrogen bonds with Asn45 as the mediator. This result illuminates the pivotal role of the conserved position 45 among family 11 endoxylanases, determining an alkaline pH optimum by asparagine residues or an acidic pH optimum by an aspartate. The asymmetric interactions of neighboring tryptophan side chains with respect to the catalytic dyad can be comprehended as a result of hydrogen bonding and aromatic stacking. Most of the chemical shift-pH profiles of the backbone amides exhibit biphasic behavior with two distinct inflection points, which correspond to the pK(a) values of the nearby acidic side chains. However, the alternation of both positive and negative slopes of individual amide titration curves is interpreted as a consequence of a simultaneous reorganization of side chain conformational space at pH approximately 6 and/or an overall change in the hydrogen network in the substrate binding cleft.     

Structure and pH sensitivity of the transmembrane segment 3 of rhodopsin

Activation of G protein-coupled receptors (GPCRs) originates in ligand-induced protein conformational changes that are transmitted to the cytosolic receptor surface. In the photoreceptor rhodopsin, and possibly other rhodopsin-like GPCRs, protonation of a carboxylic acid in the conserved E(D)RY motif at the cytosolic end of transmembrane helix 3 (TM3) is coupled to receptor activation. Here, we have investigated the structure of synthetic peptides derived from rhodopsin TM3. Polarized FTIR spectroscopy reveals a helical structure of a 31-mer TM3 peptide reconstituted into PC vesicles with a large tilt of 40-50 degrees of the helical axis relative to the membrane normal. Helical structure is also observed for the TM3 peptide in detergent micelles and depends on pH, especially in the C-terminal sequence. In addition, the fluorescence emission of the single tyrosine of the D(E)RY motif in the TM3 peptide exhibits a pronounced pH sensitivity that is abolished when Glu is replaced by Gln, demonstrating that protonation of the conserved Glu side chain affects the structure in the environment of the D(E)RY motif of TM3. The pH regulation of the C-terminal TM3 structure may be an intrinsic feature of the E(D)RY motif in other class I receptors, allowing the coupling of protonation and conformation of membrane-exposed residues in full-length GPCRs.     

Evidence for the rate of the final step in the bacteriorhodopsin photocycle being controlled by the proton release group: R134H mutant

Light absorbed by bacteriorhodopsin (bR) leads to a proton being released at the extracellular surface of the purple membrane. Structural studies as well as studies of mutants of bR indicate that several groups form a pathway for proton transfer from the Schiff base to the extracellular surface. These groups include D85, R82, E204, E194, and water molecules. Other residues may be important in tuning the initial state pK(a) values of these groups and in mediating light-induced changes of the pK(a) values. A potentially important residue is R134: it is located close to E194 and might interact electrostatically to affect the pK(a) of E194 and light-induced proton release. In this study we investigated effects of the substitution of R134 with a histidine on light-induced proton release and on the photocycle transitions associated with proton transfer. By measuring the light-induced absorption changes versus pH, we found that the R134H mutation results in an increase in the pK(a) of the proton release group in both the M (0.6 pK unit) and O (0.7 pK unit) intermediate states. This indicates the importance of R134 in tuning the pK(a) of the group that, at neutral and high pH, releases the proton upon M formation (fast proton release) and that, at low pH, releases the proton simultaneously with O decay (slow proton release). The higher pK(a) of the proton release group found in R134H correlates with the slowing of the rate of the O --> bR transition at low pH and probably is the cause of this slowing. The pH dependence of the fraction of the O intermediate is altered in R134H compared to the WT but is similar to that in the E194D mutant: a very small amount of O is present at neutral pH, but the fraction of O increases greatly upon decreasing the pH. These results provide further support for the hypothesis that the O --> bR transition is controlled by the rate of deprotonation of the proton release group. These data also provide further evidence for the importance of the R134-E194 interaction in modulating proton release from D85 after light has led to its being protonated.     

Protein-protein binding is often associated with changes in protonation state

pK(a) values of ionizable residues have been calculated using the PROPKA method and structures of 75 protein-protein complexes and their corresponding free forms. These pK(a) values were used to compute changes in protonation state of individual residues, net changes in protonation state of the complex relative to the uncomplexed proteins, and the correction to a binding energy calculated assuming standard protonation states at pH 7. For each complex, two different structures for the uncomplexed form of the proteins were used: the X-ray structures determined for the proteins in the absence of the other protein and the individual protein structures taken from the structure of the complex (referred to as unbound and bound structures, respectively). In 28 and 77% of the cases considered here, protein-protein binding is accompanied by a complete (>95%) or significant (>50%) change in protonation state of at least one residue using unbound structures. Furthermore, in 36 and 61% of the cases, protein-protein binding is accompanied by a complete or significant net change in protonation state of the complex relative to the separated monomers. Using bound structures, the corresponding values are 12, 51, 20, and 48%. Comparison to experimental data suggest that using unbound and bound structures lead to over- and underestimation of binding-induced protonation state changes, respectively. Thus, we conclude that protein-protein binding is often associated with changes in protonation state of amino acid residues and with changes in the net protonation state of the proteins. The pH-dependent correction to the binding energy contributes at least one order of magnitude to the binding constant in 45 and 23%, using unbound and bound structures, respectively.     

Functional role of the "ionic lock"--an interhelical hydrogen-bond network in family A heptahelical receptors

Activation of family A G-protein-coupled receptors involves a rearrangement of a conserved interhelical cytoplasmic hydrogen bond network between the E(D)RY motif on transmembrane helix 3 (H3) and residues on H6, which is commonly termed the cytoplasmic “ionic lock.” Glu134^3.49 of the E(D)RY motif also forms an intrahelical salt bridge with neighboring Arg135^3.50 in the dark-state crystal structure of rhodopsin. We examined the roles of Glu134^3.49 and Arg135^3.50 on H3 and Glu247^6.30 and Glu249^6.32 on H6 on the activation of rhodopsin using Fourier transform infrared spectroscopy of wild-type and mutant pigments reconstituted into lipid membranes. Activation of rhodopsin is pH-dependent with proton uptake during the transition from the inactive Meta I to the active Meta II state. Glu134^3.49 of the ERY motif is identified as the proton-accepting group, using the Fourier transform infrared protonation signature and the absence of a pH dependence of activation in the E134Q mutant. Neutralization of Arg135^3.50 similarly leads to pH-independent receptor activation, but with structural alterations in the Meta II state. Neutralization of Glu247^6.30 and Glu249^6.32 on H6, which are involved in interhelical interactions with H3 and H7, respectively, led to a shift toward Meta II in the E247Q and E249Q mutants while retaining the pH sensitivity of the equilibrium. Disruption of the interhelical interaction of Glu247^6.30 and Glu249^6.32 on H6 with H3 and H7 plays its role during receptor activation, but neutralization of the intrahelical salt bridge between Glu134^3.49 and Arg135^3.50 is considerably more critical for shifting the photoproduct equilibrium to the active conformation. These conclusions are discussed in the context of recent structural data of the β₂-adrenergic receptor.     

pK(a) calculations along a bacteriorhodopsin molecular dynamics trajectory

Electrostatic calculations of pK(a-values) are reported along a 400 ps molecular dynamics trajectory of bacteriorhodopsin. The sensitivity of calculated pK(a) values to a number of structural factors and factors related to the modelling of the electrostatics are also studied. The results are very sensitive to the choice of internal dielectric constant of the protein (in the interval 2-4). Moreover it is important to include internal water molecules and to average over a long enough portion ( approximately 100 ps) of an equilibrium molecular dynamics trajectory. The internal waters are necessary to get an ion-counter ion complex with the Schiff base and Arg 82 protonated and the aspartic groups (85 and 212) deprotonated. The fluctuations along the MD-trajectory do not change the protonation state of internal residues at neutral pH. However, at other pH values the averaging along a trajectory maybe crucial to get correct protonation states. A relationship is found between the arginine group 82, the aspartic group 85 and the glutamate group 204. Glu 204 is protonated in the ground state but the pK(a) value decreases towards deprotonation when the chromophore isomerizes into the cis state.     

On the Hill plot of NMR data for titration of protein residues

The Hill plots of NMR titration data for protein residues disclose more clearly than the usual titration curves the presence of multiple weak perturbations originating from other titratable groups, and should be used whenever the conventional curve fitting is poor. For a quantitative interpretation, we derive here expressions for the Hill equation and the Hill coefficient when the titration of the observed group is perturbed by more than one titratable group. When the generalized Hill equation is fitted to the data, values of the interaction parameters between the observed group and the others are extracted provided that there are no mutual interactions between the latter groups. The method is applied to the titration data of two histidyl residues of l-arginine phosphotransferase (E.C. 2.7.3.3.) in the transition state analogue complex (enzyme-Mg²⁺-ADP-NOsk3/–l-Arg). From the Hill plots, interactions with three titratable groups are disclosed for both residues, and the fitting with the Hill equation reveals that they experience perturbations from the same three groups. Microscopic pK values are obtained for all the involved groups, indicating large changes (up to 3 pH units) upon protonation of the interacting groups. As compared to the conventional fitting procedure, the use and fitting of Hill plots yields from NMR data more information on the neighbourhood of enzyme residues and on the changes intervening therein through the steps involved in the catalysis.     

Protonation state of Asp30 exerts crucial influence over surface loop rearrangements responsible for NO release in nitrophorin 4

pK(a) values of ionizable residues were calculated for the crystal structures describing the pH and NO binding dependant conformations of nitrophorin 4, a pH sensitive NO carrier heme protein. Comparison of resultant H-bonding patterns allowed the identification of the amino acids that take part in signaling pH change. We carried out MD simulations to show that the protonation state of Asp30, buried in the closed conformation, is crucial for maintaining the tight packed conformation of the closed form of the complex - presenting a model for the functional decrease of NO binding affinity of nitrophorins at physiological pH.     

The ionization state of D37 in E. coli porin OmpF and the nature of conductance fluctuations in D37 mutants

Understanding the flow of ions through E. coli porin outer membrane protein F (OmpF) requires knowledge of the charge state of all titratable residues located along the permeation pathway. Earlier theoretical studies proved successful in the calculation of the pK values of most residues. The (apparent) pK of Asp37 (D37), on the other hand, appeared rather sensitive to the (unknown) protein dielectric used. We addressed the protonation state of D37 experimentally by replacing D37 with a (neutral) valine. This D37V mutant expressed reduced cation selectivity, in agreement with the view that D37 in wild-type (WT) OmpF is fully ionized, i.e., deprotonated. The introduction of a (positively charged) arginine at position 37 evoked current fluctuations. Similar behavior was observed in the D37K mutant and the cysteine mutants D37C-MTSEA and D37C-MTSET. Nontitratable [2-(trimethylammonium)ethyl]-methanethiosulfonate (MTSET) carries a permanent and pH-independent charge of 1e, implying that the fluctuations of the D37C-MTSET mutant do not represent (de)protonation reactions of MTSET. We therefore conclude that these fluctuations reflect transitions between conformational substates evoked by structural instabilities due to the positive charge at that particular position in the pore lumen. Based on the similarities between D37C-MTSET fluctuations and those seen in the other mutants, notably D37K, the underlying mechanism of these fluctuations may be (essentially) the same in all four mutants studied.     

Deactivation and proton transfer in light-induced metarhodopsin II/metarhodopsin III conversion: a time-resolved fourier transform infrared spectroscopic study

Vertebrate rhodopsin shares with other retinal proteins the 11-cis-retinal chromophore and the light-induced 11-cis/trans isomerization triggering its activation pathway. However, only in rhodopsin the retinylidene Schiff base bond to the apoprotein is eventually hydrolyzed, making a complex regeneration pathway necessary. Metabolic regeneration cannot be short-cut, and light absorption in the active metarhodopsin (Meta) II intermediate causes anti/syn isomerization around the retinylidene linkage rather than reversed trans/cis isomerization. A new deactivating pathway is thereby triggered, which ends in the Meta III "retinal storage" product. Using time-resolved Fourier transform infrared spectroscopy, we show that the identified steps of receptor activation, including Schiff base deprotonation, protein structural changes, and proton uptake by the apoprotein, are all reversed. However, Schiff base reprotonation is much faster than the activating deprotonation, whereas the protein structural changes are slower. The final proton release occurs with pK approximately 4.5, similar to the pK of a free Glu residue and to the pK at which the isolated opsin apoprotein becomes active. A forced deprotonation, equivalent to the forced protonation in the activating pathway, which occurs against the unfavorable pH of the medium, is not observed. This explains properties of the final Meta III product, which displays much higher residual activity and is less stable than rhodopsin arising from regeneration with 11-cis-retinal. We propose that the anti/syn conversion can only induce a fast reorientation and distance change of the Schiff base but fails to build up the full set of dark ground state constraints, presumably involving the Glu(134)/Arg(135) cluster.     

Influence of a lipid interface on protein dynamics in a fungal lipase.

Lipases catalyze lipolytic reactions and for optimal activity they require a lipid interface. To study the effect of a lipid aggregate on the behavior of the enzyme at the interfacial plane and how the aggregate influences an attached substrate or product molecule in time and space, we have performed molecular dynamics simulations. The simulations were performed over 1 to 2 ns using explicit SPC water. The interaction energies between protein and lipid are mainly due to van der Waals contributions reflecting the hydrophobic nature of the lipid molecules. Estimations of the protonation state of titratable residues indicated that the negative charge on the fatty acid is stabilized by interactions with the titratable residues Tyr-28, His-143, and His-257. In the presence of a lipid patch, the active site lid opens wider than observed in the corresponding simulations in an aqueous environment. In that lid conformation, the hydrophobic residues Ile-85, Ile-89, and Leu-92 are embedded in the lipid patch. The behavior of the substrate or product molecule is sensitive to the environment. Entering and leaving of substrate molecules could be observed in presence of the lipid patch, whereas the product forms strong hydrogen bonds with Ser-82, Ser-144, and Trp-88, suggesting that the formation of hydrogen bonds may be an important contribution to the mechanism by which product inhibition might take place.     

pK(a) calculations for class C beta-lactamases: the role of Tyr-150

The Poisson-Boltzmann method was used to compute the pK(a) values of titratable residues in a set of class C beta-lactamases. In these calculations, the pK(a) of the phenolic group of residue Tyr150 is the only one to stand out with an abnormally low value of 8.3, more than one pK(a) unit lower than the measured reference value for tyrosine in solution. Other important residues of the catalytic pocket, such as the conserved Lys67, Lys315, His314, and Glu272 (hydrogen-bonded to the ammonium group of Lys315), display normal protonation states at neutral pH. pK(a) values were also computed in catalytically impaired beta-lactamase mutants. Comparisons between the relative k(cat) values and the Tyr150 pK(a) value in these mutants revealed a striking correlation. In active enzymes, this pK(a) value is always lower than the solution reference value while it is close to normal in inactive enzymes. These results thus support the hypothesis that the phenolate form of Tyr150 is responsible for the activation of the nucleophilic serine. The possible roles of Lys67 and Lys315 during catalysis are also discussed.     

Oxidation of the non-heme iron complex in photosystem II

In photosystem II (PSII), the redox properties of the non-heme iron complex (Fe complex) are sensitive to the redox state of quinones (Q(A/)(B)), which may relate to the electron/proton transfer. We calculated the redox potentials for one-electron oxidation of the Fe complex in PSII [E(m)(Fe)] based on the reference value E(m)(Fe) = +400 mV at pH 7 in the Q(A)(0)Q(B)(0) state, considering the protein environment in atomic detail and the associated changes in protonation pattern. Our model yields the pH dependence of E(m)(Fe) with -60 mV/pH as observed in experimental redox titration. We observed significant deprotonation at D1-Glu244 in the hydrophilic loop region upon Fe complex oxidation. The calculated pK(a) value for D1-Glu244 depends on the Fe complex redox state, yielding a pK(a) of 7.5 and 5.5 for Fe(2+) and Fe(3+), respectively. To account for the pH dependence of E(m)(Fe), a model involving not only D1-Glu244 but also the other titratable residues (five Glu in the D-de loops and six basic residues near the Fe complex) seems to be needed, implying the existence of a network of residues serving as an internal proton reservoir. Reduction of Q(A/B) yields +302 mV and +268 mV for E(m)(Fe) in the Q(A)(-)Q(B)(0) and Q(A)(0)Q(B)(-) states, respectively. Upon formation of the Q(A)(0)Q(B)(-) state, D1-His252 becomes protonated. Forming Fe(3+)Q(B)H(2) by a proton-coupled electron transfer process from the initial state Fe(2+)Q(B)(-) results in deprotonation of D1-His252. The two EPR signals observed at g = 1.82 and g = 1.9 in the Fe(2+)Q(A)(-) state of PSII may be attributed to D1-His252 with variable and fixed protonation, respectively.     

Effect of divalent cations and cytosine protonation on thermodynamic properties of intermolecular DNA double and triple helices

The contribution of divalent cations and cytosine protonation to conformation and stability of duplex and triplex formation were intensively investigated and characterized by ultraviolet (UV), circular dichroism (CD), differential scanning calorimetry (DSC), and electrophoresis mobility shift assay (EMSA). CD spectra showed that the divalent cations investigated would not significantly distort nucleotide geometry, while UV and DSC melting experiments revealed that the cation binding abilities to duplexes and triplexes were clearly dependent on the types of cations under near physiological conditions. The calorimetric enthalpies were generally underestimated relative to the corresponding van't Hoff enthalpies for Hoogsteen and Watson-Crick transitions, but free energy changes derived from the DSC measurements were in good agreement with those derived from the UV measurements. The adjacent placing of the C(+) x G.C triplets in triplexes lowered the stabilities of not only Hoogsteen base-pairing but also Watson-Crick base-pairing. The protonation contribution of the given cytosine residues might depend on the local and global structure of the protonated cytosine complex. A rigid structural targeted-strand would favor the protonation of cytosine residues. The apparent pK(a) values for parallel duplex and triplex investigated were determined to be 6.4 and 7.6, respectively, which are considerably heightened by 2.1 and 3.3 pH unit as compared to the intrinsic pK(a) value of the free cytosine residues.     

Role of Asp193 in chromophore-protein interaction of pharaonis phoborhodopsin (sensory rhodopsin II)

Pharaonis phoborhodopsin (ppR; also pharaonis sensory rhodopsin II, psRII) is a receptor of the negative phototaxis of Natronobacterium pharaonis. By spectroscopic titration of D193N and D193E mutants, the pK(a) of the Schiff base was evaluated. Asp193 corresponds to Glu204 of bacteriorhodopsin (bR). The pK(a) of the Schiff base (SBH(+)) of D193N was approximately 10.1-10.0 (at XH(+)) and approximately 11.4-11.6 (at X) depending on the protonation state of a certain residue (designated by X) and independent of Cl(-), whereas those of the wild type and D193E were >12. The pK(a) values of XH(+) were approximately 11.8-11.2 at the state of SB, 10.5 at SBH(+) state in the presence of Cl(-), and 9.6 at SBH(+) without Cl(-). These imply the presence of a long-range interaction in the extracellular channel. Asp193 was suggested to be deprotonated in the present dodecyl-maltoside (DDM) solubilized wild-type ppR, which is contrary to Glu204 of bR. In the absence of salts, the irreversible denaturation of D193N (but not the wild type and D193E) occurred via a metastable state, into which the addition of Cl(-) reversed the intact pigment. This suggests that the negative charge at residue 193, which can be substituted by Cl(-), is necessary to maintain the proper conformation in the DDM-solubilized ppR.