Revisiting gap locations in amino acid sequence alignments and a proposal for a method to improve them by introducing solvent accessibility

PC_ali: a tool for improved multiple alignments and evolutionary inference based on a hybrid protein sequence and structure similarity score

MOTIVATION

Evolutionary inference depends crucially on the quality of multiple sequence alignments (MSA), which is problematic for distantly related proteins. Since protein structure is more conserved than sequence, it seems natural to use structure alignments for distant homologs. However, structure alignments may not be suitable for inferring evolutionary relationships.

RESULTS

Here we examined four protein similarity measures that depend on sequence and structure (fraction of aligned residues, sequence identity, fraction of superimposed residues, and contact overlap), finding that they are intimately correlated but none of them provides a complete and unbiased picture of conservation in proteins. Therefore, we propose the new hybrid protein sequence and structure similarity score PC_sim based on their main principal component. The corresponding divergence measure PC_div shows the strongest correlation with divergences obtained from individual similarities, suggesting that it infers accurate evolutionary divergences. We developed the program PC_ali that constructs protein MSAs either de novo or modifying an input MSA, using a similarity matrix based on PC_sim. The program constructs a starting MSA based on the maximal cliques of the graph of these PAs and it refines it through progressive alignments along the tree reconstructed with PC_div. Compared with eight state-of-the-art multiple structure or sequence alignment tools, PC_ali achieves higher or equal aligned fraction and structural scores, sequence identity higher than structure aligners although lower than sequence aligners, highest score PC_sim, and highest similarity with the MSAs produced by other tools and with the reference MSA Balibase.

AVAILABILITY AND IMPLEMENTATION

https://github.com/ugobas/PC_ali.

Computational study of the impact of nucleotide variations on highly conserved proteins: In the case of actin

Sequencing of individual human genomes enables studying relationship among nucleotide variations, amino acid substitutions, effect on protein structures and diseases. Many studies have found general tendencies, for instance, that pathogenic variations tend to be found in the buried regions of the protein structures, that benign variations tend to be found on the surface of the proteins, and that variations on evolutionary conserved residues tend to be pathogenic. These tendencies were deduced from globular proteins with standard evolutionary changes in amino acid sequences. In this study, we investigated the variation distribution on actin, one of the highly conserved proteins. Many nucleotide variations and three-dimensional structures of actin have been registered in databases. By combining those data, we found that variations buried inside the protein were rather benign and variations on the surface of the protein were pathogenic. This idiosyncratic distribution of the variation impact is likely ascribed to the extensive use of the surface of the protein for protein-protein interactions in actin.

Multiple mutations in RNA polymerase β-subunit gene (rpoB) in Streptomyces incarnatus NRRL8089 enhance production of antiviral antibiotic sinefungin: modeling rif cluster region by density functional theory

Streptomyces incarnatus NRRL8089 produces the antiviral, antifungal, antiprotozoal nucleoside antibiotic sinefungin. To enhance sinefungin production, multiple mutations were introduced to the rpoB gene encoding RNA polymerase (RNAP) β-subunit at the target residues, D447, S453, H457, and R460. Sparse regression analysis using elastic-net lasso-ridge penalties on previously reported H457X mutations identified a numeric parameter set, which suggested that H457R/Y/F may cause production enhancement. H457R/R460C mutation successfully enhanced the sinefungin production by 3-fold, while other groups of mutations, such as D447G/R460C or D447G/H457Y, made moderate or even negative effects. To identify why the rif cluster residues have diverse effects on sinefungin production, an RNAP/DNA/mRNA complex model was constructed by homology modeling and molecular dynamics simulation. The 4 residues were located near the mRNA strand. Density functional theory-based calculation suggested that D447, H457, and R460 are in direct contact with ribonucleotide, and partially positive charges are induced by negatively charged chain of mRNA.

Genomic and Transcriptomic Analyses of Bioluminescence Genes in the Enope Squid Watasenia scintillans

