A (15)N CPMG relaxation dispersion experiment more resistant to resonance offset and pulse imperfection

Optimization of TROSY- and anti-TROSY-based 15N CPMG relaxation dispersion experiments through phase cycling

CPMG relaxation dispersion studies of biomolecular dynamics on the μs-ms timescale can provide detailed kinetic, thermodynamic, and structural insights into function. Frequently, the 15N spin serves as the probe of choice, as uniform incorporation of the 15N isotope is facile and cost-effective, and the interpretation of the resulting data is often relatively straightforward. In conventional CPMG relaxation dispersion experiments the application of CPMG pulses with constant radiofrequency (RF) phase can lead to artifactual dispersion profiles that result from off-resonance effects, RF field inhomogeneity, and pulse miscalibration. The development of CPMG experiments with the [0013]-phase cycle has significantly reduced the impact of pulse imperfections over a greater bandwidth of frequency offsets in comparison to constant phase experiments. Application of 15N-TROSY-based CPMG schemes to studies of the dynamics of large molecules is necessary for high sensitivity, yet the correct incorporation of the [0013]-phase cycle is non-trivial. Here we present TROSY- and anti-TROSY-based 15N CPMG experiments with the [0013]-phase cycling scheme and demonstrate, through comprehensive numerical simulations and experimental validation, enhanced resistance to pulse imperfections relative to traditional schemes utilizing constant phase CPMG pulses. Notably, exchange parameters derived from the new experiments are in good agreement with those obtained using other, more established, 15N-based CPMG approaches.

Design of amyloidogenic peptide traps

Segments of proteins with high β-strand propensity can self-associate to form amyloid fibrils implicated in many diseases. We describe a general approach to bind such segments in β-strand and β-hairpin conformations using de novo designed scaffolds that contain deep peptide-binding clefts. The designs bind their cognate peptides in vitro with nanomolar affinities. The crystal structure of a designed protein-peptide complex is close to the design model, and NMR characterization reveals how the peptide-binding cleft is protected in the apo state. We use the approach to design binders to the amyloid-forming proteins transthyretin, tau, serum amyloid A1 and amyloid β1-42 (Aβ42). The Aβ binders block the assembly of Aβ fibrils as effectively as the most potent of the clinically tested antibodies to date and protect cells from toxic Aβ42 species.

Thermodynamic Coupling of the tandem RRM domains of hnRNP A1 underlie its Pleiotropic RNA Binding Functions

The functional properties of RNA-binding proteins (RBPs) require allosteric regulation through inter-domain communication. Despite the foundational importance of allostery to biological regulation, almost no studies have been conducted to describe the biophysical nature by which inter-domain communication manifests in RBPs. Here, we show through high-pressure studies with hnRNP A1 that inter-domain communication is vital for the unique stability of its N- terminal domain containing a tandem of RNA Recognition Motifs (RRMs). Despite high sequence similarity and nearly identical tertiary structures, the two RRMs exhibit drastically different stability under pressure. RRM2 unfolds completely under high-pressure as an individual domain, but when appended to RRM1, it remains stable. Variants in which inter-domain communication is disrupted between the tandem RRMs show a large decrease in stability under pressure. Carrying these mutations over to the full-length protein for in vivo experiments revealed that the mutations affected the ability of the disordered C-terminus to engage in protein-protein interactions and more importantly, they also influenced the RNA binding capacity. Collectively, this work reveals that thermodynamic coupling between the tandem RRMs of hnRNP A1 accounts for its allosteric regulatory functions.

Exploiting conformational dynamics to modulate the function of designed proteins

With the recent success in calculating protein structures from amino acid sequences using artificial intelligence-based algorithms, an important next step is to decipher how dynamics is encoded by the primary protein sequence so as to better predict function. Such dynamics information is critical for protein design, where strategies could then focus not only on sequences that fold into particular structures that perform a given task, but would also include low-lying excited protein states that could influence the function of the designed protein. Herein, we illustrate the importance of dynamics in modulating the function of C34, a designed α/β protein that captures β-strands of target ligands and is a member of a family of proteins designed to sequester β-strands and β hairpins of aggregation-prone molecules that lead to a variety of pathologies. Using a strategy to "see" regions of apo C34 that are invisible to NMR spectroscopy as a result of pervasive conformational exchange, as well as a mutagenesis approach whereby C34 molecules are stabilized into a single conformer, we determine the structures of the predominant conformations that are sampled by C34 and show that these attenuate the affinity for cognate peptide. Subsequently, the observed motion is exploited to develop an allosterically regulated peptide binder whose binding affinity can be controlled through the addition of a second molecule. Our study emphasizes the unique role that NMR can play in directing the design process and in the construction of new molecules with more complex functionality.

