Large domain fluctuations on 50-ns timescale enable catalytic activity in phosphoglycerate kinase

Metastable States in the Hinge-Bending Landscape of an Enzyme in an Atomistic Cytoplasm Simulation

Many enzymes undergo major conformational changes to function in cells, particularly when they bind to more than one substrate. We quantify the large-amplitude hinge-bending landscape of human phosphoglycerate kinase (PGK) in a human cytoplasm. Approximately 70 μs of all-atom simulations, upon coarse graining, reveal three metastable states of PGK with different hinge angle distributions and additional substates. The "open" state was more populated than the "semi-open" or "closed" states. In addition to free energies and barriers within the landscape, we characterized the average transition state passage time of ≈0.3 μs and reversible substrate and product binding. Human PGK in a dilute solution simulation shows a transition directly from the open to closed states, in agreement with previous SAXS experiments, suggesting that the cell-like model environment promotes stability of the human PGK semi-open state. Yeast PGK also sampled three metastable states within the cytoplasm model, with the closed state favored in our simulation.

Signature of functional enzyme dynamics in quasielastic neutron scattering spectra: The case of phosphoglycerate kinase

We present an analysis of high-resolution quasi-elastic neutron scattering spectra of phosphoglycerate kinase which elucidates the influence of the enzymatic activity on the dynamics of the protein. We show that in the active state the inter-domain motions are amplified and the intra-domain asymptotic power-law relaxation ∝t-α is accelerated, with a reduced coefficient α. Employing an energy landscape picture of protein dynamics, this observation can be translated into a widening of the distribution of energy barriers separating conformational substates of the protein.

Interchain Hydrodynamic Interaction and Internal Friction of Polyelectrolytes

Polyelectrolytes (PE) are polymeric macromolecules in aqueous solutions characterized by their chain topology and intrinsic charge in a neutralizing fluid. Structure and dynamics are related to several characteristic screening length scales determined by electrostatic, excluded volume, and hydrodynamic interactions. We examine PE dynamics in dilute to semidilute conditions using dynamic light scattering, neutron spinecho spectroscopy, and pulse field gradient NMR spectroscopy. We connect macroscopic diffusion to segmental chain dynamics, revealing a decoupling of local chain dynamics from interchain interactions. Collective diffusion is described within a colloidal picture, including electrostatic and hydrodynamic interactions. Chain dynamics is characterized by the classical Zimm model of a neutral chain retarded by internal friction. We observe that hydrodynamic interactions are not fully screened between chains and that the internal friction within the chain increases with an increase in ion condensation on the chain.

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.

Neutron Scattering Techniques and Complementary Methods for Structural and Functional Studies of Biological Macromolecules and Large Macromolecular Complexes

Abstract

The review describes the application of small-angle scattering (SAS) of neutrons and complementary methods to study the structures of biomacromolecules. Here we cover SAS techniques, such as the contrast variation, the neutron spin-echo, and the solution of direct and inverse problems of three-dimensional reconstruction of the structures of macromolecules from SAS spectra by means of molecular modeling. A special section is devoted to specific objects of research, such as supramolecular complexes, influenza virus nucleoprotein, and chromatin.

Structure and Dynamics of Ribonuclease A during Thermal Unfolding: The Failure of the Zimm Model

Disordered regions as found in intrinsically disordered proteins (IDP) or during protein folding define response time to stimuli and protein folding times. Neutron spin-echo spectroscopy is a powerful tool to directly access the collective motions of the unfolded chain to enlighten the physical origin of basic conformational relaxation. During the thermal unfolding of native ribonuclease A, we examine the structure and dynamics of the disordered state within a two-state transition model using polymer models, including internal friction, to describe the chain dynamics. The presence of four disulfide bonds alters the disordered configuration to a more compact configuration compared to a Gaussian chain that is defined by the additional links, as demonstrated by coarse-grained simulation. The dynamics of the disordered chain is described by Zimm dynamics with internal friction (ZIF) between neighboring amino acids. Relaxation times are dominated by mode-independent internal friction. Internal friction relaxation times show an Arrhenius-like behavior with an activation energy of 33 kJ/mol. The Zimm dynamics is dominated by internal friction and suggest that the characteristic motions correspond to overdamped elastic modes similar to the motions observed for folded proteins but within a pool of disordered configurations spanning the configurational space. For IDP, internal friction dominates while solvent friction and hydrodynamic interactions are smaller corrections.

