Collapse and search dynamics of apomyoglobin folding revealed by submillisecond observations of alpha-helical content and compactness

Non-Native Structures of Apomyoglobin and Apoleghemoglobin in Folding Intermediates Related to the Protein Misfolding

Protein folding is essential for a polypeptide chain to acquire its proper structure and function. Globins are a superfamily of ubiquitous heme-binding α-helical proteins whose function is principally to regulate oxygen homoeostasis. In this review, we explore the hierarchical helical formation in the globin proteins apomyoglobin and leghemoglobin, and we discuss the existence of non-native and misfolded structures occurring during the course of folding to its native state. This review summarizes the research aimed at characterizing and comparing the equilibrium and kinetic intermediates, as well as delineating the complete folding pathway at a molecular level, in order to answer the following questions: "What is the mechanism of misfolding via a folding intermediate? Does the non-native structure stabilize the contemporary intermediate structure? Does the non-native structure induce slower folding?" The role of the non-native structures in the folding intermediate related to misfolding is also discussed.

Genetically encoded non-canonical amino acids reveal asynchronous dark reversion of chromophore, backbone and side-chains in EL222

Photoreceptors containing the light-oxygen-voltage (LOV) domain elicit biological responses upon excitation of their flavin mononucleotide (FMN) chromophore by blue light. The mechanism and kinetics of dark-state recovery are not well understood. Here we incorporated the non-canonical amino acid p-cyanophenylalanine (CNF) by genetic code expansion technology at forty-five positions of the bacterial transcription factor EL222. Screening of light-induced changes in infrared (IR) absorption frequency, electric field and hydration of the nitrile groups identified residues CNF31 and CNF35 as reporters of monomer/oligomer and caged/decaged equilibria, respectively. Time-resolved multi-probe UV/Visible and IR spectroscopy experiments of the lit-to-dark transition revealed four dynamical events. Predominantly, rearrangements around the A'α helix interface (CNF31 and CNF35) precede FMN-cysteinyl adduct scission, folding of α-helices (amide bands), and relaxation of residue CNF151. This study illustrates the importance of characterizing all parts of a protein and suggests a key role for the N-terminal A'α extension of the LOV domain in controlling EL222 photocycle length. This article is protected by copyright. All rights reserved.

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.

Advances in Mixer Design and Detection Methods for Kinetics Studies of Macromolecular Folding and Binding on the Microsecond Time Scale

Many important biological processes such as protein folding and ligand binding are too fast to be fully resolved using conventional stopped-flow techniques. Although advances in mixer design and detection methods have provided access to the microsecond time regime, there is room for improvement in terms of temporal resolution and sensitivity. To address this need, we developed a continuous-flow mixing instrument with a dead time of 12 to 27 µs (depending on solution viscosity) and enhanced sensitivity, sufficient for monitoring tryptophan or tyrosine fluorescence changes at fluorophore concentrations as low as 1 µM. Relying on commercially available laser microfabrication services, we obtained an integrated mixer/flow-cell assembly on a quartz chip, based on a cross-channel configuration with channel dimensions and geometry designed to minimize backpressure. By gradually increasing the width of the observation channel downstream from the mixing region, we are able to monitor a reaction progress time window ranging from ~10 µs out to ~3 ms. By combining a solid-state UV laser with a Galvano-mirror scanning strategy, we achieved highly efficient and uniform fluorescence excitation along the flow channel. Examples of applications, including refolding of acid-denatured cytochrome c triggered by a pH jump and binding of a peptide ligand to a PDZ domain, demonstrate the capability of the technique to resolve fluorescence changes down to the 10 µs time regime on modest amounts of reagents.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

Resolving Site-Specific Heterogeneity of the Unfolded State under Folding Conditions

Understanding the properties of the unfolded state under folding conditions is of fundamental importance for gaining mechanistic insight into folding as well as misfolding reactions. Toward achieving this objective, the folding reaction of a small protein, monellin, has been resolved structurally and temporally, with the use of the multisite time-resolved FRET methodology. The present study establishes that the initial polypeptide chain collapse is not only heterogeneous but also structurally asymmetric and nonuniform. The population-averaged size for the segments spanning parts of the β-sheet decreases much more than that for the α-helix. Multisite measurements enabled specific and nonspecific components of the initial chain collapse to be discerned. The expanded and compact intermediate subensembles have the properties of a nonspecifically collapsed (hence, random-coil-like) and specifically collapsed (hence, globular) polymer, respectively. During subsequent folding, both the subensembles underwent contraction to varying extents at the four monitored segments, which was close to gradual in nature. The expanded intermediate subensemble exhibited an additional very slow contraction, suggestive of the presence of non-native interactions that result in a higher effective viscosity slowing down intrachain motions under folding conditions.

Further thermo-stabilization of thermophilic rhodopsin from Thermus thermophilus JL-18 through engineering in extramembrane regions

