Alignment of protein structures in the presence of domain motions

Crystal structure of Leptospira LSS_01692 reveals a dimeric structure and induces inflammatory responses through Toll-like receptor 2-dependent NF-κB and MAPK signal transduction pathways

Leptospirosis is a commonly overlooked zoonotic disease that occurs in tropical and subtropical regions. Recent studies have divided the Leptospira spp. into three groups based on virulence, including pathogenic, intermediate, and saprophytic species. Pathogenic species express a protein family with leucine-rich repeat (LRR) domains, which are less expressed or absent in nonpathogenic species, highlighting the importance of this protein family in leptospirosis. However, the role of LRR domain proteins in the pathogenesis of leptospirosis is still unknown and requires further investigation. In this study, the 3D structure of LSS_01692 (rLRR38) was obtained using X-ray crystallography at a resolution of 3.2 Å. The results showed that rLRR38 forms a typical horseshoe structure with 11 α-helices and 11 β-sheets and an antiparallel dimeric structure. The interactions of rLRR38 with extracellular matrix and cell surface receptors were evaluated using ELISA and single-molecule atomic force microscopy. The results showed that rLRR38 interacted with fibronectin, collagen IV, and Toll-like receptor 2 (TLR2). Incubating HK2 cells with rLRR38 induced two downstream inflammation responses (IL-6 and MCP-1) in the TLR2 signal transduction pathway. The TLR2-TLR1 complex showed the most significant upregulation effects under rLRR38 treatment. Inhibitors also significantly inhibited nuclear factor κB and mitogen-activated protein kinases signals transduction under rLRR38 stimulation. In conclusion, rLRR38 was determined to be a novel LRR domain protein in 3D structure and demonstrated as a TLR2-binding protein that induces inflammatory responses. These structural and functional studies provide a deeper understanding of the pathogenesis of leptospirosis.

Identification of dual active sites in Caenorhabditis elegans GANA-1 protein: an ortholog of the human α-GAL a and α-NAGA enzymes

Fabry disease (FD) is caused by a defective α-galactosidase A (α-GAL A) enzyme responsible for breaking down globotriaosylceramide (Gb₃). To develop affordable therapeutics, more effort is needed to obtain insights into the underlying mechanism of FD and understanding human α-GAL A structure and function in related animal models. We adopted C. elegans as a model to elucidate the sequence and 3D structure of its GANA-1 enzyme and compared it to human α-GAL A. We constructed GANA-1 3D structure by homology modelling and validated the quality of the predicted GANA-1 structure, followed by computational docking of human ligands. The GANA-1 protein shared sequence similarities up to 42.1% with the human α-GAL A in silico and had dual active sites. GANA-1 homology modelling showed that 11 out of 13 amino acids in the first active site of GANA-1 protein overlapped with the human α-GAL A active site, indicating the prospect for substrate cross-reaction. Computational molecular docking using human ligands like Gb₃ (first pocket)_, 4-nitrophenyl-α-D-galactopyranoside (second pocket), α-galactose (second pocket), and N-acetyl-D-galactosamine (second pocket) showed negative binding energy. This revealed that the ligands were able to bind within both GANA-1 active sites, mimicking the human α-GAL A and α-NAGA enzymes. We identified human compounds with adequate docking scores, predicting robust interactions with the GANA-1 active site. Our data suggested that the C. elegans GANA-1 enzyme may possess structural and functional similarities to human α-GAL A, including an intrinsic capability to metabolize Gb₃ deposits.Communicated by Ramaswamy H. Sarma.

Structure and function of the Ecto-Nucleotide Pyrophosphatase-Phosphodiesterase (ENPP) family: tidying up diversity

Ecto-nucleotide pyrophosphatase/phosphodiesterase (ENPP) family members (ENPP1–7) have been implicated in key biological and pathophysiological processes, including nucleotide and phospholipid signaling, bone mineralization, fibrotic diseases, and tumor-associated immune cell infiltration. ENPPs are single-pass transmembrane ecto-enzymes, with notable exceptions of ENPP2 (Autotaxin) and ENNP6, which are secreted and glycosylphosphatidylinositol (GPI)-anchored, respectively. ENNP1 and ENNP2 are the best characterized and functionally the most interesting members. Here, we review the structural features of ENPP1–7 to understand how they evolved to accommodate specific substrates and mediate different biological activities. ENPPs are defined by a conserved phosphodiesterase (PDE) domain. In ENPP1–3, the PDE domain is flanked by two N-terminal somatomedin B-like domains and a C-terminal inactive nuclease domain that confers structural stability, whereas ENPP4–7 only possess the PDE domain. Structural differences in the substrate-binding site endow each protein with unique characteristics. Thus, ENPP1, ENPP3, ENPP4, and ENPP5 hydrolyze nucleotides, whereas ENPP2, ENPP6, and ENNP7 evolved as phospholipases through adaptions in the catalytic domain. These adaptations explain the different biological and pathophysiological functions of individual members. Understanding the ENPP members as a whole advances our insights into common mechanisms, highlights their functional diversity, and helps to explore new biological roles.

N-terminal residues are crucial for quaternary structure and active site conformation for the phosphoserine aminotransferase from enteric human parasite E. histolytica

Phosphoserine aminotransferase (PSAT) is a pyridoxal-5'phosphate (PLP)-dependent enzyme that catalyzes the second reversible step in the phosphoserine biosynthetic pathway producing serine. The crystal structure of E. histolytica PSAT (EhPSAT) complexed with PLP was elucidated at 3.0 Å resolution and the structures of its mutants, EhPSAT_Δ45 and EhPSAT_Δ4, at 1.8 and 2.4 Å resolution respectively. Deletion of 45 N-terminal residues (EhPSAT_Δ45) resulted in an inactive protein, the structure showed a dimeric arrangement drastically different from that of the wild-type protein, with the two monomers translated and rotated by almost 180° with respect to each other; causing a rearrangement of the active site to which PLP was unable to bind. Deletion of first N-terminal 15 (EhPSAT_Δ15) and four 11th to 14th residues (EhPSAT_Δ4) yielded up to 98% and 90% decrease in the activity respectively. Absence of aldimine linkage between PLP-Lys in the crystal structure of EhPSAT_Δ4 mutant explains for such decrease in activity and describes the importance of these N-terminal residues. Furthermore, a halide-binding site was found in close proximity to the active site. A stretch of six amino acids (146-NNTIYG-151) only conserved in the Entamoeba genus, contributes to halide binding may explain that the halide inhibition could be specific to Entamoeba.

A Single Mutation Traps a Half-Sites Reactive Enzyme in Midstream, Explaining Asymmetry in Hydride Transfer