Watasenia scintillans, a sparkling enope squid, has bioluminescence organs to illuminate its body with its own luciferase activity. To clarify the molecular mechanism underlying its scintillation, we analysed high-throughput sequencing data acquired previously and obtained draft genome sequences accomplished with comparative genomic data among the cephalopods. The genome mapped by transcriptome data showed that (1) RNA editing contributed to transcriptome variation of lineage specific genes, such as W. scintillans luciferase, and (2) two types of luciferase enzymes were characterized with reasonable 3D models docked to a luciferin molecule. We report two different types of luciferase in one organism and possibly related to variety of colour types in the W. scintillans fluorescent organs.

Sequence alignment using machine learning for accurate template-based protein structure prediction

MOTIVATION

Template-based modeling, the process of predicting the tertiary structure of a protein by using homologous protein structures, is useful if good templates can be found. Although modern homology detection methods can find remote homologs with high sensitivity, the accuracy of template-based models generated from homology-detection-based alignments is often lower than that from ideal alignments.

RESULTS

In this study, we propose a new method that generates pairwise sequence alignments for more accurate template-based modeling. The proposed method trains a machine learning model using the structural alignment of known homologs. It is difficult to directly predict sequence alignments using machine learning. Thus, when calculating sequence alignments, instead of a fixed substitution matrix, this method dynamically predicts a substitution score from the trained model. We evaluate our method by carefully splitting the training and test datasets and comparing the predicted structure's accuracy with that of state-of-the-art methods. Our method generates more accurate tertiary structure models than those produced from alignments obtained by other methods.

AVAILABILITY AND IMPLEMENTATION

https://github.com/shuichiro-makigaki/exmachina.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TAZ encodes tafazzin, a transacylase essential for cardiolipin formation and central to the etiology of Barth syndrome

Tafazzin, which is encoded by the TAZ gene, catalyzes transacylation to form mature cardiolipin and shows preference for the transfer of a linoleic acid (LA) group from phosphatidylcholine (PC) to monolysocardiolipin (MLCL) with influence from mitochondrial membrane curvature. The protein contains domains and motifs involved in targeting, anchoring, and an active site for transacylase activity. Tafazzin activity affects many aspects of mitochondrial structure and function, including that of the electron transport chain, fission-fusion, as well as apoptotic signaling. TAZ mutations are implicated in Barth syndrome, an underdiagnosed and devastating disease that primarily affects male pediatric patients with a broad spectrum of disease pathologies that impact the cardiovascular, neuromuscular, metabolic, and hematologic systems.

SGK regulates pH increase and cyclin B-Cdk1 activation to resume meiosis in starfish ovarian oocytes

Tight regulation of intracellular pH (pH_i) is essential for biological processes. Fully grown oocytes, having a large nucleus called the germinal vesicle, arrest at meiotic prophase I. Upon hormonal stimulus, oocytes resume meiosis to become fertilizable. At this time, the pH_i increases via Na⁺/H⁺ exchanger activity, although the regulation and function of this change remain obscure. Here, we show that in starfish oocytes, serum- and glucocorticoid-regulated kinase (SGK) is activated via PI3K/TORC2/PDK1 signaling after hormonal stimulus and that SGK is required for this pH_i increase and cyclin B-Cdk1 activation. When we clamped the pH_i at 6.7, corresponding to the pH_i of unstimulated ovarian oocytes, hormonal stimulation induced cyclin B-Cdk1 activation; thereafter, oocytes failed in actin-dependent chromosome transport and spindle assembly after germinal vesicle breakdown. Thus, this SGK-dependent pH_i increase is likely a prerequisite for these events in ovarian oocytes. We propose a model that SGK drives meiotic resumption via concomitant regulation of the pH_i and cell cycle machinery.

Identification of a High-Frequency Somatic NLRC4 Mutation as a Cause of Autoinflammation by Pluripotent Cell-Based Phenotype Dissection

OBJECTIVE

To elucidate the genetic background of a patient with neonatal-onset multisystem inflammatory disease (NOMID) with no NLRP3 mutation.

METHODS

A Japanese male child diagnosed as having NOMID was studied. The patient did not have any NLRP3 mutation, even as low-frequency mosaicism. We performed whole-exome sequencing on the patient and his parents. Induced pluripotent stem cells (iPSCs) were established from the patient's fibroblasts. The iPSCs were then differentiated into monocyte lineage to evaluate the cytokine profile.

