Chapter 4 Molecular Dynamics Simulation of Lipid–Protein Interactions

Modeling the molecular fingerprint of protein-lipid interactions of MLKL on complex bilayers

Lipids, the structural part of membranes, play important roles in biological functions. However, our understanding of their implication in key cellular processes such as cell division and protein-lipid interaction is just emerging. This is the case for molecular interactions in mechanisms of cell death, where the role of lipids for protein localization and subsequent membrane permeabilization is key. For example, during the last stage of necroptosis, the mixed lineage kinase domain-like (MLKL) protein translocates and, eventually, permeabilizes the plasma membrane (PM). This process results in the leakage of cellular content, inducing an inflammatory response in the microenvironment that is conducive to oncogenesis and metastasis, among other pathologies that exhibit inflammatory activity. This work presents insights from long all-atom molecular dynamics (MD) simulations of complex membrane models for the PM of mammalian cells with an MLKL protein monomer. Our results show that the binding of the protein is initially driven by the electrostatic interactions of positively charged residues. The protein bound conformation modulates lipid recruitment to the binding site, which changes the local lipid environment recruiting PIP lipids and cholesterol, generating a unique fingerprint. These results increase our knowledge of protein-lipid interactions at the membrane interface in the context of molecular mechanisms of the necroptotic pathway, currently under investigation as a potential treatment target in cancer and inflamatory diseases.

Lipid-mediated prestin organization in outer hair cell membranes and its implications in sound amplification

Prestin is a high-density motor protein in the outer hair cells (OHCs), whose conformational response to acoustic signals alters the shape of the cell, thereby playing a major role in sound amplification by the cochlea. Despite recent structures, prestin's intimate interactions with the membrane, which are central to its function remained unresolved. Here, employing a large set (collectively, more than 0.5 ms) of coarse-grained molecular dynamics simulations, we demonstrate the impact of prestin's lipid-protein interactions on its organization at densities relevant to the OHCs and its effectiveness in reshaping OHCs. Prestin causes anisotropic membrane deformation, which mediates a preferential membrane organization of prestin where deformation patterns by neighboring copies are aligned constructively. The resulting reduced membrane rigidity is hypothesized to maximize the impact of prestin on OHC reshaping. These results demonstrate a clear case of protein-protein cooperative communication in membrane, purely mediated by interactions with lipids.

Lipid Configurations from Molecular Dynamics Simulations

The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na+ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the 31P chemical shift anisotropy and 14N and CD3 NMR quadrupole splittings.

Molecular Dynamics Simulations of the Host Defense Peptide Temporin L and Its Q3K Derivative: An Atomic Level View from Aggregation in Water to Bilayer Perturbation

Temporin L (TempL) is a 13 residue Host Defense Peptide (HDP) isolated from the skin of frogs. It has a strong affinity for lipopolysaccharides (LPS), which is related to its high activity against Gram-negative bacteria and also to its strong tendency to neutralize the pro-inflammatory response caused by LPS release from inactivated bacteria. A designed analog with the Q3K substitution shows an enhancement in both these activities. In the present paper, Molecular Dynamics (MD) simulations have been used to investigate the origin of these improved properties. To this end, we have studied the behavior of the peptides both in water solution and in the presence of LPS lipid-A bilayers, demonstrating that the main effect through which the Q3K substitution improves the peptide activities is the destabilization of peptide aggregates in water.

Heterogeneous dielectric generalized Born model with a van der Waals term provides improved association energetics of membrane-embedded transmembrane helices

The heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is extended by the addition of a van der Waals dispersion term to better describe the nonpolar components of the free energy of solvation. The new model, termed HDGBvdW, improves the energy estimates in the hydrophobic interior of the membrane, where polar and charged species are rarely found and nonpolar interactions become significant. The implicit van der Waals term for the membrane environment extends the model from Gallicchio et al. (J. Comput. Chem. 2004, 25, 479) by combining separate contributions from each of the membrane components. The HDGBvdW model is validated with a series of test cases ranging from membrane insertion and pair association profiles of amino acid side chain analogs and transmembrane helices. Overall, the HDGBvdW model leads to increased agreement with explicit membrane simulation results and experimental data. © 2016 Wiley Periodicals, Inc.

Fluorescence spectroscopy and molecular dynamics simulations in studies on the mechanism of membrane destabilization by antimicrobial peptides

Since their initial discovery, 30 years ago, antimicrobial peptides (AMPs) have been intensely investigated as a possible solution to the increasing problem of drug-resistant bacteria. The interaction of antimicrobial peptides with the cellular membrane of bacteria is the key step of their mechanism of action. Fluorescence spectroscopy can provide several structural details on peptide-membrane systems, such as partition free energy, aggregation state, peptide position and orientation in the bilayer, and the effects of the peptides on the membrane order. However, these "low-resolution" structural data are hardly sufficient to define the structural requirements for the pore formation process. Molecular dynamics simulations, on the other hand, provide atomic-level information on the structure and dynamics of the peptide-membrane system, but they need to be validated experimentally. In this review we summarize the information that can be obtained by both approaches, highlighting their versatility and complementarity, suggesting that their synergistic application could lead to a new level of insight into the mechanism of membrane destabilization by AMPs.

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

Molecular dynamics (MD) simulation is a standard tool used to assess the motion of biomolecules at atomic resolution. It requires a so-called "force field" that allows the evaluation of an empirical energy from the 3D coordinates of the atoms in the system. In this chapter, the application of MD simulations to membrane proteins and peptides is described with a particular emphasis on the comparison of MD results to experimental data. Such a comparison can be used either for (1) validating the results of a simulation, (2) interpreting an experiment at the atomic level, or (3) calibrating the force field. This last step is particularly important for the use of MD as a predictive tool. As an illustration, a comparison of (2)H NMR experiments to MD simulations of a transmembrane peptide is presented and discussed.

Congruency between biophysical data from multiple platforms and molecular dynamics simulation of the double-super helix model of nascent high-density lipoprotein

The predicted structure and molecular trajectories from >80 ns molecular dynamics simulation of the solvated Double-Super Helix (DSH) model of nascent high-density lipoprotein (HDL) were determined and compared with experimental data on reconstituted nascent HDL obtained from multiple biophysical platforms, including small angle neutron scattering (SANS) with contrast variation, hydrogen-deuterium exchange tandem mass spectrometry (H/D-MS/MS), nuclear magnetic resonance spectroscopy (NMR), cross-linking tandem mass spectrometry (MS/MS), fluorescence resonance energy transfer (FRET), electron spin resonance spectroscopy (ESR), and electron microscopy. In general, biophysical constraints experimentally derived from the multiple platforms agree with the same quantities evaluated using the simulation trajectory. Notably, key structural features postulated for the recent DSH model of nascent HDL are retained during the simulation, including (1) the superhelical conformation of the antiparallel apolipoprotein A1 (apoA1) chains, (2) the lipid micellar-pseudolamellar organization, and (3) the solvent-exposed Solar Flare loops, proposed sites of interaction with LCAT (lecithin cholesteryl acyltransferase). Analysis of salt bridge persistence during simulation provides insights into structural features of apoA1 that forms the backbone of the lipoprotein. The combination of molecular dynamics simulation and experimental data from a broad range of biophysical platforms serves as a powerful approach to studying large macromolecular assemblies such as lipoproteins. This application to nascent HDL validates the DSH model proposed earlier and suggests new structural details of nascent HDL.