Protease-sensitive synthetic prions.

Prions arise when the cellular prion protein (PrP(C)) undergoes a self-propagating conformational change; the resulting infectious conformer is designated PrP(Sc). Frequently, PrP(Sc) is protease-resistant but protease-sensitive (s) prions have been isolated in humans and other animals. We report here that protease-sensitive, synthetic prions were generated in vitro during polymerization of recombinant (rec) PrP into amyloid fibers. In 22 independent experiments, recPrP amyloid preparations, but not recPrP monomers or oligomers, transmitted disease to transgenic mice (n = 164), denoted Tg9949 mice, that overexpress N-terminally truncated PrP. Tg9949 control mice (n = 174) did not spontaneously generate prions although they were prone to late-onset spontaneous neurological dysfunction. When synthetic prion isolates from infected Tg9949 mice were serially transmitted in the same line of mice, they exhibited sPrP(Sc) and caused neurodegeneration. Interestingly, these protease-sensitive prions did not shorten the life span of Tg9949 mice despite causing extensive neurodegeneration. We inoculated three synthetic prion isolates into Tg4053 mice that overexpress full-length PrP; Tg4053 mice are not prone to developing spontaneous neurological dysfunction. The synthetic prion isolates caused disease in 600-750 days in Tg4053 mice, which exhibited sPrP(Sc). These novel synthetic prions demonstrate that conformational changes in wild-type PrP can produce mouse prions composed exclusively of sPrP(Sc).

Natural and synthetic prion structure from X-ray fiber diffraction.

A conformational isoform of the mammalian prion protein (PrP(Sc)) is the sole component of the infectious pathogen that causes the prion diseases. We have obtained X-ray fiber diffraction patterns from infectious prions that show cross-beta diffraction: meridional intensity at 4.8 A resolution, indicating the presence of beta strands running approximately at right angles to the filament axis and characteristic of amyloid structure. Some of the patterns also indicated the presence of a repeating unit along the fiber axis, corresponding to four beta-strands. We found that recombinant (rec) PrP amyloid differs substantially from highly infectious brain-derived prions, both in structure as demonstrated by the diffraction data, and in heterogeneity as shown by electron microscopy. In addition to the strong 4.8 A meridional reflection, the recPrP amyloid diffraction is characterized by strong equatorial intensity at approximately 10.5 A, absent from brain-derived prions, and indicating the presence of stacked beta-sheets. Synthetic prions recovered from transgenic mice inoculated with recPrP amyloid displayed structural characteristics and homogeneity similar to those of naturally occurring prions. The relationship between the structural differences and prion infectivity is uncertain, but might be explained by any of several hypotheses: only a minority of recPrP amyloid possesses a replication-competent conformation, the majority of recPrP amyloid has to undergo a conformational maturation to acquire replication competency, or inhibitory forms of recPrP amyloid interfere with replication during the initial transmission.

Therapeutic approaches to protein-misfolding diseases.

Several sporadic and genetic diseases are caused by protein misfolding. These include cystic fibrosis and other devastating diseases of childhood as well as Alzheimer's, Parkinson's and other debilitating maladies of the elderly. A unified view of the molecular and cellular pathogenesis of these conditions has led to the search for chemical chaperones that can slow, arrest or revert disease progression. Molecules are now emerging that link our biophysical insights with our therapeutic aspirations.

Cell division modulates prion accumulation in cultured cells.

The phenotypic effect of prions on host cells is influenced by the physical properties of the prion strain and its level of accumulation. In mammalian cell cultures, prion accumulation is determined by the interplay between de novo prion formation, catabolism, cell division, and horizontal cell-to-cell transmission. Understanding this dynamic enables the analytical modeling of protein-based heritability and infectivity. Here, we quantitatively measured these competing effects in a subline of neuroblastoma (N2a) cells and propose a concordant reaction mechanism to explain the kinetics of prion propagation. Our results show that cell division leads to a predictable reduction in steady-state prion levels but not to complete clearance. Scrapie-infected N2a cells were capable of accumulating different steady-state levels of prions, dictated partly by the rate of cell division. We also show that prions in this subline of N2a cells are transmitted primarily from mother to daughter cells, rather than horizontal cell-to-cell transmission. We quantitatively modeled our kinetic results based on a mechanism that assumes a subpopulation of prions is capable of self-catalysis, and the levels of this subpopulation reach saturation in fully infected cells. Our results suggest that the apparent effectiveness of antiprion compounds in culture may be strongly influenced by the growth phase of the target cells.

