5-Enolpyruvylshikimate-3-phosphate synthase: determination of the protonation state of active site residues by the semiempirical method

EPSP synthase (EPSPS) catalyzes the addition of shikimate-3-phosphate (S3P) and phosphoenolpyruvate (PEP) to form a tetrahedral intermediate (TI) that is converted to 5-enolpyruvylshikimate-3-phosphate (EPSP) and inorganic phosphate. A semiempirical molecular modeling study of the EPSPS active site containing the TI was implemented for the assignment of the protonation states of four basic residues, Lys22, Lys340, His385, and Lys411, based on the evaluation of 16 different protonation states and comparison of the resulting energy minimized heavy atoms coordinates with available X-ray crystallographic data of the D313A mutant of EPSPS. The results, employing both gas phase and continuum solvent models, are indicative that after the TI formation the histidine residue is most probably in neutral form (N^ε-protonated) and the lysine residues are in protonated form, which suggests that none of the presently proposed assignments of aminoacid residues involved in the reaction mechanism could be completely correct. The protonated state of Lys22 in the presence of the TI supports the proposal that this residue is a general acid catalyst for TI breakdown. Modeling of the native enzyme active site suggests that Asp313 residue has only minor effects on the definition of the TI position inside the active site. Hydrogen-bonds distances suggest that, in order to act as a base, Asp313 needs the intermediacy of a hydroxyl group of the TI for effecting the attack on the TI methyl group in the elimination step leading to EPSP, as suggested previously in the literature.     

Direct Determination of Protonation States of Histidine Residues in a 2 Å Neutron Structure of Deoxy-Human Normal Adult Hemoglobin and Implications for the Bohr Effect

We have investigated the protonation states of histidine residues (potential Bohr groups) in the deoxy form (T state) of human hemoglobin by direct determination of hydrogen (deuterium) positions with the neutron protein crystallography technique. The reversible binding of protons is key to the allosteric regulation of human hemoglobin. The protonation states of 35 of the 38 His residues were directly determined from neutron scattering omit maps, with 3 of the remaining residues being disordered. Protonation states of 5 equivalent His residues—αHis20, αHis50, αHis89, βHis143, and βHis146—differ between the symmetry-related globin subunits. The distal His residues, αHis58 and βHis63, are protonated in the α1β1 heterodimer and are neutral in α2β2. Buried residue αHis103 is found to be protonated in both subunits. These distal and buried residues have the potential to act as Bohr groups. The observed protonation states of His residues are compared to changes in their pK_a values during the transition from the T to the R state and the results provide some new insights into our understanding of the molecular mechanism of the Bohr effect.     

Neutron structure of the T26H mutant of T4 phage lysozyme provides insight into the catalytic activity of the mutant enzyme and how it differs from that of wild type

T4 phage lysozyme is an inverting glycoside hydrolase that degrades the murein of bacterial cell walls by cleaving the β‐1,4‐glycosidic bond. The substitution of the catalytic Thr26 residue to a histidine converts the wild type from an inverting to a retaining enzyme, which implies that the original general acid Glu11 can also act as an acid/base catalyst in the hydrolysis. Here, we have determined the neutron structure of the perdeuterated T26H mutant to clarify the protonation states of Glu11 and the substituted His26, which are key in the retaining reaction. The 2.09‐Å resolution structure shows that the imidazole group of His26 is in its singly protonated form in the active site, suggesting that the deprotonated N?2 atom of His26 can attack the anomeric carbon of bound substrate as a nucleophile. The carboxyl group of Glu11 is partially protonated and interacts with the unusual neutral state of the guanidine moiety of Arg145, as well as two heavy water molecules. Considering that one of the water‐binding sites has the potential to be occupied by a hydronium ion, the bulk solvent could be the source for the protonation of Glu11. The respective protonation states of Glu11 and His26 are consistent with the bond lengths determined by an unrestrained refinement of the high‐resolution X‐ray structure of T26H at 1.04‐Å resolution. The detail structural information, including the coordinates of the deuterium atoms in the active site, provides insight into the distinctively different catalytic activities of the mutant and wild type enzymes.     

