A solution NMR approach to the measurement of amphiphile immersion depth and orientation in membrane model systems

Solvent paramagnetic relaxation enhancement as a versatile method for studying structure and dynamics of biomolecular systems

Solvent paramagnetic relaxation enhancement (sPRE) is a versatile nuclear magnetic resonance (NMR)-based method that allows characterization of the structure and dynamics of biomolecular systems through providing quantitative experimental information on solvent accessibility of NMR-active nuclei. Addition of soluble paramagnetic probes to the solution of a biomolecule leads to paramagnetic relaxation enhancement in a concentration-dependent manner. Here we review recent progress in the sPRE-based characterization of structural and dynamic properties of biomolecules and their complexes, and aim to deliver a comprehensive illustration of a growing number of applications of the method to various biological systems. We discuss the physical principles of sPRE measurements and provide an overview of available co-solute paramagnetic probes. We then explore how sPRE, in combination with complementary biophysical techniques, can further advance biomolecular structure determination, identification of interaction surfaces within protein complexes, and probing of conformational changes and low-population transient states, as well as deliver insights into weak, nonspecific, and transient interactions between proteins and co-solutes. In addition, we present examples of how the incorporation of solvent paramagnetic probes can improve the sensitivity of NMR experiments and discuss the prospects of applying sPRE to NMR metabolomics, drug discovery, and the study of intrinsically disordered proteins.

Targeted Delivery of Adamantylated Peptidoglycan Immunomodulators in Lipid Nanocarriers: NMR Shows That Cargo Fragments Are Available on the Surface

We present an in-depth investigation of the membrane interactions of peptidoglycan (PGN)-based immune adjuvants designed for lipid-based delivery systems using NMR spectroscopy. The derivatives contain a cargo peptidoglycan (PGN) dipeptide fragment and an adamantyl group, which serves as an anchor to the lipid bilayer. Furthermore, derivatives with a mannose group that can actively target cell surface receptors on immune cells are also studied. We showed that the targeting mannose group and the cargo PGN fragment are both available on the lipid bilayer surface, thereby enabling interactions with cognate receptors. We found that the nonmannosylated compounds are incorporated stronger into the lipid assemblies than the mannosylated ones, but the latter compounds penetrate deeper in the bilayer. This might be explained by stronger electrostatic interactions available for zwitterionic nonmannosylated derivatives as opposed to the compounds in which the charged N-terminus is capped by mannose groups. The higher incorporation efficiency of the nonmannosylated compounds correlated with a larger relative enhancement in immune stimulation activities upon lipid incorporation compared to that of the derivatives with the mannose group. The chirality of the adamantyl group also influenced the incorporation efficiency, which in turn correlated with membrane-associated conformations that affect possible intermolecular interactions with lipid molecules. These findings will help in improving the development of PGN-based immune adjuvants suitable for delivery in lipid nanoparticles.

Biological membranes in EV biogenesis, stability, uptake, and cargo transfer: an ISEV position paper arising from the ISEV membranes and EVs workshop

Paracrine and endocrine roles have increasingly been ascribed to extracellular vesicles (EVs) generated by multicellular organisms. Central to the biogenesis, content, and function of EVs are their delimiting lipid bilayer membranes. To evaluate research progress on membranes and EVs, the International Society for Extracellular Vesicles (ISEV) conducted a workshop in March 2018 in Baltimore, Maryland, USA, bringing together key opinion leaders and hands-on researchers who were selected on the basis of submitted applications. The workshop was accompanied by two scientific surveys and covered four broad topics: EV biogenesis and release; EV uptake and fusion; technologies and strategies used to study EV membranes; and EV transfer and functional assays. In this ISEV position paper, we synthesize the results of the workshop and the related surveys to outline important outstanding questions about EV membranes and describe areas of consensus. The workshop discussions and survey responses reveal that while much progress has been made in the field, there are still several concepts that divide opinion. Good consensus exists in some areas, including particular aspects of EV biogenesis, uptake and downstream signalling. Areas with little to no consensus include EV storage and stability, as well as whether and how EVs fuse with target cells. Further research is needed in these key areas, as a better understanding of membrane biology will contribute substantially towards advancing the field of extracellular vesicles.

