Lipid Dynamics Studied by Calculation of 31P Solid-State NMR Spectra Using Ensembles from Molecular Dynamics Simulations

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Single-crystal NMR spectroscopy

Single-crystal (SC) NMR spectroscopy is a solid-state NMR method that has been used since the early days of NMR to study the magnitude and orientation of tensorial nuclear spin interactions in solids. This review first presents the field of SC NMR instrumentation, then provides a survey of software for analysis of SC NMR data, and finally it highlights selected applications of SC NMR in various fields of research. The aim of the last part is not to provide a complete review of all SC NMR literature but to provide examples that demonstrate interesting applications of SC NMR.

Structure of Phospholipid Mixed Micelles (Bicelles) Studied by Small-Angle X-ray Scattering

Mixed phospholipid micelles (bicelles) are widely applied in nuclear magnetic resonance (NMR) studies of membrane proteins in solution, as they can solubilize these proteins and provide a membrane-like environment. In this work, the structure of bicelles of dihexanoyl phosphatidyl choline (DHPC) and dimyristoyl phosphatidyl choline (DMPC) at different ratios was determined by small-angle X-ray scattering (SAXS) at 37 °C. Samples with concentrations as applied for NMR measurements with 28 wt % lipids were diluted to avoid concentration effects in the SAXS data. The DMPC/DHPC ratio within the bicelles was kept constant by diluting with solutions of finite DHPC concentrations, where the concentration of free DHPC is the same as in the original solution. Absolute-scale modeling of the SAXS data using molecular and concentration constraints reveals a relatively complex set of morphologies of the lipid aggregates as a function of the molar ratio Q of DMPC to DHPC. At Q = 0 (pure DHPC lipids), oblate core-shell micelles are present. At Q = 0.5, the bicelles have a tablet-shaped core-shell cylindrical form with an ellipsoidal cross section. For Q = 1, 2, 3.2, and 4, the bicelles have a rectangular cuboidal structure with a core and a shell, for which the overall length and width increase with Q. At Q = ∞ (pure DMPC), there is coexistence between multilamellar structures and free bilayers. For Q = 1-4, the hydrocarbon core is relatively narrow and the headgroup thickness on the flat areas is larger than that of, respectively, pure DHPC and DMPC, suggesting some mixing of DHPC into these areas and staggering of the molecules. This is further supported by comparisons of the ratio of the areas of rim and flat parts and estimates of the composition of the flat areas.

Nanodisc-to-Nanofiber Transition of Noncovalent Peptide-Phospholipid Assemblies

We report a novel molecular architecture of peptide-phospholipid coassemblies. The amphiphilic peptide Ac-18A-NH₂ (18A), which was designed to mimic apolipoprotein α-helices, has been shown to form nanodisc structures with phospholipid bilayers. We show that an 18A peptide cysteine substitution at residue 11, 18A[A11C], forms fibrous assemblies with 1-palmitoyl-2-oleoyl-phosphatidylcholine at a lipid-to-peptide (L/P) molar ratio of 1, a fiber diameter of 10-20 nm, and a length of more than 1 μm. Furthermore, 18A[A11C] can form nanodiscs with these lipid bilayers at L/P ratios of 4-6. The peptide adopts α-helical structures in both the nanodisc and nanofiber assemblies, although the α-helical bundle structures were evident only in the nanofibers, and the phospholipids of the nanofibers were not lamellar. Fluorescence spectroscopic analysis revealed that the peptide and lipid molecules in the nanofibers exhibited different solvent accessibility and hydrophobicity from those of the nanodiscs. Furthermore, the cysteine substitution at residue 11 did not result in disulfide bond formation, although it was responsible for the nanofiber formation, suggesting that this free sulfhydryl group has an important functional role. Alternatively, the disulfide dimer of 18A[A11C] preferentially formed nanodiscs, even at an L/P ratio of 1. Interconversions of these discoidal and fibrous assemblies were induced by the stepwise addition of free 18A[A11C] or liposomes into the solution. Furthermore, these structural transitions could also be induced by the introduction of oxidative and reductive stresses to the assemblies. Our results demonstrate that heteromolecular lipid-peptide complexes represent a novel approach to the construction of controllable and functional nanoscale assemblies.

