Detecting folding motifs and similarities in protein structures

Structure of mitochondrial transcription termination factor 3 reveals a novel nucleic acid-binding domain

In mammalian cells, a family of mitochondrial transcription termination factors (MTERFs) regulates mitochondrial gene expression. MTERF family members share a approximately 270 residues long MTERF-domain required for DNA binding and transcription regulation. However, the structure of this widely conserved domain is unknown. Here, we show that the MTERF-domain of human MTERF3 forms a half-doughnut-shaped right-handed superhelix. The superhelix is built from alpha-helical tandem repeats that display a novel triangular three-helix motif. This repeat motif, which we denote the MTERF-motif, is a conserved structural element present in proteins from metazoans, plants, and protozoans. Furthermore, a narrow, strongly positively charged nucleic acid-binding path is found in the middle of the concave side of the half-doughnut. This arrangement suggests a half clamp nucleic acid-binding mode for MTERF-domains.

Crystal structure of the electron carrier domain of the reaction center cytochrome c(z) subunit from green photosynthetic bacterium Chlorobium tepidum

In green sulfur photosynthetic bacteria, the cytochrome c(z) (cyt c(z)) subunit in the reaction center complex mediates electron transfer mainly from menaquinol/cytochrome c oxidoreductase to the special pair (P840) of the reaction center. The cyt c(z) subunit consists of an N-terminal transmembrane domain and a C-terminal soluble domain that binds a single heme group. The periplasmic soluble domain has been proposed to be highly mobile and to fluctuate between oxidoreductase and P840 during photosynthetic electron transfer. We have determined the crystal structure of the oxidized form of the C-terminal functional domain of the cyt c(z) subunit (C-cyt c(z)) from thermophilic green sulfur bacterium Chlorobium tepidum at 1.3-A resolution. The overall fold of C-cyt c(z) consists of four alpha-helices and is similar to that of class I cytochrome c proteins despite the low similarity in their amino acid sequences. The N-terminal structure of C-cyt c(z) supports the swinging mechanism previously proposed in relation with electron transfer, and the surface properties provide useful information on possible interaction sites with its electron transfer partners. Several characteristic features are observed for the heme environment: These include orientation of the axial ligands with respect to the heme plane, surface-exposed area of the heme, positions of water molecules, and hydrogen-bond network involving heme propionate groups. These structural features are essential for elucidating the mechanism for regulating the redox state of cyt c(z).

The Mrs1 splicing factor binds the bI3 group I intron at each of two tetraloop-receptor motifs

Most large ribozymes require protein cofactors in order to function efficiently. The yeast mitochondrial bI3 group I intron requires two proteins for efficient splicing, Mrs1 and the bI3 maturase. Mrs1 has evolved from DNA junction resolvases to function as an RNA cofactor for at least two group I introns; however, the RNA binding site and the mechanism by which Mrs1 facilitates splicing were unknown. Here we use high-throughput RNA structure analysis to show that Mrs1 binds a ubiquitous RNA tertiary structure motif, the GNRA tetraloop-receptor interaction, at two sites in the bI3 RNA. Mrs1 also interacts at similar tetraloop-receptor elements, as well as other structures, in the self-folding Azoarcus group I intron and in the RNase P enzyme. Thus, Mrs1 recognizes general features found in the tetraloop-receptor motif. Identification of the two Mrs1 binding sites now makes it possible to create a model of the complete six-component bI3 ribonucleoprotein. All protein cofactors bind at the periphery of the RNA such that every long-range RNA tertiary interaction is stabilized by protein binding, involving either Mrs1 or the bI3 maturase. This work emphasizes the strong evolutionary pressure to bolster RNA tertiary structure with RNA-binding interactions as seen in the ribosome, spliceosome, and other large RNA machines.

Analyses of subnanometer resolution cryo-EM density maps

Today, electron cryomicroscopy (cryo-EM) can routinely achieve subnanometer resolutions of complex macromolecular assemblies. From a density map, one can extract key structural and functional information using a variety of computational analysis tools. At subnanometer resolution, these tools make it possible to isolate individual subunits, identify secondary structures, and accurately fit atomic models. With several cryo-EM studies achieving resolutions beyond 5Å, computational modeling and feature recognition tools have been employed to construct backbone and atomic models of the protein components directly from a density map. In this chapter, we describe several common classes of computational tools that can be used to analyze and model subnanometer resolution reconstructions from cryo-EM. A general protocol for analyzing subnanometer resolution density maps is presented along with a full description of steps used in analyzing the 4.3Å resolution structure of Mm-cpn.

BclxL changes conformation upon binding to wild-type but not mutant p53 DNA binding domain

p53 can induce apoptosis through mitochondrial membrane permeabilization by interaction of its DNA binding region with the anti-apoptotic proteins BclxL and Bcl2. However, little is known about the action of p53 at the mitochondria in molecular detail. By using NMR spectroscopy and fluorescence polarization we characterized the binding of wild-type and mutant p53 DNA binding domains to BclxL and show that the wild-type p53 DNA binding domain leads to structural changes in the BH3 binding region of BclxL, whereas mutants fail to induce such effects due to reduced affinity. This was probed by induced chemical shift and residual dipolar coupling data. These data imply that p53 partly achieves its pro-apoptotic function at the mitochondria by facilitating interaction between BclxL and BH3-only proteins in an allosteric mode of action. Furthermore, we characterize for the first time the binding behavior of Pifithrin-mu, a specific small molecule inhibitor of the p53-BclxL interaction, and present a structural model of the protein-ligand complex. A rather unusual behavior is revealed whereby Pifithrin-mu binds to both sides of the protein-protein complex. These data should facilitate the rational design of more potent specific BclxL-p53 inhibitors.

Exploring ribozyme conformational changes with X-ray crystallography

Relating three-dimensional fold to function is a central challenge in RNA structural biology. Toward this goal, X-ray crystallography has long been considered the "gold standard" for structure determinations at atomic resolution, although NMR spectroscopy has become a powerhouse in this arena as well. In the area of dynamics, NMR remains the dominant technique to probe the magnitude and timescales of molecular motion. Although the latter area remains largely unassailable by conventional crystallographic methods, inroads have been made on proteins using Laue radiation on timescales of ms to ns. Proposed 'fourth generation' radiation sources, such as free-electron X-ray lasers, promise ps- to fs-timescale resolution, and credible evidence is emerging that supports the feasibility of single molecule imaging. At present however, the preponderance of RNA structural information has been derived from timescale and motion insensitive crystallographic techniques. Importantly, developments in computing, automation and high-flux synchrotron sources have propelled the rapidity of 'conventional' RNA crystal structure determinations to timeframes of hours once a suitable set of phases is obtained. With a sufficient number of crystal structures, it is possible to create a structural ensemble that can provide insight into global and local molecular motion characteristics that are relevant to biological function. Here we describe techniques to explore conformational changes in the hairpin ribozyme, a representative non-protein-coding RNA catalyst. The approaches discussed include: (i) construct choice and design using prior knowledge to improve X-ray diffraction; (ii) recognition of long-range conformational changes and (iii) use of single-base or single-atom changes to create ensembles. The methods are broadly applicable to other RNA systems.

