Stochastic protein folding simulation in the three-dimensional HP-model

We present results from three-dimensional protein folding simulations in the HP-model on ten benchmark problems. The simulations are executed by a simulated annealing-based algorithm with a time-dependent cooling schedule. The neighbourhood relation is determined by the pull-move set. The results provide experimental evidence that the maximum depth D of local minima of the underlying energy landscape can be upper bounded by D < n(2/3). The local search procedure employs the stopping criterion (In/delta)(D/gamma) where m is an estimation of the average number of neighbouring conformations, gamma relates to the mean of non-zero differences of the objective function for neighbouring conformations, and 1-delta is the confidence that a minimum conformation has been found. The bound complies with the results obtained for the ten benchmark problems. (c) 2008 Elsevier Ltd. All rights reserved.

Population-based local search for protein folding simulation in the MJ energy model and cubic lattices

We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however. PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves. (C) 2009 Elsevier Ltd. All rights reserved.

Approximating the Set of Local Minima in Partial RNA Folding Landscapes

Motivation: We study a stochastic method for approximating the set of local minima in partial RNA folding landscapes associated with a bounded-distance neighbourhood of folding conformations. The conformations are limited to RNA secondary structures without pseudoknots. The method aims at exploring partial energy landscapes pL induced by folding simulations and their underlying neighbourhood relations. It combines an approximation of the number of local optima devised by Garnier and Kallel (2002) with a run-time estimation for identifying sets of local optima established by Reeves and Eremeev (2004).

Results: The method is tested on nine sequences of length between 50 nt and 400 nt, which allows us to compare the results with data generated by RNAsubopt and subsequent barrier tree calculations. On the nine sequences, the method captures on average 92% of local minima with settings designed for a target of 95%. The run-time of the heuristic can be estimated by O(n2D?ln?), where n is the sequence length, ? is the number of local minima in the partial landscape pL under consideration and D is the maximum number of steepest descent steps in attraction basins associated with pL.

General, PDB-based collective variables for protein folding

New, automated forms of data-analysis are required in order to understand the high-dimensional trajectories that are obtained from molecular dynamics simulations on proteins. Dimensionality reduction algorithms are particularly appealing in this regard as they allow one to construct unbiased, low-dimensional representations of the trajectory using only the information encoded in the trajectory. The downside of this approach is that different sets of coordinates are required for each different chemical systems under study precisely because the coordinates are constructed using information from the trajectory. In this paper we show how one can resolve this problem by using the sketch-map algorithm that we recently proposed to construct a low-dimensional representation of the structures contained in the protein data bank (PDB). We show that the resulting coordinates are as useful for analysing trajectory data as coordinates constructed using landmark configurations taken from the trajectory and that these coordinates can thus be used for understanding protein folding across a range of systems.

Enhancing the Manufacturing Knowledge of Undergraduate Engineering Students: A Case Study of a Design-Build-Test Challenge Involving Folding Bicycles

Many engineers currently in professional practice will have gained a degree level qualification which involved studying a curriculum heavy with mathematics and engineering science. While this knowledge is vital to the engineering design process so also is manufacturing knowledge, if the resulting designs are to be both technically and commercially viable.
The methodology advanced by the CDIO Initiative aims to improve engineering education by teaching in the context of Conceiving, Designing, Implementing and Operating products, processes or systems. A key element of this approach is the use of Design-Built-Test (DBT) projects as the core of an integrated curriculum. This approach facilitates the development of professional skills as well as the application of technical knowledge and skills developed in other parts of the degree programme. This approach also changes the role of lecturer to that of facilitator / coach in an active learning environment in which students gain concrete experiences that support their development.
The case study herein describes Mechanical Engineering undergraduate student involvement in the manufacture and assembly of concept and functional prototypes of a folding bicycle.

The conserved Helix C region in the superfamily of interferon-a/interleukin-10-related cytokines corresponds to a high-affinity binding site for the HSP70 chaperone DnaK.

HSP70 chaperones mediate protein folding by ATP-dependent interaction with short linear peptide segments that are exposed on unfolded proteins. The mode of action of the Escherichia coli homolog DnaK is representative of all HSP70 chaperones, including the endoplasmic reticulum variant BiP/GRP78. DnaK has been shown to be effective in assisting refolding of a wide variety of prokaryotic and eukaryotic proteins, including the -helical homodimeric secretory cytokine interferon- (IFN-). We screened solid-phase peptide libraries from human and mouse IFN- to identify DnaK-binding sites. Conserved DnaK-binding sites were identified in the N-terminal half of helix B and in the C-terminal half of helix C, both of which are located at the IFN- dimer interface. Soluble peptides derived from helices B and C bound DnaK with high affinity in competition assays. No DnaK-binding sites were found in the loops connecting the -helices. The helix C DnaK-binding site appears to be conserved in most members of the superfamily of interleukin (IL)-10-related cytokines that comprises, apart from IL-10 and IFN-, a series of recently discovered small secretory proteins, including IL-19, IL-20, IL-22/IL-TIF, IL-24/MDA-7 (melanoma differentiation-associated gene), IL-26/AK155, and a number of viral IL-10 homologs. These cytokines belong to a relatively small group of homodimeric proteins with highly interdigitated interfaces that exhibit the strongly hydrophobic character of the interior core of a single-chain folded domain. We propose that binding of DnaK to helix C in the superfamily of IL-10-related cytokines may constitute the hallmark of a novel conserved regulatory mechanism in which HSP70-like chaperones assist in the formation of a hydrophobic dimeric "folding" interface.