RESULTS

We established multiple iPSC clones from a patient with NOMID and incidentally found that the phenotypes of monocytes from iPSC clones were heterogeneous and could be grouped into disease and normal phenotypes. Because each iPSC clone was derived from a single somatic cell, we hypothesized that the patient had somatic mosaicism of an interleukin-1β-related gene. Whole-exome sequencing of both representative iPSC clones and the patient's blood revealed a novel heterozygous NLRC4 mutation, p.T177A (c.529A>G), as a specific mutation in diseased iPSC clones. Knockout of the NLRC4 gene using the clustered regularly interspaced short palindromic repeat/Cas9 system in a mutant iPSC clone abrogated the pathogenic phenotype.

CONCLUSION

Our findings indicate that the patient has somatic mosaicism of a novel NLRC4 mutation. To our knowledge, this is the first case showing that somatic mutation of NLRC4 causes autoinflammatory symptoms compatible with NOMID. The present study demonstrates the significance of prospective genetic screening combined with iPSC-based phenotype dissection for individualized diagnoses.

Parallel evolution of non-homologous isofunctional enzymes in methionine biosynthesis

Experimental validation of enzyme function is crucial for genome interpretation, but it remains challenging because it cannot be scaled up to accommodate the constant accumulation of genome sequences. We tackled this issue for the MetA and MetX enzyme families, phylogenetically unrelated families of acyl-L-homoserine transferases involved in L-methionine biosynthesis. Members of these families are prone to incorrect annotation because MetX and MetA enzymes are assumed to always use acetyl-CoA and succinyl-CoA, respectively. We determined the enzymatic activities of 100 enzymes from diverse species, and interpreted the results by structural classification of active sites based on protein structure modeling. We predict that >60% of the 10,000 sequences from these families currently present in databases are incorrectly annotated, and suggest that acetyl-CoA was originally the sole substrate of these isofunctional enzymes, which evolved to use exclusively succinyl-CoA in the most recent bacteria. We also uncovered a divergent subgroup of MetX enzymes in fungi that participate only in L-cysteine biosynthesis as O-succinyl-L-serine transferases.

A hydrophobic spine stabilizes a surface-exposed α-helix according to analysis of the solvent-accessible surface area

Background

Most of hydrophilic and hydrophobic residues are thought to be exposed and buried in proteins, respectively. In contrast to the majority of the existing studies on protein folding characteristics using protein structures, in this study, our aim was to design predictors for estimating relative solvent accessibility (RSA) of amino acid residues to discover protein folding characteristics from sequences.

Methods

The proposed 20 real-value RSA predictors were designed on the basis of the support vector regression method with a set of informative physicochemical properties (PCPs) obtained by means of an optimal feature selection algorithm. Then, molecular dynamics simulations were performed for validating the knowledge discovered by analysis of the selected PCPs.

Results

The RSA predictors had the mean absolute error of 14.11% and a correlation coefficient of 0.69, better than the existing predictors. The hydrophilic-residue predictors preferred PCPs of buried amino acid residues to PCPs of exposed ones as prediction features. A hydrophobic spine composed of exposed hydrophobic residues of an α-helix was discovered by analyzing the PCPs of RSA predictors corresponding to hydrophobic residues. For example, the results of a molecular dynamics simulation of wild-type sequences and their mutants showed that proteins 1MOF and 2WRP_H16I (Protein Data Bank IDs), which have a perfectly hydrophobic spine, have more stable structures than 1MOF_I54D and 2WRP do (which do not have a perfectly hydrophobic spine).

Conclusions

We identified informative PCPs to design high-performance RSA predictors and to analyze these PCPs for identification of novel protein folding characteristics. A hydrophobic spine in a protein can help to stabilize exposed α-helices.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1368-z) contains supplementary material, which is available to authorized users.

Novel AVPR2 mutation causing partial nephrogenic diabetes insipidus in a Japanese family

BACKGROUND

X-linked recessive congenital nephrogenic diabetes insipidus (NDI) is caused by mutations of the arginine vasopressin type 2 receptor gene (AVPR2). More than 200 mutations of the AVPR2 gene with complete NDI have been reported although only 15 mutations with partial NDI has been reported to date.

METHODS

We herein report a Japanese kindred with partial NDI. The proband is an 8-year-old boy who was referred to our hospital for nocturnal enuresis. Water deprivation test and hypertonic saline test suggested partial renal antidiuretic hormone arginine vasopressin (AVP) resistance.

RESULTS

Analysis of genomic DNA revealed a novel missense mutation (p.L161P) in the patient. The patient's mother was heterozygous for the mutation. Three-dimensional (3-D) modeling study showed that L161P possibly destabilizes the transmembrane domain of the V2 receptor, resulting in its misfolding or mislocalization.

CONCLUSIONS

Distinguishing partial NDI from nocturnal enuresis is important. A clinical clue for diagnosis of partial NDI is an incompatibly high level of AVP despite normal serum osmolality.

Molecular evolution of gas cavity in [NiFeSe] hydrogenases resurrected in silico

Oxygen tolerance of selenium-containing [NiFeSe] hydrogenases (Hases) is attributable to the high reducing power of the selenocysteine residue, which sustains the bimetallic Ni–Fe catalytic center in the large subunit. Genes encoding [NiFeSe] Hases are inherited by few sulphate-reducing δ-proteobacteria globally distributed under various anoxic conditions. Ancestral sequences of [NiFeSe] Hases were elucidated and their three-dimensional structures were recreated in silico using homology modelling and molecular dynamic simulation, which suggested that deep gas channels gradually developed in [NiFeSe] Hases under absolute anaerobic conditions, whereas the enzyme remained as a sealed edifice under environmental conditions of a higher oxygen exposure risk. The development of a gas cavity appears to be driven by non-synonymous mutations, which cause subtle conformational changes locally and distantly, even including highly conserved sequence regions.

Case study on the evolution of hetero-oligomer interfaces based on the differences in paralogous proteins

We addressed the evolutionary trace of hetero-oligomer interfaces by comparing the structures of paralogous proteins; one of them is a monomer or homo-oligomer and the other is a hetero-oligomer. We found different trends in amino acid conservation pattern and hydrophobicity between homo-oligomer and hetero-oligomer. The degree of amino acid conservation in the interface of homo-oligomer has no obvious difference from that in the surface, whereas the degree of conservation is much higher in the interface of hetero-oligomer. The interface of homo-oligomer has a few very conserved residue positions, whereas the residue conservation in the interface of hetero-oligomer tends to be higher. In addition, the interface of hetero-oligomer has a tendency of being more hydrophobic compared with the one in homo-oligomer. We conjecture that these differences are related to the inherent symmetry in homo-oligomers that cannot exist in hetero-oligomers. Paucity of the structural data precludes statistical tests of these tendencies, yet the trend can be applied to the prediction of the interface of hetero-oligomer. We obtained putative interfaces of the subunits in CPSF (cleavage and polyadenylation specificity factor), one of the human pre-mRNA 3′-processing complexes. The locations of predicted interface residues were consistent with the known experimental data.

Structural and functional analyses of Barth syndrome-causing mutations and alternative splicing in the tafazzin acyltransferase domain

Tafazzin is a mitochondrial phospholipid transacylase, and its mutations cause Barth syndrome (BTHS). Human tafazzin gene produces four distinct alternatively spliced transcripts. To understand the molecular mechanisms of tafazzin deficiency, we performed an atomic resolution analysis of the influence of the BTHS mutations and of alternative splicing on the structure and function of tafazzin. From the three-dimensional (3D) homology modeling of tafazzin, we identified candidate amino acid residues that contribute to cardiolipin binding and to mitochondrial membrane associations that facilitate acyl-transfer reactions. Primate specific exon 5, which is alternatively spliced, is predicted to correspond to an intrinsically unstructured region in the protein. We proposed that this region should change the substrate-binding affinity and/or contribute to primate-specific molecular interactions. Exon 7, another alternatively spliced exon, encodes a region forming a part of the putative substrate-binding cleft, suggesting that the gene products lacking exon 7 will lose their substrate-binding ability. We demonstrate a clear localization of the BTHS mutations at residues responsible for membrane association, substrate binding, and the conformational stability of tafazzin. These findings provide new insights into the function of defective tafazzin and the pathogenesis of BTHS at the level of protein 3D structure and the evolution of alternatively spliced exons in primates.

