Molecular dynamics simulations predict a tilted orientation for the helical region of dynorphin A(1-17) in dimyristoylphosphatidylcholine bilayers

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.

Functional Characterization of Spinocerebellar Ataxia Associated Dynorphin A Mutant Peptides

Mutations in the prodynorphin gene (PDYN) are associated with the development of spinocerebellar ataxia type 23 (SCA23). Pathogenic missense mutations are localized predominantly in the PDYN region coding for the dynorphin A (DynA) neuropeptide and lead to persistently elevated mutant peptide levels with neurotoxic properties. The main DynA target in the central nervous system is the kappa opioid receptor (KOR), a member of the G-protein coupled receptor family, which can elicit signaling cascades mediated by G-protein dissociation as well as β-arrestin recruitment. To date, a thorough analysis of the functional profile for the pathogenic SCA23 DynA mutants at KOR is still missing. To elucidate the role of DynA mutants, we used a combination of assays to investigate the differential activation of G-protein subunits and β-arrestin. In addition, we applied molecular modelling techniques to provide a rationale for the underlying mechanism. Our results demonstrate that DynA mutations, associated with a severe ataxic phenotype, decrease potency of KOR activation, both for G-protein dissociation as well as β-arrestin recruitment. Molecular modelling suggests that this loss of function is due to disruption of critical interactions between DynA and the receptor. In conclusion, this study advances our understanding of KOR signal transduction upon DynA wild type or mutant peptide binding.

Insight into the interactions, residue snorkeling, and membrane disordering potency of a single antimicrobial peptide into different lipid bilayers

Pardaxin, with a bend-helix-bend-helix structure, is a membrane-active antimicrobial peptide that its membrane activity depends on the lipid bilayer composition. Herein, all-atom molecular dynamics (MD) simulations were performed to provide further molecular insight into the interactions, structural dynamics, orientation behavior, and cationic residues snorkeling of pardaxin in the DMPC, DPPC, POPC, POPG, POPG/POPE (3:1), and POPG/POPE (1:3) lipid bilayers. The results showed that the C-terminal helix of the peptide was maintained in all six types of the model-bilayers and pardaxin was tilted into the DMPC, DPPC, and POPG/POPE mixed bilayers more than the POPC and POPG bilayers. As well as, the structure of zwitterionic membranes was more affected by the peptide than the anionic bilayers. Taken together, the study demonstrated that the cationic residues of pardaxin snorkeled toward the interface of lipid bilayers and all phenylalanine residues of the peptide played important roles in the peptide-membrane interactions. We hope that this work will provide a better understanding of the interactions of antimicrobial peptides with the membranes.

Altered secondary structure of Dynorphin A associates with loss of opioid signalling and NMDA-mediated excitotoxicity in SCA23

Spinocerebellar ataxia type 23 (SCA23) is caused by missense mutations in prodynorphin, encoding the precursor protein for the opioid neuropeptides α-neoendorphin, Dynorphin (Dyn) A and Dyn B, leading to neurotoxic elevated mutant Dyn A levels. Dyn A acts on opioid receptors to reduce pain in the spinal cord, but its cerebellar function remains largely unknown. Increased concentration of or prolonged exposure to Dyn A is neurotoxic and these deleterious effects are very likely caused by an N-methyl-d-aspartate-mediated non-opioid mechanism as Dyn A peptides were shown to bind NMDA receptors and potentiate their glutamate-evoked currents. In the present study, we investigated the cellular mechanisms underlying SCA23-mutant Dyn A neurotoxicity. We show that SCA23 mutations in the Dyn A-coding region disrupted peptide secondary structure leading to a loss of the N-terminal α-helix associated with decreased κ-opioid receptor affinity. Additionally, the altered secondary structure led to increased peptide stability of R6W and R9C Dyn A, as these peptides showed marked degradation resistance, which coincided with decreased peptide solubility. Notably, L5S Dyn A displayed increased degradation and no aggregation. R6W and wt Dyn A peptides were most toxic to primary cerebellar neurons. For R6W Dyn A, this is likely because of a switch from opioid to NMDA- receptor signalling, while for wt Dyn A, this switch was not observed. We propose that the pathology of SCA23 results from converging mechanisms of loss of opioid-mediated neuroprotection and NMDA-mediated excitotoxicity.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Structure and orientation of bovine lactoferrampin in the mimetic bacterial membrane as revealed by solid-state NMR and molecular dynamics simulation

