Structural and dynamic characterization of partially folded states of apomyoglobin and implications for protein folding

Mechanism of Protein Aggregation Inhibition by Arginine: Blockage of Anionic Side Chains Favors Unproductive Encounter Complexes

Aggregation refers to the assembly of proteins into nonphysiological higher order structures. While amyloid has been studied extensively, much less is known about amorphous aggregation, a process that interferes with protein expression and storage. Free arginine (Arg+) is a widely used aggregation inhibitor, but its mechanism remains elusive. Focusing on myoglobin (Mb), we recently applied atomistic molecular dynamics (MD) simulations for gaining detailed insights into amorphous aggregation (Ng J. Phys. Chem. B 2021, 125, 13099). Building on that approach, the current work for the first time demonstrates that MD simulations can directly elucidate aggregation inhibition mechanisms. Comparative simulations with and without Arg+ reproduced the experimental finding that Arg+ significantly decreased the Mb aggregation propensity. Our data reveal that, without Arg+, protein-protein encounter complexes readily form salt bridges and hydrophobic contacts, culminating in firmly linked dimeric aggregation nuclei. Arg+ promotes the dissociation of encounter complexes. These "unproductive" encounter complexes are favored because Arg+ binding to D- and E- lowers the tendency of these anionic residues to form interprotein salt bridges. Side chain blockage is mediated largely by the guanidinium group of Arg+, which binds carboxylates through H-bond-reinforced ionic contacts. Our MD data revealed Arg+ self-association into a dynamic quasi-infinite network, but we found no evidence that this self-association is important for protein aggregation inhibition. Instead, aggregation inhibition by Arg+ is similar to that mediated by free guanidinium ions. The computational strategy used here should be suitable for the rational design of aggregation inhibitors with enhanced potency.

Aqueous alkaline phosphate facilitates the non-exchangeable deuteration of peptides and proteins

The incorporation of deuterium into peptides and proteins holds broad applications across various fields, such as drug development and structural characterization. Nevertheless, current methods for peptide/protein deuteration often target exchangeable labile sites or require harsh conditions for stable modification. In this study, we present a late-stage approach utilizing an alkaline phosphate solution to achieve deuteration of non-exchangeable backbone sites of peptides and proteins. The specific deuteration regions are identified through ultraviolet photodissociation (UVPD) and mass spectrometry analysis. This deuteration strategy demonstrates site and structure selectivity, with a notable affinity for labeling the α-helix regions of myoglobin. The deuterium method is particularly suitable for peptides and proteins that remain stable under high pH conditions.

Influence of Amino Acid Substitutions in ApoMb on Different Stages of Unfolding of Amyloids

To date, most research on amyloid aggregation has focused on describing the structure of amyloids and the kinetics of their formation, while the conformational stability of fibrils remains insufficiently explored. The aim of this work was to investigate the effect of amino acid substitutions on the stability of apomyoglobin (ApoMb) amyloids. A study of the amyloid unfolding of ApoMb and its six mutant variants by urea has been carried out. Changes in the structural features of aggregates during unfolding were recorded by far-UV CD and native electrophoresis. It was shown that during the initial stage of denaturation, amyloids' secondary structure partially unfolds. Then, the fibrils undergo dissociation and form intermediate aggregates weighing approximately 1 MDa, which at the last stage of unfolding decompose into 18 kDa monomeric unfolded molecules. The results of unfolding transitions suggest that the stability of the studied amyloids relative to the intermediate aggregates and of the latter relative to unfolded monomers is higher for ApoMb variants with substitutions that increase the hydrophobicity of the residues. The results presented provide a new insight into the mechanism of stabilization of protein aggregates and can serve as a base for further investigations of the amyloids' stability.

From Immunogenic Peptides to Intrinsically Disordered Proteins

It is hard to evaluate the role of individual mentors in the genesis of important ideas. In the case of our realization that proteins do not have to be stably folded to be functional, the influence of Richard Lerner and our collaborative work in the 1980s on the conformations of immunogenic peptides provided a base level of thinking about the nature of polypeptides in water solutions that led us to formulate and develop our ideas on the importance of intrinsic disorder in proteins. This review describes how the insights gained into the behavior of peptides led directly to the realization that proteins were not only capable of being functional while disordered, but also that disorder provided a distinct functional advantage in many important cellular processes.

Pre-Molten, Wet, and Dry Molten Globules en Route to the Functional State of Proteins

Transitions between the unfolded and native states of the ordered globular proteins are accompanied by the accumulation of several intermediates, such as pre-molten globules, wet molten globules, and dry molten globules. Structurally equivalent conformations can serve as native functional states of intrinsically disordered proteins. This overview captures the characteristics and importance of these molten globules in both structured and intrinsically disordered proteins. It also discusses examples of engineered molten globules. The formation of these intermediates under conditions of macromolecular crowding and their interactions with nanomaterials are also reviewed.

Transient On- and Off-Pathway Protein Folding Intermediate States Characterized with NMR Relaxation Dispersion