In Escherichia coli thymidylate synthase (EcTS), rate-determining hydride transfer from the cofactor 5,10-methylene-5,6,7,8-tetrahydrofolate to the intermediate 5-methylene-2'-deoxyuridine 5'-monophosphate occurs by hydrogen tunneling, requiring precise alignment of reactants and a closed binding cavity, sealed by the C-terminal carboxyl group. Mutations that destabilize the closed conformation of the binding cavity allow small molecules such as β-mercaptoethanol (β-ME) to enter the active site and compete with hydride for addition to the 5-methylene group of the intermediate. The C-terminal deletion mutant of EcTS produced the β-ME adduct in proportions that varied dramatically with cofactor concentration, from 50% at low cofactor concentrations to 0% at saturating cofactor conditions, suggesting communication between active sites. We report the 2.4 Å X-ray structure of the C-terminal deletion mutant of E. coli TS in complex with a substrate and a cofactor analogue, CB3717. The structure is asymmetric, with reactants aligned in a manner consistent with hydride transfer in only one active site. In the second site, CB3717 has shifted to a site where the normal cofactor would be unlikely to form 5-methylene-2'-deoxyuridine 5'-monophosphate, consistent with no formation of the β-ME adduct. The structure shows how the binding of the cofactor at one site triggers hydride transfer and borrows needed stabilization from substrate binding at the second site. It indicates pathways through the dimer interface that contribute to allostery relevant to half-sites reactivity.

NAHAL-Flex: A Numerical and Alphabetical Hinge Detection Algorithm for Flexible Protein Structure Alignment

Flexible proteins are proteins that have conformational changes in their structures. Protein flexibility analysis is critical for classifying and understanding protein functionality. For that analysis, the hinge areas where proteins show flexibility must be detected. To detect the location of the hinges, previous methods have utilized the three-dimensional (3D) structure of proteins, which is highly computational. To reduce the computational complexity, this study proposes a novel text-based method using structural alphabets (SAs) for detecting the hinge position, called NAHAL-Flex. Protein structures were encoded to a particular type of SA called the protein folding shape code (PFSC), which remains unaffected by location, scale, and rotation. The flexible regions of the proteins are the only places in which letter sequences can be distorted. With this knowledge, it is possible to find the longest alignment path of two letter sequences using a dynamic programming (DP) algorithm. Then, the proposed method looks for regions where the alphabet sequence is distorted to find the most probable hinge positions. In order to reduce the number of hinge positions, a genetic algorithm (GA) was utilized to find the best candidate hinge points. To evaluate the method's effectiveness, four different flexible and rigid protein databases, including two small datasets and two large datasets, were utilized. For the small dataset, the NAHAL-Flex method was comparable to state-of-the-art structural flexible alignment methods. The result for the large datasets show that NAHAL-Flex outperforms some well-known alignment methods, e.g., DaliLite, Matt, DeepAlign, and TM-align; the speed of NAHAL-Flex was faster and its result was more accurate than the other methods.

Crystal structure of lipoate-bound lipoate ligase 1, LipL1, from Plasmodium falciparum

Plasmodium falciparum lipoate protein ligase 1 (PfLipL1) is an ATP-dependent ligase that belongs to the biotin/lipoate A/B protein ligase family (PFAM PF03099). PfLipL1 is the only known canonical lipoate ligase in Pf and functions as a redox switch between two lipoylation routes in the parasite mitochondrion. Here, we report the crystal structure of a deletion construct of PfLipL1 (PfLipL1_Δ243-279 ) bound to lipoate, and validate the lipoylation activity of this construct in both an in vitro lipoylation assay and a cell-based lipoylation assay. This characterization represents the first step in understanding the redox dependence of the lipoylation mechanism in malaria parasites. Proteins 2017; 85:1777-1783. © 2017 Wiley Periodicals, Inc.

Structural insight into β-Clamp and its interaction with DNA Ligase in Helicobacter pylori

Helicobacter pylori, a gram-negative and microaerophilic bacterium, is the major cause of chronic gastritis, gastric ulcers and gastric cancer. Owing to its central role, DNA replication machinery has emerged as a prime target for the development of antimicrobial drugs. Here, we report 2Å structure of β-clamp from H. pylori (Hpβ-clamp), which is one of the critical components of DNA polymerase III. Despite of similarity in the overall fold of eubacterial β-clamp structures, some distinct features in DNA interacting loops exists that have not been reported previously. The in silico prediction identified the potential binders of β-clamp such as alpha subunit of DNA pol III and DNA ligase with identification of β-clamp binding regions in them and validated by SPR studies. Hpβ-clamp interacts with DNA ligase in micromolar binding affinity. Moreover, we have successfully determined the co-crystal structure of β-clamp with peptide from DNA ligase (not reported earlier in prokaryotes) revealing the region from ligase that interacts with β-clamp.

Calculating and scoring high quality multiple flexible protein structure alignments

MOTIVATION

Calculating multiple protein structure alignments (MSAs) is important for understanding functional and evolutionary relationships between protein families, and for modeling protein structures by homology. While incorporating backbone flexibility promises to circumvent many of the limitations of rigid MSA algorithms, very few flexible MSA algorithms exist today. This article describes several novel improvements to the Kpax algorithm which allow high quality flexible MSAs to be calculated. This article also introduces a new Gaussian-based MSA quality measure called 'M-score', which circumvents the pitfalls of RMSD-based quality measures.

RESULTS

As well as calculating flexible MSAs, the new version of Kpax can also score MSAs from other aligners and from previously aligned reference datasets. Results are presented for a large-scale evaluation of the Homstrad, SABmark and SISY benchmark sets using Kpax and Matt as examples of state-of-the-art flexible aligners and 3DCOMB as an example of a state-of-the-art rigid aligner. These results demonstrate the utility of the M-score as a measure of MSA quality and show that high quality MSAs may be achieved when structural flexibility is properly taken into account.

AVAILABILITY AND IMPLEMENTATION

Kpax 5.0 may be downloaded for academic use at http://kpax.loria.fr/

CONTACT

dave.ritchie@inria.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structural investigation and inhibitory response of halide on phosphoserine aminotransferase from Trichomonas vaginalis

BACKGROUND

Phosphoserine aminotransferase (PSAT) catalyses the second reversible step of the phosphoserine biosynthetic pathway in Trichomonas vaginalis, which is crucial for the synthesis of serine and cysteine.

METHODS

PSAT from T. vaginalis (TvPSAT) was analysed using X-ray crystallography, enzyme kinetics, and molecular dynamics simulations.

RESULTS

