Spectroscopic studies of protein folding: linear and nonlinear methods

Raman Spectra of Amino Acids and Peptides from Machine Learning Polarizabilities

Raman spectroscopy is an important tool in the study of vibrational properties and composition of molecules, peptides, and even proteins. Raman spectra can be simulated based on the change of the electronic polarizability with vibrations, which can nowadays be efficiently obtained via machine learning models trained on first-principles data. However, the transferability of the models trained on small molecules to larger structures is unclear, and direct training on large structures is prohibitively expensive. In this work, we first train two machine learning models to predict the polarizabilities of all 20 amino acids. Both models are carefully benchmarked and compared to density functional theory (DFT) calculations, with the neural network method being found to offer better transferability. By combination of machine learning models with classical force field molecular dynamics, Raman spectra of all amino acids are also obtained and investigated, showing good agreement with experiments. The models are further extended to small peptides. We find that adding structures containing peptide bonds to the training set greatly improves predictions, even for peptides not included in training sets.

A TeZla micromixer for interrogating the early and broad folding landscape of G-quadruplex via multistage velocity descending

Biomacromolecular folding kinetics involves fast folding events and broad timescales. Current techniques face limitations in either the required time resolution or the observation window. In this study, we developed the TeZla micromixer, integrating Tesla and Zigzag microstructures with a multistage velocity descending strategy. TeZla achieves a significant short mixing dead time (40 µs) and a wide time window covering four orders of magnitude (up to 300 ms). Using this unique micromixer, we explored the folding landscape of c-Myc G4 and its noncanonical-G4 derivatives with different loop lengths or G-vacancy sites. Our findings revealed that c-Myc can bypass folding intermediates and directly adopt a G4 structure in the cation-deficient buffer. Moreover, we found that the loop length and specific G-vacancy site could affect the folding pathway and significantly slow down the folding rates. These results were also cross-validated with real-time NMR and circular dichroism. In conclusion, TeZla represents a versatile tool for studying biomolecular folding kinetics, and our findings may ultimately contribute to the design of drugs targeting G4 structures.

Inconsistent hydrogen bond-mediated vibrational coupling of amide I

Using infrared spectroscopy and density functional theory (DFT) calculations, we scrutinized an amide (dimethylformamide) as a "model" compound to interpret the interactions of amide 1 with different phenol derivatives (para-chlorophenol (PCP) and para-cresol (CP)) as "model guest molecules". We established the involvement of amide I in vibrational coupling with symmetric and asymmetric C[double bond, length as m-dash]C modes of different phenolic derivatives and how their coupling was dependent upon different guest aromatic phenolic compounds. Interestingly, substitution of phenol perturbed the pattern of vibrational coupling with amide I. The symmetric and asymmetric C[double bond, length as m-dash]C modes of PC were coupled significantly with amide 1. For PCP, the symmetric C[double bond, length as m-dash]C mode coupled significantly, but the asymmetric mode coupled negligibly, with amide I. Here, we reveal the nature of vibrational coupling based on the structure of a guest molecule hydrogen-bonded with amide I. Our conclusions could be valuable for depiction of the unusual dynamics of coupled amide-I modes as well as the dependency of vibrational coupling on altered factors.

An "Onion-like" Model of Protein Unfolding: Collective versus Site Specific Approaches

Approximating protein unfolding by an all-or-none cooperative event is a convenient assumption that can provide precious global information on protein stability. It is however quickly emerging that the scenario is far more complex and that global denaturation curves often hide a rich heterogeneity of states that are largely probe dependent. In this review, we revisit the importance of gaining site-specific information on the unfolding process. We focus on nuclear magnetic resonance, as this is the main technique able to provide site-specific information. We review historical and most modern approaches that have allowed an appreciable advancement of the field of protein folding. We also demonstrate how unfolding is a reporter dependent event, suggesting the outmost importance of selecting the reporter carefully.

Potential of Raman spectroscopic techniques to study proteins

Proteins are large, complex molecules responsible for various biological processes. However, protein misfolding may lead to various life-threatening diseases. Therefore, it is vital to understand the shape and structure of proteins. Despite numerous techniques, a mechanistic understanding of the protein folding process is still unclear. Therefore, new techniques are continually being explored. In the present article, we have discussed the importance of Raman spectroscopy, Raman Optical Activity (ROA) and various other advancements in Raman spectroscopy to understand protein structure and conformational changes based on the review of our earlier work and recent literature. A Raman spectrum of a protein provides unique signatures for various secondary structures like helices, beta-sheets, turns, random structures, etc., and various amino acid residues such as tyrosine, tryptophan, and phenylalanine. We have shown how Raman spectra can differentiate between bovine serum albumin (BSA) and lysozyme protein based on their difference in sequence and structure (primary, secondary and tertiary). Although it is challenging to elucidate the structure of a protein using a Raman spectrum alone, Raman spectra can be used to differentiate small changes in conformations of proteins such as BSA during melting. Various new advancements in technique and data analyses in Raman spectroscopic studies of proteins have been discussed. The last part of the review focuses on the importance of the ROA spectrum to understand additional features about proteins. The ROA spectrum is rich in information about the protein backbone due to its rigidity compared to its side chains. Furthermore, the ROA spectra of lysozyme and BSA have been presented to show how ROA provides extra information about the solvent properties of proteins.

Tuning the electronic transition energy of indole via substitution: application to identify tryptophan-based chromophores that absorb and emit visible light

