Use of a novel Hill-climbing genetic algorithm in protein folding simulations

We have developed a novel Hill-climbing genetic algorithm (GA) for simulation of protein folding. The program (written in C) builds a set of Cartesian points to represent an unfolded polypeptide's backbone. The dihedral angles determining the chain's configuration are stored in an array of chromosome structures that is copied and then mutated. The fitness of the mutated chain's configuration is determined by its radius of gyration. A four-helix bundle was used to optimise simulation conditions, and the program was compared with other, larger, genetic algorithms on a variety of structures. The program ran 50% faster than other GA programs. Overall, tests on 100 non-redundant structures gave comparable results to other genetic algorithms, with the Hill-climbing program running from between 20 and 50% faster. Examples including crambin, cytochrome c, cytochrome B and hemerythrin gave good secondary structure fits with overall alpha carbon atom rms deviations of between 5 and 5.6 Angstrom with an optimised hydrophobic term in the fitness function. (C) 2003 Elsevier Ltd. All rights reserved.

Use of essential and molecular dynamics to study gamma B-crystallin unfolding after non-enzymic post-translational modifications

Essential and Molecular Dynamics (ED/MD) have been used to model the conformational changes of a protein implicated in a conformational disease-cataract, the largest cause of blindness in the world-after non-enzymic post-translational modification. Cyanate modification did not significantly alter flexibility, while the Schiff's base adduct produced a more flexible N-terminal domain, and intra-secondary structure regions, than either the cyanate adduct or the native structure. Glycation also increased linker flexibility and disrupted the charge network. A number of post-translational adducts showed structural disruption around Cys15 and increased linker flexibility; this may be important in subsequent protein aggregation. Our modelling results are in accord with experimental evidence, and show that ED/MD is a useful tool in modelling conformational changes in proteins implicated in disease processes. (C) 2003 Published by Elsevier Ltd.

Interactions of bovine serum albumin with ethylene oxide/butylene oxide copolymers in aqueous solution

The interactions of bovine serum albumin (BSA) with three ethylene oxide/butylene oxide (E/B) copolymers having different block lengths and varying molecular architectures is examined in this study in aqueous solutions. Dynamic light scattering (DLS) indicates the absence of BSA-polymer binding in micellar systems of copolymers with lengthy hydrophilic blocks. On the contrary, stable protein-polyrner aggregates were observed in the case of E18B10 block copolymer. Results from DLS and SAXS suggest the dissociation of E/B copolymer micelles in the presence of protein and the absorption of polymer chains to BSA surface. At high protein loadings, bound BSA adopts a more compact conformation in solution. The secondary structure of the protein remains essentially unaffected even at high polymer concentrations. Raman spectroscopy was used to give insight to the configurations of the bound molecules in concentrated solutions. In the vicinity of the critical gel concentration of E18B10 introduction of BSA can dramatically modify the phase diagram, inducing a gel-sol-gel transition. The overall picture of the interaction diagram of the E18B10-BSA reflects the shrinkage of the suspended particles due to destabilization of micelles induced by BSA and the gelator nature of the globular protein. SAXS and rheology were used to further characterize the structure and flow behavior of the polymer-protein hybrid gels and sols.

The adsorbed conformation of globular proteins at the air/water interface

External reflection FTIR spectroscopy and surface pressure measurements were used to compare conformational changes in the adsorbed structures of three globular proteins at the air/water interface. Of the three proteins studied, lysozyme, bovine serum albumin and P-lactoglobulin, lysozyme was unique in its behaviour. Lysozyme adsorption was slow, taking approximately 2.5 h to reach a surface pressure plateau (from a 0.07 mM solution), and led to significant structural change. The FTIR spectra revealed that lysozyme formed a highly networked adsorbed layer of unfolded protein with high antiparallel beta-sheet content and that these changes occurred rapidly (within 10 min). This non-native secondary structure is analogous to that of a 3D heat-set protein gel, suggesting that the adsorbed protein formed a highly networked interfacial layer. Albumin and P-lactoglobulin adsorbed rapidly (reaching a plateau within 10 min) and with little chance to their native secondary structure.

