Design and synthesis of new enzymes based on the lactate dehydrogenase framework

Phage display identification of age-associated TNFα-mediated cardiac oxidative induction

Age-associated alterations in the actions of tumor necrosis factor-α (TNFα) in the heart with impaired cardioprotective pathways and enhanced apoptotic induction may contribute to the increased severity of cardiovascular pathology in older persons. To identify the molecular events mediating these changes in the microvasculature of the aging rodent heart, the biochemical properties of in vivo phage-display cyclic peptide cardiac biopanning were studied. Analysis of individual amino acid positions revealed that the center of the peptide motif (amino acid position 4) had a significantly higher frequency of aromatic amino acid side chains in phage homing to the old hearts compared with young controls (18 mo old, 11% vs. 3 mo old, 3%, P < 0.05). This subset of phage motifs revealed an age-associated homology with oxidoreductase enzymes (homology: 18 mo, 7/7; 3 mo, 0/2), suggesting the substrates and/or binding sites of these enzymes are increased in the aging hearts. Immunostaining for the oxidoreductase substrate 4-hydroxy-2-nonenal (HNE), a cardiotoxic lipid peroxidation product, demonstrated a twofold higher density of HNE(+) cells in PBS-treated hearts of old mice (18 mo) compared with young controls (3 mo) (18 mo, 3.2 ± 2.8 vs. 3 mo, 1.0 ± 0.9 cells/HPF, P < 0.05). Moreover, intracardiac injection of TNFα resulted in a significantly greater increase in HNE staining in the old hearts (18 mo, 16.9 ± 13.8 vs. 3 mo, 9.1 ± 6.0 cells/HPF, P < 0.05). Overall, these studies demonstrate that aging-associated alterations in TNFα-mediated pathways with induction of reactive oxidative species and changes in vascular surface binding sites may contribute mechanistically to the increased cardiovascular pathology of the aging heart.

Identification of Babesia bovis l-lactate dehydrogenase as a potential chemotherapeutical target against bovine babesiosis

In this study, we characterized a novel Babesia bovis cDNA clone obtained by immunoscreening the cDNA expression phage library with B. bovis-infected bovine serum. The genetic analyses showed that it contained an open reading frame of 993 bp, which was considered to encode B. bovis L-lactate dehydrogenase (BbLDH: E.C. 1.1.1.27) because of the strikingly high amino acid identities of its gene product to the LDHs of Plasmodium falciparum and Toxoplasma gondii. Immunological analyses with the anti-recombinant BbLDH mouse serum showed that 36 kDa of the native BbLDH was expressed not only in the cytoplasm of intra- and extraerythrocytic parasites but also along the membrane of infected erythrocytes. The kinetic properties of recombinant BbLDH proved a certain enzymatic activity of LDH, and the activity was significantly inhibited by the addition of gossypol, a competitive inhibitor of protozoan LDHs. Moreover, 100 microM of the gossypol irretrievably arrested the in vitro growth of B. bovis. The results demonstrated that BbLDH provides a suitable drug target for the design of novel babesial chemotherapeutics.

Stochastic roadmap simulation: an efficient representation and algorithm for analyzing molecular motion

Classic molecular motion simulation techniques, such as Monte Carlo (MC) simulation, generate motion pathways one at a time and spend most of their time in the local minima of the energy landscape defined over a molecular conformation space. Their high computational cost prevents them from being used to compute ensemble properties (properties requiring the analysis of many pathways). This paper introduces stochastic roadmap simulation (SRS) as a new computational approach for exploring the kinetics of molecular motion by simultaneously examining multiple pathways. These pathways are compactly encoded in a graph, which is constructed by sampling a molecular conformation space at random. This computation, which does not trace any particular pathway explicitly, circumvents the local-minima problem. Each edge in the graph represents a potential transition of the molecule and is associated with a probability indicating the likelihood of this transition. By viewing the graph as a Markov chain, ensemble properties can be efficiently computed over the entire molecular energy landscape. Furthermore, SRS converges to the same distribution as MC simulation. SRS is applied to two biological problems: computing the probability of folding, an important order parameter that measures the "kinetic distance" of a protein's conformation from its native state; and estimating the expected time to escape from a ligand-protein binding site. Comparison with MC simulations on protein folding shows that SRS produces arguably more accurate results, while reducing computation time by several orders of magnitude. Computational studies on ligand-protein binding also demonstrate SRS as a promising approach to study ligand-protein interactions.

Protein adaptations to temperature and pressure: complementary roles of adaptive changes in amino acid sequence and internal milieu

Retention of required structural and functional properties of proteins in species adapted to different temperatures and pressures is achieved through variation in amino acid sequence and accumulation of small organic solutes that stabilize protein traits. Conservation of ligand binding and catalytic rate can be achieved by minor differences in sequence. For orthologs of lactate dehydrogenase-A (A(4)-LDH) temperature adaptation may involve only a single amino acid substitution. Adaptation involves changes in conformational mobility of regions of A(4)-LDH that undergo movement during ligand binding, movements that are rate-limiting to catalysis. A model that integrates adaptations in sequence and intracellular milieu is developed on the basis of conformational microstates. Although orthologs of different thermally adapted species vary in stability, at physiological temperatures it is hypothesized that a similar ensemble of conformational microstates exists for all orthologs. Organic solutes stabilize this ensemble of microstates. Differences among orthologs in responses to organic solutes at a common temperature lead to similar responses at normal body temperatures. Because protein stability increases at high protein concentrations, intrinsic stabilities of proteins may reflect the protein concentrations of the cellular compartments in which they occur. Protein-stabilizing solutes like trimethylamine-N-oxide (TMAO) conserve protein function and structure at elevated hydrostatic pressures.

Computer simulations aimed at structure prediction of supersecondary motifs in proteins

It is well established that protein structures are more conserved than protein sequences. One-third of all known protein structures can be classified into ten protein folds, which themselves are composed mainly of alpha-helical hairpin, beta hairpin, and betaalphabeta supersecondary structural elements. In this study, we explore the ability of a recent Monte Carlo-based procedure to generate the 3D structures of eight polypeptides that correspond to units of supersecondary structure and three-stranded antiparallel beta sheet. Starting from extended or misfolded compact conformations, all Monte Carlo simulations show significant success in predicting the native topology using a simplified chain representation and an energy model optimized on other structures. Preliminary results on model peptides from nucleotide binding proteins suggest that this simple protein folding model can help clarify the relation between sequence and topology.

Intrinsic versus extrinsic stabilization of enzymes: the interaction of solutes and temperature on A4-lactate dehydrogenase orthologs from warm-adapted and cold-adapted marine fishes

We examined the effects of temperature and stabilizing solutes on A4-lactate dehydrogenase (A4-LDH) from warm- and cold-adapted fishes, to determine how extrinsic stabilizers affect orthologs with different intrinsic stabilities. Conformational changes during substrate binding are rate-limiting for A4-LDH, thus stabilization due to intrinsic or extrinsic factors leads to decreased activity. A4-LDH from a warm-temperate goby (Gillichthys mirabilis), which has lower values for kcat and the Michaelis constant for pyruvate ( K m PYR), was intrinsically more stable than the orthologs of the cold-adapted Antarctic notothenioids Parachaenichthys charcoti and Chionodraco rastrospinosus, as shown by a higher apparent transition ('melting') temperature (Tm(APP)). We used four solutes, glycerol, sucrose, trimethylamine-N-oxide and poly(ethylene glycol) 8000, which stabilize proteins through different modes of preferential exclusion, to study temperature-solute interactions of the three orthologs. Changes in Tm(APP) were similar for all orthologs in each solute tested, but the catalytic rate of G. mirabilis A4-LDH was decreased most by solutes and increased most by temperature. In contrast, the K m PYR values of the Antarctic orthologs were more affected than that of the goby by both solutes and temperature. We conclude that (a) preferential exclusion of solutes functions within the native state of A4-LDH to favor conformational microstates with minimal surface area; (b) the varied effects of the different solutes on the kinetic properties are due to the interaction between this nonspecific stabilization and the differing intrinsic stabilities of the orthologs; (c) the catalytic rates of A4-LDH orthologs are equally affected by stabilizing solutes, if measurements are made at physiologically appropriate temperatures; and (d) global stability and localized flexibility of these A4-LDH orthologs may evolve independently.

Selective active site inhibitors of human lactate dehydrogenases A4, B4, and C4

Human lactate dehydrogenases (LDH-A4, -B4, and -C4) are highly homologous with 84-89% sequence similarities and 69-75% amino acid identities. Active site residues are especially conserved. Gossypol, a natural product from cotton seed, is a non-selective competitive inhibitor of NADH binding to LDH, with K(i) values of 1.9, 1.4, and 4.2 microM for LDH-A4, -B4, and -C4, respectively. However, derivatives of gossypol and structural analogs of gossypol in the substituted 2,3-dihydroxy-1-naphthoic acid family exhibited markedly greater selectivity and, in many cases, greater potency. For gossypol derivatives, greater than 35-fold selectivity was observed. For dihydroxynaphthoic acids with substituents at the 4- and 7-positions, greater than 200-fold selectivity was observed. Inhibition was consistently competitive with the binding of NADH, with dissociation constants as low as 30 nM. By comparison, a series of N-substituted oxamic acids, which are competitive inhibitors of the binding of pyruvate to LDH, exhibited very modest selectivity. These results suggest that substituted dihydroxynaphthoic acids are good lead compounds for the development of selective LDH inhibitors. Selective inhibitors of LDH-C4 targeted to the dinucleotide fold may hold promise as male antifertility drugs. Selective inhibitors of LDH-A4 and -B4 may be useful for studies of lactic acidemia associated with ischemic events. More broadly, the results raise the question of the general utility of drug design targeted at the dinucleotide binding sites of dehydrogenases/reductases.

Review: Protein function at thermal extremes: balancing stability and flexibility

No organism can survive across the entire temperature range found in the biosphere, and a given species can rarely support active metabolism across more than a few tens of degrees C. Nevertheless, life can be maintained at surprisingly extreme temperatures, from below -50 to over 110 degrees C. That proteins, which are assembled with the same 20 amino acids in all species, can function well at both extremes of this range illustrates the plasticity available in the construction of these macromolecules. In studying proteins from extremophiles, researchers have found no new amino acids, covalent modifications or structural motifs that explain the ability of these molecules to function in such harsh environments. Rather, subtle redistributions of the same intramolecular interactions required for protein stabilization at moderate temperatures are sufficient to maintain structural integrity at hot or cold extremes. The key to protein function, whether in polar seas or hot springs, is the maintenance of an appropriate balance between molecular stability on the one hand and structural flexibility on the other. Stability is needed to ensure the appropriate geometry for ligand binding, as well as to avoid denaturation, while flexibility is necessary to allow catalysis at a metabolically appropriate rate. Comparisons of homologous proteins from organisms spanning a wide range of thermal habitats show that adaptive mutations, as well as stabilizing solutes, maintain a balance between these two attributes, regardless of the temperature at which the protein functions.

Structural basis for altered activity of M- and H-isozyme forms of human lactate dehydrogenase

Lactate dehydrogenase (LDH) interconverts pyruvate and lactate with concomitant interconversion of NADH and NAD(+). Although crystal structures of a variety of LDH have previously been described, a notable absence has been any of the three known human forms of this glycolytic enzyme. We have now determined the crystal structures of two isoforms of human LDH-the M form, predominantly found in muscle; and the H form, found mainly in cardiac muscle. Both structures have been crystallized as ternary complexes in the presence of the NADH cofactor and oxamate, a substrate-like inhibitor. Although each of these isoforms has different kinetic properties, the domain structure, subunit association, and active-site regions are indistinguishable between the two structures. The pK(a) that governs the K(M) for pyruvate for the two isozymes is found to differ by about 0.94 pH units, consistent with variation in pK(a) of the active-site histidine. The close similarity of these crystal structures suggests the distinctive activity of these enzyme isoforms is likely to result directly from variation of charged surface residues peripheral to the active site, a hypothesis supported by electrostatic calculations based on each structure. Proteins 2001;43:175-185.

A comparative analysis of the evolutionary patterning and mechanistic bases of lactate dehydrogenase thermal stability in porcelain crabs, genus Petrolisthes