Protein Conformational Exchanges Modulated by the Environment of Outer Membrane Vesicles

Protein function, in many cases, is strongly coupled to the dynamics and conformational equilibria of the protein. The environment surrounding proteins is critical for their dynamics and can dramatically affect the conformational equilibria and subsequently the activities of proteins. However, it is unclear how protein conformational equilibria are modulated by their crowded native environments. Here we reveal that outer membrane vesicle (OMV) environments modulate the conformational exchanges of Im7 protein at its local frustrated sites and shift the conformation toward its ground state. Further experiments show both macromolecular crowding and quinary interactions with the periplasmic components stabilize the ground state of Im7. Our study highlights the key role that the OMV environment plays in the protein conformational equilibria and subsequently the conformation-related protein functions. Furthermore, the long-lasting nuclear magnetic resonance measurement time of proteins within OMVs indicates that they could serve as a promising system for investigating protein structures and dynamics in situ via nuclear magnetic spectroscopy.

Slow conformational changes in the rigid and highly stable chymotrypsin inhibitor 2

Slow conformational changes are often directly linked to protein function. It is however less clear how such processes may perturb the overall folding stability of a protein. We previously found that the stabilizing double mutant L49I/I57V in the small protein chymotrypsin inhibitor 2 from barley led to distributed increased nanosecond and faster dynamics. Here we asked what effects the L49I and I57V substitutions, either individually or together, have on the slow conformational dynamics of CI2. We used ¹⁵ N CPMG spin relaxation dispersion experiments to measure the kinetics, thermodynamics and structural changes associated with slow conformational change in CI2. These changes result in an excited state that is populated to 4.3% at 1 °C. As the temperature is increased the population of the excited state decreases. Structural changes in the excited state are associated with residues that interact with water molecules that have well defined positions and are found at these positions in all crystal structures of CI2. The substitutions in CI2 have only little effect on the structure of the excited state whereas the stability of the excited state to some extent follows the stability of the main state. The minor state is thus most populated for the most stable CI2 variant and least populated for the least stable variant. We hypothesize that the interactions between the substituted residues and the well-ordered water molecules links subtle structural changes around the substituted residues to the region in the protein that experience slow conformational changes. This article is protected by copyright. All rights reserved.

Effects of ligand binding on dynamics of fatty acid binding protein and interactions with membranes

Intracellular transport of fatty acids involves binding of ligands to their carrier fatty acid binding proteins (FABPs) and interactions of ligand-free and -bound FABPs with membranes. Previous studies focused on ligand-free FABPs. Here, our amide hydrogen exchange data showed that oleic acid binding to human intestinal FABP (hIFABP) stabilizes the protein, most likely through enhancing the hydrogen-bonding network, and induces rearrangement of sidechains even far away from the ligand binding site. Using NMR relaxation techniques, we found that the ligand binding affects not only conformational exchanges between major and minor states but also the affinity of hIFABP to nanodiscs. Analyses of the relaxation and amide exchange data suggested that two minor native-like states existing in both ligand-free and -bound hIFABPs originate from global "breathing" motions, while one minor native-like state comes from local motions. The amide hydrogen exchange data also indicated that helix αII undergoes local unfolding through which ligands can exit from the binding cavity.

Conformational exchange of fatty acid binding protein induced by protein-nanodisc interactions