The crystal structure of TvPSAT was determined to 2.15Å resolution, and is the first protozoan PSAT structure to be reported. The active site of TvPSAT structure was found to be in a closed conformation, and at the active site PLP formed an internal aldimine linkage to Lys 202. In TvPSAT, Val 340 near the active site while it is Arg in most other members of the PSAT family, might be responsible in closing the active site. Kinetic studies yielded Km values of 54 μM and 202 μM for TvPSAT with OPLS and AKG, respectively. Only iodine inhibited the TvPSAT activity while smaller halides could not inhibit.

CONCLUSION

Results from the structure, comparative molecular dynamics simulations, and the inhibition studies suggest that iodine is the only halide that can bind TvPSAT strongly and may thus inhibit the activity of TvPSAT. The long loop between β8 and α8 at the opening of the TvPSAT active site cleft compared to other PSATs, suggests that this loop may help control the access of substrates to the TvPSAT active site and thus influences the enzyme kinetics.

GENERAL SIGNIFICANCE

Our structural and functional studies have improved our understanding of how PSAT helps this organism persists in the environment.

The Structural and Functional Characterization of Mammalian ADP-dependent Glucokinase

The enzyme-catalyzed phosphorylation of glucose to glucose-6-phosphate is a reaction central to the metabolism of all life. ADP-dependent glucokinase (ADPGK) catalyzes glucose-6-phosphate production, utilizing ADP as a phosphoryl donor in contrast to the more well characterized ATP-requiring hexokinases. ADPGK is found in Archaea and metazoa; in Archaea, ADPGK participates in a glycolytic role, but a function in most eukaryotic cell types remains unknown. We have determined structures of the eukaryotic ADPGK revealing a ribokinase-like tertiary fold similar to archaeal orthologues but with significant differences in some secondary structural elements. Both the unliganded and the AMP-bound ADPGK structures are in the "open" conformation. The structures reveal the presence of a disulfide bond between conserved cysteines that is positioned at the nucleotide-binding loop of eukaryotic ADPGK. The AMP-bound ADPGK structure defines the nucleotide-binding site with one of the disulfide bond cysteines coordinating the AMP with its main chain atoms, a nucleotide-binding motif that appears unique to eukaryotic ADPGKs. Key amino acids at the active site are structurally conserved between mammalian and archaeal ADPGK, and site-directed mutagenesis has confirmed residues essential for enzymatic activity. ADPGK is substrate inhibited by high glucose concentration and shows high specificity for glucose, with no activity for other sugars, as determined by NMR spectroscopy, including 2-deoxyglucose, the glucose analogue used for tumor detection by positron emission tomography.

Crystal Structure of the Human tRNA m(1)A58 Methyltransferase-tRNA(3)(Lys) Complex: Refolding of Substrate tRNA Allows Access to the Methylation Target

Human tRNA3(Lys) is the primer for reverse transcription of HIV; the 3' end is complementary to the primer-binding site on HIV RNA. The complementarity ends at the 18th base, A58, which in tRNA3(Lys) is modified to remove Watson-Crick pairing. Motivated to test the role of the modification in terminating the primer-binding sequence and thus limiting run-on transcription, we asked how the modification of RNA could be accomplished. tRNA m(1)A58 methyltransferase (m(1)A58 MTase) methylates N1 of A58, which is buried in the TΨC-loop of tRNA, from cofactor S-adenosyl-L-methionine. This conserved tRNA modification is essential for stability of initiator tRNA in Saccharomyces cerevisiae. Reported here, three structures of human tRNA m(1)A58 MTase in complex with human tRNA3(Lys) and the product S-adenosyl-L-homocysteine show a dimer of heterodimers in which each heterodimer comprises a catalytic chain, Trm61, and a homologous but noncatalytic chain, Trm6, repurposed as a tRNA-binding subunit that acts in trans; tRNAs bind across the dimer interface such that Trm6 from the opposing heterodimer brings A58 into the active site of Trm61. T-loop and D-loop are splayed apart showing how A58, normally buried in tRNA, becomes accessible for modification. This result has broad impact on our understanding of the mechanisms of modifying internal sites in folded tRNA. The structures serve as templates for design of inhibitors that could be used to test tRNA m(1)A58 MTase's impact on retroviral priming and transcription.

Crystal structures and kinetics of Type III 3-phosphoglycerate dehydrogenase reveal catalysis by lysine

D-Phosphoglycerate dehydrogenase (PGDH) catalyzes the first committed step of the phosphorylated serine biosynthesis pathway. Here, we report for the first time, the crystal structures of Type IIIK PGDH from Entamoeba histolytica in the apo form, as well as in complexes with substrate (3-phosphoglyceric acid) and cofactor (NAD(+) ) to 2.45, 1.8 and 2.2 Å resolution, respectively. Comparison of the apo structure with the substrate-bound structure shows that the substrate-binding domain is rotated by ~ 20° to close the active-site cleft. The cofactor-bound structure also shows a closed-cleft conformation, in which NAD(+) is bound to the nucleotide-binding domain and a formate ion occupies the substrate-binding site. Superposition of the substrate- and cofactor-bound structures represents a snapshot of the enzyme in the active form, where C2 of the substrate and C4N of the cofactor are 2.2 Å apart, and the amino group of Lys263 is close enough to the substrate to remove the proton from the hydroxyl group of PGA, indicating the role of Lys in the catalysis. Mutation of Lys263 to Ala yields just 0.8% of the specific activity of the wild-type enzyme, revealing that Lys263 indeed plays an integral role in the catalytic activity. The detectable activity of the mutant, however, indicates that after 20° rotation of the substrate-binding domain, the resulting positions of the substrate and cofactor are sufficiently close to make a productive reaction.

Conformation-independent structural comparison of macromolecules with ProSMART

The identification and exploration of (dis)similarities between macromolecular structures can help to gain biological insight, for instance when visualizing or quantifying the response of a protein to ligand binding. Obtaining a residue alignment between compared structures is often a prerequisite for such comparative analysis. If the conformational change of the protein is dramatic, conventional alignment methods may struggle to provide an intuitive solution for straightforward analysis. To make such analyses more accessible, the Procrustes Structural Matching Alignment and Restraints Tool (ProSMART) has been developed, which achieves a conformation-independent structural alignment, as well as providing such additional functionalities as the generation of restraints for use in the refinement of macromolecular models. Sensible comparison of protein (or DNA/RNA) structures in the presence of conformational changes is achieved by enforcing neither chain nor domain rigidity. The visualization of results is facilitated by popular molecular-graphics software such as CCP4mg and PyMOL, providing intuitive feedback regarding structural conservation and subtle dissimilarities between close homologues that can otherwise be hard to identify. Automatically generated colour schemes corresponding to various residue-based scores are provided, which allow the assessment of the conservation of backbone and side-chain conformations relative to the local coordinate frame. Structural comparison tools such as ProSMART can help to break the complexity that accompanies the constantly growing pool of structural data into a more readily accessible form, potentially offering biological insight or influencing subsequent experiments.

