Probing columnar discotic liquid crystals by EPR spectroscopy with a rigid-core nitroxide spin probe

Predicting and interpreting EPR spectra of POPC lipid bilayers with transmembrane α-helical peptides from all-atom molecular dynamics simulations

This study reports a large-scale all-atom MD simulation of POPC lipid bilayers in the presence of different concentrations of the transmembrane peptide acetyl-K2(LA)12K2-amide ((LA)12) and doped with 5-PC paramagnetic spin probes used in EPR studies. We apply a combined MD-EPR simulation methodology for the prediction of EPR spectra directly and completely from MD trajectories. This approach serves three major purposes. Firstly, comparing predicted EPR spectra with experimental ones, which are highly sensitive to motions, provides an ultimate test bed for the force fields currently employed for modeling lipid bilayer systems with embedded proteins or peptides. Secondly, simulations of EPR spectra directly from the atomistic MD models simplify the interpretation of the EPR line shapes and their changes induced by the presence of peptides in the lipid bilayer. These changes are directly linked to the dynamics and order of spin probes and POPC host molecules. Lastly and importantly, we demonstrate how the MD-EPR methodology can be employed to test the validity and limitations of the widely used approach for the estimation of the order parameter of lipids directly from the EPR experimental line shapes.

Direct Prediction of EPR Spectra from Lipid Bilayers: Understanding Structure and Dynamics in Biological Membranes

Of the many biophysical techniques now being brought to bear on studies of membranes, electron paramagnetic resonance (EPR) of nitroxide spin probes was the first to provide information about both mobility and ordering in lipid membranes. Here, we report the first prediction of variable temperature EPR spectra of model lipid bilayers in the presence and absence of cholesterol from the results of large scale fully atomistic molecular dynamics (MD) simulations. Three types of structurally different spin probes were employed in order to study different parts of the bilayer. Our results demonstrate very good agreement with experiment and thus confirm the accuracy of the latest lipid force fields. The atomic resolution of the simulations allows the interpretation of the molecular motions and interactions in terms of their impact on the sensitive EPR line shapes. Direct versus indirect effects of cholesterol on the dynamics of spin probes are analysed. Given the complexity of structural organisation in lipid bilayers, the advantage of using a combined MD-EPR simulation approach is two-fold. Firstly, prediction of EPR line shapes directly from MD trajectories of actual phospholipid structures allows unambiguous interpretation of EPR spectra of biological membranes in terms of complex motions. Secondly, such an approach provides an ultimate test bed for the up-to-date MD simulation models employed in the studies of biological membranes, an area that currently attracts great attention.

Prediction of EPR Spectra of Lyotropic Liquid Crystals using a Combination of Molecular Dynamics Simulations and the Model-Free Approach

We report the first application of fully atomistic molecular dynamics (MD) simulations to the prediction of the motional electron paramagnetic resonance (EPR) spectra of lyotropic liquid crystals in different aggregation states doped with a paramagnetic spin probe. The purpose of this study is twofold. First, given that EPR spectra are highly sensitive to the motions and order of the spin probes doped within lyotropic aggregates, simulation of EPR line shapes from the results of MD modelling provides an ultimate test bed for the force fields currently employed to model such systems. Second, the EPR line shapes are simulated using the motional parameters extracted from MD trajectories using the Model-Free (MF) approach. Thus a combined MD-EPR methodology allowed us to test directly the validity of the application of the MF approach to systems with multi-component molecular motions. All-atom MD simulations using the General AMBER Force Field (GAFF) have been performed on sodium dodecyl sulfate (SDS) and dodecyltrimethylammonium chloride (DTAC) liquid crystals. The resulting MD trajectories were used to predict and interpret the EPR spectra of pre-micellar, micellar, rod and lamellar aggregates. The predicted EPR spectra demonstrate good agreement with most of experimental line shapes thus confirming the validity of both the force fields employed and the MF approach for the studied systems. At the same time simulation results confirm that GAFF tends to overestimate the packing and the order of the carbonyl chains of the surfactant molecules.