A comprehensive and non-redundant database of protein domain movements

Comparative transcriptome profile reveals insight into the antibacterial immunity mechanism of the loach (Misgurnus anguillicaudatus) fed with soybean fermented broth during lipopolysaccharide (LPS) exposure

Loach (Misgurnus anguillicaudatus) is a common freshwater commercial fish species in China. The meat of this fish is a good source of protein and other nutrients that are needed for human health. Aquaculture challenges such as diseases and pest susceptibility, excessive density, and nutritional deficiency result in low production of loach rather than increased demand. Due to a lack of knowledge about the immune system of loaches, we carried out this study to better understand its antibacterial molecular mechanism. Here, we performed RNA sequencing from liver tissue obtained from soya bean-fermented fed loach after subjecting it to the LPS challenge. The results revealed a total of 18,399 differentially expressed genes (DEGs) in the LPS-treated and control groups. There were 7482 DEGs that were upregulated and 10,917 DEGs were downregulated. The enrichment analysis of DEGs revealed that the majority of DEGs were found to be abundant in the pathways of DNA replication, spliceosome, nucleotide exception repair, cell cycle, and Herpes simplex virus 1 infection. Furthermore, qRT-PCR analysis of 21 selected DEGs demonstrated that the transcriptomic data is extremely reliable. Overall, this study provides insight into the molecular features and control mechanisms of genes that affect loach growth. The availability of this information will also contribute to the enhancement of the breeding and protection of loach resources.

The chromosome-level genome assembly of Fraxinus americana provides insights into the evolution of Oleaceae plants

White ash (Fraxinus americana linn.) originates from the southeastern United States. It is a tall and fast-growing tree species with strong salt-alkali resistance and cold tolerance, making it an important reforestation species and widely planted worldwide. Here, we completed the chromosome-level reference genome assembly of F. americana based on Illumina, PacBio, and Hi-C reads, with a genome size of 878.98 Mb, an N50 of 3.27 Mb, and a heterozygosity rate of 0.3 %. Based on de novo prediction, transcriptome prediction, and homology-based protein prediction, we obtained 39,538 genes. Approximately 843.21 Mb of the assembly genome was composed of 37,928 annotated protein-coding genes, with a gene function annotation rate of 95.93 %. 99.94 % of the overlap clusters (877.44 Mb) were anchored to 23 chromosomes. Synteny analysis of F. americana and other Oleaceae plants showed that F. americana underwent frequent chromosome rearrangements. The amplification of the Ale transposons effectively promoted the genome size of F. americana. Compared with other Oleaceae plants, the Glutathione S-transferase (GST) gene family in the F. americana genome has undergone significant expansion, which may help F. americana cope with adverse natural environments. Furthermore, we found that key enzyme-coding gene families related to lignin biosynthesis were expanded and highly expressed in F. americana leaves. These key genes drive lignin synthesis and benefit F. americana in fast-growing, as well as resisting biotic and abiotic stress. Overall, the F. americana genome assembly provides insights into the evolution of Oleaceae plants and provides abundant resources for breeding and germplasm conservation of white ash.

Kinetic and structural characterization of carboxyspermidine dehydrogenase of polyamine biosynthesis

Polyamines are positively charged alkylamines ubiquitous among eukaryotes, prokaryotes, and archaea. Humans obtain polyamines through dietary intake, metabolic production, or uptake of polyamines made by gut microbes. The polyamine biosynthetic pathway used by most gut microbes differs from that used by human cells. This alternative pathway employs carboxyspermidine dehydrogenase (CASDH), an enzyme with limited characterization. Here, we solved a 1.94 Å X-ray crystal structure of Bacteroides fragilis CASDH by molecular replacement. BfCASDH is composed of three domains with a fold similar to saccharopine dehydrogenase but with a distinct active site arrangement. Using steady-state methods, we determined kcat and kcat/Km for BfCASDH and Clostridium leptum CASDH using putrescine, diaminopropane, aspartate semi-aldehyde, NADH, and NADPH as substrates. These data revealed evidence of cooperativity in BfCASDH. Putrescine is the likely polyamine substrate and NADPH is the coenzyme used to complete the reaction, forming carboxyspermidine as a product. These data provide the first kinetic characterization of CASDH-a key enzyme in the production of microbial polyamines.

Structural basis for recognition and methylation of p97 by METTL21D, a valosin-containing protein lysine methyltransferase

p97 is a human AAA+ (ATPase associated with diverse cellular activities, also known as valosin-containing protein [VCP]) ATPase, which is involved in diverse cellular processes such as membrane fusion and proteolysis. Lysine-specific methyltransferase of p97 (METTL21D) was identified as a class I methyltransferase that catalyzes the trimethylation of Lys315 of p97, a so-called VCP lysine methyltransferase (VCPKMT). Interestingly, VCPKMT disassembles a single hexamer ring consisting of p97-D1 domain and methylates Lys315 residue. Herein, the structures of S-adenosyl-L-methionine-bound VCPKMT and S-adenosyl-L-homocysteine-bound VCPKMT in complex with p97 N/D1 (N21-Q458) were reported at a resolution of 1.8 Å and 2.8 Å, respectively. The structures revealed the molecular details for the recognition and methylation of monomeric p97 by VCPKMT. Using biochemical analysis, we also investigated whether the methylation of full-length p97 could be sufficiently enhanced through cooperation between VCPKMT and the C terminus of alveolar soft part sarcoma locus (ASPL). Our study provides the groundwork for future structural and mechanistic studies of p97 and inhibitors.

Evolutionary conservation of protein dynamics: insights from all-atom molecular dynamics simulations of 'peptidase' domain of Spt16