Polarity of Hydrated Phosphatidylcholine Headgroups

The headgroup (H) stratum (sometimes called the polar region) of membrane bilayers is a relevant yet poorly understood solvation phase for small molecules and macromolecules interacting with the membranes. Solvation of compounds in bilayer strata is characterized experimentally by wide- and small-angle X-ray scattering, neutron diffraction, and various NMR techniques. The quantification is tedious and only available for a limited set of small molecules. Our recently published model of liposome partitioning of small molecules shows that solvation of compounds in the H-stratum of fluid phosphatidylcholine (PC) bilayers correlates well with their solvation in hydrated diacetyl phosphatidylcholine (DAcPC), and solvation in the core (C) depends in a similar way on that in n-hexadecane. These two correlations became a basis for a model describing the location of compounds in the H- and C-strata and at the connecting interface as a nonlinear function of the fragment solvation characteristics of the compounds. In this study, refractivity of hydrated DAcPC phases with varying water contents was measured and polarity was determined using the steady-state fluorescence of indole and Nile Red. The results were compared with the published data obtained by other techniques for PC bilayers in liposomes or on solid supports. The demonstrated qualitative agreement, as well as the polarity and refractivity dependencies on the DAcPC concentration, supports the suitability of hydrated DAcPC as the H-stratum surrogate. Interestingly, depending on hydrations typical for the H-strata of fluid PC bilayers, the dielectric constant could decrease significantly from 31.0 to 7.3 for 16 and 8 water molecules per headgroup, respectively. Although additional experiments are needed for confirmation, this observation could help set proper dielectric constant magnitudes in continuum-based computational models of accumulation and crossing of the PC bilayers with varying hydration levels thanks to the temperature or the structure of fatty acid chains.

"Snorkelling" vs. "diving" in mixed micelles probed by means of a molecular bathymeter

A photoactive bathymeter based on a carboxylic acid moiety covalently linked to a signalling methoxynaphthalene (MNP) fluorophore has been designed to prove the concept of "snorkelling" vs. "diving" in mixed micelles (MM). The carboxylic acid "floats" on the MM surface, while the MNP unit sinks deep in MM. The rate constants of MNP fluorescence quenching by iodide, which remains basically in water, consistently decrease with increasing spacer length, revealing different regions. This is associated with the distance MNP should "dive" in MM to achieve protection from aqueous reactants. Unequivocal proof of the exergonic photoinduced electron transfer was obtained from the UV-visible spectral signature of I₃^- upon steady-state photolysis. The applicability of the bathymeter was examined upon testing a family of MNP derivatives. The obtained results were validated by comparison with different lipophilicity tests: (i) a modified version of the K_ow partition coefficient and (ii) the retention factor on thin layer chromatography. This concept could potentially be extended to test drugs or pharmacophores exhibiting any photoactive moiety.

Sensitivity enhancement for membrane proteins reconstituted in parallel and perpendicular oriented bicelles obtained by using repetitive cross-polarization and membrane-incorporated free radicals

Multidimensional separated local-field and spin-exchange experiments employed by oriented-sample solid-state NMR are essential for structure determination and spectroscopic assignment of membrane proteins reconstituted in macroscopically aligned lipid bilayers. However, these experiments typically require a large number of scans in order to establish interspin correlations. Here we have shown that a combination of optimized repetitive cross polarization (REP-CP) and membrane-embedded free radicals allows one to enhance the signal-to-noise ratio by factors 2.4-3.0 in the case of Pf1 coat protein reconstituted in magnetically aligned bicelles with their normals being either parallel or perpendicular to the main magnetic field. Notably, spectral resolution is not affected at the 2:1 radical-to-protein ratio. Spectroscopic assignment of Pf1 coat protein in the parallel bicelles has been established as an illustration of the method. The proposed methodology will advance applications of oriented-sample NMR technique when applied to samples containing smaller quantities of proteins and three-dimensional experiments.

