Checking the pH-induced conformational transition of prion protein by molecular dynamics simulations: effect of protonation of histidine residues

The role of acidic pH in the conversion of human prion protein to the pathogenic isoform is investigated by means of molecular dynamics simulations, focusing the attention on the effect of protonation of histidine residues on the conformational behavior of human PrPC globular domain. Our simulations reveal a significant loss of alpha-helix content under mildly acidic conditions, due to destructuration of the C-terminal part of HB (thus suggesting a possible involvement of HB into the conformational transition leading to the pathogenic isoform) and a transient lengthening of the native beta-sheet. Protonation of His-187 and His-155 seems to be crucial for the onset of the conformational rearrangement. This finding can be related to the existence of a pathogenic mutation, H187R, which is associated with GSS syndrome. Finally, the relevance of our results for the location of a Cu2+-binding pocket in the C-terminal part of the prion is discussed.     

Influence of pH on NMR structure and stability of the human prion protein globular domain

The NMR structure of the globular domain of the human prion protein (hPrP) with residues 121-230 at pH 7.0 shows the same global fold as the previously published structure determined at pH 4.5. It contains three alpha-helices, comprising residues 144-156, 174-194, and 200-228, and a short anti-parallel beta-sheet, comprising residues 128-131 and 161-164. There are slight, strictly localized, conformational changes at neutral pH when compared with acidic solution conditions: helix alpha1 is elongated at the C-terminal end with residues 153-156 forming a 310-helix, and the population of helical structure in the C-terminal two turns of helix alpha 2 is increased. The protonation of His155 and His187 presumably contributes to these structural changes. Thermal unfolding monitored by far UV CD indicates that hPrP-(121-230) is significantly more stable at neutral pH. Measurements of amide proton protection factors map local differences in protein stability within residues 154-157 at the C-terminal end of helix alpha 1 and residues 161-164 of beta-strand 2. These two segments appear to form a separate domain that at acidic pH has a larger tendency to unfold than the overall protein structure. This domain could provide a "starting point" for pH-induced unfolding and thus may be implicated in endosomic PrPC to PrPSc conformational transition resulting in transmissible spongiform encephalopathies.     

The pH-dependent structural variation of complementarity-determining region H3 in the crystal structures of the Fv fragment from an anti-dansyl monoclonal antibody.

The Fv fragment from an anti-dansyl antibody was optimally crystallized into two crystal forms having slightly different lattice dimensions at pH 5.25 and 6.75. The two crystal structures were determined and refined at high resolution at 112 K (at 1.45 A for the crystal at pH 5.25 and at 1.55 A for that at pH 6.75). In the two crystal structures, marked differences were identified in the first half of CDRH3 s having an amino acid sequence of Ile95H-Tyr96H-Tyr97H-His98H-Tyr99H-Pro1 00H-Trp100aH-Phe100bH-Ala101H- Tyr102H. NMR pH titration experiments revealed the p Kavalues of four histidine residues (His27dL, His93L, His55H and His98H) exposed to solvent. Only His98H (p Ka=6.3) completely changed its protonation state between the two crystallization conditions. In addition, the environmental structures including hydration water molecules around the four histidine residues were carefully compared. While the hydration structures around His27dL, His93L and His55H were almost invariant between the two crystal structures, those around His98Hs showed great difference in spite of the small conformational difference of His98H between the two crystal structures. These spectroscopic and crystallographic findings suggested that the change in the protonation state in His98H was responsible for the structural differences between pH 5.25 and 6.75. In addition, the most plausible binding site of the dansyl group was mapped into the present structural models with our previous NMR experimental results. The complementarity-determining regions H1, H3 and the N-terminal region in the VH domain formed the site. The side-chain of Tyr96H occupied the site and interacted with Phe27H of H1, giving a clue for the binding mode of the dansyl group in the site.     

Concerted protonation of key histidines triggers membrane interaction of the diphtheria toxin T domain

The translocation domain (T domain) of the diphtheria toxin contributes to the transfer of the catalytic domain from the cell endosome to the cytosol, where it blocks protein synthesis. Translocation is initiated when endosome acidification induces the interaction of the T domain with the membrane of the compartment. We found that the protonation of histidine side chains triggers the conformational changes required for membrane interaction. All histidines are involved in a concerted manner, but none is indispensable. However, the preponderance of each histidine varies according to the transition observed. The pair His(223)-His(257) and His(251) are the most sensitive triggers for the formation of the molten globule state in solution, whereas His(322)-His(323) and His(251) are the most sensitive triggers for membrane binding. Interestingly, the histidines are located at key positions throughout the structure of the protein, in hinges and at the interface between each of the three layers of helices forming the domain. Their protonation induces local destabilizations, disrupting the tertiary structure and favoring membrane interaction. We propose that the selection of histidine residues as triggers of membrane interaction enables the T domain to initiate translocation at the rather mild pH found in the endosome, contributing to toxin efficacy.     