Protein function is encoded in its sequence, manifested in its three-dimensional structure, and facilitated by its dynamics. Studies have suggested that protein structures with higher sequence similarity could have more similar patterns of dynamics. However, such studies of protein dynamics within and across protein families typically rely on coarse-grained models, or approximate metrics like crystallographic B-factors. This study uses µs scale molecular dynamics (MD) simulations to explore the conservation of dynamics among homologs of ∼50 kDa N-terminal module of Spt16 (Spt16N). Spt16N from Saccharomyces cerevisiae (Sc-Spt16N) and three of its homologs with 30-40% sequence identities were available in the PDB. To make our data-set more comprehensive, the crystal structure of an additional homolog (62% sequence identity with Sc-Spt16N) was solved at 1.7 Å resolution. Cumulative MD simulations of 6 µs were carried out on these Spt16N structures and on two additional protein structures with varying degrees of similarity to it. The simulations revealed that correlation in patterns of backbone fluctuations vary linearly with sequence identity. This trend could not be inferred using crystallographic B-factors. Further, normal mode analysis suggested a similar pattern of inter-domain (inter-lobe) motions not only among Spt16N homologs, but also in the M24 peptidase structure. On the other hand, MD simulation results highlighted conserved motions that were found unique for Spt16N protein, this along with electrostatics trends shed light on functional aspects of Spt16N.Communicated by Ramaswamy H. Sarma.

Tertiary and Quaternary Structure Organization in GMP Synthetases: Implications for Catalysis

Glutamine amidotransferases, enzymes that transfer nitrogen from Gln to various cellular metabolites, are modular, with the amidotransferase (GATase) domain hydrolyzing Gln, generating ammonia and the acceptor domain catalyzing the addition of nitrogen onto its cognate substrate. GMP synthetase (GMPS), an enzyme in the de novo purine nucleotide biosynthetic pathway, is a glutamine amidotransferase that catalyzes the synthesis of GMP from XMP. The reaction involves activation of XMP though adenylation by ATP in the ATP pyrophosphatase (ATPPase) active site, followed by channeling and attack of NH₃ generated in the GATase pocket. This complex chemistry entails co-ordination of activity across the active sites, allosteric activation of the GATase domain to modulate Gln hydrolysis and channeling of ammonia from the GATase to the acceptor active site. Functional GMPS dimers associate through the dimerization domain. The crystal structure of the Gln-bound complex of Plasmodium falciparum GMPS (PfGMPS) for the first time revealed large-scale domain rotation to be associated with catalysis and leading to the juxtaposition of two otherwise spatially distal cysteinyl (C113/C337) residues. In this manuscript, we report on an unusual structural variation in the crystal structure of the C89A/C113A PfGMPS double mutant, wherein a larger degree of domain rotation has led to the dissociation of the dimeric structure. Furthermore, we report a hitherto overlooked signature motif tightly related to catalysis.

Structural insights into the mechanism of archaellar rotational switching

Signal transduction via phosphorylated CheY towards the flagellum and the archaellum involves a conserved mechanism of CheY phosphorylation and subsequent conformational changes within CheY. This mechanism is conserved among bacteria and archaea, despite substantial differences in the composition and architecture of archaellum and flagellum, respectively. Phosphorylated CheY has higher affinity towards the bacterial C-ring and its binding leads to conformational changes in the flagellar motor and subsequent rotational switching of the flagellum. In archaea, the adaptor protein CheF resides at the cytoplasmic face of the archaeal C-ring formed by the proteins ArlCDE and interacts with phosphorylated CheY. While the mechanism of CheY binding to the C-ring is well-studied in bacteria, the role of CheF in archaea remains enigmatic and mechanistic insights are absent. Here, we have determined the atomic structures of CheF alone and in complex with activated CheY by X-ray crystallography. CheF forms an elongated dimer with a twisted architecture. We show that CheY binds to the C-terminal tail domain of CheF leading to slight conformational changes within CheF. Our structural, biochemical and genetic analyses reveal the mechanistic basis for CheY binding to CheF and allow us to propose a model for rotational switching of the archaellum.

Signal transduction via phosphorylated CheY is conserved in bacteria and archaea. In this study, the authors employ structural biochemistry combined with cell biology to delineate the mechanism of CheY recognition by the adaptor protein CheF.

Enhanced Conformational Sampling with an Adaptive Coarse-Grained Elastic Network Model Using Short-Time All-Atom Molecular Dynamics

Compared to all-atom molecular dynamics (AA-MD) simulations, coarse-grained (CG) MD simulations can significantly reduce calculation costs. However, existing CG-MD methods are unsuitable for sampling structures that depart significantly from the initial structure without any biased force. In this study, we developed a new adaptive CG elastic network model (ENM), in which the dynamic cross-correlation coefficient based on short-time AA-MD of at most ns order is considered. By applying Bayesian optimization to search for a suitable parameter among the vast parameter space of adaptive CG-ENM, we succeeded in reducing the searching cost to approximately 10% of those for random sampling and exhaustive sampling. To evaluate the performance of adaptive CG-ENM, we applied the new methodology to adenylate kinase (ADK) and glutamine binding protein (GBP) in the apo state. The results showed that the structural ensembles explored by adaptive CG-ENM could be considerably more diverse than those by conventional ENMs with enhanced sampling such as temperature replica exchange MD and long-time AA-MD of 1 μs. In particular, some of the structures sampled by adaptive ENM are relatively close to the holo-type structures of ADK and GBP. Furthermore, as a challenging task, to demonstrate the advantages of the CG model with lower calculation cost, we applied our new methodology to a larger biomolecule, integrin (αV) in the inactive state. Then, we sampled various structural ensembles, including extended structures that are apparently different from inactive ones.

D3PM: a comprehensive database for protein motions ranging from residue to domain

Background

Knowledge of protein motions is significant to understand its functions. While currently available databases for protein motions are mostly focused on overall domain motions, little attention is paid on local residue motions. Albeit with relatively small scale, the local residue motions, especially those residues in binding pockets, may play crucial roles in protein functioning and ligands binding.

Results

A comprehensive protein motion database, namely D3PM, was constructed in this study to facilitate the analysis of protein motions. The protein motions in the D3PM range from overall structural changes of macromolecule to local flip motions of binding pocket residues. Currently, the D3PM has collected 7679 proteins with overall motions and 3513 proteins with pocket residue motions. The motion patterns are classified into 4 types of overall structural changes and 5 types of pocket residue motions. Impressively, we found that less than 15% of protein pairs have obvious overall conformational adaptations induced by ligand binding, while more than 50% of protein pairs have significant structural changes in ligand binding sites, indicating that ligand-induced conformational changes are drastic and mainly confined around ligand binding sites. Based on the residue preference in binding pocket, we classified amino acids into “pocketphilic” and “pocketphobic” residues, which should be helpful for pocket prediction and drug design.