The kinetic properties of orthologous homologs (orthologs) of enzymes are typically correlated with environmental temperatures in species adapted to different thermal regimes, but correlations between adaptation temperature and enzyme thermal stability are less clear. Although the thermal stability of a protein is related chiefly to its primary structure (including post-translational modification), thermal stability can also be altered by extrinsic factors present in the intracellular milieu. Here, we present a comparative analysis of the thermal stability of lactate dehydrogenase (LDH) orthologs from 22 congeneric species of porcelain crab (genera Petrolisthes and Allopetrolisthes) from a broad range of thermal habitats. Interspecific diversity of LDH stability is high: temperatures required for a 50 % loss of activity in 10 min ranged from 65 to 75.5 degrees C, corresponding to half-lives of less than 1 min to more than 3 h at 70 degrees C. Although stability is positively correlated with maximal habitat temperature in some sister taxa, phylogenetic comparative analysis incorporating all 22 species does not indicate that the interspecific diversity of LDH stability represents an adaptive response to current thermal habitats. Examination of the mechanistic bases of LDH stabilization indicates that differences in stability are related both to properties of the LDH molecule itself (intrinsic stability) and to the effects of extrinsic protein(s). Intrinsic differences were shown by the unfolding of structure during heating, as measured by circular dichroism spectroscopy. Stabilizing effects of extrinsic proteins are implied by the results of cellular fractionation experiments that removed low-molecular-mass solutes and proteins from the muscle homogenates. We conclude that the overall structural stability and functional properties of proteins can evolve independently and that in vivo protein-protein interactions can provide another means to regulate protein stability selectively.

Generating ensemble averages for small proteins from extended conformations by Monte Carlo simulations

Since it is not feasible to determine the structure of every protein by experiment, algorithms delivering the folded conformation of a protein solely from its amino acid sequence are desirable. Here the diffusion-process controlled-Monte Carlo approach has been applied to generating ensemble averages for three small proteins with 31, 36, and 46 residues. Starting from extended conformations and using an energy model that was developed on other protein models, the simulations find nativelike structures deviating by 3 A rms from the experimental structures for the main chain atoms. The balance between long-range and short-range interactions is discussed briefly in the context of stability and folding.

Probing the affinity and specificity of yeast alcohol dehydrogenase I for coenzymes

Yeast (Saccharomyces cerevisiae) alcohol dehydrogenase I (SceADH) binds NAD+ and NADH less tightly and turns over substrates more rapidly than does horse (Equus caballus) liver alcohol dehydrogenase E isoenzyme (EcaADH), and neither enzyme uses NADP efficiently. Amino acid residues in the proposed adenylate binding pocket of SceADH were substituted in attempts to improve affinity for coenzymes or reactivity with NADP. Substitutions in SceADH (Gly202Ile or Ser246Ile) with the corresponding residues in the adenine binding site of the homologous EcaADH have modest effects on coenzyme binding and other kinetic constants, but the Ser246Ile substitution decreases turnover numbers by 350-fold. The Ser176Phe substitution (also near adenine site) significantly decreases affinity for coenzymes and turnover numbers. In the consensus nucleotide-binding betaalphabeta fold sequence, SceADH has two alanine residues (177-GAAGGLG-183) instead of the Leu200 in EcaADH (199-GLGGVG-204); the Ala178-Ala179 to Leu substitution significantly decreases affinity for coenzymes and turnover numbers. Some NADP-dependent enzymes have an Ala corresponding to Gly183 in SceADH; the Gly183Ala substitution significantly decreases affinity for coenzymes and turnover numbers. NADP-dependent enzymes usually have a neutral residue instead of the Asp (Asp201 in SceADH) that interacts with the hydroxyl groups of the adenosine ribose, along with a basic residue (at position 202 or 203) to stabilize the 2'-phosphate of NADP. The Gly203Arg change in SceADH does not significantly affect the kinetics. The Gly183Ala or Gly203Arg substitutions do not enable SceADH to use NADP+ as coenzyme. SceADH with the single Asp201Gly or double Asp201Gly:Gly203Arg substitutions have similar, low activity with NADP+. The results suggest that several of the amino acid residues participate in coenzyme binding and that conversion of specificity for coenzyme requires multiple substitutions.

Structural basis for cold adaptation. Sequence, biochemical properties, and crystal structure of malate dehydrogenase from a psychrophile Aquaspirillium arcticum

Aquaspillium arcticum is a psychrophilic bacterium that was isolated from arctic sediment and grows optimally at 4 degrees C. We have cloned, purified, and characterized malate dehydrogenase from A. arcticum (Aa MDH). We also have determined the crystal structures of apo-Aa MDH, Aa MDH.NADH binary complex, and Aa MDH.NAD.oxaloacetate ternary complex at 1.9-, 2.1-, and 2.5-A resolutions, respectively. The Aa MDH sequence is most closely related to the sequence of a thermophilic MDH from Thermus flavus (Tf MDH), showing 61% sequence identity and over 90% sequence similarity. Stability studies show that Aa MDH has a half-life of 10 min at 55 degrees C, whereas Tf MDH is fully active at 90 degrees C for 1 h. Aa MDH shows 2-3-fold higher catalytic efficiency compared with a mesophilic or a thermophilic MDH at the temperature range 4-10 degrees C. Structural comparison of Aa MDH and Tf MDH suggests that the increased relative flexibility of active site residues, favorable surface charge distribution for substrate and cofactor, and the reduced intersubunit ion pair interactions may be the major factors for the efficient catalytic activity of Aa MDH at low temperatures.

A streptavidin mutant with altered ligand-binding specificity

