Positioning and Stabilization of Dynorphin Peptides in Membrane Bilayers: the Mechanistic Role of Aromatic and Basic Residues Revealed from Comparative MD Simulations

Probing electrophysiological activity of amphiphilic Dynorphin A in planar neutral membranes reveals both ion channel-like activity and neuropeptide translocation

Dynorphin A (DynA) is an endogenous neuropeptide that besides acting as a ligand of the κ-opioid receptor, presents some non-opioid pathophysiological properties associated to its ability to induce cell permeability similarly to cell-penetrating peptides (CPPs). Here, we use electrophysiology experiments to show that amphiphilic DynA generates aqueous pores in neutral membranes similar to those reported previously in charged membranes, but we also find other events thermodynamically incompatible with voltage-driven ion channel activity (i.e. non-zero currents with no applied voltage in symmetric salt conditions, reversal potentials that exceed the theoretical limit for a given salt concentration gradient). By comparison with current traces generated by other amphiphilic molecule known to spontaneously cross membranes, we hypothesize that DynA could directly translocate across neutral bilayers, a feature never observed in charged membranes following the same electrophysiological protocol. Our findings suggest that DynA interaction with the cellular membrane is modulated by the lipid charge distribution, enabling either passive ionic transport via membrane remodeling and pore formation or by peptide direct internalization independent of cellular transduction pathways.

The World of GPCR dimers - Mapping dopamine receptor D2 homodimers in different activation states and configuration arrangements

G protein-coupled receptors (GPCRs) are known to dimerize, but the molecular and structural basis of GPCR dimers is not well understood. In this study, we developed a computational framework to generate models of symmetric and asymmetric GPCR dimers using different monomer activation states and identified their most likely interfaces with molecular details. We chose the dopamine receptor D2 (D2R) homodimer as a case study because of its biological relevance and the availability of structural information. Our results showed that transmembrane domains 4 and 5 (TM4 and TM5) are mostly found at the dimer interface of the D2R dimer and that these interfaces have a subset of key residues that are mostly nonpolar from TM4 and TM5, which was in line with experimental studies. In addition, TM2 and TM3 appear to be relevant for D2R dimers. In some cases, the inactive configuration is unaffected by the partnered protomer, whereas in others, the active protomer adopts the properties of an inactive receptor. Additionally, the β-arrestin configuration displayed the properties of an active receptor in the absence of an agonist, suggesting that a switch to another meta-state during dimerization occurred. Our findings are consistent with the experimental data, and this method can be adapted to study heterodimers and potentially extended to include additional proteins such as G proteins or β-arrestins. In summary, this approach provides insight into the impact of the conformational status of partnered protomers on the overall quaternary GPCR macromolecular structure and dynamics.

Theoretical insight into dihydrogen activation with β-diketiminato ligand supported Group 13 and 14 elements: mechanism and activity difference

Mechanisms and activity difference for dihydrogen activation with the BDI-supported Group 13 and 14 elements are comparatively studied.

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells' membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane's mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Membrane driven spatial organization of GPCRs

Spatial organization of G-protein coupled receptors (GPCRs) into dimers and higher order oligomers has been demonstrated in vitro and in vivo. The pharmacological readout was shown to depend on the specific interfaces, but why particular regions of the GPCR structure are involved, and how ligand-determined states change them remains unknown. Here we show why protein-membrane hydrophobic matching is attained upon oligomerization at specific interfaces from an analysis of coarse-grained molecular dynamics simulations of the spontaneous diffusion-interaction of the prototypical beta2-adrenergic (β₂AR) receptors in a POPC lipid bilayer. The energy penalty from mismatch is significantly reduced in the spontaneously emerging oligomeric arrays, making the spatial organization of the GPCRs dependent on the pattern of mismatch in the monomer. This mismatch pattern is very different for β₂AR compared to the highly homologous and structurally similar β₁AR, consonant with experimentally observed oligomerization patterns of β₂AR and β₁AR. The results provide a mechanistic understanding of the structural context of oligomerization.

The cost of living in the membrane: a case study of hydrophobic mismatch for the multi-segment protein LeuT

Many observations of the role of the membrane in the function and organization of transmembrane (TM) proteins have been explained in terms of hydrophobic mismatch between the membrane and the inserted protein. For a quantitative investigation of this mechanism in the lipid-protein interactions of functionally relevant conformations adopted by a multi-TM segment protein, the bacterial leucine transporter (LeuT), we employed a novel method, Continuum-Molecular Dynamics (CTMD), that quantifies the energetics of hydrophobic mismatch by combining the elastic continuum theory of membrane deformations with an atomistic level description of the radially asymmetric membrane-protein interface from MD simulations. LeuT has been serving as a model for structure-function studies of the mammalian neurotransmitter:sodium symporters (NSSs), such as the dopamine and serotonin transporters, which are the subject of intense research in the field of neurotransmission. The membrane models in which LeuT was embedded for these studies were composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid, or 3:1 mixture of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) lipids. The results show that deformation of the host membrane alone is not sufficient to alleviate the hydrophobic mismatch at specific residues of LeuT. The calculations reveal significant membrane thinning and water penetration due to the specific local polar environment produced by the charged K288 of TM7 in LeuT, that is membrane-facing deep inside the hydrophobic milieu of the membrane. This significant perturbation is shown to result in unfavorable polar-hydrophobic interactions at neighboring hydrophobic residues in TM1a and TM7. We show that all the effects attributed to the K288 residue (membrane thinning, water penetration, and the unfavorable polar-hydrophobic interactions at TM1a and TM7), are abolished in calculations with the K288A mutant. The involvement of hydrophobic mismatch is somewhat different in the functionally distinct conformations (outward-open, occluded, inward-open) of LeuT, and the differences are shown to connect to structural elements (e.g., TM1a) known to play key roles in transport. This finding suggests a mechanistic hypothesis for the enhanced transport activity observed for the K288A mutant, suggesting that the unfavorable hydrophobic-hydrophilic interactions hinder the motion of TM1a in the functionally relevant conformational transition to the inward-open state. Various extents of such unfavorable interactions, involving exposure to the lipid environment of adjacent hydrophobic and polar residues, are common in multi-segment transmembrane proteins, and must be considered to affect functionally relevant conformational transitions.

Quantitative modeling of membrane deformations by multihelical membrane proteins: application to G-protein coupled receptors

The interpretation of experimental observations of the dependence of membrane protein function on the properties of the lipid membrane environment calls for a consideration of the energy cost of protein-bilayer interactions, including the protein-bilayer hydrophobic mismatch. We present a novel (to our knowledge) multiscale computational approach for quantifying the hydrophobic mismatch-driven remodeling of membrane bilayers by multihelical membrane proteins. The method accounts for both the membrane remodeling energy and the energy contribution from any partial (incomplete) alleviation of the hydrophobic mismatch by membrane remodeling. Overcoming previous limitations, it allows for radially asymmetric bilayer deformations produced by multihelical proteins, and takes into account the irregular membrane-protein boundaries. The approach is illustrated by application to two G-protein coupled receptors: rhodopsin in bilayers of different thickness, and the serotonin 5-HT(2A) receptor bound to pharmacologically different ligands. Analysis of the results identifies the residual exposure that is not alleviated by bilayer adaptation, and its quantification at specific transmembrane segments is shown to predict favorable contact interfaces in oligomeric arrays. In addition, our results suggest how distinct ligand-induced conformations of G-protein coupled receptors may elicit different functional responses through differential effects on the membrane environment.

Cholesterol modulates the membrane effects and spatial organization of membrane-penetrating ligands for G-protein coupled receptors

