Mutations that significantly change the stability, flexibility and quaternary structure of the l-lactate dehydrogenase from Bacillus megaterium

Structure, Function, and Thermodynamics of Lactate Dehydrogenases from Humans and the Malaria Parasite P. falciparum

Temperature adaptation is ubiquitous among all living organisms, yet the molecular basis for this process remains poorly understood. It can be assumed that for parasite-host systems, the same enzymes found in both organisms respond to the same selection factor (human body temperature) with similar structural changes. Herein, we report the existence of a reversible temperature-dependent structural transition for the glycolytic enzyme lactate dehydrogenase (LDH) from the malaria parasite Plasmodium falciparum (pfLDH) and human heart (hhLDH) occurring in the temperature range of human fever. This transition is observed for LDHs from psychrophiles, mesophiles, and moderate thermophiles in their operating temperature range. Thermodynamic analysis reveals unique thermodynamic signatures of the LDH-substrate complexes defining a specific temperature range to which human LDH is adapted and parasite LDH is not, despite their common mesophilic nature. The results of spectroscopic analysis combined with the available crystallographic data reveal the existence of an active center within pfLDH that imparts psychrophilic structural properties to the enzyme. This center consists of two pockets, one formed by the five amino acids (5AA insert) within the substrate specificity loop and the other by the active site, that mutually regulate one another in response to temperature and induce structural and functional changes in the Michaelis complex. Our findings pave the way toward a new strategy for malaria treatments and drug design using therapeutic agents that inactivate malarial LDH selectively at a specific temperature range of the cyclic malaria paroxysm.

Interaction of extracellular domain 2 of the human retina-specific ATP-binding cassette transporter (ABCA4) with all-trans-retinal

The retina-specific ATP-binding cassette (ABC) transporter, ABCA4, is essential for transport of all-trans-retinal from the rod outer segment discs in the retina and is associated with a broad range of inherited retinal diseases, including Stargardt disease, autosomal recessive cone rod dystrophy, and fundus flavimaculatus. A unique feature of the ABCA subfamily of ABC transporters is the presence of highly conserved, long extracellular loops or domains (ECDs) with unknown function. The high degree of sequence conservation and mapped disease-associated mutations in these domains suggests an important physiological significance. Conformational analysis using CD spectroscopy of purified, recombinant ECD2 protein demonstrated that it has an ordered and stable structure composed of 27 +/- 3% alpha-helix, 20 +/- 3% beta-pleated sheet, and 53 +/- 3% coil. Significant conformational changes were observed in disease-associated mutant proteins. Using intrinsic tryptophan fluorescence emission spectrum of ECD2 polypeptide and fluorescence anisotropy, we have demonstrated that this domain specifically interacts with all-trans-retinal. Furthermore, the retinal interaction appeared preferential for the all-trans-isomer and was directly measurable through fluorescence anisotropy analysis. Our results demonstrate that the three macular degeneration-associated mutations lead to significant changes in the secondary structure of the ECD2 domain of ABCA4, as well as in its interaction with all-trans-retinal.

A novel computer algorithm improves antibody epitope prediction using affinity-selected mimotopes: a case study using monoclonal antibodies against the West Nile virus E protein

Understanding antibody function is often enhanced by knowledge of the specific binding epitope. Here, we describe a computer algorithm that permits epitope prediction based on a collection of random peptide epitopes (mimotopes) isolated by antibody affinity purification. We applied this methodology to the prediction of epitopes for five monoclonal antibodies against the West Nile virus (WNV) E protein, two of which exhibit therapeutic activity in vivo. This strategy was validated by comparison of our results with existing F(ab)-E protein crystal structures and mutational analysis by yeast surface display. We demonstrate that by combining the results of the mimotope method with our data from mutational analysis, epitopes could be predicted with greater certainty. The two methods displayed great complementarity as the mutational analysis facilitated epitope prediction when the results with the mimotope method were equivocal and the mimotope method revealed a broader number of residues within the epitope than the mutational analysis. Our results demonstrate that the combination of these two prediction strategies provides a robust platform for epitope characterization.