The earliest events in the folding of a protein are in general poorly understood. We used NMR R2 relaxation dispersion experiments to study transient local collapse events in the unfolded-state (U) conformational ensemble of apomyoglobin (apoMb). Local residual secondary structure (seen in regions corresponding to the A, D, E, and H helices of the folded protein) is largely unchanged over the pH range of 2.3-2.75, yet a significant pH-dependent increase in the conformational exchange contribution to the R2 relaxation rate (Rex) indicates that transient intramolecular contacts occur on a microsecond to millisecond time scale at pH 2.75. A comparison of 15N and 13CO relaxation dispersion data at pH 2.75 for residues in the A, B, G, and H regions, which participate in the earliest folding intermediates, indicates that chain collapse and secondary structure formation are rapid and concomitant. Increasingly stabilizing conditions (lower temperature, higher pH) result in the observation of a relaxation dispersion in the C, CD, and E regions of the protein, which are known to fold at later stages. Mutation of Trp14 in the A-helix region to Ala eliminates conformational exchange throughout the protein, and the mutation of hydrophobic residues in other regions results in the selective inhibition of conformational exchange in the B, G, or H regions. The R2 dispersion data for WT apoMb at pH 2.75 and 10 °C are best fit to a four-state model ABGH ⇆ AGH ⇆ U ⇆ ABCD that includes on-pathway (AGH and ABGH) and off-pathway (ABCD) transiently folded states, both of which are required to explain the behavior of the mutant proteins. The off-pathway intermediate is destabilized at higher temperatures. Our analysis provides insights into the earliest stages of apoMb folding where the collapsing polypeptide chain samples both productive and nonproductive states with stabilized secondary structure.

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

Conformational Dynamics of Intrinsically Disordered Proteins Regulate Biomolecular Condensate Chemistry

Motions in biomolecules are critical for biochemical reactions. In cells, many biochemical reactions are executed inside of biomolecular condensates formed by ultradynamic intrinsically disordered proteins. A deep understanding of the conformational dynamics of intrinsically disordered proteins in biomolecular condensates is therefore of utmost importance but is complicated by diverse obstacles. Here we review emerging data on the motions of intrinsically disordered proteins inside of liquidlike condensates. We discuss how liquid–liquid phase separation modulates internal motions across a wide range of time and length scales. We further highlight the importance of intermolecular interactions that not only drive liquid–liquid phase separation but appear as key determinants for changes in biomolecular motions and the aging of condensates in human diseases. The review provides a framework for future studies to reveal the conformational dynamics of intrinsically disordered proteins in the regulation of biomolecular condensate chemistry.

Exploring the folding energy landscapes of heme proteins using a hybrid AWSEM-heme model

Heme is an active center in many proteins. Here we explore computationally the role of heme in protein folding and protein structure. We model heme proteins using a hybrid model employing the AWSEM Hamiltonian, a coarse-grained forcefield for the protein chain along with AMBER, an all-atom forcefield for the heme. We carefully designed transferable force fields that model the interactions between the protein and the heme. The types of protein-ligand interactions in the hybrid model include thioester covalent bonds, coordinated covalent bonds, hydrogen bonds, and electrostatics. We explore the influence of different types of hemes (heme b and heme c) on folding and structure prediction. Including both types of heme improves the quality of protein structure predictions. The free energy landscape shows that both types of heme can act as nucleation sites for protein folding and stabilize the protein folded state. In binding the heme, coordinated covalent bonds and thioester covalent bonds for heme c drive the heme toward the native pocket. The electrostatics also facilitates the search for the binding site.

Mechanism of Thermal Protein Aggregation: Experiments and Molecular Dynamics Simulations on the High-Temperature Behavior of Myoglobin

Proteins that encounter unfavorable solvent conditions are prone to aggregation, a phenomenon that remains poorly understood. This work focuses on myoglobin (Mb) as a model protein. Upon heating, Mb produces amorphous aggregates. Thermal unfolding experiments at low concentration (where aggregation is negligible), along with centrifugation assays, imply that Mb aggregation proceeds via globally unfolded conformers. This contrasts studies on other proteins that emphasized the role of partially folded structures as aggregate precursors. Molecular dynamics (MD) simulations were performed to gain insights into the mechanism by which heat-unfolded Mb molecules associate with one another. A prerequisite for these simulations was the development of a method for generating monomeric starting structures. Periodic boundary condition artifacts necessitated the implementation of a partially immobilized water layer lining the walls of the simulation box. Aggregation simulations were performed at 370 K to track the assembly of monomeric Mb into pentameric species. Binding events were preceded by multiple unsuccessful encounters. Even after association, protein-protein contacts remained in flux. Binding was mediated by hydrophobic contacts, along with salt bridges that involved hydrophobically embedded Lys residues. Overall, this work illustrates that atomistic MD simulations are well suited for garnering insights into protein aggregation mechanisms.

Component of Cannabis, Cannabidiol, as a Possible Drug against the Cytotoxicity of Aβ(31-35) and Aβ(25-35) Peptides: An Investigation by Molecular Dynamics and Well-Tempered Metadynamics Simulations

In this work cannabidiol (CBD) was investigated as a possible drug against the cytotoxicity of Aβ(31-35) and Aβ(25-35) peptides with the help of atomistic molecular dynamics (MD) and well-tempered metadynamics simulations. Four interrelated mechanisms of possible actions of CBD are proposed from our computations. This implies that one mechanism can be a cause or/and a consequence of another. CBD is able to decrease the aggregation of peptides at certain concentrations of compounds in water. This particular action is more prominent for Aβ(25-35), since originally Aβ(31-35) did not exhibit aggregation properties in aqueous solutions. Interactions of CBD with the peptides affect secondary structures of the latter ones. Clusters of CBD are seen as possible adsorbents of Aβ(31-35) and Aβ(25-35) since peptides are tending to aggregate around them. And last but not least, CBD exhibits binding to MET35. All four mechanisms of actions can possibly inhibit the Aβ-cytotoxicity as discussed in this paper. Moreover, the amount of water also played a role in peptide clustering: with a growing concentration of peptides in water without a drug, the aggregation of both Aβ(31-35) and Aβ(25-35) increased. The number of hydrogen bonds between peptides and water was significantly higher for simulations with Aβ(25-35) at the higher concentration of peptides, while for Aβ(31-35) that difference was rather insignificant. The presence of CBD did not substantially affect the number of hydrogen bonds in the simulated systems.

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.