Complex formation of bovine serum albumin with a poly(ethylene glycol) lipid conjugate

In this work, we report the formation of complexes by self-assembly of bovine serum albumin (BSA) with a poly(ethylene glycol) lipid conjugate (PEG(2000)-PE) in phosphate saline buffer solution (pH 7.4). Three different sets of samples have been studied. The BSA concentration remained fixed (1, 0.01, or 0.001 wt % BSA) within each set of samples, while the PEG(2000)-PE concentration was varied. Dynamic light scattering (DLS), rheology, and small-angle X-ray scattering (SAXS) were used to study samples with 1 wt % BSA. DLS showed that BSA/PEG(2000)-PE aggregates have a size intermediate between a BSA monomer and a PEG(2000)-PE micelle. Rheology suggested that BSA/PEG(2000)-PE complexes might be surrounded by a relatively compact PEG-lipid shell, while SAXS results showed that depletion forces do not take an important role in the stabilization of the complexes. Samples containing 0.01 wt % BSA were studied by circular dichroism (CD) and ultraviolet fluorescence spectroscopy (UV). UV results showed that at low concentrations of PEG-lipid, PEG(2000)-PE binds to tryptophan (Trp) groups in BSA, while at high concentrations of PEG-lipid the Trp groups are exposed to water. CD results showed that changes in Trp environment take place with a minimal variation of the BSA secondary structure elements. Finally, samples containing 0.001 wt % BSA were studied by zeta-potential experiments. Results showed that steric interactions might play an important role in the stabilization of the BSA/PEG(2000)-PE complexes.

Prediction of novel and analogous folds using fragment assembly and fold recognition

A number of new and newly improved methods for predicting protein structure developed by the Jones–University College London group were used to make predictions for the CASP6 experiment. Structures were predicted with a combination of fold recognition methods (mGenTHREADER, nFOLD, and THREADER) and a substantially enhanced version of FRAGFOLD, our fragment assembly method. Attempts at automatic domain parsing were made using DomPred and DomSSEA, which are based on a secondary structure parsing algorithm and additionally for DomPred, a simple local sequence alignment scoring function. Disorder prediction was carried out using a new SVM-based version of DISOPRED. Attempts were also made at domain docking and “microdomain” folding in order to build complete chain models for some targets.

The DISOPRED server for the prediction of protein disorder

Dynamically disordered regions appear to be relatively abundant in eukaryotic proteomes. The DISOPRED server allows users to submit a protein sequence, and returns a probability estimate of each residue in the sequence being disordered. The results are sent in both plain text and graphical formats, and the server can also supply predictions of secondary structure to provide further structural information.

Improvement of the GenTHREADER method for genomic fold recognition

Motivation: In order to enhance genome annotation, the fully automatic fold recognition method GenTHREADER has been improved and benchmarked. The previous version of GenTHREADER consisted of a simple neural network which was trained to combine sequence alignment score, length information and energy potentials derived from threading into a single score representing the relationship between two proteins, as designated by CATH. The improved version incorporates PSI-BLAST searches, which have been jumpstarted with structural alignment profiles from FSSP, and now also makes use of PSIPRED predicted secondary structure and bi-directional scoring in order to calculate the final alignment score. Pairwise potentials and solvation potentials are calculated from the given sequence alignment which are then used as inputs to a multi-layer, feed-forward neural network, along with the alignment score, alignment length and sequence length. The neural network has also been expanded to accommodate the secondary structure element alignment (SSEA) score as an extra input and it is now trained to learn the FSSP Z-score as a measurement of similarity between two proteins. Results: The improvements made to GenTHREADER increase the number of remote homologues that can be detected with a low error rate, implying higher reliability of score, whilst also increasing the quality of the models produced. We find that up to five times as many true positives can be detected with low error rate per query. Total MaxSub score is doubled at low false positive rates using the improved method.

What are the baselines for protein fold recognition?