Fluorescent amino acids (FAAs) offer significant advantages over fluorescent proteins in applications where the fluorophore size needs to be limited or minimized. A long-sought goal in biological spectroscopy/microcopy is to develop visible FAAs by modifying the indole ring of tryptophan. Herein, we examine the absorption spectra of a library of 4-substituted indoles and find that the frequency of the absorption maximum correlates linearly with the global electrophilicity index of the substituent. This finding permits us to identify two promising candidates, 4-formyltryptophan (4CHO-Trp) and 4-nitrotryptophan (4NO₂-Trp), both of which can be excited by visible light. Further fluorescence measurements indicate that while 4CHO-indole (and 4CHO-Trp) emits cyan fluorescence with a reasonably large quantum yield (ca. 0.22 in ethanol), 4NO₂-indole is essentially non-fluorescent, suggesting that 4CHO-Trp (4NO₂-Trp) could be useful as a fluorescence reporter (quencher). In addition, we present a simple method for synthesizing 4CHO-Trp.

Infrared spectroscopy probes ion binding geometries

Infrared (IR) spectroscopy is a well-established technique for probing the structure, behavior, and surroundings of molecules in their native environments. Its characteristics-most specifically high structural sensitivity, ready applicability to aqueous samples, and broad availability-make it a valuable enzymological technique, particularly for the interrogation of ion binding sites. While IR spectroscopy of the "garden variety" (steady state at room temperature with wild-type proteins) is versatile and powerful in its own right, the combination of IR spectroscopy with specialized experimental schemes for leveraging ultrafast time resolution, protein labeling, and other enhancements further extends this utility. This book chapter provides the fundamental physical background and literature context essential for harnessing IR spectroscopy in the general context of enzymology with specific focus on interrogation of ion binding. Studies of lanthanide ions binding to calmodulin are highlighted as illustrative examples of this process. Appropriate sample preparation, data collection, and spectral interpretation are discussed from a detail-oriented and practical perspective with the goal of facilitating the reader's rapid progression from reading words in a book to collecting and analyzing their own data in the lab.

Two-dimensional IR spectroscopy reveals a hidden Fermi resonance band in the azido stretch spectrum of β-azidoalanine

Azido stretch modes in a variety of azido-derivatized nonnatural amino acids and nucleotides have been used as a site-specific infrared (IR) probe for monitoring changes in their conformations and local electrostatic environments. The vibrational bands of azide probes are often accompanied by complex line shapes with shoulder peaks, which may arise either from incomplete background subtraction, Fermi resonance, or multiple conformers. The isotope substitution in the infrared probe has thus been introduced to remove Fermi resonances without causing a significant perturbation to the structure. Here, we synthesized and labeled the mid-N atoms of aliphatic azide derivatives with 15N to study the effects of isotope labelling on their vibrational properties. The FT-IR spectra of the aliphatic azide with asymmetric lineshape became a single symmetric band upon isotope substitution, which might be an indication of the removal of the hidden Fermi resonance from the system. We also noticed that the 2D-IR spectrum of unlabeled aliphatic azide has cross-peaks, even though it is not apparently identifiable. The 1D slice spectra obtained from the 2D-IR spectra reveal the existence of a hidden Fermi resonance peak. Furthermore, we show that this weak Fermi resonance does not produce discernible oscillatory beating patterns in the IR pump-probe spectrum, which has been used as evidence of the Fermi resonance. Therefore, we confirm that isotope labelling combined with 2D-IR spectroscopy is the most efficient and incisive way to identify the origin of small shoulder peaks in the linear and nonlinear vibrational spectra of various IR probe molecules.

Real-time tracking of protein unfolding with time-resolved x-ray solution scattering

The correct folding of proteins is of paramount importance for their function, and protein misfolding is believed to be the primary cause of a wide range of diseases. Protein folding has been investigated with time-averaged methods and time-resolved spectroscopy, but observing the structural dynamics of the unfolding process in real-time is challenging. Here, we demonstrate an approach to directly reveal the structural changes in the unfolding reaction. We use nano- to millisecond time-resolved x-ray solution scattering to probe the unfolding of apomyoglobin. The unfolding reaction was triggered using a temperature jump, which was induced by a nanosecond laser pulse. We demonstrate a new strategy to interpret time-resolved x-ray solution scattering data, which evaluates ensembles of structures obtained from molecular dynamics simulations. We find that apomyoglobin passes three states when unfolding, which we characterize as native, molten globule, and unfolded. The molten globule dominates the population under the conditions investigated herein, whereas native and unfolded structures primarily contribute before the laser jump and 30 μs after it, respectively. The molten globule retains much of the native structure but shows a dynamic pattern of inter-residue contacts. Our study demonstrates a new strategy to directly observe structural changes over the cause of the unfolding reaction, providing time- and spatially resolved atomic details of the folding mechanism of globular proteins.

Ultrafast pH-jump two-dimensional infrared spectroscopy

We present a pH-jump two-dimensional infrared (2D IR) spectrometer to probe pH-dependent conformational changes from nanoseconds to milliseconds. The design incorporates a nanosecond 355 nm source into a pulse-shaper-based 2D IR spectrometer to trigger dissociation of a caged proton prior to probing subsequent conformational changes with femtosecond 2D IR spectroscopy. We observe a blue shift in the amide I mode (C═O stretch) of diglycine induced by protonation of the terminal amine. This method combines the bond-specific structural sensitivity of ultrafast 2D IR with triggered conformational dynamics, providing structural access to multiscale biomolecular transformations such as protein folding.

Vibrational Coupling on Stepwise Hydrogen Bond Formation of Amide I

