Difference between revisions of "Structure-based mutation analysis GLA"
m (→D244H (Mutation 7)) |
(→R363C (Mutation 10)) |
||
(43 intermediate revisions by the same user not shown) | |||
Line 18: | Line 18: | ||
==GROMACS== |
==GROMACS== |
||
[http://www.gromacs.org/ GROMACS] is mostly a package to perform molecular dynamics, but it also provides energy calculations. |
[http://www.gromacs.org/ GROMACS] is mostly a package to perform molecular dynamics, but it also provides energy calculations. |
||
− | For the mutations we used the forcefield AMBER03 and for the wildtype AMBER03, AMBERGS and CHARMM27. |
+ | For the mutations we used the forcefield [http://www.gromacs.org/Documentation/Terminology/Force_Fields/AMBER AMBER03] and for the wildtype [http://www.gromacs.org/Documentation/Terminology/Force_Fields/AMBER AMBER03], [http://www.gromacs.org/Documentation/Terminology/Force_Fields/AMBER AMBERGS] and [http://www.gromacs.org/Documentation/Terminology/Force_Fields/CHARMM CHARMM27]. |
Additionally to the energy calculation task we did a runtime analysis with values from nsteps=10 to nsteps=1500. The results are shown in the results section of this task. |
Additionally to the energy calculation task we did a runtime analysis with values from nsteps=10 to nsteps=1500. The results are shown in the results section of this task. |
||
According to the task description we created an MDP file with the following content: |
According to the task description we created an MDP file with the following content: |
||
Line 94: | Line 94: | ||
|constraints || Which constraints should be used |
|constraints || Which constraints should be used |
||
|} |
|} |
||
− | Within the GROMACS work step we used the script fetchpdb. It checks if the given input is a valid PDB entry. If the check was successful it downloads the PDB file, extracts it and removes the packed version. |
||
+ | Within the GROMACS work step we used the script fetchpdb. It checks if the given input is a valid PDB entry. If the check was successful it downloads the PDB file, extracts it and removes the packed version. |
||
− | =Results= |
||
+ | |||
+ | =Results & Discussion= |
||
==Structure Selection== |
==Structure Selection== |
||
There are [http://www.uniprot.org/uniprot/P06280#section_x-ref several structure files available] for our protein: |
There are [http://www.uniprot.org/uniprot/P06280#section_x-ref several structure files available] for our protein: |
||
Line 149: | Line 150: | ||
[[File:GLA_structure_mutation_mapped_mutations.png|thumb|400px|center|Figure 1: Representation of the protein α-galactosidase A. The residues which will be mutated are colored red. Asp170 and A231 are part of the active site and colored cyan. The ligand is colored green.]] |
[[File:GLA_structure_mutation_mapped_mutations.png|thumb|400px|center|Figure 1: Representation of the protein α-galactosidase A. The residues which will be mutated are colored red. Asp170 and A231 are part of the active site and colored cyan. The ligand is colored green.]] |
||
− | |||
− | * close to active site? |
||
− | * comparison of side chain conformation |
||
− | ** clashes? |
||
− | * polar interactions |
||
− | * hydrophobic/hydrophilic? |
||
− | * surface |
||
===M42T (Mutation 1)=== |
===M42T (Mutation 1)=== |
||
Line 172: | Line 166: | ||
===S65T (Mutation 2)=== |
===S65T (Mutation 2)=== |
||
− | The side chain |
+ | The side chain conformations of the wildtype and the mutated residues are shown in figure 6 to 8. The hydroxyl and methyl group of threonine point towards totally different direction in the tool-based conformation (see figure 6 and 8), but there are no clashes with other residues. The wildtype S65 forms five hydrogen bonds with surrounding residues (C63, K67, L68, F69). The mutated residue T65 forms also five hydrogen bonds, but one of them is with E66 instead of K67. |
Serine is hydrophilic and the residue is exposed on the surface, but threonine is also a hydrophilic residue. There is no remarkable change of the surface of the protein. |
Serine is hydrophilic and the residue is exposed on the surface, but threonine is also a hydrophilic residue. There is no remarkable change of the surface of the protein. |
||
Line 185: | Line 179: | ||
Since the tool-based side chain conformation of the mutated residue (see figure 11) is pretty similar to the side chain conformation of the wildtype (see figure 9), there are no clashes in the mutated structure. The angle of the hydroxyl group is slightly different in the manual side chain conformation (see figure 10). |
Since the tool-based side chain conformation of the mutated residue (see figure 11) is pretty similar to the side chain conformation of the wildtype (see figure 9), there are no clashes in the mutated structure. The angle of the hydroxyl group is slightly different in the manual side chain conformation (see figure 10). |
||
− | The |
+ | The carboxyl group in the backbone of I117 forms two hydrogen bonds with L120 and A121. These bonds are also abundant in the mutated structure. Serine is a hydrophilic amino acid, but this is no problem, because it is part of the surface. Overall there is no remarkable difference in the surface of the protein due to the mutation. |
{| border="0" style="text-align:center; border-spacing:0" cellpadding="2" cellspacing="3" align = "center" |
{| border="0" style="text-align:center; border-spacing:0" cellpadding="2" cellspacing="3" align = "center" |
||
Line 227: | Line 221: | ||
Since the structure of proline is a special case, the comparison of the side chain conformation is not trivial. Threonine has a higher space requirement than proline. The difference between the manual and tool-based side chain conformation is a variation of one angle in the side chain. There are two small collisions with surrounding residues (I219 and N228) due to the higher space requirements of threonine (see figure 25). |
Since the structure of proline is a special case, the comparison of the side chain conformation is not trivial. Threonine has a higher space requirement than proline. The difference between the manual and tool-based side chain conformation is a variation of one angle in the side chain. There are two small collisions with surrounding residues (I219 and N228) due to the higher space requirements of threonine (see figure 25). |
||
− | P205 forms two hydrogen bond with M208 and W209. The bond with M208 is also formed by T205, but the hydrogen bond to W209 is lost. Therefor a hydrogen bond to N228 is established. The substitution of proline to threonine does not have any influence on the surface of the protein. |
+ | P205 forms two hydrogen bond with M208 and W209. The bond with M208 is also formed by the mutated residue T205, but the hydrogen bond to W209 is lost. Therefor a hydrogen bond to N228 is established. The substitution of proline to threonine does not have any influence on the surface of the protein. |
{| border="0" style="text-align:center; border-spacing:0" cellpadding="2" cellspacing="3" align = "center" |
{| border="0" style="text-align:center; border-spacing:0" cellpadding="2" cellspacing="3" align = "center" |
||
Line 237: | Line 231: | ||
===D244H (Mutation 7)=== |
===D244H (Mutation 7)=== |
||
− | The side chain of the wildtype has a very flat angle to the alpha-helix (see figure 26), whereas the manual and tool-based side chain conformations show an almost orthogonal orientation (see figure 27 and 28). Hence there are no collisions with other residues. |
+ | The side chain of the wildtype has a very flat angle to the alpha-helix (see figure 26), whereas the manual and tool-based side chain conformations show an almost orthogonal orientation (see figure 27 and 28). Hence there are no collisions with other residues. The manual and tool-based conformations are also very similar to each other. |
− | D244 forms four hydrogen bonds. The carboxyl group of the side chain forms two bonds with arginine at position 356. The remaining two bonds are formed by the backbone of D244 with K240 and F248. These two hydrogen bonds are also established by the mutated residue H244 and the side chain forms a third bond S247. |
+ | D244 forms four hydrogen bonds. The carboxyl group of the side chain forms two bonds with arginine at position 356. The remaining two bonds are formed by the backbone of D244 with K240 and F248. These two hydrogen bonds are also established by the mutated residue H244 and the side chain forms a third bond with S247. |
This position is exposed on the surface of the protein and aspartic acid and histidine are hydrophilic amino acids. The orthogonal angle of histidine to the alpha-helix adds a bulge and a small hole to the surface of the protein (see figure 29). |
This position is exposed on the surface of the protein and aspartic acid and histidine are hydrophilic amino acids. The orthogonal angle of histidine to the alpha-helix adds a bulge and a small hole to the surface of the protein (see figure 29). |
||
Line 264: | Line 258: | ||
===Q321E (Mutation 9)=== |
===Q321E (Mutation 9)=== |
||
− | * close to active site? |
||
− | * comparison of side chain conformation |
||
− | ** clashes? |
||
− | * polar interactions |
||
− | * hydrophobic/hydrophilic? |
||
− | * surface |
||
− | |||
Since the structures of glutamine and glutamic acid are so similar, the side chain conformation shows almost no differences (see figure 33, 34 and 35). There is only one angle in the side chain which is slightly different that does not lead to a collision with other residues. |
Since the structures of glutamine and glutamic acid are so similar, the side chain conformation shows almost no differences (see figure 33, 34 and 35). There is only one angle in the side chain which is slightly different that does not lead to a collision with other residues. |
||
Line 284: | Line 271: | ||
===R363C (Mutation 10)=== |
===R363C (Mutation 10)=== |
||
+ | The side chain confirmations of this position are shown in figures 36 to 38. The conformation between the manual and the tool-based ones is very similar. Arginine has a much higher space requirment, hence there is no collusion with other residues due to the substitution to cysteine. |
||
− | * close to active site? |
||
+ | |||
− | * comparison of side chain conformation |
||
+ | The side chain of R363 forms five hydrogen bonds with D335, N336 and N355. These bonds are all lost due to the substitution to cysteine. The backbone of R363 establishs another two hydrogen bonds with I407 which are preserved in the mutation C363. So overall there is a loss of five hydrogen bonds. |
||
− | ** clashes? |
||
+ | |||
− | * polar interactions |
||
+ | Arginine is a hydrophilic amino acid and R363 is exposed on the surface. The substitution to the hydrophobic amino acid cysteine could be a problem. The mutation also introduces a big hole in the surface (see figure 39). |
||
− | * hydrophobic/hydrophilic? |
||
+ | |||
− | * surface |
||
{| border="0" style="text-align:center; border-spacing:0" cellpadding="2" cellspacing="3" align = "center" |
{| border="0" style="text-align:center; border-spacing:0" cellpadding="2" cellspacing="3" align = "center" |
||
− | | [[Image:Fabry_Disease_1R47_res363_wt.png|thumb| |
+ | | [[Image:Fabry_Disease_1R47_res363_wt.png|thumb|150px|Figure 36: Close-up of arginine (wildtype) at position 363 in the protein GLA.]] |
− | | [[Image:Fabry_Disease_1R47_res363_mut.png|thumb| |
+ | | [[Image:Fabry_Disease_1R47_res363_mut.png|thumb|150px|Figure 37: Close-up of cysteine (mutated) at position 363 in the protein GLA. The side chain conformation was done manually.]] |
− | | [[Image:Fabry_Disease_1R47_res363_tool.png|thumb| |
+ | | [[Image:Fabry_Disease_1R47_res363_tool.png|thumb|150px|Figure 38: Close-up of cysteine (mutated) at position 363 in the protein GLA. The side chain conformation was done by SCRWL.]] |
+ | |[[Image:Fabry_Disease_mut10_surface.png|thumb|400px|Figure 39: Surface representation of the protein GLA. (A) The wildtype residue arginine at position 363 is colored green. (B) The mutated residue cysteine at position 363 is colored red. The mutation leads to a hole in the surface of the protein.]] |
||
|} |
|} |
||
==Energy Comparison== |
==Energy Comparison== |
||
− | The results of the energy comparison are presented in the table below. |
+ | The results of the energy comparison are presented in the table below. Due to the fact that the result of minimise for mutation 8 clearly differed from the other results, the run was repeated with the outcome from the first run as input. Thus, there is the number 8.2 in the table. |
+ | |||
− | Due to the fact that the result of the first run of the eighth mutation clearly differed from the other results, the run was repeated with the outcome from the first run as input. Thus, there is the number 8.2. This observation shows that minimise has a decreased tolerance for clashes in comparison to the other tools. Their results for the eighth run are not outstanding and seem not to be affected by the fact that a proline was inserted into a helix. Furthermore, their results seem to be almost equally with respect to some variance. Only the comparison of the FoldX results of the mutations with the wildtype show, that the inserted mutations have a huge influence on the energy of the protein. |
||
{| class="wikitable" border="1" style="text-align:center; border-spacing:0;" |
{| class="wikitable" border="1" style="text-align:center; border-spacing:0;" |
||
|- |
|- |
||
− | ! Number |
+ | ! rowspan=2 scope=col width=60 | Number |
− | ! AA-Position |
+ | ! rowspan=2 scope=col width=80 | AA-Position |
+ | ! rowspan=2 scope=col width=150 | Amino acid change |
||
− | ! Codon change |
||
+ | ! rowspan=2 scope=col width=80 | SCWRL4 |
||
− | ! Amino acid change |
||
+ | ! colspan= 2 | FoldX |
||
− | ! SCWRL4 |
||
+ | ! colspan= 2 | Minimise |
||
− | ! FoldX |
||
+ | ! colspan= 2 | Gromacs |
||
− | ! FoldX - Diff |
||
− | ! Minimise |
||
− | ! Minimise - Diff |
||
|- |
|- |
||
+ | ! scope=col width=70 | Energy |
||
− | |colspan="4"| WT||-||-20.93 ||- ||-20481.23 ||- |
||
+ | ! scope=col width=80 | Difference |
||
+ | ! scope=col width=70 | Energy |
||
+ | ! scope=col width=80 | Difference |
||
+ | ! scope=col width=70 | Energy |
||
+ | ! scope=col width=80 | Difference |
||
|- |
|- |
||
− | | |
+ | |colspan="3"| Wildtype ||-||-20.93 ||- ||-20481.23 ||- || -91307.7 || |
|- |
|- |
||
− | | |
+ | | 1 || 42 || Met -> Thr ||343.25||157.29 ||-178.22 ||-20324.41 ||-156.82 || -90528.4 || -779.3 |
|- |
|- |
||
− | | |
+ | | 2 || 65 || Ser -> Thr ||327.798||152.87 ||-173.8 ||-20339.34 ||-141.89 || -90481.9 || -825.8 |
|- |
|- |
||
− | | |
+ | | 3 || 117 || Ile -> Ser ||333.027||157.97 ||-178.9 ||-20353.47 ||-127.76 || -90654 || -653.7 |
|- |
|- |
||
− | | |
+ | | 4 || 143 || Ala -> Thr ||333.944||154.40 || -175.33||-20339.32 ||-141.91 || -90541 || -766.7 |
|- |
|- |
||
− | | |
+ | | 5 || 186 || His -> Arg ||323.717||154.57 ||-175.5 ||-20321.32 ||-159.91 || -91011.7 || -296 |
|- |
|- |
||
− | | |
+ | | 6 || 205 || Pro -> Thr ||340.619||155.96 ||-176.89 ||-20345.87 ||-135.36 || -90782.2 || -524.8 |
|- |
|- |
||
− | | |
+ | | 7 || 244 || Asp -> His ||333.594||152.08 ||-173.01 ||-20393.12 ||-88.11 || -90232.9 || -1074.8 |
|- |
|- |
||
− | |8 |
+ | | 8 || 283 || Gln -> Pro ||332.631||159.91 ||-180.84 ||-8027.71 ||-12453.52 || -87316 || -3991.7 |
|- |
|- |
||
− | | |
+ | |8.2||-||-||-||-||-||-19134.48||-1346.95 || - || - |
|- |
|- |
||
− | | |
+ | | 9 || 321 || Gln -> Glu ||332.853||160.95 ||-181.88 ||-20246.98 ||-234.25 || -90090.3 || -1217.4 |
+ | |- |
||
+ | | 10 || 363 || Arg -> Cys ||330.56||150.50 ||-171.43 ||-20295.77 ||-185.46 || -89721.8 || -1585.9 |
||
|- |
|- |
||
|} |
|} |
||
+ | |||
+ | |||
+ | ===FoldX=== |
||
+ | There is no mutation that stands out with a very high or very low energy. Every energy difference ranges from approximately -170 to -185. Hence it is very tough to make a clear distinction between neutral and non-neutral mutations. |
||
+ | |||
+ | ===Minimise=== |
||
+ | The results of Minimise cover a much wider range than the results of FoldX. This is mainly due to the result of mutation 8 which shows a significantly lower result than the others. Mutation 9 and 10 should also be mentioned as mutations with a very high difference to the wildtype. Mutation 7 has the lowest difference and the other mutations are somewhere in between. |
||
+ | |||
+ | ===Gromacs=== |
||
+ | Mutation 5 is the mutation with the lowest difference to the wildtype. Once again, mutation 8 has the highest difference. The next least stable structures are the ones of mutation 10, 9 and 7. Mutation 5 has the lowest difference |
||
+ | |||
+ | ==Individual Evaluation of the Mutations== |
||
+ | We evaluate the results of the visual examination and the energy comparison in the following sections individual for each mutation and try to predict the influence of the mutation on the protein. |
||
+ | |||
+ | ===M42T (Mutation 1)=== |
||
+ | The structural influence of the substitution from methionine to threonine at this position can be reduced to the change of the surface in this case, since this is the only aspect in which a remarkable difference occurs. Another aspect is the loss of disulfide bonds due to the mutation which could lead to a loss of stability of the protein. The [http://www.uniprot.org/uniprot/P06280#section_features UniProt entry] does not contain a disulfid bond at this position. The calculated energy differences show a slight tendency towards a non-neutral mutation. Because of this and the remarkable change of the proteins surface, we assume that his mutation is non-neutral. |
||
+ | |||
+ | Prediction: a '''non-neutral''' mutation. |
||
+ | |||
+ | ===S65T (Mutation 2)=== |
||
+ | The side chain conformation of the mutated residue is slightly different, but there a no clashes with other residues and the hydrogen bonds are preserved or substituted, respectively. The calculated energy differences are among the lowest and since there is no significant structual change, we assume a neutral mutation. |
||
+ | |||
+ | Prediction: a '''neutral''' mutation. |
||
+ | |||
+ | ===I117S (Mutation 3)=== |
||
+ | The substitution of isoleucine to serine does not lead to a remarkable change in any of the visual examined aspects. The energy differences of Gromacs and Minimise are on the lower end of all mutations. Because there is no aspect that indicates a influence on the functionality of the protein, we predict a neutral mutation. |
||
+ | |||
+ | Prediction: a '''neutral''' mutation. |
||
+ | |||
+ | ===A143T (Mutation 4)=== |
||
+ | This mutation has a huge influecne on the structure of the protein. There are two aspects that show a remarkable change. First, the mutations leads to a clash with another residue. Second, there is a change of the surface of the protein which is the major aspect. Since this mutation is nearby the active site, the ligand is unlikely to be able to bind at the active site and this will probably have a strong influence on the functionality of the protein. The energy comparison does not show such a strong tendency to a non-neutral mutation, but the observed changes do not lead neccessarily to a change of the energy. Nevertheless, this mutation is very likely to be a non-neutral mutation. |
||
+ | |||
+ | Prediction: a '''non-neutral''' mutation. |
||
+ | |||
+ | ===H186R (Mutation 5)=== |
||
+ | Even though the substitution to arginine leads to clash with another residue at this position, the energy values does not indicate an extraordinary result in relation to the other mutations. Quite the contrary, the value of Gromacs is one of the lowest. Since the collision is very small, the hydrogen bonds of this position are presvered and the calculated energy differences are quite low, we assume that this mutation is neutral. |
||
+ | |||
+ | Prediction: a '''neutral''' mutation. |
||
+ | |||
+ | ===P205T (Mutation 6)=== |
||
+ | Even though the substitution of proline to threonine leads to collisions with two other residues, the energy comparison does not indicate a less stable protein in relation to the energies of the other mutations. This mutation has the third or second lowest difference in the results of Minimise or Gromacs, respectively. Since the energy comparison indicates a neutral mutation and the collisions are very slight, we assume that this mutation is neutral. |
||
+ | |||
+ | Prediction: a '''neutral''' mutation. |
||
+ | |||
+ | ===D244H (Mutation 7)=== |
||
+ | Overall there is a loss of hydrogen bond due to this mutation. The energy results are quite controversial, because it has one of the highest differences in the results of Gromacs, but the lowest in the results of Minimise. It is also on the lower end in the results of the FoldX, but its result are not that meaningful. |
||
+ | |||
+ | The mutation also changes the surface of the protein, but it is tough to evaluate the influence on the functionality of the protein. Since the energy results are leaning towards a neutral mutation, we assume that this mutation is neutral. |
||
+ | |||
+ | Prediction: a '''neutral''' mutation. |
||
+ | |||
+ | ===Q283P (Mutation 8)=== |
||
+ | It is very likely that the substitution to proline leads to a partial loss of the alpha-helix in this area. This is supported by the results of the energy calculations which show the significantly highest difference of all mutations. We predict that this mutation is non-neutral, since the structual influence of this mutation is probably very high. |
||
+ | |||
+ | Prediction: a '''non-neutral''' mutation. |
||
+ | |||
+ | ===Q321E (Mutation 9)=== |
||
+ | Even though the glutamic acid is very similar to glutamine and the mutation has the same conformation, the engery results strongly indicate a non-neutral mutation. The difference of this mutation is among the highest for all thre programs. The two missing hydrogen bonds are probably partially responsible for this. Because the calculated energies show such a clear result, we predict a non-neutral mutation. |
||
+ | |||
+ | Preidction: a '''non-neutral''' mutation. |
||
+ | |||
+ | ===R363C (Mutation 10)=== |
||
+ | The mutation leads to a loss of five hydrogen bonds. This could result in a loss of stability of the protein and it could explain the results of the energy calculation which show high differences (Minimise and Gromacs). The substitution also introduces remarkable changes on the surface. |
||
+ | |||
+ | Since the structual examination and the energy calculations indicate a strong influence on the functionality of the protein, we assume that this mutation is non-neutral. |
||
+ | |||
+ | Prediction: a '''non-neutral''' mutation. |
||
+ | |||
+ | ==Overall Evaluation== |
||
+ | ===Overview=== |
||
+ | The correct assignment (solution) and the predictions of the ten mutations are listed in the table below. |
||
+ | |||
+ | We applied the following coloring scheme: |
||
+ | * green: the prediction was correct. |
||
+ | * red: the prediction was wrong. |
||
+ | |||
+ | |||
+ | {| class="wikitable" border="1" style="text-align:center; border-spacing:0;" |
||
+ | |- |
||
+ | ! scope=col width=80 | Number || scope=col width=80 | Mutation || scope=col width=80 | Solution || scope=col width=80 | Prediction |
||
+ | |- |
||
+ | |1 || M42T || non-neutral || style="background: #9ACD32;" | non-neutral |
||
+ | |- |
||
+ | |2 || S65T || neutral || style="background: #9ACD32;" | neutral |
||
+ | |- |
||
+ | |3 || I117S || neutral || style="background: #9ACD32;" | neutral |
||
+ | |- |
||
+ | |4 || A143T || non-neutral || style="background: #9ACD32;" | non-neutral |
||
+ | |- |
||
+ | |5 || H186R || neutral || style="background: #9ACD32;" | neutral |
||
+ | |- |
||
+ | |6 || P205T || non-neutral || style="background: #FF4500;" | neutral |
||
+ | |- |
||
+ | |7 || D244H || non-neutral || style="background: #FF4500;" | neutral |
||
+ | |- |
||
+ | |8 || Q238P || non-neutral || style="background: #9ACD32;" | non-neutral |
||
+ | |- |
||
+ | |9 || Q321E || non-neutral || style="background: #9ACD32;" | non-neutral |
||
+ | |- |
||
+ | |10 || R363C || non-neutral || style="background: #9ACD32;" | non-neutral |
||
+ | |} |
||
+ | |||
+ | |||
+ | '''Mutation 6: P205T''' |
||
+ | |||
+ | This mutation seemed to be non-neutral because of a collision with another residue, but the energy calculations clearly indicated a neutral mutation. Hence, we may have overweighed the energy calculation of the programs. |
||
+ | |||
+ | |||
+ | '''Mutation 7: D244H''' |
||
+ | |||
+ | The solution is quite surprising, because the energy difference was the lowest in the results of Minimise. Otherwise we would have probably predicted a non-neutral mutation. |
||
+ | |||
+ | ===Discussion of Structure-based Mutation Analysis=== |
||
+ | At first, we thought it would be easier to make a prediction about the mutation than the sequence-based analysis, because we figured that we would be able to see the structual influence of the mutation more clearly. On the one hand, it was tough to evaluate the influence of the individual changes like change of the surface or loss of hydrogen bonds. On the other hand, it was also hard how to interpret and weigh the results of the energy calculations. For instance, the results of FoldX were all pretty similar and close together. We could have assumed that these results all indicate a non-neutral mutation or that they all suggest a neutral mutation. We could have even tried to distinguish between the small differences. So we considered every aspect and tried to identify a trend in the results to make a prediction about the mutation. It may be appropiate to combine the sequence- and structure-based mutation analysis to improve the confidence of the prediction. |
||
==Gromacs== |
==Gromacs== |
||
[[Image:Gla_structure_based_Time_plot.png|thumb|Figure 11: nstep vs. Elapsed Time in Gromacs.]] |
[[Image:Gla_structure_based_Time_plot.png|thumb|Figure 11: nstep vs. Elapsed Time in Gromacs.]] |
||
+ | This section shows the results of Gromacs in more detail. |
||
===Wildtype=== |
===Wildtype=== |
Latest revision as of 00:58, 8 September 2011
by Benjamin Drexler and Fabian Grandke
Contents
Introduction
In this task we analyse the structure of our protein to find out what effects the point mutations have. Therefore we created mutated structures and compared them to the wild-type protein. Several tools based on different methods have been used to achieve that aim. We used the mutations that we have chosen in the previous task.
Methods
In the first step of this task we had to find available protein structures for our protein and to decide which one would be the best for our detailed analysis. We set several cut-offs to exclude improper structures. The following tools have been used to perform the energy calulations. They were used as described in the task description.
SCWRL
SCWRL was initially developed by Dunbrack et al. in 1997. We use SCWRL4<ref name=dunb>G. G. Krivov, M. V. Shapovalov, and R. L. Dunbrack, Jr. Improved prediction of protein side-chain conformations with SCWRL4. Proteins (2009)</ref> which was published in 2009. The program takes a PDB file and a sequence file as input. By usage of a rotamer library, collision detection, and a residue interaction graph the optimal side-chain conformation is calculated, based on the backbone and the mutated sequence given in the input files. The output is a PDB file containing the conformation and the total minimal energy of the graph in STDOUT.
FoldX
FoldX was developed by Serrano et al. in 2002<ref name=serr>Guerois R, Nielsen JE, Serrano L., Predicting Changes in the Stability of Proteins and Protein Complexes: A Study of More Than 1000 Mutation. Journal of Molecular Biology (2002)</ref>. We used version FoldX 3.0 beta 4. The program provides the calculation of determination of energy effects of point mutations. It provides different run modes, but basically it takes a PDB file as input calculates several single energies(e.g. Van der Waals, Electrostatics, ...) and returns the single energies together with the total energy as output.
Minimise
Before this tool from the virtual box was used we had to remove the hydrogens and waters from the PDB file with the script repairPDB. Afterwards we were able to compare the energies differences between the wildtype and the mutated protein.
GROMACS
GROMACS is mostly a package to perform molecular dynamics, but it also provides energy calculations.
For the mutations we used the forcefield AMBER03 and for the wildtype AMBER03, AMBERGS and CHARMM27.
Additionally to the energy calculation task we did a runtime analysis with values from nsteps=10 to nsteps=1500. The results are shown in the results section of this task.
According to the task description we created an MDP file with the following content:
title = PBSA minimization in vacuum
cpp = /usr/bin/cpp
define = -DFLEXIBLE -DPOSRES
implicit_solvent = GBSA
integrator = steep
emtol = 1.0
nsteps = 500
nstenergy = 1
energygrps = System
ns_type = grid
coulombtype = cut-off
rcoulomb = 1.0
rvdw = 1.0
constraints = none
pbc = no
Keyword | Describtion<ref name=manual>Gromacs Manual</ref> |
---|---|
General | |
title | Name of Project |
cpp | Location of c-preprocessor |
Preprocessing | |
define | Defines to pass to the preprocessor; -DFLEXIBLE:include flexible water in stead of rigid water into your topology; -DPOSRES: include posre.itp into your topology, used for position restraints |
Implicit Solvent | |
implicit_solvent | Simulation with implicit solvent using the Generalized Born formalism |
Run Control | |
integrator | Steepest descent algorithm for energy minimization |
nsteps | Maximum number of steps to integrate or minimize |
Energy minimization | |
emtol | Rhe minimization is converged when the maximum force is smaller than this value |
Output | |
nstenergy | Frequency to write energies to energy file |
Tables | |
energygrps | Group(s) to write to energy file |
Neighbor searching | |
ns_type | Type of neighbor searching |
pbc | Remove the periodicity (make molecule whole again) |
Electrostatics | |
coulombtype | Type of coulomb energy |
rcoulomb | Distance for the Coulomb cut-off |
VDW | |
rvdw | distance for the LJ or Buckingham cut-off |
Bonds | |
constraints | Which constraints should be used |
Within the GROMACS work step we used the script fetchpdb. It checks if the given input is a valid PDB entry. If the check was successful it downloads the PDB file, extracts it and removes the packed version.
Results & Discussion
Structure Selection
There are several structure files available for our protein:
PDB ID | Resolution [Å] | ph-Value | R-Factor | Coverage [%] | Missing Residues |
---|---|---|---|---|---|
1R46 | 3.25 | 8.0 | 0.262 | 99.7 | 422-429 |
1R47 | 3.45 | 8.0 | 0.285 | 99.5 | 422-429 |
3GXN | 3.01 | NULL | 0.239 | 88.08 | 422-429 |
3GXP | 2.20 | NULL | 0.204 | 81.9 | 422-429 |
3GXT | 2.70 | NULL | 0.245 | 97.29 | 422-429 |
3HG2 | 2.30 | 4.6 | 0.178 | 97.32 | 422-429 |
3HG3 | 1.90 | 6.5 | 0.167 | 98.64 | 427-435 |
3HG4 | 2.30 | 4.6 | 0.166 | 99.86 | 422-429 |
3HG5 | 2.30 | 4.6 | 0.192 | 100 | 422-429 |
3LX9 | 2.04 | 6.5 | 0.178 | 98.92 | 423-435 |
3LXA | 3.04 | 6.5 | 0.216 | 99.52 | 427-435 |
3LXB | 2.85 | 6.5 | 0.227 | 99.3 | 427-435 |
3LXC | 2.35 | 6.5 | 0.186 | 98.31 | 423-435 |
We set two cutoffs to decide which structures are excluded:
- ph-value: < 6.5
- resolution: > 2.7
After we applied the cutoffs to our set of structures three were left (exclusion factors are colored red in the table). One of them was slightly better than the other ones so we decided to use 3HG3 (worse values are colored gray in the table). Additionally 3GH3 has the best overall resolution and R-factor (colored green). As the missing residues are very similar for all structures they are not further taken into account.
Visual Examination of the Mutations
Figure 1 shows the protein α-galactosidase A and the residues which will be mutated. In the following sections, we compare the side chain conformation of the mutated residues and discuss the influence of the mutation. Aspects will be, inter alia, loss of polar interactions and clashes with other residues.
SCWRL was used to model the side chain conformation of the mutated residue and we use the term tool-based to describe this side chain conformation. The side chain conformation which was done according to this tutorial is referred to as manual side chain conformation.
M42T (Mutation 1)
Figures 2 to 4 show the side chain conformation of the residue 42 in α-galactosidase A. The only difference between the manual and the tool-based side chain conformation is a variation in the conformation of the hydroxyl group of threonine (see figure 3 and 4). The tool-based side chain conformation does not lead to any clashes with the surrounding residues.
The wildtype M42 has two hydrogen bonds with E87 and Y88. Since these hydrogen bonds are formed by the carboxylgroup of the backbone, they are also abundant in the mutation, but T42 also forms a hydrogen bond with G85.
A part of T42 is exposed to the surface (see figure 5B). This is no problem, since threonine is slightly hydrophilic. But the mutation introduces a small hole into the surface of GLA (see figure 5A and 5B).
S65T (Mutation 2)
The side chain conformations of the wildtype and the mutated residues are shown in figure 6 to 8. The hydroxyl and methyl group of threonine point towards totally different direction in the tool-based conformation (see figure 6 and 8), but there are no clashes with other residues. The wildtype S65 forms five hydrogen bonds with surrounding residues (C63, K67, L68, F69). The mutated residue T65 forms also five hydrogen bonds, but one of them is with E66 instead of K67.
Serine is hydrophilic and the residue is exposed on the surface, but threonine is also a hydrophilic residue. There is no remarkable change of the surface of the protein.
I117S (Mutation 3)
Since the tool-based side chain conformation of the mutated residue (see figure 11) is pretty similar to the side chain conformation of the wildtype (see figure 9), there are no clashes in the mutated structure. The angle of the hydroxyl group is slightly different in the manual side chain conformation (see figure 10).
The carboxyl group in the backbone of I117 forms two hydrogen bonds with L120 and A121. These bonds are also abundant in the mutated structure. Serine is a hydrophilic amino acid, but this is no problem, because it is part of the surface. Overall there is no remarkable difference in the surface of the protein due to the mutation.
A143T (Mutation 4)
This is the only mutation which is close to the active site of the protein (see figure 12). The general orientation of the side chain is very similar between the wildtype, manual and tool confirmation (see figure 13, 14 and 15). The hydroxyl group of threonine leads to an increase of the space requirements with the tool-based side chain conformation. Hence there is a clash between the hydroxyl group of T143 and the residue D93 (see figure 16).
A143 forms two hydrogen bonds with the carboxyl group of its backbone to D92 and T141. These hydrogen bonds are also formed by T143, but a third hydrogen bond is established between the oxygen of the hydroxyl group of T143 and D92.
T143 is part of the surface, so it is no problem that this position becomes hydrophilic. The mutation changes the surface and since this mutation is nearby the active site, this will probably have an influence on the work rate of the protein (see figure 17). The entry to the active site becomes more narrow due to the mutation. Hence it will be less likely or even impossible that the ligand binds in the active site.
H186R (Mutation 5)
Figures 13 to 15 show the side chain conformations of the residue. The general orientation of the side chain is similar. Because of this, the residue points towards the exterior of the protein. There is small collision with the residue D153 which is located in alpha-helix nearby (see figure 16).
H186 forms three hydrogen bonds with its backbone to D182, L189 and A190, which are also formed by the mutated residue R186. Both amino acids, histidine and arginine, are hydrophobic even though the position itself is exposed on the surface. The change on the surface is not noteworthy.
P205T (Mutation 6)
Since the structure of proline is a special case, the comparison of the side chain conformation is not trivial. Threonine has a higher space requirement than proline. The difference between the manual and tool-based side chain conformation is a variation of one angle in the side chain. There are two small collisions with surrounding residues (I219 and N228) due to the higher space requirements of threonine (see figure 25).
P205 forms two hydrogen bond with M208 and W209. The bond with M208 is also formed by the mutated residue T205, but the hydrogen bond to W209 is lost. Therefor a hydrogen bond to N228 is established. The substitution of proline to threonine does not have any influence on the surface of the protein.
D244H (Mutation 7)
The side chain of the wildtype has a very flat angle to the alpha-helix (see figure 26), whereas the manual and tool-based side chain conformations show an almost orthogonal orientation (see figure 27 and 28). Hence there are no collisions with other residues. The manual and tool-based conformations are also very similar to each other.
D244 forms four hydrogen bonds. The carboxyl group of the side chain forms two bonds with arginine at position 356. The remaining two bonds are formed by the backbone of D244 with K240 and F248. These two hydrogen bonds are also established by the mutated residue H244 and the side chain forms a third bond with S247.
This position is exposed on the surface of the protein and aspartic acid and histidine are hydrophilic amino acids. The orthogonal angle of histidine to the alpha-helix adds a bulge and a small hole to the surface of the protein (see figure 29).
Q283P (Mutation 8)
Figures 30 to 32 show the side chain conformation at position 283. The manual and the tool-based side chain conformation are very similar (see figure 31 and 32). This is probably due to the fact that proline is a very rigid amino acid. Glutamine has much higher space requirements. Hence the substitution to proline does not lead to a collision with other residues.
The hydroxyl group of the side chain of Q238 forms two hydrogen bonds with M267 and V268. The backbone establishs another four hydrogen bonds with Q278, Q279, L286 and W287. These four residues are located in the same alpha-helix. The mutated residue P283 forms only two of these four hydrogen bonds with L286 and W287. This might result in partial loss of the alpha-helix which is no surprise, since proline is considered to be a helix breaker<ref name=proline_helix>Gunasekaran et al., "Stereochemical punctuation marks in protein structures: glycine and proline containing helix stop signals.". J Mol Biol. 1998 Feb 6. PubMed</ref>.
Q283 is part of the core, so the substitution does not change the surface.
Q321E (Mutation 9)
Since the structures of glutamine and glutamic acid are so similar, the side chain conformation shows almost no differences (see figure 33, 34 and 35). There is only one angle in the side chain which is slightly different that does not lead to a collision with other residues.
The oxygen in the side chain of Q321 forms a hydrogen bond with T39, the amine forms two hydrogen bonds with P40 and N320. The backbone of Q321 establishs a hydrogen bond with I317 which is also formed by the mutation E321. The side chain of the mutation E321 only establishs one hydrogen bond with T39. So there are two bonds missing due to the mutation.
Even though this position is exposed on the surface, there are no remarkable changes due to the substitution.
R363C (Mutation 10)
The side chain confirmations of this position are shown in figures 36 to 38. The conformation between the manual and the tool-based ones is very similar. Arginine has a much higher space requirment, hence there is no collusion with other residues due to the substitution to cysteine.
The side chain of R363 forms five hydrogen bonds with D335, N336 and N355. These bonds are all lost due to the substitution to cysteine. The backbone of R363 establishs another two hydrogen bonds with I407 which are preserved in the mutation C363. So overall there is a loss of five hydrogen bonds.
Arginine is a hydrophilic amino acid and R363 is exposed on the surface. The substitution to the hydrophobic amino acid cysteine could be a problem. The mutation also introduces a big hole in the surface (see figure 39).
Energy Comparison
The results of the energy comparison are presented in the table below. Due to the fact that the result of minimise for mutation 8 clearly differed from the other results, the run was repeated with the outcome from the first run as input. Thus, there is the number 8.2 in the table.
Number | AA-Position | Amino acid change | SCWRL4 | FoldX | Minimise | Gromacs | |||
---|---|---|---|---|---|---|---|---|---|
Energy | Difference | Energy | Difference | Energy | Difference | ||||
Wildtype | - | -20.93 | - | -20481.23 | - | -91307.7 | |||
1 | 42 | Met -> Thr | 343.25 | 157.29 | -178.22 | -20324.41 | -156.82 | -90528.4 | -779.3 |
2 | 65 | Ser -> Thr | 327.798 | 152.87 | -173.8 | -20339.34 | -141.89 | -90481.9 | -825.8 |
3 | 117 | Ile -> Ser | 333.027 | 157.97 | -178.9 | -20353.47 | -127.76 | -90654 | -653.7 |
4 | 143 | Ala -> Thr | 333.944 | 154.40 | -175.33 | -20339.32 | -141.91 | -90541 | -766.7 |
5 | 186 | His -> Arg | 323.717 | 154.57 | -175.5 | -20321.32 | -159.91 | -91011.7 | -296 |
6 | 205 | Pro -> Thr | 340.619 | 155.96 | -176.89 | -20345.87 | -135.36 | -90782.2 | -524.8 |
7 | 244 | Asp -> His | 333.594 | 152.08 | -173.01 | -20393.12 | -88.11 | -90232.9 | -1074.8 |
8 | 283 | Gln -> Pro | 332.631 | 159.91 | -180.84 | -8027.71 | -12453.52 | -87316 | -3991.7 |
8.2 | - | - | - | - | - | -19134.48 | -1346.95 | - | - |
9 | 321 | Gln -> Glu | 332.853 | 160.95 | -181.88 | -20246.98 | -234.25 | -90090.3 | -1217.4 |
10 | 363 | Arg -> Cys | 330.56 | 150.50 | -171.43 | -20295.77 | -185.46 | -89721.8 | -1585.9 |
FoldX
There is no mutation that stands out with a very high or very low energy. Every energy difference ranges from approximately -170 to -185. Hence it is very tough to make a clear distinction between neutral and non-neutral mutations.
Minimise
The results of Minimise cover a much wider range than the results of FoldX. This is mainly due to the result of mutation 8 which shows a significantly lower result than the others. Mutation 9 and 10 should also be mentioned as mutations with a very high difference to the wildtype. Mutation 7 has the lowest difference and the other mutations are somewhere in between.
Gromacs
Mutation 5 is the mutation with the lowest difference to the wildtype. Once again, mutation 8 has the highest difference. The next least stable structures are the ones of mutation 10, 9 and 7. Mutation 5 has the lowest difference
Individual Evaluation of the Mutations
We evaluate the results of the visual examination and the energy comparison in the following sections individual for each mutation and try to predict the influence of the mutation on the protein.
M42T (Mutation 1)
The structural influence of the substitution from methionine to threonine at this position can be reduced to the change of the surface in this case, since this is the only aspect in which a remarkable difference occurs. Another aspect is the loss of disulfide bonds due to the mutation which could lead to a loss of stability of the protein. The UniProt entry does not contain a disulfid bond at this position. The calculated energy differences show a slight tendency towards a non-neutral mutation. Because of this and the remarkable change of the proteins surface, we assume that his mutation is non-neutral.
Prediction: a non-neutral mutation.
S65T (Mutation 2)
The side chain conformation of the mutated residue is slightly different, but there a no clashes with other residues and the hydrogen bonds are preserved or substituted, respectively. The calculated energy differences are among the lowest and since there is no significant structual change, we assume a neutral mutation.
Prediction: a neutral mutation.
I117S (Mutation 3)
The substitution of isoleucine to serine does not lead to a remarkable change in any of the visual examined aspects. The energy differences of Gromacs and Minimise are on the lower end of all mutations. Because there is no aspect that indicates a influence on the functionality of the protein, we predict a neutral mutation.
Prediction: a neutral mutation.
A143T (Mutation 4)
This mutation has a huge influecne on the structure of the protein. There are two aspects that show a remarkable change. First, the mutations leads to a clash with another residue. Second, there is a change of the surface of the protein which is the major aspect. Since this mutation is nearby the active site, the ligand is unlikely to be able to bind at the active site and this will probably have a strong influence on the functionality of the protein. The energy comparison does not show such a strong tendency to a non-neutral mutation, but the observed changes do not lead neccessarily to a change of the energy. Nevertheless, this mutation is very likely to be a non-neutral mutation.
Prediction: a non-neutral mutation.
H186R (Mutation 5)
Even though the substitution to arginine leads to clash with another residue at this position, the energy values does not indicate an extraordinary result in relation to the other mutations. Quite the contrary, the value of Gromacs is one of the lowest. Since the collision is very small, the hydrogen bonds of this position are presvered and the calculated energy differences are quite low, we assume that this mutation is neutral.
Prediction: a neutral mutation.
P205T (Mutation 6)
Even though the substitution of proline to threonine leads to collisions with two other residues, the energy comparison does not indicate a less stable protein in relation to the energies of the other mutations. This mutation has the third or second lowest difference in the results of Minimise or Gromacs, respectively. Since the energy comparison indicates a neutral mutation and the collisions are very slight, we assume that this mutation is neutral.
Prediction: a neutral mutation.
D244H (Mutation 7)
Overall there is a loss of hydrogen bond due to this mutation. The energy results are quite controversial, because it has one of the highest differences in the results of Gromacs, but the lowest in the results of Minimise. It is also on the lower end in the results of the FoldX, but its result are not that meaningful.
The mutation also changes the surface of the protein, but it is tough to evaluate the influence on the functionality of the protein. Since the energy results are leaning towards a neutral mutation, we assume that this mutation is neutral.
Prediction: a neutral mutation.
Q283P (Mutation 8)
It is very likely that the substitution to proline leads to a partial loss of the alpha-helix in this area. This is supported by the results of the energy calculations which show the significantly highest difference of all mutations. We predict that this mutation is non-neutral, since the structual influence of this mutation is probably very high.
Prediction: a non-neutral mutation.
Q321E (Mutation 9)
Even though the glutamic acid is very similar to glutamine and the mutation has the same conformation, the engery results strongly indicate a non-neutral mutation. The difference of this mutation is among the highest for all thre programs. The two missing hydrogen bonds are probably partially responsible for this. Because the calculated energies show such a clear result, we predict a non-neutral mutation.
Preidction: a non-neutral mutation.
R363C (Mutation 10)
The mutation leads to a loss of five hydrogen bonds. This could result in a loss of stability of the protein and it could explain the results of the energy calculation which show high differences (Minimise and Gromacs). The substitution also introduces remarkable changes on the surface.
Since the structual examination and the energy calculations indicate a strong influence on the functionality of the protein, we assume that this mutation is non-neutral.
Prediction: a non-neutral mutation.
Overall Evaluation
Overview
The correct assignment (solution) and the predictions of the ten mutations are listed in the table below.
We applied the following coloring scheme:
- green: the prediction was correct.
- red: the prediction was wrong.
Number | Mutation | Solution | Prediction |
---|---|---|---|
1 | M42T | non-neutral | non-neutral |
2 | S65T | neutral | neutral |
3 | I117S | neutral | neutral |
4 | A143T | non-neutral | non-neutral |
5 | H186R | neutral | neutral |
6 | P205T | non-neutral | neutral |
7 | D244H | non-neutral | neutral |
8 | Q238P | non-neutral | non-neutral |
9 | Q321E | non-neutral | non-neutral |
10 | R363C | non-neutral | non-neutral |
Mutation 6: P205T
This mutation seemed to be non-neutral because of a collision with another residue, but the energy calculations clearly indicated a neutral mutation. Hence, we may have overweighed the energy calculation of the programs.
Mutation 7: D244H
The solution is quite surprising, because the energy difference was the lowest in the results of Minimise. Otherwise we would have probably predicted a non-neutral mutation.
Discussion of Structure-based Mutation Analysis
At first, we thought it would be easier to make a prediction about the mutation than the sequence-based analysis, because we figured that we would be able to see the structual influence of the mutation more clearly. On the one hand, it was tough to evaluate the influence of the individual changes like change of the surface or loss of hydrogen bonds. On the other hand, it was also hard how to interpret and weigh the results of the energy calculations. For instance, the results of FoldX were all pretty similar and close together. We could have assumed that these results all indicate a non-neutral mutation or that they all suggest a neutral mutation. We could have even tried to distinguish between the small differences. So we considered every aspect and tried to identify a trend in the results to make a prediction about the mutation. It may be appropiate to combine the sequence- and structure-based mutation analysis to improve the confidence of the prediction.
Gromacs
This section shows the results of Gromacs in more detail.
Wildtype
Force Field | Average | Error Estimat | RMSD | Tot-Drift (kJ/mol) |
---|---|---|---|---|
Bond | ||||
AMBERGS | 1826.99 | 420 | 4409.39 | -2499.37 |
AMBER03 | 1639.74 | 410 | 4358.68 | -2424.42 |
CHARMM27 | 2908.14 | 350 | 4779.8 | -2033.44 |
Angle | ||||
AMBERGS | 5496.47 | 74 | 476.18 | 408.548 |
AMBER03 | 5324.13 | 72 | 469.75 | 369.24 |
CHARMM27 | 7975.2 | 86 | 798.12 | 432.901 |
Potential | ||||
AMBERGS | -114713 | 1200 | 5648.79 | -7915.46 |
AMBER03 | -91307.7 | 1200 | 5559.82 | -7839.05 |
CHARMM27 | 136.699 | 32 | 64.3892 | 227.896 |
Mutations
Force Field | Average | Error Estimat | RMSD | Tot/Drift |
---|---|---|---|---|
Bond | ||||
1 | 1815.39 | 570 | 5166.85 | -3384.48 |
2 | 1862.77 | 610 | 5331.85 | -3618.04 |
3 | 1773.13 | 520 | 4937.34 | -3068.93 |
4 | 1828.63 | 580 | 5229.18 | -3479.09 |
5 | 1870.95 | 610 | 5361.67 | -3713.22 |
6 | 1816.6 | 550 | 5091.81 | -3303.34 |
7 | 1819.7 | 570 | 5173.34 | -3397.07 |
8 | 2992.15 | 1700 | -nan | -10631.8 |
9 | 2083.16 | 830 | -nan | -4913.82 |
10 | 1867.42 | 620 | 5390.82 | -3693.03 |
Angle | ||||
1 | 5183.95 | 85 | 360.959 | 550.303 |
2 | 5195.33 | 80 | 364.473 | 515.645 |
3 | 5196.5 | 89 | 353.256 | 586.473 |
4 | 5175.59 | 85 | 364.496 | 547.465 |
5 | 5113.99 | 81 | 365.511 | 526.244 |
6 | 5200.44 | 85 | 356.964 | 553.934 |
7 | 5261.77 | 87 | 365.202 | 565.196 |
8 | 5178.73 | 76 | -nan | 215.036 |
9 | 5201.95 | 76 | -nan | 442.141 |
10 | 5174.48 | 88 | 375.775 | 555.294 |
Potential | ||||
1 | -90528.4 | 1600 | 7234.09 | -10149.1 |
2 | -90481.9 | 1600 | 7442.03 | -10340 |
3 | -90654 | 1500 | 6928.73 | -9614.54 |
4 | -90541 | 1600 | 7311.04 | -10343.7 |
5 | -91011.7 | 1600 | 7484.45 | -10592.5 |
6 | -90782.2 | 1600 | 7226.99 | -10188.5 |
7 | -90232.9 | 1600 | 7236.24 | -10198 |
8 | -87316 | 3600 | -nan | -23670.3 |
9 | -90090.3 | 1900 | -nan | -12335.3 |
10 | -89721.8 | 1700 | 7523.88 | -10750.1 |
Mutation | Plot |
---|---|
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 |
References
<references />