Determination of the amino acid sequence requirements for catalysis by the highly proficient orotidine monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) catalyzes the decarboxylation of orotidine 5'-monophosphate to uridine 5'-monophosphate during pyrimidine nucleotide biosynthesis. This enzyme is one of the most proficient known, exhibiting a rate enhancement of over 17 orders of magnitude over the uncatalyzed rate. An interesting question is whether the high proficiency of ODCase is associated with a highly optimized sequence of active site residues. This question was addressed by randomizing 24 residue positions in and around the active site of the E. coli ODCase (pyrF) by site-directed mutagenesis. The libraries of mutants were selected for function from a multicopy plasmid or by single-copy replacement at the pyrF locus on the E. coli chromosome. Stringent sequence requirements for function were found for the mutants expressed from the chromosomal pyrF locus. Six positions were not tolerant of substitutions and several others accepted very limited substitutions. In contrast, all positions could be substituted to some extent when the library mutants were expressed from a multicopy plasmid. For the conserved quartet of charged residues Lys44-Asp71-Lys73-Asp76, a cysteine substitution was found to provide function at positions 71 and 76. A lower pK(a) for both cysteine mutants supports a mechanism whereby the thiolate group of cysteine substitutes for the negatively charged aspartate side chain. The partial function mutants such as D71C and D76C exhibit reduced catalytic efficiency relative to wild type but nevertheless provide a rate enhancement of 15 orders of magnitude over the uncatalyzed rate indicating the catalytic proficiency of the enzyme is robust and tolerant of mutation.     

Role of enzyme-ribofuranosyl contacts in the ground state and transition state for orotidine 5'-phosphate decarboxylase: a role for substrate destabilization?

The crystal structure of yeast orotidine 5'-monophosphate decarboxylase (ODCase) complexed with the inhibitor 6-hydroxyuridine 5'-phosphate (BMP) reveals the presence of a series of strong interactions between enzyme residues and functional groups of this ligand. Enzyme contacts with the phosphoribofuranosyl moiety of orotidine 5'-phosphate (OMP) have been shown to contribute at least 16.6 kcal/mol of intrinsic binding free energy to the stabilization of the transition state for the reaction catalyzed by yeast ODCase. In addition to these enzyme-ligand contacts, active site residues contributed by both subunits of the dimeric enzyme are positioned to form hydrogen bonds with the 2'- and 3'-OH groups of the ligand's ribosyl moiety. These involve Thr-100 of one subunit and Asp-37 of the opposite subunit, respectively. To evaluate the contributions of these ribofuranosyl contacts to ground state and transition state stabilization, Thr-100 and Asp-37 were each mutated to alanine. Elimination of the enzyme's capacity to contact individual ribosyl OH groups reduced the k(cat)/K(m) value of the T100A enzyme by 60-fold and that of the D37A enzyme by 300-fold. Removal of the 2'-OH group from the substrate OMP decreased the binding affinity by less than a factor of 10, but decreased k(cat) by more that 2 orders of magnitude. Upon removal of the complementary hydroxymethyl group from the enzyme, little further reduction in k(cat)/K(m) for 2'-deoxyOMP was observed. To assess the contribution made by contacts involving both ribosyl hydroxyl groups at once, the ability of the D37A mutant enzyme to decarboxylate 2'-deoxyOMP was measured. The value of k(cat)/K(m) for this enzyme-substrate pair was 170 M(-1) s(-1), representing a decrease of more than 7.6 kcal/mol of binding free energy in the transition state. To the extent that electrostatic repulsion in the ground state can be tested by these simple alterations, the results do not lend obvious support to the view that electrostatic destabilization in the ground state enzyme-substrate complex plays a major role in catalysis.     

Continuum electrostatic model for the binding of cytochrome c2 to the photosynthetic reaction center from Rhodobacter sphaeroides

