Molecular Dynamics and EPR spectroscopy on lipid bilayers: new approaches to study biological membranes

Lead Research Organisation: University of East Anglia
Department Name: Chemistry

Abstract

Lipid bilayers are the main building blocks of biological membranes. They play a key part in many important biological mechanisms in membranes such as providing living cells with energy, organising and regulating enzyme activities, facilitating the transduction of information and even the supply of substrates for biosynthesis and for signalling molecules. It is widely accepted that membranes do not form homogeneous fluid lipid phases but, in contrast, lipids are organised into phase separated dynamical domains depending on various conditions.
Knowledge of molecular interactions, thermodynamics and system composition effects are crucial for understanding the role which different lipids play in vital life processes in biological membranes. This knowledge is also important for the design of drug delivery systems based on liposomes (artificial vesicles composed of lipid bilayers). An example would be "trigger release liposomes" where temperature sensitive liposomes could be engineered in a way to have phase separated domains to release their content upon trigger.

Of the 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 mobility and ordering in lipid membranes and lipid bilayer systems. Spin probes, specially designed chemical agents that carry a stable unpaired electron, can be introduced within complex partially ordered molecular systems in order to report on the order and dynamics of surrounding molecules. They can probe different depths / parts of the bilayer and also be attached to embedded peptides and proteins. Because an electron has a magnetic moment it can interact with an external magnetic field. EPR measures this interaction in the form of spectral line shape. The orientation of the spin label to the magnetic field has a dramatic effect on this line shape and therefore molecular mobility, dynamics and distribution can be studied. EPR is a technique that acts as a snapshot of very fast molecular motions and can resolve molecular re-orientational dynamics of the introduced spin probe over times shorter than a billionth of a second. However, the analysis of the rich and complex in information EPR lineshape requires full computer simulation. Current approaches rely on simplified parametrised models of motion and require fitting of EPR spectra with multiple adjustable parameters. Such approaches in many cases do not provide an unambiguous interpretation of the spectra preventing definite conclusions about motion and order in multi-component lipid bilayers to be reached.

The last decade has seen radical improvement in the molecular modelling of complex molecular and bio-molecular systems including lipid bilayers using Molecular Dynamics (MD) simulation techniques. MD simulations are now much faster and more accurate allowing researchers to predict complex molecular phenomena using actual structures.
This project will bring together MD and EPR and will attempt for the first time simulation of EPR spectra of biological membranes directly from the results of MD. The advantage of such an approach is twofold. Firstly, it will provide the improvement and will facilitate the interpretation of EPR of biological membranes. Secondly, our MD-EPR methodology will serve as a test bed for advanced computational models for lipid bilayers simulations.

We will use the unique combination of expertises from UEA in both EPR and atomistic MD simulations of spin labelled bio-molecules and coarse-grained simulations on large scale systems provided by Durham.
We will use a novel MD-EPR methodology to address the key problems of understanding molecular interactions, thermodynamics and system composition effects on the formation and dynamics of lipid domains, the organisation and dynamics of lipids around trans-membrane proteins, and the role of cholesterol as a lipid bilayer stabiliser.

Planned Impact

Lipid bilayers represent one of the most important classes of supramolecular assemblies for life forms on earth. Our novel MD-EPR analysis tool will, for the first time, bring together state-of-the-art MD modelling of lipid bilayer systems with EPR in order to investigate dynamics and molecular organisation in lipid membranes. The knowledge generated by this project will contribute to our understanding of molecular interaction mechanisms which are responsible for the dynamics and self-assembly of phospholipids at the lipid bilayers and in protein-lipid complexes. It will also provide critical insights into the role cholesterol plays in these processes. Since lipid bilayers play significant roles in many important bio-molecular processes (in membrane transports, signal transduction, apoptosis) such knowledge is vital for advances in biological and medical research in general. Thus over a longer period this will contribute to the enhancement of the quality of life and health.
Understanding intermolecular interactions in the lipid membranes is also vital for designing novel lipid based materials with desired properties for pharmaceutical applications. Thus in a short term pharmaceutical companies working on the development of lipid-based drug delivery systems such as "trigger release liposomes" would benefit from the knowledge generated by our research.
The outcomes of the project would also be of relevance to industrial companies working on the rational design of hybrid nanomaterials containing lipids for biomedical and technological applications (lipid bilayer coated nanomaterials such as noble metal and silica nanoparticles, semiconductor quantum dots and carbon nanotubes).
Another outcome of the project will be novel advanced software for the prediction of EPR spectra from MD simulations of lipid bilayers with spin probes. Our simulation program will be able to predict directly from MD simulations the effects of hydrophobicity and oxygen permittivity on the EPR data. These new functionalities will be integrated within the general MD-EPR simulation suite and will be available from a dedicated website. This will contribute to the dissemination of our results and will increase the international impact of our research.
Since the project will contribute to fundamental understanding of structure and self-assembly of biological membranes, which are essential elements of all living cells, it is important to EPSRC's key theme Healthcare Technologies. It is also directly related to research area "Biophysics and Soft matter Physics" as well as has relevance to the 'Physical Sciences Grand Challenges' strategic priority, in particular 'Understanding the physics of life'.
Finally, EPSRC has recently made a large strategic investment into the development of new functional materials with lipid bilayer membranes by the UK groups based in London, Nottingham, Leeds, Durham and Cambridge (Programme grant EP/J017566/1). Our MD and EPR studies will complement the application of other physical techniques to study lipid bilayers and the work of research groups in the UK.

