Team:TU-Eindhoven/Modeling Results

iGEM TU Eindhoven

For determining the mutations necessary for the creation of orthogonal pairs our written protocol was used. Our starting PDB was the PDB submitted by C. Ottman 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 P and Q are CT52 proteins.

Computational Alanine Scanning

After relaxing and backrubbing the PDB several of the best scoring PDBs were submitted to the Robetta server[1] 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. Ottman 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: The Results of the computational alanine scanning represented as a heat map with each row containing a PDB resulting from the backrub simulation and each column a residue represented by the ID of that residue in the original PDB. Values of the heat map are ΔΔG values in Rosetta Energy Units (R.E.U.)

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.

Point Mutant Scanning
Figure 2: A schematic representation of the definition of an orthogonal pair: a mutated CT52 only interacts with a mutated T14-3-3 and the wildtype only interacts with the wildtype

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.

Figure 3: The effect of mutations depicted on heat maps, displayed in the rows are the relevant residues of the CT52 and in the columns th e a ŵ iŶ o a Đ id to ǁh iĐ h th e resid u e is ŵ u ta ted , ǀa lu es d isp la yed a re th e Δ Δ G ǀa lu es iŶ R .E .U . A) effeĐt o f d esta ď iliziŶ g ŵ u tations in the CT52 on the wildtype T14-3-3 mutated CT52 binding interface. B) effect of destabilizing mutations in the T14-3-3 on the mutated T14-3-3 wildtype CT52 binding interface. Displayed mutations are not the mutations in the T14-3-3, but the mutations in the CT52 for which these T14-3-3 mutations compensate. C) The energy increase/decrease of the mutated pair compared to the wildtype pair.
Choosing Mutations

Based on the results represented in figure 3A-C, 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.

Determining the properties of the mutated proteins

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. link to model (WIP)