Discrete Hf18 Metal-oxo Cluster as a Heterogeneous Nanozyme for Site-Specific Proteolysis

The selective hydrolysis of proteins by non-enzymatic catalysis is difficult to achieve, yet it is crucial for applications in biotechnology and proteomics. Herein, we report that discrete hafnium metal-oxo cluster [Hf₁₈ O₁₀ (OH)₂₆ (SO₄ )₁₃ ⋅(H₂ O)₃₃ ] (Hf₁₈ ), which is centred by the same hexamer motif found in many MOFs, acts as a heterogeneous catalyst for the efficient hydrolysis of horse heart myoglobin (HHM) in low buffer concentrations. Among 154 amino acids present in the sequence of HHM, strictly selective cleavage at only 6 solvent accessible aspartate residues was observed. Mechanistic experiments suggest that the hydrolytic activity is likely derived from the actuation of Hf^IV Lewis acidic sites and the Brønsted acidic surface of Hf₁₈ . X-ray scattering and ESI-MS revealed that Hf₁₈ is completely insoluble in these conditions, confirming the HHM hydrolysis is caused by a heterogeneous reaction of the solid Hf₁₈ cluster, and not from smaller, soluble Hf species that could leach into solution.

Lessons Learned from 50 Years of Hemoglobin Research: Unstirred and Cell-Free Layers, Electrostatics, Baseball Gloves, and Molten Globules

Significance: Over the past 50 years, the mechanisms for O₂ storage and transport have been determined quantitatively on distance scales from millimeters to tenths of nanometers and timescales from seconds to picoseconds. Recent Advances: In this review, I have described four key conclusions from work done by my group and our close colleagues. (i) O₂ uptake by mammalian red cells is limited by diffusion through unstirred water layers adjacent to the cell surface and across cell-free layers adjacent to vessel walls. (ii) In most vertebrates, hemoglobins (Hbs) and myoglobins (Mbs), the distal histidine at the E7 helical position donates a strong hydrogen bond to bound O₂, which selectively enhances O₂ affinity, prevents carbon monoxide poisoning, and markedly slows autoxidation. (iii) O₂ binding to mammalian Hbs and Mbs occurs by migration of the ligand through a channel created by upward rotation of the His(E7) side chain, capture in the empty space of the distal pocket, and then coordination with the ferroprotoporphyrin IX (heme) iron atom. (iv) The assembly of Mbs and Hbs occurs by formation of molten globule intermediates, in which the N- and C-terminal helices have almost fully formed secondary structures, but the heme pockets are disordered and followed by high-affinity binding of heme. Critical Issues: These conclusions indicate that there are often compromises between O₂ transport function, holoprotein stability, and the efficiency of assembly. Future Directions: However, the biochemical mechanisms underlying these conclusions provide the framework for understanding globin evolution in greater detail and for engineering more efficient and stable globins.

Life in Phases: Intra- and Inter- Molecular Phase Transitions in Protein Solutions

Proteins, these evolutionarily-edited biological polymers, are able to undergo intramolecular and intermolecular phase transitions. Spontaneous intramolecular phase transitions define the folding of globular proteins, whereas binding-induced, intra- and inter- molecular phase transitions play a crucial role in the functionality of many intrinsically-disordered proteins. On the other hand, intermolecular phase transitions are the behind-the-scenes players in a diverse set of macrosystemic phenomena taking place in protein solutions, such as new phase nucleation in bulk, on the interface, and on the impurities, protein crystallization, protein aggregation, the formation of amyloid fibrils, and intermolecular liquid-liquid or liquid-gel phase transitions associated with the biogenesis of membraneless organelles in the cells. This review is dedicated to the systematic analysis of the phase behavior of protein molecules and their ensembles, and provides a description of the major physical principles governing intramolecular and intermolecular phase transitions in protein solutions.

Segmental Motions of Proteins under Non-native States Evaluated Using Quasielastic Neutron Scattering

Characterization of the dynamics of disordered polypeptide chains is required to elucidate the behavior of intrinsically disordered proteins and proteins under non-native states related to the folding process. Here we develop a method using quasielastic neutron scattering, combined with small-angle X-ray scattering and dynamic light scattering, to evaluate segmental motions of proteins as well as diffusion of the entire molecules and local side-chain motions. We apply this method to RNase A under the unfolded and molten-globule (MG) states. The diffusion coefficients arising from the segmental motions are evaluated and found to be different between the unfolded and MG states. The values obtained here are consistent with those obtained using the fluorescence-based techniques. These results demonstrate not only feasibility of this method but also usefulness to characterize the behavior of proteins under various disordered states.

Perspective: the essential role of NMR in the discovery and characterization of intrinsically disordered proteins

The 2019 ISMAR Prize recognized NMR studies of disordered proteins. Here we provide a highly personal perspective on the discovery of intrinsically disordered proteins and the development and application of NMR methods to characterize their conformational ensembles, dynamics, and interactions.

Synergistic Structural Information from Covalent Labeling and Hydrogen-Deuterium Exchange Mass Spectrometry for Protein-Ligand Interactions