Molecular modeling of the misfolded insulin subunit and amyloid fibril.

Insulin, a small hormone protein comprising 51 residues in two disulfide-linked polypeptide chains, adopts a predominantly alpha-helical conformation in its native state. It readily undergoes protein misfolding and aggregates into amyloid fibrils under a variety of conditions. Insulin is a unique model system in which to study protein fibrillization, since its three disulfide bridges are retained in the fibrillar state and thus limit the conformational space available to the polypeptide chains during misfolding and fibrillization. Taking into account this unique conformational restriction, we modeled possible monomeric subunits of the insulin amyloid fibrils using beta-solenoid folds, namely, the beta-helix and beta-roll. Both models agreed with currently available biophysical data. We performed molecular dynamics simulations, which allowed some limited insights into the relative structural stability, suggesting that the beta-roll subunit model may be more stable than the beta-helix subunit model. We also constructed beta-solenoid-based insulin fibril models and conducted fiber diffraction simulation to identify plausible fibril architectures of insulin amyloid. A comparison of simulated fiber diffraction patterns of the fibril models to the experimental insulin x-ray fiber diffraction data suggests that the model fibers composed of six twisted beta-roll protofilaments provide the most reasonable fit to available experimental diffraction patterns and previous biophysical studies.

Analysis of the sequence and structural features of the left-handed beta-helical fold.

The left-handed parallel beta-helix (LbetaH) is a structurally repetitive, highly regular, and symmetrical fold formed by coiling of elongated beta-sheets into helical "rungs." This canonical fold has recently received interest as a possible solution to the fibril structure of amyloid and as a building block of self-assembled nanotubular structures. In light of this interest, we aimed to understand the structural requirements of the LbetaH fold. We first sought to determine the sequence characteristics of the repeats by analyzing known structures to identify positional preferences of specific residues types. We then used molecular dynamics simulations to demonstrate the stabilizing effect of successive rungs and the hydrophobic core of the LbetaH. We show that a two-rung structure is the minimally stable LbetaH structure. In addition, we defined the structure-based sequence preference of the LbetaH and undertook a genome-wide sequence search to determine the prevalence of this unique protein fold. This profile-based LbetaH search algorithm predicted a large fraction of LbetaH proteins from microbial origins. However, the relative number of predicted LbetaH proteins per specie was approximately equal across the genomes from prokaryotes to eukaryotes.

Automatic extraction of protein point mutations using a graph bigram association.

Protein point mutations are an essential component of the evolutionary and experimental analysis of protein structure and function. While many manually curated databases attempt to index point mutations, most experimentally generated point mutations and the biological impacts of the changes are described in the peer-reviewed published literature. We describe an application, Mutation GraB (Graph Bigram), that identifies, extracts, and verifies point mutations from biomedical literature. The principal problem of point mutation extraction is to link the point mutation with its associated protein and organism of origin. Our algorithm uses a graph-based bigram traversal to identify these relevant associations and exploits the Swiss-Prot protein database to verify this information. The graph bigram method is different from other models for point mutation extraction in that it incorporates frequency and positional data of all terms in an article to drive the point mutation-protein association. Our method was tested on 589 articles describing point mutations from the G protein-coupled receptor (GPCR), tyrosine kinase, and ion channel protein families. We evaluated our graph bigram metric against a word-proximity metric for term association on datasets of full-text literature in these three different protein families. Our testing shows that the graph bigram metric achieves a higher F-measure for the GPCRs (0.79 versus 0.76), protein tyrosine kinases (0.72 versus 0.69), and ion channel transporters (0.76 versus 0.74). Importantly, in situations where more than one protein can be assigned to a point mutation and disambiguation is required, the graph bigram metric achieves a precision of 0.84 compared with the word distance metric precision of 0.73. We believe the graph bigram search metric to be a significant improvement over previous search metrics for point mutation extraction and to be applicable to text-mining application requiring the association of words.

Directed evolution of an anti-prion protein scFv fragment to an affinity of 1 pM and its structural interpretation.

Bovine spongiform encephalopathy (BSE) is a fatal neurodegenerative prion disease affecting cattle that is transmissible to humans, manifesting as a variant of Creutzfeldt-Jakob disease (vCJD) likely following the consumption of meat contaminated with BSE prions. High-affinity antibodies are a prerequisite for the development of simple, highly sensitive and non-invasive diagnostic tests that are able to detect even small amounts of the disease-associated PrP conformer (PrP(Sc)). We describe here the affinity maturation of a single-chain Fv antibody fragment with a binding affinity of 1 pM to a peptide derived from the unstructured region of bovine PrP (BoPrP (90-105)). This is the tightest peptide-binding antibody reported to date and may find useful application in diagnostics, especially when PrP(Sc) is pretreated by denaturation and/or proteolysis for peptide-like presentation. Several rounds of directed evolution and off-rate selection with ribosome display were performed using an antibody library generated from a single PrP binder with error-prone PCR and DNA-shuffling. As the correct determinations of affinities in this range are not straightforward, competition biosensor techniques and KinExA methods were both applied and compared. Structural interpretation of the affinity improvement was performed based on the crystal structure of the original prion binder in complex with the BoPrP (95-104) peptide by modeling the corresponding mutations.

Structure-activity relationship study of prion inhibition by 2-aminopyridine-3,5-dicarbonitrile-based compounds: parallel synthesis, bioactivity, and in vitro pharmacokinetics.

2-Aminopyridine-3,5-dicarbonitrile compounds were previously identified as mimetics of dominant-negative prion protein mutants and inhibit prion replication in cultured cells. Here, we report findings from a comprehensive structure-activity relationship study of the 6-aminopyridine-3,5-dicarbonitrile scaffold. We identify compounds with significantly improved bioactivity (approximately 40-fold) against replication of the infectious prion isoform (PrPSc) and suitable pharmacokinetic profiles to warrant evaluation in animal models of prion disease.

Electron crystallography of the scrapie prion protein complexed with heavy metals.

The insolubility of the disease-causing isoform of the prion protein (PrP(Sc)) has prevented studies of its three-dimensional structure at atomic resolution. Electron crystallography of two-dimensional crystals of N-terminally truncated PrP(Sc) (PrP 27-30) and a miniprion (PrP(Sc)106) provided the first insights at intermediate resolution on the molecular architecture of the prion. Here, we report on the structure of PrP 27-30 and PrP(Sc)106 negatively stained with heavy metals. The interactions of the heavy metals with the crystal lattice were governed by tertiary and quaternary structural elements of the protein as well as the charge and size of the heavy metal salts. Staining with molybdate anions revealed three prominent densities near the center of the trimer that forms the unit cell, coinciding with the location of the beta-helix that was proposed for the structure of PrP(Sc). Differential staining also confirmed the location of the internal deletion of PrP(Sc)106 at or near these densities.

Time-controlled transcardiac perfusion cross-linking for the study of protein interactions in complex tissues.

Because of their sensitivity to solubilizing detergents, membrane protein assemblies are difficult to study. We describe a protocol that covalently conserves protein interactions through time-controlled transcardiac perfusion cross-linking (tcTPC) before disruption of tissue integrity. To validate tcTPC for identifying protein-protein interactions, we established that tcTPC allowed stringent immunoaffinity purification of the gamma-secretase complex in high salt concentrations and detergents and was compatible with mass spectrometric identification of cross-linked aph-1, presenilin-1 and nicastrin. We then applied tcTPC to identify more than 20 proteins residing in the vicinity of the cellular prion protein (PrPC), suggesting that PrP is embedded in specialized membrane regions with a subset of molecules that, like PrP, use a glycosylphosphatidylinositol anchor for membrane attachment. Many of these proteins have been implicated in cell adhesion/neuritic outgrowth, and harbor immunoglobulin C2 and fibronectin type III-like motifs.

GPCRDB information system for G protein-coupled receptors.

The GPCRDB is a molecular class-specific information system that collects, combines, validates and disseminates heterogeneous data on G protein-coupled receptors (GPCRs). The database stores data on sequences, ligand binding constants and mutations. The system also provides computationally derived data such as sequence alignments, homology models, and a series of query and visualization tools. The GPCRDB is updated automatically once every 4-5 months and is freely accessible at http://www.gpcr.org/7tm/.

New local potential useful for genome annotation and 3D modeling.

A new potential energy function representing the conformational preferences of sequentially local regions of a protein backbone is presented. This potential is derived from secondary structure probabilities such as those produced by neural network-based prediction methods. The potential is applied to the problem of remote homolog identification, in combination with a distance-dependent inter-residue potential and position-based scoring matrices. This fold recognition jury is implemented in a Java application called JThread. These methods are benchmarked on several test sets, including one released entirely after development and parameterization of JThread. In benchmark tests to identify known folds structurally similar to (but not identical with) the native structure of a sequence, JThread performs significantly better than PSI-BLAST, with 10% more structures identified correctly as the most likely structural match in a fold library, and 20% more structures correctly narrowed down to a set of five possible candidates. JThread also improves the average sequence alignment accuracy significantly, from 53% to 62% of residues aligned correctly. Reliable fold assignments and alignments are identified, making the method useful for genome annotation. JThread is applied to predicted open reading frames (ORFs) from the genomes of Mycoplasma genitalium and Drosophila melanogaster, identifying 20 new structural annotations in the former and 801 in the latter.

Co-evolutionary analysis reveals insights into protein-protein interactions.

Protein-protein interactions play crucial roles in biological processes. Experimental methods have been developed to survey the proteome for interacting partners and some computational approaches have been developed to extend the impact of these experimental methods. Computational methods are routinely applied to newly discovered genes to infer protein function and plausible protein-protein interactions. Here, we develop and extend a quantitative method that identifies interacting proteins based upon the correlated behavior of the evolutionary histories of protein ligands and their receptors. We have studied six families of ligand-receptor pairs including: the syntaxin/Unc-18 family, the GPCR/G-alpha's, the TGF-beta/TGF-beta receptor system, the immunity/colicin domain collection from bacteria, the chemokine/chemokine receptors, and the VEGF/VEGF receptor family. For correlation scores above a defined threshold, we were able to find an average of 79% of all known binding partners. We then applied this method to find plausible binding partners for proteins with uncharacterized binding specificities in the syntaxin/Unc-18 protein and TGF-beta/TGF-beta receptor families. Analysis of the results shows that co-evolutionary analysis of interacting protein families can reduce the search space for identifying binding partners by not only finding binding partners for uncharacterized proteins but also recognizing potentially new binding partners for previously characterized proteins. We believe that correlated evolutionary histories provide a route to exploit the wealth of whole genome sequences and recent systematic proteomic results to extend the impact of these studies and focus experimental efforts to categorize physiologically or pathologically relevant protein-protein interactions.

Differences between the prion protein and its homolog Doppel: a partially structured state with implications for scrapie formation.

The key event in the pathogenesis of prion diseases is a conformational change in the prion protein (PrP). Models for conversion of PrP(C) into PrP(Sc) typically implicate an, as yet, unidentified intermediate. In an attempt to identify such an intermediate, we used native-state hydrogen exchange monitored with NMR. Although we were unable to detect an intermediate directly, we observed substantial protection above that expected based upon measurements of the global stability of PrP (>2 kcal mol(-1) super protection). This super protection implicates either structure in the denatured state or the presence of an intermediate. Similar experiments with Doppel, a homolog of PrP that does not form infectious prions, failed to demonstrate such super protection. This suggests that the partially structured state of PrP encompassing portions of the B and C helices, may be a significant factor in the ability of PrP to convert from PrP(C) to PrP(Sc).

Toward the synthesis of artificial proteins: the discovery of an amphiphilic helical peptoid assembly.

While nature exploits folded biopolymers to achieve molecular recognition and catalysis, comparable abiological heteropolymer systems have been difficult to create. We synthesized and identified abiological peptoid heteroploymers capable of binding a dye. Using combinatorial synthesis, we constructed a library of 3400 amphiphilic 15-mer peptoids on an ultra-high-capacity beaded support. Individual macrobeads, each containing a single peptoid sequence, were arrayed into plates, cleaved, and screened in aqueous solution to locate dye binding heteropolymer assemblies. Resynthesis and characterization demonstrated the formation of defined helical assemblies as judged by size-exclusion chromatography, circular dichroism, and analytical ultracentrifugation. Inspired by nature's process of sequence variation and natural selection, we identified rare abiological sequence-specific heteropolymers that begin to mimic the structure and functional properties of their biological counterparts.

Structure-activity relationships of the antimalarial agent artemisinin. 6. The development of predictive in vitro potency models using CoMFA and HQSAR methodologies.