Fatty acid binding proteins (FABPs) can facilitate the transfer of long-chain fatty acids between intracellular membranes across considerable distances. The transfer process involves fatty acids, their donor membrane and acceptor membrane, and FABPs, implying that potential protein-membrane interactions exist. Despite intensive studies on FABP-membrane interactions, the interaction mode remains elusive, and the protein-membrane association and dissociation rates are inconsistent. In this study, we used nanodiscs (NDs) as mimetic membranes to investigate FABP-membrane interactions. Our NMR experiments showed that human intestinal FABP interacts weakly with both negatively charged and neutral membranes, but it prefers the negatively charged one. Through simultaneous analysis of NMR relaxation in the rotating-frame (R_1ρ), relaxation dispersion, chemical exchange saturation transfer, and dark-state exchange saturation transfer data, we estimated the affinity of the protein to negatively charged NDs, the dissociation rate, and apparent association rate. We further showed that the protein in the ND-bound state adopts a conformation different from the native structure and the second helix is very likely involved in interactions with NDs. We also found a membrane-induced FABP conformational state that exists only in the presence of NDs. This state is native-like, different from other conformational states in structure, unbound to NDs, and in dynamic equilibrium with the ND-bound state.

Dengue virus protease activity modulated by dynamics of protease cofactor

The viral protease domain (NS3pro) of dengue virus is essential for virus replication, and its cofactor NS2B is indispensable for the proteolytic function. Although several NS3pro-NS2B complex structures have been obtained, the dynamic property of the complex remains poorly understood. Using NMR relaxation techniques, here we found that NS3pro-NS2B exists in both closed and open conformations that are in dynamic equilibrium on a submillisecond timescale in aqueous solution. Our structural information indicates that the C-terminal region of NS2B is disordered in the minor open conformation but folded in the major closed conformation. Using mutagenesis, we showed that the closed-open conformational equilibrium can be shifted by changing NS2B stability. Moreover, we revealed that the proteolytic activity of NS3pro-NS2B correlates well with the population of the closed conformation. Our results suggest that the closed-open conformational equilibrium can be used by both nature and humanity to control the replication of dengue virus.

Protein Dynamics Influence the Enzymatic Activity of Phospholipase A/Acyltransferases 3 and 4

Phospholipase A/acyltransferase 3 (PLAAT3) and PLAAT4 are enzymes involved in the synthesis of bioactive lipids. Despite sequential and structural similarities, the two enzymes differ in activity and specificity. The relation between the activity and dynamics of the N-terminal domains of PLAAT3 and PLAAT4 was studied. PLAAT3 has a much higher melting temperature and exhibits less nanosecond and millisecond dynamics in the active site, in particular in loop L2(B6), as shown by NMR spectroscopy and molecular dynamics calculations. Swapping the L2(B6) loops between the two PLAAT enzymes results in strongly increased phospholipase activity in PLAAT3 but no reduction in PLAAT4 activity, indicating that this loop contributes to the low activity of PLAAT3. The results show that, despite structural similarity, protein dynamics differ substantially between the PLAAT variants, which can help to explain the activity and specificity differences.

Ligand Entry into Fatty Acid Binding Protein via Local Unfolding Instead of Gap Widening

Fatty acid binding proteins play an important role in the transportation of fatty acids. Despite intensive studies, how fatty acids enter the protein cavity for binding is still controversial. Here, a gap-closed variant of human intestinal fatty acid binding protein was generated by mutagenesis, in which the gap is locked by a disulfide bridge. According to its structure determined here by NMR, this variant has no obvious openings as the ligand entrance and the gap cannot be widened by internal dynamics. Nevertheless, it still takes up fatty acids and other ligands. NMR relaxation dispersion, chemical exchange saturation transfer, and hydrogen-deuterium exchange experiments show that the variant exists in a major native state, two minor native-like states, and two locally unfolded states in aqueous solution. Local unfolding of either βB-βD or helix 2 can generate an opening large enough for ligands to enter the protein cavity, but only the fast local unfolding of helix 2 is relevant to the ligand entry process.

Revisiting 1HN CPMG relaxation dispersion experiments: a simple modification can eliminate large artifacts

