A mutation at the interface between domains causes rearrangement of domains in 3-isopropylmalate dehydrogenase

Elucidating the preference of dimeric over monomeric form for thermal stability of Thermus thermophilus isopropylmalate dehydrogenase: A molecular dynamics perspective

An oligomer usually refers to a macromolecular complex formed by non-covalent interactions of monomers. Several thermophilic proteins are oligomers. The significance of oligomerization of individual proteins for stability at higher temperature is of prime importance for understanding evolution and increasing industrial productivity. The functional form of Thermus thermophilius isopropylmalate dehydrogenase (IPMDH), a widely studied protein to understand the factors affecting the thermal stability of a protein is a dimer, a simplest oligomer. To decipher the relationship between the effects of oligomerization on thermal stability of a protein, we have applied all-atom molecular mechanics approach by analyzing how temperature effects dynamics of a subunit in the presence and absence of another subunit in dimeric (SS) and monomeric forms (SA), respectively, before its denaturation begins. Comparing the difference in overall dynamic structural aspects at two different temperatures, 300 K and 337 K. Analysis of root mean square deviation (RMSD), root mean square fluctuations (RMSF) and Cα-Cα distance with an increase in temperature from 300 K to 337 K for a total of 0.2 μs reveals higher thermal stability of the dimer as compared to monomer. In contrast to dimeric form, the monomer is relatively stable at 300 K but cannot withstand the structural stability at 337 K leading to loosening of intramolecular interactions with maximum fluctuation at B23-B24 within a subunit. Energetic and structural properties indicate that B24-B24' is the major contributor to maintaining subunit-subunit interaction at 337 K. Correlation between the favorable interaction energy (IE) with the minimal perturbance in Cα atoms of domain 2 in a subunit in the presence of another subunit enhances the rigidity of the domain with subunit-subunit interaction. Overall, the study indicates that the dimeric over monomeric form enhances the protein's thermal stability and not all major subunit interacting regions contribute equally in maintaining the former.

Deciphering the Dynamics of Non-Covalent Interactions Affecting Thermal Stability of a Protein: Molecular Dynamics Study on Point Mutant of Thermus thermophilus Isopropylmalate Dehydrogenase

Thermus thermophilius isopropylmalate dehydrogenase catalyzes oxidative decarboxylation and dehydrogenation of isopropylmalate. Substitution of leucine to alanine at position 172 enhances the thermal stability among the known point mutants. Exploring the dynamic properties of non-covalent interactions such as saltbridges, hydrogen bonds and hydrophobic interactions to explain thermal stability of a protein is interesting in its own right. In this study dynamic changes in the non-covalent interactions are studied to decipher the deterministic features of thermal stability of a protein considering a case study of a point mutant in Thermus thermophilus isopropylmalate dehydrogenase. A total of four molecular dynamic simulations of 0.2 μs were carried out on wild type and mutant's functional dimers at 300 K and 337 K. Higher thermal stability of the mutant as compared to wild type is revealed by root mean square deviation, root mean square fluctuations and Cα-Cα distance with an increase in temperature from 300 K to 337 K. Most of the regions of wild type fluctuate higher than the corresponding regions of mutant with an increase in temperature. Cα-Cα distance analysis suggests that long distance networks are significantly affected in wild type as compared to the mutant. Short lived contacts are higher in wild type, while long lived contacts are lost at 337 K. The mutant forms less hydrogen bonds with water as compared to wild type at 337 K. In contrast to wild type, the mutant shows significant increase in unique saltbridges, hydrogen bonds and hydrophobic contacts at 337 K. The current study indicates that there is a strong inter-dependence of thermal stability on the way in which non-covalent interactions reorganize, and it is rewarding to explore this connection in single mutant studies.

Epistasis effects of multiple ancestral-consensus amino acid substitutions on the thermal stability of glycerol kinase from Cellulomonas sp. NT3060

Thermostable variants of the Cellulomonas sp. NT3060 glycerol kinase have been constructed by through the introduction of ancestral-consensus mutations. We produced seven mutants, each having an ancestral-consensus amino acid residue that might be present in the common ancestors of both bacteria and of archaea, and that appeared most frequently at the position of 17 glycerol kinase sequences in the multiple sequence alignment. The thermal stabilities of the resulting mutants were assessed by determining their melting temperatures (Tm), which was defined as the temperature at which 50% of the initial catalytic activity is lost after 15 min of incubation, as well as when the half-life of the catalytic activity occurs at a temperature of 60°C (t1/2). Three mutants showed increased stabilities compared to the wild-type protein. We then produced five more mutants with multiple amino acid substitutions. Some of the resulting mutants showed thermal stabilities much greater than those expected given the stabilities of the respective mutants with single mutations. Therefore, the effects of mutations are not always simply additive and some amino acid substitutions, which do not affect or only slightly improve stability when individually introduced into the protein, show substantial stabilizing effects in combination with other mutations.

Induced allostery in the directed evolution of an enantioselective Baeyer-Villiger monooxygenase

The molecular basis of allosteric effects, known to be caused by an effector docking to an enzyme at a site distal from the binding pocket, has been studied recently by applying directed evolution. Here, we utilize laboratory evolution in a different way, namely to induce allostery by introducing appropriate distal mutations that cause domain movements with concomitant reshaping of the binding pocket in the absence of an effector. To test this concept, the thermostable Baeyer-Villiger monooxygenase, phenylacetone monooxygenase (PAMO), was chosen as the enzyme to be employed in asymmetric Baeyer-Villiger reactions of substrates that are not accepted by the wild type. By using the known X-ray structure of PAMO, a decision was made regarding an appropriate site at which saturation mutagenesis is most likely to generate mutants capable of inducing allostery without any effector compound being present. After screening only 400 transformants, a double mutant was discovered that catalyzes the asymmetric oxidative kinetic resolution of a set of structurally different 2-substituted cyclohexanone derivatives as well as the desymmetrization of three different 4-substituted cyclohexanones, all with high enantioselectivity. Molecular dynamics (MD) simulations and covariance maps unveiled the origin of increased substrate scope as being due to allostery. Large domain movements occur that expose and reshape the binding pocket. This type of focused library production, aimed at inducing significant allosteric effects, is a viable alternative to traditional approaches to "designed" directed evolution that address the binding site directly.