Cloning of an epoxide hydrolase-encoding gene from Aspergillus niger M200, overexpression in E. coli, and modification of activity and enantioselectivity of the enzyme by protein engineering

The gene encoding an epoxide hydrolase from Aspergillus niger M200 has been cloned and its sequence determined. The gene is interrupted by seven introns, one exon being only nine nucleotides long. The non-coding 5'- and 3'-regions of the mRNA are composed of 47 and 76 nucleotides, respectively. Overexpression of the fungal epoxide hydrolase in E. coli TOP10 has led to a 15-fold increase in specific activity (compared to the wild-type strain). Saturation mutagenesis at codon 217 resulted in the discovery of nine enzyme variants showing in several cases profound differences in activity and enantioselectivity towards various epoxides when compared to the data of the wild-type enzyme. The site 217 is located at the entrance of the tunnel that provides the substrate with access to the active site. The exchange of Ala at this position for Cys has led to a doubled enantioselectivity (E-value of 5.0) towards benzyl glycidyl ether. The same substitution resulted in a threefold-enhanced activity of the enzyme towards allyl glycidyl ether and styrene oxide without affecting enantioselectivity. The variant A217L showed an enhanced enantioselectivity towards tert-butyl glycidyl ether reaching an E-value of 100 (from 60 for the wild-type enzyme). Replacement of A217 by Val has led to higher activity towards allyl glycidyl ether by a factor of six. The substitutions Ala-->Glu and Ala-->Gln increased the enantioselectivity towards allyl glycidyl ether and styrene oxide by over 50% to E-values of 10 and 16, respectively. The study underlines that single amino acid exchanges in the substrate tunnel region can lead to significant improvements in enantioselectivity and activity of the epoxide hydrolase from A. niger M200.

A Rapid, Three-Step Process for the Preformulation of aRecombinant Ricin Toxin A-Chain Vaccine

A systematic, three-step approach was employed to develop a stable, optimized formulation of ricin toxin A-chain V76M/Y80A (rRTA) for use as a vaccine against ricin toxicity. The method first uses spectroscopic techniques to evaluate the stability of rRTA as a function of temperature and pH. To synthesize the data, empirical phase diagrams are generated to display the conditions under which the protein maintains particular conformational states. Following identification of optimal pH conditions, light scattering and fluorescence assays are employed to screen a wide variety of compounds for their abilities to stabilize rRTA. Once stabilizers were identified, the ability of rRTA to adsorb to aluminum salt adjuvants was evaluated. Desorption of the protein from the adjuvant was also analyzed. Using this approach, the optimal formulation conditions for rRTA were determined to be pH 6.0 utilizing glycerol as a stabilizer and Alhydrogel as an adjuvant. Such an approach has the potential to significantly reduce the time it takes to get vaccines into clinical testing.

A systematic approach to stabilizing EBA-175 RII-NG for use as a malaria vaccine

Region II of the erythrocyte-binding antigen (EBA-175 RII) has been identified as a promising target for a malaria vaccine. A systematic approach to identify optimal preformulation conditions of a non-glycosylated (NG) antigen, EBA-175 RII-NG, has been developed. This approach consists of development of an empirical temperature/pH phase diagram, high throughput stabilizer screening and aluminum salt adjuvant adsorption studies. Using these physical methods, we developed a stable formulation for EBA-175 RII-NG at pH 6.0 with sucrose and Brij 35 as stabilizers and Adju-Phos as an adjuvant. This approach should be generally applicable to guiding the development of stable vaccine formulations.

Linkage between fructose 1,6-bisphosphate binding and the dimer-tetramer equilibrium of Escherichia coli glycerol kinase: critical behavior arising from change of ligand stoichiometry

