Indolicidin binding induces thinning of a lipid bilayer

Thermodynamics of Micelle Formation and Membrane Fusion Modulate Antimicrobial Lipopeptide Activity

Antimicrobial lipopeptides (AMLPs) are antimicrobial drug candidates that preferentially target microbial membranes. One class of AMLPs, composed of cationic tetrapeptides attached to an acyl chain, have minimal inhibitory concentrations in the micromolar range against a range of bacteria and fungi. Previously, we used coarse-grained molecular dynamics simulations and free energy methods to study the thermodynamics of their interaction with membranes in their monomeric state. Here, we extended the study to the biologically relevant micellar state, using, to our knowledge, a novel reaction coordinate based on hydrophobic contacts. Using umbrella sampling along this reaction coordinate, we identified the critical transition states when micelles insert into membranes. The results indicate that the binding of these AMLP micelles to membranes is thermodynamically favorable, but in contrast to the monomeric case, there are significant free energy barriers. The height of these free energy barriers depends on the membrane composition, suggesting that the AMLPs' ability to selectively target bacterial membranes may be as much kinetic as thermodynamic. This mechanism highlights the importance of considering oligomeric state in solution as criterion when optimizing peptides or lipopeptides as antibiotic leads.

Retinal Conformation Changes Rhodopsin's Dynamic Ensemble

G protein-coupled receptors are vital membrane proteins that allosterically transduce biomolecular signals across the cell membrane. However, the process by which ligand binding induces protein conformation changes is not well understood biophysically. Rhodopsin, the mammalian dim-light receptor, is a unique test case for understanding these processes because of its switch-like activity; the ligand, retinal, is bound throughout the activation cycle, switching from inverse agonist to agonist after absorbing a photon. By contrast, the ligand-free opsin is outside the activation cycle and may behave differently. We find that retinal influences rhodopsin dynamics using an ensemble of all-atom molecular dynamics simulations that in aggregate contain 100 μs of sampling. Active retinal destabilizes the inactive state of the receptor, whereas the active ensemble was more structurally homogenous. By contrast, simulations of an active-like receptor without retinal present were much more heterogeneous than those containing retinal. These results suggest allosteric processes are more complicated than a ligand inducing protein conformational changes or simply capturing a shifted ensemble as outlined in classic models of allostery.

Simulation-Based Approaches for Determining Membrane Permeability of Small Compounds

Predicting the rate of nonfacilitated permeation of solutes across lipid bilayers is important to drug design, toxicology, and signaling. These rates can be estimated using molecular dynamics simulations combined with the inhomogeneous solubility-diffusion model, which requires calculation of the potential of mean force and position-dependent diffusivity of the solute along the transmembrane axis. In this paper, we assess the efficiency and accuracy of several methods for the calculation of the permeability of a model DMPC bilayer to urea, benzoic acid, and codeine. We compare umbrella sampling, replica exchange umbrella sampling, adaptive biasing force, and multiple-walker adaptive biasing force for the calculation of the transmembrane PMF. No definitive advantage for any of these methods in their ability to predict the membrane permeability coefficient Pm was found, provided that a sufficiently long equilibration is performed. For diffusivities, a Bayesian inference method was compared to a generalized Langevin method, both being sensitive to chosen parameters and the slow relaxation of membrane defects. Agreement within 1.5 log units of the computed Pm with experiment is found for all permeants and methods. Remaining discrepancies can likely be attributed to limitations of the force field as well as slowly relaxing collective movements within the lipid environment. Numerical calculations based on model profiles show that Pm can be reliably estimated from only a few data points, leading to recommendations for calculating Pm from simulations.

Sampling errors in free energy simulations of small molecules in lipid bilayers

Free energy simulations are a powerful tool for evaluating the interactions of molecular solutes with lipid bilayers as mimetics of cellular membranes. However, these simulations are frequently hindered by systematic sampling errors. This review highlights recent progress in computing free energy profiles for inserting molecular solutes into lipid bilayers. Particular emphasis is placed on a systematic analysis of the free energy profiles, identifying the sources of sampling errors that reduce computational efficiency, and highlighting methodological advances that may alleviate sampling deficiencies. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Archetypal tryptophan-rich antimicrobial peptides: properties and applications

Drug-resistant microorganisms ('superbugs') present a serious challenge to the success of antimicrobial treatments. Subsequently, there is a crucial need for novel bio-control agents. Many antimicrobial peptides (AMPs) show a broad-spectrum activity against bacteria, fungi or viruses and are strong candidates to complement or substitute current antimicrobial agents. Some AMPs are also effective against protozoa or cancer cells. The tryptophan (Trp)-rich peptides (TRPs) are a subset of AMPs that display potent antimicrobial activity, credited to the unique biochemical properties of tryptophan that allow it to insert into biological membranes. Further, many Trp-rich AMPs cross bacterial membranes without compromising their integrity and act intracellularly, suggesting interactions with nucleic acids and enzymes. In this work, we overview some archetypal TRPs derived from natural sources, i.e., indolicidin, tritrpticin and lactoferricin, summarising their biochemical properties, structures, antimicrobial activities, mechanistic studies and potential applications.

