Computational identification of self-inhibitory peptides from envelope proteins

Fusion process is known to be the initial step of viral infection and hence targeting the entry process is a promising strategy to design antiviral therapy. The self-inhibitory peptides derived from the enveloped (E) proteins function to inhibit the proteinprotein interactions in the membrane fusion step mediated by the viral E protein. Thus, they have the potential to be developed into effective antiviral therapy. Herein, we have developed a Monte Carlo-based computational method with the aim to identify and optimize potential peptide hits from the E proteins. The stability of the peptides, which indicates their potential to bind in situ to the E proteins, was evaluated by two different scoring functions, dipolar distance-scaled, finite, ideal-gas reference state and residue-specific all-atom probability discriminatory function. The method was applied to a-helical Class I HIV-1 gp41, beta-sheet Class II Dengue virus (DENV) type 2 E proteins, as well as Class III Herpes Simplex virus-1 (HSV-1) glycoprotein, a E protein with a mixture of a-helix and beta-sheet structural fold. The peptide hits identified are in line with the druggable regions where the self-inhibitory peptide inhibitors for the three classes of viral fusion proteins were derived. Several novel peptides were identified from either the hydrophobic regions or the functionally important regions on Class II DENV-2 E protein and Class III HSV-1 gB. They have potential to disrupt the proteinprotein interaction in the fusion process and may serve as starting points for the development of novel inhibitors for viral E proteins. 

A conserved spider silk domain acts as a molecular switch that controls fibre assembly

A huge variety of proteins are able to form fibrillar structures(1), especially at high protein concentrations. Hence, it is surprising that spider silk proteins can be stored in a soluble form at high concentrations and transformed into extremely stable fibres on demand(2,3). Silk proteins are reminiscent of amphiphilic block copolymers containing stretches of polyalanine and glycine-rich polar elements forming a repetitive core flanked by highly conserved non-repetitive amino-terminal(4,5) and carboxy-terminal(6) domains. The N-terminal domain comprises a secretion signal, but further functions remain unassigned. The C-terminal domain was implicated in the control of solubility and fibre formation(7) initiated by changes in ionic composition(8,9) and mechanical stimuli known to align the repetitive sequence elements and promote beta-sheet formation(10-14). However, despite recent structural data(15), little is known about this remarkable behaviour in molecular detail. Here we present the solution structure of the C-terminal domain of a spider dragline silk protein and provide evidence that the structural state of this domain is essential for controlled switching between the storage and assembly forms of silk proteins. In addition, the C-terminal domain also has a role in the alignment of secondary structural features formed by the repetitive elements in the backbone of spider silk proteins, which is known to be important for the mechanical properties of the fibre.

Protein Folding Simulations: Confinement, External Fields and Sequence Design

