Optimizing the Michaelis complex of trimethylamine dehydrogenase: identification of interactions that perturb the ionization of substrate and facilitate catalysis with trimethylamine base.

Recent evidence from isotope studies supports the view that catalysis by trimethylamine dehydrogenase (TMADH) proceeds from a Michaelis complex involving trimethylamine base and not, as thought previously, trimethylammonium cation. In native TMADH reduction of the flavin by substrate (perdeuterated trimethylamine) is influenced by two ionizations in the Michaelis complex with pK(a) values of 6.5 and 8.4; maximal activity is realized in the alkaline region. The latter ionization has been attributed to residue His-172 and, more recently, the former to the ionization of substrate itself. In the Michaelis complex, the ionization of substrate (pK(a) approximately 6.5 for perdeuterated substrate) is perturbed by approximately -3.3 to -3.6 pH units compared with that of free trimethylamine (pK(a) = 9.8) and free perdeuterated trimethylamine (pK(a) = 10.1), respectively, thus stabilizing trimethylamine base by approximately 2 kJ mol(-1). We show, by targeted mutagenesis and stopped-flow studies that this reduction of the pK(a) is a consequence of electronic interaction with residues Tyr-60 and His-172, thus these two residues are key for optimizing catalysis in the physiological pH range. We also show that residue Tyr-174, the remaining ionizable group in the active site that we have not targeted previously by mutagenesis, is not implicated in the pH dependence of flavin reduction. Formation of a Michaelis complex with trimethylamine base is consistent with a mechanism of amine oxidation that we advanced in our previous computational and kinetic studies which involves nucleophilic attack by the substrate nitrogen atom on the electrophilic C4a atom of the flavin isoalloxazine ring. Stabilization of trimethylamine base in the Michaelis complex over that in free solution is key to optimizing catalysis at physiological pH in TMADH, and may be of general importance in the mechanism of other amine dehydrogenases that require the unprotonated form of the substrate for catalysis.     

Effects of mutations of the active site arginine residues in 4-oxalocrotonate tautomerase on the pKa values of active site residues and on the pH dependence of catalysis.

The unusually low pK(a) value of the general base catalyst Pro-1 (pK(a) = 6.4) in 4-oxalocrotonate tautomerase (4-OT) has been ascribed to both a low dielectric constant at the active site and the proximity of the cationic residues Arg-11 and Arg-39 [Stivers, J. T., Abeygunawardana, C., Mildvan, A. S., Hajipour, G., and Whitman, C. P. (1996) Biochemistry 35, 814-823]. In addition, the pH-rate profiles in that study showed an unidentified protonated group essential for catalysis with a pK(a) of 9.0. To address these issues, the pK(a) values of the active site Pro-1 and lower limit pK(a) values of arginine residues were determined by direct (15)N NMR pH titrations. The pK(a) values of Pro-1 and of the essential acid group were determined independently from pH-rate profiles of the kinetic parameters of 4-OT in arginine mutants of 4-OT and compared with those of wild type. The chemical shifts of all of the Arg Nepsilon resonances in wild-type 4-OT and in the R11A and R39Q mutants were found to be independent of pH over the range 4.9-9.7, indicating that no arginine is responsible for the kinetically determined pK(a) of 9.0 for an acidic group in free 4-OT. With the R11A mutant, where k(cat)/K(m) was reduced by a factor of 10(2.9), the pK(a) of Pro-1 was not significantly altered from that of the wild-type enzyme (pK(a) = 6.4 +/- 0.2) as revealed by both direct (15)N NMR titration (pK(a) = 6.3 +/- 0.1) and the pH dependence of k(cat)/K(m) (pK(a) = 6.4 +/- 0.2). The pH-rate profiles of both k(cat)/K(m) and k(cat) for the reaction of the R11A mutant with the dicarboxylate substrate, 2-hydroxymuconate, showed humps, i.e., sharply defined maxima followed by nonzero plateaus. The humps disappeared in the reaction with the monocarboxylate substrate, 2-hydroxy-2,4-pentadienoate, indicating that, unlike the wild-type enzyme which reacts only with the dianionic form of the dicarboxylic substrate, the R11A mutant reacts with both the 6-COOH and 6-COO(-) forms, with the 6-COOH form being 12-fold more active. This reversal in the preferred ionization state of the 6-carboxyl group of the substrate that occurs upon mutation of Arg-11 to Ala provides strong evidence that Arg-11 interacts with the 6-carboxylate of the substrate. In the R39Q mutant, where k(cat)/K(m) was reduced by a factor of 10(3), the kinetically determined pK(a) value for Pro-1 was 4.6 +/- 0.2, while the ionization of Pro-1 showed negative cooperativity with an apparent pK(a) of 7.1 +/- 0.1 determined by 1D (15)N NMR. From the Hill coefficient of 0.54, it can be shown that the apparent pK(a) value of 7.1 could result most simply from the averaging of two limiting pK(a) values of 4.6 and 8.2. Mutation of Arg-39, by altering the structure of the beta-hairpin which covers the active site, could result in an increase in the solvent exposure of Pro-1, raising its upper limit pK(a) value to 8.2. In the R39A mutant, the kinetically determined pK(a) of Pro-1 was also low, 5.0 +/- 0.2, indicating that in both the R39Q and R39A mutants, only the sites with low pK(a) values were kinetically operative. With the fully active R61A mutant, the kinetically determined pK(a) of Pro-1 (pK(a) = 6.5 +/- 0.2) agreed with that of wild-type 4-OT. It is concluded that the unusually low pK(a) of Pro-1 shows little contribution from electrostatic effects of the nearby cationic Arg-11, Arg-39, and Arg-61 residues but results primarily from a site of low local dielectric constant.     