Free Energy Calculations for the Peripheral Binding of Proteins/Peptides to an Anionic Membrane. 1. Implicit Membrane Models

The binding of peptides and proteins to the surface of complex lipid membranes is important in many biological processes such as cell signaling and membrane remodeling. Computational studies can aid experiments by identifying physical interactions and structural motifs that determine the binding affinity and specificity. However, previous studies focused on either qualitative behaviors of protein/membrane interactions or the binding affinity of small peptides. Motivated by this observation, we set out to develop computational protocols for bimolecular binding to charged membranes that are applicable to both peptides and large proteins. In this work, we explore a method based on an implicit membrane/solvent model (generalized Born with a simple switching in combination with the Gouy-Chapman-Stern model for a charged interface), which we expect to lead to useful results when the binding does not implicate significant membrane deformation and local demixing of lipids. We show that the binding free energy can be efficiently computed following a thermodynamic cycle similar to protein-ligand binding calculations, especially when a Bennett acceptance ratio based protocol is used to consider both the membrane bound and solution conformational ensembles. Test calculations on a series of peptides show that our computational approach leads to binding affinities in encouraging agreement with experimental data, including for the challenging example of the bringing of flexible MARCKS-ED peptides to membranes. The calculations highlight that for a membrane with a significant fraction of anionic lipids, it is essential to include the effect of ion adsorption using the Stern model, which significantly modifies the effective surface charge. This implicit membrane model based computational protocol helps lay the groundwork for more systematic analysis of protein/peptide binding to membranes of complex shape and composition.

Penetration of HIV-1 Tat47-57 into PC/PE Bilayers Assessed by MD Simulation and X-ray Scattering

The interactions of the basic, cell-penetrating region (Y47GRKKRRQRRR57) of the HIV-1 Tat protein with dioleoylphosphatidylcholine (DOPC) bilayers were previously assessed by comparing experimental X-ray diffuse scattering with atomistic molecular dynamics simulations. Here, we extend this investigation by evaluating the influence of phosphatidylethanolamine (PE) lipids. Using experimental bilayer form factors derivedfrom X-ray diffuse scattering data as a guide, our simulations indicate that Tat peptides localize close to the carbonyl-glycerol group in the headgroup region of bilayers composed of either DOPC or DOPC:DOPE (1:1) lipid. Our results also suggest that Tat peptides may more frequently insert into the hydrophobic core of bilayers composed of PC:PE (1:1) lipids than into bilayers composed entirely of PC lipids. PE lipids may facilitate peptide translocation across a lipid bilayer by stabilizing intermediate states in which hydrated peptides span the bilayer.

Thermodynamics of antimicrobial lipopeptide binding to membranes: origins of affinity and selectivity

The development of novel antibiotic drugs is one of the most pressing biomedical problems due to the increasing number of antibiotic-resistant pathogens. Antimicrobial peptides and lipopeptides are a promising category of candidates, but the molecular origins of their antimembrane activity is unclear. Here we explore a series of recently developed antimicrobial lipopeptides, using coarse-grained molecular-dynamics simulations and free energy methods to uncover the thermodynamics governing their binding to membranes. Specifically, we quantify C16-KGGK's binding affinity to the two types of membrane by umbrella sampling. We also examined the origin of C16-KGGK's selectivity for bacterial versus mammalian membranes by systematically varying the peptide sequence and salt concentration. Our data showed that the C16 hydrophobic tail is the main contributor to its affinity to lipid membrane, whereas the peptide portion is mainly responsible for its selectivity. Furthermore, the electrostatic interaction between the cationic peptide and anionic bacterial membrane plays a significant role in the selectivity.

Action of the multifunctional peptide BP100 on native biomembranes examined by solid-state NMR

Membrane composition is a key factor that regulates the destructive activity of antimicrobial peptides and the non-leaky permeation of cell penetrating peptides in vivo. Hence, the choice of model membrane is a crucial aspect in NMR studies and should reflect the biological situation as closely as possible. Here, we explore the structure and dynamics of the short multifunctional peptide BP100 using a multinuclear solid-state NMR approach. The membrane alignment and mobility of this 11 amino acid peptide was studied in various synthetic lipid bilayers with different net charge, fluidity, and thickness, as well as in native biomembranes harvested from prokaryotic and eukaryotic cells. (19)F-NMR provided the high sensitivity and lack of natural abundance background that are necessary to observe a labelled peptide even in protoplast membranes from Micrococcus luteus and in erythrocyte ghosts. Six selectively (19)F-labeled BP100 analogues gave remarkably similar spectra in all of the macroscopically oriented membrane systems, which were studied under quasi-native conditions of ambient temperature and full hydration. This similarity suggests that BP100 has the same surface-bound helical structure and high mobility in the different biomembranes and model membranes alike, independent of charge, thickness or cholesterol content of the system. (31)P-NMR spectra of the phospholipid components did not indicate any bilayer perturbation, so the formation of toroidal wormholes or micellarization can be excluded as a mechanism of its antimicrobial or cell penetrating action. However, (2)H-NMR analysis of the acyl chain order parameter profiles showed that BP100 leads to considerable membrane thinning and thereby local destabilization.

