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.

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.

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.

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.

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.

Arsenate sensitivity in ericoid and ectomycorrhizal fungi

Isolates of the endomycorrhizal fungus Hymenoscyphus ericae and the ectomycorrhizal fungus Hebeloma crustuliniforme from soils uncontaminated with AsO₄/^3-, were compared with regard to their sensitivity to AsO₄/^3- in solution culture. When grown in liquid media amended with a range of AsO₄/^3- concentrations, H. ericae demonstrated reduced sensitivity to AsO₄/^3- compared to H. crustuliniforme. The concentrations causing 50% inhibition of growth (EC50) were 1.33 mol/m³ and 0.33 mol/m³, respectively, for H. ericae and H. crustuliniforme. The compound AsO₄/^3- is a PO₄/^3- analogue for the plasmalemma PO₄/^3- transporter. The presence of PO₄/^3- in the media at high concentrations ameliorated the toxic effects of AsO₄/^3- in both the ericoid and ectomycorrhizal fungi. This could be due to both suppression of the PO₄/^3- transporter under high phosporus status and competition of PO₄/^3- with AsO₄/^3- for the transport protein. The kinetics of AsO₄/^3- influx in H. ericae and H. crustuliniforme were also investigated. Hymenoscyphus ericae demonstrated a high K(m) value, 0.071 mol/m³, consistent with values obtained for AsO₄/^3- -tolerant plants. We suggest that the high K(m) value may be a mechanism used by H. ericae to express reduced AsO₄/^3- sensitivity. The ecological significance of this reduced sensitivity is also discussed.

Role of the AP2 beta-appendage hub in recruiting partners for clathrin-coated vesicle assembly

Adaptor protein complex 2 alpha and beta-appendage domains act as hubs for the assembly of accessory protein networks involved in clathrin-coated vesicle formation. We identify a large repertoire of beta-appendage interactors by mass spectrometry. These interact with two distinct ligand interaction sites on the beta-appendage (the "top" and "side" sites) that bind motifs distinct from those previously identified on the alpha-appendage. We solved the structure of the beta-appendage with a peptide from the accessory protein Eps15 bound to the side site and with a peptide from the accessory cargo adaptor beta-arrestin bound to the top site. We show that accessory proteins can bind simultaneously to multiple appendages, allowing these to cooperate in enhancing ligand avidities that appear to be irreversible in vitro. We now propose that clathrin, which interacts with the beta-appendage, achieves ligand displacement in vivo by self-polymerisation as the coated pit matures. This changes the interaction environment from liquid-phase, affinity-driven interactions, to interactions driven by solid-phase stability ("matricity"). Accessory proteins that interact solely with the appendages are thereby displaced to areas of the coated pit where clathrin has not yet polymerised. However, proteins such as beta-arrestin (non-visual arrestin) and autosomal recessive hypercholesterolemia protein, which have direct clathrin interactions, will remain in the coated pits with their interacting receptors.

Studying N-linked glycosylation of receptor tyrosine kinases

Metabolic alterations have been identified as a frequent event in cancer. This is often associated with increased flux through glycolysis, and also a secondary pathway to glycolysis, hexosamine biosynthetic pathway (HBP). HBP provides substrate for N-linked glycosylation, which occurs in the endoplasmic reticulum and the Golgi apparatus. N-linked glycosylation supports protein folding and correct sorting of proteins to plasma membrane and secretion. This process generates complex glycoforms, which can be recognized by other proteins and glycosylation of receptor tyrosine kinases (RTK) can also regulate their plasma-membrane retention time. Of special interest for experimental biologists, plants produce proteins, termed lectins, which bind with high specificity to glyco-conjugates. For the purposes of molecular biology, plant lectins can be conjugated to different moieties, such as agarose beads, which enable precipitation of specifically glycosylated proteins. In this chapter, we describe in detail how to perform pull-down experiments with commercially available lectins to identify changes in the glycosylation of RTKs.

Choline Kinase Alpha as an Androgen Receptor Chaperone and Prostate Cancer Therapeutic Target

BACKGROUND: The androgen receptor (AR) is a major drug target in prostate cancer (PCa). We profiled the AR-regulated kinome to identify clinically relevant and druggable effectors of AR signaling.

METHODS: Using genome-wide approaches, we interrogated all AR regulated kinases. Among these, choline kinase alpha (CHKA) expression was evaluated in benign (n = 195), prostatic intraepithelial neoplasia (PIN) (n = 153) and prostate cancer (PCa) lesions (n = 359). We interrogated how CHKA regulates AR signaling using biochemical assays and investigated androgen regulation of CHKA expression in men with PCa, both untreated (n = 20) and treated with an androgen biosynthesis inhibitor degarelix (n = 27). We studied the effect of CHKA inhibition on the PCa transcriptome using RNA sequencing and tested the effect of CHKA inhibition on cell growth, clonogenic survival and invasion. Tumor xenografts (n = 6 per group) were generated in mice using genetically engineered prostate cancer cells with inducible CHKA knockdown. Data were analyzed with χ(2) tests, Cox regression analysis, and Kaplan-Meier methods. All statistical tests were two-sided.

RESULTS: CHKA expression was shown to be androgen regulated in cell lines, xenografts, and human tissue (log fold change from 6.75 to 6.59, P = .002) and was positively associated with tumor stage. CHKA binds directly to the ligand-binding domain (LBD) of AR, enhancing its stability. As such, CHKA is the first kinase identified as an AR chaperone. Inhibition of CHKA repressed the AR transcriptional program including pathways enriched for regulation of protein folding, decreased AR protein levels, and inhibited the growth of PCa cell lines, human PCa explants, and tumor xenografts.