Despite the key roles of proteins and nucleic acids in biology, understanding their labile structures and hydrogen bond interactions with guest molecules has posed a critical challenge to the scientific community. In this report, I have used dimethylformamide as a model amide to account for amide hydrogen bond interactions of protein. To quantify hydrogen bond conformation and the structural change, I have monitored the amide I infrared (IR) stretching frequencies while varying the pK_a of phenol derivatives. For all phenol derivatives, amide I has formed one hydrogen bond and two hydrogen bond conformation. It has been observed that the formation constant for one hydrogen bond is higher than that of two hydrogen bonds for all phenol derivatives. During the formation of hydrogen bond with amide I, IR absorbance of C═C transition is enhanced for all phenol derivatives. Enhancement of the IR absorbance of the C═C transition indicates hydrogen bond-assisted vibrational coupling between the amide I and phenol ring transition. The relative coupling constant is estimated to be higher for single hydrogen-bonded conformer than the double hydrogen-bonded conformer. This is an intriguing result as the frequency difference between the two coupled transitions predicts otherwise. Using IR absorption spectroscopy, a delicate interplay between hydrogen bonding conformations and intermolecular vibrational coupling between amide I and H-bond donor phenol molecules has been shown. This study can be used as a point of reference for understanding the structural information of proteins, peptides, and nucleosides having hydrogen bond interaction with any drug or ligand molecules. My results as well provide an insight into the vibrational coupling of carbonyl and C═C transition of nucleobases.

Application of millisecond time-resolved solid state NMR to the kinetics and mechanism of melittin self-assembly

Common experimental approaches for characterizing structural conversion processes such as protein folding and self-assembly do not report on all aspects of the evolution from an initial state to the final state. Here, we demonstrate an approach that is based on rapid mixing, freeze-trapping, and low-temperature solid-state NMR (ssNMR) with signal enhancements from dynamic nuclear polarization (DNP). Experiments on the folding and tetramerization of the 26-residue peptide melittin following a rapid pH jump show that multiple aspects of molecular structure can be followed with millisecond time resolution, including secondary structure at specific isotopically labeled sites, intramolecular and intermolecular contacts between specific pairs of labeled residues, and overall structural order. DNP-enhanced ssNMR data reveal that conversion of conformationally disordered melittin monomers at low pH to α-helical conformations at neutral pH occurs on nearly the same timescale as formation of antiparallel melittin dimers, about 6 to 9 ms for 0.3 mM melittin at 24 °C in aqueous solution containing 20% (vol/vol) glycerol and 75 mM sodium phosphate. Although stopped-flow fluorescence data suggest that melittin tetramers form quickly after dimerization, ssNMR spectra show that full structural order within melittin tetramers develops more slowly, in ∼60 ms. Time-resolved ssNMR is likely to find many applications to biomolecular structural conversion processes, including early stages of amyloid formation, viral capsid formation, and protein-protein recognition.

Non-Additive Effects of Binding Site Mutations in Calmodulin

Despite decades of research on ion-sensing proteins, gaps persist in the understanding of ion binding affinity and selectivity even in well-studied proteins such as calmodulin. Site-directed mutagenesis is a powerful and popular tool for addressing outstanding questions about biological ion binding and is employed to selectively deactivate binding sites and insert chromophores at advantageous positions within ion binding structures. However, even apparently nonperturbative mutations can distort the binding dynamics they are employed to measure. We use Fourier transform infrared (FTIR) and ultrafast two-dimensional infrared (2D IR) spectroscopy of the carboxylate asymmetric stretching mode in calmodulin as a mutation- and label-independent probe of the conformational perturbations induced in calmodulin's binding sites by two classes of mutation, tryptophan insertion and carboxylate side-chain deletion, commonly used to study ion binding in proteins. Our results show that these mutations not only affect ion binding but also induce changes in calmodulin's conformational landscape along coordinates not probed by vibrational spectroscopy, remaining invisible without additional perturbation of binding site structure. Comparison of FTIR line shapes with 2D IR diagonal slices provides a clear example of how nonlinear spectroscopy produces well-resolved line shapes, refining otherwise featureless spectral envelopes into more informative vibrational spectra of proteins.

Reversible control of enzyme-inhibitor interactions with light illumination using a photoresponsive surfactant

The effects of a photoresponsive surfactant and light illumination on the complex formed between ribonuclease A (RNase A) and a protein ribonuclease inhibitor (RI) have been investigated to develop a light-based technique to reactivate an enzyme through surfactant-induced dissociation of the enzyme-inhibitor complex. The photoresponsive surfactant undergoes a photoisomerization from the relatively hydrophobic trans isomer under visible light to the relatively hydrophilic cis isomer upon UV illumination, providing a means to reversibly control protein-inhibitor interactions. In the absence of surfactant, RI binds tightly to RNase A through noncovalent interactions, which inhibits the enzyme activity. Upon addition of the surfactant under visible light, RNase A is reactivated, regaining ~75% of the native activity in the absence of RI. In the presence of the surfactant under UV light, however, the enzyme remains inhibited. Fluorescence spectroscopy, dynamic light scattering, and circular dichroism spectroscopy reveal that RI dramatically unfolds upon addition of the trans form of the surfactant, while RNase A does not undergo noticeable structural changes under the same conditions. This indicates that RNase A reactivation occurs through dissociation of the enzyme-inhibitor complex arising from surfactant-induced unfolding of the inhibitor. As a result, photoresponsive surfactant and light illumination can be used as a novel light-based technique to dissociate enzyme-inhibitor complexes and, thus, reactivate an inhibited enzyme.

Site-Specific Peptide Probes Detect Buried Water in a Lipid Membrane

