The origin of esterase activity of Parkinson's disease causative factor DJ-1 implied by evolutionary trace analysis of its prokaryotic homolog HchA.

DJ-1, a causative gene for hereditary recessive Parkinsonism, is evolutionarily conserved across eukaryotes and prokaryotes. Structural analyses of DJ-1 and its homologs suggested the 106th Cys is a nucleophilic cysteine functioning as the catalytic center of hydratase or hydrolase activity. Indeed, DJ-1 and its homologs can convert highly electrophilic α-oxoaldehydes such as methylglyoxal into α-hydroxy acids as hydratase in vitro, and oxidation-dependent ester hydrolase (esterase) activity has also been reported for DJ-1. The mechanism underlying such plural activities, however, has not been fully characterized. To address this knowledge gap, we conducted a series of biochemical assays assessing the enzymatic activity of DJ-1 and its homologs. We found no evidence for esterase activity in any of the Escherichia coli DJ-1 homologs. Furthermore, contrary to previous reports, we found that oxidation inactivated rather than facilitated DJ-1 esterase activity. The E. coli DJ-1 homolog HchA possesses phenylglyoxalase and methylglyoxalase activities but lacks esterase activity. Since evolutionary trace analysis identified the 186th H as a candidate residue involved in functional differentiation between HchA and DJ-1, we focused on H186 of HchA and found that an esterase activity was acquired by H186A mutation. Introduction of reverse mutations into the equivalent position in DJ-1 (A107H) selectively eliminated its esterase activity without compromising α-oxoaldehyde hydratase activity. The obtained results suggest that differences in the amino acid sequences near the active site contributed to acquisition of esterase activity in vitro and provide an important clue to the origin and significance of DJ-1 esterase activity.

Atg15 is a vacuolar phospholipase that disintegrates organelle membranes.

Atg15 (autophagy-related 15) is a vacuolar phospholipase essential for the degradation of cytoplasm-to-vacuole targeting (Cvt) bodies and autophagic bodies, hereinafter referred to as intravacuolar/intralysosomal autophagic compartments (IACs), but it remains unknown if Atg15 directly disrupts IAC membranes. Here, we show that the recombinant Chaetomium thermophilum Atg15 lipase domain (CtAtg15(73-475)) possesses phospholipase activity. The activity of CtAtg15(73-475) was markedly elevated by limited digestion. We inserted the human rhinovirus 3C protease recognition sequence and found that cleavage between S159 and V160 was important to activate CtAtg15(73-475). Our molecular dynamics simulation suggested that the cleavage facilitated conformational change around the active center of CtAtg15, resulting in an exposed state. We confirmed that CtAtg15 could disintegrate S. cerevisiae IAC in vivo. Further, both mitochondria and IAC of S. cerevisiae were disintegrated by CtAtg15. This study suggests Atg15 plays a role in disrupting any organelle membranes delivered to vacuoles by autophagy.

A multipoint guidance mechanism for ß-barrel folding on the SAM complex.

Mitochondrial ß-barrel proteins are essential for the transport of metabolites, ions and proteins. The sorting and assembly machinery (SAM) mediates their folding and membrane insertion. We report the cryo-electron microscopy structure of the yeast SAM complex carrying an early eukaryotic ß-barrel folding intermediate. The lateral gate of Sam50 is wide open and pairs with the last ß-strand (ß-signal) of the substrate-the 19-ß-stranded Tom40 precursor-to form a hybrid barrel in the membrane plane. The Tom40 barrel grows and curves, guided by an extended bridge with Sam50. Tom40's first ß-segment (ß1) penetrates into the nascent barrel, interacting with its inner wall. The Tom40 amino-terminal segment then displaces ß1 to promote its pairing with Tom40's last ß-strand to complete barrel formation with the assistance of Sam37's dynamic α-protrusion. Our study thus reveals a multipoint guidance mechanism for mitochondrial ß-barrel folding.

Phosprof: pathway analysis database of drug response based on phosphorylation activity measurements.

