PubMed Health⌕ Search

Biomedical subjects

Annick Thomas

Publications and source records attributed to Annick Thomas.

At least 19 recordsLinked to original sources

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.

Algorithms↗

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.

Adult↗

Juxtamembrane protein segments that contribute to recruitment of cholesterol into domains.

We investigated the properties of several peptides with sequences related to LWYIK, a segment found in the gp41 protein of HIV and believed to play a role in sequestering this protein to a cholesterol-rich domain in the membrane. This segment fulfills the requirements to be classified as a CRAC motif that has been suggested to predict those proteins that will partition into cholesterol-rich regions of the membrane. All of the peptides were studied with the terminal amino and carboxyl groups blocked, i.e., as N-acetyl-peptide-amides. Effects of cholesterol on the intensity of W emission generally parallel DSC evidence of sequestration of cholesterol. Modeling studies indicate that all of these peptides tend to partition with their mass center at the membrane interface at the level of the hydroxyl of cholesterol. Interaction with cholesterol is dual: van der Waals interactions between mainly hydrophobic surfaces and electrostatic stabilization of the cholesterol OH group. Thus, both experiments and modeling studies indicate that the preference of CRAC motifs for cholesterol-rich domains might be related to a membrane interfacial preference of the motif, to a capacity to wrap and block the cholesterol polar OH group by H-bond interactions, and to a capacity for peptide aromatic side chains to stack with cholesterol. These results were supported by studies of single mutations in the gp41 protein of HIV-1, in which L(679) is replaced with I. Despite the similarity of the properties of these amino acid residues, this single substitution resulted in a marked attenuation of the ability of JC53-BL HeLa-based HIV-1 indicator cells to form syncytia.

Amino Acid Sequence↗

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.

Arabidopsis Proteins↗

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.

Amino Acid Sequence↗

Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure.

We analyzed structural features of 11,038 direct atomic contacts (either electrostatic, H-bonds, hydrophobic, or other van der Waals interactions) extracted from 139 protein-DNA and 49 protein-RNA nonhomologous complexes from the Protein Data Bank (PDB). Globally, H-bonds are the most frequent interactions (approximately 50%), followed by van der Waals, hydrophobic, and electrostatic interactions. From the protein viewpoint, hydrophilic amino acids are over-represented in the interaction databases: Positively charged amino acids mainly contact nucleic acid phosphate groups but can also interact with base edges. From the nucleotide point of view, DNA and RNA behave differently: Most protein-DNA interactions involve phosphate atoms, while protein-RNA interactions involve more frequently base edge and ribose atoms. The increased participation of DNA phosphate involves H-bonds rather than salt bridges. A statistical analysis was performed to find the occurrence of amino acid-nucleotide pairs most different from chance. These pairs were analyzed individually. Finally, we studied the conformation of DNA in the interaction sites. Despite the prevalence of B-DNA in the database, our results suggest that A-DNA is favored in the interaction sites.

Amino Acids↗

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).

Amino Acid Sequence↗

Role of the lid hydrophobicity pattern in pancreatic lipase activity.

Pancreatic lipase is a soluble globular protein that must undergo structural modifications before it can hydrolyze oil droplets coated with bile salts. The binding of colipase and movement of the lipase lid open access to the active site. Mechanisms triggering lid mobility are unclear. The *KNILSQIVDIDGI* fragment of the lid of the human pancreatic lipase is predicted by molecular modeling to be a tilted peptide. Tilted peptides are hydrophobicity motifs involved in membrane fusion and more globally in perturbations of hydrophobic/hydrophilic interfaces. Analysis of this lid fragment predicts no clear consensus of secondary structure that suggests that its structure is not strongly sequence determined and could vary with environment. Point mutations were designed to modify the hydrophobicity profile of the [240-252] fragment and their consequences on the lipase-mediated catalysis were tested. Two mutants, in which the tilted peptide motif was lost, also have poor activity on bile salt-coated oil droplets and cannot be reactivated by colipase. Conversely, one mutant in which a different tilted peptide is created retains colipase dependence. These results suggest that the tilted hydrophobicity pattern of the [240-252] fragment is neither important for colipase binding to lipase, nor for interfacial binding but is important to trigger the maximal catalytic efficiency of lipase in the presence of bile salt.

Algorithms↗

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.

Amino Acid Sequence↗

Partial atomic charges of amino acids in proteins.

Using a semiempirical quantum mechanical procedure (FCPAC) we have calculated the partial atomic charges of amino acids from 494 high-resolution protein structures. To analyze the influence of the protein's environment, we considered each residue under two conditions: either as the center of a tripeptide with PDB structure geometry (free) or as the center of 13-16 amino acid clusters extracted from the PDB structure (buried). The partial atomic charges from residues in helices and in sheets were separated. The FCPAC partial atomic charges of the Cbeta and Calpha of most residues correlate with their helix propensity, positively for Cbeta and negatively for Calpha (r2 = 0.76 and 0.6, respectively). The main consequence of burying residues in proteins is the polarization of the backbone C=O bond, which is more pronounced in helices than in sheets. The average shift of the oxygen partial charges that results from burying is -0.120 in helix and -0.084 in sheet with the charge of the proton as unit. Linear correlations are found between the average NMR chemical shifts and the average FCPAC partial charges of Calpha (r2 = 0.8-0.85), N (r3 = 0.67-0.72), and Cbeta (r2 = 0.62) atoms. Correlations for helix and beta-sheet FCPAC partial charges show parallel regressions, suggesting that the charge variations due to burying in proteins differentiate between the dihedral angle effects and the polarization of backbone atoms.

