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

Mapping the amide-I vibrations of model dipeptides with secondary structure sensitivity and amino acid residue specificity, and its application to amyloid β-peptide in aqueous solution

Vibrational spectroscopy has been known as particularly well-suited for deciphering the polypeptide's structure. To decode structural information encoded in IR spectra, we developed amide-I frequency maps on the basis of model dipeptides to correlate the amide-I frequency of interest to the combination of the calculated secondary structure dependent amide-I frequency by using DFT method and the electrostatic potentials that projected onto the amide unit from the micro-environment within molecular mechanics force field. The constructed maps were applied to model dipeptides and amyloid β-peptide fragment (Aβ_25-35). The dipeptide specified map (DS map) and the hybrid map (HYB map) predicted amide-I bands of Aβ_25-35 in solution satisfactorily reproduce experimental observation, and indicate the preference of forming β-sheet and random coil structure for Aβ_25-35 in D₂O just as the results of cluster analysis suggested. These maps with secondary structural sensitivity and amino acid residue specificity open up a way for the interpretation of amide-I vibrations and show their potentials in the understanding of molecular structure of polypeptides in solution.

Probing the Differential Dynamics of the Monomeric and Dimeric Insulin from Amide-I IR Spectroscopy

The monomer-dimer equilibrium for insulin is one of the essential steps in forming the receptor-binding competent monomeric form of the hormone. Despite this importance, the thermodynamic stability, in particular for modified insulins, is quite poorly understood, in part, due to experimental difficulties. This work explores one- and two-dimensional infrared spectroscopy in the range of the amide-I band for the hydrated monomeric and dimeric wild-type hormone. It is found that for the monomer the frequency fluctuation correlation function (FFCF) and the one-dimensional infrared spectra are position sensitive. The spectra of the -CO probes at the dimerization interface (residues Phe24, Phe25, and Tyr26) split and indicate an asymmetry despite the overall (formal) point symmetry of the dimer structure. Also, the decay times of the FFCF for the same -CO probe in the monomer and the dimer can differ by up to 1 order of magnitude, for example, for residue PheB24, which is solvent exposed for the monomer but at the interface for the dimer. The spectroscopic shifts correlate approximately with the average number of hydration waters and the magnitude of the FFCF at time zero. However, this correlation is only qualitative due to the heterogeneous and highly dynamical environment. Based on density functional theory calculations, the dominant contribution for solvent-exposed -CO is found to arise from the surrounding water (∼75%), whereas the protein environment contributes considerably less. The results suggest that infrared spectroscopy is a positionally sensitive probe of insulin dimerization, in particular in conjunction with isotopic labeling of the probe.

Development of Vibrational Frequency Maps for Nucleobases

Vibrational spectroscopy provides a powerful tool to probe the structure and dynamics of nucleic acids because specific normal modes, particularly the base carbonyl stretch modes, are highly sensitive to the hydrogen bonding patterns and stacking configurations in these biomolecules. In this work, we develop vibrational frequency maps for the C═O and C═C stretches in nucleobases that allow the calculations of their site frequencies directly from molecular dynamics simulations. We assess the frequency maps by applying them to nucleobase derivatives in aqueous solutions and nucleosides in organic solvents and demonstrate that the predicted infrared spectra are in good agreement with experimental measurements. The frequency maps can be readily used to model the linear and nonlinear vibrational spectroscopy of nucleic acids and elucidate the molecular origin of the experimentally observed spectral features.

A neural network protocol for electronic excitations of N-methylacetamide

UV absorption is widely used for characterizing proteins structures. The mapping of UV spectra to atomic structure of proteins relies on expensive theoretical simulations, circumventing the heavy computational cost which involves repeated quantum-mechanical simulations of excited-state properties of many fluctuating protein geometries, which has been a long-time challenge. Here we show that a neural network machine-learning technique can predict electronic absorption spectra of N-methylacetamide (NMA), which is a widely used model system for the peptide bond. Using ground-state geometric parameters and charge information as descriptors, we employed a neural network to predict transition energies, ground-state, and transition dipole moments of many molecular-dynamics conformations at different temperatures, in agreement with time-dependent density-functional theory calculations. The neural network simulations are nearly 3,000× faster than comparable quantum calculations. Machine learning should provide a cost-effective tool for simulating optical properties of proteins.

Glutamine Side Chain 13C═18O as a Nonperturbative IR Probe of Amyloid Fibril Hydration and Assembly

Infrared (IR) spectroscopy has provided considerable insight into the structures, dynamics, and formation mechanisms of amyloid fibrils. IR probes, such as main chain ¹³C═¹⁸O, have been widely employed to obtain site-specific structural information, yet only secondary structures and strand-to-strand arrangements can be probed. Very few nonperturbative IR probes are available to report on the side-chain conformation and environments, which are critical to determining sheet-to-sheet arrangements in steric zippers within amyloids. Polar residues, such as glutamine, contribute significantly to the stability of amyloids and thus are frequently found in core regions of amyloid peptides/proteins. Furthermore, polyglutamine (polyQ) repeats form toxic aggregates in several neurodegenerative diseases. Here we report the synthesis and application of a new nonperturbative IR probe-glutamine side chain ¹³C═¹⁸O. We use side chain ¹³C═¹⁸O labeling and isotope dilution to detect the presence of intermolecularly hydrogen-bonded arrays of glutamine side chains (Gln ladders) in amyloid-forming peptides. Moreover, the line width of the ¹³C═¹⁸O peak is highly sensitive to its local hydration environment. The IR data from side chain labeling allows us to unambiguously determine the sheet-to-sheet arrangement in a short amyloid-forming peptide, GNNQQNY, providing insight that was otherwise inaccessible through main chain labeling. With several different fibril samples, we also show the versatility of this IR probe in studying the structures and aggregation kinetics of amyloids. Finally, we demonstrate the capability of modeling amyloid structures with IR data using the integrative modeling platform (IMP) and the potential of integrating IR with other biophysical methods for more accurate structural modeling. Together, we believe that side chain ¹³C═¹⁸O will complement main chain isotope labeling in future IR studies of amyloids and integrative modeling using IR data will significantly expand the power of IR spectroscopy to elucidate amyloid assemblies.

Engineering Proteins at Interfaces: From Complementary Characterization to Material Surfaces with Designed Functions

Once materials come into contact with a biological fluid containing proteins, proteins are generally—whether desired or not—attracted by the material's surface and adsorb onto it. The aim of this Review is to give an overview of the most commonly used characterization methods employed to gain a better understanding of the adsorption processes on either planar or curved surfaces. We continue to illustrate the benefit of combining different methods to different surface geometries of the material. The thus obtained insight ideally paves the way for engineering functional materials that interact with proteins in a predetermined manner.

Molecular simulations of self-assembling bio-inspired supramolecular systems and their connection to experiments

In bionanotechnology, the field of creating functional materials consisting of bio-inspired molecules, the function and shape of a nanostructure only appear through the assembly of many small molecules together. The large number of building blocks required to define a nanostructure combined with the many degrees of freedom in packing small molecules has long precluded molecular simulations, but recent advances in computational hardware as well as software have made classical simulations available to this strongly expanding field. Here, we review the state of the art in simulations of self-assembling bio-inspired supramolecular systems. We will first discuss progress in force fields, simulation protocols and enhanced sampling techniques using recent examples. Secondly, we will focus on efforts to enable the comparison of experimentally accessible observables and computational results. Experimental quantities that can be measured by microscopy, spectroscopy and scattering can be linked to simulation output either directly or indirectly, via quantum mechanical or semi-empirical techniques. Overall, we aim to provide an overview of the various computational approaches to understand not only the molecular architecture of nanostructures, but also the mechanism of their formation.

A Direct, Quantitative Connection between Molecular Dynamics Simulations and Vibrational Probe Line Shapes

A quantitative connection between molecular dynamics simulations and vibrational spectroscopy of probe-labeled systems would enable direct translation of experimental data into structural and dynamical information. To constitute this connection, all-atom molecular dynamics (MD) simulations were performed for two SCN probe sites (solvent-exposed and buried) in a calmodulin-target peptide complex. Two frequency calculation approaches with substantial nonelectrostatic components, a quantum mechanics/molecular mechanics (QM/MM)-based technique and a solvatochromic fragment potential (SolEFP) approach, were used to simulate the infrared probe line shapes. While QM/MM results disagreed with experiment, SolEFP results matched experimental frequencies and line shapes and revealed the physical and dynamic bases for the observed spectroscopic behavior. The main determinant of the CN probe frequency is the exchange repulsion between the probe and its local structural neighbors, and there is a clear dynamic explanation for the relatively broad probe line shape observed at the "buried" probe site. This methodology should be widely applicable to vibrational probes in many environments.

Direct Anionic Effect on Water Structure and Indirect Anionic Effect on Peptide Backbone Hydration State Revealed by Thin-Layer Infrared Spectroscopy

In this work, the anionic effect on water structure and on the peptide backbone and water interaction was investigated directly in aqueous solution using thin-layer transmission infrared spectroscopy. The chaotropic anions were found to weaken the water hydrogen-bonding strength and red shift the HOH bending frequency, while the kosmotropic anions were found to strengthen the water hydrogen-bonding network and blue shift the HOH bending frequency. The kosmotropes, especially F^-, blue shift the vibrational frequencies of both amide II and amide III bands of N-methylacetamide (NMA), indicating NMA is in the "salting-in" state; while the chaotropes (Cl^-, NO₃^-, Br^-, I^-, and SCN^-) red shift the frequencies of the two normal modes, indicating NMA is in the "salting-out" state. Furthermore, the changes of the vibrational frequencies of the HOH bending, amide II and III bands were found to generally follow the Hofmeister anionic series. Our results suggest that hydrated anion influences the peptide backbone mainly through the N-H group, but a weak and indirect effect through the amide C═O group also contributes. Thus, these amide modes can be used as vibrational measures of anionic influences on peptide backbone's hydration state. Our work also suggests that deuteration of the amide unit decreases the sensitivity of the amide II and III vibrational modes in this regard.

Strategy for Modeling the Electrostatic Responses of the Spectroscopic Properties of Proteins

For better understanding and more efficient use of the spectroscopic probes (vibrational and NMR) of the local electrostatic situations inside proteins, appropriate modeling of the properties of those probes is essential. The present study is devoted to examining the strategy for constructing such models. A more well-founded derivation than the ones in previous studies is given in constructing the models. Theoretical analyses are conducted on two representative example cases related to proteins, i.e., the peptide group of the main chains and the CO and NO ligands to the Fe²⁺ ion of heme, with careful treatment of the behavior of electrons in the electrostatic responses and with verification of consistency with observable quantities. It is shown that, for the stretching frequencies and NMR chemical shifts, it is possible to construct reasonable electrostatic interaction models that encompass the situations of hydration and uniform electric field environment and thus are applicable also to the cases of nonuniform electrostatic situations, which are highly expected for inside of proteins.

Transient two-dimensional vibrational spectroscopy of an operating molecular machine

Synthetic molecular machines are promising building blocks for future nanoscopic devices. However, the details of their mechanical behaviour are in many cases still largely unknown. A deeper understanding of mechanics at the molecular level is essential for the design and construction of complex nanodevices. Here, we show that transient two-dimensional infrared (T2DIR) spectroscopy makes it possible to monitor the conformational changes of a translational molecular machine during its operation. Translation of a macrocyclic ring from one station to another on a molecular thread is initiated by a UV pulse. The arrival of the shuttling macrocycle at the final station is visible from a newly appearing cross peak between these two moieties. To eliminate spectral congestion in the T2DIR spectra, we use a subtraction method applicable to many other complex molecular systems. The T2DIR spectra indicate that the macrocycle adopts a boat-like conformation at the final station, which contrasts with the chair-like conformation at the initial station.

A deeper understanding of the mechanics of molecular machines is limited by the fast motions which are in the nanosecond or picosecond timescale. Here the authors present a real-time observation of structural changes in a rotaxane-based molecular shuttle by transient two-dimensional infrared spectroscopy.

Revealing the Solvation Structure and Dynamics of Carbonate Electrolytes in Lithium-Ion Batteries by Two-Dimensional Infrared Spectrum Modeling

Carbonate electrolytes in lithium-ion batteries play a crucial role in conducting lithium ions between two electrodes. Mixed solvent electrolytes consisting of linear and cyclic carbonates are commonly used in commercial lithium-ion batteries. To understand how the linear and cyclic carbonates introduce different solvation structures and dynamics, we performed molecular dynamics simulations of two representative electrolyte systems containing either linear or cyclic carbonate solvents. We then modeled their two-dimensional infrared (2DIR) spectra of the carbonyl stretching mode of these carbonate molecules. We found that the chemical exchange process involving formation and dissociation of lithium-ion/carbonate complexes is responsible for the growth of 2DIR cross peaks with increasing waiting time. In addition, we also found that cyclic carbonates introduce faster dynamics of dissociation and formation of lithium-ion/carbonate complexes than linear carbonates. These findings provide new insights into understanding the lithium-ion mobility and its interplay with solvation structure and ultrafast dynamics in carbonate electrolytes used in lithium-ion batteries.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

The dynamics of peptide-water interactions in dialanine: An ultrafast amide I 2D IR and computational spectroscopy study

We present a joint experimental and computational study of the dynamic interactions of dialanine (Ala-Ala) with water, comparing the results of ultrafast 2D IR and infrared transient absorption spectroscopy of its amide I vibration with spectra modeled from molecular dynamics (MD) simulations. The experimental data are analyzed to describe vibrational frequency fluctuations, vibrational energy relaxation, and chemical exchange processes. The origin of these processes in the same underlying fluctuating forces allows a common description in terms of the fluctuations and conformational dynamics of the peptide and associated solvent. By comparing computational spectroscopy from MD simulations with multiple force fields and water models, we describe how the dynamics of water hydrogen bond fluctuations and switching processes act as a source of friction that governs the dephasing and vibrational relaxation, and provide a description of coupled water and peptide motions that give rise to spectroscopic exchange processes.

Watching Proteins Wiggle: Mapping Structures with Two-Dimensional Infrared Spectroscopy

Proteins exhibit structural fluctuations over decades of time scales. From the picosecond side chain motions to aggregates that form over the course of minutes, characterizing protein structure over these vast lengths of time is important to understanding their function. In the past 15 years, two-dimensional infrared spectroscopy (2D IR) has been established as a versatile tool that can uniquely probe proteins structures on many time scales. In this review, we present some of the basic principles behind 2D IR and show how they have, and can, impact the field of protein biophysics. We highlight experiments in which 2D IR spectroscopy has provided structural and dynamical data that would be difficult to obtain with more standard structural biology techniques. We also highlight technological developments in 2D IR that continue to expand the scope of scientific problems that can be accessed in the biomedical sciences.

Structural and Spectroscopic Properties of Assemblies of Self-Replicating Peptide Macrocycles