Expression, purification, and kinetic characterization of full-length human fibroblast activation protein

Human fibroblast activation protein (FAP), an integral membrane serine protease, was produced in insect cells as a hexa-His-tagged protein using a recombinant baculovirus expression system. Two isoforms of FAP, glycosylated and nonglycosylated, were identified by Western blotting using an anti-His-tag antibody and separated by lectin chromatography. The glycosylated FAP was purified to near homogeneity using immobilized metal affinity chromatography and was shown to have both postprolyl dipeptidyl peptidase and postgelatinase activities. In contrast, the nonglycosylated isoform demonstrated no detectable gelatinase activity by either zymography or a fluorescence-based gelatinase activity assay. The kinetic parameters of the dipeptidyl peptidase activity for glycosylated FAP were determined using dipeptide Ala-Pro-7-amino-trifluoromethyl-coumarin as the substrate. The k(cat) is 2.0 s(-1) and k(cat)/K(m) is 1.0 x 10(4) M(-1) s(-1) at pH 8.5. The pH dependence of k(cat) reveals two ionization groups with pK(a1) of 7.0 and pK(a2) of 11.0. The pH profile of k(cat)/K(m) yields similar results with pK(a1) 6.2 and pK(a2) 11.0. The neutral pK(a1) is associated with His at the active site. The basic pK(a2) might be contributed from an ionization group that is not involved directly in catalysis, instead associated with the stability of the active site structure.     

Spectrophotometric titration of phenolic groups of pepsin.

The ionization of tyrosine residues in diazotized pepsin under various solvent conditions was studied. All tyrosyl residues of the protein titrated normally with a pK of 10.02 in 6 M guanidine hydrochloride solution. On the other hand, two stages in the phenolic group titration curve were observed for the inactivated protein in the absence of guanidine hydrochloride; only about 10 tyrosine residues ionized reversibly up to pH 11, above which titration was irreversible. The irreversible titration zone corresponds to the pH range 11--13 in which unfolding, leading to the random coil state, was shown to occur by circular dichroism and viscosity measurements. The number of tyrosine residues exposed in the native and alkali-denatured (pH 7.5) states of diazotized protein were also studied by solvent perturbation techniques; 10 and 12 groups are exposed in the native and denatured states, respectively.     