Conclusion

D3PM is a comprehensive database about protein motions ranging from residue to domain, which should be useful for exploring diverse protein motions and for understanding protein function and drug design. The D3PM is available on www.d3pharma.com/D3PM/index.php.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04595-0.

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.

"Protein" no longer means what it used to

Every biologist knows that the word protein describes a group of macromolecules essential to sustain life on Earth. As biologists, we are invariably trained under a protein paradigm established since the early twentieth century. However, in recent years, the term protein unveiled itself as an euphemism to describe the overwhelming heterogeneity of these compounds. Most of our current studies are targeted on carefully selected subsets of proteins, but we tend to think and write about these as representative of the whole population. Here we discuss how seeking for universal definitions and general rules in any arbitrarily segmented study would be misleading about the conclusions. Of course, it is not our purpose to discourage the use of the word protein. Instead, we suggest to embrace the extended universe of proteins to reach a deeper understanding of their full potential, realizing that the term encompasses a group of molecules very heterogeneous in terms of size, shape, chemistry and functions, i.e. the term protein no longer means what it used to.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Order from disorder in the sarcomere: FATZ forms a fuzzy but tight complex and phase-separated condensates with α-actinin

In sarcomeres, α-actinin cross-links actin filaments and anchors them to the Z-disk. FATZ (filamin-, α-actinin-, and telethonin-binding protein of the Z-disk) proteins interact with α-actinin and other core Z-disk proteins, contributing to myofibril assembly and maintenance. Here, we report the first structure and its cellular validation of α-actinin-2 in complex with a Z-disk partner, FATZ-1, which is best described as a conformational ensemble. We show that FATZ-1 forms a tight fuzzy complex with α-actinin-2 and propose an interaction mechanism via main molecular recognition elements and secondary binding sites. The obtained integrative model reveals a polar architecture of the complex which, in combination with FATZ-1 multivalent scaffold function, might organize interaction partners and stabilize α-actinin-2 preferential orientation in Z-disk. Last, we uncover FATZ-1 ability to phase-separate and form biomolecular condensates with α-actinin-2, raising the question whether FATZ proteins can create an interaction hub for Z-disk proteins through membraneless compartmentalization during myofibrillogenesis.

Edge expansion parallel cascade selection molecular dynamics simulation for investigating large-amplitude collective motions of proteins

We propose edge expansion parallel cascade selection molecular dynamics (eePaCS-MD) as an efficient adaptive conformational sampling method to investigate the large-amplitude motions of proteins without prior knowledge of the conformational transitions. In this method, multiple independent MD simulations are iteratively conducted from initial structures randomly selected from the vertices of a multi-dimensional principal component subspace. This subspace is defined by an ensemble of protein conformations sampled during previous cycles of eePaCS-MD. The edges and vertices of the conformational subspace are determined by solving the "convex hull problem." The sampling efficiency of eePaCS-MD is achieved by intensively repeating MD simulations from the vertex structures, which increases the probability of rare event occurrence to explore new large-amplitude collective motions. The conformational sampling efficiency of eePaCS-MD was assessed by investigating the open-close transitions of glutamine binding protein, maltose/maltodextrin binding protein, and adenylate kinase and comparing the results to those obtained using related methods. In all cases, the open-close transitions were simulated in ∼10 ns of simulation time or less, offering 1-3 orders of magnitude shorter simulation time compared to conventional MD. Furthermore, we show that the combination of eePaCS-MD and accelerated MD can further enhance conformational sampling efficiency, which reduced the total computational cost of observing the open-close transitions by at most 36%.

Exploring Large Domain Motions in Proteins Using Atomistic Molecular Dynamics with Enhanced Conformational Sampling

Conformational transitions in multidomain proteins are essential for biological functions. The Apo conformations are typically open and flexible, while the Holo states form more compact conformations stabilized by protein-ligand interactions. Unfortunately, the atomically detailed mechanisms for such open-closed conformational changes are difficult to be accessed experimentally as well as computationally. To simulate the transitions using atomistic molecular dynamics (MD) simulations, efficient conformational sampling algorithms are required. In this work, we propose a new approach based on generalized replica-exchange with solute tempering (gREST) for exploring the open-closed conformational changes in multidomain proteins. Wherein, selected surface charged residues in a target protein are defined as the solute region in gREST simulation and the solute temperatures are different in replicas and exchanged between them to enhance the domain motions. This approach is called gREST selected surface charged residues (gREST_SSCR) and is applied to the Apo and Holo states of ribose binding protein (RBP) in solution. The conformational spaces sampled with gREST_SSCR are much wider than those with the conventional MD, sampling open-closed conformational changes while maintaining RBP domains’ stability. The free-energy landscapes of RBP in the Apo and Holo states are drawn along with twist and hinge angles of the two moving domains. The inter-domain salt-bridges that are not observed in the experimental structures are also important in the intermediate states during the conformational changes.

Comparative structural and mechanistic studies of 4-hydroxy-tetrahydrodipicolinate reductases from Mycobacterium tuberculosis and Vibrio vulnificus

BACKGROUND

The products of the lysine biosynthesis pathway, meso-diaminopimelate and lysine, are essential for bacterial survival. This paper focuses on the structural and mechanistic characterization of 4-hydroxy-tetrahydrodipicolinate reductase (DapB), which is one of the enzymes from the lysine biosynthesis pathway. DapB catalyzes the conversion of (2S, 4S)-4-hydroxy-2,3,4,5-tetrahydrodipicolinate (HTPA) to 2,3,4,5-tetrahydrodipicolinate in an NADH/NADPH dependent reaction. Genes coding for DapBs were identified as essential for many pathogenic bacteria, and therefore DapB is an interesting new target for the development of antibiotics.

METHODS

We have combined experimental and computational approaches to provide novel insights into mechanism of the DapB catalyzed reaction.

RESULTS

Structures of DapBs originating from Mycobacterium tuberculosis and Vibrio vulnificus in complexes with NAD⁺, NADP⁺, as well as with inhibitors, were determined and described. The structures determined by us, as well as currently available structures of DapBs from other bacterial species, were compared and used to elucidate a mechanism of reaction catalyzed by this group of enzymes. Several different computational methods were used to provide a detailed description of a plausible reaction mechanism.