It is known that a hyperthermostable protein tolerable at temperatures over 100°C can be designed from a soluble globular protein by introducing mutations. To expand the applicability of this technology to membrane proteins, here we report a further thermo-stabilization of the thermophilic rhodopsin from Thermus thermophilus JL-18 as a model membrane protein. Ten single mutations in the extramembrane regions were designed based on a computational prediction of folding free-energy differences upon mutation. Experimental characterizations using the UV-visible spectroscopy and the differential scanning calorimetry revealed that four of ten mutations were thermo-stabilizing: V79K, T114D, A115P, and A116E. The mutation-structure relationship of the TR constructs was analyzed using molecular dynamics simulations at 300 K and at 1800 K that aimed simulating structures in the native and in the random-coil states, respectively. The native-state simulation exhibited an ion-pair formation of the stabilizing V79K mutant as it was designed, and suggested a mutation-induced structural change of the most stabilizing T114D mutant. On the other hand, the random-coil-state simulation revealed a higher structural fluctuation of the destabilizing mutant S8D when compared to the wild type, suggesting that the higher entropy in the random-coil state deteriorated the thermal stability. The present thermo-stabilization design in the extramembrane regions based on the free-energy calculation and the subsequent evaluation by the molecular dynamics may be useful to improve the production of membrane proteins for structural studies.

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.

The Molten Globule, and Two-State vs. Non-Two-State Folding of Globular Proteins

From experimental studies of protein folding, it is now clear that there are two types of folding behavior, i.e., two-state folding and non-two-state folding, and understanding the relationships between these apparently different folding behaviors is essential for fully elucidating the molecular mechanisms of protein folding. This article describes how the presence of the two types of folding behavior has been confirmed experimentally, and discusses the relationships between the two-state and the non-two-state folding reactions, on the basis of available data on the correlations of the folding rate constant with various structure-based properties, which are determined primarily by the backbone topology of proteins. Finally, a two-stage hierarchical model is proposed as a general mechanism of protein folding. In this model, protein folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the molten globule state as a folding intermediate. The two-state folding is thus merely a simplified version of the hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

Protein tertiary structure and the myoglobin phase diagram

We develop an effective theory approach to investigate the phase properties of globular proteins. Instead of interactions between individual atoms or localized interaction centers, the approach builds directly on the tertiary structure of a protein. As an example we construct the phase diagram of (apo)myoglobin with temperature (T) and acidity (pH) as the thermodynamical variables. We describe how myoglobin unfolds from the native folded state to a random coil when temperature and acidity increase. We confirm the presence of two molten globule folding intermediates, and we predict an abrupt transition between the two when acidity changes. When temperature further increases we find that the abrupt transition line between the two molten globule states terminates at a tricritical point, where the helical structures fade away. Our results also suggest that the ligand entry and exit is driven by large scale collective motions that destabilize the myoglobin F-helix.

A microfluidic flow-focusing device for low sample consumption serial synchrotron crystallography experiments in liquid flow

Serial synchrotron crystallography allows low X-ray dose, room-temperature crystal structures of proteins to be determined from a population of microcrystals. Protein production and crystallization is a non-trivial procedure and it is essential to have X-ray-compatible sample environments that keep sample consumption low and the crystals in their native environment. This article presents a fast and optimized manufacturing route to metal-polyimide microfluidic flow-focusing devices which allow for the collection of X-ray diffraction data in flow. The flow-focusing conditions allow for sample consumption to be significantly decreased, while also opening up the possibility of more complex experiments such as rapid mixing for time-resolved serial crystallography. This high-repetition-rate experiment allows for full datasets to be obtained quickly (∼1 h) from crystal slurries in liquid flow. The X-ray compatible microfluidic chips are easily manufacturable, reliable and durable and require sample-flow rates on the order of only 30 µl h^-1.

PFDB: A standardized protein folding database with temperature correction

We constructed a standardized protein folding kinetics database (PFDB) in which the logarithmic rate constants of all listed proteins are calculated at the standard temperature (25 °C). A temperature correction based on the Eyring–Kramers equation was introduced for proteins whose folding kinetics were originally measured at temperatures other than 25 °C. We verified the temperature correction by comparing the logarithmic rate constants predicted and experimentally observed at 25 °C for 14 different proteins, and the results demonstrated improvement of the quality of the database. PFDB consists of 141 (89 two-state and 52 non-two-state) single-domain globular proteins, which has the largest number among the currently available databases of protein folding kinetics. PFDB is thus intended to be used as a standard for developing and testing future predictive and theoretical studies of protein folding. PFDB can be accessed from the following link: http://lee.kias.re.kr/~bala/PFDB.

Complex Folding Landscape of Apomyoglobin at Acidic pH Revealed by Ultrafast Kinetic Analysis of Core Mutants

