Synthesis,Characterization and Applications of a Perdeuterated Amphipol

Amphipols are short amphipathic polymers that can substitute for detergents at the hydrophobic surface of membrane proteins (MPs), keeping them soluble in the absence of detergents while stabilizing them. The most widely used amphipol, known as A8-35, is comprised of a polyacrylic acid (PAA) main chain grafted with octylamine and isopropylamine. Among its many applications, A8-35 has proven particularly useful for solution-state NMR studies of MPs, for which it can be desirable to eliminate signals originating from the protons of the surfactant. In the present work, we describe the synthesis and properties of perdeuterated A8-35 (perDAPol). Perdeuterated PAA was obtained by radical polymerization of deuterated acrylic acid. It was subsequently grafted with deuterated amines, yielding perDAPol. The number-average molar mass of hydrogenated and perDAPol, ~4 and ~5 kDa, respectively, was deduced from that of their PAA precursors, determined by size exclusion chromatography in tetrahydrofuran following permethylation. Electrospray ionization–ion mobility spectrometry–mass spectrometry measurements show the molar mass and distribution of the two APols to be very similar. Upon neutron scattering, the contrast match point of perDAPol is found to be ~120 % D₂O. In ¹H-¹H nuclear overhauser effect NMR spectra, its contribution is reduced to ~6 % of that of hydrogenated A8-35, making it suitable for extended uses in NMR spectroscopy. PerDAPol ought to also be of use for inelastic neutron scattering studies of the dynamics of APol-trapped MPs, as well as small-angle neutron scattering and analytical ultracentrifugation.     

Folding and stability of outer membrane protein A (OmpA) from Escherichia coli in an amphipathic polymer, amphipol A8-35

Amphipols are a class of amphipathic polymers designed to maintain membrane proteins in aqueous solutions in the absence of detergents. Denatured β-barrel membrane proteins, like outer membrane proteins OmpA from Escherichia coli and FomA from Fusobacterium nucleatum, can be folded by dilution of the denaturant urea in the presence of amphipol A8-35. Here, the folding kinetics and stability of OmpA in A8-35 have been investigated. Folding is well described by two parallel first-order processes, whose half-times, ~5 and ~70 min, respectively, are independent of A8-35 concentration. The faster process contributed ~55–64 % to OmpA folding. Folding into A8-35 was faster than into dioleoylphosphatidylcholine bilayers and complete at ratios as low as ~0.17 g/g A8-35/OmpA, corresponding to ~1–2 A8-35 molecules per OmpA. Activation energies were determined from the temperature dependence of folding kinetics, monitored both by electrophoresis, which reports on the formation of stable OmpA tertiary structure, and by fluorescence spectroscopy, which reflects changes in the environment of tryptophan side chains. The two methods yielded consistent estimates, namely ~5–9 kJ/mol for the fast process and ~29–37 kJ/mol for the slow one, which is lower than is observed for OmpA folding into dioleoylphosphatidylcholine bilayers. Folding and unfolding titrations with urea demonstrated that OmpA folding into A8-35 is reversible and that amphipol-refolded OmpA is thermodynamically stable at room temperature. Comparison of activation energies for folding and unfolding in A8-35 versus detergent indicates that stabilization of A8-35-trapped OmpA against denaturation by urea is a kinetic, not a thermodynamic phenomenon.     

Molecular Dynamics Simulations of a Membrane Protein/Amphipol Complex

Amphipathic polymers known as “amphipols” provide a highly stabilizing environment for handling membrane proteins in aqueous solutions. A8-35, an amphipol with a polyacrylate backbone and hydrophobic grafts, has been extensively characterized and widely employed for structural and functional studies of membrane proteins using biochemical and biophysical approaches. Given the sensitivity of membrane proteins to their environment, it is important to examine what effects amphipols may have on the structure and dynamics of the proteins they complex. Here we present the first molecular dynamics study of an amphipol-stabilized membrane protein, using Escherichia coli OmpX as a model. We begin by describing the structure of the complexes formed by supplementing OmpX with increasing amounts of A8-35, in order to determine how the amphipol interacts with the transmembrane and extramembrane surfaces of the protein. We then compare the dynamics of the protein in either A8-35, a detergent, or a lipid bilayer. We find that protein dynamics on all accessible length scales is restrained by A8-35, which provides a basis to understanding some of the stabilizing and functional effects of amphipols that have been experimentally observed.     