Amino Acids↗

Aromatic side-chain interactions in proteins: near- and far-sequence Tyr-X pairs.

In the present study, an extensive analysis of the aromatic Tyr-X interactions is performed on a data set of 593 PDB structures, X being Phe, His, Tyr, and Trp. The nonredundant Tyr-X pairs (2645) were retained and separated by both the residue distance in the sequence and the secondary structures they bridge. Similar to the Phe-X and His-X pairs, the far-sequence Tyr-X pairs (X partner > five apart in the sequence: 74%) show comparable secondary structures and conformers for either type of X partner, in contrast with the near-sequence Tyr-X pairs (26%). As the Phe-X pairs, the near-sequence Tyr-X pairs stabilize secondary structures, mainly the alpha- helices (positions 1, 3, and 4) and the beta-strands (position 2). Like the Phe-X and His-X pairs, most far-sequence Tyr-X pairs (34%) bridge beta-strands and only 11% bridge helices. As for the Phe-X and the His-X pairs, the X partners of the far-sequence Tyr-X pairs are frequently "above" the tyrosine ring with tilted and normal rings, whereas the X partner of the near-sequence Tyr-X pairs gradually moves from the "aside" to the "above" location, together with a progressive decrease of normal and increase of parallel rings, respectively. Unlike the His-X pairs, the interactions of the hetroatom in Tyr-X pairs are only favored with a sequence position +4 and over, owing to the spatial accessibility of the heteroatom.

Amino Acids, Aromatic↗

Functional analysis of the cell division protein FtsW of Escherichia coli.

Site-directed mutagenesis experiments combined with fluorescence microscopy shed light on the role of Escherichia coli FtsW, a membrane protein belonging to the SEDS family that is involved in peptidoglycan assembly during cell elongation, division, and sporulation. This essential cell division protein has 10 transmembrane segments (TMSs). It is a late recruit to the division site and is required for subsequent recruitment of penicillin-binding protein 3 (PBP3) catalyzing peptide cross-linking. The results allow identification of several domains of the protein with distinct functions. The localization of PBP3 to the septum was found to be dependent on the periplasmic loop located between TMSs 9 and 10. The E240-A249 amphiphilic peptide in the periplasmic loop between TMSs 7 and 8 appears to be a key element in the functioning of FtsW in the septal peptidoglycan assembly machineries. The intracellular loop (containing the R166-F178 amphiphilic peptide) between TMSs 4 and 5 and Gly 311 in TMS 8 are important components of the amino acid sequence-folding information.

Bacterial Proteins↗

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.

Amino Acids, Aromatic↗

Piracetam inhibits the lipid-destabilising effect of the amyloid peptide Abeta C-terminal fragment.

Amyloid peptide (Abeta) is a 40/42-residue proteolytic fragment of a precursor protein (APP), implicated in the pathogenesis of Alzheimer's disease. The hypothesis that interactions between Abeta aggregates and neuronal membranes play an important role in toxicity has gained some acceptance. Previously, we showed that the C-terminal domain (e.g. amino acids 29-42) of Abeta induces membrane permeabilisation and fusion, an effect which is related to the appearance of non-bilayer structures. Conformational studies showed that this peptide has properties similar to those of the fusion peptide of viral proteins i.e. a tilted penetration into membranes. Since piracetam interacts with lipids and has beneficial effects on several symptoms of Alzheimer's disease, we investigated in model membranes the ability of piracetam to hinder the destabilising effect of the Abeta 29-42 peptide. Using fluorescence studies and 31P and 2H NMR spectroscopy, we have shown that piracetam was able to significantly decrease the fusogenic and destabilising effect of Abeta 29-42, in a concentration-dependent manner. While the peptide induced lipid disorganisation and subsequent negative curvature at the membrane-water interface, the conformational analysis showed that piracetam, when preincubated with lipids, coats the phospholipid headgroups. Calculations suggest that this prevents appearance of the peptide-induced curvature. In addition, insertion of molecules with an inverted cone shape, like piracetam, into the outer membrane leaflet should make the formation of such structures energetically less favourable and therefore decrease the likelihood of membrane fusion.

Amyloid beta-Peptides↗

Insertion of X-ray structures of proteins in membranes.

Few structures of membrane proteins are known and their relationships with the membrane are unclear. In a previous report, 20 X-ray structures of transmembrane proteins were analyzed in silico for their orientation in a 36A-thick membrane [J. Mol. Graph. Model. 20 (2001) 235]. In this paper, we use the same approach to analyze how the insertion of the X-ray structures varies with the bilayer thickness. The protein structures are kept constant and, at each membrane thickness, the protein is allowed to tilt and rotate in order to accommodate at their best. The conditions are said to be optimal when the energy of insertion is minimal. The results show that most helix bundles require thicker membranes than porin barrels. Moreover, in a few instances, the ideal membrane thickness is unrealistic with respect to natural membranes supporting that the X-ray structure requires adaptation to stabilize in membrane. For instance, the squalene cyclase could adapt by bending the side chains of its ring of lysine and arginine in order to increase the hydrophobic surface in contact with membranes. We analyzed the distribution of amino acids in the water, interface and acyl chain layers of the membrane and compared with the literature.

Carrier Proteins↗

Analysis of accessible surface of residues in proteins.

We analyzed the total, hydrophobic, and hydrophilic accessible surfaces (ASAs) of residues from a nonredundant bank of 587 3D structure proteins. In an extended fold, residues are classified into three families with respect to their hydrophobicity balance. As expected, residues lose part of their solvent-accessible surface with folding but the three groups remain. The decrease of accessibility is more pronounced for hydrophobic than hydrophilic residues. Amazingly, Lysine is the residue with the largest hydrophobic accessible surface in folded structures. Our analysis points out a clear difference between the mean (other studies) and median (this study) ASA values of hydrophobic residues, which should be taken into consideration for future investigations on a protein-accessible surface, in order to improve predictions requiring ASA values. The different secondary structures correspond to different accessibility of residues. Random coils, turns, and beta-structures (outside beta-sheets) are the most accessible folds, with an average of 30% accessibility. The helical residues are about 20% accessible, and the difference between the hydrophobic and the hydrophilic residues illustrates the amphipathy of many helices. Residues from beta-sheets are the most inaccessible to solvent (10% accessible). Hence, beta-sheets are the most appropriate structures to shield the hydrophobic parts of residues from water. We also show that there is an equal balance between the hydrophobic and the hydrophilic accessible surfaces of the 3D protein surfaces irrespective of the protein size. This results in a patchwork surface of hydrophobic and hydrophilic areas, which could be important for protein interactions and/or activity.

Amino Acids↗

Revisiting the Ramachandran plot: hard-sphere repulsion, electrostatics, and H-bonding in the alpha-helix.

What determines the shape of the allowed regions in the Ramachandran plot? Although Ramachandran explained these regions in terms of 1-4 hard-sphere repulsions, there are discrepancies with the data where, in particular, the alphaR, alphaL, and beta-strand regions are diagonal. The alphaR-region also varies along the alpha-helix where it is constrained at the center and the amino terminus but diffuse at the carboxyl terminus. By analyzing a high-resolution database of protein structures, we find that certain 1-4 hard-sphere repulsions in the standard steric map of Ramachandran do not affect the statistical distributions. By ignoring these steric clashes (NH(i+1) and O(i-1)C), we identify a revised set of steric clashes (CbetaO, O(i-1)N(i+1), CbetaN(i+1), O(i-1)Cbeta, and O(i-1)O) that produce a better match with the data. We also find that the strictly forbidden region in the Ramachandran plot is excluded by multiple steric clashes, whereas the outlier region is excluded by only one significant steric clash. However, steric clashes alone do not account for the diagonal regions. Using electrostatics to analyze the conformational dependence of specific interatomic interactions, we find that the diagonal shape of the alphaR and alphaL-regions also depends on the optimization of the NH(i+1) and O(i-1)C interactions, and the diagonal beta-strand region is due to the alignment of the CO and NH dipoles. Finally, we reproduce the variation of the Ramachandran plot along the alpha-helix in a simple model that uses only H-bonding constraints. This allows us to rationalize the difference between the amino terminus and the carboxyl terminus of the alpha-helix in terms of backbone entropy.

Amino Acids↗

Aromatic side-chain interactions in proteins. I. Main structural features.

In a data set of 593 nonhomologous proteins from the PDB, we have analyzed the pairing of phenylalanine, tyrosine, tryptophan, and histidine residues with their closest aromatic partner. The frequency distribution of the shortest interatomic distance of partners is bimodal with a sharp peak at approximately 3.8 A and a wider one at a longer distance. Only the 3.8 A peak corresponds to direct ring-ring interactions thus aromatic pairs. The aromatic pairs were separated into two classes, near-sequence pairs and far-sequence pairs. Near sequence pairs stabilize local structure, and far-sequence pairs stabilize tertiary structure. Far-sequence pairs (74% of all pairs) mainly bridge two beta-strands, followed by pairs that bridge a beta-strand and a helix, and pairs that bridge a beta-strand and a random coil structure. Pairs that bridge helices are rare. The secondary structure of the near-sequence pairs depends on the partner distance in the sequence. When the partners are 1, 3, or 4 residues apart in the sequence, pairs are mostly found in helical structures. When the partners are two apart, pairs are mostly found in the same beta-strand. Analysis of the frequency of near sequence pairs supports the hypothesis that aromatic pairing occurs after, rather than before, the formation of secondary structures.

Amino Acids, Aromatic↗