Escherichia coli glycerol kinase (EC 2.7.1.30; ATP-glycerol 3-phosphotransferase) is inhibited allosterically by fructose 1,6-bisphosphate (FBP), and this inhibition is a primary mechanism by which glucose controls glycerol utilization in vivo. Earlier work indicates that glycerol kinase displays a dimer-tetramer equilibrium in solution, FBP shifts the equilibrium toward the tetramer, and tetramer formation is required for FBP inhibition. However, equilibrium constants for FBP binding and dimer-tetramer assembly that describe the linkage between these processes are unknown. Here, decreased fluorescence anisotropy of extrinsic fluorophores fluorescein and 2',7'-difluorofluorescein due to homo fluorescence resonance energy transfer (homo-FRET) is used to quantitate tetramer assembly and FBP binding. Glycerol kinase is labeled with extrinsic fluorophores covalently attached to an engineered surface cysteine residue under conditions that prevent labeling of native cysteine residues. Tryptic peptide mapping and MALDI-MS verify labeling at the engineered site only. Initial velocity studies show the labeling does not alter the catalytic properties or FBP inhibition. The steady-state fluorescence anisotropy of enzyme with a labeling stoichiometry of approximately 0.1 mol of fluorophore/mol of subunit is not sensitive to increased protein concentration or binding of FBP, indicating the absence of homo-FRET. However, steady-state fluorescence anisotropy of enzyme with a labeling stoichiometry of approximately 0.4 mol of fluorophore/mol of subunit decreases with increasing protein concentration, which is consistent with depolarization due to homo-FRET. The protein concentration dependence of the decreased fluorescence anisotropy is described by a dimer-tetramer equilibrium with an apparent dissociation constant of 61 +/- 7 nM (subunits) at pH 7.0 and 25 degrees C. FBP binds to both the dimer and tetramer of glycerol kinase, and the FBP concentration dependence of the apparent dissociation constant for the dimer-tetramer equilibrium shows critical behavior. The apparent dissociation constant decreases and then increases with increasing FBP concentration, reaching a minimum at about 20 mM FBP. Critical behavior is seen also in the FBP dependence of the inhibition. The critical behavior arises because tetramer dissociation increases FBP stoichiometry from two sites per tetramer to four half-sites per two dimers. The phenomenological description of the coupling between tetramer assembly and FBP binding shows antagonistic binding of FBP to the two sites on the tetramer, indicating that the strong positive cooperativity observed for FBP inhibition of catalytic activity (Hill coefficient approximately 1.5) is due to the approximately 4000-fold higher affinity of the tetramer for FBP rather than to positive coupling between the two FBP sites.

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.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.

Protein thermal stability: insights from atomic displacement parameters (B values)

