CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

A detailed computational investigation on the structural and spectroscopic properties of propolisbenzofuran B

This investigation deals with some structural and spectroscopic aspects of propolisbenzofuran B molecule as one of the most important bioactive molecules which exists in the bee propolis composition. FT-IR vibrational analysis carried-out at B3LYP/6-311++G(d,p) level of the theory. ¹H and ¹³C NMR chemical shift have been predicted with GIAO method. TD- DFT calculations have been established to predict the UV- Vis spectral analysis for propolisbenzofuran B molecule. The detailed structural analysis such as electronic characterization, HOMO and LUMO, DOS plot, Molecular Electronic Potential (MEP), Natural Bond Orbital (NBO) are performed and discussed for studied molecule.

Structural heterogeneity in a parent ground-state structure of AnPixJg2 revealed by theory and spectroscopy

We investigated the red absorbing, dark stable state (Pr state) of the second GAF domain of the cyanobacteriochrome AnPixJ (AnPixJg2) by a molecular dynamics simulation of 1 μs duration. Our results reveal two distinct conformational isoforms of the chromophore, from which only one was known from crystallographic experiments. The interconversion between both isoforms is accompanied by alterations in the hydrogen bond pattern between the chromophore and the protein and the solvation structure of the chromophore binding pocket. The existence of sub-states in the Pr form of AnPixJg2 is supported by the results from experimental ¹³C MAS NMR spectroscopy. Our finding is consistent with the observation of structural heterogeneity in other cyanobacteriochromes and phytochromes.

Ab-Initio Molecular Dynamics Simulations and Calculations of Spectroscopic Parameters in Hydrogen-Bonding Liquids in Confinement (Project 8)

We investigate the effect of several nanoscale confinements on structural and dynamical properties of liquid water and binary aqueous mixtures. By means of molecular dynamics simulations based on density functional theory and atomistic force fields. Our main focus is on the dependence on the structure and the hydrogen-bonding-network of the liquids near the confinement interface at atomistic resolution. As a complementary aspect, spatially resolved profiles of the proton NMR chemical shift values are used to quantify the local strength of the hydrogen-bond-network.

Improved Quantum Chemical NMR Chemical Shift Prediction of Metabolites in Aqueous Solution toward the Validation of Unknowns

A quantum-chemistry based protocol, termed MOSS-DFT, is presented for the prediction of ¹³C and ¹H NMR chemical shifts of a wide range of organic molecules in aqueous solution, including metabolites. Molecular motif-specific linear scaling parameters are reported for five different density functional theory (DFT) methods (B97-2/pcS-1, B97-2/pcS-2, B97-2/pcS-3, B3LYP/pcS-2, and BLYP/pcS-2), which were applied to a large set of 176 metabolite molecules. The chemical shift root-mean-square deviations (RMSD) for the best method, B97-2/pcS-3, are 1.93 and 0.154 ppm for ¹³C and ¹H chemical shifts, respectively. Excellent results have been obtained for chemical shifts of methyl and aromatic ¹³C and ¹H that are not directly bonded to a heteroatom (O, N, S, or P) with RMSD values of 1.15/0.079 and 1.31/0.118 ppm, respectively. This study not only demonstrates how NMR chemical shift in aqueous environment can be improved over the commonly used global linear scaling approach, but also allows for motif-specific error estimates, which are useful for an improved chemical shift-based verification of metabolite candidates of metabolomics samples containing unknown components.

Phycocyanobilin in solution – a solvent triggered molecular switch

The chromophore phycocyanobilin changes its spectroscopic behaviour upon solvent change. Our calculations trace this effect back to conformational switching, induced by changes in the hydrogen bonding pattern.

Conformational Sampling by Ab Initio Molecular Dynamics Simulations Improves NMR Chemical Shift Predictions

Car-Parrinello molecular dynamics simulations were performed for N-methyl acetamide as a small test system for amide groups in protein backbones, and NMR chemical shifts were calculated based on the generated ensemble. If conformational sampling and explicit solvent molecules are taken into account, excellent agreement between the calculated and experimental chemical shifts is obtained. These results represent a landmark improvement over calculations based on classical molecular dynamics (MD) simulations especially for amide protons, which are predicted too high-field shifted based on the latter ensembles. We were able to show that the better results are caused by the solute-solvents interactions forming shorter hydrogen bonds as well as by the internal degrees of freedom of the solute. Inspired by these results, we propose our approach as a new tool for the validation of force fields due to its power of identifying the structural reasons for discrepancies between the experimental and calculated data.