Amino acid residues involved in the catalytic mechanism of NAD-dependent glutamate dehydrogenase from Halobacterium salinarum

The pH dependence of kinetic parameters for a competitive inhibitor (glutarate) was determined in order to obtain information on the chemical mechanism for NAD-dependent glutamate dehydrogenase from Halobacterium salinarum. The maximum velocity is pH dependent, decreasing at low pHs giving a pK value of 7.19+/-0.13, while the V/K for l-glutamate at 30 degrees C decreases at low and high pHs, yielding pK values of 7.9+/-0.2 and 9.8+/-0.2, respectively. The glutarate pKis profile decreases at high pHs, yielding a pK of 9. 59+/-0.09 at 30 degrees C. The values of ionization heat calculated from the change in pK with temperature are: 1.19 x 10(4), 5.7 x 10(3), 7 x 10(3), 6.6 x 10(3) cal mol-1, for the residues involved. All these data suggest that the groups required for catalysis and/or binding are lysine, histidine and tyrosine. The enzyme shows a time-dependent loss in glutamate oxidation activity when incubated with diethyl pyrocarbonate (DEPC). Inactivation follows pseudo-first-order kinetics with a second-order rate constant of 53 M-1min-1. The pKa of the titratable group was pK1=6.6+/-0.6. Inactivation with ethyl acetimidate also shows pseudo-first-order kinetics as well as inactivation with TNM yielding second-order constants of 1.2 M-1min-1 and 2.8 M-1min-1, and pKas of 8.36 and 9.0, respectively. The proposed mechanism involves hydrogen binding of each of the two carboxylic groups to tyrosyl residues; histidine interacts with one of the N-hydrogens of the l-glutamate amino group. We also corroborate the presence of a conservative lysine that has a remarkable ability to coordinate a water molecule that would act as general base.     

Chemical modification of cholesterol-hydroxylating cytochrome P450 from bovine adrenal cortex mitochondria by acetylimidazole

The state and role of tyrosine residues in adrenocortical cytochrome P-450scc was investigated. Spectrophotometric titration experiments showed the existence of two types of tyrosine residues in the hemoprotein molecule, i.e., exposed (pK 9.25) and buried (pK 10.75) ones. Chemical modification of the exposed tyrosine residues with N-acetylimidazole was carried out. The data obtained indicate the involvement of these residues in the interaction with adrenodoxin.     

Application of the pH-jump method to the titration of tyrosine residues in bovine alpha-lactalbumin.

A stopped-flow technique has been developed for the zero-time spectrophotometric titration of tyrosine residues in the purely native or in the purely alkaline denatured state of alpha-lactalbumin that undergoes an alkaline conformational transition in the pH region of tyrosine ionization. The progressive absorption change at 298 nm caused by a pH jump from neutral pH is shown to result from the change in ionization of the tyrosine residues brought about by a first-order process of the conformational transition. Extrapolation to zero time gives the titration curve for purely native alpha-lactalbumin. Similarly, the pH jump from highly alkaline pH gives the titration curve for the purely alkaline denatured protein. The method should be generally applicable to other proteins that contain tyrosines. Analysis of the titration curves suggests that the four tyrosines in native alpha-lactalbumin have pK values of 10.5, 11.8, 11.8, and 12.7, respectively. After the alkaline transconformation, all of them become titrated normally with a pK value of 10.3. A comparison of these results with the ionization behavior of tyrosines in hen egg white and human lysozymes is presented and discussed in terms of differences in the sequences of the proteins.     

Assignment of 1H NMR resonances of histidine and other aromatic residues in met-, cyano-, oxy-, and (carbon monoxy)myoglobins