Mechanochemical enzymes and protein machines as hydrodynamic force dipoles: the active dimer model

Mechanochemically active enzymes change their shapes within every turnover cycle. Therefore, they induce circulating flows in the solvent around them and behave as oscillating hydrodynamic force dipoles. Because of non-equilibrium fluctuating flows collectively generated by the enzymes, mixing in the solution and diffusion of passive particles within it are expected to get enhanced. Here, we investigate the intensity and statistical properties of such force dipoles in the minimal active dimer model of a mechanochemical enzyme. In the framework of this model, novel estimates for hydrodynamic collective effects in solution and in lipid bilayers under rapid rotational diffusion are derived, and available experimental and computational data is examined.

How fast are the motions of tertiary-structure elements in proteins?

Protein motions occur on multiple time and distance scales. Large-scale motions of protein tertiary-structure elements, i.e., domains, are particularly intriguing as they are essential for the catalytic activity of many enzymes and for the functional cycles of protein machines and motors. Theoretical estimates suggest that domain motions should be very fast, occurring on the nanosecond or microsecond time scales. Indeed, free-energy barriers for domain motions are likely to involve salt bridges, which can break in microseconds. Experimental methods that can directly probe domain motions on fast time scales have appeared only in recent years. This Perspective discusses briefly some of these techniques, including nuclear magnetic resonance and single-molecule fluorescence spectroscopies. We introduce a few recent studies that demonstrate ultrafast domain motions and discuss their potential roles. Particularly surprising is the observation of tertiary-structure element dynamics that are much faster than the functional cycles in some protein machines. These swift motions can be rationalized on a case-by-case basis. For example, fast domain closure in multi-substrate enzymes may be utilized to optimize relative substrate orientation. Whether a large mismatch in time scales of conformational dynamics vs functional cycles is a general design principle in proteins remains to be determined.

Effect of red blood cell shape changes on haemoglobin interactions and dynamics: a neutron scattering study

By using a combination of experimental neutron scattering techniques, it is possible to obtain a statistical perspective on red blood cell (RBC) shape in suspensions, and the inter-relationship with protein interactions and dynamics inside the confinement of the cell membrane. In this study, we examined the ultrastructure of RBC and protein-protein interactions of haemoglobin (Hb) in them using ultra-small-angle neutron scattering and small-angle neutron scattering (SANS). In addition, we used the neutron backscattering method to access Hb motion on the ns time scale and Å length scale. Quasi-elastic neutron scattering (QENS) experiments were performed to measure diffusive motion of Hb in RBCs and in an RBC lysate. By using QENS, we probed both internal Hb dynamics and global protein diffusion, on the accessible time scale and length scale by QENS. Shape changes of RBCs and variation of intracellular Hb concentration were induced by addition of the Na⁺-selective ionophore monensin and the K⁺-selective one, valinomycin. The experimental SANS and QENS results are discussed within the framework of crowded protein solutions, where free motion of Hb is obstructed by mutual interactions.

Reduced Internal Friction by Osmolyte Interaction in Intrinsically Disordered Myelin Basic Protein

Urea is a strong denaturing osmolyte that disrupts noncovalent bonds in proteins. Here, we present a small-angle neutron scattering (SANS) and neutron spin-echo spectroscopy (NSE) study on the structure and dynamics of the intrinsically disordered myelin basic protein (MBP) denatured by urea. SANS results show that urea-denatured MBP is more compact than ideal polymers, while its secondary structure content is entirely lost. NSE experiments reveal concomitantly an increase of the relaxation time and of the amplitude of internal motions in urea-denatured MBP as compared to native MBP. If interpreted in terms of the Zimm model including internal friction (ZIF), the internal friction parameter decreased by a factor of 6.5. Urea seems to not only smooth local energy barriers, reducing internal friction on a local scale, but also significantly reduces the overall depth of the global energy landscape. This leads to a nearly complete loss of restoring forces beyond entropic forces and in turn allows for larger motional amplitudes. Obviously, the noncovalent H-bonds are largely eliminated, driving the unfolded protein to be more similar to a synthetic polymer.