Self-replication at the molecular level is often seen as essential to the early origins of life. Recently a mechanism of self-replication has been discovered in which replicator self-assembly drives the process. We have studied one of the examples of such self-assembling self-replicating molecules to a high level of structural detail using a combination of computational and spectroscopic techniques. Molecular Dynamics simulations of self-assembled stacks of peptide-derived replicators provide insights into the structural characteristics of the system and serve as the basis for semiempirical calculations of the UV-vis, circular dichroism (CD) and infrared (IR) absorption spectra that reflect the chiral organization and peptide secondary structure of the stacks. Two proposed structural models are tested by comparing calculated spectra to experimental data from electron microscopy, CD and IR spectroscopy, resulting in a better insight into the specific supramolecular interactions that lead to self-replication. Specifically, we find a cooperative self-assembly process in which β-sheet formation leads to well-organized structures, while also the aromatic core of the macrocycles plays an important role in the stability of the resulting fibers.

Electrostatic frequency maps for amide-I mode of β-peptide: Comparison of molecular mechanics force field and DFT calculations

The spectroscopy of amide-I vibrations has been widely utilized for the understanding of dynamical structure of polypeptides. For the modeling of amide-I spectra, two frequency maps were built for β-peptide analogue (N-ethylpropionamide, NEPA) in a number of solvents within different schemes (molecular mechanics force field based, GM map; DFT calculation based, GD map), respectively. The electrostatic potentials on the amide unit that originated from solvents and peptide backbone were correlated to the amide-I frequency shift from gas phase to solution phase during map parameterization. GM map is easier to construct with negligible computational cost since the frequency calculations for the samples are purely based on force field, while GD map utilizes sophisticated DFT calculations on the representative solute-solvent clusters and brings insight into the electronic structures of solvated NEPA and its chemical environments. The results show that the maps' predicted amide-I frequencies present solvation environmental sensitivities and exhibit their specific characters with respect to the map protocols, and the obtained vibrational parameters are in satisfactory agreement with experimental amide-I spectra of NEPA in solution phase. Although different theoretical schemes based maps have their advantages and disadvantages, the present maps show their potentials in interpreting the amide-I spectra for β-peptides, respectively.

Amyloid Beta Peptide Folding in Reverse Micelles

Previously published experimental studies have suggested that when the 40-residue amyloid beta peptide is encapsulated in a reverse micelle, it folds into a structure that may nucleate amyloid fibril formation (Yeung, P. S.-W.; Axelsen, P. H. J. Am. Chem. Soc. 2012, 134, 6061 ). The factors that induce the formation of this structure have now been identified in a multi-microsecond simulation of the same reverse micelle system that was studied experimentally. Key features of the polypeptide-micelle interaction include the anchoring of a hydrophobic residue cluster into gaps in the reverse micelle surface, the formation of a beta turn at the anchor point that brings N- and C-terminal segments of the polypeptide into proximity, high ionic strength that promotes intramolecular hydrogen bond formation, and deformation of the reverse micelle surface to facilitate interactions with the surface along the entire length of the polypeptide. Together, these features cause the simulation-derived vibrational spectrum to red shift in a manner that reproduces the red-shift previously reported experimentally. On the basis of these findings, a new mechanism is proposed whereby membranes nucleate fibril formation and facilitate the in-register alignment of polypeptide strands that is characteristic of amyloid fibrils.

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.

Interplay between Hydrogen Bonding and Vibrational Coupling in Liquid N-Methylacetamide

Intrinsically disordered proteins play an important role in biology, and unraveling their labile structure presents a vital challenge. However, the dynamical structure of such proteins thwarts their study by standard techniques such as X-ray diffraction and NMR spectroscopy. Here, we use a neat liquid composed of N-methylacetamide molecules as a model system to elucidate dynamical and structural properties similar to those one can expect to see in intrinsically disordered proteins. To examine the structural dynamics in the neat liquid, we combine molecular dynamics, response-function-based spectral simulations, and two-dimensional polarization-resolved infrared spectroscopy in the amide I (CO stretch) region. The two-dimensional spectra reveal a delicate interplay between hydrogen bonding and intermolecular vibrational coupling effects, observed through a fast anisotropy decay. The present study constitutes a general platform for understanding the structure and dynamics of highly disordered proteins.

