Flexibility of aromatic residues in acetylcholinesterase

From Proteopedia

Flexibility of aromatic residues in the active-gorge of AChE

The high aromatic content of the deep and narrow active-site gorge of acetylcholinesterase (AChE) is a remarkable feature of this enzyme.There are lined along the gorge of Torpedo californica AChE (TcAChE), F120, F288, F290, F330, F331, W84, W233, W279, W432, Y70, Y121, Y130, Y334, and Y442. The side-chain conformational analyses based on the multuple available crystal structures and molecular dyanmics (MD) simulation trajectories show that the degree of flexibility of these 14 aromatic side chains is diverse. While those of are both very flexible, the side-chain conformations of appear to be fixed. Residues located on, or adjacent to the , viz. , display different flexibilities in the MD simulations and in the crystal structures.

Group I: residues with fixed side-chain conformation

Fig. 1. The χ₁/χ₂ plot of F120.

F120 as well as W233, W432, Y70, Y121, F331, F288, and F290 fall into group I, the category of residues with fixed side-chain conformations in both crystal structures and MD simulation trajectory. Fig. 1 shows that both the experimental and MD data are concentrated in a specific region. In this figure, the grey dots are derived from a 20-ns MD trajectory of native TcAChE and each one represents a pair of χ₁ and χ₂ angles calculated based on snapshot structures of TcAChE extracted from the MD trajectory at 1 ps intervals. The plot contains 20,000 such dots in total. The red pentacle is the crystal structure (pdb code 1ea5) used for the MD simulation. The dark triangles present pairs of c1 and c2 angles calculated based on 89 crystal structures of AChE deposited in the PDB.

Group II: residues with flexible side chains

Fig. 2. The χ₁/χ₂ plots of W279. Details in the plots are as in Fig. 1 except that the grey areas in the right plot are the favorable regions for the side-chain of Trp predicted by PROCHECK.

The side-chain flexibility of W279 is quite high in both data sets. The MD simulation of the native enzyme could produce most of side-chain conformations revealed by the crystal structures of complexes, suggesting that not the induced-fit but rather the preexisting equilibrium dynamics are involved in the conformational changes of W279 with regard to binding of these ligands. Moreover, it is striking that a standard χ₁/χ₂ rotomer plot of tryptophan would miss several of these regions both explored in MD simulation and crystal structures, such as group f (Figure 2).Rational drug design, based on structural information, might profit by taking into account conformational heterogeneity observed in MD simulations.

Fig. 3. a) The χ₁/χ₂ plot of F330. Details are as in Fig. 1. b)The χ₁/χ₂ of all Phe residues in TcAChE (pdb code 1ea5) overlayed on the favored regions for these residues derived by the program PROCHECK.

F330 is known for to have high side-chain flexibility as revealed by earlier experimental and simulation studies. The same conclusion was obtained in the current study. The only puzzle is that in some native structures F330 adopts an unfavored side-chain orientation according to PROCHECK. The re-examination of the electron density of these ‘strange’ native structures indicated that there may be a PEG bound in the active-site gorge of the structures in which the side-chain conformation of F330 is unfavorable.

Group III: Residues with different degree of side-chain flexibility in MD simulations and crytal structures

Fig. 4. The χ₁/χ₂ plots of W84. MD trajectory of the TcAChE monomer is shown without (a) or with (b) main-chain fixed during the simulation.

W84 as well as Y130, Y442, and Y334 behave differently in crystal structures and in MD trajectories. Comparison of the two data sets revealed that not only the side-chains of the four residues but also the main-chain of the omega-loop (C67-C94) on which the four residues are located or are adjacent to are different. In all the crystal structures, the conformation of the omega-loop is quite fixed because of the crystal packing. In contrast, this loop is very flexible in the MD simulation. A 2nd MD simulation, with main-chain fixed, which mimics the effect of the crystal packing, however, produced side-chain conformations similar to those in the crystal structures (Fig. 4b).

Conclusions

The side-chain flexibility analyses of the 14 aromatic residues based on the multiple crystal structures, together with the MD simulation trajectories, could benefit structure-based drug design for AChE, and understanding of the dynamic properties of the active-site gorge as well as ligand trafficking within it. The good agreement between conformational sampling of the crystal structures and the MD simulations suggests that the latter is a valid method to explore the conformational landscape of the residues as well as of the whole protein.

Additional Resources

For additional information, see: Alzheimer's Disease

References

• [1] Yechun Xu, Jacques Philippe Colletier, Hualiang Jiang, Israel Silman, Joel L. Sussman, Martin Weik. Induced-fit of preexisting equilibrium dynamics? Lessons from protein crystallography and MD simulations on acetylcholinesterase and implications for structure-based drug design. Prot. Sci. 2008, 17, 601-605.
• [2] Yechun Xu, Jacques Philippe Colletier, Martin Weik, Hualiang Jiang, John Moult, Israel Silman, Joel L. Sussman. Flexibility of aromatic residues in the active-site gorge of acetylcholinesterase: X-ray vs MD. Biophys. J. 2008, 95, 2500-2511.

Ache_apo_14aromaticresidues