The biotin-binding site of streptavidin was modified to alter its ligand-binding specificity. In natural streptavidin, the side chains of N23 and S27 make two of the three hydrogen bonds with the ureido oxygen of biotin. These two residues were mutated to severely weaken biotin binding while attempting to maintain the affinity for two biotin analogs, 2-iminobiotin and diaminobiotin. Redesigning of the biotin-binding site used the difference in local electrostatic charge distribution between biotin and these biotin analogs. Free energy calculations predicted that the introduction of a negative charge at the position of S27 plus the mutation N23A should disrupt two of the three hydrogen bonds between natural streptavidin and the ureido oxygen of biotin. In contrast, the imino hydrogen of 2-iminobiotin should form a hydrogen bond with the side chain of an acidic amino acid at position 27. This should reduce the biotin-binding affinity by approximately eight orders of magnitude, while leaving the affinities for these biotin analogs virtually unaffected. In good agreement with these predictions, a streptavidin mutant with the N23A and S27D substitutions binds 2-iminobiotin with an affinity (Ka) of 1 x 10(6) M-1, two orders of magnitude higher than that for biotin (1 x 10(4) M-1). In contrast, the binding affinity of this streptavidin mutant for diaminobiotin (2.7 x 10(4) M-1) was lower than predicted (2.9 x 10(5) M-1), suggesting the position of the diaminobiotin in the biotin-binding site was not accurately determined by modeling.

Hot spots in cold adaptation: localized increases in conformational flexibility in lactate dehydrogenase A4 orthologs of Antarctic notothenioid fishes

To elucidate mechanisms of enzymatic adaptation to extreme cold, we determined kinetic properties, thermal stabilities, and deduced amino acid sequences of lactate dehydrogenase A4 (A4-LDH) from nine Antarctic (-1.86 to 1 degree C) and three South American (4 to 10 degree C) notothenioid teleosts. Higher Michaelis-Menten constants (Km) and catalytic rate constants (kcat) distinguish orthologs of Antarctic from those of South American species, but no relationship exists between adaptation temperature and the rate at which activity is lost because of heat denaturation. In all species, active site residues are conserved fully, and differences in kcat and Km are caused by substitutions elsewhere in the molecule. Within geographic groups, identical kinetic properties are generated by different substitutions. By combining our data with A4-LDH sequences for other vertebrates and information on roles played by localized conformational changes in setting kcat, we conclude that notothenioid A4-LDHs have adapted to cold temperatures by increases in flexibility in small areas of the molecule that affect the mobility of adjacent active-site structures. Using these findings, we propose a model that explains linked temperature-adaptive variation in Km and kcat. Changes in sequence that increase flexibility of regions of the enzyme involved in catalytic conformational changes may reduce energy (enthalpy) barriers to these rate-governing shifts in conformation and, thereby, increase kcat. However, at a common temperature of measurement, the higher configurational entropy of a cold-adapted enzyme may foster conformations that bind ligands poorly, leading to high Km values relative to warm-adapted orthologs.

Charge and solvation effects in anion recognition centers: an inquiry exploiting reactive arginines

Following a long-standing suggestion of Riordan et al. [Riordan, J. F., McElvany, K. D., and Borders, C. L., Jr. (1977) Science 195, 884-885], we sought to exploit chemically activated arginines as probes to characterize the microenvironmental effects in enzymes that mediate the recognition of anionic substrates. A micellar simulation study establishes that octylguanidine (OGn) becomes chemically activated upon incorporation into both cetyltrimethylammonium bromide (CTAB) and Triton X-100 micelles and that the activations correlate with the pKa diminutions induced in its guanidinium group by the effects of electrostatic or nonelectrostatic nature as reflected in the results of pH and salt titration experiments. Next, a protein modification study establishes that the modifiable arginines in a number of enzymes also have diminished pKa's, again due to effects of electrostatic or nonelectrostatic nature as reflected in the results of pH and salt titration experiments. Warwicker's finite difference Poisson--Boltzmann algorithm [Warwicker, J. (1992) J. Mol. Biol. 223, 247-257] is applied to several of the enzymes with available crystal structure coordinates, and indeed, their chemically activated arginines are found to be in an electrostatic microenvironment that can diminish their pKa's, with the magnitudes of these diminutions matching closely the diminutions measured experimentally. Finally, the chemically activated arginines are examined with respect to their atomic atmosphere and are thus found to occur in a local microenvironment that would facilitate their roles as anion anchors. Thus, electrostatic and solvation effects are found to be critical determinants of the arginine role as an anion anchor.

Protein folds and functions

BACKGROUND

The recent rapid increase in the number of available three-dimensional protein structures has further highlighted the necessity to understand the relationship between biological function and structure. Using structural classification schemes such as SCOP, CATH and DALI, it is now possible to explore global relationships between protein fold and function, something which was previously impractical.

RESULTS

Using a relational database of CATH data we have generated fold distributions for arbitrary selections of proteins automatically. These distributions have been examined in the light of protein function and bound ligand. Different enzyme classes are not clearly reflected in distributions of protein class and architecture, whereas the type of bound ligand has a much more dramatic effect.

CONCLUSIONS

The availability of structural classification data has enabled this novel overview analysis. We conclude that function at the top level of the EC number enzyme classification is not related to fold, as only a very few specific residues are actually responsible for enzyme activity. Conversely, the fold is much more closely related to ligand type.

Analysis of the structure and substrate binding of Phormidium lapideum alanine dehydrogenase

The structure of the hexameric L-alanine dehydrogenase from Phormidium lapideum reveals that the subunit is constructed from two domains, each having the common dinucleotide binding fold. Despite there being no sequence similarity, the fold of alanine dehydrogenase is closely related to that of the family of D-2-hydroxyacid dehydrogenases, with a similar location of the active site, suggesting that these enzymes are related by divergent evolution. L-alanine dehydrogenase and the 2-hydroxyacid dehydrogenases also use equivalent functional groups to promote substrate recognition and catalysis. However, they are arranged differently on the enzyme surface, which has the effect of directing opposite faces of the keto acid to the dinucleotide in each case, forcing a change in absolute configuration of the product.

Different types of interactions involving cysteine sulfhydryl group in proteins