The first crystal structures of a family 19 class IV chitinase: the enzyme from Norway spruce

Chitinases help plants defend themselves against fungal attack, and play roles in other processes, including development. The catalytic modules of most plant chitinases belong to glycoside hydrolase family 19. We report here x-ray structures of such a module from a Norway spruce enzyme, the first for any family 19 class IV chitinase. The bi-lobed structure has a wide cleft lined by conserved residues; the most interesting for catalysis are Glu113, the proton donor, and Glu122, believed to be a general base that activate a catalytic water molecule. Comparisons to class I and II enzymes show that loop deletions in the class IV proteins make the catalytic cleft shorter and wider; from modeling studies, it is predicted that only three N-acetylglucosamine-binding subsites exist in class IV. Further, the structural comparisons suggest that the family 19 enzymes become more closed on substrate binding. Attempts to solve the structure of the complete protein including the associated chitin-binding module failed, however, modeling studies based on close relatives indicate that the binding module recognizes at most three N-acetylglucosamine units. The combined results suggest that the class IV enzymes are optimized for shorter substrates than the class I and II enzymes, or alternatively, that they are better suited for action on substrates where only small regions of chitin chain are accessible. Intact spruce chitinase is shown to possess antifungal activity, which requires the binding module; removing this module had no effect on measured chitinase activity.

Nucleoside monophosphate complex structures of the endonuclease domain from the influenza virus polymerase PA subunit reveal the substrate binding site inside the catalytic center

Highly pathogenic influenza virus strains currently in circulation pose a significant risk of a global pandemic. Following the reported crystal structure of the endonuclease domain from the avian influenza virus polymerase PA subunit, here we report the results of a systematic X-ray crystallographic analysis of its complex with adenosine, uridine, and thymidine nucleoside monophosphates (NMPs). Electron density corresponding to the monophosphate moiety of each nucleotide was apparent in each NMP complex and bound to the catalytic metal. A hydrophobic site was found to contribute to nucleoside binding. The NMP complex structures should represent the conformation of the bound product after nuclease cleavage. Moreover, one solvent molecule was found to occupy an equivalent position to the second reported Mn(2+) ion, where it mediates the interaction between bound NMPs and the N-terminal PA domain in the presence of the Mg(2+) ion. The results presented here indicate a possible cleavage mechanism and identify a distinct nucleotide binding pocket. The identification of this binding pocket opens a new avenue for anti-influenza drug discovery, targeting the cap-dependent endonuclease, in response to the worldwide threat of influenza.

Structure of Escherichia coli malate dehydrogenase at 1.45 A resolution

The structure of apo malate dehydrogenase from Escherichia coli has been determined to 1.45 A resolution. The crystals belonged to space group C2, with unit-cell parameters a = 146.0, b = 52.0, c = 168.9 A, beta = 102.2 degrees. The structure was determined with the molecular-replacement pipeline program BALBES and was refined to a final R factor of 18.6% (R(free) = 21.4%). The final model has two dimers in the asymmetric unit. In each dimer one monomer contains the active-site loop in the open conformation, whereas in the opposing monomer the active-site loop is disordered.

Structure-function analyses of isochorismate-pyruvate lyase from Pseudomonas aeruginosa suggest differing catalytic mechanisms for the two pericyclic reactions of this bifunctional enzyme

The isochorismate-pyruvate lyase from Pseudomonas aeruginosa (PchB) catalyzes two pericyclic reactions in a single active site. PchB physiologically produces salicylate and pyruvate from isochorismate for ultimate incorporation of the salicylate into the siderophore pyochelin. PchB also produces prephenate from chorismate, most likely due to structural homology to the Escherchia coli chorismate mutase. The molecular basis of catalysis among enzymatic pericyclic reactions is a matter of debate, one view holding that catalysis may be derived from electrostatic transition state stabilization and the opposing view that catalysis is derived from the generation of a reactive substrate conformation. Mutant forms of PchB were generated by site-directed mutagenesis at the site (K42) hypothesized to be key for electrostatic transition state stabilization (K42A, K42Q, K42E, and K42H). The loop containing K42 is mobile, and a mutant to slow loop dynamics was also designed (A43P). Finally, a previously characterized mutation (I87T) was also produced. Circular dichroism was used to assess the overall effect on secondary structure as a result of the mutations, and X-ray crystallographic structures are reported for K42A with salicylate and pyruvate bound and for apo-I87T. The data illustrate that the active site architecture is maintained in K42A-PchB, which indicates that differences in activity are not caused by secondary structural changes or by differences in active site loop conformation but rather by the chemical nature of this key residue. In contrast, the I87T structure demonstrates considerable mobility, suggesting that loop dynamics and conformational plasticity may be important for efficient catalysis. Finally, the mutational effects on k(cat) provide evidence that the two activities of PchB are not covariant and that a single hypothesis may not provide a sufficient explanation for catalysis.

Structural plasticity of malaria dihydroorotate dehydrogenase allows selective binding of diverse chemical scaffolds

Malaria remains a major global health burden and current drug therapies are compromised by resistance. Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) was validated as a new drug target through the identification of potent and selective triazolopyrimidine-based DHODH inhibitors with anti-malarial activity in vivo. Here we report x-ray structure determination of PfDHODH bound to three inhibitors from this series, representing the first of the enzyme bound to malaria specific inhibitors. We demonstrate that conformational flexibility results in an unexpected binding mode identifying a new hydrophobic pocket on the enzyme. Importantly this plasticity allows PfDHODH to bind inhibitors from different chemical classes and to accommodate inhibitor modifications during lead optimization, increasing the value of PfDHODH as a drug target. A second discovery, based on small molecule crystallography, is that the triazolopyrimidines populate a resonance form that promotes charge separation. These intrinsic dipoles allow formation of energetically favorable H-bond interactions with the enzyme. The importance of delocalization to binding affinity was supported by site-directed mutagenesis and the demonstration that triazolopyrimidine analogs that lack this intrinsic dipole are inactive. Finally, the PfDHODH-triazolopyrimidine bound structures provide considerable new insight into species-selective inhibitor binding in this enzyme family. Together, these studies will directly impact efforts to exploit PfDHODH for the development of anti-malarial chemotherapy.

Vimentin coil 1A-A molecular switch involved in the initiation of filament elongation

Interestingly, our previously published structure of the coil 1A fragment of the human intermediate filament protein vimentin turned out to be a monomeric alpha-helical coil instead of the expected dimeric coiled coil. However, the 39-amino-acid-long helix had an intrinsic curvature compatible with a coiled coil. We have now designed four mutants of vimentin coil 1A, modifying key a and d positions in the heptad repeat pattern, with the aim of investigating the molecular criteria that are needed to stabilize a dimeric coiled-coil structure. We have analysed the biophysical properties of the mutants by circular dichroism spectroscopy, analytical ultracentrifugation and X-ray crystallography. All four mutants exhibited an increased stability over the wild type as indicated by a rise in the melting temperature (T(m)). At a concentration of 0.1 mg/ml, the T(m) of the peptide with the single point mutation Y117L increased dramatically by 46 degrees C compared with the wild-type peptide. In general, the introduction of a single stabilizing point mutation at an a or a d position did induce the formation of a stable dimer as demonstrated by sedimentation equilibrium experiments. The dimeric oligomerisation state of the Y117L peptide was furthermore confirmed by X-ray crystallography, which yielded a structure with a genuine coiled-coil geometry. Most notably, when this mutation was introduced into full-length vimentin, filament assembly was completely arrested at the unit-length filament (ULF) level, both in vitro and in cDNA-transfected cultured cells. Therefore, the low propensity of the wild-type coil 1A to form a stable two-stranded coiled coil is most likely a prerequisite for the end-to-end annealing of ULFs into filaments. Accordingly, the coil 1A domains might "switch" from a dimeric alpha-helical coiled coil into a more open structure, thus mediating, within the ULFs, the conformational rearrangements of the tetrameric subunits that are needed for the intermediate filament elongation reaction.

Structure of the capsid amino-terminal domain from the betaretrovirus, Jaagsiekte sheep retrovirus

Jaagsiekte sheep retrovirus is a betaretrovirus and the causative agent of pulmonary adenocarcinoma, a transmissible lung tumour of sheep. Here we report the crystal structure of the capsid amino-terminal domain and examine the self-association properties of Jaagsiekte sheep retrovirus capsid. We find that the structure is remarkably similar to the amino-terminal domain of the alpharetrovirus, avian leukosis virus, revealing a previously undetected evolutionary similarity. Examination of capsid self-association suggests a mode of assembly not driven by the strong capsid carboxy-terminal domain interactions that characterise capsid assembly in the lentiviruses. Based on these data, we propose this structure provides a model for the capsid of betaretroviruses including the HML-2 family of endogenous human betaretroviruses.

Structure of YraM, a protein essential for growth of Haemophilus influenzae

Nontypeable Haemophilus influenzae is an obligate human parasite that often causes middle ear infections in children and exacerbates chronic obstructive pulmonary disorder, the fourth leading cause of death in the United States. There are no effective vaccines available for this strain. The lipoprotein YraM (gene HI1655) was identified as essential for the growth and viability of H. influenzae but its function is unknown. Sequence comparisons showed that YraM is a fusion of two protein modules. We grew crystals of the carboxyl-terminal module of YraM comprising residues 257-573 (YraM-C), phased the diffraction data by the multiwavelength anomalous diffraction technique, and refined the model to a crystallographic R-factor of 0.16 (R(free) = 0.19) with data to 1.35 A resolution. The two-domain structure of YraM-C adopts a fold similar to that observed for the open, unliganded forms of several periplasmic binding proteins (PBPs) involved in bacterial active transport. Sequence alignments of YraM homologues from other Gram-negative species showed that the most conserved residues of YraM-C cluster between the two domains in the location where other PBPs bind their cognate ligand. Modeling of YraM-C into a closed conformation similar to the leucine-bound form of the Leu/Ile/Val-binding protein (LIVBP) shows a putative binding pocket larger than the leucine-binding site in LIVBP. The pocket has both polar and nonpolar surfaces, with the latter located in the same area where a leucine side chain binds to LIVBP. We discuss possible biological functions of YraM considering its predicted location in the outer membrane, a novel place for such a binding protein.

The crystal structure of auracyanin A at 1.85 A resolution: the structures and functions of auracyanins A and B, two almost identical "blue" copper proteins, in the photosynthetic bacterium Chloroflexus aurantiacus

Auracyanins A and B are two closely similar "blue" copper proteins produced by the filamentous anoxygenic phototrophic bacterium Chloroflexus aurantiacus. Both proteins have a water-soluble 140-residue globular domain, which is preceded in the sequence by an N-terminal tail. The globular domains of auracyanins A and B have sequences that are 38% identical. The sequences of the N-terminal tails, on the other hand, are distinctly different, suggesting that auracyanins A and B occupy different membrane sites and have different functions. The crystal structure of auracyanin A has been solved and refined at 1.85 A resolution. The polypeptide fold is similar to that of auracyanin B (Bond et al. in J Mol Biol 306:47-67, 2001), but the distribution of charged and polar residues on the molecular surface is different. The Cu-site dimensions of the two auracyanins are identical. This is unexpected, since auracyanin A has a shorter polypeptide loop between two of the Cu-binding residues, and the two proteins have significantly different EPR, UV-visible and resonance Raman spectra. The genes for the globular domains of auracyanins A and B have been cloned in a bacterial expression system, enabling purification of large quantities of protein. It is shown that auracyanin A is expressed only when C. aurantiacus cells are grown in light, whereas auracyanin B is expressed under dark as well as light conditions. The inference is that auracyanin A has a function in photosynthesis, and that auracyanin B has a function in aerobic respiration.

Model-building strategies for low-resolution X-ray crystallographic data

Interpretation of low-resolution X-ray crystallographic data can prove to be a difficult task. The challenges faced in electron-density interpretation, the strategies that have been employed to overcome them and developments to automate the process are reviewed.

The interpretation of low-resolution X-ray crystallographic data proves to be challenging even for the most experienced crystallographer. Ambiguity in the electron-density map makes main-chain tracing and side-chain assignment difficult. However, the number of structures solved at resolutions poorer than 3.5 Å is growing rapidly and the structures are often of high biological interest and importance. Here, the challenges faced in electron-density interpretation, the strategies that have been employed to overcome them and developments to automate the process are reviewed. The methods employed in model generation from electron microscopy, which share many of the same challenges in providing high-confidence models of macromolecular structures and assemblies, are also considered.

Structure and proposed activity of a member of the VapBC family of toxin-antitoxin systems. VapBC-5 from Mycobacterium tuberculosis

In prokaryotes, cognate toxin-antitoxin pairs have long been known, but no three-dimensional structure has been available for any given complex from Mycobacterium tuberculosis. Here we report the crystal structure and activity of a member of the VapBC family of complexes from M. tuberculosis. The toxin VapC-5 is a compact, 150 residues, two domain alpha/beta protein. Bent around the toxin is the VapB-5 antitoxin, a 33-residue alpha-helix. Assays suggest that the toxin is an Mg-enabled endoribonuclease, inhibited by the antitoxin. The lack of DNase activity is consistent with earlier suggestions that the complex represses its own operon. Furthermore, analysis of the interactions in the binding of the antitoxin to the toxin suggest that exquisite control is required to protect the bacteria cell from toxic VapC-5.

Novel dimerization mode of the human Bcl-2 family protein Bak, a mitochondrial apoptosis regulator