Probing the Effects of Gating on the Ion Occupancy of the K+ Channel Selectivity Filter Using Two-Dimensional Infrared Spectroscopy

The interplay between the intracellular gate and the selectivity filter underlies the structural basis for gating in potassium ion channels. Using a combination of protein semisynthesis, two-dimensional infrared (2D IR) spectroscopy, and molecular dynamics (MD) simulations, we probe the ion occupancy at the S1 binding site in the constricted state of the selectivity filter of the KcsA channel when the intracellular gate is open and closed. The 2D IR spectra resolve two features, whose relative intensities depend on the state of the intracellular gate. By matching the experiment to calculated 2D IR spectra of structures predicted by MD simulations, we identify the two features as corresponding to states with S1 occupied or unoccupied by K⁺. We learn that S1 is >70% occupied when the intracellular gate is closed and <15% occupied when the gate is open. Comparison of MD trajectories show that opening of the intracellular gate causes a structural change in the selectivity filter, which leads to a change in the ion occupancy. This work reveals the complexity of the conformational landscape of the K⁺ channel selectivity filter and its dependence on the state of the intracellular gate.

Perspective: Quantum mechanical methods in biochemistry and biophysics

In this perspective article, I discuss several research topics relevant to quantum mechanical (QM) methods in biophysical and biochemical applications. Due to the immense complexity of biological problems, the key is to develop methods that are able to strike the proper balance of computational efficiency and accuracy for the problem of interest. Therefore, in addition to the development of novel ab initio and density functional theory based QM methods for the study of reactive events that involve complex motifs such as transition metal clusters in metalloenzymes, it is equally important to develop inexpensive QM methods and advanced classical or quantal force fields to describe different physicochemical properties of biomolecules and their behaviors in complex environments. Maintaining a solid connection of these more approximate methods with rigorous QM methods is essential to their transferability and robustness. Comparison to diverse experimental observables helps validate computational models and mechanistic hypotheses as well as driving further development of computational methodologies.

Solvent-Independent Anharmonicity for Carbonyl Oscillators

The physical origins of vibrational frequency shifts have been extensively studied in order to understand noncovalent intermolecular interactions in the condensed phase. In the case of carbonyls, vibrational solvatochromism, MD simulations, and vibrational Stark spectroscopy suggest that the frequency shifts observed in simple solvents arise predominately from the environment's electric field due to the vibrational Stark effect. This is contrary to many previously invoked descriptions of vibrational frequency shifts, such as bond polarization, whereby the bond's force constant and/or partial nuclear charges are altered due to the environment, often illustrated in terms of favored resonance structures. Here we test these hypotheses using vibrational solvatochromism as measured using 2D IR to assess the solvent dependence of the bond anharmonicity. These results indicate that the carbonyl bond's anharmonicity is independent of solvent as tested using hexanes, DMSO, and D₂O and is supported by simulated 2D spectra. In support of the linear vibrational Stark effect, these 2D IR measurements are consistent with the assertion that the Stark tuning rate is unperturbed by the electric field generated by both hydrogen and non-hydrogen bonding environments and further extends the general applicability of carbonyl probes for studying intermolecular interactions.

Refining Disordered Peptide Ensembles with Computational Amide I Spectroscopy: Application to Elastin-Like Peptides