The crystal structure of ferritin from Chlorobium tepidum reveals a new conformation of the 4-fold channel for this protein family

Ferritins are ubiquitous iron-storage proteins found in all kingdoms of life. They share a common architecture made of 24 subunits of five α-helices. The recombinant Chlorobium tepidum ferritin (rCtFtn) is a structurally interesting protein since sequence alignments with other ferritins show that this protein has a significantly extended C-terminus, which possesses 12 histidine residues as well as several aspartate and glutamic acid residues that are potential metal ion binding residues. We show that the macromolecular assembly of rCtFtn exhibits a cage-like hollow shell consisting of 24 monomers that are related by 4-3-2 symmetry; similar to the assembly of other ferritins. In all ferritins of known structure the short fifth α-helix adopts an acute angle with respect to the four-helix bundle. However, the crystal structure of the rCtFtn presented here shows that this helix adopts a new conformation defining a new assembly of the 4-fold channel of rCtFtn. This conformation allows the arrangement of the C-terminal region into the inner cavity of the protein shell. Furthermore, two Fe(III) ions were found in each ferroxidase center of rCtFtn, with an average FeA-FeB distance of 3 Å; corresponding to a diferric μ-oxo/hydroxo species. This is the first ferritin crystal structure with an isolated di-iron center in an iron-storage ferritin. The crystal structure of rCtFtn and the biochemical results presented here, suggests that rCtFtn presents similar biochemical properties reported for other members of this protein family albeit with distinct structural plasticity.

Parallel implementation of 3D protein structure similarity searches using a GPU and the CUDA

Searching for similar 3D protein structures is one of the primary processes employed in the field of structural bioinformatics. However, the computational complexity of this process means that it is constantly necessary to search for new methods that can perform such a process faster and more efficiently. Finding molecular substructures that complex protein structures have in common is still a challenging task, especially when entire databases containing tens or even hundreds of thousands of protein structures must be scanned. Graphics processing units (GPUs) and general purpose graphics processing units (GPGPUs) can perform many time-consuming and computationally demanding processes much more quickly than a classical CPU can. In this paper, we describe the GPU-based implementation of the CASSERT algorithm for 3D protein structure similarity searching. This algorithm is based on the two-phase alignment of protein structures when matching fragments of the compared proteins. The GPU (GeForce GTX 560Ti: 384 cores, 2GB RAM) implementation of CASSERT (“GPU-CASSERT”) parallelizes both alignment phases and yields an average 180-fold increase in speed over its CPU-based, single-core implementation on an Intel Xeon E5620 (2.40GHz, 4 cores). In this paper, we show that massive parallelization of the 3D structure similarity search process on many-core GPU devices can reduce the execution time of the process, allowing it to be performed in real time. GPU-CASSERT is available at: http://zti.polsl.pl/dmrozek/science/gpucassert/cassert.htm.

Nucleotide-induced asymmetry within ATPase activator ring drives σ54-RNAP interaction and ATP hydrolysis

It is largely unknown how the typical homomeric ring geometry of ATPases associated with various cellular activities enables them to perform mechanical work. Small-angle solution X-ray scattering, crystallography, and electron microscopy (EM) reconstructions revealed that partial ATP occupancy caused the heptameric closed ring of the bacterial enhancer-binding protein (bEBP) NtrC1 to rearrange into a hexameric split ring of striking asymmetry. The highly conserved and functionally crucial GAFTGA loops responsible for interacting with σ54-RNA polymerase formed a spiral staircase. We propose that splitting of the ensemble directs ATP hydrolysis within the oligomer, and the ring's asymmetry guides interaction between ATPase and the complex of σ54 and promoter DNA. Similarity between the structure of the transcriptional activator NtrC1 and those of distantly related helicases Rho and E1 reveals a general mechanism in homomeric ATPases whereby complex allostery within the ring geometry forms asymmetric functional states that allow these biological motors to exert directional forces on their target macromolecules.

A new family of proteins related to the HEAT-like repeat DNA glycosylases with affinity for branched DNA structures

The recently discovered HEAT-like repeat (HLR) DNA glycosylase superfamily is widely distributed in all domains of life. The present bioinformatics and phylogenetic analysis shows that HLR DNA glycosylase superfamily members in the genus Bacillus form three subfamilies: AlkC, AlkD and AlkF/AlkG. The crystal structure of AlkF shows structural similarity with the DNA glycosylases AlkC and AlkD, however neither AlkF nor AlkG display any DNA glycosylase activity. Instead, both proteins have affinity to branched DNA structures such as three-way and Holliday junctions. A unique β-hairpin in the AlkF/AlkG subfamily is most likely inserted into the DNA major groove, and could be a structural determinant regulating DNA substrate affinity. We conclude that AlkF and AlkG represent a new family of HLR proteins with affinity for branched DNA structures.

Utility of the bacteriophage RB69 polymerase gp43 as a surrogate enzyme for herpesvirus orthologs

Viral polymerases are important targets in drug discovery and development efforts. Most antiviral compounds that are currently approved for treatment of infection with members of the herpesviridae family were shown to inhibit the viral DNA polymerase. However, biochemical studies that shed light on mechanisms of drug action and resistance are hampered primarily due to technical problems associated with enzyme expression and purification. In contrast, the orthologous bacteriophage RB69 polymerase gp43 has been crystallized in various forms and therefore serves as a model system that provides a better understanding of structure–function relationships of polymerases that belong the type B family. This review aims to discuss strengths, limitations, and opportunities of the phage surrogate with emphasis placed on its utility in the discovery and development of anti-herpetic drugs.

Antigenic switching of hepatitis B virus by alternative dimerization of the capsid protein

Chronic hepatitis B virus (HBV) infection afflicts millions worldwide with cirrhosis and liver cancer. HBV e-antigen (HBeAg), a clinical marker for disease severity, is a nonparticulate variant of the protein (core antigen, HBcAg) that forms the building-blocks of capsids. HBeAg is not required for virion production, but is implicated in establishing immune tolerance and chronic infection. Here, we report the crystal structure of HBeAg, which clarifies how the short N-terminal propeptide of HBeAg induces a radically altered mode of dimerization relative to HBcAg (∼140° rotation), locked into place through formation of intramolecular disulfide bridges. This structural switch precludes capsid assembly and engenders a distinct antigenic repertoire, explaining why the two antigens are cross-reactive at the T cell level (through sequence identity) but not at the B cell level (through conformation). The structure offers insight into how HBeAg may establish immune tolerance for HBcAg while evading its robust immunogenicity.

Formatt: Correcting protein multiple structural alignments by incorporating sequence alignment

BACKGROUND

The quality of multiple protein structure alignments are usually computed and assessed based on geometric functions of the coordinates of the backbone atoms from the protein chains. These purely geometric methods do not utilize directly protein sequence similarity, and in fact, determining the proper way to incorporate sequence similarity measures into the construction and assessment of protein multiple structure alignments has proved surprisingly difficult.

RESULTS

We present Formatt, a multiple structure alignment based on the Matt purely geometric multiple structure alignment program, that also takes into account sequence similarity when constructing alignments. We show that Formatt outperforms Matt and other popular structure alignment programs on the popular HOMSTRAD benchmark. For the SABMark twilight zone benchmark set that captures more remote homology, Formatt and Matt outperform other programs; depending on choice of embedded sequence aligner, Formatt produces either better sequence and structural alignments with a smaller core size than Matt, or similarly sized alignments with better sequence similarity, for a small cost in average RMSD.

CONCLUSIONS

Considering sequence information as well as purely geometric information seems to improve quality of multiple structure alignments, though defining what constitutes the best alignment when sequence and structural measures would suggest different alignments remains a difficult open question.

Crystallization, dehydration and experimental phasing of WbdD, a bifunctional kinase and methyltransferase from Escherichia coli O9a

WbdD is a bifunctional kinase/methyltransferase that is responsible for regulation of lipopolysaccharide O antigen polysaccharide chain length in Escherichia coli serotype O9a. Solving the crystal structure of this protein proved to be a challenge because the available crystals belonging to space group I23 only diffracted to low resolution (>95% of the crystals diffracted to resolution lower than 4 Å and most only to 8 Å) and were non-isomorphous, with changes in unit-cell dimensions of greater than 10%. Data from a serendipitously found single native crystal that diffracted to 3.0 Å resolution were non-isomorphous with a lower (3.5 Å) resolution selenomethionine data set. Here, a strategy for improving poor (3.5 Å resolution) initial phases by density modification and cross-crystal averaging with an additional 4.2 Å resolution data set to build a crude model of WbdD is desribed. Using this crude model as a mask to cut out the 3.5 Å resolution electron density yielded a successful molecular-replacement solution of the 3.0 Å resolution data set. The resulting map was used to build a complete model of WbdD. The hydration status of individual crystals appears to underpin the variable diffraction quality of WbdD crystals. After the initial structure had been solved, methods to control the hydration status of WbdD were developed and it was thus possible to routinely obtain high-resolution diffraction (to better than 2.5 Å resolution). This novel and facile crystal-dehydration protocol may be useful for similar challenging situations.

Evolutionary inaccuracy of pairwise structural alignments

Motivation: Structural alignment methods are widely used to generate gold standard alignments for improving multiple sequence alignments and transferring functional annotations, as well as for assigning structural distances between proteins. However, the correctness of the alignments generated by these methods is difficult to assess objectively since little is known about the exact evolutionary history of most proteins. Since homology is an equivalence relation, an upper bound on alignment quality can be found by assessing the consistency of alignments. Measuring the consistency of current methods of structure alignment and determining the causes of inconsistencies can, therefore, provide information on the quality of current methods and suggest possibilities for further improvement.

Results: We analyze the self-consistency of seven widely-used structural alignment methods (SAP, TM-align, Fr-TM-align, MAMMOTH, DALI, CE and FATCAT) on a diverse, non-redundant set of 1863 domains from the SCOP database and demonstrate that even for relatively similar proteins the degree of inconsistency of the alignments on a residue level is high (30%). We further show that levels of consistency vary substantially between methods, with two methods (SAP and Fr-TM-align) producing more consistent alignments than the rest. Inconsistency is found to be higher near gaps and for proteins of low structural complexity, as well as for helices. The ability of the methods to identify good structural alignments is also assessed using geometric measures, for which FATCAT (flexible mode) is found to be the best performer despite being highly inconsistent. We conclude that there is substantial scope for improving the consistency of structural alignment methods.

Contact:msadows@nimr.mrc.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Structural insight into activation mechanism of Toxoplasma gondii nucleoside triphosphate diphosphohydrolases by disulfide reduction

The intracellular parasite Toxoplasma gondii produces two nucleoside triphosphate diphosphohydrolases (NTPDase1 and -3). These tetrameric, cysteine-rich enzymes require activation by reductive cleavage of a hitherto unknown disulfide bond. Despite a 97% sequence identity, both isozymes differ largely in their ability to hydrolyze ATP and ADP. Here, we present crystal structures of inactive NTPDase3 as an apo form and in complex with the product AMP to resolutions of 2.0 and 2.2 Å, respectively. We find that the enzyme is present in an open conformation that precludes productive substrate binding and catalysis. The cysteine bridge 258-268 is identified to be responsible for locking of activity. Crystal structures of constitutively active variants of NTPDase1 and -3 generated by mutation of Cys(258)-Cys(268) show that opening of the regulatory cysteine bridge induces a pronounced contraction of the whole tetramer. This is accompanied by a 12° domain closure motion resulting in the correct arrangement of all active site residues. A complex structure of activated NTPDase3 with a non-hydrolyzable ATP analog and the cofactor Mg(2+) to a resolution of 2.85 Å indicates that catalytic differences between the NTPDases are primarily dictated by differences in positioning of the adenine base caused by substitution of Arg(492) and Glu(493) in NTPDase1 by glycines in NTPDase3.

Computational exploration of structural hypotheses for an additional sequence in a mammalian mitochondrial protein

BACKGROUND

Proteins involved in mammalian mitochondrial translation, when compared to analogous bacterial proteins, frequently have additional sequence regions whose structural or functional roles are not always clear. For example, an additional short insert sequence in the bovine mitochondrial initiation factor 2 (IF2(mt)) seems sufficient to fulfill the added role of eubacterial initiation factor IF1. Prior to our recent cryo-EM study that showed IF2(mt) to structurally occupy both the IF1 and IF2 binding sites, the spatial separation of these sites, and the short length of the insert sequence, posed ambiguity in whether it could perform the role of IF1 through occupation of the IF1 binding site on the ribosome.

RESULTS