The resolved 1H NMR resonances of the aromatic region in the 270-MHz NMR spectrum of sperm whale, horse, and pig metmyoglobin (metMb) have been assigned, including the observable H-2 and H-4 histidine resonances, the tryptophan H-2 resonances, and upfield-shifted resonances from one tyrosine residue. The use of different Mb species, carboxymethylation, and matching of pK values allows the assignment of the H-4 resonances, which agree in only three cases out of seven with scalar-correlated two-dimensional NMR spectroscopy assignments by others. The conversion to hydroxymyoglobin at high pH involves rearrangements throughout the molecule and is observed by many assigned residues. In sperm whale ferric cyanomyoglobin, nine H-2 and eight H-4 histidine resonances have been assigned, including the His-97 H-2 resonance and tyrosine resonances from residues 103 and 146. The hyperfine-shifted resonances from heme and near-heme protons observe a shift with a pK = 5.3 +/- 0.3 (probably due to deprotonation of His-97, pK = 5.6) and another shift at pK = 10.8 +/- 0.3. The spectrum of high-spin ferrous sperm whale deoxymyoglobin is very similar to that of metMb, which allows the assignment of seven surface histidine H-2 and H-4 resonances and also resonances from the two tryptophan residues and one tyrosine. In diamagnetic sperm whale (carbon monoxy)myoglobin (COMb), 10 His H-2 and 11 His H-4 resonances are observed, and 8 H-2 and 9 H-4 resonances are assigned, including His-64 H-4, the distal histidine. This important resonance is not observed in sperm whale oxymyoglobin, which in general shows very similar titration curves to COMb. Histidine-36 shows unusual titration behavior in the paramagnetic derivatives but normal behavior in the diamagnetic derivatives, which is discussed in the accompanying paper [Bradbury, J. H., & Carver, J. A. (1984) Biochemistry (following paper in this issue)].     

Experimental pK(a) values of buried residues: analysis with continuum methods and role of water penetration

Lys-66 and Glu-66, buried in the hydrophobic interior of staphylococcal nuclease by mutagenesis, titrate with pK(a) values of 5.7 and 8.8, respectively (Dwyer et al., Biophys. J. 79:1610-1620; García-Moreno E. et al., Biophys. Chem. 64:211-224). Continuum calculations with static structures reproduced the pK(a) values when the protein interior was treated with a dielectric constant (epsilon(in)) of 10. This high apparent polarizability can be rationalized in the case of Glu-66 in terms of internal water molecules, visible in crystallographic structures, hydrogen bonded to Glu-66. The water molecules are absent in structures with Lys-66; the high polarizability cannot be reconciled with the hydrophobic environment surrounding Lys-66. Equilibrium thermodynamic experiments showed that the Lys-66 mutant remained folded and native-like after ionization of the buried lysine. The high polarizability must therefore reflect water penetration, minor local structural rearrangement, or both. When in pK(a) calculations with continuum methods, the internal water molecules were treated explicitly, and allowed to relax in the field of the buried charged group, the pK(a) values of buried residues were reproduced with epsilon(in) in the range 4-5. The calculations show that internal waters can modulate pK(a) values of buried residues effectively, and they support the hypothesis that the buried Lys-66 is in contact with internal waters even though these are not seen crystallographically. When only the one or two innermost water molecules were treated explicitly, epsilon(in) of 5-7 reproduced the pK(a) values. These values of epsilon(in) > 4 imply that some conformational reorganization occurs concomitant with the ionization of the buried groups.     

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a).     

Tissue factor alters the pK(a) values of catalytically important factor VIIa residues