The characterization of intrinsically disordered protein (IDP) ensembles is complicated both by inherent heterogeneity and by the fact that many common experimental techniques function poorly when applied to IDPs. For this reason, the development of alternative structural tools for probing IDP ensembles has attracted considerable attention. Here we describe our recent work in developing experimental and computational tools for characterizing IDP ensembles using Amide I (backbone carbonyl stretch) vibrational spectroscopy. In this approach, the infrared (IR) absorption frequencies of isotope-labeled amide bonds probe their local electrostatic environments and structures. Empirical frequency maps allow us to use this spectroscopic data as a direct experimental test of atomistic structural models. We apply these methods to a family of short elastin-like peptides (ELPs), fragments of the elastin protein based around the Pro-Gly turn motif characteristic of the elastomeric segments of the full protein. Using a maximum entropy analysis that applies constraints from experimental spectra to weighting predicted spectra from molecular dynamics (MD) ensembles, we find that peptides with Ala or Val side chains preceding the Pro-Gly turn unit exhibit a stronger tendency toward extended structures than do Gly-Pro-Gly motifs, suggesting an important role for steric interactions in tuning the molecular properties of elastin.

Vibrational solvatochromism. III. Rigorous treatment of the dispersion interaction contribution

A rigorous first principles theory of vibrational solvatochromism including the intermolecular dispersion interaction, which is based on the effective fragment potential method, is developed. The present theory is an extended version of our previous vibrational solvatochromism model that took into account the Coulomb, exchange-repulsion, and induction interactions. We show that the frequency shifts of the amide I mode of N-methylacetamide in H2O and CDCl3, when combined with molecular dynamics simulations, can be quantitatively reproduced by the theory, which indicates that the dispersion interaction contribution to the vibrational frequency shift is not always negligibly small. Nonetheless, the reason that the purely Coulombic interaction model for vibrational solvatochromism works well for describing amide I mode frequency shifts in polar solvents is because the electrostatic contribution is strong and highly sensitive to the relative orientation of surrounding solvent molecules, which is in stark contrast with polarization, dispersion, and exchange-repulsion contributions. It is believed that the theory presented and discussed here will be of great use in quantitatively describing vibrational solvatochromism and electrochromism of infrared probes in not just polar solvent environments but also in biopolymers such as proteins.

Instantaneous ion configurations in the K+ ion channel selectivity filter revealed by 2D IR spectroscopy

Potassium channels are responsible for the selective permeation of K+ ions across cell membranes. K+ ions permeate in single file through the selectivity filter, a narrow pore lined by backbone carbonyls that compose four K+ binding sites. Here, we report on the two-dimensional infrared (2D IR) spectra of a semisynthetic KcsA channel with site-specific heavy (13C18O) isotope labels in the selectivity filter. The ultrafast time resolution of 2D IR spectroscopy provides an instantaneous snapshot of the multi-ion configurations and structural distributions that occur spontaneously in the filter. Two elongated features are resolved, revealing the statistical weighting of two structural conformations. The spectra are reproduced by molecular dynamics simulations of structures with water separating two K+ ions in the binding sites, ruling out configurations with ions occupying adjacent sites.

Vibrational Stark Effects of Carbonyl Probes Applied to Reinterpret IR and Raman Data for Enzyme Inhibitors in Terms of Electric Fields at the Active Site

IR and Raman frequency shifts have been reported for numerous probes of enzyme transition states, leading to diverse interpretations. In the case of the model enzyme ketosteroid isomerase (KSI), we have argued that IR spectral shifts for a carbonyl probe at the active site can provide a connection between the active site electric field and the activation free energy (Fried et al. Science 2014, 346, 1510-1514). Here we generalize this approach to a much broader set of carbonyl probes (e.g., oxoesters, thioesters, and amides), first establishing the sensitivity of each probe to an electric field using vibrational Stark spectroscopy, vibrational solvatochromism, and MD simulations, and then applying these results to reinterpret data already in the literature for enzymes such as 4-chlorobenzoyl-CoA dehalogenase and serine proteases. These results demonstrate that the vibrational Stark effect provides a general framework for estimating the electrostatic contribution to the catalytic rate and may provide a metric for the design or modification of enzymes. Opportunities and limitations of the approach are also described.

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.

Communication: Quantitative multi-site frequency maps for amide I vibrational spectroscopy