Distinct evolutionary rate in the eye field transcription factors found by estimation of ancestral protein structure

Eye-field transcription factors (EFTFs) are a set of genes that compose a regulatory network for eye development in animals, which are highly conserved among various animal phyla. To investigate the processes of conservation and diversification of the transcription factors for eye development, we examined the structural changes in the EFTF proteins by estimating the ancestral sequences with the available genome information. Among the different types of EFTFs, we selected otx2, tbx3, rx1, pax6, six3/6, lhx2 and nr2e1 because they are highly conserved in bilaterian animals. We searched the genome sequences of representative animal phyla for EFTF protein sequences. With deduced ancestral sequences and three-dimensional structures of EFTFs, we traced the evolutionary changes in amino acid residues and found that the DNA-binding domains were always more conserved than other regions, and that the other regions showed distinct evolutionary rates. The EFTF rx1, which resides at the pivotal part of the EFTF network, had a faster evolutionary rate than the others. These results indicated that the evolutionary rates of each protein in the EFTF network, which were expected to be consistent with each other to maintain the interactions in the network, were not constant among or within the factors, but rather, varied to a significant extent.

Two novel HSD3B2 missense mutations with diverse residual enzymatic activities for Δ5-steroids

CONTEXT

Classical 3β-hydroxysteroid dehydrogenase (3β-HSD) deficiency (3β-HSDD) is caused by loss-of-function mutations in the HSD3B2 gene encoding type II 3β-HSD, which has a key role in steroid biosynthesis, converting Δ5-steroids to Δ4-steroids in adrenal glands and gonads.

PATIENT

A patient (46, XX) was found to have elevated 17-hydroxyprogesterone (17-OHP) [203 nmol/l (normal range: 2·94 ± 0·9 nmol/l)] by newborn screening. Endocrinological examination revealed dramatically increased Δ5-steroids [e.g. 17-OH pregnenolone: 910 nmol/l (normal range: 12·6 ± 10·5 nmol/l)]. The patient had virilization of external genitalia with labial fusion, suggesting classical 3β-HSDD.

METHODS AND RESULTS

Consistent with the endocrinological data, the patient was a compound heterozygote for two novel missense mutations (p.Y190C and p.S218P) that were identified in HSD3B2. Both Y190 and S218 are conserved among mammals. The mutant proteins had severely impaired residual enzymatic activity in vitro, although both mutants retained higher activity for 17-OH pregnenolone than for the other Δ5-steroids. In a three-dimensional model of the enzyme based on the known structures of similar proteins, both mutations were located extremely close to the predicted substrate-binding pocket. This suggests that the mutations can cause a local conformational change in the substrate-binding pocket, leading to alterations of the binding affinities for Δ5-steroids.

CONCLUSIONS

We identified two novel missense mutations of HSD3B2 that resulted in unbalanced residual enzymatic activities for Δ5-steroids. As a potential novel mechanism, we propose that the mutations, which differently affect the activity towards different substrates, the effects of these mutations provide novel insights into the pathophysiology of 3β-HSDD.

X-ray structure analysis and characterization of AFUEI, an elastase inhibitor from Aspergillus fumigatus

Elastase from Aspergillus sp. is an important factor for aspergillosis. AFUEI is an inhibitor of the elastase derived from Aspergillus fumigatus. AFUEI is a member of the I78 inhibitor family and has a high inhibitory activity against elastases of Aspergillus fumigatus and Aspergillus flavus, human neutrophil elastase and bovine chymotrypsin, but does not inhibit bovine trypsin. Here we report the crystal structure of AFUEI in two crystal forms. AFUEI is a wedge-shaped protein composed of an extended loop and a scaffold protein core. The structure of AFUEI shows remarkable similarity to serine protease inhibitors of the potato inhibitor I family, although they are classified into different inhibitor families. A structural comparison with the potato I family inhibitors suggests that the extended loop of AFUEI corresponds to the binding loop of the potato inhibitor I family, and AFUEI inhibits its cognate proteases through the same mechanism as the potato I family inhibitors.

A bicyclic 1-deoxygalactonojirimycin derivative as a novel pharmacological chaperone for GM1 gangliosidosis

Lysosomal β-galactosidase (β-Gal) deficiency causes a group of disorders that include neuronopathic GM1 gangliosidosis and non-neuronopathic Morquio B disease. We have previously proposed the use of small molecule ligands of β-Gal as pharmacological chaperones (PCs) for the treatment of GM1 gangliosidosis brain pathology. Although it is still under development, PC therapy has yielded promising preclinical results in several lysosomal diseases. In this study, we evaluated the effect of bicyclic 1-deoxygalactonojirimycin (DGJ) derivative of the sp(2)-iminosugar type, namely 5N,6S-(N'-butyliminomethylidene)-6-thio-1- deoxygalactonojirimycin (6S-NBI-DGJ), as a novel PC for human mutant β-Gal. In vitro, 6S-NBI-DGJ had the ability to inhibit the activity of human β-Gal in a competitive manner and was able to protect this enzyme from heat-induced degradation. Computational analysis supported that the rigid glycone bicyclic core of 6S-NBI-DGJ binds to the active site of the enzyme, with the aglycone N'-butyl substituent, in a precise E-orientation, located at a hydrophobic region nearby. Chaperone potential profiling indicated significant increases of enzyme activity in 24 of 88 β-Gal mutants, including four common mutations. Finally, oral administration of 6S-NBI-DGJ ameliorated the brain pathology of GM1 gangliosidosis model mice. These results suggest that 6S-NBI-DGJ is a novel PC that may be effective on a broad range of β-Gal mutants.

Structure of N-formylglycinamide ribonucleotide amidotransferase II (PurL) from Thermus thermophilus HB8

The crystal structure of PurL from Thermus thermophilus HB8 (TtPurL; TTHA1519) was determined in complex with an adenine nucleotide, PO(4)(3-) and Mg(2+) at 2.35 Å resolution. TtPurL consists of 29 α-helices and 28 β-strands, and one loop is disordered. TtPurL consists of four domains, A1, A2, B1 and B2, and the structures of the A1-B1 and A2-B2 domains were almost identical to each other. Although the sequence identity between TtPurL and PurL from Thermotoga maritima (TmPurL) is higher than that between TtPurL and the PurL domain of the large PurL from Salmonella typhimurium (StPurL), the secondary structure of TtPurL is much more similar to that of StPurL than to that of TmPurL.

Substrate recognition mechanism and substrate-dependent conformational changes of an ROK family glucokinase from Streptomyces griseus

Carbon catabolite repression (CCR) is a widespread phenomenon in many bacteria that is defined as the repression of catabolic enzyme activities for an unfavorable carbon source by the presence of a preferable carbon source. In Streptomyces, secondary metabolite production often is negatively affected by the carbon source, indicating the involvement of CCR in secondary metabolism. Although the CCR mechanism in Streptomyces still is unclear, glucokinase is presumably a central player in CCR. SgGlkA, a glucokinase from S. griseus, belongs to the ROK family glucokinases, which have two consensus sequence motifs (1 and 2). Here, we report the crystal structures of apo-SgGlkA, SgGlkA in complex with glucose, and SgGlkA in complex with glucose and adenylyl imidodiphosphate (AMPPNP), which are the first structures of an ROK family glucokinase. SgGlkA is divided into a small α/β domain and a large α+β domain, and it forms a dimer-of-dimer tetrameric configuration. SgGlkA binds a β-anomer of glucose between the two domains, and His157 in consensus sequence 1 plays an important role in the glucose-binding mechanism and anomer specificity of SgGlkA. In the structures of SgGlkA, His157 forms an HC3-type zinc finger motif with three cysteine residues in consensus sequence 2 to bind a zinc ion, and it forms two hydrogen bonds with the C1 and C2 hydroxyls of glucose. When the three structures are compared, the structure of SgGlkA is found to be modified by the binding of substrates. The substrate-dependent conformational changes of SgGlkA may be related to the CCR mechanism in Streptomyces.