JPred4: a protein secondary structure prediction server

Purriato is a conserved small open reading frame gene that interacts with the CASA pathway to regulate muscle homeostasis and epithelial tissue growth in Drosophila

Recent advances in proteogenomic techniques and bioinformatic pipelines have permitted the detection of thousands of translated small Open Reading Frames (smORFs), which contain less than 100 codons, in eukaryotic genomes. Hundreds of these actively translated smORFs display conserved sequence, structure and evolutionary signatures indicating that the translated peptides could fulfil important biological roles. Despite their abundance, only tens of smORF genes have been fully characterised; these act mainly as regulators of canonical proteins involved in essential cellular processes. Importantly, some of these smORFs display conserved functions with their mutations being associated with pathogenesis. Thus, investigating smORF roles in Drosophila will not only expand our understanding of their functions but it may have an impact in human health. Here we describe the function of a novel and essential Drosophila smORF gene named purriato (prto). prto belongs to an ancient gene family whose members have expanded throughout the Protostomia clade. prto encodes a transmembrane peptide which is localized in endo-lysosomes and perinuclear and plasma membranes. prto is dynamically expressed in mesodermal tissues and imaginal discs. Targeted prto knockdown (KD) in these organs results in changes in nuclear morphology and endo-lysosomal distributions correlating with the loss of sarcomeric homeostasis in muscles and reduction of mitosis in wing discs. Consequently, prto KD mutants display severe reduction of motility, and shorter wings. Finally, our genetic interaction experiments show that prto function is closely associated to the CASA pathway, a conserved mechanism involved in turnover of mis-folded proteins and linked to muscle dystrophies and neurodegenerative diseases. Thus, this study shows the relevance of smORFs in regulating important cellular functions and supports the systematic characterisation of this class of genes to understand their functions and evolution.

Intracytoplasmic-membrane development in alphaproteobacteria involves the homolog of the mitochondrial crista-developing protein Mic60

Mitochondrial cristae expand the surface area of respiratory membranes and ultimately allow for the evolutionary scaling of respiration with cell volume across eukaryotes. The discovery of Mic60 homologs among alphaproteobacteria, the closest extant relatives of mitochondria, suggested that cristae might have evolved from bacterial intracytoplasmic membranes (ICMs). Here, we investigated the predicted structure and function of alphaproteobacterial Mic60, and a protein encoded by an adjacent gene Orf52, in two distantly related purple alphaproteobacteria, Rhodobacter sphaeroides and Rhodopseudomonas palustris. In addition, we assessed the potential physical interactors of Mic60 and Orf52 in R. sphaeroides. We show that the three α helices of mitochondrial Mic60's mitofilin domain, as well as its adjacent membrane-binding amphipathic helix, are present in alphaproteobacterial Mic60. The disruption of Mic60 and Orf52 caused photoheterotrophic growth defects, which are most severe under low light conditions, and both their disruption and overexpression led to enlarged ICMs in both studied alphaproteobacteria. We also found that alphaproteobacterial Mic60 physically interacts with BamA, the homolog of Sam50, one of the main physical interactors of eukaryotic Mic60. This interaction, responsible for making contact sites at mitochondrial envelopes, has been conserved in modern alphaproteobacteria despite more than a billion years of evolutionary divergence. Our results suggest a role for Mic60 in photosynthetic ICM development and contact site formation at alphaproteobacterial envelopes. Overall, we provide support for the hypothesis that mitochondrial cristae evolved from alphaproteobacterial ICMs and have therefore improved our understanding of the nature of the mitochondrial ancestor.

Machine learning for the identification of respiratory viral attachment machinery from sequences data

At the outset of an emergent viral respiratory pandemic, sequence data is among the first molecular information available. As viral attachment machinery is a key target for therapeutic and prophylactic interventions, rapid identification of viral "spike" proteins from sequence can significantly accelerate the development of medical countermeasures. For six families of respiratory viruses, covering the vast majority of airborne and droplet-transmitted diseases, host cell entry is mediated by the binding of viral surface glycoproteins that interact with a host cell receptor. In this report it is shown that sequence data for an unknown virus belonging to one of the six families above provides sufficient information to identify the protein(s) responsible for viral attachment. Random forest models that take as input a set of respiratory viral sequences can classify the protein as "spike" vs. non-spike based on predicted secondary structure elements alone (with 97.3% correctly classified) or in combination with N-glycosylation related features (with 97.0% correctly classified). Models were validated through 10-fold cross-validation, bootstrapping on a class-balanced set, and an out-of-sample extra-familial validation set. Surprisingly, we showed that secondary structural elements and N-glycosylation features were sufficient for model generation. The ability to rapidly identify viral attachment machinery directly from sequence data holds the potential to accelerate the design of medical countermeasures for future pandemics. Furthermore, this approach may be extendable for the identification of other potential viral targets and for viral sequence annotation in general in the future.

IGPRED-MultiTask: A Deep Learning Model to Predict Protein Secondary Structure, Torsion Angles and Solvent Accessibility

Protein secondary structure, solvent accessibility and torsion angle predictions are preliminary steps to predict 3D structure of a protein. Deep learning approaches have achieved significant improvements in predicting various features of protein structure. In this study, IGPRED-Multitask, a deep learning model with multi task learning architecture based on deep inception network, graph convolutional network and a bidirectional long short-term memory is proposed. Moreover, hyper-parameters of the model are fine-tuned using Bayesian optimization, which is faster and more effective than grid search. The same benchmark test data sets as in the OPUS-TASS paper including TEST2016, TEST2018, CASP12, CASP13, CASPFM, HARD68, CAMEO93, CAMEO93_HARD, as well as the train and validation sets, are used for fair comparison with the literature. Statistically significant improvements are observed in secondary structure prediction on 4 datasets, in phi angle prediction on 2 datasets and in psi angel prediction on 3 datasets compared to the state-of-the-art methods. For solvent accessibility prediction, TEST2016 and TEST2018 datasets are used only to assess the performance of the proposed model.

A tetramerization domain in prokaryotic and eukaryotic transcription regulators homologous to p53

Transcriptional regulation usually requires the action of several proteins that either repress or activate a promotor of an open reading frame. These proteins can counteract each other, thus allowing tight regulation of the transcription of the corresponding genes, where tight repression is often linked to DNA looping or cross-linking. Here, the tetramerization domain of the bacterial gene repressor Rco from Bacillus subtilis plasmid pLS20 (RcopLS20) has been identified and its structure is shown to share high similarity to the tetramerization domain of the well known p53 family of human tumor suppressors, despite lacking clear sequence homology. In RcopLS20, this tetramerization domain is responsible for inducing DNA looping, a process that involves multiple tetramers. In accordance, it is shown that RcopLS20 can form octamers. This domain was named TetDloop and its occurrence was identified in other Bacillus species. The TetDloop fold was also found in the structure of a transcriptional repressor from Salmonella phage SPC32H. It is proposed that the TetDloop fold has evolved through divergent evolution and that the TetDloop originates from a common ancestor predating the occurrence of multicellular life.

Insertion sequences related to ISAjo2 target pdif and dif sites and belong to a new IS family, the IS1202 family