The present study probes how well computational structure prediction methods can a priori address hypothesized roles of such additional sequences by creating quasi-atomic models of IF2(mt) using bacterial IF2 cryo-EM densities (that lack the insert sequences). How such initial IF2(mt) predictions differ from the observed IF2(mt) cryo-EM map and how they can be suitably improved using further sequence analysis and flexible fitting are analyzed.

CONCLUSIONS

By hypothesizing that the insert sequence occupies the IF1 binding site, continuous IF2(mt) models that occupy both the IF2 and IF1 binding sites can be predicted computationally. These models can be improved by flexible fitting into the IF2(mt) cryo-EM map to get reasonable quasi-atomic IF2(mt) models, but the exact orientation of the insert structure may not be reproduced. Specific eukaryotic insert sequence conservation characteristics can be used to predict alternate IF2(mt) models that have minor secondary structure rearrangements but fewer unusually extended linker regions. Computational structure prediction methods can thus be combined with medium-resolution cryo-EM maps to explore structure-function hypotheses for additional sequence regions and to guide further biochemical experiments, especially in mammalian systems where high-resolution structures are difficult to determine.

Crystal structure of plasminogen activator inhibitor-1 in an active conformation with normal thermodynamic stability

The serpin plasminogen activator inhibitor-1 (PAI-1) is a crucial regulator in fibrinolysis and tissue remodeling. PAI-1 has been associated with several pathological conditions and is a validated prognostic marker in human cancers. However, structural information about the native inhibitory form of PAI-1 has been elusive because of its inherent conformational instability and rapid conversion to a latent, inactive structure. Here we report the crystal structure of PAI-1 W175F at 2.3 Å resolution as the first model of the metastable native molecule. Structural comparison with a quadruple mutant (14-1B) previously used as representative of the active state uncovered key differences. The most striking differences occur near the region that houses three of the four mutations in the 14-1B PAI-1 structure. Prominent changes are localized within a loop connecting β-strand 3A with the F helix, in which a previously observed 3(10)-helix is absent in the new structure. Notably these structural changes are found near the binding site for the cofactor vitronectin. Because vitronectin is the only known physiological regulator of PAI-1 that slows down the latency conversion, the structure of this region is important. Furthermore, the previously identified chloride-binding site close to the F-helix is absent from the present structure and likely to be artifactual, because of its dependence on the 14-1B mutations. Instead we found a different chlorine-binding site that is likely to be present in wild type PAI-1 and that more satisfactorily accounts for the chlorine stabilizing effect on PAI-1.

Engagement of arginine finger to ATP triggers large conformational changes in NtrC1 AAA+ ATPase for remodeling bacterial RNA polymerase

The NtrC-like AAA+ ATPases control virulence and other important bacterial activities through delivering mechanical work to σ54-RNA polymerase to activate transcription from σ54-dependent genes. We report the first crystal structure for such an ATPase, NtrC1 of Aquifex aeolicus, in which the catalytic arginine engages the γ-phosphate of ATP. Comparing the new structure with those previously known for apo and ADP-bound states supports a rigid-body displacement model that is consistent with large-scale conformational changes observed by low-resolution methods. First, the arginine finger induces rigid-body roll, extending surface loops above the plane of the ATPase ring to bind σ54. Second, ATP hydrolysis permits Pi release and retraction of the arginine with a reversed roll, remodeling σ54-RNAP. This model provides a fresh perspective on how ATPase subunits interact within the ring-ensemble to promote transcription, directing attention to structural changes on the arginine-finger side of an ATP-bound interface.

Insertion domain within mammalian mitochondrial translation initiation factor 2 serves the role of eubacterial initiation factor 1

Mitochondria have their own translational machineries for the synthesis of thirteen polypeptide chains that are components of the complexes that participate in the process of oxidative phosphorylation (or ATP generation). Translation initiation in mammalian mitochondria requires two initiation factors, IF2(mt) and IF3(mt), instead of the three that are present in eubacteria. The mammalian IF2(mt) possesses a unique 37 amino acid insertion domain, which is known to be important for the formation of the translation initiation complex. We have obtained a three-dimensional cryoelectron microscopic map of the mammalian IF2(mt) in complex with initiator fMet-tRNA(iMet) and the eubacterial ribosome. We find that the 37 amino acid insertion domain interacts with the same binding site on the ribosome that would be occupied by the eubacterial initiation factor IF1, which is absent in mitochondria. Our finding suggests that the insertion domain of IF2(mt) mimics the function of eubacterial IF1, by blocking the ribosomal aminoacyl-tRNA binding site (A site) at the initiation step.

A mathematical framework for protein structure comparison

Comparison of protein structures is important for revealing the evolutionary relationship among proteins, predicting protein functions and predicting protein structures. Many methods have been developed in the past to align two or multiple protein structures. Despite the importance of this problem, rigorous mathematical or statistical frameworks have seldom been pursued for general protein structure comparison. One notable issue in this field is that with many different distances used to measure the similarity between protein structures, none of them are proper distances when protein structures of different sequences are compared. Statistical approaches based on those non-proper distances or similarity scores as random variables are thus not mathematically rigorous. In this work, we develop a mathematical framework for protein structure comparison by treating protein structures as three-dimensional curves. Using an elastic Riemannian metric on spaces of curves, geodesic distance, a proper distance on spaces of curves, can be computed for any two protein structures. In this framework, protein structures can be treated as random variables on the shape manifold, and means and covariance can be computed for populations of protein structures. Furthermore, these moments can be used to build Gaussian-type probability distributions of protein structures for use in hypothesis testing. The covariance of a population of protein structures can reveal the population-specific variations and be helpful in improving structure classification. With curves representing protein structures, the matching is performed using elastic shape analysis of curves, which can effectively model conformational changes and insertions/deletions. We show that our method performs comparably with commonly used methods in protein structure classification on a large manually annotated data set.

Protein structure comparison is important for understanding the evolutionary relationships among proteins, predicting protein functions, and predicting protein structures. Despite its importance, there have been no rigorous mathematical or statistical frameworks for protein structure comparison. One notable issue in this field is that with many different similarity measures used in comparing protein structures, none of them are proper distances when protein structures of different sequences are compared. In this study, we develop a mathematical framework for protein structure comparison by treating protein structures as three dimensional curves. A formal distance, geodesic distance, can be computed for any two protein structures. In this framework, population-specific variations within protein families can be characterized through building probability distributions for structures of protein families. The mean and covariance computed from groups of protein structures can also help to improve the classifications of protein structures. With curves representing protein structures, the matching is performed using elastic shape analysis of curves, which can effectively model conformational changes and insertions/deletions.

X-ray-induced deterioration of disulfide bridges at atomic resolution

