We report a crystallographic and computational analysis of two mutant forms of the Alicyclobacillus acidocaldarius thioredoxin (BacTrx) done in order to evaluate the contribution of two specific amino acids to the thermostability of BacTrx. Our results suggest that the thermostability of BacTrx may be modulated by mutations affecting the overall electrostatic energy of the protein.In the past few decades, many studies have aimed at unraveling the relationship between protein structure and thermostability. Sequence and structural comparisons between functionally homologous proteins from mesophilic and thermophilic organisms have pointed out that high thermal stability is usually related to multiple factors (8,11,16,17,19,20,25,37,39,41,42,(45)(46)(47)(48). Recently, a series of findings have suggested that thermostability can be also achieved by single ad hoc mutations, possibly allowing the formation of ion pairs or ion networks (26, 44) and/or optimization of long-range coulombic interactions on the protein surface (15,23,36,40). A striking example of this phenomenon is given by previous studies of two mutant forms, Lys18Gly and Arg82Glu, of the thermophilic Alicyclobacillus acidocaldarius thioredoxin (BacTrx) (5, 28) and of a Glu85Arg mutant of the mesophilic Escherichia coli thioredoxin (Trx) (22). Molecular dynamic simulations (33-35) allowed the identification of two surface residues able to modulate the thermostability of both proteins. In fact, replacement of Lys18 or Arg82 of BacTrx with the corresponding residues of Trx (Gly21 and Glu85, respectively) led to a dramatic decrease in thermostability (34). Alternatively, the reverse replacement of Glu85 with one Arg residue made Trx remarkably more thermostable (35).Here we report a crystallographic analysis of the Lys18Gly (K18G) and Arg82Glu (R82E) mutant forms of BacTrx that was done in order to ascertain the structural consequences of these two mutations. Moreover, on the basis of the experimental structures, we performed computational studies on K18G, R82E, and the Glu85Arg (E85R) mutant form of Trx, both at ab initio quantum mechanical and molecular mechanic levels, to evaluate both the energetic contribution of a hydrogen bond highlighted by the X-ray experiments and the effect of the mutations on the overall electrostatic energy of the proteins.Crystallization and data collection. Crystals of the K18G and R82E mutant forms were grown by the hanging-drop method at 22°C. Two microliters of a 16-mg/ml protein solution was mixed with 2 l of reservoir solution. K18G crystals grew in 0.1 M HEPES (pH 8)-25% polyethylene glycol 8000-0.2 M calcium acetate. R82E crystals grew from a solution containing 0.1 M cacodylate (pH 6.5)-25% polyethylene glycol 8000-0.1 M zinc acetate.Diffraction data at 1.90 Å resolution for both K18G and R82E were collected at the EMBL Beamline X11 of the Deutsches Elektronen-Synchrotron, Hamburg, Germany, at 100K. For data collection, the crystals were cryoprotected by addition of 10% glycerol to the reservoir solution and flash-fro...