A reinterpretation of neutron scattering experiments on a lipidated Ras peptide using replica exchange molecular dynamics

The Ras family of proteins plays crucial roles in a variety of cell signaling networks where they have the function of a molecular switch. Their particular medical relevance arises from mutations in these proteins that are implicated in ~30% of human cancers. The various Ras proteins exhibit a high degree of homology in their soluble domains but extremely high variability in the membrane anchoring regions that are crucial for protein function and are the focus of this study. We have employed replica exchange molecular dynamics computer simulations to study a doubly lipidated heptapeptide, corresponding to the C-terminus of the human N-Ras protein, incorporated into a dimyristoylphosphatidylcholine lipid bilayer. This same system has previously been investigated experimentally utilizing a number of techniques, including neutron scattering. Here we present results of well converged simulations that describe the subtle changes in scattering density in terms of the location of the peptide and its lipid modifications and in terms of changes in phospholipid density arising from the incorporation of the peptide into the membrane bilayer. The detailed picture that emerges from the combination of experimental and computational data exemplifies the power of combining isotopic substitution neutron scattering with atomistic molecular dynamics simulation. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Subunit Stoichiometry of the Clostridium botulinum Type A Neurotoxin Complex Determined Using Denaturing Capillary Electrophoresis

A denaturing capillary electrophoresis method was developed to evaluate the subunit stoichiometry of the Clostridium botulinum type A neurotoxin complex. The results indicate that the neurotoxin complex contains single copies of the 150 kDa neurotoxin and the non-toxic non-hemagglutinating subunits as well as 5–6 HA17, 4–5 HA23, 3–4 HA48, and 8–9 HA34 subunits. The calculated molecular mass for a complex with this stoichiometry is between 880 and 1,000 kDa. The molecular mass of the intact complex was determined using size-exclusion HPLC (SE-HPLC) and SE-HPLC in conjunction with multi-angle laser light scattering detection. Based on a comparison to a mixture of standard proteins, SE-HPLC analysis yielded a molecular mass of 880 kDa while light scattering analysis indicated a weight average molecular mass of 925 ± 45 kDa. The close agreement between the molecular mass values determined by the three approaches supports the subunit stoichiometry proposed for the C. botulinum type A neurotoxin complex.     

Fundamental and biotechnological applications of neutron scattering measurements for macromolecular dynamics

To explore macromolecular dynamics on the picosecond timescale, we used neutron spectroscopy. First, molecular dynamics were analyzed for the hyperthermophile malate dehydrogenase from Methanococcus jannaschii and a mesophilic homologue, the lactate dehydrogenase from Oryctolagus cunniculus muscle. Hyperthermophiles have elaborate molecular mechanisms of adaptation to extremely high temperature. Using a novel elastic neutron scattering approach that provides independent measurements of the global flexibility and of the structural resilience (rigidity), we have demonstrated that macromolecular dynamics represents one of these molecular mechanisms of thermoadaptation. The flexibilities were found to be similar for both enzymes at their optimal activity temperature and the resilience is higher for the hyperthermophilic protein. Secondly, macromolecular motions were examined in a native and immobilized dihydrofolate reductase (DHFR) from Escherichia coli. The immobilized mesophilic enzyme has increased stability and decreased activity, so that its properties are changed to resemble those of the thermophilic enzyme. Are these changes reflected in dynamical behavior? For this study, we performed quasielastic neutron scattering measurements to probe the protein motions. The residence time is 7.95 ps for the native DHFR and 20.36 ps for the immobilized DHFR. The average height of the potential barrier to local motions is therefore increased in the immobilized DHFR, with a difference in activation energy equal to 0.54 kcal/mol, which is, using the theoretical rate equation, of the same order than expected from calculation.     

Low-temperature molecular dynamics simulations of horse heart cytochrome c and comparison with inelastic neutron scattering data