Structure-based prediction of drug distribution across the headgroup and core strata of a phospholipid bilayer using surrogate phases

Solvation of drugs in the core (C) and headgroup (H) strata of phospholipid bilayers affects their physiological transport rates and accumulation. These characteristics, especially a complete drug distribution profile across the bilayer strata, are tedious to obtain experimentally, to the point that even simplified preferred locations are only available for a few dozen compounds. Recently, we showed that the partition coefficient (P) values in the system of hydrated diacetyl phosphatidylcholine (DAcPC) and n-hexadecane (C16), as surrogates of the H- and C-strata of the bilayer composed of the most abundant mammalian phospholipid, PC, agree well with the preferred bilayer location of compounds. High P values are typical for lipophiles accumulating in the core, and low P values are characteristic of cephalophiles preferring the headgroups. This simple pattern does not hold for most compounds, which usually have more even distribution and may also accumulate at the H/C interface. To model complete distribution, the correlates of solvation energies are needed for each drug state in the bilayer: (1) for the H-stratum it is the DAcPC/W P value, calculated as the ratio of the C16/W and C16/DAcPC (W for water) P values; (2) for the C-stratum, the C16/W P value; (3) for the H/C interface, the P values for all plausible molecular poses are characterized using the fragment DAcPC/W and C16/W solvation parameters for the parts of the molecule embedded in the H- and C-strata, respectively. The correlates, each scaled by two Collander coefficients, were used in a nonlinear, mass-balance based model of intrabilayer distribution, which was applied to the easily measurable overall P values of compounds in the DMPC (M = myristoyl) bilayers and monolayers as the dependent variables. The calibrated model for 107 neutral compounds explains 94% of experimental variance, achieves similar cross-validation levels, and agrees well with the nontrivial, experimentally determined bilayer locations for 27 compounds. The resulting structure-based prediction system for intrabilayer distribution will facilitate more realistic modeling of passive transport and drug interactions with those integral membrane proteins, which have the binding sites located in the bilayer, such as some enzymes, influx and efflux transporters, and receptors. If only overall bilayer accumulation is of interest, the 1-octanol/W P values suffice to model the studied set.

Measuring membrane penetration with depth-dependent fluorescence quenching: distribution analysis is coming of age

Depth-dependent fluorescence quenching by lipid-attached quenchers (e.g., bromine atoms and doxyl groups) is an important tool for determining the penetration of proteins and peptides into lipid bilayers. Extracting quantitative information and accurate calculations of the depth of the fluorophore are complicated by thermal disorder, resulting in broad distributions of the transverse positions of both quenchers and fluorophores. Twenty-one years ago a methodology called distribution analysis (DA) was introduced, based on the emerging view of the complexity of the transverse organization of lipid bilayer structure. The method is aimed at extracting quantitative information on membrane penetration, such as position and width of fluorophore's distribution along the depth coordinate and its exposure to the lipid phase. Here we review recent progress in refining the DA method and illustrate its applications to protein-membrane interactions. We demonstrate how basic assumptions of the DA approach can be validated using molecular dynamics simulations and how the precision of depth determination is improved by applying a new protocol based on a combination of steady-state and time-resolved fluorescence quenching. Using the example of the MPER fragment of the membrane-spanning domain of the HIV-1 gp41 fusion protein, we illustrate how DA applications and computer simulations can be used together to reveal the molecular organization of a protein-membrane complex. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Sensitivity enhancement and contrasting information provided by free radicals in oriented-sample NMR of bicelle-reconstituted membrane proteins