Blood coagulation is triggered when the serine protease factor VIIa (fVIIa) binds to cell surface tissue factor (TF) to form the active enzyme-cofactor complex. TF binding to fVIIa allosterically augments the enzymatic activity of fVIIa toward macromolecular substrates and small peptidyl substrates. The mechanism of this enhancement remains unclear. Our previous studies have indicated that soluble TF (sTF; residues 1-219) alters the pH dependence of fVIIa amidolytic activity (Neuenschwander et al. (1993) Thromb. Haemostasis 70, 970), indicating an effect of TF on critical ionizations within the fVIIa active center. The pKa values and identities of these ionizable groups are unknown. To gain additional insight into this effect, we have performed a detailed study of the pH dependence of fVIIa amidolytic activity. Kinetic constants of Chromozym t-PA (MeSO(2)-D-Phe-Gly-Arg-pNA) hydrolysis at various pH values were determined for fVIIa alone and in complex with sTF. The pH dependence of both enzymes was adequately represented using a diprotic model. For fVIIa alone, two ionizations were observed in the free enzyme (pK(E1) = 7.46 and pK(E2) = 8.67), with at least a single ionization apparent in the Michaelis complex (pK(ES1) similar 7.62). For the fVIIa-TF complex, the pK(a) of one of the two important ionizations in the free enzyme was shifted to a more basic value (pK(E1) = 7.57 and pK(E2) = 9.27), and the ionization in the Michaelis complex was possibly shifted to a more acidic pH (pK(ES1) = 6.93). When these results are compared to those obtained for other well-studied serine proteases, K(E1) and K(ES1) are presumed to represent the ionization of the overall catalytic triad in the absence and presence of substrate, respectively, while K(E2) is presumed to represent ionization of the alpha-amino group of Ile(153). Taken together, these results would suggest that sTF binding to fVIIa alters the chemical environment of the fVIIa active site by protecting Ile(153) from deprotonation in the free enzyme while deprotecting the catalytic triad as a whole when in the Michaelis complex.     

Interaction of morphocycline with beryllium ions]

Behaviour of morphocycline (H5R) and its complex with beryllium ions in acqueous solutions was studied fluorimetrically. The ionization constants of H5R were estimated at pH 1.5-15 according to the data of fluorimetric determination with respect to OH-group: pK1 6.40, pK2 8.25, pK3 10.65, pK4 11.30. Two constants characterizing the deprotonization process with respect to the carbonylic group and nitrogen were also estimated: pK01--1.0 (greater than C = 0) and pK02 4.75 (--N=). Formation of an intensively fluorescing compound less than [Be3(OH)3(H2O2)5]2HR greater than 2+ was observed at pH 6.0-7.0. The cause of such intensive fluorescence was deformation of ion [Be3(OH)3(H2O)6]3+ because of its exclusion from the coordinating sphere of one molecule of water during the complex formation and decreasing of level H comes from II as compared to the morphocycline level II comes from n. A procedure for detecting morphocycline in the blood of humans and animals was developed.     

The pH-dependence of pepsin-catalysed reactions

1. The pH-dependence of the pepsin-catalysed hydrolysis of three peptide substrates was studied by using a method for the continuous monitoring of the formation of ninhydrin-positive products. 2. Two peptide acid substrates, N-acetyl-l-phenylalanyl-l-phenylalanine and N-acetyl-l-phenylalanyl-l-phenylalanyl-glycine, show apparent pK(a) values of 1.1 and 3.5 in the plots of k(0)/K(m) versus pH. By contrast a neutral substrate, N-acetyl-l-phenylalanyl-l-phenylalanine amide, shows apparent pK(a) values of 1.0 and 4.7. 3. Together with the data of the preceding paper (Knowles, Sharp & Greenwell, 1969), these results are taken to indicate that the rate of pepsin-catalysed hydrolysis is controlled by the ionization of two groups, which on the free enzyme have apparent pK(a) values of 1.0 and 4.7. It is apparent that the anions of peptide acid substrates are not perceptibly bound to the enzyme, resulting in apparent pK(a) values of 3.5 for the dependence of k(0)/K(m) for these materials.     

A possible method for characterizing the secondary structure of ribonucleic acids

