A binary complex of malate dehydrogenase from the thermophilic bacterium Thermus flavus (tMDH) with NADH has been crystallized from poly(ethylene glycol) 3500, pH 8.5, yielding diffraction-quality crystals in space group P2(1)2(1)2(1). The structure was solved at 1.9-A resolution using molecular replacement and refined to an R factor of 15.8% with good geometry. The primary sequence of tMDH is 55% identical to that of cytoplasmic malate dehydrogenase (cMDH) [Birktoft, J. J., Rhodes, G., & Banaszak, L. J. (1989) Biochemistry 28, 6065-6081], and overall their three-dimensional structures are very similar. Like cMDH, tMDH crystallized as a dimer with one coenzyme bound per subunit. The coenzyme binds in the extended conformation, and most of the interactions with enzyme are similar to those in cMDH. In tMDH, small local conformational changes are caused by the replacement of a glutamic acid for the aspartic acid involved in hydrogen bonding to the adenine ribose of NADH. Comparison of tMDH with cMDH reveals that both tMDH subunits more closely resemble the B subunit of cMDH which therefore is the more likely representative of the solution conformation. While cMDH is inactivated at temperatures above about 50 degrees C, tMDH is fully active at 90 degrees C. On the basis of the X-ray crystal structure, a number of factors have been identified which are likely to contribute to the relative thermostability of tMDH compared to cMDH. The most striking of the differences involves the introduction of four ion pairs per monomer. All of these ion pairs are solvent-accessible. Three of these ion pairs are located in the dimer interface, Glu27-Lys31, Glu57-Lys168, and Glu57-Arg229, and one ion pair, Glu275-Arg149, is at the domain interface within each subunit. Additionally, we observe incorporation of additional alanines into alpha-helices of tMDH and, in one instance, incorporation of an aspartate that functions as a counterchange to an alpha-helix dipole. The possible contributions of these and other factors to protein thermostability in tMDH are discussed.

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.