Various types of interactions involving the sulfhydryl group of free cysteine residues have been analyzed using known protein structures. In a hydrogen bond the -SH group is more amenable to donating its proton to a carbonyl group, rather than acting as a proton acceptor. It rarely interacts with a carboxylate group, and is a poor ligand to bind an anionic substrate. It is quite prone to make contacts that are definitely non-hydrogen bond type. In the S...C=O interaction the S atom is placed on the face of an amide group (mostly from the main-chain, but there are cases from the side-chain also) close to the C atom. Cases of S...N interaction, where the S atom is on top of the N atom of another residue (both main-, as well as side-chains, including the guanidinium group) are also observed. A considerable number of Cys residues have aromatic residues as neighbors, and here too, the preferred mode of interaction is along the face. The intra-residue S...C=O interaction constrains the main-chain and side-chain torsion angles (psi and chi1), whereas the inter-residue interactions are non-local and stabilize the tertiary structure. The S...C=O interaction may have a role in lowering the pKa values of the Cys residues in enzyme active sites.

Flexible ligand docking using conformational ensembles

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.

Substrate and cofactor specificity and selective inhibition of lactate dehydrogenase from the malarial parasite P. falciparum

Lactate dehydrogenase from the malarial parasite Plasmodium falciparum has many amino acid residues that are unique compared to any other known lactate dehydrogenase. This includes residues that define the substrate and cofactor binding sites. Nevertheless, parasite lactate dehydrogenase exhibits high specificity for pyruvic acid, even more restricted than the specificity of human lactate dehydrogenases M4 and H4. Parasite lactate dehydrogenase exhibits high catalytic efficiency in the reduction of pyruvate, kcat/Km = 9.0 x 10(8) min(-1) M(-1). Parasite lactate dehydrogenase also exhibits similar cofactor specificity to the human isoforms in the oxidation of L-lactate with NAD+ and with a series of NAD+ analogs, suggesting a similar cofactor binding environment in spite of the numerous amino acid differences. Parasite lactate dehydrogenase exhibits an enhanced kcat with the analog 3-acetylpyridine adenine dinucleotide (APAD+) whereas the human isoforms exhibit a lower kcat. This differential response to APAD+ provides the kinetic basis for the enzyme-based detection of malarial parasites. A series of inhibitors structurally related to the natural product gossypol were shown to be competitive inhibitors of the binding of NADH. Slight changes in structure produced marked changes in selectivity of inhibition of lactate dehydrogenase. 7-p-Trifluoromethylbenzyl-8-deoxyhemigossylic acid inhibited parasite lactate dehydrogenase, Ki = 0.2 microM, which was 65- and 400-fold tighter binding compared to the M4 and H4 isoforms of human lactate dehydrogenase. The results suggest that the cofactor site of parasite lactate dehydrogenase may be a potential target for structure-based drug design.

Comparison of protein surfaces using a genetic algorithm

A genetic algorithm (GA) is described which is used to compare the solvent-accessible surfaces of two proteins or fragments of proteins, represented by a dot surface calculated using the Connolly algorithm. The GA is used to move one surface relative to the other to locate the most similar surface region between the two. The matching process is enhanced by the use of the surface normals and shape terms provided by the Connolly program and also by a simple hydrogen-bonding descriptor and an additional shape descriptor. The algorithm has been tested in applications ranging from the comparison of small surface patches to the comparison of whole protein surfaces, and it has performed correctly in all cases. Examples of the matches are given and a quantitative analysis of the quality of the matches is performed. A number of possible future enhancements to the program are described which would allow the GA to be used for more complex surface comparisons.

The contribution of electrostatic and van der Waals interactions to the stereospecificity of the reaction catalyzed by lactate dehydrogenase

Continuum electrostatic calculations in conjunction with molecular dynamics simulations have been used to investigate the source of the stereospecificity in the hydride transfer reaction catalyzed by lactate dehydrogenase (LDH). These studies show that favorable electrostatic interactions between the carboxamide group of the reduced nicotinamide adenine dinucleotide coenzyme and protein residues of the active site of LDH can account for much if not all of the stereospecificity of the LDH-catalyzed reaction, with A-side hydride transfer more than 10(7) times greater than B-side transfer. Unfavorable steric interactions within the binding complex for B-side transfer are not found.

Toxoplasma gondii expresses two distinct lactate dehydrogenase homologous genes during its life cycle in intermediate hosts

Two Toxoplasma gondii genes were characterized that are differentially expressed during the parasite's life cycle. The genes named LDH1 and LDH2, respectively, encode polypeptides similar to the enzyme lactate dehydrogenase (LDH; L-lactate:NAD+ oxidoreductase, EC 1.1.1.27) from a variety of organisms. They show 64.0% nucleotide identity in the coding region and both have an intron at the same relative position. The deduced amino acid sequences of LDH1 and LDH2 share 71.1% identity. LDH1 and LDH2 are most similar to an LDH of Plasmodium falciparum (46.5% and 48.5% amino acid identities, respectively). The mRNA of LDH2 was only detected in the bradyzoite stage, while the mRNA of LDH1 was detected in both the bradyzoite and tachyzoite stages. However, by isoelectric focusing and immunoblot analysis, only one LDH isoform was found to be expressed in each stage. Furthermore, the expression of a reporter gene carrying chloramphenicol acetyltransferase (CAT) coding sequence and the putative LDH2 promoter sequence was significantly up-regulated by growing parasites in tissue culture in media with alkaline pH (pH 8.2, a condition known to induce the expression of bradyzoite-specific antigens), while the expression of a CAT reporter construct carrying the putative LDH1 promoter sequence was down-regulated by similar treatment. These results indicate that LDH expression is developmentally regulated in T. gondii and suggest a possible correlation between stage conversion and alteration in carbohydrate or energy metabolism in this parasite.

Removal of substrate inhibition in a lactate dehydrogenase from human muscle by a single residue change

High concentrations of ketoacid substrates inhibit most natural hydroxyacid dehydrogenases due to the formation of an abortive enzyme-NAD+-ketoacid complex. It was postulated that this substrate inhibition could be eliminated from lactate dehydrogenases if the rate of NAD+ dissociation could be increased. An analysis of the crystal structure of mammalian LDHs showed that the amide of the nicotinamide cofactor formed a water-bridged hydrogen bond to S163. The LDH of Plasmodium falciparum is not inhibited by its substrate and, uniquely, in this enzyme the serine is replaced by a leucine. In the S163L mutant of human LDH-M4 pyruvate inhibition is, indeed, abolished and the enzyme retains high activity. However, the major contribution to this effect comes from a weakening of the interaction of pyruvate with the enzyme-coenzyme complex.