The E(280)/E(260) ratio was found to be suitable for following the ionization of cytosine residues of polynucleotides on the basis of studies with model compounds such as oligoguanylic acid, oligocytidylic acid, a complex formed between polyadenylic acid and polyuridylic acid, and a copolymer of guanylic acid and cytidylic acid, provided that changes in secondary structure were taken into account. The pK of cytosine residues of a polynucleotide in the amorphous form was found to be 4.70 at 25 degrees in 0.1m-sodium phosphate on the basis of titration at 75-85 degrees and on the assumption that the heat of ionization was the same as the value (5.2kcal./mole) found for CMP. In contrast, the pK of cytosine residues in the double-helical form of DNA was found to be about 3.25. These observations were utilized in estimating the fraction of cytosine residues in helical segments of ribosomal RNA, a copolymer of guanylic acid and cytidylic acid, and a copolymer of adenylic acid, guanylic acid, uridylic acid and cytidylic acid. The ionization of guanine and uracil residues was estimated from changes in the E(270)/E(260) ratio and E(230)/E(260) ratio respectively. In the amorphous form of RNA both residues had the same pK, whereas in the double-helical form ionization was suppressed. The fraction of guanine and uracil residues in amorphous segments may be estimated from the titration curves. The difference in the denaturation spectrum of adenine--uracil and guanine--cytosine base pairs at 280mmu was enhanced in acidic solutions whereas E(260) was hardly affected. Hence a comparison of the increments in E(280) and E(260) obtained on increasing the temperature at constant pH may be used to distinguish the melting ranges of helical domains differing in nucleotide composition. In alkaline solutions comparison of the increments in E(260) and E(270) yields similar information. In acidic solutions the fraction of cytosine residues involved in helical secondary structure, the degree of ionization of cytosine residues and the fraction of adenine--uracil base pairs denatured may be estimated from DeltaE(265) and DeltaE(280). In alkaline solutions the fractions of guanine and uracil residues involved in secondary structure and the degrees of ionization of these residues may be estimated from DeltaE(230), DeltaE(245), DeltaE(260) and DeltaE(280).     

Quantitative structure toxicity relationships for catechols in isolated rat hepatocytes