Binding of ruthenium(III) anti-tumor drugs to human lactoferrin probed by high resolution X-ray crystallographic structure analyses

The binding to human lactoferrin of three Ru(III) complexes with anti-tumor activity has been investigated by X-ray crystallography in order to gain insights into how such complexes might be carried during transferrin-mediated delivery to cells. The complexes, HIm[RuIm₂Cl₄], HInd[RuInd₂Cl₄] and (HInd)₂ [RuIndCl₅], where Im?=?imidazole and Ind?=?indazole, were diffused into crystals of apo-lactoferrin (apoLf). X-ray diffraction data were collected to 2.6?Å, 2.2?Å and 2.4?Å respectively. The binding sites for the Ru complexes were determined from difference Fouriers, in comparison with native apoLf; the two indazole-apoLf complexes were also refined crystallographically to final R factors of 0.202 (for 8.0 to 2.3?Å data) and 0.192 (for 8.0 to 2.4?Å data) respectively. Two types of binding site were identified, a high-affinity site at His 253 in the open N-lobe iron-binding cleft of apoLf (and by analogy a similar one at His 597 in the C-lobe), and lower-affinity sites at surface-exposed His residues, primarily His 590 and His 654. The exogenous heterocyclic ligands remain bound to Ru, at least at the His 253 site, and modelling suggests that the nature and number of these ligands may determine whether the closed structure that is required for receptor binding could be formed or not. The results also highlight the importance of His residues for binding such complexes and the value of heavy atom binding studies from crystallographic analyses for identifying non-specific binding sites on proteins.     

Electron acceptors of bacterial photosynthetic reaction centers II. H+ binding coupled to secondary electron transfer in the quinone acceptor complex

The photoreduction of ubiquinone in the electron acceptor complex (Q₁Q₁₁) of photosynthetic reaction centers from Rhodopseudomonas sphaeroides, R26, was studied in a series of short, saturating flashes. The specific involvement of H⁺ in the reduction was revealed by the pH dependence of the electron transfer events and by net H⁺ binding during the formation of ubiquinol, which requires two turnovers of the photochemical act. On the first flash Q₁₁ receives an electron via Q₁ to form a stable ubisemiquinone anion (Q?^?₁₁); the second flash generates Q?^?₁. At low pH the two semiquinones rapidly disproportionate with the uptake of 2 H⁺, to produce Q₁₁H₂. This yields out-of-phase binary oscillations for the formation of anionic semiquinone and for H⁺ uptake. Above pH 6 there is a progressive increase in H⁺ binding on the first flash and an equivalent decrease in binding on the second flash until, at about pH 9.5, the extent of H⁺ binding is the same on all flashes. The semiquinone oscillations, however, are undiminished up to pH 9. It is suggested that a non-chromophoric, acid-base group undergoes a pK shift in response to the appearance of the anionic semiquinone and that this group is the site of protonation on the first flash. The acid-base group, which may be in the reaction center protein, appears to be subsequently involved in the protonation events leading to fully reduced ubiquinol. The other proton in the two electron reduction of ubiquinone is always taken up on the second flash and is bound directly to Q?^?₁₁. At pH values above 8.0, it is rate limiting for the disproportionation and the kinetics, which are diffusion controlled, are properly responsive to the prevailing pH. Below pH 8, however, a further step in the reaction mechanism was shown to be rate limiting for both H⁺ binding electron transfer following the second flash.     

Protonated Form: The Potent Form of Potassium-Competitive Acid Blockers

