Human butyrylcholinesterase polymorphism: Molecular modeling

Prolonged apnoea following injection of ester-containing myoralaxants was first described in 1953. Because a large part of administered succinylcholine is shortly hydrolyzed by plasma butyrylcholinesterase (BChE) under normal conditions, prolonged apnoea was attributed to deficiency in BChE. It was found that BChE deficiency was due to genetic variations. Human BChE gene shows a large polyallelism. About 75 natural mutations of the BCHE gene have been documented so far [1]. Most of them cause alteration in BChE activity through point mutation effect on catalytic activity. Frame shifts and stop codons may also affect expression, or cause truncations in the sequence.

Recently, two novel BChE "silent" variants, Val204Asp [2] and Ala34Val [3], causing prolonged neuromuscular block after administration of mivacurium, were discovered. Mutations were genetically and kinetically characterized. The aim of the current study was to understand how these mutations determine "silent" phenotype.

Molecular dynamics studies were carried out with NAMD 2.9 software at the Lomonosov supercomputer. Charmm 36 force field was used, periodical boundary conditions, 1 atm pressure, 298 K. 100 ns molecular dynamics runs were performed for the wild-type BChE and its mutants Val204Asp and Ala34Val.

Unlike wild-type BChE, which retained its operative catalytic triad through the whole MD simulation, the catalytic triad of mutants was disrupted, making chemical step impossible. Val204Asp mutation leads to reorganization of hydrogen bonding network around the catalytic triad, which in turn increases the distance between catalytic residue main chains. Mutation Ala34Val, located on the protein surface, leads to increased fluctuations in the Ω-loop and subsequent disruption of the gorge structure, including disruption of the catalytic triad and formation of new hydrogen bonds involving catalytic center residues.

Comparative study of the "silent" Ala328Asp mutant and the catalytically active mutant Ala328Cys shows that MD approach can discriminate between the differential effects of point mutations at a same position.