Under mildly acidic conditions (pH 4-4.5) apomyoglobin (apoMb) adopts a partially structured equilibrium state ( M-state) that structurally resembles a kinetic intermediate encountered at a late stage of folding to the native structure at neutral pH. We have previously reported that the M-state is formed rapidly (<1 ms) via a multistate process and thus offers a unique opportunity for exploring early stages of folding by both experimental and computational techniques. In order to gain structural insight into intermediates and barriers at the residue level, we studied the folding/unfolding kinetics of 12 apoMb mutants at pH 4.2 using fluorescence-detected ultrafast mixing techniques. Global analysis of the submillisecond folding/unfolding kinetics vs urea concentration for each variant, based on a sequential four-state mechanism ( U ⇔ I ⇔ L ⇔ M), allowed us to determine elementary rate constants and their dependence on urea concentration for most transitions. Comparison of the free energy diagrams constructed from the kinetic data of the mutants with that of wild-type apoMb yielded quantitative information on the effects of mutations on the free energy (ΔΔ G) of both intermediates and the first two kinetic barriers encountered during folding. Truncation of conserved aliphatic side chains on helices A, G, and H gives rise to a stepwise increase in ΔΔ G as the protein advances from U toward M, consistent with progressive stabilization of native-like contacts within the primary core of apoMb. Helix-helix contacts in the primary core contribute little to the first folding barrier ( U ⇔ I) and thus are not required for folding initiation but are critical for the stability of the late intermediate, L, and the M-state. Alanine substitution of hydrophobic residues at more peripheral helix-helix contact sites of the native structure, which are still absent or unstable in the M-state, shows both positive (destabilizing) and negative (stabilizing) ΔΔ G, indicating that non-native contacts are formed initially and weakened or lost as a result of subsequent structural rearrangement steps.

Cyclic olefin copolymer as an X-ray compatible material for microfluidic devices

The combination of microfluidics and X-ray methods attracts a lot of attention from researchers as it brings together the high controllability of microfluidic sample environments and the small length scales probed by X-rays. In particular, the fields of biophysics and biology have benefited enormously from such approaches. We introduce a straightforward fabrication method for X-ray compatible microfluidic devices made solely from cyclic olefin copolymers. We benchmark the performance of the devices against other devices including more commonly used Kapton windows and obtain data of equal quality using small angle X-ray scattering. An advantage of the devices presented here is that no gluing between interfaces is necessary, rendering the production very reliable. As a biophysical application, we investigate the early time points of the assembly of vimentin intermediate filament proteins into higher-order structures. This weakly scattering protein system leads to high quality data in the new devices, thus opening up the way for numerous future applications.

How Does Your Protein Fold? Elucidating the Apomyoglobin Folding Pathway

Although each type of protein fold and in some cases individual proteins within a fold classification can have very different mechanisms of folding, the underlying biophysical and biochemical principles that operate to cause a linear polypeptide chain to fold into a globular structure must be the same. In an aqueous solution, the protein takes up the thermodynamically most stable structure, but the pathway along which the polypeptide proceeds in order to reach that structure is a function of the amino acid sequence, which must be the final determining factor, not only in shaping the final folded structure, but in dictating the folding pathway. A number of groups have focused on a single protein or group of proteins, to determine in detail the factors that influence the rate and mechanism of folding in a defined system, with the hope that hypothesis-driven experiments can elucidate the underlying principles governing the folding process. Our research group has focused on the folding of the globin family of proteins, and in particular on the monomeric protein apomyoglobin. Apomyoglobin (apoMb) folds relatively slowly (∼2 s) via an ensemble of obligatory intermediates that form rapidly after the initiation of folding. The folding pathway can be dissected using rapid-mixing techniques, which can probe processes in the millisecond time range. Stopped-flow measurements detected by circular dichroism (CD) or fluorescence spectroscopy give information on the rates of folding events. Quench-flow experiments utilize the differential rates of hydrogen-deuterium exchange of amide protons protected in parts of the structure that are folded early; protection of amides can be detected by mass spectrometry or proton nuclear magnetic resonance spectroscopy (NMR). In addition, apoMb forms an intermediate at equilibrium at pH ∼ 4, which is sufficiently stable for it to be structurally characterized by solution methods such as CD, fluorescence and NMR spectroscopies, and the conformational ensembles formed in the presence of denaturing agents and low pH can be characterized as models for the unfolded states of the protein. Newer NMR techniques such as measurement of residual dipolar couplings in the various partly folded states, and relaxation dispersion measurements to probe invisible states present at low concentrations, have contributed to providing a detailed picture of the apomyoglobin folding pathway. The research summarized in this Account was aimed at characterizing and comparing the equilibrium and kinetic intermediates both structurally and dynamically, as well as delineating the complete folding pathway at a residue-specific level, in order to answer the question: "What is it about the amino acid sequence that causes each molecule in the unfolded protein ensemble to start folding, and, once started, to proceed towards the formation of the correctly folded three-dimensional structure?"

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.