Transmembrane peptides contain polar residues in the interior of the membrane, which may alter the electrostatic environment and favor hydration in the otherwise nonpolar environment of the membrane core. Here, we demonstrate a general, nonperturbative strategy to probe hydration of the peptide backbone at specific depths within the bilayer using a combination of site-specific isotope labels, ultrafast two-dimensional infrared spectroscopy, and spectral modeling based on molecular dynamics simulations. Our results show that the amphiphilic pH-low insertion peptide supports a highly heterogeneous environment, with significant backbone hydration of nonpolar residues neighboring charged residues. For example, a leucine residue located as far as 1 nm into the hydrophobic bulk reports hydrogen-bonded populations as high as ∼20%. These findings indicate that the polar nature of these residues may facilitate the transport of water molecules into the hydrophobic core of the membrane.

Exposing the Nucleation Site in α-Helix Folding: A Joint Experimental and Simulation Study

One of the fundamental events in protein folding is α-helix formation, which involves sequential development of a series of helical hydrogen bonds between the backbone C═O group of residues i and the -NH group of residues i + 4. While we now know a great deal about α-helix folding dynamics, a key question that remains to be answered is where the productive helical nucleation event occurs. Statistically, a helical nucleus (or the first helical hydrogen-bond) can form anywhere within the peptide sequence in question; however, the one that leads to productive folding may only form at a preferred location. This consideration is based on the fact that the α-helical structure is inherently asymmetric, due to the specific alignment of the helical hydrogen bonds. While this hypothesis is plausible, validating it is challenging because there is not an experimental observable that can be used to directly pinpoint the location of the productive nucleation process. Therefore, in this study we combine several techniques, including peptide cross-linking, laser-induced temperature-jump infrared spectroscopy, and molecular dynamics simulations, to tackle this challenge. Taken together, our experimental and simulation results support an α-helix folding mechanism wherein the productive nucleus is formed at the N-terminus, which propagates toward the C-terminal end of the peptide to yield the folded structure. In addition, our results show that incorporation of a cross-linker can lead to formation of differently folded conformations, underscoring the need for all-atom simulations to quantitatively assess the proposed cross-linking design.

4-Cyanoindole-2'-deoxyribonucleoside as a Dual Fluorescence and Infrared Probe of DNA Structure and Dynamics

Unnatural nucleosides possessing unique spectroscopic properties that mimic natural nucleobases in both size and chemical structure are ideally suited for spectroscopic measurements of DNA/RNA structure and dynamics in a site-specific manner. However, such unnatural nucleosides are scarce, which prompts us to explore the utility of a recently found unnatural nucleoside, 4-cyanoindole-2'-deoxyribonucleoside (4CNI-NS), as a site-specific spectroscopic probe of DNA. A recent study revealed that 4CNI-NS is a universal nucleobase that maintains the high fluorescence quantum yield of 4-cyanoindole and that among the four natural nucleobases, only guanine can significantly quench its fluorescence. Herein, we further show that the C≡N stretching frequency of 4CNI-NS is sensitive to the local environment, making it a useful site-specific infrared probe of oligonucleotides. In addition, we demonstrate that the fluorescence-quencher pair formed by 4CNI-NS and guanine can be used to quantitatively assess the binding affinity of a single-stranded DNA to the protein system of interest via fluorescence spectroscopy, among other applications. We believe that this fluorescence binding assay is especially useful as its potentiality allows high-throughput screening of DNA⁻protein interactions.

Linear and Nonlinear Infrared Spectroscopies Reveal Detailed Solute-Solvent Dynamic Interactions of a Nitrosyl Ruthenium Complex in Solution

In this work, the solvation of a nitrosyl ruthenium complex, [(CH₃)₄N][RuCl₃(qn)(NO)] (with qn = deprotonated 8-hydroxyquinoline), which is a potential NO-releasing molecule in the bio-environment, was studied in two bio-friendly solvents, namely deuterated dimethyl sulfoxide (dDMSO) and water (D₂O). A blue-shifted NO stretching frequency was observed in water with respect to that in dDMSO, which was believed to be due to ligand-solvent hydrogen-bonding interactions, one N═O···D and particularly three Ru-Cl···D, that show competing effects on the NO bond length. The dynamic differences of the NO stretch in these two solvents were further revealed by transient pump-probe IR and two-dimensional IR results: faster vibrational relaxation and faster spectral diffusion (SD) were observed in D₂O, confirming stronger solvent-solute interaction and also faster solvent structural dynamics in D₂O than in DMSO. Further, a significant non-decaying residual in the SD dynamics was observed in D₂O but not in DMSO, suggesting the formation of a stable solvation shell in water due to strong multi-site ligand-solvent hydrogen-bonding interactions, which is in agreement with the observed blue-shifted NO stretching frequency. This work demonstrates that small solvent molecules such as water can form a relatively rigid solvation shell for certain transition metal complexes due to cooperative ligand-solvent interactions and show slower dynamics.

Equilibrium versus Nonequilibrium Peptide Dynamics: Insights into Transient 2D IR Spectroscopy

Over the past two decades, two-dimensional infrared (2D IR) spectroscopy has evolved from the theoretical underpinnings of nonlinear spectroscopy as a means of investigating detailed molecular structure on an ultrafast time scale. The combined time and spectral resolution over which spectra can be collected on complex molecular systems has led to the precise structural resolution of dynamic species that have previously been impossible to directly observe through traditional methods. The adoption of 2D IR spectroscopy for the study of protein folding and peptide interactions has provided key details of how small changes in conformations can exert major influences on the activities of these complex molecular systems. Traditional 2D IR experiments are limited to molecules under equilibrium conditions, where small motions and fluctuations of these larger molecules often still lead to functionality. Utilizing techniques that allow the rapid initiation of chemical or structural changes in conjunction with 2D IR spectroscopy, i.e., transient 2D IR, a vast dynamic range becomes available to the spectroscopist uncovering structural content far from equilibrium. Furthermore, this allows the observation of reaction pathways of these macromolecules under quasi- and nonequilibrium conditions.

