Comparison of the solution conformations of a cell-adhesive peptide LBE and its reverse sequence EBL

T-cell adhesion is mediated by an ICAM-1/LFA-1 interaction; this interaction plays a crucial role in T-cell activation during immune response. LBE peptide, which is derived from the beta-subunit of LFA-1, has been shown to inhibit ICAM-1/LFA-1-mediated T-cell adhesion. In this work, we studied the solution conformations of LBE peptide and its reverse sequence (EBL) by NMR, CD and molecular dynamics simulations. Reverse peptides have been used as controls in biological studies. The effect of reversing the sequence of LBE to EBL peptides on their respective conformations is important in understanding their biological properties in vitro or in vivo. The NMR studies for these peptides were carried out in water and in TFE/water solvent systems. In 40% TFE/water, both peptides exhibited helical conformation. CD studies suggested that the LBE exhibits 30% helical conformation, while the EBL exhibits 20% helical conformation. From the NMR and MD simulation studies, it was evident that the peptides exhibited a stable helical conformation; a stable helical structure was found at Leu6 to Leu15 for LBE and at Gly9 to Leu17 for EBL. The helical conformations of LBE and EBL may be in equilibrium with other possible conformers; the other conformers contain loop and turn structures. Both peptides bind to divalent cations because the LBE is derived from the cation-binding region of the LFA-1. This study shows that reversing the peptide sequence did not alter the secondary structure of the corresponding sequence. Hence, caution must be exercised when using reverse peptides as controls in biological studies. This report will improve our ability to design a better inhibitor of ICAM-1/LFA-1 interaction.     

Determination of the primary structure of homologous proteins by sequence analysis of peptide mixtures

We propose a rapid method to determine the primary structure of a protein knowing the sequence of a homologous protein. The method consists in submitting both the reduced and alkylated proteins to an enzymatic or chemical hydrolysis and performing the sequence analysis of the peptide mixtures. The assessment of the unknown sequence and the degree of identity of the two proteins are reached by comparing the two sequence analyses. The sequences of all the possible peptides present in the two mixtures are reconstructed and the differences in the two sequences are determined. If necessary, the differences can be confirmed by performing a mass spectrometric analysis of the two mixtures. We used this procedure on two homologous proteins of known sequence to furnish an application example of the method.     

DNA conformations and their sequence preferences

The geometry of the phosphodiester backbone was analyzed for 7739 dinucleotides from 447 selected crystal structures of naked and complexed DNA. Ten torsion angles of a near-dinucleotide unit have been studied by combining Fourier averaging and clustering. Besides the known variants of the A-, B- and Z-DNA forms, we have also identified combined A + B backbone-deformed conformers, e.g. with α/γ switches, and a few conformers with a syn orientation of bases occurring e.g. in G-quadruplex structures. A plethora of A- and B-like conformers show a close relationship between the A- and B-form double helices. A comparison of the populations of the conformers occurring in naked and complexed DNA has revealed a significant broadening of the DNA conformational space in the complexes, but the conformers still remain within the limits defined by the A- and B- forms. Possible sequence preferences, important for sequence-dependent recognition, have been assessed for the main A and B conformers by means of statistical goodness-of-fit tests. The structural properties of the backbone in quadruplexes, junctions and histone-core particles are discussed in further detail.     

Clines in clock genes: fine-tuning circadian rhythms to the environment

The dissection of the circadian clock into its molecular components represents the most striking and well-studied example of a gene regulatory network underlying a complex behavioural trait. By contrast, the evolutionary analysis of the clock has developed more slowly. Here we review studies that have surveyed intraspecific clock gene variation over large geographical areas and have discovered latitudinal clines in gene frequencies. Such spatial patterns traditionally suggest that natural selection shapes genetic variation, but it is equally possible that population history, or a mixture of demography and selection, could contribute to the clines. We discuss how population genetics, together with functional assays, can illuminate these possible cases of natural selection in Drosophila clock genes.     

DynaPred: a structure and sequence based method for the prediction of MHC class I binding peptide sequences and conformations