Elucidating structure and topology of membrane proteins (MPs) is essential for unveiling functionality of these important biological constituents. Oriented-sample solid-state NMR (OS-NMR) is capable of providing such information on MPs under nearly physiological conditions. However, two dimensional OS-NMR experiments can take several days to complete due to long longitudinal relaxation times combined with the large number of scans to achieve sufficient signal sensitivity in biological samples. Here, free radicals 5-DOXYL stearic acid, TEMPOL, and CAT-1 were added to uniformly (15)N-labeled Pf1 coat protein reconstituted in DMPC/DHPC bicelles, and their effect on the longitudinal relaxation times (T1Z) was investigated. The dramatically shortened T1Z's allowed for the signal gain per unit time to be used for either: (i) up to a threefold reduction of the total experimental time at 99% magnetization recovery or (ii) obtaining up to 74% signal enhancement between the control and radical samples during constant experimental time at "optimal" relaxation delays. In addition, through OS-NMR and high-field EPR studies, free radicals were able to provide positional constraints in the bicelle system, which provide a description of the location of each residue in Pf1 coat protein within the bicellar membranes. This information can be useful in the determination of oligomerization states and immersion depths of larger membrane proteins.

Studying the structure and dynamics of biomolecules by using soluble paramagnetic probes

Characterisation of the structure and dynamics of large biomolecules and biomolecular complexes by NMR spectroscopy is hampered by increasing overlap and severe broadening of NMR signals. As a consequence, the number of available NMR spectroscopy data is often sparse and new approaches to provide complementary NMR spectroscopy data are needed. Paramagnetic relaxation enhancements (PREs) obtained from inert and soluble paramagnetic probes (solvent PREs) provide detailed quantitative information about the solvent accessibility of NMR-active nuclei. Solvent PREs can be easily measured without modification of the biomolecule; are sensitive to molecular structure and dynamics; and are therefore becoming increasingly powerful for the study of biomolecules, such as proteins, nucleic acids, ligands and their complexes in solution. In this Minireview, we give an overview of the available solvent PRE probes and discuss their applications for structural and dynamic characterisation of biomolecules and biomolecular complexes.

Molecular dynamics simulations of depth distribution of spin-labeled phospholipids within lipid bilayer

Spin-labeled lipids are commonly used as fluorescence quenchers in studies of membrane penetration of dye-labeled proteins and peptides using depth-dependent quenching. Accurate calculations of depth of the fluorophore rely on the use of several spin labels placed in the membrane at various positions. The depth of the quenchers (spin probes) has to be determined independently; however, experimental determination of transverse distributions of spin probe depths is difficult. In this Article, we use molecular dynamics (MD) simulations to study the membrane behavior and depth distributions of spin-labeled phospholipids in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. To probe different depths within the bilayer, a series containing five Doxyl-labeled lipids (n-Doxyl PC) has been studied, in which a spin moiety was covalently attached to nth carbon atoms (where n = 5, 7, 10, 12, and 14) of the sn-2 stearoyl chain of the host phospholipid. Our results demonstrate that the chain-attached spin labels are broadly distributed across the model membrane and their environment is characterized by a high degree of mobility and structural heterogeneity. Despite the high thermal disorder, the depth distributions of the Doxyl labels were found to correlate well with their attachment positions, indicating that the distribution of the spin label within the model membrane is dictated by the depth of the nth lipid carbon atom and not by intrinsic properties of the label. In contrast, a much broader and heterogeneous distribution was observed for a headgroup-attached Tempo spin label of Tempo-PC lipids. MD simulations reveal that, due to the hydrophobic nature, a Tempo moiety favors partitioning from the headgroup region deeper into the membrane. Depending on the concentration of Tempo-PC lipids, the probable depth of the Tempo moiety could span a range from 14.4 to 18.2 Å from the membrane center. Comparison of the MD-estimated immersion depths of Tempo and n-Doxyl labels with their suggested experimental depth positions allows us to review critically the possible sources of error in depth-dependent fluorescence quenching studies.

Topology and immersion depth of an integral membrane protein by paramagnetic rates from dissolved oxygen

