Site-specific reflex response of ubiquitin to loop insertions

Predicting the structural effects of insertions in proteins by homology modeling remains a challenge. To investigate the molecular basis for conformational adaptations to insertions, ten mutants of ubiquitin were generated by introducing five different inserts, varying from five to 11 residues in size, at two different sites. Most insertion sequences were derived from homologous positions in structurally homologous ubiquitin-like proteins; to test sequence specificity, insertions were made into both homologous and non-homologous sites in ubiquitin. Structural inferences from NMR data suggest that each insertion site shows a reflex response to insertions: the sequence of the insertion has much less impact on structural adaptations than does the site of the insertion. Further, each site responds to insertions in a unique but consistent manner. For a given insertion site, different inserted sequences give rise to different stabilities, but the relationship between stability and sequence is not yet clear. However, the change in stability is similar for all insertions in a given site.     

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.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

Nucleotide sequences at the boundaries between gene and insertion regions in the rDNA of Drosophila melanogaster

Ribosomal RNA genes interrupted by type 1 insertions of 1 kb and 0.5 kb have been sequenced through the insertion region and compared with an uninterrupted gene. The 0.5 kb insertion is flanked by a duplication of a 14 bp segment that is present once in the uninterrupted gene; the 1 kb insertion is flanked by a duplication of 11 of these 14 bp. Short insertions are identical in their entire length to downstream regions of long insertions. No internal repeats occur in the insertion. The presence of target site duplications suggests that type 1 insertions arose by the introduction of transposable elements into rDNA. Short sequence homologies between the upstream ends of the insertions and the 28S' boundaries of the rRNA coding region suggest that short type 1 insertions may have arisen by recombination from longer insertions.We have sequenced both boundaries of two molecules containing type 2 insertions and the upstream boundary of a third; the points of interruption at the upstream boundary (28S' site) differ from each other in steps of 2 bp. Between the boundary in the 0.5 kb type 1 insertion and the type 2 boundaries there are distances of 74, 76, and 78 bp. At the downstream boundary (28S' site) the two sequenced type 2 insertions are identical. The rRNA coding region of one molecule extends across the insertion without deletion or duplication, but a 2 bp deletion in the RNA coding region is present in the second molecule. Stretches of 13 or 22 adenine residues occur at the downstream (28S') end of the two type 2 insertions.     

Comparison of the refined crystal structures of liganded and unliganded chicken, yeast and trypanosomal triosephosphate isomerase.

The refined crystal structures of chicken, yeast and trypanosomal triosephosphate isomerase (TIM) have been compared. TIM is known to exist in an "open" (unliganded) and "closed" (liganded) conformation. For chicken TIM only the refined open structure is available, whereas for yeast TIM and trypanosomal TIM refined structures of both the open and the closed structure have been used for this study. Comparison of these structures shows that the open structures of chicken TIM, yeast TIM and trypanosomal TIM are essentially identical. Also it is shown that the closed structures of yeast TIM and trypanosomal TIM are essentially identical. The conformational difference between the open and closed structures concerns a major shift (7 A) in loop-6. Minor shifts are observed in the two adjacent loops, loop-5 (1 A) and loop-7 (1 A). The pairwise comparison of the three different TIM barrels shows that the 105C alpha atoms of the core superimpose within 0.9 A. The sequences of these three TIMs have a pairwise sequence identity of approximately 50%. The residues that line the active site are 100% conserved. The residues interacting with each other across the dimer interface show extensive variability, but the direct hydrogen bonds between the two subunits are well conserved. The orientation of the two monomers with respect to each other is almost identical in the three different TIM structures. There are 56 (22%) conserved residues out of approximately 250 residues in 13 sequences. The functions of most of these conserved residues can be understood from the available open and closed structures of the three different TIMs. Some of these residues are quite far from the active site. For example, at a distance of 19 A from the active site there is a conserved saltbridge interaction between residues at the C-terminal ends of alpha-helix-6 and alpha-helix-7. This anchoring contrasts with the large conformational flexibility of loop-6 and loop-7 near the N termini of these helices. The flexibility of loop-6 is facilitated by a conserved large empty cavity near the N terminus of alpha-helix-6, which exists only in the open conformation.     

Accommodation of insertion mutations on the surface and in the interior of staphylococcal nuclease.

Alignment of homologous amino acid sequences reveals that insertion mutations are fairly common in evolution. Hitherto, the structural consequences of insertion mutations on the surface and in the interior of proteins of known structures have received little attention. We report here the high-resolution X-ray crystal structures of 2 site-directed insertion mutants of staphylococcal nuclease. The structure of the first insertion mutant, in which 2 glycine residues were inserted on the protein surface in the amino-terminal beta-strand, has been solved to 1.70 A resolution and refined to a crystallographic R value of 0.182. The inserted residues are accommodated in a special 3-residue beta-bulge. A bridging water molecule in the newly created cavity satisfies the hydrogen bonding requirements of the beta-sheet by forming a bifurcated hydrogen bond to 1 beta-strand, and a single hydrogen bond to the other beta-strand. The second insertion mutant contains a single leucine residue inserted at the end of the third beta-strand. The structure was solved to 2.0 A resolution and refined to a final R value of 0.196. The insertion is accommodated in a register shift that changes the conformation of the flexible loop portion of the molecule, relaxing and widening the omega turn. This structural alteration results in changes in position and coordination of a bound calcium ion important for catalysis. These structures illustrate important differences in how amino acid insertions are accommodated: as localized bulges, and as extensive register shifts.     

Crystal structure of phosphoserine phosphatase from Methanococcus jannaschii, a hyperthermophile, at 1.8 A resolution

BACKGROUND: D-Serine is a co-agonist of the N-methyl-D-aspartate subtype of glutamate receptors, a major neurotransmitter receptor family in mammalian nervous systems. D-Serine is converted from L-serine, 90% of which is the product of the enzyme phosphoserine phosphatase (PSP). PSP from M. jannaschii (MJ) shares significant sequence homology with human PSP. PSPs and P-type ATPases are members of the haloacid dehalogenase (HAD)-like hydrolase family, and all members share three conserved sequence motifs. PSP and P-type ATPases utilize a common mechanism that involves Mg(2+)-dependent phosphorylation and autodephosphorylation at an aspartyl side chain in the active site. The strong resemblance in sequence and mechanism implies structural similarity among these enzymes. RESULTS: The PSP crystal structure resembles the NAD(P) binding Rossmann fold with a large insertion of a four-helix-bundle domain and a beta hairpin. Three known conserved sequence motifs are arranged next to each other in space and outline the active site. A phosphate and a magnesium ion are bound to the active site. The active site is within a closed environment between the core alpha/beta domain and the four-helix-bundle domain. CONCLUSIONS: The crystal structure of MJ PSP was determined at 1.8 A resolution. Critical residues were assigned based on the active site structure and ligand binding geometry. The PSP structure is in a closed conformation that may resemble the phosphoserine bound state or the state after autodephosphorylation. Compared to a P-type ATPase (Ca(2+)-ATPase) structure, which is in an open state, this PSP structure appears also to be a good model for the closed conformation of P-type ATPase.     

Transfer RNA Genes Are Genomic Targets for De Novo Transposition of the Yeast Retrotransposon Ty3

Insertions of the yeast element Ty3 resulting from induced retrotransposition were characterized in order to identify the genomic targets of transposition. The DNA sequences of the junctions between Ty3 and flanking DNA were determined for two insertions of an unmarked element. Each insertion was at position -17 from the 5' end of a tRNA-coding sequence. Ninety-one independent insertions of a marked Ty3 element were studied by Southern blot analysis. Pairs of independent insertions into seven genomic loci accounted for 14 of these insertions. The DNA sequence flanking the insertion site was determined for at least one member of each pair of integrated elements. In each case, insertion was at position -16 or -17 relative to the 5' end of one of seven different tRNA genes. This proportion of genomic loci used twice for Ty3 integration is consistent with that predicted by a Poisson distribution for a number of genomic targets roughly equivalent to the estimated number of yeast tRNA genes. In addition, insertions upstream of the same tRNA gene in one case were at different positions, but in all cases were in the same orientation. Thus, genomic insertions of Ty3 in a particular orientation are apparently specified by the target, while the actual position of the insertion relative to the tRNA-coding sequence can vary slightly.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

