Wild type MgCYP51 and a selection of MgCYP51 variants were modelled. Amino acid alterations currently observed in the field were modelled in order to elucidate the contribution of single and multiple alterations on the resistance and activity profile of MgCYP51. We used the PocketFinder (http://www.modelling.leeds.ac.uk/pocketfinder/
) algorithm to calculate the predicted size of the pocket volumes of our 14 models and the root mean squared deviation (RMSD) as a measure of the overall change in the conformation of our models compared to the wild type model. Together these give a measurement of the effect an alteration or combination of alterations is likely to cause on the structure of MgCYP51 (). In addition, we have used molecular modelling to simulate the docking of epoxiconazole, prochloraz, tebuconazole and triadimenol in wild type and mutant models to gain an insight into the mechanism behind the associated azole resistance. We have considered the effect that individual alterations may contribute and the overall profiles expected in variants carrying multiple changes currently observed in the field.
Binding pocket volumes of modelled M. graminicola CYP51 mutants, and their RMSD values compared to wild type, ordered according to cavity volume.
Wild type M. graminicola CYP51
CYP51 is predicted to possess a typical P450 fold with the generally conserved structural core formed by helices E, I, J, K, and L around the haem prosthetic group (), similar to the recently characterised human CYP51 
, and the access to the haem cavity guarded by the loop between helices B and C. The predicted volume of the wild type protein is 1685 Å3
, including the space occupied by the haem. A total of 86 residues were identified as lining the pocket of the wild type protein (text S1
), including 5 residues (V136 (on the B–C loop), I381 (K–L loop), Y459, G460 and Y461) observed as sites of alteration leading to azole resistance in M. graminicola 
. shows the location in the wild type protein structure of the residues subject to alteration. V136 and Y137 reside on the access channel end of the binding pocket, while I381, Y459, G460 and Y461 are at the haem end. L50 and A379 do not directly border the binding pocket, but form part of secondary structures that are immediately adjacent to the cavity. S188 and N513 are highly exposed on the outside of the protein, far removed from the binding regions.
The numbers of residues bordering the binding pocket and predicted to be within 3 Å of epoxiconazole, prochloraz, tebuconazole and triadimenol are shown in (the residues predicted to be within 3.0–3.5 Å and 3.5–4.5 Å of docked azoles are listed in table S1
). The docking of azoles to the wild type MgCYP51 reveals a pattern of residues immediately proximal to the azole docking site known to be subject to mutation which results in azole resistance (V136, Y137, I381, Y459, G460, and Y461). All these residues are implicated, for instance, by the interactions of triadimenol docking (closer than 4.5 Å). By this analysis, L50, S188, A379 and N513 are confirmed to be located away from the binding region. This modelling is consistent with the use of azoles leading to the emergence of the key resistance mutations at Y137, I381V and Y459.
Impact of alterations
The effect of individual alterations upon CYP51 structure was assessed by modelling the altered forms of the protein. Test sequences were modified accordingly and input into the modelling pipeline. details the changes in homologues adopted by the automated modelling pipeline and changes in sequence homology, which lead to the structural changes observed in the modelled structures. The proximity of residues to each of the azoles for which binding was simulated was assessed and compared with the wild type protein.
For the purpose of discussion of the impacts of alterations, the variants are divided into three groups; i) Alterations that do not affect azole binding, ii) Alterations affecting the size and accessibility of the binding pocket, iii) Alterations affecting the localisation of residues interacting with azoles.
i) Alterations that do not affect azole binding
The inclusion of L50S on its own causes only a slight increase in cavity volume (). The L50S substitution appears in the majority of subsequent combinations of multiple mutations and its location away from the cysteine pocket (L/S50 is not within 4.5 Å of any of the docked azoles in any of our 14 models (table S1
) and its marginal impact on protein conformation indicate that it may indeed represent a natural variant in its own right which is not associated with azole resistance. The L50S alteration also does not affect activity as shown by expression in yeast 
L50S S188N N513K
L50S S188N N513K is a natural variant found in wild type populations in New Zealand sensitive to all azoles (Cools et al., unpublished). Although the model shows a significant increase in the volume of the haem cavity (), and generally fewer residues in close proximity to the azoles, key residues for interaction, namely Y137 and I381 are as close or closer to the docked azoles than in the wild type protein (table S1
), being within 3.5 Å of all azoles, an arrangement consistent with general azole sensitivity. V136 is within 4.5 Å of triadimenol, tebuconazole and epoxiconazole, but not prochloraz, also like the wild type protein.
Leroux et al., (2007) have suggested that the mutations recorded at positions 50, 188, 379 and 513 do not seem to be related to DMI resistance, but suggest that M. graminicola
populations might be constituted of two independent entities that represent sympatric species 
. They suggest that S188N and N513K (along with L50S) are always seen together and may have originated from the same haplotype rather than having been selected through conferring azole resistance 
. Brunner et al 
have proposed an alternative hypothesis, that rapid evolution of the CYP51 gene has led to two distinct populations of molecules, one based on the accumulation of a series of point mutations within the ancestral wild-type template and another that is highly diverged as a result of a few intragenic recombination events that brought together the different point mutations from different regions of the gene.
This genotype may exhibit general sensitivity to the above azoles and might present an increased selective pressure for mutation at the 136 and 381 positions. This is supported by the observation that the L50S variant model brings I381 into closer interaction with triadimenol (from within 4.5 to <3.0 Å), tebuconazole and epoxiconazole (now within 4.5 Å) and remains within 3.0 Å of prochloraz (table S1
). There is a reduction in the overall number of residues within 4.5 Å of triadimenol (13 to 9) and prochloraz (11 to 10), but an increase in those within 4.5 Å of tebuconazole (12 to 14). There are also more polar residues within 4.5 Å of tebuconazole and epoxiconazole ().
ii) Alterations affecting the size and accessibility of the binding pocket
The predicted structure of the L50S I381V Y461H variant displays a notable constriction of the binding pocket, while the cavity of L50S S188N N513K is more than double the size of the wild type pocket, suggesting that the slight increase in cavity volume brought about by L50S on its own is modified considerably by the other alterations with which it is combined.
Variants including I381V
Several of the multiple alterations modelled are seen in current populations. The structural models will be discussed in relation to the resistance profiles observed and compared to wild type structure to account for resistance, beginning with those that include the influential I381V alteration.
Multiple alterations including I381V display some sensitivity to prochloraz. L50S I381V Y461H is a common natural variant, which causes resistance to epoxiconazole and slightly decreased sensitivity to prochloraz (). The predicted cavity volume of this mutant protein model is notably smaller than wild type at 1217 Å3 and therefore the number of residues within interaction range of docked azoles is increased ().
EC50 and resistance factor values of CYP51 variants and structural changes conferring resistance.
It is the incorporation of I381V that causes the decrease in the cavity volume. The single I381V substitution reduces the cavity volume to 1438 Å3
, and similarly causes a relative decrease in the cavity volume in other mutants where it appears (). Interestingly, L50S I381V and I381V have never been observed in populations in isolation, and do not complement CYP51 function in a heterologous system 
, indicating that the conformational change also results in decreased activity. Leroux and colleagues have observed high resistance of this laboratory variant to triadimenol, tebuconazole and epoxiconazole 
I381V is usually found in populations with changes in the 459–461 region 
, indicating that CYP51 function is lost with I381V and this mutation is only selected after other alterations have occurred. Indeed, it has been proposed that a recombination event brought these two mutations together 
. It has also been shown in the heterologous system that the combination of I381V and Y459–Y461 rescues the lethality of I381V alone 
I381/V381 is within 4.5 Å of docked azoles in 36 of our 56 docking simulation models (table S2
) and so is clearly under substantial evolutionary pressure. The I381V substitution results in the entire 459–461 section and Y137 being worked closer to the docked azole (, table S1
) – substantial structural changes which may account for the lack of CYP51 activity seen and a consequent lack of selective pressure without the presence of other mutations.
L50S I381V Y461H is a common variant, resistant to tebuconazole and sensitive to prochloraz. No strict pattern of individual residues within interaction range of docked azoles can account for resistance, in fact Y137 is brought to within 3 Å of all azoles including prochloraz and the proximity of V136 with triadimenol, tebuconazole and epoxiconazole is maintained. It may be that the substitution of tyrosine by histidine at position 461 in itself may cause azole resistance, though this alone does not explain why L50S I381V Y461H causes further decreases in sensitivity to tebuconazole and epoxiconazole compared to L50S Y461H. Proximity to V381, G460 and S315 also appears to be related to observed resistance. All are out of range of tebuconazole, while V381 is within 3 Å and G460 is within 3.5 Å of prochloraz. This fits in with the observed EC50 and resistance factor values (), which are markedly higher for tebuconazole than for prochloraz.
Therefore, here it appears that it is the general constriction of the cavity combined with removal of structure at the 460–461 end that may be responsible for the observed resistance to tebuconazole and epoxiconazole, while prochloraz gives a different profile, lying closely to V381 and G460. This protein also conferred impaired growth in a heterologous system 
, indicating that the conformational change also results in decreased activity.
Our modelling suggests that the main effect of the I381V alteration appears to have been to reduce the interactions of triadimenol, tebuconazole and epoxiconazole by constriction of the haem cavity at the K and L helix end, and subsequent removal from interaction range of key tyrosine and polar residues elsewhere around the cavity, rather than removal of the 381 position itself. In the wild type model, I381 is within 3 Å of prochloraz and within 4.5 Å of triadimenol, but out of range of epoxiconazole and tebuconazole. In the L50S model, I381 is generally around 1.5 Å closer to the docked azoles, being within 3 Å of prochloraz and triadimenol and within 4.5 Å of epoxiconazole and tebuconazole (table S1
). I381V has not been observed in the field without L50S. Even following multiple alterations including I381V, the 381 position remains in close range of prochloraz, though not the other azoles, which appears to be the strongest factor in the retention of sensitivity to prochloraz.
The ΔY459/G460 deletion
Alterations of residues 459 and 460 of MgCYP51 occur frequently. Docking of azoles in our models reveals these residues are within 4.5 Å in 34 of the 40 docking models that do not include the ΔY459/G460 deletion, indicating that these amino acids are important in azole binding. The residues 459–461 lie within a region (438–463) which is specific to fungal CYP51s. It is this region which renders fungal models of MgCYP51 that use only a single homologue as a template unreliable. Cañas-Gutiérrez et al. used the Myocbacterium tuberculosis
CYP51 structure as a template for an M. fijiensis
model and found the area between position 442–472 (their alignment) resulted in an energetic difference between the model and the scaffold 
. In contrast, our method is not constrained to a single homologue and our results below show how important this region is in relation to azole resistance, particularly in the ΔY459/G460 mutant. Indeed, this is the first fungal CYP51 model that provides structural evidence for the importance of Y459–Y461 for azole binding.
The incorporation of the Y459/G460 deletion is consistently associated with large increases in cavity volume. A deletion at Y459/G460 is predicted to increase the cavity volume to 4755 Å3 and in combination with other mutations also increases the cavity volume (L50S S188N ΔY459/G460 N513K (5725 Å3) compared to L50S S188N N513K (3754 Å3). This increase in volume is to be expected, as the deletion results in the removal of that whole section of beta turn from the vicinity of the haem pocket (). The removal of this region out of the haem cavity to the outside of the protein is a feature of all mutants carrying the ΔY459/G460. Variants carrying the deletion also display some of the highest RMSD values compared to wild type.
Azole docking in mutant proteins.
ΔY459/G460 consequently brings about a dramatic effect on the proximity of key residues to each of the azoles docked. Y137, Y461, and of course the deleted Y459 and G460, are entirely removed from the binding pocket whilst the lysine residues K133, K148 and K149 move into the pocket (). These changes are brought about by gross rearrangement of the two sections involved, namely that containing K133, Y137, K148 and K149 and the beta turn section 450–470 containing residues 459–461 (). In B. graminis f. sp. hordei
the combination of the Y136F and K147Q substitution displays additive effects towards triadimenol resistance 
The loss of interaction with Y137 and Y459 balanced by the introduction of the two lysine residues within range for interaction is in good keeping with the observation of medium level resistance to all the azoles tested 
. Interestingly this alteration has been predicted to have arisen in two individual events by Brunner et al. 
, further demonstrating the importance of this deletion in the development of higher level azole resistance. This is the first model to show the considerable effect of deletion of 459/460.
L50S S188N ΔY459/G460 N513K and L50S S188N I381V ΔY459/G460 N513K
The largest change in pocket volume was found in the natural variant L50S S188N ΔY459/G460 N513K (5724 Å3). All the mutated residues are withdrawn from the azole binding region, except for I381 which remains within 3.5 Å and 3 Å of triadimenol and prochloraz respectively. Q313 is brought within range of tebuconazole, and H314 within reach of all the azoles. Y123 is within 3 Å of triadimenol, tebuconazole and epoxiconazole and within 4.5 Å of prochloraz. This variant demonstrates substantial resistance to epoxiconazole and tebuconazole and more modest resistance to prochloraz (), which suggests that proximity to I381 is a factor in determining the extent of resistance to a given azole, but particularly prochloraz.
For the natural variant, L50S S188N I381V ΔY459/G460 N513K, the model showed the haem cavity volume increased to 5204 Å3
and other than V381 (within 3.5 Å of triadimenol and prochloraz), all mutated residues are removed from the binding pocket. However Y123 is within 3 Å of all azoles - this is the only polar residue found in interaction range of all docked azoles. Isolates carrying this variant are resistant to tebuconazole and sensitive to prochloraz 
, our models suggesting again that the proximity of the 381 position to prochloraz is the determining factor in the differential sensitivities observed.
L50S S188N ΔY459/G460 N513K and L50S S188N I381V ΔY459/G460 N513K may reduce sensitivity to some azoles in a similar way since Y137 and lysine residues (K133, K148 and K149) are all more than 4.5 Å from docked azoles.
iii) Alterations affecting the localisation of residues interacting with azoles
The Y137F substitution brings about a substantial increase in the size of the pocket to over twice that of the wild type (3769 Å3
, ). This is the largest change in cavity volume brought about by any of the single alterations modelled. However, Y137F does not exert its effect by reducing the number of polar residues in proximity to the bound azoles () and it also does not remove the altered residue from the pocket. In the model F137 is closer to the docked azoles than Y137 (within 4.5 Å of all azoles and within <3 Å of triadimenol, table S1
) and is pushed into an obstructive position prohibiting the binding of triadimenol (). The docking of the other azoles within the haem cavity is not constrained by the occlusion of the pocket around the 137 position; rather it is the loss of the hydroxyl group incurred by the substitution of tyrosine with phenylalanine which may result in the minor decrease in sensitivity to tebuconazole, epoxiconazole and prochloraz observed. The conformational change caused by the substitution also serves to take Y459 away (around 2 Å) from the azole binding position. Y459 is within 3 Å of all docked azoles in the wild type, but in the mutant is only within 4.5 Å of triadimenol and is moved out of interaction range (>4.5 Å) of tebuconazole, epoxiconazole and prochloraz (table S1
These modelling observations are consistent with resistance recorded in M. graminicola
and other fungi with mutations of equivalent residues. The mutations Y132H in Candida albicans 
, Y136F in M. fijiensis 
, Y136F in Erysiphe gramainis f.sp. hordei 
and Y136F in Uncinula necator 
and have all been associated with azole resistance. Y137F in M. graminicola
is specifically associated with triadimenol resistance 
. Strains carrying the Y137F substitution are believed to have emerged due to triadimenol usage 
and although the mutation was common in the early 1990s it is now rare in most European populations 
, having been replaced by mutations in the 459–461 region. In addition the large conformational change predicted by our model is further supported by the fact that the activity of MgCYP51 Y137F is greatly reduced, to around 10% of that of the wild type 
Y459D, G460D, Y461H and Y461S alterations
Although single substitutions at positions 459, 460 and 461 all lead to an increase in cavity volume to varying extents (), indeed G460D leads to a more than doubling of the cavity volume (to 3501 Å3), the main effects of these single substitutions are with regard to the localisation of azole-interacting residues.
With the addition of the substitution Y459C in L50S Y459C, a rare variant, Y137 is withdrawn from the pocket (>4.5 Å from all docked azoles) while I381 and Y461 remain close to the azole. Y461 is predicted to be brought closer (<3.0 Å) than in L50S or the wild type. By substitution of Y459 to cysteine, the azole proximity at that position widens when triadimaenol and tebuconazole are tested (<3 to <3.5 Å and <3 to <4.5 Å respectively).
Our modelling of the Y459D substitution shows that V136 and Y137 are withdrawn from the azole binding pocket (>4.5 Å from all docked azoles) whilst I381 is brought closer to bound triadimenol (<3 Å) and tebuconazole (<4.5 Å).
In contrast the G460D substitution shows that azoles maintain close contact with both Y137 and Y459, however, V136 is removed from interaction range (>4.5 Å). The mutated D460 residue itself is pushed into close proximity with triadimenol, tebuconazole and epoxiconazole (<3 Å).
With the Y461H substitution, the most notable difference is at the site of the mutated residue where the azole interaction with Y461 is lost due to its change to histidine. Close proximity is maintained between triadimenol and I381, indeed markedly closer proximity (<3.0 Å compared to <4.5 Å) than for the wild type (table S1
). The inclusion of Y461H in L50S Y461H results in residue Y123 being closer to all the docked azoles (<3 Å) whilst Y137 and Y459 also remain close.
L50S Y461S unsurprisingly provides a similar profile to L50S Y461H - since serine has a smaller side chain than tyrosine, S461 is further from the tested azoles than Y461 >4.5 Å for all docked azoles). However, unlike Y461H, Y137 is withdrawn from the azole binding region (>4.5 Å from all docked azoles).
In summary, substitutions at positions 459–461 cause azole resistance by moving residues V136 and or Y137 further from the docked azoles, whereas the biggest impact on azole binding of alteration at Y461 is likely to be due to the loss of the tyrosine residue. An inference from our structural modelling is that the agricultural azoles impose uniformly intense selective pressure at residues between 459–461 forcing the organism toward the adoption of an evolutionary strategy of more drastic deletion rather than single substitutions in that region. This idea is supported by the fact that ΔY459/G460 is present in the most prevalent genotypes recently observed.
L50S V136A Y461H
V136A has not been found alone in M. graminicola
. Similar to I381V, when introduced alone and expressed in yeast it is lethal and can be partially rescued by combining with Y459–Y461. However, L50S V136A Y461H is a quite common variant that causes resistance to prochloraz and it has been postulated that the V136A substitution has been selected by the use of this imidazole 
. Although V136 is not within 4.5 Å of prochloraz in any of our mutants, our model agrees with this finding. Rather than prochloraz resistance being directly via a change in the interacting residue it is proposed that a conformational change results in Y137 and Y123 moving further from docked prochloraz that reduces affinity.
The haem cavity volume of the natural variant, L50S V136A Y461H, is only slightly increased to 1971 Å3 but more residues are within 3 Å of docked azoles than in the wild type. In particular, Y137 and residue 461 are closer to all docked azoles, though for prochloraz, Y137 is more than 4 Å away. However, V136 proximity with triadimenol, tebuconazole and epoxiconazole is lost. This variant causes resistance to prochloraz, far greater than with Y461H alone. Compared to L50S Y461H, the main difference between the two model structures is that I381 is about 1 Å further away in L50S V136A Y461H.
V136 is calculated as being closer than 4.5 Å to the docked azole in 15 of our 56 azole docking models, but not for prochloraz binding (table S1
). In L50S V136A Y461H, residues 459–461 are closer to the docked azoles than in L50S Y461H.
With L50S V136A S188N ΔY459/G460 N513K, the haem cavity volume is increased to 4355 Å3 causing a great reduction in the number of residues within 3 Å of docked azoles (). K148 is within 3 Å of tebuconazole, but further away (3.5–4.5 Å) from the other azoles, which correlates with the observed sensitivity to tebuconazole and conversely the observed resistance to epoxiconazole and prochloraz (, ).
Inclusion of A379G
The inclusion of A379G in L50S S188N A379G I381V ΔY459/G460 N513K decreases the cavity volume only slightly to 5060 Å3
from 5204 Å3
(compared to L50S S188N I381V ΔY459/G460 N513K, ). In this variant, Y137 is around 4 Å from epoxiconazole and tebuconazole (), which gives some scope for resistance with only the weakest form of H-bond being possible. Indeed, both V136 and Y137 are within interaction range of triadimenol, tebuconazole () and epoxiconazole, but the conspicuous proximity of V381 to prochloraz alone is repeated. V381 is around 2 Å closer to prochloraz than the other azoles (). Stammler et al., (2008) 
have found that this construct is highly resistant against epoxiconazole and tebuconazole, but sensitive to prochloraz (). It is interesting that the incorporation of the A379G alteration has little impact upon the binding capacity of the nearby 381 position (table S1
), which can be associated with interaction with prochloraz and consequent sensitivity (). The more compact structure of the valine side chain might allow richer hydrophobic interactions with prochloraz than the wild type isoleucine side chain accounting for the sensitivity seen to prochloraz in this strain.
The presence of I381V in combination with other alterations increases following the introduction of tebuconazole and decreases after prochloraz, with no difference following epoxiconazole 
. This sensitivity to prochloraz correlates with proximity to Y123 and particularly, V381. These findings support a notable role for Y123 and V381 in maintaining sensitivity to prochloraz, and suggest that despite their proximity, K133 and V136 are not important in the interactions with triadimenol, tebuconazole and epoxiconazole.
Recent strains and compensatory sensitivity for tebuconazole and prochloraz
In the recent isolates, it is interesting to note that the I381V and A379G alterations increase resistance factors to tebuconazole but lower them for the imidazole prochloraz, whereas the V136A alteration increased resistance to prochloraz 
. It appears therefore that recent alterations serve to convey resistance to one of either tebuconazole or prochloraz, which is compensated by sensitivity to the other. This suggests that continued administration of combinations of azoles is feasible, providing that alterations are restricted to combinations of the currently altered residue set.
The final column of provides a summary of the main structural mechanisms underlying the resistance profiles observed for the main resistance variants. Three distinct mechanisms emerge, namely obstruction of binding or loss of interaction due to the replacement of a specific residue; constriction of the binding cavity by incorporation of the I381V alteration, resulting in the cavity being inaccessible to larger azoles such as tebuconazole; extensive increase in the volume of the binding cavity, leading to isolation of the azole from the key interacting residues.
Alterations at other sites
A consistent feature of all the models of azole docking in this study is the proximity of A311. It is the only residue predicted to be closer than 4.5 Å to epoxiconazole, triadimenol, tebuconazole or prochloraz. The equivalent residue has previously been identified in M. fijiensis
as being mutated (A313G) though a link to azole resistance could not be confirmed 
. An alignment of M. graminicola
and M. fijiensis
CYP51 sequences with Mycobacterium tuberculosis
reveals that the residue lies in the substrate signal recognition sequence 4 (SRS4) and is also located within the I helix 
. Substitutions at this residue (A311G) have also been identified in M. graminicola
, and indeed used to be more common (Fraaije et al., unpublished). A mutation affecting A311 is likely to further reduce DMI sensitivity. The modelling also implicates the neighbouring residue G312, though to a lesser extent (table S2
), and the G312A substitution has also been found previously. D134 is a residue in close proximity to bound triadimenol and epoxiconazole. Interestingly the D134G substitution used to be quite rare but now appears on the increase (Cools et al., unpublished). This illustrates the capacity of the modelling approach to not only rationalise the effects of alterations that have been characterised in the laboratory but also to identify other residues that may susceptible to alteration in the past or in the future and to test theoretical combinations in a predictive manner.
We have modelled the Mycosphaerella graminicola CYP51 enzyme along with 13 altered sequences containing one or a combination of amino acid alterations which have been associated with azole resistance. The multi-homologue models were generated in a wholly automated fashion, blind of prior functional inference and unbiased by homologue selection. The structural models support observations of resistance in field populations indicating that we have created a robust model which may be used for prediction of resistance to azoles and possible new mutations that may arise.
Our molecular modelling provides a structural rationale for the emergence of azole resistance in M. graminicola and suggests an element of compensation between atomic rearrangement and cavity extension in the later variants has occurred in the evolutionary response of MgCYP51 to different fungicide exposures. Extending the volume of the haem cavity while limiting structural rearrangement would appear to be a feasible evolutionary solution that allows the accommodation of the larger azole molecules without their engagement in interaction, while limiting rearrangement of specific side chains important for the general function of the enzyme that is essential for growth.
Our molecular modelling of M. graminicola CYP51 and its variants provides a robust structure-function rationale for the binding of azoles including the importance of Y137 in the binding interactions and valuable insights of mechanisms conferring resistance. Examination of the wild type binding pocket implicates several residues in close proximity to the positions of the bound azole molecules that are subject to alteration, presenting a picture of the various ligands being spatially constrained by specific regions bordering the pocket that become the focus of genetic mutation events and those following mutations result in subsequent resistance to that group of related compounds. The most striking of these is the previously structurally uncharacterised 459–461 region, which is the focus of several different single site substitutions and a 2 amino acid deletion, at 459–460 in various resistant strains. In structural terms, the deletion is strikingly effective, resulting in the removal of the entire beta turn section from the vicinity of the binding pocket. The models account for several observed functional effects of specific mutations, including the loss of triadimenol sensitivity in isolates possessing the Y137F variant, the lower sensitivity to tebuconazole of I381V variants, increased resistance to prochloraz of V136A variants together with sensitivity to triadimenol and tebuconazole; and finally, resistance to tebuconazole and decreased sensitivity to epoxiconazole of the A379G I381V ΔY459/G460 combination coupled to sensitivity to prochloraz.
We provide detailed descriptions for how different combinations of mutations affect binding of four different azole fungicides based on molecular models of CYP51. All these findings support the idea that the different binding specificities of different azoles to particular CYP51 variants has contributed to observed amount of polymorphism in the CYP51 gene, as a result of the mosaic of different fungicides that have been applied in wheat fields of Europe.
The structure – function relationships underpinning the observation that resistance of new combinations of mutations to new azole fungicides may increase sensitivity to older azole fungicides might inform effective rotation of fungicides or could be exploited to create new fungicide mixtures.
The measurement of the volume of the haem cavity was a useful metric for assessment of the scope for general resistance to this group of compounds by virtue of their accommodation without the potential for key binding interactions. It appears that this combination of methods provides the basis to a reliable in silico predictive approach for assessing the probability of particular variants exhibiting resistance to particular azole fungicides that would be of value in predicting and managing the use of fungicides in the control of this important wheat pathogen. This approach to understanding binding of azoles can be used as a screen to identify new fungicidal compounds and also to predict which azole classes will be affected by any particular combination of CYP51 mutations.
In terms of studying protein evolution more generally, our analyses have provided much insight into the role of specific amino acids in binding azoles in M. graminicola, and how conformational changes brought about by the alteration of particular residues, often proximal to the ligand, can result in modified binding capacity of other key residues. In doing so, we have established a promising in silico protocol for the assessment of alterations and their effects in proteins where evolutionary stress is introduced by the extensive application of a closely related class of agents targeting a specific protein binding domain.