Molecular dynamics simulations of T-20 HIV fusion inhibitor interacting with model membranes

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Molecular dynamics simulations of T-2410 and T-2429 HIV fusion inhibitors interacting with model membranes: Insight into peptide behavior, structure and dynamics

T-2410 and T-2429 are HIV fusion inhibitor peptides (FI) designed to present a higher efficiency even against HIV strains that developed resistance against other FIs. Similar peptides were shown to interact with model membranes both in the liquid disordered phase and in the liquid ordered state. Those results indicated that such interaction is important to function and could be correlated with their effectiveness. Extensive molecular dynamics simulations were carried out to investigate the interactions between both T-2410 and T-2429 with bilayers of pure 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and a mixture of POPC/cholesterol (Chol) (1:1). It was observed that both peptides interact strongly with both membrane systems, especially with the POPC/Chol systems, where these peptides show the highest number of H-bonds observed so far. T-2410 and T-2429 showed higher extent of interaction with bilayers when compared to T-20 or T-1249 in previous studies. This is most notable in POPC/Chol membranes where, although able to form H-bonds with Chol, they do so to a lesser extent than T-1249 does, the latter being the only FI peptide so far that was observed to form H-bonds with Chol. This behavior suggests that interaction of FI peptides with rigid Chol rich membranes may not be as dependent from peptide/Chol H-bond formation as previous results of T-1249 behavior led to believe. As in other similar peptides, the higher ability to interact with membranes shown by T-2410 and T2429 is probably correlated with its higher inhibitory efficiency.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Computer-Aided Approaches for Targeting HIVgp41

Virus-cell fusion is the primary means by which the human immunodeficiency virus-1 (HIV) delivers its genetic material into the human T-cell host. Fusion is mediated in large part by the viral glycoprotein 41 (gp41) which advances through four distinct conformational states: (i) native, (ii) pre-hairpin intermediate, (iii) fusion active (fusogenic), and (iv) post-fusion. The pre-hairpin intermediate is a particularly attractive step for therapeutic intervention given that gp41 N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR) domains are transiently exposed prior to the formation of a six-helix bundle required for fusion. Most peptide-based inhibitors, including the FDA-approved drug T20, target the intermediate and there are significant efforts to develop small molecule alternatives. Here, we review current approaches to studying interactions of inhibitors with gp41 with an emphasis on atomic-level computer modeling methods including molecular dynamics, free energy analysis, and docking. Atomistic modeling yields a unique level of structural and energetic detail, complementary to experimental approaches, which will be important for the design of improved next generation anti-HIV drugs.

Computational analysis of local membrane properties

In the field of biomolecular simulations, dynamics of phospholipid membranes is of special interest. A number of proteins, including channels, transporters, receptors and short peptides are embedded in lipid bilayers and tightly interact with phospholipids. While the experimental measurements report on the spatial and/or temporal average membrane properties, simulation results are not restricted to the average properties. In the current study, we present a collection of methods for an efficient local membrane property calculation, comprising bilayer thickness, area per lipid, deuterium order parameters, Gaussian and mean curvature. The local membrane property calculation allows for a direct mapping of the membrane features, which subsequently can be used for further analysis and visualization of the processes of interest. The main features of the described methods are highlighted in a number of membrane systems, namely: a pure dimyristoyl-phosphatidyl-choline (DMPC) bilayer, a fusion peptide interacting with a membrane, voltage-dependent anion channel protein embedded in a DMPC bilayer, cholesterol enriched bilayer and a coarse grained simulation of a curved palmitoyl-oleoyl-phosphatidyl-choline lipid membrane. The local membrane property analysis proves to provide an intuitive and detailed view on the observables that are otherwise interpreted as averaged bilayer properties.

Effect of amphipathic HIV fusion inhibitor peptides on POPC and POPC/cholesterol membrane properties: a molecular simulation study

T-20 and T-1249 fusion inhibitor peptides were shown to interact with 1-palmitoyl-2-oleyl-phosphatidylcholine (POPC) (liquid disordered, ld) and POPC/cholesterol (1:1) (POPC/Chol) (liquid ordered, lo) bilayers, and they do so to different extents. Although they both possess a tryptophan-rich domain (TRD), T-20 lacks a pocket binding domain (PBD), which is present in T-1249. It has been postulated that the PBD domain enhances FI interaction with HIV gp41 protein and with model membranes. Interaction of these fusion inhibitor peptides with both the cell membrane and the viral envelope membrane is important for function, i.e., inhibition of the fusion process. We address this problem with a molecular dynamics approach focusing on lipid properties, trying to ascertain the consequences and the differences in the interaction of T-20 and T-1249 with ld and lo model membranes. T-20 and T-1249 interactions with model membranes are shown to have measurable and different effects on bilayer structural and dynamical parameters. T-1249’s adsorption to the membrane surface has generally a stronger influence in the measured parameters. The presence of both binding domains in T-1249 appears to be paramount to its stronger interaction, and is shown to have a definite importance in membrane properties upon peptide adsorption.

Molecular dynamics simulation of HIV fusion inhibitor T-1249: insights on peptide-lipid interaction

T-1249 is a peptide that inhibits the fusion of HIV envelope with the target cell membrane. Recent results indicate that T-1249, as in the case of related inhibitor peptide T-20 (enfuvirtide), interacts with membranes, more extensively in the bilayer liquid disordered phase than in the liquid ordered state, which could be linked to its effectiveness. Extensive molecular dynamics simulations (100 ns) were carried out to investigate the interaction between T-1249 and bilayers of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and POPC/cholesterol (1 : 1). It was observed that T-1249 interacts to different extents with both membrane systems and that peptide interaction with the bilayer surface has a local effect on membrane structure. Formation of hydrogen bonding between certain peptide residues and several acceptor and donor groups in the bilayer molecules was observed. T-1249 showed higher extent of interaction with bilayers when compared to T-20. This is most notable in POPC/Chol membranes, owing to more peptide residues acting as H bond donors and acceptors between the peptide and the bilayer lipids, including H-bonds formed with cholesterol. This behavior is at variance with that of T-20, which forms no H bonds with cholesterol. This higher ability to interact with membranes is probably correlated with its higher inhibitory efficiency.