Interactions of Bcl-2 family proteins play a regulatory role in mitochondrial apoptosis. The pro-apoptotic protein Bak resides in the outer mitochondrial membrane, and the formation of Bak homo- or heterodimers is involved in the regulation of apoptosis. The previously reported structure of the human Bak protein (residues Glu16-Gly186) revealed that a zinc ion was coordinated with two pairs of Asp160 and His164 residues from the symmetry-related molecules. This zinc-dependent homodimer was regarded as an anti-apoptotic dimer. In the present study, we determined the crystal structure of the human Bak residues Ser23-Asn185 at 2.5A, and found a distinct type of homodimerization through Cys166 disulfide bridging between the symmetry-related molecules. In the two modes of homodimerization, the molecular interfaces are completely different. In the membrane-targeted model of the S-S bridged dimer, the BH3 motifs are too close to the membrane to interact directly with the anti-apoptotic relatives, such as Bcl-x(L). Therefore, the Bak dimer structure reported here may represent a pro-apoptotic mode under oxidized conditions.

Ligand-induced conformational changes in a thermophilic ribose-binding protein

BACKGROUND

Members of the periplasmic binding protein (PBP) superfamily are involved in transport and signaling processes in both prokaryotes and eukaryotes. Biological responses are typically mediated by ligand-induced conformational changes in which the binding event is coupled to a hinge-bending motion that brings together two domains in a closed form. In all PBP-mediated biological processes, downstream partners recognize the closed form of the protein. This motion has also been exploited in protein engineering experiments to construct biosensors that transduce ligand binding to a variety of physical signals. Understanding the mechanistic details of PBP conformational changes, both global (hinge bending, twisting, shear movements) and local (rotamer changes, backbone motion), therefore is not only important for understanding their biological function but also for protein engineering experiments.

RESULTS

Here we present biochemical characterization and crystal structure determination of the periplasmic ribose-binding protein (RBP) from the hyperthermophile Thermotoga maritima in its ribose-bound and unliganded state. The T. maritima RBP (tmRBP) has 39% sequence identity and is considerably more resistant to thermal denaturation (app Tm value is 108 degrees C) than the mesophilic Escherichia coli homolog (ecRBP) (app Tm value is 56 degrees C). Polar ligand interactions and ligand-induced global conformational changes are conserved among ecRBP and tmRBP; however local structural rearrangements involving side-chain motions in the ligand-binding site are not conserved.

CONCLUSION

Although the large-scale ligand-induced changes are mediated through similar regions, and are produced by similar backbone movements in tmRBP and ecRBP, the small-scale ligand-induced structural rearrangements differentiate the mesophile and thermophile. This suggests there are mechanistic differences in the manner by which these two proteins bind their ligands and are an example of how two structurally similar proteins utilize different mechanisms to form a ligand-bound state.

The first structure of a glycoside hydrolase family 61 member, Cel61B from Hypocrea jecorina, at 1.6 A resolution

The glycoside hydrolase (GH) family 61 is a long-recognized, but still recondite, class of proteins, with little known about the activity, mechanism or function of its more than 70 members. The best-studied GH family 61 member, Cel61A of the filamentous fungus Hypocrea jecorina, is known to be an endoglucanase, but it is not clear if this represents the main activity or function of this family in vivo. We present here the first structure for this family, that of Cel61B from H. jecorina. The best-quality crystals were formed in the presence of nickel, and the crystal structure was solved to 1.6 A resolution using a single-wavelength anomalous dispersion method with nickel as the source of anomalous scatter. Cel61B lacks a carbohydrate-binding module and is a single-domain protein that folds into a twisted beta-sandwich. A structure-aided sequence alignment of all GH family 61 proteins identified a highly conserved group of residues on the surface of Cel61B. Within this patch of mostly polar amino acids was a site occupied by the intramolecular nickel hexacoordinately bound in the solved structure. In the Cel61B structure, there is no easily identifiable carbohydrate-binding cleft or pocket or catalytic center of the types normally seen in GHs. A structural comparison search showed that the known structure most similar to Cel61B is that of CBP21 from the Gram-negative soil bacterium Serratia marcescens, a member of the carbohydrate-binding module family 33 proteins. A polar surface patch highly conserved in that structural family has been identified in CBP21 and shown to be involved in chitin binding and in the protein's enhancement of chitinase activities. The analysis of the Cel61B structure is discussed in light of our continuing research to better understand the activities and function of GH family 61.

Nucleotide-binding flexibility in ultrahigh-resolution structures of the SRP GTPase Ffh

Two structures of the nucleotide-bound NG domain of Ffh, the GTPase subunit of the bacterial signal recognition particle (SRP), have been determined at ultrahigh resolution in similar crystal forms. One is GDP-bound and one is GMPPCP-bound. The asymmetric unit of each structure contains two protein monomers, each of which exhibits differences in nucleotide-binding conformation and occupancy. The GDP-bound Ffh NG exhibits two binding conformations in one monomer but not the other and the GMPPCP-bound protein exhibits full occupancy of the nucleotide in one monomer but only partial occupancy in the other. Thus, under the same solution conditions, each crystal reveals multiple binding states that suggest that even when nucleotide is bound its position in the Ffh NG active site is dynamic. Some differences in the positioning of the bound nucleotide may arise from differences in the crystal-packing environment and specific factors that have been identified include the relative positions of the N and G domains, small conformational changes in the P-loop, the positions of waters buried within the active site and shifts in the closing loop that packs against the guanine base. However, ;loose' binding may have biological significance in promoting facile nucleotide exchange and providing a mechanism for priming the SRP GTPase prior to its activation in its complex with the SRP receptor.

TOPS++FATCAT: fast flexible structural alignment using constraints derived from TOPS+ Strings Model

Background

Protein structure analysis and comparison are major challenges in structural bioinformatics. Despite the existence of many tools and algorithms, very few of them have managed to capture the intuitive understanding of protein structures developed in structural biology, especially in the context of rapid database searches. Such intuitions could help speed up similarity searches and make it easier to understand the results of such analyses.

Results

We developed a TOPS++FATCAT algorithm that uses an intuitive description of the proteins' structures as captured in the popular TOPS diagrams to limit the search space of the aligned fragment pairs (AFPs) in the flexible alignment of protein structures performed by the FATCAT algorithm. The TOPS++FATCAT algorithm is faster than FATCAT by more than an order of magnitude with a minimal cost in classification and alignment accuracy. For beta-rich proteins its accuracy is better than FATCAT, because the TOPS+ strings models contains important information of the parallel and anti-parallel hydrogen-bond patterns between the beta-strand SSEs (Secondary Structural Elements). We show that the TOPS++FATCAT errors, rare as they are, can be clearly linked to oversimplifications of the TOPS diagrams and can be corrected by the development of more precise secondary structure element definitions.

Software Availability

The benchmark analysis results and the compressed archive of the TOPS++FATCAT program for Linux platform can be downloaded from the following web site:

Conclusion

TOPS++FATCAT provides FATCAT accuracy and insights into protein structural changes at a speed comparable to sequence alignments, opening up a possibility of interactive protein structure similarity searches.

Vaccinia virus proteins A52 and B14 Share a Bcl-2-like fold but have evolved to inhibit NF-kappaB rather than apoptosis

Vaccinia virus (VACV), the prototype poxvirus, encodes numerous proteins that modulate the host response to infection. Two such proteins, B14 and A52, act inside infected cells to inhibit activation of NF-kappaB, thereby blocking the production of pro-inflammatory cytokines. We have solved the crystal structures of A52 and B14 at 1.9 A and 2.7 A resolution, respectively. Strikingly, both these proteins adopt a Bcl-2-like fold despite sharing no significant sequence similarity with other viral or cellular Bcl-2-like proteins. Unlike cellular and viral Bcl-2-like proteins described previously, A52 and B14 lack a surface groove for binding BH3 peptides from pro-apoptotic Bcl-2-like proteins and they do not modulate apoptosis. Structure-based phylogenetic analysis of 32 cellular and viral Bcl-2-like protein structures reveals that A52 and B14 are more closely related to each other and to VACV N1 and myxoma virus M11 than they are to other viral or cellular Bcl-2-like proteins. This suggests that a progenitor poxvirus acquired a gene encoding a Bcl-2-like protein and, over the course of evolution, gene duplication events have allowed the virus to exploit this Bcl-2 scaffold for interfering with distinct host signalling pathways.

The first structure of dipeptidyl-peptidase III provides insight into the catalytic mechanism and mode of substrate binding

Dipeptidyl-peptidases III (DPP III) are zinc-dependent enzymes that specifically cleave the first two amino acids from the N terminus of different length peptides. In mammals, DPP III is associated with important physiological functions and is a potential biomarker for certain types of cancer. Here, we present the 1.95-A crystal structure of yeast DPP III representing the prototype for the M49 family of metallopeptidases. It shows a novel fold with two domains forming a wide cleft containing the catalytic metal ion. DPP III exhibits no overall similarity to other metallopeptidases, such as thermolysin and neprilysin, but zinc coordination and catalytically important residues are structurally conserved. Substrate recognition is accomplished by a binding site for the N terminus of the peptide at an appropriate distance from the metal center and by a series of conserved arginine residues anchoring the C termini of different length substrates.

Structure-based function prediction: approaches and applications

The ever increasing number of protein structures determined by structural genomic projects has spurred much interest in the development of methods for structure-based function prediction. Existing methods can be roughly classified in two groups: some use a comparative approach looking for the presence of structural motifs possibly associated with a known biochemical function. Other methods try to identify functional patches on the surface of a protein using only its physicochemical characteristics. This review will cover both kinds of approaches to structure-based function prediction as well as their use in real-world cases. The main issues and limitations in using protein structure to predict function will also be discussed. These are mainly: the assessment of the statistical significance of structural similarities and the extent to which these methods depend on the accuracy and availability of structural data.

Insight from TonB hybrid proteins into the mechanism of iron transport through the outer membrane

We created hybrid proteins to study the functions of TonB. We first fused the portion of Escherichia coli tonB that encodes the C-terminal 69 amino acids (amino acids 170 to 239) of TonB downstream from E. coli malE (MalE-TonB69C). Production of MalE-TonB69C in tonB(+) bacteria inhibited siderophore transport. After overexpression and purification of the fusion protein on an amylose column, we proteolytically released the TonB C terminus and characterized it. Fluorescence spectra positioned its sole tryptophan (W213) in a weakly polar site in the protein interior, shielded from quenchers. Affinity chromatography showed the binding of the TonB C-domain to other proteins: immobilized TonB-dependent (FepA and colicin B) and TonB-independent (FepADelta3-17, OmpA, and lysozyme) proteins adsorbed MalE-TonB69C, revealing a general affinity of the C terminus for other proteins. Additional constructions fused full-length TonB upstream or downstream of green fluorescent protein (GFP). TonB-GFP constructs had partial functionality but no fluorescence; GFP-TonB fusion proteins were functional and fluorescent. The activity of the latter constructs, which localized GFP in the cytoplasm and TonB in the cell envelope, indicate that the TonB N terminus remains in the inner membrane during its biological function. Finally, sequence analyses revealed homology in the TonB C terminus to E. coli YcfS, a proline-rich protein that contains the lysin (LysM) peptidoglycan-binding motif. LysM structural mimicry occurs in two positions of the dimeric TonB C-domain, and experiments confirmed that it physically binds to the murein sacculus. Together, these findings infer that the TonB N terminus remains associated with the inner membrane, while the downstream region bridges the cell envelope from the affinity of the C terminus for peptidoglycan. This architecture suggests a membrane surveillance model of action, in which TonB finds occupied receptor proteins by surveying the underside of peptidoglycan-associated outer membrane proteins.

Three-dimensional structure of an intact glycoside hydrolase family 15 glucoamylase from Hypocrea jecorina

The three-dimensional structure of a complete Hypocrea jecorina glucoamylase has been determined at 1.8 A resolution. The presented structure model includes the catalytic and starch binding domains and traces the course of the 37-residue linker segment. While the structures of other fungal and yeast glucoamylase catalytic and starch binding domains have been determined separately, this is the first intact structure that allows visualization of the juxtaposition of the starch binding domain relative to the catalytic domain. The detailed interactions we see between the catalytic and starch binding domains are confirmed in a second independent structure determination of the enzyme in a second crystal form. This second structure model exhibits an identical conformation compared to the first structure model, which suggests that the H. jecorina glucoamylase structure we report is independent of crystal lattice contact restraints and represents the three-dimensional structure found in solution. The proposed starch binding regions for the starch binding domain are aligned with the catalytic domain in the three-dimensional structure in a manner that supports the hypothesis that the starch binding domain serves to target the glucoamylase at sites where the starch granular matrix is disrupted and where the enzyme might most effectively function.

The crystal structures of oxidized forms of human peroxiredoxin 5 with an intramolecular disulfide bond confirm the proposed enzymatic mechanism for atypical 2-Cys peroxiredoxins

Peroxiredoxin 5 (PRDX5) belongs to the PRDX superfamily of thiol-dependent peroxidases able to reduce hydrogen peroxide, alkyl hydroperoxides and peroxynitrite. PRDX5 is classified in the atypical 2-Cys subfamily of PRDXs. In this subfamily, the oxidized form of the enzyme is characterized by the presence of an intramolecular disulfide bridge between the peroxidatic and the resolving cysteine residues. We report here three crystal forms in which this intramolecular disulfide bond is indeed observed. The structures are characterized by the expected local unfolding of the peroxidatic loop, but also by the unfolding of the resolving loop. A new type of interface between PRDX molecules is described. The three crystal forms were not oxidized in the same way and the influence of the oxidizing conditions is discussed.

Towards atomic resolution structural determination by single-particle cryo-electron microscopy

