Transthyretin (TTR) is a vertebrate-exclusive transport protein that plays a key role in binding and distributing thyroid hormones. However, its evolutionary origin lies in the duplication of the gene that encoding the enzyme 5-hydroxyisourate hydrolase (HIUase), which is involved in uric acid metabolism. Unlike TTR, HIUase is ubiquitous in both prokaryotes and eukaryotes, with the exception of hominids. Both HIUase and TTR subfamilies form homotetramers, possessing an internal charged cavity between the two dimer pairs. Based on their high degree of structural similarity, we hypothesized that specific in silico substitutions would enable the interconversion between these protein functions. Using an evolution-based approach, we engineered two putative protein sequences, where correlated locally conserved positions from one subfamily representative sequence were substituted by the other, and vice versa. Applying computational modeling techniques, the best models were refined, validated, and their cavity volumes, three-dimensional geometries, propensity to aggregation and electrostatic potentials were analyzed. Molecular dynamics simulations were performed with the reference proteins and the engineered mutants in the bound and unbound states. We demonstrate that the volumes and geometries differ from one another, due to size and physicochemical differences between their ligands. The bound state mutant complexes are stable, and the enzymatic assay demonstrated active new enzymes. Our work suggests that the evolution-based protein engineering approach used has residue-specific resolution to identify locally conserved residues in the sequence of evolutionarily related proteins, such as HIUase and TTR.
Keywords: 5-hydroxyisourate hydrolase; Transthyretin; molecular modeling; protein design strategies; protein engineering; protein evolution.