Conformational changes during proteolytic processing of a picornavirus capsid proteins

We have used synthetic peptide antibodies to probe conformational changes that occur during the cleavage cascade which generates the capsid proteins of a picornavirus. The initial translation product of 97 kDa, the precursor of all four structural proteins, is cleaved to form a 63 kDa fragment which, we show, has significantly different folding characteristics to both its larger parent and its products. We demonstrate that proteolytic cleavages as distant as 520 residues from epitopes confer sufficiently large conformational changes as to render them unrecognisable. To our knowledge, this is the first demonstration of this phenomenon in the picornavirus system.

Inhibition of Secretion of Interleukin (IL)-12/IL-23 Family Cytokines by 4-Trifluoromethyl-celecoxib Is Coupled to Degradation via the Endoplasmic Reticulum Stress Protein HERP

Interleukin-12 (IL-12), p80, and IL-23 are structurally related cytokines sharing a p40 subunit. We have recently demonstrated that celecoxib and its COX-2-independent analogue 4-trifluoromethyl-celecoxib (TFM-C) inhibit secretion but not transcription of IL-12 (p35/p40) and p80 (p40/p40). This is associated with a mechanism involving altered cytokine-chaperone interaction in the endoplasmic reticulum (ER). In the present study, we found that celecoxib and TFM-C also block secretion of IL-23 (p40/p19 heterodimers). Given the putative ER-centric mode of these compounds, we performed a comprehensive RTPCR analysis of 23 ER-resident chaperones/foldases and associated co-factors. This revealed that TFM-C induced 1.5-3-fold transcriptional up-regulation of calreticulin, GRP78, GRP94, GRP170, ERp72, ERp57, ERdj4, and ERp29. However, more significantly, a 7-fold up-regulation of homocysteine-inducible ER protein (HERP) was observed. HERP is part of a high molecular mass protein complex involved in ER-associated protein degradation (ERAD). Using co-immunoprecipitation assays, we show that TFM-C induces protein interaction of p80 and IL-23 with HERP. Both HERP siRNA knockdown and HERP overexpression coupled to cycloheximide chase assays revealed that HERP is necessary for degradation of intracellularly retained p80 by TFM-C. Thus, our data suggest that targeting cytokine folding in the ER by small molecule drugs could be therapeutically exploited to alleviate in appropriate inflammation in autoimmune conditions.

A topological algorithm for identification of structural domains of proteins

Background: Identification of the structural domains of proteins is important for our understanding of the organizational principles and mechanisms of protein folding, and for insights into protein function and evolution. Algorithmic methods of dissecting protein of known structure into domains developed so far are based on an examination of multiple geometrical, physical and topological features. Successful as many of these approaches are, they employ a lot of heuristics, and it is not clear whether they illuminate any deep underlying principles of protein domain organization. Other well-performing domain dissection methods rely on comparative sequence analysis. These methods are applicable to sequences with known and unknown structure alike, and their success highlights a fundamental principle of protein modularity, but this does not directly improve our understanding of protein spatial structure.

Structural and Molecular Biology of Type I Galactosemia: Disease-associated Mutations

Type I galactosemia results from reduced galactose 1-phosphate uridylyltransferase (GALT) activity. Signs of disease include damage to the eyes, brain, liver, and ovaries. However, the exact nature and severity of the pathology depends on the mutation(s) in the patient's genes and his/her environment. Considerable enzymological and structural knowledge has been accumulated and this provides a basis to explain, at a biochemical level, impairment in the enzyme in the more than 230 disease-associated variants, which have been described. The most common variant, Q188R, occurs close to the active site and the dimer interface. The substitution probably disrupts both UDP-sugar binding and homodimer stability. Other alterations, for example K285N, occur close to the surface of the enzyme and most likely affect the folding and stability of the enzyme. There are a number of unanswered questions in the field, which require resolution. These include the possibility that the main enzymes of galactose metabolism form a supramolecular complex and the need for a high resolution crystal structure of human GALT. (C) 2011 IUBMB IUBMB Life, 63(11): 949-954, 2011

Effects of Charge Location on the Absorptions and Lifetimes of Protonated Tyrosine Peptides in Vacuo