How a short pore forming peptide spans the lipid membrane

Many antimicrobial peptides function by forming pores in the plasma membrane of the target cells. Intriguingly, some of these peptides are very short, and thus, it is not known how they can span the membrane, or whether other mechanisms of cell disruption are dominant. Here, the conformation and orientation of the 14-residue peptaibol SPF-5506-A₄ (SPF) are investigated in lipid environments by atomistic and coarse grained molecular dynamics (MD) simulations, circular dichroism, and nuclear magnetic resonance (NMR) experiments. The MD simulations show that SPF is inserted spontaneously in a transmembrane orientation in both 1,2-dimyristoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers resulting in thinning of the bilayers near the peptides, which drives the peptide aggregation. Furthermore, the backbone conformation of the peptide in the bilayer bound state is different from that of the NMR model solved in small bicelles. These results demonstrate that mutual adaption between the peptides and the membrane is likely to be important for pore formation.

Bicelles and Other Membrane Mimics: Comparison of Structure, Properties, and Dynamics from MD Simulations

The increased interest in studying membrane proteins has led to the development of new membrane mimics such as bicelles and nanodiscs. However, only limited knowledge is available of how these membrane mimics are affected by embedded proteins and how well they mimic a lipid bilayer. Herein, we present molecular dynamics simulations to elucidate structural and dynamic properties of small bicelles and compare them to a large alignable bicelle, a small nanodisc, and a lipid bilayer. Properties such as lipid packing and properties related to embedding both an α-helical peptide and a transmembrane protein are investigated. The small bicelles are found to be very dynamic and mainly assume a prolate shape substantiating that small bicelles cannot be regarded as well-defined disclike structures. However, addition of a peptide results in an increased tendency to form disc-shaped bicelles. The small bicelles and the nanodiscs show increased peptide solvation and difference in peptide orientation compared to embedding in a bilayer. The large bicelle imitated a bilayer well with respect to both curvature and peptide solvation, although peripheral binding of short tailed lipids to the embedded proteins is observed, which could hinder ligand binding or multimer formation.

Solid-state NMR methods for oriented membrane proteins

Oriented-sample solid-state NMR represents one of few experimental methods capable of characterising the membrane-bound conformation of proteins in the cell membrane. Since the technique was developed 25 years ago, the technique has been applied to study the structure of helix bundle membrane proteins and antimicrobial peptides, characterise protein-lipid interactions, and derive information on dynamics of the membrane anchoring of membrane proteins. We will review the major developments in various aspects of oriented-sample solid-state NMR, including sample-preparation methods, pulse sequences, theory required to interpret the experiments, perspectives for and guidelines to new experiments, and a number of representative applications.

Centerband-only analysis of rotor-unsynchronized spin echo for measurement of lipid (31) P chemical shift anisotropy

Structural diversity and molecular flexibility of phospholipids are essential for biological membranes to play key roles in numerous cellular processes. Uncovering the behavior of individual lipids in membrane dynamics is crucial for understanding the molecular mechanisms underlying biological functions of cell membranes. In this paper, we introduce a simple method to investigate dynamics of lipid molecules in multi-component systems by measuring the (31) P chemical shift anisotropy (CSA) under magic angle spinning (MAS) conditions. For achieving both signal separation and CSA determination, we utilized a centerband-only analysis of rotor-unsynchronized spin echo (COARSE). This analysis is based on the curve fitting of periodic modulation of centerband intensity along the interpulse delay time in rotor-unsynchronized spin-echo experiments. The utility of COARSE was examined by using phospholipid vesicles, a three-component lipid raft model system, and archaeal purple membranes. We found that the apparent advantages of this method are high resolution and high sensitivity given by the moderate MAS speed and the one-dimensional acquisition with short spin-echo delays. COARSE provides an alternative method for CSA measurement that is effective in the investigation of lipid polymorphologies.