CONCLUSIONS

This is the first report presenting the detailed mechanism of reaction catalyzed by DapB.

GENERAL SIGNIFICANCE

Structural data in combination with information on the reaction mechanism provide a background for development of DapB inhibitors, including transition-state analogues.

The roles of SDHAF2 and dicarboxylate in covalent flavinylation of SDHA, the human complex II flavoprotein

Mitochondrial complex II, also known as succinate dehydrogenase (SDH), is an integral-membrane heterotetramer (SDHABCD) that links two essential energy-producing processes, the tricarboxylic acid (TCA) cycle and oxidative phosphorylation. A significant amount of information is available on the structure and function of mature complex II from a range of organisms. However, there is a gap in our understanding of how the enzyme assembles into a functional complex, and disease-associated complex II insufficiency may result from incorrect function of the mature enzyme or from assembly defects. Here, we investigate the assembly of human complex II by combining a biochemical reconstructionist approach with structural studies. We report an X-ray structure of human SDHA and its dedicated assembly factor SDHAF2. Importantly, we also identify a small molecule dicarboxylate that acts as an essential cofactor in this process and works in synergy with SDHAF2 to properly orient the flavin and capping domains of SDHA. This reorganizes the active site, which is located at the interface of these domains, and adjusts the pK_a of SDHA^R451 so that covalent attachment of the flavin adenine dinucleotide (FAD) cofactor is supported. We analyze the impact of disease-associated SDHA mutations on assembly and identify four distinct conformational forms of the complex II flavoprotein that we assign to roles in assembly and catalysis.

Structure of P46, an immunodominant surface protein from Mycoplasma hyopneumoniae: interaction with a monoclonal antibody

Mycoplasma hyopneumoniae is a prokaryotic pathogen that colonizes the respiratory ciliated epithelial cells in swine. Infected animals suffer respiratory lesions, causing major economic losses in the porcine industry. Characterization of the immunodominant membrane-associated proteins from M. hyopneumoniae may be instrumental in the development of new therapeutic approaches. Here, the crystal structure of P46, one of the main surface-antigen proteins, from M. hyopneumoniae is presented and shows N- and C-terminal α/β domains connected by a hinge. The structures solved in this work include a ligand-free open form of P46 (3.1 Å resolution) and two ligand-bound structures of P46 with maltose (2.5 Å resolution) and xylose (3.5 Å resolution) in open and closed conformations, respectively. The ligand-binding site is buried in the cleft between the domains at the hinge region. The two domains of P46 can rotate with respect to each other, giving open or closed alternative conformations. In agreement with this structural information, sequence analyses show similarities to substrate-binding members of the ABC transporter superfamily, with P46 facing the extracellular side as a functional subunit. In the structure with xylose, P46 was also bound to a high-affinity (K_d = 29 nM) Fab fragment from a monoclonal antibody, allowing the characterization of a structural epitope in P46 that exclusively involves residues from the C-terminal domain. The Fab structure in the complex with P46 shows only small conformational rearrangements in the six complementarity-determining regions (CDRs) with respect to the unbound Fab (the structure of which is also determined in this work at 1.95 Å resolution). The structural information that is now available should contribute to a better understanding of sugar nutrient intake by M. hyopneumoniae. This information will also allow the design of protocols and strategies for the generation of new vaccines against this important swine pathogen.

Investigation of sequence features of hinge-bending regions in proteins with domain movements using kernel logistic regression

Background

Hinge-bending movements in proteins comprising two or more domains form a large class of functional movements. Hinge-bending regions demarcate protein domains and collectively control the domain movement. Consequently, the ability to recognise sequence features of hinge-bending regions and to be able to predict them from sequence alone would benefit various areas of protein research. For example, an understanding of how the sequence features of these regions relate to dynamic properties in multi-domain proteins would aid in the rational design of linkers in therapeutic fusion proteins.

Results

The DynDom database of protein domain movements comprises sequences annotated to indicate whether the amino acid residue is located within a hinge-bending region or within an intradomain region. Using statistical methods and Kernel Logistic Regression (KLR) models, this data was used to determine sequence features that favour or disfavour hinge-bending regions. This is a difficult classification problem as the number of negative cases (intradomain residues) is much larger than the number of positive cases (hinge residues). The statistical methods and the KLR models both show that cysteine has the lowest propensity for hinge-bending regions and proline has the highest, even though it is the most rigid amino acid. As hinge-bending regions have been previously shown to occur frequently at the terminal regions of the secondary structures, the propensity for proline at these regions is likely due to its tendency to break secondary structures. The KLR models also indicate that isoleucine may act as a domain-capping residue. We have found that a quadratic KLR model outperforms a linear KLR model and that improvement in performance occurs up to very long window lengths (eighty residues) indicating long-range correlations.

Conclusion

In contrast to the only other approach that focused solely on interdomain hinge-bending regions, the method provides a modest and statistically significant improvement over a random classifier. An explanation of the KLR results is that in the prediction of hinge-bending regions a long-range correlation is at play between a small number amino acids that either favour or disfavour hinge-bending regions. The resulting sequence-based prediction tool, HingeSeek, is available to run through a webserver at hingeseek.cmp.uea.ac.uk.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

Morphing and docking visualisation of biomolecular structures using Multi-Dimensional Scaling

Protein structures are often solved at atomic resolution in two states defining a functional movement but intervening conformations are usually unknown. Morphing methods generate intervening conformations between two known structures. When viewed as an animation using molecular graphics, a smooth, direct morph enables the eye to track changes in structure that might be otherwise missed. We present a morphing method that aims to linearly interpolate interatomic distances and which uses SMACOF (Scaling by MAjorisation of COmplicated Function) and multigrid techniques with a cut-off distance based weighting that optimizes the MolProbity score of intervening structures. The all-atom morphs are smooth, move directly between the two structures, and are shown, in general, to pass closer to a set of known intermediates than those generated using other methods. The techniques are also used for docking by putting the unbound structures in a "near-approach pose" and then morphing to the bound complex. The resulting GPU-accelerated tools are available on a webserver, Morphit_Pro, at http://morphit-pro.cmp.uea.ac.uk/ and more than 5000 domains movements available at the DynDom website can now be viewed as morphs http://morphit-pro.cmp.uea.ac.uk/dyndom/.

