Molecular dynamics of membrane peptides and proteins: principles and comparison to experimental data

Methodologies for the analysis of instantaneous lipid diffusion in MD simulations of large membrane systems

Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically comprise millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library.

Simultaneous determination of interfacial molarities of amide bonds, carboxylate groups, and water by chemical trapping in micelles of amphiphiles containing peptide bond models

Chemical trapping is a powerful approach for obtaining experimental estimates of interfacial molarities of weakly basic nucleophiles in the interfacial regions of amphiphile aggregates. Here, we demonstrate that the chemical probe 4-hexadecyl-2,6-dimethylbenzenediazonium ion (16-ArN(2)(+)) reacts competitively with interfacial water, with the amide carbonyl followed by cleavage of the headgroups from the tail at the amide oxygen, and with the terminal carboxylate groups in micelles of two N-acyl amino-acid amphiphiles, sodium N-lauroylsarcosinate (SLS) and sodium N-lauroylglycinate (SLG), simple peptide bond model amphiphiles. Interfacial molarities (in moles per liter of interfacial volume) of these three groups were obtained from product yields, assuming that selectivity toward a particular nucleophile compared to water is the same in an aqueous reference solution and in the interfacial region. Interfacial carboxylate group molarities are ~1.5 M in both SLS and SLG micelles, but the concentration of the amide carbonyl for SLS micelles is ~4.6-5 times less (ca. 0.7 M) than that of SLG micelles (~3 M). The proton on the secondary N of SLG helps solubilize the amide bond in the aqueous region, but the methyl on the tertiary N of SLS helps solubilize the amide bond in the micellar core, reducing its reaction with 16-ArN(2)(+). Application of chemical trapping to proteins in membrane mimetic interfaces should provide insight into the topology of the protein within the interface because trapping of the amide carbonyl and cleavage at the C-N bond occurs only within the interface, and fragment characterization marks those peptide bonds located within the interface.

The Transmembrane Helix Tilt May Be Determined by the Balance between Precession Entropy and Lipid Perturbation

Hydrophobic helical peptides interact with lipid bilayers in various modes, determined by the match between the length of the helix’s hydrophobic core and the thickness of the hydrocarbon region of the bilayer. For example, long helices may tilt with respect to the membrane normal to bury their hydrophobic cores in the membrane, and the lipid bilayer may stretch to match the helix length. Recent molecular dynamics simulations and potential of mean force calculations have shown that some TM helices whose lengths are equal to, or even shorter than, the bilayer thickness may also tilt. The tilt is driven by a gain in the helix precession entropy, which compensates for the free energy penalty resulting from membrane deformation. Using this free energy balance, we derived theoretically an equation of state, describing the dependence of the tilt on the helix length and membrane thickness. To this end, we conducted coarse-grained Monte Carlo simulations of the interaction of helices of various lengths with lipid bilayers of various thicknesses, reproducing and expanding the previous molecular dynamics simulations. Insight from the simulations facilitated the derivation of the theoretical model. The tilt angles calculated using the theoretical model agree well with our simulations and with previous calculations and measurements.

A solution NMR view of protein dynamics in the biological membrane

Structure determination of membrane-associated proteins (MPs) represents a frontier of structural biology that is characterized by unique challenges in sample preparation and data acquisition. No less important is our ability to study the dynamics of MPs, since MP flexibility and characteristic motions often make sizeable contributions to their function. This review focuses on solution state NMR methods to characterize dynamics of MPs in the membrane environment. NMR approaches to study molecular motions on a wide range of time-scales and their application to membrane proteins are described. Studies of polytopic and bitopic MPs demonstrating the power of such methods to characterize the dynamic behavior of MPs and their interaction with the membrane-mimicking surroundings are presented. Attempts are made to place the dynamic conclusions into a biological context. The importance and limitations of such investigations guarantee that further developments in this field will be actively pursued.