Novichkova_2019_Chem.Biol.Interact_13ChEPon_309_108699

The crystal structures of truncated forms of cholinesterases provide good models for assessing the role of non-covalent interactions in dimer assembly in the absence of cross-linking disulfide bonds. These structures identify the four-helix bundle that serves as the interface for formation of acetylcholinesterase and butyrylcholinesterase dimers. Here we performed a theoretical comparison of the structural and energetic factors governing dimerization. This included identification of inter-subunit and intra-subunit hydrogen bonds and hydrophobic interactions, evaluation of solvent-accessible surfaces, and estimation of electrostatic contributions to dimerization. To reveal the contribution to dimerization of individual amino acids within the contact area, free energy perturbation alanine screening was performed. Markov state modelling shows that the loop between the alpha13 and alpha14 helices in BChE is unstable, and occupies 4 macro-states. The order of magnitude of mean first passage times between these macrostates is ~10(-8)s. Replica exchange molecular dynamics umbrella sampling calculations revealed that the free energy of human BChE dimerization is -15.5kcal/mol, while that for human AChE is -26.4kcal/mol. Thus, the C-terminally truncated human butyrylcholinesterase dimer is substantially less stable than that of human acetylcholinesterase. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:CHEMBIOINT:1.