Potassium-competitive acid blockers (P-CABs) are highly safe and active drugs targeting H⁺,K⁺-ATPase to cure acid-related gastric diseases. In this study, we for the first time investigate the interaction mechanism between the protonated form of P-CABs and human H⁺,K⁺-ATPase using homology modeling, molecular docking, molecular dynamics and binding free energy calculation methods. The results explain why P-CABs have higher activities with higher pKa values or at lower pH. With positive charge, the protonated forms of P-CABs have more competitive advantage to block potassium ion into luminal channel and to bind with H⁺,K⁺-ATPase via electrostatic interactions. The binding affinity of the protonated form is more favorable than that of the neutral P-CABs. In particular, Asp139 should be a very important binding site for the protonated form of P-CABs through hydrogen bonds and electrostatic interactions. These findings could promote the rational design of novel P-CABs.     

A goose-type lysozyme from ostrich (Struthio camelus) egg white: multiple roles of His101 in its enzymatic reaction

A goose-type lysozyme from ostrich egg white (OEL) was produced by Escherichia coli expression system, and the role of His101 of OEL in the enzymatic reaction was investigated by NMR spectroscopy, thermal unfolding, and theoretical modeling of the enzymatic hydrolysis of hexa-N-acetylchitohexaose, (GlcNAc)₆. Although the binding of tri-N-acetylchitotriose, (GlcNAc)₃, to OEL perturbed several backbone resonances in the ¹H–¹⁵N HSQC spectrum, the chemical shift of the backbone resonance of His101 was not significantly affected. However, apparent pK_a values of His101 and Lys102 determined from the pH titration curves of the backbone chemical shifts were markedly shifted by (GlcNAc)₃ binding. Thermal unfolding experiments and modeling study of (GlcNAc)₆ hydrolysis using a His101-mutated OEL (H101A-OEL) revealed that the His101 mutation affected not only sugar residue affinities at subsites ?3 and ?2 but also the rate constant for bond cleavage. His101 appears to play multiple roles in the substrate binding and the catalytic reaction.     

Histidine protonation states are key in the LigI catalytic reaction mechanism

Lignin is one of the world's most abundant organic polymers, and 2-pyrone-4,6-dicarboxylate lactonase (LigI) catalyzes the hydrolysis of 2-pyrone-4,6-dicarboxylate (PDC) in the degradation of lignin. The pH has profound effects on enzyme catalysis and therefore we studied this in the context of LigI. We found that changes of the pH mostly affects surface residues, while the residues at the active site are more subject to changes of the surrounding microenvironment. In accordance with this, a high pH facilitates the deprotonation of the substrate. Detailed free energy calculations by the empirical valence bond (EVB) approach revealed that the overall hydrolysis reaction is more likely when the three active site histidines (His31, His33 and His180) are protonated at the ? site, however, protonation at the δ site may be favored during specific steps of the reaction. Our studies have uncovered the determinant role of the protonation state of the active site residues His31, His33 and His180 in the hydrolysis of PDC.     

N-NMR characterization of His residues in and around the active site of FeSOD

We have exploited ¹⁵N-NMR to observe histidine (His) side chains in and around the active site of Fe-containing superoxide dismutase (FeSOD). In the oxidized state, we observe all the non-ligand His side chains and in the reduced state we can account for all the signals in the imidazole spectral region in terms of the non-ligand His′, paramagnetically displaced signals from two backbone amides, and the side chain of glutamine 69 (Gln69). We also observe signals from the His′ that ligate Fe^II. These confirm that neither the Q69H nor the Q69E mutation strongly affects the Fe^II electronic structure, despite the 250 mV and > 660 mV increases in E_m they produce, respectively. In the Q69H mutant, we observe two new signals attributable to the His introduced into the active site in place of Gln69. One corresponds to a protonated N and the other is strongly paramagnetically shifted, to 500 ppm. The strong paramagnetic effects support the existence of an H-bond between His69 and the solvent molecule coordinated to Fe^II, as proposed based on crystallography. Based on previous information that His69 is neutral, we infer that the shifted N is not protonated. Therefore, we propose that this N represents a site of H-bond acceptance from coordinated solvent, representing a reversal of the polarity of this H-bond from that in WT (wild-type) FeSOD protein. We also present evidence that substrate analogs bind to Fe^IISOD outside the Fe^II coordination sphere, affecting Gln69 but without direct involvement of His30.     