Hydrogen-deuterium exchange (HDX) mass spectrometry (MS) and covalent labeling (CL) MS are typically considered to be complementary methods for protein structural analysis, because one probes the protein backbone, while the other probes side chains. For protein-ligand interactions, we demonstrate in this work that the two labeling techniques can provide synergistic structural information about protein-ligand binding when reagents like diethylpyrocarbonate (DEPC) are used for CL because of the differences in the reaction rates of DEPC and HDX. Using three model protein-ligand systems, we show that the slower time scale for DEPC labeling makes it only sensitive to changes in solvent accessibility and insensitive to changes in protein structural fluctuations, whereas HDX is sensitive to changes in both solvent accessibility and structural fluctuations. When used together, the two methods more clearly reveal binding sites and ligand-induced changes to structural fluctuations that are distant from the binding site, which is more comprehensive information than either technique alone can provide. We predict that these two methods will find widespread usage together for more deeply understanding protein-ligand interactions.

Misfolding of a Single Disulfide Bonded Globular Protein into a Low-Solubility Species Conformationally and Biophysically Distinct from the Native One

In practice and despite Anfinsen’s dogma, the refolding of recombinant multiple SS-bonded proteins is famously difficult because misfolded species with non-native SS-bonds appear upon the oxidization of their cysteine residues. On the other hand, single SS-bond proteins are thought to be simple to refold because their cysteines have only one SS-bond partner. Here, we report that dengue 4 envelope protein domain 3 (DEN4 ED3), a single SS-bonded protein can be irreversibly trapped into a misfolded species through the formation of its sole intramolecular SS-bond. The misfolded species had a much lower solubility than the native one at pHs higher than about 7, and circular dichroism measurements clearly indicated that its secondary structure content was different from the native species. Furthermore, the peaks in the Heteronuclear Single Quantum Correlation spectroscopy (HSQC) spectrum of DEN4 ED3 from the supernatant fraction were sharp and well dispersed, reflecting the beta-sheeted native structure, whereas the spectrum of the precipitated fraction showed broad signals clustered near its center suggesting no or little structure and a strong tendency to aggregate. The two species had distinct biophysical properties and could interconvert into each other only by cleaving and reforming the SS-bond, strongly suggesting that they are topologically different. This phenomenon can potentially happen with any single SS-bonded protein, and our observation emphasizes the need for assessing the conformation and biophysical properties of bacterially produced therapeutic proteins in addition to their chemical purities.

Charging and supercharging of proteins for mass spectrometry: recent insights into the mechanisms of electrospray ionization

Molecular dynamics simulations have uncovered mechanistic details of the protein ESI process under various experimental conditions.

Initial Protein Unfolding Events in Ubiquitin, Cytochrome c and Myoglobin Are Revealed with the Use of 213 nm UVPD Coupled to IM-MS

The initial stages of protein unfolding may reflect the stability of the entire fold and can also reveal which parts of a protein can be perturbed, without restructuring the rest. In this work, we couple UVPD with activated ion mobility mass spectrometry to measure how three model proteins start to unfold. Ubiquitin, cytochrome c and myoglobin ions produced via nESI from salty solutions are subjected to UV irradiation pre-mobility separation; experiments are conducted with a range of source conditions which alter the conformation of the precursor ion as shown by the drift time profiles. For all three proteins, the compact structures result in less fragmentation than more extended structures which emerge following progressive in-source activation. Cleavage sites are found to differ between conformational ensembles, for example, for the dominant charge state of cytochrome c [M + 7H]⁷⁺, cleavage at Phe10, Thr19 and Val20 was only observed in activating conditions whilst cleavage at Ala43 is dramatically enhanced. Mapping the photo-cleaved fragments onto crystallographic structures provides insight into the local structural changes that occur as protein unfolding progresses, which is coupled to global restructuring observed in the drift time profiles. Graphical Abstract.

Interrupted Pressure-Jump NMR Experiments Reveal Resonances of On-Pathway Protein Folding Intermediate

Previous pressure-jump NMR experiments on a pressure-sensitized double mutant of ubiquitin showed evidence that its folding occurs via two parallel, comparably efficient pathways: a single barrier and a two-barrier pathway. An interrupted folding NMR experiment is introduced, where for a brief period the pressure is dropped to atmospheric conditions (1 bar), followed by a jump back to high pressure for signal detection. Conventional, forward sampling of the indirect dimension during the low-pressure period correlates the ¹⁵N or ¹³C' chemical shifts of the unfolded protein at 1 bar to the ¹H frequencies of both the unfolded and folded proteins at high pressure. Remarkably, sampling the data of the same experiment in the reverse direction yields the frequencies of proteins present at the end of the low-pressure interval, which include unfolded, intermediate, and folded species. Although the folding intermediate ¹⁵N shifts differ strongly from natively folded protein, its ¹³C' chemical shifts, which are more sensitive probes for secondary structure, closely match those of the folded protein and indicate that the folding intermediate must have a structure that is quite similar to the native state.

Impact of cellular health conditions on the protein folding state in mammalian cells

By using in-cell NMR experiments, we have demonstrated that the protein folding state in cells is significantly influenced by the cellular health conditions. hAK1 was denatured in cells under stressful culture conditions, while it remained functional and properly folded in cells continuously supplied with a fresh medium.

sw ApoMb Amyloid Aggregation under Nondenaturing Conditions: The Role of Native Structure Stability

Investigation of the molecular mechanisms underlying amyloid-related human diseases attracts close attention. These diseases, the number of which currently is above 40, are characterized by formation of peptide or protein aggregates containing a cross-β structure. Most of the amyloidogenesis mechanisms described so far are based on experimental studies of aggregation of short peptides, intrinsically disordered proteins, or proteins under denaturing conditions, and studies of amyloid aggregate formations by structured globular proteins under conditions close to physiological ones are still in the initial stage. We investigated the effect of amino acid substitutions on propensity of the completely helical protein sperm whale apomyoglobin (sw ApoMb) for amyloid formation from its structured state in the absence of denaturing agents. Stability and aggregation of mutated sw ApoMb were studied using circular dichroism, Fourier transform infrared spectroscopy, x-ray diffraction, native electrophoresis, and electron microscopy techniques. Here, we demonstrate that stability of the protein native state determines both protein aggregation propensity and structural peculiarities of formed aggregates. Specifically, structurally stable mutants show low aggregation propensity and moderately destabilized sw ApoMb variants form amyloids, whereas their strongly destabilized mutants form both amyloids and nonamyloid aggregates.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

Monitoring Hydrogen Exchange During Protein Folding by Fast Pressure Jump NMR Spectroscopy

A method is introduced that permits direct observation of the rates at which backbone amide hydrogens become protected from solvent exchange after rapidly dropping the hydrostatic pressure inside the NMR sample cell from denaturing (2.5 kbar) to native (1 bar) conditions. The method is demonstrated for a pressure-sensitized ubiquitin variant that contains two Val to Ala mutations. Increased protection against hydrogen exchange with solvent is monitored as a function of time during the folding process. Results for 53 backbone amides show narrow clustering with protection occurring with a time constant of ca. 85 ms, but slower protection is observed around a reverse turn near the C-terminus of the protein. Remarkably, the native NMR spectrum returns with this slower time constant of ca. 150 ms, indicating that the almost fully folded protein retains molten globule characteristics with severe NMR line broadening until the final hydrogen bonds are formed. Prior to crossing the transition state barrier, hydrogen exchange protection factors are close to unity, but with slightly elevated values in the β₁-β₂ hairpin, previously shown to be already lowly populated in the urea-denatured state.

Characterization of intermediate state of myoglobin in the presence of PEG 10 under physiological conditions

The macromolecular crowding progressively has been gaining prominence in recent years as it acts as a sword with double-edge on protein stability and folding, i.e., showing assorted results of having both stabilizing and destabilizing effects. We studied the effects of different concentrations of polyethylene glycol (PEG-10) on structure and stability of myoglobin. The tertiary structure of myoglobin was found to be perturbed in the presence of polyethylene glycol, however there was insignificant change in the secondary structure. It was observed that polyethylene glycol induces molten globule state in myoglobin, where the intermediate state holds hydrophobic patches and larger hydrodynamic volume as compared to the native protein. In addition, isothermal titration calorimetry (ITC) showed strong binding between myoglobin and polyethylene glycol, at the physiological pH. We hypothesize that polyethylene glycol induces molten globule conformation in myoglobin by interacting with heme group of myoglobin. We caution that the binding of protein with crowder and other soft interactions need to be gravely well thought-out when studying macromolecular crowding. In our case, destabilizing protein-crowder interactions could compete and overcome the stabilizing exclusion volume effect.

Folding of apomyoglobin: Analysis of transient intermediate structure during refolding using quick hydrogen deuterium exchange and NMR

The structures of apomyoglobin folding intermediates have been widely analyzed using physical chemistry methods including fluorescence, circular dichroism, small angle X-ray scattering, NMR, mass spectrometry, and rapid mixing. So far, at least two intermediates (on sub-millisecond- and millisecond-scales) have been demonstrated for apomyoglobin folding. The combination of pH-pulse labeling and NMR is a useful tool for analyzing the kinetic intermediates at the atomic level. Its use has revealed that the latter-phase kinetic intermediate of apomyoglobin (6 ms) was composed of helices A, B, G and H, whereas the equilibrium intermediate, called the pH 4 molten-globule intermediate, was composed mainly of helices A, G and H. The improved strategy for the analysis of the kinetic intermediate was developed to include (1) the dimethyl sulfoxide method, (2) data processing with the various labeling times, and (3) a new in-house mixer. Particularly, the rapid mixing revealed that helices A and G were significantly more protected at the earlier stage (400 µs) of the intermediate (former-phase intermediate) than the other helices. Mutation studies, where each hydrophobic residue was replaced with an alanine in helices A, B, E, F, G and H, indicated that both non-native and native-like structures exist in the latter-phase folding intermediate. The N-terminal part of helix B is a weak point in the intermediate, and the docking of helix E residues to the core of the A, B, G and H helices was interrupted by a premature helix B, resulting in the accumulation of the intermediate composed of helices A, B, G and H. The prediction-based protein engineering produced important mutants: Helix F in a P88K/A90L/S92K/A94L mutant folded in the latter-phase intermediate, although helix F in the wild type does not fold even at the native state. Furthermore, in the L11G/W14G/A70L/G73W mutant, helix A did not fold but helix E did, which is similar to what was observed in the kinetic intermediate of apoleghemoglobin. Thus, this protein engineering resulted in a changed structure for the apomyoglobin folding intermediate.

The Neurite Outgrowth Inhibitory Nogo-A-Δ20 Region Is an Intrinsically Disordered Segment Harbouring Three Stretches with Helical Propensity

