Assessing the acid-base and conformational properties of histidine residues in human prion protein (125-228) by means of pK(a) calculations and molecular dynamics simulations

A thorough study of the acid-base behavior of the four histidines and the other titratable residues of the structured domain of human prion protein (125-228) is presented. By using multi-tautomer electrostatic calculations, average titration curves have been built for all titratable residues, using the whole bundles of NMR structures determined at pH 4.5 and 7.0. According to our results, (1) only histidine residues are likely to be involved in the first steps of the pH-driven conformational transition of prion protein; (2) the pK(a)'s of His140 and His177 are approximately 7.0, whereas those of His155 and His187 are < 5.5. 10-ns long molecular dynamics simulations have been performed on five different models, corresponding to the most significant combinations of histidine protonation states. A critical comparison between the available NMR structures and our computational results (1) confirms that His155 and His187 are the residues whose protonation is involved in the conformational rearrangement of huPrP in mildly acidic condition, and (2) shows how their protonation leads to the destructuration of the C-terminal part of HB and to the loss of the last turn of HA that represent the crucial microscopic steps of the rearrangement.     

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.     

Acid-induced equilibrium folding intermediate of human platelet profilin

The acid-induced unfolding of human platelet profilin (HPP) can be minimally modeled as a three-state process. Equilibrium unfolding studies have been performed on human platelet profilin1 (HPP) and monitored by far-UV circular dichroism, tryptophan fluorescence, ANS binding, and NMR spectroscopy. Far-UV CD measurements obtained by acid titration demonstrate that HPP unfolds via a three-state mechanism (N --> I --> U), with a highly populated intermediate between pH 4 and 5. Approximately 80% of native helical secondary structural content remains at pH 4, as indicated by monitoring the CD signal at 222 nm. The stability (DeltaGH2O) of the native conformation at pH 7.0 (obtained by monitoring the change in tryptophan signal as a function of urea concentration) is 5.56 +/- 0.51 kcal mol-1; however, the DeltaGH2O for the intermediate species at pH 4 is 2.01 +/- 0.47 kcal mol-1. The calculated m-values for the pH 7.0 and pH 4.0 species were 1.64 +/- 0.15 and 1.34 +/- 0.17 kcal mol-1 M-1, respectively, which is an indication that the native and intermediate species are similarly compact. Additionally, translational diffusion measurements obtained by NMR spectroscopy and ANS binding studies are consistent with a globular and compact conformation at both pH 7.0 and 4.0. The pKa values for the two histidine (His) residues located on helix 4 of HPP were determined to be 5.6 and 5.7 pH units. These pKa values coincide with the midpoint of the far-UV CD acid titration curve and suggest that the protonation of one or both His residues may play a role in the formation of the unfolding intermediate. Stable intermediate species populate the 2D 1H-15N HSQC NMR spectra between pH 4 and 5. A number of backbone and side-chain resonances show significant perturbations relative to the native spectrum; however, considerable nativelike tertiary contacts remain. Interestingly, the residues on HPP that are significantly altered at low pH coincide with segments of the G-actin binding surface and poly-l-proline binding interface. The earlier reports that a decrease in pH below 6.0 induces structural alterations in profilin, favoring dissociation of the profilin-actin complex, corresponds with the structural alterations observed in the partially unfolded species. Our findings suggest that a novel mechanism for pH induced disruption of the profilin-G-actin complex involve a nativelike unfolding intermediate of profilin.     

Fast events in protein folding: structural volume changes accompanying the early events in the N-->I transition of apomyoglobin induced by ultrafast pH jump

Ultrafast, laser-induced pH jump with time-resolved photoacoustic detection has been used to investigate the early protonation steps leading to the formation of the compact acid intermediate (I) of apomyoglobin (ApoMb). When ApoMb is in its native state (N) at pH 7.0, rapid acidification induced by a laser pulse leads to two parallel protonation processes. One reaction can be attributed to the binding of protons to the imidazole rings of His24 and His119. Reaction with imidazole leads to an unusually large contraction of -82 +/- 3 ml/mol, an enthalpy change of 8 +/- 1 kcal/mol, and an apparent bimolecular rate constant of (0.77 +/- 0.03) x 10(10) M(-1) s(-1). Our experiments evidence a rate-limiting step for this process at high ApoMb concentrations, characterized by a value of (0. 60 +/- 0.07) x 10(6) s(-1). The second protonation reaction at pH 7. 0 can be attributed to neutralization of carboxylate groups and is accompanied by an apparent expansion of 3.4 +/- 0.2 ml/mol, occurring with an apparent bimolecular rate constant of (1.25 +/- 0.02) x 10(11) M(-1) s(-1), and a reaction enthalpy of about 2 kcal/mol. The activation energy for the processes associated with the protonation of His24 and His119 is 16.2 +/- 0.9 kcal/mol, whereas that for the neutralization of carboxylates is 9.2 +/- 0.9 kcal/mol. At pH 4.5 ApoMb is in a partially unfolded state (I) and rapid acidification experiments evidence only the process assigned to carboxylate protonation. The unusually large contraction and the high energetic barrier observed at pH 7.0 for the protonation of the His residues suggests that the formation of the compact acid intermediate involves a rate-limiting step after protonation.     

