Membrane binding of lipidated Ras peptides and proteins--the structural point of view

Biological membranes are interesting interfaces, at which important biological processes occur. In addition to integral membrane proteins, a number of proteins bind to the membrane surface and associate with it. Posttranslational lipid modification is one important mechanism, by which soluble molecules develop a propensity towards the membrane and reversibly bind to it. Membrane binding by insertion of hydrophobic lipid moieties is relevant for up to 10% of all cellular proteins. A particular interesting lipid-modified protein is the small GTPase Ras, which plays a key role in cellular signal transduction. Until recently, the structural basis for membrane binding of Ras was not well-defined. However, with the advent of new synthesis techniques and the advancement of several biophysical methods, a number of structural and dynamical features about membrane binding of Ras proteins have been revealed. This review will summarize the chemical biology of Ras and discuss in more detail the biophysical and structural features of the membrane bound C-terminus of the protein.     

Spectral fitting for signal assignment and structural analysis of uniformly <Superscript>13</Superscript>C-labeled solid proteins by simulated annealing based on chemical shifts and spin dynamics

We describe an approach for the signal assignment and structural analysis with a suite of two-dimensional (13)C-(13)C magic-angle-spinning solid-state NMR spectra of uniformly (13)C-labeled peptides and proteins. We directly fit the calculated spectra to experimental ones by simulated annealing in restrained molecular dynamics program CNS as a function of atomic coordinates. The spectra are calculated from the conformation dependent chemical shift obtained with SHIFTX and the cross-peak intensities computed for recoupled dipolar interactions. This method was applied to a membrane-bound 14-residue peptide, mastoparan-X. The obtained C', C(alpha) and C(beta) chemical shifts agreed with those reported previously at the precisions of 0.2, 0.7 and 0.4 ppm, respectively. This spectral fitting program also provides backbone dihedral angles with a precision of about 50 degrees from the spectra even with resonance overlaps. The restraints on the angles were improved by applying protein database program TALOS to the obtained chemical shifts. The peptide structure provided by these restraints was consistent with the reported structure at the backbone RMSD of about 1 A.     

Membrane interactions and dynamics of a 21-mer cytotoxic peptide: A solid-state NMR study

We have investigated the membrane interactions and dynamics of a 21-mer cytotoxic model peptide that acts as an ion channel by solid-state NMR spectroscopy. To shed light on its mechanism of membrane perturbation, ³¹P and ²H NMR experiments were performed on 21-mer peptide-containing bicelles. ³¹P NMR results indicate that the 21-mer peptide stabilizes the bicelle structure and orientation in the magnetic field and perturbs the lipid polar head group conformation. On the other hand, ²H NMR spectra reveal that the 21-mer peptide orders the lipid acyl chains upon binding. ¹⁵N NMR experiments performed in DMPC bilayers stacked between glass plates also reveal that the 21-mer peptide remains at the bilayer surface. ¹⁵N NMR experiments in perpendicular DMPC bicelles indicate that the 21-mer peptide does not show a circular orientational distribution in the bicelle planar region. Finally, ¹³C NMR experiments were used to study the 21-mer peptide dynamics in DMPC multilamellar vesicles. By analyzing the ¹³CO spinning sidebands, the results show that the 21-mer peptide is immobilized upon membrane binding. In light of these results, we propose a model of membrane interaction for the 21-mer peptide where it lies at the bilayer surface and perturbs the lipid head group conformation.     

Lessons from computer simulations of Ras proteins in solution and in membrane

Background

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

Scope of review

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

Major conclusions

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

General significance

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.     

Ca-induced stimulation of the membrane binding of Escherichia coli SecA and its association with signal peptides of secretory proteins