Coordination to lanthanide ions distorts binding site conformation in calmodulin

The Ca²⁺-sensing protein calmodulin (CaM) is a popular model of biological ion binding since it is both experimentally tractable and essential to survival in all eukaryotic cells. CaM modulates hundreds of target proteins and is sensitive to complex patterns of Ca²⁺ exposure, indicating that it functions as a sophisticated dynamic transducer rather than a simple on/off switch. Many details of this transduction function are not well understood. Fourier transform infrared (FTIR) spectroscopy, ultrafast 2D infrared (2D IR) spectroscopy, and electronic structure calculations were used to probe interactions between bound metal ions (Ca²⁺ and several trivalent lanthanide ions) and the carboxylate groups in CaM's EF-hand ion-coordinating sites. Since Tb³⁺ is commonly used as a luminescent Ca²⁺ analog in studies of protein-ion binding, it is important to characterize distinctions between the coordination of Ca²⁺ and the lanthanides in CaM. Although functional assays indicate that Tb³⁺ fully activates many Ca²⁺-dependent proteins, our FTIR spectra indicate that Tb³⁺, La³⁺, and Lu³⁺ disrupt the bidentate coordination geometry characteristic of the CaM binding sites' strongly conserved position 12 glutamate residue. The 2D IR spectra indicate that, relative to the Ca²⁺-bound form, lanthanide-bound CaM exhibits greater conformational flexibility and larger structural fluctuations within its binding sites. Time-dependent 2D IR lineshapes indicate that binding sites in Ca²⁺-CaM occupy well-defined configurations, whereas binding sites in lanthanide-bound-CaM are more disordered. Overall, the results show that binding to lanthanide ions significantly alters the conformation and dynamics of CaM's binding sites.

Polymers and biopolymers at interfaces

This review updates recent progress in the understanding of the behaviour of polymers at surfaces and interfaces, highlighting examples in the areas of wetting, dewetting, crystallization, and 'smart' materials. Recent developments in analysis tools have yielded a large increase in the study of biological systems, and some of these will also be discussed, focussing on areas where surfaces are important. These areas include molecular binding events and protein adsorption as well as the mapping of the surfaces of cells. Important techniques commonly used for the analysis of surfaces and interfaces are discussed separately to aid the understanding of their application.

Mass spectrometric characterization of protein structures and protein complexes in condensed and gas phase

Proteins are essential for almost all physiological processes of life. They serve a myriad of functions which are as varied as their unique amino acid sequences and their corresponding three-dimensional structures. To fulfill their tasks, most proteins depend on stable physical associations, in the form of protein complexes that evolved between themselves and other proteins. In solution (condensed phase), proteins and/or protein complexes are in constant energy exchange with the surrounding solvent. Albeit methods to describe in-solution thermodynamic properties of proteins and of protein complexes are well established and broadly applied, they do not provide a broad enough access to life-science experimentalists to study all their proteins' properties at leisure. This leaves great desire to add novel methods to the analytical biochemist's toolbox. The development of electrospray ionization created the opportunity to characterize protein higher order structures and protein complexes rather elegantly by simultaneously lessening the need of sophisticated sample preparation steps. Electrospray mass spectrometry enabled us to translate proteins and protein complexes very efficiently into the gas phase under mild conditions, retaining both, intact protein complexes, and gross protein structures upon phase transition. Moreover, in the environment of the mass spectrometer (gas phase, in vacuo), analyte molecules are free of interactions with surrounding solvent molecules and, therefore, the energy of inter- and intramolecular forces can be studied independently from interference of the solvating environment. Provided that gas phase methods can give information which is relevant for understanding in-solution processes, gas phase protein structure studies and/or investigations on the characterization of protein complexes has rapidly gained more and more attention from the bioanalytical scientific community. Recent reports have shown that electrospray mass spectrometry provides direct access to six prime protein complex properties: stabilities, compositions, binding surfaces (epitopes), disassembly processes, stoichiometries, and thermodynamic parameters.

Meandering Down the Energy Landscape of Protein Folding: Are We There Yet?

As judged by a single publication metric, the activity in the protein folding field has been declining over the past 5 years, after enjoying a decade-long growth. Does this development indicate that the field is sunsetting or is this decline only temporary? Upon surveying a small territory of its landscape, we find that the protein folding field is still quite active and many important findings have emerged from recent experimental studies. However, it is also clear that only continued development of new techniques and methods, especially those enabling dissection of the fine details and features of the protein folding energy landscape, will fuel this old field to move forward.

Dual time-resolved temperature-jump fluorescence and infrared spectroscopy for the study of fast protein dynamics

