Design, expression, and purification of a Flaviviridae polymerase using a high-throughput approach to facilitate crystal structure determination

Bovine viral diarrhea virus (BVDV) nonstructural protein 5B is an RNA-dependent RNA polymerase, essential for viral replication. Initial attempts to crystallize a soluble form of the 695-residue BVDV polymerase did not produce any crystals. Limited proteolysis, homology modeling, and mutagenesis data were used to aid the design of polymerase constructs that might crystallize more readily. Limited proteolysis of the polymerase with trypsin identified a domain boundary within the protein. Homology modeling of the polymerase, based on the structure of hepatitis C virus polymerase, indicated that the two polymerases share a 23% identical "core," although overall sequence identity is low. Eighty-four expression clones of the BVDV polymerase were designed by fine-sampling of chain termini at the boundaries of domain and of active truncated forms of the polymerase. The resulting constructs were expressed in Escherichia coli and purified using high-throughput methods. Soluble truncated proteins were subjected to crystallization trials in a 96-well format, and two of these proteins were successfully crystallized.     

A dihydropyridine receptor alpha1s loop region critical for skeletal muscle contraction is intrinsically unstructured and binds to a SPRY domain of the type 1 ryanodine receptor

The II-III loop of the dihydropyridine receptor (DHPR) alpha(1s) subunit is a modulator of the ryanodine receptor (RyR1) Ca(2+) release channel in vitro and is essential for skeletal muscle contraction in vivo. Despite its importance, the structure of this loop has not been reported. We have investigated its structure using a suite of NMR techniques which revealed that the DHPR II-III loop is an intrinsically unstructured protein (IUP) and as such belongs to a burgeoning structural class of functionally important proteins. The loop does not possess a stable tertiary fold: it is highly flexible, with a strong N-terminal helix followed by nascent helical/turn elements and unstructured segments. Its residual structure is loosely globular with the N and C termini in close proximity. The unstructured nature of the II-III loop may allow it to easily modify its interaction with RyR1 following a surface action potential and thus initiate rapid Ca(2+) release and contraction. The in vitro binding partner for the II-III was investigated. The II-III loop interacts with the second of three structurally distinct SPRY domains in RyR1, whose function is unknown. This interaction occurs through two preformed N-terminal alpha-helical regions and a C-terminal hydrophobic element. The A peptide corresponding to the helical N-terminal region is a common probe of RyR function and binds to the same SPRY domain as the full II-III loop. Thus the second SPRY domain is an in vitro binding site for the II-III loop. The possible in vivo role of this region is discussed.     

A "loop entropy reduction" phage-display selection for folded amino acid sequences

As a step toward selecting folded proteins from libraries of randomized sequences, we have designed a 'loop entropy reduction'-based phage-display method. The basic premise is that insertion of a long disordered sequence into a loop of a host protein will substantially destabilize the host because of the entropic cost of closing a loop in a disordered chain. If the inserted sequence spontaneously folds into a stable structure with the N and C termini close in space, however, this entropic cost is diminished. The host protein function can, therefore, be used to select folded inserted sequences without relying on specific properties of the inserted sequence. This principle is tested using the IgG binding domain of protein L and the lck SH2 domain as host proteins. The results indicate that the loop entropy reduction screen is capable of discriminating folded from unfolded sequences when the proper host protein and insertion point are chosen.     

The structure of the mammalian signal recognition particle (SRP) receptor as prototype for the interaction of small GTPases with Longin domains

