Microarray analyses reveal strain-specific antibody responses to Plasmodium falciparum apical membrane antigen 1 variants following natural infection and vaccination.

Vaccines based on Plasmodium falciparum apical membrane antigen 1 (AMA1) have failed due to extensive polymorphism in AMA1. To assess the strain-specificity of antibody responses to malaria infection and AMA1 vaccination, we designed protein and peptide microarrays representing hundreds of unique AMA1 variants. Following clinical malaria episodes, children had short-lived, sequence-independent increases in average whole-protein seroreactivity, as well as strain-specific responses to peptides representing diverse epitopes. Vaccination resulted in dramatically increased seroreactivity to all 263 AMA1 whole-protein variants. High-density peptide analysis revealed that vaccinated children had increases in seroreactivity to four distinct epitopes that exceeded responses to natural infection. A single amino acid change was critical to seroreactivity to peptides in a region of AMA1 associated with strain-specific vaccine efficacy. Antibody measurements using whole antigens may be biased towards conserved, immunodominant epitopes. Peptide microarrays may help to identify immunogenic epitopes, define correlates of vaccine protection, and measure strain-specific vaccine-induced antibodies.

Comprehensive Profiling of the Rheumatoid Arthritis Antibody Repertoire.

OBJECTIVE: Autoantibodies against citrullinated proteins are found in 64-89% of rheumatoid arthritis (RA) patients, with 88-99% specificity. This study was undertaken to create an unbiased, comprehensive profile of serum antibodies against the human proteome, including the citrullinome and the homocitrullinome, in RA patients, using a high-density peptide array. METHODS: Our high-density peptide array, consisting of >4.6 million peptides, contained the entire annotated human proteome. The 20,246 proteins were represented as overlapping 16-mer peptides. In addition to native peptides, citrullinated and homocitrullinated peptides were included, as substitutions for arginine and lysine, and provided a comprehensive screen against all possible epitopes. Twenty-six serum samples (from 8 controls and 18 RA patients) were profiled on the high-density peptide array. Using RA-specific epitopes, we constructed an 8-epitope diagnostic biomarker on a Gyrolab xPlore instrument with a cohort of 92 serum samples (from 29 controls and 63 RA patients). The diagnostic biomarker was further validated with an independent cohort of 181 serum samples (from 54 controls and 127 RA patients). RESULTS: In the initial cohort the diagnostic performance of the 8-epitope biomarker yielded 96.6% specificity and 92.1% sensitivity. The overall diagnostic performance in the validation cohort was 94.4% specificity and 85% sensitivity. In both cohorts, the performance of the 8-epitope diagnostic biomarker compared favorably against the Abnova cyclic citrullinated peptide 2 (CCP2) assay. Using data from the peptide array, we identified novel RA-specific epitopes and formed the basis of a new RA diagnostic assay. CONCLUSION: Comprehensive antibody profiling using a high-density peptide array not only identified novel RA-specific epitopes but also allowed us to construct a novel diagnostic biomarker that is as specific as and more sensitive than the Abnova CCP2 assay.

Stepwise Evolution Improves Identification of Diverse Peptides Binding to a Protein Target.

Considerable efforts have been made to develop technologies for selection of peptidic molecules that act as substrates or binders to a protein of interest. Here we demonstrate the combination of rational peptide array library design, parallel screening and stepwise evolution, to discover novel peptide hotspots. These hotspots can be systematically evolved to create high-affinity, high-specificity binding peptides to a protein target in a reproducible and digitally controlled process. The method can be applied to synthesize both linear and cyclic peptides, as well as peptides composed of natural and non-natural amino acid analogs, thereby enabling screens in a much diverse chemical space. We apply this method to stepwise evolve peptide binders to streptavidin, a protein studied for over two decades and report novel peptides that mimic key interactions of biotin to streptavidin.

Effectiveness and limitations of local structural entropy optimization in the thermal stabilization of mesophilic and thermophilic adenylate kinases.