Time-resolved temperature-jump (T-jump) coupled with fluorescence and infrared (IR) spectroscopy is a powerful technique for monitoring protein dynamics. Although IR spectroscopy of the polypeptide amide I mode is more technically challenging, it offers complementary information because it directly probes changes in the protein backbone, whereas, fluorescence spectroscopy is sensitive to the environment of specific side chains. With the advent of widely tunable quantum cascade lasers (QCL) it is possible to efficiently probe multiple IR frequencies with high sensitivity and reproducibility. Here we describe a dual time-resolved T-jump fluorescence and IR spectrometer and its application to study protein folding dynamics. A Q-switched Ho:YAG laser provides the T-jump source for both time-resolved IR and fluorescence spectroscopy, which are probed by a QCL and Ti:Sapphire laser, respectively. The Ho:YAG laser simultaneously pumps the time-resolved IR and fluorescence spectrometers. The instrument has high sensitivity, with an IR absorbance detection limit of <0.2mOD and a fluorescence sensitivity of 2% of the overall fluorescence intensity. Using a computer controlled QCL to rapidly tune the IR frequency it is possible to create a T-jump induced difference spectrum from 50ns to 0.5ms. This study demonstrates the power of the dual time-resolved T-jump fluorescence and IR spectroscopy to resolve complex folding mechanisms by complementary IR absorbance and fluorescence measurements of protein dynamics.

Microscopic nucleation and propagation rates of an alanine-based α-helix

An infrared temperature-jump (T-jump) study by Huang et al. (Proc. Natl. Acad. Sci. U. S. A., 2002, 99, 2788-2793) showed that the conformational relaxation kinetics of an alanine-based α-helical peptide depend not only on the final temperature (Tf) but also on the initial temperature (Ti) when Tf is fixed. Their finding indicates that the folding free energy landscape of this peptide is non-two-state like, allowing for the population of conformational ensembles with different helical lengths and relaxation times in the temperature range of the experiment. Because α-helix folding involves two fundamental events, nucleation and propagation, the results of Huang et al. thus present a unique opportunity to determine their rate constants - a long-sought goal in the study of the helix-coil transition dynamics. Herein, we capitalize on this notion and develop a coarse-grained kinetic model to globally fit the thermal unfolding curve and T-jump kinetic traces of this peptide. Using this strategy, we are able to explicitly determine the microscopic rate constants of the kinetic steps encountered in the nucleation and propagation processes. Our results reveal that the time taken to form one α-helical turn (i.e., an α-helical segment with one helical hydrogen bond) is about 315 ns, whereas the time taken to elongate this nucleus by one residue (or backbone unit) is 5.9 ns, depending on the position of the residue.

Probing the Acid-Induced Packing Structure Changes of the Molten Globule Domains of a Protein near Equilibrium Unfolding

Using simultaneously scanning small-angle X-ray scattering (SAXS) and UV-vis absorption with integrated online size exclusion chromatography, supplemental with molecular dynamics simulations, we unveil the long-postulated global structure evolution of a model multidomain protein bovine serum albumin (BSA) during acid-induced unfolding. Our results differentiate three global packing structures of the three molten globule domains of BSA, forming three intermediates I₁, I₂, and E along the unfolding pathway. The I₁-I₂ transition, overlooked in all previous studies, involves mainly coordinated reorientations across interconnected molten globule subdomains, and the transition activates a critical pivot domain opening of the protein for entering into the E form, with an unexpectedly large unfolding free energy change of -9.5 kcal mol^-1, extracted based on the observed packing structural changes. The revealed local packing flexibility and rigidity of the molten globule domains in the E form elucidate how collective motions of the molten globule domains profoundly influence the folding-unfolding pathway of a multidomain protein.

Simple method to introduce an ester infrared probe into proteins

The ester carbonyl stretching vibration has recently been shown to be a sensitive and convenient infrared (IR) probe of protein electrostatics due to the linear dependence of its frequency on local electric field. While an ester moiety can be easily incorporated into peptides via solid-phase synthesis, currently there is no method available to site-specifically incorporate it into a large protein. Herein, we show that it is possible to use a cysteine alkylation reaction to achieve this goal and demonstrate the feasibility of this simple method by successfully incorporating a methyl ester group (CH₂ COOCH₃ ) into a model peptide (YGGCGG), two amyloid-forming peptides derived from the insulin B chain and Aβ, and bovine serum albumin (BSA). IR results obtained with those peptide and protein systems further confirm the utility of this vibrational probe in monitoring, for example, the structural integrity of amyloid fibrils and ligand binding-induced changes in protein local hydration status.

Kinetics of peptide folding in lipid membranes

Despite our extensive understanding of water-soluble protein folding kinetics, much less is known about the folding dynamics and mechanisms of membrane proteins. However, recent studies have shown that for relatively simple systems, such as peptides that form a transmembrane α-helix, helical dimer, or helix-turn-helix, it is possible to assess the kinetics of several important steps, including peptide binding to the membrane from aqueous solution, peptide folding on the membrane surface, helix insertion into the membrane, and helix-helix association inside the membrane. Herein, we provide a brief review of these studies and also suggest new initiation and probing methods that could lead to improved temporal and structural resolution in future experiments.

Infrared and Fluorescence Assessment of Protein Dynamics: From Folding to Function

While folding or performing functions, a protein can sample a rich set of conformational space. However, experimentally capturing all of the important motions with sufficient detail to allow a mechanistic description of their dynamics is nontrivial since such conformational events often occur over a wide range of time and length scales. Therefore, many methods have been employed to assess protein conformational dynamics, and depending on the nature of the conformational transition in question, some may be more advantageous than others. Herein, we describe our recent efforts, and also those of others, wherever appropriate, to use infrared- and fluorescence-based techniques to interrogate protein folding and functional dynamics. Specifically, we focus on discussing how to use extrinsic spectroscopic probes to enhance the structural resolution of these techniques and how to exploit various cross-linking strategies to acquire dynamic and mechanistic information that was previously difficult to attain.

Direct measurement of the tryptophan-mediated photocleavage kinetics of a protein disulfide bond