The eukaryotic signal recognition particle (SRP) and its receptor (SR) play a central role in co-translational targeting of secretory and membrane proteins to the endoplasmic reticulum. The SR is a heterodimeric complex assembled by the two GTPases SRalpha and SRbeta, which is membrane-anchored. Here we present the 2.45-A structure of mammalian SRbeta in its Mg2+ GTP-bound state in complex with the minimal binding domain of SRalpha termed SRX. SRbeta is a member of the Ras-GTPase superfamily closely related to Arf and Sar1, while SRX belongs to the SNARE-like superfamily with a fold also known as longin domain. SRX binds to the P loop and the switch regions of SRbeta-GTP. The binding mode and structural similarity with other GTPase-effector complexes suggests a co-GAP (GTPase-activating protein) function for SRX. Comparison with the homologous yeast structure and other longin domains reveals a conserved adjustable hydrophobic surface within SRX which is of central importance for the SRbeta-GTP:SRX interface. A helix swap in SRX results in the formation of a dimer in the crystal structure. Based on structural conservation we present the SRbeta-GTP:SRX structure as a prototype for conserved interactions in a variety of GTPase regulated targeting events occurring at endomembranes.     

Crystal structure of vesicular stomatitis virus matrix protein

The vesicular stomatitis virus (VSV) matrix protein (M) interacts with cellular membranes, self-associates and plays a major role in virus assembly and budding. We present the crystallographic structure, determined at 1.96 A resolution, of a soluble thermolysin resistant core of VSV M. The fold is a new fold shared by the other vesiculovirus matrix proteins. The structure accounts for the loss of stability of M temperature-sensitive mutants deficient in budding, and reveals a flexible loop protruding from the globular core that is important for self-assembly. Membrane floatation shows that, together with the M lysine-rich N-terminal peptide, a second domain of the protein is involved in membrane binding. Indeed, the structure reveals a hydrophobic surface located close to the hydrophobic loop and surrounded by conserved basic residues that may constitute this domain. Lastly, comparison of the negative-stranded virus matrix proteins with retrovirus Gag proteins suggests that the flexible link between their major membrane binding domain and the rest of the structure is a common feature shared by these proteins involved in budding and virus assembly.     

Characterization of the PR domain of RIZ1 histone methyltransferase

RIZ1 (PRDM2) and PRDI-BF1 (PRDM1) are involved in B cell differentiation and the development of B cell lymphomas. These proteins are expressed in two forms that differ by the presence or absence of a PR domain. The protein product that retains the PR domain is anti-tumorigenic while the product that lacks the PR domain is oncogenic and over-expressed in tumor cells. The conserved PR domain is homologous to the SET domain from a family of histone methyltransferases. RIZ1 is also a histone methyltransferase and methylates lysine 9 in histone H3. This activity has been mapped to the PR domain. In the present study, deuterium exchange mass spectrometry was used to define the structural boundaries of the RIZ1 PR domain and to map sites of missense mutations that occur in human cancers and reduce methyltransferase activity. Flexible segments were selectively deleted to produce protein products that crystallize for structural studies. Segments at the carboxyl terminus of the PR domain that are involved in methylation of H3 were shown to be flexible, similar to SET domains, suggesting that the PR and SET methyltransferases may belong to an emerging class of proteins that contain mobile functional regions.     

Crystal structure of the Shank PDZ-ligand complex reveals a class I PDZ interaction and a novel PDZ-PDZ dimerization

The Shank/proline-rich synapse-associated protein family of multidomain proteins is known to play an important role in the organization of synaptic multiprotein complexes. For instance, the Shank PDZ domain binds to the C termini of guanylate kinase-associated proteins, which in turn interact with the guanylate kinase domain of postsynaptic density-95 scaffolding proteins. Here we describe the crystal structures of Shank1 PDZ in its peptide free form and in complex with the C-terminal hexapeptide (EAQTRL) of guanylate kinase-associated protein (GKAP1a) determined at 1.8- and 2.25-A resolutions, respectively. The structure shows the typical class I PDZ interaction of PDZ-peptide complex with the consensus sequence -X-(Thr/Ser)-X-Leu. In addition, Asp-634 within the Shank1 PDZ domain recognizes the positively charged Arg at -1 position and hydrogen bonds, and salt bridges between Arg-607 and the side chains of the ligand at -3 and -5 positions contribute further to the recognition of the peptide ligand. Remarkably, whether free or complexed, Shank1 PDZ domains form dimers with a conserved beta B/beta C loop and N-terminal beta A strands, suggesting a novel model of PDZ-PDZ homodimerization. This implies that antiparallel dimerization through the N-terminal beta A strands could be a common configuration among PDZ dimers. Within the dimeric structure, the two-peptide binding sites are arranged so that the N termini of the bound peptide ligands are in close proximity and oriented toward the 2-fold axis of the dimer. This configuration may provide a means of facilitating dimeric organization of PDZ-target assemblies.     