A quantitative study of organic phosphate-amine interaction by 31P nuclear magnetic resonance

Interactions between the phosphate group of 4-deoxypyridoxine 5′-phosphate and different protonated amines were quantitatively measured by means of {³¹P}-¹H nuclear magnetic double resonance technique combined with pD titration. An interaction of the phosphate group with added amine resulted in a measurable difference in the ³¹P chemical shift of these phosphate-containing samples with and without amine [Δδ(³¹P)]. Basic amino acids and biogenic amines had significant measurable Δδ(³¹P) values. No interactions were observed for acidic or neutral α, β and γ-amino acids.     

Ribose 5-phosphate isomerase type B from Trypanosoma cruzi: kinetic properties and site-directed mutagenesis reveal information about the reaction mechanism

Trypanosoma cruzi, the human parasite that causes Chagas disease, contains a functional pentose phosphate pathway, probably essential for protection against oxidative stress and also for R5P (ribose 5-phosphate) production for nucleotide synthesis. The haploid genome of the CL Brener clone of the parasite contains one gene coding for a Type B Rpi (ribose 5-phosphate isomerase), but genes encoding Type A Rpis, most frequent in eukaryotes, seem to be absent. The RpiB enzyme was expressed in Escherichia coli as a poly-His tagged active dimeric protein, which catalyses the reversible isomerization of R5P to Ru5P (ribulose 5-phosphate) with Km values of 4 mM (R5P) and 1.4 mM (Ru5P). 4-phospho-D-erythronohydroxamic acid, an analogue to the reaction intermediate when the Rpi acts via a mechanism involving the formation of a 1,2-cis-enediol, inhibited the enzyme competitively, with an IC50 value of 0.7 mM and a Ki of 1.2 mM. Site-directed mutagenesis allowed the demonstration of a role for His102, but not for His138, in the opening of the ribose furanosic ring. A major role in catalysis was confirmed for Cys69, since the C69A mutant was inactive in both forward and reverse directions of the reaction. The present paper contributes to the know-ledge of the mechanism of the Rpi reaction; in addition, the absence of RpiBs in the genomes of higher animals makes this enzyme a possible target for chemotherapy of Chagas disease.     

Role of catalytic residues in the formation of a tetrahedral adduct in the acylation reaction of bovine β-trypsin: A molecular orbital study

In the acylation reaction of serine proteases the effect of amino acid residues on the geometrical change of the catalytic site from Michaelis to tetrahedral state was studied by using ab initio molecular orbital calculations. Amino acid residues in the catalytic site and the peptide substrate were calculated as a quantum mechanical region, and all the other amino acid residues and the calcium ion were included in the calculation as the electrostatic effects. The effects of Asp102, Asp194, N-terminus and the oxyanion binding site are large. The oxyanion binding site directly stabilizes the tetrahedral substrate. Asp102 stabilizes the enzyme intermediate, interacting with the protonated His57 residue. In order to elucidate the roles of Asp102 and the oxyanion binding site, energy decomposition analyses were done for the intermolecular interactions. The contribution of Asp102 and the oxyanion binding site to the decrease of energy in the geometrical change is due to the electrostatic effect. The energies of the proton shuttle from Ser195 O^γ to the leaving group of the substrate were calculated for amide and ester substrate models.     

Modelling C3 photosynthesis: A sensitivity analysis of the photosynthetic carbon-reduction cycle