Disulfide cleavage is one of the major causes underlying ultraviolet (UV) light-induced protein damage. While previous studies have provided strong evidence to support the notion that this process is mediated by photo-induced electron transfer from the excited state of an aromatic residue (e.g., tryptophan) to the disulfide bond, many mechanistic details are still lacking. For example, we do not know how quickly this process occurs in a protein environment. Herein, we design an experiment, which uses the unfolding kinetics of a protein as an observable, to directly assess the kinetics and mechanism of photo-induced disulfide cleavage. Our results show that this disulfide bond cleavage event takes place in ∼2 μs via a mechanism involving electron transfer from the triplet state of a tryptophan (Trp) residue to the disulfide bond. Furthermore, we find that one of the photoproducts of this reaction, a Trp-SR adduct, is formed locally, thus preventing the protein from re-cross-linking. Taken together, these findings suggest that a Trp-disulfide pair could be used as a photo-trigger to initiate protein folding dynamics and control the biological activities of disulfide-containing peptides.

Probing the Folding-Unfolding Transition of a Thermophilic Protein, MTH1880

The folding mechanism of typical proteins has been studied widely, while our understanding of the origin of the high stability of thermophilic proteins is still elusive. Of particular interest is how an atypical thermophilic protein with a novel fold maintains its structure and stability under extreme conditions. Folding-unfolding transitions of MTH1880, a thermophilic protein from Methanobacterium thermoautotrophicum, induced by heat, urea, and GdnHCl, were investigated using spectroscopic techniques including circular dichorism, fluorescence, NMR combined with molecular dynamics (MD) simulations. Our results suggest that MTH1880 undergoes a two-state N to D transition and it is extremely stable against temperature and denaturants. The reversibility of refolding was confirmed by spectroscopic methods and size exclusion chromatography. We found that the hyper-stability of the thermophilic MTH1880 protein originates from an extensive network of both electrostatic and hydrophobic interactions coordinated by the central β-sheet. Spectroscopic measurements, in combination with computational simulations, have helped to clarify the thermodynamic and structural basis for hyper-stability of the novel thermophilic protein MTH1880.

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.

How Sensitive is the Amide I Vibration of the Polypeptide Backbone to Electric Fields?

Site-selective isotopic labeling of amide carbonyls offers a nonperturbative means to introduce a localized infrared probe into proteins. Although this strategy has been widely used to investigate various biological questions, the dependence of the underlying amide I vibrational frequency on electric fields (or Stark tuning rate) has not been fully determined, which prevents it from being used in a quantitative manner in certain applications. Herein, through the use of experiments and molecular dynamics simulations, the Stark tuning rate of the amide I vibration of an isotopically labeled backbone carbonyl in a transmembrane α-helix is determined to be approximately 1.4 cm(-1) /(MV/cm). This result provides a quantitative basis for using this vibrational model to assess local electric fields in proteins, among other applications. For instance, by using this value, we are able to show that the backbone region of a dipeptide has a surprisingly low dielectric constant.

Modeling Vibrational Spectra of Ester Carbonyl Stretch in Water and DMSO Based on Molecular Dynamics Simulation

On the basis of molecular dynamics simulation, we model the ester carbonyl stretch FTIR signals of methyl acetate in D2O and DMSO. An ab initio map is constructed at the B3LYP/6-311++G** level to relate the carbonyl stretch frequency to the external electric field. Using this map, fluctuating Hamiltonian of the carbonyl stretch is constructed from the MD simulation trajectory. The IR spectra calculated based on this Hamiltonian are found to be in good agreement with the experiment. For methyl acetate in D2O, hydrogen bonding on alkoxy oxygen causes a blue shift of frequency, while that on carbonyl oxygen causes a red shift. Two peaks observed in FTIR signals originate from the balance of these two effects. Furthermore, in both D2O and DMSO solutions, correlations are found between the instantaneous electric field on C═O and the frequencies. Broader line width of the signal in D2O suggests a more inhomogeneous electric field distribution due to the complicated hydrogen-bonding environment.

Biomolecular Crowding Arising from Small Molecules, Molecular Constraints, Surface Packing, and Nano-Confinement

The effect of macromolecular crowding on the structure, dynamics, and reactivity of biomolecules is well established and the relevant research has been extensively reviewed. Herein, we focus our discussion on crowding effects arising from small cosolvent molecules and densely packed surface conditions. In addition, we highlight recent efforts that capitalize on the excluded volume effect for various tailored biochemical and biophysical applications. Specifically, we discuss how a targeted increase in local mass density can be exploited to gain insight into the folding dynamics of the protein of interest and how confinement via reverse micelles can be used to study a range of biophysical questions, from protein hydration dynamics to amyloid formation.

Interaction Networks in Protein Folding via Atomic-Resolution Experiments and Long-Time-Scale Molecular Dynamics Simulations

The integration of atomic-resolution experimental and computational methods offers the potential for elucidating key aspects of protein folding that are not revealed by either approach alone. Here, we combine equilibrium NMR measurements of thermal unfolding and long molecular dynamics simulations to investigate the folding of gpW, a protein with two-state-like, fast folding dynamics and cooperative equilibrium unfolding behavior. Experiments and simulations expose a remarkably complex pattern of structural changes that occur at the atomic level and from which the detailed network of residue-residue couplings associated with cooperative folding emerges. Such thermodynamic residue-residue couplings appear to be linked to the order of mechanistically significant events that take place during the folding process. Our results on gpW indicate that the methods employed in this study are likely to prove broadly applicable to the fine analysis of folding mechanisms in fast folding proteins.

Sensing pH via p-cyanophenylalanine fluorescence: Application to determine peptide pKa and membrane penetration kinetics