NMR dynamic studies suggest that allosteric activation regulates ligand binding in chicken liver bile acid-binding protein

Apo chicken liver bile acid-binding protein has been structurally characterized by NMR. The dynamic behavior of the protein in its apo- and holo-forms, complexed with chenodeoxycholate, has been determined via (15)N relaxation and steady state heteronuclear (15)N((1)H) nuclear Overhauser effect measurements. The dynamic parameters were obtained at two pH values (5.6 and 7.0) for the apoprotein and at pH 7.0 for the holoprotein, using the model free approach. Relaxation studies, performed at three different magnetic fields, revealed a substantial conformational flexibility on the microsecond to millisecond time scales, mainly localized in the C-terminal face of the beta-barrel. The observed dynamics are primarily caused by the protonation/deprotonation of a buried histidine residue, His(98), located on this flexible face. A network of polar buried side chains, defining a spine going from the E to J strand, is likely to provide the long range connectivity needed to communicate motion from His(98) to the EF loop region. NMR data are accompanied by molecular dynamics simulations, suggesting that His(98) protonation equilibrium is the triggering event for the modulation of a functionally important motion, i.e. the opening/closing at the protein open end, whereas ligand binding stabilizes one of the preexisting conformations (the open form). The results presented here, complemented with an analysis of proteins belonging to the intracellular lipid-binding protein family, are consistent with a model of allosteric activation governing the binding mechanism. The functional role of this mechanism is thoroughly discussed within the framework of the mechanism for the enterohepatic circulation of bile acids.     

High resolution NMR studies of histidine-substituted and histidine-perturbed hemoglobin variants. Histidine assignments, electrostatic interactions at the protein surface, and implications for hemoglobin S polymerization

The aromatic region of the proton NMR spectrum of human adult hemoglobin (HbA) contains resonances from at least 11 titratable histidine residues. Assignments for five beta chain histidines have previously been proposed. In order to further characterize the aromatic spectra of HbA we studied 11 histidine-substituted and -perturbed hemoglobin variants in oxy and deoxy states and at different pH values by 400 MHz NMR spectroscopy. We propose assignments for the resonances corresponding to the C2 protons of His alpha 20, His alpha 72, His alpha 112, and His beta 77 in oxy and deoxy spectra and of His beta 97 and His beta 117 in deoxy spectra. Our assignments for His beta 2 and His beta 117 in the oxy state agree with those previously reported for the CO form, but in the deoxy state our spectra suggest a different assignment. Studies with Hb variants in which a histidine is perturbed by a neighboring substitution suggest additional assignments for His alpha 50 and His alpha 89 and demonstrate a strong dependence of the imidazole ring pK on hydrogen bond interactions and on the net charge of neighboring residues. Some of the newly proposed assignments of histidine resonances are used to discuss specific intermolecular interactions implicating His alpha 20, His beta 77, and His beta 117 in deoxy HbS polymers.     

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.     

Effect of 2-fluorohistidine labeling of the anthrax protective antigen on stability, pore formation, and translocation

The action of anthrax toxin relies in part upon the ability of the protective antigen (PA) moiety to form a heptameric pore in the endosomal membrane, providing a portal for entry of the enzymic moieties of the toxin into the cytosol. Pore formation is dependent on a conformational change in the heptameric prepore that occurs in the neutral to mildly acidic pH range, and it has been hypothesized that protonation of one or more histidine residues triggers this transition. To test this hypothesis, we used biosynthetic methods to incorporate the unnatural amino acid analogue 2-fluorohistidine (2-FHis) into PA. 2-FHis is isosteric with histidine but resists protonation at physiological pH values due to a dramatically reduced side-chain pKa ( approximately 1). We found that 2-FHis-labeled PA was biologically inactive, as judged by its inability to deliver a model intracellular effector, LFN-DTA, to the cytosol of CHO-K1 cells. However, whereas 2-FHis blocked a conformational transition in the full-length PA83 protein in the pH 5-6 range, the pH dependence of prepore-to-pore conversion of (PA63)7 was unchanged from the wild-type protein, implying that this conversion is not dependent on His protonation. Consistent with this result, the labeled, trypsin-activated PA was able to permeabilize liposomes to K+ and retained pore-forming activity in planar phospholipid bilayers. The pores in planar bilayers were incapable, however, of translocating a model ligand in response to a transmembrane pH gradient or elevated voltage. The results indicate that protonation of residues other than His, presumably Glu and/or Asp side chains, triggers pore formation in vitro, but His residues are nonetheless important for PA functioning in vivo.     