Armadillo: domain boundary prediction by amino acid composition

The identification and annotation of protein domains provides a critical step in the accurate determination of molecular function. Both computational and experimental methods of protein structure determination may be deterred by large multi-domain proteins or flexible linker regions. Knowledge of domains and their boundaries may reduce the experimental cost of protein structure determination by allowing researchers to work on a set of smaller and possibly more successful alternatives. Current domain prediction methods often rely on sequence similarity to conserved domains and as such are poorly suited to detect domain structure in poorly conserved or orphan proteins. We present here a simple computational method to identify protein domain linkers and their boundaries from sequence information alone. Our domain predictor, Armadillo (http://armadillo.blueprint.org), uses any amino acid index to convert a protein sequence to a smoothed numeric profile from which domains and domain boundaries may be predicted. We derived an amino acid index called the domain linker propensity index (DLI) from the amino acid composition of domain linkers using a non-redundant structure dataset. The index indicates that Pro and Gly show a propensity for linker residues while small hydrophobic residues do not. Armadillo predicts domain linker boundaries from Z-score distributions and obtains 35% sensitivity with DLI in a two-domain, single-linker dataset (within +/-20 residues from linker). The combination of DLI and an entropy-based amino acid index increases the overall Armadillo sensitivity to 56% for two domain proteins. Moreover, Armadillo achieves 37% sensitivity for multi-domain proteins, surpassing most other prediction methods. Armadillo provides a simple, but effective method by which prediction of domain boundaries can be obtained with reasonable sensitivity. Armadillo should prove to be a valuable tool for rapidly delineating protein domains in poorly conserved proteins or those with no sequence neighbors. As a first-line predictor, domain meta-predictors could yield improved results with Armadillo predictions.     

Lysine methylation as a routine rescue strategy for protein crystallization

Crystallization remains a critical step in X-ray structure determination. Because it is not generally possible to rationally predict crystallization conditions, commercial screens have been developed which sample a wide range of crystallization space. While this approach has proved successful in many cases, a significant number of proteins fail to crystallize despite being soluble and monodispersed. It is established that chemical modification can facilitate the crystallization of otherwise intractable proteins. Here we describe a method for the reductive methylation of lysine residues which is simple, inexpensive, and efficient, and report on its application to ten proteins. We describe the effect of methylation on the physico-chemical properties of these proteins, and show that it led to diffraction-quality crystals from four proteins and structures for three that had hitherto proved refractory to crystallization. The method is suited to both low- and high-throughput laboratories.     

XANNpred: neural nets that predict the propensity of a protein to yield diffraction-quality crystals

Production of diffracting crystals is a critical step in determining the three-dimensional structure of a protein by X-ray crystallography. Computational techniques to rank proteins by their propensity to yield diffraction-quality crystals can improve efficiency in obtaining structural data by guiding both protein selection and construct design. XANNpred comprises a pair of artificial neural networks that each predict the propensity of a selected protein sequence to produce diffraction-quality crystals by current structural biology techniques. Blind tests show XANNpred has accuracy and Matthews correlation values ranging from 75% to 81% and 0.50 to 0.63 respectively; values of area under the receiver operator characteristic (ROC) curve range from 0.81 to 0.88. On blind test data XANNpred outperforms the other available algorithms XtalPred, PXS, OB-Score, and ParCrys. XANNpred also guides construct design by presenting graphs of predicted propensity for diffraction-quality crystals against residue sequence position. The XANNpred-SG algorithm is likely to be most useful to target selection in structural genomics consortia, while the XANNpred-PDB algorithm is more suited to the general structural biology community. XANNpred predictions that include sliding window graphs are freely available from http://www.compbio.dundee.ac.uk/xannpred     

The structural basis of mode of activation and functional diversity: a case study with HtrA family of serine proteases

HtrA (High temperature requirement protease A) proteins that are primarily involved in protein quality control belong to a family of serine proteases conserved from bacteria to humans. HtrAs are oligomeric proteins that share a common trimeric pyramidal architecture where each monomer comprises a serine protease domain and one or two PDZ domains. Although the overall structural integrity is well maintained and they exhibit similar mechanism of activation, subtle conformational changes and structural plasticity especially in the flexible loop regions and domain interfaces lead to differences in their active site conformation and hence in their specificity and functions.     

Fine-tuning of protein domain boundary by minimizing potential coiled coil regions

Structural determination of individual protein domains isolated from multidomain proteins is a common approach in the post-genomic era. Novel and thus uncharacterized domains liberated from intact proteins often self-associate due to incorrectly defined domain boundaries. Self-association results in missing signals, poor signal dispersion and a low signal-to-noise ratio in ¹H–¹⁵N HSQC spectra. We have found that a putative, non-canonical coiled coil region close to a domain boundary can cause transient hydrophobic self-association and monomer–dimer equilibrium in solution. Here we propose a rational method to predict putative coiled coil regions adjacent to the globular core domain using the program COILS. Except for the amino acid sequence, no preexisting knowledge concerning the domain is required. A small number of mutant proteins with a minimized coiled coil region have been rationally designed and tested. The engineered domains exhibit decreased self-association as assessed by ¹H–¹⁵N HSQC spectra with improved peak dispersion and sharper cross peaks. Two successful examples of isolating novel N-terminal domains from AAA-ATPases are demonstrated. Our method is useful for the experimental determination of domain boundaries suited for structural genomics studies. Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Scanning mutagenesis of regions in the Galpha protein Gpa1 that are predicted to interact with yeast mating pheromone receptors.

The mechanism by which receptors activate heterotrimeric G proteins was examined by scanning mutagenesis of the Saccharomyces cerevisiae pheromone-responsive Galpha protein (Gpa1). The juxtaposition of high-resolution structures for rhodopsin and its cognate G protein transducin predicted that at least six regions of Galpha are in close proximity to the receptor. Mutagenesis was targeted to residues in these domains in Gpa1, which included four loop regions (beta2-beta3, alpha2-beta4, alpha3-beta5, and alpha4-beta6) as well as the N and C termini. The mutants displayed a range of phenotypes from nonsignaling to constitutive activation of the pheromone pathway. The constitutive activity of some mutants could be explained by decreased production of Gpa1, which permits unregulated signaling by Gbetagamma. However, the constitutive activity caused by the F344C and E335C mutations in the alpha2-beta4 loop and F378C in the alpha3-beta5 loop was not due to decreased protein levels, and was apparently due to defects in sequestering Gbetagamma. The strongest loss of the function mutant, which was not detectably induced by a pheromone, was caused by a K314C substitution in the beta2-beta3 loop. Several other mutations caused weak signaling phenotypes. Altogether, these results suggest that residues in different interface regions of Galpha contribute to activation of signaling.     

The N-terminal domain (IF2N) of bacterial translation initiation factor IF2 is connected to the conserved C-terminal domains by a flexible linker

Bacterial translation initiation factor IF2 is a multidomain protein that is an essential component of a system for ensuring that protein synthesis begins at the correct codon within a messenger RNA. Full-length IF2 from Escherichia coli and seven fragments of the protein were expressed, purified, and characterized using nuclear magnetic resonance (NMR) and circular dichroism (CD) methods. Interestingly, resonances of the 6 kD IF2N domain located at the extreme N terminus of IF2 can be clearly identified within the NMR spectra of the full-length 97-kD protein. (15)N NMR relaxation rate data indicate that (1) the IF2N domain is internally well ordered and tumbles in solution in a manner that is independent of the other domains of the IF2 protein, and (2) the IF2N domain is connected to the C-terminal regions of IF2 by a flexible linker. Chemical shifts of resonances within the isolated IF2N domain do not significantly differ from those of the corresponding residues within the context of the full-length 97-kD protein, indicating that IF2N is a structurally independent unit that does not strongly interact with other regions of IF2. CD and NMR data together provide evidence that Domains I-III of IF2 have unstructured and flexible regions as well as substantial helical content; CD data indicate that the helical content of these regions decreases significantly at temperatures above 35 degrees C. The features of structurally well-ordered N- and C-terminal domains connected by a flexible linker with significant helical content are reminiscent of another translation initiation factor, IF3.     

A novel site on the Galpha -protein that recognizes heptahelical receptors

Specific domains of the G-protein alpha subunit have been shown to control coupling to heptahelical receptors. The extreme N and C termini and a region between alpha4 and alpha5 helices of the G-protein alpha subunit are known to determine selective interaction with the receptors. The metabotropic glutamate receptor 2 activated both mouse Galpha(15) and its human homologue Galpha(16), whereas metabotropic glutamate receptor 8 activated Galpha(15) only. The extreme C-terminal 20 amino acid residues are identical between the Galpha(15) and Galpha(16) and are therefore unlikely to be involved in coupling selectivity. Our data reveal two regions on Galpha(16) that inhibit its coupling to metabotropic glutamate receptor 8. On a three-dimensional model, both regions are found in a close proximity to the extreme C terminus of Galpha(16). One module comprises alpha4 helix, alpha4-beta6 loop (L9 Loop), beta6 sheet, and alpha5 helix. The other, not described previously, is located within the loop that links the N-terminal alpha helix to the beta1 strand of the Ras-like domain of the alpha subunit. Coupling of Galpha(16) protein to the metabotropic glutamate receptor 8 is partially modulated by each module alone, whereas both modules are needed to eliminate the coupling fully.     

Structure, dynamics and binding characteristics of the second PDZ domain of PTP-BL

The PDZ domains of the protein tyrosine phosphatase PTP-BL mediate interactions by binding to specific amino acid sequences in target proteins. The solution structure of the second PDZ domain of PTP-BL, PDZ2, displays a compact fold with six β strands and two α-helices. A unique feature of this domain compared to the canonical PDZ fold is an extended flexible loop at the base of the binding pocket, termed L1, that folds back onto the protein backbone, a feature that is shared by both the murine and human orthologues. The structure of PDZ2 differs significantly from the orthologous human structure. A comparison of structural quality indicators clearly demonstrates that the PDZ2 ensemble is statistically more reasonable than that of the human orthologue. The analysis of ¹⁵N relaxation data for PDZ2 shows a normal pattern, with more rigid secondary structures and more flexible loop structures. Close to the binding pocket, Leu85 and Thr88 display greater mobility when compared to surrounding residues. Peptide binding studies demonstrated a lack of interaction between murine PDZ2 and the C terminus of the murine Fas/CD95 receptor, suggesting that the Fas/CD95 receptor is not an in vivo target for PDZ2. In addition, PDZ2 specifically binds the C termini of both human Fas/CD95 receptor and the RIL protein, despite RIL containing a non-canonical PDZ-interacting sequence of E-x-V. A model of PDZ2 with the RIL peptide reveals that the PDZ2 binding pocket is able to accommodate the bulkier side-chain of glutamic acid while maintaining crucial protein to peptide hydrogen bond interactions.     

Structure and dynamics of ASC2, a pyrin domain-only protein that regulates inflammatory signaling

Pyrin domain (PYD)-containing proteins are key components of pathways that regulate inflammation, apoptosis, and cytokine processing. Their importance is further evidenced by the consequences of mutations in these proteins that give rise to autoimmune and hyperinflammatory syndromes. PYDs, like other members of the death domain (DD) superfamily, are postulated to mediate homotypic interactions that assemble and regulate the activity of signaling complexes. However, PYDs are presently the least well characterized of all four DD subfamilies. Here we report the three-dimensional structure and dynamic properties of ASC2, a PYD-only protein that functions as a modulator of multidomain PYD-containing proteins involved in NF-kappaB and caspase-1 activation. ASC2 adopts a six-helix bundle structure with a prominent loop, comprising 13 amino acid residues, between helices two and three. This loop represents a divergent feature of PYDs from other domains with the DD fold. Detailed analysis of backbone 15N NMR relaxation data using both the Lipari-Szabo model-free and reduced spectral density function formalisms revealed no evidence of contiguous stretches of polypeptide chain with dramatically increased internal motion, except at the extreme N and C termini. Some mobility in the fast, picosecond to nanosecond timescale, was seen in helix 3 and the preceding alpha2-alpha3 loop, in stark contrast to the complete disorder seen in the corresponding region of the NALP1 PYD. Our results suggest that extensive conformational flexibility in helix 3 and the alpha2-alpha3 loop is not a general feature of pyrin domains. Further, a transition from complete disorder to order of the alpha2-alpha3 loop upon binding, as suggested for NALP1, is unlikely to be a common attribute of pyrin domain interactions.     

Membrane topology of gamma-secretase component PEN-2

PEN-2 is an integral membrane protein that is a necessary component of the gamma-secretase complex, which is central in the pathogenesis of Alzheimer's disease and is also required for Notch signaling. In the absence of PEN-2, Notch signaling fails to guide normal development in Caenorhabditis elegans, and amyloid beta peptide is not generated from the amyloid precursor protein. Human PEN-2 is a 101-amino acid protein containing two putative transmembrane domains. To understand its interaction with other gamma-secretase components, it is important to know the membrane topology of each member of the complex. To characterize the membrane topology of PEN-2, we introduced single amino acid changes in each of the three hydrophilic regions of PEN-2 to generate N-linked glycosylation sites. We found that the N-linked glycosylation sites present in the N- and C-terminal domains of PEN-2 were utilized, whereas a site in the hydrophilic "loop" region connecting the two transmembrane domains was not. The addition of a carbohydrate structure in the N-terminal domain of PEN-2 prevented association with presenilin 1, whereas glycosylation in the C-terminal region of PEN-2 did not, suggesting that the N-terminal domain is important for interactions with presenilin 1. Immunofluorescence microscopy with selective permeabilization of the plasma membrane of cells expressing epitope-tagged forms of PEN-2 confirmed the lumenal location of both the N and C termini. A protease protection assay also demonstrated that the loop domain of PEN-2 is cytosolic. Thus, PEN-2 spans the membrane twice, with the N and C termini facing the lumen of the endoplasmic reticulum.     

Crystal contact-free conformation of an intrinsically flexible loop in protein crystal: Tim21 as the case study

BackgroundIn protein crystals, flexible loops are frequently deformed by crystal contacts, whereas in solution, the large motions result in the poor convergence of such flexible loops in NMR structure determinations. We need an experimental technique to characterize the structural and dynamic properties of intrinsically flexible loops of protein molecules.MethodsWe designed an intended crystal contact-free space (CCFS) in protein crystals, and arranged the flexible loop of interest in the CCFS. The yeast Tim 21 protein was chosen as the model protein, because one of the loops (loop 2) is distorted by crystal contacts in the conventional crystal.ResultsYeast Tim21 was fused to the MBP protein by a rigid α-helical linker. The space created between the two proteins was used as the CCFS. The linker length provides adjustable freedom to arrange loop 2 in the CCFS. We re-determined the NMR structure of yeast Tim21, and conducted MD simulations for comparison. Multidimensional scaling was used to visualize the conformational similarity of loop 2. We found that the crystal contact-free conformation of loop 2 is located close to the center of the ensembles of the loop 2 conformations in the NMR and MD structures.ConclusionsLoop 2 of yeast Tim21 in the CCFS adopts a representative, dominant conformation in solution.General significanceNo single powerful technique is available for the characterization of flexible structures in protein molecules. NMR analyses and MD simulations provide useful, but incomplete information. CCFS crystallography offers a third route to this goal.