Several insertion sequences (IS) found in various Acinetobacter species exhibit target specificity. They are found, in the same orientation, 5 bp from the XerC binding site of the pdif sites associated with dif modules in Acinetobacter plasmids, and searches revealed they are also found near chromosomal dif sites of Acinetobacter species. These IS are 1.5 kb long, bounded by 24-26 bp imperfect terminal inverted repeats (TIRs) and encode a large transposase of 441-457 aa. They generate 5 bp target site duplications (TSDs). Structural predictions of the ISAjo2 transposase, TnpAjo2, modelled on TnsB of Tn7 revealed two N-terminal HTH domains followed by an RNaseH fold (DDE domain), a β barrel and a C-terminal domain. Similar to Tn7, the outer IS ends are 5'-TGT and ACA-3', and an additional Tnp binding site, corresponding to the internal portion of the IR, is found near each end. However, the Acinetobacter IS do not encode further proteins related to those required by Tn7 for targeted transposition, and the transposase may interact directly with XerC bound to a dif-like site. We propose that these IS, currently in the IS1202 group in the not characterized yet (NCY) category in ISFinder, are part of a distinct IS1202 family. Other IS listed as in the IS1202 group encode transposases related to TnpAjo2 (25-56 % amino acid identity) and have similar TIRs but fall into three groups based on the TSD length (3-5, >15, 0 bp). Those with 3-5 bp TSDs may also target dif-like sites but targets were not found for the other groups.

New biochemistry in the Rhodanese-phosphatase superfamily: emerging roles in diverse metabolic processes, nucleic acid modifications, and biological conflicts