Recent advances in cryo-electron microscopy and single-particle reconstruction (collectively referred to as 'cryoEM') have made it possible to determine the three-dimensional (3D) structures of several macromolecular complexes at near-atomic resolution ( approximately 3.8-4.5A). These achievements were accomplished by overcoming the challenges in sample handling, instrumentation, image processing, and model building. At near-atomic resolution, many detailed structural features can be resolved, such as the turns and deep grooves of helices, strand separation in beta sheets, and densities for loops and bulky amino acid side chains. Such structural data of the cytoplasmic polyhedrosis virus (CPV), the Epsilon 15 bacteriophage and the GroEL complex have provided valuable constraints for atomic model building using integrative tools, thus significantly enhancing the value of the cryoEM structures. The CPV structure revealed a drastic conformational change from a helix to a beta hairpin associated with RNA packaging and replication, coupling of RNA processing and release, and the long sought-after polyhedrin-binding domain. These latest advances in single-particle cryoEM provide exciting opportunities for the 3D structural determination of viruses and macromolecular complexes that are either too large or too heterogeneous to be investigated by conventional X-ray crystallography or nuclear magnetic resonance (NMR) methods.

The backbone structure of the thermophilic Thermoanaerobacter tengcongensis ribose binding protein is essentially identical to its mesophilic E. coli homolog

Background

Comparison of experimentally determined mesophilic and thermophilic homologous protein structures is an important tool for understanding the mechanisms that contribute to thermal stability. Of particular interest are pairs of homologous structures that are structurally very similar, but differ significantly in thermal stability.

Results

We report the X-ray crystal structure of a Thermoanaerobacter tengcongensis ribose binding protein (tteRBP) determined to 1.9 Å resolution. We find that tteRBP is significantly more stable (^appT_m value ~102°C) than the mesophilic Escherichia coli ribose binding protein (ecRBP) (^appT_m value ~56°C). The tteRBP has essentially the identical backbone conformation (0.41 Å RMSD of 235/271 C_α positions and 0.65 Å RMSD of 270/271 C_α positions) as ecRBP. Classification of the amino acid substitutions as a function of structure therefore allows the identification of amino acids which potentially contribute to the observed thermal stability of tteRBP in the absence of large structural heterogeneities.

Conclusion

The near identity of backbone structures of this pair of proteins entails that the significant differences in their thermal stabilities are encoded exclusively by the identity of the amino acid side-chains. Furthermore, the degree of sequence divergence is strongly correlated with structure; with a high degree of conservation in the core progressing to increased diversity in the boundary and surface regions. Different factors that may possibly contribute to thermal stability appear to be differentially encoded in each of these regions of the protein. The tteRBP/ecRBP pair therefore offers an opportunity to dissect contributions to thermal stability by side-chains alone in the absence of large structural differences.

Structural implications of Siglec-5-mediated sialoglycan recognition

Sialic acid (Sia) Ig-like binding lectins are important mediators of recognition and signaling events among myeloid cells. To investigate the molecular mechanism underlying sialic acid Ig-like lectin (Siglec) functions, we determined the crystal structure of the two N-terminal extracellular domains of human myeloid cell inhibitory receptor Siglec-5 (CD170) and its complexes with two sialylated carbohydrates. The native structure revealed an unusual conformation of the CC' ligand specificity loop and a unique interdomain disulfide bond. The alpha(2,3)- and alpha(2,6)-sialyllactose complexed structures showed a conserved Sia recognition motif that involves both Arg124 and a portion of the G-strand in the V-set domain forming beta-sheet-like hydrogen bonds with the glycerol side chain of the Sia. Only few protein contacts to the subterminal sugars are observed and mediated by the highly variable GG' linker and CC' loop. These structural observations, in conjunction with surface plasmon resonance binding assays, provide mechanistic insights into linkage-dependent Siglec carbohydrate recognition and suggest that Siglec-5 and other CD33-related Siglec receptors are more promiscuous in sialoglycan recognition than previously understood.

Structural analysis of the synaptic protein neuroligin and its beta-neurexin complex: determinants for folding and cell adhesion

The neuroligins are postsynaptic cell adhesion proteins whose associations with presynaptic neurexins participate in synaptogenesis. Mutations in the neuroligin and neurexin genes appear to be associated with autism and mental retardation. The crystal structure of a neuroligin reveals features not found in its catalytically active relatives, such as the fully hydrophobic interface forming the functional neuroligin dimer; the conformations of surface loops surrounding the vestigial active center; the location of determinants that are critical for folding and processing; and the absence of a macromolecular dipole and presence of an electronegative, hydrophilic surface for neurexin binding. The structure of a beta-neurexin-neuroligin complex reveals the precise orientation of the bound neurexin and, despite a limited resolution, provides substantial information on the Ca2+-dependent interactions network involved in trans-synaptic neurexin-neuroligin association. These structures exemplify how an alpha/beta-hydrolase fold varies in surface topography to confer adhesion properties and provide templates for analyzing abnormal processing or recognition events associated with autism.

Classification of protein folds

The diversity and complexity of bioinformatics tools currently available for protein sequence analysis can make it difficult to know where to begin when presented with a new sequence. In this article, we present a protocol outlining one approach to sequence analysis that should give as comprehensive a picture as possible as to the likely structure and function of a protein given the limits of available tools. We also provide worked examples showing how these tools can have an impact on the understanding of protein function prior to experimental studies.

Crystal structures of mammalian glutamine synthetases illustrate substrate-induced conformational changes and provide opportunities for drug and herbicide design

Glutamine synthetase (GS) catalyzes the ligation of glutamate and ammonia to form glutamine, with concomitant hydrolysis of ATP. In mammals, the activity eliminates cytotoxic ammonia, at the same time converting neurotoxic glutamate to harmless glutamine; there are a number of links between changes in GS activity and neurodegenerative disorders, such as Alzheimer's disease. In plants, because of its importance in the assimilation and re-assimilation of ammonia, the enzyme is a target of some herbicides. GS is also a central component of bacterial nitrogen metabolism and a potential drug target. Previous studies had investigated the structures of bacterial and plant GSs. In the present publication, we report the first structures of mammalian GSs. The apo form of the canine enzyme was solved by molecular replacement and refined at a resolution of 3 A. Two structures of human glutamine synthetase represent complexes with: a) phosphate, ADP, and manganese, and b) a phosphorylated form of the inhibitor methionine sulfoximine, ADP and manganese; these structures were refined to resolutions of 2.05 A and 2.6 A, respectively. Loop movements near the active site generate more closed forms of the eukaryotic enzymes when substrates are bound; the largest changes are associated with the binding of the nucleotide. Comparisons with earlier structures provide a basis for the design of drugs that are specifically directed at either human or bacterial enzymes. The site of binding the amino acid substrate is highly conserved in bacterial and eukaryotic GSs, whereas the nucleotide binding site varies to a much larger degree. Thus, the latter site offers the best target for specific drug design. Differences between mammalian and plant enzymes are much more subtle, suggesting that herbicides targeting GS must be designed with caution.

Structure of CrmE, a virus-encoded tumour necrosis factor receptor

Vaccinia virus (VACV), the smallpox vaccine, encodes many proteins that subvert the host immune response. One of these, cytokine response modifier E (CrmE), is secreted by infected cells and protects these cells from apoptotic challenge by tumour necrosis factor alpha (TNFalpha). We have expressed recombinant CrmE from VACV strain Lister in Escherichia coli, shown that the purified protein is monomeric in solution and competent to bind TNFalpha, and solved the structure to 2.0 A resolution. This is the first structure of a virus-encoded tumour necrosis factor receptor (TNFR). CrmE shares significant sequence similarity with mammalian type 2 TNF receptors (TNFSFR1B, p75; TNFR type 2). The structure confirms that CrmE adopts the canonical TNFR fold but only one of the two "ligand-binding" loops of TNFRSF1A is conserved in CrmE, suggesting a mechanism for the higher affinity of poxvirus TNFRs for TNFalpha over lymphotoxin-alpha. The roles of dimerisation and pre-ligand-assembly domains (PLADs) in poxvirus and mammalian TNFR activity are discussed.

Crystal structures of a family 19 chitinase from Brassica juncea show flexibility of binding cleft loops

Brassica juncea chitinase is an endo-acting, pathogenesis-related protein that is classified into glycoside hydrolase family 19, with highest homology (50-60%) in its catalytic domain to class I plant chitinases. Here we report X-ray structures of the chitinase catalytic domain from wild-type (apo, as well as with chloride ions bound) and a Glu234Ala mutant enzyme, solved by molecular replacement and refined at 1.53, 1.8 and 1.7 A resolution, respectively. Confirming our earlier mutagenesis studies, the active-site residues are identified as Glu212 and Glu234. Glu212 is believed to be the catalytic acid in the reaction, whereas Glu234 is thought to have a dual role, both activating a water molecule in its attack on the anomeric carbon, and stabilizing the charged intermediate. The molecules in the various structures differ significantly in the conformation of a number of loops that border the active-site cleft. The differences suggest an opening and closing of the enzyme during the catalytic cycle. Chitin is expected to dock first near Glu212, which will protonate it. Conformational changes then bring Glu234 closer, allowing it to assist in the following steps. These observations provide important insights into catalysis in family 19 chitinases.

Crystal Structures of a Poplar Thioredoxin Peroxidase that Exhibits the Structure of Glutathione Peroxidases: Insights into Redox-driven Conformational Changes

Glutathione peroxidases (GPXs) are a group of enzymes that regulate the levels of reactive oxygen species in cells and tissues, and protect them against oxidative damage. Contrary to most of their counterparts in animal cells, the higher plant GPX homologues identified so far possess cysteine instead of selenocysteine in their active site. Interestingly, the plant GPXs are not dependent on glutathione but rather on thioredoxin as their in vitro electron donor. We have determined the crystal structures of the reduced and oxidized form of Populus trichocarpaxdeltoides GPX5 (PtGPX5), using a selenomethionine derivative. PtGPX5 exhibits an overall structure similar to that of the known animal GPXs. PtGPX5 crystallized in the assumed physiological dimeric form, displaying a pseudo ten-stranded beta sheet core. Comparison of both redox structures indicates that a drastic conformational change is necessary to bring the two distant cysteine residues together to form an intramolecular disulfide bond. In addition, a computer model of a complex of PtGPX5 and its in vitro recycling partner thioredoxin h1 is proposed on the basis of the crystal packing of the oxidized form enzyme. A possible role of PtGPX5 as a heavy-metal sink is also discussed.

Crystal structure of bovine lipoyltransferase in complex with lipoyl-AMP

Lipoic acid is an essential cofactor of the alpha-ketoacid dehydrogenase complexes and the glycine cleavage system. It is covalently attached to a specific lysine residue of the subunit of the complexes. The bovine lipoyltransferase (bLT) catalyzes the lipoic acid attachment reaction using lipoyl-AMP as a substrate, forming a lipoylated protein and AMP. To gain insights into the reaction mechanism at the atomic level, we have determined the crystal structure of bLT at 2.10 A resolution. Unexpectedly, the purified recombinant bLT contains endogenous lipoyl-AMP. The structure of bLT consists of N-terminal and C-terminal domains, and lipoyl-AMP is bound to the active site in the N-terminal domain, adopting a U-shaped conformation. The lipoyl moiety is buried in the hydrophobic pocket, forming van der Waals interactions, and the AMP moiety forms numerous hydrogen bonds with bLT in another tunnel-like cavity. These interactions work together to expose the C10 atom of lipoyl-AMP to the surface of the bLT molecule. The carbonyl oxygen atom of lipoyl-AMP interacts with the invariant Lys135. The interaction might stimulate the positive charge of the C10 atom of lipoyl-AMP, and consequently facilitate the nucleophilic attack by the lysine residue of the lipoate-acceptor protein, accompanying the bond cleavage between the carbonyl group and the phosphate group. We discuss the structural differences between bLT and the lipoate-protein ligase A from Escherichia coli and Thermoplasma acidophilum. We further demonstrate that bLT in mitochondria also contains endogenous lipoylmononucleotide, being ready for the lipoylation of apoproteins.

X-ray structure of the T. aquaticus FtsY:GDP complex suggests functional roles for the C-terminal helix of the SRP GTPases

FtsY and Ffh are structurally similar prokaryotic Signal Recognition Particle GTPases that play an essential role in the Signal Recognition Particle (SRP)-mediated cotranslational targeting of proteins to the membrane. The two GTPases assemble in a GTP-dependent manner to form a heterodimeric SRP targeting complex. We report here the 2.1 A X-ray structure of FtsY from T. aquaticus bound to GDP. The structure of the monomeric protein reveals, unexpectedly, canonical binding interactions for GDP. A comparison of the structures of the monomeric and complexed FtsY NG GTPase domain suggests that it undergoes a conformational change similar to that of Ffh NG during the assembly of the symmetric heterodimeric complex. However, in contrast to Ffh, in which the C-terminal helix shifts independently of the other subdomains, the C-terminal helix and N domain of T. aquaticus FtsY together behave as a rigid body during assembly, suggesting distinct mechanisms by which the interactions of the NG domain "module" are regulated in the context of the two SRP GTPases.

Characterization of a C-C bond hydrolase from Sphingomonas wittichii RW1 with novel specificities towards polychlorinated biphenyl metabolites