Overall and site-specific X-ray-induced damage to porcine pancreatic elastase was studied at atomic resolution at temperatures of 100 and 15 K. The experiments confirmed that irradiation causes small movements of protein domains and bound water molecules in protein crystals. These structural changes occur not only at 100 K but also at temperatures as low as 15 K. An investigation of the deterioration of disulfide bridges demonstrated the following. (i) A decrease in the occupancy of S(γ) atoms and the appearance of new cysteine rotamers occur simultaneously. (ii) The occupancy decrease is observed for all S(γ) atoms, while new rotamers arise for some of the cysteine residues; the appearance of new conformations correlates with the accessibility to solvent. (iii) The sum of the occupancies of the initial and new conformations of a cysteine residue is approximately equal to the occupancy of the second cysteine residue in the bridge. (iv) The most pronounced changes occur at doses below 1.4 × 10(7) Gy, with only small changes occurring at higher doses. Comparison of the radiation-induced changes in an elastase crystal at 100 and 15 K suggested that the dose needed to induce a similar level of deterioration of the disulfide bonds and atomic displacements at 15 K to those seen at 100 K is more than two times higher.

Torsional network model: normal modes in torsion angle space better correlate with conformation changes in proteins

We introduce the torsional network model (TNM), an elastic network model whose degrees of freedom are the torsion angles of the protein backbone. Normal modes of the TNM displace backbone atoms including C(β) maintaining their covalent geometry. For many proteins, low frequency TNM modes are localized in torsion space yet collective in Cartesian space, reminiscent of hinge motions. A smaller number of TNM modes than anisotropic network model modes are enough to represent experimentally observed conformation changes. We observed significant correlation between the contribution of each normal mode to equilibrium fluctuations and to conformation changes, and defined the excess correlation with respect to a simple neutral model. The stronger this excess correlation, the lower the predicted free energy barrier of the conformation change and the fewer modes contribute to the change.

deconSTRUCT: general purpose protein database search on the substructure level

deconSTRUCT webserver offers an interface to a protein database search engine, usable for a general purpose detection of similar protein (sub)structures. Initially, it deconstructs the query structure into its secondary structure elements (SSEs) and reassembles the match to the target by requiring a (tunable) degree of similarity in the direction and sequential order of SSEs. Hierarchical organization and judicious use of the information about protein structure enables deconSTRUCT to achieve the sensitivity and specificity of the established search engines at orders of magnitude increased speed, without tying up irretrievably the substructure information in the form of a hash. In a post-processing step, a match on the level of the backbone atoms is constructed. The results presented to the user consist of the list of the matched SSEs, the transformation matrix for rigid superposition of the structures and several ways of visualization, both downloadable and implemented as a web-browser plug-in. The server is available at http://epsf.bmad.bii.a-star.edu.sg/struct_server.html.

Accuracy analysis of multiple structure alignments

Protein structure alignment methods are essential for many different challenges in protein science, such as the determination of relations between proteins in the fold space or the analysis and prediction of their biological function. A number of different pairwise and multiple structure alignment (MStA) programs have been developed and provided to the community. Prior knowledge of the expected alignment accuracy is desirable for the user of such tools. To retrieve an estimate of the performance of current structure alignment methods, we compiled a test suite taken from literature and the SISYPHUS database consisting of proteins that are difficult to align. Subsequently, different MStA programs were evaluated regarding alignment correctness and general limitations. The analysis shows that there are large differences in the success between the methods in terms of applicability and correctness. The latter ranges from 44 to 75% correct core positions. Taking only the best method result per test case this number increases to 84%. We conclude that the methods available are applicable to difficult cases, but also that there is still room for improvements in both, practicability and alignment correctness. An approach that combines the currently available methods supported by a proper score would be useful. Until then, a user should not rely on just a single program.

BlockMaster: partitioning protein kinase structures using normal-mode analysis

Protein kinases are key signaling enzymes which are dysregulated in many health disorders and therefore represent major targets of extensive drug discovery efforts. Their regulation in the cell is exerted via various mechanisms, including control of the 3D conformation of their catalytic domains. We developed a procedure, BlockMaster, for partitioning protein structures into semirigid blocks and flexible regions based on residue-residue correlations calculated from normal modes. BlockMaster provided correct partitioning into domains and subdomains of several test set proteins for which documented expert annotation of subdomains exists. When applied to representative structures of protein kinases, BlockMaster identified semirigid blocks within the traditional N-terminal and C-terminal lobes of the kinase domain. In general, the block regions had elevated helical content and reduced, but significant, coil content compared to the nonblock (flexible) regions. The specificity-determining regions, previously used to derive inhibitory peptides, were found to be more flexible in the tyrosine kinases than in serine/threonine kinases. Two blocks were identified which spanned both lobes. The first, which we termed the "pivot" block, included the alphaC-beta4 loop in the N-terminal lobe and part of the activation loop in the C-terminal lobe and appeared in both the active and inactive conformations of the kinases. The second, which we termed the "loop" block, differed between the active and inactive conformations. In the structures of active kinases, this block included part of the activation loop in the C-terminal lobe and the alphaC helix in the N-terminal lobe, representing a known interaction that stabilizes the active conformation. In the inactive structures, this block included G loop residues instead of the alphaC residues. This novel inactive "loop" block may stabilize the inactive conformation and thus downregulate kinase activity.

Pushing structural information into the yeast interactome by high-throughput protein docking experiments

The last several years have seen the consolidation of high-throughput proteomics initiatives to identify and characterize protein interactions and macromolecular complexes in model organisms. In particular, more that 10,000 high-confidence protein-protein interactions have been described between the roughly 6,000 proteins encoded in the budding yeast genome (Saccharomyces cerevisiae). However, unfortunately, high-resolution three-dimensional structures are only available for less than one hundred of these interacting pairs. Here, we expand this structural information on yeast protein interactions by running the first-ever high-throughput docking experiment with some of the best state-of-the-art methodologies, according to our benchmarks. To increase the coverage of the interaction space, we also explore the possibility of using homology models of varying quality in the docking experiments, instead of experimental structures, and assess how it would affect the global performance of the methods. In total, we have applied the docking procedure to 217 experimental structures and 1,023 homology models, providing putative structural models for over 3,000 protein-protein interactions in the yeast interactome. Finally, we analyze in detail the structural models obtained for the interaction between SAM1-anthranilate synthase complex and the MET30-RNA polymerase III to illustrate how our predictions can be straightforwardly used by the scientific community. The results of our experiment will be integrated into the general 3D-Repertoire pipeline, a European initiative to solve the structures of as many as possible protein complexes in yeast at the best possible resolution. All docking results are available at http://gatealoy.pcb.ub.es/HT_docking/.

