For determining the mutations necessary for the creation of orthogonal pairs our written protocol was used. Our starting Protein DataBase file (PDB) was the PDB submitted by C. Ottmann et al. containing a complex of a T14-3-3 primer, fusicoccin, and 2 CT52 (PDB ID: 2O98) in which chain A and B are T14-3-3 monomers and chain P and Q are CT52 proteins.
After subjecting the PDB to relax and backrub simulations in Rosetta several of the best scoring PDBs (according to Rosetta) were submitted to the Robetta server1 for computational alanine scanning, together with a mutations list containing the residues that are relevant for the binding of the T14-3-3-FC-CT52 complex based on a paper by C. Ottmann et al.2
The alanine scans gave us insight in which residues were the most critical for the binding interactions and therefore the best options for mutation testing.
Figure 1 displays the Computational Alanine scan in the form of a heat map. It shows that the most important residues in the T14-3-3 are 56 63 136 181 182 189 229 233, and for CT52 932 935 936 939 943 945 947 948 950 951 952 953 954 955 956.
We create an orthogonal pair (see figure 2) by first destabilizing the binding interface by introducing a mutation in the CT52, since CT52 is a more flexible protein and thus more prone to changes in secondary structure than T14-3-3.
In the next step we restabilize the binding interface by introducing compensating mutations in the T14-3-3. In the first step all relevant residues found with the Computational alanine scanning are mutated to 19 different amino acids (all except cysteine) and then evaluated on their destabilizing effect (see figure 3A). Higher destabilization means a better orthogonality between the mutated CT52 and the wildtype T14-3-3.
In the second step all residues in T14-3-3 that are within 6 Å of the mutated residue are one by one mutated into 19 mutants to test the restabilizing effects of these mutations. Some of these mutations were combined into two point mutations to increase the restabilizing effect. This step ensures that the mutated CT52 binds to the mutated T14-3-3. Energy increase/decrease of the mutated pair compared to the wildtype pair are shown in figure 3C. The third step is to determine the orthogonality, that is the energy increase, between the wildtype CT52 and the mutated T14-3-3, this is shown in figure 3B.
Based on the results represented in
figure 3A-CFor example, we first looked at heatmap A, the darker the colour, the stronger the destabilizing effect and thus the better the orthogonality between the mutated CT52 and the wildtype T14-3-3: if you look at the residue 947 on the y-axis you can see several dark colors, one of which is at mutation F on the x-axis. Next, we looked at the same coordinates on heatmap B (so residue 947 mutation F) to see whether those same coordinates also have a dark color in this heatmap. It did, which means the destabilizing effect of the mutated T14-3-3 and wildtype CT-52 is also quite strong. Lastly we looked at heatmap C, for the binding to still be strong enough to actually bind, the free Gibbs energy should preferably be lower than the wildtype-wildtype binding, this means the color in the heatmap should be blue for the chosen coordinates. For the coordinates residue 947 mutation F, the color is indeed slightly blue. For our application we need to be able to regulate the interaction with fusicoccin, so the CT52 shouldn’t bind strongly to the scaffold without fusicoccin present, this means that the Gibbs free energy should not decrease too much, which is the case when the color is dark blue, that this happens. This is the reason that, for example, the mutation I947W wasn’t chosen.
, several potential mutations were chosen to test in wet lab experiments, see
table 1.
The chosen mutations have high ΔΔG changes in figure 3A and 3B, and are blue (meaning it has a decrease in ΔΔG compared to the wildtype). A lower energy in our mutated pair is preferred over a higher energy, not only because a lower energy means a more stable binding, but also because Rosetta’s point mutant scan application tends to be a bit off in its predictions and thus a high energy could be even more unfavorable than predicted. For our application, however, the decrease in energy should not be too large, as this could decrease the regulability of our complex.
In figure 4 these mutations are presented in a crystal structure of the T14-3-3/CT52 complex.
In order to determine the properties of our mutated complexes, most importantly the binding affinity/dissociation constant and the cooperativity (which is defined as the increase/decrease in binding affinity of the second CT52 because of the first bound CT52), a model based on mass-action and Michaelis-Menten kinetics was developed to convert the activity measured in the readout to these properties.