Molecular dynamics (MD) simulation combined with inelastic neutron scattering can provide information about the thermal dynamics of proteins, especially the low-frequency vibrational modes responsible for large movement of some parts of protein molecules. We performed several 30-ns MD simulations of cytochrome c (Cyt c) in a water box for temperatures ranging from 110 to 300 K and compared the results with those from experimental inelastic neutron scattering. The low-frequency vibrational modes were obtained via dynamic structure factors, S(Q, ω), obtained both from inelastic neutron scattering experiments and calculated from MD simulations for Cyt c in the same range of temperatures. The well known thermal transition in structural movements of Cyt c is clearly seen in MD simulations; it is, however, confined to unstructured fragments of loops Ω₁ and Ω₂; movement of structured loop Ω₃ and both helical ends of the protein is resistant to thermal disturbance. Calculated and experimental S(Q, ω) plots are in qualitative agreement for low temperatures whereas above 200 K a boson peak vanishes from the calculated plots. This may be a result of loss of crystal structure by the protein–water system compared with the protein crystal.     

Temperature dependence of the generalized frequency distribution of water molecules: comparison of experiments and molecular dynamics simulations

The results of inelastic neutron scattering experiments on water in the temperature interval 300–623 K along the coexistence curve are compared with data obtained from molecular dynamics simulations. In general, a good agreement between experiments and calculations is observed and it serves as a satisfactory test of the potential models employed. The temperature dependence of the generalized frequency distribution of water molecules obtained by both experiment and computer simulation demonstrates the accordance with the temperature evolution of the water structure, extracted from neutron and X-ray diffraction measurements.     

Structure of the K-turn U4 RNA: a combined NMR and SANS study

K-turn motifs are universal RNA structural elements providing a binding platform for proteins in several cellular contexts. Their characteristic is a sharp kink in the phosphate backbone that puts the two helical stems of the protein-bound RNA at an angle of 60°. However, to date no high-resolution structure of a naked K-turn motif is available. Here, we present the first structural investigation at atomic resolution of an unbound K-turn RNA (the spliceosomal U4-Kt RNA) by a combination of NMR and small-angle neutron scattering data. With this study, we wish to address the question whether the K-turn structural motif assumes the sharply kinked conformation in the absence of protein binders and divalent cations. Previous studies have addressed this question by fluorescence resonance energy transfer, biochemical assays and molecular dynamics simulations, suggesting that the K-turn RNAs exist in equilibrium between a kinked conformation, which is competent for protein binding, and a more extended conformation, with the population distribution depending on the concentration of divalent cations. Our data shows that the U4-Kt RNA predominantly assumes the more extended conformation in the absence of proteins and divalent cations. The internal loop region is well structured but adopts a different conformation from the one observed in complex with proteins. Our data suggests that the K-turn consensus sequence does not per se code for the kinked conformation; instead the sharp backbone kink requires to be stabilized by protein binders.     

Small angle neutron scattering on an absolute intensity scale and the internal surface of diatom frustules from three species of differing morphologies

The internal nanostructure of the diatoms Cyclotella meneghiniana, Seminavis robusta and Achnanthes subsessilis was investigated using small angle neutron scattering (SANS) to examine thin biosilica samples, consisting of isotropic (powder) from their isolated cell walls. The interpretation of SANS data was assisted by several other measurements. The N₂ adsorption, interpreted within the Branuer–Emmet–Teller isotherm, yielded the specific surface area of the material. Fourier transform infrared (FTIR) and Raman spectroscopy indicates that the isolated material is amorphous silica with small amounts of organic cell wall materials acting as a filling material between the silica particles. A two-phase (air and amorphous silica) model was used to interpret small angle neutron scattering data. After correction for instrumental resolution, the measurements on two SANS instruments covered an extended range of scattering vectors 0.0011 nm^?1 < q < 5.6 nm^?1, giving an almost continuous SANS curve over a range of scattering vectors, q, on an absolute scale of intensity for each sample. Each of the samples gave a characteristic scattering curve where log (intensity) versus log (q) has a ?4 dependence, with other features superimposed. In the high-q regime, departure from this behaviour was observed at a length-scales equivalent to the proposed unitary silica particle. The limiting Porod scattering law was used to determine the specific area per unit of volume of each sample illuminated by the neutron beam. The Porod behaviour, and divergence from this behaviour, is discussed in terms of various structural features and the proposed mechanisms for the bio-assembly of unitary silica particles in frustules.     