Structural changes of homodimers in the PDB

Protein complexes are involved in various biological phenomena. These complexes are intrinsically flexible, and structural changes are essential to their functions. To perform a large-scale automated analysis of the structural changes of complexes, we combined two original methods. An application, SCPC, compares two structures of protein complexes and decides the match of binding mode. Another application, Motion Tree, identifies rigid-body motions in various sizes and magnitude from the two structural complexes with the same binding mode. This approach was applied to all available homodimers in the Protein Data Bank (PDB). We defined two complex-specific motions: interface motion and subunit-spanning motion. In the former, each subunit of a complex constitutes a rigid body, and the relative movement between subunits occurs at the interface. In the latter, structural parts from distinct subunits constitute a rigid body, providing the relative movement spanning subunits. All structural changes were classified and examined. It was revealed that the complex-specific motions were common in the homodimers, detected in around 40% of families. The dimeric interfaces were likely to be small and flat for interface motion, while large and rugged for subunit-spanning motion. Interface motion was accompanied by a drastic change in contacts at the interface, while the change in the subunit-spanning motion was moderate. These results indicate that the interface properties of homodimers correlated with the type of complex-specific motion. The study demonstrates that the pipeline of SCPC and Motion Tree is useful for the massive analysis of structural change of protein complexes.

Molecular dynamics simulation of bacterial flagella

The bacterial flagellum is a biological nanomachine for the locomotion of bacteria, and is seen in organisms like Salmonella and Escherichia coli. The flagellum consists of tens of thousands of protein molecules and more than 30 different kinds of proteins. The basal body of the flagellum contains a protein export apparatus and a rotary motor that is powered by ion motive force across the cytoplasmic membrane. The filament functions as a propeller whose helicity is controlled by the direction of the torque. The hook that connects the motor and filament acts as a universal joint, transmitting torque generated by the motor to different directions. This report describes the use of molecular dynamics to study the bacterial flagellum. Molecular dynamics simulation is a powerful method that permits the investigation, at atomic resolution, of the molecular mechanisms of biomolecular systems containing many proteins and solvent. When applied to the flagellum, these studies successfully unveiled the polymorphic supercoiling and transportation mechanism of the filament, the universal joint mechanism of the hook, the ion transfer mechanism of the motor stator, the flexible nature of the transport apparatus proteins, and activation of proteins involved in chemotaxis.

Characterizing Protein Dynamics with Integrative Use of Bulk and Single-Molecule Techniques

A protein dynamically samples multiple conformations, and the conformational dynamics enables protein function. Most biophysical measurements are ensemble-based, with the observables averaged over all members of the ensemble. Though attainable, the decomposition of the observables to the constituent conformational states can be computationally expensive and ambiguous. Here we show that the incorporation of single-molecule fluorescence resonance energy transfer (smFRET) data resolves the ambiguity and affords protein ensemble structures that are more precise and accurate. Using K63-linked diubiquitin, we characterize the dynamic domain arrangements of the model system, with the use of chemical cross-linking coupled with mass spectrometry (CXMS), small-angle X-ray scattering (SAXS), and smFRET techniques. CXMS allows the modeling of protein conformational states that are alternatives to the crystal structure. SAXS provides ensemble-averaged low-resolution shape information. Importantly, smFRET affords state-specific populations, and the FRET distances validate the ensemble structures obtained by refining against CXMS and SAXS restraints. Together, the integrative use of bulk and single-molecule techniques affords better insight into protein dynamics and shall be widely implemented in structural biology.

Redox Switch for the Inhibited State of Yeast Glycogen Synthase Mimics Regulation by Phosphorylation

Glycogen synthase (GS) is the rate limiting enzyme in the synthesis of glycogen. Eukaryotic GS is negatively regulated by covalent phosphorylation and allosterically activated by glucose-6-phosphate (G-6-P). To gain structural insights into the inhibited state of the enzyme, we solved the crystal structure of yGsy2-R589A/R592A to a resolution of 3.3 Å. The double mutant has an activity ratio similar to the phosphorylated enzyme and also retains the ability to be activated by G-6-P. When compared to the 2.88 Å structure of the wild-type G-6-P activated enzyme, the crystal structure of the low-activity mutant showed that the N-terminal domain of the inhibited state is tightly held against the dimer-related interface thereby hindering acceptor access to the catalytic cleft. On the basis of these two structural observations, we developed a reversible redox regulatory feature in yeast GS by substituting cysteine residues for two highly conserved arginine residues. When oxidized, the cysteine mutant enzyme exhibits activity levels similar to the phosphorylated enzyme but cannot be activated by G-6-P. Upon reduction, the cysteine mutant enzyme regains normal activity levels and regulatory response to G-6-P activation.

ChSeq: A database of chameleon sequences

Chameleon sequences (ChSeqs) refer to sequence strings of identical amino acids that can adopt different conformations in protein structures. Researchers have detected and studied ChSeqs to understand the interplay between local and global interactions in protein structure formation. The different secondary structures adopted by one ChSeq challenge sequence-based secondary structure predictors. With increasing numbers of available Protein Data Bank structures, we here identify a large set of ChSeqs ranging from 6 to 10 residues in length. The homologous ChSeqs discovered highlight the structural plasticity involved in biological function. When compared with previous studies, the set of unrelated ChSeqs found represents an about 20-fold increase in the number of detected sequences, as well as an increase in the longest ChSeq length from 8 to 10 residues. We applied secondary structure predictors on our ChSeqs and found that methods based on a sequence profile outperformed methods based on a single sequence. For the unrelated ChSeqs, the evolutionary information provided by the sequence profile typically allows successful prediction of the prevailing secondary structure adopted in each protein family. Our dataset will facilitate future studies of ChSeqs, as well as interpretations of the interplay between local and nonlocal interactions. A user-friendly web interface for this ChSeq database is available at prodata.swmed.edu/chseq.

Quantitative method for the assignment of hinge and shear mechanism in protein domain movements

Motivation: A popular method for classification of protein domain movements apportions them into two main types: those with a ‘hinge’ mechanism and those with a ‘shear’ mechanism. The intuitive assignment of domain movements to these classes has limited the number of domain movements that can be classified in this way. Furthermore, whether intended or not, the term ‘shear’ is often interpreted to mean a relative translation of the domains.

