Coupled protein domain motion in Taq polymerase revealed by neutron spin-echo spectroscopy

Nanoscale dynamics of the cadherin-catenin complex bound to vinculin revealed by neutron spin echo spectroscopy

We report a neutron spin echo (NSE) study of the nanoscale dynamics of the cell-cell adhesion cadherin-catenin complex bound to vinculin. Our measurements and theoretical physics analyses of the NSE data reveal that the dynamics of full-length α-catenin, β-catenin, and vinculin residing in the cadherin-catenin-vinculin complex become activated, involving nanoscale motions in this complex. The cadherin-catenin complex is the central component of the cell-cell adherens junction (AJ) and is fundamental to embryogenesis, tissue wound healing, neuronal plasticity, cancer metastasis, and cardiovascular health and disease. A highly dynamic cadherin-catenin-vinculin complex provides the molecular dynamics basis for the flexibility and elasticity that are necessary for the AJs to function as force transducers. Our theoretical physics analysis provides a way to elucidate these driving nanoscale motions within the complex without requiring large-scale numerical simulations, providing insights not accessible by other techniques. We propose a three-way "motorman" entropic spring model for the dynamic cadherin-catenin-vinculin complex, which allows the complex to function as a flexible and elastic force transducer.

An efficient approach for overproduction of DNA polymerase from Pyrococcus furiosus using an optimized autoinduction system in Escherichia coli

High fidelity DNA polymerase from Pyrococcus furiosus (Pfupol) is an attractive alternative to the highly popular DNA polymerase from Thermus aquaticus. Because this enzyme is in great demand for biotechnological applications, optimizing Pfupol production is essential to supplying the industry's expanding demand. T7-induced promoter expression in Escherichia coli expression systems is used to express recombinant Pfupol; however, this method is not cost-effective. Here, we have effectively developed an optimized process for the autoinduction approach of Pfupol expression in a defined medium. To better examine Pfupol's activities, its purified fraction was used. A 71 mg/L of pure Pfupol was effectively produced, resulting in a 2.6-fold increase in protein yield when glucose, glycerol, and lactose were added in a defined medium at concentrations of 0.05%, 1%, and 0.6%, respectively, and the condition for production in a 5 L bioreactor was as follow: 200 rpm, 3 vvm, and 10% inoculant. Furthermore, the protein exhibited 1445 U/mg of specific activity when synthesized in its active state. This work presents a high level of Pfupol production, which makes it an economically viable and practically useful approach.

Codon optimization of a gene encoding DNA polymerase from Pyrococcus furiosus and its expression in Escherichia coli

BACKGROUND

DNA polymerase is an essential component in PCR assay for DNA synthesis. Improving DNA polymerase with characteristics indispensable for a powerful assay is crucial because it can be used in wide-range applications. Derived from Pyrococcus furiosus, Pfu DNA polymerase (Pfu pol) is one of the excellent polymerases due to its high fidelity. Therefore, we aimed to develop Pfu pol from a synthetic gene with codon optimization to increase its protein yield in Escherichia coli.

RESULTS

Recombinant Pfu pol was successfully expressed and purified with a two-step purification process using nickel affinity chromatography, followed by anion exchange chromatography. Subsequently, the purified Pfu pol was confirmed by Western blot analysis, resulting in a molecular weight of approximately 90 kDa. In the final purification process, we successfully obtained a large amount of purified enzyme (26.8 mg/L). Furthermore, the purified Pfu pol showed its functionality and efficiency when tested for DNA amplification using the standard PCR.

CONCLUSIONS

Overall, a high-level expression of recombinant Pfu pol was achieved by employing our approach in the present study. In the future, our findings will be useful for studies on synthesizing recombinant DNA polymerase in E. coli expression system.

How neutron scattering techniques benefit investigating structures and dynamics of monoclonal antibody