Bovine lactoferrampin (LFampinB) is a newly discovered antimicrobial peptide found in the N1-domain of bovine lactoferrin (268-284), and consists of 17 amino-acid residues. It is important to determine the orientation and structure of LFampinB in bacterial membranes to reveal the antimicrobial mechanism. We therefore performed (13)C and (31)P NMR, (13)C-(31)P rotational echo double resonance (REDOR), potassium ion-selective electrode, and quartz-crystal microbalance measurements for LFampinB with mimetic bacterial membrane and molecular-dynamics simulation in acidic membrane. (31)P NMR results indicated that LFampinB caused a defect in mimetic bacterial membranes. Ion-selective electrode measurements showed that ion leakage occurred for the mimetic bacterial membrane containing cardiolipin. Quartz-crystal microbalance measurements revealed that LFampinB had greater affinity to acidic phospholipids than that to neutral phospholipids. (13)C DD-MAS and static NMR spectra showed that LFampinB formed an α-helix in the N-terminus region and tilted 45° to the bilayer normal. REDOR dephasing patterns between carbonyl carbon nucleus in LFampinB and phosphorus nuclei in lipid phosphate groups were measured by (13)C-(31)P REDOR and the results revealed that LFampinB is located in the interfacial region of the membrane. Molecular-dynamics simulation showed the tilt angle to be 42° and the rotation angle to be 92.5° for Leu(3), which are in excellent agreement with the experimental values.

Aurein 2.3 functionality is supported by oblique orientated α-helical formation

In this study, an amphibian antimicrobial peptide, aurein 2.3, was predicted to use oblique orientated α-helix formation in its mechanism of membrane destabilisation. Molecular dynamic (MD) simulations and circular dichroism (CD) experimental data suggested that aurein 2.3 exists in solution as unstructured monomers and folds to form predominantly α-helical structures in the presence of a dimyristoylphosphatidylcholine membrane. MD showed that the peptide was highly surface active, which supported monolayer data where the peptide induced surface pressure changes>34 mNm(-1). In the presence of a lipid membrane MD simulations suggested that under hydrophobic mismatch the peptide is seen to insert via oblique orientation with a phenylalanine residue (PHE3) playing a key role in the membrane interaction. There is evidence of snorkelling leucine residues leading to further membrane disruption and supporting the high level of lysis observed using calcein release assays (76%). Simulations performed at higher peptide/lipid ratio show peptide cooperativity is key to increased efficiency leading to pore-formation.

Molecular dynamics simulation of Bombolitin II in the dipalmitoylphosphatidylcholine membrane bilayer

The orientation behavior of Bombolitin II (BLT2) in the dipalmitoylphosphatidylcholine membrane bilayer was investigated by using molecular-dynamics simulation. During the 20-ns simulation, the BLT2 began to tilt and finally reached the angle of 51° from the membrane-normal. The structure of the peptide formed the amphipathic α-helical structure during the entire simulation time. The peptide tilts with its hydrophobic side faced to the hydrophobic core of the bilayer. We analyzed the mechanism of the tilting behavior of the peptide associated with the membrane in detail. The analysis showed that the hydrogen-bond interaction and the electrostatic interaction were found to exist between Lys(12) and a lipid molecule. These interactions are considered to work as an important factor in tilting the peptide to the membrane-normal.

COMPUTER SIMULATIONS OF MEMBRANE-LYTIC PEPTIDES: PERSPECTIVES IN DRUG DESIGN

Structure activity relationships were investigated for membrane-lytic peptides (MLP) Ltc1 and Ltc2a from the latarcin family. The peptides were studied via long-term molecular dynamics (MD) simulations in different membrane environments (detergent micelles, mixed lipid bilayers mimiking eukaryotic and bacterial membranes). The calculated structure of Ltc2a in sodium dodecyl sulfate micelle agrees well with the data obtained by 1H-NMR spectroscopy. This validates the applied modeling approach. The binding mode of MLPs is governed by several factors: (i) the membrane surface curvature; (ii) the conformational plasticity and hydrophobic organization of the peptide, which depend on the arrangement of charged, non-polar and helix-breaking residues in the amino acid sequence. In contrast to Ltc1, insertion of Ltc2a into model membranes induces significant changes in dynamic behavior of lipids in the contact region. Such a prominent membrane destabilization correlates with high membrane-lytic activity of Ltc2a. In all cases the "membrane response" has a local character and is caused by formation of specific peptide-lipid contacts. Results of MD simulations of Ltc2a in model membranes were used to develop a number of its analogs with predefined activity.