Previously, it was found that Ca²⁺ stimulates the intrinsic Escherichia coli SecA ATPase activity [Kim et al., FEBS Lett. 493 (2001) 12-16]. Now, we suggest that Ca²⁺ is required for efficient interaction of SecA with membranes and the signal peptide of ribose-binding protein. When the amount of external Ca²⁺ was enhanced, the amounts of membrane-bound SecA and its lipid/ATPase activity increased. In the presence of entrapped Ca²⁺ in liposomes, the binding was also stimulated in a Ca²⁺ concentration-dependent manner. The effect of Ca²⁺ on the functional regulation of SecA was also evident in the presence of the signal peptides of secretory proteins, which the interaction of SecA with the signal peptide increased with increasing Ca²⁺ concentration in the presence of membranes. However, other divalent cations including Mg²⁺, Mn²⁺, and Zn²⁺ had inhibitory or no effect, suggesting a specific role of Ca²⁺ in SecA interaction with lipid bilayers and signal peptides.     

The Solution Structure of a Cyclic Analog of Neuropeptide Y with High Y<Subscript>1</Subscript> Receptor Affinity by NMR,CD and MD Simulations

The conformation of a cyclic analog of neuropeptide Y [Tyr¹--Lys--Gly--Arg--cyclo^5/8-(Glu⁵--Tyr--Ile--Lys⁸)--Leu--Ile¹⁰--Thr--Arg--Pro--Arg--Tyr¹⁵--NH₂; cEK-NPY] with high Y₁ receptor affinity was studied using ¹H, ¹³C and ¹⁵N 2D-NMR and CD in three diverse media-viz. DMSO-d₆, water (pH 4.0) and 50% hexafluoroacetone (HFA). The conformation of cEK-NPY was interpreted based on chemical shift (¹H, ¹³C and ¹⁵N), temperature coefficients of the NH chemical shifts, ³J_NHα coupling constants and the pattern of intra and inter-residue NOE’s and the CD spectrum. In both DMSO and water, there is a preponderance of a β-strand structure, while HFA promotes an α-helical structure, which is discontinuous in the mid-region of the peptide, due to the constraints of the lactam ring. The solution structures were generated using Restrained Molecular Dynamics simulations and further refined by Mardigras to R factors between 0.55 and 0.65. The role of its conformations in its biological activity is discussed.     

The effect of calcium on the properties of charged phospholipid bilayers

We have performed molecular dynamics simulations to investigate the structure and dynamics of charged bilayers as well as the distribution of counterions at the bilayer interface. For this, we have considered the negatively charged di-myristoyl-phosphatidyl-glycerol (DMPG) and di-myristoyl-phosphatidyl-serine (DMPS) bilayers as well as a protonated di-myristoyl-phosphatidyl-serine (DMPSH) bilayer. We were particularly interested in calcium ions due to their important role in biological systems. Simulations performed in the presence of calcium ions (DMPG, DMPS) or sodium ions (DMPS) were run for 45-60 ns. Simulation results for DMPG are compared with fluorescence measurements. The average areas per molecule were 47.4 ± 0.5 Å² (DMPG with calcium), 47.3 ± 0.5 Å² (DMPS with calcium), 51.3 ± 1.0 Å² (DMPS with sodium) and 45.3 ± 0.5 Å² (DMPSH). The structure of the negatively charged lipids is significantly affected by the counterions, where calcium ions have a more pronounced effect than sodium ions. Calcium ions were found to be tightly bound to the anionic groups of the lipid molecules and as such appear to constitute an integral part of the membrane interface on nanoseconds time scales. In contrast to sodium ions, calcium ions are localised in a narrow (∼ 10 Å) band around the phosphate group. The interaction of calcium with the lipid molecules enhances the molecular packing of the PG and PS lipids. This observation is in good agreement with emission spectra of the membrane partitioning probe Laurdan in DMPG multilamellar vesicles that indicate an increase in the ordering of the DMPG bilayer due to the presence of calcium. Our results indicate that calcium ions, which often function as a second messengers in living cells have a pronounced effect on membrane structures, which may have implications during signal transduction events.     

Molecular modeling study on the dynamical structural features of human smoothened receptor and binding mechanism of antagonist LY2940680 by metadynamics simulation and free energy calculation

Background

The smoothened (SMO) receptor, one of the Class F G protein coupled receptors (GPCRs), is an essential component of the canonical hedgehog signaling pathway which plays a key role in the regulation of embryonic development in animals. The function of the SMO receptor can be modulated by small-molecule agonists and antagonists, some of which are potential antitumour agents. Understanding the binding mode of an antagonist in the SMO receptor is crucial for the rational design of new antitumour agents.

Methods

Molecular dynamics (MD) simulation and dynamical network analysis are used to study the dynamical structural features of SMO receptor. Metadynamics simulation and free energy calculation are employed to explore the binding mechanism between the antagonist and SMO receptor.

Results

The MD simulation results and dynamical network analysis show that the conserved KTXXXW motif in helix VIII has strong interaction with helix I. The α-helical extension of transmembrane 6 (TM6) is detected as part of the ligand-binding pocket and dissociation pathway of the antagonist. The metadynamics simulation results illustrate the binding mechanism of the antagonist in the pocket of SMO receptor, and free energy calculation shows the antagonist needs to overcome about 38 kcal/mol of energy barrier to leave the binding pocket of SMO receptor.

Conclusions

The unusually long TM6 plays an important role on the binding behavior of the antagonist in the pocket of SMO receptor.

General significance

The results can not only profile the binding mechanism between the antagonist and Class F GPCRs, but also supply the useful information for the rational design of a more potential small molecule antagonist bound to SMO receptor.     

The binding of cytochrome c to neuroglobin: a docking and surface plasmon resonance study

It has recently been proposed that the role of neuroglobin in the protection of neurons from ischaemia induced cell death requires the formation of a transient complex with cytochrome c. No such complex has yet been isolated. Here, we present the results of soft docking calculations, which indicate one major binding site for cytochrome c to neuroglobin. The results yield a plausible structure for the most likely complex structure in which the hemes of each protein are in close contact. NMR analysis identifies the formation of a weak complex in which the heme group of cytochrome c is involved. surface plasmon resonance studies provide a value of 45muM for the equilibrium constant for cytochrome c binding to neuroglobin, which increases significantly as the ionic strength of the solution increases. The temperature dependence of the binding constant indicates that the complex formation is associated with a small unfavourable enthalpy change (1.9kcalmol(-1)) and a moderately large, favourable entropy change (14.8calmol(-1)deg(-1)). The sensitivity of the binding constant to the presence of salt suggests that the complex formation involves electrostatic interactions.     

Biological and structural characteristics of the binding peptides from the sporozoite proteins essential for cell traversal (SPECT)-1 and -2

The sporozoite microneme proteins essential for cell traversal, SPECT-1 and SPECT-2, are considered attractive pre-erythrocytic immune targets due to the key role they play in crossing of the malaria parasite across the dermis and the liver sinusoidal wall, prior to invasion of hepatocytes. In this study, the sequences of SPECT-1 and SPECT-2 were mapped using 20 mer-long synthetic peptides to identify high-activity binding peptides (HABPs) to HeLa cells. 17 HABPs with enzyme sensitive bindings to HeLa cells were identified: 3 predominantly α-helical in SPECT-1, and 10 α-helical and 4 β-turns/random coils in SPECT-2. Immunofluorescence assays (IFA) with antibodies raised in rabbits against chemically synthesized B-cell epitopes suggests the presence of these two proteins in the micronemes and in sporozoite membrane. ¹H NMR studies showed that HABPs located in the membrane-attack complex/perforin (MACPF) domain of SPECT-2 share high similarity with the 3D structure of C8α. Altogether, the results highlight the potential of including HABPs from SPECT-1 and SPECT-2 as components of a fully effective multistage, multiepitopic, minimal subunit-based synthetic vaccine against Plasmodium falciparum malaria.     

The effect of structural parameters and positive charge distance on the interaction free energy of antimicrobial peptides with membrane surface

