Stabilization provided by neighboring strands is critical for the mechanical stability of proteins

Single-molecule force spectroscopy studies and steered molecular dynamics simulations have revealed that protein topology and pulling geometry play important roles in determining the mechanical stability of proteins. Most studies have focused on local interactions that are associated with the force-bearing β-strands. Interactions mediated by neighboring strands are often overlooked. Here we use Top7 and barstar as model systems to illustrate the critical importance of the stabilization effect provided by neighboring β-strands on the mechanical stability. Using single-molecule atomic force microscopy, we showed that Top7 and barstar, which have similar topology in their force-bearing region, exhibit vastly different mechanical-stability characteristics. Top7 is mechanically stable and unfolds at ∼150 pN, whereas barstar is mechanically labile and unfolds largely below 50 pN. Steered molecular dynamics simulations revealed that stretching force peels one force-bearing strand away from barstar to trigger unfolding, whereas Top7 unfolds via a substructure-sliding mechanism. This previously overlooked stabilization effect from neighboring β-strands is likely to be a general mechanism in protein mechanics and can serve as a guideline for the de novo design of proteins with significant mechanical stability and novel protein topology.     

The effect of temperature on mechanical resistance of the native and intermediate states of I27

We investigated the effect of temperature on the mechanical unfolding of I27 from human cardiac titin, employing a custom-built temperature control device for single-molecule atomic force microscopy measurement. A sawtooth pattern was observed in the force curves where each force peak reports on the unfolding of an I27 domain. In early unfolding events, we observed a “hump-like” deviation due to the detachment of β-strand A from each I27 domain. The force at which the humps appear was ∼130 pN and showed no temperature dependence, at least in the temperature range of 2°C-30°C. The hump structure was successfully analyzed with a two-state worm-like chain model, and the Gibbs free energy difference of the detachment reaction was estimated to be 11.6 ± 0.58 kcal/mol and found to be temperature independent. By contrast, upon lowering the temperature, the mean unfolding force from the partly unfolded intermediate state was found to markedly increase and the unfolding force distribution to broaden significantly, suggesting that the distance (x_u) between the folded and transition states in the energy landscape along the pulling direction is decreased. These results suggest that the local structure of β-strand A are stabilized by topologically simple local hydrogen-bond network and that the temperature does not affect the detachment reaction thermodynamically and kinetically, whereas the interaction between the β-strands A′ and G, which is a critical region for its mechanical stability, is strongly dependent on the temperature.     

Complex Stability of Single Proteins Explored by Forced Unfolding Experiments

In the last decade atomic force microscopy has been used to measure the mechanical stability of single proteins. These force spectroscopy experiments have shown that many water-soluble and membrane proteins unfold via one or more intermediates. Recently, Li and co-workers found a linear correlation between the unfolding force of the native state and the intermediate in fibronectin, which they suggested indicated the presence of a molecular memory or multiple unfolding pathways (1). Here, we apply two independent methods in combination with Monte Carlo simulations to analyze the unfolding of α-helices E and D of bacteriorhodopsin (BR). We show that correlation analysis of unfolding forces is very sensitive to errors in force calibration of the instrument. In contrast, a comparison of relative forces provides a robust measure for the stability of unfolding intermediates. The proposed approach detects three energetically different states of α-helices E and D in trimeric BR. These states are not observed for monomeric BR and indicate that substantial information is hidden in forced unfolding experiments of single proteins.     

Connecting Thermal and Mechanical Protein (Un)folding Landscapes

Molecular dynamics simulations supplement single-molecule pulling experiments by providing the possibility of examining the full free energy landscape using many coordinates. Here, we use an all-atom structure-based model to study the force and temperature dependence of the unfolding of the protein filamin by applying force at both termini. The unfolding time-force relation τ(F) indicates that the force-induced unfolding behavior of filamin can be characterized into three regimes: barrier-limited low- and intermediate-force regimes, and a barrierless high-force regime. Slope changes of τ(F) separate the three regimes. We show that the behavior of τ(F) can be understood from a two-dimensional free energy landscape projected onto the extension X and the fraction of native contacts Q. In the low-force regime, the unfolding rate is roughly force-independent due to the small (even negative) separation in X between the native ensemble and transition state ensemble (TSE). In the intermediate-force regime, force sufficiently separates the TSE from the native ensemble such that τ(F) roughly follows an exponential relation. This regime is typically explored by pulling experiments. While X may fail to resolve the TSE due to overlap with the unfolded ensemble just below the folding temperature, the overlap is minimal at lower temperatures where experiments are likely to be conducted. The TSE becomes increasingly structured with force, whereas the average order of structural events during unfolding remains roughly unchanged. The high-force regime is characterized by barrierless unfolding, and the unfolding time approaches a limit of ∼10 μs for the highest forces we studied. Finally, a combination of X and Q is shown to be a good reaction coordinate for almost the entire force range.     

Coarse-Grained Simulations of Topology-Dependent Mechanisms of Protein Unfolding and Translocation Mediated by ClpY ATPase Nanomachines

Clp ATPases are powerful ring shaped nanomachines which participate in the degradation pathway of the protein quality control system, coupling the energy from ATP hydrolysis to threading substrate proteins (SP) through their narrow central pore. Repetitive cycles of sequential intra-ring ATP hydrolysis events induce axial excursions of diaphragm-forming central pore loops that effect the application of mechanical forces onto SPs to promote unfolding and translocation. We perform Langevin dynamics simulations of a coarse-grained model of the ClpY ATPase-SP system to elucidate the molecular details of unfolding and translocation of an α/β model protein. We contrast this mechanism with our previous studies which used an all-α SP. We find conserved aspects of unfolding and translocation mechanisms by allosteric ClpY, including unfolding initiated at the tagged C-terminus and translocation via a power stroke mechanism. Topology-specific aspects include the time scales, the rate limiting steps in the degradation pathway, the effect of force directionality, and the translocase efficacy. Mechanisms of ClpY-assisted unfolding and translocation are distinct from those resulting from non-allosteric mechanical pulling. Bulk unfolding simulations, which mimic Atomic Force Microscopy-type pulling, reveal multiple unfolding pathways initiated at the C-terminus, N-terminus, or simultaneously from both termini. In a non-allosteric ClpY ATPase pore, mechanical pulling with constant velocity yields larger effective forces for SP unfolding, while pulling with constant force results in simultaneous unfolding and translocation.     

Repetitive pulling catalyzes co-translocational unfolding of barnase during import through a mitochondrial pore

We present a computational study of barnase unfolding during import into mitochondria through a model translocon. In contrast to thermal (or chemical) unfolding, the major intermediates of co-translocational unfolding are mainly mediated by non-native interactions accompanying the protein configurations induced by pulling forces. These energy contributions, combined with backbone topological constraints imposed by the model pore, result in milestones along the unfolding pathways which are significantly different not only from those experienced during thermal (or chemical) denaturation, but also from those observed in single-molecule pulling by both ends without pore constraints. Two on-pathway major translocation intermediates trapped in long-lived states by significantly high unfolding barriers are identified. A fraction of these pathways can, however, skip such local kinetic traps and result in extremely fast translocations, leading to a dramatic kinetic partitioning spanning approximately four orders of magnitude. The fraction of fast translocation events is shown to increase upon switching the pull off and on, when compared to pulling at constant force. This suggests a "catalytic" mechanism by which the mitochondrial import machinery regulates this partitioning by repetitively pulling in cycles. A number of mutation sites that alter the kinetic "flow" of the unfolding trajectories are suggested and tested.     

Dynameomics: A Consensus View of the Protein Unfolding/Folding Transition State Ensemble across a Diverse Set of Protein Folds

The Dynameomics project aims to simulate a representative sample of all globular protein metafolds under both native and unfolding conditions. We have identified protein unfolding transition state (TS) ensembles from multiple molecular dynamics simulations of high-temperature unfolding in 183 structurally distinct proteins. These data can be used to study individual proteins and individual protein metafolds and to mine for TS structural features common across all proteins. Separating the TS structures into four different fold classes (all proteins, all-α, all-β, and mixed α/β and α + β) resulted in no significant difference in the overall protein properties. The residues with the most contacts in the native state lost the most contacts in the TS ensemble. On average, residues beginning in an α-helix maintained more structure in the TS ensemble than did residues starting in β-strands or any other conformation. The metafolds studied here represent 67% of all known protein structures, and this is, to our knowledge, the largest, most comprehensive study of the protein folding/unfolding TS ensemble to date. One might have expected broad distributions in the average global properties of the TS relative to the native state, indicating variability in the amount of structure present in the TS. Instead, the average global properties converged with low standard deviations across metafolds, suggesting that there are general rules governing the structure and properties of the TS.     

