Abstract:Precisely quantifying the magnitude and direction of electric fields in proteins has long been an outstanding challenge in understanding biological functions. Nitrile vibrational Stark effect probes have been shown to be minimally disruptive to the protein structure and can be better direct reporters of local electrostatic field in the native state of a protein than other measures such as pKa shifts of titratable residues. However, interpretations of the connection between measured vibrational energy and electric field rely on the accurate molecular understanding of interactions of the nitrile group and its environment, particularly from hydrogen bonding. In this work, we compared the extent of hydrogen bonding calculated in two common force fields, the fixed charge force field Amber03 and polarizable force field AMOEBA, at 10 locations of cyanocysteine (CNC) in staphylococcal nuclease (SNase) against the experimental nitrile absorption frequency in terms of full width at half-maximum (FWHM) and frequency temperature line slope (FTLS). We observed that the number of hydrogen bonds correlated well in AMOEBA trajectories with respect to both the FWHM (r = 0.88) and the FTLS (r = -0.85), whereas the correlation of Amber03 trajectories was less reliable because the Amber03 force field predicted more hydrogen bonds in some mutants. Moreover, we demonstrated that contributions from the interactions between CNC and nearby water molecules were significant in AMOEBA trajectories but were not predicted by Amber03. We conclude that although the nitrile absorption peak shape could be qualitatively predicted by the fixed charge Amber03 force field, the detailed electrostatic environment measured by the nitrile probe in terms of the extent of hydrogen bonding could only be accurately observed in the AMOEBA trajectories, where the permanent dipole, quadrupole, and dipole-induced-dipole polarizable interactions were all taken into account. The significance of this finding to the goal of accurately predicting electric fields in complex biomolecular environments is discussed.

Geometry and Excitation Energy Fluctuations of NMA in Aqueous Solution with CHARMM, AMBER, OPLS, and GROMOS Force Fields: Implications for Protein Ultraviolet Spectra Simulation.

Ultraviolet Spectroscopy of Protein Backbone Transitions in Aqueous Solution: Combined QM and MM Simulations.

Multipolar Force Fields for Amide-I Spectroscopy from Conformational Dynamics of the Alanine-Trimer

A neural network protocol for electronic excitations of N-methylacetamide.

Comprehensive analysis of nitrile probe IR shifts and intensities in proteins: experiment and critical evaluation of simulations

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model.

Molecular Dynamics Simulations of N-methylacetamide (NMA) in Water by the ABEEM/MM Model

AFNMR: Automated Fragmentation Quantum Mechanical Calculation of NMR Chemical Shifts for Biomolecules

Hybrid molecular dynamics-quantum mechanics simulations of solute spectral properties in the condensed phase: evaluation of simulation parameters

Fragment Quantum Chemical Approach to Geometry Optimization and Vibrational Spectrum Calculation of Proteins

Lack of the ESIPT Band of Aromatic Ortho-Aminoaldehyde Derivatives Triggered by N-H Vibration

Vibrational Spectroscopy of N3- in the Gas and Condensed Phase

AMOEBA Force Field Predicts Accurate Hydrogen Bond Counts of Nitriles in SNase by Revealing Water-Protein Interaction in Vibrational Absorption Frequencies

Anharmonic Vibrational Analysis of Biomolecules and Solvated Molecules Using Hybrid QM/MM Computations

A neural network protocol for electronic excitations of <i>N</i> -methylacetamide

First-Principles Simulation of Amide and Aromatic Side Chain Ultraviolet Spectroscopy of A Cyclic Dipeptide

Two-Dimensional Electronic Correlation Spectroscopy of the N Pi* and Pi Pi* Protein Backbone Transitions: A Simulation Study

Vibrational circular dichroism spectroscopy with a classical polarizable force field: alanine in the gas and condensed phases

X-ray Absorption Spectra for Aqueous Ammonia and Ammonium: Quantum Mechanical versus Molecular Mechanical Embedding Schemes

Time-dependent density functional investigation on electronic spectra of 4′-N-dimethylamino-3-hydroxyflavone