Protein structural plasticity exemplified by insertion and deletion mutants in T4 lysozyme

High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST

The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.

A review of visualisations of protein fold networks and their relationship with sequence and function

Proteins form arguably the most significant link between genotype and phenotype. Understanding the relationship between protein sequence and structure, and applying this knowledge to predict function, is difficult. One way to investigate these relationships is by considering the space of protein folds and how one might move from fold to fold through similarity, or potential evolutionary relationships. The many individual characterisations of fold space presented in the literature can tell us a lot about how well the current Protein Data Bank represents protein fold space, how convergence and divergence may affect protein evolution, how proteins affect the whole of which they are part, and how proteins themselves function. A synthesis of these different approaches and viewpoints seems the most likely way to further our knowledge of protein structure evolution and thus, facilitate improved protein structure design and prediction.

Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

An atlas of amyloid aggregation: the impact of substitutions, insertions, deletions and truncations on amyloid beta fibril nucleation

Multiplexed assays of variant effects (MAVEs) guide clinical variant interpretation and reveal disease mechanisms. To date, MAVEs have focussed on a single mutation type-amino acid (AA) substitutions-despite the diversity of coding variants that cause disease. Here we use Deep Indel Mutagenesis (DIM) to generate a comprehensive atlas of diverse variant effects for a disease protein, the amyloid beta (Aβ) peptide that aggregates in Alzheimer's disease (AD) and is mutated in familial AD (fAD). The atlas identifies known fAD mutations and reveals that many variants beyond substitutions accelerate Aβ aggregation and are likely to be pathogenic. Truncations, substitutions, insertions, single- and internal multi-AA deletions differ in their propensity to enhance or impair aggregation, but likely pathogenic variants from all classes are highly enriched in the polar N-terminal region of Aβ. This comparative atlas highlights the importance of including diverse mutation types in MAVEs and provides important mechanistic insights into amyloid nucleation.

Deep RNA Sequencing Revealed Fusion Junctional Heterogeneity May Predict Crizotinib Treatment Efficacy in ALK-Rearranged NSCLC

INTRODUCTION

Gene fusion variants in ALK-rearranged NSCLC may predict patient outcomes, but previous results have been inconclusive. Fusion isoforms coexisting in the same tumor may affect the efficacy of targeted therapy, but they have not been investigated.

METHODS

Patients with ALK-rearranged NSCLC who received crizotinib treatments were recruited. Precrizotinib tumor tissues were analyzed by the anchored multiplex polymerase chain reaction for targeted RNA sequencing. Kaplan-Meier and Cox regression were used to compare overall and progression-free survivals.

RESULTS

