Lushchekina_2016_Chem.Biol.Interact_259_223

Conformational dynamics of wild-type human butyrylcholinesterase (BChE), two mutants of residue Ala328, the catalytically active Ala328Cys, and the catalytically inactive (silent) Ala328Asp, and their interactions with butyrylcholine were studied. The aim was to understand the molecular mechanisms by which point mutations may lead to silent BChE variant or alter catalytic activity. Importance of BChE natural variants is due to medical consequences, i.e. prolonged apnea, following administration of the myorelaxant esters, succinylcholine and mivacurium. Comparison of molecular dynamics (MD) simulations for the three model systems showed that: 1) the active mutant Ala328Cys mutant has some changes in configuration of catalytic residues, which do not prevent binding of butyrylcholine to the active site; 2) in the naturally-occurring silent variant Ala328Asp, the Asp328 carboxylate may either form a salt bridge with Lys339 or a H-bond with His438. In the first case, the Omega-loop swings off the gorge, disrupting the pi-cation binding site and the catalytic triad. In the second case, binding of cationic substrates in the catalytic center is also impaired. MD simulations carried out in 0.15 M NaCl, close to physiological ionic strength conditions, favored the second situation. It was seen that Asp328 forms a H-bond with the catalytic triad His438, which in turn disrupts the catalytic machinery. Therefore, we concluded that the Ala328Asp variant is not catalytically active because of that dramatic event. Computational results, consistent with in vitro biochemical data and clinical observations, validate our MD approach.