Bent DNA is needed for recombinational enhancer activity in the site-specific recombination system Cin of bacteriophage P1. The role of FIS protein

A series of recombinational enhancer mutants was constructed by manipulating the ClaI site between the two FIS binding sites of the Hin enhancer. These mutants include insertions from two to 12 base-pairs and two deletions of one or two base-pairs. Recombinational enhancer activity was found only with four mutants carrying either a four base-pair substitution, ten base-pair insertions or a one base-pair deletion, respectively; two other ten base-pair insertion mutants, however, were inactive, although FIS protein binding was unaffected. So, besides binding of FIS protein to its specific sites within the enhancer sequence and the correct helical positioning of these sites on the DNA, another criterion for enhancer activity must be fulfilled. DNA bending assays identify this requirement as a change of the enhancer DNA conformation, which FIS protein is able to induce and to stabilize. This conformational change of the DNA can be blocked by mutations in the central segment between the two FIS binding sites of the Hin enhancer. This sequence has special functions for the recombinational enhancer activity.     

Characterization of the novel protein kinase activity present in the R1 subunit of herpes simplex virus ribonucleotide reductase.

We have compared the protein kinase activities of the R1 subunits from herpes simplex virus types 1 (HSV-1) and 2 (HSV-2) ribonucleotide reductase following expression in Escherichia coli. Autophosphorylation activity was observed when kinase assays were performed with immunoprecipitated R1 or proteins purified to homogeneity, and the activity was stimulated by the basic protein protamine. Transphosphorylation of histones or calmodulin by purified or immunoprecipitated HSV-1 and HSV-2 R1 was not observed, and our results suggest that the activities of these two proteins are similar. We further characterized the protein kinase activity of HSV-1 R1 by producing insertion and deletion mutants constructed with a plasmid expressing R1 amino acids 1 to 449. C-terminal deletion analysis identified the catalytic core of the enzyme as comprising residues 1 to 292, and this polypeptide will be useful for structural determinations by X-ray crystallography. Insertion of a 4-amino-acid sequence at sites within the protein kinase domain identified regions essential for activity; insertions at residues 22 and 112 completely inactivated activity, and an insertion at residue 136 reduced activity sixfold. Similar insertions at residues 257, 262, 292, and 343 had no effect on activity. The ATP analog 5'-fluorosulfonylbenzoyladenosine, which covalently modifies conventional eukaryotic kinases at an essential lysine residue within the active site, did label HSV R1, but this labelling occurred outside the N-terminal domain. These data indicate that the HSV R1 kinase is novel and distinct from other eukaryotic protein kinases.     

Multimodal Mechanism of Action for the Cdc34 Acidic Loop: A CASE STUDY FOR WHY UBIQUITIN-CONJUGATING ENZYMES HAVE LOOPS AND TAILS*

Together with ubiquitin ligases (E3), ubiquitin-conjugating enzymes (E2) are charged with the essential task of synthesizing ubiquitin chains onto protein substrates. Some 75% of the known E2s in the human proteome contain unique insertions in their primary sequences, yet it is largely unclear what effect these insertions impart on the ubiquitination reaction. Cdc34 is an important E2 with prominent roles in cell cycle regulation and signal transduction. The amino acid sequence of Cdc34 contains an insertion distal to the active site that is absent in most other E2s, yet this acidic loop (named for its four invariably conserved acidic residues) is critical for Cdc34 function both in vitro and in vivo. Here we have investigated how the acidic loop in human Cdc34 promotes ubiquitination, identifying two key molecular events during which the acidic loop exerts its influence. First, the acidic loop promotes the interaction between Cdc34 and its ubiquitin ligase partner, SCF. Second, two glutamic acid residues located on the distal side of the loop collaborate with an invariably conserved histidine on the proximal side of the loop to suppress the pK_a of an ionizing species on ubiquitin or Cdc34 which greatly contributes to Cdc34 catalysis. These results demonstrate that insertions can guide E2s to their physiologically relevant ubiquitin ligases as well as provide essential modalities that promote catalysis.     

A collection of sequenced and mapped Ds transposon insertion sites in Arabidopsis thaliana

Insertional mutagenesis is a powerful tool for generating knockout mutations that facilitate associating biological functions with as yet uncharacterized open reading frames (ORFs) identified by genomic sequencing or represented in EST databases. We have generated a collection of Dissociation (Ds) transposon lines with insertions on all 5 Arabidopsis chromosomes. Here we report the insertion sites in 260 independent single-transposon lines, derived from four different Ds donor sites. We amplified and determined the genomic sequence flanking each transposon, then mapped its insertion site by identity of the flanking sequences to the corresponding sequence in the Arabidopsis genome database. This constitutes the largest collection of sequence-mapped Ds insertion sites unbiased by selection against the donor site. Insertion site clusters have been identified around three of the four donor sites on chromosomes 1 and 5, as well as near the nucleolus organizers on chromosomes 2 and 4. The distribution of insertions between ORFs and intergenic sequences is roughly proportional to the ratio of genic to intergenic sequence. Within ORFs, insertions cluster near the translational start codon, although we have not detected insertion site selectivity at the nucleotide sequence level. A searchable database of insertion site sequences for the 260 transposon insertion sites is available at http://sgio2.biotec.psu.edu/sr. This and other collections of Arabidopsis lines with sequence-identified transposon insertion sites are a valuable genetic resource for functional genomics studies because the transposon location is precisely known, the transposon can be remobilized to generate revertants, and the Ds insertion can be used to initiate further local mutagenesis.     

IS1 insertion generates duplication of a nine base pair sequence at its target site.

Three independent integrations of the E. coli insertion sequence, IS1, into the gal operon have been analyzed. DNA sequences of portions of the wild-type galT gene which act as the target sites for these insertions, as well as the corresponding gal/IS1 junctions, are reported. Two features are particularly noteworthy. First, similar sequences appearing in inverted orientation consitute the ends of IS1: 18 of the terminal 23 base pairs at each end are identical. Second, in all three insertions, a 9 base pair segment found once in the wild-type sequence at the site of insertion is duplicated and appears in the same orientation at each end of the inserted element. The sequence of this 9 base pair repeat is different for each insertion analyzed. No homology between the inverted repeat sequences at the ends of IS1 and the sequences of the target sites is observed. Models for the mechanism of IS1 insertion are proposed.     

Roles of the species-specific subdomain and the N-terminal peptide of Toxoplasma gondii ferredoxin-NADP+ reductase in ferredoxin binding

The plant-type ferredoxin/ferredoxin-NADP(+) reductase (Fd/FNR) redox system found in parasites of the phylum Apicomplexa has been proposed as a target for novel drugs used against life-threatening diseases such as malaria and toxoplasmosis. Like many proteins from these protists, apicomplexan FNRs are characterized by the presence of unique peptide insertions of variable length and yet unknown function. Since three-dimensional data are not available for any of the parasite FNRs, we used limited proteolysis to carry out an extensive study of the conformation of Toxoplasma gondii FNR. This led to identification of 11 peptide bonds susceptible to the action of four different proteases. Cleavage sites are clustered in four regions of the enzyme, which include two of its three species-specific insertions. Such regions are thus predicted to form flexible surface loops. The protein substrate Fd protected FNR against cleavage both at its N-terminal peptide and at its largest sequence insertion (28 residues). Deletion by protein engineering of the species-specific subdomain containing the latter insertion resulted in an enzyme form that, although catalytically active, displayed a 10-fold decreased affinity for Fd. In contrast, removal of the first 15 residues of the enzyme unexpectedly enhanced its interaction with Fd. Thus, two flexible polypeptide regions of T. gondii FNR are involved in Fd interaction but have opposite roles in modulating the binding affinity for the protein ligand. In this respect, T. gondii FNR differs from plant FNRs, where the N-terminal peptide contributes to the stabilization of their complex with Fd.     

Structural insight into dynamic bypass of the major cisplatin‐DNA adduct by Y‐family polymerase Dpo4