Publications

10 25 50
publication icon
Catte A (2016) Molecular electrometer and binding of cations to phospholipid bilayers. in Physical chemistry chemical physics : PCCP

publication icon
Catte A (2018) Direct Prediction of EPR Spectra from Lipid Bilayers: Understanding Structure and Dynamics in Biological Membranes. in Chemphyschem : a European journal of chemical physics and physical chemistry

publication icon
Jayasooriya UA (2016) Rate of Molecular Transfer of Allyl Alcohol across an AOT Surfactant Layer Using Muon Spin Spectroscopy. in Langmuir : the ACS journal of surfaces and colloids

 
Description It is widely recognised that lipid bilayers and their various phases play a significant role in many vital biological processes in cell membranes and that understanding their behaviour at a molecular level is crucial. The main objective of this award was to bring together theory and experiment in order to shed light on the molecular organisation and dynamics in domain forming mixtures of lipid bilayers at different compositions and conditions and in the presence of transmembrane peptides. This was achieved by combining state-of-the-art Molecular Dynamics (MD) simulations, Electron Paramagnetic Resonance (EPR) spectroscopy and a novel methodology for prediction of EPR spectra from MD trajectories.

We have reported 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 MD simulations. Three types of structurally different nitroxide spin probes were employed in order to study different parts of the bilayer. EPR spectra were predicted directly from long MD trajectories using our MD-EPR simulation methodology. The results demonstrate very good agreement with the measured spectra confirming the accuracy of the latest lipid force fields (Slipids). 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. The advantage of our direct MD-EPR simulation approach over previous ones is that it explicitly accounts for the complexity of both local and global motions of both the probes and the host phospholipids. We demonstrated direct and indirect effects of cholesterol on the motions and order of spin probes and the role that hydrogen bonds play in such interactions. Analysis indicated strong linear regression between re-orientational motions of the nitroxide group in spin probe molecules and associated parts in phospholipid and cholesterol molecules upon varying the temperature. In addition, we have carried out the study of the effects on the structure and mobility of both lipids and spin probes arising from the presence of cardiolipin in lipid membranes.

We have performed the first prediction of the EPR spectra from MD simulations of lipid bilayers and their binary and ternary mixtures in the presence of transmembrane peptides. Two cases were considered, namely, spin labelled transmembrane peptides in lipid bilayers and pure peptides in lipid bilayers doped with spin probes. Our results showed very good agreement with experiment confirming the accuracy of the currently employed MD models. Analysis demonstrated that in the presence of a peptide the acyl chains of the nearby phospholipid molecules become much disordered and disarrange significantly. This has strong impact on the lateral and rotational diffusion of both the lipids and the probes.

Given the complexity of structural organisation in lipid bilayers the advantage of using the 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 multi-component motions, avoiding the uncertainty arising from the fitting of spectra using multiple adjustable parameters employed in previous studies. Secondly, such a synergistic MD-EPR 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.
We have extended the application of the developed MD-EPR methodology to other bio-molecular systems including its first application to spin labelled duplex and single-strand DNA fragments.

In addition, we have applied molecular modelling to study the processes of peptide aggregation and pore formation in lipid bilayers and vesicles by the highly cationic antimicrobial peptide (AMP), Chrysophsin-3, using coarse-grained MD (CG-MD) simulations and potential of mean force calculations. AMPs are small cationic proteins that are able to destabilise a lipid bilayer structure through one or more modes of action. We showed that aggregation of peptides at the surface, leads to pronounced deformation of lipid bilayers, leading in turn to lipid protrusions for peptide : ligand ratios > 1:12. In addition, aggregation of Chrysophsin-3 peptides within the centre of a lipid bilayer leads to spontaneous formation of pores and aggregates, consistent with previously reported experimental data.

