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.     

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.     

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.     

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 structural analysis of protein–protein interactions by NMR spectroscopy

A comprehensive understanding of protein–protein interactions is an important next step in our quest to understand how the information contained in a genome is put into action. Although a number of experimental techniques can report on the existence of a protein– protein interaction, very few can provide detailed structural information. NMR spectroscopy is one of these, and in recent years several complementary NMR approaches, including residual dipolar couplings and the use of paramagnetic effects, have been developed that can provide insight into the structure of protein–protein complexes. In this article, we review these approaches and comment on their strengths and weaknesses.     

βVI Turns in peptides and proteins: A model peptide mimicry

To investigate the role of proline in defining β turn conformations within cyclic hexa- and pentapeptides we synthesized and determined the conformations of a series of L - and D -proline-containing peptides by means of 2D NMR spectroscopy and restrained molecular dynamics simulations. Due to cis/trans isomerism the L -proline peptides adopt at least two different conformations that are analyzed and compared to the structures of the corresponding D -proline peptides. The cis conformations of the compounds cyclo(-Pro-Ala-Ala-Pro-Ala-Ala-), cyclo(-Arg-Gly-Asp-Phe-Pro-Gly-), cyclo(-Arg-Gly-Asp-Phe-Pro-Ala-), cyclo(-Pro-Ala-Ala-Ala-Ala--), and cyclo(-Pro-Ala-Pro-Ala-Ala-) form uncommon βVI turns that mimic the turn geometries found in crystallographically refined protein structures at such a detailed level that even preferred side chain orientations are reproduced. The ratios of the cis/trans isomers are analyzed in terms of the steric demand of the proline-following residue. The conformational details derived from this study illustrate the importance of the examination of small model compounds derived from protein loop regions, especially if bioactive recognition sequences, such as RGD (Arg-Gly-Asp), are incorporated. © 1993 Wiley-Liss, Inc.     

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.     

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.     

Ligand binding and homology modelling of insect odorant‐binding proteins

This review describes the main characteristics of odorant‐binding proteins (OBPs) for homology modelling and presents a summary of structure prediction studies on insect OBPs, along with the steps involved and some limitations and improvements. The technique involves a computing approach to model protein structures and is based on a comparison between a target (unknown structure) and one or more templates (experimentally determined structures). As targets for structure prediction, OBPs are considered to play a functional role for recognition, desorption, scavenging, protection and transportation of hydrophobic molecules (odourants) across an aqueous environment (lymph) to olfactory receptor neurones (ORNs) located in sensilla, the main olfactory units of insect antennae. Lepidopteran pheromone‐binding proteins, a subgroup of OBPs, are characterized by remarkable structural features, in which high sequence identities (approximately 30%) among these OBPs and a large number of available templates can facilitate the prediction of precise homology models. Approximately 30 studies have been performed on insect OBPs using homology modelling as a tool to predict their structures. Although some of the studies have assessed ligand‐binding affinity using structural information and biochemical measurements, few have performed docking and molecular dynamic (MD) simulations as a virtual method to predict best ligands. Docking and MD simulations are discussed in the context of discovery of novel semiochemicals (super‐ligands) using homology modelling to conceive further strategies in insect management.     

Analyzing large‐scale structural change in proteins: Comparison of principal component projection and sammon mapping

Effective analysis of large-scale conformational transitions in macromolecules requires transforming them into a lower dimensional representation that captures the dominant motions. Herein, we apply and compare two different dimensionality reduction techniques, namely, principal component analysis (PCA), a linear method, and Sammon mapping, which is nonlinear. The two methods are used to analyze four different protein transition pathways of varying complexity, obtained by using either the conjugate peak refinement method or constrained molecular dynamics. For the return-stroke in myosin, both Sammon mapping and PCA show that the conformational change is dominated by a simple rotation of a rigid body. Also, in the case of the T-->R transition in hemoglobin, both methods are able to identify the two main quaternary transition events. In contrast, in the cases of the unfolding transition of staphylococcal nuclease or the signaling switch of Ras p21, which are both more complex conformational transitions, only Sammon mapping is able to identify the distinct phases of motion.     

Probing the binding mechanism of capecitabine to human serum albumin using spectrometric methods,molecular modeling,and chemometrics approach

