Protocols for metallo- and serine-β-lactamase free energy predictions: insights from cross-class inhibitors

Jonna Jasmin Güven,Papu Kalita,Marko Hanževački,Adrian Mulholland,Antonia Mey

DOI: https://doi.org/10.26434/chemrxiv-2024-9ksxr

2024-09-23

Abstract:While relative binding free energy (RBFE) calculations using alchemical methods are routinely carried out for many pharmaceutically relevant protein targets, challenges remain. For example, open-source tools do not support the easy setup and simulation of metalloproteins, particularly when ligands directly coordinate to the metal site. Here, we evaluate the performance of RBFE methods for KPC-2, a serine-β-lactamase (SBL), and two non-bonded metal parameter setups for VIM-2, a metallo-β-lactamase (MBL) with two active site zinc ions. We tested two different ways of modeling the ligand-zinc interactions. First, a restraint-based approach, in which FF14SB zinc parameters are combined with harmonic restraints between the zincs and their coordinating residues. The second approach uses an upgraded Amber force field (UAFF) for zinc-metalloproteins with adjusted partial charges and non-bonded terms of zinc-coordinating residues. Molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) simulations show that the crystallographically observed zinc coordination is not retained in MM simulations with either zinc parameter set for a series of known phosphonic acid-based inhibitors bound to VIM-2. These phosphonic acid-based inhibitors exhibit known cross-class affinity for SBLs and MBLs and serve as a benchmark for RBFE calculations for VIM-2, after validation with KPC-2. The KPC-2 free energy of binding estimates are within expected literature accuracies for the lig- and series with a mean absolute error of 0.43 (0.25, 0.66) kcal/mol and a Pearson’s correlation coefficient of 0.93 (0.86, 0.98). For VIM-2, the UAFF approach has improved correlation from 0.48 (−0.02, 0.86) to 0.66 (−0.09, 0.94), compared to the restraint approach. The presented strategies for handling ligands coordinating to metal sites highlight that simple metal parameter models can provide some predictive free energy estimate for metalloprotein–ligand systems, but leave room for improvement in their ease of use, modeling of coordination sites, and as a result, their accuracy.

Chemistry

What problem does this paper attempt to address?

The paper attempts to address the issue of how to accurately predict the binding affinity between ligands and enzymes using the Relative Binding Free Energy (RBFE) method in metallo-β-lactamase and serine-β-lactamase. Specifically: 1. **Technical Challenges**: Current open-source tools face difficulties when dealing with metalloproteins, especially in cases where the ligand directly interacts with the metal site. This is particularly challenging for the parameterization of zinc ions in metallo-β-lactamase. 2. **Research Objectives**: To evaluate two different modeling approaches: one based on a constraint method and the other using an improved Amber force field (UAFF), to better simulate the interaction between ligands and metal sites, and to validate the accuracy of these methods in predicting binding free energies. 3. **Practical Applications**: By developing effective computational methods to aid in the design of new inhibitors, particularly in combating antibiotic-resistant bacteria, the study aims to improve the efficiency and accuracy of the drug discovery process. In summary, the paper aims to overcome the technical challenges of current computational methods in handling metalloproteins, thereby providing a reliable theoretical foundation for designing broad-spectrum β-lactamase inhibitors.

Protocols for metallo- and serine-β-lactamase free energy predictions: insights from cross-class inhibitors

A multiscale approach to predict the binding mode of metallo beta‐lactamase inhibitors

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations

Non-bonded force field parameters from MBIS partitioning of the molecular electron density improve binding affinity predictions of the T4-lysozyme double mutant

Meta-Analysis Reveals That Absolute Binding Free-Energy Calculations Approach Chemical Accuracy

Predictions Of Binding Of A Diverse Set Of Ligands To Gelatinase-A By A Combination Of Molecular Dynamics And Continuum Solvent Models

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Predicting antibiotic resistance in complex protein targets using alchemical free energy methods

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Comparison of Methodologies for Absolute Binding Free Energy Calculations of Ligands to Intrinsically Disordered Proteins

Fast and Effective Prediction of the Absolute Binding Free Energies of Covalent Inhibitors of SARS-CoV-2 Main Protease and 20S Proteasome

Towards chemical accuracy for alchemical free energy calculations with hybrid physics-based machine learning / molecular mechanics potentials

Calculation of Absolute Protein-Ligand Binding Affinity Using Path and Endpoint Approaches

Equilibrium and Non-equilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors

An evaluation of machine learning/molecular mechanics end-state corrections with mechanical embedding to calculate relative protein-ligand binding free energies

Nonbonded Force Field Parameters from MBIS Partitioning of the Molecular Electron Density Improve Binding Affinity Predictions of the T4-Lysozyme Double Mutant

Performance and Analysis of the Alchemical Transfer Method for Binding Free Energy Predictions of Diverse Ligands

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery