Refinement of Peptide Conformational Ensembles by 2D IR Spectroscopy: Application to Ala‒Ala‒Ala

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─Quo Vadis?

Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.

Progress in infrared spectroscopy as an efficient tool for predicting protein secondary structure

Infrared (IR) spectroscopy is a highly sensitive technique that provides complete information on chemical compositions. The IR spectra of proteins or peptides give rise to nine characteristic IR absorption bands. The amide I bands are the most prominent and sensitive vibrational bands and widely used to predict protein secondary structures. The interference of H₂O absorbance is the greatest challenge for IR protein secondary structure prediction. Much effort has been made to reduce/eliminate the interference of H₂O, simplify operation steps, and increase prediction accuracy. Progress in sampling and equipment has rendered the Fourier transform infrared (FTIR) technique suitable for determining the protein secondary structure in broader concentration ranges, greatly simplifying the operating steps. This review highlights the recent progress in sample preparation, data analysis, and equipment development of FTIR in A/T mode, with a focus on recent applications of FTIR spectroscopy in the prediction of protein secondary structure. This review also provides a brief introduction of the progress in ATR-FTIR for predicting protein secondary structure and discusses some combined IR methods, such as AFM-based IR spectroscopy, that are used to analyze protein structural dynamics and protein aggregation.

Conformation and Metal Cation Binding of Zwitterionic Alanine Tripeptide in Saline Solutions by Infrared Vibrational Spectroscopy and Molecular Dynamics Simulations

In this work, linear infrared (IR) spectroscopy and molecular dynamics (MD) simulations were used to examine the interaction of different metal cations (Na⁺, Ca²⁺, Mg²⁺, and Zn²⁺) with backbone (amide C═O) and C-terminal carboxylate (COO^-) groups in zwitterionic alanine tripeptide (Ala3) in aqueous solutions with varying saline concentrations. Circular dichroism spectra and MD results suggest that Ala3 is predominantly in polyproline-II (PP_II) conformation, whose amide-I and asymmetric carboxylate stretching IR vibration signatures are also supported by quantum-chemistry calculations. The zwitterionic form of Ala3 separates the two amide-I modes in frequency, which are weakly coupled modes, as revealed by two-dimensional IR measurement, and can be used to probe backbone-cation interactions at different scenarios (near charged or neutral chemical groups respectively). Cation concentration-dependent IR frequency red shifts in the amide-I mode are seen for both amide-I modes, whereas blue shifts are also seen in the amide-I mode far from the NH₃⁺ group. The observed spectral changes are discussed from the perspective of the salting-in and salting-out abilities of the cations. In addition, all the metal cations studied here (except Zn²⁺) can specifically coordinate to the COO^- group in bidentate and pseudo-bridging forms simultaneously. For Zn²⁺, only the pseudo-bridging form exists. Our results shed light on the macroscopic protein salting-in and salting-out phenomena from the perspective of key chemical bonds in peptides.

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

The dynamics and spectroscopy of N-methyl-acetamide (NMA) and trialanine in solution are characterized from molecular dynamics simulations using different energy functions, including a conventional point charge (PC)-based force field, one based on a multipolar (MTP) representation of the electrostatics, and a semiempirical DFT method. For the 1D infrared spectra, the frequency splitting between the two amide-I groups is 10 cm^-1 from the PC, 13 cm^-1 from the MTP, and 47 cm^-1 from self-consistent charge density functional tight-binding (SCC-DFTB) simulations, compared with 25 cm^-1 from experiment. The frequency trajectory required for the frequency fluctuation correlation function (FFCF) is determined from individual normal mode (INM) and full normal mode (FNM) analyses of the amide-I vibrations. The spectroscopy, time-zero magnitude of the FFCF C(t = 0), and the static component Δ₀² from simulations using MTP and analysis based on FNM are all consistent with experiments for (Ala)₃. Contrary to this, for the analysis excluding mode-mode coupling (INM), the FFCF decays to zero too rapidly and for simulations with a PC-based force field, the Δ₀² is too small by a factor of two compared with experiments. Simulations with SCC-DFTB agree better with experiment for these observables than those from PC-based simulations. The conformational ensemble sampled from simulations using PCs is consistent with the literature (including P_II, β, α_R, and α_L), whereas that covered by the MTP-based simulations is dominated by P_II with some contributions from β and α_R. This agrees with and confirms recently reported Bayesian-refined populations based on 1D infrared experiments. FNM analysis together with a MTP representation provides a meaningful model to correctly describe the dynamics of hydrated trialanine.

