Compaction during Protein Folding Studied by Real-Time NMR Diffusion Experiments

Translational diffusion of unfolded and intrinsically disordered proteins

Translational (or self-diffusion) coefficient in dilute solution is inversely proportional to the size of a diffusing molecule, and hence self-diffusion coefficient measurements have been applied to determine the effective hydrodynamic radii for a range of native and nonnative protein conformations. In particular, translational diffusion coefficient measurements are useful to estimate the hydrodynamic radius of natively (or intrinsically) disordered proteins in solution, and, thereby, probe the compactness of a protein as well as its change when environmental parameters such as temperature, solution pH, or protein concentration are varied. The situation becomes more complicated in concentrated solutions. In this review, we discuss the translational diffusion of disordered proteins in dilute and crowded solutions, focusing primarily on the information provided by pulsed-field gradient NMR technique, and draw analogies to well-structured globular proteins and synthetic polymers.

Hydrophobic Collapse in N-Methylacetamide-Water Mixtures

Aqueous N-methylacetamide solutions were investigated by polarization-resolved pump-probe and 2D infrared spectroscopy (2D IR), using the amide I mode as a reporter. The 2D IR results are compared with molecular dynamics simulations and spectral calculations to gain insight into the molecular structures in the mixture. N-Methylacetamide and water molecules tend to form clusters with "frozen" amide I dynamics. This is driven by a hydrophobic collapse as the methyl groups of the N-methylacetamide molecules cluster in the presence of water. Since the studied system can be considered as a simplified model for the backbone of proteins, the present study forms a convenient basis for understanding the structural and vibrational dynamics in proteins. It is particularly interesting to find out that a hydrophobic collapse as the one driving protein folding is observed in such a simple system.

Theoretical and Experimental Investigation of the Translational Diffusion of Proteins in the Vicinity of Temperature-Induced Unfolding Transition

Translational diffusion is the most fundamental form of transport in chemical and biological systems. The diffusion coefficient is highly sensitive to changes in the size of the diffusing species; hence, it provides important information on the variety of macromolecular processes, such as self-assembly or folding-unfolding. Here, we investigate the behavior of the diffusion coefficient of a macromolecule in the vicinity of heat-induced transition from folded to unfolded state. We derive the equation that describes the diffusion coefficient of the macromolecule in the vicinity of the transition and use it to fit the experimental data from pulsed-field-gradient nuclear magnetic resonance (PFG NMR) experiments acquired for two globular proteins, lysozyme and RNase A, undergoing temperature-induced unfolding. A very good qualitative agreement between the theoretically derived diffusion coefficient and experimental data is observed.

Probing the non-native H helix translocation in apomyoglobin folding intermediates

Apomyoglobin folds via sequential helical intermediates that are formed by rapid collapse of the A, B, G, and H helix regions. An equilibrium molten globule with a similar structure is formed near pH 4. Previous studies suggested that the folding intermediates are kinetically trapped states in which folding is impeded by non-native packing of the G and H helices. Fluorescence spectra of mutant proteins in which cysteine residues were introduced at several positions in the G and H helices show differential quenching of W14 fluorescence, providing direct evidence of translocation of the H helix relative to helices A and G in both the kinetic and equilibrium intermediates. Förster resonance energy transfer measurements show that a 5-({2-[(acetyl)amino]ethyl}amino)naphthalene-1-sulfonic acid acceptor coupled to K140C (helix H) is closer to Trp14 (helix A) in the equilibrium molten globule than in the native state, by a distance that is consistent with sliding of the H helix in an N-terminal direction by approximately one helical turn. Formation of an S108C–L135C disulfide prevents H helix translocation in the equilibrium molten globule by locking the G and H helices into their native register. By enforcing nativelike packing of the A, G, and H helices, the disulfide resolves local energetic frustration and facilitates transient docking of the E helix region onto the hydrophobic core but has only a small effect on the refolding rate. The apomyoglobin folding landscape is highly rugged, with several energetic bottlenecks that frustrate folding; relief of any one of the major identified bottlenecks is insufficient to speed progression to the transition state.

Simulating the peptide folding kinetic related spectra based on the Markov State Model

Optical spectroscopic tools are used to monitor protein folding/unfolding dynamics after a fast triggering such as the laser induced temperature jump. These techniques provide new opportunities for comparison between theory and simulations and atom-level understanding protein folding mechanism. However, the direct comparison still face two main challenges: a gap between folding relevant timescales (microseconds or above) and length of molecular dynamics simulations (typically tens to hundreds of nanoseconds), and difficulty in directly calculating spectroscopic observables from simulation configurations. Markov State Model (MSM) approach is one of the most powerful means which can increase simulations timescale up to microsecond or even millisecond. We address progress on modeling infrared and fluorescence spectroscopic signals of temperature jump induced unfolding dynamics for a few small proteins. The harmoniousness between experiment and theoretical can both improve our understanding of protein folding mechanisms and provide direct validation of those theoretical models.

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

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

