The elastic network interpolation (ENI)1 is a computationally efficient and physically realistic method to generate conformational transition intermediates between two forms of a given protein. However it can be asked whether these calculated conformations provide good representatives for these intermediates. In this study, we use ENI to generate conformational transition intermediates between the open form and the closed forms of adenylate kinase (AK). Based on Cα-only intermediates, we construct atomic intermediates by grafting all the atoms of known AK structures onto the Cα atoms and then perform CHARMM energy minimization to remove steric conflicts and optimize these intermediate structures. We compare the energy profiles for all intermediates from both the CHARMM force-field and from knowledge-based energy functions. We find that the CHARMM energies can successfully capture the two energy minima representing the open AK and closed AK forms, while the energies computed from the knowledge-based energy functions can detect the local energy minimum representing the closed AK form and show some general features of the transition pathway with a somewhat similar energy profile as the CHARMM energies. The combinatorial extension (CE) structural alignment2 and the k-means clustering algorithm are then used to show that known PDB structures closely resemble computed intermediates along the transition pathway.
Available at: http://works.bepress.com/robert-jernigan/58/
This is the peer reviewed version of the following article: Feng, Yaping, Lei Yang, Andrzej Kloczkowski, and Robert L. Jernigan. "The energy profiles of atomic conformational transition intermediates of adenylate kinase." Proteins: Structure, Function, and Bioinformatics 77, no. 3 (2009): 551-558, which has been published in final form at DOI: 10.1002/prot.22467. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.