Many attempts have been made to find hints explaining the relationship between physicochemical and structural properties of antimicrobial peptides (AMPs) which are relevant to their antimicrobial activities. We here found that there is a difference in the percentages of hydrophobic, hydrophilic, and charged residues between AMPs killing both bacteria and fungi (Group A) and AMPs that only kill bacteria (Group B). The percentage of charged residues in Group A AMPs is highly elevated, while in Group B the percentage of hydrophobic residues is increased. This result suggests a sequence-based mechanism of selectivity for AMPs. Moreover, we examined how the distance between basic residues affects the interaction free energy of AMPs with the membrane surface, since most of the known AMPs act by membrane perturbation. We measured the average distance between basic residues throughout the 3D structure of AMPs by defining D_pr parameter and calculated the interaction free energy for 10 AMPs that interacted with the DPPC membrane using molecular dynamics simulation. We found that the changes of the interaction free energy correlates with the change of D_pr by a linear regression coefficient of r^2?=?.47 and a cubic regression coefficient of r^2?=?.70.     

Membrane curvature and surface area per lipid affect the conformation and oligomeric state of HIV-1 fusion peptide: A combined FTIR and MD simulation study

Fourier-transformed infrared spectroscopy (FTIR) and molecular dynamics (MD) simulation results are presented to support our hypothesis that the conformation and the oligomeric state of the HIV-1 gp41 fusion domain or fusion peptide (gp41-FP) are determined by the membrane surface area per lipid (APL), which is affected by the membrane curvature. FTIR of the gp41-FP in the Aerosol-OT (AOT) reversed micellar system showed that as APL decreases from ∼ 50 to 35 Å² by varying the AOT/water ratio, the FP changes from the monomeric α-helical to the oligomeric β-sheet structure. MD simulations in POPE lipid bilayer systems showed that as the APL decreases by applying a negative surface tension, helical monomers start to unfold into turn-like structures. Furthermore, an increase in the applied lateral pressure during nonequilibrium MD simulations favored the formation of β-sheet structure. These results provide better insight into the relationship between the structures of the gp41-FP and the membrane, which is essential in understanding the membrane fusion process. The implication of the results of this work on what is the fusogenic structure of the HIV-1 FP is discussed.     

The effects of oxidised phospholipids and cholesterol on the biophysical properties of POPC bilayers

Oxidation of unsaturated membrane phospholipids by oxidative stress is associated with inflammation, infection, numerous diseases and neurodegenerative disorders. Lipid oxidation is observed in experimental samples when the parent lipid is exposed to oxidative stressors. The effect of phospholipid oxidation on the properties of biological membranes are still being explored, while low concentrations (0.1–2.0?mol%) of oxidised phospholipids are associated with disease states [1]. Previous computational studies have focused on the effect of high concentrations (~50?mol%) of oxidised phospholipids on binary lipid bilayers. This work systematically characterises the effect of lower concentrations (~10?mol%) of two oxidised lipid species, PoxnoPC (1-palmitoyl-2-(9′-oxo-nonanoyl)-sn-glycero-3-phosphocholine) or PazePC (1-palmitoyl-2-azelaoyl-sn-glycero-3-phosphocholine), on POPC/cholesterol and pure POPC bilayers. During μs atomistic simulations in pure POPC bilayers, PoxnoPC and PazePC reoriented their oxidised sn-2 acyl chains towards the solution, and PazePC adopted an extended conformation. The addition of 20?mol% cholesterol not only modulated the fluidity of the bilayers; it also modulated the flexibility of the PoxnoPC oxidised sn-2 tail, reducing bilayer disorder. In contrast, the addition of cholesterol had little effect on bilayers containing PazePC. Our studies show that the effect of oxidised lipids on the biophysical properties of a multicomponent bilayer cannot be intuitively extrapolated from a binary lipid system.     

The solution structure of ParD, the antidote of the ParDE toxin antitoxin module, provides the structural basis for DNA and toxin binding