Results: Numbers of occurrences of four different types of residue contact changes between domains were optimally combined by logistic regression using the training set of domain movements intuitively classified as hinge and shear to produce a predictor for hinge and shear. This predictor was applied to give a 10-fold increase in the number of examples over the number previously available with a high degree of precision. It is shown that overall a relative translation of domains is rare, and that there is no difference between hinge and shear mechanisms in this respect. However, the shear set contains significantly more examples of domains having a relative twisting movement than the hinge set. The angle of rotation is also shown to be a good discriminator between the two mechanisms.

Availability and implementation: Results are free to browse at http://www.cmp.uea.ac.uk/dyndom/interface/.

Contact:sjh@cmp.uea.ac.uk.

Supplementary information:Supplementary data are available at Bioinformatics online.

Structure of the trehalose-6-phosphate phosphatase from Brugia malayi reveals key design principles for anthelmintic drugs

Parasitic nematodes are responsible for devastating illnesses that plague many of the world's poorest populations indigenous to the tropical areas of developing nations. Among these diseases is lymphatic filariasis, a major cause of permanent and long-term disability. Proteins essential to nematodes that do not have mammalian counterparts represent targets for therapeutic inhibitor discovery. One promising target is trehalose-6-phosphate phosphatase (T6PP) from Brugia malayi. In the model nematode Caenorhabditis elegans, T6PP is essential for survival due to the toxic effect(s) of the accumulation of trehalose 6-phosphate. T6PP has also been shown to be essential in Mycobacterium tuberculosis. We determined the X-ray crystal structure of T6PP from B. malayi. The protein structure revealed a stabilizing N-terminal MIT-like domain and a catalytic C-terminal C2B-type HAD phosphatase fold. Structure-guided mutagenesis, combined with kinetic analyses using a designed competitive inhibitor, trehalose 6-sulfate, identified five residues important for binding and catalysis. This structure-function analysis along with computational mapping provided the basis for the proposed model of the T6PP-trehalose 6-phosphate complex. The model indicates a substrate-binding mode wherein shape complementarity and van der Waals interactions drive recognition. The mode of binding is in sharp contrast to the homolog sucrose-6-phosphate phosphatase where extensive hydrogen-bond interactions are made to the substrate. Together these results suggest that high-affinity inhibitors will be bi-dentate, taking advantage of substrate-like binding to the phosphoryl-binding pocket while simultaneously utilizing non-native binding to the trehalose pocket. The conservation of the key residues that enforce the shape of the substrate pocket in T6PP enzymes suggest that development of broad-range anthelmintic and antibacterial therapeutics employing this platform may be possible.

Here, we describe the structure of trehalose-6-phosphate phosphatase (T6PP) from Brugia malayi. This enzyme is essential to the organism; deletion of the gene encoding T6PP results in toxic accumulation of trehalose 6-phosphate. Structure-guided mutagenesis coupled with kinetic analyses revealed residues important for binding and catalysis. The model for substrate binding suggests a binding mode in which shape complementarity plays a major role. Conservation of binding residues among T6PP orthologs present in pathogenic nematodes and bacteria favors T6PP as a suitable target for broad-range anthelmintic and antibacterial drug design.

Structural basis for 2'-phosphate incorporation into glycogen by glycogen synthase

Glycogen is a glucose polymer that contains minor amounts of covalently attached phosphate. Hyperphosphorylation is deleterious to glycogen structure and can lead to Lafora disease. Recently, it was demonstrated that glycogen synthase catalyzes glucose-phosphate transfer in addition to its characteristic glucose transfer reaction. Glucose-1,2-cyclic-phosphate (GCP) was proposed to be formed from UDP-Glc breakdown and subsequently transferred, thus providing a source of phosphate found in glycogen. To gain further insight into the molecular basis for glucose-phosphate transfer, two structures of yeast glycogen synthase were determined; a 3.0-Å resolution structure of the complex with UMP/GCP and a 2.8-Å resolution structure of the complex with UDP/glucose. Structural superposition of the complexes revealed that the bound ligands and most active site residues are positioned similarly, consistent with the use of a common transfer mechanism for both reactions. The N-terminal domain of the UDP-glucose complex was found to be 13.3° more closed compared with a UDP complex. However, the UMP · GCP complex was 4.8° less closed than the glucose complex, which may explain the low efficiency of GCP transfer. Modeling of either α- or β-glucose or a mixture of both anomers can account for the observed electron density of the UDP-glucose complex. NMR studies of UDP-Glc hydrolysis by yeast glycogen synthase were used to verify the stereochemistry of the product, and they also showed synchronous GCP accumulation. The similarities in the active sites of glycogen synthase and glycogen phosphorylase support the idea of a common catalytic mechanism in GT-B enzymes independent of the specific reaction catalyzed.

Classification of domain movements in proteins using dynamic contact graphs

A new method for the classification of domain movements in proteins is described and applied to 1822 pairs of structures from the Protein Data Bank that represent a domain movement in two-domain proteins. The method is based on changes in contacts between residues from the two domains in moving from one conformation to the other. We argue that there are five types of elemental contact changes and that these relate to five model domain movements called: “free”, “open-closed”, “anchored”, “sliding-twist”, and “see-saw.” A directed graph is introduced called the “Dynamic Contact Graph” which represents the contact changes in a domain movement. In many cases a graph, or part of a graph, provides a clear visual metaphor for the movement it represents and is a motif that can be easily recognised. The Dynamic Contact Graphs are often comprised of disconnected subgraphs indicating independent regions which may play different roles in the domain movement. The Dynamic Contact Graph for each domain movement is decomposed into elemental Dynamic Contact Graphs, those that represent elemental contact changes, allowing us to count the number of instances of each type of elemental contact change in the domain movement. This naturally leads to sixteen classes into which the 1822 domain movements are classified.

Protein conformational switches: from nature to design