Local structural entropy (LSE) is a descriptor for the extent of conformational heterogeneity in short protein sequences that is computed from structural information derived from the Protein Data Bank. Reducing the LSE of a protein sequence by introducing amino acid mutations can result in fewer conformational states and thus a more stable structure, indicating that LSE optimization can be used as a protein stabilization method. Here, we describe a series of LSE optimization experiments designed to stabilize mesophilic and thermophilic adenylate kinases (AKs) and report crystal structures of LSE-optimized AK variants. In the mesophilic AK, thermal stabilization by LSE reduction was effective but limited. Structural analyses of the LSE-optimized mesophilic AK variants revealed a strong correlation between LSE and the apolar buried surface area. Additional mutations designed to introduce noncovalent interactions between distant regions of the polypeptide resulted in further stabilization. Unexpectedly, optimizing the LSE of the thermophilic AK resulted in a decrease in thermal stability. This destabilization was reduced when charged residues were excluded from the possible substitutions during LSE optimization. These observations suggest that stabilization by LSE reduction may result from the optimization of local hydrophobic contacts. The limitations of this process are likely due to ignorance of other interactions that bridge distant regions in a given amino acid sequence. Our results illustrate the effectiveness and limitations of LSE optimization as a protein stabilization strategy and highlight the importance and complementarity of local conformational stability and global interactions in protein thermal stability.

Optimal design of thermally stable proteins.

MOTIVATION: For many biotechnological purposes, it is desirable to redesign proteins to be more structurally and functionally stable at higher temperatures. For example, chemical reactions are intrinsically faster at higher temperatures, so using enzymes that are stable at higher temperatures would lead to more efficient industrial processes. We describe an innovative and computationally efficient method called Improved Configurational Entropy (ICE), which can be used to redesign a protein to be more thermally stable (i.e. stable at high temperatures). This can be accomplished by systematically modifying the amino acid sequence via local structural entropy (LSE) minimization. The minimization problem is modeled as a shortest path problem in an acyclic graph with nonnegative weights and is solved efficiently using Dijkstra's method.

Structure of human J-type co-chaperone HscB reveals a tetracysteine metal-binding domain.

Iron-sulfur proteins play indispensable roles in a broad range of biochemical processes. The biogenesis of iron-sulfur proteins is a complex process that has become a subject of extensive research. The final step of iron-sulfur protein assembly involves transfer of an iron-sulfur cluster from a cluster-donor to a cluster-acceptor protein. This process is facilitated by a specialized chaperone system, which consists of a molecular chaperone from the Hsc70 family and a co-chaperone of the J-domain family. The 3.0 A crystal structure of a human mitochondrial J-type co-chaperone HscB revealed an L-shaped protein that resembles Escherichia coli HscB. The important difference between the two homologs is the presence of an auxiliary metal-binding domain at the N terminus of human HscB that coordinates a metal via the tetracysteine consensus motif CWXCX(9-13)FCXXCXXXQ. The domain is found in HscB homologs from animals and plants as well as in magnetotactic bacteria. The metal-binding site of the domain is structurally similar to that of rubredoxin and several zinc finger proteins containing rubredoxin-like knuckles. The normal mode analysis of HscB revealed that this L-shaped protein preferentially undergoes a scissors-like motion that correlates well with the conformational changes of human HscB observed in the crystals.

Bioinformatic method for protein thermal stabilization by structural entropy optimization.

Engineering proteins for higher thermal stability is an important and difficult challenge. We describe a bioinformatic method incorporating sequence alignments to redesign proteins to be more stable through optimization of local structural entropy. Using this method, improved configurational entropy (ICE), we were able to design more stable variants of a mesophilic adenylate kinase with only the sequence information of one psychrophilic homologue. The redesigned proteins display considerable increases in their thermal stabilities while still retaining catalytic activity. ICE does not require a three-dimensional structure or a large number of homologous sequences, indicating a broad applicability of this method. Our results also highlight the importance of entropy in the stability of protein structures.

Structure and dynamics of gamma-SNAP: insight into flexibility of proteins from the SNAP family.

Soluble N-ethylmaleimide-sensitive factor attachment protein gamma (gamma-SNAP) is a member of an eukaryotic protein family involved in intracellular membrane trafficking. The X-ray structure of Brachydanio rerio gamma-SNAP was determined to 2.6 A and revealed an all-helical protein comprised of an extended twisted-sheet of helical hairpins with a helical-bundle domain on its carboxy-terminal end. Structural and conformational differences between multiple observed gamma-SNAP molecules and Sec17, a SNAP family protein from yeast, are analyzed. Conformational variation in gamma-SNAP molecules is matched with great precision by the two lowest frequency normal modes of the structure. Comparison of the lowest-frequency modes from gamma-SNAP and Sec17 indicated that the structures share preferred directions of flexibility, corresponding to bending and twisting of the twisted sheet motif. We discuss possible consequences related to the flexibility of the SNAP proteins for the mechanism of the 20S complex disassembly during the SNAP receptors recycling.