Pump Slice Amplitudes: A Simple and Robust Method for Connecting Two-Dimensional Infrared and Fourier Transform Infrared Spectra

Ultrafast two-dimensional infrared (2D IR) spectroscopy and Fourier transform infrared (FTIR) spectroscopy are often performed in tandem, with FTIR typically used to interpret and provide hypotheses for 2D IR experiments. Comparisons between 2D IR and FTIR spectra can also be used to examine the structure and orientation in systems of coupled vibrational chromophores. The most common method for comparing 2D IR and FTIR lineshapes, the diagonal slice method, contains significant artifacts when applied to oscillators with low anharmonicities. Here, we introduce a new technique, the pump slice amplitude (PSA) method, for relating 2D IR lineshapes to FTIR lineshapes and compare PSAs against diagonal slices using theoretical and experimental spectra. We find that PSAs are significantly more similar to FTIR lineshapes than diagonal slices in systems with low anharmonicity.

Substituent Effects on the Vibrational Properties of the CN Stretch Mode of Aromatic Nitriles: IR Probes Useful for Time-resolved IR Spectroscopy

Developing ideal IR probes is essential to understand the structure and dynamics of biomolecules with time-resolved IR spectroscopies and imaging techniques. Especially, nitrile (CN) group has recently been proposed to serve as IR probes of the local environment of proteins. Herein, we investigated the effect of a substituent on the vibrational properties of the benzonitrile. The electron-donating and withdrawing character of p-substituent on benzonitrile are expected to modulate the vibrational frequency, molar extinction coefficient, and vibrational lifetime of CN probe. FT-IR revealed the positive correlation between electron-donating character and the molar extinction coefficient of CN stretch mode. Infrared pump-probe (IR-PP) measurements showed that the vibrational lifetime of CN stretch mode exhibits a relatively weak correlation with the electron-donating strength. Among the investigated samples, 4-dimethylamino benzonitrile with the strongest electron-donating strength shows enhanced absorption and extended vibrational lifetime. Utilizing substituent effects will be a practical strategy to improve the performance of the IR probe.

Computational IR Spectroscopy of Insulin Dimer Structure and Conformational Heterogeneity

We have investigated the structure and conformational dynamics of insulin dimer using a Markov state model (MSM) built from extensive unbiased atomistic molecular dynamics simulations and performed infrared spectral simulations of the insulin MSM to describe how structural variation within the dimer can be experimentally resolved. Our model reveals two significant conformations to the dimer: a dominant native state consistent with other experimental structures of the dimer and a twisted state with a structure that appears to reflect a ∼55° clockwise rotation of the native dimer interface. The twisted state primarily influences the contacts involving the C-terminus of insulin's B chain, shifting the registry of its intermolecular hydrogen bonds and reorganizing its side-chain packing. The MSM kinetics predict that these configurations exchange on a 14 μs time scale, largely passing through two Markov states with a solvated dimer interface. Computational amide I spectroscopy of site-specifically ¹³C¹⁸O labeled amides indicates that the native and twisted conformation can be distinguished through a series of single and dual labels involving the B24F, B25F, and B26Y residues. Additional structural heterogeneity and disorder is observed within the native and twisted states, and amide I spectroscopy can also be used to gain insight into this variation. This study will provide important interpretive tools for IR spectroscopic investigations of insulin structure and transient IR kinetics experiments studying the conformational dynamics of insulin dimer.

Effect of the Hydration Shell on the Carbonyl Vibration in the Ala-Leu-Ala-Leu Peptide