The protein-tyrosine/dual-specificity phosphatases and rhodanese domains constitute a sprawling superfamily of Rossmannoid domains that use a conserved active site with a cysteine to catalyze a range of phosphate-transfer, thiotransfer, selenotransfer and redox activities. While these enzymes have been extensively studied in the context of protein/lipid head group dephosphorylation and various thiotransfer reactions, their overall diversity and catalytic potential remain poorly understood. Using comparative genomics and sequence/structure analysis, we comprehensively investigate and develop a natural classification for this superfamily. As a result, we identified several novel clades, both those which retain the catalytic cysteine and those where a distinct active site has emerged in the same location (e.g. diphthine synthase-like methylases and RNA 2' OH ribosyl phosphate transferases). We also present evidence that the superfamily has a wider range of catalytic capabilities than previously known, including a set of parallel activities operating on various sugar/sugar alcohol groups in the context of NAD+-derivatives and RNA termini, and potential phosphate transfer activities involving sugars and nucleotides. We show that such activities are particularly expanded in the RapZ-C-DUF488-DUF4326 clade, defined here for the first time. Some enzymes from this clade are predicted to catalyze novel DNA-end processing activities as part of nucleic-acid-modifying systems that are likely to function in biological conflicts between viruses and their hosts.

Distinguishing features of fold-switching proteins

Though many folded proteins assume one stable structure that performs one function, a small-but-increasing number remodel their secondary and tertiary structures and change their functions in response to cellular stimuli. These fold-switching proteins regulate biological processes and are associated with autoimmune dysfunction, severe acute respiratory syndrome coronavirus-2 infection, and more. Despite their biological importance, it is difficult to computationally predict fold switching. With the aim of advancing computational prediction and experimental characterization of fold switchers, this review discusses several features that distinguish fold-switching proteins from their single-fold and intrinsically disordered counterparts. First, the isolated structures of fold switchers are less stable and more heterogeneous than single folders but more stable and less heterogeneous than intrinsically disordered proteins (IDPs). Second, the sequences of single fold, fold switching, and intrinsically disordered proteins can evolve at distinct rates. Third, proteins from these three classes are best predicted using different computational techniques. Finally, late-breaking results suggest that single folders, fold switchers, and IDPs have distinct patterns of residue-residue coevolution. The review closes by discussing high-throughput and medium-throughput experimental approaches that might be used to identify new fold-switching proteins.

Computational prediction of extracellular loops of the Por39 outer membrane porin of Rhodospirillum rubrum suitable for epitope surface display

Outer membrane porins from Gram-negative bacteria are established vehicles for the production of vaccines. Typically, one or more of the extracellular loops of a porin are replaced by a peptide encoding a foreign epitope, and recombinant porin is then used as a vaccine. However, many host strains are potentially pathogenic, and also produce toxic lipopolysaccharide (LPS), both of which are undesirable for safety reasons. In contrast, the outer membrane porins from photosynthetic, purple bacteria have no known human pathology and produce only weakly toxic LPS. The purple bacterium Rhodospirillum rubrum is well-suited for large-scale biotechnology, and expresses a major porin, Por39, which is a candidate for a vaccine platform. Unfortunately, the atomic structure of Por39 could not be determined so far, and Por39 shows only a weak homology to other porins of known structure, making the assignment of external loops difficult. Here, we construct a knowledge-based model of Por39 using secondary structure constraints from both the low sequence homology to the 2POR porin from Rhodobacter capsulatus, for which the X-ray structure is known, as well as those obtained using secondary structure prediction packages. The secondary structure predictions were used to constrain a three-dimensional model created using the I-TASSER package. The modelling procedure was validated by predicting the structure of 2POR using the same strategy, but excluding the 2POR X-ray structure from the I-TASSER database. The final Por39 model allows three external loops to be defined precisely, and could also be used to obtain an initial model for the closely related Por41 using molecular modelling. These structures provide a good starting point for the insertion of epitopes with vaccine potential.

Metagenome diversity illuminates origins of pathogen effectors

Recent metagenome assembled genome (MAG) analyses have profoundly impacted Rickettsiology systematics. Discovery of basal lineages (Mitibacteraceae and Athabascaceae) with predicted extracellular lifestyles reveals an evolutionary timepoint for the transition to host dependency, which occurred independent of mitochondrial evolution. Notably, these basal rickettsiae carry the Rickettsiales vir homolog (rvh) type IV secretion system (T4SS) and purportedly use rvh to kill congener microbes rather than parasitize host cells as described for derived rickettsial pathogens. MAG analysis also substantially increased diversity for genus Rickettsia and delineated a basal lineage (Tisiphia) that stands to inform on the rise of human pathogens from protist and invertebrate endosymbionts. Herein, we probed Rickettsiales MAG and genomic diversity for the distribution of Rickettsia rvh effectors to ascertain their origins. A sparse distribution of most Rickettsia rvh effectors outside of Rickettsiaceae lineages indicates unique rvh evolution from basal extracellular species and other rickettsial families. Remarkably, nearly every effector was found in multiple divergent forms with variable architectures, illuminating profound roles for gene duplication and recombination in shaping effector repertoires in Rickettsia pathogens. Lateral gene transfer plays a prominent role shaping the rvh effector landscape, as evinced by the discover of many effectors on plasmids and conjugative transposons, as well as pervasive effector gene exchange between Rickettsia and Legionella species. Our study exemplifies how MAGs can provide incredible insight on the origins of pathogen effectors and how their architectural modifications become tailored to eukaryotic host cell biology.

Genome-wide screening of the RNase T2 gene family and functional analyses in jujube (Ziziphus jujuba Mill.)

BACKGROUND

Ribonuclease (RNase T2) plays crucial roles in plant evolution and breeding. However, there have been few studies on the RNase T2 gene family in Ziziphus jujuba Mill., one of important dried fruit tree species. Recently, the released sequences of the reference genome of jujube provide a good chance to perform genome-wide identification and characterization of ZjRNase gene family in the jujube.

RESULTS

In this study, we identified four members of RNase T2 in jujube distributed on three chromosomes and unassembled chromosomes. They all contained two conserved sites (CASI and CASII). Analysis of the phylogenetic relationships revealed that the RNase T2 genes in jujube could be divided into two groups: ZjRNase1 and ZjRNase2 belonged to class I, while ZjRNase3 and ZjRNase4 belonged to class II. Only ZjRNase1 and ZjRNase2 expression were shown by the jujube fruit transcriptome analysis. So ZjRNase1 and ZjRNase2 were selected functional verification by overexpression transformation of Arabidopsis. The overexpression of these two genes led to an approximately 50% reduction in seed number, which deserve further attention. Moreover, the leaves of the ZjRNase1 overexpression transgenic lines were curled and twisted. Overexpression of ZjRNase2 resulted in shortened and crisp siliques and the production of trichomes, and no seeds were produced.

CONCLUSION

In summary, these findings will provide new insights into the molecular mechanisms of low number of hybrid seeds in jujube and a reference for the future molecular breeding of jujube.

WG-ICRN: Protein 8-state secondary structure prediction based on Wasserstein generative adversarial networks and residual networks with Inception modules

Protein secondary structure is the basis of studying the tertiary structure of proteins, drug design and development, and the 8-state protein secondary structure can provide more adequate protein information than the 3-state structure. Therefore, this paper proposes a novel method WG-ICRN for predicting protein 8-state secondary structures. First, we use the Wasserstein generative adversarial network (WGAN) to extract protein features in the position-specific scoring matrix (PSSM). The extracted features are combined with PSSM into a new feature set of WG-data, which contains richer feature information. Then, we use the residual network (ICRN) with Inception to further extract the features in WG-data and complete the prediction. Compared with the residual network, ICRN can reduce parameter calculations and increase the width of feature extraction to obtain more feature information. We evaluated the prediction performance of the model using six datasets. The experimental results show that the WGAN has excellent feature extraction capabilities, and ICRN can further improve network performance and improve prediction accuracy. Compared with four popular models, WG-ICRN achieves better prediction performance.

Ca2+-Dependent and -Independent Calmodulin Binding to the Cytoplasmic Loop of Gap Junction Connexins

Ca²⁺/calmodulin (Ca²⁺/CaM) interaction with connexins (Cx) is well-established; however, the mechanistic basis of regulation of gap junction function by Ca²⁺/CaM is not fully understood. Ca²⁺/CaM is predicted to bind to a domain in the C-terminal portion of the intracellular loop (CL2) in the vast majority of Cx isoforms and for a number of Cx-s this prediction has proved correct. In this study, we investigate and characterise both Ca²⁺/CaM and apo-CaM binding to selected representatives of each of the α, β and γ connexin family to develop a better mechanistic understanding of CaM effects on gap junction function. The affinity and kinetics Ca²⁺/CaM and apo-CaM interactions of CL2 peptides of β-Cx32, γ-Cx35, α-Cx43, α-Cx45 and α-Cx57 were investigated. All five Cx CL2 peptides were found to have high affinity for Ca²⁺/CaM with dissociation constants (K_d(+Ca)) from 20 to 150 nM. The limiting rate of binding and the rates of dissociation covered a broad range. In addition, we obtained evidence for high affinity Ca²⁺-independent interaction of all five peptides with CaM, consistent with CaM remaining anchored to gap junctions in resting cells. However, for the α-Cx45 and α-Cx57 CL2 peptides, Ca²⁺-dependent association at resting [Ca²⁺] of 50-100 nM is indicated in these complexes as one of the CaM Ca²⁺ binding sites displays high affinity with K_d of 70 and 30 nM for Ca²⁺, respectively. Furthermore, complex conformational changes were observed in peptide-apo-CaM complexes with the structure of CaM compacted or stretched by the peptide in a concentration dependent manner suggesting that the CL2 domain may undergo helix-to-coil transition and/or forms bundles, which may be relevant in the hexameric gap junction. We demonstrate inhibition of gap junction permeability by Ca²⁺/CaM in a dose dependent manner, further cementing Ca²⁺/CaM as a regulator of gap junction function. The motion of a stretched CaM-CL2 complex compacting upon Ca²⁺ binding may bring about the Ca²⁺/CaM block of the gap junction pore by a push and pull action on the CL2 C-terminal hydrophobic residues of transmembrane domain 3 (TM3) in and out of the membrane.

The deficiency of poly-β-1,6-N-acetyl-glucosamine deacetylase trigger A. baumannii to convert to biofilm-independent colistin-tolerant cells

Acinetobacter baumannii is a nosocomial pathogen that can be resistant to antibiotics by rapidly modulating its anti-drug mechanisms. The multidrug-resistant A. baumannii has been considered one of the most threatening pathogens to our society. Biofilm formation and persistent cells within the biofilm matrix are recognized as intractable problems, especially in hospital-acquired infections. Poly-β-1,6-N-acetyl-glucosamine (PNAG) is one of the important building blocks in A. baumannii's biofilm. Here, we discover a protein phosphoryl-regulation on PNAG deacetylase, AbPgaB1, in which residue Ser411 was phosphorylated. The phosphoryl-regulation on AbPgaB1 modulates the product turnover rate in which deacetylated PNAG is produced and reflected in biofilm production. We further uncovered the PgaB deficient A. baumannii strain shows the lowest level of biofilm production but has a high minimal inhibition concentration to antibiotic colistin and tetracycline. Based on bactericidal post-antibiotic effects and time-dependent killing assays with antibacterial drugs, we claim that the PgaB-deficient A. baumannii converts to colistin-tolerant cells. This study utilizes a biofilm-independent colistin-tolerant model of A. baumannii to further investigate its characteristics and mechanisms to better understand clinical outcomes.

Uphill energy transfer mechanism for photosynthesis in an Antarctic alga

Prasiola crispa, an aerial green alga, forms layered colonies under the severe terrestrial conditions of Antarctica. Since only far-red light is available at a deep layer of the colony, P. crispa has evolved a molecular system for photosystem II (PSII) excitation using far-red light with uphill energy transfer. However, the molecular basis underlying this system remains elusive. Here, we purified a light-harvesting chlorophyll (Chl)-binding protein complex from P. crispa (Pc-frLHC) that excites PSII with far-red light and revealed its ring-shaped structure with undecameric 11-fold symmetry at 3.13 Å resolution. The primary structure suggests that Pc-frLHC evolved from LHCI rather than LHCII. The circular arrangement of the Pc-frLHC subunits is unique among eukaryote LHCs and forms unprecedented Chl pentamers at every subunit‒subunit interface near the excitation energy exit sites. The Chl pentamers probably contribute to far-red light absorption. Pc-frLHC's unique Chl arrangement likely promotes PSII excitation with entropy-driven uphill excitation energy transfer.

Ensemble deep learning models for protein secondary structure prediction using bidirectional temporal convolution and bidirectional long short-term memory

Protein secondary structure prediction (PSSP) is a challenging task in computational biology. However, existing models with deep architectures are not sufficient and comprehensive for deep long-range feature extraction of long sequences. This paper proposes a novel deep learning model to improve Protein secondary structure prediction. In the model, our proposed bidirectional temporal convolutional network (BTCN) can extract the bidirectional deep local dependencies in protein sequences segmented by the sliding window technique, the bidirectional long short-term memory (BLSTM) network can extract the global interactions between residues, and our proposed multi-scale bidirectional temporal convolutional network (MSBTCN) can further capture the bidirectional multi-scale long-range features of residues while preserving the hidden layer information more comprehensively. In particular, we also propose that fusing the features of 3-state and 8-state Protein secondary structure prediction can further improve the prediction accuracy. Moreover, we also propose and compare multiple novel deep models by combining bidirectional long short-term memory with temporal convolutional network (TCN), reverse temporal convolutional network (RTCN), multi-scale temporal convolutional network (multi-scale bidirectional temporal convolutional network), bidirectional temporal convolutional network and multi-scale bidirectional temporal convolutional network, respectively. Furthermore, we demonstrate that the reverse prediction of secondary structure outperforms the forward prediction, suggesting that amino acids at later positions have a greater impact on secondary structure recognition. Experimental results on benchmark datasets including CASP10, CASP11, CASP12, CASP13, CASP14, and CB513 show that our methods achieve better prediction performance compared to five state-of-the-art methods.

Screening and identification of potential MERS-CoV papain-like protease (PLpro) inhibitors; Steady-state kinetic and Molecular dynamic studies

MERS-CoV belongs to the coronavirus group. Recent years have seen a rash of coronavirus epidemics. In June 2012, MERS-CoV was discovered in the Kingdom of Saudi Arabia, with 2,591 MERSA cases confirmed by lab tests by the end of August 2022 and 894 deaths at a case-fatality ratio (CFR) of 34.5% documented worldwide. Saudi Arabia reported the majority of these cases, with 2,184 cases and 813 deaths (CFR: 37.2%), necessitating a thorough understanding of the molecular machinery of MERS-CoV. To develop antiviral medicines, illustrative investigation of the protein in coronavirus subunits are required to increase our understanding of the subject. In this study, recombinant expression and purification of MERS-CoV (PLpro), a primary goal for the development of 22 new inhibitors, were completed using a high throughput screening methodology that employed fragment-based libraries in conjunction with structure-based virtual screening. Compounds 2, 7, and 20, showed significant biological activity. Moreover, a docking analysis revealed that the three compounds had favorable binding mood and binding free energy. Molecular dynamic simulation demonstrated the stability of compound 2 (2-((Benzimidazol-2-yl) thio)-1-arylethan-1-ones) the strongest inhibitory activity against the PLpro enzyme. In addition, disubstitutions at the meta and para locations are the only substitutions that may boost the inhibitory action against PLpro. Compound 2 was chosen as a MERS-CoV PLpro inhibitor after passing absorption, distribution, metabolism, and excretion studies; however, further investigations are required.

Systems and strategies for plant protein expression

At least a quarter of the protein-encoding genes in plant genomes are predicted to encode enzymes for which no physiological function is known. Determining functions for these uncharacterized enzymes is key to understanding plant metabolism. Functional characterization typically requires expression and purification of recombinant enzymes to be used in enzyme assays and/or for protein structure elucidation studies. Here, we describe several practical considerations used to improve the heterologous expression and purification of Arabidopsis thaliana and Zea mays NAD(P)HX dehydratase (NAXD) and NAD(P)HX epimerase (NAXE), two enzymes that are involved in repair of chemically damaged NAD(P)H cofactors. We provide protocols for transit peptide prediction and construct design, expression in Escherichia coli, and purification of NAXD and NAXE. Many of these strategies are generally applicable to the purification of any plant protein.

Evolution of Rev7 interactions in eukaryotic TLS DNA polymerase Polζ

Translesion synthesis (TLS) DNA polymerase Polζ is crucial for the bypass replication over sites of DNA damage. The Rev7 subunit of Polζ is a HORMA (Hop1, Rev7, Mad2) protein that facilitates recruitment of Polζ to the replication fork via interactions with the catalytic subunit Rev3 and the translesion synthesis scaffold protein Rev1. Human Rev7 (hRev7) interacts with two Rev7-binding motifs (RBMs) of hRev3 by a mechanism conserved among HORMA proteins whereby the safety-belt loop of hRev7 closes on the top of the ligand. The two copies of hRev7 tethered by the two hRev3-RBMs form a symmetric head-to-head dimer through the canonical HORMA dimerization interface. Recent cryo-EM structures reveal that Saccharomyces cerevisiae Polζ (scPolζ) also includes two copies of scRev7 bound to distinct regions of scRev3. Surprisingly, the HORMA dimerization interface is not conserved in scRev7, with the two scRev7 protomers forming an asymmetric head-to-tail dimer with a much smaller interface than the hRev7 dimer. Here, we validated the two adjacent RBM motifs in scRev3, which bind scRev7 with affinities that differ by two orders of magnitude and confirmed the 2:1 stoichiometry of the scRev7:Rev3 complex in solution. However, our biophysical studies reveal that scRev7 does not form dimers in solution either on its own accord or when tethered by the two RBMs in scRev3. These findings imply that the scRev7 dimer observed in the cryo-EM structures is induced by scRev7 interactions with other Polζ subunits and that Rev7 homodimerization via the HORMA interface is a mechanism that emerged later in evolution.

An Interplay between Mitochondrial and ER Targeting of a Bacterial Signal Peptide in Plants

Protein targeting is essential in eukaryotic cells to maintain cell function and organelle identity. Signal peptides are a major type of targeting sequences containing a tripartite structure, which is conserved across all domains in life. They are frequently included in recombinant protein design in plants to increase yields by directing them to the endoplasmic reticulum (ER) or apoplast. The processing of bacterial signal peptides by plant cells is not well understood but could aid in the design of efficient heterologous expression systems. Here we analysed the signal peptide of the enzyme PmoB from methanotrophic bacteria. In plant cells, the PmoB signal peptide targeted proteins to both mitochondria and the ER. This dual localisation was still observed in a mutated version of the signal peptide sequence with enhanced mitochondrial targeting efficiency. Mitochondrial targeting was shown to be dependent on a hydrophobic region involved in transport to the ER. We, therefore, suggest that the dual localisation could be due to an ER-SURF pathway recently characterised in yeast. This work thus sheds light on the processing of bacterial signal peptides by plant cells and proposes a novel pathway for mitochondrial targeting in plants.

FapA is an Intrinsically Disordered Chaperone for Pseudomonas Functional Amyloid FapC

Bacterial functional amyloids contribute to biofilm development by bacteria and provide protection from the immune system and prevent antibiotic treatment. Strategies to target amyloid formation and interrupt biofilm formation have attracted recent interest due to their antimicrobial potential. Functional amyloid in Pseudomonas (Fap) includes FapC as the major component of the fibril while FapB is a minor component suggested to function as a nucleator of FapC. The system also includes the small periplasmic protein FapA, which has been shown to regulate fibril composition and morphology. The interplay between these three components is central in Fap fibril biogenesis. Here we present a comprehensive biophysical and spectroscopy analysis of FapA, FapB and FapC and provide insight into their molecular interactions. We show that all three proteins are primarily disordered with some regions with structural propensities for α-helix and β-sheet. FapA inhibits FapC fibrillation by targeting the nucleation step, whereas for FapB the elongation step is modulated. Furthermore, FapA alters the morphology of FapC (more than FapB) fibrils. Complex formation is observed between FapA and FapC, but not between FapA and FapB, and likely involves the N-terminus of FapA. We conclude that FapA is an intrinsically disordered chaperone for FapC that guards against fibrillation within the periplasm. This new understanding of a natural protective mechanism of Pseudomonas against amyloid formations can serve as inspiration for strategies blocking biofilm formation in infections.

Comparative and expression analyses of AP2/ERF genes reveal copy number expansion and potential functions of ERF genes in Solanaceae

BACKGROUND

The AP2/ERF gene family is a superfamily of transcription factors that are important in the response of plants to abiotic stress and development. However, comprehensive research of the AP2/ERF genes in the Solanaceae family is lacking.

RESULTS

Here, we updated the annotation of AP2/ERF genes in the genomes of eight Solanaceae species, as well as Arabidopsis thaliana and Oryza sativa. We identified 2,195 AP2/ERF genes, of which 368 (17%) were newly identified. Based on phylogenetic analyses, we observed expansion of the copy number of these genes, especially those belonging to specific Ethylene-Responsive Factor (ERF) subgroups of the Solanaceae. From the results of chromosomal location and synteny analyses, we identified that the AP2/ERF genes of the pepper (Capsicum annuum), the tomato (Solanum lycopersicum), and the potato (Solanum tuberosum) belonging to ERF subgroups form a tandem array and most of them are species-specific without orthologs in other species, which has led to differentiation of AP2/ERF gene repertory among Solanaceae. We suggest that these genes mainly emerged through recent gene duplication after the divergence of these species. Transcriptome analyses showed that the genes have a putative function in the response of the pepper and tomato to abiotic stress, especially those in ERF subgroups.

CONCLUSIONS

Our findings will provide comprehensive information on AP2/ERF genes and insights into the structural, evolutionary, and functional understanding of the role of these genes in the Solanaceae.

The molecular mechanism of sialic acid transport mediated by Sialin

Malfunction of the sialic acid transporter caused by various genetic mutations in the SLC17A5 gene encoding Sialin leads to a spectrum of neurodegenerative conditions called free sialic acid storage disorders. Unfortunately, how Sialin transports sialic acid/proton (H+) and how pathogenic mutations impair its function are poorly defined. Here, we present the structure of human Sialin in an inward-facing partially open conformation determined by cryo-electron microscopy, representing the first high-resolution structure of any human SLC17 member. Our analysis reveals two unique features in Sialin: (i) The H+ coupling/sensing requires two highly conserved Glu residues (E171 and E175) instead of one (E175) as implied in previous studies; and (ii) the normal function of Sialin requires the stabilization of a cytosolic helix, which has not been noticed in the literature. By mapping known pathogenic mutations, we provide mechanistic explanations for corresponding functional defects. We propose a structure-based mechanism for sialic acid transport mediated by Sialin.

Investigation of the Structure of Full-Length Tau Proteins with Coarse-Grained and All-Atom Molecular Dynamics Simulations

Tau proteins not only have many important biological functions but also are associated with several neurodegenerative diseases, such as Parkinson's disease and Alzheimer's disease (AD). However, it is still a challenge to identify the atomic structure of full-length tau proteins due to their lengthy and disordered characteristics and the factor that there are no crystal structures of full-length tau proteins available. We performed multi- and large-scale molecular dynamics simulations of the full-length tau monomer (the 2N4R isoform and 441 residues) in aqueous solution under biological conditions with coarse-grained and all-atom force fields. The obtained atomic structures produced radii of gyration and chemical shifts that are in excellent agreement with those of experiment. The generated monomer structure ensemble would be very useful for further studying the oligomerization mechanism and discovering tau oligomerization inhibitors, which are important events in AD drug development.

A novel toxin-antitoxin system swpAB alters gene expression patterns and reduces virulence expression in enterohemorrhagic Escherichia coli

Toxin-antitoxin (TA) systems are found widely among many bacteria, including enterohemorrhagic Escherichia coli (EHEC), but their functions are still poorly understood. In this study, we identified and characterized a novel TA system belonging to the relBE family, classified as a type II TA system, found in EHEC. The protein encoded by the toxin gene is homologous to RelE ribonuclease. Using various conditions for increasing the toxin activity, high-level induction of a toxin gene, and repression of an antitoxin gene in wild-type EHEC, we showed that the TA system, named swpAB (switching of gene expression profile), is involved in selective repression of a set of genes, including some virulence genes, and in the reduction of adherence capacity, rather than in suppression of bacterial growth. A detailed analysis of the profiles of RNA levels along sequences at 15 min after high expression of swpA revealed that two virulence genes, espA and tir, were direct targets of the SwpA toxin. These results suggested that the swpAB system can alter gene expression patterns and change bacterial physiological activity without affecting bacterial growth.

Diversity and prevalence of type VI secretion system effectors in clinical Pseudomonas aeruginosa isolates

Pseudomonas aeruginosa is an opportunistic pathogen and a major driver of morbidity and mortality in people with Cystic Fibrosis (CF). The Type VI secretion system (T6SS) is a molecular nanomachine that translocates effectors across the bacterial membrane into target cells or the extracellular environment enabling intermicrobial interaction. P. aeruginosa encodes three T6SS clusters, the H1-, H2- and H3-T6SS, and numerous orphan islands. Genetic diversity of T6SS-associated effectors in P. aeruginosa has been noted in reference strains but has yet to be explored in clinical isolates. Here, we perform a comprehensive bioinformatic analysis of the pangenome and T6SS effector genes in 52 high-quality clinical P. aeruginosa genomes isolated from CF patients and housed in the Personalised Approach to P. aeruginosa strain repository. We confirm that the clinical CF isolate pangenome is open and principally made up of accessory and unique genes that may provide strain-specific advantages. We observed genetic variability in some effector/immunity encoding genes and show that several well-characterised vgrG and PAAR islands are absent from numerous isolates. Our analysis shows clear evidence of disruption to T6SS genomic loci through transposon, prophage, and mobile genetic element insertions. We identified an orphan vgrG island in P. aeruginosa strain PAK and five clinical isolates using in silico analysis which we denote vgrG7, predicting a gene within this cluster to encode a Tle2 lipase family effector. Close comparison of T6SS loci in clinical isolates compared to reference P. aeruginosa strain PAO1 revealed the presence of genes encoding eight new T6SS effectors with the following putative functions: cytidine deaminase, lipase, metallopeptidase, NADase, and pyocin. Finally, the prevalence of characterised and putative T6SS effectors were assessed in 532 publicly available P. aeruginosa genomes, which suggests the existence of accessory effectors. Our in silico study of the P. aeruginosa T6SS exposes a level of genetic diversity at T6SS genomic loci not seen to date within P. aeruginosa, particularly in CF isolates. As understanding the effector repertoire is key to identifying the targets of T6SSs and its efficacy, this comprehensive analysis provides a path for future experimental characterisation of these mediators of intermicrobial competition and host manipulation.

Theoretical 3D Modeling of NLRP3 Inflammasome Complex

The NOD-like receptor pyrin domain containing 3 (NLRP3) is a multidomain protein that plays a key role in innate immune response. Structures of NLRP3 in different conformational states and bound to cognate partners are available. In this chapter we present an approach to model the oligomeric structure of NLRP3 by homology modeling using multiple templates, symmetry, and refinement. The overall process presented here represents advanced exercise in structural modeling that provides unique insights into the biological role and activation of NLRP3 oligomer. Finally, the same approach can be easily adapted to the rest of the members of the NLRP family.

Pneumococcal capsule expression is controlled through a conserved, distal cis-regulatory element during infection

Streptococcus pneumoniae (the pneumococcus) is the major cause of bacterial pneumonia in the US and worldwide. Studies have shown that the differing chemical make-up between serotypes of its most important virulence factor, the capsule, can dictate disease severity. Here we demonstrate that control of capsule synthesis is also critical for infection and facilitated by two broadly conserved transcription factors, SpxR and CpsR, through a distal cis-regulatory element we name the 37-CE. Strikingly, changing only three nucleotides within this sequence is sufficient to render pneumococcus avirulent. Using in vivo and in vitro approaches, we present a model where SpxR interacts as a unique trimeric quaternary structure with the 37-CE to enable capsule repression in the airways. Considering its dramatic effect on infection, variation of the 37-CE between serotypes suggests this molecular switch could be a critical contributing factor to this pathogen's serotype-specific disease outcomes.

ABC-transporter CFTR folds with high fidelity through a modular, stepwise pathway

The question how proteins fold is especially pointed for large multi-domain, multi-spanning membrane proteins with complex topologies. We have uncovered the sequence of events that encompass proper folding of the ABC transporter CFTR in live cells by combining kinetic radiolabeling with protease-susceptibility assays. We found that CFTR folds in two clearly distinct stages. The first, co-translational, stage involves folding of the 2 transmembrane domains TMD1 and TMD2, plus one nucleotide-binding domain, NBD1. The second stage is a simultaneous, post-translational increase in protease resistance for both TMDs and NBD2, caused by assembly of these domains onto NBD1. Our assays probe every 2-3 residues (on average) in CFTR. This in-depth analysis at amino-acid level allows detailed analysis of domain folding and importantly also the next level: assembly of the domains into native, folded CFTR. Defects and changes brought about by medicines, chaperones, or mutations also are amenable to analysis. We here show that the well-known disease-causing mutation F508del, which established cystic fibrosis as protein-folding disease, caused co-translational misfolding of NBD1 but not TMD1 nor TMD2 in stage 1, leading to absence of stage-2 folding. Corrector drugs rescued stage 2 without rescuing NBD1. Likewise, the DxD motif in NBD1 that was identified to be required for export of CFTR from the ER we found to be required already upstream of export as CFTR mutated in this motif phenocopies F508del CFTR. The highly modular and stepwise folding process of such a large, complex protein explains the relatively high fidelity and correctability of its folding.

MLProbs: A Data-Centric Pipeline for Better Multiple Sequence Alignment

In this paper, we explore using the data-centric approach to tackle the Multiple Sequence Alignment (MSA) construction problem. Unlike the algorithm-centric approach, which reduces the construction problem to a combinatorial optimization problem based on an abstract mathematical model, the data-centric approach explores using classification models trained from existing benchmark data to guide the construction. We identified two simple classifications to help us choose a better alignment tool and determine whether and how much to carry out realignment. We show that shallow machine-learning algorithms suffice to train sensitive models for these classifications. Based on these models, we implemented a new multiple sequence alignment pipeline, called MLProbs. Compared with 10 other popular alignment tools over four benchmark databases (namely, BAliBASE, OXBench, OXBench-X and SABMark), MLProbs consistently gives the highest TC score. More importantly, MLProbs shows non-trivial improvement for protein families with low similarity; in particular, when evaluated against the 1,356 protein families with similarity ≤ 50%, MLProbs achieves a TC score of 56.93, while the next best three tools are in the range of [55.41, 55.91] (increased by more than 1.8%). We also compared the performance of MLProbs and other MSA tools in two real-life applications - Phylogenetic Tree Construction Analysis and Protein Secondary Structure Prediction - and MLProbs also had the best performance. In our study, we used only shallow machine-learning algorithms to train our models. It would be interesting to study whether deep-learning methods can help make further improvements, so we suggest some possible research directions in the conclusion section.

CryoEM structure of a post-assembly MS-ring reveals plasticity in stoichiometry and conformation

The flagellar motor supports bacterial chemotaxis, a process that allows bacteria to move in response to their environment. A central feature of this motor is the MS-ring, which is composed entirely of repeats of the FliF subunit. This MS-ring is critical for the assembly and stability of the flagellar switch and the entire flagellum. Despite multiple independent cryoEM structures of the MS-ring, there remains a debate about the stoichiometry and organization of the ring-building motifs (RBMs). Here, we report the cryoEM structure of a Salmonella MS-ring that was purified from the assembled flagellar switch complex (MSC-ring). We term this the 'post-assembly' state. Using 2D class averages, we show that under these conditions, the post-assembly MS-ring can contain 32, 33, or 34 FliF subunits, with 33 being the most common. RBM3 has a single location with C32, C33, or C34 symmetry. RBM2 is found in two locations with RBM2inner having C21 or C22 symmetry and an RBM2outer-RBM1 having C11 symmetry. Comparison to previously reported structures identifies several differences. Most strikingly, we find that the membrane domain forms 11 regions of discrete density at the base of the structure rather than a contiguous ring, although density could not be unambiguously interpreted. We further find density in some previously unresolved areas, and we assigned amino acids to those regions. Finally, we find differences in interdomain angles in RBM3 that affect the diameter of the ring. Together, these investigations support a model of the flagellum with structural plasticity, which may be important for flagellar assembly and function.

GLB-3: A resilient, cysteine-rich, membrane-tethered globin expressed in the reproductive and nervous system of Caenorhabditis elegans

The popular genetic model organism Caenorhabditis elegans (C. elegans) encodes 34 globins, whereby the few that are well-characterized show divergent properties besides the typical oxygen carrier function. Here, we present a biophysical characterization and expression analysis of C. elegans globin-3 (GLB-3). GLB-3 is predicted to exist in two isoforms and is expressed in the reproductive and nervous system. Knockout of this globin causes a 99% reduction in fertility and reduced motility. Spectroscopic analysis reveals that GLB-3 exists as a bis-histidyl-ligated low-spin form in both the ferrous and ferric heme form. A function in binding of diatomic gases is excluded on the basis of the slow CO-binding kinetics. Unlike other globins, GLB-3 is also not capable of reacting with H₂O₂, H₂S, and nitrite. Intriguingly, not only does GLB-3 contain a high number of cysteine residues, it is also highly stable under harsh conditions (pH = 2 and high concentrations of H₂O₂). The resilience diminishes when the N- and C-terminal extensions are removed. Redox potentiometric measurements reveal a slightly positive redox potential (+8 ± 19 mV vs. SHE), suggesting that the heme iron may be able to oxidize cysteines. Electron paramagnetic resonance shows that formation of an intramolecular disulphide bridge, involving Cys70, affects the heme-pocket region. The results suggest an involvement of the globin in (cysteine) redox chemistry.

Structures of RecBCD in complex with phage-encoded inhibitor proteins reveal distinctive strategies for evasion of a bacterial immunity hub

Following infection of bacterial cells, bacteriophage modulate double-stranded DNA break repair pathways to protect themselves from host immunity systems and prioritise their own recombinases. Here, we present biochemical and structural analysis of two phage proteins, gp5.9 and Abc2, which target the DNA break resection complex RecBCD. These exemplify two contrasting mechanisms for control of DNA break repair in which the RecBCD complex is either inhibited or co-opted for the benefit of the invading phage. Gp5.9 completely inhibits RecBCD by preventing it from binding to DNA. The RecBCD-gp5.9 structure shows that gp5.9 acts by substrate mimicry, binding predominantly to the RecB arm domain and competing sterically for the DNA binding site. Gp5.9 adopts a parallel coiled-coil architecture that is unprecedented for a natural DNA mimic protein. In contrast, binding of Abc2 does not substantially affect the biochemical activities of isolated RecBCD. The RecBCD-Abc2 structure shows that Abc2 binds to the Chi-recognition domains of the RecC subunit in a position that might enable it to mediate the loading of phage recombinases onto its single-stranded DNA products.

Pan-kinome of Legionella expanded by a bioinformatics survey

The pathogenic Legionella bacteria are notorious for delivering numerous effector proteins into the host cell with the aim of disturbing and hijacking cellular processes for their benefit. Despite intensive studies, many effectors remain uncharacterized. Motivated by the richness of Legionella effector repertoires and their oftentimes atypical biochemistry, also by several known atypical Legionella effector kinases and pseudokinases discovered recently, we undertook an in silico survey and exploration of the pan-kinome of the Legionella genus, i.e., the union of the kinomes of individual species. In this study, we discovered 13 novel (pseudo)kinase families (all are potential effectors) with the use of non-standard bioinformatic approaches. Together with 16 known families, we present a catalog of effector and non-effector protein kinase-like families within Legionella, available at http://bioinfo.sggw.edu.pl/kintaro/ . We analyze and discuss the likely functional roles of the novel predicted kinases. Notably, some of the kinase families are also present in other bacterial taxa, including other pathogens, often phylogenetically very distant from Legionella. This work highlights Nature's ingeniousness in the pathogen-host arms race and offers a useful resource for the study of infection mechanisms.

Membrane curvature sensing and stabilization by the autophagic LC3 lipidation machinery

How the highly curved phagophore membrane is stabilized during autophagy initiation is a major open question in autophagosome biogenesis. Here, we use in vitro reconstitution on membrane nanotubes and molecular dynamics simulations to investigate how core autophagy proteins in the LC3 (Microtubule-associated proteins 1A/1B light chain 3) lipidation cascade interact with curved membranes, providing insight into their possible roles in regulating membrane shape during autophagosome biogenesis. ATG12(Autophagy-related 12)-ATG5-ATG16L1 was up to 100-fold enriched on highly curved nanotubes relative to flat membranes. At high surface density, ATG12-ATG5-ATG16L1 binding increased the curvature of the nanotubes. While WIPI2 (WD repeat domain phosphoinositide-interacting protein 2) binding directs membrane recruitment, the amphipathic helix α2 of ATG16L1 is responsible for curvature sensitivity. Molecular dynamics simulations revealed that helix α2 of ATG16L1 inserts shallowly into the membrane, explaining its curvature-sensitive binding to the membrane. These observations show how the binding of the ATG12-ATG5-ATG16L1 complex to the early phagophore rim could stabilize membrane curvature and facilitate autophagosome growth.

In silico investigation on structure-function relationship of members belonging to the human SLC52 transporter family

Riboflavin is an essential water-soluble vitamin that needs to be provided through the diet because of the conversion into flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN), important cofactors in hundreds of flavoenzymes. The adsorption and distribution of riboflavin is mediated by transmembrane transporters of the SLC52 family, namely RFVT1-3, whose mutations are mainly associated with two diseases, MADD and the Brown-Vialetto-Van Laere syndrome. Interest in RFVTs as pharmacological targets has increased in the last few years due to their overexpression in several cancer cells, which can be exploited both by blocking the uptake of riboflavin into the cancerous cells, and by performing cancer targeted delivery of drugs with a high affinity for RFVTs. In this work, we propose three-dimensional structural models for all three human riboflavin transporters obtained by state-of-the-art artificial intelligence-based methods, which were then further refined with molecular dynamics simulations. Furthermore, two of the most notable mutations concerning RFVT2 and RFVT3 (W31S and N21S, respectively) were investigated studying the interactions between the wild-type and mutated transporters with riboflavin.