Ionization of His 55 at the dimer interface of dynein light-chain LC8 is coupled to dimer dissociation

LC8 is a highly conserved light-chain subunit of cytoplasmic dynein that interacts with a wide variety of cellular proteins and is presumed to play a fundamental role in dynein assembly and cargo recruitment and in the assembly of protein complexes unrelated to dynein. LC8 is a dimer at physiological pH but dissociates to a folded monomer at pH < 4.8. We have suggested that acid-induced dimer dissociation is due to protonation of His 55, which is stacked against His 55' and completely buried in the dimer interface. In this work, we show that the pH-induced dissociation is reversible and indeed governed by the ionization state of His 55. Mutagenesis of His 55 to Lys results in a monomer in the pH range of 3-8, while the mutation to Ala results in a dimer in the same pH range. Mutations that disrupt intermolecular hydrogen bonds between Tyr 65 and Lys 44' and His 55 and Thr 67' do not change the association state of the dimer. Titration curves for His 55 and the two other histidines, His 72 and 68, were determined by (13)C-(1)H NMR for H55K and for WT-LC8 in the monomeric and dimeric states. The pK(a) values of His 72 and His 68 are 6 in the WT dimer and 6.2-6.5 in monomeric H55K, while the pK(a) of His 55 is about 4.5 in the WT dimer. These results indicate that deprotonation of His 55 is linked to dimer formation and that mutation of His 55 to a small neutral residue or to a positively charged residue uncouples the protonation and dissociation processes.     

Two independent histidines, one in human prolactin and one in its receptor, are critical for pH-dependent receptor recognition and activation

Human prolactin (hPRL), a member of the family of hematopoietic cytokines, functions as both an endocrine hormone and autocrine/paracrine growth factor. We have previously demonstrated that recognition of the hPRL·receptor depends strongly on solution acidity over the physiologic range from pH 6 to pH 8. The hPRL·receptor binding interface contains four histidines whose protonation is hypothesized to regulate pH-dependent receptor recognition. Here, we systematically dissect its molecular origin by characterizing the consequences of His to Ala mutations on pH-dependent receptor binding kinetics, site-specific histidine protonation, and high resolution structures of the intermolecular interface. Thermodynamic modeling of the pH dependence to receptor binding affinity reveals large changes in site-specific protonation constants for a majority of interface histidines upon complexation. Removal of individual His imidazoles reduces these perturbations in protonation constants, which is most likely explained by the introduction of solvent-filled, buried cavities in the crystallographic structures without inducing significant conformational rearrangements.     

Indirect DNA Readout on the Protein Side: Coupling between Histidine Protonation, Global Structural Cooperativity, Dynamics, and DNA Binding of the Human Papillomavirus Type 16 E2C Domain

DNA sequence recognition by the homodimeric C-terminal domain of the human papillomavirus type 16 E2 protein (E2C) is known to involve both direct readout and DNA-dependent indirect readout mechanisms, while protein-dependent indirect readout has been deduced but not directly observed. We have investigated coupling between specific DNA binding and the dynamics of the unusual E2C fold, using pH as an external variable. Nuclear magnetic resonance and isothermal titration calorimetry show that pH titration of His318 in the complex interface and His288 in the core of the domain is coupled to both binding and the dynamics of the β-barrel core of E2C, with a tradeoff between dimer stability and function. Specific DNA binding is, in turn, coupled to the slow dynamics and amide hydrogen exchange in the entire β-barrel, reaching residues far apart from the DNA recognition elements but not affecting the two helices of each monomer. The changes are largest in the dimerization interface, suggesting that the E2C β-barrel acts as a hinge that regulates the relative position of the DNA recognition helices. In conclusion, the cooperative dynamics of the human papillomavirus type 16 E2C β-barrel is coupled to sequence recognition in a protein-dependent indirect readout mechanism. The patterns of residue substitution in genital papillomaviruses support the importance of the protonation states of His288 and His318 and suggest that protein-dependent indirect readout and histidine pH titration may regulate DNA binding in the cell.     