Cloning, sequencing and functional expression of a DNA encoding pig cytosolic malate dehydrogenase: purification and characterization of the recombinant enzyme

Using the polymerase chain reaction, DNA encoding cytosolic malate dehydrogenase (cMDH) has been cloned from a pig heart cDNA library. Large amounts of the enzyme (30 mg per litre of original culture) have been produced in Escherichia coli using an inducible expression vector (pKK223-3) in which the 5'-non-coding region of the gene was replaced with the tac promoter. The complete nucleotide sequence of the DNA is reported for the first time. The recombinant cMDH purified was shown to be identical to the native enzyme according to: chromatographic behaviour, isoelectric point, N-terminal amino acid sequence, and physiochemical and catalytic properties.

Arginine to tryptophan substitution in the active site of a human lactate dehydrogenase variant--LDHB GUA1: postulated effects on subunit structure and catalysis

A variant of lactate dehydrogenase (LDHB GUA1) was previously identified among the Guaymi Indians of Panama and Costa Rica. The LDHB GUA1 variant is enzymatically inactive; however, the variant subunits alter the electrophoretic mobility of the tetramers that include active LDHA and LDHB subunits. The kinetic properties of the tetrameric enzyme, comprised of the inactive B plus active A subunits, are similar to properties of the heterotetramers with active B subunits, except for the reduced specific activity. We have determined that a single C.G to T.A transition changes an Arg to a Trp at amino acid residue 106. This substitution explains the increase in net negative charge observed by protein electrophoresis. This Arg 106 residue is absolutely conserved throughout evolution. Published high-resolution crystal structures of LDH reveal that this residue is within the hinge of a loop that closes over the active site of the subunit upon binding of substrate and cofactor and also has a direct role in catalysis. Computer modeling of the variant enzyme suggests that replacement of this Arg residue with a Trp does not induce significant change in the structure of the active site. However, this substitution would result in disruption of enzyme activity through the inability of the uncharged tryptophan side-chain to polarize the substrate carbonyl bond. This would explain the loss of the catalytic function with maintenance of normal kinetic properties in the heterotetramers containing the variant subunits. The ability to maintain normal, tissue-specific kinetic properties could explain the absence of clinical manifestations in the homozygous LDHB GUA1 individuals.

Hydration in drug design. 1. Multiple hydrogen-bonding features of water molecules in mediating protein-ligand interactions

Water is known to play an important role in the recognition and stabilization of the interaction between a ligand and its site. This has important implications for drug design. Analyses of 19 high-resolution crystal structures of protein-ligand complexes reveal the multiple hydrogen-bonding feature of water molecules mediating protein-ligand interactions. Most of the water molecules (nearly 80%) involved in bridging the protein and the ligand can make three or more hydrogen bonds when distance and bond angles are used as criteria to define hydrogen-bonding interactions. Isotropic B-factors have been used to take into account the mobility of water molecules. The water molecules at binding sites bridge the protein and ligand, and interact with other water molecules to form a complex network of interconnecting hydrogen bonds. Some water molecules at the site do not directly bridge between the protein and the ligand, but may contribute indirectly to the stability of the complex by holding bridging water molecules in the right position through a network of hydrogen bonds. These water networks are probably crucial for the stability of the protein-ligand complex and are important for any site-directed drug design strategies.

Hydration in drug design. 2. Influence of local site surface shape on water binding

If water molecules are strongly bound at a protein-ligand interface, they are unlikely to be displaced during ligand binding. Such water molecules can change the shape of the ligand binding site and thus affect strategies for drug design. To understand the nature of water binding, and factors influencing it, water molecules at the ligand binding sites of 26 high-resolution protein-ligand complexes have been examined here. Water molecules bound in deep grooves and cavities between the protein and the ligand are located in the indentations on the protein-site surface, but not in the indentations on the ligand surface. The majority of the water molecules bound in deep indentations on the protein-site surface make multiple polar contacts with the protein surface. This may indicate a strong binding of water molecules in deep indentations on protein-site surfaces. The local shape of the site surface may influence the binding of water molecules that mediate protein-ligand interactions.

Hydration in drug design. 3. Conserved water molecules at the ligand-binding sites of homologous proteins

Water molecules are known to play an important rôle in mediating protein-ligand interactions. If water molecules are conserved at the ligand-binding sites of homologous proteins, such a finding may suggest the structural importance of water molecules in ligand binding. Structurally conserved water molecules change the conventional definition of 'binding sites' by changing the shape and complementarity of these sites. Such conserved water molecules can be important for site-directed ligand/drug design. Therefore, five different sets of homologous protein/protein-ligand complexes have been examined to identify the conserved water molecules at the ligand-binding sites. Our analysis reveals that there are as many as 16 conserved water molecules at the FAD binding site of glutathione reductase between the crystal structures obtained from human and E. coli. In the remaining four sets of high-resolution crystal structures, 2-4 water molecules have been found to be conserved at the ligand-binding sites. The majority of these conserved water molecules are either bound in deep grooves at the protein-ligand interface or completely buried in cavities between the protein and the ligand. All these water molecules, conserved between the protein/protein-ligand complexes from different species, have identical or similar apolar and polar interactions in a given set. The site residues interacting with the conserved water molecules at the ligand-binding sites have been found to be highly conserved among proteins from different species; they are more conserved compared to the other site residues interacting with the ligand. These water molecules, in general, make multiple polar contacts with protein-site residues.

NAD-binding domains of dehydrogenases

The nicotinamide adenine dinucleotide (NAD)-binding domains of dehydrogenases, containing a conserved double beta-alpha-beta-alpha-beta motif, are common structural feature of many enzymes that bind NAD, nicotinamide adenine dinucleotide phosphate (NADP) and related cofactors. Features of this folding pattern that create a natural binding site for such molecules have been described. The domain continues to appear in many structures, in the form of a common core with different peripheral additions or variations. Other structures that bind NAD and related molecules use entirely different topologies, although, in many, a phosphate group appears at the N terminus of an alpha helix. Ferredoxin reductase seems to show convergent evolution, containing a single beta-alpha-beta motif that is similar both in its structure and in its interactions with the ligand to a region in dehydrogenases.

