The deacylation of substituted benzol-alpha-chymotrypsins.

An extensive comparison of deacylation rates of mono- and disubstituted benzoyl-α-chymotrypsins indicates that no steric effects on rate or apparent pK_a of deacylation are detectable within this series. Some anomalous effects on deacylation rate appear to be associated with fluoro- and nitro-substituents in particular positions on the ring and may be attributable to specific interactions at the enzyme active site. The extensive series of structurally similar acyl-enzymes prepared has allowed a thorough analysis of the effect of acyl group pK_a on the apparent pK_a of deacylation. The data indicates that polar effects on the apparent pK_a are probably negligible. Rho for the deacylation reaction is in good agreement with model reactions for an imidazole general base-catalyzed model reaction.     

Determination of 1-aryl-4-propylpiperazine pKa values: The substituent on aryl modulates basicity

In order to design a potential drug, it is important to know its pK_a because the protonation state of the molecule will be critical for ligand–receptor interaction and for the pharmacokinetic of the molecule. pK_a values of a series of 1-(substitutedphenyl)-4-propylpiperazines were measured to study how the presence of a substituent on the phenyl ring modulates the basicity of N-4 nitrogen. pK_a values indicated that the position of the substituent was crucial. In general, the introduction of the substituent in ortho-position of the phenyl ring increased the basicity of the molecule. This effect appeared to be related to steric and conformational effects and not to the electronic properties of the substituent. On the other hand, meta- and para-substituted derivatives showed a slight decrease of pK_a that was qualitatively consistent with the electronic properties of the substituent.     

Quantifying the electronic cis effect of phosphine, arsine and stibine ligands by use of rhodium(I) Vaska-type complexes

The cis effects of phosphine, arsine and stibine ligands have been evaluated by measuring the IR stretching frequency in dichloromethane of the carbonyl ligand in a series of Rh(I) Vaska-type complexes, trans-[RhCl(CO)(L)₂]. These data were correlated with those obtained by Tolman for the electronic trans influences in the [Ni(L)(CO)₃] complexes. The electronic contribution, χ_Fc, of ferrocenyl was determined as 0.8 from these plots by evaluating PPh₂Fc as ligand. In order to accommodate arsine and stibine ligands an additional correction term, to compensate for differences in the donor atom, was added to Tolman’s equation for calculation of the Tolman electronic parameter of phosphine ligands. In the resulting equation: ν(CO_Ni)=2056.1+∑_i=1³χ_i+C_L values for C_L of C_P=0, C_As=−1.5 and C_Sb=−3.1 are suggested for phosphine, arsine and stibine ligands, respectively. The crystal and molecular structures of trans-[RhCl(CO)(PPh₂Fc)₂] · 2C₆H₆, trans-[RhCl(CO){P(NMe₂)₃}₂] and trans-[RhCl(CO)(AsPh₃)₂] are reported. The Tolman cone angles for PPh₂Fc and P(NMe₂)₃ were determined as 169° and 166°, while the effective cone angles for PPh₂Fc, P(NMe₂)₃ and AsPh₃ were determined as 171°, 168° and 147°, respectively.     

Kinetics and mechanism of the anation reactions of aquocobaloximes containing bulky phosphine ligands

The anation reactions of phosphino-aquocobaloximes containing bulky P donor ligands with HSO₃⁻ lead through a multistep mechanism to the formation of SO₃Co(DH)₂SO₃³⁻; the first step leads to the water substitution. The successive steps involve the hydrolysis of the phosphine and the ligation of HSO₃⁻, the dissociation of phosphine being the rate determining step.The steric bulk of the P ligand appears to be the factor determining the phosphine lability.     

How enzymes harness highly unfavorable proton transfer reactions

Acid–base reactions that are exceedingly unfavorable under standard conditions can be catalytically important at enzyme active sites. For example, in triose phosphate isomerase, a glutamate side chain (nominal pK _a ≈ 4 in solution) can in fact deprotonate a CH group that is vicinal to a carbonyl (pK _a ≈ 18 in solution). This is true because of three distinct interactions: (a) ground state pK _a shifts due to environment polarity and electrostatics; (b) dramatic increases in effective molarity due to optimization of proximity and orientation; and (c) transition state pK _a shifts due to binding interactions and the formation of strong low barrier hydrogen bonds. In this report, we review the literature showing that the sum of these three effects supplies more than enough free energy to push forward proton transfer reactions that under standard conditions are exceedingly nonspontaneous and slow.     

Redox potentials of trinuclear μ-oxo ruthenium acetate clusters with N-heterocyclic ligands