The present Thesis looks at the problem of protein folding using Monte Carlo and Langevin simulations, three topics in protein folding have been studied: 1) the effect of confining potential barriers, 2) the effect of a static external field and 3) the design of amino acid sequences which fold in a short time and which have a stable native state (global minimum). Regarding the first topic, we studied the confinement of a small protein of 16 amino acids known as 1NJ0 (PDB code) which has a beta-sheet structure as a native state. The confinement of proteins occurs frequently in the cell environment. Some molecules called Chaperones, present in the cytoplasm, capture the unfolded proteins in their interior and avoid the formation of aggregates and misfolded proteins. This mechanism of confinement mediated by Chaperones is not yet well understood. In the present work we considered two kinds of potential barriers which try to mimic the confinement induced by a Chaperon molecule. The first kind of potential was a purely repulsive barrier whose only effect is to create a cavity where the protein folds up correctly. The second kind of potential was a barrier which includes both attractive and repulsive effects. We performed Wang-Landau simulations to calculate the thermodynamical properties of 1NJ0. From the free energy landscape plot we found that 1NJ0 has two intermediate states in the bulk (without confinement) which are clearly separated from the native and the unfolded states. For the case of the purely repulsive barrier we found that the intermediate states get closer to each other in the free energy landscape plot and eventually they collapse into a single intermediate state. The unfolded state is more compact, compared to that in the bulk, as the size of the barrier decreases. For an attractive barrier modifications of the states (native, unfolded and intermediates) are observed depending on the degree of attraction between the protein and the walls of the barrier. The strength of the attraction is measured by the parameter $\epsilon$. A purely repulsive barrier is obtained for $\epsilon=0$ and a purely attractive barrier for $\epsilon=1$. The states are changed slightly for magnitudes of the attraction up to $\epsilon=0.4$. The disappearance of the intermediate states of 1NJ0 is already observed for $\epsilon =0.6$. A very high attractive barrier ($\epsilon \sim 1.0$) produces a completely denatured state. In the second topic of this Thesis we dealt with the interaction of a protein with an external electric field. We demonstrated by means of computer simulations, specifically by using the Wang-Landau algorithm, that the folded, unfolded, and intermediate states can be modified by means of a field. We have found that an external field can induce several modifications in the thermodynamics of these states: for relatively low magnitudes of the field ($<2.06 \times 10^8$ V/m) no major changes in the states are observed. However, for higher magnitudes than ($6.19 \times 10^8$ V/m) one observes the appearance of a new native state which exhibits a helix-like structure. In contrast, the original native state is a $\beta$-sheet structure. In the new native state all the dipoles in the backbone structure are aligned parallel to the field. The design of amino acid sequences constitutes the third topic of the present work. We have tested the Rate of Convergence criterion proposed by D. Gridnev and M. Garcia ({\it work unpublished}). We applied it to the study of off-lattice models. The Rate of Convergence criterion is used to decide if a certain sequence will fold up correctly within a relatively short time. Before the present work, the common way to decide if a certain sequence was a good/bad folder was by performing the whole dynamics until the sequence got its native state (if it existed), or by studying the curvature of the potential energy surface. There are some difficulties in the last two approaches. In the first approach, performing the complete dynamics for hundreds of sequences is a rather challenging task because of the CPU time needed. In the second approach, calculating the curvature of the potential energy surface is possible only for very smooth surfaces. The Rate of Convergence criterion seems to avoid the previous difficulties. With this criterion one does not need to perform the complete dynamics to find the good and bad sequences. Also, the criterion does not depend on the kind of force field used and therefore it can be used even for very rugged energy surfaces.

Polimorfismos del gen Butirilcolinesterasa responsables de reacciones adversas en pacientes consumidores de “cocaína”

La butirilcolinesterasa humana (BChE; EC 3.1.1.8) es una enzima polimórfica sintetizada en el hígado y en el tejido adiposo, ampliamente distribuida en el organismo y encargada de hidrolizar algunos ésteres de colina como la procaína, ésteres alifáticos como el ácido acetilsalicílico, fármacos como la metilprednisolona, el mivacurium y la succinilcolina y drogas de uso y/o abuso como la heroína y la cocaína. Es codificada por el gen BCHE (OMIM 177400), habiéndose identificado más de 100 variantes, algunas no estudiadas plenamente, además de la forma más frecuente, llamada usual o silvestre. Diferentes polimorfismos del gen BCHE se han relacionado con la síntesis de enzimas con niveles variados de actividad catalítica. Las bases moleculares de algunas de esas variantes genéticas han sido reportadas, entre las que se encuentra las variantes Atípica (A), fluoruro-resistente del tipo 1 y 2 (F-1 y F-2), silente (S), Kalow (K), James (J) y Hammersmith (H). En este estudio, en un grupo de pacientes se aplicó el instrumento validado Lifetime Severity Index for Cocaine Use Disorder (LSI-C) para evaluar la gravedad del consumo de “cocaína” a lo largo de la vida. Además, se determinaron Polimorfismos de Nucleótido Simple (SNPs) en el gen BCHE conocidos como responsables de reacciones adversas en pacientes consumidores de “cocaína” mediante secuenciación del gen y se predijo el efecto delos SNPs sobre la función y la estructura de la proteína, mediante el uso de herramientas bio-informáticas. El instrumento LSI-C ofreció resultados en cuatro dimensiones: consumo a lo largo de la vida, consumo reciente, dependencia psicológica e intento de abandono del consumo. Los estudios de análisis molecular permitieron observar dos SNPs codificantes (cSNPs) no sinónimos en el 27.3% de la muestra, c.293A>G (p.Asp98Gly) y c.1699G>A (p.Ala567Thr), localizados en los exones 2 y 4, que corresponden, desde el punto de vista funcional, a la variante Atípica (A) [dbSNP: rs1799807] y a la variante Kalow (K) [dbSNP: rs1803274] de la enzima BChE, respectivamente. Los estudios de predicción In silico establecieron para el SNP p.Asp98Gly un carácter patogénico, mientras que para el SNP p.Ala567Thr, mostraron un comportamiento neutro. El análisis de los resultados permite proponer la existencia de una relación entre polimorfismos o variantes genéticas responsables de una baja actividad catalítica y/o baja concentración plasmática de la enzima BChE y algunas de las reacciones adversas ocurridas en pacientes consumidores de cocaína.