Structure Determination of Functional Membrane Proteins using Small-Angle Neutron Scattering (SANS) with Small, Mixed-Lipid Liposomes: Native Beef Heart Mitochondrial Cytochrome c Oxidase Forms Dimers

The low-resolution three-dimensional structure of purified native beef heart mitochondrial cytochrome c oxidase (COX) in asolectin unilamellar liposomes has been measured by small-angle neutron scattering under the conditions where the protein remains fully functional. From a neutron scattering perspective, the use of mixed-lipid liposomes provided for a more homogeneous matrix than can be achieved using a single lipid. As a result, the measurements were able to be performed under conditions where the liposome scattering was essentially eliminated (contrast-matched conditions). The protein structure in the membrane was modeled as a simple parallelepiped with side lengths of (59 × 70 × 120) Å with uncertainties, respectively, (11, 12, 20 Å). The molecular mass calculated for a typical protein with this volume is estimated to be (410 ± 124) kDa, which indicates the mass of a COX dimer. The longest dimension has some uncertainty due to intermolecular scattering contributing to the data. Nevertheless, that length was estimated using an average protein density and the known dimer molecular mass. Using the same cross sectional dimensions for the structure, the length is estimated to be 120 Å. However, the measured scattering curve of the dimer in the liposome differs significantly from that calculated from the X-ray structure of the dimer in a crystal of mixed micelles (PDB 3AG1). The calculated SANS scattering from the crystal structure was fit with a parallelepiped, measuring (59 × 101 × 129) Å with fitting uncertainties, respectively, (2, 3, 3 Å). Our results suggest that COX is a functional dimer when reconstituted into mixed-lipid liposomes.     

Structure and composition of influenza virus. A small-angle neutron scattering study

A detailed analysis is presented of the small-angle neutron scattering curves of homogeneous solutions of influenza B virus, both intact and after treatment with bromelain, which removes the external glycoprotein spikes. The two sets of data are consistent with the following low-resolution structure: the virus particles are spherical, about 1200 A in diameter and of Mr about 180 X 10(6). The lipid bilayer is centred at a radius of 425 A, is 40 A to 50 A thick and constitutes 25% to 28% of the virus mass. The surface glycoproteins, predominantly haemagglutinin, contribute 40% to 46% of the total mass. Surprisingly little protein is found in the interior of the virus. It is suggested that the reason for this is that many particles do not contain the full complement of ribonucleoprotein complexes. These results are in good agreement with recent scanning transmission electron microscopic measurements of molecular mass and cryo-electron microscopic observations of the same preparations. Appendix 1 describes a new method of deriving spherical shell models from contrast variation neutron scattering data on viruses, in which scattering curves from all measured contrasts are used simultaneously. There is also a discussion of the assumptions and limitations implicit in the structural interpretation of such models, with emphasis on viruses containing lipid bilayers. Appendix 2 examines the effect on the scattering curves of various arrangements of the surface glycoproteins.     

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.     

Structural dynamics of protein backbone φ angles: extended molecular dynamics simulations versus experimental 3 J scalar couplings

³ J scalar couplings report on the conformational averaging of backbone φ angles in peptides and proteins, and therefore represent a potentially powerful tool for studying the details of both structure and dynamics in solution. We have compared an extensive experimental dataset with J-couplings predicted from unrestrained molecular dynamics simulation using enhanced sampling available from accelerated molecular dynamics or using long timescale trajectories (200 ns). The dynamic fluctuations predicted to be present along the backbone, in agreement with residual dipolar coupling analysis, are compatible with the experimental ³ J scalar couplings providing a slightly better reproduction of these experimental parameters than a high-resolution static structure.     

Application of Amphipols for Structure–Functional Analysis of TRP Channels