Dynamic structure of bombolitin II bound to lipid bilayers as revealed by solid-state NMR and molecular-dynamics simulation

Bombolitin II (BLT2) is one of the hemolytic heptadecapeptides originally isolated from the venom of a bumblebee. Structure and orientation of BLT2 bound to 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) membranes were determined by solid-state (31)P and (13)C NMR spectroscopy. (31)P NMR spectra showed that BLT2-DPPC membranes were disrupted into small particles below the gel-to-liquid crystalline phase transition temperature (T(c)) and fused to form a magnetically oriented vesicle system where the membrane surface is parallel to the magnetic fields above the T(c). (13)C NMR spectra of site-specifically (13)C-labeled BLT2 at the carbonyl carbons were observed and the chemical shift anisotropies were analyzed to determine the dynamic structure of BLT2 bound to the magnetically oriented vesicle system. It was revealed that the membrane-bound BLT2 adopted an α-helical structure, rotating around the membrane normal with the tilt angle of the helical axis at 33°. Interatomic distances obtained from rotational-echo double-resonance experiments further showed that BLT2 adopted a straight α-helical structure. Molecular dynamics simulation performed in the BLT2-DPPC membrane system showed that the BLT2 formed a straight α-helix and that the C-terminus was inserted into the membrane. The α-helical axis is tilted 30° to the membrane normal, which is almost the same as the value obtained from solid-state NMR. These results suggest that the membrane disruption induced by BLT2 is attributed to insertion of BLT2 into the lipid bilayers.

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.

Is there a difference in metal ion-based inhibition between members of thionin family: Molecular dynamics simulation study

Thionins have a considerable potential as antimicrobial compounds although their application may be restricted by metal ion-based inhibition of membrane permeabilizing activity. We previously reported the properties associated with the proposed mechanism of metal ion-based inhibition of beta-purothionin. In this study, we investigated the effects of metal ions on alpha-hordothionin which differs from beta-purothionin by eight out of 45 residues. Three of the differing residues are thought to be involved in the mechanism of metal ion-based inhibition in beta-purothionin. The structure and dynamics of alpha-hordothionin were explored using unconstrained molecular dynamics (MD) simulations in explicit water as a function of metal ions. Although the global fold is almost identical to that of beta-purothionin, alpha-hordothionin displays reduced fluctuating motions. Moreover, alpha-hordothionin is more resistant to the presence of metal ions than beta-purothionin. Mg(+2) ions do not affect alpha-hordothionin, whereas K(+) ions induce perturbations in the alpha2 helix, modify dynamics and electrostatic properties. Nevertheless, these changes are considerably smaller than those in beta-purothionin. The proposed mechanism of metal ion-based inhibition involves the hydrogen bonding network of Arg5-Arg30-Gly27, which regulates dynamic unfolding of the alpha2 C-end which is similar to beta-purothionin response. The key residues responsible for the increased resistance for alpha-hordothionin are Gly27 and Gly42 which replace Asn27 and Asp42 involved into the mechanism of metal ion-based inhibition in beta-purothionin. Comparison of MD simulations of alpha-hordothionin with beta-purothionin reveals dynamic properties which we believe are intrinsic properties of thionins with four disulphide bonds.

On the thermodynamic stability of a charged arginine side chain in a transmembrane helix