Of the 51 studied subjects, EML4-ALK variant types v1, v2, v3, and others were detected in 23 (45.1%), five (9.8%), 19 (37.3%), and four patients (7.8%), respectively. Multiple EML4-ALK RNA isoforms were detected in 24 tumors (47.1%), and single isoform in 27 (52.9%). Most of the v3 tumors (16 of 19) harbored both v3a and v3b RNA isoforms. Multiple isoforms were also detected in eight non-v3 tumors (33.3% of all 24 multiple isoforms; five v1, two v5', and one v2). Compared with patients with single isoform, those with multiple isoforms had worse progression-free (hazard ratio and 95% confidence interval: 2.45 [1.06-5.69]) and overall (hazard ratio [95% confidence interval]: 3.74 [1.26-11.13]) survivals after adjusting for potential confounders including variant type. Using the patient-derived H2228 cells known to express v3a and v3b, our single-cell polymerase chain reaction detected either v3a or v3b in most single cells. Treatment of H2228 cells by three ALK inhibitors revealed increased ratios of v3a-to-v3b expression over time.

CONCLUSIONS

Intratumoral EML4-ALK isoforms may predict the efficacy of targeted therapy in ALK-rearranged NSCLC. Temporal changes of intratumoral fusion isoforms may result from differential selection pressures that a drug might have on one isoform over another. Larger studies on fusion heterogeneity using RNA sequencing are warranted.

Structure elements can be predicted using the contact volume among protein residues

Previously, the structure elements of dihydrofolate reductase (DHFR) were determined using comprehen­sive Ala-insertion mutation analysis, which is assumed to be a kind of protein “building blocks.” It is hypo­thesized that our comprehension of the structure elements could lead to understanding how an amino acid sequence dictates its tertiary structure. However, the comprehensive Ala-insertion mutation analysis is a time- and cost-consuming process and only a set of the DHFR structure elements have been reported so far. Therefore, developing a computational method to predict structure elements is an urgent necessity. We focused on intramolecular residue–residue contacts to predict the structure elements. We introduced a simple and effective parameter: the overlapped contact volume (CV) among the residues and calculated the CV along the DHFR sequence using the crystal structure. Our results indicate that the CV profile can recapitulate its precipitate ratio profile, which was used to define the structure elements in the Ala-insertion mutation analysis. The CV profile allowed us to predict structure elements like the experimentally determined structure elements. The strong correlation between the CV and precipitate ratio profiles indicates the importance of the intramolecular residue–residue contact in maintaining the tertiary structure. Additionally, the CVs between the structure elements are considerably more than those between a structure element and a linker or two linkers, indicating that the structure elements play a funda­mental role in increasing the intramolecular adhesion. Thus, we propose that the structure elements can be considered a type of “building blocks” that maintain and dictate the tertiary structures of proteins.

Rosetta-Enabled Structural Prediction of Permissive Loop Insertion Sites in Proteins

While loop motifs frequently play a major role in protein function, our understanding of how to rationally engineer proteins with novel loop domains remains limited. In the absence of rational approaches, the incorporation of loop domains often destabilizes proteins, thereby requiring massive screening and selection to identify sites that can accommodate loop insertion. We developed a computational strategy for rapidly scanning the entire structure of a scaffold protein to determine the impact of loop insertion at all possible amino acid positions. This approach is based on the Rosetta kinematic loop modeling protocol and was demonstrated by identifying sites in lipase that were permissive to insertion of the LAP peptide. Interestingly, the identification of permissive sites was dependent on the contribution of the residues in the near-loop environment on the Rosetta score and did not correlate with conventional structural features (e.g., B-factors). As evidence of this, several insertion sites (e.g., following residues 17, 47-49, and 108), which were predicted and confirmed to be permissive, interrupted helices, while others (e.g., following residues 43, 67, 116, 119, and 121), which are situated in loop regions, were nonpermissive. This approach was further shown to be predictive for β-glucosidase and human phosphatase and tensin homologue (PTEN), and to facilitate the engineering of insertion sites through in silico mutagenesis. By enabling the design of loop-containing protein libraries with high probabilities of soluble expression, this approach has broad implications in many areas of protein engineering, including antibody design, improving enzyme activity, and protein modification.

Cell-cell and virus-cell fusion assay-based analyses of alanine insertion mutants in the distal α9 portion of the JRFL gp41 subunit from HIV-1

Membrane fusion is the first essential step in HIV-1 replication. The gp41 subunit of HIV-1 envelope protein (Env), a class I fusion protein, achieves membrane fusion by forming a structure called a six-helix bundle composed of N- and C-terminal heptad repeats. We have recently shown that the distal portion of the α9 helix in the C-terminal heptad repeat of X4-tropic HXB2 Env plays a critical role in the late-stage membrane fusion and viral infection. Here, we used R5-tropic JRFL Env and constructed six alanine insertion mutants, 641+A to 646+A, in the further distal portion of α9 where several glutamine residues are conserved (the number corresponds to the position of the inserted alanine in JRFL Env). 644+A showed the most severe defect in syncytia formation. Decreased fusion pore formation activity, revealed by a dual split protein assay, was observed in all mutants except 641+A. Sequence analysis and substitution of inserted 644A with Gln revealed that the glutamine residue at position 644 that forms complex hydrogen-bond networks with other polar residues on the surface of the six-helix bundle is critical for cell-cell fusion. We also developed a split NanoLuc® (Nluc) reporter-based assay specific to the virus-cell membrane fusion step to analyze several of the mutants. Interestingly syncytia-competent mutants failed to display Nluc activities. In addition to defective fusion activity, a reduction of Env incorporation into virions may further contribute to differences in cell-cell and virus-cell fusions.

Multistep mutational transformation of a protein fold through structural intermediates

New protein folds may evolve from existing folds through metamorphic evolution involving a dramatic switch in structure. To mimic pathways by which amino acid sequence changes could induce a change in fold, we designed two folded hybrids of Xfaso 1 and Pfl 6, a pair of homologous Cro protein sequences with ~40% identity but different folds (all-α vs. α + β, respectively). Each hybrid, XPH1 or XPH2, is 85% identical in sequence to its parent, Xfaso 1 or Pfl 6, respectively; 55% identical to its noncognate parent; and ~70% identical to the other hybrid. XPH1 and XPH2 also feature a designed hybrid chameleon sequence corresponding to the C-terminal region, which switched from α-helical to β-sheet structure during Cro evolution. We report solution nuclear magnetic resonance (NMR) structures of XPH1 and XPH2 at 0.3 Å and 0.5 Å backbone root mean square deviation (RMSD), respectively. XPH1 retains a global fold generally similar to Xfaso 1, and XPH2 retains a fold similar to Pfl 6, as measured by TM-align scores (~0.7), DALI Z-scores (7-9), and backbone RMSD (2-3 Å RMSD for the most ordered regions). However, these scores also indicate significant deviations in structure. Most notably, XPH1 and XPH2 have different, and intermediate, secondary structure content relative to Xfaso 1 and Pfl 6. The multistep progression in sequence, from Xfaso 1 to XPH1 to XPH2 to Pfl 6, thus involves both abrupt and gradual changes in folding pattern. The plasticity of some protein folds may allow for "polymetamorphic" evolution through intermediate structures.

Competitive species interactions constrain abiotic adaptation in a bacterial soil community

Studies of abiotic adaptation often consider single species in isolation, yet natural communities contain many coexisting species which could limit or promote abiotic adaptation. Here we show, using soil bacterial communities, that evolving in the presence of a competitor constrained abiotic adaptation. Specifically, Pseudomonas fluorescens evolved alone was fitter than P. fluorescens evolved alongside Pseudomonas putida, when P. putida was absent. Genome analyses indicated this was due to mutation of the acetate scavenger actP, which occurred exclusively, and almost universally, in single‐species‐evolved clones. actP disruption was associated with increased growth in soil compared with wild‐type actP, but this benefit was abolished when P. putida was present, suggesting a role for carbon scavenging transporters in species interactions, possibly through nutrient competition. Our results show that competitive species interactions can limit the evolutionary response to abiotic selection, because the fitness benefits of abiotic adaptive mutations were negated in more complex communities.

Insertion of antihypertensive peptides in acidic subunit from amaranth 11S induces contrasting effects in stability

The insertion of peptides is a biotechnology tool widely used to improve the nutraceutical properties of proteins. Because the effect of these insertions in protein stability and function is difficult to predict, it should be determined experimentally. In this study, we created two variants of amarantin acidic subunit and analyzed them along with other four proteins reported previously. We measured their response against two destabilizing agents: temperature and urea. The six proteins presented the insertion of antihypertensive peptides (VYVYVYVY or RIPP) in the variable regions of the protein. We observed that their effect strongly depended on the site of the insertion. The insertion in the variable region I stabilized the protein both thermally and chemically, but it affected the inhibitory activity of the angiotensin-converting enzyme in vitro. In contrast, insertions in other three regions were severely destabilizing, producing molten globules. Our findings reveal that the insertion of bioactive peptides in variable regions of a protein can increase or decrease the protein's thermal and chemical stability and that these conformational changes may also alter its final activity.

Six-helix bundle completion in the distal C-terminal heptad repeat region of gp41 is required for efficient human immunodeficiency virus type 1 infection

Background

The native pre-fusion structure of gp120/gp41 complex of human immunodeficiency virus type 1 was recently revealed. In the model, the helices of gp41 (α6, α7, α8, and α9) form a four-helix collar underneath trimeric gp120. Gp41 is a class I fusion protein and mediates membrane fusion by forming a post-fusion structure called the six-helix bundle (6HB). The comparison of the pre- and post-fusion structures revealed the large conformational changes in gp41 during the antiparallel packing of the N- and C-terminal heptad repeats (NHRs and CHRs) in membrane fusion. Several mutagenesis studies of gp41 performed in the past were interpreted based on 6HB, the only available structure at that time. To obtain an insight about the current pre-fusion structural model and conformational changes during membrane fusion, alanine insertion mutagenesis of the NHR, CHR and connecting loop regions of HXB2 gp41 was performed. The effects of mutations on biosynthesis and membrane fusion were analyzed by immunoblotting and fusion assays, respectively. The extent of membrane fusion was evaluated by split luciferase-based pore formation and syncytia formation assays, respectively.

Results

Consistent with the current structural model, drastic negative effects of mutations on biosynthesis and membrane fusion were observed for NHR, loop, and proximal regions of CHR (up to amino acid position 643). The insertions in α9 after it leaves the four-helix collar were tolerable for biosynthesis. These CHR mutants showed varying effects on membrane fusion. Insertion at position 644 or 645 resulted in poor pore and syncytia formation. Efficient pore and syncytia formation almost similar to that of the wild type was observed for insertion at position 647, 648 or 649. However, recovery of virus infectivity was only observed for the insertions beyond position 648.

Conclusions

The mutagenesis data for HXB2 gp41 is in agreement with the recent pre-fusion structure model. The virus infection data suggested that fusion pores sufficiently large enough for the release of the virus genome complex are formed after the completion of 6HB beyond position 648.

Electronic supplementary material

The online version of this article (10.1186/s12977-018-0410-9) contains supplementary material, which is available to authorized users.

Computational prediction of the tolerance to amino-acid deletion in green-fluorescent protein

Proteins evolve through two primary mechanisms: substitution, where mutations alter a protein's amino-acid sequence, and insertions and deletions (indels), where amino acids are either added to or removed from the sequence. Protein structure has been shown to influence the rate at which substitutions accumulate across sites in proteins, but whether structure similarly constrains the occurrence of indels has not been rigorously studied. Here, we investigate the extent to which structural properties known to covary with protein evolutionary rates might also predict protein tolerance to indels. Specifically, we analyze a publicly available dataset of single-amino-acid deletion mutations in enhanced green fluorescent protein (eGFP) to assess how well the functional effect of deletions can be predicted from protein structure. We find that weighted contact number (WCN), which measures how densely packed a residue is within the protein's three-dimensional structure, provides the best single predictor for whether eGFP will tolerate a given deletion. We additionally find that using protein design to explicitly model deletions results in improved predictions of functional status when combined with other structural predictors. Our work suggests that structure plays fundamental role in constraining deletions at sites in proteins, and further that similar biophysical constraints influence both substitutions and deletions. This study therefore provides a solid foundation for future work to examine how protein structure influences tolerance of more complex indel events, such as insertions or large deletions.

Atomistic picture of conformational exchange in a T4 lysozyme cavity mutant: an experiment-guided molecular dynamics study

Despite the importance of dynamics to protein function there is little information about the states that are formed as the protein explores its conformational landscape or about the mechanism by which transitions between the different states occur. Here we used a combined NMR spin relaxation and unbiased molecular dynamics (MD) approach to investigate the exchange process by which a cavity in an L99A mutant of T4 lysozyme (T4L 99A) interconverts between an empty and occupied form that involves repositioning of an aromatic residue, Phe114. Although structures of the end-states of the exchange process are available, insight into the mechanism by which the transition takes place cannot be obtained from experiment and the timescales involved are too slow to address using brute force MD. Using spin relaxation NMR methods, we have identified a triple-mutant of T4L that undergoes the same exchange process as T4L L99A but where the minor state lifetime has decreased significantly so that the spontaneous conformational transition to the major state can be studied using all-atom MD simulations. The simulation trajectories were analyzed using Markov state models and the energy landscape so obtained is in good agreement with expectations based on NMR studies. Notably there is no large-scale perturbation of the structure during the transition, multiple intermediates are formed between the two similar exchanging conformations and the free energy barrier between these two well-folded, compact forms is small (6kBT), only slightly larger than for processes considered to be barrierless.

Comparative analysis of genome-wide Mlo gene family in Cajanus cajan and Phaseolus vulgaris

The Mlo gene was discovered in barley because the mutant 'mlo' allele conferred broad-spectrum, non-race-specific resistance to powdery mildew caused by Blumeria graminis f. sp. hordei. The Mlo genes also play important roles in growth and development of plants, and in responses to biotic and abiotic stresses. The Mlo gene family has been characterized in several crop species, but only a single legume species, soybean (Glycine max L.), has been investigated so far. The present report describes in silico identification of 18 CcMlo and 20 PvMlo genes in the important legume crops Cajanus cajan (L.) Millsp. and Phaseolus vulgaris L., respectively. In silico analysis of gene organization, protein properties and conserved domains revealed that the C. cajan and P. vulgaris Mlo gene paralogs are more divergent from each other than from their orthologous pairs. The comparative phylogenetic analysis classified CcMlo and PvMlo genes into three major clades. A comparative analysis of CcMlo and PvMlo proteins with the G. max Mlo proteins indicated close association of one CcMlo, one PvMlo with two GmMlo genes, indicating that there was no further expansion of the Mlo gene family after the separation of these species. Thus, most of the diploid species of eudicots might be expected to contain 15-20 Mlo genes. The genes CcMlo12 and 14, and PvMlo11 and 12 are predicted to participate in powdery mildew resistance. If this prediction were verified, these genes could be targeted by TILLING or CRISPR to isolate powdery mildew resistant mutants.

Structural and dynamic changes associated with beneficial engineered single-amino-acid deletion mutations in enhanced green fluorescent protein

Single-amino-acid deletions are a common part of the natural evolutionary landscape but are rarely sampled during protein engineering owing to limited and prejudiced molecular understanding of mutations that shorten the protein backbone. Single-amino-acid deletion variants of enhanced green fluorescent protein (EGFP) have been identified by directed evolution with the beneficial effect of imparting increased cellular fluorescence. Biophysical characterization revealed that increased functional protein production and not changes to the fluorescence parameters was the mechanism that was likely to be responsible. The structure EGFP(D190Δ) containing a deletion within a loop revealed propagated changes only after the deleted residue. The structure of EGFP(A227Δ) revealed that a `flipping' mechanism was used to adjust for residue deletion at the end of a β-strand, with amino acids C-terminal to the deletion site repositioning to take the place of the deleted amino acid. In both variants new networks of short-range and long-range interactions are generated while maintaining the integrity of the hydrophobic core. Both deletion variants also displayed significant local and long-range changes in dynamics, as evident by changes in B factors compared with EGFP. Rather than being detrimental, deletion mutations can introduce beneficial structural effects through altering core protein properties, folding and dynamics, as well as function.

Random single amino acid deletion sampling unveils structural tolerance and the benefits of helical registry shift on GFP folding and structure

Altering a protein’s backbone through amino acid deletion is a common evolutionary mutational mechanism, but is generally ignored during protein engineering primarily because its effect on the folding-structure-function relationship is difficult to predict. Using directed evolution, enhanced green fluorescent protein (EGFP) was observed to tolerate residue deletion across the breadth of the protein, particularly within short and long loops, helical elements, and at the termini of strands. A variant with G4 removed from a helix (EGFP^G4Δ) conferred significantly higher cellular fluorescence. Folding analysis revealed that EGFP^G4Δ retained more structure upon unfolding and refolded with almost 100% efficiency but at the expense of thermodynamic stability. The EGFP^G4Δ structure revealed that G4 deletion caused a beneficial helical registry shift resulting in a new polar interaction network, which potentially stabilizes a cis proline peptide bond and links secondary structure elements. Thus, deletion mutations and registry shifts can enhance proteins through structural rearrangements not possible by substitution mutations alone.

•

Using directed evolution, the impact of amino acid deletion on EGFP is explored

•

Loops, helices, and strand termini are especially tolerant to amino acid deletion

•

A deletion mutant that enhances cellular production and fluorescence is identified

•

Structure reveals that a helical registry shift creates a new polar network

A role for indels in the evolution of Cro protein folds

Insertions and deletions in protein sequences, or indels, can disrupt structure and may result in changes in protein folds during evolution or in association with alternative splicing. Pfl 6 and Xfaso 1 are two proteins in the Cro family that share a common ancestor but have different folds. Sequence alignments of the two proteins show two gaps, one at the N terminus, where the sequence of Xfaso 1 is two residues shorter, and one near the center of the sequence, where the sequence of Pfl 6 is five residues shorter. To test the potential importance of indels in Cro protein evolution, we generated hybrid variants of Pfl 6 and Xfaso 1 with indels in one or both regions, chosen according to several plausible sequence alignments. All but one deletion variant completely unfolded both proteins, showing that a longer N-terminal sequence was critical for Pfl 6 folding and a longer central region sequence was critical for Xfaso 1 folding. By contrast, Xfaso 1 tolerated a longer N-terminal sequence with little destabilization, and Pfl 6 tolerated central region insertions, albeit with substantial effects on thermal stability and some perturbation of the surrounding structure. None of the mutations appeared to convert one stable fold into the other. On the basis of this two-protein comparison, short insertion and deletion mutations probably played a role in evolutionary fold change in the Cro family, but were also not the only factors.

Crystallization and preliminary X-ray crystallographic analysis of the oxysterol-binding protein Osh3 from Saccharomyces cerevisiae

Oxysterol-binding protein (OSBP) related proteins (ORPs) are conserved from yeast to humans and are implicated in regulation of sterol homeostasis and in signal transduction pathways. Osh3 of Saccharomyces cerevisiae is a pleckstrin-homology (PH) domain-containing ORP member that regulates phosphoinositide metabolism at endoplasmic reticulum-plasma membrane contact sites. The N-terminal PH domain of Osh3 was purified and crystallized as a lysozyme fusion and the resulting crystal diffracted to 2.3 Å resolution. The crystal belonged to the monoclinic space group C2, with unit-cell parameters a=98.03, b=91.31, c=84.13 Å, β=81.41°. With two molecules in the asymmetric unit, the Matthews coefficient was 3.13 Å3 Da(-1). Initial attempts to solve the structure by molecular-replacement techniques using T4 lysozyme as a search model were successful. The C-terminal OSBP-related domain (OBD) of Osh3 was crystallized by the vapour-diffusion method and the resulting crystal diffracted to 1.5 Å resolution. The crystal was orthorhombic, belonging to space group P2(1)2(1)2(1), with unit-cell parameters a=41.57, b=87.52, c=100.58 Å. With one molecule in the asymmetric unit, the Matthews coefficient was 2.01 Å3 Da(-1). Initial attempts to solve the structure by the single-wavelength anomalous dispersion technique using bromine were successful.

Evolutionary pressures on simple sequence repeats in prokaryotic coding regions

Simple sequence repeats (SSRs) are indel mutational hotspots in genomes. In prokaryotes, SSR loci can cause phase variation, a microbial survival strategy that relies on stochastic, reversible on–off switching of gene activity. By analyzing multiple strains of 42 fully sequenced prokaryotic species, we measure the relative variability and density distribution of SSRs in coding regions. We demonstrate that repeat type strongly influences indel mutation rates, and that the most mutable types are most strongly avoided across genomes. We thoroughly characterize SSR density and variability as a function of N→C position along protein sequences. Using codon-shuffling algorithms that preserve amino acid sequence, we assess evolutionary pressures on SSRs. We find that coding sequences suppress repeats in the middle of proteins, and enrich repeats near termini, yielding U-shaped SSR density curves. We show that for many species this characteristic shape can be attributed to purely biophysical constraints of protein structure. In multiple cases, however, particularly in certain pathogenic bacteria, we observe over enrichment of SSRs near protein N-termini significantly beyond expectation based on structural constraints. This increases the probability that frameshifts result in non-functional proteins, revealing that these species may evolutionarily tune SSR positions in coding regions to facilitate phase variation.

Lessons from the lysozyme of phage T4

An overview is presented of some of the major insights that have come from studies of the structure, stability, and folding of T4 phage lysozyme. A major purpose of this review is to provide the reader with a complete tabulation of all of the variants that have been characterized, including melting temperatures, crystallographic data, Protein Data Bank access codes, and references to the original literature. The greatest increase in melting temperature (T(m)) for any point mutant is 5.1 degrees C for the mutant Ser 117 --> Val. This is achieved in part not only by hydrophobic stabilization but also by eliminating an unusually short hydrogen bond of 2.48 A that apparently has an unfavorable van der Waals contact. Increases in T(m) of more than 3-4 degrees C for point mutants are rare, whereas several different types of destabilizing substitutions decrease T(m) by 20 degrees C or thereabouts. The energetic cost of cavity creation and its relation to the hydrophobic effect, derived from early studies of "large-to-small" mutants in the core of T4 lysozyme, has recently been strongly supported by related studies of the intrinsic membrane protein bacteriorhodopsin. The L99A cavity in the C-terminal domain of the protein, which readily binds benzene and many other ligands, has been the subject of extensive study. Crystallographic evidence, together with recent NMR analysis, suggest that these ligands are admitted by a conformational change involving Helix F and its neighbors. A total of 43 nonisomorphous crystal forms of different monomeric lysozyme mutants were obtained plus three more for synthetically-engineered dimers. Among the 43 space groups, P2(1)2(1)2(1) and P2(1) were observed most frequently, consistent with the prediction of Wukovitz and Yeates.

Systematic analysis of short internal indels and their impact on protein folding

Background

Protein sequence insertions/deletions (indels) can be introduced during evolution or through alternative splicing (AS). Alternative splicing is an important biological phenomenon and is considered as the major means of expanding structural and functional diversity in eukaryotes. Knowledge of the structural changes due to indels is critical to our understanding of the evolution of protein structure and function. In addition, it can help us probe the evolution of alternative splicing and the diversity of functional isoforms. However, little is known about the effects of indels, in particular the ones involving core secondary structures, on the folding of protein structures. The long term goal of our study is to accurately predict the protein AS isoform structures. As a first step towards this goal, we performed a systematic analysis on the structural changes caused by short internal indels through mining highly homologous proteins in Protein Data Bank (PDB).

Results

We compiled a non-redundant dataset of short internal indels (2-40 amino acids) from highly homologous protein pairs and analyzed the sequence and structural features of the indels. We found that about one third of indel residues are in disordered state and majority of the residues are exposed to solvent, suggesting that these indels are generally located on the surface of proteins. Though naturally occurring indels are fewer than engineered ones in the dataset, there are no statistically significant differences in terms of amino acid frequencies and secondary structure types between the "Natural" indels and "All" indels in the dataset. Structural comparisons show that all the protein pairs with short internal indels in the dataset preserve the structural folds and about 85% of protein pairs have global RMSDs (root mean square deviations) of 2Å or less, suggesting that protein structures tend to be conserved and can tolerate short insertions and deletions. A few pairs with high RMSDs are results of relative domain positions of the proteins, probably due to the intrinsically dynamic nature of the proteins.

Conclusions

The analysis demonstrated that protein structures have the "plasticity" to tolerate short indels. This study can provide valuable guides in modeling protein AS isoform structures and homologous proteins with indels through placing the indels at the right locations since the accuracy of sequence alignments dictate model qualities in homology modeling.

A survey of proteins encoded by non-synonymous single nucleotide polymorphisms reveals a significant fraction with altered stability and activity

On average, each human gene has approximately four SNPs (single nucleotide polymorphisms) in the coding region, half of which are nsSNPs (non-synonymous SNPs) or missense SNPs. Current attention is focused on those that are known to perturb function and are strongly linked to disease. However, the vast majority of SNPs have not been investigated for the possibility of causing disease. We set out to assess the fraction of nsSNPs that encode proteins that have altered stability and activity, for this class of variants would be candidates to perturb cellular function. We tested the thermostability and, where possible, the catalytic activity for the most common variant (wild-type) and minor variants (total of 46 SNPs) for 16 human enzymes for which the three-dimensional structures were known. There were significant differences in the stability of almost half of the variants (48%) compared with their wild-type counterparts. The catalytic efficiency of approx. 14 variants was significantly altered, including several variants of human PKM2 (pyruvate kinase muscle 2). Two PKM2 variants, S437Y and E28K, also exhibited changes in their allosteric regulation compared with the wild-type enzyme. The high proportion of nsSNPs that affect protein stability and function, albeit subtly, underscores the need for experimental analysis of the diverse human proteome.

Amino acid insertion reveals a necessary three-helical intermediate in the folding pathway of the colicin E7 immunity protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ϕ-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.

Structural and thermodynamic characterization of T4 lysozyme mutants and the contribution of internal cavities to pressure denaturation

Using small-angle X-ray scattering (SAXS) and tryptophan fluorescence spectroscopy, we have identified multiple compact denatured states of a series of T4 lysozyme mutants that are stabilized by high pressures. Recent studies imply that the mechanism of pressure denaturation is the penetration of water into the protein rather than the transfer of hydrophobic residues into water. To investigate water penetration and the volume change associated with pressure denaturation, we studied the solution behavior of four T4 lysozyme mutants having different cavity volumes at low and neutral pH up to a pressure of 400 MPa (0.1 MPa = 0.9869 atm). At low pH, L99A T4 lysozyme expanded from a compact folded state to a partially unfolded state with a corresponding change in radius of gyration from 17 to 32 A. The volume change upon denaturation correlated well with the total cavity volume, indicating that all of the molecule's major cavities are hydrated with pressure. As a direct comparison to high-pressure crystal structures of L99A T4 lysozyme solved at neutral pH [Collins, M. D., Hummer, G., Quillin, M. L., Matthews, B. W., and Gruner, S. M. (2005) Proc. Natl. Acad. Sci. U.S.A. 102, 16668-16671], pressure denaturation of L99A and the structurally similar L99G/E108V mutant was studied at neutral pH. The pressure-denatured state at neutral pH is even more compact than at low pH, and the small volume changes associated with denaturation suggest that the preferential filling of large cavities is responsible for the compactness of the pressure-denatured state. These results confirm that pressure denaturation is characteristically distinct from thermal or chemical denaturation.

Development of a family of redox-sensitive green fluorescent protein indicators for use in relatively oxidizing subcellular environments

Green fluorescent protein (GFP) indicators were previously developed that rapidly and quantitatively respond to changes in the thiol/disulfide equilibrium within subcellular compartments. In these indicators, surface-exposed cysteines residues were introduced so as to form a labile redox-active disulfide that in turn controls the emission properties of the internal chromophore. The biosensors have been shown to be effective reporters of the thiol/disulfide status within reducing compartments such as the mitochondria and cytosol for several cell types. However, due to the high thermodynamic stability of the introduced disulfide bond, the indicators are not useful for quantitative analysis within more oxidizing compartments such as the endoplasmic reticulum. Here we report the development of a new family of GFP-based redox indicators (roGFP1-iX) in which the thermodynamic stability of the disulfide is substantially lowered by insertion of a single amino acid into the main chain, adjacent to cysteine 147. The insertions result in indicators with midpoint potentials of -229 to -246 mV and are thus better suited for study of relatively oxidizing subcellular compartments. Atomic resolution crystallographic analyses suggest that two important factors act to destabilize the disulfide linkage in roGFP1-iX. In the oxidized state, an unusual non-proline cis-peptide bond adjacent to one of the cysteines introduces geometric strain into the system, while in the reduced state, a dramatic loop opening lowers the effective concentration of the reacting species.

Predicting the magnitude of the reflex response to insertions in ubiquitin

The ability to predict the structural response of a protein to an insertion would be a significant advance for the fields of homology modeling and protein design. However, the effects of insertions on protein conformation are not well understood. Previous work has demonstrated that for two loops in ubiquitin, the primary determinant of the structural adaptation to insertions is the insertion site rather than the sequence of the insertion; this phenomenon was termed the reflex response of loops to insertions. We report herein the analysis of ubiquitin mutants with insertions in two other loops. This study demonstrates that the insertion site is the primary determinant of the response to insertions for these two new loops as well, which further supports the reflex response hypothesis. We also attempted to predict the relative magnitudes of the responses at each site but were unsuccessful. Using the additional data collected in this work, we have refined our predictive hypothesis.

Investigating protein structural plasticity by surveying the consequence of an amino acid deletion from TEM-1 beta-lactamase

While the deletion of an amino acid is a common mutation observed in nature, it is generally thought to be disruptive to protein structure. Using a directed evolution approach, we find that the enzyme TEM-1 beta-lactamase was broadly tolerant to the deletion mutations sampled. Circa 73% of the variants analysed retained activity towards ampicillin, with deletion mutations observed in helices and strands as well as regions important for structure and function. Several deletion variants had enhanced activity towards ceftazidime compared to the wild-type TEM-1 demonstrating that removal of an amino acid can have a beneficial outcome.

A posteriori design of crystal contacts to improve the X-ray diffraction properties of a small RNA enzyme

The hairpin ribozyme is a small catalytic RNA comprising two helix-loop-helix domains linked by a four-way helical junction (4WJ). In its most basic form, each domain can be formed independently and reconstituted without a 4WJ to yield an active enzyme. The production of such minimal junctionless hairpin ribozymes is achievable by chemical synthesis, which has allowed structures to be determined for numerous nucleotide variants. However, abasic and other destabilizing core modifications hinder crystallization. This investigation describes the use of a dangling 5'-U to form an intermolecular U.U mismatch, as well as the use of synthetic linkers to tether the loop A and B domains, including (i) a three-carbon propyl linker (C3L) and (ii) a nine-atom triethylene glycol linker (S9L). Both linker constructs demonstrated similar enzymatic activity, but S9L constructs yielded crystals that diffracted to 2.65 A resolution or better. In contrast, C3L variants diffracted to 3.35 A and exhibited a 15 A expansion of the c axis. Crystal packing of the C3L construct showed a paucity of 6(1) contacts, which comprise numerous backbone to 2'-OH hydrogen bonds in junctionless and S9L complexes. Significantly, the crystal packing in minimal structures mimics stabilizing features observed in the 4WJ hairpin ribozyme structure. The results demonstrate how knowledge-based design can be used to improve diffraction and overcome otherwise destabilizing defects.

Plasticity of amyloid fibrils

In experiments designed to characterize the basis of amyloid fibril stability through mutational analysis of the Abeta (1-40) molecule, fibrils exhibit consistent, significant structural malleability. In these results, and in other properties, amyloid fibrils appear to more resemble plastic materials generated from synthetic polymers than globular proteins. Thus, like synthetic polymers and plastics, amyloid fibrils exhibit both polymorphism, the ability of one polypeptide to form aggregates of different morphologies, and isomorphism, the ability of different polypeptides to grow into a fibrillar amyloid morphology. This view links amyloid with the prehistorical and 20th century use of proteins as starting materials to make films, fibers, and plastics, and with the classic protein fiber stretching experiments of the Astbury group. Viewing amyloids from the point of view of the polymer chemist may shed new light on a number of issues, such as the role of protofibrils in the mechanism of amyloid formation, the biological potency of fibrils, and the prospects for discovering inhibitors of amyloid fibril formation.

Active TEM-1 beta-lactamase mutants with random peptides inserted in three contiguous surface loops

Engineering of alternative binding sites on the surface of an enzyme while preserving the enzymatic activity would offer new opportunities for controlling the activity by binding of non-natural ligands. Loops and turns are the natural substructures in which binding sites might be engineered with this purpose. We have genetically inserted random peptide sequences into three relatively rigid and contiguous loops of the TEM-1 beta-lactamase and assessed the tolerance to insertion by the percentage of active mutants. Our results indicate that tolerance to insertion could not be correlated to tolerance to mutagenesis. A turn between two beta-strands bordering the active site was observed to be tolerant to random mutagenesis but not to insertions. Two rigid loops comprising rather well-conserved amino acid residues tolerated insertions, although with some constraints. Insertions between the N-terminal helix and the first beta-strand generated active libraries if cysteine residues were included at both ends of the insert, suggesting the requirement for a stabilizing disulfide bridge. Random sequences were relatively well accommodated within the loop connecting the final beta-strand to the C-terminal helix, particularly if the wild-type residue was retained at one of the loops' end. This suggests two strategies for increasing the percentage of active mutants in insertion libraries. The amino acid distribution in the engineered loops was analyzed and found to be less biased against hydrophobic residues than in natural medium-sized loops. The combination of these activity-selected libraries generated a huge library containing active hybrid enzymes with all three loops modified.

Sequential reorganization of beta-sheet topology by insertion of a single strand

Insertions, duplications, and deletions of sequence segments are thought to be major evolutionary mechanisms that increase the structural and functional diversity of proteins. Alternative splicing, for example, is an intracellular editing mechanism that is thought to generate isoforms for 30%-50% of all human genes. Whereas the inserted sequences usually display only minor structural rearrangements at the insertion site, recent observations indicate that they may also cause more dramatic structural displacements of adjacent structures. In the present study we test how artificially inserted sequences change the structure of the beta-sheet region in T4 lysozyme. Copies of two different beta-strands were inserted into two different loops of the beta-sheet, and the structures were determined. Not surprisingly, one insert "loops out" at its insertion site and forms a new small beta-hairpin structure. Unexpectedly, however, the second insertion leads to displacement of adjacent strands and a sequential reorganization of the beta-sheet topology. Even though the insertions were performed at two different sites, looping out occurred at the C-terminal end of the same beta-strand. Reasons as to why a non-native sequence would be recruited to replace that which occurs in the native protein are discussed. Our results illustrate how sequence insertions can facilitate protein evolution through both local and nonlocal changes in structure.

The Purkinje cell degeneration 5J mutation is a single amino acid insertion that destabilizes Nna1 protein

In the mouse, Purkinje cell degeneration (pcd) is a recessive mutation characterized by degeneration of cerebellar Purkinje cells, retinal photoreceptors, olfactory bulb mitral neurons, and certain thalamic neurons, and is accompanied by defective spermatogenesis. Previous studies of pcd have led to the identification of Nna1 as the causal gene; however, how loss of Nna1 function results in neurodegeneration remains unresolved. One useful approach for establishing which functional domains of a protein underlie a recessive phenotype has been to determine the genetic basis of the various alleles at the locus of interest. Because none of the pcd alleles analyzed at the time of the identification of Nna1 provided insight into the molecular basis of Nna1 loss-of-function, we obtained a recent pcd remutation--pcd5J, and after determining that its phenotype is comparable to existing pcd severe alleles, we sought its genetic basis by sequencing Nna1. In this article we report that pcd5J results from the insertion of a single GAC triplet encoding an aspartic acid residue at position 775 of Nna1. Although this insertion does not affect Nna1 expression at the RNA level, Nna1pcd-5J protein expression is markedly decreased. Pulse-chase experiments reveal that the aspartic acid insertion dramatically destabilizes Nna1pcd-5J protein, accounting for the observation that pcd5J is a severe allele. The presence of a readily detectable genetic mutation in pcd5J confirms that Nna1 loss-of-function alone underlies the broad pcd phenotype and will facilitate further studies of how Nna1 loss-of-function produces neurodegeneration and defective spermatogenesis in pcd mice.

Engineering Soluble Monomeric Streptavidin with Reversible Biotin Binding Capability

Monomeric streptavidin with reversible biotin binding capability has many potential applications. Because a complete biotin binding site in each streptavidin subunit requires the contribution of tryptophan 120 from a neighboring subunit, monomerization of the natural tetrameric streptavidin can generate streptavidin with reduced biotin binding affinity. Three residues, valine 55, threonine 76, and valine 125, were changed to either arginine or threonine to create electrostatic repulsion and steric hindrance at the interfaces. The double mutation (T76R,V125R) was highly effective to monomerize streptavidin. Because interfacial hydrophobic residues are exposed to solvent once tetrameric streptavidin is converted to the monomeric state, a quadruple mutein (T76R,V125R,V55T,L109T) was developed. The first two mutations are for monomerization, whereas the last two mutations aim to improve hydrophilicity at the interface to minimize aggregation. Monomerization was confirmed by four different approaches including gel filtration, dynamic light scattering, sensitivity to proteinase K, and chemical cross-linking. The quadruple mutein remained in the monomeric state at a concentration greater than 2 mg/ml. Its kinetic parameters for interaction with biotin suggest excellent reversible biotin binding capability, which enables the mutein to be easily purified on the biotin-agarose matrix. Another mutein (D61A,W120K) was developed based on two mutations that have been shown to be effective in monomerizing avidin. This streptavidin mutein was oligomeric in nature. This illustrates the importance in selecting the appropriate residues and approaches for effective monomerization of streptavidin.

A comprehensive and non-redundant database of protein domain movements

MOTIVATION

The current DynDom database of protein domain motions is a user-created database that suffers from selectivity and redundancy. The aim of the analysis presented here was to overcome both these limitations and to produce both a comprehensive and a non-redundant description of domain movements from structures stored in the current protein data bank.

RESULTS

A multi-step procedure is applied that starts with grouping proteins in the structural databank into families based on sequence similarity. Multiple sequence alignment, conformational clustering and a dimensional clustering method based on the Gram-Schmidt algorithm are applied to members of each family to remove dynamic redundancy in their domain movements. Representative domain movements are described in terms of domains, hinge axes and hinge-bending residues using the DynDom program. The results show that within an average family of 11.5 members, there are on average only 1.31 different domain movements indicating a high redundancy in the movements these structures represent. This verifies earlier findings that domain movements are usually highly controlled. Despite the removal of this considerable redundancy, the process has resulted in double the number of domain movements stored in the user-created database. The data are organized in a relational database with a web-interface.

AVAILABILITY

The database can be browsed and searched at http://www.cmp.uea.ac.uk/dyndom

CONTACT

sjh@cmp.uea.ac.uk.

Recent advances in our understanding of protein conformational stability from a pharmaceutical perspective

The marginal conformational stability of proteins has made them in some cases less than ideal candidates for pharmaceutical agents. Recent progress in our understanding of protein structure and stability has provided the opportunity to design the desired degree of stability into protein drug candidates. Modifications such as the optimisation of interior side-chain packing, the introduction of new ion-pairs, as well as the design of stabilising disulfide bridges and ligand binding sites, all offer the opportunity to produce proteins with enhanced stability properties.

Use of sequence duplication to engineer a ligand-triggered, long-distance molecular switch in T4 lysozyme

We have designed a molecular switch in a T4 lysozyme construct that controls a large-scale translation of a duplicated helix. As shown by crystal structures of the construct with the switch on and off, the conformational change is triggered by the binding of a ligand (guanidinium ion) to a site that in the wild-type protein was occupied by the guanidino head group of an Arg. In the design template, a duplicated helix is flanked by two loop regions of different stabilities. In the "on" state, the N-terminal loop is weakly structured, whereas the C-terminal loop has a well defined conformation that is stabilized by means of nonbonded interactions with the Arg head group. The truncation of the Arg to Ala destabilizes this loop and switches the protein to the "off" state, in which the duplicated helix is translocated approximately 20 A. Guanidinium binding restores the key interactions, restabilizes the C-terminal loop, and restores the "on" state. Thus, the presence of an external ligand, which is unrelated to the catalytic activity of the enzyme, triggers the inserted helix to translate 20 A away from the binding site. The results illustrate a proposed mechanism for protein evolution in which sequence duplication followed by point mutation can lead to the establishment of new function.

A Closed Binding Pocket and Global Destabilization Modify the Binding Properties of an Alternatively Spliced Form of the Second PDZ Domain of PTP-BL

PTP-BL is a large phosphatase that is implicated in cellular processes as diverse as cytokinesis, actin-cytoskeletal rearrangement, and apoptosis. Five PDZ domains mediate its cellular role by binding to the C termini of target proteins, forming multiprotein complexes. The second PDZ domain (PDZ2) binds to the C termini of the tumor suppressor protein APC and the LIM domain-containing protein RIL; however, in one splice variant, PDZ2as, a 5 residue insertion abrogates this binding. The insert causes distinct structural and dynamical changes in the alternatively spliced PDZ2: enlarging the L1 loop between beta2 and beta3, both lengthening and changing the orientation of the alpha2 helix, giving the base of the binding pocket less flexibility to accommodate ligands, and destabilizing the entire domain. These changes render the binding pocket incapable of binding C termini, possibly having implications in the functional role of PTP-BL.

Induction of substrate specificity shifts by placement of alanine insertions within the consensus amphipathic region of the Escherichia coli GABA (gamma-aminobutyric acid) transporter encoded by gabP

The Escherichia coli GABA (gamma-aminobutyric acid) permease GabP is a prototypical APC (amine/polyamine/choline) super-family transporter that has a CAR (consensus amphipathic region) containing multiple specificity determinants, ostensibly organized on two helical surfaces, one hydrophobic [SHS (sensitive hydrophobic surface)] and the other hydrophilic [SPS (sensitive polar surface)]. To gauge the functional effects of placing alanine insertions at close intervals across the entire GabP CAR, 64 insertion variants were constructed. Insertions, particularly those in the SHS and the SPS, were highly detrimental to steady-state [(3)H]GABA accumulation. TSR (transport specificity ratio) analysis, employing [(3)H]nipecotic acid and [(14)C]GABA, showed that certain alanine insertions were associated with a specificity shift (i.e. a change in k (cat)/ K (m)). An insertion (INS Ala-269) located N-terminal to the SHS increased specificity for [(3)H]nipecotic acid relative to [(14)C]GABA, whereas an insertion (INS Ala-321) located C-terminal to the SPS had the opposite effect. Overall, the results are consistent with a working hypothesis that the GabP CAR contains extensive functional surfaces that may be manipulated by insertion mutagenesis to alter the specificity ( k (cat)/ K (m)) phenotype. The thermodynamic basis of TSR analysis provides generality, suggesting that amino acid insertions could affect specificity in many other transporters, particularly those such as the E. coli phenylalanine permease PheP [Pi, Chow and Pittard (2002) J. Bacteriol. 184, 5842-5847] that have a functionally significant CAR-like domain.

Long-distance conformational changes in a protein engineered by modulated sequence duplication

There are few, if any, known instances in which a biological signal is transmitted via a large conformational change through the body of a protein. We describe here a mutant of T4 lysozyme that was engineered to permit structural change at a distance. The design uses a tandem sequence repeat that makes it possible to transmit large-scale structural changes from one end of an alpha-helix to the other over a distance of 17-25 A. The method should be of general applicability and may make it possible to introduce a mutation at one site in a protein that will induce large-scale changes in the structure at a spatially remote site.

SNPs on human chromosomes 21 and 22 -- analysis in terms of protein features and pseudogenes

SNPs are useful for genome-wide mapping and the study of disease genes. Previous studies have focused on SNPs in specific genes or SNPs pooled from a variety of different sources. Here, a systematic approach to the analysis of SNPs in relation to various features on a genome-wide scale, with emphasis on protein features and pseudogenes, is presented. We have performed a comprehensive analysis of 39,408 SNPs on human chromosomes 21 and 22 from the SNP consortium (TSC) database, where SNPs are obtained by random sequencing using consistent and uniform methods. Our study indicates that the occurrence of SNPs is lowest in exons and higher in repeats, introns and pseudogenes. Moreover, in comparing genes and pseudogenes, we find that the SNP density is higher in pseudogenes and the ratio of nonsynonymous to synonymous changes is also much higher. These observations may be explained by the increased rate of SNP accumulation in pseudogenes, which presumably are not under selective pressure. We have also performed secondary structure prediction on all coding regions and found that there is no preferential distribution of SNPs in a -helices, b -sheets or coils. This could imply that protein structures, in general, can tolerate a wide degree of substitutions. Tables relating to our results are available from http://genecensus.org/pseudogene.

K(+)-independent gastric H(+),K(+)-atpase activity. Dissociation of K(+)-independent dephosphorylation and preference for the E1 conformation by combined mutagenesis of transmembrane glutamate residues

Several mutations of residues Glu(795) and Glu(820) present in M5 and M6 of the catalytic subunit of gastric H(+),K(+)-ATPase have resulted in a K(+)-independent, SCH 28080-sensitive ATPase activity, caused by a high spontaneous dephosphorylation rate. The mutants with this property also have a preference for the E(1) conformation. This paper investigates the question of whether these two phenomena are coupled. This possibility was studied by combining mutations in residue Glu(343), present in M4, with those in residues 795 and 820. When in combined mutants Glu and/or Gln residues were present at positions 343, 795, and 820, the residue at position 820 dominated the behavior: a Glu giving K(+)-activated ATPase activity and an E(2) preference and a Gln giving K(+)-independent ATPase activity and an E(1) preference. With an Asp at position 343, the enzyme could be phosphorylated, but the dephosphorylation was blocked, independent of the presence of either a Glu or a Gln at positions 795 and 820. However, in these mutants, the direction of the E(2) <--> E(1) equilibrium was still dominated by the 820 residue: a Glu giving E(2) and a Gln giving E(1). This indicates that the preference for the E(1) conformation of the E820Q mutation is independent of an active dephosphorylation process.

Structural consequences of an amino acid deletion in the B1 domain of protein G

We describe the NMR structure of a deletion mutant of the B1 IgG-binding domain from Group G Streptococcus. The deletion occurs within the last beta-strand of the protein, where it may potentially have a deleterious effect on the stability of the protein if the protein were not able to conformationally adjust to the perturbation. In particular, the deletion changes the registry of the final three residues in the sheet, forcing a polar Thr to be buried in the interior of the protein and exposing a hydrophobic Val to solvent. The deletion could also potentially create a large cavity in the beta-sheet and force the alpha- and gamma-carboxylates of the C-terminal Glu residue into a partially buried region of the sheet. The structure of the mutant illustrates how the conformation of the protein adjusts to the deletion, thereby mitigating some of the potentially deleterious consequences. Although the elements of secondary structure are retained between the mutant and the wt domain, there are multiple small adjustments in the segments connecting secondary structure elements. In particular, a hydrogen bond between the Glu57 carboxylates and two main chain amides is introduced that alters the conformation in the loop connecting the helix to strand 3. In addition, to minimize hydrophobic surface exposure, the turn connecting strands 1 and 2 folds toward the core so that the molecular volume is decreased.

Exploring steric constraints on protein mutations using MAGE/PROBE

When planning a mutation to test some hypothesis, one crucial question is whether the new side chain is compatible with the existing structure; only if it is compatible can the interpretation of mutational results be straightforward. This paper presents a simple way of using the sensitive geometry of all-atom contacts (including hydrogens) to answer that question. The interactive MAGE/PROBE system lets the biologist explore conformational space for the mutant side chain, with an interactively updated kinemage display of its all-atom contacts to the original structure. The Autobondrot function in PROBE systematically explores that same conformational space, outputting contact scores at each point, which are then contoured and displayed. These procedures are applied here in two types of test cases, with known mutant structures. In ricin A chain, the ability of a neighboring glutamate to rescue activity of an active-site mutant is modeled successfully. In T4 lysozyme, six mutations to Leu are analyzed within the wild-type background structure, and their Autobondrot score maps correctly predict whether or not their surroundings must shift significantly in the actual mutant structures; interactive examination of contacts for the conformations involved explains which clashes are relieved by the motions. These programs are easy to use, are available free for UNIX or Microsoft Windows operating systems, and should be of significant help in choosing good mutation experiments or in understanding puzzling results.

Analysis of protein dimerization and ligand binding of orphan receptor HNF4alpha

Hepatocyte nuclear factor 4alpha (HNF4alpha) (NR2A1), an orphan member of the nuclear receptor superfamily, binds DNA exclusively as a homodimer even though it is very similar in amino acid sequence to retinoid X receptor alpha (RXRalpha), which heterodimerizes readily with other receptors. Here, experimental analysis of residues involved in protein dimerization and studies on a reported ligand for HNF4alpha are combined with a structural model of the HNF4alpha ligand-binding domain (LBD) (residues 137 to 384). When K300 (in helix 9) and E327 (in helix 10) of HNF4alpha1 were converted to the analogous residues in RXRalpha (E390 and K417, respectively) the resulting construct did not heterodimerize with the wild-type HNF4alpha, although it was still able to form homodimers and bind DNA. Furthermore, the double mutant did not heterodimerize with RXR or RAR but was still able to dimerize in solution with an HNF4alpha construct truncated at amino acid residue 268. This suggests that the charge compatibility between helices 9 and 10 is necessary, but not sufficient, to determine dimerization partners, and that additional residues in the HNF4alpha LBD are also important in dimerization. The structural model of the HNF4alpha LBD and an amino acid sequence alignment of helices 9 and 10 in various HNF4 and other receptor genes indicates that a K(X)(26)E motif can be used to identify HNF4 genes from other organisms and that a (E/D(X)(26-29)K/R) motif can be used to predict heterodimerization of many, but not all, receptors with RXR. In vitro analysis of another HNF4alpha mutant construct indicates that helix 10 also plays a structural role in the conformational integrity of HNF4alpha. The structural model and experimental analysis indicate that fatty acyl CoA thioesters, the proposed HNF4alpha ligands, are not good candidates for a traditional ligand for HNF4alpha. Finally, these results provide insight into the mechanism of action of naturally occurring mutations in the human HNF4alpha gene found in patients with maturity onset diabetes of the young 1 (MODY1).