MOTIVATION: The binding of endogenous antigenic peptides to MHC class I molecules is an important step during the immunologic response of a host against a pathogen. Thus, various sequence- and structure-based prediction methods have been proposed for this purpose. The sequence-based methods are computationally efficient, but are hampered by the need of sufficient experimental data and do not provide a structural interpretation of their results. The structural methods are data-independent, but are quite time-consuming and thus not suited for screening of whole genomes. Here, we present a new method, which performs sequence-based prediction by incorporating information obtained from molecular modeling. This allows us to perform large databases screening and to provide structural information of the results. RESULTS: We developed a SVM-trained, quantitative matrix-based method for the prediction of MHC class I binding peptides, in which the features of the scoring matrix are energy terms retrieved from molecular dynamics simulations. At the same time we used the equilibrated structures obtained from the same simulations in a simple and efficient docking procedure. Our method consists of two steps: First, we predict potential binders from sequence data alone and second, we construct protein-peptide complexes for the predicted binders. So far, we tested our approach on the HLA-A0201 allele. We constructed two prediction models, using local, position-dependent (DynaPred(POS)) and global, position-independent (DynaPred) features. The former model outperformed the two sequence-based methods used in our evaluation; the latter shows a much higher generalizability towards other alleles than the position-dependent models. The constructed peptide structures can be refined within seconds to structures with an average backbone RMSD of 1.53 A from the corresponding experimental structures.     

An Escherichia coli twin-arginine signal peptide switches between helical and unstructured conformations depending on the hydrophobicity of the environment.

The Tat system catalyzes the transport of folded globular proteins across the bacterial plasma membrane and the chloroplast thylakoid. It recognizes cleavable signal peptides containing a critical twin-arginine motif but little is known of the overall structure of these peptides. In this report, we have analyzed the secondary structure of the SufI signal peptide, together with those of two nonfunctional variants in which the region around the twin-arginine, RRQFI, is replaced by KKQFI or RRQAA. Circular dichroism studies show that the SufI peptide exists as an unstructured peptide in aqueous solvent with essentially no stable secondary structure. In membrane-mimetic environments such as SDS micelles or water/trifluoroethanol, however, the peptide adopts a structure containing up to about 40% alpha-helical content. Secondary structure predictions and molecular modelling programs strongly suggest that the helical region begins at, or close to, the twin-arginine motif. Studies on the thermal stability of the helix demonstrate a sharp transition between the unstructured and helical states, suggesting that the peptide exists in one of two distinct states. The two nonfunctional peptides exhibit almost identical spectra and properties to the wild-type SufI peptide, indicating that it is the arginine sidechains, and not their contribution to the helical structure, that are critical in this class of peptide.     

Conformational study of the preferred conformations of the peptide sequence VP3(110-121) of HAV by circular dichroism and molecular mechanics

A conformational analysis of the fragment 110–121 of VP3 coating protein of the hepatitis A virus was carried out using circular dichroism spectroscopy and computational studies. The latter studies indicate the tendency of the peptide to adopt hairpin-type structures. Circular dichroism experiments indicate that, in spite of the fact that the isolated peptide exhibits no structure under different experimental conditions, negatively charged liposomes induce a secondary structure that agrees with the results of the computational study.     

Conformational study of the preferred conformations of the peptide sequence VP3(110-121) of HAV by circular dichroism and molecular mechanics

Summary A conformational analysis of the fragment 110–121 of VP3 coating protein of the hepatitis A virus was carried out using circular dichroism spectroscopy and computational studies. The latter studies indicate the tendency of the peptide to adopt hairpin-type structures. Circular dichroism experiments indicate that, in spite of the fact that the isolated peptide exhibits no structure under different experimental conditions, negatively charged liposomes induce a secondary structure that agrees with the results of the computational study.     

Short peptide amyloid organization: stabilities and conformations of the islet amyloid peptide NFGAIL

Experimentally, short peptides have been shown to form amyloids similar to those of their parent proteins. Consequently, they present useful systems for studies of amyloid conformation. Here we simulate extensively the NFGAIL peptide, derived from the human islet amyloid polypeptide (residues 22-27). We simulate different possible strand/sheet organizations, from dimers to nonamers. Our simulations indicate that the most stable conformation is an antiparallel strand orientation within the sheets and parallel between sheets. Consistent with the alanine mutagenesis, we find that the driving force is the hydrophobic effect. Whereas the NFGAIL forms stable oligomers, the NAGAIL oligomer is unstable, and disintegrates very quickly after the beginning of the simulation. The simulations further identify a minimal seed size. Combined with our previous simulations of the prion-derived AGAAAAGA peptide, AAAAAAAA, and the Alzheimer Abeta fragments 16-22, 24-36, 16-35, and 10-35, and the solid-state NMR data for Abeta fragments 16-22, 10-35, and 1-40, some insight into the length and the sequence matching effects may be obtained.     