Thermal unfolding of myoglobin in the Landau-Ginzburg-Wilson approach

The Landau-Ginzburg-Wilson paradigm is applied to model the low-temperature crystallographic Cα backbone structure of sperm whale myoglobin. The Glauber protocol is employed to simulate its response to an increase in ambient temperature. The myoglobin is found to unfold from its native state by a succession of α-helical intermediates, fully in line with the observed folding and unfolding patterns in denaturation experiments. In particular, a molten globule intermediate is identified with experimentally correct attributes. A detailed, experimentally testable contact map is constructed to characterize the specifics of the unfolding pathway, including the formation of long-range interactions. The results reveal how the unfolding process of a protein is driven by the interplay between, and a successive melting of, its modular secondary structure components.

Recent advances in X-ray compatible microfluidics for applications in soft materials and life sciences

The increasingly narrow and brilliant beams at X-ray facilities reduce the requirements for both sample volume and data acquisition time. This creates new possibilities for the types and number of sample conditions that can be examined but simultaneously increases the demands in terms of sample preparation. Microfluidic-based sample preparation techniques have emerged as elegant alternatives that can be integrated directly into the experimental X-ray setup remedying several shortcomings of more traditional methods. We review the use of microfluidic devices in conjunction with X-ray measurements at synchrotron facilities in the context of 1) mapping large parameter spaces, 2) performing time resolved studies of mixing-induced kinetics, and 3) manipulating/processing samples in ways which are more demanding or not accessible on the macroscale. The review covers the past 15 years and focuses on applications where synchrotron data collection is performed in situ, i.e. directly on the microfluidic platform or on a sample jet from the microfluidic device. Considerations such as the choice of materials and microfluidic designs are addressed. The combination of microfluidic devices and measurements at large scale X-ray facilities is still emerging and far from mature, but it definitely offers an exciting array of new possibilities.

Direct monitoring of calcium-triggered phase transitions in cubosomes using small-angle X-ray scattering combined with microfluidics

This article introduces a simple microfluidic device that can be combined with synchrotron small-angle X-ray scattering (SAXS) for monitoring dynamic structural transitions. The microfluidic device is a thiol–ene-based system equipped with 125 µm-thick polystyrene windows, which are suitable for X-ray experiments. The device was prepared by soft lithography using elastomeric molds followed by a simple UV-initiated curing step to polymerize the chip material and simultaneously seal the device with the polystyrene windows. The microfluidic device was successfully used to explore the dynamics of the structural transitions of phytantriol/dioleoylphosphatidylglycerol-based cubosomes on exposure to a buffer containing calcium ions. The resulting SAXS data were resolved in the time frame between 0.5 and 5.5 s, and a calcium-triggered structural transition from an internal inverted-type cubic phase of symmetryIm3mto an internal inverted-type cubic phase of symmetryPn3mwas detected. The combination of microfluidics with X-ray techniques opens the door to the investigation of early dynamic structural transitions, which is not possible with conventional techniques such as glass flow cells. The combination of microfluidics with X-ray techniques can be used for investigating protein unfolding, for monitoring the formation of nanoparticles in real time, and for other biomedical and pharmaceutical investigations.

Statistical Mechanical Model for pH-Induced Protein Folding: Application to Apomyoglobin

Despite the major role of pH in protein folding and stability, a quantitative understanding of the pH-induced protein folding mechanism remains elusive. Two conventional models, the Monod-Wyman-Changeux and Linderstrøm-Lang smeared charge models, respectively, have been used to analyze the formation/disruption of specific native structures and fluctuating non-native states. However, there are only a few models that can represent the overall kinetic events of folding/unfolding independent of the properties of relevant molecular species, which has hampered the efforts to systematically analyze pH-induced folding. Here, we constructed a statistical mechanical model that incorporates the protonation mechanism of conventional models along with a combined manual search and least-squares fitting procedure, which was used to investigate the folding of horse apomyoglobin over a wide pH range (2.2-6.7), with a time window ranging from ∼40 μs to ∼100 s, using continuous-/stopped-flow fluorescence at 8 °C. Quantitative analysis assuming a five-state sequential scheme indicated that (1) pH-induced folding/unfolding is represented by both specific binding and Coulombic interactions; (2) kinetic folding/unfolding intermediates share kinetic mechanisms with the equilibrium intermediate, indicating their equivalence; and (3) native-like properties are acquired successively during folding by intermediates and in transition states. This model could also be applied to a variety of association/dissociation processes.

Sulfmyoglobin Conformational Change: A Role in the Decrease of Oxy-Myoglobin Functionality