A novel post-developmental role of the Hox genes underlies normal adult behavior

The molecular mechanisms underlying the stability of mature neurons and neural circuits are poorly understood. Here we explore this problem and discover that the Hox genes are a component of the genetic program that maintains normal neural function in adult Drosophila. We show that post-developmental downregulation of the Hox gene Ultrabithorax (Ubx) in adult neurons leads to substantial anomalies in flight. Mapping the cellular basis of these effects reveals that Ubx is required within a subset of dopaminergic neurons, and cell circuitry analyses in combination with optogenetics allow us to link these dopaminergic neurons to flight control. Functional imaging experiments show that Ubx is necessary for normal dopaminergic activity, and neuron-specific RNA-sequencing defines two previously uncharacterized ion channel-encoding genes as potential mediators of Ubx behavioral roles. Our study thus reveals a novel role of the Hox system in controlling adult behavior and neural function. Based on the broad evolutionary conservation of the Hox system across distantly related animal phyla, we predict that the Hox genes might play neurophysiological roles in adult forms of other species, including humans.

Cryo-EM structure of ssDNA bacteriophage ΦCjT23 provides insight into early virus evolution

The origin of viruses remains an open question. While lack of detectable sequence similarity hampers the analysis of distantly related viruses, structural biology investigations of conserved capsid protein structures facilitate the study of distant evolutionary relationships. Here we characterize the lipid-containing ssDNA temperate bacteriophage ΦCjT23, which infects Flavobacterium sp. (Bacteroidetes). We report ΦCjT23-like sequences in the genome of strains belonging to several Flavobacterium species. The virion structure determined by cryogenic electron microscopy reveals similarities to members of the viral kingdom Bamfordvirae that currently consists solely of dsDNA viruses with a major capsid protein composed of two upright β-sandwiches. The minimalistic structure of ΦCjT23 suggests that this phage serves as a model for the last common ancestor between ssDNA and dsDNA viruses in the Bamfordvirae. Both ΦCjT23 and the related phage FLiP infect Flavobacterium species found in several environments, suggesting that these types of viruses have a global distribution and a shared evolutionary origin. Detailed comparisons to related, more complex viruses not only expand our knowledge about this group of viruses but also provide a rare glimpse into early virus evolution.