Proteins are the main perpetrators of most biological processes. However, they seldom act alone, and most cellular functions are, in fact, carried out by large macromolecular complexes and regulated through intricate protein-protein interaction networks. Consequently, large efforts have been devoted to unveil protein interrelationships in a high-throughput manner, and the last several years have seen the consecution of the first interactome drafts for several model organisms. Unfortunately, these studies only reveal whether two proteins interact, but not the molecular bases of these interactions. A full comprehension of how proteins bind and form complexes can only come from high-resolution, three-dimensional (3D) structures, since they provide the key quasi-atomic details necessary to understand how the individual components in a complex or pathway are assembled and coordinated to function as a molecular unit. Here, we use protein docking experiments, in a high-throughput manner, to predict the 3D structure of over 3,000 interactions in yeast, which will be used to complement the complex structures obtained within the 3D-Repertoire pan-European initiative (http://www.3drepertoire.org).

The SALAMI protein structure search server

Protein structures often show similarities to another which would not be seen at the sequence level. Given the coordinates of a protein chain, the SALAMI server atwww.zbh.uni-hamburg.de/salami will search the protein data bank and return a set of similar structures without using sequence information. The results page lists the related proteins, details of the sequence and structure similarity and implied sequence alignments. Via a simple structure viewer, one can view superpositions of query and library structures and finally download superimposed coordinates. The alignment method is very tolerant of large gaps and insertions, and tends to produce slightly longer alignments than other similar programs.

Cross-over between discrete and continuous protein structure space: insights into automatic classification and networks of protein structures

Structural classifications of proteins assume the existence of the fold, which is an intrinsic equivalence class of protein domains. Here, we test in which conditions such an equivalence class is compatible with objective similarity measures. We base our analysis on the transitive property of the equivalence relationship, requiring that similarity of A with B and B with C implies that A and C are also similar. Divergent gene evolution leads us to expect that the transitive property should approximately hold. However, if protein domains are a combination of recurrent short polypeptide fragments, as proposed by several authors, then similarity of partial fragments may violate the transitive property, favouring the continuous view of the protein structure space. We propose a measure to quantify the violations of the transitive property when a clustering algorithm joins elements into clusters, and we find out that such violations present a well defined and detectable cross-over point, from an approximately transitive regime at high structure similarity to a regime with large transitivity violations and large differences in length at low similarity. We argue that protein structure space is discrete and hierarchic classification is justified up to this cross-over point, whereas at lower similarities the structure space is continuous and it should be represented as a network. We have tested the qualitative behaviour of this measure, varying all the choices involved in the automatic classification procedure, i.e., domain decomposition, alignment algorithm, similarity score, and clustering algorithm, and we have found out that this behaviour is quite robust. The final classification depends on the chosen algorithms. We used the values of the clustering coefficient and the transitivity violations to select the optimal choices among those that we tested. Interestingly, this criterion also favours the agreement between automatic and expert classifications. As a domain set, we have selected a consensus set of 2,890 domains decomposed very similarly in SCOP and CATH. As an alignment algorithm, we used a global version of MAMMOTH developed in our group, which is both rapid and accurate. As a similarity measure, we used the size-normalized contact overlap, and as a clustering algorithm, we used average linkage. The resulting automatic classification at the cross-over point was more consistent than expert ones with respect to the structure similarity measure, with 86% of the clusters corresponding to subsets of either SCOP or CATH superfamilies and fewer than 5% containing domains in distinct folds according to both SCOP and CATH. Almost 15% of SCOP superfamilies and 10% of CATH superfamilies were split, consistent with the notion of fold change in protein evolution. These results were qualitatively robust for all choices that we tested, although we did not try to use alignment algorithms developed by other groups. Folds defined in SCOP and CATH would be completely joined in the regime of large transitivity violations where clustering is more arbitrary. Consistently, the agreement between SCOP and CATH at fold level was lower than their agreement with the automatic classification obtained using as a clustering algorithm, respectively, average linkage (for SCOP) or single linkage (for CATH). The networks representing significant evolutionary and structural relationships between clusters beyond the cross-over point may allow us to perform evolutionary, structural, or functional analyses beyond the limits of classification schemes. These networks and the underlying clusters are available at http://ub.cbm.uam.es/research/ProtNet.php

Making order of the fast-growing information on proteins is essential for gaining evolutionary and functional knowledge. The most successful approaches to this task are based on classifications of protein structures, such as SCOP and CATH, which assume a discrete view of the protein structure space as a collection of separated equivalence classes (folds). However, several authors proposed that protein domains should be regarded as assemblies of polypeptide fragments, which implies that the protein–structure space is continuous. Here, we assess these views of domain space through the concept of transitivity; i.e., we test whether structure similarity of A with B and B with C implies that A and C are similar, as required for consistent classification. We find that the domain space is approximately transitive and discrete at high similarity and continuous at low similarity, where transitivity is severely violated. Comparing our classification at the cross-over similarity with CATH and SCOP, we find that they join proteins at low similarity where classification is inconsistent. Part of this discrepancy is due to structural divergence of homologous domains, which are forced to be in a single cluster in CATH and SCOP. Structural and evolutionary relationships between consistent clusters are represented as a network in our approach, going beyond current protein classification schemes. We conjecture that our results are related to a change of evolutionary regime, from uniparental divergent evolution for highly related domains to assembly of large fragments for which the classical tree representation is unsuitable.

X-ray-radiation-induced cooperative atomic movements in protein

X-rays interact with biological matter and cause damage. Proteins and other macromolecules are damaged primarily by ionizing X-ray photons and secondarily by reactive radiolytic chemical species. In particular, protein molecules are damaged during X-ray diffraction experiments with protein crystals, which is, in many cases, a serious hindrance to structure solution. The local X-ray-induced structural changes of the protein molecule have been studied using a number of model systems. However, it is still not well understood whether these local chemical changes lead to global structural changes in protein and what the mechanism is. We present experimental evidence at atomic resolution indicating the movement of large parts of the protein globule together with bound water molecules in the early stages of radiation damage to the protein crystal. The data were obtained from a crystal cryocooled to approximately 100 K and diffracting to 1 A. The movement of the protein structural elements occurs simultaneously with the decarboxylation of several glutamate and aspartate residues that mediate contacts between moving protein structural elements and with the rearrangement of the water network. The analysis of the anisotropy of atomic displacement parameters reveals that the observed atomic movements occur at different rates in different unit cells of the crystal. Thus, the examination of the cooperative atomic movement enables us to better understand how radiation-induced local chemical and structural changes of the protein molecule eventually lead to disorder in protein crystals.