Involvement of the C-terminal tail in the activity of Drosophila alcohol dehydrogenase. Evaluation of truncated proteins constructed by site-directed mutagenesis

Drosophila alcohol dehydrogenase belongs to the heterogeneous family of short-chain dehydrogenases/reductases, which does not include the well characterized mammalian alcohol dehydrogenases. Although it is clear that the main biological role of this enzyme is in alcohol oxidation, in the absence of the three-dimensional conformation only partial information on the protein regions involved in the active site, and the coenzyme and substrate interacting cavities is available. Two segments have already been identified, a coenzyme-binding segment at the N-terminus, and the reactive Tyr152 and Lys156 residues. Limited proteolytic assays had suggested the involvement of the 13 C-terminal amino acids in the function of the enzyme. By site-directed mutagenesis, we have constructed eight different truncated mutant enzymes and expressed them in Escherichia coli. The purified mutant enzymes have been recovered and characterized using monoclonal antibodies. Kinetic analysis and stability assays have been performed, and clearly demonstrate the contribution of the last 13 amino acids to the activity. We hypothesize that the C-terminal tail constitutes an essential region for maintaining the hydrophobicity of the catalytic pocket needed for binding of the substrate.

Redesign of the substrate specificity of Escherichia coli aspartate aminotransferase to that of Escherichia coli tyrosine aminotransferase by homology modeling and site-directed mutagenesis

Although several high-resolution X-ray crystallographic structures have been determined for Escherichia coli aspartate aminotransferase (eAATase), efforts to crystallize E. coli tyrosine aminotransferase (eTATase) have been unsuccessful. Sequence alignment analyses of eTATase and eAATase show 43% sequence identity and 72% sequence similarity, allowing for conservative substitutions. The high similarity of the two sequences indicates that both enzymes must have similar secondary and tertiary structures. Six active site residues of eAATase were targeted by homology modeling as being important for aromatic amino acid reactivity with eTATase. Two of these positions (Thr 109 and Asn 297) are invariant in all known aspartate aminotransferase enzymes, but differ in eTATase (Ser 109 and Ser 297). The other four positions (Val 39, Lys 41, Thr 47, and Asn 69) line the active site pocket of eAATase and are replaced by amino acids with more hydrophobic side chains in eTATase (Leu 39, Tyr 41, Ile 47, and Leu 69). These six positions in eAATase were mutated by site-directed mutagenesis to the corresponding amino acids found in eTATase in an attempt to redesign the substrate specificity of eAATase to that of eTATase. Five combinations of the individual mutations were obtained from mutagenesis reactions. The redesigned eAATase mutant containing all six mutations (Hex) displays second-order rate constants for the transamination of aspartate and phenylalanine that are within an order of magnitude of those observed for eTATase. Thus, the reactivity of eAATase with phenylalanine was increased by over three orders of magnitude without sacrificing the high transamination activity with aspartate observed for both enzymes.(ABSTRACT TRUNCATED AT 250 WORDS)

Novel enzymes as reagents

Enzymes are now used in a wide range of analytical methods, primarily for the measurement of substrates and as labels in immunoassays. Enzymes from microbial sources are becoming the preferred choice because of ease of extraction, catabolic activities and opportunities for enhancing yields. Protein engineering of bacterial enzymes will also play a role in future enzyme design with respect to improvements in specificity, reaction kinetics and stability.

Pediococcus acidilactici ldhD gene: cloning, nucleotide sequence, and transcriptional analysis

The gene encoding D-lactate dehydrogenase was isolated on a 2.9-kb insert from a library of Pediococcus acidilactici DNA by complementation for growth under anaerobiosis of an Escherichia coli lactate dehydrogenase and pyruvate-formate lyase double mutant. The nucleotide sequence of ldhD encodes a protein of 331 amino acids (predicted molecular mass of 37,210 Da) which shows similarity to the family of D-2-hydroxyacid dehydrogenases. The enzyme encoded by the cloned fragment is equally active on pyruvate and hydroxypyruvate, indicating that the enzyme has both D-lactate and D-glycerate dehydrogenase activities. Three other open reading frames were found in the 2.9-kb insert, one of which (rpsB) is highly similar to bacterial genes coding for ribosomal protein S2. Northern (RNA) blotting analyses indicated the presence of a 2-kb dicistronic transcript of ldhD (a metabolic gene) and rpsB (a putative ribosomal protein gene) together with a 1-kb monocistronic rpsB mRNA. These transcripts are abundant in the early phase of exponential growth but steadily fade away to disappear in the stationary phase. Primer extension analysis identified two distinct promoters driving either cotranscription of ldhD and rpsB or transcription of rpsB alone.

A single amino acid mutation enhances the thermal stability of Escherichia coli malate dehydrogenase

The stability of wild-type Escherichia coli malate dehydrogenase was compared with a mutant form of the enzyme with the amino acid residue at position 102 changed from arginine to glutamine. The mutation occurs on the underside of a mobile loop which closes over the active-site cleft on formation of the enzyme/cofactor/substrate ternary complex. The mutant enzyme is kinetically compromised while the wild-type enzyme is highly specific for oxaloacetate. The mutant enzyme was shown to be more resistant to irreversible thermal denaturation by thermal inactivation experiments and high-sensitivity differential scanning calorimetry than the wild-type enzyme. In contrast, resistance of both enzymes to reversible unfolding in guanidinium chloride was similar. Circular dichroic spectropolarimetry shows the secondary structures of the enzymes are similar but there is a demonstrable difference in tertiary structure. From the position of the mutation, it is conjectured that the substitution on a mobile surface loop results in partial closure of the loop and greater resistance to thermal inactivation of the mutant enzyme. However, molecular modelling combined with circular dichroic spectropolarimetry indicate that the mutation may have a more widespread effect on the structure than simply partial closure of the mobile surface loop as the environment of distant tyrosine residues is altered. Resistance of the wild-type enzyme to thermal inactivation can be increased by cofactor addition, which may have the effect of partial closure of the mobile surface loop, but has little effect on the mutant enzyme.

Solvent isotope effects and the nature of electrophilic catalysis in the action of the lactate dehydrogenase of Bacillus stearothermophilus