Delineating the conformational dynamics of intermediate structures on the unfolding pathway of β-lactoglobulin in aqueous urea and dimethyl sulfoxide

Abstract

The funnel shaped energy landscape model of the protein folding suggests that progression of folding proceeds through multiple pathways, having the multiple intermediates which leads to multidimensional free-energy surface. Herein, we applied all-atom MD simulation to conduct a comparative study on the structure of β-lactoglobulin (β-LgA) in aqueous mixture of 8?M urea and 8?M dimethyl sulfoxide (DMSO), at different temperatures. The cumulative results of multiple simulations suggest a common unfolding pathway of β-LgA, occurred through the stable and meta-stable intermediates (I), in both urea and DMSO. However, the free-energy landscape (FEL) analyses show that the structural transitions of I-states are energetically different. In urea, FEL shows distinct ensemble of intermediates, I1 and I2, separated by the energy barrier of ～3.0?kcal mol^?1. Similarly, we find the population of two distinct I1 and I2 states in DMSO, however, the I1 appeared transiently around ～30–35?ns and is short-lived. But, the I2 ensemble is observed structurally compact and long-lived (～50–150?ns) as compared to unfolding in urea. Furthermore, the I1 and I2 are separated through a high energy barrier of ～6.0?kcal mol^?1. Thus, our results provide the structural insights of intermediates which essentially bear the signature of a different unfolding pathway of β-LgA in urea and DMSO.

Abbreviations β-LgA β-lactoglobulin

DMSO dimethyl sulfoxide

FEL free-energy landscape

GdmCl guanidinium chloride

I intermediate state

MG molten globule state

PME particle mesh Ewald

Q fraction of native contacts

RMSD root mean square deviation

RMSF root mean square fluctuation

R_g radius of gyration

SASA solvent Accessible Surface Area

scSASA the side chain SASA

Trp tryptophan

Communicated by Ramaswamy H. Sarma     

Forced-unfolding and force-quench refolding of RNA hairpins

Nanomanipulation of individual RNA molecules, using laser optical tweezers, has made it possible to infer the major features of their energy landscape. Time-dependent mechanical unfolding trajectories, measured at a constant stretching force (f_S) of simple RNA structures (hairpins and three-helix junctions) sandwiched between RNA/DNA hybrid handles show that they unfold in a reversible all-or-none manner. To provide a molecular interpretation of the experiments we use a general coarse-grained off-lattice Gō-like model, in which each nucleotide is represented using three interaction sites. Using the coarse-grained model we have explored forced-unfolding of RNA hairpin as a function of f_S and the loading rate (r_f). The simulations and theoretical analysis have been done both with and without the handles that are explicitly modeled by semiflexible polymer chains. The mechanisms and timescales for denaturation by temperature jump and mechanical unfolding are vastly different. The directed perturbation of the native state by f_S results in a sequential unfolding of the hairpin starting from their ends, whereas thermal denaturation occurs stochastically. From the dependence of the unfolding rates on r_f and f_S we show that the position of the unfolding transition state is not a constant but moves dramatically as either r_f or f_S is changed. The transition-state movements are interpreted by adopting the Hammond postulate for forced-unfolding. Forced-unfolding simulations of RNA, with handles attached to the two ends, show that the value of the unfolding force increases (especially at high pulling speeds) as the length of the handles increases. The pathways for refolding of RNA from stretched initial conformation, upon quenching f_S to the quench force f_Q, are highly heterogeneous. The refolding times, upon force-quench, are at least an order-of-magnitude greater than those obtained by temperature-quench. The long f_Q-dependent refolding times starting from fully stretched states are analyzed using a model that accounts for the microscopic steps in the rate-limiting step, which involves the trans to gauche transitions of the dihedral angles in the GAAA tetraloop. The simulations with explicit molecular model for the handles show that the dynamics of force-quench refolding is strongly dependent on the interplay of their contour length and persistence length and the RNA persistence length. Using the generality of our results, we also make a number of precise experimentally testable predictions.     

Mechanical unfolding of RNA: from hairpins to structures with internal multiloops

Mechanical unfolding of RNA structures, ranging from hairpins to ribozymes, using laser optical tweezer experiments have begun to reveal the features of the energy landscape that cannot be easily explored using conventional experiments. Upon application of constant force (f), RNA hairpins undergo cooperative transitions from folded to unfolded states whereas subdomains of ribozymes unravel one at a time. Here, we use a self-organized polymer model and Brownian dynamics simulations to probe mechanical unfolding at constant force and constant-loading rate of four RNA structures of varying complexity. For simple hairpins, such as P5GA, application of constant force or constant loading rate results in bistable cooperative transitions between folded and unfolded states without populating any intermediates. The transition state location (DeltaxFTS) changes dramatically as the loading rate is varied. At loading rates comparable to those used in laser optical tweezer experiments, the hairpin is plastic, with DeltaxFTS being midway between folded and unfolded states; whereas at high loading rates, DeltaxFTS moves close to the folded state, i.e., RNA is brittle. For the 29-nucleotide TAR RNA with the three-nucleotide bulge, unfolding occurs in a nearly two-state manner with an occasional pause in a high free energy metastable state. Forced unfolding of the 55 nucleotides of the Hepatitis IRES domain IIa, which has a distorted L-shaped structure, results in well-populated stable intermediates. The most stable force-stabilized intermediate represents straightening of the L-shaped structure. For these structures, the unfolding pathways can be predicted using the contact map of the native structures. Unfolding of a RNA motif with internal multiloop, namely, the 109-nucleotide prohead RNA that is part of the 29 DNA packaging motor, at constant value of rf occurs with three distinct rips that represent unraveling of the paired helices. The rips represent kinetic barriers to unfolding. Our work shows 1), the response of RNA to force is largely determined by the native structure; and 2), only by probing mechanical unfolding over a wide range of forces can the underlying energy landscape be fully explored.     

Mechanical Unfolding of Acylphosphatase Studied by Single-Molecule Force Spectroscopy and MD Simulations

Single-molecule manipulation methods provide a powerful means to study protein transitions. Here we combined single-molecule force spectroscopy and steered molecular-dynamics simulations to study the mechanical properties and unfolding behavior of the small enzyme acylphosphatase (AcP). We find that mechanical unfolding of AcP occurs at relatively low forces in an all-or-none fashion and is decelerated in the presence of a ligand, as observed in solution measurements. The prominent energy barrier for the transition is separated from the native state by a distance that is unusually long for α/β proteins. Unfolding is initiated at the C-terminal strand (β_T) that lies at one edge of the β-sheet of AcP, followed by unraveling of the strand located at the other. The central strand of the sheet and the two helices in the protein unfold last. Ligand binding counteracts unfolding by stabilizing contacts between an arginine residue (Arg-23) and the catalytic loop, as well as with β_T of AcP, which renders the force-bearing units of the protein resistant to force. This stabilizing effect may also account for the decelerated unfolding of ligand-bound AcP in the absence of force.     

Force-Regulated Refolding of the Mechanosensory Domain in the Platelet Glycoprotein Ib-IX Complex

In platelets, the glycoprotein (GP) Ib-IX receptor complex senses blood shear flow and transmits the mechanical signals into platelets. Recently, we have discovered a juxtamembrane mechanosensory domain (MSD) within the GPIbα subunit of GPIb-IX. Mechanical unfolding of the MSD activates GPIb-IX signaling into platelets, leading to their activation and clearance. Using optical tweezer-based single-molecule force measurement, we herein report a systematic biomechanical characterization of the MSD in its native, full-length receptor complex and a recombinant, unglycosylated MSD in isolation. The native MSD unfolds at a resting rate of 9 × 10^?3 s^?1. Upon exposure to pulling forces, MSD unfolding accelerates exponentially over a force scale of 2.0 pN. Importantly, the unfolded MSD can refold with or without applied forces. The unstressed refolding rate of MSD is ～17 s^?1 and slows exponentially over a force scale of 3.7 pN. Our measurements confirm that the MSD is relatively unstable, with a folding free energy of 7.5 k_BT. Because MSD refolding may turn off GPIb-IX’s mechanosensory signals, our results provide a mechanism for the requirement of a continuous pulling force of >15 pN to fully activate GPIb-IX.     