Phylogenetic Analyses of Microbial Hydrolytic Dehalogenases Reveal Polyphyletic Origin

Hydrolytic dehalogenases form an important class of dehalogenases that include haloacid dehalogenase, haloalkane dehalogenase, haloacetate dehalogenase, and atrazine chlorohydrolase. These enzymes are involved in biodegradation of various environmental pollutants and therefore it is important to understand their phylogeny. In the present study, it was found that the enzymes haloalkane and haloacetate dehalogenases share a common ancestry with enzymes such as carboxyesterase, epoxide hydrolase, and lipases, which can be traced to ancestral α/β hydrolase fold enzyme. Haloacid dehalogenases and atrazine chlorohydrolases have probabaly evolved from ancestral enzymes with phosphatase and deaminases activity, respectively. These findings were supported by the similarities in the secondary structure, key catalytic motifs and placement of catalytic residues. The phylogeny of haloalkane dehalogenases and haloacid dehalogenases differs from 16S rRNA gene phylogeny, suggesting spread through horizontal gene transfer. Hydrolytic dehalogenases are polyphyletic and do not share a common evolutionay history, the functional similarities are due to convergent evolution. The present study also identifies key functional residues, mutating which, can help in generating better enzymes for clean up of the persistent environmental pollutants using enzymatic bioremediation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-022-01043-8.