Protein phosphorylation plays a fundamental role in many cellular processes. Proteins are phosphorylated by kinases, which have been studied as drug targets for the treatment of various diseases, particularly cancer. Because kinases have multiple roles in interconnected molecular pathways, their specific regulation is required to enhance beneficial and reduce adversarial effects of drugs. Using our previously developed platform, we measured phosphorylation profiles of MCF7 and K562 cells treated with 94 clinical drugs. These phosphorylation profiles can provide insights into pathway activities and biological functions. Here, we introduce Phosprof, a novel database of drug response based on phosphorylation activity. Phosprof is able to present up- or downregulated phosphorylated signature proteins on pathway maps, significant pathways on the hierarchal tree in signal transduction and commonly perturbed pathways affected by the selected drugs. It also serves as a useful web interface for new or known drug profile search based on their molecular similarity with the 94 drugs. Phosprof can be helpful for further investigation of drug responses in terms of phosphorylation by utilizing the various approved drugs whose target phenotypes are known. DATABASE URL: https://phosprof.medals.jp/.

Imputation-free reconstructions of three-dimensional chromosome architectures in human diploid single-cells using allele-specified contacts.

Single-cell Hi-C analysis of diploid human cells is difficult because of the lack of dense chromosome contact information and the presence of homologous chromosomes with very similar nucleotide sequences. Thus here, we propose a new algorithm to reconstruct the three-dimensional (3D) chromosomal architectures from the Hi-C dataset of single diploid human cells using allele-specific single-nucleotide variations (SNVs). We modified our recurrence plot-based algorithm, which is suitable for the estimation of the 3D chromosome structure from sparse Hi-C datasets, by newly incorporating a function of discriminating SNVs specific to each homologous chromosome. Here, we eventually regard a contact map as a recurrence plot. Importantly, the proposed method does not require any imputation for ambiguous segment information, but could efficiently reconstruct 3D chromosomal structures in single human diploid cells at a 1-Mb resolution. Datasets of segments without allele-specific SNVs, which were considered to be of little value, can also be used to validate the estimated chromosome structure. Introducing an additional mathematical measure called a refinement further improved the resolution to 40-kb or 100-kb. The reconstruction data supported the notion that human chromosomes form chromosomal territories and take fractal structures where the dimension for the underlying chromosome structure is a non-integer value.

Unfolding is the driving force for mitochondrial import and degradation of the Parkinson's disease-related protein DJ-1.

Diverse genes associated with familial Parkinson's disease (familial Parkinsonism) have been implicated in mitochondrial quality control. One such gene, PARK7 encodes the protein DJ-1, pathogenic mutations of which trigger its translocation from the cytosol to the mitochondrial matrix. The translocation of steady-state cytosolic proteins like DJ-1 to the mitochondrial matrix upon missense mutations is rare, and the underlying mechanism remains to be elucidated. Here, we show that the protein unfolding associated with various DJ-1 mutations drives its import into the mitochondrial matrix. Increasing the structural stability of these DJ-1 mutants restores cytosolic localization. Mechanistically, we show that a reduction in the structural stability of DJ-1 exposes a cryptic N-terminal mitochondrial-targeting signal (MTS), including Leu10, which promotes DJ-1 import into the mitochondrial matrix for subsequent degradation. Our work describes a novel cellular mechanism for targeting a destabilized cytosolic protein to the mitochondria for degradation.

MDContactCom: a tool to identify differences of protein molecular dynamics from two MD simulation trajectories in terms of interresidue contacts.

SUMMARY: Comparing results from multiple MD simulations performed under different conditions is essential during the initial stages of analysis. We propose a tool called MD Contact Comparison (MDContactCom) that compares residue-residue contact fluctuations of two MD trajectories, quantifies the differences, identifies sites that exhibit large differences and visualizes those sites on the protein structure. Using this method, it is possible to identify sites affected by varying simulation conditions and reveal the path of propagation of the effect even when differences between the 3D structure of the molecule and the fluctuation RMSF of each residue is unclear. MDContactCom can monitor differences in complex protein dynamics between two MD trajectories and identify candidate sites to be analyzed in more detail. As such, MDContactCom is a versatile software package for analyzing most MD simulations. AVAILABILITY AND IMPLEMENTATION: MDContactCom is freely available for download on GitLab. The software is implemented in Python3. https://gitlab.com/chiemotono/mdcontactcom. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Porin Associates with Tom22 to Regulate the Mitochondrial Protein Gate Assembly.

Mitochondria import nearly all of their resident proteins from the cytosol, and the TOM complex functions as their entry gate. The TOM complex undergoes a dynamic conversion between the majority population of a three-channel gateway ("trimer") and the minor population that lacks Tom22 and has only two Tom40 channels ("dimer"). Here, we found that the porin Por1 acts as a sink to bind newly imported Tom22. This Por1 association thereby modulates Tom22 integration into the TOM complex, guaranteeing formation of the functional trimeric TOM complex. Por1 sequestration of Tom22 dissociated from the trimeric TOM complex also enhances the dimeric TOM complex, which is preferable for the import of TIM40/MIA-dependent proteins into mitochondria. Furthermore, Por1 appears to contribute to cell-cycle-dependent variation of the functional trimeric TOM complex by chaperoning monomeric Tom22, which arises from the cell-cycle-controlled variation of phosphorylated Tom6.

SAHG, a comprehensive database of predicted structures of all human proteins.

Most proteins from higher organisms are known to be multi-domain proteins and contain substantial numbers of intrinsically disordered (ID) regions. To analyse such protein sequences, those from human for instance, we developed a special protein-structure-prediction pipeline and accumulated the products in the Structure Atlas of Human Genome (SAHG) database at http://bird.cbrc.jp/sahg. With the pipeline, human proteins were examined by local alignment methods (BLAST, PSI-BLAST and Smith-Waterman profile-profile alignment), global-local alignment methods (FORTE) and prediction tools for ID regions (POODLE-S) and homology modeling (MODELLER). Conformational changes of protein models upon ligand-binding were predicted by simultaneous modeling using templates of apo and holo forms. When there were no suitable templates for holo forms and the apo models were accurate, we prepared holo models using prediction methods for ligand-binding (eF-seek) and conformational change (the elastic network model and the linear response theory). Models are displayed as animated images. As of July 2010, SAHG contains 42,581 protein-domain models in approximately 24,900 unique human protein sequences from the RefSeq database. Annotation of models with functional information and links to other databases such as EzCatDB, InterPro or HPRD are also provided to facilitate understanding the protein structure-function relationships.

Thermodynamic and kinetic determinants of Thermotoga maritima cold shock protein stability: a structural and dynamic analysis.

The cold shock protein (CSP) from hyperthermophile Thermotoga maritima (TmCSP) is only marginally stable (DeltaG(T(opt)) = 0.3 kcal/mol) at 353 K, the optimum environmental temperature (T(opt)) for T. maritima. In comparison, homologous CSPs from E. coli (DeltaG(T(opt)) = 2.2 kcal/mol) and B. subtilis (DeltaG(T(opt)) = 1.5 kcal/mol) are at least five times more stable at 310 K, the T(opt) for the mesophiles. Yet at the room temperature, TmCSP is more stable (DeltaG(T(R)) = 4.7 kcal/mol) than its homologues (DeltaG(T(R)) = 3.0 kcal/mol for E. coli CSP and DeltaG(T(R)) = 2.1 kcal/mol for B. subtilis CSP). This unique observation suggests that kinetic, rather than thermodynamic, barriers toward unfolding might help TmCSP native structure at high temperatures. Consistently, the unfolding rate of TmCSP is considerably slower than its homologues. High temperature (600 K) complete unfolding molecular dynamics (MD) simulations of TmCSP support our hypothesis and reveal an unfolding scheme unique to TmCSP. For all the studied homologues of TmCSP, the unfolding process first starts at the C-terminal region and N-terminal region unfolds in the end. But for TmCSP, both the terminals resist unfolding for consistently longer simulation times and, in the end, unfold simultaneously. In TmCSP, the C-terminal region is better fortified and has better interactions with the N-terminal region due to the charged residues, R2, E47, E49, H61, K63, and E66, being in spatial vicinity. The electrostatic interactions among these residues are unique to TmCSP. Consistently, the room temperature MD simulations show that TmCSP is more rigid at its N- and C-termini as compared to its homologues from E. coli, B. subtilis, and B. caldolyticus.

Large-scale identification and characterization of alternative splicing variants of human gene transcripts using 56,419 completely sequenced and manually annotated full-length cDNAs.