Functional recovery from central neurotrauma, such as spinal cord injury, is limited by myelin-associated inhibitory proteins. The most prominent example, Nogo-A, imposes an inhibitory cue for nerve fibre growth via two independent domains: Nogo-A-Δ20 (residues 544-725 of the rat Nogo-A sequence) and Nogo-66 (residues 1026-1091). Inhibitory signalling from these domains causes a collapse of the neuronal growth cone via individual receptor complexes, centred around sphingosine 1-phosphate receptor 2 (S1PR2) for Nogo-A-Δ20 and Nogo receptor 1 (NgR1) for Nogo-66. Whereas the helical conformation of Nogo-66 has been studied extensively, only little structural information is available for the Nogo-A-Δ20 region. We used nuclear magnetic resonance (NMR) spectroscopy to assess potential residual structural propensities of the intrinsically disordered Nogo-A-Δ20. Using triple resonance experiments, we were able to assign 94% of the non-proline backbone residues. While secondary structure analysis and relaxation measurements highlighted the intrinsically disordered character of Nogo-A-Δ20, three stretches comprising residues 561EAIQESL567, 639EAMNVALKALGT650, and 693SNYSEIAK700 form transient α-helical structures. Interestingly, 561EAIQESL567 is situated directly adjacent to one of the most conserved regions of Nogo-A-Δ20 that contains a binding motif for β1-integrin. Likewise, 639EAMNVALKALGT650 partially overlaps with the epitope recognized by 11C7, a Nogo-A-neutralizing antibody that promotes functional recovery from spinal cord injury. Diffusion measurements by pulse-field gradient NMR spectroscopy suggest concentration- and oxidation state-dependent dimerisation of Nogo-A-Δ20. Surprisingly, NMR and isothermal titration calorimetry (ITC) data could not validate previously shown binding of extracellular loops of S1PR2 to Nogo-A-Δ20.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

Probing the Time Scale of FPOP (Fast Photochemical Oxidation of Proteins): Radical Reactions Extend Over Tens of Milliseconds

Hydroxyl radical (⋅OH) labeling with mass spectrometry detection reports on protein conformations and interactions. Fast photochemical oxidation of proteins (FPOP) involves ⋅OH production via H2O2 photolysis by UV laser pulses inside a flow tube. The experiments are conducted in the presence of a scavenger (usually glutamine) that shortens the ⋅OH lifetime. The literature claims that FPOP takes place within 1 μs. This ultrafast time scale implies that FPOP should be immune to labeling-induced artifacts that may be encountered with other techniques. Surprisingly, the FPOP time scale has never been validated in direct kinetic measurements. Here we employ flash photolysis for probing oxidation processes under typical FPOP conditions. Bleaching of the reporter dye cyanine-5 (Cy5) served as readout of the time-dependent radical milieu. Surprisingly, Cy5 oxidation extends over tens of milliseconds. This time range is four orders of magnitude longer than expected from the FPOP literature. We demonstrate that the glutamine scavenger generates metastable secondary radicals in the FPOP solution, and that these radicals lengthen the time frame of Cy5 oxidation. Cy5 and similar dyes are widely used for monitoring the radical dose experienced by proteins in solution. The measured Cy5 kinetics thus strongly suggest that protein oxidation in FPOP extends over a much longer time window than previously thought (i.e., many milliseconds instead of one microsecond). The optical approach developed here should be suitable for assessing the performance of future FPOP-like techniques with improved temporal labeling characteristics. Graphical Abstract ᅟ.

Under-folded proteins: Conformational ensembles and their roles in protein folding, function, and pathogenesis

For decades, protein function was intimately linked to the presence of a unique, aperiodic crystal‐like structure in a functional protein. The two only places for conformational ensembles of under‐folded (or partially folded) protein forms in this picture were either the end points of the protein denaturation processes or transiently populated folding intermediates. Recent years witnessed dramatic change in this perception and conformational ensembles, which the under‐folded proteins are, have moved from the shadow. Accumulated to date data suggest that a protein can exist in at least three global forms–functional and folded, functional and intrinsically disordered (nonfolded), and nonfunctional and misfolded/aggregated. Under‐folded protein states are crucial for each of these forms, serving as important folding intermediates of ordered proteins, or as functional states of intrinsically disordered proteins (IDPs) and IDP regions (IDPRs), or as pathology triggers of misfolded proteins. Based on these observations, conformational ensembles of under‐folded proteins can be classified as transient (folding and misfolding intermediates) and permanent (IDPs and stable misfolded proteins). Permanently under‐folded proteins can further be split into intentionally designed (IDPs and IDPRs) and unintentionally designed (misfolded proteins). Although intrinsic flexibility, dynamics, and pliability are crucial for all under‐folded proteins, the different categories of under‐foldedness are differently encoded in protein amino acid sequences. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 870–887, 2013.

Effects of heme binding on myoglobin folding: Coarse grained molecular simulations

Many proteins contain cofactors, such as heme, ATP and metal ions. Binding of cofactors is not only essential for their biological functions, but also can reshape the intrinsic energy landscape of protein molecules and modulate the folding and stability. However, the molecular mechanism of cofactor coupled protein folding is not well understood. In this work, we study the cofactor coupled folding of myoglobin, which is a typical cofactor (heme) containing protein, by performing molecular dynamics simulations with a structure-based protein model developed based on the energy landscape theory. We showed that the heme binding increases the stability of the myoglobin. More importantly, the heme binding tends to increase the protein folding cooperativity, and switch the folding process from a “three-state” mechanism to a “two-state” mechanism. We also showed that the folding pathways of the myoglobin can be modulated by the heme binding. By performing comparative simulations, we revealed that the above effects of heme binding are resulted from the heme induced folding of F-helix, which is otherwise unstructured at apo state, and the heme mediated contacting interactions around the heme binding site. The simulation results are consistent with available experimental data, and provide insights into the molecular mechanism of the effects of cofactor binding on the protein folding and stability.

Picosecond to nanosecond dynamics provide a source of conformational entropy for protein folding

Myoglobin can be trapped in fully folded structures, partially folded molten globules, and unfolded states under stable equilibrium conditions.

Nuclear Magnetic Resonance Characterization of the Type III Secretion System Tip Chaperone Protein PcrG of Pseudomonas aeruginosa