Probing the pH-dependent structural features of alpha-KTx12.1, a potassium channel blocker from the scorpion Tityus serrulatus

Potassium channels are widespread in living cells and are involved in many diseases. The scorpion toxin alpha-KTx(12.1) interacts with various K(+) channels, suggesting its capacity to match diverse channel pores. It is recognized that tissue injuries may affect the pH at toxins site of action, thereby modulating both protein conformation and activity. To better understand its molecular mechanism of action, we studied alpha-KTx(12.1) using pH as a tool to explore its plasticity and NMR in combination with MD calculations to detect it. The toxin solution structure consists of an alpha-helix and a triple-stranded beta-sheet stabilized by four disulfide bridges. The NMR results show, in addition, that His28 possesses an unusually low pK(a) of 5.2. The best set of protein conformers is obtained at pH 4.5, while at pH 7.0, the reduced number of NOEs resulting from a faster hydrogen exchange does not allow to reach a good structural convergence. Nonetheless, MD calculations show that the toxin structure does not vary significantly in that pH range, while conformational changes and modifications of the surface charge distribution occur when His28 is fully protonated. Moreover, essential dynamics analysis reveals variations in the toxin's coherent motions. In conclusion, His28, with its low pK(a) value, provides alpha-KTx(12.1) with the ability to preserve its active conformation over a wide pH interval, thus expanding the range of cellular conditions where the toxin can fully exhibit its activity. Overall, the results further underline the role of histidine as a natural controller of proteins' functionality.     

Constant pH molecular dynamics (CpHMD) and molecular docking studies of CquiOBP1 pH-induced ligand releasing mechanism

The odorant binding protein of Culex quinquefasciatus (CquiOBP1), expressed on the insect antenna, is crucial for the investigation of trapping baited with oviposition semi-chemicals and controlling mosquito populations. The acidic titratable residues pKa prediction and the ligand binding poses investigation in two systems (pH 7 and pH 5) are studied by constant pH molecular dynamics (CpHMD) and molecular docking methods. Research results reveal that the change of the protonation states would disrupt some important H-bonds, such as Asp 66-Asp 70, Glu 105-Asn 102, etc. The cleavage of these H-bonds leads to the movement of the relative position of hydrophobic tunnel, N- and C- termini loops and pH-sensing triad (His23-Tyr54-Val125) in acid solution. Ligand MOP has lower affinity and shows different binding poses to protein CquiOBP1 at pH 5. This ligand may be released from another tunnel between helices α3 and α4 in acidic environment. However, it would bind to the protein with high affinity in neutral environment. This work could provide more penetrating understanding of the pH-induced ligand-releasing mechanism.     

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.     

Electrostatic interactions in the acid denaturation of alpha-lactalbumin determined by NMR.

alpha-Lactalbumin (alpha-LA) undergoes a pH-dependent unfolding from the native state to a partially unfolded state (the molten globule state). To understand the role of electrostatic interactions in protein denaturation, NMR and CD pH titration experiments are performed on guinea pig alpha-LA. Variation of pH over the range of 7.0 to 2.0 simultaneously leads to the acid denaturation of the protein and the titration of individual ionizable groups. The pH titrations are interpreted in the context of these coupled events, and indicate that acid denaturation in alpha-LA is a cooperative event that is triggered by the protonation of two ionizable residues. Our NMR results suggest that the critical electrostatic interactions that contribute to the denaturation of alpha-LA are concentrated in the calcium binding region of the protein.     

Increase in the conformational flexibility of beta 2-microglobulin upon copper binding: a possible role for copper in dialysis-related amyloidosis