DHPC strongly affects the structure and oligomerization propensity of Alzheimer's Aβ(1-40) peptide

Alzheimer's disease (AD) is thought to depend on the deleterious action of amyloid fibrils or oligomers derived from β-amyloid (Aβ) peptide. Out of various known Aβ alloforms, the 40-residue peptide Aβ(1-40) occurs at highest concentrations inside the brains of AD patients. Its aggregation properties critically depend on lipids, and it was thus proposed that lipids could play a major role in AD. To better understand their possible effects on the structure of Aβ and on the ability of this peptide to form potentially detrimental amyloid structures, we here analyze the interactions between Aβ(1-40) and 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC). DHPC has served, due to its controlled properties, as a major model system for studying general lipid properties. Here, we show that DHPC concentrations of 8 mM or higher exert dramatic effects on the conformation of soluble Aβ(1-40) peptide and induce the formation of β-sheet structure at high levels. By contrast, we find that DHPC concentrations well below the critical micelle concentration present no discernible effect on the conformation of soluble Aβ, although they substantially affect the peptide's oligomerization and fibrillation kinetics. These data imply that subtle lipid-peptide interactions suffice in controlling the overall aggregation properties and drastically accelerate, or delay, the fibrillation kinetics of Aβ peptide in near-physiological buffer solutions.

Structural classification of the amide I sites of a beta-hairpin with isotope label 2DIR spectroscopy

We present a theoretical study of the possibility to use isotope label two-dimensional infrared (2DIR) spectroscopy to obtain site specific structural information in trpzip2. This small beta-hairpin peptide was designed as a model system for studying protein folding in beta-sheet structures. In order to unravel the folding mechanism, the surroundings of local sites should be characterized, which in principle is possible by using 2DIR in combination with isotope labeling. This requires a classification that correlates local structures to two-dimensional spectra. To this end, we provide the first spectral simulation of the isotope label spectra of all the amide I sites in trpzip2. We find that the anti-diagonal width of the 2DIR peak associated with a labelled site is a good measure of solvent exposure and the key parameter to distinguish between solvent exposed and internal sites. The diagonal widths are not particularly sensitive to this, but they do reveal the presence of slowly interconverting turn structures.

Hydration dynamics in a partially denatured ensemble of the globular protein human alpha-lactalbumin investigated with molecular dynamics simulations

Atomistic molecular dynamics simulations are used to probe changes in the nature and subnanosecond dynamical behavior of solvation waters that accompany partial denaturation of the globular protein, human alpha-lactalbumin. A simulated ensemble of subcompact conformers, similar to the molten globule state of human alpha-lactalbumin, demonstrates a marginal increase in the amount of surface solvation relative to the native state. This increase is accompanied by subtle but distinct enhancement in surface water dynamics, less favorable protein-water interactions, and a marginal decrease in the anomalous behavior of solvation water dynamics. The extent of solvent influx is not proportional to the increased surface area, and the partially denatured conformers are less uniformly solvated compared to their native counterpart. The observed solvation in partially denatured conformers is lesser in extent compared to earlier experimental estimates in molten globule states, and is consistent with more recent descriptions based on nuclear magnetic relaxation dispersion studies.

Conserved folding pathways of alpha-lactalbumin and lysozyme revealed by kinetic CD, fluorescence, NMR, and interrupted refolding experiments

In this report, it is shown by a combination of stopped-flow CD, fluorescence, and time-resolved NMR studies that the Ca(2+)-induced refolding of bovine alpha-lactalbumin (BLA) at constant denaturant concentration (4 M urea) exhibits triple-exponential kinetics. In order to distinguish between parallel folding pathways and a strictly sequential formation of the native state, interrupted refolding experiments were conducted. We show here that the Ca(2+)-induced refolding of BLA involves parallel pathways and the transient formation of a folding intermediate on the millisecond timescale. Our data furthermore suggest that the two structurally homologous proteins BLA and hen egg white lysozyme share a common folding mechanism. We provide evidence that the guiding role of long-range interactions in the unfolded state of lysozyme in mediating intersubdomain interactions during folding is replaced in the case of BLA by the Ca(2+)-binding site. Time-resolved NMR spectroscopy, in combination with fast ion release from caged compounds, enables the measurement of complex protein folding kinetics at protein concentrations as low as 100 microM and the concomitant detection of conformational transitions with rate constants of up to 8 s(-1).