Histidine pK(a) shifts and changes of tautomeric states induced by the binding of gallium-protoporphyrin IX in the hemophore HasA(SM)

The HasA(SM) hemophore, secreted by Serratia marcescens, binds free or hemoprotein bound heme with high affinity and delivers it to a specific outer membrane receptor, HasR. In HasA(SM), heme is held by two loops and coordinated to iron by two residues, His 32 and Tyr 75. A third residue His 83 was shown recently to play a crucial role in heme ligation. To address the mechanistic issues of the heme capture and release processes, the histidine protonation states were studied in both apo- and holo-forms of HasA(SM) in solution. Holo-HasA(SM) was formed with gallium-protoporphyrin IX (GaPPIX), giving rise to a diamagnetic protein. By use of heteronuclear correlation NMR spectroscopy, the imidazole side-chain (15)N and (1)H resonances of the six HasA(SM) histidines were assigned and their pKa values and predominant tautomeric states according to pH were determined. We show that protonation states of the heme pocket histidines can modulate the nucleophilic character of the two axial ligands and, consequently, control the heme binding. In particular, the essential role of the His 83 is emphasized according to its direct interaction with Tyr 75.     

Role of individual histidines in the pH-dependent global stability of human chloride intracellular channel 1

Chloride intracellular channel proteins exist in both a soluble cytosolic form and a membrane-bound form. The mechanism of conversion between the two forms is not properly understood, although one of the contributing factors is believed to be the variation in pH between the cytosol (~7.4) and the membrane (~5.5). We systematically mutated each of the three histidine residues in CLIC1 to an alanine at position 74 and a phenylalanine at positions 185 and 207. We examined the effect of the histidine-mediated pH dependence on the structure and global stability of CLIC1. None of the mutations were found to alter the global structure of the protein. However, the stability of H74A-CLIC1 and H185F-CLIC1, as calculated from the equilibrium unfolding data, is no longer dependent on pH because similar trends are observed at pH 7.0 and 5.5. The crystal structures show that the mutations result in changes in the local hydrogen bond coordination. Because the mutant total free energy change upon unfolding is not different from that of the wild type at pH 7.0, despite the presence of intermediates that are not seen in the wild type, we propose that it may be the stability of the intermediate state rather than the native state that is dependent on pH. On the basis of the lower stability of the intermediate in the H74A and H185F mutants compared to that of the wild type, we conclude that both His74 and His185 are involved in triggering the pH changes to the conformational stability of wild-type CLIC1 via their protonation, which stabilizes the intermediate state.     

The charge-relay system of serine proteinases: proton magnetic resonance titration studies of the four histidines of porcine trypsin.

Peaks corresponding to the C₍₂₎-protons of all four histidine residues of porcine β-trypsin were resolved in 250 MHz nuclear magnetic resonance spectra after deuteration of the slowly exchangeable N-H groups (whose resonances obscure the histidine peaks) by reversible unfolding of the protein in D₂O. One of the four peaks was assigned to the charge-relay histidine in the active site of trypsin (His(57) in the bovine chymotrypsinogen numbering system). Whereas the three other histidine C₍₂₎-peaks exhibited normal titration curves with single pK′ values of 7.20, 6.71 and 6.67, the peak assigned to His(57) had an abnormal titration curve showing two protonation steps in the pH range from 1 to 9. The first protonation with a pH′_mid of 5.0 is rapid on the nuclear magnetic resonance time-scale; the second with a pH′_mid of 4.5 is slow and apparently involves conformational transitions between two states having lifetimes of approximately 18 ms.In the complex between porcine β-trypsin and bovine pancreatic trypsin inhibitor (Kunitz) His(57) was found to be insensitive to pH over the range from 4 to 9 and its chemical shift resembles that of His(57) in the singly protonated charge relay of free trypsin. This result provides direct evidence that the trypsin charge relay acts as a proton acceptor in the initial catalytic step which leads to the formation of a tetrahedral complex. In the presence of equimolar bovine pancreatic trypsin inhibitor (Kunitz) the pH'_mid of the conformational transition that affects the charge-relay histidine is lowered from 4.5 to approximately 3.5.     

pK(a) calculations of calbindin D(9k): effects of Ca(2+) binding, protein dielectric constant, and ionic strength