Biological membranes consist of bilayer arrangements of lipids forming a hydrophobic core that presents a physical barrier to all polar and charged molecules. This long-held notion has recently been challenged by biological translocon-based experiments that report small apparent free energies to insert charged side chains near the center of a transmembrane (TM) helix. We have carried out fully atomistic simulations to provide the free-energy profile for moving a TM helix containing a protonated Arg side chain across a lipid bilayer. Our results reveal the fundamental thermodynamics governing the stability of charged side chains in membranes and the microscopic interactions involved. Despite local membrane deformations, where large amounts of water and lipid head groups are pulled into the bilayer to interact with Arg, the free-energy barrier is 17 kcal/mol. We provide a rationale for the differences in our microscopic free energies and cell biological experiments using free-energy calculations that indicate that a protonated Arg at the central residue of a TM helix of the Leader peptidase might reside close to the interface and not at the membrane center. Our findings have implications for the gating mechanisms of voltage-gated ion channels, suggesting that movements of protonated Arg residues through the membrane will be prohibited.

Conformational analysis of dynorphin A (1-13) using hydrogen-deuterium exchange and tandem mass spectrometry

Trifluoroethanol (TFE)-induced conformational changes in dynorphin A (1-13) were investigated using charge-state distribution (CSD) and hydrogen-deuterium exchange (HDX), combined with electrospray ionization (ESI) mass spectrometry (MS). Individual amino acids involved in secondary structural elements were identified by collision-induced dissociation-tandem mass spectrometry (MS/MS). It was observed that dynorphin A (1-13) largely exists in an unfolded conformation and a folded structure in increasing concentrations of TFE. In 50% TFE, it forms an alpha-helix that encompasses residues 1-9 and remains flexible from residues 10 to 13.

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.

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.

Role of lysine residues in membrane anchoring of saposin C

Molecular dynamics (MD) simulations of the N-terminal region of saposin C, containing amino acid residues 4-20 (saposin C4-20), were performed over 2.5 ns in 1,2-dioleoyl-sn-glycero-3-phosphoserine (DOPS) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) monolayers. The simulations revealed several strong specific interactions of lysine 13 (Lys13) and lysine 17 (Lys17) in saposin C4-20 with the anionic phospholipids, which are required for membrane anchoring of the peptide. Membrane anchoring of saposin C4-20 facilitates saposin C-induced liposomal membrane fusion. Substitutions of Lys13 or Lys17 with alanine or glutamic acid led to a substantial loss of saposin C's fusogenicity. However, arginine replacement of Lys13 or Lys17 caused a partial loss of saposin C's fusogenic activity. The membrane anchoring of saposin C was altered in the presence of 0.4 M sodium chloride. Differential salt effects on Lys-mutant saposin Cs were observed using Trp fluorescence analysis. Low salt concentration had a more significant impact on Lys-mutant saposin C with a negatively charged amino acid residue replacement than those mutants with a positively charged or neutral residue replacement. These results indicate that positively charged amino acids at positions 13 and 17 are required for the fusogenic function of saposin C. In addition, the side-chain structure of lysine is crucial to the precise membrane anchoring which is necessary for the total fusion activity of saposin C. The MD simulations and vesicle size measurements of lysine-mutant saposins confirm the importance of the two lysine residues in saposin C4-20 for saposin C-induced fusion of negatively charged phospholipid membranes.

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.

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.

Membrane localization and flexibility of a lipidated ras peptide studied by molecular dynamics simulations

Lipid-modified membrane-binding proteins are essential in signal transduction events of the cell, a typical example being the GTPase ras. Recently, membrane binding of a doubly lipid-modified heptapeptide from the C-terminus of the human N-ras protein was studied by spectroscopic techniques. It was found that membrane binding is mainly due to lipid chain insertion, but it is also favored by interactions between apolar side chains and the hydrophobic region of the membrane. Here, 10 explicit solvent molecular dynamics simulations for a total time of about 150 ns are used to investigate the atomic details of the peptide-membrane association. The 16:0 peptide lipid chains are more mobile than the 14:0 phospholipid chains, which is in agreement with (2)H NMR experiments. Peptide-lipid and peptide-solvent interactions, backbone and side-chain distributions, as well as the effects of lipidated peptide insertion onto the structure, and dynamics of a 1,2-dimyristoylglycero-3-phosphocholine bilayer are described. The simulation results validate the structural model proposed by the analysis of spectroscopic data and highlight the main aspects of the insertion mechanism. The peptide in the membrane is rather rigid over the simulation time scale of about 10 ns, but different partially extended conformations devoid of backbone hydrogen bonds are observed in different trajectories.

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.