A key pathological event in dialysis-related amyloidosis is the fibril formation of beta(2)-microglobulin (beta 2-m). Because beta 2-m does not form fibrils in vitro, except under acidic conditions, predisposing factors that may drive fibril formation at physiological pH have been the focus of much attention. One factor that may be implicated is Cu(2+) binding, which destabilizes the native state of beta 2-m and thus stabilizes the amyloid precursor. To address the Cu(2+)-induced destabilization of beta 2-m at the atomic level, we studied changes in the conformational dynamics of beta 2-m upon Cu(2+) binding. Titration of beta 2-m with Cu(2+) monitored by heteronuclear NMR showed that three out of four histidines (His13, His31, and His51) are involved in the binding at pH 7.0. (1)H-(15)N heteronuclear NOE suggested increased backbone dynamics for the residues Val49 to Ser55, implying that the Cu(2+) binding at His51 increased the local dynamics of beta-strand D. Hydrogen/deuterium exchange of amide protons showed increased flexibility of the core residues upon Cu(2+) binding. Taken together, it is likely that Cu(2+) binding increases the pico- to nanosecond fluctuation of the beta-strand D on which His51 exists, which is propagated to the core of the molecule, thus promoting the global and slow fluctuations. This may contribute to the overall destabilization of the molecule, increasing the equilibrium population of the amyloidogenic intermediate.     

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.     

Kinetics and mechanism of the acid transition of the active site in plastocyanin

Exchange on the microsecond time scale between the protonated and deprotonated forms of His92 in the copper site of reduced plastocyanin from the cyanobacteria Anabaena variabilis was monitored using 15N NMR relaxation measurements. On the basis of the dependence of the kinetics on pH and phosphate buffer concentration, we propose a two-step model for the protonation of the copper site in agreement with previous crystallographic studies. It is shown that the proton transfer is the rate-limiting step in the reaction at low buffer concentrations, whereas at high buffer concentrations, another step becomes rate-limiting. We suggest that the latter step is a concerted dissociation of His92 from the Cu(I) ion and a 180 degrees rotation of the imidazole ring, which precede the protonation. The first-order rate constant for the dissociation of His92 from the Cu(I) ion is estimated to be 2.4 x 10(4) s(-1). Also, a cooperative effect of the protonation of the remote His61 on the protonation of His92 and the redox properties of the protein was investigated by substituting His61 with asparagine. The mutation causes a modest change in both the pKa value of His92 and the redox potential of the protein.     

Peptide binding to HLA-DP proteins at pH 5.0 and pH 7.0: a quantitative molecular docking study

ABSTRACT: BACKGROUND: HLA-DPs are class II MHC proteins mediating immune responses to many diseases. Peptides bind MHC class II proteins in the acidic environment within endosomes. Acidic pH markedly elevates association rate constants but dissociation rates are almost unchanged in the pH range 5.0 - 7.0. This pH-driven effect can be explained by the protonation/deprotonation states of Histidine, whose imidazole has a pKa of 6.0. At pH 5.0, imidazole ring is protonated, making Histidine positively charged and very hydrophilic, while at pH 7.0 imidazole is unprotonated, making Histidine less hydrophilic. We develop here a method to predict peptide binding to the four most frequent HLA-DP proteins: DP1, DP41, DP42 and DP5, using a molecular docking protocol. Dockings to virtual combinatorial peptide libraries were performed at pH 5.0 and pH 7.0. RESULTS: The X-ray structure of the peptide - HLA-DP2 protein complex was used as a starting template to model by homology the structure of the four DP proteins. The resulting models were used to produce virtual combinatorial peptide libraries constructed using the single amino acid substitution (SAAS) principle. Peptides were docked into the DP binding site using AutoDock at pH 5.0 and pH 7.0. The resulting scores were normalized and used to generate Docking Score-based Quantitative Matrices (DS-QMs). The predictive ability of these QMs was tested using an external test set of 484 known DP binders. They were also compared to existing servers for DP binding prediction. The models derived at pH 5.0 predict better than those derived at pH 7.0 and showed significantly improved predictions for three of the four DP proteins, when compared to the existing servers. They are able to recognize 50% of the known binders in the top 5% of predicted peptides. CONCLUSIONS: The higher predictive ability of DS-QMs derived at pH 5.0 may be rationalised by the additional hydrogen bond formed between the backbone carbonyl oxygen belonging to the peptide position before p1 (p-1) and the protonated -nitrogen of His79beta. Additionally, protonated His residues are well accepted at most of the peptide binding core positions which is in a good agreement with the overall negatively charged peptide binding site of most MHC proteins.     

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.     

Protein-Protein Docking with Dynamic Residue Protonation States

