Backbone Fluctuations Research Articles

Arginine-phosphate salt bridges between histones and DNA: intermolecular actuators that control nucleosome architecture.

Structural bioinformatics and van der Waals density functional theory are combined to investigate the mechanochemical impact of a major class of histone-DNA interactions, namely, the formation of salt bridges between arginine residues in histones and phosphate groups on the DNA backbone. Principal component analysis reveals that the configurational fluctuations of the sugar-phosphate backbone display sequence-specific directionality and variability, and clustering of nucleosome crystal structures identifies two major salt-bridge configurations: a monodentate form in which the arginine end-group guanidinium only forms one hydrogen bond with the phosphate, and a bidentate form in which it forms two. Density functional theory calculations highlight that the combination of sequence, denticity, and salt-bridge positioning enables the histones to apply a tunable mechanochemical stress to the DNA via precise and specific activation of backbone deformations. The results suggest that selection for specific placements of van der Waals contacts, with high-precision control of the spatial distribution of intermolecular forces, may serve as an underlying evolutionary design principle for the structure and function of nucleosomes, a conjecture that is corroborated by previous experimental studies.

PH dependence of conformational fluctuations of the protein backbone.

Proton binding equilibria (pK(a) values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK(a) values measured for individual ionizable residues with NMR spectroscopy are actually population-weighted averages of the pK(a) in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK(a) values of surface Glu and Asp residues are affected by pH-sensitive fluctuations of the backbone between folded and locally unfolded conformations. (15)N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps-ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub-global unfolding near carboxylic groups. This study has implications for structure-based pKa calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH-sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co-operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK(a) values of ionizable residues in proteins.

Probabilistic Determination of Native State Ensembles of Proteins

The motions of biological macromolecules are tightly coupled to their functions. However, while the study of fast motions has become increasingly feasible in recent years, the study of slower, biologically important motions remains difficult. Here, we present a method to construct native state ensembles of proteins by the combination of physical force fields and experimental data through modern statistical methodology. As an example, we use NMR residual dipolar couplings to determine a native state ensemble of the extensively studied third immunoglobulin binding domain of protein G (GB3). The ensemble accurately describes both local and nonlocal backbone fluctuations as judged by its reproduction of complementary experimental data. While it is difficult to assess precise time-scales of the observed motions, our results suggest that it is possible to construct realistic conformational ensembles of biomolecules very efficiently. The approach may allow for a dramatic reduction in the computational as well as experimental resources needed to obtain accurate conformational ensembles of biological macromolecules in a statistically sound manner.

Internally labeled Cy3/Cy5 DNA constructs show greatly enhanced photo-stability in single-molecule FRET experiments

DNA constructs labeled with cyanine fluorescent dyes are important substrates for single-molecule (sm) studies of the functional activity of protein–DNA complexes. We previously studied the local DNA backbone fluctuations of replication fork and primer–template DNA constructs labeled with Cy3/Cy5 donor–acceptor Förster resonance energy transfer (FRET) chromophore pairs and showed that, contrary to dyes linked ‘externally’ to the bases with flexible tethers, direct ‘internal’ (and rigid) insertion of the chromophores into the sugar-phosphate backbones resulted in DNA constructs that could be used to study intrinsic and protein-induced DNA backbone fluctuations by both smFRET and sm Fluorescent Linear Dichroism (smFLD). Here we show that these rigidly inserted Cy3/Cy5 chromophores also exhibit two additional useful properties, showing both high photo-stability and minimal effects on the local thermodynamic stability of the DNA constructs. The increased photo-stability of the internal labels significantly reduces the proportion of false positive smFRET conversion ‘background’ signals, thereby simplifying interpretations of both smFRET and smFLD experiments, while the decreased effects of the internal probes on local thermodynamic stability also make fluctuations sensed by these probes more representative of the unperturbed DNA structure. We suggest that internal probe labeling may be useful in studies of many DNA–protein interaction systems.

Conformational instability of human prion protein upon residue modification: a molecular dynamics simulation study.

Technical strategies like amino acid substitution and residue modification have been widely used to characterize the importance of key amino acids and the role that each residue plays in the structural and functional properties of protein molecules. However, there is no systematic approach to assess the impact of the substituted/modified amino acids on the conformational dynamics of proteins. In this investigation to clarify the effects of residue modifications on the structural dynamics of human prion protein (PrP), a comparative molecular dynamics simulation study on the native and the amino acid-substituted analog at position 208 of PrP has been performed. It is believed that Arginine to Histidine mutation at position 208 is responsible for the structural transition of the native form of human prion protein to the pathogenic isoform causing Creutzfeldt-Jakob disease (CJD). So, three 10 ns molecular dynamics simulations on three model constructs have been performed. Simulation results indicated considerable differences of conformational fluctuations for Alanine substituted construct (PrPALA) and the analog form (PrPSB) comprising the neutralized state of the Arginine residue at position 208 of the human prion protein. According to our data, substitution of the Arginine residue by the uncharged state of this residue induces some reversible structural alterations in the intrinsically flexible loop area including residues 167-171 of PrP. Thus, deprotonation of Arg(208) is a weak perturbation to the structural fluctuations of the protein backbone and the resulting construct behaves almost identical as its native form. Otherwise, Alanine substitution at position 208 imposed an irreversible impact on the secondary and tertiary structure of the protein, which leads to conformational instabilities in the remote hot region comprising residues 190-195 of the C-terminal part of helix 2. Based on the results, it could be deduced that the observed conformational transitions upon Arg(208) to His point mutation, which is the main reason for CJD, may be mainly related to the structural instabilities due to the induced-conformational changes that caused alterations in local/spatial arrangements of the force distributions in the backbone of the human prion protein.

The role of topology and thermal backbone fluctuations on sacrificial bond efficacy in mechanical metalloproteins

Sacrificial bonding is a ubiquitous cross-linking strategy for increasing toughness that is found throughout nature in various biological materials such as bone, wood, silk and mussel byssal threads. However, the molecular mechanism of sacrificial bonding remains only poorly understood. Molecular modeling possesses a strong potential to provide insights into the behavior of these cross-links. Here we use Monte Carlo simulations to investigate the effect of reversible sacrificial binding sites on the mechanical properties of single linear polymer chains based on load-bearing metalloproteins found in the mussel byssus. It is shown that the topology of the bonds determines the position and spacing of sacrificial force peaks, while the height of these peaks is intimately tied to the magnitude of thermal fluctuations in the chain that are dependent on effective chain length. These results bear important implications for understanding natural systems and for the generation of strong and ductile biomimetic polymers.

Technological advances in site-directed spin labeling of proteins

Molecular flexibility over a wide time range is of central importance to the function of many proteins, both soluble and membrane. Revealing the modes of flexibility, their amplitudes, and time scales under physiological conditions is the challenge for spectroscopic methods, one of which is site-directed spin labeling EPR (SDSL-EPR). Here we provide an overview of some recent technological advances in SDSL-EPR related to investigation of structure, structural heterogeneity, and dynamics of proteins. These include new classes of spin labels, advances in measurement of long range distances and distance distributions, methods for identifying backbone and conformational fluctuations, and new strategies for determining the kinetics of protein motion.

On the role of thermal backbone fluctuations in myoglobin ligand gate dynamics

We construct an energy function that describes the crystallographic structure of sperm whale myoglobin backbone. As a model in our construction, we use the Protein Data Bank entry 1ABS that has been measured at liquid helium temperature. Consequently, the thermal B-factor fluctuations are very small, which is an advantage in our construction. The energy function that we utilize resembles that of the discrete nonlinear Schrödinger equation. Likewise, ours supports topological solitons as local minimum energy configurations. We describe the 1ABS backbone in terms of topological solitons with a precision that deviates from 1ABS by an average root-mean-square distance, which is less than the experimentally observed Debye-Waller B-factor fluctuation distance. We then subject the topological multi-soliton solution to extensive numerical heating and cooling experiments, over a very wide range of temperatures. We concentrate in particular to temperatures above 300 K and below the Θ-point unfolding temperature, which is around 348 K. We confirm that the behavior of the topological multi-soliton is fully consistent with Anfinsen's thermodynamic principle, up to very high temperatures. We observe that the structure responds to an increase of temperature consistently in a very similar manner. This enables us to characterize the onset of thermally induced conformational changes in terms of three distinct backbone ligand gates. One of the gates is made of the helix F and the helix E. The two other gates are chosen similarly, when open they provide a direct access route for a ligand to reach the heme. We find that out of the three gates we investigate, the one which is formed by helices B and G is the most sensitive to thermally induced conformational changes. Our approach provides a novel perspective to the important problem of ligand entry and exit.

Structure and dynamics of an imidazoline nitroxide side chain with strongly hindered internal motion in proteins

A disulfide-linked imidazoline nitroxide side chain (V1) has a similar and highly constrained internal motion at diverse topological sites in a protein, unlike that for the disulfide-linked pyrroline nitroxide side chain (R1) widely used in site directed spin labeling EPR. Crystal structures of V1 at two positions in a helix of T4 Lysozyme and quantum mechanical calculations suggest the source of the constraints as intra-side chain interactions of the disulfide sulfur atoms with both the protein backbone and the 3-nitrogen in the imidazoline ring. These interactions apparently limit the conformation of the side chain to one of only three possible rotamers, two of which are observed in the crystal structure. An inter-spin distance measurement in frozen solution using double electron–electron resonance (DEER) gives a value essentially identical to that determined from the crystal structure of the protein containing two copies of V1, indicating that lattice forces do not dictate the rotamers observed. Collectively, the results suggest the possibility of predetermining a unique rotamer of V1 in helical structures. In general, the reduced rotameric space of V1 compared to R1 should simplify interpretation of inter-spin distance information in terms of protein structure, while the highly constrained internal motion is expected to extend the dynamic range for characterizing large amplitude nanosecond backbone fluctuations.

Backbone Conformational Fluctuations as Determinants of pKa Values of Surface Residues in Proteins

The pKa values of surface ionizable residues in proteins are usually similar to the pKa values of ionizable residues in water. This is still difficult to reproduce using structure-based electrostatics calculations with static structures, which tend to exaggerate the shifts in pKa values because the magnitude of Coulomb effects is exaggerated and the balance between Coulomb and dehydration effects is not reproduced correctly. These problems can be minimized by using artificially high protein dielectric constants, but sometimes the discrepancy persists even when the protein interior is treated with the dielectric constant of water. This suggests that pKa values are affected by the conformational dynamics of the protein, which are not reflected in the crystal structure. Molecular dynamics or Monte Carlo simulations can be used to attempt to reproduce these dynamic effects implicitly, but there are no useful data for testing this approach directly. To examine the role of fluctuations of the backbone in determining pKa values, NMR spectroscopy was used to measure the pKa values of all 20 Asp and Glu residues in variants of staphylococcal nuclease with Gly substituted at select locations. These substitutions were intended to enhance backbone fluctuations without affecting the overall protein structure. Some Gly substitutions detectably shifted the pKa values of some carboxylic groups. Crystal structures of these variants showed no significant changes relative to the structure of the reference protein. Calculations with continuum methods using these crystal structures do not reproduce the measured shifts in pKa values. NMR spectroscopy studies suggest Gly substitutions affect the dynamics of the protein. Efforts are underway to better characterize these changes. Our results suggest that the high apparent polarizability of proteins might be due to subtle structural fluctuations that are difficult to reproduce computationally.

Impact of deglycosylation and thermal stress on conformational stability of a full length murine igG2a monoclonal antibody: Observations from molecular dynamics simulations

With the rise of antibody based therapeutics as successful medicines, there is an emerging need to understand the fundamental antibody conformational dynamics and its implications towards stability of these medicines. Both deglycosylation and thermal stress have been shown to cause conformational destabilization and aggregation in monoclonal antibodies. Here, we study instabilities caused by deglycosylation and by elevated temperature (400 K) by performing molecular dynamic simulations on a full length murine IgG2a mAb whose crystal structure is available in the Protein Data bank. C(α)-atom root mean square deviation and backbone root mean square fluctuation calculations show that deglycosylation perturbs quaternary and tertiary structures in the C(H) 2 domains. In contrast, thermal stress pervades throughout the antibody structure and both Fabs and Fc regions are destabilized. The thermal stress applied in this study was not sufficient to cause large scale unfolding within the simulation time and most amino acid residues showed similar average solvent accessible surface area and secondary structural conformations in all trajectories. C(H) 3 domains were the most successful at resisting the conformational destabilization. The simulations helped identify aggregation prone regions, which may initiate cross-β motif formation upon deglycosylation and upon applying thermal stress. Deglycosylation leads to increased backbone fluctuations and solvent exposure of a highly conserved APR located in the edge β-strand A of the C(H) 2 domains. Aggregation upon thermal stress is most likely initiated by two APRs that overlap with the complementarity determining regions. This study has important implications for rational design of antibody based therapeutics that are resistant towards aggregation.

Probing Sequence-Specific DNA Flexibility in A-Tracts and Pyrimidine-Purine Steps by Nuclear Magnetic Resonance 13C Relaxation and Molecular Dynamics Simulations

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.

Exploring structural and optical properties of fluorescent proteins by squeezing: modeling high-pressure effects on the mStrawberry and mCherry red fluorescent proteins.

Molecular dynamics calculations of pressure effects on mStrawberry and mCherry fluorescent proteins are reported. The simulations reveal that mStrawberry has much floppier structure at atmospheric pressure, as evidenced by larger backbone fluctuations and the coexistence of two conformers that differ by Ser146 orientation. Consequently, pressure increase has a larger effect on mStrawberry, making its structure more rigid and reducing the population of one of the conformers. The most significant effect of pressure increase is in the hydrogen-bonding network between the chromophore and the nearby residues. The quantum-mechanics/molecular mechanics calculations of excitation energies in mStrawberry explain the observed blue shift and identify Lys70 as the residue that has the most pronounced effect on the spectra. The results suggest that pressure increase causes an initial increase of fluorescence yield only for relatively floppy fluorescent proteins, whereas the fluorescent proteins that have more rigid structures have quantum yields close to their maximum. The results suggest that a low quantum yield in fluorescent proteins is dynamic in nature and depends on the range of thermal motions of the chromophore and fluctuations in the H-bonding network rather than on their average structure.

Mapping Molecular Flexibility of Proteins with Site-Directed Spin Labeling: A Case Study of Myoglobin

Site-directed spin labeling (SDSL) has potential for mapping protein flexibility under physiological conditions. The purpose of the present study was to explore this potential using 38 singly spin-labeled mutants of myoglobin distributed throughout the sequence. Correlation of the EPR spectra with protein structure provides new evidence that the site-dependent variation in line shape, and hence motion of the spin label, is due largely to differences in mobility of the helical backbone in the ns time range. Fluctuations between conformational substates, typically in the μs-ms time range, are slow on the EPR time scale, and the spectra provide a snapshot of conformational equilibria frozen in time as revealed by multiple components in the spectra. A recent study showed that osmolyte perturbation can positively identify conformational exchange as the origin of multicomponent spectra (López et al. (2009), Protein Sci. 18, 1637). In the present study, this new strategy is employed in combination with line shape analysis and pulsed-EPR interspin distance measurements to investigate the conformation and flexibility of myoglobin in three folded and partially folded states. The regions identified to be in conformational exchange in the three forms agree remarkably well with those assigned by NMR, but the faster time scale of EPR allows characterization of localized states not detected in NMR. Collectively, the results suggest that SDSL-EPR and osmolyte perturbation provide a facile means for mapping the amplitude of fast backbone fluctuations and for detecting sequences in slow conformational exchange in folded and partially folded protein sequences.

Effect of Backbone Flexibility on Charge Transfer Rates in Peptide Nucleic Acid Duplexes

Charge transfer (CT) properties are compared between peptide nucleic acid structures with an aminoethylglycine backbone (aeg-PNA) and those with a γ-methylated backbone (γ-PNA). The common aeg-PNA is an achiral molecule with a flexible structure, whereas γ-PNA is a chiral molecule with a significantly more rigid structure than aeg-PNA. Electrochemical measurements show that the CT rate constant through an aeg-PNA bridging unit is twice the CT rate constant through a γ-PNA bridging unit. Theoretical calculations of PNA electronic properties, which are based on a molecular dynamics structural ensemble, reveal that the difference in the CT rate constant results from the difference in the extent of backbone fluctuations of aeg- and γ-PNA. In particular, fluctuations of the backbone affect the local electric field that broadens the energy levels of the PNA nucleobases. The greater flexibility of the aeg-PNA gives rise to more broadening, and a more frequent appearance of high-CT rate conformations than in γ-PNA.

A Delicate Interplay of Structure, Dynamics, and Thermodynamics for Function: A High Pressure NMR Study of Outer Surface Protein A

Outer surface protein A (OspA) is a crucial protein in the infection of Borrelia burgdorferi causing Lyme disease. We studied conformational fluctuations of OspA with high-pressure 15N/1H two-dimensional NMR along with high-pressure fluorescence spectroscopy. We found evidence within folded, native OspA for rapid local fluctuations of the polypeptide backbone in the nonglobular single layer β-sheet connecting the N- and C-terminal domains with τ << ms, which may give the two domains certain independence in mobility and thermodynamic stability. Furthermore, we found that folded, native OspA is in equilibrium (τ >> ms) with a minor conformer I, which is almost fully disordered and hydrated for the entire C-terminal part of the polypeptide chain from β8 to the C-terminus. Conformer I is characterized with ΔG0 = 32 ± 9 kJ/mol and ΔV0 = −140 ± 40 mL/mol, populating only ∼0.001% at 40°C at 0.1 MPa, pH 5.9. Because in the folded conformer the receptor binding epitope of OspA is buried in the C-terminal domain, its transition into conformer I under in vivo conditions may be critical for the infection of B. burgdorferi. The formation and stability of the peculiar conformer I are apparently supported by a large packing defect or cavity located in the C-terminal domain.

The High Dielectric Constant of Staphylococcal Nuclease Is Encoded in Its Structural Architecture

The pK(a) values of Lys-66, Glu-66, and Asp-66 buried in the interior of the staphylococcal nuclease Δ+PHS variant were reported to be shifted by as much as 5 pK(a) units from their normal values. Reproducing the pK(a) of these buried ionizable residues using continuum electrostatic calculations required the use of a high protein dielectric constant of 10 or higher. The apparent high dielectric constant has been rationalized as a consequence of a local structural reorganization or increased fluctuations in the microenvironment of the mutation site (Chimenti, M. S., et al. J. Mol. Biol. 2011, 405, 361-377). We have calculated the dielectric constant of Δ+PHS and the Lys-66, Asp-66, and Glu-66 mutants from first principles using the Kirkwood-Fröhlich equation and discovered that staphylococcal nuclease has a naturally high dielectric constant ranging from 20 to 30. This high dielectric constant does not change significantly with the mutation of residue 66 or with the ionization of the mutated residues. Calculation of the spatial dependence of the dielectric constant for the microenvironment of residue-66 produces a value of about 10, which matches well with the apparent dielectric constant needed to reproduce pK(a) values from continuum electrostatic calculations. Our results suggest an alternative explanation that the high dielectric constant of staphylococcal nuclease is a property resulting from the intrinsic backbone fluctuations originating from its structural architecture.

Dynamic structure of the polytheonamide B channel studied by normal mode analysis

Polytheonamide B is a linear 48-residue peptide with alternating d- and l-amino acids. It forms a single β6.3-helix structure and acts as a cation-selective ion channel. We performed normal mode analysis to investigate the dynamical properties of the polytheonamide B channel. We calculated root-mean-square (RMS) displacements of all backbone atoms, RMS fluctuations of the backbone dihedral angles φ and ψ, and correlation factors for the Cα atom and for dihedral angle fluctuations. The analysis revealed that polytheonamide B has elastic rod properties in very low-frequency modes and that the librational motions of backbone amide planes have the modes with relatively low frequencies, which is relevant to the function of polytheonamide B. The librational motions occur almost independently and are weakly anti-correlated with those of the neighbouring hydrogen-bonded amide planes. Backbone fluctuations show that polytheonamide B is roughly uniformly flexible over the entire helix. We also compared our results with those of gramicidin A, an analogue of polytheonamide B and found some common features in the dynamics of the ‘channel wall’. It was suggested that polytheonamide B utilises similar mechanisms to gramicidin A for conducting cations across the membrane.

Increased Hydrophobicity and Decreased Backbone Flexibility Explain the Lower Solubility of a Cataract-Linked Mutant of γD-Crystallin

A number of point mutations in γD-crystallin are associated with human cataract. The Pro23-to-Thr (P23T) mutation is perhaps the most common, is geographically widespread, and presents itself in a variety of phenotypes. It is therefore important to understand the molecular basis of lens opacity due to this mutation. In our earlier studies, we noted that P23T shows retrograde and sharply lowered solubility, most likely due to the emergence of hydrophobic patches involved in protein aggregation. Binding of 4,4′-dianilino-1,1′-binaphthyl-5,5′-disulfonate (Bis-ANS) dye (a probe commonly used for detecting surface hydrophobicity) competed with aggregation, suggesting that the residues involved in Bis-ANS binding are also involved in protein aggregation. Here, using NMR spectroscopy in conjunction with Bis-ANS binding, we identify three residues (Y16, D21, and Y50) in P23T that are involved in binding the dye. Furthermore, using 15N NMR relaxation experiments, we show that, in the mutant protein, backbone fluctuations are restricted to the picosecond-to-nanosecond and microsecond timescales relative to the wild type. Our present studies specify the residues involved in these two pivotal characteristics of the mutant protein, namely increased surface hydrophobicity and restricted mobility of the protein backbone, which can explain the nucleation and further propagation of protein aggregates. Thus, we have now identified the residues in the P23T mutant that give rise to novel hydrophobic surfaces, as well as those regions of the protein backbone where fluctuations in different timescales are restricted, providing a comprehensive understanding of how lens opacity could result from this mutation.

Binding Energetics of Ferredoxin‐NADP+ Reductase with Ferredoxin and Its Relation to Function

To obtain insight into the motional features of proteins for enzymatic function, we studied binding reactions between ferredoxin-NADP(+) reductase (FNR) and ferredoxin (Fd) using isothermal titration calorimetry and NMR-based magnetic relaxation and hydrogen/deuterium exchange (HD(ex)). Fd-FNR binding was accompanied by endothermic reactions and driven by the entropy gain. Component-wise analysis of the net entropy change revealed that increases in the conformational entropy of the Fd-FNR complex contributed largely to stabilizing the complex. Intriguingly, analyses of magnetic relaxation and HD(ex) rates with X-ray B factor implied that Fd binding led to both structural stiffening and softening of FNR. Enhanced FNR backbone fluctuations suggest favorable contributions to the net conformational entropy. Fd-bound FNR further showed that relatively large-scale motions of the C terminus, a gatekeeper for interactions of NADP(+) (H), were quenched in the closed form, thereby facilitating exit of NADP(+) (H). This can provide a first dynamic structure-based explanation for the negative cooperativity between Fd and NADP(+) (H) via FNR.