This work is focused at understanding the interaction of H2S with Myoglobin (Mb), in particular the Sulfmyoglobin (SMb) product, whose physiological role is controversial and not well understood. The scattering curves, Guinier, Kratky, Porod and P(r) plots were analyzed for oxy-Mb and oxy-Hemoglobin I (oxyHbI) in the absence and presence of H2S, using Small and Wide Angle X-ray Scattering (SAXS/WAXS) technique. Three dimensional models were also generated from the SAXS/WAXS data. The results show that SMb formation, produced by oxyMb and H2S interaction, induces a change in the protein conformation where its envelope has a very small cleft and the protein is more flexible, less rigid and compact. Based on the direct relationship between Mb's structural conformation and its functionality, we suggest that the conformational change observed upon SMb formation plays a contribution to the protein decrease in O2 affinity and, therefore, on its functionality.

Microsecond Rearrangements of Hydrophobic Clusters in an Initially Collapsed Globule Prime Structure Formation during the Folding of a Small Protein

Determining how polypeptide chain collapse initiates structure formation during protein folding is a long standing goal. It has been challenging to characterize experimentally the dynamics of the polypeptide chain, which lead to the formation of a compact kinetic molten globule (MG) in about a millisecond. In this study, the sub-millisecond events that occur early during the folding of monellin from the guanidine hydrochloride-unfolded state have been characterized using multiple fluorescence and fluorescence resonance energy transfer probes. The kinetic MG is shown to form in a noncooperative manner from the unfolded (U) state as a result of at least three different processes happening during the first millisecond of folding. Initial chain compaction completes within the first 37μs, and further compaction occurs only after structure formation commences at a few milliseconds of folding. The transient nonnative and native-like hydrophobic clusters with side chains of certain residues buried form during the initial chain collapse and the nonnative clusters quickly disassemble. Subsequently, partial chain desolvation occurs, leading to the formation of a kinetic MG. The initial chain compaction and subsequent chain rearrangement appear to be barrierless processes. The two structural rearrangements within the collapsed globule appear to prime the protein for the actual folding transition.

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.

Where the complex things are: single molecule and ensemble spectroscopic investigations of protein folding dynamics

Progress in our understanding of the simple folding dynamics of small proteins and the complex dynamics of large proteins is reviewed. Recent characterizations of the folding transition path of small proteins revealed a simple dynamics explainable by the native centric model. In contrast, the accumulated data showed the substates containing residual structures in the unfolded state and partially populated intermediates, causing complexity in the early folding dynamics of small proteins. The size of the unfolded proteins in the absence of denaturants is likely expanded but still controversial. The steady progress in the observation of folding of large proteins has clarified the rapid formation of long-range contacts that seem inconsistent with the native centric model, suggesting that the folding strategy of large proteins is distinct from that of small proteins.

279(Val→Phe) Polymorphism of lipoprotein-associated phospholipase A(2) resulted in changes of folding kinetics and recognition to substrate

INTRODUCTION

PLA2G7 encodes Lp-PLA2 having role in the formation of atherosclerotic plaques by catalyzing its substrate, phosphatydilcholine (PC), to be pro-inflammatory substances. The increased risk for coronary artery disease (CAD) in Asian population has been related with this enzyme. 279(Val→Phe) variant was reported to have a protective role against CAD due to, in part, secretion defect or loss of enzymatic function. Therefore, We study folding kinetics and enzyme-substrate interaction in 279(Val→Phe) by using clinical and computational biology approach.

METHODS

Polymorphisms were detected by genotyping among 103 acute myocardial infarction patients and 37 controls. Folding Lp-PLA2 was simulated using GROMACS software by assessing helicity, hydrogen bond formation and stability. The interactions of Lp-PLA2 and its substrate were simulated using Pyrx software followed by molecular dynamics simulation using YASARA software.

RESULT

Polymorphism of 279(Val→Phe) was represented by the change of nucleotide from G to T of 994th PLA2G7 gene. The folding simulation suggested a decreased percentage of α-helix, hydrogen bond formation, hydrogen bond stability and hydrophobicity in 279(Val→Phe). The PC did not interact with active site of 279(Val→Phe) as paradoxically observed in 279 valine. 279(Val→Phe) polymorphism is likely to cause unstable binding to the substrate and decrease the enzymatic activity as observed in molecular dynamics simulations. The results of our computational biology study supported a protected effect of 279(Val→Phe) Polymorphism showed by the odd ratio for MI of 0.22 (CI 95% 0.035-1.37) in this study.

CONCLUSION

279(Val→Phe) Polymorphism of Lp-PLA2 may lead to decrease the enzymatic activity via changes of folding kinetics and recognition to its substrate.

Sequence- and Temperature-Dependent Properties of Unfolded and Disordered Proteins from Atomistic Simulations