Calbindin is a small (75 residues) helix-loop-helix ("EF-hand") calcium-binding protein belonging to the calmodulin superfamily. It binds two Ca(2+) ions. Continuum electrostatics in combination with the boundary element method was employed for the calculation of the acid-dissociation constants K(a) (pK(a) = -log K(a)) values of all titratable residues in the protein. The objectives were to determine quantitatively the effects of divalent ion binding and small ion-induced structural changes on predicted pK(a)'s. Computations were carried out for the apo and holo form of calbindin, for which both X-ray and NMR structures were available. Comparison was made with several sets of experimental pK(a) values determined by NMR spectroscopy. Different choices of the dielectric constant (ranging from 4 to 78.5) for calbindin and variations in ionic strength (from 0 to 0.3 M) were investigated in a systematic fashion. Removal of the two bound Ca(2+) ions increases the pK(a) values of all residues if no conformational changes were allowed. If conformational differences between the apo and holo were accounted for, shifts in either direction were observed. Titrating groups that are directly involved in Ca(2+) binding (Asp and Glu) required a dielectric constant of 78.5 for the holo structure to obtain a reasonable estimate of their pK(a)'s. For the apo structure, passable values for the pK(a)'s of these ligating groups could be determined if the structure was allowed to relax upon ion removal.     

NMR detection of multiple transitions to low-populated states in azurin

Transitions to conformational states with very low populations were detected for the reduced blue copper protein azurin from Pseudomonas aeruginosa by applying constant relaxation time CPMG measurements to the backbone (15)N nuclei at three magnetic fields (11.7, 14.1, and 18.8 T) and three temperatures (25.7, 35.4, and 44.8 degrees C). Two exchange processes with different rate constants could be discriminated despite populations of the excited states below 1% and spatial neighborhood of the two processes. The group of (15)N nuclei involved in the faster process exhibits at 44.8 degrees C a forward rate constant of 11.7+/-2.4 s(-1) and a population of the exited state of 0.39+/-0.07%. They surround the aromatic ring of histidine 35 whose protonation state is coupled to the flipping of a neighboring peptide plane. For the slower process, the forward rate constant and population of the exited state at 44.8 degrees C are 4.1+/-0.1 s(-1) and 0.45+/-0.02%, respectively. The residues involved cluster nearby the copper ion, which is separated from the protonation site of histidine 35 by about 8 A, indicating conformational rearrangements involving the copper coordinating loops. The dependence of the equilibrium constant on the temperature is consistent with an enthalpy-dominated transition around the copper, but an entropy-controlled transition near histidine 35. The detection by nuclear magnetic resonance of millisecond to second conformational transitions near the copper ion suggests a low energy-cost rearrangement of the copper-binding site that may be necessary for efficient electron transfer.     

Histidine pK(a) values for the N-lobe of human transferrin: effect of substitution of binding site Asp by Ser (D63S)

The pK(a) values have been determined for eight of the nine histidine residues and the amino terminus of the N-lobe of human apo-transferrin (hTF/2N), and for seven of the nine histidine residues and the amino terminus of the protein Asp63Ser hTF/2N containing a mutation of the Fe(3+)-ligand Asp63 to Ser63. Calculations suggested that substitution of aspartate by serine would result in decreases of the pK(a) values of most of the histidine residues in the protein. This was found to be the case experimentally, and allowed assignment of the varepsilonCH resonance of His249. For the wild-type protein, the His residue with a pK(a) of 7.40 was assigned as His249, whereas for the mutant, no observable His residue had a pK(a) value higher than 6.9. The protonated form of His249 appears to be stabilised by interactions with Asp63, and the high pK(a) value may be critical for ensuring the release of iron at endosomal pH (5.5). The mutation lowered the apparent binding constant of hTF/2N for the synergistic anion oxalate from log K 4.0 to log K 3.3. (1)H NMR spectral changes induced by Ga(3+) binding to the mutant are compared to those observed for the wild-type protein.     

Heteronuclear NMR and crystallographic studies of wild-type and H187Q Escherichia coli uracil DNA glycosylase: electrophilic catalysis of uracil expulsion by a neutral histidine 187.