Lung infection with Pseudomonas aeruginosa is the leading cause of death among cystic fibrosis patients. To initiate infection, P. aeruginosa assembles a protein nanomachine, the type III secretion system (T3SS), to inject bacterial proteins directly into target host cells. An important regulator of the P. aeruginosa T3SS is the chaperone protein PcrG, which forms a complex with the tip protein, PcrV. In addition to its role as a chaperone to the tip protein, PcrG also regulates protein secretion. PcrG homologues are also important in the T3SS of other pathogens such as Yersinia pestis, the causative agent of bubonic plague. The atomic structure of PcrG or any member of the family of tip protein chaperones is currently unknown. Here, we show by circular dichroism and nuclear magnetic resonance (NMR) spectroscopy that PcrG lacks a tertiary structure. However, it is not completely disordered but contains secondary structures dominated by two long α-helices from residue 16 to 41 and from residue 55 to 76. The helices of PcrG are partially formed, have similar backbone dynamics, and are flexible. NMR titrations show that the entire length of PcrG residues from position 9 to 76 is involved in binding to PcrV. PcrG adds to the growing list of partially folded or unstructured proteins with important roles in type III secretion.

High-Pressure EPR and Site-Directed Spin Labeling for Mapping Molecular Flexibility in Proteins

High hydrostatic pressure is a powerful probe of protein conformational flexibility. Pressurization reveals regions of elevated compressibility, and thus flexibility, within individual conformational states, but also shifts conformational equilibria such that "invisible" excited states become accessible for spectroscopic characterization. The central aim of this chapter is to describe recently developed instrumentation and methodologies that enable high-pressure site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) experiments on proteins and to demonstrate the information content of these experiments by highlighting specific recent applications. A brief introduction to the thermodynamics of proteins under pressure is presented first, followed by a discussion of the principles underlying SDSL-EPR detection of pressure effects in proteins, and the suitability of SDSL-EPR for this purpose in terms of timescale and ability to characterize conformational heterogeneity. Instrumentation and practical considerations for variable-pressure continuous wave EPR and pressure-resolved double electron-electron resonance (PR DEER) experiments are reviewed, and finally illustrations of data analysis using recent applications are presented. Although high-pressure SDSL-EPR is in its infancy, the recent applications presented highlight the considerable potential of the method to (1) identify compressible (flexible) regions in a folded protein; (2) determine thermodynamic parameters that relate conformational states in equilibrium; (3) populate and characterize excited states of proteins undetected at atmospheric pressure; (4) reveal the structural heterogeneity of conformational ensembles and provide distance constraints on the global structure of pressure-populated states with PR DEER.

The LcrG Tip Chaperone Protein of the Yersinia pestis Type III Secretion System Is Partially Folded

The type III secretion system (T3SS) is essential in the pathogenesis of Yersinia pestis, the causative agent of plague. A small protein, LcrG, functions as a chaperone to the tip protein LcrV, and the LcrG-LcrV interaction is important in regulating protein secretion through the T3SS. The atomic structure of the LcrG family is currently unknown. However, because of its predicted helical propensity, many have suggested that the LcrG family forms a coiled-coil structure. Here, we show by NMR and CD spectroscopy that LcrG lacks a tertiary structure and it consists of three partially folded α-helices spanning residues 7-38, 41-46, and 58-73. NMR titrations of LcrG with LcrV show that the entire length of a truncated LcrG (residues 7-73) is involved in binding to LcrV. However, there is regional variation in how LcrG binds to LcrV. The C-terminal region of a truncated LcrG (residues 52-73) shows tight binding interaction with LcrV while the N-terminal region (residues 7-51) shows weaker interaction with LcrV. This suggests that there are at least two binding events when LcrG binds to LcrV. Biological assays and mutagenesis indicate that the C-terminal region of LcrG (residues 52-73) is important in blocking protein secretion through the T3SS. Our results reveal structural and mechanistic insights into the atomic conformation of LcrG and how it binds to LcrV.

Evidence for a Shared Mechanism in the Formation of Urea-Induced Kinetic and Equilibrium Intermediates of Horse Apomyoglobin from Ultrarapid Mixing Experiments

In this study, the equivalence of the kinetic mechanisms of the formation of urea-induced kinetic folding intermediates and non-native equilibrium states was investigated in apomyoglobin. Despite having similar structural properties, equilibrium and kinetic intermediates accumulate under different conditions and via different mechanisms, and it remains unknown whether their formation involves shared or distinct kinetic mechanisms. To investigate the potential mechanisms of formation, the refolding and unfolding kinetics of horse apomyoglobin were measured by continuous- and stopped-flow fluorescence over a time range from approximately 100 μs to 10 s, along with equilibrium unfolding transitions, as a function of urea concentration at pH 6.0 and 8°C. The formation of a kinetic intermediate was observed over a wider range of urea concentrations (0–2.2 M) than the formation of the native state (0–1.6 M). Additionally, the kinetic intermediate remained populated as the predominant equilibrium state under conditions where the native and unfolded states were unstable (at ~0.7–2 M urea). A continuous shift from the kinetic to the equilibrium intermediate was observed as urea concentrations increased from 0 M to ~2 M, which indicates that these states share a common kinetic folding mechanism. This finding supports the conclusion that these intermediates are equivalent. Our results in turn suggest that the regions of the protein that resist denaturant perturbations form during the earlier stages of folding, which further supports the structural equivalence of transient and equilibrium intermediates. An additional folding intermediate accumulated within ~140 μs of refolding and an unfolding intermediate accumulated in <1 ms of unfolding. Finally, by using quantitative modeling, we showed that a five-state sequential scheme appropriately describes the folding mechanism of horse apomyoglobin.

Apoglobin Stability Is the Major Factor Governing both Cell-free and in Vivo Expression of Holomyoglobin