The vibrational spectrum of the Ala-Leu-Ala-Leu peptide in solution, computed from first-principles simulations, shows a prominent band in the amide I region that is assigned to stretching of carbonyl groups. Close inspection reveals combined but slightly different contributions by the three carbonyl groups of the peptide. The shift in their exact vibrational signature is in agreement with the different probabilities of these groups to form hydrogen bonds with the solvent. The central carbonyl group has a hydrogen bond probability intermediate to the other two groups due to interchanges between different hydrogen-bonded states. Analysis of the interaction energies of individual water molecules with that group shows that shifts in its frequency are directly related to the interactions with the water molecules in the first hydration shell. The interaction strength is well correlated with the hydrogen bond distance and hydrogen bond angle, though there is no perfect match, allowing geometrical criteria for hydrogen bonds to be used as long as the sampling is sufficient to consider averages. The hydrogen bond state of a carbonyl group can therefore serve as an indicator of the solvent's effect on the vibrational frequency.

Specific and non-specific interactions between metal cations and zwitterionic alanine tripeptide in saline solutions reported by the symmetric carboxylate stretching and amide-II vibrations

The "specific" interaction between metal cations (Na⁺, Ca²⁺, Mg²⁺, and Zn²⁺) and the charged COO^- group, and the "non-specific" interaction between these cations and the peptide backbone of a zwitterionic trialanine (Ala3) in aqueous solutions were examined in detail, using linear infrared (IR) absorptions of the COO^- symmetric stretching and the amide-II (mainly the C-N stretching) modes as IR probes. Different IR spectral changes in peak positions and intensities of the two IR probes clearly demonstrate their sensitivities to nearby cation distributions in distance and population. Quantum chemistry calculations and molecular dynamics simulations were used to describe the cation-peptide interaction picture. These combined results suggest that Na⁺ and Ca²⁺ tend to bind to the COO^- group in the bidentate form, while Mg²⁺ and Zn²⁺ tend to bind to the COO^- group in the pseudo-bridging form. The results also show that while all three divalent cations indirectly interact with the peptide backbone with large population, Ca²⁺ and Mg²⁺ can be sometimes distributed very close to the backbone. Such a non-specific cation interaction can be moderately sensed by the C-N stretching of the amide-II mode when cations approach the polar amide C[double bond, length as m-dash]O group, and is also influenced by the NH₃⁺ charge group located at the N-terminus. The results suggest that the experimentally observed complication of the Hofmeister cation series shall be understood as a combined specific and non-specific cation-peptide interactions.

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.

Insulin Dissociates by Diverse Mechanisms of Coupled Unfolding and Unbinding

The protein hormone insulin exists in various oligomeric forms, and a key step in binding its cellular receptor is dissociation of the dimer. This dissociation process and its corresponding association process have come to serve as paradigms of coupled (un)folding and (un)binding more generally. Despite its fundamental and practical importance, the mechanism of insulin dimer dissociation remains poorly understood. Here, we use molecular dynamics simulations, leveraging recent developments in umbrella sampling, to characterize the energetic and structural features of dissociation in unprecedented detail. We find that the dissociation is inherently multipathway with limiting behaviors corresponding to conformational selection and induced fit, the two prototypical mechanisms of coupled folding and binding. Along one limiting path, the dissociation leads to detachment of the C-terminal segment of the insulin B chain from the protein core, a feature believed to be essential for receptor binding. We simulate IR spectroscopy experiments to aid in interpreting current experiments and identify sites where isotopic labeling can be most effective for distinguishing the contributions of the limiting mechanisms.

Efficient Ensemble Refinement by Reweighting

Ensemble refinement produces structural ensembles of flexible and dynamic biomolecules by integrating experimental data and molecular simulations. Here we present two efficient numerical methods to solve the computationally challenging maximum-entropy problem arising from a Bayesian formulation of ensemble refinement. Recasting the resulting constrained weight optimization problem into an unconstrained form enables the use of gradient-based algorithms. In two complementary formulations that differ in their dimensionality, we optimize either the log-weights directly or the generalized forces appearing in the explicit analytical form of the solution. We first demonstrate the robustness, accuracy, and efficiency of the two methods using synthetic data. We then use NMR J-couplings to reweight an all-atom molecular dynamics simulation ensemble of the disordered peptide Ala-5 simulated with the AMBER99SB*-ildn-q force field. After reweighting, we find a consistent increase in the population of the polyproline-II conformations and a decrease of α-helical-like conformations. Ensemble refinement makes it possible to infer detailed structural models for biomolecules exhibiting significant dynamics, such as intrinsically disordered proteins, by combining input from experiment and simulation in a balanced manner.