RESUMO
The increase in the resistance of mutant strains of Neisseria gonorrhoeae to the antibiotic ceftriaxone is pronounced in the decrease in the second-order acylation rate constant, k2/KS, by penicillin-binding protein 2 (PBP2). These changes can be caused by both the decrease in the acylation rate constant, k2, and the weakening of the binding affinity, i.e., an increase in the substrate constant, KS. A501X mutations in PBP2 affect second-order acylation rate constants. The PBP2A501V variant exhibits a higher k2/KS value, whereas for PBP2A501R and PBP2A501P variants, these values are lower. We performed molecular dynamic simulations with both classical and QM/MM potentials to model both acylation energy profiles and conformational dynamics of four PBP2 variants to explain the origin of k2/KS changes. The acylation reaction occurs in two elementary steps, specifically, a nucleophilic attack by the oxygen atom of the Ser310 residue and C-N bond cleavage in the ß-lactam ring accompanied by the elimination of the leaving group of ceftriaxone. The energy barrier of the first step increases for PBP2 variants with a decrease in the observed k2/KS value. Submicrosecond classic molecular dynamic trajectories with subsequent cluster analysis reveal that the conformation of the ß3-ß4 loop switches from open to closed and its flexibility decreases for PBP2 variants with a lower k2/KS value. Thus, the experimentally observed decrease in the k2/KS in A501X variants of PBP2 occurs due to both the decrease in the acylation rate constant, k2, and the increase in KS.