The factors contributing to the thermal stability of proteins from thermophilic origins are matters of intense debate and investigation. Thermophilic proteins are thought to possess better packed interiors than their mesophilic counterparts, leading to lesser overall flexibility and a corresponding reduction in surface-to-volume ratio. These observations prompted an analysis of B values reported in high-resolution X-ray crystal structures of mesophilic and thermophilic proteins. In this analysis, the following aspects were addressed: (1) frequency distribution of normalized B values (B' factors) over all the proteins and for individual amino acids; (2) amino acid compositions in high B value regions of polypeptide chains; (3) variation in the B values from core to the surface of proteins in terms of their radius of gyration; and (4) degree of dispersion of normalized B values in spheres around the Calpha atoms. The analysis revealed that (1) Ser and Thr have lesser flexibility in thermophiles than in mesophiles, (2) the proportion of Glu and Lys in high B value regions of thermophiles is higher and that of Ser and Thr is lower and (3) the dispersion of B values within spheres at Calpha atoms is similar in mesophiles and thermophiles. These observations reflect plausible differences in the dynamics of thermophilic and mesophilic proteins and suggest amino acid substitutions that are likely to change thermal stability.

Directed evolution of a thermostable esterase

We have used in vitro evolution to probe the relationship between stability and activity in a mesophilic esterase. Previous studies of these properties in homologous enzymes evolved for function at different temperatures have suggested that stability at high temperatures is incompatible with high catalytic activity at low temperatures through mutually exclusive demands on enzyme flexibility. Six generations of random mutagenesis, recombination, and screening stabilized Bacillus subtilis p-nitrobenzyl esterase significantly (>14 degreesC increase in Tm) without compromising its catalytic activity at lower temperatures. Furthermore, analysis of the stabilities and activities of large numbers of random mutants indicates that these properties are not inversely correlated. Although enhanced thermostability does not necessarily come at the cost of activity, the process by which the molecule adapts is important. Mutations that increase thermostability while maintaining low-temperature activity are very rare. Unless both properties are constrained (by natural selection or screening) the evolution of one by the accumulation of single amino acid substitutions typically comes at the cost of the other, regardless of whether the two properties are inversely correlated or not correlated at all.

Extremely thermostable L(+)-lactate dehydrogenase from Thermotoga maritima: cloning, characterization, and crystallization of the recombinant enzyme in its tetrameric and octameric state

L(+)-lactate dehydrogenase (LDH; E.C.1.1.1.27) from the hyperthermophilic bacterium Thermotoga maritima has been shown to represent the most stable LDH isolated so far (Wrba A, Jaenicke R, Huber R, Stetter KO, 1990, Eur J Biochem 188:195-201). In order to obtain the enzyme in amounts sufficient for physical characterization, and to analyze the molecular basis of its intrinsic stability, the gene was cloned and expressed functionally in Escherichia coli. Growth of the cells and purification of the enzyme were performed aerobically at 26 degrees C, i.e., ca. 60 degrees below the optimal growth temperature of Thermotoga. Two enzyme species with LDH activity were purified to homogeneity. Crystals of the enzyme obtained at 4 degrees C show satisfactory diffraction suitable for X-ray analysis up to a resolution of 2.8 A. As shown by gel-permeation chromatography, chemical crosslinking, light scattering, analytical ultracentrifugation, and electron microscopy, the two LDH species represent homotetramers and homooctamers (i.e., dimers of tetramers), with a common subunit molecular mass of 35 kDa. The spectroscopic characteristics (UV absorption, fluorescence emission, near- and far-UV CD) of the two species are indistinguishable. The calculated alpha-helix content is 45%, in accordance with the result of homology modeling. Compared to the tetrameric enzyme, the octamer exhibits reduced specific activity, whereas KM is unalatered. The extreme intrinsic stability of the protein is reflected by its unaltered catalytic activity over 4 h at 85 degrees C; irreversible thermal denaturation becomes significant at approximately 95 degrees C. The anomalous resistance toward chemical denaturation using guanidinium chloride and urea confirms this observation. Both the high optimal temperature and the pH optimum of the catalytic activity correspond to the growth conditions of T. maritima in its natural habitat.

Determinants of enzyme thermostability observed in the molecular structure of Thermus aquaticus D-glyceraldehyde-3-phosphate dehydrogenase at 25 Angstroms Resolution

The crystal structure of holo D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from the extreme thermophile Thermus aquaticus has been solved at 2.5 Angstroms resolution. To study the determinants of thermostability, we compare our structure to four other GAPDHs. Salt links, hydrogen bonds, buried surface area, packing density, surface to volume ratio, and stabilization of alpha-helices and beta-turns are analyzed. We find a strong correlation between thermostability and the number of hydrogen bonds between charged side chains and neutral partners. These charged-neutral hydrogen bonds provide electrostatic stabilization without the heavy desolvation penalty of salt links. The stability of thermophilic GAPDHs is also correlated with the number of intrasubunit salt links and total hydrogen bonds. Charged residues, therefore, play a dual role in stabilization by participating not only in salt links but also in hydrogen bonds with a neutral partner. Hydrophobic effects allow for discrimination between thermophiles and psychrophiles, but not within the GAPDH thermophiles. There is, however, an association between thermostability and decreasing enzyme surface to volume ratio. Finally, we describe several interactions present in both our GAPDH and a hyperthermophilic GAPDH that are absent in the less thermostable GAPDHs. These include a four-residue salt link network, a hydrogen bond near the active site, an intersubunit salt link, and several buried Ile residues.

Engineering thermostability: lessons from thermophilic proteins

As several groups begin to tap the rich pickings found in the Archaea--a vast kingdom that stretches the concept of life as we know it--the structures of proteins from hyperthermophiles are being elucidated. Certain features are beginning to emerge, such as compactness and hydrophobic clustering, but the ability to engineer these features into temperature-intolerant proteins is still some way off.

The crystal structure of citrate synthase from the thermophilic archaeon, Thermoplasma acidophilum

BACKGROUND

The Archaea constitute a phylogenetically distinct, evolutionary domain and comprise organisms that live under environmental extremes of temperature, salinity and/or anaerobicity. Different members of the thermophilic Archaea tolerate temperatures in the range 55-110 degrees C, and the comparison of the structures of their enzymes with the structurally homogolous enzymes of mesophilic organisms (optimum growth temperature range 15-45 degrees C) may provide important information on the structural basis of protein thermostability. We have chosen citrate synthase, the first enzyme of the citric acid cycle, as a model enzyme for such studies.

RESULTS

We have determined the crystal structure of Thermoplasma acidophilum citrate synthase to 2.5 A and have compared it with the citrate synthase from pig heart, with which it shares a high degree of structural homology, but little sequence identity (20%).

CONCLUSIONS

The three-dimensional structural comparison of thermophilic and mesophilic citrate synthases has permitted catalytic and substrate-binding residues to be tentatively assigned in the archaeal, thermophilic enzyme, and has identified structural features that may be responsible for its thermostability.

Influence of a mutation in the transmembrane domain of the p185c-erbB2 oncogene-encoded protein studied by molecular dynamics simulations

The c-erbB2 proto-oncogene encodes for a protein of 185kDa (p185) which becomes transforming upon the Val-->Glu transmembrane amino acid substitution. The transforming ability seems to be due to a substitution-resulting constitutive activation of the tyrosine kinase cytosolic domain of the protein. These observations prompted us to evaluate the structural and dynamical behavior of the transmembrane region of the wild and transforming p185 protein in order to understand the role of this region in the transduction mechanism. 160 ps molecular dynamics simulations in vacuo have been performed on two peptides corresponding to the sequence [651-679] of p185c-erbB2 protein and its transforming mutant Val659-->Glu659. These two sequences include the transmembrane domain and are initially postulated to be in an alpha-helix conformation. Noticeable differences in the flexibility of the two peptides are shown. The nontransforming sequence seems rather flexible and several conformational changes are detected at the junction of the mutation point [658-659] and at position Val665-Val666 during the 160 ps simulations. On the contrary, no transitions were observed for the mutated sequence which adopts a stable alpha-helix conformation. This difference in flexibility could be hypothesized as a factor involved in the regulation of the tyrosine kinase activity of p185c-erbB2.

Protein engineering for unusual environments

The ability to use proteins in unusual or non-natural environments greatly expands their potential applications in biotechnology. Because natural selection has neither maximized the stability of proteins nor optimized them to function under unusual conditions, there is considerable room for their improvement by protein engineering. Significant advances reported within the past year include a dramatic demonstration of a protein's ability to tolerate changes in its amino acid sequence, large increases in protein stability, and the use of random mutagenesis to obtain novel enzymatic properties. Approaches using random or site-directed mutagenesis have been successful in generating proteins able to function in an extended range of environments.