Deuterium oxide at atom fractions of deuterium from 0.0 to 0.97 has an effect of less than 20% on the kinetic term kcat/KmB (believed to reflect the transition state for the hydride-transfer step) for the reduction of pyruvic acid by NADH at 55 degrees C, with catalysis by the tetrameric form of the lactate dehydrogenase of Bacillus stearothermophilus. This observation suggests that the hydride-transfer event is not assisted by protonic bridging to the carbonyl group being reduced. The results are consistent with protonic bridging only if an opposing isotope effect is present, for example from a generalized conformation or solvation change. The results are consistent with other forms of electrophilic catalysis.

Activity and fluorescence changes of lactate dehydrogenase induced by guanidine hydrochloride in reverse micelles

Denaturants activate several multimeric enzymes in reverse micelles [Garza-Ramos, G., Darszon, A., Tuena de Gómez-Puyou, M. & Gómez-Puyou, A. (1992) Eur. J. Biochem. 205, 509-517]. Here, the effect on activity and intrinsic fluorescence of pig heart lactate dehydrogenase (LDH) in reverse micelles [formed with 0.2 M cetyltrimethylammonium bromide in octane/hexanol (8.6:1, by vol.)] was explored at various water and guanidine hydrochloride (Gdn/HCl) concentrations. Emission fluorescence spectra of LDH in aqueous media and in micelles were similar. As in all aqueous media, 1.0 M Gdn/HCl in the water phase of reverse micelles produced fluorescence quenching and a blue shift of the maximal emission. In 5.0 M Gdn/HCl, instead of the red shift and significant quenching seen in water, the maximum emission further shifted to the blue and was only slightly quenched. Gdn/HCl titrations of activity and fluorescence changes of LDH in micelles with different water contents showed that at Wo ([H2O]/[surfactant]) of 6.6, 8.3, or 12.5, increasing concentrations of Gdn/HCl up to 0.6 M produced small changes in fluorescence, whereas activity increased several-fold. At higher denaturant concentrations, activity decreased with significant fluorescence changes. In reverse micelles with 1 M Gdn/HCl, Vmax but not Km of LDH decreased with time. Under these conditions, there was progressive quenching of LDH fluorescence. The results show that in reverse micelles different Gdn/HCl concentrations induce variations in activity with or without alterations of the intrinsic fluorescence of LDH. The results also indicate that in reverse micelles, concentrations of Gdn/HCl below 1.0 M cause an enhancement of protein flexibility; this is accompanied by a marked increase in activity without important changes in intrinsic fluorescence. 1.0 M Gdn/HCl produces perturbations of inter-subunit contacts that lead to fluorescence quenching and loss of catalytic activity, probably as consequence of dimerization of tetrameric LDH.

Engineering surface loops of proteins--a preferred strategy for obtaining new enzyme function

A prerequisite for the rational redesign of enzymes is that altering amino acids in an attempt to obtain new biological function does not unexpectedly alter the globular, natural framework of the native protein on which the design is being executed. The results of combinatorial-mutagenesis strategies suggest that random variation of amino acid sequence is most easily tolerated at the solvent-exposed surfaces of a protein. This review analyses effective redesigns of Bacillus stearothermophilus lactate dehydrogenase (bsLDH), in which all residue variations are at solvent-exposed surfaces. The majority of these variations were located within surface loops, which interconnect stable secondary structures traversing the globular core of the protein.

Substitution of the amino acid at position 102 with polar and aromatic residues influences substrate specificity of lactate dehydrogenase

The Gln residue at amino acid position 102 of Bacillus stearothermophilus lactate dehydrogenase was replaced with Ser, Thr, Tyr, or Phe to investigate the effect on substrate recognition. The Q102S and Q102T mutant enzymes were found to have a broader range of substrate specificity (measured by kcat/Km) than the wild-type enzyme. However, it is evident that either Ser or Thr at position 102 are of a size able to accommodate a wide variety of substrates in the active site and substrate specificity appears to rely largely on size discrimination in these mutants. The Q102F and Q102Y mutant enzymes have low catalytic efficiency and do not show this relaxed substrate specificity. However, their activities are restored by the presence of an aromatic substrate. All of the enzymes have a very low catalytic efficiency with branched chain aliphatic substrates.

Dissecting the contributions of a specific side-chain interaction to folding and catalysis of Bacillus stearothermophilus lactate dehydrogenase

X-ray crystallography predicts hydrogen-bonding interactions between the side chains of Thr198 and two other amino acid residues, Glu194 (adjacent to the catalytic His195) and Ser318 (on the alpha-H helix which rearranges on substrate binding). In order to investigate the contribution of this conserved amino acid residue, Thr198, two mutants of Bacillus stearothermophilus lactate dehydrogenase were created (Val198 and Ile198). The steady-state kinetic parameters for both mutant enzymes were very similar with increased substrate Km and reduced kcat when compared with the wild-type enzyme. The mutation Val198 allowed non-productive binding of pyruvate to the unprotonated form of His195. Steady-state kinetic parameters determined for the Val198 mutant enzyme in high solvent viscosity suggested both an altered rate-limiting step in catalysis and implicated Thr198 in allosteric activation by the effector fructose 1,6-bisphosphate (Fru1,6P2). A shift in the Fru1,6P2 activation constant for the Val198 mutant enzyme suggested that Thr198 stabilises the catalytically competent (Fru1,6P2-activated) form of the enzyme by 6.6 kJ/mol. However, Thr198 was not important for maintaining the thermal stability of the Fru1,6P2-activated form. Equilibrium unfolding in guanidinium chloride indicated that Thr198 contributes 17.2 kJ/mol subunits towards the tertiary structural stability. The results emphasise the importance of the side chain-hydroxyl group of Thr198 which is required for (a) productive substrate binding, (b) allosteric activation and (c) protein conformational stability. The characteristics of the B. stearothermophilus lactate dehydrogenase mutations reported here were significantly different from those of the same mutations made in the corresponding position of the analogous enzyme Thermus flavus malate dehydrogenase [Nishiyama, M., Shimada, K., Horinouchi, S., & Beppu, T. (1991) J. Biol. Chem. 266, 14294-14299].