Dataset Information

Discerning the catalytic mechanism of Staphylococcus aureus sortase A with QM/MM free energy calculations.

ABSTRACT: Sortases are key virulence factors in Gram-positive bacteria. These enzymes embed surface proteins in the cell wall through a transpeptidation reaction that involves recognizing a penta-peptide "sorting signal" in a target protein, cleaving it, and covalently attaching it to a second substrate that is later inserted into the cell wall. Although well studied, several aspects of the mechanism by which sortases perform these functions remains unclear. In particular, experiments have revealed two potential sorting signal binding motifs: a "Threonine-Out" (Thr-Out) structure in which the catalytically critical threonine residues protrudes into solution, and a "Threonine-In" (Thr-In) configuration in which this residue inserts into the binding site. To determine which of these is the biologically relevant state, we have performed a series of conventional and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations of the Staphylococcus aureus sortase A (SrtA) enzyme bound to a sorting signal substrate. Through the use of multi-dimensional metadynamics, our simulations were able to both map the acylation mechanism of SrtA in the Thr-In and Thr-Out states, as well as determine the free energy minima and barriers along these reactions. Results indicate that in both states the catalytic mechanisms are similar, however the free energy barriers are lower in the Thr-In configuration, suggesting that Thr-In is the catalytically relevant state. This has important implications for advancing our understanding of the mechanisms of sortase enzymes, as well we for future structure based drug design efforts aimed at inhibiting sortase function in vivo.

SUBMITTER: Shrestha P

PROVIDER: S-EPMC4917473 | biostudies-literature | 2016 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Discerning the catalytic mechanism of Staphylococcus aureus sortase A with QM/MM free energy calculations.

Shrestha Pooja P Wereszczynski Jeff J

Journal of molecular graphics & modelling 20160427

Sortases are key virulence factors in Gram-positive bacteria. These enzymes embed surface proteins in the cell wall through a transpeptidation reaction that involves recognizing a penta-peptide "sorting signal" in a target protein, cleaving it, and covalently attaching it to a second substrate that is later inserted into the cell wall. Although well studied, several aspects of the mechanism by which sortases perform these functions remains unclear. In particular, experiments have revealed two po ...[more]

PMID: 27172839

Similar Datasets

Project description:Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly, the high computational cost involved in performing proper configurational sampling and, in particular, in obtaining ab initio QM/MM (QM(ai)/MM) free-energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work, we point out the subtle efficiency problems associated with this approach and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. 1995, 99, 17516)], which we refer to as paradynamics (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free-energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free-energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely, the S(N)2 reaction between a chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost-effective than the metadynamics approach (at least in its current implementation). We also discuss the general features of the PD approach in terms of its ability to explore complex systems and clarify that it is not a specialized approach limited to only accelerating QM(ai)/MM calculations with proper sampling, but rather can be used in a wide variety of applications. In fact, we point out that the use of a reference (CG) potential coupled with its PD refinement, as well as our renormalization approach, provides very general and powerful strategies that can be used very effectively to explore any property that has been studied by the MTD approach.

Project description:We use the PBE0/6-31G* density functional method to perform ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations under periodic boundary conditions with rigorous electrostatics using the ambient potential composite Ewald method in order to test the convergence of MM ? QM/MM free energy corrections for the prediction of 17 small-molecule solvation free energies and eight ligand binding free energies to T4 lysozyme. The "indirect" thermodynamic cycle for calculating free energies is used to explore whether a series of reference potentials improve the statistical quality of the predictions. Specifically, we construct a series of reference potentials that optimize a molecular mechanical (MM) force field's parameters to reproduce the ab initio QM/MM forces from a QM/MM simulation. The optimizations form a systematic progression of successively expanded parameters that include bond, angle, dihedral, and charge parameters. For each reference potential, we calculate benchmark quality reference values for the MM ? QM/MM correction by performing the mixed MM and QM/MM Hamiltonians at 11 intermediate states, each for 200 ps. We then compare forward and reverse application of Zwanzig's relation, thermodynamic integration (TI), and Bennett's acceptance ratio (BAR) methods as a function of reference potential, simulation time, and the number of simulated intermediate states. We find that Zwanzig's equation is inadequate unless a large number of intermediate states are explicitly simulated. The TI and BAR mean signed errors are very small even when only the end-state simulations are considered, and the standard deviations of the TI and BAR errors are decreased by choosing a reference potential that optimizes the bond and angle parameters. We find a robust approach for the data sets of fairly rigid molecules considered here is to use bond + angle reference potential together with the end-state-only BAR analysis. This requires QM/MM simulations to be performed in order to generate reference data to parametrize the bond + angle reference potential, and then this same simulation serves a dual purpose as the full QM/MM end state. The convergence of the results with respect to time suggests that computational resources may be used more efficiently by running multiple simulations for no more than 50 ps, rather than running one long simulation.

Dataset Information

Discerning the catalytic mechanism of Staphylococcus aureus sortase A with QM/MM free energy calculations.

Publications

Discerning the catalytic mechanism of Staphylococcus aureus sortase A with QM/MM free energy calculations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets