Spectral signatures of heterogeneous protein ensembles revealed by MD Simulations of 2DIR spectra

The secondary structure of diatom silaffin peptide R5 determined by two-dimensional infrared spectroscopy

Diatoms, unicellular marine organisms, harness short peptide repeats of the protein silaffin to transform silicic acid into biosilica nanoparticles. This process has been a white whale for material scientists due to its potential in biomimetic applications, ranging from medical to microelectronic fields. Replicating diatom biosilicification will depend on a thorough understanding of the silaffin peptide structure during the reaction, yet existing models in the literature offer conflicting views on peptide folding during silicification. In our study, we employed two-dimensional infrared spectroscopy (2DIR) within the amide I region to determine the secondary structure of the silaffin repeat unit 5 (R5), both pre- and post-interaction with silica. The 2DIR experiments are complemented by molecular dynamics (MD) simulations of pure R5 reacting with silicate. Subsequently, theoretical 2DIR spectra calculated from these MD trajectories allowed us to compare calculated spectra with experimental data, and to determine the diverse structural poses of R5. Our findings indicate that unbound R5 predominantly forms β-strand structures alongside various atypical secondary structures. Post-silicification, there's a noticeable shift: a decrease in β-strands coupled with an increase in turn-type and bend-type configurations. We theorize that this structural transformation stems from silicate embedding within R5's hydrogen-bond network, prompting the peptide backbone to contract and adapt around the biosilica precursors.

Theoretical basis for interpreting heterodyne chirality-selective sum frequency generation spectra of water

Chirality-selective vibrational sum frequency generation (chiral SFG) spectroscopy has emerged as a powerful technique for the study of biomolecular hydration water due to its sensitivity to the induced chirality of the first hydration shell. Thus far, water O-H vibrational bands in phase-resolved heterodyne chiral SFG spectra have been fit using one Lorentzian function per vibrational band, and the resulting fit has been used to infer the underlying frequency distribution. Here, we show that this approach may not correctly reveal the structure and dynamics of hydration water. Our analysis illustrates that the chiral SFG responses of symmetric and asymmetric O-H stretch modes of water have opposite phase and equal magnitude and are separated in energy by intramolecular vibrational coupling and a heterogeneous environment. The sum of the symmetric and asymmetric responses implies that an O-H stretch in a heterodyne chiral SFG spectrum should appear as two peaks with opposite phase and equal amplitude. Using pairs of Lorentzian functions to fit water O-H stretch vibrational bands, we improve spectral fitting of previously acquired experimental spectra of model β-sheet proteins and reduce the number of free parameters. The fitting allows us to estimate the vibrational frequency distribution and thus reveals the molecular interactions of water in hydration shells of biomolecules directly from chiral SFG spectra.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

Interactions of potato-derived and human recombinant 5-lipoxygenase with sec-O-glucosylhamaudol by multi-spectroscopy and molecular docking

5-lipoxygenase (5-LOX) was a key enzyme involved in many inflammatory diseases. Sec-O-glucosylhamaudol (SOG) was a chromone found in Saposhnikovia divaricata (Turcz.) Schischk (S. divaricate). The potato-derived 5-LOX (p-5-LOX) and human recombinant 5-LOX (h-5-LOX) were selected as model protein due to their simple usability and high stability in this study. Thus, the binding interactions of p-5-LOX and h-5-LOX with SOG were investigated by multi-spectroscopy and molecular docking. As a result, the fluorescence intensities of the two 5-LOX were quenched statically by SOG. However, the binding ability of SOG to h-5-LOX was higher than that of p-5-LOX at the same temperature. The results of multi-spectroscopy revealed that the conformation and micro-environment of the two 5-LOX proteins were changed after binding with SOG. Fluorescence assay and molecular docking indicated that hydrogen bond and electrostatic gravitation were the main forces between the two 5-LOX and SOG. Our results here suggested that SOG may exert anti-inflammatory effect by inhibiting 5-LOX activity.

AIM: A Mapping Program for Infrared Spectroscopy of Proteins

Here, we present a new analysis program, AIM, that allows extracting the vibrational amide-I Hamiltonian using molecular dynamics trajectories for protein infrared spectroscopy modeling. The constructed Hamiltonians can be used as input for spectral calculations allowing the calculation of infrared absorption spectra, vibrational circular dichroism, and two-dimensional infrared spectra. These spectroscopies allow the study of the structure and dynamics of proteins. We will explain the essence of how AIM works and give examples of the information and spectra that can be obtained with the program using the Trypsin Inhibitor as an example. AIM is freely available from GitHub, and the package contains a demonstration allowing easy introduction to the use of the program.

Temperature-dependent local conformations and conformational distributions of cyanine dimer labeled single-stranded-double-stranded DNA junctions by 2D fluorescence spectroscopy

DNA replication and the related processes of genome expression require binding, assembly, and function of protein complexes at and near single-stranded (ss)-double-stranded (ds) DNA junctions. These central protein-DNA interactions are likely influenced by thermally induced conformational fluctuations of the DNA scaffold across an unknown distribution of functionally relevant states to provide regulatory proteins access to properly conformed DNA binding sites. Thus, characterizing the nature of conformational fluctuations and the associated structural disorder at ss-dsDNA junctions is critical for understanding the molecular mechanisms of these central biological processes. Here, we describe spectroscopic studies of model ss-dsDNA fork constructs that contain dimers of "internally labeled" cyanine (iCy3) chromophore probes that have been rigidly inserted within the sugar-phosphate backbones of the DNA strands. Our combined analyses of absorbance, circular dichroism, and two-dimensional fluorescence spectroscopy permit us to characterize the local conformational parameters and conformational distributions. We find that the DNA sugar-phosphate backbones undergo abrupt successive changes in their local conformations-initially from a right-handed and ordered DNA state to a disordered splayed-open structure and then to a disordered left-handed conformation-as the dimer probes are moved across the ss-dsDNA junction. Our results suggest that the sugar-phosphate backbones at and near ss-dsDNA junctions adopt specific position-dependent local conformations and exhibit varying extents of conformational disorder that deviate widely from the Watson-Crick structure. We suggest that some of these conformations can function as secondary-structure motifs for interaction with protein complexes that bind to and assemble at these sites.

A computational insight into the relationship between side chain IR line shapes and local environment in fibril-like structures

Infrared spectroscopy is a widely used technique to characterize protein structures and protein mediated processes. While the amide I band provides information on proteins' secondary structure, amino acid side chains are used as infrared probes for the investigation of protein reactions and local properties. In this paper, we use a hybrid quantum mechanical/classical molecular dynamical approach based on the perturbed matrix method to compute the infrared band due to the C=O stretching mode of amide-containing side chains. We calculate, at first, the infrared band of zwitterionic glutamine in water and obtain results in very good agreement with the experimental data. Then, we compute the signal arising from glutamine side chains in a microcrystal of the yeast prion Sup35-derived peptide, GNNQQNY, with a fibrillar structure. The infrared bands obtained by selective isotopic labeling of the two glutamine residues, Q4 and Q5, of each peptide were experimentally used to investigate the local hydration in the fibrillar microcrystal. The experimental spectra of the two glutamine residues, which experience different hydration environments, feature different spectral signals that are well reproduced by the corresponding calculated spectra. In addition, the analysis of the simulated spectra clarifies the molecular origin of the experimentally observed spectroscopic differences that arise from the different local electric field experienced by the two glutamine residues, which is, in turn, determined by a different hydrogen bonding pattern.

Does liquid-liquid phase separation drive peptide folding?

Proline–arginine (PR) dipeptide repeats have been shown to undergo liquid–liquid phase separation and are an example of a growing number of intrinsically disordered proteins that can assemble into membraneless organelles. These structures have been posited as nucleation sites for pathogenic protein aggregation. As such, a better understanding of the effects that the increased local concentration and volumetric crowding within droplets have on peptide secondary structure is necessary. Herein we use Fourier transform infrared (FTIR) and two-dimensional infrared (2DIR) spectroscopy to show that formation of droplets by PR₂₀ accompanies changes in the amide-I spectra consistent with folding into poly-proline helical structures.

Two-dimensional infrared spectroscopy reveals folding of an intrinsically disordered peptide when sequestered into a model “membrane-less” organelle.

Bottom-Up Approach to Assess the Molecular Structure of Aqueous Poly(N-Isopropylacrylamide) at Room Temperature via Infrared Spectroscopy

The structure of poly(N-isopropylacrylamide) (PNIPAM) in solution is still an unresolved topic. Here, the PNIPAM structure in water was investigated using a bottom-up approach, involving the monomer, dimer, and trimer, and a combination of infrared (IR) spectroscopies as well as molecular dynamics simulations. The experiments show that the monomer and oligomers exhibit a broad and asymmetric amide I band with two underlying transitions, while PNIPAM presents the same major transitions and a minor one. Analysis of the 2D IR spectra and theoretical modeling of the amide I band indicates that the two transitions of the monomer do not have the same molecular origin as the oligomers and the polymer. In the monomer, the two bands originate from the ultrafast rotation of its ethyl group, which leads to different solvation structures for the various rotational conformers. In the case of the oligomers, the asymmetry and splitting of the amide I band is caused by the vibrational coupling among adjacent amide side chains. Moreover, it is deduced from the simulations that the oligomers have three distinct backbone conformations for neighboring amides. In particular, two of the backbone conformations have a closed and compact structure, while in the third, the backbone is open and elongated. The bottom-up approach allowed us to infer that such backbone conformations exist in PNIPAM as well. Consequently, the two major amide I transitions of the polymer are also assigned to split amide I transitions resulting from the vibrationally coupled nearest-neighboring amides. In contrast, the additional minor transition observed in PNIPAM is assigned to unsolvated amide units of the polymer. The proposed molecular model successfully describes that PNIPAM amide I band changes with temperature in terms of its molecular structure. This new model strongly suggests that PNIPAM does not have a completely random backbone structure, but has distinct backbone conformers between neighboring amides.

Studies on the structure-activity relationship and interaction mechanism of flavonoids and xanthine oxidase through enzyme kinetics, spectroscopy methods and molecular simulations

In this study, some flavonoids were screened as potent xanthine oxidase (XO) inhibitors in vitro. Flavonoid 9 was demonstrated to exhibit the inhibitory activity through a ping-pong mechanism. Further structure-activity relationship revealed that different structural elements had greatly influenced the inhibition effect on XO and underlined the requirement of hydroxyl groups at C5 and C4' of flavonoid type I. Moreover, some bioactive flavonoids could efficiently quench the intrinsic fluorescence of XO by either static or static-dynamic mixed mechanism. The synchronous fluorescence, ANS-binding fluorescence, Fourier transform infrared spectra and circular dichroism suggested that active flavonoids could bind to the active center of XO, prevent the entrance of substrate, and induce the rearrangement and conformation change of its secondary structures, ultimately resulting in the significant inhibition effect. Additionally, molecular docking further confirmed these conclusions and highlighted the great importance of hydrophobic interactions and hydrogen bonds for the formation of stable complex conformation.

Modeling amino-acid side chain infrared spectra: the case of carboxylic residues

Infrared (IR) spectroscopy is commonly utilized for the investigation of protein structures and protein-mediated processes.

Integrated multi-spectroscopic and molecular modelling techniques to probe the interaction mechanism between salvianolic acid A and α‑glucosidase

α-Glucosidase (AG) is an important drug target for the treatment of type 2 diabetes mellitus in humans due to the potential effect of down regulating glucose absorption in patients. In our previous study, salvianolic acid A (SAA) was found to exhibit potent AG inhibitory activity, whereas the interaction mechanism was still ambiguous. Herein, the interaction mechanism of SAA and AG was investigated by multi-spectroscopic methods along with molecular docking. As a result, it was found that SAA reversibly inhibited AG in a competitive manner with IC₅₀ of 16.44 ± 0.18 μM, and the inhibition belonged to a multi-phase kinetics process with a first-order reaction. The intrinsic fluorescence of AG could be strongly quenched by SAA through a static quenching mechanism. The negative Gibbs free energy change and positive values of enthalpy and entropy change revealed that the binding of SAA to AG was spontaneous and dominated mainly by hydrophobic interactions, and only a single binding site was determined for them. Analysis of synchronous fluorescence, ANS-binding fluorescence, circular dichroism and Fourier transform infrared spectra suggested that the binding of SAA to AG induced rearrangement and conformational changes of the enzyme. Besides, further molecular modelling validated that SAA could bind to the active domain and prevent the entrance of substrate, resulting in the inhibition of AG activity. These findings provide new insights into understanding the interaction mechanism of SAA on AG.

Effect of Salvia miltiorrhiza on acetylcholinesterase: Enzyme kinetics and interaction mechanism merging with molecular docking analysis

Acetylcholinesterase (AchE) serves as an important target for Alzheimer's disease. Salvia miltiorrhiza has been used to treat cardiovascular disease for hundreds of years. However, the interaction between S. miltiorrhiza and AchE is still inadequate. Herein, an integrated method including molecular docking and experimental studies was employed to investigate the interaction. Consequently, some components were screened as potent AchE inhibitors by in silico and in vitro. Among them, miltirone (MT) and salvianolic acid A (SAA) reversibly inhibited AchE in a mixed-competitive manner. Fluorescence data revealed that SAA and salvianolic acid C (SAC) strongly quenched the intrinsic fluorescence of AchE through a static quenching mechanism, and the binding was spontaneous and dominated by hydrophobic interaction inferred by the thermodynamic parameters. The synchronous and ANS-binding fluorescence spectra suggested that SAA and SAC could bind to the enzyme and induce its conformation changes of secondary structures, which was further confirmed by Fourier transform infrared spectra. Meanwhile, molecular docking presented the probable binding modes of inhibitors to AchE and highlighted the key role of hydrophobic interaction and hydrogen bonds for the stability of docking complex. These findings put more insights into understanding the interaction of S. miltiorrhiza chemicals and AchE, as well as Alzheimer's disease.

Investigation of the interaction between salvianolic acid C and xanthine oxidase: Insights from experimental studies merging with molecular docking methods

Xanthine oxidase (XO) has emerged as an important target for gout. In our previous study, salvianolic acid C (SAC) was found to show potent XO inhibitory activity, whereas the interaction mechanism was still not clear. Herein, an integrated approach consisting of enzyme kinetics, multi-spectroscopic methods and molecular docking was employed to investigate the interaction between SAC and XO. Consequently, SAC exhibited a rapid and mixed-type inhibition of XO with IC₅₀ of 5.84 ± 0.18 μM. The fluorescence data confirmed that SAC presented a strong fluorescence quenching effect through a static quenching procedure. The values of enthalpy change, entropy change and Gibbs free energy change indicated that their binding was spontaneous and driven mainly by hydrophobic interactions. Analysis of synchronous fluorescence, circular dichroism and fourier transform infrared spectra demonstrated that SAC induced conformational changes of the enzyme. Besides, further molecular docking revealed that SAC occupied the catalytic center resulting in the inhibition of XO activity. This study provides a comprehensive understanding on the interaction mechanism of SAC on XO.

Insight into the internal structure of amyloid-β oligomers by isotope-edited Fourier transform infrared spectroscopy

Isotope-edited infrared spectroscopy reveals the structural unit of amyloid-β oligomers.

Physiological Calcium Concentrations Slow Dynamics at the Lipid-Water Interface

Phospholipids can interact strongly with ions at physiological concentrations, and these interactions can alter membrane properties. Here, we describe the effects of calcium ions on the dynamics in phospholipid membranes. We used a combination of time-resolved ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations. We found that millimolar Ca²⁺ concentrations lead to slower fluctuations in the local environment at the lipid-water interface of membranes with phosphatidylserine. The effect was only observed in bilayers containing anionic phosphatidylserine; membranes composed of only zwitterionic phosphatidylcholine did not experience a slowdown. Local water dynamics were measured using the ester groups as label-free probes and were found to be up to 50% slower with 2.5 mM Ca²⁺. Molecular dynamics simulations show that Ca²⁺ primarily binds to the carboxylate group of phosphatidylserines. These findings have implications for apoptotic and diseased cells in which phosphatidylserine is exposed to extracellular calcium and for the biophysical effects of divalent cations on lipid bilayers.

Characterization of the interaction between human serum albumin and diazinon via spectroscopic and molecular docking methods

Human serum albumin (HSA) is a soluble blood protein which binds to small molecules (such as drugs and toxins) and transfers them within the blood circulation. In this research, the interaction of diazinon, as a toxic organophosphate, with HSA was investigated. Various biophysical methods such as fluorescence, ultraviolet–visible (UV-vis), Fourier transform infrared spectroscopy, and molecular docking were utilized to characterize the binding properties of diazinon to HSA under physiological-like condition. The UV-vis spectroscopy showed that the absorption increased and the fluorescence intensity of HSA decreased regularly with regard to the gradual increases of the concentrations of diazinon. Due to the binding constant of ( ka = 3.367 × 10+4 M−1), the α-helix structure for the first day and 35 days of incubation were obtained 66.09–55.4% and 59.99–46.48%, respectively, and their amounts in other secondary structures (β-sheet, β-anti, and random (r) coils) were increased. The molecular docking revealed a good binding site in HSA (Trp-214) for diazinon which was related to the considerable alterations in HSA secondary and tertiary structures. There is a close relationship between the secondary structure of protein and its biological activity and after 35 days of incubation, the high toxic concentrations of diazinon can make HSA to be partially unfolded and lose its structure.

A Free Energy Barrier Caused by the Refolding of an Oligomeric Intermediate Controls the Lag Time of Amyloid Formation by hIAPP

Transiently populated oligomers formed en route to amyloid fibrils may constitute the most toxic aggregates associated with many amyloid-associated diseases. Most nucleation theories used to describe amyloid aggregation predict low oligomer concentrations and do not take into account free energy costs that may be associated with structural rearrangements between the oligomer and fiber states. We have used isotope labeling and two-dimensional infrared spectroscopy to spectrally resolve an oligomeric intermediate during the aggregation of the human islet amyloid protein (hIAPP or amylin), the protein associated with type II diabetes. A structural rearrangement includes the F₂₃G₂₄A₂₅I₂₆L₂₇ region of hIAPP, which starts from a random coil structure, evolves into ordered β-sheet oligomers containing at least 5 strands, and then partially disorders in the fibril structure. The supercritical concentration is measured to be between 150 and 250 μM, which is the thermodynamic parameter that sets the free energy of the oligomers. A 3-state kinetic model fits the experimental data, but only if it includes a concentration independent free energy barrier >3 kcal/mol that represents the free energy cost of refolding the oligomeric intermediate into the structure of the amyloid fibril; i.e., "oligomer activation" is required. The barrier creates a transition state in the free energy landscape that slows fibril formation and creates a stable population of oligomers during the lag phase, even at concentrations below the supercritical concentration. Largely missing in current kinetic models is a link between structure and kinetics. Our experiments and modeling provide evidence that protein structural rearrangements during aggregation impact the populations and kinetics of toxic oligomeric species.

Computing infrared spectra of proteins using the exciton model

The ability to compute from first principles the infrared spectrum of a protein in solution phase representing a biological system would provide a useful connection to atomistic models of protein structure and dynamics. Indeed, such calculations are a vital complement to 2DIR experimental measurements, allowing the observed signals to be interpreted in terms of detailed structural and dynamical information. In this article, we have studied nine structurally and spectroscopically well-characterized proteins, representing a range of structural types. We have simulated the equilibrium conformational dynamics in an explicit point charge water model. Using the resulting trajectories based on MD simulations, we have computed the one and two dimensional infrared spectra in the Amide I region, using an exciton approach, in which a local mode basis of carbonyl stretches is considered. The role of solvent in shifting the Amide I band (by 30 to 50 cm-1 ) is clearly evident. Similarly, the conformational dynamics contribute to the broadening of peaks in the spectrum. The inhomogeneous broadening in both the 1D and 2D spectra reflects the significant conformational diversity observed in the simulations. Through the computed 2D cross-peak spectra, we show how different pulse schemes can provide additional information on the coupled vibrations. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

Assessing Spectral Simulation Protocols for the Amide I Band of Proteins

We present a benchmark study of spectral simulation protocols for the amide I band of proteins. The amide I band is widely used in infrared spectroscopy of proteins due to the large signal intensity, high sensitivity to hydrogen bonding, and secondary structural motifs. This band has, thus, proven valuable in many studies of protein structure-function relationships. We benchmark spectral simulation protocols using two common force fields in combination with several electrostatic mappings and coupling models. The results are validated against experimental linear absorption and two-dimensional infrared spectroscopy for three well-studied proteins. We find two-dimensional infrared spectroscopy to be much more sensitive to the simulation protocol than linear absorption and report on the best simulation protocols. The findings demonstrate that there is still room for ideas to improve the existing models for the amide I band of proteins.

Application of two-dimensional infrared spectroscopy to benchmark models for the amide I band of proteins

In this paper, we present a novel benchmarking method for validating the modelling of vibrational spectra for the amide I region of proteins. We use the linear absorption spectra and two-dimensional infrared spectra of four experimentally well-studied proteins as a reference and test nine combinations of molecular dynamics force fields, vibrational frequency mappings, and coupling models. We find that two-dimensional infrared spectra provide a much stronger test of the models than linear absorption does. The best modelling approach in the present study still leaves significant room for future improvement. The presented benchmarking scheme, thus, provides a way of validating future protocols for modelling the amide I band in proteins.

Isotope-enriched protein standards for computational amide I spectroscopy

We present a systematic isotope labeling study of the protein G mutant NuG2b as a step toward the production of reliable, structurally stable, experimental standards for amide I infrared spectroscopic simulations. By introducing isotope enriched amino acids into a minimal growth medium during bacterial expression, we induce uniform labeling of the amide bonds following specific amino acids, avoiding the need for chemical peptide synthesis. We use experimental data to test several common amide I frequency maps and explore the influence of various factors on map performance. Comparison of the predicted absorption frequencies for the four maps tested with empirical assignments to our experimental spectra yields a root-mean-square error of 6-12 cm(-1), with outliers of at least 12 cm(-1) in all models. This means that the predictions may be useful for predicting general trends such as changes in hydrogen bonding configuration; however, for finer structural constraints or absolute frequency assignments, the models are unreliable. The results indicate the need for careful testing of existing literature maps and shed light on possible next steps for the development of quantitative spectral maps.

Structural disorder of folded proteins: isotope-edited 2D IR spectroscopy and Markov state modeling

The conformational heterogeneity of the N-terminal domain of the ribosomal protein L9 (NTL91-39) in its folded state is investigated using isotope-edited two-dimensional infrared spectroscopy. Backbone carbonyls are isotope-labeled ((13)C=(18)O) at five selected positions (V3, V9, V9G13, G16, and G24) to provide a set of localized spectroscopic probes of the structure and solvent exposure at these positions. Structural interpretation of the amide I line shapes is enabled by spectral simulations carried out on structures extracted from a recent Markov state model. The V3 label spectrum indicates that the β-sheet contacts between strands I and II are well folded with minimal disorder. The V9 and V9G13 label spectra, which directly probe the hydrogen-bond contacts across the β-turn, show significant disorder, indicating that molecular dynamics simulations tend to overstabilize ideally folded β-turn structures in NTL91-39. In addition, G24-label spectra provide evidence for a partially disordered α-helix backbone that participates in hydrogen bonding with the surrounding water.

Amide I Vibrational Properties Affected by Hydrogen Bonding Out-of-Plane of the Peptide Group

The amide I vibrational properties of a peptide-water complex in various intermolecular configurations are analyzed theoretically to see whether a water molecule with a weak out-of-plane hydrogen bond really induces a large low-frequency shift. It is shown that the frequency shift strongly depends on the C═O···H angle, with a larger low-frequency shift as the C═O···H becomes more bent, suggesting that the so-called hydrated helix with a rather low amide I frequency has an additional water molecule located out-of-plane of the peptide group as compared with a typical α-helix. The infrared intensity also depends on the angular position of water. A new model parameter set (that can be combined with molecular dynamics) is developed for a more correct representation of the hydration-induced frequency shift. The question regarding the scalar and vectorial nature of the molecular properties related to the frequency shift is also discussed.

A molecular interpretation of 2D IR protein folding experiments with Markov state models

The folding mechanism of the N-terminal domain of ribosomal protein L9 (NTL91-39) is studied using temperature-jump (T-jump) amide I' two-dimensional infrared (2D IR) spectroscopy in combination with spectral simulations based on a Markov state model (MSM) built from millisecond-long molecular dynamics trajectories. The results provide evidence for a compact well-structured folded state and a heterogeneous fast-exchanging denatured state ensemble exhibiting residual secondary structure. The folding rate of 26.4 μs(-1) (at 80°C), extracted from the T-jump response of NTL91-39, compares favorably with the 18 μs(-1) obtained from the MSM. Structural decomposition of the MSM and analysis along the folding coordinate indicates that helix-formation nucleates the global folding. Simulated difference spectra, corresponding to the global folding transition of the MSM, are in qualitative agreement with measured T-jump 2D IR spectra. The experiments demonstrate the use of T-jump 2D IR spectroscopy as a valuable tool for studying protein folding, with direct connections to simulations. The results suggest that in addition to predicting the correct native structure and folding time constant, molecular dynamics simulations carried out with modern force fields provide an accurate description of folding mechanisms in small proteins.

Electrostatic frequency shifts in amide I vibrational spectra: direct parameterization against experiment

The interpretation of protein amide I infrared spectra has been greatly assisted by the observation that the vibrational frequency of a peptide unit reports on its local electrostatic environment. However, the interpretation of spectra remains largely qualitative due to a lack of direct quantitative connections between computational models and experimental data. Here, we present an empirical parameterization of an electrostatic amide I frequency map derived from the infrared absorption spectra of 28 dipeptides. The observed frequency shifts are analyzed in terms of the local electrostatic potential, field, and field gradient, evaluated at sites near the amide bond in molecular dynamics simulations. We find that the frequency shifts observed in experiment correlate very well with the electric field in the direction of the C=O bond evaluated at the position of the amide oxygen atom. A linear best-fit mapping between observed frequencies and electric field yield sample standard deviations of 2.8 and 3.7 cm(-1) for the CHARMM27 and OPLS-AA force fields, respectively, and maximum deviations (within our data set) of 9 cm(-1). These results are discussed in the broader context of amide I vibrational models and the effort to produce quantitative agreement between simulated and experimental absorption spectra.

Microsecond folding experiments and simulations: a match is made

For the past two decades, protein folding experiments have been speeding up from the second or millisecond time scale to the microsecond time scale, and full-atom simulations have been extended from the nanosecond to the microsecond and even millisecond time scale. Where the two meet, it is now possible to compare results directly, allowing force fields to be validated and refined, and allowing experimental data to be interpreted in atomistic detail. In this perspective we compare recent experiments and simulations on the microsecond time scale, pointing out the progress that has been made in determining native structures from physics-based simulations, refining experiments and simulations to provide more quantitative underlying mechanisms, and tackling the problems of multiple reaction coordinates, downhill folding, and complex underlying structure of unfolded or misfolded states.

Solvent and conformation dependence of amide I vibrations in peptides and proteins containing proline

We present a mixed quantum-classical model for studying the amide I vibrational dynamics (predominantly CO stretching) in peptides and proteins containing proline. There are existing models developed for determining frequencies of and couplings between the secondary amide units. However, these are not applicable to proline because this amino acid has a tertiary amide unit. Therefore, a new parametrization is required for infrared-spectroscopic studies of proteins that contain proline, such as collagen, the most abundant protein in humans and animals. Here, we construct the electrostatic and dihedral maps accounting for solvent and conformation effects on frequency and coupling for the proline unit. We examine the quality and the applicability of these maps by carrying out spectral simulations of a number of peptides with proline in D(2)O and compare with experimental observations.

Vibrational spectroscopic determination of local solvent electric field, solute-solvent electrostatic interaction energy, and their fluctuation amplitudes

IR probes have been extensively used to monitor local electrostatic and solvation dynamics. Particularly, their vibrational frequencies are highly sensitive to local solvent electric field around an IR probe. Here, we show that the experimentally measured vibrational frequency shifts can be inversely used to determine local electric potential distribution and solute-solvent electrostatic interaction energy. In addition, the upper limits of their fluctuation amplitudes are estimated by using the vibrational bandwidths. Applying this method to fully deuterated N-methylacetamide (NMA) in D(2)O and examining the solvatochromic effects on the amide I' and II' mode frequencies, we found that the solvent electric potential difference between O(═C) and D(-N) atoms of the peptide bond is about 5.4 V, and thus, the approximate solvent electric field produced by surrounding water molecules on the NMA is 172 MV/cm on average if the molecular geometry is taken into account. The solute-solvent electrostatic interaction energy is estimated to be -137 kJ/mol, by considering electric dipole-electric field interaction. Furthermore, their root-mean-square fluctuation amplitudes are as large as 1.6 V, 52 MV/cm, and 41 kJ/mol, respectively. We found that the water electric potential on a peptide bond is spatially nonhomogeneous and that the fluctuation in the electrostatic peptide-water interaction energy is about 10 times larger than the thermal energy at room temperature. This indicates that the peptide-solvent interactions are indeed important for the activation of chemical reactions in aqueous solution.

Mechanism and conformational studies of farrerol binding to bovine serum albumin by spectroscopic methods

The mechanism and conformational changes of farrerol binding to bovine serum albumin (BSA) were studied by spectroscopic methods including fluorescence quenching technique, UV-vis absorption, circular dichroism (CD) spectroscopy and Fourier transform infrared (FT-IR) spectroscopy under simulative physiological conditions. The results of fluorescence titration revealed that farrerol could strongly quench the intrinsic fluorescence of BSA through a static quenching procedure. The thermodynamic parameters enthalpy change and entropy change for the binding were calculated to be -29.92 kJ mol(-1) and 5.06 J mol(-1) K(-1) according to the van't Hoff equation, which suggested that the both hydrophobic interactions and hydrogen bonds play major role in the binding of farrerol to BSA. The binding distance r deduced from the efficiency of energy transfer was 3.11 nm for farrerol-BSA system. The displacement experiments of site markers and the results of fluorescence anisotropy showed that warfarin and farrerol shared a common binding site I corresponding to the subdomain IIA of BSA. Furthermore, the studies of synchronous fluorescence, CD and FT-IR spectroscopy showed that the binding of farrerol to BSA induced conformational changes in BSA.

Spectroscopic analysis of protein Fe–NO complexes

The toxic free radical NO (nitric oxide) has diverse biological roles in eukaryotes and bacteria, being involved in signalling, vasodilation, blood clotting and immunity, and as an intermediate in microbial denitrification. The predominant biological mechanism of detecting NO is through the formation of iron nitrosyl complexes, although this is a deleterious process for other iron-containing enzymes. We have previously applied techniques such as UV–visible and EPR spectroscopy to the analysis of protein Fe–NO complex formation in order to study how NO controls the activity of the bacterial transcriptional regulators NorR and NsrR. These studies have analysed NO-dependent biological activity both in vitro and in vivo using diverse biochemical, molecular and spectroscopic methods. Recently, we have applied ultrafast 2D-IR (two-dimensional IR) spectroscopy to the analysis of NO–protein interactions using Mb (myoglobin) and Cc (cytochrome c) as model haem proteins. The ultrafast fluctuations of Cc and Mb show marked differences, indicating altered flexibility of the haem pockets. We have extended this analysis to bacterial catalase enzymes that are known to play a role in the nitrosative stress response by detoxifying peroxynitrite. The first 2D-IR analysis of haem nitrosylation and perspectives for the future are discussed.