Correlating protein function and stability through the analysis of single amino acid substitutions

Background

Mutations resulting in the disruption of protein function are the underlying causes of many genetic diseases. Some mutations affect the number of expressed proteins while others alter the activity on a per-molecule basis. Single amino acid substitutions as caused by non-synonymous Single Nucleotide Polymorphisms (nsSNPs) often disrupt function by altering protein structure and/or stability, but can also wreak havoc by directly impacting functional binding sites. Given the experimental three-dimensional (3D) structure of a protein, we can try to differentiate between the "effect on structure/stability" and the "effect on binding". However, experimental 3D structures are available for only 1% of all known proteins; the magnitude of stability change caused by a given mutation is more widely available.

Results

Here, we analyze to which extent the functional effect of a mutation can be predicted from the effect on protein stability. We find that simple sequence-based methods succeed in predicting functional effects of nsSNPs. In fact, such methods consistently outperform approaches that predict functional change through the application of binary thresholds to stability change. We also observed that if stability is affected, functional change is easier to predict than when stability is not affected.

Conclusion

Our results confirmed that stability change is somehow related to function change. However, we also show that the knowledge of stability changes in no way suffices to predict functional changes and that many function changing mutations have no effect on stability.

Cold-adaptation mechanism of mutant enzymes of 3-isopropylmalate dehydrogenase from Thermus thermophilus

Random mutagenesis of Thermus thermophilus 3-isopropylmalate dehydrogenase revealed that a substitution of Val126Met in a hinge region caused a marked increase in specific activity, particularly at low temperatures, although the site is far from the binding residues for 3-isopropylmalate and NAD. To understand the molecular mechanism, residue 126 was substituted with one of eight other residues, Gly, Ala, Ser, Thr, Glu, Leu, Ile or Phe. Circular dichroism analyses revealed a decreased thermal stability of the mutants (Delta T ((1/2))= 0-13 degrees C), indicating structural perturbations caused by steric conflict with surrounding residues having larger side chains. Kinetic parameters, k(cat) and K(m) values for isopropylmalate and NAD, were also affected by the mutation, but the resulting k(cat)/K(m) values were similar to that of the wild-type enzyme, suggesting that the change in the catalytic property is caused by the change in free-energy level of the Michaelis complex state relative to that of the initial state. The kinetic parameters and activation enthalpy change (Delta H (double dagger)) showed good correlation with the van der Waals volume of residue 126. These results suggested that the artificial cold adaptation (enhancement of k(cat) value at low temperatures) resulted from the destabilization of the ternary complex caused by the increase in the volume of the residue at position 126.

Crystal structures of mutants of Thermus thermophilus IPMDH adapted to low temperatures

Random mutagenesis on thermophilic 3-isopropylmalate dehydrogenases (IPMDH; EC 1.1.1.85) produced mutant enzymes which adapt to low temperatures. These mutants had higher activity at lower temperatures than the wild-type enzyme without losing high thermostability. Here we report three structures of the mutants of Thermus thermophilus IPMDH determined by X-ray diffraction which was adapted to a low-temperature environment. Two of them have unstable coenzyme binding states and the other one probably has a stable substrate binding state. The present research suggests that the adaptation is correlated with the binding of either coenzyme or the substrate.

Crystal structures of 3-isopropylmalate dehydrogenases with mutations at the C-terminus: crystallographic analyses of structure-stability relationships

Thermal stability of the Thermus thermophilus isopropylmalate dehydrogenase enzyme was substantially lost upon the deletion of three residues from the C-terminus. However, the stability was partly recovered by the addition of two, four and seven amino acid residues (called HD177, HD708 and HD711, respectively) to the C-terminal region of the truncated enzyme. Three structures of these mutant enzymes were determined by an X-ray diffraction method. All protein crystals belong to space group P2(1) and their structures were solved by a standard molecular replacement method where the original dimer structure of the A172L mutant was used as a search model. Thermal stability of these mutant enzymes is discussed based on the 3D structure with special attention to the width of the active-site groove and the minor groove, distortion of beta-sheet pillar structure and size of cavity in the domain-domain interface around the C-terminus. Our previous studies revealed that the thermal stability of isopropylmalate dehydrogenase increases when the active-site cleft is closed (the closed form). In the present study it is shown that the active-site cleft can be regulated by open-close movement of the minor groove located at the opposite side to the active-site groove on the same subunit, through a paperclip-like motion.

Effect of polar side chains at position 172 on thermal stability of 3-isopropylmalate dehydrogenase from Thermus thermophilus

To understand the role of the amino acid residue at position 172 in the conformational stability, four mutant enzymes of Thermus thermophilus 3-isopropylmalate dehydrogenase in which Ala172 was replaced with Asp, Glu, Asn, and Gln were prepared by site-directed mutagenesis. Three mutants were more stable than the wild-type enzyme. No significant change in catalytic properties was found in the mutant enzymes. The molecular modeling studies suggested that the enhanced thermostability of the mutant enzymes resulted from the formation of extra electrostatic interactions and/or improvement of hydrophobic packing of the interior core.