Carr-Purcell-Meiboom-Gill relaxation dispersion experiments are commonly used to probe biomolecular dynamics on the millisecond timescale. The simplest experiment involves using backbone ¹⁵N spins as probes of motion and pulse sequences are now available for providing accurate dispersion profiles in this case. In contrast, ¹H-based experiments recorded on fully protonated samples are less common because of difficulties associated with homonuclear scalar couplings that can result in transfer of magnetization between coupled spins, leading to significant artifacts. Herein we examine a version of the ¹H^N CPMG experiment that has been used in our laboratory where a pair of CPMG pulse trains comprising non-selective, high power ¹H refocusing pulses sandwich an amide selective pulse that serves to refocus scalar-coupled evolution by the end of the train. The origin of the artifacts in our original scheme is explained and a new, significantly improved sequence is presented. The utility of the new experiment is demonstrated by obtaining flat ¹H^N dispersion profiles in a protonated protein system that is not expected to undergo millisecond timescale dynamics, and subsequently by measuring profiles on a cavity mutant of T4 lysozyme that exchanges between a pair of distinct states, establishing that high quality data can be generated even for fully protonated samples.

The Role of Protein Thermodynamics and Primary Structure in Fibrillogenesis of Variable Domains from Immunoglobulin Light Chains

Immunoglobulin light-chain amyloidosis is a protein aggregation disease that leads to proteinaceous deposits in a variety of organs in the body and, if untreated, ultimately results in death. The mechanisms by which light-chain aggregation occurs are not well understood. Here we have used solution NMR spectroscopy and biophysical studies to probe immunoglobulin variable domain λV6-57 V_L aggregation, a process that appears to drive the degenerative phenotypes in amyloidosis patients. Our results establish that aggregation proceeds via the unfolded state. We identify, through NMR relaxation experiments recorded on the unfolded domain ensemble, a series of hotspots that could be involved in the initial phases of aggregate formation. Mutational analysis of these hotspots reveals that the region that includes K16-R24 is particularly aggregation prone. Notably, this region includes the site of the R24G substitution, a mutation that is found in variable domains of λ light-chain deposits in 25% of patients. The R24G λV6-57 V_L domain aggregates more rapidly than would be expected on the basis of thermodynamic stability alone, while substitutions in many of the aggregation-prone regions significantly slow down fibril formation.

A methyl 1H double quantum CPMG experiment to study protein conformational exchange

Protein conformational changes play crucial roles in enabling function. The Carr-Purcell-Meiboom-Gill (CPMG) experiment forms the basis for studying such dynamics when they involve the interconversion between highly populated and sparsely formed states, the latter having lifetimes ranging from ~ 0.5 to ~ 5 ms. Among the suite of experiments that have been developed are those that exploit methyl group probes by recording methyl ¹H single quantum (Tugarinov and Kay in J Am Chem Soc 129:9514-9521, 2007) and triple quantum (Yuwen et al. in Angew Chem Int Ed Engl 55:11490-11494, 2016) relaxation dispersion profiles. Here we build upon these by developing a third experiment in which methyl ¹H double quantum coherences evolve during a CPMG relaxation element. By fitting single, double, and triple quantum datasets, akin to recording the single quantum dataset at static magnetic fields of B_o, 2B_o and 3B_o, we show that accurate exchange values can be obtained even in cases where exchange rates exceed 10,000 s^-1. The utility of the double quantum experiment is demonstrated with a pair of cavity mutants of T4 lysozyme (T4L) with ground and excited states interchanged and with exchange rates differing by fourfold (~ 900 s^-1 and ~ 3600 s^-1), as well as with a fast-folding domain where the unfolded state lifetime is ~ 80 µs.

Residue selective 15N CEST and CPMG experiments for studies of millisecond timescale protein dynamics

Proteins are intrinsically dynamic molecules and undergo exchanges among multiple conformations to perform biological functions. The CPMG relaxation dispersion and CEST experiments are two important solution NMR techniques for characterizing the conformational exchange processes on the millisecond timescale. Traditional pseudo 3D ¹⁵N CEST and CPMG experiments have certain limitations in their applications. For example, both experiments have low sensitivity for broadened resonances, and the process of optimizing sample conditions and experimental parameters are often time consuming. To overcome these limitations, we herein present a new set of residue selective ¹⁵N CEST and CPMG pulse sequences by employing the Hartmann-Hahn cross-polarization transfer of magnetization in both 1D and 2D schemes. Combined with frequency labeling in the indirect dimension using only a small number of increments, the pulse sequences in the 2D scheme can be applied on resonances in overlapped regions of the ¹H-¹⁵N HSQC spectrum. The pulse sequences were further applied on several proteins, demonstrating their advantages over the traditional CEST and CPMG experiments under specific circumstances.

