Modeling of non-covalent complexes of the cell-penetrating peptide CADY and its siRNA cargo.

CADY is a cell-penetrating peptide spontaneously making non-covalent complexes with Short interfering RNAs (siRNAs) in water. Neither the structure of CADY nor that of the complexes is resolved. We have calculated and analyzed 3D models of CADY and of the non-covalent CADY-siRNA complexes in order to understand their formation and stabilization. Data from the ab initio calculations and molecular dynamics support that, in agreement with the experimental data, CADY is a polymorphic peptide partly helical. Taking into consideration the polymorphism of CADY, we calculated and compared several complexes with peptide/siRNA ratios of up to 40. Four complexes were run by using molecular dynamics. The initial binding of CADYs is essentially due to the electrostatic interactions of the arginines with siRNA phosphates. Due to a repetitive arginine motif (XLWR(K)) in CADY and to the numerous phosphate moieties in the siRNA, CADYs can adopt multiple positions at the siRNA surface leading to numerous possibilities of complexes. Nevertheless, several complex properties are common: an average of 14±1 CADYs is required to saturate a siRNA as compared to the 12±2 CADYs experimentally described. The 40 CADYs/siRNA that is the optimal ratio for vector stability always corresponds to two layers of CADYs per siRNA. When siRNA is covered by the first layer of CADYs, the peptides still bind despite the electrostatic repulsion. The peptide cage is stabilized by hydrophobic CADY-CADY contacts thanks to CADY polymorphism. The analysis demonstrates that the hydrophobicity, the presence of several positive charges and the disorder of CADY are mandatory to make stable the CADY-siRNA complexes.

Characterization of SLC26A9 in patients with CF-like lung disease.

Diffuse bronchiectasis is a common problem in respiratory clinics. We hypothesized that mutations in the solute carrier 26A9 (SLC26A9) gene, encoding for a chloride (Cl(-)) transporter mainly expressed in lungs, may lead to defects in mucociliary clearance. We describe two missense variants in the SLC26A9 gene in heterozygote patients presenting with diffuse idiopathic bronchiectasis : p.Arg575Trp, identified in a patient also heterozygote for p.Phe508del in the CFTR gene; and p.Val486Ile. Expression of both mutants in Xenopus laevis oocytes abolished SLC26A9-mediated Cl(-) conductance without decreasing protein membrane expression. Coexpression of CFTR with SLC26A9-p.Val486Ile resulted in a significant increase in the Cl(-) current induced by PKA stimulation, similar to that obtained in oocytes expressing CFTR and SLC26A9-WT. In contrast, coexpression of CFTR with SLC26A9-p.Arg575Trp inhibited SLC26A9-enhanced CFTR activation upon PKA. Further structure-function analyses led us to propose a site encompassing Arg575 in the SLC26A9-STAS domain for CFTR-SLC26A9 interaction. We hypothesize that SLC26A9-p.Arg575Trp prevented SLC26A9-mediated functional activation of CFTR by altering SLC26A9-CFTR interaction. Although we cannot confirm that these mutations by themselves are deleterious, we propose that they trigger the pathogenic role of a single CFTR mutation and provide insight into a novel mechanism of Cl(-) transport alteration across the respiratory mucosa, based on functional inhibition of CFTR.

In silico predictions of 3D structures of linear and cyclic peptides with natural and non-proteinogenic residues.

We extended the use of Peplook, an in silico procedure for the prediction of three-dimensional (3D) models of linear peptides to the prediction of 3D models of cyclic peptides and thanks to the ab initio calculation procedure, to the calculation of peptides with non-proteinogenic amino acids. Indeed, such peptides cannot be predicted by homology or threading. We compare the calculated models with NMR and X-ray models and for the cyclic peptides, with models predicted by other in silico procedures (Pep-Fold and I-Tasser). For cyclic peptides, on a set of 38 peptides, average root mean square deviation of backbone atoms (BB-RMSD) was 3.8 and 4.1 Å for Peplook and Pep-Fold, respectively. The best results are obtained with I-Tasser (2.5 Å) although evaluations were biased by the fact that the resolved Protein Data Bank models could be used as template by the server. Peplook and Pep-Fold give similar results, better for short (up to 20 residues) than for longer peptides. For peptides with non-proteinogenic residues, performances of Peplook are sound with an average BB-RMSD of 3.6 Å for 'non-natural peptides' and 3.4 Å for peptides combining non-proteinogenic residues and cyclic structure. These results open interesting possibilities for the design of peptidic drugs.

The role of proline in the membrane re-entrant helix of caveolin-1.

Caveolin-1 has a segment of hydrophobic amino acids comprising approximately residues 103-122. We have performed an in silico analysis of the conformational preference of this segment of caveolin-1 using PepLook. We find that there is one main group of stable conformations corresponding to a hydrophobic U bent model that would not traverse the membrane. Furthermore, the calculations predict that substituting the Pro(110) residue with an Ala will change the conformation to a straight hydrophobic helix that would traverse the membrane. We have expressed the P110A mutant of caveolin-1, with a FLAG tag at the N terminus, in HEK 293 cells. We evaluate the topology of the proteins with confocal immunofluorescence microscopy in these cells. We find that FLAG tag at the N terminus of the wild type caveolin-1 is not reactive with antibodies unless the cell membrane is permeabilized with detergent. This indicates that in these cells, the hydrophobic segment of this protein is not transmembrane but takes up a bent conformation, making the protein monotopic. In contrast, the FLAG tag at the N terminus of the P110A mutant is equally exposed to antibodies, before and after membrane permeabilization. We also find that the P110A mutation causes a large reduction of endocytosis of caveolae, cellular lipid accumulation, and lipid droplet formulation. In addition, we find that this mutation markedly reduces the ability of caveolin-1 to form structures with the characteristic morphology of caveolae or to partition into the detergent-resistant membranes of these cells. Thus, the single Pro residue in the membrane-inserting segment of caveolin-1 plays an important role in both the membrane topology and localization of the protein as well as its functions.

Standardized evaluation of protein stability.

We compare mean force potential values of a large series of PDB models of proteins and peptides and find that, either as monomers or polymers, proteins longer than 200-250 residues have equivalent MFP values that are averaged to -65+/-3 kcal/aa. This value is named the standard or stability value. The standard value is reached irrespective of sequences and 3D folds. Peptides are too short to follow the rule and frequently exist as populations of conformers; one exception is peptides in amyloid fibrils. Fibrils surpass the standard value in accordance with their uppermost stability. In parallel, we calculate median MFP values of amino acids in stably folded PDB models of proteins: median values vary from -25 for Gly to -115 kcal/aa for Trp. These median values are used to score primary sequences of proteins: all sequences converge to a mean value of -63.5+/-2.5 kcal/aa, i.e., only 1.5 kcal less than the folded model standard. Sequences from unfolded proteins have lower values. This supports the conclusion that sequences carry in an important message and more specifically that diversity of amino acids in sequences is mandatory for stability. We also use the median amino acid MFP to score residue stability in 3D folds. This demonstrates that 3D folds are compromises between fragments of high and fragments of low scores and that functional residues are often but not always in the extreme score values. The approach opens to possibilities of evaluating any 3D model and of detecting functional residues and should help in conducting mutation assays.

Realistic modeling approaches of structure-function properties of CPPs in non-covalent complexes.

Transfers of cargoes into cells by means of carrier peptides are multi-steps biological phenomenon the mechanisms of which are unclear. We here discuss bases of realistic in silico molecular modeling approaches of the formation of non-covalent complexes considering CPPs and cargo diversities.

Cholesterol interaction with proteins that partition into membrane domains: an overview.

Biological membranes are complex structures composed largely of proteins and lipids. These components have very different structural and physical properties and consequently they do not form a single homogeneous mixture. Rather components of the mixture are more enriched in some regions than in others. This can be demonstrated with simple lipid mixtures that spontaneously segregate components so as to form different lipid phases that are immiscible with one another. The segregation of molecular components of biological membranes also involves proteins. One driving force that would promote the segregation of membrane components is the preferential interaction between a protein and certain lipid components. Among the varied lipid components of mammalian membranes, the structure and physical properties of cholesterol is quite different from that of other major membrane lipids. It would therefore be expected that in many cases proteins would have very different energies of interaction with cholesterol vs. those of other membrane lipids. This would be sufficient to cause segregation of components in membranes. The factors that facilitate the interaction of proteins with cholesterol are varied and are not yet completely understood. However, there are certain groups that are present in some proteins that facilitate interaction of the protein with cholesterol. These groups include saturated acyl chains of lipidated proteins, as well as certain amino acid sequences. Although there is some understanding as to why these particular groups favour interaction with cholesterol, our knowledge of these molecular features is not sufficiently developed to allow for the design of agents that will modify such binding.

Insight into the cellular uptake mechanism of a secondary amphipathic cell-penetrating peptide for siRNA delivery.

Delivery of siRNA remains a major limitation to their clinical application, and several technologies have been proposed to improve their cellular uptake. We recently described a peptide-based nanoparticle system for efficient delivery of siRNA into primary cell lines: CADY. CADY is a secondary amphipathic peptide that forms stable complexes with siRNA and improves their cellular uptake independently of the endosomal pathway. In the present work, we have combined molecular modeling, spectroscopy, and membrane interaction approaches in order to gain further insight into CADY/siRNA particle mechanism of interaction with biological membrane. We demonstrate that CADY forms stable complexes with siRNA and binds phospholipids tightly, mainly through electrostatic interactions. Binding to siRNA or phospholipids triggers a conformational transition of CADY from an unfolded state to an alpha-helical structure, thereby stabilizing CADY/siRNA complexes and improving their interactions with cell membranes. Therefore, we propose that CADY cellular membrane interaction is driven by its structural polymorphism which enables stabilization of both electrostatic and hydrophobic contacts with surface membrane proteoglycan and phospholipids.

Structural polymorphism of two CPP: an important parameter of activity.

Despite numerous investigations, the important structural features of Cell Penetrating Peptides (CPPs) remain unclear as demonstrated by the difficulties encountered in designing new molecules. In this study, we focused our interest on Penetratin and Transportan and several of their variants. Penetratin W48F and Penetratin W48F/W56F exhibit a reduced and a complete lack of cellular uptake, respectively; TP07 and TP10 present a similar cellular uptake as Transportan and TP08, TP13 and TP15 display no or weak internalization capacity. We applied the algorithmic method named PepLook to analyze the peptide polymorphism. The study reveals common conformational characteristics for the CPPs and their permeable variants: they all are polymorphic. Negative, non permeable, mutants share the opposite feature since they are monomorphic. Finally, we support the hypothesis that structural polymorphism may be crucial since it provides peptides with the possibility of adapting their conformation to medium hydrophobicity and or to partner diversity.

Large changes in the CRAC segment of gp41 of HIV do not destroy fusion activity if the segment interacts with cholesterol.

The membrane-proximal external region (MPER) of the gp41 fusion protein of HIV is highly conserved among isolates of this virus and is considered a target for vaccine development. This region also appears to play a role in membrane fusion as well as localization of the virus to cholesterol-rich domains in membranes. The carboxyl terminus of MPER has the sequence LWYIK and appears to have an important role in cholesterol interactions. We have tested how amino acid substitutions that would affect the conformational flexibility of this segment could alter its interaction with cholesterol. We studied a family of peptides (all peptides as N-acetyl-peptide amides) with P, G, or A substituting for W and I of the LWYIK sequence. The peptide having the greatest effect on cholesterol distribution in membranes was the most flexible one, LGYGK. The corresponding mutation in gp41 resulted in a protein retaining 72% of the fusion activity of the wild-type protein. Two other peptides were synthesized, also containing two Gly residues, GWGIK and LWGIG, and did not have the ability to sequester cholesterol as efficiently as LGYGK did. Making the corresponding mutants of gp41 showed that these other two double Gly substitutions resulted in proteins that were much less fusogenic, although they were equally well expressed at the cell surface. The study demonstrates that drastic changes can be made in the LWYIK segment with the retention of a significant fraction of the fusogenic activity, as long as the mutant proteins interact with cholesterol.

Computational study of colipase interaction with lipid droplets and bile salt micelles.

Colipase is a key element in the lipase-catalyzed hydrolysis of dietary lipids. Although devoid of enzymatic activity, colipase promotes the pancreatic lipase activity in physiological intestinal conditions by anchoring the enzyme at the surface of lipid droplets. Analysis of structures of NMR colipase models and simulations of their interactions with various lipid aggregates, lipid droplet, and bile salt micelle, were carried out to determine and to map the lipid binding sites on colipase. We show that the micelle and the oil droplet bind to the same side of colipase 3D structure, mainly the hydrophobic fingers. Moreover, it appears that, although colipase has a single direction of interaction with a lipid interface, it does not bind in a specific way but rather oscillates between different positions. Indeed, different NMR models of colipase insert different fragments of sequence in the interface, either simultaneously or independently. This supports the idea that colipase finger plasticity may be crucial to adapt the lipase activity to different lipid aggregates.

The N-terminal juxtamembranous domain of KCNQ1 is critical for channel surface expression: implications in the Romano-Ward LQT1 syndrome.

N-terminal mutations in the KCNQ1 channel are frequently linked to fatal arrhythmias in newborn children and adolescents but the cellular mechanisms involved in this dramatic issue remain, however, to be discovered. Here, we analyzed the trafficking of a series of N-terminal truncation mutants and identified a critical trafficking motif of KCNQ1. This determinant is located in the juxtamembranous region preceding the first transmembrane domain of the protein. Three mutations (Y111C, L114P and P117L) implicated in inherited Romano-Ward LQT1 syndrome, are embedded within this domain. Reexpression studies in both COS-7 cells and cardiomyocytes showed that the mutant proteins fail to exit the endoplasmic reticulum. KCNQ1 subunits harboring Y111C or L114P exert a dominant negative effect on the wild-type KCNQ1 subunit by preventing plasma membrane trafficking of heteromultimeric channels. The P117L mutation had a less pronounced effect on the trafficking of heteromultimeric channels but altered the kinetics of the current. Furthermore, we showed that the trafficking determinant in KCNQ1 is structurally and functionally conserved in other KCNQ channels and constitutes a critical trafficking determinant of the KCNQ channel family. Computed structural predictions correlated the potential structural changes introduced by the mutations with impaired protein trafficking. In conclusion, our studies unveiled a new role of the N-terminus of KCNQ channels in their trafficking and its implication in severe forms of LQT1 syndrome.

Enrichment of anhydrous milk fat in polyunsaturated fatty acid residues from linseed and rapeseed oils through enzymatic interesterification.

Lipozyme TL IM was used in a solvent-free batch and microaqueous system for enzymatic interesterification of anhydrous milkfat (AMF) with linseed oil (LO) in binary blends and with rapeseed oil (RO) in one ternary blend. The aim was to obtain and characterize physicochemically fats enriched with unsaturated C 18 fatty acids (oleic, linoleic, and, especially, linolenic acids) from natural vegetable oils. Binary blends of AMF/LO 100/0, 90/10, 80/20, 70/30, and 60/40 (w/w) were interesterified. The change in triacylglycerol (TAG) profiles showed that quasi-equilibrium conditions were reached after 4-6 h of reaction. Free fatty acid contents <1%. The decrease in solid fat content and in dropping point temperature obtained with increasing content of LO and interesterification resulted in good plastic properties for the products originating from the blends 70/30 and 60/40. This was confirmed by textural measurements. Melting profiles determined by differential scanning calorimetry showed complete disappearance of low-melting TAGs from LO and the formation of intermediary species with a lower melting temperature. Oxidative stability of the interesterified products was diminished with increasing LO content, resulting in low oxidation induction times. A ternary blend composed of AMF/RO/LO 70/20/10 gave satisfactory rheological and oxidative properties, fulfilling the requirements for a marketable spread and, moreover, offering increased potential health benefits due to the enriched content in polyunsaturated fatty acid residues.

Rational design of complementary peptides to the betaAmyloid 29-42 fusion peptide: an application of PepDesign.

Peptides in solution currently exist under several conformations; an equilibrium which varies with solvent polarity. Despite or because of this structure versatility, peptides can be selective biological tools: they can adapt to a target, vary conformation with solvents and so on. These capacities are crucial for cargo carriers. One promising way of using peptides in biotechnologies is to decipher their medium-sequence-structure-function relationships and one approach is molecular modelling. Only few "in silico" methods of peptide design are described in the literature. Most are used in support of experimental screening of peptide libraries. However, the way they are made does not teach us much for future researches. In this paper, we describe an "in silico" method (PepDesign) which starts by analysing the native interaction of a peptide with a target molecule in order to define which points are important. From there, a modelling protocol for the design of 'better' peptides is set. The PepDesign procedure calculates new peptides fulfilling the hypothesis, tests the conformational space of these peptides in interaction with the target by angular dynamics and goes up to the selection of the best peptide based on the analysis of complex structure properties. Experimental biological assays are finally used to test the selected peptides, hence to validate the approach. Applications of PepDesign are wide because the procedure will remain similar irrespective of the target which can be a protein, a drug or a nucleic acid. In this paper, we describe the design of peptides which binds to the fusogenic helical form of the C-terminal domain of the Abeta peptide (Abeta29-42).

Prediction of peptide structure: how far are we?

Rational design of peptides is a challenge, which would benefit from a better knowledge of the rules of sequence-structure-function relationships. Peptide structures can be approached by spectroscopy and NMR techniques but data from these approaches too frequently diverge. Structures can also be calculated in silico from primary sequence information using three algorithms: Pepstr, Robetta, and PepLook. The most recent algorithm, PepLook introduces indexes for evaluating structural polymorphism and stability. For peptides with converging experimental data, calculated structures from PepLook and, to a lesser extent from Pepstr, are close to NMR models. The PepLook index for polymorphism is low and the index for stability points out possible binding sites. For peptides with divergent experimental data, calculated and NMR structures can be similar or, can be different. These differences are apparently due to polymorphism and to different conditions of structure assays and calculations. The PepLook index for polymorphism maps the fragments encoding disorder. This should provide new means for the rational design of peptides.

Tilted peptides: the history.

Nature has selected peptide motifs for protein functions. It is clear that specific sequence motifs can identify families of enzymes. These sequence motifs are one dimensional signatures and nature has also developed two dimension motifs which cannot be read in the one dimension of sequence language but can be detected in the three dimensional properties of a secondary structure. One of such motifs is tilted peptides. They do not correspond to any consensus of sequence but correspond to a consensus motif where hydrophobicity balance is used as a functional device. In the nineteen eighties, the first tilted peptide was deciphered from the sequence of a virus fusion protein by molecular modelling. It was described as a protein fragment hydrophobic enough to insert into a membrane but too short to span it. The fragment exhibited an asymmetric distribution of hydrophobicity along the helix long axis and hence, was unable to lie parallel or perpendicular to a membrane surface but adopted an orientation in between. Hydrophobicity motif was a very new and very challenging concept and tilted peptides were rapidly found to be involved in several mechanisms of virus fusion. They were also found to be involved in protein secretion and future studies could establish their involvement in the destabilization of 3D protein structures and in the alpha to beta transconformations, which drive the generation of amyloid deposits.

In vitro characterization of the homogalacturonan-binding domain of the wall-associated kinase WAK1 using site-directed mutagenesis.

Wall-associated kinase 1--WAK1 is a transmembrane protein containing a cytoplasmic Ser/Thr kinase domain and an extracellular domain in contact with the pectin fraction of the plant cell wall in Arabidopsis thaliana (L.) HEYNH. In a previous paper [Decreux, A., Messiaen, J., 2005. Wall-associated kinase WAK1 interacts with cell wall pectins in a calcium-induced conformation. Plant Cell Physiol. 46, 268-278], we showed that a recombinant peptide expressed in yeast corresponding to amino acids 67-254 of the extracellular domain of WAK1 specifically interacts with commercial non-methylesterified homogalacturonic acid, purified homogalacturonans from Arabidopsis and oligogalacturonides in a calcium-induced conformation. In this report, we used a receptor binding domain sequence-based prediction method to identify four putative binding sites in the extracellular domain of WAK1, in which cationic amino acids were selected for substitution by site-directed mutagenesis. Interaction studies between mutated forms of WAK1 and homogalacturonans allowed us to identify and confirm at least five specific amino acids involved in the interaction with homogalacturonan dimers and multimers. The presence of this homogalacturonan-binding domain within the extracellular domain of WAK1 is discussed in terms of cell wall architecture and signal transduction.

Fusogenic Alzheimer's peptide fragment Abeta (29-42) in interaction with lipid bilayers: secondary structure, dynamics, and specific interaction with phosphatidyl ethanolamine polar heads as revealed by solid-state NMR.

The interaction of the native Alzheimer's peptide C-terminal fragment Abeta (29-42), and two mutants (G33A and G37A) with neutral lipid bilayers made of POPC and POPE in a 9:1 molar ratio was investigated by solid-state NMR. This fragment and the lipid composition were selected because they represent the minimum requirement for the fusogenic activity of the Alzheimer's peptide. The chemical shifts of alanine methyl isotropic carbon were determined by MAS NMR, and they clearly demonstrated that the major form of the peptide equilibrated in membrane is not in a helical conformation. (2)H NMR, performed with acyl chain deuterated POPC, demonstrated that there is no perturbation of the acyl chain's dynamics and of the lipid phase transition temperature. (2)H NMR, performed with alanine methyl-deuterated peptide demonstrated that the peptide itself has a limited mobility below and above the lipid phase transition temperature (molecular order parameter equal to 0.94). MAS (31)P NMR revealed a specific interaction with POPE polar head as seen by the enhancement of POPE phosphorus nuclei T(2) relaxation. All these results are in favor of a beta-sheet oligomeric association of the peptide at the bilayer interface, preferentially recruiting phosphatidyl ethanolamine polar heads.

Aromatic side-chain interactions in proteins. Near- and far-sequence His-X pairs.

Several studies have analysed aromatic interactions, involving mostly phenylalanine, tyrosine and tryptophan. Only a few studies have considered histidine as an interacting aromatic residue. An extensive analysis of aromatic His-X interactions is performed here on a data set of 593 PDB structures: 68% of the histidine are involved in aromatic pairs and 1271 non-redundant His-X pairs were analysed. Thirty percent of these pairs involve an aromatic partner less than 6 apart in the sequence. These near-sequence pairs correspond to conformations which stabilise secondary structures, mainly alpha-helices when the residues are 4 apart and beta-strands when they are 2 apart in the sequence. The partners of the other His-X pairs (887, 70%) are more than 5 apart in the sequence. Of these far-sequence pairs, 35% bridge beta strands and only 9% helices. The near-sequence pairs are sterically constrained as supported by conformer distribution. The X partners of far-sequence His-X pairs are mainly "above" the histidine ring with tilted and normal rings, corresponding to a "T shape; face to edge" orientation. Phenylalanine, the only aromatic residue with no heteroatom, is a disfavoured partner, whereas histidine is the preferred one. Heteroatom-heteroatom interactions are favoured in near-sequence as well as in far-sequence His-His, His-Trp and His-Tyr pairs.