Development and validation of transferable amide I vibrational frequency maps for peptides

Infrared (IR) spectroscopy of the amide I band has been widely utilized for the analysis of peptides and proteins. Theoretical modeling of IR spectra of proteins requires an accurate and efficient description of the amide I frequencies. In this paper, amide I frequency maps for protein backbone and side chain groups are developed from experimental spectra and vibrational lifetimes of N-methylacetamide and acetamide in different solvents. The frequency maps, along with established nearest-neighbor frequency shift and coupling schemes, are then applied to a variety of peptides in aqueous solution and reproduce experimental spectra well. The frequency maps are designed to be transferable to different environments; therefore, they can be used for heterogeneous systems, such as membrane proteins.

Molecular theory and simulation of coherence transfer in metal carbonyls and its signature on multidimensional infrared spectra

We present a general and comprehensive theoretical and computational framework for modeling ultrafast multidimensional infrared spectra of a vibrational excitonic system in liquid solution. Within this framework, we describe the dynamics of the system in terms of a quantum master equation that can account for population relaxation, dephasing, coherence-to-coherence transfer, and coherence-to-population transfer. A unique feature of our approach is that, in principle, it does not rely on any adjustable fitting parameters. More specifically, the anharmonic vibrational Hamiltonian is derived from ab initio electronic structure theory, and the system-bath coupling is expressed explicitly in terms of liquid degrees of freedom whose dynamics can be obtained via molecular dynamics simulations. The applicability of the new approach is demonstrated by employing it to model the recently observed signatures of coherence transfer in the two-dimensional spectra of dimanganese decacarbonyl in liquid cyclohexane. The results agree well with experiment and shed new light on the nature of the molecular interactions and dynamics underlying the spectra and the interplay between dark and bright states, their level of degeneracy, and the nature of their interactions with the solvent.

Modeling the vibrational dynamics and nonlinear infrared spectra of coupled amide I and II modes in peptides

The amide vibrational modes play an important role in energy transport and relaxation in polypeptides and proteins and provide us with spectral markers for structure and structural dynamics of these macromolecules. Here, we present a detailed model to describe the dynamic properties of the amide I and amide II modes and the resulting linear and nonlinear spectra. These two modes have large oscillator strengths, and their mutual coupling plays an important role in their relaxation. Using first-principles calculations of NMA-d(7) and a dipeptide in a fluctuating bath described by molecular dynamics simulations, we model the frequencies of the local vibrations as well as the coupling between them. Both the coherent couplings and the fluctuations induced by contact with their environment are taken into account. We apply the resulting model of interacting fluctuating oscillators to study the collective vibrations and the partially coherent transport of vibrational energy through a model α-helix. We find that the instantaneous vibrations are delocalized over a few (up to four) amide units, while the coherences in the helix survive for 0.5-1 ps, leading to coherent transport on a similar time scale.

Melting of a beta-hairpin peptide using isotope-edited 2D IR spectroscopy and simulations

Isotope-edited two-dimensional infrared spectroscopy has been used to characterize the conformational heterogeneity of the beta-hairpin peptide TrpZip2 (TZ2) across its thermal unfolding transition. Four isotopologues were synthesized to probe hydrogen bonding and solvent exposure of the beta-turn (K8), the N-terminus (S1), and the midstrand region (T10 and T3T10). Isotope-shifts, 2D lineshapes, and other spectral changes to the amide I 2D IR spectra of labeled TZ2 isotopologues were observed as a function of temperature. Data were interpreted on the basis of structure-based spectroscopic modeling of conformers obtained from extensive molecular dynamics simulations. The K8 spectra reveal two unique turn geometries, the type I' beta-turn observed in the NMR structure, and a less populated disordered or bulged loop. The data indicate that structures at low temperature resemble the folded NMR structure with typical cross-strand hydrogen bonds, although with a subpopulation of misformed turns. As the temperature is raised from 25 to 85 degrees C, the fraction of population with a type I' turn increases, but the termini also fray. Hydrogen bonding contacts in the midstrand region remain at all temperatures although with increasing thermal disorder. Our data show no evidence of an extended chain or random coil state for the TZ2 peptide at any temperature. The methods demonstrated here offer an approach to characterizing conformational variation within the folded or unfolded states of proteins and peptides.