CSGT-DFT calculation of 13C and 15N NMR shielding of the backbone amide group, 13Cα, and 13Cβ in ω-Conotoxin GVIA

Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes

NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients (r²) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.

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.

Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 3. Conformational Sampling and Explicit Solvents Model

Fragment-based quantum chemical calculations are able to accurately calculate NMR chemical shifts even for very large molecules like proteins. But even with systematic optimization of the level of theory and basis sets as well as the use of implicit solvents models, some nuclei like polar protons and nitrogens suffer from poor predictions. Two properties of the real system, strongly influencing the experimental chemical shifts but almost always neglected in the calculations, will be discussed here in great detail: (1) conformational averaging and (2) interactions with first-shell solvent molecules. Classical molecular dynamics simulations in explicit water were carried out for obtaining a representative ensemble including the arrangement of neighboring solvent molecules, which was then subjected to quantum chemical calculations. We could demonstrate with the small test system N-methyl acetamide (NMA) that the calculated chemical shifts show immense variations of up to 6 ppm and 50 ppm for protons and nitrogens, respectively, depending on the snapshot taken from a classical molecular dynamics simulation. Applying the same approach to the HA2 domain of the influenza virus glycoprotein hemagglutinin, a 32-amino-acid-long polypeptide, and comparing averaged values to the experiment, chemical shifts of nonpolar protons and carbon atoms in proteins were calculated with unprecedented accuracy. Additionally, the mean absolute error could be reduced by a factor of 2.43 for polar protons, and reasonable correlations were obtained for nitrogen and carbonyl carbon in contrast to all other studies published so far.

Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 2. Level of Theory, Basis Set, and Solvents Model Dependence

It has been demonstrated that the fragmentation scheme of our adjustable density matrix assembler (ADMA) approach for the quantum chemical calculations of very large systems is well-suited to calculate NMR chemical shifts of proteins [ Frank et al. Proteins2011, 79, 2189-2202 ]. The systematic investigation performed here on the influences of the level of theory, basis set size, inclusion or exclusion of an implicit solvent model, and the use of partial charges to describe additional parts of the macromolecule on the accuracy of NMR chemical shifts demonstrates that using a valence triple-ζ basis set leads to large improvement compared to the results given in the previous publication. Additionally, moving from the B3LYP to the mPW1PW91 density functional and including partial charges and implicit solvents gave the best results with mean absolute errors of 0.44 ppm for hydrogen atoms excluding H(N) atoms and between 1.53 and 3.44 ppm for carbon atoms depending on the size and also on the accuracy of the protein structure. Polar hydrogen and nitrogen atoms are more difficult to predict. For the first, explicit hydrogen bonds to the solvents need to be included and, for the latter, going beyond DFT to post-Hartree-Fock methods like MP2 is probably required. Even if empirical methods like SHIFTX+ show similar performance, our calculations give for the first time very reliable chemical shifts that can also be used for complexes of proteins with small-molecule ligands or DNA/RNA. Therefore, taking advantage of its ab initio nature, our approach opens new fields of application that would otherwise be largely inaccessible due to insufficient availability of data for empirical parametrization.

Backbone conformational preferences and pseudorotational ring puckering of 1-aminocyclopentane-1-carboxylic acid

We have used quantum mechanical calculations at the B3LYP/6-311G(d,p) level to determine the conformational preferences of the N-acetyl-N'-methylamide derivative of 1-aminocyclopentane-1-carboxylic acid in the gas phase, chloroform solution, and water solution. The backbone conformation of this dipeptide has been described through the dihedral angles varphi and psi, while the pseudorotational phase angle was used to define the conformation of the cyclopentane ring. Results indicate that the backbone flexibility of this amino acid is restricted by the cyclic nature of the side chain, the relative stability of the different conformations depending on the polarity of the environment. The potential energy of the pseudorotation was also studied as a function of the backbone conformation. Interestingly, the conformation of the cyclic side chain depends on the backbone arrangement. Furthermore, the number of pseudorotational states accessible at room temperature is high in all the investigated environments, especially in aqueous solution. Finally, a set of force-field parameters for classical molecular mechanics calculations was developed for the investigated amino acid. Molecular dynamics simulations in both chloroform and aqueous solutions were performed to demonstrate the reliability of such parameters.