Contribution of rotational diffusion to pulsed field gradient diffusion measurements

NMR diffusion experiments employing pulsed field gradients are well established as sensitive probes of the displacement of individual nuclear spins in a sample. Conventionally such measurements are used as a measure of translational diffusion, but here we demonstrate that under certain conditions rotational motion will contribute very significantly to the experimental data. This situation occurs when at least one spatial dimension of the species under study exceeds the root mean square displacement associated with translational diffusion, and leads to anomalously large apparent diffusion coefficients when conventional analytical procedures are employed. We show that in such a situation the effective diffusion coefficient is a function of the duration of the diffusion delay used, and that this dependence provides a means of characterizing the dimensions of the species under investigation.

Protein folding and unfolding studied at atomic resolution by fast two-dimensional NMR spectroscopy

Atom-resolved real-time studies of kinetic processes in proteins have been hampered in the past by the lack of experimental techniques that yield sufficient temporal and atomic resolution. Here we present band-selective optimized flip-angle short transient (SOFAST) real-time 2D NMR spectroscopy, a method that allows simultaneous observation of reaction kinetics for a large number of nuclear sites along the polypeptide chain of a protein with an unprecedented time resolution of a few seconds. SOFAST real-time 2D NMR spectroscopy combines fast NMR data acquisition techniques with rapid sample mixing inside the NMR magnet to initiate the kinetic event. We demonstrate the use of SOFAST real-time 2D NMR to monitor the conformational transition of alpha-lactalbumin from a molten globular to the native state for a large number of amide sites along the polypeptide chain. The kinetic behavior observed for the disappearance of the molten globule and the appearance of the native state is monoexponential and uniform along the polypeptide chain. This observation confirms previous findings that a single transition state ensemble controls folding of alpha-lactalbumin from the molten globule to the native state. In a second application, the spontaneous unfolding of native ubiquitin under nondenaturing conditions is characterized by amide hydrogen exchange rate constants measured at high pH by using SOFAST real-time 2D NMR. Our data reveal that ubiquitin unfolds in a gradual manner with distinct unfolding regimes.

Diffusion NMR spectroscopy in supramolecular and combinatorial chemistry: an old parameter--new insights

Intermolecular interactions in solution play an important role in molecular recognition, which lies at the heart of supramolecular and combinatorial chemistry. Diffusion NMR spectroscopy gives information over such interactions and has become the method of choice for simultaneously measuring diffusion coefficients of multicomponent systems. The diffusion coefficient reflects the effective size and shape of a molecular species. Applications of this technique include the estimation of association constants and mapping the intermolecular interactions in multicomponent systems as well as investigating aggregation, ion pairing, encapsulation, and the size and structure of labile systems. Diffusion NMR spectroscopy can also be used to virtually separate mixtures and screen for specific ligands of different receptors, and may assist in finding lead compounds.

Diffusion NMR spectroscopy: folding and aggregation of domains in p53

Protein interactions and aggregation phenomena are probably amongst the most ubiquitous types of interactions in biological systems; they play a key role in many cellular processes. The ability to identify weak intermolecular interactions is a unique feature of NMR spectroscopy. In recent years, pulsed-field gradient NMR spectroscopy has become a convenient method to study molecular diffusion in solution. Since the diffusion coefficient of a certain molecule under given conditions correlates with its effective molecular weight, size, and shape, it is evident that diffusion can be used to map intermolecular interactions or aggregation events. Complex models can be derived from comparison of experimental diffusion data with those predicted by hydrodynamic simulations. In this review, we will give an introduction to pulsed-field gradient NMR spectroscopy and the hydrodynamic properties of proteins and peptides. Furthermore, we show examples for applying these techniques to a helical peptide and its hydrophobic oligomerization, as well as to the dimerization behavior and folding of p53.

Identification of partially disordered peptide intermediates through residue-specific NMR diffusion measurements

It is technically challenging to detect low-population partially disordered species that are in equilibrium with the folded and unfolded states. Residue-specific translational diffusion experiments measured by pulsed field gradient NMR have been used to detect the presence, and define the conformation, of such equilibrium intermediates. The experiment is demonstrated for equilibrium solutions of related triple helical peptides that model a small region of type I collagen with and without a mutation known to cause osteogenesis imperfecta. The data show that residue-specific diffusion coefficients of an interconverting trimer to monomer system can allow discrimination between a simple two-state model and more complex multistate models involving partially disordered intermediates.

Elucidation of the Protein Folding Landscape by NMR

