Hybrid RHF/MP2 geometry optimizations with the Effective Fragment Molecular Orbital Method

The frozen domain effective fragment molecular orbital method is extended to allow for the treatment of a single fragment at the MP2 level of theory. The approach is applied to the conversion of chorismate to prephenate by chorismate mutase, where the substrate is treated at the MP2 level of theory while the rest of the system is treated at the RHF level. MP2 geometry optimization is found to lower the barrier by up to 3.5 kcal/mol compared to RHF optimzations and ONIOM energy refinement and leads to a smoother convergence with respect to the basis set for the reaction profile. For double zeta basis sets the increase in CPU time relative to RHF is roughly a factor of two.Comment: 11 pages, 3 figure

Interface of the polarizable continuum model of solvation with semi-empirical methods in the GAMESS program

An interface between semi-empirical methods and the polarized continuum model (PCM) of solvation successfully implemented into GAMESS following the approach by Chudinov et al (Chem. Phys. 1992, 160, 41). The interface includes energy gradients and is parallelized. For large molecules such as ubiquitin a reasonable speedup (up to a factor of six) is observed for up to 16 cores. The SCF convergence is greatly improved by PCM for proteins compared to the gas phase

Protein structure validation and refinement using amide proton chemical shifts derived from quantum mechanics

We present the ProCS method for the rapid and accurate prediction of protein backbone amide proton chemical shifts - sensitive probes of the geometry of key hydrogen bonds that determine protein structure. ProCS is parameterized against quantum mechanical (QM) calculations and reproduces high level QM results obtained for a small protein with an RMSD of 0.25 ppm (r = 0.94). ProCS is interfaced with the PHAISTOS protein simulation program and is used to infer statistical protein ensembles that reflect experimentally measured amide proton chemical shift values. Such chemical shift-based structural refinements, starting from high-resolution X-ray structures of Protein G, ubiquitin, and SMN Tudor Domain, result in average chemical shifts, hydrogen bond geometries, and trans-hydrogen bond (h3JNC') spin-spin coupling constants that are in excellent agreement with experiment. We show that the structural sensitivity of the QM-based amide proton chemical shift predictions is needed to refine protein structures to this agreement. The ProCS method thus offers a powerful new tool for refining the structures of hydrogen bonding networks to high accuracy with many potential applications such as protein flexibility in ligand binding.Comment: PLOS ONE accepted, Nov 201

Metabolic fingerprinting of <i>Lactobacillus paracasei</i>:the optimal quenching strategy

BACKGROUND: Quenching in cold buffered methanol at −40 °C has long been the preferred method for sub-second inactivation of cell metabolism during metabolic fingerprinting. However, methanol is known to cause intracellular metabolite leakage of microbial cells, making the distinction between intra- and extracellular metabolites in microbial systems challenging. In this paper, we tested three quenching protocols proposed for microbial cultures: fast filtration, cold buffered methanol and cold glycerol saline. RESULTS: Our results clearly showed that cold glycerol saline quenching resulted in the best recovery of intracellular metabolites in Lactobacillus paracasei subsp. paracasei (L. paracasei). Membrane integrity assayed by propidium iodide revealed that approximately 10 % of the L. paracasei cell membranes were damaged by contact with the cold buffered methanol solution, whilst cold glycerol saline quenching led to minimal cell damage. Due to the nature of the L. paracasei culture, fast filtration took several minutes, which is far from ideal for metabolites with high intracellular turnover rates. CONCLUSION: The implementation of a reliable, reproducible quenching method is essential within the metabolomics community. Cold glycerol saline prevented leakage of intracellular metabolites, and, thus, allowed more accurate determinations of intracellular metabolite levels