Expression levels in animal muscle tissues and in Escherichia coli vary widely for naturally occurring mammalian myoglobins (Mb). To explore this variation, we developed an in vitro transcription and wheat germ extract-based translation assay to examine quantitatively the factors that govern expression of holoMb. We constructed a library of naturally occurring Mbs from two terrestrial and four deep-diving aquatic mammals and three distal histidine mutants designed to enhance apoglobin stability but decrease hemin affinity. A strong linear correlation is observed between cell-free expression levels of holo-metMb variants and their corresponding apoglobin stabilities, which were measured independently by guanidine HCl-induced unfolding titrations using purified proteins. In contrast, there is little dependence of expression on hemin affinity. Our results confirm quantitatively that deep diving mammals have highly stable Mbs that express to higher levels in animal myocytes, E. coli, and the wheat germ cell-free system than Mbs from terrestrial mammals. Our theoretical analyses show that the rate of aggregation of unfolded apoMb is very large, and as a result, the key factor for high level expression of holoMb, and presumably other heme proteins, is an ultra high fraction of folded, native apoglobin that is capable of rapidly binding hemin. This fraction is determined by the overall equilibrium folding constant and not hemin affinity. These results also demonstrate that the cell-free transcription/translation system can be used as a high throughput platform to screen for apoglobin stability without the need to generate large amounts of protein for in vitro unfolding measurements.

Spider wrapping silk fibre architecture arising from its modular soluble protein precursor

Spiders store spidroins in their silk glands as high concentration aqueous solutions, spinning these dopes into fibres with outstanding mechanical properties. Aciniform (or wrapping) silk is the toughest spider silk and is devoid of the short amino acid sequence motifs characteristic of the other spidroins. Using solution-state NMR spectroscopy, we demonstrate that the 200 amino acid Argiope trifasciata AcSp1 repeat unit contrasts with previously characterized spidroins, adopting a globular 5-helix bundle flanked by intrinsically disordered N- and C-terminal tails. Split-intein-mediated segmental NMR-active isotope-enrichment allowed unambiguous demonstration of modular and malleable "beads-on-a-string" concatemeric behaviour. Concatemers form fibres upon manual drawing with silk-like morphology and mechanical properties, alongside secondary structuring and orientation consistent with native AcSp1 fibres. AcSp1 structural stability varies locally, with the fifth helix denaturing most readily. The structural transition of aciniform spidroin from a mostly α-helical dope to a mixed α-helix/β-sheet-containing fibre can be directly related to spidroin architecture and stability.

Kinetic intermediates of holo- and apo-myoglobin studied using HDX-TIMS-MS and molecular dynamic simulations

In the present work, the kinetic intermediates of holo- and apo-myoglobin were studied by correlating the ion-neutral collision cross section and time resolved H/D back exchange rate simultaneously in a trapped ion mobility spectrometer coupled to a mass spectrometer (HDX-TIMS-MS). The high mobility resolution of the TIMS cell permitted the observation of multiple IMS bands and complementary molecular dynamics simulations resulted in the assignment of candidate structures for each experimental condition studied (e.g., holo [M + 8H](+8)-[M + 9H](+9) and apo [M + 9H](+9)-[M + 19H](+19)). Inspection of the kinetic intermediates suggests that the tertiary structure of apomyoglobin unfolds quickly upon the loss of the Fe protoporphyrin IX that stabilizes the interactions between the A, G, and H helices. In the absence of the porphyrin heme, the apomyoglobin unfolds to Xn kinetic intermediates that vary in the extent of unfolding as a result of the observed charge state.

Structural characterization of holo- and apo-myoglobin in the gas phase by ultraviolet photodissociation mass spectrometry

Ultraviolet photodissociation (UVPD) mass spectrometry is employed to investigate the structure of holo-myoglobin as well as its apo form transferred to the gas phase by native electrospray. UVPD provided insight into the stability of native structural elements of holo-myoglobin. The fragmentation yields from UVPD showed the greatest overall correlation with B-factors generated from the crystal structure of apo-myoglobin, particularly for the more disordered loop regions. Solvent accessibility measurements also showed some correlation with the UVPD fragmentation of holo-myoglobin. Comparison of UVPD of holo- and apo-myoglobin revealed similarities in fragmentation yields, particularly for the lower charge states (8 and 9+). Both holo- and apo-myoglobin exhibited low fragmentation yields for the AGH helical core, whereas regions known to interact with the heme show suppressed fragmentation for holo-myoglobin. The fragment yields from HCD showed the lowest correlation with B-factor values and rather reflected preferential charge-directed backbone cleavages.

Structural characterization of MG and pre-MG states of proteins by MD simulations, NMR, and other techniques

Almost all proteins fold via a number of partially structured intermediates such as molten globule (MG) and pre-molten globule states. Understanding the structure of these intermediates at atomic level is often a challenge, as these states are observed under extreme conditions of pH, temperature, and chemical denaturants. Furthermore, several other processes such as chemical modification, site-directed mutagenesis (or point mutation), and cleavage of covalent bond of natural proteins often lead to MG like partially unfolded conformation. However, the dynamic nature of proteins in these states makes them unsuitable for most structure determination at atomic level. Intermediate states studied so far have been characterized mostly by circular dichroism, fluorescence, viscosity, dynamic light scattering measurements, dye binding, infrared techniques, molecular dynamics simulations, etc. There is a limited amount of structural data available on these intermediate states by nuclear magnetic resonance (NMR) and hence there is a need to characterize these states at the molecular level. In this review, we present characterization of equilibrium intermediates by biophysical techniques with special reference to NMR.