Nuclear magnetic resonance structure of the nucleic acid-binding domain of severe acute respiratory syndrome coronavirus nonstructural protein 3

The nuclear magnetic resonance (NMR) structure of a globular domain of residues 1071 to 1178 within the previously annotated nucleic acid-binding region (NAB) of severe acute respiratory syndrome coronavirus nonstructural protein 3 (nsp3) has been determined, and N- and C-terminally adjoining polypeptide segments of 37 and 25 residues, respectively, have been shown to form flexibly extended linkers to the preceding globular domain and to the following, as yet uncharacterized domain. This extension of the structural coverage of nsp3 was obtained from NMR studies with an nsp3 construct comprising residues 1066 to 1181 [ nsp3(1066-1181)] and the constructs nsp3(1066-1203) and nsp3(1035-1181). A search of the protein structure database indicates that the globular domain of the NAB represents a new fold, with a parallel four-strand beta-sheet holding two alpha-helices of three and four turns that are oriented antiparallel to the beta-strands. Two antiparallel two-strand beta-sheets and two 3(10)-helices are anchored against the surface of this barrel-like molecular core. Chemical shift changes upon the addition of single-stranded RNAs (ssRNAs) identified a group of residues that form a positively charged patch on the protein surface as the binding site responsible for the previously reported affinity for nucleic acids. This binding site is similar to the ssRNA-binding site of the sterile alpha motif domain of the Saccharomyces cerevisiae Vts1p protein, although the two proteins do not share a common globular fold.

Ribonucleocapsid formation of severe acute respiratory syndrome coronavirus through molecular action of the N-terminal domain of N protein

Conserved among all coronaviruses are four structural proteins: the matrix (M), small envelope (E), and spike (S) proteins that are embedded in the viral membrane and the nucleocapsid phosphoprotein (N), which exists in a ribonucleoprotein complex in the lumen. The N-terminal domain of coronaviral N proteins (N-NTD) provides a scaffold for RNA binding, while the C-terminal domain (N-CTD) mainly acts as oligomerization modules during assembly. The C terminus of the N protein anchors it to the viral membrane by associating with M protein. We characterized the structures of N-NTD from severe acute respiratory syndrome coronavirus (SARS-CoV) in two crystal forms, at 1.17 A (monoclinic) and at 1.85 A (cubic), respectively, resolved by molecular replacement using the homologous avian infectious bronchitis virus (IBV) structure. Flexible loops in the solution structure of SARS-CoV N-NTD are now shown to be well ordered around the beta-sheet core. The functionally important positively charged beta-hairpin protrudes out of the core, is oriented similarly to that in the IBV N-NTD, and is involved in crystal packing in the monoclinic form. In the cubic form, the monomers form trimeric units that stack in a helical array. Comparison of crystal packing of SARS-CoV and IBV N-NTDs suggests a common mode of RNA recognition, but they probably associate differently in vivo during the formation of the ribonucleoprotein complex. Electrostatic potential distribution on the surface of homology models of related coronaviral N-NTDs suggests that they use different modes of both RNA recognition and oligomeric assembly, perhaps explaining why their nucleocapsids have different morphologies.

Amino terminal interaction in the prion protein identified using fusion to green fluorescent protein

In contrast to the well-characterized carboxyl domain, the amino terminal half of the mature cellular prion protein has no defined structure. Here, following fusion of mouse prion protein fragments to green fluorescence protein as a reporter of protein stability, we report extreme variability in fluorescence level that is dependent on the prion fragment expressed. In particular, exposure of the extreme amino terminus in the context of a truncated prion protein molecule led to rapid degradation, whereas the loss of only six amino terminal residues rescued high level fluorescence. Study of the precise endpoints and residue identity associated with high fluorescence suggested a domain within the amino terminal half of the molecule defined by a long-range intramolecular interaction between 23KKRPKP28 and 143DWED146 and dependent upon the anti-parallel beta-sheet ending at residue 169 and normally associated with the structurally defined carboxyl terminal domain. This previously unreported interaction may be significant for understanding prion bioactivity and for structural studies aimed at the complete prion structure.