The ligands of certain G-protein coupled receptors (GPCRs) are membrane soluble and reach their target from the lipid bilayer. Lipid composition and dynamics will therefore modulate the activity of these receptors, but specific roles of lipid components, including the ubiquitous cholesterol (Chol), are not clear. We have probed the organization and dynamics of such a lipid-bilayer-penetrating ligand, the endogenous ligand for the κ-opioid receptor (KOR) dynorphin A (1-17) (DynA), using molecular dynamics (MD) simulations of DynA in cholesterol-depleted and cholesterol-enriched model membranes. DynA is found to penetrate deep inside fluid dimyristoylphosphatidylcholine (DMPC) bilayers, and resides with its N-terminal helix at ∼6 Å away from the bilayer midplane, in a tilted orientation, at an ∼50° angle with respect to the membrane normal. In contrast, DynA inside DMPC/Chol membranes with 20% cholesterol (DMPC/Chol) is situated with its helical segment ∼5 A higher, i.e., closer to the lipid/water interface and in a relatively vertical orientation. The DMPC membrane shows greater thinning around the insertion and permits a stronger influx of water inside the hydrocarbon core than the DMPC/Chol membranes. Relating these results to data about key GPCR residues that have been implicated in interactions with membrane-inserting GPCR ligands, we conclude that the position of DynA in DMPC/Chol, but not in pure DMPC, correlates with generally proposed GPCR ligand entry pathways. Our predictions provide a possible mechanistic explanation as to why DynA binding to KOR, and the subsequent activation of the receptor, is facilitated in cholesterol-enriched environments. A quantitative description of DynA-induced membrane deformations is obtained with a continuum theory of membrane deformations (CTMD) that is based on hydrophobic matching. Comparison with the MD data reveals the significance of the lipid tail packing energy contribution in the DMPC/Chol mixtures in predicting equilibrium membrane shape around DynA. On this basis, specific corrections are introduced to this energy term within the CTMD framework, thereby extending the applicability of the CTMD framework to lipid raft mixtures and their interactions with GPCR proteins and their ligands.

Force field dependence of phospholipid headgroup and acyl chain properties: comparative molecular dynamics simulations of DMPC bilayers

The reliability of molecular simulations largely depends on the quality of the empirical force field parameters. Force fields used in lipid simulations continue to be improved to enhance the agreement with experiments for a number of different properties. In this work, we have carried out molecular dynamics simulations of neat DMPC bilayers using united-atom Berger force field and three versions of all-atom CHARMM force fields. Three different systems consisting of 48, 72, and 96 lipids were studied. Both particle mesh Ewald (PME) and spherical cut-off schemes were used to evaluate the long-range electrostatic interactions. In total, 21 simulations were carried out and analyzed to find out the dependence of lipid properties on force fields, system size, and schemes to calculate long-range interactions. The acyl chain order parameters calculated from Berger and the recent versions of CHARMM simulations have shown generally good agreement with the experimental results. However, both sets of force fields deviate significantly from the experimentally observed P-C dipolar coupling values for the carbon atoms that link the choline and glycerol groups with the phosphate groups. Significant differences are also observed in several headgroup parameters between CHARMM and Berger simulations. Our results demonstrate that when changes were introduced to improve CHARMM force field using PME scheme, all the headgroup parameters have not been reoptimized. The headgroup properties are likely to play a significant role in lipid-lipid, protein-lipid, and ligand-lipid interactions and hence headgroup parameters in phospholipids require refinement for both Berger and CHARMM force fields.

Dynamic models of G-protein coupled receptor dimers: indications of asymmetry in the rhodopsin dimer from molecular dynamics simulations in a POPC bilayer