Over the past several decades, great progresses have been made for the pharmaceutical industry of monoclonal antibody (mAb). More and more mAb products were approved for human therapeutics. This review describes the state of art of utilizing neutron scattering to investigate mAbs, in the aspects of structures, dynamics, physicochemical stability, functionality, etc. Firstly, brief histories of mAbs and neutron scattering, as well as some basic knowledges and principles of neutron scattering were introduced. Then specific examples were demonstrated. For the structure and structural evolution investigation of in dilute and concentrated mAbs solution, in situ small angle neutron scattering (SANS) was frequently utilized. Neutron reflectometry (NR) is powerful to probe the absorption behaviors of mAbs on various surfaces and interfaces. While for dynamic investigation, quasi-elastic scattering techniques such as neutron spin echo (NSE) demonstrate the capabilities. With this review, how to utilize and take advantages of neutron scattering on investigating structures and dynamics of mAbs were demonstrated and discussed.

Nanosecond structural dynamics of intrinsically disordered β-casein micelles by neutron spectroscopy

β-casein undergoes a reversible endothermic self-association, forming protein micelles of limited size. In its functional state, a single β-casein monomer is unfolded, which creates a high structural flexibility, which is supposed to play a major role in preventing the precipitation of calcium phosphate particles. We characterize the structural flexibility in terms of nanosecond molecular motions, depending on the temperature by quasielastic neutron scattering. Our major questions are: Does the self-association reduce the chain flexibility? How does the dynamic spectrum of disordered caseins differ from a compactly globular protein? How does the dynamic spectrum of β-casein in solution differ from that of a protein in hydrated powder states? We report on two relaxation processes on a nanosecond and a sub-nanosecond timescale for β-casein in solution. Both processes are analyzed by Brownian oscillator model, by which the spring constant can be defined in the isotropic parabolic potential. The slower process, which is analyzed by neutron spin echo, seems a characteristic feature of the unfolded structure. It requires bulk solvent and is not seen in hydrated protein powders. The faster process, which is analyzed by neutron backscattering, has a smaller amplitude and requires hydration water, which is also observed with folded proteins in the hydrated state. The self-association had no significant influence on internal relaxation, and thus, a β-casein protein monomer flexibility is preserved in the micelle. We derive spring constants of the faster and slower motions of β-caseins in solution and compared them with those of some proteins in various states (folded or hydrated powder).

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.

Activated nanoscale actin-binding domain motion in the catenin-cadherin complex revealed by neutron spin echo spectroscopy

As the core component of the adherens junction in cell-cell adhesion, the cadherin-catenin complex transduces mechanical tension between neighboring cells. Structural studies have shown that the cadherin-catenin complex exists as an ensemble of flexible conformations, with the actin-binding domain (ABD) of α-catenin adopting a variety of configurations. Here, we have determined the nanoscale protein domain dynamics of the cadherin-catenin complex using neutron spin echo spectroscopy (NSE), selective deuteration, and theoretical physics analyses. NSE reveals that, in the cadherin-catenin complex, the motion of the entire ABD becomes activated on nanosecond to submicrosecond timescales. By contrast, in the α-catenin homodimer, only the smaller disordered C-terminal tail of ABD is moving. Molecular dynamics (MD) simulations also show increased mobility of ABD in the cadherin-catenin complex, compared to the α-catenin homodimer. Biased MD simulations further reveal that the applied external forces promote the transition of ABD in the cadherin-catenin complex from an ensemble of diverse conformational states to specific states that resemble the actin-bound structure. The activated motion and an ensemble of flexible configurations of the mechanosensory ABD suggest the formation of an entropic trap in the cadherin-catenin complex, serving as negative allosteric regulation that impedes the complex from binding to actin under zero force. Mechanical tension facilitates the reduction in dynamics and narrows the conformational ensemble of ABD to specific configurations that are well suited to bind F-actin. Our results provide a protein dynamics and entropic explanation for the observed force-sensitive binding behavior of a mechanosensitive protein complex.

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

A novel approach for high-level expression and purification of GST-fused highly thermostable Taq DNA polymerase in Escherichia coli

Polymerases are enzymes that synthesize long chains or polymers of nucleic acids including DNA or RNA from nucleotides. They assemble nucleic acids by copying a DNA or RNA template strand using base-pairing interactions. One of the polymerase enzymes, Taq DNA polymerase, originally isolated from Thermus aquaticus (Taq) is a widely used enzyme in molecular biology so far. The thermostable properties of this enzyme have contributed majorly to the specificity, automation, and efficacy of the polymerase chain reaction (PCR), making it a powerful tool for today's molecular biology researches across the globe. The purification of Taq DNA polymerase from the native host results in low yield, more labor and time consumption. Therefore, many studies have been previously conducted to obtain this enzyme using alternative hosts. So far, all the existing methodologies are more laborious, time-consuming and require heavy expense. We used a novel approach to purify the enzyme with relatively high efficiency, yield and minimum time consumption using Escherichia coli (E. coli) as an alternative host. We cloned a 2500 base pair Taq DNA polymerase gene into pGEX-4T-1 vector, containing a GST-tag, downstream of tac promoter and overexpressed it using isopropyl β-d-1-thiogalactopyranoside (IPTG) as an inducer. The enzyme was efficiently purified using novel chromatography approaches and was used in routine PCR assays in our laboratory. Our findings suggest a novel approach to facilitate the availability of polymerases for molecular and diagnostic studies. In the future, it may be used for the purification of other recombinant peptides or proteins used in structural biology and proteomics-based researches.

Transition between protein-like and polymer-like dynamic behavior: Internal friction in unfolded apomyoglobin depends on denaturing conditions

Equilibrium dynamics of different folding intermediates and denatured states is strongly connected to the exploration of the conformational space on the nanosecond time scale and might have implications in understanding protein folding. For the first time, the same protein system apomyoglobin has been investigated using neutron spin-echo spectroscopy in different states: native-like, partially folded (molten globule) and completely unfolded, following two different unfolding paths: using acid or guanidinium chloride (GdmCl). While the internal dynamics of the native-like state can be understood using normal mode analysis based on high resolution structural information of myoglobin, for the unfolded and even for the molten globule states, models from polymer science are employed. The Zimm model accurately describes the slowly-relaxing, expanded GdmCl-denaturated state, ignoring the individuality of the different aminoacid side chain. The dynamics of the acid unfolded and molten globule state are similar in the framework of the Zimm model with internal friction, where the chains still interact and hinder each other: the first Zimm relaxation time is as large as the internal friction time. Transient formation of secondary structure elements in the acid unfolded and presence of α-helices in the molten globule state lead to internal friction to a similar extent.

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.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

α-Catenin Structure and Nanoscale Dynamics in Solution and in Complex with F-Actin

As a core component of the adherens junction, α-catenin stabilizes the cadherin/catenin complexes to the actin cytoskeleton for the mechanical coupling of cell-cell adhesion. α-catenin also modulates actin dynamics, cell polarity, and cell-migration functions that are independent of the adherens junction. We have determined the solution structures of the α-catenin monomer and dimer using in-line size-exclusion chromatography small-angle X-ray scattering, as well as the structure of α-catenin dimer in complex to F-actin filament using selective deuteration and contrast-matching small angle neutron scattering. We further present the first observation, to our knowledge, of the nanoscale dynamics of α-catenin by neutron spin-echo spectroscopy, which explicitly reveals the mobile regions of α-catenin that are crucial for binding to F-actin. In solution, the α-catenin monomer is more expanded than either protomer shown in the crystal structure dimer, with the vinculin-binding M fragment and the actin-binding domain being able to adopt different configurations. The α-catenin dimer in solution is also significantly more expanded than the dimer crystal structure, with fewer interdomain and intersubunit contacts than the crystal structure. When in complex to F-actin, the α-catenin dimer has an even more open and extended conformation than in solution, with the actin-binding domain further separated from the main body of the dimer. The α-catenin-assembled F-actin bundle develops into an ordered filament packing arrangement at increasing α-catenin/F-actin molar ratios. Together, the structural and dynamic studies reveal that α-catenin possesses dynamic molecular conformations that prime this protein to function as a mechanosensor protein.

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.

Anesthetic Alterations of Collective Terahertz Oscillations in Tubulin Correlate with Clinical Potency: Implications for Anesthetic Action and Post-Operative Cognitive Dysfunction

Anesthesia blocks consciousness and memory while sparing non-conscious brain activities. While the exact mechanisms of anesthetic action are unknown, the Meyer-Overton correlation provides a link between anesthetic potency and solubility in a lipid-like, non-polar medium. Anesthetic action is also related to an anesthetic's hydrophobicity, permanent dipole, and polarizability, and is accepted to occur in lipid-like, non-polar regions within brain proteins. Generally the protein target for anesthetics is assumed to be neuronal membrane receptors and ion channels, however new evidence points to critical effects on intra-neuronal microtubules, a target of interest due to their potential role in post-operative cognitive dysfunction (POCD). Here we use binding site predictions on tubulin, the protein subunit of microtubules, with molecular docking simulations, quantum chemistry calculations, and theoretical modeling of collective dipole interactions in tubulin to investigate the effect of a group of gases including anesthetics, non-anesthetics, and anesthetic/convulsants on tubulin dynamics. We found that these gases alter collective terahertz dipole oscillations in a manner that is correlated with their anesthetic potency. Understanding anesthetic action may help reveal brain mechanisms underlying consciousness, and minimize POCD in the choice and development of anesthetics used during surgeries for patients suffering from neurodegenerative conditions with compromised cytoskeletal microtubules.

Controllable Activation of Nanoscale Dynamics in a Disordered Protein Alters Binding Kinetics

The phosphorylation of specific residues in a flexible disordered activation loop yields precise control of signal transduction. One paradigm is the phosphorylation of S339/S340 in the intrinsically disordered tail of the multi-domain scaffolding protein NHERF1, which affects the intracellular localization and trafficking of NHERF1 assembled signaling complexes. Using neutron spin echo spectroscopy (NSE), we show salt-concentration-dependent excitation of nanoscale motion at the tip of the C-terminal tail in the phosphomimic S339D/S340D mutant. The "tip of the whip" that is unleashed is near the S339/S340 phosphorylation site and flanks the hydrophobic Ezrin-binding motif. The kinetic association rate constant of the binding of the S339D/S340D mutant to the FERM domain of Ezrin is sensitive to buffer salt concentration, correlating with the excited nanoscale dynamics. The results suggest that electrostatics modulates the activation of nanoscale dynamics of an intrinsically disordered protein, controlling the binding kinetics of signaling partners. NSE can pinpoint the nanoscale dynamics changes in a highly specific manner.

Intermediate scattering function for macromolecules in solutions probed by neutron spin echo

The neutron-spin-echo method (NSE) is a powerful technique for studying internal dynamics of macromolecules in solutions because it can simultaneously probe length and time scales comparable to intramolecular density fluctuations of macromolecules. Recently, there has been increased, strong interest in studying protein internal motions using NSE. The coherent intermediate scattering function (ISF) measured by NSE depends on internal, rotational, and translational motions of macromolecules in solutions. It is thus critical, but highly nontrivial, to separate the internal motion from other motions in order to properly understand protein internal dynamics. Even though many experiments are performed at relatively high concentrations, current theories of calculating the ISF of concentrated protein solutions are either inaccurate or flawed by incorrect assumptions for realistic protein systems with anisotropic shapes. Here, a theoretical framework is developed to establish the quantitative relationship of different motions included in the ISF. This theory based on the dynamic decoupling approximation is applicable to a wide range of protein concentrations, including dilute cases. It is also, in general, useful for studying many other types of macromolecule systems studied by NSE.

Modulation of hemoglobin dynamics by an allosteric effector

Hemoglobin (Hb) is an extensively studied paradigm of proteins that alter their function in response to allosteric effectors. Models of its action have been used as prototypes for structure‐function relationships in many proteins, and models for the molecular basis of its function have been deeply studied and extensively argued. Recent reports suggest that dynamics may play an important role in its function. Relatively little is known about the slow, correlated motions of hemoglobin subunits in various structural states because experimental and computational strategies for their characterization are challenging. Allosteric effectors such as inositol hexaphosphate (IHP) bind to both deoxy‐Hb and HbCO, albeit at different sites, leading to a lowered oxygen affinity. The manner in which these effectors impact oxygen binding is unclear and may involve changes in structure, dynamics or both. Here we use neutron spin echo measurements accompanied by wide‐angle X‐ray scattering to show that binding of IHP to HbCO results in an increase in the rate of coordinated motions of Hb subunits relative to one another with little if any change in large scale structure. This increase of large‐scale dynamics seems to be coupled with a decrease in the average magnitude of higher frequency modes of individual residues. These observations indicate that enhanced dynamic motions contribute to the functional changes induced by IHP and suggest that they may be responsible for the lowered oxygen affinity triggered by these effectors.

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.

Visualizing the nanoscale: protein internal dynamics and neutron spin echo spectroscopy

The most complex molecular machines are proteins found within cells. Protein dynamics, in particular dynamics on nanoscales, presents us with a novel paradigm for cell signaling: the idea that proteins and protein complexes can communicate directly within themselves to effect long-range information transfer, via coupled domains and correlated residue clusters. This idea has been little explored, in large part because of a paucity of experimental techniques that can address the necessary questions. Here we review recent progress in developing a promising new approach, neutron spin echo spectroscopy.

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.

Extended Law of Corresponding States Applied to Solvent Isotope Effect on a Globular Protein

Investigating proteins with techniques such as NMR or neutron scattering frequently requires the partial or complete substitution of D2O for H2O as a solvent, often tacitly assuming that such a solvent substitution does not significantly alter the properties of the protein. Here, we report a systematic investigation of the solvent isotope effect on the phase diagram of the lens protein γB-crystallin in aqueous solution as a model system exhibiting liquid-liquid phase separation. We demonstrate that the observed strong variation of the critical temperature Tc can be described by the extended law of corresponding states for all H2O/D2O ratios, where scaling of the temperature by Tc or the reduced second virial coefficient accurately reproduces the binodal, spinodal, and osmotic compressibility. These findings highlight the impact of H2O/D2O substitution on γB-crystallin properties and warrant further investigations into the universality of this phenomenon and its underlying mechanisms.

Innovative qPCR using interfacial effects to enable low threshold cycle detection and inhibition relief

Molecular diagnostics offers quick access to information but fails to operate at a speed required for clinical decision-making. Our novel methodology, droplet-on-thermocouple silhouette real-time polymerase chain reaction (DOTS qPCR), uses interfacial effects for droplet actuation, inhibition relief, and amplification sensing. DOTS qPCR has sample-to-answer times as short as 3 min 30 s. In infective endocarditis diagnosis, DOTS qPCR demonstrates reproducibility, differentiation of antibiotic susceptibility, subpicogram limit of detection, and thermocycling speeds of up to 28 s/cycle in the presence of tissue contaminants. Langmuir and Gibbs adsorption isotherms are used to describe the decreasing interfacial tension upon amplification. Moreover, a log-linear relationship with low threshold cycles is presented for real-time quantification by imaging the droplet-on-thermocouple silhouette with a smartphone. DOTS qPCR resolves several limitations of commercially available real-time PCR systems, which rely on fluorescence detection, have substantially higher threshold cycles, and require expensive optical components and extensive sample preparation. Due to the advantages of low threshold cycle detection, we anticipate extending this technology to biological research applications such as single cell, single nucleus, and single DNA molecule analyses. Our work is the first demonstrated use of interfacial effects for sensing reaction progress, and it will enable point-of-care molecular diagnosis of infections.

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.

Nanoscale protein domain motion and long-range allostery in signaling proteins- a view from neutron spin echo sprectroscopy

Many cellular proteins are multi-domain proteins. Coupled domain-domain interactions in these multidomain proteins are important for the allosteric relay of signals in the cellular signaling networks. We have initiated the application of neutron spin echo spectroscopy to the study of nanoscale protein domain motions on submicrosecond time scales and on nanometer length scale. Our NSE experiments reveal the activation of protein domain motions over a long distance of over more than 100 Å in a multidomain scaffolding protein NHERF1 upon binding to another protein Ezrin. Such activation of nanoscale protein domains motions is correlated with the allosteric assembly of multi-protein complexes by NHERF1 and Ezrin. Here, we summarize the theoretical framework that we have developed, which uses simple concepts from nonequilibrium statistical mechanics to interpret the NSE data, and employs a mobility tensor to describe nanoscale protein domain motion. Extracting nanoscale protein domain motion from the NSE does not require elaborate molecular dynamics simulations, or complex fits to rotational motion, or elastic network models. The approach is thus more robust than multiparameter techniques that require untestable assumptions. We also demonstrate that an experimental scheme of selective deuteration of a protein subunit in a complex can highlight and amplify specific domain dynamics from the abundant global translational and rotational motions in a protein. We expect NSE to provide a unique tool to determine nanoscale protein dynamics for the understanding of protein functions, such as how signals are propagated in a protein over a long distance to a distal domain.

A structural mechanism for calcium transporter headpiece closure

To characterize the conformational dynamics of sarcoplasmic reticulum (SR) calcium pump (SERCA) we performed molecular dynamics simulations beginning with several different high-resolution structures. We quantified differences in structural disorder and dynamics for an open conformation of SERCA versus closed structures and observed that dynamic motions of SERCA cytoplasmic domains decreased with decreasing domain–domain separation distance. The results are useful for interpretation of recent intramolecular Förster resonance energy transfer (FRET) distance measurements obtained for SERCA fused to fluorescent protein tags. Those previous physical measurements revealed several discrete structural substates and suggested open conformations of SERCA are more dynamic than compact conformations. The present simulations support this hypothesis and provide additional details of SERCA molecular mechanisms. Specifically, all-atoms simulations revealed large-scale translational and rotational motions of the SERCA N-domain relative to the A- and P-domains during the transition from an open to a closed headpiece conformation over the course of a 400 ns trajectory. The open-to-closed structural transition was accompanied by a disorder-to-order transition mediated by an initial interaction of an N-domain loop (Nβ5-β6, residues 426–436) with residues 133–139 of the A-domain. Mutation of three negatively charged N-domain loop residues abolished the disorder-to-order transition and prevented the initial domain–domain interaction and subsequent closure of the cytoplasmic headpiece. Coarse-grained molecular dynamics simulations were in harmony with all-atoms simulations and physical measurements and revealed a close communication between fluorescent protein tags and the domain to which they were fused. The data indicate that previous intramolecular FRET distance measurements report SERCA structure changes with high fidelity and suggest a structural mechanism that facilitates the closure of the SERCA cytoplasmic headpiece.

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.

Dynamical Threshold of Diluteness of Soft Colloids

Soft colloids are hybrids between linear polymers and hard colloids. Their solutions exhibit rich phase phenomenon due to their unique microstructure. In scaling theories, a geometrically defined overlap concentration c* is used to identify the concentration regimes of their solutions characterized with distinct conformational properties. Previous experiments showed that the average size of soft colloids remains invariant below c* and varies characteristically above it. This observation reveals the causality between the conformational evolution and the physical overlap between neighboring particles. Using neutron scattering, we demonstrate that the competition between the interparticle translational diffusion and intraparticle internal dynamics leads to significant conformational evolution below c*. Substantial structural dehydration and slowing-down of internal dynamics are both observed before physical overlap develops. Well below c*, a new threshold of diluteness c_D* emerges as the crossover between the characteristic times associated with these two relaxation processes. Below this dynamically defined c_D*, the two relaxation processes are essentially uncoupled, and therefore, the majority of the soft colloids retain their unperturbed conformational dimensions. Our observation demonstrates the importance of incorporating dynamical degrees of freedom in defining the threshold of diluteness for this important class of soft matter.

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.

Diffusion-based microfluidic PCR for "one-pot" analysis of cells

Genetic analysis starting with cell samples often requires multi-step processing including cell lysis, DNA isolation/purification, and polymerase chain reaction (PCR) based assays. When conducted on a microfluidic platform, the compatibility among various steps often demands a complicated procedure and a complex device structure. Here we present a microfluidic device that permits a "one-pot" strategy for multi-step PCR analysis starting from cells. Taking advantage of the diffusivity difference, we replace the smaller molecules in the reaction chamber by diffusion while retaining DNA molecules inside. This simple scheme effectively removes reagents from the previous step to avoid interference and thus permits multi-step processing in the same reaction chamber. Our approach shows high efficiency for PCR and potential for a wide range of genetic analysis including assays based on single cells.

Molecular conformation of the full-length tumor suppressor NF2/Merlin--a small-angle neutron scattering study

The tumor suppressor protein Merlin inhibits cell proliferation upon establishing cell-cell contacts. Because Merlin has high level of sequence similarity to the Ezrin-Radixin-Moesin family of proteins, the structural model of Ezrin-Radixin-Moesin protein autoinhibition and cycling between closed/resting and open/active conformational states is often employed to explain Merlin function. However, recent biochemical studies suggest alternative molecular models of Merlin function. Here, we have determined the low-resolution molecular structure and binding activity of Merlin and a Merlin(S518D) mutant that mimics the inactivating phosphorylation at S518 using small-angle neutron scattering and binding experiments. Small-angle neutron scattering shows that, in solution, both Merlin and Merlin(S518D) adopt a closed conformation, but binding experiments indicate that a significant fraction of either Merlin or Merlin(S518D) is capable of binding to the target protein NHERF1. Upon binding to the phosphatidylinositol 4,5-bisphosphate lipid, the wild-type Merlin adopts a more open conformation than in solution, but Merlin(S518D) remains in a closed conformation. This study supports a rheostat model of Merlin in NHERF1 binding and contributes to resolving a controversy about the molecular conformation and binding activity of Merlin.

Conformational dynamics of Thermus aquaticus DNA polymerase I during catalysis

Despite the fact that DNA polymerases have been investigated for many years and are commonly used as tools in a number of molecular biology assays, many details of the kinetic mechanism they use to catalyze DNA synthesis remain unclear. Structural and kinetic studies have characterized a rapid, pre-catalytic open-to-close conformational change of the Finger domain during nucleotide binding for many DNA polymerases including Thermus aquaticus DNA polymerase I (Taq Pol), a thermostable enzyme commonly used for DNA amplification in PCR. However, little has been performed to characterize the motions of other structural domains of Taq Pol or any other DNA polymerase during catalysis. Here, we used stopped-flow Förster resonance energy transfer to investigate the conformational dynamics of all five structural domains of the full-length Taq Pol relative to the DNA substrate during nucleotide binding and incorporation. Our study provides evidence for a rapid conformational change step induced by dNTP binding and a subsequent global conformational transition involving all domains of Taq Pol during catalysis. Additionally, our study shows that the rate of the global transition was greatly increased with the truncated form of Taq Pol lacking the N-terminal domain. Finally, we utilized a mutant of Taq Pol containing a de novo disulfide bond to demonstrate that limiting protein conformational flexibility greatly reduced the polymerization activity of Taq Pol.

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.

Nanoporous scaffold for DNA polymerase: pore-size optimisation of mesoporous silica for DNA amplification

The simple and selective immobilisation of a thermostable DNA polymerase on mesoporous silicas was achieved, and DNA amplification activity was retained under the pore-size regulation of the mesoporous silicas.

Nanoscale protein dynamics: a new frontier for neutron spin echo spectroscopy

Recent studies show that neutron spin echo spectroscopy (NSE) can reveal long-range protein domain motions on nanometer lengthscales and on nanosecond to microsecond timescales. This unique capability of NSE provides new opportunities to understand protein dynamics and functions, such as how binding signals are propagated in a protein to distal sites. Here we review our applications of NSE to the study of nanoscale protein domain motions in a set of cell signaling proteins. We summarize the theoretical framework we have developed, which allows one to interpret the NSE data (Biophys. J. 99, 3473 (2010) and Proc. Natl. Acad. Sci. USA 102, 17646 (2005)). Our theoretical framework uses simple concepts from nonequilibrium statistical mechanics, and does not require elaborate molecular dynamics simulations, complex fits to rotational motion, or elastic network models. It is thus more robust than multiparameter techniques that require untestable assumptions. We also demonstrate our experimental scheme involving deuterium labeling of a protein domain or a subunit in a protein complex. We show that our selective deuteration scheme can highlight and resolve specific domain dynamics from the abundant global translational and rotational motions in a protein. Our approach thus clears significant hurdles to the application of NSE for the study of protein dynamics in solution.

Multiple pathways promote dynamical coupling between catalytic domains in Escherichia coli prolyl-tRNA synthetase