Sphingomonas wittichii RW1 degrades chlorinated dibenzofurans and dibenzo-p-dioxins via meta cleavage. We used inverse PCR to amplify dxnB2, a gene encoding one of three meta-cleavage product (MCP) hydrolases identified in the organism that are homologues of BphD involved in biphenyl catabolism. Purified DxnB2 catalyzed the hydrolysis of 8-OH 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoate (HOPDA) approximately six times faster than for HOPDA at saturating substrate concentrations. Moreover, the specificity of DxnB2 for HOPDA (k(cat)/K(m) = 1.2 x 10(7) M(-1) s(-1)) was about half that of the BphDs of Burkholderia xenovorans LB400 and Rhodococcus globerulus P6, two potent polychlorinated biphenyl (PCB)-degrading strains. Interestingly, DxnB2 transformed 3-Cl and 4-OH HOPDAs, compounds that inhibit the BphDs and limit PCB degradation. DxnB2 had a higher specificity for 9-Cl HOPDA than for HOPDA but a lower specificity for 8-Cl HOPDA (k(cat)/K(m) = 1.7 x 10(6) M(-1) s(-1)), the chlorinated analog of 8-OH HOPDA produced during dibenzofuran catabolism. Phylogenetic analyses based on structure-guided sequence alignment revealed that DxnB2 belongs to a previously unrecognized class of MCP hydrolases, evolutionarily divergent from the BphDs although the physiological substrates of both enzyme types are HOPDAs. However, both classes of enzymes have mainly small hydrophobic residues lining the subsite that binds the C-6 phenyl of HOPDA, in contrast to the bulky hydrophobic residues (Phe106, Phe135, Trp150, and Phe197) found in the class II enzymes that prefer substrates possessing a C-6 alkyl. Thr196 and/or Asn203 appears to be an important determinant of specificity for DxnB2, potentially forming hydrogen bonds with the 8-OH substituent. This study demonstrates that the substrate specificities of evolutionarily divergent hydrolases may be useful for degrading mixtures of pollutants, such as PCBs.

Two lutein molecules in LHCII have different conformations and functions: Insights into the molecular mechanism of thermal dissipation in plants

When LHCII forms aggregates, the internal conformational changes will result in chlorophyll fluorescence quenching. Uncovering the molecular mechanism of this phenomenon will help us to understand how plants dissipate the excess excitation energy through non-photochemical quenching (NPQ) process. The crystal structure of spinach and pea LHCII have been published, and recently, we solved another crystal structure of LHCII from cucumber at 2.66A resolution. Here we present the first direct structural evidence indicating that the two lutein(Lut) molecules bound in each LHCII monomer have different conformations, Lut621 has a more twisted conformation than that of Lut620. The intimate interaction between the Lut620 and Chla612/Chla611 dimer leads to form a hetero-trimer, which is considered to be a potential quenching site. We also discovered that the dehydration of the LHCII crystals resulted in a notable shrinkage of the crystal unit cell dimensions which was accompanied by a red-shift of the fluorescence emission spectra of the crystals. These phenomena suggest the changes in the crystal packing during dehydration might be the cause of internal conformational changes within LHCII. We proposed a conformational change related NPQ model based on the structure analysis.

Identification of secondary structure elements in intermediate-resolution density maps

An increasing number of structural studies of large macromolecular complexes, both in X-ray crystallography and cryo-electron microscopy, have resulted in intermediate-resolution (5-10 A) density maps. Despite being limited in resolution, significant structural and functional information may be extractable from these maps. To aid in the analysis and annotation of these complexes, we have developed SSEhunter, a tool for the quantitative detection of alpha helices and beta sheets. Based on density skeletonization, local geometry calculations, and a template-based search, SSEhunter has been tested and validated on a variety of simulated and authentic subnanometer-resolution density maps. The result is a robust, user-friendly approach that allows users to quickly visualize, assess, and annotate intermediate-resolution density maps. Beyond secondary structure element identification, the skeletonization algorithm in SSEhunter provides secondary structure topology, which is potentially useful in leading to structural models of individual molecular components directly from the density.

Structure of the T109S mutant of Escherichia coli dihydroorotase complexed with the inhibitor 5-fluoroorotate: catalytic activity is reflected by the crystal form

Crystals of a single-point mutant (T109S) of Escherichia coli dihydroorotase (DHOase) with diminished activity grown in the presence of L-dihydroorotate (L-DHO) are tetragonal, with a monomer in the asymmetric unit. These crystals are extremely unstable and disintegrate shortly after formation, which is followed by the growth of orthorhombic crystals from the remnants of the tetragonal crystals or at new nucleation sites. Orthorhombic crystals, for which a structure has previously been reported [Thoden et al. (2001), Biochemistry, 40, 6989-6997; Lee et al. (2005), J. Mol. Biol. 348, 523-533], contain a dimer of DHOase in the asymmetric unit; the active site of one monomer contains the substrate N-carbamyl-L-aspartate (L-CA-asp) and the active site of the other monomer contains the product of the reaction, L-DHO. In the subunit with L-DHO in the active site, a surface loop (residues 105-115) is 'open'. In the other subunit, with L-CA-asp in the active site, the loop folds inwards, forming specific hydrogen bonds from the loop to the L-CA-asp. The tetragonal crystal form can be stabilized by crystallization in the presence of the inhibitor 5-fluoroorotate (FOA), a product (L-DHO) mimic. Crystals of the complex of T109S DHOase with FOA are tetragonal, space group P4(1)2(1)2, with unit-cell parameters a = b = 72.6, c = 176.1 A. The structure has been refined to R and R(free) values of 0.218 and 0.257, despite severe anisotropy of the diffraction. In this structure, the flexible loops are both in the 'open' conformation, which is consistent with FOA, like L-DHO, binding at both sites. The behaviour of the T109S mutant crystals of DHOase in the presence of L-DHO is explained by initial binding of L-DHO to both subunits, followed by slow conversion to L-CA-asp, with consequent movement of the flexible loop and dissolution of the crystals. Orthorhombic crystals are then able to grow in the presence of L-DHO and L-CA-asp.

EMatch: an efficient method for aligning atomic resolution subunits into intermediate-resolution cryo-EM maps of large macromolecular assemblies

Structural analysis of biological machines is essential for inferring their function and mechanism. Nevertheless, owing to their large size and instability, deciphering the atomic structure of macromolecular assemblies is still considered as a challenging task that cannot keep up with the rapid advances in the protein-identification process. In contrast, structural data at lower resolution is becoming more and more available owing to recent advances in cryo-electron microscopy (cryo-EM) techniques. Once a cryo-EM map is acquired, one of the basic questions asked is what are the folds of the components in the assembly and what is their configuration. Here, a novel knowledge-based computational method, named EMatch, towards tackling this task for cryo-EM maps at 6-10 A resolution is presented. The method recognizes and locates possible atomic resolution structural homologues of protein domains in the assembly. The strengths of EMatch are demonstrated on a cryo-EM map of native GroEL at 6 A resolution.

EMatch: discovery of high resolution structural homologues of protein domains in intermediate resolution cryo-EM maps

Cryo-EM has become an increasingly powerful technique for elucidating the structure, dynamics, and function of large flexible macromolecule assemblies that cannot be determined at atomic resolution. However, due to the relatively low resolution of cryo-EM data, a major challenge is to identify components of complexes appearing in cryo-EM maps. Here, we describe EMatch, a novel integrated approach for recognizing structural homologues of protein domains present in a 6-10 A resolution cryo-EM map and constructing a quasi-atomic structural model of their assembly. The method is highly efficient and has been successfully validated on various simulated data. The strength of the method is demonstrated by a domain assembly of an experimental cryo-EM map of native GroEL at 6 A resolution.