We expand the spectroscopic utility of a well-known infrared and fluorescence probe, p-cyanophenylalanine, by showing that it can also serve as a pH sensor. This new application is based on the notion that the fluorescence quantum yield of this unnatural amino acid, when placed at or near the N-terminal end of a polypeptide, depends on the protonation status of the N-terminal amino group of the peptide. Using this pH sensor, we are able to determine the N-terminal pKa values of nine tripeptides and also the membrane penetration kinetics of a cell-penetrating peptide. Taken together, these examples demonstrate the applicability of using this unnatural amino acid fluorophore to study pH-dependent biological processes or events that accompany a pH change.

Tightening up the structure, lighting up the pathway: Application of molecular constraints and light to manipulate protein folding, self-assembly and function

Chemical cross-linking provides an effective avenue to reduce the conformational entropy of polypeptide chains and hence has become a popular method to induce or force structural formation in peptides and proteins. Recently, other types of molecular constraints, especially photoresponsive linkers and functional groups, have also found increased use in a wide variety of applications. Herein, we provide a concise review of using various forms of molecular strategies to constrain proteins, thereby stabilizing their native states, gaining insight into their folding mechanisms, and/or providing a handle to trigger a conformational process of interest with light. The applications discussed here cover a wide range of topics, ranging from delineating the details of the protein folding energy landscape to controlling protein assembly and function.

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.

Novel microscale approaches for easy, rapid determination of protein stability in academic and commercial settings

Chemical denaturant titrations can be used to accurately determine protein stability. However, data acquisition is typically labour intensive, has low throughput and is difficult to automate. These factors, combined with high protein consumption, have limited the adoption of chemical denaturant titrations in commercial settings. Thermal denaturation assays can be automated, sometimes with very high throughput. However, thermal denaturation assays are incompatible with proteins that aggregate at high temperatures and large extrapolation of stability parameters to physiological temperatures can introduce significant uncertainties. We used capillary-based instruments to measure chemical denaturant titrations by intrinsic fluorescence and microscale thermophoresis. This allowed higher throughput, consumed several hundred-fold less protein than conventional, cuvette-based methods yet maintained the high quality of the conventional approaches. We also established efficient strategies for automated, direct determination of protein stability at a range of temperatures via chemical denaturation, which has utility for characterising stability for proteins that are difficult to purify in high yield. This approach may also have merit for proteins that irreversibly denature or aggregate in classical thermal denaturation assays. We also developed procedures for affinity ranking of protein–ligand interactions from ligand-induced changes in chemical denaturation data, and proved the principle for this by correctly ranking the affinity of previously unreported peptide–PDZ domain interactions. The increased throughput, automation and low protein consumption of protein stability determinations afforded by using capillary-based methods to measure denaturant titrations, can help to revolutionise protein research. We believe that the strategies reported are likely to find wide applications in academia, biotherapeutic formulation and drug discovery programmes.

•

Chemical denaturant titrations are slow, lengthy and consume lots of protein sample.

•

We developed very fast methods for measuring chemical denaturant titrations.

•

Automated titrations can be measured in minutes using only μl sample volumes.

•

Ligand interactions rapidly screened and ranked via changes in protein stability

•

These advances make chemical denaturation more suitable for commercial research.

Experimental validation of the role of trifluoroethanol as a nanocrowder

Trifluoroethanol (TFE) is commonly used to induce protein secondary structure, especially α-helix formation. Due to its amphiphilic nature, however, TFE can also self-associate to form micellelike, nanometer-sized clusters. Herein, we hypothesize that such clusters can act as nanocrowders to increase protein folding rates via the excluded volume effect. To test this hypothesis, we measure the conformational relaxation kinetics of an intrinsically disordered protein, the phosphorylated kinase inducible domain (pKID), which forms a helix–turn–helix in TFE solutions. We find that the conformational relaxation rate of pKID displays a rather complex dependence on TFE percentage (v/v): while it first decreases between 0 and 5%, between 5 and 15% the rate increases and then remains relatively unchanged between 15 and 30% and finally decreases again at higher percentages (i.e., 50%). This trend coincides with the fact that TFE clustering is maximized in the range of 15–30%, thus providing validation of our hypothesis. Another line of supporting evidence comes from the observation that the relaxation rate of a monomeric helical peptide, which due to its predominantly local interactions in the folded state is less affected by crowding, does not show a similar TFE dependence.

Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium-binding protein exhibits a broad range of structural transitions, including domain opening and folding-unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well-known, it has been challenging to examine the dynamics of transition between these structures in real-time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C-terminal domain in response to Ca(2+) binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation-sensitive tyrosine residues, opening a path to time-resolved studies of CaM dynamics using Raman spectroscopy.

Ester carbonyl vibration as a sensitive probe of protein local electric field

The ability to quantify the local electrostatic environment of proteins and protein/peptide assemblies is key to gaining a microscopic understanding of many biological interactions and processes. Herein, we show that the ester carbonyl stretching vibration of two non-natural amino acids, L-aspartic acid 4-methyl ester and L-glutamic acid 5-methyl ester, is a convenient and sensitive probe in this regard, since its frequency correlates linearly with the local electrostatic field for both hydrogen-bonding and non-hydrogen-bonding environments. We expect that the resultant frequency-electric-field map will find use in various applications. Furthermore, we show that, when situated in a non-hydrogen-bonding environment, this probe can also be used to measure the local dielectric constant (ε). For example, its application to amyloid fibrils formed by Aβ(16-22) revealed that the interior of such β-sheet assemblies has an ε value of approximately 5.6.