An accurate method for predicting the amide I vibrational spectrum of a given protein structure has been sought for many years. Significant progress has been made recently by sampling structures from molecular dynamics simulations and mapping local electrostatic variables onto the frequencies of individual amide bonds. Agreement with experiment, however, has remained largely qualitative. Previously, we used dipeptide fragments and isotope-labeled constructs of the protein G mimic NuG2b as experimental standards for developing and testing amide I frequency maps. Here, we combine these datasets to test different frequency-map models and develop a novel method to produce an optimized four-site potential (4P) map based on the CHARMM27 force field. Together with a charge correction for glycine residues, the optimized map accurately describes both experimental datasets, with average frequency errors of 2-3 cm(-1). This 4P map is shown to be convertible to a three-site field map which provides equivalent performance, highlighting the viability of both field- and potential-based maps for amide I spectral modeling. The use of multiple sampling points for local electrostatics is found to be essential for accurate map performance.

Roles of the scalar and vector components of the solvation effects on the vibrational properties of hydrogen- or halogen-bond accepting stretching modes

Solvation-induced vibrational frequency shifts and infrared (IR) intensity changes of the hydrogen- or halogen-bond accepting stretching modes, especially their dependence on the angular position of the hydrogen- or halogen-bond donating molecule, are examined theoretically. Calculations are carried out for some modes of hydrogen- or halogen-bonding molecular complexes, including the S[double bond, length as m-dash]O stretch of dimethyl sulfoxide-(13)C2H2O, the C[triple bond, length as m-dash]N stretch of acetonitrileH2O, and the amide I' mode of the N-methylacetamide-d1BrNC 1 : 1 complex. It is shown that, in all the example cases dealt with in this study, the frequency shift depends rather strongly on the hydrogen- or halogen-bond angle (e.g., S[double bond, length as m-dash]OH angle), with a larger low-frequency shift as the hydrogen or halogen bond becomes more bent, indicating the generality of the results obtained for the amide I' mode of the N-methylacetamide-d1(2)H2O 1 : 1 complex in a previous study. Contrary to our vague expectation, the frequency shift is not well correlated to the hydrogen- or halogen-bond distance or strength, but nevertheless, it is well reproduced by an electrostatic interaction model if it is carefully constructed by considering the scalar and vector components separately in a reasonable way. On the basis of this electrostatic interaction model, the reason why our vague expectation is not realized is clarified, and a unified understanding is achieved on the hydration-induced high-frequency shift of the C[triple bond, length as m-dash]N stretch and the low-frequency shifts of the S[double bond, length as m-dash]O stretch and amide I'. With regard to the IR intensity, it is shown that, in some of the example cases, it also has rather strong angular position dependence. The mechanism of the IR intensity changes is estimated by analyzing the dipole derivative vector, especially its angular relation with the hydrogen or halogen bond.

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.

General Applicable Frequency Map for the Amide-I Mode in β-Peptides

In this work, a general applicable amide-I vibrational frequency map (GA map) for β-peptides in a number of common solvents was constructed, based on a peptide derivative, N-ethylpropionamide (NEPA). The map utilizes force fields at the ab initio computational level to accurately describe molecular structure and solute-solvent interactions, and also force fields at the molecular mechanics level to take into account long-range solute-solvent interactions. The results indicate that the GA map works reasonably for mapping the vibrational frequencies of the amide-I local-modes for β-peptides, holding promises for understanding the complicated infrared spectra of the amide-I mode in β-polypeptides.

Fragment quantum chemical approach to geometry optimization and vibrational spectrum calculation of proteins

Geometry optimization and vibrational spectra (infrared and Raman spectra) calculations of proteins are carried out by a quantum chemical approach using the EE-GMFCC (electrostatically embedded generalized molecular fractionation with conjugate caps) method (J. Phys. Chem. A, 2013, 117, 7149).