We use all-atom molecular simulation with explicit solvent to study the properties of selected intrinsically disordered proteins and unfolded states of foldable proteins, which include chain dimensions and shape, secondary structure propensity, solvent accessible surface area, and contact formation. We find that the qualitative scaling behavior of the chains matches expectations from theory under ambient conditions. In particular, unfolded globular proteins tend to be more collapsed under the same conditions than charged disordered sequences of the same length. However, inclusion of explicit solvent in addition naturally captures temperature-dependent solvation effects, which results in an initial collapse of the chains as temperature is increased, in qualitative agreement with experiment. There is a universal origin to the collapse, revealed in the change of hydration of individual residues as a function of temperature: namely, that the initial collapse is driven by unfavorable solvation free energy of individual residues, which in turn has a strong temperature dependence. We also observe that in unfolded globular proteins, increased temperature also initially favors formation of native-like (rather than non-native-like) structure. Our results help to establish how sequence encodes the degree of intrinsic disorder or order as well as its response to changes in environmental conditions.

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.

Calmodulin inhibitors from natural sources: an update

Calmodulin (CaM) plays a central role in regulating a myriad of cellular functions in physiological and pathophysiological processes, thus representing an important drug target. In previous reviews, our group has reported relevant information regarding natural anti-CaM compounds up to 2009. Natural sources continue to provide a diverse and unique reservoir of CaM inhibitors for drug and research tool discovery. This review provides an update of natural products with reported CaM inhibitory properties, which includes around 70 natural products and some synthetic analogues, belonging to different structural classes. Most of these natural inhibitors were isolated from fungi and plants and belong to the stilbenoid, polyketide, alkaloid, and peptide structural classes. These products were discovered mainly using a fluorescence-based method on rationally designed biosensors, which are highly specific, low-cost, and selective and have short reaction times. The effect of several antimitotic drugs on Ca(2+)-hCaM is also described.

Dimer domain swapping versus monomer folding in apo-myoglobin studied by molecular simulations

Using a coarse-grained symmetrized Go model, we performed a series of folding simulations of two apo-myoglobin molecules restrained at a high density, addressing competition of formation of a domain-swapped dimer with folding to two monomer structures.

In situ microfluidic dialysis for biological small-angle X-ray scattering

Owing to the demand for low sample consumption and automated sample changing capabilities at synchrotron small-angle X-ray (solution) scattering (SAXS) beamlines, X-ray microfluidics is receiving continuously increasing attention. Here, a remote-controlled microfluidic device is presented for simultaneous SAXS and ultraviolet absorption measurements during protein dialysis, integrated directly on a SAXS beamline. Microfluidic dialysis can be used for monitoring structural changes in response to buffer exchange or, as demonstrated, protein concentration. By collecting X-ray data during the concentration procedure, the risk of inducing protein aggregation due to excessive concentration and storage is eliminated, resulting in reduced sample consumption and improved data quality. The proof of concept demonstrates the effect of halted or continuous flow in the microfluidic device. No sample aggregation was induced by the concentration process at the levels achieved in these experiments. Simulations of fluid dynamics and transport properties within the device strongly suggest that aggregates, and possibly even higher-order oligomers, are preferentially retained by the device, resulting in incidental sample purification. Hence, this versatile microfluidic device enables investigation of experimentally induced structural changes under dynamically controllable sample conditions.

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.

Ultrafast hydrogen exchange reveals specific structural events during the initial stages of folding of cytochrome c

Many proteins undergo a sharp decrease in chain dimensions during early stages of folding, prior to the rate-limiting step in folding. However, it remains unclear whether compact states are the result of specific folding events or a general hydrophobic collapse of the poly peptide chain driven by the change in solvent conditions. To address this fundamental question, we extended the temporal resolution of NMR-detected H/D exchange labeling experiments into the microsecond regime by adopting a microfluidics approach. By observing the competition between H/D exchange and folding as a function of labeling pH, coupled with direct measurement of exchange rates in the unfolded state, we were able to monitor hydrogen-bond formation for over 50 individual backbone NH groups within the initial 140 microseconds of folding of horse cytochrome c. Clusters of solvent-shielded amide protons were observed in two α-helical segments in the C-terminal half of the protein, while the N-terminal helix remained largely unstructured, suggesting that proximity in the primary structure is a major factor in promoting helix formation and association at early stages of folding, while the entropically more costly long-range contacts between the N- and C-terminal helices are established only during later stages. Our findings clearly indicate that the initial chain condensation in cytochrome c is driven by specific interactions among a subset of α-helical segments rather than a general hydrophobic collapse.

Average conformations determined from PRE data provide high-resolution maps of transient tertiary interactions in disordered proteins