Electrostatic interactions are important for protein-protein association. In this study, we examined the electrostatic interactions between two proteins, cytochrome c(2) (cyt c(2)) and the reaction center (RC) from the photosynthetic bacterium Rhodobacter sphaeroides, that function in intermolecular electron transfer in photosynthesis. Electrostatic contributions to the binding energy for the cyt c(2)-RC complex were calculated using continuum electrostatic methods based on the recent cocrystal structure [Axelrod, H. L., et al. (2002) J. Mol. Biol. 319, 501-515]. Calculated changes in binding energy due to mutations of charged interface residues agreed with experimental results for a protein dielectric constant epsilon(in) of 10. However, the electrostatic contribution to the binding energy for the complex was close to zero due to unfavorable desolvation energies that compensate for the favorable Coulomb attraction. The electrostatic energy calculated as a function of displacement of the cyt c(2) from the bound position showed a shallow minimum at a position near but displaced from the cocrystal configuration. These results show that although electrostatic steering is present, other short-range interactions must be present to contribute to the binding energy and to determine the structure of the complex. Calculations made to model the experimental data on association rates indicate a solvent-separated transition state for binding in which the cyt c(2) is displaced approximately 8 A above its position in the bound complex. These results are consistent with a two-step model for protein association: electrostatic docking of the cyt c(2) followed by desolvation to form short-range van der Waals contacts for rapid electron transfer.     

Role of the electrostatic loop charged residues in Cu,Zn superoxide dismutase.

We have expressed and characterized a mutant of Xenopus laevis Cu,Zn superoxide dismutase in which four highly conserved charged residues belonging to the electrostatic loop have been replaced by neutral side chains: Lys120 --> Leu, Asp130 --> Gln, Glu131 --> Gln, and Lys134 --> Thr. At low ionic strength, the mutant enzyme is one of the fastest superoxide dismutases ever assayed (k = 6.7 x 10(9) M(-1) s(-1), at pH 7 and mu = 0.02 M). Brownian dynamics simulations give rise to identical enzyme-substrate association rates for both wild-type and mutant enzymes, ruling out the possibility that enhancement of the activity is due to pure electrostatic factors. Comparative analysis of the experimental catalytic rate of the quadruple and single mutants reveals the nonadditivity of the mutation effects, indicating that the hyperefficiency of the mutant is due to a decrease of the energy barrier and/or to an alternative pathway for the diffusion of superoxide within the active site channel. At physiological ionic strength the catalytic rate of the mutant at neutral pH is similar to that of the wild-type enzyme as it is to the catalytic rate pH dependence. Moreover, mutation effects are additive. These results show that, at physiological salt conditions, electrostatic loop charged residues do not influence the diffusion pathway of the substrate and, if concomitantly neutralized, are not essential for high catalytic efficiency of the enzyme, pointing out the role of the metal cluster and of the invariant Arg141 in determining the local electrostatic forces facilitating the diffusion of the substrate towards the active site.     

Computational and mutational analysis of human glutaredoxin (thioltransferase): probing the molecular basis of the low pKa of cysteine 22 and its role in catalysis

Human glutaredoxin (GRx), also known as thioltransferase, is a 12 kDa thiol-disulfide oxidoreductase that is highly selective for reduction of glutathione-containing mixed disulfides. The apparent pK(a) for the active site Cys22 residue is approximately 3.5. Previously we observed that the catalytic enhancement by glutaredoxin could be ascribed fully to the difference between the pK(a) of its Cys22 thiol moiety and the pK(a) of the product thiol, each acting as a leaving group in the enzymatic and nonenzymatic reactions, respectively [Srinivasan et al. (1997), Biochemistry 36, 3199-3206]. Continuum electrostatic calculations suggest that the low pK(a) of Cys22 results primarily from stabilization of the thiolate anion by a specific ion-pairing with the positively charged Lys19 residue, although hydrogen bonding interactions with Thr21 also appear to contribute. Variants of Lys19 were considered to further assess the predicted role of Lys19 on the pK(a) of Cys22. The variants K19Q and K19L were generated by molecular modeling, and the pK(a) value for Cys22 was calculated for each variant. For K19Q, the predicted Cys22 pK(a) is 7.3, while the predicted value is 8.3 for K19L. The effects of the mutations on the interaction energy between the adducted glutathionyl moiety and GRx were roughly estimated from the van der Waals and electrostatic energies between the glutathionyl moiety and proximal protein residues in a mixed disulfide adduct of GRx and glutathione, i.e., the GRx-SSG intermediate. The values for the K19 mutants differed by only a small amount compared to those for the wild type enzyme intermediate. Together, the computational analysis predicted that the mutant enzymes would have markedly reduced catalytic rates while retaining the glutathionyl specificity displayed by the wild type enzyme. Accordingly, we constructed and characterized the K19L and K19Q mutants of two forms of the GRx enzyme. Each of the mutants retained glutathionyl specificity as predicted and displayed diminution in activity, but the decreases in activity were not to the extent predicted by the theoretical calculations. Changes in the respective Cys22-thiol pK(a) values of the mutant enzymes, as shown by pH profiles for iodoacetamide inactivation of the respective enzymes, clearly revealed that the K19-C22 ion pair cannot fully account for the low pK(a) of the Cys22 thiol. Additional contributions to stabilization of the Cys22 thiolate are likely donated by Thr21 and the N-terminal partial positive charge of the neighboring alpha-helix.     