Neutron scattering for STRUCTURAL BIOLOGY: Modern neutron sources illuminate the complex functions of living systems

Modern neutron sources illuminate the complex functions of living systems.

Localised contacts lead to nanosecond hinge motions in dimeric bovine serum albumin

Domain motions in proteins are crucial for biological function. In the present manuscript, we present a neutron spin-echo spectroscopy (NSE) study of native bovine serum albumin (BSA) in solution. NSE allows to probe both global and internal dynamics of the BSA monomer and dimer equilibrium that is formed in solution. Using a model independent approach, we were able to identify an internal dynamic process in BSA that is visible in addition to global rigid-body diffusion of the BSA monomer and dimer mixture. The observed internal protein motion is characterised by a relaxation time of 43 ns. The overdamped Brownian oscillator was considered as an alternative analytical theory that was able to describe the internal process as first-order approximation. More detailed information on the physical nature of the internal protein motion was extracted from the q-dependent internal diffusion coefficients ΔDeff(q) that were detected by NSE in addition to global rigid-body translational and rotational diffusion. The ΔDeff(q) were interpreted using normal mode analysis based on the available crystal structures of the BSA monomer and dimer as structural test models. Normal mode analysis demonstrates that the observed internal dynamic process can be attributed to bending motion of the BSA dimer. The native BSA monomer does not show any internal dynamics on the time- and length-scales probed by NSE. An intermolecular disulphide bridge or a direct structural contact between the BSA monomers forms a localised link acting as a molecular hinge in the BSA dimer. The effect of that hinge on the observed motion of BSA in the used dimeric structural model is discussed in terms of normal modes in a molecular picture.

Heat shock-induced chaperoning by Hsp70 is enabled in-cell

Recent work has shown that weak protein-protein interactions are susceptible to the cellular milieu. One case in point is the binding of heat shock proteins (Hsps) to substrate proteins in cells under stress. Upregulation of the Hsp70 chaperone machinery at elevated temperature was discovered in the 1960s, and more recent studies have shown that ATPase activity in one Hsp70 domain is essential for control of substrate binding by the other Hsp70 domain. Although there are several denaturant-based assays of Hsp70 activity, reports of ATP-dependent binding of Hsp70 to a globular protein substrate under heat shock are scarce. Here we show that binding of heat-inducible Hsp70 to phosphoglycerate kinase (PGK) is remarkably different in vitro compared to in-cell. We use fluorescent-labeled mHsp70 and ePGK, and begin by showing that mHsp70 passes the standard β-galactosidase assay, and that it does not self-aggregate until 50°C in presence of ATP. Yet during denaturant refolding or during in vitro heat shock, mHsp70 shows only ATP-independent non-specific sticking to ePGK, as evidenced by nearly identical results with an ATPase activity-deficient K71M mutant of Hsp70 as a control. Addition of Hsp40 (co-factor) or Ficoll (crowder) does not reduce non-specific sticking, but cell lysate does. Therefore, Hsp70 does not act as an ATP-dependent chaperone on its substrate PGK in vitro. In contrast, we observe only specific ATP-dependent binding of mHsp70 to ePGK in mammalian cells, when compared to the inactive Hsp70 K71M mutant. We hypothesize that enhanced in-cell activity is not due to an unknown co-factor, but simply to a favorable shift in binding equilibrium caused by the combination of crowding and osmolyte/macromolecular interactions present in the cell. One candidate mechanism for such a favorable shift in binding equilibrium is the proven ability of Hsp70 to bind near-native states of substrate proteins in vitro. We show evidence for early onset of binding in-cell. Our results suggest that Hsp70 binds PGK preemptively, prior to its full unfolding transition, thus stabilizing it against further unfolding. We propose a "preemptive holdase" mechanism for Hsp70-substrate binding. Given our result for PGK, more proteins than one might think based on in vitro assays may be chaperoned by Hsp70 in vivo. The cellular environment thus plays an important role in maintaining proper Hsp70 function.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Small-angle X-ray scattering study of the kinetics of light-dark transition in a LOV protein