ParD is the antidote of the plasmid-encoded toxin-antitoxin (TA) system ParD-ParE. These modules rely on differential stabilities of a highly expressed but labile antidote and a stable toxin expressed from one operon. Consequently, loss of the coding plasmid results in loss of the protective antidote and poisoning of the cell. The antidote protein usually also exhibits an autoregulatory function of the operon. In this paper, we present the solution structure of ParD. The repressor activity of ParD is mediated by the N-terminal half of the protein, which adopts a ribbon-helix-helix (RHH) fold. The C-terminal half of the protein is unstructured in the absence of its cognate binding partner ParE. Based on homology with other RHH proteins, we present a model of the ParD-DNA interaction, with the antiparallel beta-strand being inserted into the major groove of DNA. The fusion of the N-terminal DNA-binding RHH motif to the toxin-binding unstructured C-terminal domain is discussed in its evolutionary context.     

The structures of T87I phosphono-CheY and T87I/Y106W phosphono-CheY help to explain their binding affinities to the FliM and CheZ peptides

CheY is a response regulator in bacterial chemotaxis. Escherichia coli CheY mutants T87I and T87I/Y106W CheY are phosphorylatable on Asp57 but unable to generate clockwise rotation of the flagella. To understand this phenotype in terms of structure, stable analogs of the two CheY-P mutants were synthesized: T87I phosphono-CheY and T87I phosphono-CheY. Dissociation constants for peptides derived from flagellar motor protein FliM and phosphatase CheZ were determined for phosphono-CheY and the two mutants. The peptides bind phosphono-CheY almost as strongly as CheY-P; however, they do not bind T87I phosphono-CheY or T87I/Y106W phosphono-CheY, implying that the mutant proteins cannot bind FliM or CheZ tightly in vivo. The structures of T87I phosphono-CheY and T87I/Y106W phosphono-CheY were solved to resolutions of 1.8 and 2.4 Å, respectively. The increased bulk of I87 forces the side-chain of Y106 or W106, into a more solvent-accessible conformation, which occludes the peptide-binding site.     

The first laminin G-like domain of protein S is essential for binding and activation of Tyro3 receptor and intracellular signalling

The homologous proteins Gas6 and protein S (ProS1) are both natural ligands for the TAM (Tyro3, Axl, MerTK) receptor tyrosine kinases. ProS1 selectively activates Tyro3; however, the precise molecular interface of the ProS1-Tyro3 contact has not been characterised. We used a set of chimeric proteins in which each of the C-terminal laminin G-like (LG) domains of ProS1 were swapped with those of Gas6, as well as a set of ProS1 mutants with novel added glycosylations within LG1. Alongside wildtype ProS1, only the chimera containing ProS1 LG1 domain stimulated Tyro3 and Erk phosphorylation in human cancer cells, as determined by Western blot. In contrast, Gas6 and chimeras containing minimally the Gas6 LG1 domain stimulated Axl and Akt phosphorylation. We performed in silico homology modelling and molecular docking analysis to construct and evaluate structural models of both ProS1-Tyro3 and Gas6-Axl ligand-receptor interactions. These analyses revealed a contact between the ProS1 LG1 domain and the first immunoglobulin domain of Tyro3, which was similar to the Gas6-Axl interaction, and involved long-range electrostatic interactions that were further stabilised by hydrophobic and polar contacts. The mutant ProS1 proteins, which had added glycosylations within LG1 but which were all outside of the modelled contact region, all activated Tyro3 in cells with no hindrance. In conclusion, we show that the LG1 domain of ProS1 is necessary for activation of the Tyro3 receptor, involving protein-protein interaction interfaces that are homologous to those of the Gas6-Axl interaction.     

The system of galactans from Cryptonemia crenulata (Halymeniaceae, Halymeniales) and the structure of two major fractions. Kinetic studies on the alkaline cyclization of the unusual diad G2S-->D(L)6S