In the last decade it has become evident that disordered states of proteins play important physiological and pathological roles and that the transient tertiary interactions often present in these systems can play a role in their biological activity. The structural characterization of such states has so far largely relied on ensemble representations, which in principle account for both their local and global structural features. However, these approaches are inherently of low resolution due to the large number of degrees of freedom of conformational ensembles and to the sparse nature of the experimental data used to determine them. Here, we overcome these limitations by showing that tertiary interactions in disordered states can be mapped at high resolution by fitting paramagnetic relaxation enhancement data to a small number of conformations, which can be as low as one. This result opens up the possibility of determining the topology of cooperatively collapsed and hidden folded states when these are present in the vast conformational landscape accessible to disordered states of proteins. As a first application, we study the long-range tertiary interactions of acid-unfolded apomyoglobin from experimentally measured paramagnetic relaxation enhancement data.

A simple three-dimensional-focusing, continuous-flow mixer for the study of fast protein dynamics

We present a simple, yet flexible microfluidic mixer with a demonstrated mixing time as short as 80 μs that is widely accessible because it is made of commercially available parts. To simplify the study of fast protein dynamics, we have developed an inexpensive continuous-flow microfluidic mixer, requiring no specialized equipment or techniques. The mixer uses three-dimensional, hydrodynamic focusing of a protein sample stream by a surrounding sheath solution to achieve rapid diffusional mixing between the sample and sheath. Mixing initiates the reaction of interest. Reactions can be spatially observed by fluorescence or absorbance spectroscopy. We characterized the pixel-to-time calibration and diffusional mixing experimentally. We achieved a mixing time as short as 80 μs. We studied the kinetics of horse apomyoglobin (apoMb) unfolding from the intermediate (I) state to its completely unfolded (U) state, induced by a pH jump from the initial pH of 4.5 in the sample stream to a final pH of 2.0 in the sheath solution. The reaction time was probed using the fluorescence of 1-anilinonaphthalene-8-sulfonate (1,8-ANS) bound to the folded protein. We observed unfolding of apoMb within 760 μs, without populating additional intermediate states under these conditions. We also studied the reaction kinetics of the conversion of pyruvate to lactate catalyzed by lactate dehydrogenase using the intrinsic tryptophan emission of the enzyme. We observe sub-millisecond kinetics that we attribute to Michaelis complex formation and loop domain closure. These results demonstrate the utility of the three-dimensional focusing mixer for biophysical studies of protein dynamics.

Nonuniform chain collapse during early stages of staphylococcal nuclease folding detected by fluorescence resonance energy transfer and ultrarapid mixing methods

The development of tertiary structure during folding of staphylococcal nuclease (SNase) was studied by time-resolved fluorescence resonance energy transfer measured using continuous- and stopped-flow techniques. Variants of this two-domain protein containing intradomain and interdomain fluorescence donor/acceptor pairs (Trp and Cys-linked fluorophore or quencher) were prepared to probe the intradomain and interdomain structural evolution accompanying SNase folding. The intra-domain donor/acceptor pairs are within the β-barrel domain (Trp27/Cys64 and Trp27/Cys97) and the interdomain pair is between the α-helical domain and the β-barrel domain (Trp140/Cys64). Time-resolved energy transfer efficiency accompanying folding and unfolding at different urea concentrations was measured over a time range from 30 μs to ≈ 10 s. Information on average donor/acceptor distances at different stages of the folding process was obtained by using a quantitative kinetic modeling approach. The average distance for the donor/acceptor pairs in the β-barrel domain decreases to nearly native values whereas that of the interdomain donor/acceptor pairs remains unchanged in the earliest intermediate (<500 μs of refolding). This indicates a rapid nonuniform collapse resulting in an ensemble of heterogeneous conformations in which the central region of the β-barrel domain is well developed while the C-terminal α-helical domain remains disordered. The distance between Trp140 and Cys64 decreases to native values on the 100-ms time scale, indicating that the α-helical domain docks onto the preformed β-barrel at a late stage of the folding. In addition, the unfolded state is found to be more compact under native conditions, suggesting that changes in solvent conditions may induce a nonspecific hydrophobic collapse.

Submillisecond protein folding events monitored by rapid mixing and mass spectrometry-based oxidative labeling

Kinetic measurements can provide insights into protein folding mechanisms. However, the initial (submillisecond) stages of folding still represent a formidable analytical challenge. A number of ultrarapid triggering techniques have been available for some time, but coupling of these techniques with detection methods that are capable of providing detailed structural information has proven to be difficult. The current work addresses this issue by combining submillisecond mixing with laser-induced oxidative labeling. Apomyoglobin (aMb) serves as a model system for our measurements. Exposure of the protein to a brief pulse of hydroxyl radical (·OH) at different time points during folding introduces covalent modifications at solvent accessible side chains. The extent of labeling is monitored using mass spectrometry-based peptide mapping, providing spatially resolved measurements of changes in solvent accessibility. The submillisecond mixer used here improves the time resolution by a factor of 50 compared to earlier ·OH labeling experiments from our laboratory. Data obtained in this way indicate that early aMb folding events are driven by both local and sequence-remote docking of hydrophobic side chains. Assembly of a partially formed A(E)G(H) scaffold after 0.2 ms is followed by stepwise consolidation that ultimately yields the native state. Major conformational changes go to completion within 0.1 s. The technique introduced here is capable of providing in-depth structural information on very short time scales that have thus far been dominated by low resolution (global) spectroscopic probes. By employing submillisecond mixing in conjunction with slower mixing techniques, it is possible to observe complete folding pathways, from fractions of a millisecond all the way to minutes.