Light, oxygen, voltage (LOV) photoreceptors consist of conserved photo-responsive domains in bacteria, archaea, plants and fungi, and detect blue-light via a flavin cofactor. We investigated the blue-light induced conformational transition of the dimeric photoreceptor PpSB1-LOV-R66I from Pseudomonas putida in solution by using small-angle X-ray scattering (SAXS). SAXS experiments of the fully populated light- and dark-states under steady-state conditions revealed significant structural differences between the two states that are in agreement with the known structures determined by crystallography. We followed the transition from the light- to the dark-state by using SAXS measurements in real-time. A two-state model based on the light- and dark-state conformations could describe the measured time-course SAXS data with a relaxation time τREC of ~ 34 to 35 min being larger than the recovery time found with UV/vis spectroscopy. Unlike the flavin chromophore-based UV/vis method that is sensitive to the local chromophore environment in flavoproteins, SAXS-based assay depends on protein conformational changes and provides with an alternative to measure the recovery kinetics.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

Solution scattering approaches to dynamical ordering in biomolecular systems

Clarification of solution structure and its modulation in proteins and protein complexes is crucially important to understand dynamical ordering in macromolecular systems. Small-angle x-ray scattering (SAXS) and small-angle neutron scattering (SANS) are among the most powerful techniques to derive structural information. Recent progress in sample preparation, instruments and software analysis is opening up a new era for small-angle scattering. In this review, recent progress and trends of SAXS and SANS are introduced from the point of view of instrumentation and analysis, touching on general features and standard methods of small-angle scattering. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Thermal activation of 'allosteric-like' large-scale motions in a eukaryotic Lactate Dehydrogenase

Conformational changes occurring during the enzymatic turnover are essential for the regulation of protein functionality. Individuating the protein regions involved in these changes and the associated mechanical modes is still a challenge at both experimental and theoretical levels. We present here a detailed investigation of the thermal activation of the functional modes and conformational changes in a eukaryotic Lactate Dehydrogenase enzyme (LDH). Neutron Spin Echo spectroscopy and Molecular Dynamics simulations were used to uncover the characteristic length- and timescales of the LDH nanoscale motions in the apo state. The modes involving the catalytic loop and the mobile region around the binding site are activated at room temperature, and match the allosteric reorganisation of bacterial LDHs. In a temperature window of about 15 degrees, these modes render the protein flexible enough and capable of reorganising the active site toward reactive configurations. On the other hand an excess of thermal excitation leads to the distortion of the protein matrix with a possible anti-catalytic effect. Thus, the temperature activates eukaryotic LDHs via the same conformational changes observed in the allosteric bacterial LDHs. Our investigation provides an extended molecular picture of eukaryotic LDH's conformational landscape that enriches the static view based on crystallographic studies alone.

Structure and domain dynamics of human lactoferrin in solution and the influence of Fe(III)-ion ligand binding

Background

Human lactoferrin is an iron-binding protein of the innate immune system consisting of two connected lobes, each with a binding site located in a cleft. The clefts in each lobe undergo a hinge movement from open to close when Fe³⁺ is present in the solution and can be bound. The binding mechanism was assumed to relate on thermal domain fluctuations of the cleft domains prior to binding. We used Small Angle Neutron Scattering and Neutron Spin Echo Spectroscopy to determine the lactoferrin structure and domain dynamics in solution.

Results

When Fe³⁺ is present in solution interparticle interactions change from repulsive to attractive in conjunction with emerging metas aggregates, which are not observed without Fe³⁺. The protein form factor shows the expected change due to lobe closing if Fe³⁺ is present. The dominating motions of internal domain dynamics with relaxation times in the 30–50 ns range show strong bending and stretching modes with a steric suppressed torsion, but are almost independent of the cleft conformation. Thermally driven cleft closing motions of relevant amplitude are not observed if the cleft is open.

Conclusion

The Fe³⁺ binding mechanism is not related to thermal equilibrium fluctuations closing the cleft. A likely explanation may be that upon entering the cleft the iron ion first binds weakly which destabilizes and softens the hinge region and enables large fluctuations that then close the cleft resulting in the final formation of the stable iron binding site and, at the same time, stable closed conformation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13628-016-0032-3) contains supplementary material, which is available to authorized users.

Determination of functional collective motions in a protein at atomic resolution using coherent neutron scattering

Protein function often depends on global, collective internal motions. However, the simultaneous quantitative experimental determination of the forms, amplitudes, and time scales of these motions has remained elusive. We demonstrate that a complete description of these large-scale dynamic modes can be obtained using coherent neutron-scattering experiments on perdeuterated samples. With this approach, a microscopic relationship between the structure, dynamics, and function in a protein, cytochrome P450cam, is established. The approach developed here should be of general applicability to protein systems.

Quasielastic neutron scattering in biology: Theory and applications

Neutrons scatter quasielastically from stochastic, diffusive processes, such as overdamped vibrations, localized diffusion and transitions between energy minima. In biological systems, such as proteins and membranes, these relaxation processes are of considerable physical interest. We review here recent methodological advances and applications of quasielastic neutron scattering (QENS) in biology, concentrating on the role of molecular dynamics simulation in generating data with which neutron profiles can be unambiguously interpreted. We examine the use of massively-parallel computers in calculating scattering functions, and the application of Markov state modeling. The decomposition of MD-derived neutron dynamic susceptibilities is described, and the use of this in combination with NMR spectroscopy. We discuss dynamics at very long times, including approximations to the infinite time mean-square displacement and nonequilibrium aspects of single-protein dynamics. Finally, we examine how neutron scattering and MD can be combined to provide information on lipid nanodomains. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Atomistic details of protein dynamics and the role of hydration water

BACKGROUND

The importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited.

SCOOP OF REVIEW

Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes.

MAJOR CONCLUSIONS

Internal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10-100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

GENERAL SIGNIFICANCE

The proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Fast antibody fragment motion: flexible linkers act as entropic spring

A flexible linker region between three fragments allows antibodies to adjust their binding sites to an antigen or receptor. Using Neutron Spin Echo Spectroscopy we observed fragment motion on a timescale of 7 ns with motional amplitudes of about 1 nm relative to each other. The mechanistic complexity of the linker region can be described by a spring model with Brownian motion of the fragments in a harmonic potential. Displacements, timescale, friction and force constant of the underlying dynamics are accessed. The force constant exhibits a similar strength to an entropic spring, with friction of the fragment matching the unbound state. The observed fast motions are fluctuations in pre-existing equilibrium configurations. The Brownian motion of domains in a harmonic potential is the appropriate model to examine functional hinge motions dependent on the structural topology and highlights the role of internal forces and friction to function.

Structural and mechanistic insight into substrate binding from the conformational dynamics in apo and substrate-bound DapE enzyme

Conformational dynamics induced by substrate binding in DapE enzyme.

Structure and dynamics of a compact state of a multidomain protein, the mercuric ion reductase

The functional efficacy of colocalized, linked protein domains is dependent on linker flexibility and system compaction. However, the detailed characterization of these properties in aqueous solution presents an enduring challenge. Here, we employ a novel, to our knowledge, combination of complementary techniques, including small-angle neutron scattering, neutron spin-echo spectroscopy, and all-atom molecular dynamics and coarse-grained simulation, to identify and characterize in detail the structure and dynamics of a compact form of mercuric ion reductase (MerA), an enzyme central to bacterial mercury resistance. MerA possesses metallochaperone-like N-terminal domains (NmerA) tethered to its catalytic core domain by linkers. The NmerA domains are found to interact principally through electrostatic interactions with the core, leashed by the linkers so as to subdiffuse on the surface over an area close to the core C-terminal Hg(II)-binding cysteines. How this compact, dynamical arrangement may facilitate delivery of Hg(II) from NmerA to the core domain is discussed.

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

Hierarchical Conformational Analysis of Native Lysozyme Based on Sub-Millisecond Molecular Dynamics Simulations

Hierarchical organization of free energy landscape (FEL) for native globular proteins has been widely accepted by the biophysics community. However, FEL of native proteins is usually projected onto one or a few dimensions. Here we generated collectively 0.2 milli-second molecular dynamics simulation trajectories in explicit solvent for hen egg white lysozyme (HEWL), and carried out detailed conformational analysis based on backbone torsional degrees of freedom (DOF). Our results demonstrated that at micro-second and coarser temporal resolutions, FEL of HEWL exhibits hub-like topology with crystal structures occupying the dominant structural ensemble that serves as the hub of conformational transitions. However, at 100ns and finer temporal resolutions, conformational substates of HEWL exhibit network-like topology, crystal structures are associated with kinetic traps that are important but not dominant ensembles. Backbone torsional state transitions on time scales ranging from nanoseconds to beyond microseconds were found to be associated with various types of molecular interactions. Even at nanoseconds temporal resolution, the number of conformational substates that are of statistical significance is quite limited. These observations suggest that detailed analysis of conformational substates at multiple temporal resolutions is both important and feasible. Transition state ensembles among various conformational substates at microsecond temporal resolution were observed to be considerably disordered. Life times of these transition state ensembles are found to be nearly independent of the time scales of the participating torsional DOFs.

Slow internal protein dynamics in solution

Large-scale domain dynamics in proteins are found when flexible linkers or hinges connect domains. The related conformational changes are often related to the function of the protein,for example by arranging the active center after substrate binding or allowing transport and release of products. The adaptation of a specific active structure is referred to as ‘induced fit’ and is challenged by models such as ‘conformational sampling’. Newer models about protein unction include some flexibility within the protein structure or even internal dynamics of the protein. As larger domains contribute to the configurational changes, the timescale of the involved motions is slowed down. The role of slow domain dynamics is being increasingly recognized as essential to understanding the function of proteins. Neutron spin echospectroscopy (NSE) is a technique that is able to access the related timescales from 0.1 up to several hundred nanoseconds and simultaneously covers the length scale relevant for protein domain movements of several nanometers distance between domains. Here we focus on these large-scale domain fluctuations and show how the structure and dynamics of proteins can be assessed by small-angle neutron scattering and NSE.

Observation of small cluster formation in concentrated monoclonal antibody solutions and its implications to solution viscosity

Monoclonal antibodies (mAbs) are a major class of biopharmaceuticals. It is hypothesized that some concentrated mAb solutions exhibit formation of a solution phase consisting of reversibly self-associated aggregates (or reversible clusters), which is speculated to be responsible for their distinct solution properties. Here, we report direct observation of reversible clusters in concentrated solutions of mAbs using neutron spin echo. Specifically, a stable mAb solution is studied across a transition from dispersed monomers in dilute solution to clustered states at more concentrated conditions, where clusters of a preferred size are observed. Once mAb clusters have formed, their size, in contrast to that observed in typical globular protein solutions, is observed to remain nearly constant over a wide range of concentrations. Our results not only conclusively establish a clear relationship between the undesirable high viscosity of some mAb solutions and the formation of reversible clusters with extended open structures, but also directly observe self-assembled mAb protein clusters of preferred small finite size similar to that in micelle formation that dominate the properties of concentrated mAb solutions.

Conformational state distributions and catalytically relevant dynamics of a hinge-bending enzyme studied by single-molecule FRET and a coarse-grained simulation

Over the last few decades, a view has emerged showing that multidomain enzymes are biological machines evolved to harness stochastic kicks of solvent particles into highly directional functional motions. These intrinsic motions are structurally encoded, and Nature makes use of them to catalyze chemical reactions by means of ligand-induced conformational changes and states redistribution. Such mechanisms align reactive groups for efficient chemistry and stabilize conformers most proficient for catalysis. By combining single-molecule Förster resonance energy transfer measurements with normal mode analysis and coarse-grained mesoscopic simulations, we obtained results for a hinge-bending enzyme, namely phosphoglycerate kinase (PGK), which support and extend these ideas. From single-molecule Förster resonance energy transfer, we obtained insight into the distribution of conformational states and the dynamical properties of the domains. The simulations allowed for the characterization of interdomain motions of a compact state of PGK. The data show that PGK is intrinsically a highly dynamic system sampling a wealth of conformations on timescales ranging from nanoseconds to milliseconds and above. Functional motions encoded in the fold are performed by the PGK domains already in its ligand-free form, and substrate binding is not required to enable them. Compared to other multidomain proteins, these motions are rather fast and presumably not rate-limiting in the enzymatic reaction. Ligand binding slightly readjusts the orientation of the domains and feasibly locks the protein motions along a preferential direction. In addition, the functionally relevant compact state is stabilized by the substrates, and acts as a prestate to reach active conformations by means of Brownian motions.

Coherent neutron scattering and collective dynamics in the protein, GFP

Collective dynamics are considered to be one of the major properties of soft materials, including biological macromolecules. We present coherent neutron scattering studies of the low-frequency vibrations, the so-called boson peak, in fully deuterated green fluorescent protein (GFP). Our analysis revealed unexpectedly low coherence of the atomic motions in GFP. This result implies a low amount of in-phase collective motion of the secondary structural units contributing to the boson peak vibrations and fast conformational fluctuations on the picosecond timescale. These observations are in contrast to earlier studies of polymers and glass-forming systems, and suggest that random or out-of-phase motions of the β-strands contribute greater than two-thirds of the intensity to the low-frequency vibrational spectra of GFP.

de Gennes narrowing describes the relative motion of protein domains

The relative motion of structural domains is essential for the biological function of many proteins. Here, by analyzing neutron scattering data and performing molecular dynamics simulations, we find that interdomain motion in several proteins obeys the principle of de Gennes narrowing, in which the wave vector dependence of the interdomain diffusion coefficient is inversely proportional to the interdomain structure factor. Thus, the rate of interdomain motion is inversely proportional to the probability distribution of the spatial configurations of domains.

Protein dynamics are influenced by the order of ligand binding to an antibiotic resistance enzyme

The aminoglycoside N3 acetyltransferase-IIIb (AAC) is responsible for conferring bacterial resistance to a variety of aminoglycoside antibiotics. Nuclear magnetic resonance spectroscopy and dynamic light scattering analyses revealed a surprising result; the dynamics of the ternary complex between AAC and its two ligands, an antibiotic and coenzyme A, are dependent upon the order in which the ligands are bound. Additionally, two structurally similar aminoglycosides, neomycin and paromomycin, induce strikingly different dynamic properties when they are in their ternary complexes. To the best of our knowledge, this is the first example of a system in which two identically productive pathways of forming a simple ternary complex yield significant differences in dynamic properties. These observations emphasize the importance of the sequence of events in achieving optimal protein-ligand interactions and demonstrate that even a minor difference in molecular structure can have a profound effect on biochemical processes.

Modeling of solvent flow effects in enzyme catalysis under physiological conditions