Self assembly of a model amphiphilic phenylalanine peptide/polyethylene glycol block copolymer in aqueous solution

There has been great interest recently in peptide amphiphiles and block copolymers containing biomimetic peptide sequences due to applications in bionanotechnology. We investigate the self-assembly of the peptide-PEG amphiphile FFFF-PEG5000 containing the hydrophobic sequence of four phenylalanine residues conjugated to PEG of molar mass 5000. This serves as a simple model peptide amphiphile. At very low concentration, association of hydrophobic aromatic phenylalanine residues occurs, as revealed by circular dichroism and UV/vis fluorescence experiments. A critical aggregation concentration associated with the formation of hydrophobic domains is determined through pyrene fluorescence assays. At higher concentration, defined beta-sheets develop as revealed by FTIR spectroscopy and X-ray diffraction. Transmission electron microscopy reveals self-assembled straight fibril structures. These are much shorter than those observed for amyloid peptides, the finite length may be set by the end cap energy due to the hydrophobicity of phenylalanine. The combination of these techniques points to different aggregation processes depending on concentration. Hydrophobic association into irregular aggregates occurs at low concentration, well-developed beta-sheets only developing at higher concentration. Drying of FFFF-PEG5000 solutions leads to crystallization of PEG, as confirmed by polarized optical microscopy (POM), FTIR and X-ray diffraction (XRD). PEG crystallization does not disrupt local beta-sheet structure (as indicated by FTIR and XRD). However on longer lengthscales the beta-sheet fibrillar structure is perturbed because spheruilites from PEG crystallization are observed by POM. (C) 2009 Elsevier B.V. All rights reserved.

The role of protecting groups in the formation of organogels through a nano-fibrillar network formed by self-assembling terminally protected tripeptides

A series of eight synthetic self-assembling terminally blocked tripeptides have been studied for gelation. Some of them form gels in various aromatic solvents including benzene, toluene, xylene, and chlorobenzene. It has been found that the protecting groups play an important role in the formation of organogels. It has been observed that, if the C-terminal has been changed from methyl ester to ethyl ester the gelation property does not change significantly (keeping the N-terminal protecting group same), while the change of the protecting group from ethyl ester to isopropyl ester completely abolishes the gelation property. Similarly, keeping the identical C-terminal protecting group (methyl ester) the results of the gelation study indicate that the substitution of N-terminal protection Boc-(tert-butyloxycarbonyl) to Cbz-(benzyloxycarbonyl) does change the gelation property insignificantly, while the change from Boc- to pivaloyl (Piv-) or acetyl (Ac-) group completely eliminates the gelation property. Morphological studies of the dried gels of two of the peptides indicate the presence of an entangled nano-fibrillar network that might be responsible for gelation. FTIR studies of the gels demonstrate that an intermolecular hydrogen bonding network is formed during gelation. Results of X-ray powder diffraction studies for these gelator peptides in different states (dried gels, gel, and bulk solids) reflected that the structure in the wet gel is distinctly different from the dried gel and solid state structures. Single crystal X-ray diffraction studies of a non-gelator peptide, which is structurally similar to the gelator molecules reveal that the peptide forms an antiparallel beta-sheet structure in crystals. (c) 2007 Elsevier Ltd. All rights reserved.

Nematic and columnar ordering of a PEG-peptide conjugate in aqueous solution

The self-assembly in aqueous solution of a PEG-peptide conjugate is studied by spectroscopy, electron microscopy, rheology and small-angle Xray and neutron scattering (SAXS and SANS). The peptide fragment, FFKLVFF is based on fragment KLVFF of the amyloid beta-peptide, A beta(16-20), extended by two hydrophobic phenylalanine units. This is conjugated to PEG which confers water solubility and leads to distinct self-assembled structures. Small-angle scattering reveals the formation of cylindrical fibrils comprising a peptide core and PEG corona. This constrained structure leads to a model parallel beta-sheet self-assembled structure with a radial arrangement of beta sheets. Oil increasing concentration, successively nematic and hexagonal columnar phases are formed. The flow-induced alignment of both structures was studied in situ by SANS using a Couette cell. Shear-induced alignment is responsible for the shear thinning behaviour observed by dynamic shear rheometry. Incomplete recovery of moduli after cessation of shear is consistent with the observation from SANS of retained orientation in the sample.