The nature of the putative general acid His187 in the reaction catalyzed by Escherichia coli uracil DNA glycosylase (UDG) was investigated using X-ray crystallography and NMR spectroscopy. The crystal structures of H187Q UDG, and its complex with uracil, have been solved at 1.40 and 1.60 A resolution, respectively. The structures are essentially identical to those of the wild-type enzyme, except that the side chain of Gln187 is turned away from the uracil base and cannot interact with uracil O2. This result provides a structural basis for the similar kinetic properties of the H187Q and H187A enzymes. The ionization state of His187 was directly addressed with (1)H-(15)N NMR experiments optimized for histidine ring spin systems, which established that His187 is neutral in the catalytically active state of the enzyme (pK(a) <5.5). These NMR experiments also show that His187 is held in the N(epsilon)()2-H tautomeric form, consistent with the crystallographic observation of a 2.9 A hydrogen bond from the backbone nitrogen of Ser189 to the ring N(delta)()1 of His187. The energetic cost of breaking this hydrogen bond may contribute significantly to the low pK(a) of His187. Thus, the traditional view that a cationic His187 donates a proton to uracil O2 is incorrect. Rather, we propose a concerted mechanism involving general base catalysis by Asp64 and electrophilic stabilization of the developing enolate on uracil O2 by a neutral His187.     

The pH-induced release of iron from transferrin investigated with a continuum electrostatic model.

A reduction in pH induces the release of iron from transferrin in a process that involves a conformational change in the protein from a closed to an open form. Experimental evidence suggests that there must be changes in the protonation states of certain, as yet not clearly identified, residues in the protein accompanying this conformational change. Such changes in protonation states of residues and the consequent changes in electrostatic interactions are assumed to play a large part in the mechanism of release of iron from transferrin. Using the x-ray crystal structures of human ferri- and apo-lactoferrin, we calculated the pKa values of the titratable residues in both the closed (iron-loaded) and open (iron-free) conformations with a continuum electrostatic model. With the knowledge of a residue's pKa value, its most probable protonation state at any specified pH may be determined. The preliminary results presented here are in good agreement with the experimental observation that the binding of ferric iron and the synergistic anion bicarbonate/carbonate results in the release of approximately three H+ ions. It is suggested that the release of these three H+ ions may be accounted for, in most part, by the deprotonation of the bicarbonate and residues Tyr-92, Lys-243, Lys-282, and Lys-285 together with the protonation of residues Asp-217 and Lys-277.     

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.     

Conformation and stability properties of B17: I. Analytical investigations using circular dichroism

Structural characterization of B17, the 17?% N-terminal domain of apo B, was carried out using circular dichroic (CD) spectroscopy, where secondary and tertiary structures were studied as a function of temperature and pH. Mild acidic conditions that correlate with histidine protonation invoked a change in the α-helix and random coil contents of the protein, with no apparent change in the β-sheet structural content. Specific changes in the structure of the protein that occur in response to temperature were also investigated to understand the stability and conformational changes of B17. Far- and near-UV CDs were used to probe the thermal changes in the protein. The protonation of some histidine residues was attributed to underlie the increase in the helical content of the protein.     

High-resolution solution structure of Bacillus subtilis IIAglc

The high-resolution solution structure of the phosphocarrier protein IIA^glc from Bacillus subtilis is determined using 3D and 4D heteronuclear NMR methods. B. subtilis IIA^glc contains 162 amino acid residues and is one of the larger proteins for which high-resolution solution structure has been determined by NMR methods. The structures have been calculated from a total of 2,232 conformational constraints. Comparison with the X-ray crystal structure indicates that the overall fold is the same in solution and in crystalline environments, although some local structural differences are observed. These occur largely in turns and loops, and mostly correspond to regions with high-temperature factors in the crystal structure. The N-terminus of IIA^glc is disordered in solution. The active site is located in a concave region of the protein surface. The histidine, which accepts the phosphoryl group (His 83), interacts with a neighboring histidine (His 68) and is surrounded by hydrophobic residues. Proteins 31:258–270, 1998. © 1998 Wiley-Liss, Inc.     

Solution structure of Methylophilus methylotrophus cytochrome c": insights into the structural basis of haem-ligand detachment

Cytochrome c" from Methylophilus methylotrophus is a monohaem protein with 124 amino acid residues. The iron has two histidine ligands in the oxidised form, one of which detaches and picks up a proton when the protein is reduced. Thus, both forms are paramagnetic. The structure of the oxidised form in solution, determined from NMR data is presented. The family of structures has an average backbone rmsd value of 0.53 A, and a heavy atom rmsd value of 0.95 A, within a target function range of 32 %. This structure is related to class I cytochromes with an additional helix at the N terminus. The haem-binding site occurs in a domain essentially lacking secondary structure motifs and the axial histidinyl residues were found in an unusual near perpendicular orientation. Moreover, a disulfide bridge is present, an uncommon structural feature among c-type cytochromes. The disulfide bridge, linking cysteine residues 96 and 104, forms a loop that confers rigidity and is essential to the detachment of the axial histidine (His95) as demonstrated by chemical disruption of the S-S bond. A route for protonation of the distal histidine involving haem propionate 17 is proposed and discussed in the light of available models for complex membrane proton pumps.