Relaxation and anharmonic couplings of the O-H stretching vibration of asymmetric strongly hydrogen-bonded complexes

Population Dynamics of Stretching Excitations of p-Azido-phenylalanine Incorporated in Calmodulin-Peptide Complexes

We genetically incorporated the unnatural amino acid p-azido-phenylalanine (AzF) into the ubiquitous Ca²⁺ sensor protein calmodulin (CaM) in complex with different peptides to explore the response of the azido stretching line shape to varying binding motifs with femtosecond infrared spectroscopy. The dynamic response of the azido stretching mode varies in different CaM-peptide complexes. We model these dynamics as coherent excitations of Fermi resonances and extract a lifetime of the azido stretching vibration of about 1 ps. The resulting model parameters are commensurate with the linear infrared absorption lineshapes which suggests that the conformation-sensitive vibrational lineshape could be composed of Fermi resonances that differ between the protein-peptide complexes.

Composition-Dependent Hydrogen-Bonding Motifs and Dynamics in Brønsted Acid-Base Mixtures

In recent years the interaction of organophosphates and imines, which is at the core of Brønsted acid organocatalysis, has been established to be based on strong ionic hydrogen bonds. Yet, besides the formation of homodimers consisting of two acid molecules and heterodimers consisting of one acid and one base, also multimeric molecular aggregates are formed in solution. These multimeric aggregates consist of one base and several acid molecules. The details of the intermolecular bonding in such aggregates, however, have remained elusive. To characterize composition-dependent bonding and bonding dynamics in these aggregates, we use linear and nonlinear infrared (IR) spectroscopy at varying molar ratios of diphenyl phosphoric acid and quinaldine. We identify the individual aggregate species, giving rise to the structured, strong, and very broad infrared absorptions, which span more than 1000 cm^–1. Linear infrared spectra and density functional theory calculations of the proton transfer potential show that doubly ionic intermolecular hydrogen bonds between the acid and the base lead to absorptions which peak at ∼2040 cm^–1. The contribution of singly ionic hydrogen bonds between an acid anion and an acid molecule is observed at higher frequencies. As common to such strong hydrogen bonds, ultrafast IR spectroscopy reveals rapid, ∼ 100 fs, dissipation of energy from the proton transfer coordinate. Yet, the full dissipation of the excess energy occurs on a ∼0.8–1.1 ps time scale, which becomes longer when multimers dominate. Our results thus demonstrate the coupling and collectivity of the hydrogen bonds within these complexes, which enable efficient energy transfer.

Effect of intramolecular hydrogen-bond formation on the molecular conformation of amino acids

The molecular conformation of the carboxyl group can be crucial for its chemical properties and intermolecular interactions, especially in complex molecular environments such as polypeptides. Here, we study the conformational behaviour of the model amino acid N-acetylproline in solution at room temperature with two-dimensional infrared spectroscopy. We find that the carboxyl group of N-acetylproline adopts two distinct conformations, syn- and anti-. In the syn-conformer the O-H group is oriented at  ~60^∘ with respect to the C=O and in the anti-conformer the O-H is anti-parallel to the C=O. In hydrogen-bond accepting solvents such as dimethyl sulfoxide or water, we observe that, similar to simple carboxylic acids, around 20% of the -COOH groups adopt an anti-conformation. However, when N-acetylproline is dissolved in a weakly hydrogen-bond accepting solvent (acetonitrile), we observe the formation of a strong intramolecular hydrogen bond between the carboxyl group in the anti-conformation and the amide group, which stabilizes the anti-conformer, increasing its relative abundance to ~60%.

Peptide Side-COOH Groups Have Two Distinct Conformations under Biorelevant Conditions

The carboxyl (COOH) side chain groups of amino acids, such as aspartic acid, play an important role in biochemical processes, including enzymatic proton transport. In many theoretical studies, it was found that the (bio)chemical reactivity of the carboxyl group strongly depends on the conformation of this group. Interestingly, up to now there has been no experimental investigation of the geometry and the stability of different COOH conformers under biorelevant conditions. Here, we investigate the conformational isomerism of the side chain COOH group of N-acetyl aspartic acid amide using polarization-resolved two-dimensional infrared spectroscopy. We find that the carboxyl group shows two distinct near-planar conformers (syn and anti) when dissolved in water at room temperature. Both conformers are significantly populated in aqueous solution (75 ± 10% and 25 ± 10% for syn and anti, respectively). Molecular dynamics simulations show that the anti conformer interacts more strongly with water molecules than the syn conformer, explaining why this conformer is significantly present in aqueous solution.

Oscillatory Active-site Motions Correlate with Kinetic Isotope Effects in Formate Dehydrogenase

Thermal motions of enzymes have been invoked to explain the temperature dependence of kinetic isotope effects (KIE) in enzyme-catalyzed hydride transfers. Formate dehydrogenase (FDH) from Candida boidinii exhibits a temperature independent KIE that becomes temperature dependent upon mutation of hydrophobic residues in the active site. Ternary complexes of FDH that mimic the transition state structure allow investigation of how these mutations influence active-site dynamics. A combination of X-ray crystallography, two-dimensional infrared (2D IR) spectroscopy, and molecular dynamic simulations characterize the structure and dynamics of the active site. FDH exhibits oscillatory frequency fluctuations on the picosecond timescale, and the amplitude of these fluctuations correlates with the temperature dependence of the KIE. Both the kinetic and dynamic phenomena can be reproduced computationally. These results provide experimental evidence for a connection between the temperature dependence of KIEs and motions of the active site in an enzyme-catalyzed reaction consistent with activated tunneling models of the hydride transfer reaction.

Interpreting Quasi-Thermal Effects in Ultrafast Spectroscopy of Hydrogen-Bonded Systems

Vibrational excitation of molecules in the condensed phase relaxes through vibrational modes of decreasing energy to ultimately generate an equilibrium state in which the energy is distributed among low-frequency modes. In ultrafast vibrational spectroscopy, changes in the vibrational features of hydrogen-bonded NH and OH stretch modes are typically observed to persist long after these high-frequency vibrations have relaxed. Due to the resemblance to the spectral changes caused by heating the sample, these features are typically described as arising from a hot ground state. However, these spectral features appear on ultrafast time scales that are much too fast to result from a true thermal state, and significant differences between the thermal difference spectrum and the induced quasi-thermal changes in ultrafast spectroscopy are often observed. Here, we examine and directly compare the thermal and quasi-thermal responses of the hydrogen-bonded homodimer of 7-azaindole with temperature-dependent FTIR spectroscopy and ultrafast mid-IR continuum spectroscopy. We find that the thermal difference spectra contain contributions from both dissociation of the hydrogen bonds and from frequency shifts due to changes in the thermal population of low-frequency modes. The transient spectra in ultrafast vibrational spectroscopy are also found to contain two contributions: initial frequency shifts over 2.3 ± 0.11 ps associated with equilibration of the initial excitation, and frequency shifts associated with the excitation of several fingerprint modes, which decay over 21.8 ± 0.11 ps, giving rise to a quasi-thermal response caused by a distribution of fingerprint modes being excited within the sample ensemble. This resembles the thermal frequency shifts due to population changes of low-frequency modes, but not the overall thermal spectrum, which is dominated by features caused by dimer dissociation. These findings provide insight into the changes in the vibrational spectrum from different origins and are important for assigning, analyzing, and comparing features in thermal and ultrafast vibrational spectroscopy of hydrogen-bonded complexes.

Investigation of hydrogen bonds properties in the terephthalic acid crystal, using molecular dynamics method

The aim of this study was to perform calculations using the method of Car-Parrinello molecular dynamics, leading to the optimized geometry of the molecules of 1,4-benzenedicarboxylic acid (terephthalic acid) in crystals, for the hydrogen form and three variants of substitution of deuterium atoms inside a carboxyl group. Based on the results, trajectories and dipole moments were calculated, what makes possible to simulate vibrations in different systems, and to make calculation of theoretical infrared spectra and atomic power spectra. Theoretical results were compared with the experimental spectra, which verifies the correctness of the method and also was compared with the results obtained by quantum-mechanical calculations using DFT for the isolated dimer. Comparison of the spectra of different forms, allowed for in-depth analysis of the effect of isotopic substitution on the frequency of vibrations and shapes of bands, and confirm the presence of possible coupling effects and intra- and intermolecular interactions. Comparison with the DFT results for the dimer show influence of the crystal structure on the spectra.

Signatures of nanoconfinement on the linear and nonlinear vibrational spectroscopy of a model hydrogen-bonded complex dissolved in a polar solvent

The one-dimensional IR (1D-IR) absorption and IR pump-probe spectra of a hydrogen stretch in a model hydrogen-bonded complex dissolved in a polar solvent confined in spherical hydrophobic cavities of different sizes were simulated using ground-state mixed quantum-classical dynamics. Due to a thorough analysis of key properties of the complex and solvent from equilibrium trajectory data, we were able to gain insight into the microscopic details underlying the spectra. Both the 1D-IR and IR pump-probe spectra manifested the effects of confinement on the relative stabilities of the covalent and ionic forms of the complex through pronounced changes in their peak intensities and numbers. However, in contrast to the 1D-IR spectra, the time-resolved pump-probe spectra were found to be uniquely sensitive to the changes in the molecular dynamics as the cavity size is varied. In particular, it was found that the variations in the time evolutions of the peak intensities in the pump-probe spectra reflect the differences in the solvation dynamics associated with the various forms of the complex in different locations within the cavities. The ability to detect these differences underscores the advantage of using pump-probe spectroscopy for studying nanoconfined systems.

Characterizing the dynamics of functionally relevant complexes of formate dehydrogenase

The potential for femtosecond to picosecond time-scale motions to influence the rate of the intrinsic chemical step in enzyme-catalyzed reactions is a source of significant controversy. Among the central challenges in resolving this controversy is the difficulty of experimentally characterizing thermally activated motions at this time scale in functionally relevant enzyme complexes. We report a series of measurements to address this problem using two-dimensional infrared spectroscopy to characterize the time scales of active-site motions in complexes of formate dehydrogenase with the transition-state-analog inhibitor azide (N(3)(-)). We observe that the frequency-frequency time correlation functions (FFCF) for the ternary complexes with NAD(+) and NADH decay completely with slow time constants of 3.2 ps and 4.6 ps, respectively. This result suggests that in the vicinity of the transition state, the active-site enzyme structure samples a narrow and relatively rigid conformational distribution indicating that the transition-state structure is well organized for the reaction. In contrast, for the binary complex, we observe a significant static contribution to the FFCF similar to what is seen in other enzymes, indicating the presence of the slow motions that occur on time scales longer than our measurement window.

Theoretical modeling of infrared spectra of the hydrogen and deuterium bond in aspirin crystal

An extended quantum theoretical approach of the nu(X-H) IR lineshape of cyclic dimers of weakly H-bonded species is proposed. We have extended a previous approach [M.E.-A. Benmalti, P. Blaise, H.T. Flakus, O. Henri-Rousseau, Chem. Phys. 320 (2006) 267] by accounting for the anharmonicity of the slow mode which is described by a "Morse" potential in order to reproduce the polarized infrared spectra of the hydrogen and deuterium bond in acetylsalicylic acid (aspirin) crystals. From comparison of polarized IR spectra of isotopically neat and isotopically diluted aspirin crystals it resulted that centrosymmetric aspirin dimer was the bearer of the crystal main spectral properties. In this approach, the adiabatic approximation is performed for each separate H-bond bridge of the dimer and a strong non-adiabatic correction is introduced into the model via the resonant exchange between the fast mode excited states of the two moieties. Within the strong anharmonic coupling theory, according to which the X-H...Y high-frequency mode is anharmonically coupled to the H-bond bridge, this model incorporated the Davydov coupling between the excited states of the two moieties, the quantum direct and indirect dampings and the anharmonicity for the H-bond bridge. The spectral density is obtained within the linear response theory by Fourier transform of the damped autocorrelation functions. The evaluated spectra are in fairly good agreement with the experimental ones by using a minimum number of independent parameters. The effect of deuteration has been well reproduced by reducing simply the angular frequency of the fast mode and the anharmonic coupling parameter.

Exploring the molecular origins of protein dynamics in the active site of human carbonic anhydrase II

We present three-pulse vibrational echo measurements of azide ion bound to the active site Zn of human carbonic anhydrase II (HCA II) and of two separate active-site mutants Thr199 --> Ala (T199A) and Leu198 --> Phe (L198F). Because structural motions of the protein active site influence the frequency of bound ligands, the differences in the time scales of the frequency-frequency correlation functions (FFCFs) obtained from global fits to each set of data allow us to make inferences about the time scales of the active site dynamics of HCA II. Surprisingly, the deletion of a potential electrostatic interaction in results in very little change in the FFCF, but the insertion of the bulky phenylalanine ring in causes much faster dynamics. We conclude that the fast, sub-picosecond time scale in the correlation function is attributable to hydrogen bond dynamics, and the slow, apparently static contribution is due to the conformational flexibility of Zn-bound azide in the active site.

Ultrafast N-H vibrational dynamics of cyclic doubly hydrogen-bonded homo- and heterodimers

Hydrogen-bonded interfaces are essential structural elements in biology. Furthermore, they can mediate electron transport by coupling the electron to proton transfer within the interface. The specific hydrogen-bonding configuration and strength have a large impact on the proton transfer, which exchanges the hydrogen-bonded donor and acceptor species (i.e., NH...O --> N...HO). Modulations of the hydrogen-bonding environment, such as the hydrogen-bond stretch and twist modes, affect the proton-transfer dynamics. Here, we present transient grating and echo peak shift measurements of the NH stretch vibrations of four doubly hydrogen-bonded cyclic dimers in their electronic ground state. The equilibrium vibrational dynamics exhibit strong coherent modulations that we attribute to coupling of the high-frequency NH vibration to the low-frequency interdimer stretch and twist modes and not to interference between multiple Fermi resonances that dominate the substructure of the linear spectra.

Coupling of complex aromatic ring vibrations to solvent through hydrogen bonds: effect of varied on-ring and off-ring hydrogen-bonding substitutions

In this study, we examine the coupling of a complex ring vibration to solvent through hydrogen-bonding interactions. We compare phenylalanine, tyrosine, l-dopa, dopamine, norepinephrine, epinephrine, and hydroxyl-dl-dopa, a group of physiologically important small molecules that vary by single differences in H-bonding substitution. By examination of the temperature dependence of infrared absorptions of these molecules, we show that complex, many-atom vibrations can be coupled to solvent through hydrogen bonds and that the extent of that coupling is dependent on the degree of both on- and off-ring H-bonding substitution. The coupling is seen as a temperature-dependent frequency shift in infrared spectra, but the determination of the physical origin of that shift is based on additional data from temperature-dependent optical experiments and ab initio calculations. The optical experiments show that these small molecules are most sensitive to their immediate H-bonding environment rather than to bulk solvent properties. Ab initio calculations demonstrate H-bond-mediated vibrational coupling for the system of interest and also show that the overall small molecule solvent dependence is determined by a complex interplay of specific interactions and bulk solvation characteristics. Our findings indicate that a full understanding of biomolecule vibrational properties must include consideration of explicit hydrogen-bonding interactions with the surrounding microenvironment.