Helical membrane protein conformations and their environment

Evidence that membrane proteins respond conformationally and functionally to their environment is growing. Structural models, by necessity, have been characterized in preparations where the protein has been removed from its native environment. Different structural methods have used various membrane mimetics that have recently included lipid bilayers as a more native-like environment. Structural tools applied to lipid bilayer-embedded integral proteins are informing us about important generic characteristics of how membrane proteins respond to the lipid environment as compared with their response to other nonlipid environments. Here, we review the current status of the field, with specific reference to observations of some well-studied α-helical membrane proteins, as a starting point to aid the development of possible generic principles for model refinement.     

Position fine-tuning of caudal primary motoneurons in the zebrafish spinal cord

In zebrafish embryos, each myotome is typically innervated by three primary motoneurons (PMNs): the caudal primary (CaP), middle primary (MiP) and rostral primary (RoP). PMN axons first exit the spinal cord through a single exit point located at the midpoint of the overlying somite, which is formed beneath the CaP cell body and is pioneered by the CaP axon. However, the placement of CaP cell bodies with respect to corresponding somites is poorly understood. Here, we determined the early events in CaP cell positioning using neuropilin 1a (nrp1a):gfp transgenic embryos in which CaPs were specifically labeled with GFP. CaP cell bodies first exhibit an irregular pattern in presence of newly formed corresponding somites and then migrate to achieve their proper positions by axonogenesis stages. CaPs are generated in excess compared with the number of somites, and two CaPs often overlap at the same position through this process. Next, we showed that CaP cell bodies remain in the initial irregular positions after knockdown of Neuropilin1a, a component of the class III semaphorin receptor. Irregular CaP position frequently results in aberrant double exit points of motor axons, and secondary motor axons form aberrant exit points following CaP axons. Its expression pattern suggests that sema3ab regulates the CaP position. Indeed, irregular CaP positions and exit points are induced by Sema3ab knockdown, whose ectopic expression can alter the position of CaP cell bodies. Results suggest that Semaphorin-Neuropilin signaling plays an important role in position fine-tuning of CaP cell bodies to ensure proper exit points of motor axons.     

Formation of the local secondary structure of proteins: local sequence or environment

An analysis of the occurrence of tetrapeptides in 35 globular proteins for alpha-helix, beta-structure and coil was performed. We concluded that: the conformation of a short polypeptide segment cannot be determined on the basis of the knowledge of the amino acid sequence only; local structures of a protein are formed as the result of interactions within the whole structural domain of the protein as well as interactions with the environment.     

Perturbation of peptide conformations induced in anisotropic environments

Reversed-phase high performance liquid chromatography (RP HPLC) has been found to be a convenient and powerful tool for the study of the secondary structure of peptides. Here, the ability of proline to perturb the secondary structures of peptides induced at aqueous-lipid interfaces and the induced conformation of polyproline peptides were investigated by means of RP HPLC. For these studies, four different complete sets of substitution analogues of model peptides expected to have specific induced conformations were used. In the first two studies, a single lysine was “walked” through two 18-residue polyproline sequences (one N-acetylated, the other not). In the remaining two studies, a proline was “walked” through two different sequences that had been found earlier to be induced into an α-helical conformation during RP HPLC (an 18-residue polyalanine sequence and the amphipathic 14-residue sequence Ac-LLKLLKKLLKKLKK-NH₂). Sixty-eight individual analogues were synthesized for this study and the effect of the respective substitutions on retention times was determined. The results are consistent with the concept that, upon interaction with the C-18 of the stationary phase during RP HPLC, polyproline is induced into a type II helical conformation, polyalanine into an α-helical conformation, and Ac-LLKLLKKLLKKLKK-NH₂ into an amphipathic α-helical array. In an extension of this study, the antimicrobial activities of Ac-LLKLLKKLLKKLKK-NH₂ and its 18 proline substitution analogues were found to be inversely correlated with their RP HPLC retention times.     

