Determinants of Regioselectivity and Chemoselectivity in Fosfomycin Resistance Protein FosA from QM/MM Calculations

A molecular model to study FosA enzyme inhibition

Fosfomycin resistance protein (FosA) is a metalloenzyme known for catalyzing a nucleophilic addition reaction of glutathione to the epoxide ring of Fosfomycin, a broad-spectrum antibiotic used to combat Gram-positive pathogens. The reaction leads fosfomycin to lose its pharmacological effect, thus promotes antibiotic resistance. A small-molecule FosA inhibitor has been discovered. ANY1 (3-bromo-6-[3-(3-bromo-2-oxo-1H-pyrazolo[1,5-a]pyrimidin-6-yl)-4-nitro-1H-pyrazol-5-yl]-1H-pyrazolo[1,5-a]pyrimidin-2-one) is competitive with the antibiotic for binding the active site of the enzyme. Through Molecular Mechanics methods, using the AMBER force field, we carry out molecular dynamics simulations and binding free energy calculations to investigate the most important interactions between the enzyme and inhibitor. Our results were able to reproduce the trend of experimental data with R² of 77.51%. Furthermore, we have shown that electrostatic and van der Waals interactions, as well as cavitation energies, are favorable for maintaining the enzyme-inhibitor complex, while reactive field energies and non-polar interactions act in an unfavorable way for interactions between FosA and ANY1.

A QM/MM Evaluation of the Missing Step in the Reduction Mechanism of HMG-CoA by Human HMG-CoA Reductase

Statins are important drugs in the regulation of cholesterol levels in the human body that have as a primary target the enzyme β-hydroxy-β-methylglutaryl-CoA reductase (HMGR). This enzyme plays a crucial role in the mevalonate pathway, catalyzing the four-electron reduction of HMG-CoA to mevalonate. A second reduction step of this reaction mechanism has been the subject of much speculation in the literature, with different conflicting theories persisting to the present day. In this study, the different mechanistic hypotheses were evaluated with atomic-level detail through a combination of molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The obtained Gibbs free activation and Gibbs free reaction energy (15 kcal mol−1 and −40 kcal mol−1) show that this hydride step takes place with the involvement of a cationic His405 and Lys639, and a neutral Glu98, while Asp715 remains in an anionic state. The results provide an atomic-level portrait of this step, clearly demonstrating the nature and protonation state of the amino acid residues involved, the energetics associated, and the structure and charge of the key participating atoms in the several intermediate states, finally elucidating this missing step.

QM/MM Study of the Fosfomycin Resistance Mechanism Involving FosB Enzyme

Multidrug-resistant organisms contain antibiotic-modifying enzymes that facilitate resistance to a variety of antimicrobial compounds. Particularly, the fosfomycin (FOF) drug can be structurally modified by several FOF-modifying enzymes before it reaches the biological target. Among them, FosB is an enzyme that utilizes l-cysteine or bacillithiol in the presence of a divalent metal to open the epoxide ring of FOF and, consequently, inactivate the drug. Here, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations to explore the mechanism of the reaction involving FosB and FOF. The calculated free-energy profiles show that the cost to open the epoxide ring of FOF at the C2 atom is ∼3.0 kcal/mol higher than that at the C1 atom. Besides, our QM/MM MD results revealed the critical role of conformation change of Cys9 and Asn50 to release the drug from the active site. Overall, the present study provides insights into the mechanism of FOF-resistant proteins.

Computational Understanding of the Selectivities in Metalloenzymes

Metalloenzymes catalyze many different types of biological reactions with high efficiency and remarkable selectivity. The quantum chemical cluster approach and the combined quantum mechanics/molecular mechanics methods have proven very successful in the elucidation of the reaction mechanism and rationalization of selectivities in enzymes. In this review, recent progress in the computational understanding of various selectivities including chemoselectivity, regioselectivity, and stereoselectivity, in metalloenzymes, is discussed.

Pharmacological properties of oral antibiotics for the treatment of uncomplicated urinary tract infections