The trinuclear clusters of general composition [Ru₃O(OOCCH₃)₆(N-Het)₃], where N-Het=pyridine and pyrazine derivatives, exhibit a series of reversible waves in the range of −1.8 to 2.4 V versus SHE, in acetonitrile, ascribed to the successive [cluster]^{−2/−1/0/+1/+2/+3} redox couples. The redox potentials decrease with the pK_a of the N-heterocyclic ligands according to the equations E°(+3/+2)= 2.24−0.023 pK_a; E°(+2/+1)=1.34−0.029 pK_a; E°(+1/0)=0.36−0.039 pK_a and E°(0/−1)=−0.68− 0.074 pK_a. The dependence is greater at lower oxidation states, reflecting the role of π-backbonding in the complexes.     

Metal ion-biomolecule interactions. II. Methylmercuration of the deprotonated amino groups in adenine,guanine, and cytosine derivatives,and its relationship to amino group acidity

Proton nmr spectroscopic evidence is presented for methylmercury(II) binding to the deprotonated amino groups in adenosine, 9-methyladenine, guanosine, 1-methylguanosine, and cytidine under basic conditions. Except for the guanosine case, ¹H nmr spectra of the products from aqueous or ethanolic 1:1 mixtures of substrate and MeHgOH are consistent with methylmercuration of the deprotonated amino groups. Guanosine undergoes initial binding of MeHg to N₁, and a second equivalent of MeHgOH is necessary to effect amino binding. The nmr spectra of the complexed adenine derivatives suggest that different geometrical isomers exist in (CD₃)₂SO solution, reflecting the partial double bond character of the C₆N bond in these systems. Using a correlation relating the magnitude of the ¹⁹⁹Hg-¹H coupling constant (J) for MeHg-ligand complexes with the ligand pK_a (J = ?3.88 pK_a + 248.5, extending over 13 pK units, based on a variety of N and O donor ligands), estimates (± 0.3 pK unit) of the pK_as of the amino groups of the above substrates have been made. In this way, pK_a values of 15.5 (cytidine), 17.0 (adenosine and 9-methyladenine), 15.1 (guanosine), and 14.9 (1-methylguanosine) are obtained. In the cases where comparisons with literature pK_a data can be made, good agreement is found.     

Modulation of an Active-Site Cysteine pKa Allows PDI to Act as a Catalyst of both Disulfide Bond Formation and Isomerization

Protein disulfide isomerase (PDI) plays a central role in disulfide bond formation in the endoplasmic reticulum. It is implicated both in disulfide bond formation and in disulfide bond reduction and isomerization. To be an efficient catalyst of all three reactions requires complex mechanisms. These include mechanisms to modulate the pK_a values of the active-site cysteines of PDI. Here, we examined the role of arginine 120 in modulating the pK_a values of these cysteines. We find that arginine 120 plays a significant role in modulating the pK_a of the C-terminal active-site cysteine in the a domain of PDI and plays a role in determining the reactivity of the N-terminal active-site cysteine but not via direct modulation of its pK_a. Mutation of arginine 120 and the corresponding residue, arginine 461, in the a′ domain severely reduces the ability of PDI to catalyze disulfide bond formation and reduction but enhances the ability to catalyze disulfide bond isomerization due to the formation of more stable PDI-substrate mixed disulfides. These results suggest that the modulation of pK_a of the C-terminal active cysteine by the movement of the side chain of these arginine residues into the active-site locales has evolved to allow PDI to efficiently catalyze both oxidation and isomerization reactions.     

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

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

Effects of alcohol/water mixtures on the structure and reactivity of cytochrome c

The oxidation reaction of ferrocytochrome c (produced in situ by pulse radiolysis) by Fe(CN)^3?₆, was used to probe the effect of alcohol/water mixtures on the reactivity of the protein. Reduced cytochrome c is oxidized in a biphasic process. The relative contribution of each phase depended on: pH, alcohol concentration and temperature. pK_a values were derived from the kinetic data. These pK_a values were identical with the spectroscopic pK_a values determined under similar conditions by monitoring the 695 nm absorption band of the oxidized protein. The two phases of oxidation were therefore related to the oxidation of a relaxed and a nonrelaxed conformer of reduced cytochrome c produced in situ. A shift in the pK_a of ferricytochrome c and a retardation of the redox reactions of both the reduced and the oxidized protein were observed at low alcohol concentrations (up to 5 mol %). These low alcohol concentrations are known to affect the structure of water (Yaacobi, M. and Ben-Naim, A. (1973) J. Sol. Chem. 5, 425?443; Ben-Naim, A. (1967) J. Phys. Chem. 71, 4002?4007 and Ben-Naim, A. and Baer, S. (1964) Trans. Faraday Soc. 60, 1736?1741) but have only minor effects on the protein. Accordingly, the kinetic results are interpreted on the basis of involvement of water molecules in the reaction complex of cytochrome c with its redox substrates.     

Acidic elements in histamine H3 receptor antagonists

Antagonists of the human histamine H₃ receptor (hH₃R) often contain a second basic moiety, which is well known to boost affinity on this histamine receptor subtype. Here, we prepared compounds with acidic moieties of different pK_a values to figure out that the hH₃R tolerates these functionalities when added to a common pharmacophore blueprint. Depending on the acidic, electronic and steric features the designed ligands showed hH₃R affinities in the nanomolar concentration range. Additionally, selected ligands were tested but failed as dual acting hH₃R/hPPAR (human peroxisome proliferator-activated receptor) ligands.     

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.     

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.     

Heme-linked protonation of HCN, CO, NO and O2 complexes of reduced horseradish peroxidases.

Small reversible changes in the absorption spectra of HCN, CO, NO and O₂ complexes of ferrous diacetyldeuteroperoxidase A, not hitherto observed, were attributed to proton dissociation of a distal amino acid residue. From spectrophotometric titration data the pK_a was measured as 5.5 (HCN), 5.6 (ligand free), 6.0 (CO), 6.55 (NO) and 8.0 (O₂). The value of 8.0 for the pK_a of the O₂ complex was also obtained from a curve of pH dependence of proton uptake in the reaction of the ferrous enzyme with O₂. Absorption bands in the visible region were shifted to longer wavelengths in the order of CO to NO to O₂ which is the decreasing order of the energy of π¹ level of these diatomic ligands.The pK_a values for CO complexes of ferroperoxidases, isoenzymes A and (B+C) were varied with substituents at the 2 and 4 positions of deuterohemin IX, and the $\text{Δp}\text{K}{}_{\mathrm{<mtext>a</mtext>}}\text{Δp}\text{K}{}_3$ ratio was about 0.3 in both series of isoenzyme preparations, where pK₃ is a measure of basicity of pyrrole nitrogen.The present data support the previous conclusion (Yamada and Yamazaki (1974) Arch. Biochem. Biophys., 165, 728) that the pK_a for ferroperoxidases, measured from small reversible changes in the absorption spectra, represents a proton dissociation constant of a distal amino acid residue and that there is hydrogen bonding between the residue and a ligand atom directly bound to the iron atom.     

Retinal-Protein Interactions in Halorhodopsin from Natronomonas pharaonis: Binding and Retinal Thermal Isomerization Catalysis

Halorhodopsin from Natronomonas pharaonis (NpHR) is a member of the retinal protein group and serves as a light-driven chloride pump in which chloride ions are transported through the membrane following light absorption by the retinal chromophore. In this study, we examined two main issues: (1) factors controlling the binding of the retinal chromophore to the NpHR opsin and (2) the ability of the NpHR opsin to catalyze the thermal isomerization of retinal isomers. We have revealed that the reconstitution process of pharaonis HR (NpHR) pigment from its apoprotein and all-trans retinal depends on the pH, and the process has a pK_a of 5.8 ± 0.1. It was proposed that this pK_a is associated with the pK_a of the lysine residue that binds the retinal chromophore (Lys256). The pigment formation is regulated by the concentration of sodium chloride, and the maximum yield was observed at 3.7 M NaCl. The low yield of pigment in a lower concentration of NaCl (< 3 M) may be due to an altered conformation adopted by the apomembrane, which is not capable of forming the pigment. Unexpectedly and unlike the apomembrane of bacteriorhodopsin, NpHR opsin produces pigments with 11-cis retinal and 9-cis retinal owing to the thermal isomerization of these retinal isomers to all-trans retinal. The isomerization rate depends on the pH, and it is faster at a higher pH. The pK_a value of the isomerization process is similar to the pK_a of the binding process of these retinals, which suggests that Lys256 is also involved in the isomerization process. The isomerization is independent of the sodium chloride concentration. However, in the absence of sodium chloride, the apoprotein adopts such a conformation, which does not prevent the isomerization of retinal, but it prevents a covalent bond formation with the lysine residue. The rate and the thermodynamic parameter analysis of the retinal isomerization by NpHR apoprotein led to the conclusion that the apomembrane catalyzes the isomerization via a triplet mechanism.     

