Abstract
Oculocutaneous albinism type 3 (OCA3) is an autosomal recessive disorder caused by mutations in the TYRP1 gene. Tyrosinase-related protein 1 (Tyrp1) is involved in eumelanin synthesis, catalyzing the oxidation of 5,6-dihydroxyindole-2-carboxylic acid oxidase (DHICA) to 5,6-indolequinone-2-carboxylic acid (IQCA). Here, for the first time, four OCA3-causing mutations of Tyrp1, C30R, H215Y, D308N, and R326H, were investigated computationally to understand Tyrp1 protein stability and catalytic activity. Using the Tyrp1 crystal structure (PDB:5M8L), global mutagenesis was conducted to evaluate mutant protein stability. Consistent with the foldability parameter, C30R and H215Y should exhibit greater instability, and two other mutants, D308N and R326H, are expected to keep a native conformation. SDS-PAGE and Western blot analysis of the purified recombinant proteins confirmed that the foldability parameter correctly predicted the effect of mutations critical for protein stability. Further, the mutant variant structures were built and simulated for 100 ns to generate free energy landscapes and perform docking experiments. Free energy landscapes formed by Y362, N378, and T391 indicate that the binding clefts of C30R and H215Y mutants are larger than the wild-type Tyrp1. In docking simulations, the hydrogen bond and salt bridge interactions that stabilize DHICA in the active site remain similar among Tyrp1, D308N, and R326H. However, the strengths of these interactions and stability of the docked ligand may decrease proportionally to mutation severity due to the larger and less well-defined natures of the binding clefts in mutants. Mutational perturbations in mutants that are not unfolded may result in allosteric alterations to the active site, reducing the stability of protein-ligand interactions.