Classification and functional characterization of spidroin genes in a wandering spider, Pardosa pseudoannulata

Spiders impress us with their sophisticated use of silk and the stunningly distinct silk proteins (spidroins) in each spider species. Understanding how silks and spidroins function and evolve within the spider world is one profound interest to expand our knowledge on spider evolution. Spidroins are characterized with the divergent repeat core region flanked with the relatively conserved N- and C-terminus. The structure and number of the repeats contribute to the unique mechanical properties of the spidroin and the silk. Spidroins have been intensively studied in web-weaver spiders, but information regarding their diversity in wandering spiders remains scarce. Here, twenty spidroin genes were identified in the pond wolf spider, Pardosa pseudoannulata, belonging to the retrolateral tibial apophysis (RTA) clade. These spidroins were categorized into four classes, including twelve ampullate spidroin (AmpSp), four aciniform spidroin (AcSp), one tubuliform spidroin (TuSp), one pyriform spidroin (PiSp), and two spidroin-like proteins. Multiple copies of the AmpSp and AcSp genes were tandemly arranged in a cluster within the genome, and the N-terminal domains and repetitive sequences of the proximately located spidroins were highly similar, suggesting that the spidroin genes diversified via tandem duplication. Only four types of morphologically distinct silk glands were found in P. pseudoannulata, namely Ma, Mi, Ac, and Pi glands, consistent with the glandular affiliation hypothesis that spidroins co-evolved with glandular specialization to fit species-specific needs. Expression profiling revealed that the single tubuliform spidroin (TuSp) gene was highly expressed in gravid females and two AcSp genes displayed synchronous expression. Knock-down of the TuSp gene via RNAi resulted in fragile and cracked eggsacs and prolonged the female pre-oviposition period, validating its importance in spider reproduction. The genome-scale characterization and functional study of spidroin genes allows associating the presence of specific spidroins with silk utility in P. pseudoannulata and will expand our knowledge of spider evolution.

