Dataset Information

Interactions between phosphate and water in solution: a natural bond orbital based analysis in a QM/MM framework.

ABSTRACT: The natural bond orbital (NBO) and natural energy decomposition analysis (NEDA) calculations are used to analyze the interaction between mono-methyl phosphate-ester (MMP) and its solvation environment in a combined quantum mechanical/molecular mechanical (QM/MM) framework. The solute-solvent configurations are generated using a specific parametrization of the self-consistent-charge density functional tight-binding (SCC-DFTB) model for the MMP and TIP3P for water. The NBO and NEDA calculations are done with several QM/MM partitioning schemes with HF/6-31+G** as the QM level. Regardless of the size of the QM region, a notable amount of charge transfer is observed between MMP and the neighboring water molecules and the charge-transfer interactions are, in the NEDA framework, as important as the electric (electrostatic and polarization) components. This work illustrates that NBO based analyses are effective tools for probing intermolecular interactions in condensed phase systems.

SUBMITTER: Yang Y

PROVIDER: S-EPMC2596991 | biostudies-literature | 2007 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Interactions between phosphate and water in solution: a natural bond orbital based analysis in a QM/MM framework.

Yang Yang Y Cui Qiang Q

The journal of physical chemistry. B 20070329 16

The natural bond orbital (NBO) and natural energy decomposition analysis (NEDA) calculations are used to analyze the interaction between mono-methyl phosphate-ester (MMP) and its solvation environment in a combined quantum mechanical/molecular mechanical (QM/MM) framework. The solute-solvent configurations are generated using a specific parametrization of the self-consistent-charge density functional tight-binding (SCC-DFTB) model for the MMP and TIP3P for water. The NBO and NEDA calculations ar ...[more]

PMID: 17391023

Similar Datasets

Project description:Understanding the nature of the free-energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here, the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free-energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free-energy surface for the hydrolysis of phosphate monoesters, paying special attention to the comparison between the one water (1W) and two water (2W) paths for the proton-transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by nonsystematic QM/MM energy minimization, as well as by nonsystematic free-energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces that are necessary for reaching concrete conclusions. Here we report a systematic evaluation of the 2D (3D) free-energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provide an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of the methyl diphosphate (MDP) trianion is 6-9 kcal/mol higher than that the 2W path. This difference becomes slightly larger in the presence of the Mg(2+) ion because this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives a much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth pointing out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic model are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic model, which was calibrated on the observed barrier in solution and in which the TS charge distribution was similar to the that of the plateau (as was done in all of our previous EVB studies).

Project description:Many gram-positive pathogens possess external pili or fimbriae with which they adhere to host cells during the infection process. Unusual dual intramolecular isopeptide bonds between Asn and Lys side chains within the N-terminal and C-terminal domains of the pilus subunits have been observed initially in the Streptococcus pyogenes pilin subunit Spy0128 and subsequently in GBS52 from Streptococcus agalactiae, in the BcpA major pilin of Bacillus cereus and in the RrgB pilin of Streptococcus pneumoniae, among others. Within each pilin subunit, intramolecular isopeptide bonds serve to stabilize the protein. These bonds provide a means to withstand large external mechanical forces, as well as possibly assisting in supporting a conformation favored for pilin subunit polymerization via sortase transpeptidases. Genome-wide analyses of pili-containing gram-positive bacteria are known or suspected to contain isopeptide bonds in pilin subunits. For the autocatalytic formation of isopeptide cross-links, a conservation of three amino acids including Asn, Lys, and a catalytically important acidic Glu (or Asp) residue are responsible. However, the chemical mechanism of how isopeptide bonds form within pilin remains poorly understood. Although it is possible that several mechanistic paths could lead to isopeptide bond formation in pili, the requirement of a conserved glutamate and highly organized positioning of residues within the hydrophobic environment of the active site were found in numerous pilin crystal structures such as Spy0128 and RrgB. This suggests a mechanism involving direct coupling of lysine side chain amine to the asparagine carboxamide mediated by critical acid/base or hydrogen bonding interactions with the catalytic glutamate residue. From this mechanistic perspective, we used the QM/MM minimum free-energy path method to examine the reaction details of forming the isopeptide bonds with Spy0128 as a model pilin, specifically focusing on the role of the glutamate in catalysis. It was determined that the reaction mechanism likely consists of two major steps: the nucleophilic attack on C? by nitrogen in the unprotonated Lys ?-amino group and, then two concerted proton transfers occur during the formation of the intramolecular isopeptide bond to subsequently release ammonia. More importantly, within the dual active sites of Spy0128, Glu(117), and Glu(258) residues function as crucial catalysts for each isopeptide bond formation, respectively, by relaying two proton transfers. This work also suggests that domain-domain interactions within Spy0128 may modulate the reactivity of residues within each active site. Our results may hopefully shed light on the molecular mechanisms of pilin biogenesis in gram-positive bacteria.

Dataset Information

Interactions between phosphate and water in solution: a natural bond orbital based analysis in a QM/MM framework.

Publications

Interactions between phosphate and water in solution: a natural bond orbital based analysis in a QM/MM framework.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets