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.     

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

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 domain 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, nor complex fits to rotational motion, nor 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.     

Atomistic ensemble modeling and small-angle neutron scattering of intrinsically disordered protein complexes: applied to minichromosome maintenance protein

The minichromosome maintenance (MCM) proteins are thought to function as the replicative helicases in archaea and eukarya. In this work we determined the solution structure of the N-terminal portion of the MCM complex from the archaeon Methanothermobacter thermautotrophicus (N-mtMCM) in the presence and absence of DNA using small-angle neutron scattering (SANS). N-mtMCM is a multimeric protein complex that consists of 12 monomers, each of which contains three distinct domains and two unstructured regions. Using an all-atom approach incorporating modern force field and Monte Carlo methods to allow the unstructured regions of each monomer to be varied independently, we generated an ensemble of biologically relevant structures for the complex. An examination of the subsets of structures that were most consistent with the SANS data revealed that large movements between the three domains of N-mtMCM can occur in solution. Furthermore, changes in the SANS curves upon DNA binding could be correlated to the motion of a particular N-mtMCM domain. These results provide structural support to the previously reported biochemical observations that large domain motions are required for the activation of the MCM helicase in archaea and eukarya. The methods developed here for N-mtMCM solution structure modeling should be suitable for other large protein complexes with unstructured flexible regions.     

Activation of nanoscale allosteric protein domain motion revealed by neutron spin echo spectroscopy

NHERF1 is a multidomain scaffolding protein that assembles signaling complexes, and regulates the cell surface expression and endocytic recycling of a variety of membrane proteins. The ability of the two PDZ domains in NHERF1 to assemble protein complexes is allosterically modulated by the membrane-cytoskeleton linker protein ezrin, whose binding site is located as far as 110 Ångstroms away from the PDZ domains. Here, using neutron spin echo (NSE) spectroscopy, selective deuterium labeling, and theoretical analyses, we reveal the activation of interdomain motion in NHERF1 on nanometer length-scales and on submicrosecond timescales upon forming a complex with ezrin. We show that a much-simplified coarse-grained model suffices to describe interdomain motion of a multidomain protein or protein complex. We expect that future NSE experiments will benefit by exploiting our approach of selective deuteration to resolve the specific domain motions of interest from a plethora of global translational and rotational motions. Our results demonstrate that the dynamic propagation of allosteric signals to distal sites involves changes in long-range coupled domain motions on submicrosecond timescales, and that these coupled motions can be distinguished and characterized by NSE.     

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.     

Charge structure and counterion distribution in hexagonal DNA liquid crystal

A hexagonal liquid crystal of DNA fragments (double-stranded, 150 basepairs) with tetramethylammonium (TMA) counterions was investigated with small angle neutron scattering (SANS). We obtained the structure factors pertaining to the DNA and counterion density correlations with contrast matching in the water. Molecular dynamics (MD) computer simulation of a hexagonal assembly of nine DNA molecules showed that the inter-DNA distance fluctuates with a correlation time around 2 ns and a standard deviation of 8.5% of the interaxial spacing. The MD simulation also showed a minimal effect of the fluctuations in inter-DNA distance on the radial counterion density profile and significant penetration of the grooves by TMA. The radial density profile of the counterions was also obtained from a Monte Carlo (MC) computer simulation of a hexagonal array of charged rods with fixed interaxial spacing. Strong ordering of the counterions between the DNA molecules and the absence of charge fluctuations at longer wavelengths was shown by the SANS number and charge structure factors. The DNA-counterion and counterion structure factors are interpreted with the correlation functions derived from the Poisson-Boltzmann equation, MD, and MC simulation. Best agreement is observed between the experimental structure factors and the prediction based on the Poisson-Boltzmann equation and/or MC simulation. The SANS results show that TMA is too large to penetrate the grooves to a significant extent, in contrast to what is shown by MD simulation.     

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

Essential spaces defined by NMR structure ensembles and molecular dynamics simulation show significant overlap

Large concerted motions of proteins which span its “essential space,” are an important component of protein dynamics. We investigate to what extent structure ensembles generated with standard structure calculation techniques such as simulated annealing can capture these motions by comparing them to long-time molecular dynamics (MD) trajectories. The motions are analyzed by principal component analysis and compared using inner products of eigenvectors of the respective covariance matrices. Two very different systems are studied, the β-spectrin PH domain and the single-stranded DNA binding protein (ssDBP) from the filamentous phage Pf3. A comparison of the ensembles from NMR and MD shows significant overlap of the essential spaces, which in the case of ssDBP is extraordinarily high. The influence of variations in the specifications of distance restraints is investigated. We also study the influence of the selection criterion for the final structure ensemble on the definition of mobility. The results suggest a modified criterion that improves conformational sampling in terms of amplitudes of correlated motion. Proteins 31:370–382, 1998. © 1998 Wiley-Liss, Inc.     

Clathrin triskelia show evidence of molecular flexibility

The clathrin triskelion, which is a three-legged pinwheel-shaped heteropolymer, is a major component in the protein coats of certain post-Golgi and endocytic vesicles. At low pH, or at physiological pH in the presence of assembly proteins, triskelia will self-assemble to form a closed clathrin cage, or “basket”. Recent static light scattering and dynamic light scattering studies of triskelia in solution showed that an individual triskelion has an intrinsic pucker similar to, but differing from, that inferred from a high resolution cryoEM structure of a triskelion in a clathrin basket. We extend the earlier solution studies by performing small-angle neutron scattering (SANS) experiments on isolated triskelia, allowing us to examine a higher q range than that probed by static light scattering. Results of the SANS measurements are consistent with the light scattering measurements, but show a shoulder in the scattering function at intermediate q values (0.016 Å⁻¹), just beyond the Guinier regime. This feature can be accounted for by Brownian dynamics simulations based on flexible bead-spring models of a triskelion, which generate time-averaged scattering functions. Calculated scattering profiles are in good agreement with the experimental SANS profiles when the persistence length of the assumed semiflexible triskelion is close to that previously estimated from the analysis of electron micrographs.     

Measuring Material Microstructure Under Flow Using 1-2 Plane Flow-Small Angle Neutron Scattering

A new small-angle neutron scattering (SANS) sample environment optimized for studying the microstructure of complex fluids under simple shear flow is presented. The SANS shear cell consists of a concentric cylinder Couette geometry that is sealed and rotating about a horizontal axis so that the vorticity direction of the flow field is aligned with the neutron beam enabling scattering from the 1-2 plane of shear (velocity-velocity gradient, respectively). This approach is an advance over previous shear cell sample environments as there is a strong coupling between the bulk rheology and microstructural features in the 1-2 plane of shear. Flow-instabilities, such as shear banding, can also be studied by spatially resolved measurements. This is accomplished in this sample environment by using a narrow aperture for the neutron beam and scanning along the velocity gradient direction. Time resolved experiments, such as flow start-ups and large amplitude oscillatory shear flow are also possible by synchronization of the shear motion and time-resolved detection of scattered neutrons. Representative results using the methods outlined here demonstrate the useful nature of spatial resolution for measuring the microstructure of a wormlike micelle solution that exhibits shear banding, a phenomenon that can only be investigated by resolving the structure along the velocity gradient direction. Finally, potential improvements to the current design are discussed along with suggestions for supplementary experiments as motivation for future experiments on a broad range of complex fluids in a variety of shear motions.     

Molecular dynamics simulation of carboxy-myoglobin embedded in a trehalose-water matrix

We report on a molecular dynamics (MD) simulation of carboxy-myoglobin (MbCO) embedded in a water-trehalose system. The mean square fluctuations of protein atoms, calculated at different temperatures in the 100-300 K range, are compared with those from a previous MD simulation on an H2O-solvated MbCO and with experimental data from M?ssbauer spectroscopy and incoherent elastic neutron scattering on trehalose-coated MbCO. The results show that, for almost all the atomic classes, the amplitude of the nonharmonic motions stemming from the interconversion among the protein's conformational substates is reduced with respect to the H2O-solvated system, and their onset is shifted toward higher temperature. Moreover, our simulation shows that, at 300 K, the heme performs confined diffusive motions as a whole, leaving the underlying harmonic vibrations unaltered.     

Small angle neutron and X-ray scattering in structural biology: recent examples from the literature

Small angle scattering can provide unique structural information on the shape, domain organisation, and interactions of biomacromolecules in solution. Small angle neutron scattering (SANS) combined with deuterium labelling makes it possible to define the positions of specific components within a complex while small angle X-ray scattering (SAXS) provides more precise data on the overall shape. Here I review four recent publications, three of which were presented at the Neutrons in Biology meeting at the STFC Rutherford Appleton Laboratory in July 2007, that utilise SANS, SAXS, and complementary techniques to define the solution structure of large multidomain proteins and macromolecular complexes. These four papers emphasise the critical importance of sample quality and characterisation as well as the important role played by complementary techniques in building structural models based on small angle scattering data. They show the ability of SANS and SAXS in determining solution structures provides an important complementary structural technique for large, flexible, and glycosylated proteins where high resolution structural techniques, such as crystallography and NMR, cannot be applied.     

Interactions of the Anticancer Drug Tamoxifen with Lipid Membranes

Interactions of the hydrophobic anticancer drug tamoxifen (TAM) with lipid model membranes were studied using calcein-encapsulated vesicle leakage, attenuated total reflection Fourier transform infrared (FTIR) spectroscopy, small-angle neutron scattering (SANS), atomic force microscopy (AFM) based force spectroscopy, and all-atom molecular dynamics (MD) simulations. The addition of TAM enhances membrane permeability, inducing calcein to translocate from the interior to the exterior of lipid vesicles. A large decrease in the FTIR absorption band’s magnitude was observed in the hydrocarbon chain region, suggesting suppressed bond vibrational dynamics. Bilayer thickening was determined from SANS data. Force spectroscopy measurements indicate that the lipid bilayer area compressibility modulus K_A is increased by a large amount after the incorporation of TAM. MD simulations show that TAM decreases the lipid area and increases chain order parameters. Moreover, orientational and positional analyses show that TAM exhibits a highly dynamic conformation within the lipid bilayer. Our detailed experimental and computational studies of TAM interacting with model lipid membranes shed new light on membrane modulation by TAM.     

Kinetic asymmetry of subunit exchange of homooligomeric protein as revealed by deuteration-assisted small-angle neutron scattering

We developed a novel, to our knowledge, technique for real-time monitoring of subunit exchange in homooligomeric proteins, using deuteration-assisted small-angle neutron scattering (SANS), and applied it to the tetradecamer of the proteasome α7 subunit. Isotopically normal and deuterated tetradecamers exhibited identical SANS profiles in 81% D₂O solution. After mixing these solutions, the isotope sensitive SANS intensity in the low-q region gradually decreased, indicating subunit exchange, whereas the small-angle x-ray scattering profile remained unchanged confirming the structural integrity of the tetradecamer particles during the exchange. Kinetic analysis of zero-angle scattering intensity indicated that 1), only two of the 14 subunits were exchanged in each tetradecamer and 2), the exchange process involves at least two steps. This study underscores the usefulness of deuteration-assisted SANS, which can provide quantitative information not only on the molecular sizes and shapes of homooligomeric proteins, but also on their kinetic properties.     

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.     

Biophysical aspects of neutron scattering from vibrational modes of proteins

This review describes a major portion of the published work on neutron scattering experiments aimed at measuring large scale motions in proteins. The importance of these motions for enzyme function and oxygen transport is indicated. The theory applicable to each type of neutron scattering measurement is given and results are discussed with a view to biological relevance. New experiments are suggested and a comparison of neutron scattering data is made with results from other techniques such as raman scattering, infrared absorption, photolysis and molecular dynamics simulations.     

The dynamical Matryoshka model: 1. Incoherent neutron scattering functions for lipid dynamics in bilayers

Fluid lipid bilayers are the building blocks of biological membranes. Although there is a large amount of experimental data using incoherent quasi-elastic neutron scattering (QENS) techniques to study membranes, very little theoretical works have been developed to study the local dynamics of membranes. The main objective of this work is to build a theoretical framework to study and describe the local dynamics of lipids and derive analytical expressions of intermediate scattering functions (ISF) for QENS. As results, we developed the dynamical Matryoshka model which describes the local dynamics of lipid molecules in membrane layers as a nested hierarchical convolution of three motional processes: (i) individual motions described by the vibrational motions of H-atoms; (ii) internal motions including movements of the lipid backbone, head groups and tails, and (iii) molecule movements of the lipid molecule as a whole. The analytical expressions of the ISF associated with these movements are all derived. For use in analyzing the QENS experimental data, we also derived an analytical expression for the aggregate ISF of the Matryoshka model which involves an elastic term plus three inelastic terms of well-separated time scales and whose amplitudes and rates are functions of the lipid motions. And as an illustrative application, we used the aggregated ISF to analyze the experimental QENS data on a lipid sample of multilamellar bilayers of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine). It is clear from this analysis that the dynamical Matryoshka model describes very well the experimental data and allow extracting the dynamical parameters of the studied system.     

The helix-to-sheet transition of an HIV-1 fusion peptide derivative changes the mechanical properties of lipid bilayer membranes

A short sequence on the gp41 envelope protein of HIV-1 is integral to infection by the virus. Without this sequence, termed the fusion peptide (FP), the virus is far less effective at fusing with the cellular membrane. One of the interesting features of the isolated FP is that it transitions between an α-helical conformation and a β-sheet conformation in lipid bilayer membranes as a function of lipid composition and concentration, and the transition correlates with fusion. To better understand how the conformations of the FP impact lipid bilayer membranes, a variant of the FP that does not strongly promote fusion, termed gp41rk, was studied. Circular dichroism spectroscopy, dynamic light scattering, small-angle neutron scattering (SANS) and neutron spin echo spectroscopy (NSE) were used to relate the conformation of gp41rk to the structure and mechanical properties of lipid bilayer membrane vesicles composed of a 7:3 molar ratio mixture of 1,2-dimyristoyl-sn-glycero-3-phosphocholine and 1,2-dimyristoyl-sn-glycero-3-phospho-(1′-rac-glycerol). At a peptide-to-lipid ratio (P/L) of 1/200, it adopts an α-helical conformation, while gp41rk is a β-sheet at a P/L of 1/50 in the unilamellar vesicles. SANS reveals that the lipid bilayer membrane becomes thicker when gp41rk adopts a β-sheet conformation, which indicates that the high-concentration state of the peptide increases the order of the lipid acyl chains. At the same time, NSE demonstrates that the bilayer becomes more rigid, demonstrating that the β-sheet conformation, which correlates with fusion for the native FP sequence, stiffens the bilayer. The results have implications for the function of the FP.     

Protein dynamics in solution and powder measured by incoherent elastic neutron scattering: the influence of Q-range and energy resolution

Incoherent elastic neutron scattering (IENS) has been widely used to measure intramolecular atomic mean square displacements (MSDs) of proteins in powder and in solution. The instrumental energy resolution and the wave vector transfer (Q-range) determine, respectively, the time and length scales of observable motions. In order to investigate contributions of diffusive motions to MSDs measured by this method, we calculated the elastic intensity for several simple scattering functions. We showed that continuous translational diffusion contributes to MSDs in a Q-range where the energy width of the scattering function is of the order of the instrumental energy resolution. We discuss the choice of instruments adapted to focus on intramolecular motions in the presence of solvent or global macromolecular diffusion. The concepts developed are applied to interpret experimental data from H₂O- and D₂O-hydrated proteins. Finally, analogies between the Gaussian approximation in IENS and the Guinier approximation in small-angle scattering are discussed.