Response regulator PorX coordinates oligonucleotide signalling and gene expression to control the secretion of virulence factors

The PglZ family of proteins belongs to the alkaline phosphatase superfamily, which consists of metallohydrolases with limited sequence identity but similar metal-coordination architectures in otherwise divergent active sites. Proteins with a well-defined PglZ domain are ubiquitous among prokaryotes as essential components of BREX phage defence systems and two-component systems (TCSs). Whereas other members of the alkaline phosphatase superfamily are well characterized, the activity, structure and biological function of PglZ family proteins remain unclear. We therefore investigated the structure and function of PorX, an orphan response regulator of the Porphyromonas gingivalis TCS containing a putative PglZ effector domain. The crystal structure of PorX revealed a canonical receiver domain, a helical bundle, and an unprecedented PglZ domain, similar to the general organization of the phylogenetically related BREX-PglZ proteins. The PglZ domain of PorX features an active site cleft suitable for large substrates. An extensive search for substrates revealed that PorX is a phosphodiesterase that acts on cyclic and linear oligonucleotides, including signalling molecules such as cyclic oligoadenylates. These results, combined with mutagenesis, biophysical and enzymatic analysis, suggest that PorX coordinates oligonucleotide signalling pathways and indirectly regulates gene expression to control the secretion of virulence factors.