Analysis of Artifacts Caused by Pulse Imperfections in CPMG Pulse Trains in NMR Relaxation Dispersion Experiments

Nuclear magnetic resonance relaxation dispersion (rd) experiments provide kinetics and thermodynamics information of molecules undergoing conformational exchange. Rd experiments often use a Carr-Purcell-Meiboom-Gill (CPMG) pulse train equally separated by a spin-state selective inversion element (U-element). Even with measurement parameters carefully set, however, parts of 1H–15N correlations sometimes exhibit large artifacts that may hamper the subsequent analyses. We analyzed such artifacts with a combination of NMR measurements and simulation. We found that particularly the lowest CPMG frequency (νcpmg) can also introduce large artifacts into amide 1H–15N and aromatic 1H–13C correlations whose 15N/13C resonances are very close to the carrier frequencies. The simulation showed that the off-resonance effects and miscalibration of the CPMG π pulses generate artifact maxima at resonance offsets of even and odd multiples of νcpmg, respectively. We demonstrate that a method once introduced into the rd experiments for molecules having residual dipolar coupling significantly reduces artifacts. In the method the 15N/13C π pulse phase in the U-element is chosen between x and y. We show that the correctly adjusted sequence is tolerant to miscalibration of the CPMG π pulse power as large as ±10% for most amide 15N and aromatic 13C resonances of proteins.

CPMG Experiments for Protein Minor Conformer Structure Determination

CPMG relaxation dispersion NMR experiments have emerged as a powerful method to characterize protein minor states that are in exchange with a visible dominant conformation, and have lifetimes between ~0.5 and 5 milliseconds (ms) and populations greater than 0.5%. The structure of the minor state can, in favorable cases, be determined from the parameters provided by the CPMG relaxation dispersion experiments. Here, we go through the intricacies of setting up these powerful CPMG experiments.

Ligand Proton Pseudocontact Shifts Determined from Paramagnetic Relaxation Dispersion in the Limit of NMR Intermediate Exchange

Delineation of protein-ligand interaction modes is key for rational drug discovery. The availability of complex crystal structures is often limited by the aqueous solubility of the compounds, while lead-like compounds with micromolar affinities normally fall into the NMR intermediate exchange regime, in which severe line broadening to beyond the detection of interfacial resonances limits NMR applications. Here, we developed a new method to retrieve low-populated bound-state ¹H pseudocontact shifts (PCSs) using paramagnetic relaxation dispersion (RD). We evaluated using a ¹H PCS-RD approach in a BRM bromodomain lead-like inhibitor to filter molecular docking poses using multiple intermolecular structural restraints. Considering the universal presence of proton atoms in druglike compounds, our work will have wide application in structure-guided drug discovery even under an extreme condition of NMR intermediate exchange and low aqueous solubility of ligands.

Removal of slow-pulsing artifacts in in-phase 15N relaxation dispersion experiments using broadband 1H decoupling

Understanding of the molecular mechanisms of protein function requires detailed insight into the conformational landscape accessible to the protein. Conformational changes can be crucial for biological processes, such as ligand binding, protein folding, and catalysis. NMR spectroscopy is exquisitely sensitive to such dynamic changes in protein conformations. In particular, Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments are a powerful tool to investigate protein dynamics on a millisecond time scale. CPMG experiments that probe the chemical shift modulation of 15N in-phase magnetization are particularly attractive, due to their high sensitivity. These experiments require high power 1H decoupling during the CPMG period to keep the 15N magnetization in-phase. Recently, an improved version of the in-phase 15N-CPMG experiment was introduced, offering greater ease of use by employing a single 1H decoupling power for all CPMG pulsing rates. In these experiments however, incomplete decoupling of off-resonance amide 1H spins introduces an artefactual dispersion of relaxation rates, the so-called slow-pulsing artifact. Here, we analyze the slow-pulsing artifact in detail and demonstrate that it can be suppressed through the use of composite pulse decoupling (CPD). We report the performances of various CPD schemes and show that CPD decoupling based on the 90x-240y-90x element results in high-quality dispersion curves free of artifacts, even for amides with high 1H offset.

