Dataset Information

Theoretical Infrared Spectra: Quantitative Similarity Measures and Force Fields.

ABSTRACT: Infrared spectroscopy can provide significant insight into the structures and dynamics of molecules of all sizes. The information that is contained in the spectrum is, however, often not easily extracted without the aid of theoretical calculations or simulations. We present here the calculation of the infrared spectra of a database of 703 gas phase compounds with four different force fields (CGenFF, GAFF-BCC, GAFF-ESP, and OPLS) using normal-mode analysis. Modern force fields increasingly use virtual sites to describe, e.g., lone-pair electrons or the ?-holes on halogen atoms. This requires some adaptation of code to perform normal-mode analysis of such compounds, the implementation of which into the GROMACS software is briefly described as well. For the quantitative comparison of the obtained spectra with experimental reference data, we discuss the application of two different statistical correlation coefficients, Pearson and Spearman. The advantages and drawbacks of the different methods of comparison are discussed, and we find that both methods of comparison give the same overall picture, showing that present force field methods cannot match the performance of quantum chemical methods for the calculation of infrared spectra.

SUBMITTER: Henschel H

PROVIDER: S-EPMC7304875 | biostudies-literature | 2020 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Theoretical Infrared Spectra: Quantitative Similarity Measures and Force Fields.

Henschel Henning H Andersson Alfred T AT Jespers Willem W Mehdi Ghahremanpour Mohammad M van der Spoel David D

Journal of chemical theory and computation 20200423 5

Infrared spectroscopy can provide significant insight into the structures and dynamics of molecules of all sizes. The information that is contained in the spectrum is, however, often not easily extracted without the aid of theoretical calculations or simulations. We present here the calculation of the infrared spectra of a database of 703 gas phase compounds with four different force fields (CGenFF, GAFF-BCC, GAFF-ESP, and OPLS) using normal-mode analysis. Modern force fields increasingly use vi ...[more]

PMID: 32271575

Similar Datasets

Project description:Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3-7 cm(-1) (compared to the experimental value of 10 cm(-1)) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe...CO and Fe...OC) is 0.3 kcal mol(-1) and agrees well with the experimentally estimated value of 0.29 kcal mol(-1). A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

Project description:Infrared photodissociation analyses are supported by theoretical calculations that allow a trustworthy interpretation of experimental spectra of gaseous ions. B3LYP calculations are the most prominent method used to model IR spectra, as detailed in our bibliographic survey. However, this and other commonly used methods are known to provide inaccurate energy values and geometries, especially when it comes to long-range interactions, such as intramolecular H-bonds, which show increased anharmonicity. Therefore, we evaluated some of the most commonly used density functional theory methods (B3LYP, CAM-B3LYP, and M06-2X) and basis sets (6-31+G(d,p), 6-311++G(d,p), 6-311++G(3df,2pd), aug-cc-pVDZ, and aug-cc-pVTZ), including anharmonicity and dispersion corrections. The results were compared to MP2 calculations and to experimental high-frequency (2000-4000 cm-1) IR multiphotonic dissociation (IRMPD) spectra of two protonated model molecules containing intramolecular hydrogen bonds: biotin and tryptophan. M06-2X/6-31+G(d,p) was shown to be the most cost-effective level of theory, whereas CAM-B3LYP was the most efficient method to describe the van der Waals interactions. The use of the dispersion correction D3, proposed by Grimme, improved the description of O-H vibrations involved in H-bonding but worsened the description of N-H stretches. Anharmonic calculations were shown to be extremely expensive when compared to other approaches. The efficiencies of well-established scaling factors (SFs) in opposition to sample-dependent SFs were also discussed and the use of fitted SFs were shown to be the most cost-effective approach to predict IRMPD spectra. M06-2X/6-31+G(d,p) and CAM-B3LYP/aug-cc-pVDZ were also tested against the fingerprint region. Our results suggest that these methods can also be used for analysis in this lower frequency range and should be regarded as the methods of choice for cost-effective IRMPD simulations rather than the ubiquitous B3LYP method, especially when further molecular properties are needed.

Dataset Information

Theoretical Infrared Spectra: Quantitative Similarity Measures and Force Fields.

Publications

Theoretical Infrared Spectra: Quantitative Similarity Measures and Force Fields.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets