Infrared line shape of an alpha-carbon deuterium-labeled amino acid

Activation of Ammonia by a Carbene-Stabilized Dithiolene Zwitterion

A carbene-stabilized dithiolene zwitterion (3) activates ammonia, affording 4^• and 5, through both single-electron transfer (SET) and hydrogen atom transfer (HAT). Reaction products were characterized spectroscopically and by single-crystal X-ray diffraction. The mechanism of the formation of 4^• and 5 was probed by experimental and computational methods. This discovery is the first example of metal-free ammonia activation via HAT.

Direct observation of peptide hydrogel self-assembly

The characterization of self-assembling molecules presents significant experimental challenges, especially when associated with phase separation or precipitation. Transparent window infrared (IR) spectroscopy leverages site-specific probes that absorb in the “transparent window” region of the biomolecular IR spectrum. Carbon–deuterium (C–D) bonds are especially compelling transparent window probes since they are non-perturbative, can be readily introduced site selectively into peptides and proteins, and their stretch frequencies are sensitive to changes in the local molecular environment. Importantly, IR spectroscopy can be applied to a wide range of molecular samples regardless of solubility or physical state, making it an ideal technique for addressing the solubility challenges presented by self-assembling molecules. Here, we present the first continuous observation of transparent window probes following stopped-flow initiation. To demonstrate utility in a self-assembling system, we selected the MAX1 peptide hydrogel, a biocompatible material that has significant promise for use in drug delivery and medical applications. C–D labeled valine was synthetically introduced into five distinct positions of the twenty-residue MAX1 β-hairpin peptide. Consistent with current structural models, steady-state IR absorption frequencies and linewidths of C–D bonds at all labeled positions indicate that these side chains occupy a hydrophobic region of the hydrogel and that the motion of side chains located in the middle of the hairpin is more restricted than those located on the hairpin ends. Following a rapid change in ionic strength to initiate self-assembly, the peptide absorption spectra were monitored as function of time, allowing determination of site-specific time constants. We find that within the experimental resolution, MAX1 self-assembly occurs as a cooperative process. These studies suggest that stopped-flow transparent window FTIR can be extended to other time-resolved applications, such as protein folding and enzyme kinetics.

To facilitate the characterization of phase-transitioning molecules, site-specific non-perturbative infrared probes are leveraged for continuous observation of the self-assembly of fibrils in a peptide hydrogel following stopped-flow initiation.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Synthesis and Evaluation of the Sensitivity and Vibrational Lifetimes of Thiocyanate and Selenocyanate Infrared Reporters

Two novel 2'-deoxyadenosine (dA) analogues, Si₂-dA-SCN and Si₂-dA-SeCN, and two novel phenylalanine (Phe) analogues, Boc-Me-PheCH₂SCN and Boc-Me-PheCH₂SeCN, have been synthesized and the thiocyanate (SCN) and selenocyanate (SeCN) functional groups evaluated as vibrational reporters. The syntheses of Si₂-dA-SCN and Si₂-dA-SeCN were accomplished in three steps in 16% and 32% overall yields, respectively, and the syntheses of Boc-Me-PheCH₂SCN and Boc-Me-PheCH₂SeCN were completed in four steps in 8.9% and 2.3% overall yields, respectively. The SCN and SeCN stretch vibrational modes were shown to be sensitive to the local environment by frequency shifts and full-width half-maximum (fwhm) changes in response to tetrahydrofuran (THF) and THF/water solvent mixtures. The vibrational lifetimes of the Si₂-dA-SeCN (237±12 ps) and Boc-Me-PheCH₂SeCN (295±31 ps) in THF solution were determined by ultrafast infrared pump-probe spectroscopy to be 1.5 to 3 times longer than those for Si₂-dA-SCN (140±6 ps) and Boc-Me-PheCH₂SCN (102±4 ps). The longer lifetimes for the SeCN analogues were attributed to the better insulating effects of the heavier selenium atom compared to the sulfur atom. The solvent sensitivity and longer vibrational lifetimes compared to other vibrational reporters suggest that SCN and SeCN vibrational reporters are well suited to studying several dynamic processes including protein and nucleic acid hydration and conformational changes, however stability issues may require post-synthetic modification methods to incorporate these reporters into biomacromolecules.

Asymmetric environmental effects on the structure and vibrations of cis-[Pt(NH3)2Cl2] in condensed phases

We report the structural and vibrational properties of anticancer drug cisplatin (cis-[Pt(NH3)2Cl2]) in gas phase, in solid phase, and in aqueous solution using density functional theory (DFT) calculations, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics, and effective normal modes analysis. In contrast with the gas-phase case, asymmetric hydrogen bonding environments are found in both solid phase and aqueous solution. It is shown that the discrepancy of the molecular geometry between previous gas phase calculations and the X-ray crystal structure can be resolved by considering intermolecular hydrogen bonds in the calculations of solid phase. In addition, our simulations in solid phase and aqueous solution reveal that asymmetric environmental effects lead to several spectral features observed in experiments, such as the blue-shift in the N-H stretching region and the frequency splitting of NH3 symmetric deformation modes. Furthermore, a similar decoupling and localization of several vibrational modes of cisplatin is found in solid phase and aqueous solution, in comparison to those of O-H stretching modes of water molecules in liquid water [ J. Phys. Chem. Lett. 2013 , 4 ( 19 ), 3245 - 3250 ].

Carbon-deuterium bonds as non-perturbative infrared probes of protein dynamics, electrostatics, heterogeneity, and folding

Vibrational spectroscopy is uniquely able to characterize protein dynamics and microenvironmental heterogeneity because it possesses an inherently high temporal resolution and employs probes of ultimately high structural resolution-the bonds themselves. The use of carbon-deuterium (C-D) bonds as vibrational labels circumvents the spectral congestion that otherwise precludes the use of vibrational spectroscopy to proteins and makes the observation of single vibrations within a protein possible while being wholly non-perturbative. Thus, C-D probes can be used to site-specifically characterize conformational heterogeneity and thermodynamic stability. C-D probes are also uniquely useful in characterizing the electrostatic microenvironment experienced by a specific residue side chain or backbone due to its effect on the C-D absorption frequency. In this chapter we describe the experimental procedures required to use C-D bonds and FT IR spectroscopy to characterize protein dynamics, structural and electrostatic heterogeneity, ligand binding, and folding.

Computational infrared and two-dimensional infrared photon echo spectroscopy of both wild-type and double mutant myoglobin-CO proteins

The CO stretching mode of both wild-type and double mutant ( T67R / S92D ) MbCO (carbonmonoxymyoglobin) proteins is an ideal infrared (IR) probe for studying the local electrostatic environment inside the myoglobin heme pocket. Recently, to elucidate the conformational switching dynamics between two distinguishable states, extensive IR absorption, IR pump-probe, and two-dimensional (2D) IR spectroscopic studies for various mutant MbCO's have been performed by the Fayer group. They showed that the 2D IR spectroscopy of the double mutant, which has a peroxidase enzyme activity, reveals a rapid chemical exchange between two distinct states, whereas that of the wild-type does not. Despite the fact that a few simulation studies on these systems were already performed and reported, such complicated experimental results have not been fully reproduced nor described in terms of conformational state-to-state transition processes. Here, we first develop a distributed vibrational solvatochromic charge model for describing the CO stretch frequency shift reflecting local electric potential changes. Then, by carrying out molecular dynamic simulations of the two MbCO's and examining their CO frequency trajectories, it becomes possible to identify a proper reaction coordinate consisting of His64 imidazole ring rotation and its distance to the CO ligand. From the 2D surfaces of the resulting potential of mean forces, the spectroscopically distinguished A1 and A3 states of the wild-type as well as two more substates of the double mutant are identified and their vibrational frequencies and distributions are separately examined. Our simulated IR absorption and 2D IR spectra of the two MbCO's are directly compared with the previous experimental results reported by the Fayer group. The chemical exchange rate constants extracted from the two-state kinetic analyses of the simulated 2D IR spectra are in excellent agreement with the experimental values. On the basis of the quantitative agreement between the simulated spectra and experimental ones, we further examine the conformational differences in the heme pockets of the two proteins and show that the double mutation, T67R / S92D , suppresses the A1 population, restricts the imidazole ring rotation, and increases hydrogen-bond strength between the imidazole Nε-H and the oxygen atom of the CO ligand. It is believed that such delicate change of distal His64 imidazole ring dynamics induced by the double mutation may be responsible for its enhanced peroxidase catalytic activity as compared to the wild-type myoglobin.

Vibrational solvatochromism and electrochromism. II. Multipole analysis

Small infrared probe molecules have been widely used to study local electrostatic environment in solutions and proteins. Using a variety of time- and frequency-resolved vibrational spectroscopic methods, one can accurately measure the solvation-induced vibrational frequency shifts and the timescales and amplitudes of frequency fluctuations of such IR probes. Since the corresponding frequency shifts are directly related to the local electric field and its spatial derivatives of the surrounding solvent molecules or amino acids in proteins, one can extract information on local electric field around an IR probe directly from the vibrational spectroscopic results. Here, we show that, carrying out a multipole analysis of the solvatochromic frequency shift, the solvatochromic dipole contribution to the frequency shift is not always the dominant factor. In the cases of the nitrile-, thiocyanato-, and azido-derivatized molecules, the solvatochromic quadrupole contributions to the corresponding stretch mode frequency shifts are particularly large and often comparable to the solvatochromic dipole contributions. Noting that the higher multipole moment-solvent electric field interactions are short range effects in comparison to the dipole interaction, the H-bonding interaction-induced vibrational frequency shift can be caused by such short-range multipole-field interaction effects. We anticipate that the present multipole analysis method specifically developed to describe the solvatochromic vibrational frequency shifts will be useful to understand the intermolecular interaction-induced vibrational property changes and to find out a relationship between vibrational solvatochromism and electrochromism of IR probes in condensed phases.

Monitoring Intramolecular Proton Transfer with Two-Dimensional Infrared Spectroscopy: A Computational Prediction

Proton transfer processes are ubiquitous and play a vital role in a broad range of chemical and biochemical phenomena. The ability of two-dimensional infrared (2D IR) spectroscopy with a carbon-deuterium (C-D) reporter to monitor the kinetics of proton transfer in the model compound malonaldehyde was demonstrated computationally. One of the two carbonyl/enol carbon atoms in malonaldehyde was labeled with a C-D bond. The C-D stretch vibrational frequency provides ∼150 cm(-1) of sensitivity to the two tautomers of malonaldehyde. Mixed quantum mechanics/molecular mechanics simulations employing the self-consistent-charge density functional tight binding (SCC-DFTB) method were used to compute 2D IR line shapes for the C-D stretch of labeled malonaldehyde in aqueous solution. The 2D IR spectra reveal cross peaks from the chemical exchange of the proton. The kinetics for the growth of the cross-peaks (and the decay of the diagonal peaks) precisely match the proton transfer rate observed in the SCC-DFTB simulations.

Azido Homoalanine is a Useful Infrared Probe for Monitoring Local Electrostatistics and Sidechain Solvation in Proteins

The use of IR probes to monitor protein structure, deduce local electric field, and investigate the mechanism of enzyme catalysis and protein folding has attracted increasing attention. Here, the azidohomoalanine (Aha) is considered as a useful IR probe. The intricate details of the distinct effects of backbone peptide bonds and H-bonded water molecules on the azido stretch mode of the IR probe Aha were revealed by carrying out QM/MM MD simulations of two variants of the protein NTL9, NTL9-Met1Aha and NTL9-Ile4Aha and comparing the resulting simulated IR spectra with experiments.

Two-dimensional infrared spectroscopy of azido-nicotinamide adenine dinucleotide in water

Mid-IR active analogs of enzyme cofactors have the potential to be important spectroscopic reporters of enzyme active site dynamics. Azido-nicotinamide adenine dinucleotide (NAD(+)), which has been recently synthesized in our laboratory, is a mid-IR active analog of NAD(+), a ubiquitous redox cofactor in biology. In this study, we measure the frequency-frequency time correlation function for the antisymmetric stretching vibration of the azido group of azido-NAD(+) in water. Our results are consistent with previous studies of pseudohalides in water. We conclude that azido-NAD(+) is sensitive to local environmental fluctuations, which, in water, are dominated by hydrogen-bond dynamics of the water molecules around the probe. Our results demonstrate the potential of azido-NAD(+) as a vibrational probe and illustrate the potential of substituted NAD(+)-analogs as reporters of local structural dynamics that could be used for studies of protein dynamics in NAD-dependent enzymes.

Modulating Accidental Fermi Resonance: What a Difference a Neutron Makes

Vibrational reporters have shown significant promise as sensitive probes of local environments in proteins and nucleic acids. The utility of two potential vibrational probes, the cyanate and azide groups in phenyl cyanate and 3-azidopyridine, respectively, has been hindered by accidental Fermi resonance. Anharmonic coupling, between the fundamental -OCN or -N(3) asymmetric stretch vibration with a near resonant combination band, results in an extremely broad and complex absorption profile for each of these probes. A total of eight phenyl cyanate and six 3-azidopyridine isotopomers were synthesized and studied. Isotopic editing effectively modulated the accidental Fermi resonance - the absorption profiles of several isotopomers were greatly simplified while others remained complex. The origins of the observed profiles are discussed. Addition of a single neutron to the middle atom of the oscillator converted the absorption profile to essentially a single band resulting from either the cyanate or azide asymmetric stretch vibration.

A direct comparison of azide and nitrile vibrational probes

The synthesis of 2'-azido-5-cyano-2'-deoxyuridine, N(3)CNdU (1), from trityl-protected 2'-amino-2'-deoxyuridine was accomplished in four steps with a 12.5% overall yield. The IR absorption positions and profiles of the azide and nitrile group of N(3)CNdU were investigated in 14 different solvents and water/DMSO solvent mixtures. The azide probe was superior to the nitrile probe in terms of its extinction coefficient, which is 2-4 times larger. However, the nitrile IR absorbance profile is generally less complicated by accidental Fermi resonance. The IR frequencies of both probes undergo a substantial red shift upon going from water to aprotic solvents such as THF or DMSO. DFT calculations supported the hypothesis that the molecular origin of the higher observed frequency in water is primarily due to hydrogen bonds between the probes and water molecules.

Car-Parrinello simulation of hydrogen bond dynamics in sodium hydrogen bissulfate

We studied proton dynamics of a short hydrogen bond of the crystalline sodium hydrogen bissulfate, a hydrogen-bonded ferroelectric system. Our approach was based on the established Car-Parrinello molecular dynamics (CPMD) methodology, followed by an a posteriori quantization of the OH stretching motion. The latter approach is based on snapshot structures taken from CPMD trajectory, calculation of proton potentials, and solving of the vibrational Schrodinger equation for each of the snapshot potentials. The so obtained contour of the OH stretching band has the center of gravity at about 1540 cm(-1) and a half width of about 700 cm(-1), which is in qualitative agreement with the experimental infrared spectrum. The corresponding values for the deuterated form are 1092 and 600 cm(-1), respectively. The hydrogen probability densities obtained by solving the vibrational Schrodinger equation allow for the evaluation of potential of mean force along the proton transfer coordinate. We demonstrate that for the present system the free energy profile is of the single-well type and features a broad and shallow minimum near the center of the hydrogen bond, allowing for frequent and barrierless proton (or deuteron) jumps. All the calculated time-averaged geometric parameters were in reasonable agreement with the experimental neutron diffraction data. As the present methodology for quantization of proton motion is applicable to a variety of hydrogen-bonded systems, it is promising for potential use in computational enzymology.

Sum frequency generation and solid-state NMR study of the structure, orientation, and dynamics of polystyrene-adsorbed peptides

The power of combining sum frequency generation (SFG) vibrational spectroscopy and solid-state nuclear magnetic resonance (ssNMR) spectroscopy to quantify, with site specificity and atomic resolution, the orientation and dynamics of side chains in synthetic model peptides adsorbed onto polystyrene (PS) surfaces is demonstrated in this study. Although isotopic labeling has long been used in ssNMR studies to site-specifically probe the structure and dynamics of biomolecules, the potential of SFG to probe side chain orientation in isotopically labeled surface-adsorbed peptides and proteins remains largely unexplored. The 14 amino acid leucine-lysine peptide studied in this work is known to form an alpha-helical secondary structure at liquid-solid interfaces. Selective, individual deuteration of the isopropyl group in each leucine residue was used to probe the orientation and dynamics of each individual leucine side chain of LKalpha14 adsorbed onto PS. The selective isotopic labeling methods allowed SFG analysis to determine the orientations of individual side chains in adsorbed peptides. Side chain dynamics were obtained by fitting the deuterium ssNMR line shape to specific motional models. Through the combined use of SFG and ssNMR, the dynamic trends observed for individual side chains by ssNMR have been correlated with side chain orientation relative to the PS surface as determined by SFG. This combination provides a more complete and quantitative picture of the structure, orientation, and dynamics of these surface-adsorbed peptides than could be obtained if either technique were used separately.

Conformation dependence of the C(alpha)D(alpha) stretch mode in peptides. II. explicitly hydrated alanine peptide structures

Our previous studies of the potential utility of the C(alpha)D(alpha) stretch frequency, nu(CD), as a tool for determining conformation in peptide systems (Mirkin and Krimm, J Phys Chem A 2004, 108, 10923-10924; 2007, 111, 5300-5303) dealt with the spectroscopic characteristics of isolated alanine peptides with alpha(R), beta, and polyproline II structures. We have now extended these ab initio calculations to include various explicit-water environments interacting with such conformers. We find that the structure-discriminating feature of this technique is in fact enhanced as a result of the conformation-specific interactions of the bonding waters, in part due to our finding (Mirkin and Krimm, J Phys Chem B 2008, 112, 15268) that C(alpha)--D(alpha)...O(water) hydrogen bonds can be present in addition to those expected between water and the CO and NH of the peptide groups. In fact, nu(CD) is hardly affected by the latter bonding but can be shifted by up to 70 cm(-1) by the former hydrogen bonds. We also discuss the factors that will have to be considered in developing the molecular dynamics (MD) treatment needed to satisfactorily take account of the influence of outer water layers on the structure of the first-layer water molecules that hydrogen bond to the peptide backbone.

Carbon-deuterium vibrational probes of peptide conformation: alanine dipeptide and glycine dipeptide

The utility of alpha-carbon deuterium-labeled bonds (C(alpha)-D) as infrared reporters of local peptide conformation was investigated for two model dipeptide compounds: C(alpha)-D labeled alanine dipeptide (Adp-d(1)) and C(alpha)-D(2) labeled glycine dipeptide (Gdp-d(2)). These model compounds adopt structures that are analogous to the motifs found in larger peptides and proteins. For both Adp-d(1) and Gdp-d(2), we systematically mapped the entire conformational landscape in the gas phase by optimizing the geometry of the molecule with the values of phi and psi, the two dihedral angles that are typically used to characterize the backbone structure of peptides and proteins, held fixed on a uniform grid with 7.5 degrees spacing. Since the conformations were not generally stationary states in the gas phase, we then calculated anharmonic C(alpha)-D and C(alpha)-D(2) stretch transition frequencies for each structure. For Adp-d(1) the C(alpha)-D stretch frequency exhibited a maximum variability of 39.4 cm(-1) between the six stable structures identified in the gas phase. The C(alpha)-D(2) frequencies of Gdp-d(2) show an even more substantial difference between its three stable conformations: there is a 40.7 cm(-1) maximum difference in the symmetric C(alpha)-D(2) stretch frequencies and an 81.3 cm(-1) maximum difference in the asymmetric C(alpha)-D(2) stretch frequencies. Moreover, the splitting between the symmetric and asymmetric C(alpha)-D(2) stretch frequencies of Gdp-d(2) is remarkably sensitive to its conformation.