An NMR strategy to detect conformational differences in a protein complexed with highly analogous inhibitors in solution

This manuscript presents an NMR strategy to investigate conformational differences in protein-inhibitor complexes, when the inhibitors tightly bind to a protein at sub-nanomolar dissociation constants and are highly analogous to each other. Using HIV-1 protease (PR), we previously evaluated amide chemical shift differences, ΔCSPs, of PR bound to darunavir (DRV) compared to PR bound to several DRV analogue inhibitors, to investigate subtle but significant long-distance conformation changes caused by the inhibitor's chemical moiety variation [Khan, S. N., Persons, J. D. Paulsen, J. L., Guerrero, M., Schiffer, C. A., Kurt-Yilmaz, N., and Ishima, R., Biochemistry, (2018), 57, 1652-1662]. However, ΔCSPs are not ideal for investigating subtle PR-inhibitor interface differences because intrinsic differences in the electron shielding of the inhibitors affect protein ΔCSPs. NMR relaxation is also not suitable as it is not sensitive enough to detect small conformational differences in rigid regions among similar PR-inhibitor complexes. Thus, to gain insight into conformational differences at the inhibitor-protein interface, we recorded ¹⁵N-half filtered NOESY spectra of PR bound to two highly analogous inhibitors and assessed NOEs between PR amide protons and inhibitor protons, between PR amide protons and hydroxyl side chains, and between PR amide protons and water protons. We also verified the PR amide-water NOEs using 2D water-NOE/ROE experiments. Differences in water-amide proton NOE peaks, possibly due to amide-protein hydrogen bonds, were observed between subunit A and subunit B, and between the DRV-bound form and an analogous inhibitor-bound form, which may contribute to remote conformational changes.

Simultaneous determination of fast and slow dynamics in molecules using extreme CPMG relaxation dispersion experiments

Molecular dynamics play a significant role in how molecules perform their function. A critical method that provides information on dynamics, at the atomic level, is NMR-based relaxation dispersion (RD) experiments. RD experiments have been utilized for understanding multiple biological processes occurring at micro-to-millisecond time, such as enzyme catalysis, molecular recognition, ligand binding and protein folding. Here, we applied the recently developed high-power RD concept to the Carr-Purcell-Meiboom-Gill sequence (extreme CPMG; E-CPMG) for the simultaneous detection of fast and slow dynamics. Using a fast folding protein, gpW, we have shown that previously inaccessible kinetics can be accessed with the improved precision and efficiency of the measurement by using this experiment.

Determining rotational dynamics of the guanidino group of arginine side chains in proteins by carbon-detected NMR

A new NMR-based method is presented to determine the rotational dynamics around the N_ε–C_ζ bond of arginine to characterise the interactions mediated by arginine side chains.

Arginine residues are imperative for many active sites and protein-interaction interfaces. A new NMR-based method is presented to determine the rotational dynamics around the N_ε–C_ζ bond of arginine side chains. An application to a 19 kDa protein shows that the strengths of interactions involving arginine side chains can be characterised.

Structural Basis for pH-mediated Regulation of F-actin Severing by Gelsolin Domain 1

Six-domain gelsolin regulates actin structural dynamics through its abilities to sever, cap and uncap F-actin. These activities are modulated by various cellular parameters like Ca2+ and pH. Until now, only the molecular activation mechanism of gelsolin by Ca2+ has been understood relatively well. The fragment comprising the first domain and six residues from the linker region into the second domain has been shown to be similar to the full-length protein in F-actin severing activity in the absence of Ca2+ at pH 5. To understand how this gelsolin fragment is activated for F-actin severing by lowering pH, we solved its NMR structures at both pH 7.3 and 5 in the absence of Ca2+ and measured the pKa values of acidic amino acid residues and histidine residues. The overall structure and dynamics of the fragment are not affected significantly by pH. Nevertheless, local structural changes caused by protonation of His29 and Asp109 result in the activation on lowering the pH, and protonation of His151 directly effects filament binding since it resides in the gelsolin/actin interface. Mutagenesis studies support that His29, Asp109 and His151 play important roles in the pH-dependent severing activity of the gelsolin fragment.