We report the first genome-wide identification and characterization of alternative splicing in human gene transcripts based on analysis of the full-length cDNAs. Applying both manual and computational analyses for 56,419 completely sequenced and precisely annotated full-length cDNAs selected for the H-Invitational human transcriptome annotation meetings, we identified 6877 alternative splicing genes with 18 297 different alternative splicing variants. A total of 37,670 exons were involved in these alternative splicing events. The encoded protein sequences were affected in 6005 of the 6877 genes. Notably, alternative splicing affected protein motifs in 3015 genes, subcellular localizations in 2982 genes and transmembrane domains in 1348 genes. We also identified interesting patterns of alternative splicing, in which two distinct genes seemed to be bridged, nested or having overlapping protein coding sequences (CDSs) of different reading frames (multiple CDS). In these cases, completely unrelated proteins are encoded by a single locus. Genome-wide annotations of alternative splicing, relying on full-length cDNAs, should lay firm groundwork for exploring in detail the diversification of protein function, which is mediated by the fast expanding universe of alternative splicing variants.

Protein structure prediction using a variety of profile libraries and 3D verification.

This study is intended to construct a useful method for fold recognition, regardless of whether the proteins to be compared are evolutionarily related. We developed several descendants of our profile-profile comparison method to make use of known structural information for protein structure prediction. Our prediction strategy in CASP6 is simple. For every CASP6 target, we derived target-template alignments from several different versions of profile-profile comparisons. We then constructed and exhaustively evaluated 3D models based on those alignments. Subsequently, we selected proper model(s) among them. We specifically addressed the validation of our simple approach for protein structure prediction through CASP6 because the fold recognition results of CASP5 revealed areas of improvement in the selection of good models. Consequently, we applied a more stringent method for 3D model evaluation this time. All generated models were evaluated based on a structural quality score calculated by both Verify3D and Prosa2003 programs. It turns out that the prediction results of our human group were supported by the results of three servers. The pipeline that we constructed for our human group prediction and human intervention were also greatly effective in improving prediction models, but the efficacy of our scheme for 3D model evaluation was obscure.

Large-scale analysis of human alternative protein isoforms: pattern classification and correlation with subcellular localization signals.

We investigated human alternative protein isoforms of >2600 genes based on full-length cDNA clones and SwissProt. We classified the isoforms and examined their co-occurrence for each gene. Further, we investigated potential relationships between these changes and differential subcellular localization. The two most abundant patterns were the one with different C-terminal regions and the one with an internal insertion, which together account for 43% of the total. Although changes of the N-terminal region are less common than those of the C-terminal region, extension of the C-terminal region is much less common than that of the N-terminal region, probably because of the difficulty of removing stop codons in one isoform. We also found that there are some frequently used combinations of co-occurrence in alternative isoforms. We interpret this as evidence that there is some structural relationship which produces a repertoire of isoformal patterns. Finally, many terminal changes are predicted to cause differential subcellular localization, especially in targeting either peroxisomes or mitochondria. Our study sheds new light on the enrichment of the human proteome through alternative splicing and related events. Our database of alternative protein isoforms is available through the internet.

Molecular dynamics simulation of dimeric and monomeric forms of human prion protein: insight into dynamics and properties.

A central theme in prion protein research is the detection of the process that underlies the conformational transition from the normal cellular prion form (PrP(C)) to its pathogenic isoform (PrP(Sc)). Although the three-dimensional structures of monomeric and dimeric human prion protein (HuPrP) have been revealed by NMR spectroscopy and x-ray crystallography, the process underlying the conformational change from PrP(C) to PrP(Sc) and the dynamics and functions of PrP(C) remain unknown. The dimeric form is thought to play an important role in the conformational transition. In this study, we performed molecular dynamics (MD) simulations on monomeric and dimeric HuPrP at 300 K and 500 K for 10 ns to investigate the differences in the properties of the monomer and the dimer from the perspective of dynamic and structural behaviors. Simulations were also undertaken with Asp178Asn and acidic pH, which is known as a disease-associated factor. Our results indicate that the dynamics of the dimer and monomer were similar (e.g., denaturation of helices and elongation of the beta-sheet). However, additional secondary structure elements formed in the dimer might result in showing the differences in dynamics and properties between the monomer and dimer (e.g., the greater retention of dimeric than monomeric tertiary structure).