CONCLUSIONS: CHKA can act as an AR chaperone, providing, to our knowledge, the first evidence for kinases as molecular chaperones, making CHKA both a marker of tumor progression and a potential therapeutic target for PCa.

Biological phosphorus removal during high-rate, low-temperature, anaerobic digestion of wastewater.

We report, for the first time, extensive biologically-mediated phosphate removal from wastewater during high-rate anaerobic digestion (AD). A hybrid sludge bed/fixed-film (packed pumice stone) reactor was employed for low-temperature (12°C) anaerobic treatment of synthetic sewage wastewater. Successful phosphate removal from the wastewater (up to 78% of influent phosphate) was observed, mediated by biofilms in the reactor. Scanning electron microscopy and energy dispersive X-ray analysis revealed the accumulation of elemental phosphorus (~2%) within the sludge bed and fixed-film biofilms. 4’, 6-diamidino-2-phenylindole (DAPI) staining indicated phosphorus accumulation was biological in nature and mediated through the formation of intracellular inorganic polyphosphate (polyP) granules within these biofilms. DAPI staining further indicated that polyP accumulation was rarely associated with free cells. Efficient and consistent chemical oxygen demand (COD) removal was recorded, throughout the 732-day trial, at applied organic loading rates between 0.4-1.5 kg COD m-3 d-1 and hydraulic retention times of 8-24 hours, while phosphate removal efficiency ranged from 28-78% on average per phase. Analysis of protein hydrolysis kinetics and the methanogenic activity profiles of the biomass revealed the development, at 12˚C, of active hydrolytic and methanogenic populations. Temporal microbial changes were monitored using Illumina Miseq analysis of bacterial and archaeal 16S rRNA gene sequences. The dominant bacterial phyla present in the biomass at the conclusion of the trial were the Proteobacteria and Firmicutes and the dominant archaeal genus was Methanosaeta. Trichococcus and Flavobacterium populations, previously associated with low temperature protein degradation, developed in the reactor biomass. The presence of previously characterised polyphosphate accumulating organisms (PAOs) such as Rhodocyclus, Chromatiales, Actinobacter and Acinetobacter was recorded at low numbers. However, it is unknown as yet if these were responsible for the luxury polyP uptake observed in this system. The possibility of efficient phosphate removal and recovery from wastewater during AD would represent a major advance in the scope for widespread application of anaerobic wastewater treatment technologies.

Potential strategies for the eradication of multi-drug resistant Gram-negative bacterial infections

Antimicrobial resistance is one of the leading threats to society. The increasing burden of multidrug-resistant Gram-negative infection is particularly concerning as such bacteria are demonstrating resistance to nearly all currently licensed therapies. Various strategies have been hypothesized to treat multidrug-resistant Gram-negative infections including: targeting the Gram-negative outer membrane; neutralization of lipopolysaccharide; inhibition of bacterial efflux pumps and prevention of protein folding. Silver and silver nanoparticles, fusogenic liposomes and nanotubes are potential strategies for extending the activity of licensed, Gram-positive selective, antibiotics to Gram-negatives. This may serve as a strategy to fill the current void in pharmaceutical development in the short term. This review outlines the most promising strategies that could be implemented to solve the threat of multidrug-resistant Gram-negative infections

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

Structure and kinetics of phosphonopyruvate hydrolase from Variovorax sp. Pal2: New insight into the divergence of catalysis within the PEP mutase/isocitrate lyase superfamily .

Phosphonopyruvate (P-pyr) hydrolase (PPH), a member of the phosphoenolpyruvate (PEP) mutase/isocitrate lyase (PEPM/ICL) superfamily, hydrolyzes P-pyr and shares the highest sequence identity and functional similarity with PEPM. Recombinant PPH from Variovorax sp. Pal2 was expressed in Escherichia coli and purified to homogeneity. Analytical gel filtration indicated that the protein exists in solution predominantly as a tetramer. The PPH pH rate profile indicates maximal activity over a broad pH range.The steady-state kinetic constants determined for a rapid equilibrium ordered kinetic mechanism with Mg+2 binding first (Kd =140 ± 40 M), are kcat = 105 ± 2 s-1 and P-pyr Km = 5 ± 1 M. PEP (slow substrate kcat = 2 × 10-4 s-1), oxalate, and sulfopyruvate are competitive inhibitors with Ki values of 2.0 ± 0.1 mM, 17 ± 1 M, and 210 ± 10 M, respectively. Three PPH crystal structures have been determined, that of a ligand-free enzyme, the enzyme bound to Mg2+ and oxalate (inhibitor), and the enzyme bound to Mg2+ and P-pyr (substrate). The complex with the inhibitor was obtained by cocrystallization, whereas that with the substrate was obtained by briefly soaking crystals of the ligand-free enzyme with P-pyr prior to flash cooling. The PPH structure resembles that of the other members of the PEPM/ICL superfamily and is most similar to the functionally related enzyme, PEPM. Each monomer of the dimer of dimers exhibits an (/)8 barrel fold with the eighth helix swapped between two molecules of the dimer. Both P-pyr and oxalate are anchored to the active site by Mg2+. The loop capping the active site is disordered in all three structures, in contrast to PEPM, where the equivalent loop adopts an open or disordered conformation in the unbound state but sequesters the inhibitor from solvent in the bound state. Crystal packing may have favored the open conformation of PPH even when the enzyme was cocrystallized with the oxalate inhibitor. Structure alignment of PPH with other superfamily members revealed two pairs of invariant or conservatively replaced residues that anchor the flexible gating loop. The proposed PPH catalytic mechanism is analogous to that of PEPM but includes activation of a water nucleophile with the loop Thr118 residue.