In studies of membrane proteins, knowledge of protein topology can provide useful insight into both structure and function. In this work, we present a solution NMR method for the measurement the tilt angle and average immersion depth of alpha helices in membrane proteins, from analysis of the paramagnetic relaxation rate enhancements arising from dissolved oxygen. No modification to the micelle or protein is necessary, and the topology of both transmembrane and amphipathic helices are readily determined. We apply this method to the measure the topology of a monomeric mutant of phospholamban (AFA-PLN), a 52-residue membrane protein containing both an amphipathic and a transmembrane alpha helix. In dodecylphosphocholine micelles, the amphipathic helix of AFA-PLN was found to have a tilt angle of 87° ± 1° and an average immersion depth of 13.2 Å. The transmembrane helix was found to have an average immersion depth of 5.4 Å, indicating residues 41 and 42 are closest to the micelle centre. The resolution of paramagnetic relaxation rate enhancements from dissolved oxygen compares favourably to those from Ni (II), a hydrophilic paramagnetic species.

Temperature driven annealing of perforations in bicellar model membranes

Bicellar model membranes composed of 1,2-dimyristoylphosphatidylcholine (DMPC) and 1,2-dihexanoylphosphatidylcholine (DHPC), with a DMPC/DHPC molar ratio of 5, and doped with the negatively charged lipid 1,2-dimyristoylphosphatidylglycerol (DMPG), at DMPG/DMPC molar ratios of 0.02 or 0.1, were examined using small angle neutron scattering (SANS), (31)P NMR, and (1)H pulsed field gradient (PFG) diffusion NMR with the goal of understanding temperature effects on the DHPC-dependent perforations in these self-assembled membrane mimetics. Over the temperature range studied via SANS (300-330 K), these bicellar lipid mixtures exhibited a well-ordered lamellar phase. The interlamellar spacing d increased with increasing temperature, in direct contrast to the decrease in d observed upon increasing temperature with otherwise identical lipid mixtures lacking DHPC. (31)P NMR measurements on magnetically aligned bicellar mixtures of identical composition indicated a progressive migration of DHPC from regions of high curvature into planar regions with increasing temperature, and in accord with the "mixed bicelle model" (Triba, M. N.; Warschawski, D. E.; Devaux, P. E. Biophys. J.2005, 88, 1887-1901). Parallel PFG diffusion NMR measurements of transbilayer water diffusion, where the observed diffusion is dependent on the fractional surface area of lamellar perforations, showed that transbilayer water diffusion decreased with increasing temperature. A model is proposed consistent with the SANS, (31)P NMR, and PFG diffusion NMR data, wherein increasing temperature drives the progressive migration of DHPC out of high-curvature regions, consequently decreasing the fractional volume of lamellar perforations, so that water occupying these perforations redistributes into the interlamellar volume, thereby increasing the interlamellar spacing.

Anisotropic solvent model of the lipid bilayer. 2. Energetics of insertion of small molecules, peptides, and proteins in membranes

