PSI-BLAST-ISS: an intermediate sequence search tool for estimation of the position-specific alignment reliability

BACKGROUND: Protein sequence alignments have become indispensable for virtually any evolutionary, structural or functional study involving proteins. Modern sequence search and comparison methods combined with rapidly increasing sequence data often can reliably match even distantly related proteins that share little sequence similarity. However, even highly significant matches generally may have incorrectly aligned regions. Therefore when exact residue correspondence is used to transfer biological information from one aligned sequence to another, it is critical to know which alignment regions are reliable and which may contain alignment errors. RESULTS: PSI-BLAST-ISS is a standalone Unix-based tool designed to delineate reliable regions of sequence alignments as well as to suggest potential variants in unreliable regions. The region-specific reliability is assessed by producing multiple sequence alignments in different sequence contexts followed by the analysis of the consistency of alignment variants. The PSI-BLAST-ISS output enables the user to simultaneously analyze alignment reliability between query and multiple homologous sequences. In addition, PSI-BLAST-ISS can be used to detect distantly related homologous proteins. The software is freely available at: . CONCLUSION: PSI-BLAST-ISS is an effective reliability assessment tool that can be useful in applications such as comparative modelling or analysis of individual sequence regions. It favorably compares with the existing similar software both in the performance and functional features

Human HLTF mediates postreplication repair by its HIRAN domain-dependent replication fork remodelling

Defects in the ability to respond properly to an unrepaired DNA lesion blocking replication promote genomic instability and cancer. Human HLTF, implicated in error-free replication of damaged DNA and tumour suppression, exhibits a HIRAN domain, a RING domain, and a SWI/SNF domain facilitating DNA-binding, PCNA-polyubiquitin-ligase, and dsDNA-translocase activities, respectively. Here, we investigate the mechanism of HLTF action with emphasis on its HIRAN domain. We found that in cells HLTF promotes the filling-in of gaps left opposite damaged DNA during replication, and this postreplication repair function depends on its HIRAN domain. Our biochemical assays show that HIRAN domain mutant HLTF proteins retain their ubiquitin ligase, ATPase and dsDNA translocase activities but are impaired in binding to a model replication fork. These data and our structural study indicate that the HIRAN domain recruits HLTF to a stalled replication fork, and it also provides the direction for the movement of the dsDNA translocase motor domain for fork reversal. In more general terms, we suggest functional similarities between the HIRAN, the OB, the HARP2, and other domains found in certain motor proteins, which may explain why only a subset of DNA translocases can carry out fork reversal

Identification of new homologs of PD-(D/E)XK nucleases by support vector machines trained on data derived from profile–profile alignments

PD-(D/E)XK nucleases, initially represented by only Type II restriction enzymes, now comprise a large and extremely diverse superfamily of proteins. They participate in many different nucleic acids transactions including DNA degradation, recombination, repair and RNA processing. Different PD-(D/E)XK families, although sharing a structurally conserved core, typically display little or no detectable sequence similarity except for the active site motifs. This makes the identification of new superfamily members using standard homology search techniques challenging. To tackle this problem, we developed a method for the detection of PD-(D/E)XK families based on the binary classification of profile–profile alignments using support vector machines (SVMs). Using a number of both superfamily-specific and general features, SVMs were trained to identify true positive alignments of PD-(D/E)XK representatives. With this method we identified several PFAM families of uncharacterized proteins as putative new members of the PD-(D/E)XK superfamily. In addition, we assigned several unclassified restriction enzymes to the PD-(D/E)XK type. Results show that the new method is able to make confident assignments even for alignments that have statistically insignificant scores. We also implemented the method as a freely accessible web server at http://www.ibt.lt/bioinformatics/software/pdexk/

Programmable RNA Shredding by the Type III-A CRISPR-Cas System of Streptococcus thermophilus

Immunity against viruses and plasmids provided by CRISPR-Cas systems relies on a ribonucleoprotein effector complex that triggers the degradation of invasive nucleic acids (NA). Effector complexes of type I (Cascade) and II (Cas9-dual RNA) target foreign DNA. Intriguingly, the genetic evidence suggests that the type III-A Csm complex targets DNA, whereas biochemical data show that the type III-B Cmr complex cleaves RNA. Here we aimed to investigate NA specificity and mechanism of CRISPR interference for the Streptococcus thermophilus Csm (III-A) complex (StCsm). When expressed in Escherichia coli, two complexes of different stoichiometry copurified with 40 and 72 nt crRNA species, respectively. Both complexes targeted RNA and generated multiple cuts at 6 nt intervals. The Csm3 protein, present in multiple copies in both Csm complexes, acts as endoribonuclease. In the heterologous E. coli host, StCsm restricts MS2 RNA phage in a Csm3 nuclease-dependent manner. Thus, our results demonstrate that the type III-A StCsm complex guided by crRNA targets RNA and not DNA. Highlights • Streptococcus thermophilus type III-A Csm (StCsm) complex targets RNA •Multiple cuts are introduced in the target RNA at 6 nt intervals •Csm3 protein subunits are responsible for endoribonuclease activity of the complex •StCsm complex offers a programmable tool for RNA degradatio

Restriction endonuclease BpuJI specific for the 5′-CCCGT sequence is related to the archaeal Holliday junction resolvase family

Type IIS restriction endonucleases (REases) recognize asymmetric DNA sequences and cleave both DNA strands at fixed positions downstream of the recognition site. REase BpuJI recognizes the asymmetric sequence 5′-CCCGT, however it cuts at multiple sites in the vicinity of the target sequence. We show that BpuJI is a dimer, which has two DNA binding surfaces and displays optimal catalytic activity when bound to two recognition sites. BpuJI is cleaved by chymotrypsin into an N-terminal domain (NTD), which lacks catalytic activity but binds specifically to the recognition sequence as a monomer, and a C-terminal domain (CTD), which forms a dimer with non-specific nuclease activity. Fold recognition approach reveals that the CTD of BpuJI is structurally related to archaeal Holliday junction resolvases (AHJR). We demonstrate that the isolated catalytic CTD of BpuJI possesses end-directed nuclease activity and preferentially cuts 3 nt from the 3′-terminus of blunt-ended DNA. The nuclease activity of the CTD is repressed in the apo-enzyme and becomes activated upon specific DNA binding by the NTDs. This leads to a complicated pattern of specific DNA cleavage in the vicinity of the target site. Bioinformatics analysis identifies the AHJR-like domain in the putative Type III enzymes and functionally uncharacterized proteins

Detection of distant evolutionary relationships between protein families using theory of sequence profile-profile comparison

<p>Abstract</p> <p>Background</p> <p>Detection of common evolutionary origin (homology) is a primary means of inferring protein structure and function. At present, comparison of protein families represented as sequence profiles is arguably the most effective homology detection strategy. However, finding the best way to represent evolutionary information of a protein sequence family in the profile, to compare profiles and to estimate the biological significance of such comparisons, remains an active area of research.</p> <p>Results</p> <p>Here, we present a new homology detection method based on sequence profile-profile comparison. The method has a number of new features including position-dependent gap penalties and a global score system. Position-dependent gap penalties provide a more biologically relevant way to represent and align protein families as sequence profiles. The global score system enables an analytical solution of the statistical parameters needed to estimate the statistical significance of profile-profile similarities. The new method, together with other state-of-the-art profile-based methods (HHsearch, COMPASS and PSI-BLAST), is benchmarked in all-against-all comparison of a challenging set of SCOP domains that share at most 20% sequence identity. For benchmarking, we use a reference ("gold standard") free model-based evaluation framework. Evaluation results show that at the level of protein domains our method compares favorably to all other tested methods. We also provide examples of the new method outperforming structure-based similarity detection and alignment. The implementation of the new method both as a standalone software package and as a web server is available at <url>http://www.ibt.lt/bioinformatics/coma</url>.</p> <p>Conclusion</p> <p>Due to a number of developments, the new profile-profile comparison method shows an improved ability to match distantly related protein domains. Therefore, the method should be useful for annotation and homology modeling of uncharacterized proteins.</p

Structural and functional analysis of rare missense mutations in human chorionic gonadotrophin β-subunit

Heterodimeric hCG is one of the key hormones determining early pregnancy success. We have previously identified rare missense mutations in hCGβ genes with potential pathophysiological importance. The present study assessed the impact of these mutations on the structure and function of hCG by applying a combination of in silico (sequence and structure analysis, molecular dynamics) and in vitro (co-immunoprecipitation, immuno- and bioassays) approaches. The carrier status of each mutation was determined for 1086 North-Europeans [655 patients with recurrent miscarriage (RM)/431 healthy controls from Estonia, Finland and Denmark] using PCR-restriction fragment length polymorphism. The mutation CGB5 p.Val56Leu (rs72556325) was identified in a single heterozygous RM patient and caused a structural hindrance in the formation of the hCGα/β dimer. Although the amount of the mutant hCGβ assembled into secreted intact hCG was only 10% compared with the wild-type, a stronger signaling response was triggered upon binding to its receptor, thus compensating the effect of poor dimerization. The mutation CGB8 p.Pro73Arg (rs72556345) was found in five heterozygotes (three RM cases and two control individuals) and was inherited by two of seven studied live born children. The mutation caused ∼50% of secreted β-subunits to acquire an alternative conformation, but did not affect its biological activity. For the CGB8 p.Arg8Trp (rs72556341) substitution, the applied in vitro methods revealed no alterations in the assembly of intact hCG as also supported by an in silico analysis. In summary, the accumulated data indicate that only mutations with neutral or mild functional consequences might be tolerated in the major hCGβ genes CGB5 and CGB8

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70–75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70–80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.Cancer Research UK, Grant/Award Number: FC001003; Changzhou Science and Technology Bureau, Grant/Award Number: CE20200503; Department of Energy and Climate Change, Grant/Award Numbers: DE-AR001213, DE-SC0020400, DE-SC0021303; H2020 European Institute of Innovation and Technology, Grant/Award Numbers: 675728, 777536, 823830; Institut national de recherche en informatique et en automatique (INRIA), Grant/Award Number: Cordi-S; Lietuvos Mokslo Taryba, Grant/Award Numbers: S-MIP-17-60, S-MIP-21-35; Medical Research Council, Grant/Award Number: FC001003; Japan Society for the Promotion of Science KAKENHI, Grant/Award Number: JP19J00950; Ministerio de Ciencia e Innovación, Grant/Award Number: PID2019-110167RB-I00; Narodowe Centrum Nauki, Grant/Award Numbers: UMO-2017/25/B/ST4/01026, UMO-2017/26/M/ST4/00044, UMO-2017/27/B/ST4/00926; National Institute of General Medical Sciences, Grant/Award Numbers: R21GM127952, R35GM118078, RM1135136, T32GM132024; National Institutes of Health, Grant/Award Numbers: R01GM074255, R01GM078221, R01GM093123, R01GM109980, R01GM133840, R01GN123055, R01HL142301, R35GM124952, R35GM136409; National Natural Science Foundation of China, Grant/Award Number: 81603152; National Science Foundation, Grant/Award Numbers: AF1645512, CCF1943008, CMMI1825941, DBI1759277, DBI1759934, DBI1917263, DBI20036350, IIS1763246, MCB1925643; NWO, Grant/Award Number: TOP-PUNT 718.015.001; Wellcome Trust, Grant/Award Number: FC00100