Predicting residue solvent accessibility from protein sequence by considering the sequence environment

The solvent accessibility of each residue is predicted on the basis of the protein sequence. A set of 338 monomeric, non-homologous and high-resolution protein crystal structures is used as a learning set and a jackknife procedure is applied to each entry. The prediction is based on the comparison of the observed and the average values of the solvent-accessible area. It appears that the prediction accuracy is significantly improved by considering the residue types preceding and/or following the residue whose accessibility must be predicted. In contrast, the separate treatment of different secondary structural types does not improve the quality of the prediction. It is furthermore shown that the residue accessibility is much better predicted in small than in larger proteins. Such a discrepancy must be carefully considered in any algorithm for predicting residue accessibility.     

Structural compromise of disallowed conformations in peptide and protein structures

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.     

Isolation from chicken antrum, and primary amino acid sequence of a novel 36-residue peptide of the gastrin/CCK family

A peptide that cross-reacted with C-terminal gastrin/CCK antisera was isolated from chicken antral extracts by a combination of gel filtration and reversed-phase HPLC. The sequence was: Phe-Leu-Pro-His- Val-Phe-Ala-Glu-Leu-Ser-Asp-Arg-Lys-Gly-Phe-Val-Gln-Gly-Asn-Gly-Ala- Val-Glu-Ala-Leu-His-Asp-His-Phe-Tyr-Pro-Asp-Trp-Met-Asp-Phe(NH2). Aside from the C-terminal tetrapeptide and the Tyr residue, the molecule does not resemble other known forms of gastrin or CCK. The peptide was a potent stimulus of avian gastric acid but not pancreatic secretion. The results have important implications for the structure-activity and evolutionary relationships of the gastrin/CCK family.     

Synthesis and activities of cyclic thrombin-receptor-derived peptide analogues of the Ser42-Phe-Leu-Leu-Arg46 motif sequence containing d-Phe and/or d-Arg

Summary Cyclic analogues of the active thrombin receptor peptide SFLLR (TRP_42–46) containingd-Phe and/ord-Arg have been prepared by the solid-phase method, purified by reversed-phase HPLC and bioassayed in a rat smooth muscle contractile assay. Cyclization was achieved by forming an amide linkage between the-NH₂ and-COOH groups of the two leucine residues located at the N- and C-terminal positions of the linear protected precursor H₂N-Leu-Arg(Pmc)-Y-Phe-Leu-OH (Y=Gly,Acp) using 2-(1H-benzotriazole-1-yl)-1,1,3,3-tetramethyluronium hexafluoroborate borate (HBTU) or 2-(1H-benzotriazole-1-yl)-1,1,3,3-tetramethyluronium tetrafluoroborate (TBTU) as coupling reagents andN,N-diisopropylethylamine (DIPEA) in high dilution. Their structure was confirmed by fast by fast atom bombardment mass spectrometry and NMR methods. The cyclic peptides c-fLLrG, c-fLLRG, c-FLLrG and c-fLLrAcp, c-FLLrAcp so synthesized were assessed for their contractile activity in a rat gastric longitudinal muscle bioassay system which has been used previously to evaluate the biological activities of linear thrombin-receptor-derived polypeptides such as SFLLR (P5) and SFLLR-NH₂ (P5-NH₂).     

A protein sequence that can encode native structure by disfavoring alternate conformations

The linear sequence of amino acids contains all the necessary information for a protein to fold into its unique three-dimensional structure. Native protein sequences are known to accomplish this by promoting the formation of stable, kinetically accessible structures. Here we describe a Pro residue in the center of the third transmembrane helix of the cystic fibrosis transmembrane conductance regulator that promotes folding by a distinct mechanism: disfavoring the formation of a misfolded structure. The generality of this mechanism is supported by genome-wide transmembrane sequence analyses. Furthermore, the results provide an explanation for the increased frequency of Pro residues in transmembrane alpha-helices. Incorporation by nature of such 'negative folding determinants', aimed at preventing the formation of off-pathway structures, represents an additional mechanism by which folding information is encoded within the evolved sequences of proteins.