A new computational approach to calculating binding energies and spatial positions of small molecules, peptides, and proteins in the lipid bilayer has been developed. The method combines an anisotropic solvent representation of the lipid bilayer and universal solvation model, which predicts transfer energies of molecules from water to an arbitrary medium with defined polarity properties. The universal solvation model accounts for hydrophobic, van der Waals, hydrogen-bonding, and electrostatic solute-solvent interactions. The lipid bilayer is represented as a fluid anisotropic environment described by profiles of dielectric constant (ε), solvatochromic dipolarity parameter (π*), and hydrogen bonding acidity and basicity parameters (α and β). The polarity profiles were calculated using published distributions of quasi-molecular segments of lipids determined by neutron and X-ray scattering for DOPC bilayer and spin-labeling data that define concentration of water in the lipid acyl chain region. The model also accounts for the preferential solvation of charges and polar groups by water and includes the effect of the hydrophobic mismatch for transmembrane proteins. The method was tested on calculations of binding energies and preferential positions in membranes for small-molecules, peptides and peripheral membrane proteins that have been experimentally studied. The new theoretical approach was implemented in a new version (2.0) of our PPM program and applied for the large-scale calculations of spatial positions in membranes of more than 1000 peripheral and integral proteins. The results of calculations are deposited in the updated OPM database ( http://opm.phar.umich.edu ).

Solid-state NMR paramagnetic relaxation enhancement immersion depth studies in phospholipid bilayers

A new approach for determining the membrane immersion depth of a spin-labeled probe has been developed using paramagnetic relaxation enhancement (PRE) in solid-state NMR spectroscopy. A DOXYL spin label was placed at different sites of 1-palmitoyl-2-stearoyl-sn-glycero-3-phosphocholine (PSPC) phospholipid bilayers as paramagnetic moieties and the resulting enhancements of the longitudinal relaxation (T₁) times of ³¹P nuclei on the surface of the bilayers were measured by a standard inversion recovery pulse sequence. The ³¹P NMR spin-lattice relaxation times decrease steadily as the DOXYL spin label moves closer to the surface as well as the concentration of the spin-labeled lipids increase. The enhanced relaxation vs. the position and concentration of spin-labels indicate that PRE induced by the DOXYL spin label are significant to determine longer distances over the whole range of the membrane depths. When these data were combined with estimated correlation times τ(c), the r⁻⁶-weighted, time-averaged distances between the spin-labels and the ³¹P nuclei on the membrane surface were estimated. The application of using this solid-state NMR PRE approach coupled with site-directed spin labeling (SDSL) may be a powerful method for measuring membrane protein immersion depth.

Presentation of membrane-anchored glycosphingolipids determined from molecular dynamics simulations and NMR paramagnetic relaxation rate enhancement

Challenges for structural characterization of membrane-bound glycosphingolipids include their high internal dynamic motions and their physical proximity to membrane surfaces. Here we demonstrate that NMR paramagnetic relaxation enhancement can be used, alongside independent molecular dynamics simulations and an outer-sphere relaxation model, to quantitatively characterize the presentation (insertion depth and orientation relative to a membrane surface) of ganglioside GM1 in biologically relevant membrane environments. Longitudinal and transverse paramagnetic relaxation enhancement effects were measured for GM1, anchored to phospholipid bicelles, using both water-soluble and membrane-anchored paramagnetic probes, respectively. A method was developed to rapidly calculate paramagnetic relaxation enhancement effects from thousands of structures taken from a simulation of GM1 in a phospholipid bilayer. The combined computational and experimental approach yielded experimentally verified atomic-resolution 3D models of a highly plastic membrane-bound biomolecule.

The impact of window functions on NMR-based paramagnetic relaxation enhancement measurements in membrane proteins

Though challenging, solution NMR spectroscopy allows fundamental interrogation of the structure and dynamics of membrane proteins. One major technical hurdle in studies of helical membrane proteins by NMR is the difficulty of obtaining sufficient long range NOEs to determine tertiary structure. For this reason, long range distance information is sometimes sought through measurement of paramagnetic relaxation enhancements (PRE) of NMR nuclei as a function of distance from an introduced paramagnetic probe. Current PRE interpretation is based on the assumption of Lorentzian resonance lineshapes. However, in order to optimize spectral resolution, modern multidimensional NMR spectra are almost always subjected to resolution-enhancement, leading to distortions in the Lorentizian peak shape. Here it is shown that when PREs are derived using peak intensities (i.e., peak height) and linewidths from both real and simulated spectra that were produced using a wide range of apodization/window functions, that there is little variation in the distances determined (<1 A at the extremes). This indicates that the high degree of resolution enhancement required to obtain well-resolved spectra from helical membrane proteins is compatible with the use of PRE data as a source of distance restraints. While these conclusions are particularly important for helical membrane proteins, they are generally applicable to all PRE measurements made using resolution-enhanced data.