Based on the growing evidence that G-protein coupled receptors (GPCRs) form homo- and hetero-oligomers, models of GPCR signaling are now considering macromolecular assemblies rather than monomers, with the homo-dimer regarded as the minimal oligomeric arrangement required for functional coupling to the G-protein. The dynamic mechanisms of such signaling assemblies are unknown. To gain some insight into properties of GPCR dimers that may be relevant to functional mechanisms, we study their current structural prototype, rhodopsin. We have carried out nanosecond time-scale molecular dynamics (MD) simulations of a rhodopsin dimer and compared the results to the monomer simulated in the same type of bilayer membrane model composed of an equilibrated unit cell of hydrated palmitoyl-oleoyl-phosphatidyl choline (POPC). The dynamic representation of the homo-dimer reveals the location of structural changes in several regions of the monomeric subunits. These changes appear to be more pronounced at the dimerization interface that had been shown to be involved in the activation process [Proc Natl Acad Sci USA 102:17495, 2005]. The results are consistent with a model of GPCR activation that involves allosteric modulation through a single GPCR subunit per dimer.

Recognition of GPCRs by Peptide Ligands and Membrane Compartments theory: Structural Studies of Endogenous Peptide Hormones in Membrane Environment

One of the largest family of cell surface proteins, G-protein coupled receptors (GPCRs) regulate virtually all known physiological processes in mammals. With seven transmembrane segments, they respond to diverse range of extracellular stimuli and represent a major class of drug targets. Peptidergic GPCRs use endogenous peptides as ligands. To understand the mechanism of GPCR activation and rational drug design, knowledge of three-dimensional structure of receptor–ligand complex is important. The endogenous peptide hormones are often short, flexible and completely disordered in aqueous solution. According to “Membrane Compartments Theory”, the flexible peptide binds to the membrane in the first step before it recognizes its receptor and the membrane-induced conformation is postulated to bind to the receptor in the second step. Structures of several peptide hormones have been determined in membrane-mimetic medium. In these studies, micelles, reverse micelles and bicelles have been used to mimic the cell membrane environment. Recently, conformations of two peptide hormones have also been studied in receptor-bound form. Membrane environment induces stable secondary structures in flexible peptide ligands and membrane-induced peptide structures have been correlated with their bioactivity. Results of site-directed mutagenesis, spectroscopy and other experimental studies along with the conformations determined in membrane medium have been used to interpret the role of individual residues in the peptide ligand. Structural differences of membrane-bound peptides that belong to the same family but differ in selectivity are likely to explain the mechanism of receptor selectivity and specificity of the ligands. Knowledge of peptide 3D structures in membrane environment has potential applications in rational drug design.

Hallucinogen actions on 5-HT receptors reveal distinct mechanisms of activation and signaling by G protein-coupled receptors

We review the effect of some key advances in the characterization of molecular mechanisms of signaling by G protein-coupled receptors (GPCRs) on our current understanding of mechanisms of drugs of abuse. These advances are illustrated by results from our ongoing work on the actions of hallucinogens on serotonin (5-HT) receptors. We show how a combined computational and experimental approach can reveal specific modes of receptor activation underlying the difference in properties of hallucinogens compared with nonhallucinogenic congeners. These modes of activation-that can produce distinct ligand-dependent receptor states-are identified in terms of structural motifs (SM) in molecular models of the receptors, which were shown to constitute conserved functional microdomains (FM). The role of several SM/FMs in the activation mechanism of the GPCRs is presented in detail to illustrate how this mechanism can lead to ligand-dependent modes of signaling by the receptors. Novel bioinformatics tools are described that were designed to support the quantitative mathematical modeling of ligand-specific signaling pathways activated by the 5-HT receptors targeted by hallucinogens. The approaches for mathematical modeling of signaling pathways activated by 5-HT receptors are described briefly in the context of ongoing work on detailed biochemical models of 5-HT2A, and combined 5-HT2A/5-HT1A, receptor-mediated activation of the MAPK 1,2 pathway. The continuing need for increasingly more realistic representation of signaling in dynamic compartments within the cell, endowed with spatio-temporal characteristics obtained from experiment, is emphasized. Such developments are essential for attaining a quantitative understanding of how the multiple functions of a cell are coordinated and regulated, and to evaluate the specifics of the perturbations caused by the drugs of abuse that target GPCRs.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.