The project also contributed to the development of the advanced software package for the prediction and analysis of EPR spectroscopic data from the results of MD simulations.
Exploitation Route The developed MD-EPR methodology would be particularly valuable to both experimental researchers and molecular modelers who study complex membrane phenomena such as lipid domain aggregation in ternary lipid systems, including miscibility critical points on the phase diagrams and the formation of lipid rafts.
Sectors Healthcare,Other

 
Description The work has contributed to the completion of the first version of SpinMolDyn, an advanced software suite for prediction and analysis of EPR spectroscopic data from the results of Molecular dynamics (MD) simulations. SpinMolDyn is a multi-modular software that offers fast simulation times using several original advanced theoretical and computational techniques developed in the PI's group. As a result simulations of spectral line shapes literally take seconds or minutes on a common laptop which makes the software superior in terms of speed and accuracy to any other code with similar capability available to date. SpinMolDyn can be integrated with the commercially distributed software products developed by leading vendors and sold to both academic researchers and industries, including pharmaceutical companies that combine spectroscopic measurements and molecular modelling in the drugs discovery process. Such impact opportunities are currently being explored through the licensing agreement mechanisms. The novel MD-EPR simulation methodology with application to lipid membranes and other lipotropic systems has attracted strong interest among research groups working in collaboration with various industrial and pharmaceutical partners on the development of drug delivery systems. Our methodology is particularly attractive to them as it is based on the combination of highly sensitive experimental technique with atomistic MD modelling and is particularly designed to perform studies at elevated temperatures to meet physiological conditions. It is currently being tried by the pharmaceutical research groups at UEA who work on polymer and lipid-based liquid crystalline platforms for drug delivery. A review published by the PI in Liquid Crystals, 2018, 45, 2139-2157 covers the methodological developments and applications to lyotropic systems carried out during the project. It has attracted strong international interest and serves as a teaching tool for different researches and postgraduate students. According to the website of Liquid Crystals the review is among their top ten most read publications to date.
Sector Digital/Communication/Information Technologies (including Software),Education,Healthcare,Pharmaceuticals and Medical Biotechnology
Impact Types Economic

 
Description EPSRC DTP Studentship
Amount £68,000 (GBP)
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 10/2015 
End 02/2019
 
Description Research Grant
Amount £267,040 (GBP)
Funding ID EP/P007554/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 04/2017 
End 04/2019
 
Description NMR lipids blogpost project 
Organisation NMRlipids project
Sector Public 
PI Contribution PI and one of the PDRAs are the members of 'Matching lipid force fields with NMR data' (http://nmrlipids.blogspot.com), an international network that brings together modellers to find a lipid force field that matches measured NMR parameters, with the aim of improving molecular models for biological membranes. This in an open scientific collaboration project to understand the atomistic resolution structures of lipid bilayers which unites the efforts of research groups from around the world. We contribute the data from all-atom large scale MD simulations on lipid bilayers for the testing of quality of various force fields against experimental NMR spectral line shapes and other parameters.
Collaborator Contribution Generation of the data from MD simulations on lipid bilayers for the testing of quality of various force fields against high quality NMR data. The aim is to develop an all-atom force field that will match these data.
Impact As a result of this network international collaboration a research paper was published in 2016. Catte, et al., 'Molecular electrometer and binding of cations to phospholipid bilayers', Phys. Chem. Chem. Phys., 18, 32560-32569, (2016)
Start Year 2014
 
Title Advanced software for prediction and analysis of EPR spectroscopic data from the results of Molecular dynamics simulations 
Description The software is developed in two versions, one based on MATLAB code for just spectra generation and a more advanced one based on C programming combined with the Qt GUI application development framework for building a state-of-the-art modern GUI for cross-platform applications. It provides simulation of CW EPR spectra directly and completely form single dynamical trajectories of spin probe's re-orientational diffusional motions. Applicable to both Brownian Dynamics (BD) and Molecular Dynamics (MD) trajectories. 
Type Of Technology Software 
Year Produced 2011 
Impact So far the software has been used by the group members at UEA for the simulation and analysis of the EPR spectra of various systems with doped spin probes leading to the publication in high impact journals. 
 
Description Committee member of British Liquid Crystal Society (BLCS). 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Other academic audiences (collaborators, peers etc.)
Results and Impact From 2012 is serving as Committee member of British Liquid Crystal Society.

As a Committee member I have been invited to write a Conference report paper about the Liquid Crystals meeting in Durham in 2014
Year(s) Of Engagement Activity 2012,2014