Formation of fibrillar structures through self-assembly of designed peptide turns

Three tripeptides Boc-Phe-Aib-Val-OMe (1), Boc-Leu-Aib-p-NA-NO2 (2) and Boc-Pro-Aib-m-NA-NO2 (3) (Aib: alpha-aminoisobutyric acid; p- and m-NA: para- and meta-nitroaniline) have been designed by incorporating aromatic rings to study the self-assembly and fibril formation. Single crystal X-ray diffraction studies show that all the peptides adopt turn-like structures that are self-assembled through intermolecular hydrogen bonds and van der Waals interactions to create layers of beta-sheets. Solvent dependent NMR titration and CD studies show that the turn structures of the peptides also exist in the solution phase. The field emission scanning electron microscopic (FE-SEM) images of the peptides in the solid state reveal fibrillar structures of flat morphology that are formed through beta-sheet mediated self-assembly of the preorganized turn building blocks.

Self-assembly of peptide nanotubes in an organic solvent

The self-assembly of a modified fragment of the amyloid beta peptide, based on sequence A beta(16-20), KLVFF, extended to give AAKLVFF is studied in methanol. Self-assembly into peptide nanotubes is observed, as confirmed by electron microscopy and small-angle X-ray scattering. The secondary structure of the peptide is probed by FTIR and circular dichroism, and UV/visible spectroscopy provides evidence for the important role of aromatic interactions between phenylalanine residues in driving beta-sheet self-assembly. The beta-sheets wrap helically to form the nanotubes, the nanotube wall comprising four wrapped beta-sheets. At higher concentration, the peptide nanotubes form a nematic phase that exhibits spontaneous flow alignment as observed by small-angle neutron scattering.

A short water-soluble self-assembling peptide forms amyloid-like fibrils

A water-soluble tripeptide Val-Ile-Ala (VIA) 1, bearing sequence identity with the C-terminal portion of the Alzheimer A beta-peptide (A beta(40-42)), self-assembles, in crystalline form, to produce an intermolecularly hydrogen bonded supramolecular beta-sheet structure which self-associates to form straight, unbranched nanofibrils exhibiting amyloid-like behavior; in contrast, the synthetic tripeptide Ala-Val-Ile (AVI) 2 self-assembles to produce a beta-sheet structure that forms branched nanofibrils which do not show any characteristic features of amyloid-like fibrils.

Wheat glutenin proteins assemble into a nanostructure with unusual structural features

The proteins of wheat have a known propensity to aggregate into a variety of forms. We report here a novel nanostructure from wheat proteins, derived from a crude extract of high molecular weight glutenins. The structure was characterised by a significant thioflavin T (ThT) fluorescence and a fibrillar morphology by transmission electron microscopy (TEM). The ThT fluorescence and TEM data are suggestive of an amyloid structure, but the X-ray fibre diffraction data show a reflection pattern (4.02, 4.2-4.3, 4.6, 12.9,19.3 and 38.7 angstrom) inconsistent with both the classic amyloid form and the previously described beta-helix structure. The 4.6 angstrom reflection is consistent with that predicted for the amyloid inter-beta-strand, and the absence of the inter-beta-sheet distance at approximate to 10-11 angstrom is not unprecedented in amyloid-like structures. However, our observed X-ray reflection pattern has not been previously reported and suggests a novel wheat glutenin nanostructure. (C) 2008 Elsevier Ltd. All rights reserved.

Polymer conformation structure of wheat proteins and gluten subfractions revealed by ATR-FTIR

The polymer conformation structure of gluten extracted from a Polish wheat cultivar, Korweta, and gluten subtractions obtained from 2 U.K. breadmaking and biscuit flour cultivars, Hereward and Riband, was investigated using attenuated total reflectance Fourier transform infrared spectroscopy (ATR-FTIR). The results showed the conformation of proteins varied between flour, hydrated flour, and hydrated gluten. The beta-sheet structure increased progressively from flour to hydrated flour and to hydrated gluten. In hydrated gluten protein fractions comprising gliadin, soluble glutenin, and gel protein, beta-sheet structure increased progressively from soluble gliadin and glutenin to gluten and gel protein; beta-sheet content was also greater in the gel protein from the breadmaking flour Hereward than the biscuit flour Riband.