Interpreting trends in the binding of cyclic ureas to HIV-1 protease

The design of new HIV protease inhibitors requires an improved understanding of the physical basis of inhibitor/protein binding. Here, the binding affinities of seven aliphatic cyclic ureas to HIV-1 protease are calculated using a predominant states method and an implicit solvent model based upon finite difference solutions of the Poisson-Boltzmann equation. The calculations are able to reproduce the observed U-shaped trend of binding free energy as a function of aliphatic chain length. Interestingly, the decrease in affinity for the longest chains is attributable primarily to the energy cost of partly desolvating charged aspartic and arginine groups at the mouths of the active site. Even aliphatic chains too short to contact these charged groups directly are subject to considerable desolvation penalties. We are not aware of other systems where binding affinity trends have been attributed to long-ranged electrostatic desolvation of ionized groups. A generalized Born/surface area solvation model yields a much smaller change in desolvation energy with chain length and, therefore, does not reproduce the experimental binding affinity trends. This result suggests that the generalized Born model should be used with caution for complex, partly desolvated systems like protein binding sites. We also find that changing the assumed protonation state of the active site aspartyl dyad significantly affects the computed binding affinity trends. The protonation state of the aspartyl dyad in the presence of cyclic ureas is discussed in light of the observation that the monoprotonated state reproduces the experimental results best.     

Role of charged residues in the catalytic mechanism of hepatitis C virus NS3 protease: electrostatic precollision guidance and transition-state stabilization

Maturational cleavage of the hepatitis C virus polyprotein involves the viral chymotrypsin-like serine protease NS3. The substrate binding site of this enzyme is unusually flat and featureless. We here show that NS3 has a highly asymmetric charge distribution that is characterized by strong positive potentials in the vicinity of its active site and in the S5/S6 region. Using electrostatic potential calculations, we identified determinants of this positive potential, and the role of six different residues was explored by site-directed mutagenesis. Mutation of residues in the vicinity of the active site led to changes in k(cat) values of a peptide substrate indicating that basic amino acids play a role in the stabilization of the transition state. Charge neutralization in the S5/S6 region increased the K(m) values of peptide substrates in a manner that depended on the presence of negatively charged residues in the P5 and P6 positions. K(i) values of hexapeptide acids spanning P6-P1 (product inhibitors) were affected by charge neutralization in both the active site region and the S5/S6 region. Pre-steady-state kinetic data showed that the electrostatic surface potential is used by this enzyme to enhance collision rates between peptidic ligands and the active site. Calculations of the interaction energies of protease-substrate or protease-inhibitor complexes showed that electrostatic interaction energies oppose the formation of a tightly bound complex due to an unfavorable change in the desolvation energy. We propose that desolvation costs are minimized by avoiding the formation of individual ion pair interactions through the use of clusters of positively charged residues in the generation of local electrostatic potentials.     

Atomic Resolution Structure of the Orotidine 5′-Monophosphate Decarboxylase Product Complex Combined with Surface Plasmon Resonance Analysis: IMPLICATIONS FOR THE CATALYTIC MECHANISM*