Nearby charges affect the electronic energy levels of chromophores, with the extent of the effect being determined by the magnitude of the charge and degree of charge-chromophore separation. The molecular configuration dictates the charge chromophore distance. Hence, in this study, we aim to assess how the location of the charge influences the absorption of a set of model protonated and diprotonated peptide ions, and whether spectral differences are large enough to be identified. The studied ions were the dipeptide YK, the tripeptide KYK (Y = tyrosine; K = lysine) and their complexes with 18-crown-6-ether (CE). The CE targets the ammonium group by forming internal ionic hydrogen bonds and limits the folding of the peptide. In the tripeptide, the distance between the chromophore and the backbone ammonium is enlarged relative to that in the dipeptide. Experiments were performed in an electrostatic ion storage ring using a tunable laser system, and action spectra based on lifetime measurements were obtained in the range from 210 to 310 nm. The spectra are all quite similar though there seems to be some changes in the absorption band between 210 and 250 nm, while in the lower energy band all ions had a maximum absorption at similar to 275 nm. Lifetimes after photoexcitation were found to shorten upon protonation and lengthen upon CE complexation, in accordance with the increased number of degrees of freedom and an increase in activation energies for dissociation as the mobile proton model is no longer operative.

Functional analysis of the C-terminal domain of the WbaP protein that mediates initiation of O antigen synthesis in Salmonella enterica

WbaP catalyzes the transfer of galactose-1-phosphate onto undecaprenyl phosphate (Und-P). The enzyme belongs to a large family of bacterial membrane proteins required for initiation of the synthesis of O antigen lipopolysaccharide and polysaccharide capsules. Previous work in our laboratory demonstrated that the last transmembrane helix and C-terminal tail region of WbaP (WbaP(CT)) are sufficient for enzymatic activity. Here, we demonstrate the cytoplasmic location of the WbaP C-terminal tail and show that WbaPCT domain N-terminally fused to thioredoxin (TrxA-WbaP(CT)) exhibits improved protein folding and enhanced transferase activity. Alanine replacement of highly conserved charged or polar amino acids identified seven critical residues for enzyme activity in vivo and in vitro. Four of these residues are located in regions predicted to be a-helical. These regions and their secondary structure predictions are conserved in distinct WbaP family members, suggesting they may contribute to form a conserved catalytic center.

First report of the use of a saxitoxin-protein conjugate to develop a DNA aptamer to a small molecule toxin

Saxitoxin (STX) is a low molecular weight neurotoxin mainly produced by certain marine dinoflagellates that, along with its family of similarly related paralytic shellfish toxins, may cause the potentially fatal intoxication known as paralytic shellfish poisoning. Illness and fatality rates are low due to the effective monitoring programs that determine when toxins exceed the established regulatory action level and effectuate shellfish harvesting closures accordingly. Such monitoring programs rely on the ability to rapidly screen large volumes of samples. Many of the screening assays currently available employ antibodies or live animals. This research focused on developing an analytical recognition element that would eliminate the challenges associated with the limited availability of antibodies and the use of animals. Here we report the discovery of a DNA aptamer that targets STX. Concentration-dependent and selective binding of the aptamer to STX was determined using a surface plasmon resonance sensor. Not only does this work represent the first reported aptamer to STX, but also the first aptamer to any marine biotoxin. A novel strategy of using a toxin-protein conjugate for DNA aptamer selection was successfully implemented to overcome the challenges associated with aptamer selection to small molecules. Taking advantage of such an approach could lead to increased diversity and accessibility of aptamers to low molecular weight toxins, which could then be incorporated as analytical recognition elements in diagnostic assays for foodborne toxin detection. The selected STX aptamer sequence is provided here, making it available to any investigator for use in assay development for the detection of STX.

Using sketch-map coordinates to analyze and bias molecular dynamics simulations

When examining complex problems, such as the folding of proteins, coarse grained descriptions of the system drive our investigation and help us to rationalize the results. Oftentimes collective variables (CVs), derived through some chemical intuition about the process of interest, serve this purpose. Because finding these CVs is the most difficult part of any investigation, we recently developed a dimensionality reduction algorithm, sketch-map, that can be used to build a low-dimensional map of a phase space of high-dimensionality. In this paper we discuss how these machine-generated CVs can be used to accelerate the exploration of phase space and to reconstruct free-energy landscapes. To do so, we develop a formalism in which high-dimensional configurations are no longer represented by low-dimensional position vectors. Instead, for each configuration we calculate a probability distribution, which has a domain that encompasses the entirety of the low-dimensional space. To construct a biasing potential, we exploit an analogy with metadynamics and use the trajectory to adaptively construct a repulsive, history-dependent bias from the distributions that correspond to the previously visited configurations. This potential forces the system to explore more of phase space by making it desirable to adopt configurations whose distributions do not overlap with the bias. We apply this algorithm to a small model protein and succeed in reproducing the free-energy surface that we obtain from a parallel tempering calculation.