The pKa values of two histidine residues in human haemoglobin, the Bohr effect, and the dipole moments of alpha-helices

Studies of abnormal and chemically modified haemoglobins indicate that in 0.1 m-NaCl about 40% of the alkaline Bohr effect of human haemoglobin is contributed by the C-terminal histidine HC3(146)β. In deoxyhaemoglobin, the imidazole of this histidine forms a salt bridge with aspartate FG1(94)β, in oxyhaemoglobin or carbonmonoxyhaemoglobin it accepts a hydrogen bond from its own NH group instead. Kilmartin et al. (1973) showed that in 0.2 m-NaCl + 0.2 m-phosphate this change of ligation lowered the pK_a of the histidine from 8.0 in Hb³ to 7.1 in HbCO, but Russu et al. (1980) claimed that in bis-Tris buffer without added NaCl its pK_a in HbCO dropped no lower than 7.85, and that in this medium the C-terminal histidine made only a negligible contribution to the alkaline Bohr effect.We have compared the histidine resonances of HbCO A with those of three abnormal haemoglobins: HbCO Cowtown (His HC3(146)β → Leu), HbCO Wood (His FG4(97)β → Leu) and HbCO Malmø (His FG4(97)β → Gln). Our results show that the resonance assigned by Russu et al. to His HC3(146)β in fact belongs to His FG4(97)β. Although in Hb the pK_a of His HC3(146)β is 8.05 ± 0.05 independent of ionic strength, in HbCO its pK_a drops sharply with diminishing ionic strength, so that in the buffer employed by Russu et al. it has a pK_a of 6.2 and makes a contribution to the alkaline Bohr effect that is 57% larger than in the phosphate buffer employed by Kilmartin et al. (1973).In HbCO A, His FG4(97)β does not contribute to the Bohr effect, but in HbCO from which His HC3(146)β has been cleaved (HbCO des-His), His FG4(97)β is in equilibrium between two conformations with different pK_a values. This equilibrium varies with ionic strength and pH, and presumably also with degree of ligation of the haem moiety.In HbCO A, His FG4(97)β has a pK_a of 7.8 compared to the pK_a value of about 6.6 characteristic of free histidines at the surface of proteins. This high pK_a is accounted for by its interaction with the negative pole at the C terminus of helices F and FG. It corresponds to a free energy change of the same order as that observed in the interaction of histidines with carboxylate ions and confirms the strongly dipolar character of α-helices, which manifests itself even when they lie on the surface of the protein.     

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

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

Mechanistic informations from trans-effect in the anation reactions of some bis(dimethylglyoximato)Rh(III) complexes

The parallel lability trend for the anation reactions of RRh(DH)₂H₂O (RCH₃, CH₃CH₂, CF₃CH₂, ClCH₂) and RCo(DH)₂H₂O complexes suggests a dissociative activation process for the reactions of the organorhodoximes. The high lability of these complexes, arising from the stabilization of the transition state, is not entirely due to the trans-effect of the R group, but, at least partially, to the labilizing effect of the equatorial macrocycle.     

Ionization constants and thermodynamic parameters of histidine and derivatives

The pK_a values for the proton dissociation of carboxyl, imidazolium, and ammonium groups for histidine and ten of its derivatives were determined electrometrically at seven temperatures in the range 10–40°C. The ΔH and ΔS values were estimated from the temperature dependence of the dissociation constants of histidine and its derivatives. These results and the pK_a values compared in terms of inductive effect suggest an ion-dipole interaction between the protonated amino group and the unprotonated imidazole ring. The charge and the solvation effects of the neighboring groups are the main factors that determine the imidazole group pK_a in histidine and its studied derivatives. The N^τ-H tautomer is favored over the N^π-H by 1.6 kcal/mol, indicating that the inductive substituent effect at position 4 of the imidazole ring is the major component in determining this tautomeric preference.     

Origin of the pKa shift of the catalytic lysine in acetoacetate decarboxylase

The pK_a value of Lys115, the catalytic residue in acetoacetate decarboxylate, was calculated using atomic coordinates of the X-ray crystal structure with consideration of the protonation states of all titratable sites in the protein. The calculated pK_a value of Lys115 (pK_a(Lys115)) was unusually low (≈6) in agreement with the experimentally measured value. Although charged residues impact pK_a(Lys115) considerably in the native protein, the significant pK_a(Lys115) downshift in the protein with respect to aqueous solution was mainly due to loss of the solvation energy in the catalytic active site relative to bulk water.