Orotidine 5′-monophosphate decarboxylase (ODCase) accelerates the decarboxylation of its substrate by 17 orders of magnitude. One argument brought forward against steric/electrostatic repulsion causing substrate distortion at the carboxylate substituent as part of the catalysis has been the weak binding affinity of the decarboxylated product (UMP). The crystal structure of the UMP complex of ODCase at atomic resolution (1.03 Å) shows steric competition between the product UMP and the side chain of a catalytic lysine residue. Surface plasmon resonance analysis indicates that UMP binds 5 orders of magnitude more tightly to a mutant in which the interfering side chain has been removed than to wild-type ODCase. These results explain the low affinity of UMP and counter a seemingly very strong argument against a contribution of substrate distortion to the catalytic reaction mechanism of ODCase.     

Second step of hydrolytic dehalogenation in haloalkane dehalogenase investigated by QM/MM methods

Mechanistic studies on the hydrolytic dehalogenation catalyzed by haloalkane dehalogenases are of importance for environmental and industrial applications. Here, Car-Parrinello (CP) and ONIOM hybrid quantum-mechanical/molecular mechanics (QM/MM) are used investigate the second reaction step of the catalytic cycle, which comprises a general base-catalyzed hydrolysis of an ester intermediate (EI) to alcohol and free enzyme. We focus on the enzyme LinB from Sphingomonas paucimobilis UT26, for which the X-ray structure at atomic resolution is available. In agreement with previous proposals, our calculations suggest that a histidine residue (His272), polarized by glutamate (Glu132), acts as a base, accepting a proton from the catalytic water molecule and transferring it to an alcoholate ion. The reaction proceeds through a metastable tetrahedral intermediate, which shows an easily reversed reaction to the EI. In the formation of the products, the protonated aspartic acid (Asp108) can easily adopt conformation of the relaxed state found in the free enzyme. The overall free energy barrier of the reaction calculated by potential of the mean force integration using CP-QM/MM calculations is equal to 19.5 +/- 2 kcal . mol(-1). The lowering of the energy barrier of catalyzed reaction with respect to the water reaction is caused by strong stabilization of the reaction intermediate and transition state and their preorganization by electrostatic field of the enzyme.     

