Unraveling low-barrier hydrogen bonds in complex systems with a simple quantum topological criterion

Direct evidence of a low barrier hydrogen bond in the catalytic triad of a Serine protease

Serine proteases are one of the largest groups of enzymes, found in both eukaryotes and prokaryotes, and are responsible for many different functions. The detailed information about the hydrogen-bonds in the catalytic triad (Asp…His…Ser) of these enzymes is of importance in order to fully understand the mechanism of action. The aspartate of the triad is hydrogen bonded to the histidine but the exact nature of this bond has been under discussion for some time. It is either a common short ionic hydrogen bond (SIHB) or a delocalized low barrier hydrogen bond (LBHB) were the hydrogen bond is shorter. So far, the evidence for LBHB in proteins have not been conclusive. Here we show clear NMR evidence that LBHB does exist in NS3, a serine protease from Dengue. The one bond coupling constant between the hydrogen and nitrogen was shown to be only 52 Hz instead of the usual 90 Hz. This together with a 1H chemical shift of 19.93 ppm is evidence that the hydrogen bond distance between His and Asp is shorter than for SIHB. Our result clearly shows the existence of LBHB and will help in understanding the mechanism of the catalytic triad in the important group of serine proteases.

Ligand-Induced Proton Transfer and Low-Barrier Hydrogen Bond Revealed by X-ray Crystallography

Ligand binding can change the pKa of protein residues and influence enzyme catalysis. Herein, we report three ultrahigh resolution X-ray crystal structures of CTX-M β-lactamase, directly visualizing protonation state changes along the enzymatic pathway: apo protein at 0.79 Å, precovalent complex with nonelectrophilic ligand at 0.89 Å, and acylation transition state (TS) analogue at 0.84 Å. Binding of the noncovalent ligand induces a proton transfer from the catalytic Ser70 to the negatively charged Glu166, and the formation of a low-barrier hydrogen bond (LBHB) between Ser70 and Lys73, with a length of 2.53 Å and the shared hydrogen equidistant from the heteroatoms. QM/MM reaction path calculations determined the proton transfer barrier to be 1.53 kcal/mol. The LBHB is absent in the other two structures although Glu166 remains neutral in the covalent complex. Our data represents the first X-ray crystallographic example of a hydrogen engaged in an enzymatic LBHB, and demonstrates that desolvation of the active site by ligand binding can provide a protein microenvironment conducive to LBHB formation. It also suggests that LBHBs may contribute to stabilization of the TS in general acid/base catalysis together with other preorganized features of enzyme active sites. These structures reconcile previous experimental results suggesting alternatively Glu166 or Lys73 as the general base for acylation, and underline the importance of considering residue protonation state change when modeling protein-ligand interactions. Additionally, the observation of another LBHB (2.47 Å) between two conserved residues, Asp233 and Asp246, suggests that LBHBs may potentially play a special structural role in proteins.

The pathway for serial proton supply to the active site of nitrogenase: enhanced density functional modeling of the Grotthuss mechanism

Proton translocation along a chain of eight waters to the active site of nitrogenase is described in detail, using density functional simulations with a 269 atom system that includes surrounding amino acids.

Stereoselectivity through a network of non-classical CH weak interactions: a prospective study of a bicyclic organocatalytic scaffold

The preferred stereoisomeric product of this catalytic Diels–Alder reaction is in part determined by noncovalent CH⋯π interactions.

Are ionic liquids pairwise in gas phase? A cluster approach and in situ IR study

In this work, we discussed the vaporization and gas species of ionic liquids (ILs) by a cluster approach of quantum statistical thermodynamics proposed by R. Luwig (Phys. Chem. Chem. Phys., 10, 4333), which is a controversial issue up to date. Based on the different sized clusters (2-12 ion-pairs) of the condensed phase, the molar enthalpies of vaporization (ΔvapH, 298.15 K, 1bar) of four representative ILs, 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([Emim][NTf2]) 1-ethyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide ([Emmim][NTf2]) 1-ethyl-3-methylimidazolium chloride ([Emim]Cl) and ethylammonium nitrate ([EtAm][NO3]), were calculated. The predicted ΔvapH were increased remarkably; even the values of [EtAm][NO3] were larger than 700 kJ mol(-1) when the charged isolated ions were assumed to be gas species. However, the ΔvapH were close to experimental measurements when the gas species assumed to be anion-cation pairwise, indicating that the different conformational ion-pairs can coexist in the gas phase when the IL is evaporated. Particularly for the protic IL, [EtAm][NO3], even the neutral precursor molecules by proton transfer can occur in gas phase. In addition, it's found that the effect of hydrogen bonds on the vaporization cannot be negligible by comparing the ΔvapH of [Emim][NTf2] with [Emmim][NTf2]. The in situ and calculated IR spectra provided the further proof that the ions are pairwise in gas phase.

Correlation between electron localization and metal ion mutagenicity in DNA synthesis from QM/MM calculations

DNA polymerases require two divalent metal ions in the active site for catalysis. Mg(2+) has been confirmed to be the most probable cation utilized by most polymerases in vivo. Other metal ions are either potent mutagens or inhibitors. We used structural and topological analyses based on ab initio QM/MM calculations to study human DNA polymerase λ (Polλ) with different metals in the active site. Our results indicate a slightly longer O3'-Pα distance (∼3.6 Å) for most inhibitor cations compared to the natural and mutagenic metals (∼3.3-3.4 Å). Optimization with a larger basis set for the previously reported transition state (TS) structures (Cisneros et al., DNA Repair, 2008, 7, 1824.) gives barriers of 17.4 kcal mol(-1) and 15.1 kcal mol(-1) for the Mg(2+) and Mn(2+) catalyzed reactions respectively. Relying on the key relation between the topological signature of a metal cation and its selectivity within biological systems (de Courcy et al., J. Chem. Theor. Comput., 2010, 6, 1048.) we have performed electron localization function (ELF) topological analyses. These analyses show that all inhibitor and mutagenic metals considered, except Na(+), present a "split" of the outer-shell density of the metal. This "splitting" is not observed for the non-mutagenic Mg(2+) metal. Population and multipole analyses on the ELF basins reveal that the electronic dipolar and quadrupolar polarization is significantly different with Mg(2+) compared to all other cations. Our results shed light at the atomic level on the subtle differences between Mg(2+), mutagenic, and inhibitor metals in DNA polymerases. These results provide a correlation between the electronic distribution of the cations in the active site and the possible consequences on DNA synthesis.