One- and two-parameter quantitative structure toxicity relationship (QSTR) equations were obtained to describe the cytotoxicity of isolated rat hepatocytes induced by 23 catechols in which LD(50) represents the catechol concentration required to induce 50% cytotoxicity in 2 h. A QSTR equation logLD(50) (microM = - 0.464(+/-0.065) log P + 3.724(+/-0.114) (n = 20, r(2) = 0.740, s(y,x) = 0.372, P < 1 x 10(-6), outliers: 4-methoxycatechol, 3-methoxycatechol, L-dopa) was derived where logP represents octanol/water partitioning. Outliers were determined by adopting a statistical method to standardize the identification of outliers. When pK(a1), the first ionization constant, was considered as a contributing parameter a two-parameter QSTR equation was derived: logLD(50) (microM = - 0.343(+/-0.058) log P - 0.116(+/-0.041) pK(a1)+4.389 (+/-0.315) (n = 22, r(2) = 0.738, s(y,x) = 0.375, P < 0.01, outlier: 4-methoxycatechol). Replacing logP with logD(7.4), the partition coefficient at pH 7.4, improved the first correlation by limiting the outlier to 4-methoxycatechol: logLD(50) (microM)=-0.252(+/-0.039) logD(7.4)+3.168(+/-0.090) (n = 22, r(2) = 0.671, s(y,x) = 0.420, P < 1 x 10(-5). In this study, 4-methoxycatechol (readily autooxidizable) was found to be an outlier for all QSTR equations derived. These findings point to lipophilicity and pK(a1) as two important characteristics of catechols that can be used to predict their cytotoxicity towards isolated rat hepatocytes. The catechols with the higher lipophilicity/distribution coefficient, the lower degree of ionization and the higher pK(a(catechol)) were more toxic towards hepatocytes than the other catechols.     

MCCE analysis of the pKas of introduced buried acids and bases in staphylococcal nuclease

The pK(a)s of 96 acids and bases introduced into buried sites in the staphylococcal nuclease protein (SNase) were calculated using the multiconformation continuum electrostatics (MCCE) program and the results compared with experimental values. The pK(a)s are obtained by Monte Carlo sampling of coupled side chain protonation and position as a function of pH. The dependence of the results on the protein dielectric constant (ε(prot)) in the continuum electrostatics analysis and on the Lennard-Jones non-electrostatics parameters was evaluated. The pK(a)s of the introduced residues have a clear dependence on ε(prot,) whereas native ionizable residues do not. The native residues have electrostatic interactions with other residues in the protein favoring ionization, which are larger than the desolvation penalty favoring the neutral state. Increasing ε(prot) scales both terms, which for these residues leads to small changes in pK(a). The introduced residues have a larger desolvation penalty and negligible interactions with residues in the protein. For these residues, changing ε(prot) has a large influence on the calculated pK(a). An ε(prot) of 8-10 and a Lennard-Jones scaling of 0.25 is best here. The X-ray crystal structures of the mutated proteins are found to provide somewhat better results than calculations carried out on mutations made in silico. Initial relaxation of the in silico mutations by Gromacs and extensive side chain rotamer sampling within MCCE can significantly improve the match with experiment.     

Studies on the state of tyrosyl residues in a ribonuclease from seminal vesicles.

In order to study the state of tyrosyl residues in a ribouuclease from bovine semina vesicles [EC 3.1.4.22, RNase Vs1] several lines of experiments were carried out. Spectrophotometric titration of RNase Vs1 indicated that two out of 8 tyrosine residues were titrated very easily and their apparent pKa values were about 9.8. Next, about 4 residues were titrated at pH up to 13.5. The remaining 2 residues were titrated time-dependently at pH 13.5. In 8 M urea, about 6 tyrosine residues were titrated with apparent pK4 values of about 11.2 and about 2 residues were titrated time-dependently at pH 13.5. Acetylation of RNase Vs1 with N-acetylimidazole was studied at pH 7.5. In aqueous solution, about 1.1-3.5 tyrosine residues were acetylated, depending on the experimental conditions, and in 8 M urea, 5.3 tyrosine residues were modified. RNase Vs1 was nitrated with tetranitromethane at pH 7.5. In aqueous solution, about 2.5 tyrosine residues were nitrated very easily; the enzymatic activity of the modified enzymes was 130-200% of that of the native enzyme. In 8 M urea, the reactivity of the tyrosine residues increased and about 4-5.5 residues were modified. The results of chemical modification and spectrophotometric titration indicated that about two tyrosine residues in RNase Vs1 were exposed to the solvent and were more reactive to various reagents, and 3-4 tyrosine residues were less reactive. The final 2 residues were not accessible to the reagent even in the presence of urea, but were titraten at pH 13.5. The solvent perturbation difference spectrum using ethylene glycol as a perturbant indicated that about 4 tyrosine residues were perturbed. When the pH of the enzyme solution was changed from 7.0 to 1.0, the change in optical density of RNase Vs1 due to denaturation blue shift was about 1,600 at 287nm. The optical density change at 287 nm of native RNase Vs1 on exposure to 8 M urea and 6 M guanidine-HCl indicated that the environments of 2-3 and 4 tyrosine residues were changed by the addition of the denaturants, urea and guanidine-HCl, respectively. In RNase Vs1 having about four nitrotyrosine residues, the two most inaccessible tyrosine residues remained resistant to titration with alkali. On adding nucleotide, nitrated RNase Vs1 gave a difference spectrum in the ultraviolet region but not in 320-460 nm region, where nitrotyrosine residues absorb light. This may indicate that tyrosine residues located relatively near the surface of the molecule are not perturbed directly by nucleotide binding.     

Numerical method to analyze pH-profile from initial velocity data and semi-empirical procedure to identify catalytic residues of enzymes using pK and heat of ionization

A numerical computing method to estimate pK and heat of ionization directly from initial velocity data was developed for systematic analysis of rate parameters. Values of rate and thermodynamic parameters were determined together with standard deviations, without using Arrhenius and van't Hoff plots. Distributions of pK and heat of ionization of small-molecular derivatives of aliphatic carboxylates, imidazoles, thiols, and aliphatic amines as possible models for catalytic residues of enzymes were analyzed from a thermodynamic point of view. Each group was classified in a parallelogram without overlapping, except for a part of the imidazole and thiol groups. It was shown that values of pK and heat of ionization of carboxylates and histidine residues involved in the catalytic activity of enzymes determined from pH-profile experiments can be classified into similar, though somewhat extended, areas to those of the corresponding small-molecular derivatives. Identification of catalytic residues using values of pK and heat of ionization is proved to be a reliable method when the procedure is properly used.     

pK values of histidine residues in ribonuclease Sa: effect of salt and net charge

The primary goal of this study was to gain a better understanding of the effect of environment and ionic strength on the pK values of histidine residues in proteins. The salt-dependence of pK values for two histidine residues in ribonuclease Sa (RNase Sa) (pI=3.5) and a variant in which five acidic amino acids have been changed to lysine (5K) (pI=10.2) was measured and compared to pK values of model histidine-containing peptides. The pK of His53 is elevated by two pH units (pK=8.61) in RNase Sa and by nearly one pH unit (pK=7.39) in 5K at low salt relative to the pK of histidine in the model peptides (pK=6.6). The pK for His53 remains elevated in 1.5M NaCl (pK=7.89). The elevated pK for His53 is a result of screenable electrostatic interactions, particularly with Glu74, and a non-screenable hydrogen bond interaction with water. The pK of His85 in RNase Sa and 5K is slightly below the model pK at low salt and merges with this value at 1.5M NaCl. The pK of His85 reflects mainly effects of long-range Coulombic interactions that are screenable by salt. The tautomeric states of the neutral histidine residues are changed by charge reversal. The histidine pK values in RNase Sa are always higher than the pK values in the 5K variant. These results emphasize that the net charge of the protein influences the pK values of the histidine residues. Structure-based pK calculations capture the salt-dependence relatively well but are unable to predict absolute histidine pK values.     

Surface characteristics of phosphatidylglycerol phosphate from the extreme halophile Halobacterium cutirubrum compared with those of its deoxy analogue, at the air/water interface.

The relationship between area per molecule and surface pressure of monolayers of phosphatidylglycerol phosphate from extreme halophile Halobacterium cutrirubrum and its deoxy analogue, deoxyphosphatidylglycerol phosphate, spread at an air/water interface was examined. The effect of ionization of the primary and secondary acidic functions of the phosphate groups of the two lipids on surface characteristics of compression isotherms was determined by spreading monolayers on subphases with pH values ranging from below the apparent pKa of the primary ionization (pH 0) to greater than that of secondary ionization (pH 10.9). The limiting molecular area increases with decreasing pH below 2. Ionization of the primary phosphate functions of both phospholipids (with bulk pK1 values close to 4) is associated with a marked expansion of the films, as judged by values of limiting molecular area. Ionization of the secondary phosphate functions causes further expansion of the films, with the apparent pK2 of deoxyphosphatidylglycerol phosphate slightly less than that indicated for phosphatidylglycerol phosphate. Values of surface-compressibility modulus calculated from the surface characteristics of the phosphatidylglcerol phosphate monolayers showed that films spread on subphases with a pH of about the apparent pK1 of the primary phosphate functions were the least compressible. Increasing or decreasing subphase pH caused an increase in compressibility; this effect on compressibility was much less with monolayers of deoxyphosphatidylglycerol phosphate at high pH. The effect of inorganic counter-ions on monolayer characteristics of phosphatidylglycerol phosphate was examined by using subphases of NaCl concentrations varying from 0.01 to 1 M. The limiting molecular area was found to increase exponentially with respect to the subphase NaCl concentration.