C-D modes of deuterated side chain of leucine as structural reporters via dual-frequency two-dimensional infrared spectroscopy

Perdeuteration of the side chains of amino acids such as leucine results in appearance of reasonably strong absorption peaks around 2050-2220 cm(-1) that belong to the CD stretching modes and exhibit extinction coefficients of up to 120 M(-1) cm(-1). The properties of the CD stretching transitions in leucine-d(10) as structural labels are studied via the methods of two-dimensional infrared (2DIR) spectroscopy. The cross peaks caused by interactions of the CD stretching modes with amide I (Am-I), CO, and amide II (Am-II) modes are obtained by the dual-frequency 2DIR method. The CD stretching peaks in leucine-d(10) are characterized using DFT computational modeling as well as relaxation-assisted 2DIR (RA 2DIR) measurements. The RA 2DIR measurements showed different enhancements and different energy transport times (arrival times) for different CD/Am-II and CD/CO cross peaks; a correlation between the intermode distance, the arrival time, and the amplification factor is reported. We demonstrated that the CD transitions of leucine-d(10) amino acid can serve as convenient structural reporters via the dual-frequency 2DIR methods and discussed the potential of leucine-d(10) and other amino acids with deuterium-labeled side chains as probes of protein structure and dynamics.

Efforts toward the direct experimental characterization of enzyme microenvironments: tyrosine100 in dihydrofolate reductase

State secrets: Site-specific deuteration and FTIR studies reveal that Tyr100 in dihydrofolate reductase plays an important role in catalysis, with a strong electrostatic coupling occurring between Tyr100 and the charge that develops in the hydride-transfer transition state (see picture, NADP(+) purple, Tyr100 green). However, relaying correlated motions that facilitate catalysis from distal sites of the protein to the hydride donor may also be involved.

The determinants of stability and folding in evolutionarily diverged cytochromes c

Cytochrome c has served as a paradigm for the study of protein stability, folding, and molecular evolution, but it remains unclear how these aspects of the protein are related. For example, while the bovine and equine cytochromes c are known to have different stabilities, and possibly different folding mechanisms, it is not known how these differences arise from just three amino acid substitutions introduced during divergence. Using site-selectively incorporated carbon-deuterium bonds, we show that like the equine protein, bovine cytochrome c is induced to unfold by guanidine hydrochloride via a stepwise mechanism, but it does not populate an intermediate as is observed with the equine protein. The increased stability also results in more similar free energies of unfolding observed at different sites within the protein, giving the appearance of a more concerted mechanism. Furthermore, we show that the differences in stability and folding appear to result from a single amino acid substitution that stabilizes a helix by allowing for increased solvation of its N-terminus.

Structure and Raman Spectrum of Clavulanic Acid in Aqueous Solution

The calculation of the vibrational Raman spectrum of enzyme-bound beta-lactamase inhibitors may be of help to understand the mechanisms responsible for bacterial drug resistance. Here, we present a study of the solvation structure and the vibrational properties of clavulanate, an important beta-lactamase inhibitor, in aqueous solution as obtained from full quantum and hybrid empirical/quantum molecular dynamics simulations at ambient conditions. The analysis of the vibrational density of states indicates that hybrid empirical/quantum mechanical simulations are able to properly describe the vibrational levels of clavulanate in solution. In addition, we propose a computationally efficient protocol to calculate the vibrational Raman effect for large solute molecules in water, which is able to faithfully reproduce the experimentally recorded clavulanate Raman spectrum and discloses the possibility to employ hybrid simulations to assign the experimental Raman spectra of inhibitors bound to beta-lactamases.