The permanently chaperone-active small heat shock protein Hsp17 from Caenorhabditis elegans exhibits topological separation of its N-terminal regions

Small Heat shock proteins (sHsps) are a family of molecular chaperones that bind nonnative proteins in an ATP-independent manner. Caenorhabditis elegans encodes 16 different sHsps, among them Hsp17, which is evolutionarily distinct from other sHsps in the nematode. The structure and mechanism of Hsp17 and how these may differ from other sHsps remain unclear. Here, we find that Hsp17 has a distinct expression pattern, structural organization, and chaperone function. Consistent with its presence under nonstress conditions, and in contrast to many other sHsps, we determined that Hsp17 is a mono-disperse, permanently active chaperone in vitro, which interacts with hundreds of different C. elegans proteins under physiological conditions. Additionally, our cryo-EM structure of Hsp17 reveals that in the 24-mer complex, 12 N-terminal regions are involved in its chaperone function. These flexible regions are located on the outside of the spherical oligomer, whereas the other 12 N-terminal regions are engaged in stabilizing interactions in its interior. This allows the same region in Hsp17 to perform different functions depending on the topological context. Taken together, our results reveal structural and functional features that further define the structural basis of permanently active sHsps.

Molecular basis for DNA recognition by the maternal pioneer transcription factor FoxH1

Forkhead box H1 (FoxH1) is an essential maternal pioneer factor during embryonic development that binds to specific GG/GT-containing DNA target sequences. Here we have determined high-resolution structures of three FoxH1 proteins (from human, frog and fish species) and four DNAs to clarify the way in which FoxH1 binds to these sites. We found that the protein-DNA interactions extend to both the minor and major DNA grooves and are thus almost twice as extensive as those of other FOX family members. Moreover, we identified two specific amino acid changes in FoxH1 that allowed the recognition of GG/GT motifs. Consistent with the pioneer factor activity of FoxH1, we found that its affinity for nucleosomal DNA is even higher than for linear DNA fragments. The structures reported herein illustrate how FoxH1 binding to distinct DNA sites provides specificity and avoids cross-regulation by other FOX proteins that also operate during the maternal-zygotic transition and select canonical forkhead sites.

DLBLS_SS: protein secondary structure prediction using deep learning and broad learning system

Protein secondary structure prediction (PSSP) is not only beneficial to the study of protein structure and function but also to the development of drugs. As a challenging task in computational biology, experimental methods for PSSP are time-consuming and expensive. In this paper, we propose a novel PSSP model DLBLS_SS based on deep learning and broad learning system (BLS) to predict 3-state and 8-state secondary structure. We first use a bidirectional long short-term memory (BLSTM) network to extract global features in residue sequences. Then, our proposed SEBTCN based on temporal convolutional networks (TCN) and channel attention can capture bidirectional key long-range dependencies in sequences. We also use BLS to rapidly optimize fused features while further capturing local interactions between residues. We conduct extensive experiments on public test sets including CASP10, CASP11, CASP12, CASP13, CASP14 and CB513 to evaluate the performance of the model. Experimental results show that our model exhibits better 3-state and 8-state PSSP performance compared to five state-of-the-art models.

Role of DNA modifications in Mycoplasma gallisepticum

The epigenetics of bacteria, and bacteria with a reduced genome in particular, is of great interest, but is still poorly understood. Mycoplasma gallisepticum, a representative of the class Mollicutes, is an excellent model of a minimal cell because of its reduced genome size, lack of a cell wall, and primitive cell organization. In this study we investigated DNA modifications of the model object Mycoplasma gallisepticum and their roles. We identified DNA modifications and methylation motifs in M. gallisepticum S6 at the genome level using single molecule real time (SMRT) sequencing. Only the ANCNNNNCCT methylation motif was found in the M. gallisepticum S6 genome. The studied bacteria have one functional system for DNA modifications, the Type I restriction-modification (RM) system, MgaS6I. We characterized its activity, affinity, protection and epigenetic functions. We demonstrated the protective effects of this RM system. A common epigenetic signal for bacteria is the m6A modification we found, which can cause changes in DNA-protein interactions and affect the cell phenotype. Native methylation sites are underrepresented in promoter regions and located only near the -35 box of the promoter, which does not have a significant effect on gene expression in mycoplasmas. To study the epigenetics effect of m6A for genome-reduced bacteria, we constructed a series of M. gallisepticum strains expressing EGFP under promoters with the methylation motifs in their different elements. We demonstrated that m6A modifications of the promoter located only in the -10-box affected gene expression and downregulated the expression of the corresponding gene.

Insights into the role of Nup62 and Nup93 in assembling cytoplasmic ring and central transport channel of the nuclear pore complex

The nuclear pore complex (NPC) is a highly modular assembly of 34 distinct nucleoporins (Nups) to form a versatile transport channel between the nucleus and the cytoplasm. Among them, Nup62 is known as an essential component for nuclear transport, Nup93 for proper nuclear envelope assembly. These Nups constitute various NPC subcomplexes such as the central transport channel (CTC), the cytoplasmic ring (CR), and the inner ring (IR). However, how they play their roles in NPC assembly and transport activity is not clear. Here we delineated the interacting regions and conducted biochemical reconstitution and structural characterization of the mammalian CR complex to reveal its intrinsic dynamic behavior and a distinct "4"-shaped architecture resembling the CTC complex. Our in vitro reconstitution data demonstrate that the Nup62 coiled-coil domain is critical to form both Nup62^322-525 •Nup88^517-742 and Nup62^322-525•Nup88^517-742•Nup214^693-926 heterotrimers and both can bind to Nup93^1-150. We therefore propose that Nup93 acts as a "sensor" to bind to Nup62 shared heterotrimers including the Nup62•Nup54 heterotrimer of the CTC, which was not shown previously to be an interacting partner. Altogether, our biochemical study suggests that Nup62 via its coiled-coil domain is central to form compositionally distinct yet structurally similar heterotrimers and Nup93 binds these diverse heterotrimers nonselectively.

Single-Step Protocol for Isolating the Recombinant Extracellular Domain of the Luteinizing Hormone Receptor from the Ovis aries Testis

The luteinizing hormone receptor (LHR) is a glycoprotein member of the G protein-coupled receptors superfamily. It participates in corpus luteum formation and ovulation in females and acts in testosterone synthesis and spermatogenesis in males. In this study, we extracted RNA from sheep testicles and synthetized the cDNA to amplify the gene lhr-bed. This gene consists of 762 bp and encodes 273 amino acids of the extracellular domain of LHR. The lhr-bed was cloned into pJET1.2/blunt, then subcloned into pCOLD II, and finally, transformed in E. coli BL21 (DE3) cells. Because the induced rLHR-Bed protein was found in the insoluble fraction, we followed a modified purification protocol involving induction at 25 °C, subjection to denaturing conditions, and on-column refolding to increase solubility. We confirmed rLHR-Bed expression by means of Western blot and mass spectrometry analysis. It is currently known that the structure stem-loop 5′UTR on pCOLD II vector is stable at 15 °C. We predicted and obtained RNAfold stability at 25 °C. We successfully obtained the recombinant LHR extracellular domain, with protein yields of 0.2 mg/L, and purity levels of approximately 90%, by means of a single chromatographic purification step. The method described here may be used to obtain large quantities of rLHR-Bed in the future.