What constitutes a baseline level of success for protein fold recognition methods? As fold recognition benchmarks are often presented without any thought to the results that might be expected from a purely random set of predictions, an analysis of fold recognition baselines is long overdue. Given varying amounts of basic information about a protein—ranging from the length of the sequence to a knowledge of its secondary structure—to what extent can the fold be determined by intelligent guesswork? Can simple methods that make use of secondary structure information assign folds more accurately than purely random methods and could these methods be used to construct viable hierarchical classifications?

High yield production of a soluble bifidobacterial β-galactosidase (BbgIV) inE. coliDH5α with improved catalytic efficiency for the synthesis of prebiotic galactooligosaccharides

The bifidobacterial β-galactosidase (BbgIV) was produced in E. coli DH5α at 37 and 30 °C in a 5 L bioreactor under varied conditions of dissolved oxygen (dO2) and pH. The yield of soluble BbgIV was significantly (P < 0.05) increased once the dO2 dropped to 0–2% and remained at such low values during the exponential phase. Limited dO2 significantly (P < 0.05) increased the plasmid copy number and decreased the cells growth rate. Consequently, the BbgIV yield increased to its maximum (71–75 mg per g dry cell weight), which represented 20–25% of the total soluble proteins in the cells. In addition, the specific activity and catalytic efficiency of BbgIV were significantly (P < 0.05) enhanced under limited dO2 conditions. This was concomitant with a change in the enzyme secondary structure, suggesting a link between the enzyme structure and function. The knowledge generated from this work is very important for producing BbgIV as a biocatalyst for the development of a cost-effective process for the synthesis of prebiotic galactooligosaccharides from lactose.

Ruminal dry matter and nitrogen degradation in relation to condensed tannin and protein molecular structures in sainfoin (Onobrychis viciifolia) and lucerne (Medicago sativa).

Sainfoin is a temperate legume that contains condensed tannins (CT), i.e. polyphenols that are able to bind proteins and thus reduce protein degradation in the rumen. A reduction in protein degradation in the rumen can lead to a subsequent increase in amino acid flow to the small intestine. The effects of CT in the rumen and the intestine differ according to the amount and structure of CT and the nature of the protein molecular structure. The objective of the present study was to investigate the degradability in the rumen of three CT-containing sainfoin varieties and CT-free lucerne in relation to CT content and structure (mean degree of polymerization, proportion of prodelphinidins and cis-flavanol units) and protein structure (amide I and II bands, ratio of amide I-to-amide II, α-helix, β-sheet, ratio of α-helix-to-β-sheet). Protein molecular structures were identified using Fourier transform/infrared-attenuated total reflectance (FT/IR-ATR) spectroscopy. The in situ degradability of three sainfoin varieties (Ambra, Esparcette and Villahoz) was studied in 2008, during the first growth cycle at two harvest dates (P1 and P2, i.e. 5 May and 2 June, respectively) and at one date (P3) during the second growth cycle (2 June) and these were compared with a tannin-free legume, lucerne (Aubigny). Loss of dry matter (DMDeg) and nitrogen (NDeg) in polyester bags suspended in the rumen was measured using rumen-fistulated cows. The NDeg of lucerne compared with sainfoin was 0·80 v. 0·77 at P1, 0·78 v. 0·65 at P2 and 0·79 v. 0·70 at P3, respectively. NDeg was related to the rapidly disappearing fraction (‘a’) fraction (r=0·76), the rate of degradation (‘c’) (r=0·92), to the content (r=−0·81) and structure of CT. However, the relationship between NDeg and the slowly disappearing fraction (‘b’) was weak. There was a significant effect of date and species×date, for NDeg and ‘a’ fraction. The secondary protein structure varied with harvest date (species×date) and was correlated with the fraction ‘b’. Both tannin and protein structures influenced the NDeg degradation. CT content and structure were correlated to the ‘a’ fraction and to the ‘c’. Features of the protein molecular secondary structure were correlated to the ‘b’ fraction.