Complex stability of single proteins explored by forced unfolding experiments

In the last decade atomic force microscopy has been used to measure the mechanical stability of single proteins. These force spectroscopy experiments have shown that many water-soluble and membrane proteins unfold via one or more intermediates. Recently, Li and co-workers found a linear correlation between the unfolding force of the native state and the intermediate in fibronectin, which they suggested indicated the presence of a molecular memory or multiple unfolding pathways (1). Here, we apply two independent methods in combination with Monte Carlo simulations to analyze the unfolding of alpha-helices E and D of bacteriorhodopsin (BR). We show that correlation analysis of unfolding forces is very sensitive to errors in force calibration of the instrument. In contrast, a comparison of relative forces provides a robust measure for the stability of unfolding intermediates. The proposed approach detects three energetically different states of alpha-helices E and D in trimeric BR. These states are not observed for monomeric BR and indicate that substantial information is hidden in forced unfolding experiments of single proteins.     

Investigating the refolding pathway of human acidic fibroblast growth factor (hFGF-1) from the residual structure(s) obtained by denatured-state hydrogen/deuterium exchange

Human fibroblast growth factor 1 (hFGF-1) consists of 12 anti-parallel β-strands arranged into a β-trefoil architecture. We directly measured hydrogen/deuterium exchange rates on the urea-denatured hFGF-1 to obtain the information with regard to the persistent residual interaction(s) in the unfolded hFGF-1. Thirty-eight residues whose heteronuclear single quantum coherence cross-peaks can be observed after exchange show higher protections than those predicted for the same residues in a random coil conformation, suggesting the existence of residual structure(s). The urea-denaturation of hFGF-1 tested by both circular dichroism and fluorescence spectroscopy indicated that the unfolding process is a cooperative two-state process and that the residual structures observed did not originate from the existence of a partially structured intermediate. The coincident disappearance of the native heteronuclear single quantum coherence cross-peaks during the urea-denaturation process suggests that the residual structures observed contain no nativelike interactions. The protected residues (fold ons) in the urea-denatured state are mostly those that exchange slowly in the native state H/D exchange. The distribution of these fold ons in the native structure of hFGF-1 suggests that the refolding starts by collisions between the residual structures (microdomains) between the β-strands VI and VII, and between the β-strands II and III, which appear to be two independent refolding coordinates during the refolding process.     

Membrane-Protein Unfolding Intermediates Detected with Enhanced Precision Using a Zigzag Force Ramp

Precise quantification of the energetics and interactions that stabilize membrane proteins in a lipid bilayer is a long-sought goal. Toward this end, atomic force microscopy has been used to unfold individual membrane proteins embedded in their native lipid bilayer, typically by retracting the cantilever at a constant velocity. Recently, unfolding intermediates separated by as few as two amino acids were detected using focused-ion-beam-modified ultrashort cantilevers. However, unambiguously discriminating between such closely spaced states remains challenging, in part because any individual unfolding trajectory only occupies a subset of the total number of intermediates. Moreover, structural assignment of these intermediates via worm-like-chain analysis is hindered by brief dwell times compounded with thermal and instrumental noise. To overcome these issues, we moved the cantilever in a sawtooth pattern of 6–12 nm, offset by 0.25–1 nm per cycle, generating a “zigzag” force ramp of alternating positive and negative loading rates. We applied this protocol to the model membrane protein bacteriorhodopsin (bR). In contrast to conventional studies that extract bR’s photoactive retinal along with the first transmembrane helix, we unfolded bR in the presence of its retinal. To do so, we introduced a previously developed enzymatic-cleavage site between helices E and F and pulled from the top of the E helix using a site-specific, covalent attachment. The resulting zigzag unfolding trajectories occupied 40% more states per trajectory and occupied those states for longer times than traditional constant-velocity records. In total, we identified 31 intermediates during the unfolding of five helices of EF-cleaved bR. These included a previously reported, mechanically robust intermediate located between helices C and B that, with our enhanced resolution, is now shown to be two distinct states separated by three amino acids. Interestingly, another intermediate directly interacted with the retinal, an interaction confirmed by removing the retinal.     

Two Amyloid States of the Prion Protein Display Significantly Different Folding Patterns

It has been well established that a single amino acid sequence can give rise to several conformationally distinct amyloid states. The extent to which amyloid structures formed within the same sequence are different, however, remains unclear. To address this question, we studied two amyloid states (referred to as R- and S-fibrils) produced in vitro from highly purified full-length recombinant prion protein. Several biophysical techniques including X-ray diffraction, CD, Fourier transform infrared spectroscopy (FTIR), hydrogen-deuterium exchange, proteinase K digestion, and binding of a conformation-sensitive fluorescence dye revealed that R- and S-fibrils have substantially different secondary, tertiary, and quaternary structures. While both states displayed a 4. 8-Å meridional X-ray diffraction typical for amyloid cross-β-spines, they showed markedly different equatorial profiles, suggesting different folding pattern of β-strands. The experiments on hydrogen-deuterium exchange monitored by FTIR revealed that only small fractions of amide protons were protected in R- or S-fibrils, an argument for the dynamic nature of their cross-β-structure. Despite this fact, both amyloid states were found to be very stable conformationally as judged from temperature-induced denaturation monitored by FTIR and the conformation-sensitive dye. Upon heating to 80 °C, only local unfolding was revealed, while individual state-specific cross-β features were preserved. The current studies demonstrated that the two amyloid states formed by the same amino acid sequence exhibited significantly different folding patterns that presumably reflect two different architectures of cross-β-structure. Both S- and R-fibrils, however, shared high conformational stability, arguing that the energy landscape for protein folding and aggregation can contain several deep free-energy minima.     

(Un)Folding Mechanisms of the FBP28 WW Domain in Explicit Solvent Revealed by Multiple Rare Event Simulation Methods

We report a numerical study of the (un)folding routes of the truncated FBP28 WW domain at ambient conditions using a combination of four advanced rare event molecular simulation techniques. We explore the free energy landscape of the native state, the unfolded state, and possible intermediates, with replica exchange molecular dynamics. Subsequent application of bias-exchange metadynamics yields three tentative unfolding pathways at room temperature. Using these paths to initiate a transition path sampling simulation reveals the existence of two major folding routes, differing in the formation order of the two main hairpins, and in hydrophobic side-chain interactions. Having established that the hairpin strand separation distances can act as reasonable reaction coordinates, we employ metadynamics to compute the unfolding barriers and find that the barrier with the lowest free energy corresponds with the most likely pathway found by transition path sampling. The unfolding barrier at 300 K is ∼17 k_BT ≈ 42 kJ/mol, in agreement with the experimental unfolding rate constant. This work shows that combining several powerful simulation techniques provides a more complete understanding of the kinetic mechanism of protein folding.     

Kinetic traps in the folding of beta alpha-repeat proteins: CheY initially misfolds before accessing the native conformation

The βα-repeat class of proteins, represented by the (βα)₈ barrel and the α/β/α sandwich, are among the most common structural platforms in biology. Previous studies on the folding mechanisms of these motifs have revealed or suggested that the initial event involves the submillisecond formation of a kinetically trapped species that must at least partially unfold before productive folding to the respective native conformation can occur. To test the generality of these observations, CheY, a bacterial response regulator, was subjected to an extensive analysis of its folding reactions. Although earlier studies had proposed the formation of an off-pathway intermediate, the data available were not sufficient to rule out an alternative on-pathway mechanism. A global analysis of single- and double-jump kinetic data, combined with equilibrium unfolding data, was used to show that CheY folds and unfolds through two parallel channels defined by the state of isomerization of a prolyl peptide bond in the active site. Each channel involves a stable, highly structured folding intermediate whose kinetic properties are better described as the properties of an off-pathway species. Both intermediates subsequently flow through the unfolded state ensemble and adopt the native cis-prolyl isomer prior to forming the native state. Initial collapse to off-pathway folding intermediates is a common feature of the folding mechanisms of βα-repeat proteins, perhaps reflecting the favored partitioning to locally determined substructures that cannot directly access the native conformation. Productive folding requires the dissipation of these prematurely folded substructures as a prelude to forming the larger-scale transition state that leads to the native conformation. Results from Gō-modeling studies in the accompanying paper elaborate on the topological frustration in the folding free-energy landscape of CheY.     

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.