NMR is one of the few experimental methods that can provide detailed insights into the structure and dynamics of unfolded and partly folded states of proteins. Mapping the protein folding landscape is of central importance to understanding the mechanism of protein folding. In addition, it is now recognized that many proteins are intrinsically unstructured in their functional states, while partly folded states of several cellular proteins have been implicated in amyloid disease. NMR is uniquely suited to characterize the structures present in the conformational ensemble and probe the dynamics of the polypeptide chain in unfolded and partially folded protein states.

Stochastic Liouville equation simulation of multidimensional vibrational line shapes of trialanine

The line shapes detected in coherent femtosecond vibrational spectroscopies contain direct signatures of peptide conformational fluctuations through their effect on vibrational frequencies and intermode couplings. These effects are simulated in trialanine using a Green's function solution of a stochastic Liouville equation constructed for four collective bath coordinates (two Ramachandran angles affecting the mode couplings and two diagonal energies). We find that fluctuations of the Ramachandran angles which hardly affect the linear absorption can be effectively probed by two-dimensional spectra. The signal generated at k(1)+k(2)-k(3) is particularly sensitive to such fluctuations.

Rapid Sample-Mixing Technique for Transient NMR and Photo-CIDNP Spectroscopy: Applications to Real-Time Protein Folding

We describe the development and application of a novel rapid sample-mixing technique for real-time NMR (nuclear magnetic resonance) spectroscopy. The apparatus consists of an insert inside a conventional NMR tube coupled to a rapid injection syringe outside the NMR magnet. Efficient and homogeneous mixing of solutions in the NMR tube is achieved with a dead time of tens of milliseconds, without modification of the NMR probe or additional hardware inside the magnet. Provision is made for the inclusion of an optical fiber to allow in situ laser irradiation of samples, for example to generate photo-CIDNP (chemically induced dynamic nuclear polarization). An NMR water suppression method has been implemented to allow experiments in H(2)O as well as in deuterated solvents. The performance of the device has been tested and optimized by a variety of methods, including sensitive detection of residual pH gradients and the use of NMR imaging to monitor the extent of mixing in real time. The potential utility of this device, in conjunction with the sensitivity and selectivity of photo-CIDNP, is demonstrated by experiments on the protein hen lysozyme. These measurements involve the direct detection of spectra during real-time refolding, and the use of CIDNP pulse labeling to study a partially unfolded state of the protein under equilibrium conditions. Magnetization transfer from this disordered state to the well-characterized native state provides evidence for the remarkable persistence of nativelike elements of structure under conditions in which the protein is partially denatured and aggregation prone.

Residue-specific real-time NMR diffusion experiments define the association states of proteins during folding

Characterizing the association states of proteins during folding is critical for understanding the nature of protein-folding intermediates and protein-folding pathways, protein aggregation, and disease-related aggregation. To study the association states of unfolded, folded, and intermediate species during protein folding, we have introduced a novel residue-specific real-time NMR diffusion experiment. This experiment, a combination of NMR real-time folding experiments and 3D heteronuclear pulsed field gradient NMR diffusion experiments (LED-HSQC), measures hydrodynamic properties, or molecular sizes, of kinetic species directly during the folding process. Application of the residue-specific real-time NMR diffusion experiments to characterize the folding of the collagen triple helix motif shows that this experiment can be used to determine the association states of unfolded, folded, and kinetic intermediates with transient lifetimes simultaneously. The ratio of the apparent translational diffusion coefficients of the unfolded to the folded form of the triple helix is 0.59, which correlates very well with a theoretical ratio for monomer to linear trimer. The apparent diffusion coefficients of the kinetic intermediates formed during triple helix folding indicate the formation of trimer-like associates which is consistent with previously published kinetic and relaxation data. The residue-specific time dependence of apparent diffusion coefficients of monomer and trimer peaks also illustrates the ability to use diffusion data to probe the directionality of triple helix formation. NMR diffusion experiments provide a new strategy for the investigation of protein-folding mechanisms, both to understand the role of kinetic intermediates and to determine the time-dependent aggregation processes in human diseases.

Signatures of beta-peptide unfolding in two-dimensional vibrational echo spectroscopy: a simulation study

An ensemble of exciton Hamiltonians for the amide-I band of the folded and unfolded states of a helical beta-heptapeptide is generated using a molecular dynamics (MD) simulation. The correlated fluctuations of its parameters and their signatures in two-dimensional (2D) vibrational echo spectroscopy are computed. This technique uses infrared pulse sequences to provide ultrafast snapshots of molecular structural fluctuations, in analogy with multidimensional NMR. The present study demonstrates that, by combining a method of calculating the vibrational Hamiltonian from MD snapshots and the nonlinear exciton equations (NEE), it may be possible to simulate realistic multidimensional IR spectra of chemically and biologically interesting systems.