Y‐family DNA polymerases bypass Pt‐GG, the cisplatin‐DNA double‐base lesion, contributing to the cisplatin resistance in tumour cells. To reveal the mechanism, we determined three structures of the Y‐family DNA polymerase, Dpo4, in complex with Pt‐GG DNA. The crystallographic snapshots show three stages of lesion bypass: the nucleotide insertions opposite the 3′G (first insertion) and 5′G (second insertion) of Pt‐GG, and the primer extension beyond the lesion site. We observed a dynamic process, in which the lesion was converted from an open and angular conformation at the first insertion to a depressed and nearly parallel conformation at the subsequent reaction stages to fit into the active site of Dpo4. The DNA translocation‐coupled conformational change may account for additional inhibition on the second insertion reaction. The structures illustrate that Pt‐GG disturbs the replicating base pair in the active site, which reduces the catalytic efficiency and fidelity. The in vivo relevance of Dpo4‐mediated Pt‐GG bypass was addressed by a dpo‐4 knockout strain of Sulfolobus solfataricus, which exhibits enhanced sensitivity to cisplatin and proteomic alterations consistent with genomic stress.     

Functional Fluorescent Protein Insertions in Herpes Simplex Virus gB Report on gB Conformation before and after Execution of Membrane Fusion

Entry of herpes simplex virus (HSV) into a target cell requires complex interactions and conformational changes by viral glycoproteins gD, gH/gL, and gB. During viral entry, gB transitions from a prefusion to a postfusion conformation, driving fusion of the viral envelope with the host cell membrane. While the structure of postfusion gB is known, the prefusion conformation of gB remains elusive. As the prefusion conformation of gB is a critical target for neutralizing antibodies, we set out to describe its structure by making genetic insertions of fluorescent proteins (FP) throughout the gB ectodomain. We created gB constructs with FP insertions in each of the three globular domains of gB. Among 21 FP insertion constructs, we found 8 that allowed gB to remain membrane fusion competent. Due to the size of an FP, regions in gB that tolerate FP insertion must be solvent exposed. Two FP insertion mutants were cell-surface expressed but non-functional, while FP insertions located in the crown were not surface expressed. This is the first report of placing a fluorescent protein insertion within a structural domain of a functional viral fusion protein, and our results are consistent with a model of prefusion HSV gB constructed from the prefusion VSV G crystal structure. Additionally, we found that functional FP insertions from two different structural domains could be combined to create a functional form of gB labeled with both CFP and YFP. FRET was measured with this construct, and we found that when co-expressed with gH/gL, the FRET signal from gB was significantly different from the construct containing CFP alone, as well as gB found in syncytia, indicating that this construct and others of similar design are likely to be powerful tools to monitor the conformation of gB in any model system accessible to light microscopy.     

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe2?-O? + (?)NO → trHbN-Fe3?-OH? + NO??) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (?)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (?)NO and O? to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.     

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.     

Energetics of charge-charge interactions between residues adjacent in sequence

Statistical analysis of the residue separation between a pair of ionizable side chains within 4 ? of each other was performed on a set of 1560 non-homologous PDB structures. We found that the frequency of pairs of like charges (i.e., pairs consisting of acidic residues Asp and Glu or pairs consisting of basic residues Arg and Lys) is two orders of magnitude lower than the pairs of oppositely charged residues (salt-bridges). We also found that for pairs of like charges the distribution is skewed dramatically towards short residue separation (<3). On the basis of these observations, we hypothesize that at short residue separation the repulsion between charges does not contribute much to the protein stability and the effects are largely dominated by the long range charge-charge interactions with other ionizable groups in the protein molecule. To test this hypothesis, we incorporated various pairs of charged residues at position 63 and 64 of ubiquitin and compared the stabilities of these variants. We also performed calculations of the expected changes in the charge-charge interactions. A very good correlation between experimental changes in the stability of ubiquitin variants, and changes in the energy of charge-charge interactions provides support for the hypothesis that a pair of ionizable residues next to each other in sequence modulates protein stability via long range charge-charge interactions with the rest of the protein.