Capecitabine as a prodrug of 5-Fluorouracil plays an important role in the treatment of breast and gastrointestinal cancers. Herein, in view of the importance of this drug in chemotherapy, interaction mechanism between Capecitabine (CAP) and human serum albumin (HSA) as a major transport protein in the blood circulatory system has been investigated by using a combination of spectroscopic and molecular modeling methods. The fluorescence spectroscopic results revealed that capecitabine could effectively quench the intrinsic fluorescence of HSA through a static quenching mechanism. Evaluation of the thermodynamic parameters suggested that the binding process was spontaneous while hydrogen bonds and van der Waals forces played a major role in this interaction. The value of the binding constant (K_b = 1.820 × 10⁴) suggested a moderate binding affinity between CAP and HSA which implies its easy diffusion from the circulatory system to the target tissue. The efficiency of energy transfer and the binding distance between the donor (HSA) and acceptor (CAP) were determined according to forster theory of nonradiation energy transfer as 0.410 and 4.135 nm, respectively. Furthermore, UV–Vis spectroscopic results confirmed that the interaction was occurred between HSA and CAP and caused conformational and micro-environmental changes of HSA during the interaction. Multivariate curve resolution-alternating least square (MCR-ALS) methodology as an efficient chemometric tool was used to separate the overlapped spectra of the species. The MCR-ALS result was exploited to estimate the stoichiometry of interaction and to provide concentration and structural information about HSA-CAP interactions. Molecular docking studies suggested that CAP binds mainly to the subdomain IIA of HSA, which were compatible with those obtained by experimental data. Finally, molecular dynamics simulation (MD) was performed on the best docked complex by considering the permanence and flexibility of HSA-CAP complex in the binding site. MD result showed that CAP could steadily bind to HSA in the site I based on the formation of hydrogen bond and π-π stacking interaction in addition to hydrophobic force.     

The application of DOSY NMR and molecular dynamics simulations to explore the mechanism(s) of micelle binding of antimicrobial peptides containing unnatural amino acids

Anionic and zwitterionic micelles are often used as simple models for the lipids found in bacterial and mammalian cell membranes to investigate antimicrobial peptide‐lipid interactions. In our laboratory we have employed a variety of 1D, 2D, and diffusion ordered (DOSY) NMR experiments to investigate the interactions of antimicrobial peptides containing unnatural amino acids with SDS and DPC micelles. Complete assignment of the proton spectra of these peptides is prohibited by the incorporation of a high percentage of unnatural amino acids which don't contain amide protons into the backbone. However preliminary assignment of the TOCSY spectra of compound 23 in the presence of both micelles indicated multiple conformers are present as a result of binding to these micelles. Chemical Shift Indexing agreed with previously collected CD spectra that indicated on binding to SDS micelles compound 23 adopts a mixture of α‐helical structures and on binding to DPC micelles this peptide adopts a mixture of helical and β‐turn/sheet like structures. DOSY NMR experiments also indicated that the total positive charge and the relative placement of that charge at the N‐terminus or C‐terminus are important in determining the mole fraction of the peptide that will bind to the different micelles. DOSY and ¹H‐NMR experiments indicated that the length of Spacer #1 plays a major role in defining the binding conformation of these analogs with SDS micelles. Results obtained from molecular simulations studies of the binding of compounds 23 and 36 with SDS micelles were consistent with the observed NMR results. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 548–561, 2013.     

Molecular dynamics simulation of the induced‐fit binding process of DNA aptamer and L‐argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced‐fit binding mechanism for L ‐argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer‐Arm binding and the induced‐fit binding process at atomic resolution. The simulation results revealed that the Watson‐Crick pair (G10‐C16), C9, A12, and C17 bases play important roles in aptamer‐Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced‐fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen‐bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer‐ligand induced‐fit binding mechanism. The critical bases in aptamer‐ligand binding may provide a guideline in aptamer design for molecular recognition engineering.     

The interplay between the clustering of transmembrane proteins and coupling of anchored membrane proteins

Clustering of membrane proteins plays an important role in many cellular activities such as protein sorting and signal transduction. In this study, we used dissipative particle dynamics simulation method to investigate the clustering of anchored membrane proteins (AMPs) in the presence of transmembrane proteins (TMPs). First, our simulation results show that clustering of AMPs and that of TMPs are in fact interdependent, and depending on their hydrophobic length, both protein mixing and protein demixing are observed. Especially, the protein demixing occurs only when the hydrophobic mismatch of TMPs is negative while that of AMPs is positive. Second, our simulation results indicate that the clustering of TMPs also modulates the coupling of the clustering of AMPs in both leaflets. On the one hand, the coupling between AMPs in different leaflets will be strongly restrained if TMPs form protein mixing with AMPs in one leaflet and protein demixing with AMPs in the other leaflet. On the other hand, the coupling between AMPs can be enhanced or mediated by TMPs when TMPs mix with AMPs in both leaflets. Our results may have some implications on our understanding of how different types of membrane proteins cluster and provide a possible explanation of how TMPs participate in signal transduction across cellular membranes.     

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.