Protein structural variation in computational models and crystallographic data.

Normal mode analysis offers an efficient way of modeling the conformational flexibility of protein structures. We use anisotropic displacement parameters from crystallography to test the quality of prediction of both the magnitude and directionality of conformational flexibility. Normal modes from four simple elastic network model potentials and from the CHARMM force field are calculated for a data set of 83 diverse, ultrahigh-resolution crystal structures. While all five potentials provide good predictions of the magnitude of flexibility, all-atom potentials have a clear edge at prediction of directionality, and the CHARMM potential has the highest prediction quality. The low-frequency modes from different potentials are similar, but those computed from the CHARMM potential show the greatest difference from the elastic network models. The comprehensive evaluation demonstrates the costs and benefits of using normal mode potentials of varying complexity.

Structure and dynamics of UDP-glucose pyrophosphorylase from Arabidopsis thaliana with bound UDP-glucose and UTP.

The structure of the UDP-glucose pyrophosphorylase encoded by Arabidopsis thaliana gene At3g03250 has been solved to a nominal resolution of 1.86 Angstroms. In addition, the structure has been solved in the presence of the substrates/products UTP and UDP-glucose to nominal resolutions of 1.64 Angstroms and 1.85 Angstroms. The three structures revealed a catalytic domain similar to that of other nucleotidyl-glucose pyrophosphorylases with a carboxy-terminal beta-helix domain in a unique orientation. Conformational changes are observed between the native and substrate-bound complexes. The nucleotide-binding loop and the carboxy-terminal domain, including the suspected catalytically important Lys360, move in and out of the active site in a concerted fashion. TLS refinement was employed initially to model conformational heterogeneity in the UDP-glucose complex followed by the use of multiconformer refinement for the entire molecule. Normal mode analysis generated atomic displacement predictions in good agreement in magnitude and direction with the observed conformational changes and anisotropic displacement parameters generated by TLS refinement. The structures and the observed dynamic changes provide insight into the ordered mechanism of this enzyme and previously described oligomerization effects on catalytic activity.

Effect of low-complexity regions on protein structure determination.

It has been previously shown that protein sequences containing a quasi-repetitive assortment of amino acids are common in genomes and databases such as Swiss-Prot but are under-represented in the structure-based Protein Data Bank (PDB). Structural genomics groups have been using the absence of these "low-complexity" sequences for several years as a way to select proteins that have a good chance of successful structure determination. In this study, we examine the data deposited in the PDB as well as the available data from structural genomics groups in TargetDB and PepcDB to reveal interesting trends that could be taken into consideration when using low-complexity sequences as part of the target selection process.

Membrane association, mechanism of action, and structure of Arabidopsis embryonic factor 1 (FAC1).

Embryonic factor 1 (FAC1) is one of the earliest expressed plant genes and encodes an AMP deaminase (AMPD), which is also an identified herbicide target. This report identifies an N-terminal transmembrane domain in Arabidopsis FAC1, explores subcellular fractionation, and presents a 3.3-A globular catalytic domain x-ray crystal structure with a bound herbicide-based transition state inhibitor that provides the first glimpse of a complete AMPD active site. FAC1 contains an (alpha/beta)(8)-barrel characterized by loops in place of strands 5 and 6 that places it in a small subset of the amidohydrolase superfamily with imperfect folds. Unlike tetrameric animal orthologs, FAC1 is a dimer and each subunit contains an exposed Walker A motif that may be involved in the dramatic combined K(m) (25-80-fold lower) and V(max) (5-6-fold higher) activation by ATP. Normal mode analysis predicts a hinge motion that flattens basic surfaces on each monomer that flank the dimer interface, which suggests a reversible association between the FAC1 globular catalytic domain and intracellular membranes, with N-terminal transmembrane and disordered linker regions serving as the anchor and attachment to the globular catalytic domain, respectively.