Structural basis for the catalytic mechanism of a proficient enzyme: orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) catalyzes the decarboxylation of orotidine 5'-monophosphate, the last step in the de novo synthesis of uridine 5'-monophosphate. ODCase is a very proficient enzyme [Radzicka, A., and Wolfenden, R. (1995) Science 267, 90-93], enhancing the reaction rate by a factor of 10(17). This proficiency has been enigmatic, since it is achieved without metal ions or cofactors. Here we present a 2.5 A resolution structure of ODCase complexed with the inhibitor 1-(5'-phospho-beta-D-ribofuranosyl)barbituric acid. It shows a closely packed dimer composed of two alpha/beta-barrels with two shared active sites. The orientation of the orotate moiety of the substrate is unambiguously deduced from the structure, and previously proposed catalytic mechanisms involving protonation of O2 or O4 can be ruled out. The proximity of the OMP carboxylate group with Asp71 appears to be instrumental for the decarboxylation of OMP, either through charge repulsion or through the formation of a very short O.H.O hydrogen bond between the two carboxylate groups.     

Binding of small basic peptides to membranes containing acidic lipids: theoretical models and experimental results.

We measured directly the binding of Lys3, Lys5, and Lys7 to vesicles containing acidic phospholipids. When the vesicles contain 33% acidic lipids and the aqueous solution contains 100 mM monovalent salt, the standard Gibbs free energy for the binding of these peptides is 3, 5, and 7 kcal/mol, respectively. The binding energies decrease as the mol% of acidic lipids in the membrane decreases and/or as the salt concentration increases. Several lines of evidence suggest that these hydrophilic peptides do not penetrate the polar headgroup region of the membrane and that the binding is mainly due to electrostatic interactions. To calculate the binding energies from classical electrostatics, we applied the nonlinear Poisson-Boltzmann equation to atomic models of the phospholipid bilayers and the basic peptides in aqueous solution. The electrostatic free energy of interaction, which arises from both a long-range coulombic attraction between the positively charged peptide and the negatively charged lipid bilayer, and a short-range Born or image charge repulsion, is a minimum when approximately 2.5 A (i.e., one layer of water) exists between the van der Waals surfaces of the peptide and the lipid bilayer. The calculated molar association constants, K, agree well with the measured values: K is typically about 10-fold smaller than the experimental value (i.e., a difference of about 1.5 kcal/mol in the free energy of binding). The predicted dependence of K (or the binding free energies) on the ionic strength of the solution, the mol% of acidic lipids in the membrane, and the number of basic residues in the peptide agree very well with the experimental measurements. These calculations are relevant to the membrane binding of a number of important proteins that contain clusters of basic residues.     

Expression and sequence analysis of a cDNA encoding the orotidine-5'-monophosphate decarboxylase domain from Ehrlich ascites uridylate synthase

Orotidine-5'-monophosphate decarboxylase (OD-Case) catalyzes the conversion of orotidine 5'-monophosphate to UMP. In mammals, ODCase is present as part of a bifunctional protein which also contains orotate phosphoribosyltransferase; the preceding enzyme in the de novo UMP biosynthetic pathway. We have isolated a plasmid (pMEJ) which contains a cDNA for the ODCase domain of UMP synthase. Insertion of this sequence into an Escherichia coli expression vector (pUC12) has allowed for the expression of ODCase and not orotate phosphoribosyltransferase in E. coli. The molecular weight of the expressed protein is 26,000-27,300 from immunoblot analysis which corresponds closely to the molecular weight of the ODCase domain (28,500) isolated by tryptic digestion of UMP synthase. We have sequenced the cDNA insert of pMEJ and deduced the amino acid sequence. The molecular weight of the ODCase domain calculated from the amino acid sequence in 28,654. Comparison of the deduced amino acid sequence from pMEJ with that for yeast ODCase (a monofunctional protein) demonstrated that 52% of the amino acids were identical when the two sequences are compared. Furthermore, several stretches of the amino acid sequence have 80% or greater absolute homology.     

A unique catalytic and inhibitor-binding role for Lys93 of yeast orotidylate decarboxylase.

The presence of a proton-donating catalytic amino acid side chain in orotidylate decarboxylase (ODCase) was sought by site-directed mutagenesis. Replacement of yeast ODCase Lys93 with a cysteine resulted in a mutant protein (K93C) with no measurable activity, representing a decrease in activity by a factor of, at most, 2 x 10(-8) times the activity of the wild-type enzyme. Treatment of this mutant protein with 2-bromoethylamine, designed to append Cys93 to yield S-(2-aminoethyl)cysteine, restored activity by a factor of at least 5 x 10(5) over the untreated mutant protein. Activity could not be restored by treatment with other brominated reagents designed to replace the epsilon-amino of S-(2-aminoethyl)Cys93 with a different functional group. The overall architecture of the K93C protein was not significantly changed, as judged by the similar dimerization properties (in the absence of ligands) of the mutant enzyme compared to the wild-type enzyme. The binding affinity of the substrate orotidylate was not measurably changed by the mutation, indicating that Lys93 has an essential role in catalysis which is mechanistically distinguishable from substrate binding. Apparently the mutation removes an integral portion of the active site and does not drastically affect the structural or substrate binding properties. However, the affinities of the mutant protein for the competitive inhibitors 6-azauridylate (6-azaUMP) and UMP are significantly altered from the pattern seen with the wild-type enzyme. The K93C protein has an affinity for the neutral ligand UMP which is greater than that for the anionic 6-azaUMP, in clear contrast to the preference for 6-azaUMP displayed by the wild-type enzyme. Lys93 is apparently critical for catalysis of the substrate to product and for the binding of anionic inhibitors; the data are discussed in terms of previously existing models for transition-state analogue inhibitor binding and catalysis.     

Study of the kinetic and physical properties of the orotidine-5'-monophosphate decarboxylase domain from mouse UMP synthase produced in Saccharomyces cerevisiae

In mammals, the bifunctional protein UMP synthase contains the final two enzymatic activities, orotate phosphoribosyltransferase and orotidine-5'-monophosphate decarboxylase (ODCase), for de novo biosynthesis of UMP. The plasmid pMEJ contains a cDNA for the ODCase domain of mouse Ehrlich ascites UMP synthase. The cDNA from pMEJ was joined to the Saccharomyces cerevisiae iso-1-cytochrome c (CYC1) promoter and the first four CYC1 coding nucleotides in the plasmid pODCcyc. ODCase-deficient yeast cells (HF200x1) transformed with pODCcyc expressed an active ODCase domain with a specific activity of 20 nmol/min/mg in cell extracts. The expressed ODCase domain has a lower affinity for the substrate orotidine 5'-monophosphate and the inhibitor 6-azauridine 5'-monophosphate than intact UMP synthase or an ODCase domain isolated after proteolysis of homogenous UMP synthase. Sucrose density gradient sedimentation experiments showed that the expressed ODCase domain forms a dimer in the presence of ligands which bind at the catalytic site. These studies support the existence of an ODCase structural domain which contains the ODCase catalytic site and a dimerization surface of UMP synthase, but the domain may not have the regulatory site required to form the altered dimer form.     

Impact of ligand and protein desolvation on ligand binding to the S1 pocket of thrombin

In the present study, we investigate the impact of a tightly bound water molecule on ligand binding in the S1 pocket of thrombin. The S1 pocket contains a deeply buried deprotonated aspartate residue (Asp189) that is, due to its charged state, well hydrated in the uncomplexed state. We systematically studied the importance of this water molecule by evaluating a series of ligands that contains pyridine-type P1 side chains that could potentially alter the binding properties of this water molecule. All of the pyridine derivatives retain the original hydration state albeit sometimes with a slight perturbance. In order to prevent a direct H-bond formation with Asp189, and to create a permanent positive charge on the P1 side chain that is positioned adjacent to the Asp189 carboxylate anion, we methylated the pyridine nitrogen. This methylation resulted in displacement of water but was accompanied by a loss in binding affinity. Quantum chemical calculations of the ligand solvation free energy showed that the positively charged methylpyridinium derivatives suffer a large penalty of desolvation upon binding. Consequently, they have a substantially less favorable enthalpy of binding. In addition to the ligand desolvation penalty, the hydration shell around Asp189 has to be overcome, which is achieved in nearly all pyridinium derivatives. Only for the ortho derivative is a partial population of a water next to Asp189 found. Possibly, the gain of electrostatic interactions between the charged P1 side chain and Asp189 helps to compensate for the desolvation penalty. In all uncharged pyridine derivatives, the solvation shell remains next to Asp189, partly mediating interactions between ligand and protein. In the case of the para-pyridine derivative, a strongly disordered cluster of water sites is observed between ligand and Asp189.     

Molecular mechanics PBSA ligand binding energy and interaction of Efavirenz derivatives with HIV-1 reverse transcriptase

In order to evaluate the properties of several HIV-1 reverse transcripase(RT) inhibitors, Efavirenz (SUSTIVA) and a set of its derivatives (benzoxazinones) have been placed into the nonnucleoside analogue binding site of the enzyme by molecular docking. The resulting geometries were used for a molecular dynamics simulation and binding energy calculations. The enzyme-inhibitor binding energies were estimated from experimental inhibitory activities (IC90). The correlation of the predicted and experimental binding energies were satisfactory acceptable as indicated by r2 = 0.865. Based on MD simulations, the obtained results indicate that the tight association of the ligand to the HIV-1 RT binding pocket was based on hydrogen bonding between Efavirenz's N1 and the oxygen of the backbone of Lys 101, with an estimated average distance of 1.88 A. Moreover, electrostatic interaction was mainly contributed by two amino acid residues in the binding site; Lys 101 and His 235. MD simulations open the possibility to study the reaction of the flexible enzyme to those substances as well as the overall affinity.     

Effect of lysine side chain length on intra-helical glutamate--lysine ion pairing interactions

Ion-pairing interactions are important for protein stabilization. Despite the apparent electrostatic nature of these interactions, natural positively charged amino acids Lys and Arg have multiple methylenes linking the charged functionality to the backbone. Interestingly, the amino acids Lys and Orn have positively charged side chains that differ by only one methylene. However, only Lys is encoded and incorporated into proteins. To investigate the effect of side chain length of Lys on ion-pairing interactions, a series of 12 monomeric alpha-helical peptides containing potential Glu-Xaa (i, i+3), (i, i+4) and (i, i+5) (Xaa = Lys, Orn, Dab, Dap) interactions were studied by circular dichroism (CD) spectroscopy at pH 7 and 2. At pH 7, no Glu-Xaa (i, i+5) interaction was observed, regardless of the Xaa side chain length. Furthermore, only Lys was capable of supporting Glu-Xaa (i, i+3) interactions, whereas any Xaa side chain length supported Glu-Xaa (i, i+4) interactions. Side chain conformational analysis by molecular mechanics calculations showed that the side chain length of Lys enables the Glu-Xaa (i, i+3) interaction with lower energy conformations compared to residues with side chain lengths shorter than that of Lys. Furthermore, these calculated low energy conformers were consistent with conformations of intra-helical Glu-Lys salt bridges in a non-redundant protein structure database. Importantly, the CD spectra for peptides with Glu-Lys interactions did not alter significantly upon changing the pH because of a greater contribution to these interactions by forces other than electrostatics. Incorporating side chains just one methylene shorter (Orn) resulted in significant pH dependence or lack of interaction, suggesting that nature has chosen Lys to form durable interactions with negatively charged functional groups.     

Density functional models of the mechanism for decarboxylation in orotidine decarboxylase

The mechanism of orotidine 5-monophosphate decarboxylase (ODCase) has been modeled using hybrid Density Functional Theory (B3LYP functional). The main goal of the present study was to investigate if much larger quantum chemical models of the active site than previously used could shed new light on the mechanism. The models used include the five conserved amino acids expected to be the most important ones for catalysis. One result of this model is that a mechanism involving a direct cleavage of the C–C bond followed by a protonation of C6 by Lys93 appears unlikely, with a barrier for decarboxylation 20 kcal mol^–1 too high. Additional effects like electrostatic stress and ground-state destabilization have been estimated to have only a minor influence on the reaction barrier. The conclusion from the calculations is that the negative charge developing on the substrate during decarboxylation must be stabilized by a protonation of the carbonyl O2 of the substrate. For this mechanism, the addition of the catalytic amino acids decreases the reaction barrier by 25 kcal mol^–1, but full agreement with experimental results has still not been reached. Further modifications of this mechanism are discussed. Electronic supplementary material to this paper can be obtained by using the Springer LINK server located a http://dx.doi.org/10.1007/s00894-002-0080-2.Electronic Supplementary Material available.     

Transition state stabilization by the N-terminal anticodon-binding domain of lysyl-tRNA synthetase

Lysyl-tRNA synthetase from Bacillus stearothermophilus (B.s. LysRS) (EC ) catalyzes aminoacylation of tRNA(Lys) with l-lysine, in which l-lysine was first activated with ATP to yield an enzyme (lysyladenylate complex), and then the lysine molecule was transferred from the complex to tRNA(Lys). B.s. LysRS is a homodimeric enzyme with a subunit that consists of two domains, an N-terminal tRNA anticodon-binding domain (TAB-ND: Ser(1)-Pro(144)) and a C-terminal Class II-specific catalytic domain (CAT-CD: Lys(151)-Lys(493)). CAT-CD alone retained catalytic activity, although at a low level; TAB-ND alone showed no activity. Size exclusion chromatography revealed that CAT-CD exists as a dimer, whereas TAB-ND was a monomer. The formation of a complex consisting of these domains was detected with the guidance of surface plasmon resonance. In accordance with this, the addition of TAB-ND to CAT-CD significantly enhanced both the l-lysine activation and the tRNA aminoacylation reactions. Kinetic analysis showed that deletion of TAB-ND resulted in a significant destabilization of the transition state of CAT-CD in the l-lysine activation reaction but had little effect on the ground state of substrate binding. A significant role of a cross-subunit interaction in the enzyme between TAB-ND and CAT-CD was proposed for the stabilization of the transition state in the l-lysine activation reaction.