Two Sides of Quantum-Based Modeling of Enzyme-Catalyzed Reactions: Mechanistic and Electronic Structure Aspects of the Hydrolysis by Glutamate Carboxypeptidase

Evolution of Ceftriaxone Resistance of Penicillin-Binding Proteins 2 Revealed by Molecular Modeling

Penicillin-binding proteins 2 (PBP2) are critically important enzymes in the formation of the bacterial cell wall. Inhibition of PBP2 is utilized in the treatment of various diseases, including gonorrhea. Ceftriaxone is the only drug used to treat gonorrhea currently, and recent growth in PBP2 resistance to this antibiotic is a serious threat to human health. Our study reveals mechanistic aspects of the inhibition reaction of PBP2 from the wild-type FA19 strain and mutant 35/02 and H041 strains of Neisseria Gonorrhoeae by ceftriaxone. QM(PBE0-D3/6-31G**)/MM MD simulations show that the reaction mechanism for the wild-type PBP2 consists of three elementary steps including nucleophilic attack, C-N bond cleavage in the β-lactam ring and elimination of the leaving group in ceftriaxone. In PBP2 from the mutant strains, the second and third steps occur simultaneously. For all considered systems, the acylation rate is determined by the energy barrier of the first step that increases in the order of PBP2 from FA19, 35/02 and H041 strains. Dynamic behavior of ES complexes is analyzed using geometry and electron density features including Fukui electrophilicity index and Laplacian of electron density maps. It reveals that more efficient activation of the carbonyl group of the antibiotic leads to the lower energy barrier of nucleophilic attack and larger stabilization of the first reaction intermediate. Dynamical network analysis of MD trajectories explains the differences in ceftriaxone binding affinity: in PBP2 from the wild-type strain, the β₃-β₄ loop conformation facilitates substrate binding, whereas in PBP2 from the mutant strains, it exists in the conformation that is unfavorable for complex formation. Thus, we clarify that the experimentally observed decrease in the second-order rate constant of acylation (k₂/K_S) in PBP2 from the mutant strains is due to both a decrease in the acylation rate constant k₂ and an increase in the dissociation constant K_S.

Review on the QM/MM Methodologies and Their Application to Metalloproteins

The multiscaling quantum mechanics/molecular mechanics (QM/MM) approach was introduced in 1976, while the extensive acceptance of this methodology started in the 1990s. The combination of QM/MM approach with molecular dynamics (MD) simulation, otherwise known as the QM/MM/MD approach, is a powerful and promising tool for the investigation of chemical reactions’ mechanism of complex molecular systems, drug delivery, properties of molecular devices, organic electronics, etc. In the present review, the main methodologies in the multiscaling approaches, i.e., density functional theory (DFT), semiempirical methodologies (SE), MD simulations, MM, and their new advances are discussed in short. Then, a review on calculations and reactions on metalloproteins is presented, where particular attention is given to nitrogenase that catalyzes the conversion of atmospheric nitrogen molecules N₂ into NH₃ through the process known as nitrogen fixation and the FeMo-cofactor.

Role of water coordination at zinc binding site and its catalytic pathway of dizinc creatininase: insights from quantum cluster approach

Creatininase is a key enzyme of creatinine-metabolizing pathway in mammals, and has a great potential for diagnostic application. It catalyzes the reversible conversion of creatinine to creatine. Here, we investigated its reaction mechanism with density functional theory in conjunction with the quantum cluster approach. Three reaction pathways in which several possible proton transfers assisted by either His178 or a water ligand to Zn1 (Wat2) or both were considered. DFT calculations reveal, depending on Wat2 coordination mode at Zn1, two competitive ring-opening pathways where His178 playing a central role as a proton shuttle or both His178 and Wat2 serving as a dual catalytic role as a base and an acid, respectively. Three elementary steps were proposed for the reaction: the first involves nucleophilic attack by a bridging hydroxide to the substrate and forms a gem-diolate intermediate, followed by a proton transfer from the gem-diolate to His178 (His178 protonation is a required step for efficient proton transfers). Finally, the second proton transfer from the protonated His178 or Wat2 to the amide of substrate leads to the ring opening. The first proton transfer is the rate-limiting step of the whole reaction, in consistent with previous experimental and computational studies. A detailed understanding of the reaction mechanism of the creatininase enzyme family will also be helpful for developing a biosensor for kidney function.