Interactions between lipid-free apolipoprotein-AI and a lipopeptide incorporating the RGDS cell adhesion motif

The interaction of a designed bioactive lipopeptide C16-GGGRGDS, comprising a hexadecyl lipid chain attached to a functional heptapeptide, with the lipid-free apoliprotein, Apo-AI, is examined. This apolipoprotein is a major component of high density lipoprotein and it is involved in lipid metabolism and may serve as a biomarker for cardiovascular disease and Alzheimers’ disease. We find via isothermal titration calorimetry that binding between the lipopeptide and Apo-AI occurs up to a saturation condition, just above equimolar for a 10.7 μM concentration of Apo-AI. A similar value is obtained from circular dichroism spectroscopy, which probes the reduction in α-helical secondary structure of Apo-AI upon addition of C16-GGGRGDS. Electron microscopy images show a persistence of fibrillar structures due to self-assembly of C16-GGGRGDS in mixtures with Apo-AI above the saturation binding condition. A small fraction of spheroidal or possibly “nanodisc” structures was observed. Small-angle X-ray scattering (SAXS) data for Apo-AI can be fitted using a published crystal structure of the Apo-AI dimer. The SAXS data for the lipopeptide/ Apo-AI mixtures above the saturation binding conditions can be fitted to the contribution from fibrillar structures coexisting with flat discs corresponding to Apo-AI/lipopeptide aggregates.

Studies on Anopheles (Kerteszia) homunculus Komp (Diptera: Culicidae)

The present findings suggest that Anopheles (Kerteszia) homunculus may comprise more than one species. The rDNA ITS2 sequence data corroborate the presence of An. homunculus l.s. in Mata Atlantica, southern Brazil, and suggest that specimens from Trinidad may belong to an unnamed morphologically similar species. There is a need for additional studies to establish the geographical distribution of An. homunculus l.s. in continental South America and in Trinidad, especially in southern Mata Atlantica, Brazil.

Expression and purification of human respiratory syncytial virus recombinant fusion protein

The Human Respiratory Syncytial Virus (HRSV) fusion protein (F) was expressed in Escherichia call BL21A using the pET28a vector at 37 degrees C. The protein was purified from the soluble fraction using affinity resin. The structural quality of the recombinant fusion protein and the estimation of its secondary structure were obtained by circular dichroism. Structural models of the fusion protein presented 46% of the helices in agreement with the spectra by circular dichroism analysis. There are only few studies that succeeded in expressing the HRSV fusion protein in bacteria. This is a report on human fusion protein expression in E. call and structure analysis, representing a step forward in the development of fusion protein F inhibitors and the production of antibodies. (c) 2008 Elsevier Inc. All rights reserved.

Putative outer membrane proteins of Leptospira interrogans stimulate human umbilical vein endothelial cells (HUVECS) and express during infection

Cell adhesion molecules (CAMs) are surface receptors present in eukaryotic cells that mediate cell-cell or cell-extracellular matrix interactions. Vascular endothelium stimulation in vitro that lead to the upregulation of CAMs was reported for the pathogenic spirochaetes, including rLIC10365 of Leptospira interrogans. In this study, we report the cloning of LIC10507, LIC10508, LIC10509 genes of L interrogans using Escherichia coli as a host system. The rational for selecting these sequences is due to their location in L. interrogans serovar Copenhageni genome that has a potential involvement in pathogenesis. The genes encode for predicted lipoproteins with no assigned functions. The purified recombinant proteins were capable to promote the upregulation of intercellular adhesion molecule 1 (ICAM-1) and E-selectin on monolayers of human umbilical vein endothelial cells (HUVECS). In addition, the coding sequences are expressed in the renal tubules of animal during bacterial experimental infection. The proteins are probably located at the outer membrane of the bacteria since they are detected in detergent-phase of L interrogans Triton X-114 extract. Altogether our data suggest a possible involvement of these proteins during bacterial infection and provide new insights into the role of this region in the pathogenesis of Leptospira. (C) 2008 Elsevier Ltd. All rights reserved.