A model of the C ₃ photosynthetic system is developed which describes the sensitivity of the steadystate rate of carbon dioxide assimilation to changes in the activity of several enzymes of the system. The model requires measurements of the steady-state rate of carbon dioxide assimilation, the concentrations of several intermediates in the photosynthetic system, and the concentration of the active site of ribulose 1,5-bisphosphate carboxyalse/oxygenase (Rubisco). It is shown that in sunflowers (Helianthus annuus L.) at photon flux densities that are largely saturating for the rate of photosynthesis, the steady-stete rate of carbon dioxide assimilation is most sensitive to Rubisco activity and, to a lesser degree, to the activities of the stromal fructose, 6-bisphosphatase and the enzymes catalysing sucrose synthesis. The activities of sedoheptulose 1,7-bisphosphatase, ribulose 5-phosphate kinase, ATP synthase and the ADP-glucose pyrophosphorylase are calculated to have a negligible effect on the flux under the high-light conditions. The utility of this analysis in developing simpler models of photosynthesis is also discussed.Abbreviations c _i intercellular CO₂ concentration - C _infP ^supJ control coefficient for enzyme P with respect to flux J - DHAP dihydroxyacetonephosphate - E4P erythrose 4-phosphate - F6P fructose 6-phosphate - FBP fructose 1,6-bisphosphate - FBPase fructose 1,6-bisphosphatase - G3P glyceraldehyde 3-phosphate - G1P glucose 1-phosphate - G6P glucose 6-phosphate - Pi inorganic phosphate - PCR photosynthetic carbon reduction - PGA 3-phosphoglyceric acid - PPFD photosynthetically active photon flux density - R _n ^J response coefficient for effector n with respect to flux J - R5P ribose 5-phosphate - Rubisco ribulose 1,5-bisphosphate carboxylase/oxygenase - Ru5P ribulose 5-phosphate - RuBP ribulose 1,5-bisphosphate - S7P sedoheptulose 7-phosphate - SBP sedoheptulose 1,7-bisphosphate - SBPase sedoheptulose 1,7-bisphosphatase - SPS sucrose-phosphate synthase - Xu5P xylulose 5-phosphate - _n ^P elasticity coefficient for effector n with respect to the catalytic velocity of enzyme P This research was funded by an Australian Research Council grant to I.E.W. and was undertaken during a visity by K.A.M. to the James Cook University of North Queensland. The expert help of Glenys Hanley and Mick Kelly is greatly appreciated.     

Potential anti‐bacterial drug target: Structural characterization of 3,4‐dihydroxy‐2‐butanone‐4‐phosphate synthase from Salmonella typhimurium LT2

3,4‐Dihydroxy‐2‐butanone‐4‐phosphate synthase (DHBPS) encoded by ribB gene is one of the first enzymes in riboflavin biosynthesis pathway and catalyzes the conversion of ribulose‐5‐phosphate (Ru5P) to 3,4‐dihydroxy‐2‐butanone‐4‐phosphate and formate. DHBPS is an attractive target for developing anti‐bacterial drugs as this enzyme is essential for pathogens, but absent in humans. The recombinant DHBPS enzyme of Salmonella requires magnesium ion for its activity and catalyzes the formation of 3,4‐dihydroxy‐2‐butanone‐4‐phosphate from Ru5P at a rate of 199 nmol min^?1 mg^?1 with K_m value of 116 μM at 37°C. Further, we have determined the crystal structures of Salmonella DHBPS in complex with sulfate, Ru5P and sulfate‐zinc ion at a resolution of 2.80, 2.52, and 1.86 Å, respectively. Analysis of these crystal structures reveals that the acidic loop (residues 34–39) responsible for the acid‐base catalysis is disordered in the absence of substrate or metal ion at the active site. Upon binding either substrate or sulfate and metal ions, the acidic loop becomes stabilized, adopts a closed conformation and interacts with the substrate. Our structure for the first time reveals that binding of substrate Ru5P alone is sufficient for the stabilization of the acidic active site loop into a closed conformation. In addition, the Glu38 residue from the acidic active site loop undergoes a conformational change upon Ru5P binding, which helps in positioning the second metal ion that stabilizes the Ru5P and the reaction intermediates. This is the first structural report of DHBPS in complex with either substrate or metal ion from any eubacteria. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Properties of a pentulose-5-phosphate phosphoketolase from yeasts grown on xylose