Membrane environment modulates the pKa values of transmembrane helices

In this work, we apply the recently developed constant pH molecular dynamics technique to study protonation equilibria of titratable side chains in the context of simple transmembrane (TM) helices and explore the effect of pH on their configurations in membrane bilayers. We observe that, despite a significant shift toward neutral states, considerable population of different side chains stay in the charged state that give rise to pK_a values around 9.6 for Asp and Glu and 4.5 to 6 for His and Lys side chains, respectively. These charged states are highly stabilized by favorable interactions between head groups, water molecules, and the charged side chains that are facilitated by substantial changes in the configuration of the peptides. The pH dependent configurations and the measured pK_a values are in good agreement with relatively recent solid state NMR measurements. Our results presented here demonstrate that all-atom constant pH molecular dynamics can be applied to membrane proteins and peptides to obtain reliable pK_a values and pH dependent behavior for these systems.

Dynamics of lipids, cholesterol, and transmembrane α-helices from microsecond molecular dynamics simulations

Extensive all-atom molecular dynamics (∼24 μs total) allowed exploration of configurational space and calculation of lateral diffusion coefficients of the components of a protein-embedded, cholesterol-containing model bilayer. The three model membranes are composed of an ∼50/50 (by mole) dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayer and contained an α-helical transmembrane protein (HIV-1 gp41 TM). Despite the high concentration of cholesterol, normal Brownian motion was observed and the calculated diffusion coefficients (on the order of 10(-9) cm(2)/s) are consistent with experiments. Diffusion is sensitive to a variety of parameters, and a temperature difference of ∼4 K from thermostat artifacts resulted in 2-10-fold differences in diffusion coefficients and significant differences in lipid order, membrane thickness, and unit cell area. Also, the specific peptide sequence likely underlies the consistently observed faster diffusion in one leaflet. Although the simulations here present molecular dynamics (MD) an order of magnitude longer than those from previous studies, the three systems did not approach ergodicity. The distributions of cholesterol and DPPC around the peptides changed on the microsecond time scale, but not significantly enough to thoroughly explore configurational space. These simulations support conclusions of other recent microsecond MD in that even longer time scales are needed for equilibration of model membranes and simulations of more realistic cellular or viral bilayers.

Lipid clustering correlates with membrane curvature as revealed by molecular simulations of complex lipid bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.

Cell membranes play important roles in vivo both in shielding the cell interior from the surrounding environment and in cell function through lipid components of the membrane having roles in controlling protein function, cell signaling etc. We employ molecular dynamics simulations to explore the behavior of biologically realistic membrane models. Our simulations reveal nano-domain clustering of the glycolipid GM3 and to a lesser extent of the anionic lipid phosphatidylinositol 4,5-bisphophate (PIP₂). When including transmembrane proteins we are able to observe preferential interactions of known regulatory lipids (e.g. GM3, PIP₂ and cholesterol) with the proteins. Membrane curvature is shown to be coupled to the local lipid composition, suggestive of a link between lipid nano-domains and membrane geometry.

Towards fast, rigorous and efficient conformational sampling of biomolecules: Advances in accelerated molecular dynamics

BACKGROUND

Accelerated molecular dynamics (aMD) has been proven to be a powerful biasing method for enhanced sampling of biomolecular conformations on general-purpose computational platforms. Biologically important long timescale events that are beyond the reach of standard molecular dynamics can be accessed without losing the detailed atomistic description of the system in aMD. Over other biasing methods, aMD offers the advantages of tuning the level of acceleration to access the desired timescale without any advance knowledge of the reaction coordinate.

SCOPE OF REVIEW

Recent advances in the implementation of aMD and its applications to small peptides and biological macromolecules are reviewed here along with a brief account of all the aMD variants introduced in the last decade.

MAJOR CONCLUSIONS

In comparison to the original implementation of aMD, the recent variant in which all the rotatable dihedral angles are accelerated (RaMD) exhibits faster convergence rates and significant improvement in statistical accuracy of retrieved thermodynamic properties. RaMD in conjunction with accelerating diffusive degrees of freedom, i.e. dual boosting, has been rigorously tested for the most difficult conformational sampling problem, protein folding. It has been shown that RaMD with dual boosting is capable of efficiently sampling multiple folding and unfolding events in small fast folding proteins.

GENERAL SIGNIFICANCE

RaMD with the dual boost approach opens exciting possibilities for sampling multiple timescales in biomolecules. While equilibrium properties can be recovered satisfactorily from aMD-based methods, directly obtaining dynamics and kinetic rates for larger systems presents a future challenge. This article is part of a Special Issue entitled Recent developments of molecular dynamics.