A stochastic model for the dynamics of enzymatic catalysis in explicit, effective solvents under physiological conditions is presented. Analytically-computed first passage time densities of a diffusing particle in a spherical shell with absorbing boundaries are combined with densities obtained from explicit simulation to obtain the overall probability density for the total reaction cycle time of the enzymatic system. The method is used to investigate the catalytic transfer of a phosphoryl group in a phosphoglycerate kinase-ADP-bis phosphoglycerate system, one of the steps of glycolysis. The direct simulation of the enzyme-substrate binding and reaction is carried out using an elastic network model for the protein, and the solvent motions are described by multiparticle collision dynamics which incorporates hydrodynamic flow effects. Systems where solvent-enzyme coupling occurs through explicit intermolecular interactions, as well as systems where this coupling is taken into account by including the protein and substrate in the multiparticle collision step, are investigated and compared with simulations where hydrodynamic coupling is absent. It is demonstrated that the flow of solvent particles around the enzyme facilitates the large-scale hinge motion of the enzyme with bound substrates, and has a significant impact on the shape of the probability densities and average time scales of substrate binding for substrates near the enzyme, the closure of the enzyme after binding, and the overall time of completion of the cycle.

Functional domain motions in proteins on the ~1-100 ns timescale: comparison of neutron spin-echo spectroscopy of phosphoglycerate kinase with molecular-dynamics simulation

Protein function often requires large-scale domain motion. An exciting new development in the experimental characterization of domain motions in proteins is the application of neutron spin-echo spectroscopy (NSE). NSE directly probes coherent (i.e., pair correlated) scattering on the ~1-100 ns timescale. Here, we report on all-atom molecular-dynamics (MD) simulation of a protein, phosphoglycerate kinase, from which we calculate small-angle neutron scattering (SANS) and NSE scattering properties. The simulation-derived and experimental-solution SANS results are in excellent agreement. The contributions of translational and rotational whole-molecule diffusion to the simulation-derived NSE and potential problems in their estimation are examined. Principal component analysis identifies types of domain motion that dominate the internal motion's contribution to the NSE signal, with the largest being classic hinge bending. The associated free-energy profiles are quasiharmonic and the frictional properties correspond to highly overdamped motion. The amplitudes of the motions derived by MD are smaller than those derived from the experimental analysis, and possible reasons for this difference are discussed. The MD results confirm that a significant component of the NSE arises from internal dynamics. They also demonstrate that the combination of NSE with MD is potentially useful for determining the forms, potentials of mean force, and time dependence of functional domain motions in proteins.

Light-controlled protein dynamics observed with neutron spin echo measurements

A photoresponsive surfactant has been used as a means to control protein structure and dynamics with light illumination. This cationic azobenzene surfactant, azoTAB, which undergoes a reversible photoisomerization upon exposure to the appropriate wavelength of light, adopts a relatively hydrophobic, trans structure under visible light illumination and a relatively hydrophilic cis structure under UV light illumination. Small-angle neutron scattering (SANS) and neutron spin echo (NSE) spectroscopy were used to measure the tertiary structure and internal dynamics of lysozyme in the presence of the photosurfactant, respectively. The SANS-based in vitro structures indicate that under visible light the photosurfactant induces partial unfolding that principally occurs away from the active site near the hinge region connecting the α and β domains. Upon UV exposure, however, the protein refolds to a nativelike structure. At the same time, enhanced internal dynamics of lysozyme were detected with the surfactant in the trans form through NSE measurements of the Q-dependent effective diffusion coefficient (D(eff)) of the protein. In contrast, the D(eff) values of lysozyme in the presence of cis azoTAB largely agree with the rigid-body calculation as well as those measured for pure lysozyme, suggesting that the native protein is dormant on the nanosecond time and nanometer length scales. Lysozyme internal motions were modeled by assuming a protein of two (α and β domains) or three (α and β domains and the hinge region) domains connects by either soft linkers or rigid, freely rotating bonds. Protein dynamics were also tracked with Fourier transform infrared spectroscopy through hydrogen-deuterium exchange kinetics, which further demonstrated enhanced protein flexibility induced by the trans form of the surfactant relative to the native protein. Ensemble-averaged intramolecular fluorescent resonance energy transfer measurements similarly demonstrated the enhanced dynamics of lysozyme with the trans form of the photosurfactant. Previous results have shown a significant increase in protein activity in the presence of azoTAB in the trans conformation. Combined, these results provide insight into a unique light-based method of controlling protein structure, dynamics, and function and strongly support the relevance of large domain motions for the activity of proteins.