Solid-state NMR investigation of the depth of insertion of protegrin-1 in lipid bilayers using paramagnetic Mn2+

The depth of insertion of an antimicrobial peptide, protegrin-1 (PG-1), in lipid bilayers is investigated using solid-state NMR. Paramagnetic Mn(2+) ions bind to the surface of lipid bilayers and induce distance-dependent dipolar relaxation of nuclear spins. By comparing the signal dephasing of the peptide with that of the lipids, whose segmental depths of insertion are known, we determined the depths of several residues of PG-1 in 1,2 dilauryl-sn-glycero-3-phosphotidylcholine (DLPC) bilayers. We found that residues G2 at the N-terminus and F12 at the beta-turn of the peptide reside near the membrane surface, whereas L5 and V16 are embedded in the acyl chain region. The depths increase in the order of G2 < F12 < L5 < V16. These intensity-dephasing results are confirmed by direct measurement of the paramagnetically enhanced (13)C transverse relaxation rates. The relative depths indicate that PG-1 is tilted from the bilayer normal, which is consistent with independent solid-state NMR measurements of PG-1 orientation in the same lipids (Yamaguchi et al., 2001). They also indicate that PG-1 is fully immersed in the lipid bilayer. However, a quantitative mismatch between the bilayer thickness and PG-1 length suggests a local thinning of the DLPC bilayer by 8-10 A. The depth sensitivity of this Mn(2+) dephasing technique is tunable with the Mn(2+) concentration to focus on different regions of the lipid bilayer.

Molecular dynamics simulations of the transmembrane domain of the oncogenic ErbB2 receptor dimer in a DMPC bilayer

Molecular dynamics simulations of an atomic model of the transmembrane domain of the oncogenic ErbB2 receptor dimer embedded in an explicit dimyristoylphosphatidylcholine (DMPC) bilayer were performed for more than 4 ns. The oncogenic Glu mutation in the membrane spanning segment plays a major role in tyrosine kinase activity and receptor dimerization, and is thought to be partly responsible for the structure of the transmembrane domain of the active receptor. MD results show that the interactions between the two transmembrane helices are characteristic of a left-handed packing as previously demonstrated from in vacuo simulations. Moreover, MD results reveal the absence of persistent hydrogen bonds between the Glu side chains in a membrane environment, which raise the question of the ability for Glu alone to stabilize the TM domain of the ErbB2 receptor. Interestingly the formation of the alpha-pi motif in the two ErbB2 transmembrane helices confirms the concept of intrinsic sequence-induced conformational flexibility. From a careful analysis of our MD results, we suggest that the left-handed helix-helix packing could be the key to correctly orient the intracellular domain of the activated receptor dimer. The prediction of such interactions from computer simulations represents a new step towards the understanding of signaling mechanisms.

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.

Specific mode of interaction between components of model pulmonary surfactants using computer simulations

Atomistic molecular dynamics simulations and structural bioinformatics tools enable the identification of the exact mode of interaction between model pulmonary surfactant components. Two nanosecond long simulations of the N-terminal region of human surfactant protein-B (SP-B(1-25)) in dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylglycerol (DPPG) monolayers of different lipid surface densities reveal the preferential affinity of SP-B(1-25) for anionic phospholipids. In particular, arginine 12 and lysine 24 interact strongly and with high specificity with the phosphate group of the DPPG lipids, stabilizing the position, the orientation, and the secondary structure of the peptide in the monolayer. The peptide lies at an oblique angle to the interfacial plane, ranging between 47 degrees and 62 degrees, increasing with decreasing lipid surface density. In DPPC monolayers the interaction is largely determined by hydrophobic interactions. The non-specific nature of DPPC-SP-B(1-25) interactions allows for significant flexibility in the topology of the peptide in the lipid matrix. Bioinformatics tools are employed to generalize the simulation results to the sequences of SP-B(1-25) in other organisms. The importance of specific residues, and the role of the largely helical and amphiphilic nature of the peptide in the functionality of SP-B(1-25) are established. The synergy of classical mechanics tools with bioinformatics methods greatly enhances the molecular-level interpretation of pulmonary surfactant action and facilitates the development of design rules for synthetic surfactant analogues.