Protein conformational switches alter their shape upon receiving an input signal, such as ligand binding, chemical modification, or change in environment. The apparent simplicity of this transformation--which can be carried out by a molecule as small as a thousand atoms or so--belies its critical importance to the life of the cell as well as its capacity for engineering by humans. In the realm of molecular switches, proteins are unique because they are capable of performing a variety of biological functions. Switchable proteins are therefore of high interest to the fields of biology, biotechnology, and medicine. These molecules are beginning to be exploited as the core machinery behind a new generation of biosensors, functionally regulated enzymes, and "smart" biomaterials that react to their surroundings. As inspirations for these designs, researchers continue to analyze existing examples of allosteric proteins. Recent years have also witnessed the development of new methodologies for introducing conformational change into proteins that previously had none. Herein we review examples of both natural and engineered protein switches in the context of four basic modes of conformational change: rigid-body domain movement, limited structural rearrangement, global fold switching, and folding-unfolding. Our purpose is to highlight examples that can potentially serve as platforms for the design of custom switches. Accordingly, we focus on inducible conformational changes that are substantial enough to produce a functional response (e.g., in a second protein to which it is fused), yet are relatively simple, structurally well-characterized, and amenable to protein engineering efforts.

PSCDB: a database for protein structural change upon ligand binding

Proteins are flexible molecules that undergo structural changes to function. The Protein Data Bank contains multiple entries for identical proteins determined under different conditions, e.g. with and without a ligand molecule, which provides important information for understanding the structural changes related to protein functions. We gathered 839 protein structural pairs of ligand-free and ligand-bound states from monomeric or homo-dimeric proteins, and constructed the Protein Structural Change DataBase (PSCDB). In the database, we focused on whether the motions were coupled with ligand binding. As a result, the protein structural changes were classified into seven classes, i.e. coupled domain motion (59 structural changes), independent domain motion (70), coupled local motion (125), independent local motion (135), burying ligand motion (104), no significant motion (311) and other type motion (35). PSCDB provides lists of each class. On each entry page, users can view detailed information about the motion, accompanied by a morphing animation of the structural changes. PSCDB is available at http://idp1.force.cs.is.nagoya-u.ac.jp/pscdb/.

Structural analysis of a putative aminoglycoside N-acetyltransferase from Bacillus anthracis

For the last decade, worldwide efforts for the treatment of anthrax infection have focused on developing effective vaccines. Patients that are already infected are still treated traditionally using different types of standard antimicrobial agents. The most popular are antibiotics such as tetracyclines and fluoroquinolones. While aminoglycosides appear to be less effective antimicrobial agents than other antibiotics, synthetic aminoglycosides have been shown to act as potent inhibitors of anthrax lethal factor and may have potential application as antitoxins. Here, we present a structural analysis of the BA2930 protein, a putative aminoglycoside acetyltransferase, which may be a component of the bacterium's aminoglycoside resistance mechanism. The determined structures revealed details of a fold characteristic only for one other protein structure in the Protein Data Bank, namely, YokD from Bacillus subtilis. Both BA2930 and YokD are members of the Antibiotic_NAT superfamily (PF02522). Sequential and structural analyses showed that residues conserved throughout the Antibiotic_NAT superfamily are responsible for the binding of the cofactor acetyl coenzyme A. The interaction of BA2930 with cofactors was characterized by both crystallographic and binding studies.

Classification and annotation of the relationship between protein structural change and ligand binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank-one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.

Prediction of protein motions from amino acid sequence and its application to protein-protein interaction

BACKGROUND

Structural flexibility is an important characteristic of proteins because it is often associated with their function. The movement of a polypeptide segment in a protein can be broken down into two types of motions: internal and external ones. The former is deformation of the segment itself, but the latter involves only rotational and translational motions as a rigid body. Normal Model Analysis (NMA) can derive these two motions, but its application remains limited because it necessitates the gathering of complete structural information.

RESULTS

In this work, we present a novel method for predicting two kinds of protein motions in ordered structures. The prediction uses only information from the amino acid sequence. We prepared a dataset of the internal and external motions of segments in many proteins by application of NMA. Subsequently, we analyzed the relation between thermal motion assessed from X-ray crystallographic B-factor and internal/external motions calculated by NMA. Results show that attributes of amino acids related to the internal motion have different features from those related to the B-factors, although those related to the external motion are correlated strongly with the B-factors. Next, we developed a method to predict internal and external motions from amino acid sequences based on the Random Forest algorithm. The proposed method uses information associated with adjacent amino acid residues and secondary structures predicted from the amino acid sequence. The proposed method exhibited moderate correlation between predicted internal and external motions with those calculated by NMA. It has the highest prediction accuracy compared to a naïve model and three published predictors.

CONCLUSIONS

Finally, we applied the proposed method predicting the internal motion to a set of 20 proteins that undergo large conformational change upon protein-protein interaction. Results show significant overlaps between the predicted high internal motion regions and the observed conformational change regions.

The high-resolution structure of pig heart succinyl-CoA:3-oxoacid coenzyme A transferase

The enzyme succinyl-CoA:3-oxoacid coenzyme A transferase (SCOT) participates in the metabolism of ketone bodies in extrahepatic tissues. It catalyses the transfer of coenzyme A (CoA) from succinyl-CoA to acetoacetate with a classical ping-pong mechanism. There is biochemical evidence that the enzyme undergoes conformational changes during the reaction, but no domain movements have been reported in the available crystal structures. Here, a structure of pig heart SCOT refined at 1.5 A resolution is presented, showing that one of the four enzyme subunits in the crystallographic asymmetric unit has a molecule of glycerol bound in the active site; the glycerol molecule is hydrogen bonded to the conserved catalytic glutamate residue and is likely to occupy the cosubstrate-binding site. The binding of glycerol is associated with a substantial relative movement (a 13 degrees rotation) of two previously undefined domains that close around the substrate-binding site. The binding orientation of one of the cosubstrates, acetoacetate, is suggested based on the glycerol binding and the possibility that this dynamic domain movement is of functional importance is discussed.

Fast hinge detection algorithms for flexible protein structures