Artemisinin (1) is a unique sesquiterpene peroxide occurring as a constituent of Artemisia annua L. Because of the effectiveness of Artemisinin in the treatment of drug-resistant Plasmodium falciparum and its rapid clearance of cerebral malaria, development of clinically useful semisynthetic drugs for severe and complicated malaria (artemether, artesunate) was prompt. However, recent reports of fatal neurotoxicity in animals with dihydroartemisinin derivatives such as artemether have spawned a renewed effort to develop nontoxic analogues of artemisinin. In our effort to develop more potent, less neurotoxic agents for the oral treatment of drug-resistant malaria, we utilized comparative molecular field analysis (CoMFA) and hologram QSAR (HQSAR), beginning with a series of 211 artemisinin analogues with known in vitro antimalarial activity. CoMFA models were based on two conformational hypotheses: (a) that the X-ray structure of artemisinin represents the bioactive shape of the molecule or (b) that the hemin-docked conformation is the bioactive form of the drug. In addition, we examined the effect of inclusion or exclusion of racemates in the partial least squares (pls) analysis. Databases derived from the original 211 were split into chiral (n = 157), achiral (n = 34), and mixed databases (n = 191) after leaving out a test set of 20 compounds. HQSAR and CoMFA models were compared in terms of their potential to generate robust QSAR models. The r(2) and q(2) (cross-validated r(2)) were used to assess the statistical quality of our models. Another statistical parameter, the ratio of the standard error to the activity range (s/AR), was also generated. CoMFA and HQSAR models were developed having statistically excellent properties, which also possessed good predictive ability for test set compounds. The best model was obtained when racemates were excluded from QSAR analysis. Thus, CoMFA of the n = 157 database gave excellent predictions with outstanding statistical properties. HQSAR did an outstanding job in statistical analysis and also handled predictions well.

Synthesis and structure-activity relationship study of potent trypanocidal thio semicarbazone inhibitors of the trypanosomal cysteine protease cruzain.

American trypanosomiasis, or Chagas' disease, is the leading cause of heart disease in Latin America. Currently there is an urgent need to develop antitrypanosomal therapy due to the toxicity of existing agents and emerging drug resistance. A novel series of potent thio semicarbazone small-molecule inhibitors of the Trypanosoma cruzi cysteine protease cruzain have been identified. Some of these inhibitors have been shown to be trypanocidal. We initially discovered that 3'-bromopropiophenone thio semicarbazone (1i) inhibited cruzain and could cure mammalian cell cultures infected with T. cruzi. 3'-Bromopropiophenone thio semicarbazone showed no toxicity for mammalian cells at concentrations that were trypanocidal. Following this lead, more than 100 compounds were designed and synthesized. A specific structure-activity relationship (SAR) was established, and many potent analogues with IC(50) values in the low nanomolar range were identified. Eight additional analogues were trypanocidal in a cell culture assay, and this indicates that aryl thio semicarbazone is a productive scaffold for killing the parasites. Kinetic studies show that these are time-dependent inhibitors. Molecular modeling studies of the enzyme-inhibitor complex have led to a proposed mechanism of interaction as well as insight into the SAR of the thio semicarbazone series. The nonpeptide nature of this series, small size, and extremely low cost of production suggest this is a promising direction for the development of new antitrypanosome chemotherapy.

Screening of acyl hydrazide proteinase inhibitors for antiparasitic activity against Trypanosoma brucei.

The major cysteine proteinase (brucipain) of Trypanosoma brucei is a target for chemotherapy of African Sleeping Sickness. We have screened a non-peptidyl acyl hydrazide proteinase inhibitor library of 500 compounds for inhibition of brucipain. Those 21 compounds with IC(50) values of <40 microM were tested for efficacy against bloodstream forms of T. brucei in cell culture. Eight acyl hydrazides showed 50% or more inhibition of trypanosome replication at <1 microM. The trypanocidal acitivity of the most effective compounds was comparable with those of the commercial antitrypanosomal drugs suramin and diminazene aceturate. However, these acyl hydrazides exhibited varying cytotoxicity towards human HL-60 cells and therefore, only less favourable selectivity indices compared with the commercially available drugs. Nevertheless, the data support the potential of acyl hydrazides as antitrypanosomal chemotherapeutic agents for treatment of sleeping sickness.