Aminoacyl-tRNA synthetases are multidomain enzymes that catalyze covalent attachment of amino acids to their cognate tRNA. Cross-talk between functional domains is a prerequisite for this process. In this study, we investigate the molecular mechanism of site-to-site communication in Escherichia coli prolyl-tRNA synthetase (Ec ProRS). Earlier studies have demonstrated that evolutionarily conserved and/or co-evolved residues that are engaged in correlated motion are critical for the propagation of functional conformational changes from one site to another in modular proteins. Here, molecular simulation and bioinformatics-based analysis were performed to identify dynamically coupled and evolutionarily constrained residues that form contiguous pathways of residue-residue interactions between the aminoacylation and editing domains of Ec ProRS. The results of this study suggest that multiple pathways exist between these two domains to maintain the dynamic coupling essential for enzyme function. Moreover, residues in these interaction networks are generally highly conserved. Site-directed changes of on-pathway residues have a significant impact on enzyme function and dynamics, suggesting that any perturbation along these pathways disrupts the native residue-residue interactions that are required for effective communication between the two functional domains. Free energy analysis revealed that communication between residues within a pathway and cross-talk between pathways are important for coordinating functions of different domains of Ec ProRS for efficient catalysis.

Open conformation of ezrin bound to phosphatidylinositol 4,5-bisphosphate and to F-actin revealed by neutron scattering

Ezrin is a member of the ezrin-radixin-moesin family (ERM) of adapter proteins that are localized at the interface between the cell membrane and the cortical actin cytoskeleton, and they regulate a variety of cellular functions. The structure representing a dormant and closed conformation of an ERM protein has previously been determined by x-ray crystallography. Here, using contrast variation small angle neutron scattering, we reveal the structural changes of the full-length ezrin upon binding to the signaling lipid phosphatidylinositol 4,5-bisphosphate (PIP(2)) and to F-actin. Ezrin binding to F-actin requires the simultaneous binding of ezrin to PIP(2). Once bound to F-actin, the opened ezrin forms more extensive contacts with F-actin than generally depicted, suggesting a possible role of ezrin in regulating the interfacial structure and dynamics between the cell membrane and the underlying actin cytoskeleton. In addition, using gel filtration, we find that the conformational opening of ezrin in response to PIP(2) binding is cooperative, but the cooperativity is disrupted by a phospho-mimic mutation S249D in the 4.1-ezrin/radixin/moesin (FERM) domain of ezrin. Using surface plasmon resonance, we show that the S249D mutation weakens the binding affinity and changes the kinetics of 4.1-ERM to PIP(2) binding. The study provides the first structural view of the activated ezrin bound to PIP(2) and to F-actin.

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.

Role of coupled dynamics in the catalytic activity of prokaryotic-like prolyl-tRNA synthetases

Prolyl-tRNA synthetases (ProRSs) have been shown to activate both cognate and some noncognate amino acids and attach them to specific tRNA(Pro) substrates. For example, alanine, which is smaller than cognate proline, is misactivated by Escherichia coli ProRS. Mischarged Ala-tRNA(Pro) is hydrolyzed by an editing domain (INS) that is distinct from the activation domain. It was previously shown that deletion of the INS greatly reduced cognate proline activation efficiency. In this study, experimental and computational approaches were used to test the hypothesis that deletion of the INS alters the internal protein dynamics leading to reduced catalytic function. Kinetic studies with two ProRS variants, G217A and E218A, revealed decreased amino acid activation efficiency. Molecular dynamics studies showed motional coupling between the INS and protein segments containing the catalytically important proline-binding loop (PBL, residues 199-206). In particular, the complete deletion of INS, as well as mutation of G217 or E218 to alanine, exhibited significant effects on the motion of the PBL. The presence of coupled dynamics between neighboring protein segments was also observed through in silico mutations and essential dynamics analysis. Altogether, this study demonstrates that structural elements at the editing domain-activation domain interface participate in coupled motions that facilitate amino acid binding and catalysis by bacterial ProRSs, which may explain why truncated or defunct editing domains have been maintained in some systems, despite the lack of catalytic activity.