The therapeutic management of uncomplicated bacterial urinary tract infections (UTIs) is based on short-term courses of oral antibiotics. The preferred drugs are nitrofurantoin trimethoprim-sulfamethoxazole, fosfomycin trometamol, fluoroquinolones and β-lactam agents. The choice of agent for treating uncomplicated UTIs should be based on the pharmacokinetic characteristics of the molecule so that clinical benefit is optimized and the risk of antibacterial resistance is minimized. This article discusses the general pharmacokinetic-pharmacodynamic (PK/PD) aspects of antimicrobial chemotherapy, the PK/PD characteristics of oral antimicrobial agents for the treatment of uncomplicated UTIs and the pharmacological and therapeutic strategies for limiting or preventing bacterial resistance.

Structure and Dynamics of FosA-Mediated Fosfomycin Resistance in Klebsiella pneumoniae and Escherichia coli

Fosfomycin exhibits broad-spectrum antibacterial activity and is being reevaluated for the treatment of extensively drug-resistant pathogens. Its activity in Gram-negative organisms, however, can be compromised by expression of FosA, a metal-dependent transferase that catalyzes the conjugation of glutathione to fosfomycin, rendering the antibiotic inactive. In this study, we solved the crystal structures of two of the most clinically relevant FosA enzymes: plasmid-encoded FosA3 from Escherichia coli and chromosomally encoded FosA from Klebsiella pneumoniae (FosA^KP). The structure, molecular dynamics, catalytic activity, and fosfomycin resistance of FosA3 and FosA^KP were also compared to those of FosA from Pseudomonas aeruginosa (FosA^PA), for which prior crystal structures exist. E. coli TOP10 transformants expressing FosA3 and FosA^KP conferred significantly greater fosfomycin resistance (MIC, >1,024 μg/ml) than those expressing FosA^PA (MIC, 16 μg/ml), which could be explained in part by the higher catalytic efficiencies of the FosA3 and FosA^KP enzymes. Interestingly, these differences in enzyme activity could not be attributed to structural differences at their active sites. Instead, molecular dynamics simulations and hydrogen-deuterium exchange experiments with FosA^KP revealed dynamic interconnectivity between its active sites and a loop structure that extends from the active site of each monomer and traverses the dimer interface. This dimer interface loop is longer and more extended in FosA^KP and FosA3 than in FosA^PA, and kinetic analyses of FosA^KP and FosA^PA loop-swapped chimeric enzymes highlighted its importance in FosA activity. Collectively, these data yield novel insights into fosfomycin resistance that could be leveraged to develop new strategies to inhibit FosA and potentiate fosfomycin activity.

QM/MM study of the mechanism of reduction of 3-hydroxy-3-methylglutaryl coenzyme A catalyzed by human HMG-CoA reductase

Detailing with atomistic resolution the reaction mechanism of human HMG-CoA reductase (HMG-CoA-R) might provide valuable insights for the development of new cholesterol-lowering drugs.

Mechanism and selectivity of the dinuclear iron benzoyl-coenzyme A epoxidase BoxB

DFT calculations are used to elucidate the reaction mechanism and selectivity of BoxB catalyzed benzoyl-CoA epoxidation.

Benzoyl-CoA epoxidase is a dinuclear iron enzyme that catalyzes the epoxidation reaction of the aromatic ring of benzoyl-CoA with chemo-, regio- and stereo-selectivity. It has been suggested that this enzyme may also catalyze the deoxygenation reaction of epoxide, suggesting a unique bifunctionality among the diiron enzymes. We report a density functional theory study of this enzyme aimed at elucidating its mechanism and the various selectivities. The epoxidation is suggested to start with the binding of the O₂ molecule to the diferrous center to generate a diferric peroxide complex, followed by concerted O–O bond cleavage and epoxide formation. Two different pathways have been located, leading to (2S,3R)-epoxy and (2R,3S)-epoxy products, with barriers of 17.6 and 20.4 kcal mol^–1, respectively. The barrier difference is 2.8 kcal mol^–1, corresponding to a diastereomeric excess of about 99 : 1. Further isomerization from epoxide to phenol is found to have quite a high barrier, which cannot compete with the product release step. After product release into solution, fast epoxide–oxepin isomerization and racemization can take place easily, leading to a racemic mixture of (2S,3R) and (2R,3S) products. The deoxygenation of epoxide to regenerate benzoyl-CoA by a diferrous form of the enzyme proceeds via a stepwise mechanism. The C2–O bond cleavage happens first, coupled with one electron transfer from one iron center to the substrate, to form a radical intermediate, which is followed by the second C3–O bond cleavage. The first step is rate-limiting with a barrier of only 10.8 kcal mol^–1. Further experimental studies are encouraged to verify our results.