Cryptonemia crenulata biosynthesizes a family of dl-hybrid galactans that are based on the classical 3-linked beta-d-galactopyranosyl-->4-linked alpha-galactopyranosyl alternating sequence (A-units-->B-units). The dispersion of structures in these galactans is based on four factors, namely: (a) the amount and position of substituent groups as sulfate (major), pyruvic acid ketals, methoxyl and side substituents of beta-D-xylose and/or beta-D-galactose; (b) the ratio galactose/3,6-anhydrogalactose in the B-units; (c) the ratio D-/L-galactoses and 3,6-anhydrogalactoses also in the B-units and (d) the sequence of the diads in the linear backbone. Alkali treatment carried out on the major fraction produced a nearly quantitative formation of 3,6-anhydrogalactose units from precursor units (alpha-galactose 6-sulfate (major) and alpha-galactose 2,6-sulfate, minor). Kinetic studies show a rate constant, for the diad G2S-D(L) 6-S, of 1.7 x 10(4)s(-1) indicating a reaction faster than in lambda-carrageenans but slower than in porphyrans.     

Preparation and structural organisation of heteroleptic tetraphenylantimony(V) complexes comprising unidentately and bidentately coordinated O,O′-dialkyldithiophosphate groups: Multinuclear (C, P) CP/MAS NMR and single-crystal X-ray diffraction studies

O,O′-dipropyldithiophosphate and O,O′-di-iso-butyldithiophosphate (Dtph) tetraphenylantimony(V) complexes of the general formula [Sb(C₆H₅)₄{S₂P(OR)₂}] (R = C₃H₇, i-C₄H₉) were prepared and studied by means of ¹³C, ³¹P CP/MAS NMR spectroscopy and single-crystal X-ray diffraction. Distorted octahedral and trigonal bipyramidal molecular structures have been established for prepared complexes. These unexpected structural distinctions between chemically related compounds are defined by the principally different coordination modes of O,O′-dipropyldithiophosphate and O,O′-di-iso-butyldithiophosphate ligands in their molecular structures (i.e., S,S′-bidentate chelating and S-unidentately coordinated, respectively). To characterise quantitatively phosphorus sites in both species of dithiophosphate ligands, ³¹P chemical shift anisotropy parameters (δ_aniso and η) were calculated from spinning sideband manifolds in MAS NMR spectra. The ³¹P chemical shift tensors for the bidentate chelating and unidentately coordinated dithiophosphate ligands display a profoundly rhombic and nearly axially symmetric characters, respectively.     

The function and structural influence of selective relaxed constraint at functional intracellular loop3 of 5-HT1A serotonin-1 receptor family

Serotonin (5-HT) and its receptors have been involved in critical signal transduction mechanism and deregulation implicated in mood-related disorders. 5-HT activities are mediated through a family of transmembrane spanning serotonin receptors. Both within the family and species, 5-HT receptor protein sequence diversity and 7-transmembrane structural homogeneity have long been intriguing. In this study, we have analyzed the codon site constraint in 5-HT1 subclass receptors from 13 orthologous mammalian mRNA coding sequence. Further, the study was extended to computationally investigate the impact of non-synonymous sites with respect to function and structural significance through sequence homology algorithm and molecular dynamics simulation (MDS). Codon sites with significant posterior probability were observed in 5-HT_1A, 5‐HT_1B and 5-HT_1D receptor indicating variations in site constraint within the 5‐HT1 sub-class genes. In 5-HT_1A receptor, seven sites were detected at the functional intracellular loop₃ (ICL₃) with higher substitution rate through Codeml program. Sequence homology algorithm identifies that these sites were functionally tolerant within the mammals representing a selectively relaxed constraint at this domain. On the other hand, the root mean square deviation (rmsd) values from MDS suggest differences in structural conformation of ICL₃ models among the species. Specifically, the human ICL₃ model fluctuation was comparatively more stable than other species. Hence, we argue that these sites may have varying influence in G-proteins coupling and activation of effectors systems through downstream interacting accessory proteins of cell among the species. However, further experimental studies are required to elucidate the precise role and the seeming difference of these sites in 5-HT receptors between species.