Summary Phosphoketolase activity from nine yeasts grown on xylose occurred with both xylulose 5-phosphate (Xu5P) and ribulose 5-phosphate (Ru5P) as substrates. With extracts from five yeasts (Candida curvata, C. famata, Lipomyces starkeyi, Rhodotorula glutinis and Pachysolen tannophilus) activity was approximately the same with either substrate; with C. boidinii, Pichia media and Yarrowia lipolytica Ru5P was the preferred substrate; and with Rhodosporidium toruloides Xu5P was the better substrate. Partial purification of the phosphoketolase from C. famata was attempted: although activity of phosphoketolase towards Ru5P was decreased it was not eliminated and it is concluded that either (i) the phosphoketolase does have dual substrate specificity, in which case it should be referred to as a pentulose-5-phosphate phosphoketolase (Pu5PPK) or (ii) Ru5P-3-epimerase activity, which can interconvert Xu5P and Ru5P, may be closely associated with phosphoketolase activity. The Pu5PPK has a of 2.4mm for Xu5P, a pH optimum of 7.2–7.4 and a M _r of 5x10⁵ daltons. It is not sensitive to inhibition by citrate or acetyl-CoA at physiological concentrations.     

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.     

Some insights into the binding mechanism of the GABAA receptor: a combined docking and MM-GBSA study

Gamma-aminobutyric type A receptor (GABA_AR) is a member of the Cys-loop family of pentameric ligand gated ion channels (pLGICs). It has been identified as a key target for many clinical drugs. In the present study, we construct the structure of human 2α₁2β₂γ₂ GABA_AR using a homology modeling method. The structures of ten benzodiazepine type drugs and two non-benzodiazepine type drugs were then docked into the potential benzodiazepine binding site on the GABA_AR. By analyzing the docking results, the critical residues His102 (α₁), Phe77 (γ₂) and Phe100 (α₁) were identified in the binding site. To gain insight into the binding affinity, molecular dynamics (MD) simulations were performed for all the receptor–ligand complexes. We also examined single mutant GABA_AR (His102A) in complexes with the three drugs (flurazepam, eszopiclone and zolpidem) to elucidate receptor–ligand interactions. For each receptor–ligand complex (with flurazepam, eszopiclone and zolpidem), we calculated the average distance between the C_α of the mutant residue His102A (α₁) to the center of mass of the ligands. The results reveal that the distance between the C_α of the mutant residue His102A (α₁) to the center of flurazepam is larger than that between His102 (α₁) to flurazepam in the WT type complex. Molecular mechanic-generalized Born surface area (MM-GBSA)-based binding free energy calculations were performed. The binding free energy was decomposed into ligand-residue pairs to create a ligand-residue interaction spectrum. The predicted binding free energies correlated well (R ²?=?0.87) with the experimental binding free energies. Overall, the major interaction comes from a few groups around His102 (α₁), Phe77 (γ₂) and Phe100 (α₁). These groups of interaction consist of at least of 12 residues in total with a binding energy of more than 1 kcal mol^?1. The simulation study disclosed herein provides a meaningful insight into GABA_AR–ligand interactions and helps to arrive at a binding mode hypothesis with implications for drug design.     

Purification and characterization of phosphoribulokinase from wheat leaves