Protein-protein interactions depend on a host of environmental factors. Local pH conditions influence the interactions through the protonation states of the ionizable residues that can change upon binding. In this work, we present a pH-sensitive docking approach, pHDock, that can sample side-chain protonation states of five ionizable residues (Asp, Glu, His, Tyr, Lys) on-the-fly during the docking simulation. pHDock produces successful local docking funnels in approximately half (79/161) the protein complexes, including 19 cases where standard RosettaDock fails. pHDock also performs better than the two control cases comprising docking at pH 7.0 or using fixed, predetermined protonation states. On average, the top-ranked pHDock structures have lower interface RMSDs and recover more native interface residue-residue contacts and hydrogen bonds compared to RosettaDock. Addition of backbone flexibility using a computationally-generated conformational ensemble further improves native contact and hydrogen bond recovery in the top-ranked structures. Although pHDock is designed to improve docking, it also successfully predicts a large pH-dependent binding affinity change in the Fc–FcRn complex, suggesting that it can be exploited to improve affinity predictions. The approaches in the study contribute to the goal of structural simulations of whole-cell protein-protein interactions including all the environmental factors, and they can be further expanded for pH-sensitive protein design.     

Conformational and dynamics changes induced by bile acids binding to chicken liver bile acid binding protein

The correlation between protein motions and function is a central problem in protein science. Several studies have demonstrated that ligand binding and protein dynamics are strongly correlated in intracellular lipid binding proteins (iLBPs), in which the high degree of flexibility, principally occurring at the level of helix-II, CD, and EF loops (the so-called portal area), is significantly reduced upon ligand binding. We have recently investigated by NMR the dynamic properties of a member of the iLBP family, chicken liver bile acid binding protein (cL-BABP), in its apo and holo form, as a complex with two bile salts molecules. Binding was found to be regulated by a dynamic process and a conformational rearrangement was associated with this event. We report here the results of molecular dynamics (MD) simulations performed on apo and holo cL-BABP with the aim of further characterizing the protein regions involved in motion propagation and of evaluating the main molecular interactions stabilizing bound ligands. Upon binding, the root mean square fluctuation values substantially decrease for CD and EF loops while increase for the helix-loop-helix region, thus indicating that the portal area is the region mostly affected by complex formation. These results nicely correlate with backbone dynamics data derived from NMR experiments. Essential dynamics analysis of the MD trajectories indicates that the major concerted motions involve the three contiguous structural elements of the portal area, which however are dynamically coupled in different ways whether in the presence or in the absence of the ligands. Motions of the EF loop and of the helical region are part of the essential space of both apo and holo-BABP and sample a much wider conformational space in the apo form. Together with NMR results, these data support the view that, in the apo protein, the flexible EF loop visits many conformational states including those typical of the holo state and that the ligand acts stabilizing one of these pre-existing conformations. The present results, in agreement with data reported for other iLBPs, sharpen our knowledge on the binding mechanism for this protein family.     

Sequence-specific DNA binding by the glucocorticoid receptor DNA-binding domain is linked to a salt-dependent histidine protonation

We used isothermal titration calorimetry in the temperature range 21-25 degrees C to investigate the effect of pH on the calorimetric enthalpy (delta H(cal)) for sequence specific DNA-binding of the glucocorticoid receptor DNA-binding domain (GR DBD). Titrations were carried out in solutions containing 100 mM NaCl, 1 mM dithiothreitol, 5% glycerol by volume, and 20 mM Tris, Hepes, Mops, or sodium phosphate buffers at pH 7.5. A strong dependence of delta H(cal) on the buffer ionization enthalpy is observed, demonstrating that the DNA binding of the GR DBD is linked to proton uptake at these conditions. The apparent increase in the pK(a) for an amino acid side chain upon DNA binding is supported by the results of complementary titrations, where delta H(cal) shows a characteristic dependence on the solution pH. delta H(cal) is also a function of the NaCl concentration, with opposite dependencies in Tris and Hepes buffers, respectively, such that a similar delta H(cal) value is approached at 300 mM NaCl. This behavior shows that the DNA-binding induced protonation is inhibited by increased concentrations of NaCl. A comparison with structural data suggests that the protonation involves a histidine (His451) in the GR DBD, because in the complex this residue is located close to a DNA phosphate at an orientation that is consistent with a charged-charged hydrogen bond in the protonated state. NMR spectra show that His451 is not protonated in the unbound protein at pH 7.5. The pH dependence in delta H(cal) can be quantitatively described by a shift of the pK(a) of His451 from approximately 6 in the unbound state to close to 8 when bound to DNA at low salt concentration conditions. A simple model involving a binding competition between a proton and a Na(+) counterion to the GR DBD-DNA complex reproduces the qualitative features of the salt dependence.