Misfolding and amyloid aggregation of apomyoglobin

Apomyoglobin is an excellent example of a monomeric all α-helical globular protein whose folding pathway has been extensively studied and well characterized. Structural perturbation induced by denaturants or high temperature as well as amino acid substitution have been described to induce misfolding and, in some cases, aggregation. In this article, we review the molecular mechanism of the aggregation process through which a misfolded form of a mutated apomyoglobin aggregates at physiological pH and room temperature forming an amyloid fibril. The results are compared with data showing that either amyloid or aggregate formation occurs under particular denaturing conditions or upon cleavage of the residues corresponding to the C-terminal helix of apomyoglobin. The results are discussed in terms of the sequence regions that are more important than others in determining the amyloid aggregation process.

Molecular dynamics simulations of certain mutant peptide models from staphylococcal nuclease reveal that initial hydrophobic collapse associated with turn propensity drives β-hairpin folding

An important nucleation event during the folding of staphylococcal nuclease involves the formation of a β-hairpin by the sequence (21) DTVKLMYKGQPMTFR(35) . Earlier studies show that the turn sequence 'YKGQP' has an important role in the folding of this β-hairpin. To understand the active or passive nature of the turn sequence 'YKGQP' in the folding of the aforementioned β-hairpin sequence, we studied glycine mutant peptides Ac-(2) DTVKLMYGGQPMTFR(16) -NMe (K9G:15), Ac-(2) DTVKLMYKGGPMTFR(16) -NMe (Q11G:15), Ac-(2) DTVKLMYGGGPMTFR(16) -NMe (K9G/Q11G:15), and Ac-(2) DTVKLMGGGGGMTFR(16) -NMe (penta-G:15) by using molecular dynamics simulations, starting with two different unfolded states, polyproline II and extended conformational forms. Further, 5mer mutant turn peptides Ac-(2) YGGQP(6) -NMe (K3G:5), Ac-(2) YKGGP(6) -NMe (Q5G:5), Ac-(2) YGGGP(6) -NMe (K3G/Q5G:5), and Ac-(2) GGGGG(6) -NMe (penta-G:5) were also studied individually. Our results show that an initial hydrophobic collapse and loop closure occurs in all 15mer mutants, but only K9G:15 mutant forms a stable native-like β-hairpin. In the other 15mer mutants, the hydrophobic collapsed state would not proceed to β-hairpin formation. Of the different simulations performed for the penta-G:15 mutant, in only one simulation a nonnative β-hairpin conformation is sampled with highly flexible loop region ((8) GGGGG(12) ), which has no specific conformational preference as a 5mer. While the sequence 'YGGQP' in the K3G:5 simulation shows relatively higher β-turn propensity, the presence of this sequence in K9G:15 peptide seems to be driving the β-hairpin formation. Thus, these results seem to suggest that for the formation of a stable β-hairpin, the initial hydrophobic collapse is to be assisted by a turn propensity. Initial hydrophobic collapse alone is not sufficient to guide β-hairpin formation.

Sub-millisecond chain collapse of the Escherichia coli globin ApoHmpH

Myoglobins are ubiquitous proteins that play a seminal role in oxygen storage, transport, and NO metabolism. The folding mechanism of apomyoglobins from different species has been studied to a fair extent over the last two decades. However, integrated investigations of the entire process, including both the early (sub-ms) and late (ms-s) folding stages, have been missing. Here, we study the folding kinetics of the single-Trp Escherichia coli globin apoHmpH via a combination of continuous-flow microfluidic and stopped-flow approaches. A rich series of molecular events emerges, spanning a very wide temporal range covering more than 7 orders of magnitude, from sub-microseconds to tens of seconds. Variations in fluorescence intensity and spectral shifts reveal that the protein region around Trp120 undergoes a fast collapse within the 8 μs mixing time and gradually reaches a native-like conformation with a half-life of 144 μs from refolding initiation. There are no further fluorescence changes beyond ca. 800 μs, and folding proceeds much more slowly, up to 20 s, with acquisition of the missing helicity (ca. 30%), long after consolidation of core compaction. The picture that emerges is a gradual acquisition of native structure on a free-energy landscape with few large barriers. Interestingly, the single tryptophan, which lies within the main folding core of globins, senses some local structural consolidation events after establishment of native-like core polarity (i.e., likely after core dedydration). In all, this work highlights how the main core of the globin fold is capable of becoming fully native efficiently, on the sub-millisecond time scale.