Analysis of conformational changes is one of the keys to the understanding of protein functions and interactions. For the analysis, we often compare two protein structures, taking flexible regions like hinge regions into consideration. The Root Mean Square Deviation (RMSD) is the most popular measure for comparing two protein structures, but it is only for rigid structures without hinge regions. In this paper, we propose a new measure called RMSD considering hinges (RMSDh) and its variant RMSDh(k) for comparing two flexible proteins with hinge regions. We also propose novel efficient algorithms for computing them, which can detect the hinge positions at the same time. The RMSDh is suitable for cases where there is one small hinge region in each of the two target structures. The new algorithm for computing the RMSDh runs in linear time, which is the same as the time complexity for computing the RMSD and is faster than any of previous algorithms for hinge detection. The RMSDh(k) is designed for comparing structures with more than one hinge region. The RMSDh(k) measure considers at most k small hinge region, i.e., the RMSDh(k) value should be small if the two structures are similar except for at most k hinge regions. To compute the value, we propose an O(kn2)-time and O(n)-space algorithm based on a new dynamic programming technique. With the same computational time and space, we can enumerate the predicted hinge positions. We also test our algorithms against actual flexible protein structures, and show that the hinge positions can be correctly detected by our algorithms.

Reduced representation of protein structure: implications on efficiency and scope of detection of structural similarity

Background

Computational comparison of two protein structures is the starting point of many methods that build on existing knowledge, such as structure modeling (including modeling of protein complexes and conformational changes), molecular replacement, or annotation by structural similarity. In a commonly used strategy, significant effort is invested in matching two sets of atoms. In a complementary approach, a global descriptor is assigned to the overall structure, thus losing track of the substructures within.

Results

Using a small set of geometric features, we define a reduced representation of protein structure, together with an optimizing function for matching two representations, to provide a pre-filtering stage in a database search. We show that, in a straightforward implementation, the representation performs well in terms of resolution in the space of protein structures, and its ability to make new predictions.

Conclusions

Perhaps unexpectedly, a substantial discriminating power already exists at the level of main features of protein structure, such as directions of secondary structural elements, possibly constrained by their sequential order. This can be used toward efficient comparison of protein (sub)structures, allowing for various degrees of conformational flexibility within the compared pair, which in turn can be used for modeling by homology of protein structure and dynamics.

Crystal structures of the ATPase domains of four human Hsp70 isoforms: HSPA1L/Hsp70-hom, HSPA2/Hsp70-2, HSPA6/Hsp70B', and HSPA5/BiP/GRP78

The 70-kDa heat shock proteins (Hsp70) are chaperones with central roles in processes that involve polypeptide remodeling events. Hsp70 proteins consist of two major functional domains: an N-terminal nucleotide binding domain (NBD) with ATPase activity, and a C-terminal substrate binding domain (SBD). We present the first crystal structures of four human Hsp70 isoforms, those of the NBDs of HSPA1L, HSPA2, HSPA5 and HSPA6. As previously with Hsp70 family members, all four proteins crystallized in a closed cleft conformation, although a slight cleft opening through rotation of subdomain IIB was observed for the HSPA5-ADP complex. The structures presented here support the view that the NBDs of human Hsp70 function by conserved mechanisms and contribute little to isoform specificity, which instead is brought about by the SBDs and by accessory proteins.

ENHANCED VERSION

This article can also be viewed as an enhanced version in which the text of the article is integrated with interactive 3D representations and animated transitions. Please note that a web plugin is required to access this enhanced functionality. Instructions for the installation and use of the web plugin are available in Text S1.

First crystal structure of L-lysine 6-dehydrogenase as an NAD-dependent amine dehydrogenase

A gene encoding an L-lysine dehydrogenase was identified in the hyperthermophilic archaeon Pyrococcus horikoshii. The gene was overexpressed in Escherichia coli, and its product was purified and characterized. The expressed enzyme is the most thermostable L-lysine dehydrogenase yet described, with a half-life of 180 min at 100 degrees C. The product of the enzyme's catalytic activity is Delta(1)-piperideine-6-carboxylate, which makes this enzyme an L-lysine 6-dehydrogenase (EC 1.4.1.18) that catalyzes the reductive deamination of the epsilon- amino group and a type of NAD-dependent amine dehydrogenase. The three-dimensional structure of the enzyme was determined using the mercury-based multiple-wavelength anomalous dispersion method at a resolution of 2.44 A in the presence of NAD and sulfate ion. The asymmetric unit consisted of two subunits, and a crystallographic 2-fold axis generated the functional dimer. Each monomer consisted of a Rossmann fold domain and a C-terminal catalytic domain, and the fold of the catalytic domain showed similarity to that of saccharopine reductase. Notably, the structures of subunits A and B differed significantly. In subunit A, the active site contained a sulfate ion that was not seen in subunit B. Consequently, subunit A adopted a closed conformation, whereas subunit B adopted an open one. In each subunit, one NAD molecule was bound to the active site in an anti-conformation, indicating that the enzyme makes use of pro-R-specific hydride transfer between the two hydrides at C-4 of NADH (type A specificity). This is the first description of the three-dimensional structure of l-lysine 6-dehydrogenase as an NAD-dependent amine dehydrogenase.

FlexSnap: flexible non-sequential protein structure alignment

Background

Proteins have evolved subject to energetic selection pressure for stability and flexibility. Structural similarity between proteins that have gone through conformational changes can be captured effectively if flexibility is considered. Topologically unrelated proteins that preserve secondary structure packing interactions can be detected if both flexibility and Sequential permutations are considered. We propose the FlexSnap algorithm for flexible non-topological protein structural alignment.

Results

The effectiveness of FlexSnap is demonstrated by measuring the agreement of its alignments with manually curated non-sequential structural alignments. FlexSnap showed competitive results against state-of-the-art algorithms, like DALI, SARF2, MultiProt, FlexProt, and FATCAT. Moreover on the DynDom dataset, FlexSnap reported longer alignments with smaller rmsd.

Conclusions

We have introduced FlexSnap, a greedy chaining algorithm that reports both sequential and non-sequential alignments and allows twists (hinges). We assessed the quality of the FlexSnap alignments by measuring its agreements with manually curated non-sequential alignments. On the FlexProt dataset, FlexSnap was competitive to state-of-the-art flexible alignment methods. Moreover, we demonstrated the benefits of introducing hinges by showing significant improvements in the alignments reported by FlexSnap for the structure pairs for which rigid alignment methods reported alignments with either low coverage or large rmsd.

Availability

An implementation of the FlexSnap algorithm will be made available online at http://www.cs.rpi.edu/~zaki/software/flexsnap.