Homogeneous phosphoribulokinase (PRK; ATP: d-ribulose-5-phosphate 1-phosphotransferase, EC 2.7.1.19) was isolated from wheat leaves with a specific activity of 15 kat mg^-1 protein. The purification included ammonium sulfate cuts, isoelectric precipitation, and hydrophobic and affinity chromatography on pentylagarose and Blue Sepharose CL 6B, respectively. Gel filtration of the purified enzyme yielded a 83000 Da protein. Subunits of about 42000 Da were estimated from sodium dodecyl sulfate-polyacrylamide gels. Wheat leaf PRK was stable for at least four weeks when stored at 4°C. Saturation curves for ribulose 5-phosphate (Ru5P) and ATP followed Michaelis-Menten kinetics (K _m values: K _{m Ru5P}=50–80 M; K _{m ATP}=70 M). The saturation curve for MgCl₂ was sigmoidal (half-maximal velocity <0.5 mM). The affinity for Ru5P, ATP and Mg²⁺ was not affected by pH changes comparable to pH shifts in the stroma. In contrast to chloroplast fructose-bisphosphatase (Zimmermann et al. 1976, Eur. J. Biochem. 70, 361–367) the affinity for ligands remained unchanged in the dithiothreitol-activated and in the non-activated state. The activity of PRK was increasingly sensitive to inhibition by 3-phosphoglyceric acid with decreasing pH below pH 8.0.Abbreviations DTT dithiothreitol - EDTA ethylenediamine-tetraacetic acid - PRK phosphoribulokinase - Ru5P ribulose-5-phosphate - SDS-PAGE sodium dodecyl sulfate-polyacryl-amide gel electrophoresis     

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.     

Correlation of low-barrier hydrogen bonding and oxyanion binding in transition state analogue complexes of chymotrypsin

The structures of the hemiketal adducts of Ser 195 in chymotrypsin with N-acetyl-L-leucyl-L-phenylalanyl trifluoromethyl ketone (AcLF-CF3) and N-acetyl-L-phenylalanyl trifluoromethyl ketone (AcF-CF3) were determined to 1.4-1.5 A by X-ray crystallography. The structures confirm those previously reported at 1.8-2.1 A [Brady, K., Wei, A., Ringe, D., and Abeles, R. H. (1990) Biochemistry 29, 7600-7607]. The 2.6 A spacings between Ndelta1 of His 57 and Odelta1 of Asp 102 are confirmed at 1.3 A resolution, consistent with the low-barrier hydrogen bonds (LBHBs) between His 57 and Asp 102 postulated on the basis of spectroscopy and deuterium isotope effects. The X-ray crystal structure of the hemiacetal adduct between Ser 195 of chymotrypsin and N-acetyl-L-leucyl-L-phenylalanal (AcLF-CHO) has also been determined at pH 7.0. The structure is similar to the AcLF-CF3 adduct, except for the presence of two epimeric adducts in the R- and S-configurations at the hemiacetal carbons. In the (R)-hemiacetal, oxygen is hydrogen bonded to His 57, not the oxyanion site. On the basis of the downfield 1H NMR spectrum in solution, His 57 is not protonated at Nepsilon2, and there is no LBHB at pH >7.0. Because addition of AcLF-CHO to chymotrypsin neither releases nor takes up a proton from solution, it is concluded that the hemiacetal oxygen of the chymotrypsin-AcLF-CHO complex is a hydroxyl group and not attracted to the oxyanion site. The protonation states of the hemiacetal and His 57 are explained by the high basicity of the hemiacetal oxygen (pK(a) > 13.5) relative to that of His 57. The 13C NMR signal for the adduct of AcLF-13CHO with chymotrypsin is consistent with a neutral hemiacetal between pH 7 and 13. At pH <7.0, His 57 in the AcLF-CHO-hemiacetal complex of chymotrypsin undergoes protonation at Nepsilon2 of His 57, leading to a transition of the 15.1 ppm downfield signal to 17.8 ppm. The pK(a)s in the active sites of the AcLF-CF3 and AcLF-CHO adducts suggest an energy barrier of 6-7 kcal x mol(-1) against ionizations that change the electrostatic charge at the active site. However, ionizations of neutral His 57 in the AcLF-CHO-chymotrypsin adduct, or in free chymotrypsin, proceed with no apparent barrier. Protonation of His 57 is accompanied by LBHB formation, suggesting that stabilization by the LBHB overcomes the barrier to ionization. On the basis of the hydration constant for AcLF-13CHO and its inhibition constant, its K(d) is 16 microM, 8000-fold larger than the comparable value for AcLF-CF3.