Coexistence of multiple minor states of fatty acid binding protein and their functional relevance

Proteins are dynamic over a wide range of timescales, but determining the number of distinct dynamic processes and identifying functionally relevant dynamics are still challenging. Here we present the study on human intestinal fatty acid binding protein (hIFABP) using a novel analysis of ¹⁵N relaxation dispersion (RD) and chemical shift saturation transfer (CEST) experiments. Through combined analysis of the two types of experiments, we found that hIFABP exists in a four-state equilibrium in which three minor states interconvert directly with the major state. According to conversion rates from the major “closed” state to minor states, these minor states are irrelevant to the function of fatty acid transport. Based on chemical shifts of the minor states which could not be determined from RD data alone but were extracted from a combined analysis of RD and CEST data, we found that all the minor states are native-like. This conclusion is further supported by hydrogen-deuterium exchange experiments. Direct conversions between the native state and native-like intermediate states may suggest parallel multitrack unfolding/folding pathways of hIFABP. Moreover, hydrogen-deuterium exchange data indicate the existence of another locally unfolded minor state that is relevant to the fatty acid entry process.

Quantitative measurement of exchange dynamics in proteins via (13)C relaxation dispersion of (13)CHD2-labeled samples

Methyl groups have emerged as powerful probes of protein dynamics with timescales from picoseconds to seconds. Typically, studies involving high molecular weight complexes exploit (13)CH3- or (13)CHD2-labeling in otherwise highly deuterated proteins. The (13)CHD2 label offers the unique advantage of providing (13)C, (1)H and (2)H spin probes, however a disadvantage has been the lack of an experiment to record (13)C Carr-Purcell-Meiboom-Gill relaxation dispersion that monitors millisecond time-scale dynamics, implicated in a wide range of biological processes. Herein we develop an experiment that eliminates artifacts that would normally result from the scalar coupling between (13)C and (2)H spins that has limited applications in the past. The utility of the approach is established with a number of applications, including measurement of ms dynamics of a disease mutant of a 320 kDa p97 complex.

A New Method for Determining Structure Ensemble: Application to a RNA Binding Di-Domain Protein

Structure ensemble determination is the basis of understanding the structure-function relationship of a multidomain protein with weak domain-domain interactions. Paramagnetic relaxation enhancement has been proven a powerful tool in the study of structure ensembles, but there exist a number of challenges such as spin-label flexibility, domain dynamics, and overfitting. Here we propose a new (to our knowledge) method to describe structure ensembles using a minimal number of conformers. In this method, individual domains are considered rigid; the position of each spin-label conformer and the structure of each protein conformer are defined by three and six orthogonal parameters, respectively. First, the spin-label ensemble is determined by optimizing the positions and populations of spin-label conformers against intradomain paramagnetic relaxation enhancements with a genetic algorithm. Subsequently, the protein structure ensemble is optimized using a more efficient genetic algorithm-based approach and an overfitting indicator, both of which were established in this work. The method was validated using a reference ensemble with a set of conformers whose populations and structures are known. This method was also applied to study the structure ensemble of the tandem di-domain of a poly (U) binding protein. The determined ensemble was supported by small-angle x-ray scattering and nuclear magnetic resonance relaxation data. The ensemble obtained suggests an induced fit mechanism for recognition of target RNA by the protein.

Local Unfolding of Fatty Acid Binding Protein to Allow Ligand Entry for Binding

Fatty acid binding proteins are responsible for the transportation of fatty acids in biology. Despite intensive studies, the molecular mechanism of fatty acid entry to and exit from the protein cavity is still unclear. Here a cap-closed variant of human intestinal fatty acid binding protein was generated by mutagenesis, in which the helical cap is locked to the β-barrel by a disulfide linkage. Structure determination shows that this variant adopts a closed conformation, but still uptakes fatty acids. Stopped-flow experiments indicate that a rate-limiting step exists before the ligand association and this step corresponds to the conversion of the closed form to the open one. NMR relaxation dispersion and H-D exchange data demonstrate the presence of two excited states: one is native-like, but the other adopts a locally unfolded structure. Local unfolding of helix 2 generates an opening for ligands to enter the protein cavity, and thus controls the ligand association rate.