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

On the Dynamic Mechanism of Long-Flexible Fatty Acid Binding to Fatty Acid Binding Protein: Resolving the Long-Standing Debate

Fatty acids (FAs) are one of the essential energy sources for physiological processes, and they play a vital role in regulating immune and inflammatory responses, promoting cell differentiation and apoptosis, and inhibiting tumor growth. These functions are carried out by FA binding proteins (FABPs) that recognize and transport FAs. Although the crystal structure of the FA-FABPs complex has long been characterized, the mechanism behind FA binding and dissociation from FABP remains unclear. This study employed conventional MD simulations and enhanced sampling technologies to investigate the atomic-scale complexes of heart fatty acid binding proteins and stearic acid (SA). The results revealed two primary pathways for the binding or dissociation of the flexible long-chain ligand, with the orientation of the SA carboxyl head during dissociation determining the chosen path. Conformational changes in the portal region of FABP during the ligand binding/unbinding were found to be trivial, and the overturn of the ″cap″ or the unfolding of the α2 helix was not required. This study resolves the long-standing debate on the binding mechanism of SA with the long-flexible tail to FABP, which significantly improves the understanding of the transport mechanism of FABPs and the development of related therapeutic agents.

A combined computational-biophysical approach to understanding fatty acid binding to FABP7

Members of the fatty acid binding protein (FABP) family function as intracellular transporters of long-chain fatty acids and other hydrophobic molecules to different cellular compartments. Brain FABP (FABP7) exhibits ligand-directed differences in cellular transport. For example, when FABP7 binds to docosahexaenoic acid (DHA), the complex relocates to the nucleus and influences transcriptional activity, whereas FABP7 bound with monosaturated fatty acids remains in the cytosol. Preferential binding of FABP7 to polyunsaturated fatty acids like DHA has been previously observed and is thought to play a role in differential localization. However, we find that at 37°C, FABP7 does not display strong selectivity, suggesting that the conformational ensemble of FABP7 and its perturbation upon binding may be important. We use molecular dynamics simulations, NMR, and a variety of biophysical techniques to better understand the conformational ensemble of FABP7, how it is perturbed by fatty acid binding, and how this may be related to ligand-directed transport. We find that FABP7 has high degree of conformational heterogeneity that is substantially reduced upon ligand binding. We also observe substantial heterogeneity in ligand binding poses, which is consistent with our finding that ligand binding is resistant to mutations in key polar residues in the binding pocket. Our NMR experiments show that DHA binding leads to chemical shift perturbations in residues near the nuclear localization signal, which may point toward a mechanism of differential transport.

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.

Structure of an Unfolding Intermediate of an RRM Domain of ETR-3 Reveals Its Native-like Fold

Prevalence of one or more partially folded intermediates during protein unfolding with different secondary and ternary conformations has been identified as an integral character of protein unfolding. These transition-state species need to be characterized structurally for elucidation of their folding pathways. We have determined the three-dimensional structure of an intermediate state with increased conformational space sampling under urea-denaturing condition. The protein unfolds completely at 10 M urea but retains residual secondary structural propensities with restricted motion. Here, we describe the native state, observable intermediate state, and unfolded state for ETR-3 RRM-3, which has canonical RRM fold. These observations can shed more light on unfolding events for RRM-containing proteins.

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.

Atomistic Insights into the Functional Instability of the Second Helix of Fatty Acid Binding Protein

Structural dynamics of fatty acid binding proteins (FABPs), which accommodate poorly soluble ligands in the internalized binding cavities, are intimately related to their function. Recently, local unfolding of the α-helical cap in a variant of human intestinal FABP (IFABP) has been shown to correlate with the kinetics of ligand association, shedding light on the nature of the critical conformational reorganization. Yet, the physical origin and mechanism of the functionally relevant transient unfolding remain elusive. Here, we investigate the intrinsic structural instability of the second helix (αII) of IFABP in comparison with other segments of the protein using hydrogen-exchange NMR spectroscopy, microsecond molecular dynamics simulations, and enhanced sampling techniques. Although tertiary interactions positively contribute to the stability of helices in IFABP, the intrinsic unfolding tendency of αII is encoded in its primary sequence and can be described by the Lifson-Roig theory in the absence of tertiary interactions. The unfolding pathway of αII in intact proteins involves an on-pathway intermediate state that is characterized with the fraying of the last helical turn, captured by independent enhanced sampling methods. The simulations in this work, combined with hydrogen-exchange NMR data, provide new, to our knowledge, atomistic insights into the functional local unfolding of FABPs.

The Observation of Ligand-Binding-Relevant Open States of Fatty Acid Binding Protein by Molecular Dynamics Simulations and a Markov State Model

As a member of the fatty acids transporter family, the heart fatty acid binding proteins (HFABPs) are responsible for many important biological activities. The binding mechanism of fatty acid with FABP is critical to the understanding of FABP functions. The uncovering of binding-relevant intermediate states and interactions would greatly increase our knowledge of the binding process. In this work, all-atom molecular dynamics (MD) simulations were performed to characterize the structural properties of nativelike intermediate states. Based on multiple 6 μs MD simulations and Markov state model (MSM) analysis, several "open" intermediate states were observed. The transition rates between these states and the native closed state are in good agreement with the experimental measurements, which indicates that these intermediate states are binding relevant. As a common property in the open states, the partially unfolded α2 helix generates a larger portal and provides the driving force to facilitate ligand binding. On the other side, there are two kinds of open states for the ligand-binding HFABP: one has the partially unfolded α2 helix, and the other has the looser β-barrel with disjointing βD-βE strands. Our results provide atomic-level descriptions of the binding-relevant intermediate states and could improve our understanding of the binding mechanism.

Enhanced dynamics of conformationally heterogeneous T7 bacteriophage lysozyme native state attenuates its stability and activity

Proteins are dynamic in nature and exist in a set of equilibrium conformations on various timescale motions. The flexibility of proteins governs various biological functions, and therefore elucidation of such functional dynamics is essential. In this context, we have studied the structure–dynamics–stability–activity relationship of bacteriophage T7 lysozyme/endolysin (T7L) native-state ensemble in the pH range of 6–8. Our studies established that T7L native state is conformationally heterogeneous, as several residues of its C-terminal half are present in two conformations (major and minor) in the slow exchange time scale of nuclear magnetic resonance (NMR). Structural and dynamic studies suggested that the residues belonging to minor conformations do exhibit native-like structural and dynamic features. Furthermore, the NMR relaxation experiments unraveled that the native state is highly dynamic and the dynamic behavior is regulated by the pH, as the pH 6 conformation exhibited enhanced dynamics compared with pH 7 and 8. The stability measurements and cell-based activity studies on T7L indicated that the native protein at pH 6 is ∼2 kcal less stable and is ∼50% less active than those of pH 7 and 8. A comprehensive analysis of the T7L active site, unfolding initiation sites and the residues with altered dynamics outlined that the attenuation of stability and activity is a resultant of its enhanced dynamic properties, which, in turn, can be attributed to the protonation/deprotonation of its partially buried His residues. Our study on T7L structure–dynamics–activity paradigm could assist in engineering novel amidase-based endolysins with enhanced activity and stability over a broad pH range.

Structural basis and mechanism of the unfolding-induced activation of HdeA, a bacterial acid response chaperone

The role of protein structural disorder in biological functions has gained increasing attention in the past decade. The bacterial acid-resistant chaperone HdeA belongs to a group of "conditionally disordered" proteins, because it is inactive in its well-structured state and becomes activated via an order-to-disorder transition under acid stress. However, the mechanism for unfolding-induced activation remains unclear because of a lack of experimental information on the unfolded state conformation and the chaperone-client interactions. Herein, we used advanced solution NMR methods to characterize the activated-state conformation of HdeA under acidic conditions and identify its client-binding sites. We observed that the structure of activated HdeA becomes largely disordered and exposes two hydrophobic patches essential for client interactions. Furthermore, using the pH-dependent chemical exchange saturation transfer (CEST) NMR method, we identified three acid-sensitive regions that act as structural locks in regulating the exposure of the two client-binding sites during the activation process, revealing a multistep activation mechanism of HdeA's chaperone function at the atomic level. Our results highlight the role of intrinsic protein disorder in chaperone function and the self-inhibitory role of ordered structures under nonstress conditions, offering new insights for improving our understanding of protein structure-function paradigms.

Intra- and inter-protein couplings of backbone motions underlie protein thiol-disulfide exchange cascade

The thioredoxin (Trx)-coupled arsenate reductase (ArsC) is a family of enzymes that catalyzes the reduction of arsenate to arsenite in the arsenic detoxification pathway. The catalytic cycle involves a series of relayed intramolecular and intermolecular thiol-disulfide exchange reactions. Structures at different reaction stages have been determined, suggesting significant conformational fluctuations along the reaction pathway. Herein, we use two state-of-the-art NMR methods, the chemical exchange saturation transfer (CEST) and the CPMG-based relaxation dispersion (CPMG RD) experiments, to probe the conformational dynamics of B. subtilis ArsC in all reaction stages, namely the enzymatic active reduced state, the intra-molecular C10-C82 disulfide-bonded intermediate state, the inactive oxidized state, and the inter-molecular disulfide-bonded protein complex with Trx. Our results reveal highly rugged energy landscapes in the active reduced state, and suggest global collective motions in both the C10-C82 disulfide-bonded intermediate and the mixed-disulfide Trx-ArsC complex.

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.

Microsecond Timescale Protein Dynamics: a Combined Solid-State NMR Approach

Conformational exchange in proteins is a major determinant in protein functionality. In particular, the μs-ms timescale is associated with enzymatic activity and interactions between biological molecules. We show here that a comprehensive data set of R1ρ relaxation dispersion profiles employing multiple effective fields and tilt angles can be easily obtained in perdeuterated, partly back-exchanged proteins at fast magic-angle spinning and further complemented with chemical-exchange saturation transfer NMR experiments. The approach exploits complementary sources of information and enables the extraction of multiple exchange parameters for μs-ms timescale conformational exchange, most notably including the sign of the chemical shift differences between the ground and excited states.

Equilibrium folding dynamics of meACP in water, heavy water, and low concentration of urea

Many proteins fold in apparent two-state behavior, as partially folded intermediates only transiently accumulate and easily escape detection. Besides a native form and a mainly unfolded form, we captured a partially unfolded form of an acyl carrier protein from Micromonospora echinospora (meACP) in the folding/unfolding equilibrium using chemical exchange saturation transfer NMR experiments. The C-terminal region of the partially unfolded form is mainly folded and the N-terminal is unfolded. Furthermore, to understand how the folding process of meACP is influenced by solvent environments, we compared the folding dynamics of meACP in D₂O, H₂O and low concentration of urea. As the environment becomes more denaturing from D₂O to H₂O and then to urea, the unfolded state becomes increasingly populated, and the folding rate decreases. Adding a small amount of urea, which does not change solvent viscosity, has little effects on the unfolding rates, while changing H₂O to D₂O reduces the unfolding rates possibly due to the increase of solvent viscosity. The quantified solvent effects on the protein folding Gibbs energy and activation energy suggest that the transition state of folding may have a similar structure to the native state of the protein.

Conformational Flexibility Differentiates Naturally Occurring Bet v 1 Isoforms

The protein Bet v 1 represents the main cause for allergic reactions to birch pollen in Europe and North America. Structurally homologous isoforms of Bet v 1 can have different properties regarding allergic sensitization and Th2 polarization, most likely due to differential susceptibility to proteolytic cleavage. Using NMR relaxation experiments and molecular dynamics simulations, we demonstrate that the initial proteolytic cleavage sites in two naturally occurring Bet v 1 isoforms, Bet v 1.0101 (Bet v 1a) and Bet v 1.0102 (Bet v 1d), are conformationally flexible. Inaccessible cleavage sites in helices and strands are highly flexible on the microsecond-millisecond time scale, whereas those located in loops display faster nanosecond-microsecond flexibility. The data consistently show that Bet v 1.0102 is more flexible and conformationally heterogeneous than Bet v 1.0101. Moreover, NMR hydrogen-deuterium exchange measurements reveal that the backbone amides in Bet v 1.0102 are significantly more solvent exposed, in agreement with this isoform's higher susceptibility to proteolytic cleavage. The differential conformational flexibility of Bet v 1 isoforms, along with the transient exposure of inaccessible sites to the protein surface, may be linked to proteolytic susceptibility, representing a potential structure-based rationale for the observed differences in Th2 polarization and allergic sensitization.