Amphipathic polymers (amphipols), such as A8-35 and SApol, are a new tool for stabilizing integral membrane proteins in detergent-free conditions for structural and functional studies. Transient receptor potential (TRP) ion channels function as tetrameric protein complexes in a diverse range of cellular processes including sensory transduction. Mammalian TRP channels share ~20 % sequence similarity and are categorized into six subfamilies: TRPC (canonical), TRPV (vanilloid), TRPA (ankyrin), TRPM (melastatin), TRPP (polycystin), and TRPML (mucolipin). Due to the inherent difficulties in purifying eukaryotic membrane proteins, structural studies of TRP channels have been limited. Recently, A8-35 was essential in resolving the molecular architecture of the nociceptor TRPA1 and led to the determination of a high-resolution structure of the thermosensitive TRPV1 channel by cryo-EM. Newly developed maltose-neopentyl glycol (MNG) detergents have also proven to be useful in stabilizing TRP channels for structural analysis. In this review, we will discuss the impacts of amphipols and MNG detergents on structural studies of TRP channels by cryo-EM. We will compare how A8-35 and MNG detergents interact with the hydrophobic transmembrane domains of TRP channels. In addition, we will discuss what these cryo-EM studies reveal on the importance of screening different types of surfactants toward determining high-resolution structures of TRP channels.     

Self-similar dynamics of proteins under hydrostatic pressure—Computer simulations and experiments

Different experimental techniques, such as kinetic studies of ligand binding and fluorescence correlation spectroscopy, have revealed that the diffusive, internal dynamics of proteins exhibits autosimilarity on the time scale from microseconds to hours. Computer simulations have demonstrated that this type of dynamics is already established on the much shorter nanosecond time scale, which is also covered by quasielastic neutron scattering experiments. The autosimilarity of protein dynamics is reflected in long-time memory effects in the underlying diffusion processes, which lead to a non-exponential decay of the observed time correlation functions. Fractional Brownian dynamics is an empirical model which is able to capture the essential aspects of internal protein dynamics. Here we give a brief introduction into the theory and show how the model can be used to interpret neutron scattering experiments and molecular dynamics simulation of proteins in solution under hydrostatic pressure.     

Time-resolved quasielastic neutron scattering studies of native photosystems

The internal molecular dynamics of proteins plays an important role in a number of functional processes in native photosystems. Prominent examples include the photocycle of bacteriorhodopsin and electron transfer in the reaction center of plant photosystem II. In this regard, the recently developed technique of time-resolved quasielastic neutron scattering with laser excitation opens up new perspectives for the study of protein/membrane dynamics in specific functional states of even complex systems. The first direct observation of a functionally modulated protein dynamics has just recently been reported for the model system bacteriorhodopsin (Pieper et al., Phys. Rev. Lett. 100, 2008, 228103.), where a transient softening of the protein was observed on a timescale of ∼ 1 ms along with the large-scale structural change in the M-intermediate of bacteriorhodopsin. In contrast, photosystem II membrane fragments with inhibited electron transfer show a suppression of protein dynamics ∼160 μs after the actinic laser flash (Pieper and Renger, Biochemistry 48, 2009, 6111). This effect may reflect aggregation-like conformational changes capable of dissipation of excess excitation energy to prevent photodamage in the absence of Q_A→Q_B electron transfer. These findings indicate that proteins exhibit a remarkable flexibility to accommodate different functional processes. This contribution will discuss methodical aspects, challenges, and recent applications of laser-excited, time-resolved quasielastic neutron scattering.     

Frequent side chain methyl carbon‐oxygen hydrogen bonding in proteins revealed by computational and stereochemical analysis of neutron structures

The propensity of backbone Cα atoms to engage in carbon‐oxygen (CH···O) hydrogen bonding is well‐appreciated in protein structure, but side chain CH···O hydrogen bonding remains largely uncharacterized. The extent to which side chain methyl groups in proteins participate in CH···O hydrogen bonding is examined through a survey of neutron crystal structures, quantum chemistry calculations, and molecular dynamics simulations. Using these approaches, methyl groups were observed to form stabilizing CH···O hydrogen bonds within protein structure that are maintained through protein dynamics and participate in correlated motion. Collectively, these findings illustrate that side chain methyl CH···O hydrogen bonding contributes to the energetics of protein structure and folding. Proteins 2015; 83:403–410. © 2014 Wiley Periodicals, Inc.