PredictProtein--an open resource for online prediction of protein structural and functional features

Low frequency variants associated with leukocyte telomere length in the Singapore Chinese population

The role of low frequency variants associated with telomere length homeostasis in chronic diseases and mortalities is relatively understudied in the East-Asian population. Here we evaluated low frequency variants, including 1,915,154 Asian specific variants, for leukocyte telomere length (LTL) associations among 25,533 Singapore Chinese samples. Three East Asian specific variants in/near POT1, TERF1 and STN1 genes are associated with LTL (Meta-analysis P 2.49×10-14-6.94×10-10). Rs79314063, a missense variant (p.Asp410His) at POT1, shows effect 5.3 fold higher and independent of a previous common index SNP. TERF1 (rs79617270) and STN1 (rs139620151) are linked to LTL-associated common index SNPs at these loci. Rs79617270 is associated with cancer mortality [HR95%CI = 1.544 (1.173, 2.032), PAdj = 0.018] and 4.76% of the association between the rs79617270 and colon cancer is mediated through LTL. Overall, genetically determined LTL is particularly associated with lung adenocarcinoma [HR95%CI = 1.123 (1.051, 1.201), Padj = 0.007]. Ethnicity-specific low frequency variants may affect LTL homeostasis and associate with certain cancers.

Eusociality Shapes Convergent Patterns of Molecular Evolution across Mitochondrial Genomes of Snapping Shrimps

Eusociality is a highly conspicuous and ecologically impactful behavioral syndrome that has evolved independently across multiple animal lineages. So far, comparative genomic analyses of advanced sociality have been mostly limited to insects. Here, we study the only clade of animals known to exhibit eusociality in the marine realm-lineages of socially diverse snapping shrimps in the genus Synalpheus. To investigate the molecular impact of sociality, we assembled the mitochondrial genomes of eight Synalpheus species that represent three independent origins of eusociality and analyzed patterns of molecular evolution in protein-coding genes. Synonymous substitution rates are lower and potential signals of relaxed purifying selection are higher in eusocial relative to noneusocial taxa. Our results suggest that mitochondrial genome evolution was shaped by eusociality-linked traits-extended generation times and reduced effective population sizes that are hallmarks of advanced animal societies. This is the first direct evidence of eusociality impacting genome evolution in marine taxa. Our results also strongly support the idea that eusociality can shape genome evolution through profound changes in life history and demography.

Functional interpretation of ATAD3A variants in neuro-mitochondrial phenotypes

Background

ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane-anchored protein involved in diverse processes including mitochondrial dynamics, mitochondrial DNA organization, and cholesterol metabolism. Biallelic deletions (null), recessive missense variants (hypomorph), and heterozygous missense variants or duplications (antimorph) in ATAD3A lead to neurological syndromes in humans.

Methods

To expand the mutational spectrum of ATAD3A variants and to provide functional interpretation of missense alleles in trans to deletion alleles, we performed exome sequencing for identification of single nucleotide variants (SNVs) and copy number variants (CNVs) in ATAD3A in individuals with neurological and mitochondrial phenotypes. A Drosophila Atad3a Gal4 knockin-null allele was generated using CRISPR-Cas9 genome editing technology to aid the interpretation of variants.

Results

We report 13 individuals from 8 unrelated families with biallelic ATAD3A variants. The variants included four missense variants inherited in trans to loss-of-function alleles (p.(Leu77Val), p.(Phe50Leu), p.(Arg170Trp), p.(Gly236Val)), a homozygous missense variant p.(Arg327Pro), and a heterozygous non-frameshift indel p.(Lys568del). Affected individuals exhibited findings previously associated with ATAD3A pathogenic variation, including developmental delay, hypotonia, congenital cataracts, hypertrophic cardiomyopathy, and cerebellar atrophy. Drosophila studies indicated that Phe50Leu, Gly236Val, Arg327Pro, and Lys568del are severe loss-of-function alleles leading to early developmental lethality. Further, we showed that Phe50Leu, Gly236Val, and Arg327Pro cause neurogenesis defects. On the contrary, Leu77Val and Arg170Trp are partial loss-of-function alleles that cause progressive locomotion defects and whose expression leads to an increase in autophagy and mitophagy in adult muscles.

Conclusion

Our findings expand the allelic spectrum of ATAD3A variants and exemplify the use of a functional assay in Drosophila to aid variant interpretation.

Integrative omics identification, evolutionary and structural analysis of low affinity nitrate transporters in diatoms, diNPFs

Diatoms are one of the major and most diverse groups of phytoplankton, with chimeric genomes harbouring a combination of genes of bacterial, animal and plant origin. They have developed sophisticated mechanisms to face environmental variations. In marine environments, nutrients concentration shows significant temporal and spatial variability, influencing phytoplankton growth. Among nutrients, nitrogen, present at micromolar levels, is often a limiting resource. Here, we report a comprehensive characterization of the Nitrate Transporter 1/Peptide Transporter Family (NPF) in diatoms, diNPFs. NPFs are well characterized in many organisms where they recognize a broad range of substrates, ranging from short-chained di- and tri-peptides in bacteria, fungi and mammals to a wide variety of molecules including nitrate in higher plants. Scarce information is available for diNPFs. We integrated-omics, phylogenetic, structural and expression analyses, to infer information on their role in diatoms. diNPF genes diverged to produce two distinct clades with strong sequence and structural homology with either bacterial or plant NPFs, with different predicted sub-cellular localization, suggesting that the divergence resulted in functional diversification. Moreover, transcription analysis of diNPF genes under different laboratory and environmental growth conditions suggests that diNPF diversification led to genetic adaptations that might contribute to diatoms ability to flourish in diverse environmental conditions.

Designing of a chimeric protein contains StxB, intimin and EscC against toxicity and adherence of enterohemorrhagic Escherichia coli O157:H7 and evaluation of serum antibody titers against it

Enterohemorrhagic Escherichia coli (EHEC) O157:H7 strain is known as one of the major human foodborne pathogens. Lack of effective clinical treatment for human diarrheal diseases confirms the need for vaccine production against enteric bacteria such as E.coli O157:H7. Shiga-like toxin (Stx), EscC, and Intimin are the main important virulent factors of this enteric pathogen. In the present study, a comparative Omics analysis was conducted to identify most invasion EHEC antigenic factors as a potential immunogen. SEI (Stx-EscC-Intimin) trivalent chimeric protein was designed from the exposed and epitope rich part of these virulence factors. Sequence optimization, physicochemical properties, mRNA folding, three-dimensional structure and immunoinformatics data were investigated. The chimeric gene was synthesized with codon bias of E. coli. Recombinant protein was expressed and confirmed by western blot analysis. To evaluate the immunogenicity of the designed protein, the protein was administered to BALB/c mice and the serum IgG was determined by ELISA. Based on the Ramachandran plot, the validation data showed that 90.1 % of residues lie in the favored region. The high antigenicity of the multimeric protein was predicted by the immunoinformatic analysis. Epitope prediction had shown the proper distribution of linear and conformational B-cell epitopes and the competition of T-cell epitopes to bind MHC molecules too. Recombinant ESI Protein with 74.5 kDa was expressed in E. coli. Western blot analysis by anti-Stx antibody, confirmed a single band of chimeric protein. Consequently, the chimeric gene was designed and constructed after assessments. From in silico approach, the protein deduced from this cassette can be an immunogen candidate, and act against toxicity and adherence of EHEC.

Characterization of the Bubblegum acyl-CoA synthetase of Microchloropsis gaditana

The metabolic pathways of glycerolipids are well described in cells containing chloroplasts limited by a two-membrane envelope but not in cells containing plastids limited by four membranes, including heterokonts. Fatty acids (FAs) produced in the plastid, palmitic and palmitoleic acids (16:0 and 16:1), are used in the cytosol for the synthesis of glycerolipids via various routes, requiring multiple acyl-Coenzyme A (CoA) synthetases (ACS). Here, we characterized an ACS of the Bubblegum subfamily in the photosynthetic eukaryote Microchloropsis gaditana, an oleaginous heterokont used for the production of lipids for multiple applications. Genome engineering with TALE-N allowed the generation of MgACSBG point mutations, but no knockout was obtained. Point mutations triggered an overall decrease of 16:1 in lipids, a specific increase of unsaturated 18-carbon acyls in phosphatidylcholine and decrease of 20-carbon acyls in the betaine lipid diacylglyceryl-trimethyl-homoserine. The profile of acyl-CoAs highlighted a decrease in 16:1-CoA and 18:3-CoA. Structural modeling supported that mutations affect accessibility of FA to the MgACSBG reaction site. Expression in yeast defective in acyl-CoA biosynthesis further confirmed that point mutations affect ACSBG activity. Altogether, this study supports a critical role of heterokont MgACSBG in the production of 16:1-CoA and 18:3-CoA. In M. gaditana mutants, the excess saturated and monounsaturated FAs were diverted to triacylglycerol, thus suggesting strategies to improve the oil content in this microalga.

Structure of the cytoplasmic domain of SctV (SsaV) from the Salmonella SPI-2 injectisome and implications for a pH sensing mechanism

•

CryoEM of a full-length type III secretion system SctV resolves cytoplasmic but not transmembrane domains.

•

MD simulations show SctV protomers flexibly hinge.

•

Acidification expands the SctV ring by altering interprotomer interactions.

Bacterial type III secretion systems assemble the axial structures of both injectisomes and flagella. Injectisome type III secretion systems subsequently secrete effector proteins through their hollow needle into a host, requiring co-ordination. In the Salmonella enterica serovar Typhimurium SPI-2 injectisome, this switch is triggered by sensing the neutral pH of the host cytoplasm. Central to specificity switching is a nonameric SctV protein with an N-terminal transmembrane domain and a toroidal C-terminal cytoplasmic domain. A ‘gatekeeper’ complex interacts with the SctV cytoplasmic domain in a pH dependent manner, facilitating translocon secretion while repressing effector secretion through a poorly understood mechanism. To better understand the role of SctV in SPI-2 translocon-effector specificity switching, we purified full-length SctV and determined its toroidal cytoplasmic region’s structure using cryo-EM. Structural comparisons and molecular dynamics simulations revealed that the cytoplasmic torus is stabilized by its core subdomain 3, about which subdomains 2 and 4 hinge, varying the flexible outside cleft implicated in gatekeeper and substrate binding. In light of patterns of surface conservation, deprotonation, and structural motion, the location of previously identified critical residues suggest that gatekeeper binds a cleft buried between neighboring subdomain 4s. Simulations suggest that a local pH change from 5 to 7.2 stabilizes the subdomain 3 hinge and narrows the central aperture of the nonameric torus. Our results are consistent with a model of local pH sensing at SctV, where pH-dependent dynamics of SctV cytoplasmic domain affect binding of gatekeeper complex.

Proteogenomic Characterization of the Cement and Adhesive Gland of the Pelagic Gooseneck Barnacle Lepas anatifera

We focus on the stalked goose barnacle L. anatifera adhesive system, an opportunistic less selective species for the substrate, found attached to a variety of floating objects at seas. Adhesion is an adaptative character in barnacles, ensuring adequate positioning in the habitat for feeding and reproduction. The protein composition of the cement multicomplex and adhesive gland was quantitatively studied using shotgun proteomic analysis. Overall, 11,795 peptide sequences were identified in the gland and 2206 in the cement, clustered in 1689 and 217 proteinGroups, respectively. Cement specific adhesive proteins (CPs), proteases, protease inhibitors, cuticular and structural proteins, chemical cues, and many unannotated proteins were found, among others. In the cement, CPs were the most abundant (80.5%), being the bulk proteins CP100k and -52k the most expressed of all, and CP43k-like the most expressed interfacial protein. Unannotated proteins comprised 4.7% of the cement proteome, ranking several of them among the most highly expressed. Eight of these proteins showed similar physicochemical properties and amino acid composition to known CPs and classified through Principal Components Analysis (PCA) as new CPs. The importance of PCA on the identification of unannotated non-conserved adhesive proteins, whose selective pressure is on their relative amino acid abundance, was demonstrated.

Protein cleavage influences surface protein presentation in Mycoplasma pneumoniae

Mycoplasma pneumoniae is a significant cause of pneumonia and post infection sequelae affecting organ sites distant to the respiratory tract are common. It is also a model organism where extensive 'omics' studies have been conducted to gain insight into how minimal genome self-replicating organisms function. An N-terminome study undertaken here identified 4898 unique N-terminal peptides that mapped to 391 (56%) predicted M. pneumoniae proteins. True N-terminal sequences beginning with the initiating methionine (iMet) residue from the predicted Open Reading Frame (ORF) were identified for 163 proteins. Notably, almost half (317; 46%) of the ORFS derived from M. pneumoniae strain M129 are post-translationally modified, presumably by proteolytic processing, because dimethyl labelled neo-N-termini were characterised that mapped beyond the predicted N-terminus. An analysis of the N-terminome describes endoproteolytic processing events predominately targeting tryptic-like sites, though cleavages at negatively charged residues in P1' (D and E) with lysine or serine/alanine in P2' and P3' positions also occurred frequently. Surfaceome studies identified 160 proteins (23% of the proteome) to be exposed on the extracellular surface of M. pneumoniae. The two orthogonal methodologies used to characterise the surfaceome each identified the same 116 proteins, a 72% (116/160) overlap. Apart from lipoproteins, transporters, and adhesins, 93/160 (58%) of the surface proteins lack signal peptides and have well characterised, canonical functions in the cell. Of the 160 surface proteins identified, 134 were also targets of endo-proteolytic processing. These processing events are likely to have profound implications for how the host immune system recognises and responds to M. pneumoniae.

In Silico Analysis of Huntingtin Homologs in Lower Eukaryotes

Huntington’s disease is a rare neurodegenerative and autosomal dominant disorder. HD is caused by a mutation in the gene coding for huntingtin (Htt). The result is the production of a mutant Htt with an abnormally long polyglutamine repeat that leads to pathological Htt aggregates. Although the structure of human Htt has been determined, albeit at low resolution, its functions and how they are performed are largely unknown. Moreover, there is little information on the structure and function of Htt in other organisms. The comparison of Htt homologs can help to understand if there is a functional conservation of domains in the evolution of Htt in eukaryotes. In this work, through a computational approach, Htt homologs from lower eukaryotes have been analysed, identifying ordered domains and modelling their structure. Based on the structural models, a putative function for most of the domains has been predicted. A putative C. elegans Htt-like protein has also been analysed following the same approach. The results obtained support the notion that this protein is a orthologue of human Htt.

The impact of structural bioinformatics tools and resources on SARS-CoV-2 research and therapeutic strategies

SARS-CoV-2 is the causative agent of COVID-19, the ongoing global pandemic. It has posed a worldwide challenge to human health as no effective treatment is currently available to combat the disease. Its severity has led to unprecedented collaborative initiatives for therapeutic solutions against COVID-19. Studies resorting to structure-based drug design for COVID-19 are plethoric and show good promise. Structural biology provides key insights into 3D structures, critical residues/mutations in SARS-CoV-2 proteins, implicated in infectivity, molecular recognition and susceptibility to a broad range of host species. The detailed understanding of viral proteins and their complexes with host receptors and candidate epitope/lead compounds is the key to developing a structure-guided therapeutic design. Since the discovery of SARS-CoV-2, several structures of its proteins have been determined experimentally at an unprecedented speed and deposited in the Protein Data Bank. Further, specialized structural bioinformatics tools and resources have been developed for theoretical models, data on protein dynamics from computer simulations, impact of variants/mutations and molecular therapeutics. Here, we provide an overview of ongoing efforts on developing structural bioinformatics tools and resources for COVID-19 research. We also discuss the impact of these resources and structure-based studies, to understand various aspects of SARS-CoV-2 infection and therapeutic development. These include (i) understanding differences between SARS-CoV-2 and SARS-CoV, leading to increased infectivity of SARS-CoV-2, (ii) deciphering key residues in the SARS-CoV-2 involved in receptor-antibody recognition, (iii) analysis of variants in host proteins that affect host susceptibility to infection and (iv) analyses facilitating structure-based drug and vaccine design against SARS-CoV-2.

Low Diversity of Human Variation Despite Mostly Mild Functional Impact of De Novo Variants

Non-synonymous Single Nucleotide Variants (nsSNVs), resulting in single amino acid variants (SAVs), are important drivers of evolutionary adaptation across the tree of life. Humans carry on average over 10,000 SAVs per individual genome, many of which likely have little to no impact on the function of the protein they affect. Experimental evidence for protein function changes as a result of SAVs remain sparse – a situation that can be somewhat alleviated by predicting their impact using computational methods. Here, we used SNAP to examine both observed and in silico generated human variation in a set of 1,265 proteins that are consistently found across a number of diverse species. The number of SAVs that are predicted to have any functional effect on these proteins is smaller than expected, suggesting sequence/function optimization over evolutionary timescales. Additionally, we find that only a few of the yet-unobserved SAVs could drastically change the function of these proteins, while nearly a quarter would have only a mild functional effect. We observed that variants common in the human population localized to less conserved protein positions and carried mild to moderate functional effects more frequently than rare variants. As expected, rare variants carried severe effects more frequently than common variants. In line with current assumptions, we demonstrated that the change of the human reference sequence amino acid to the reference of another species (a cross-species variant) is unlikely to significantly impact protein function. However, we also observed that many cross-species variants may be weakly non-neutral for the purposes of quick adaptation to environmental changes, but may not be identified as such by current state-of-the-art methodology.

Function Analysis of MBF1, a Factor Involved in the Response to Amino Acid Starvation and Virulence in Candida albicans

Candida albicans is a commensal of human mucosae, but also one of the most common fungal pathogens of humans. Systemic infections caused by this fungus, mostly affecting immunocompromised patients, are associated to fatality rates as high as 50% despite the available treatments. In order to improve this situation, it is necessary to fully understand how C. albicans is able to cause disease and how it copes with the host defenses. Our previous studies have revealed the importance of the C. albicans gene MBF1 in virulence and ability to colonize internal organs of mammalian and insect hosts. MBF1 encodes a putative transcriptional regulator, and as such it likely has an impact in the regulation of C. albicans gene expression during host infection. Here, recent advances in RNA-seq technologies were used to obtain a detailed analysis of the impact of MBF1 on C. albicans gene expression both in vitro and during infection. MBF1 was involved in the regulation of several genes with a role in glycolysis and response to stress, particularly to nutritional stress. We also investigated whether an interaction existed between MBF1 and GCN4, a master regulator of response to starvation, and found that both genes were needed for resistance to amino acid starvation, suggesting some level of interaction between the two. Reinforcing this idea, we showed that the proteins encoded by both genes could interact. Consistent with the role of MBF1 in virulence, we also established that GCN4 was necessary for virulence in the mouse model of systemic infection as well as in the Galleria mellonella infection model.

Computationally-obtained structural insights into the molecular interactions between Pidilizumab and binding partners DLL1 and PD-1

Pidilizumab is a monoclonal antibody tested against several types of malignancies, such as lymphoma and metastatic melanoma, showing promising results. In 2016, the FDA put Pidilizumab's clinical studies on partial hold due to emerging evidence pointing to the antibody target uncertainty. Although initial studies indicated an interaction with the PD-1 checkpoint receptor, recent updates assert that Pidilizumab binds primarily to Notch ligand DLL1. However, a detailed description of which interactions coordinate antibody-antigen complex formation is lacking. Therefore, this study uses computational tools to identify molecular interactions between Pidilizumab and its reported targets PD-1 and DLL1. A docking methodology was validated and applied to determine the binding modes between modeled Pidilizumab scFvs and the two antigens. We used Molecular Dynamics (MD) simulations to verify the complexes' stability and submitted the resulting trajectory files to MM/PBSA and Principal Component Analysis. A set of different prediction tools determined scFv interface hot-spots. Whereas docking and MD simulations revealed that the antibody fragments do not interact straightforwardly with PD-1, ten scFv hot-spots, including Met93 and Leu112, mediated the interaction with the DLL1 C2 domain. The interaction triggered a conformational selection-like effect on DLL1, allowing new hydrogen bonds on the β3-β4 interface loop. The unprecedented structural data on Pidilizumab's interactions provided novel evidence that its legitimate target is the DLL1 protein and offered structural insight on how these molecules interact, shedding light on the pathways that could be affected by the use of this essential immunobiological.Communicated by Ramaswamy H. Sarma.

The In Silico Identification of Potential Members of the Ded1/DDX3 Subfamily of DEAD-Box RNA Helicases from the Protozoan Parasite Leishmania infantum and Their Analyses in Yeast

DEAD-box RNA helicases are ubiquitous proteins found in all kingdoms of life and that are associated with all processes involving RNA. Their central roles in biology make these proteins potential targets for therapeutic or prophylactic drugs. The Ded1/DDX3 subfamily of DEAD-box proteins is of particular interest because of their important role(s) in translation. In this paper, we identified and aligned the protein sequences of 28 different DEAD-box proteins from the kinetoplast-protozoan parasite Leishmania infantum, which is the cause of the visceral form of leishmaniasis that is often lethal if left untreated, and compared them with the consensus sequence derived from DEAD-box proteins in general, and from the Ded1/DDX3 subfamily in particular, from a wide variety of other organisms. We identified three potential homologs of the Ded1/DDX3 subfamily and the equivalent proteins from the related protozoan parasite Trypanosoma brucei, which is the causative agent of sleeping sickness. We subsequently tested these proteins for their ability to complement a yeast strain deleted for the essential DED1 gene. We found that the DEAD-box proteins from Trypanosomatids are highly divergent from other eukaryotes, and consequently they are suitable targets for protein-specific drugs.

Structural insights into the tropomodulin assembly at the pointed ends of actin filaments

Tropomodulins are a family of important regulators of actin dynamics at the pointed ends of actin filaments. Four isoforms of tropomodulin, Tmod1-Tmod4, are expressed in vertebrates. Binding of tropomodulin to the pointed end is dependent on tropomyosin, an actin binding protein that itself is represented in mammals by up to 40 isoforms. The understanding of the regulatory role of the tropomodulin/tropomyosin molecular diversity has been limited due to the lack of a three-dimensional structure of the tropomodulin/tropomyosin complex. In this study, we mapped tropomyosin residues interacting with two tropomyosin-binding sites of tropomodulin and generated a three-dimensional model of the tropomodulin/tropomyosin complex for each of these sites. The models were refined by molecular dynamics simulations and validated via building a self-consistent three-dimensional model of tropomodulin assembly at the pointed end. The model of the pointed-end Tmod assembly offers new insights in how Tmod binding ensures tight control over the pointed end dynamics.

Role of an ancient light-harvesting protein of PSI in light absorption and photoprotection

Diverse algae of the red lineage possess chlorophyll a-binding proteins termed LHCR, comprising the PSI light-harvesting system, which represent an ancient antenna form that evolved in red algae and was acquired through secondary endosymbiosis. However, the function and regulation of LHCR complexes remain obscure. Here we describe isolation of a Nannochloropsis oceanica LHCR mutant, named hlr1, which exhibits a greater tolerance to high-light (HL) stress compared to the wild type. We show that increased tolerance to HL of the mutant can be attributed to alterations in PSI, making it less prone to ROS production, thereby limiting oxidative damage and favoring growth in HL. HLR1 deficiency attenuates PSI light-harvesting capacity and growth of the mutant under light-limiting conditions. We conclude that HLR1, a member of a conserved and broadly distributed clade of LHCR proteins, plays a pivotal role in a dynamic balancing act between photoprotection and efficient light harvesting for photosynthesis.

LHCR proteins are ancient chlorophyll a-binding antennas that evolved in diverse algae of the red lineage. Here Lu et al. characterize a red lineage LHCR mutant and show reduced oxidative damage in high light but attenuated growth under low light, thus demonstrating how LHCR proteins impact the balance between photoprotection and light harvesting.

On the strong connection between nanoscale adhesion of Yad fimbriae and macroscale attachment of Yad-decorated bacteria to glycosylated, hydrophobic and hydrophilic surfaces

Yad fimbriae are currently viewed as versatile bacterial adhesins able to significantly mediate host or plant-pathogen recognition and contribute to the persistence of Escherichia coli in both the environment and within hosts. To date, however, the underlying adhesion process of Yad fimbriae on surfaces defined by controlled coating chemistries has not been evaluated on the relevant molecular scale. In this work, the interaction forces operational between Yad fimbriae expressed by genetically modified E. coli and self-assembled monolayers (SAM) differing in terms of charge, hydrophobicity or the nature of decorating sugar units are quantified by Single Molecule Force Spectroscopy (SMFS) on the nanoscale. It is found that the adhesion of Yad fimbriae onto probes functionalized with xylose is as strong as that measured with probes decorated with anti-Yad antibodies (ca. 80 to 300 pN). In contrast, the interactions of Yad with galactose, lactose, mannose, -OH, -NH₂, -COOH and -CH₃ terminated SAMs are clearly non-specific. Interpretation of SMFS measurements on the basis of worm-like-chain modeling for polypeptide nanomechanics further leads to the estimates of the maximal extension of Yad fimbriae upon stretching, of their persistence length and of their polydispersity. Finally, we show for the first time a strong correlation between the adhesion properties of Yad-decorated bacteria determined from conventional macroscopic counting methods and the molecular adhesion capacity of Yad fimbriae. This demonstration advocates the effort that should be made to understand on the nanoscale level the interactions between fimbriae and their cognate ligands. The results could further help the design of potential anti-adhesive molecules or surfaces to better fight against the virulence of bacterial pathogens.

DescribePROT: database of amino acid-level protein structure and function predictions

We present DescribePROT, the database of predicted amino acid-level descriptors of structure and function of proteins. DescribePROT delivers a comprehensive collection of 13 complementary descriptors predicted using 10 popular and accurate algorithms for 83 complete proteomes that cover key model organisms. The current version includes 7.8 billion predictions for close to 600 million amino acids in 1.4 million proteins. The descriptors encompass sequence conservation, position specific scoring matrix, secondary structure, solvent accessibility, intrinsic disorder, disordered linkers, signal peptides, MoRFs and interactions with proteins, DNA and RNAs. Users can search DescribePROT by the amino acid sequence and the UniProt accession number and entry name. The pre-computed results are made available instantaneously. The predictions can be accesses via an interactive graphical interface that allows simultaneous analysis of multiple descriptors and can be also downloaded in structured formats at the protein, proteome and whole database scale. The putative annotations included by DescriPROT are useful for a broad range of studies, including: investigations of protein function, applied projects focusing on therapeutics and diseases, and in the development of predictors for other protein sequence descriptors. Future releases will expand the coverage of DescribePROT. DescribePROT can be accessed at http://biomine.cs.vcu.edu/servers/DESCRIBEPROT/.

Nuclear and Nucleolar Localization of Bovine Adenovirus-3 Protein V

The L2 region of bovine adenovirus-3 (BAdV-3) encodes a Mastadenovirus genus-specific protein, designated as pV, which is important for the production of progeny viruses. Here, we demonstrate that BAdV-3 pV, expressed as 55 kDa protein, localizes to the nucleus and specifically targets nucleolus of the infected cells. Analysis of deletion mutants of pV suggested that amino acids 81–120, 190–210, and 380–389 act as multiple nuclear localization signals (NLS), which also appear to serve as the binding sites for importin α-3 protein, a member of the importin α/β nuclear import receptor pathway. Moreover, pV amino acids 21–50 and 380–389 appear to act as nucleolar localization signals (NoLs). Interestingly, amino acids 380–389 appear to act both as NLS and as NoLS. The presence of NoLS is essential for the production of infectious progeny virions, as deletion of both NoLs are lethal for the production of infectious BAdV-3. Analysis of mutant BAV.pVd1d3 (isolated in pV completing CRL cells) containing deletion/mutation of both NoLS in non-complementing CRL cells not only revealed the altered intracellular localization of mutant pV but also reduced the expression of some late proteins. However, it does not appear to affect the incorporation of viral proteins, including mutant pV, in BAV.pVd1d3 virions. Further analysis of CsCl purified BAV.pVd1d3 suggested the presence of thermo-labile virions with disrupted capsids, which appear to affect the infectivity of the progeny virions. Our results suggest that pV contains overlapping and non-overlapping NoLS/NLS. Moreover, the presence of both NoLS appear essential for the production of stable and infectious progeny BAV.pVd1d3 virions.

Protein Subcellular Localization Prediction

The elucidation of the subcellular localization of proteins is very important in order to deeply understand their functions. In fact, proteins activities are strictly correlated to the cellular compartment and microenvironment in which they are present.In recent years, several effective and reliable proteomics techniques and computational methods have been developed and implemented in order to identify the proteins subcellular localization. This process is often time-consuming and expensive, but the recent technological and bioinformatics progress allowed the development of more accurate and simple workflows to determine the localization, interactions, and functions of proteins.In the following chapter, a brief introduction on the importance of knowing subcellular localization of proteins will be presented. Then, sample preparation protocols, proteomic methods, data analysis strategies, and software for the prediction of proteins localization will be presented and discussed. Finally, the more recent and advanced spatial proteomics techniques will be shown.

Computational Methods for the Elucidation of Protein Structure and Interactions

Biologists are increasingly aware of the importance of protein structure in revealing function. The computational tools now exist which allow researchers to model unknown proteins simply on the basis of their primary sequence. However, for the non-specialist bioinformatician, there is a dazzling array of terminology, acronyms, and competing computer software available for this process. This review is intended to highlight the key stages of computational protein structure prediction, as well as explain the reasons behind some of the procedures and list some established workarounds for common pitfalls. Thereafter follows a review of five one-stop servers for start-to-finish structure prediction.

Computational insight into the three-dimensional structure of ADP ribosylation factor like protein 15, a novel susceptibility gene for rheumatoid arthritis

The ARL15 gene (ADP ribosylation factor like protein 15) encodes for an uncharacterized small GTP-binding protein. Its exact role in human physiology remains unknown, but a number of genetic association studies have recognised different variants in this gene to be statistically associated with numerous traits and complex diseases. We have previously reported a novel association of ARL15 with rheumatoid arthritis (RA) based on a genome-wide association study in a north Indian cohort. Subsequent investigations have provided leads for its involvement in RA pathophysiology, especially its potential as a novel therapeutic target. However, the absence of an experimentally determined tertiary structure for ARL15 significantly hinders the understanding of its biochemical and physiological functions, as well as development of potential lead molecules. We, therefore, aimed to derive a high quality, refined model of the three dimensional structure of human ARL15 protein using two different computational protein structure prediction methods - template-based threading and ab initio modelling. The best model each from among the five each derived from both the approaches was selected based on stringent quality assessment and refinement. Molecular dynamics simulations over long timescales revealed the ab initio model to be relatively more stable, and it marginally outperformed the template-based model in the quality assessment as well. A putative GTP-binding site was also predicted using homology for the ARL15 protein, where potential competitive inhibitors can be targeted. This high quality predicted model may provide insights to the biological role(s) of ARL15 and inform and guide further experimental, structural and biochemical characterization efforts.Communicated by Ramaswamy H. Sarma.

Conformation of the Intermediates in the Reaction Catalyzed by Protoporphyrinogen Oxidase: An In Silico Analysis

Protoporphyrinogen oxidase (PPO) is a critical enzyme across life as the last common step in the synthesis of many metalloporphyrins. The reaction mechanism of PPO was assessed in silico and the unstructured loop near the binding pocket was investigated. The substrate, intermediates, and product were docked in the catalytic domain of PPO using a modified Autodock method, introducing flexibility in the macrocycles. Sixteen PPO protein sequences across phyla were aligned and analyzed with Phyre2 and ProteinPredict to study the unstructured loop from residue 204–210 in the H. sapiens structure. Docking of the substrate, intermediates, and product all resulted in negative binding energies, though the substrate had a lower energy than the others by 40%. The α-H of C10 was found to be 1.4 angstroms closer to FAD than the β-H, explaining previous reports of the reaction occurring on the meso face of the substrate. A lack of homology in sequence or length in the unstructured loop indicates a lack of function for the protein reaction. This docking study supports a reaction mechanism proposed previously whereby all hydride abstractions occur on the C10 of the tetrapyrrole followed by tautomeric rearrangement to prepare the intermediate for the next reaction.

ARL8 Relieves SKIP Autoinhibition to Enable Coupling of Lysosomes to Kinesin-1

Long-range movement of organelles within the cytoplasm relies on coupling to microtubule motors, a process that is often mediated by adaptor proteins. In many cases, this coupling involves organelle- or adaptor-induced activation of the microtubule motors by conformational reversal of an autoinhibited state. Herein, we show that a similar regulatory mechanism operates for an adaptor protein named SKIP (also known as PLEKHM2). SKIP binds to the small guanosine triphosphatase (GTPase) ARL8 on the lysosomal membrane to couple lysosomes to the anterograde microtubule motor kinesin-1. Structure-function analyses of SKIP reveal that the C-terminal region comprising three pleckstrin homology (PH) domains interacts with the N-terminal region comprising ARL8- and kinesin-1-binding sites. This interaction inhibits coupling of lysosomes to kinesin-1 and, consequently, lysosome movement toward the cell periphery. We also find that ARL8 does not just recruit SKIP to the lysosomal membrane but also relieves SKIP autoinhibition, promoting kinesin-1-driven, anterograde lysosome transport. Finally, our analyses show that the largely disordered middle region of SKIP mediates self-association and that this self-association enhances the interaction of SKIP with kinesin-1. These findings indicate that SKIP is not just a passive connector of lysosome-bound ARL8 to kinesin-1 but is itself subject to intra- and inter-molecular interactions that regulate its function. We anticipate that similar organelle- or GTPase-induced conformational changes could regulate the activity of other kinesin adaptors.

Evolution of gentamicin and arsenite resistance acquisition in Ralstonia pickettii water isolates

Ralstonia pickettii are ubiquitous in water environments. Members of this species are frequently, but not always, resistant to both gentamicin and arsenite. Gentamicin and arsenite co-resistance and the putative molecular mechanisms were investigated. A group of 37 R. pickettii strains isolated from drinking water and hospital wastewater were characterized for gentamicin and arsenite resistance phenotypes, the number and size of plasmids, and screened for genetic elements associated with arsenite tolerance, Integrative and Conjugative Elements (ICEs), among other. The genomes of three representative strains were compared. Most gentamicin resistant (GR) isolates (32/33) were resistant to arsenite, and harbored ICE- and ars operon-related genes. These genetic elements were not detected in any of the five arsenite susceptible strains, regardless of the GR (n = 1) or gentamicin susceptibility (GS) (n = 4) phenotype. The comparison of the genomes of two GR (one resistant and one susceptible to arsenite) and one GS strains suggested that these phenotypes correspond to three phylogroups, distinguished by presence of some genes only in GR isolates, in addition to point mutations in functional genes. The presence of ICEs and ars operon-related genes suggest that arsenite resistance might have been acquired by GR lineages.

Vaccine design based on 16 epitopes of SARS-CoV-2 spike protein

The global outbreak of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) urgently requires an effective vaccine for prevention. In this study, 66 epitopes containing pentapeptides of SARS‐CoV‐2 spike protein in the IEDB database were compared with the amino acid sequence of SARS‐CoV‐2 spike protein, and 66 potentially immune‐related peptides of SARS‐CoV‐2 spike protein were obtained. Based on the single‐nucleotide polymorphisms analysis of spike protein of 1218 SARS‐CoV‐2 isolates, 52 easily mutated sites were identified and used for vaccine epitope screening. The best vaccine candidate epitopes in the 66 peptides of SARS‐CoV‐2 spike protein were screened out through mutation and immunoinformatics analysis. The best candidate epitopes were connected by different linkers in silico to obtain vaccine candidate sequences. The results showed that 16 epitopes were relatively conservative, immunological, nontoxic, and nonallergenic, could induce the secretion of cytokines, and were more likely to be exposed on the surface of the spike protein. They were both B‐ and T‐cell epitopes, and could recognize a certain number of HLA molecules and had high coverage rates in different populations. Moreover, epitopes 897‐913 were predicted to have possible cross‐immunoprotection for SARS‐CoV and SARS‐CoV‐2. The results of vaccine candidate sequences screening suggested that sequences (without linker, with linker GGGSGGG, EAAAK, GPGPG, and KK, respectively) were the best. The proteins translated by these sequences were relatively stable, with a high antigenic index and good biological activity. Our study provided vaccine candidate epitopes and sequences for the research of the SARS‐CoV‐2 vaccine.

Xanthomonas axonopodis pv. punicae depends on multiple non-TAL (Xop) T3SS effectors for its coveted growth inside the pomegranate plant through repressing the immune responses during bacterial blight development

Xanthomonas axonopodis pv. punicae (Xap), the bacterial blight pathogen of pomegranate, incurs substantial loss to yield and reduces export quality of this economically important fruit crop. During infection, the bacterium secretes six non-TAL (Xop) effectors into the pomegranate cells through a specialized type three secretion system (T3SS). Previously, we demonstrated the role of two key effectors, XopL and XopN in pathogenesis. Here, we investigate the role of rest effectors (XopC2, XopE1, XopQ and XopZ) on disease development. We generated null mutants for each individual effector and mutant bacterial suspension was infiltrated into pomegranate leaves. Compared to Xap wild, the mutant bacterial growth was reduced by 2.7-11.5 folds. The mutants produced lesser water-soaked lesions when infiltrated on leaves by 1.13-2.21 folds. Among the four effectors, XopC2 contributes highest for in planta bacterial growth and disease development. XopC2 efficiently suppressed the defense responses like callose deposition, reactive oxygen species (ROS) and the activation of immune responsive genes. Being a major contributor, we further characterize XopC2 for its subcellular localization, its protein structure and networking. XopC2 is localized to the plasma membrane of Nicotiana benthamiana like XopL and XopN. XopC2 is a 661 amino acids protein having 15 alpha and 17 beta helix. Our STRING and I-TASSER based analysis hinted that XopC2 interacts with multiple membrane localized plant proteins including transcription regulator of CCR4-NOT family, TTN of maintenance of chromosome family and serine/threonine-protein phosphatase 2A (PP2A) isoform. Based on the interaction it is predicted that XopC2 might involve in diverse functions like nuclear-transcribed mRNA catabolic process, maintenance of chromosome, hormone signaling and protein dephosphorylation activities and thereby suppress the plant immunity. Altogether, our study suggests that Xap largely depends on three non-TAL (Xop) effectors, including XopC2, XopL and XopN, to modulate pomegranate PTI for its unrestricted proliferation during bacterial blight development.

PssJ Is a Terminal Galactosyltransferase Involved in the Assembly of the Exopolysaccharide Subunit in Rhizobium Leguminosarum bv. Trifolii

Rhizobium leguminosarum bv. trifolii produces exopolysaccharide (EPS) composed of glucose, glucuronic acid, and galactose residues at a molar ratio 5:2:1. A majority of genes involved in the synthesis, modification, and export of exopolysaccharide are located in the chromosomal Pss-I region. In the present study, a ΔpssJ deletion mutant was constructed and shown to produce EPS lacking terminal galactose in the side chain of the octasaccharide subunit. The lack of galactose did not block EPS subunit translocation and polymerization. The in trans delivery of the pssJ gene restored the production of galactose-containing exopolysaccharide. The mutant was compromised in several physiological traits, e.g., motility and biofilm production. An impact of the pssJ mutation and changed EPS structure on the symbiotic performance was observed as improper signaling at the stage of molecular recognition, leading to formation of a significant number of non-infected empty nodules. Terminal galactosyltransferase PssJ was shown to display a structure typical for the GT-A class of glycosyltransferases and interact with other GTs and Wzx/Wzy system proteins. The latter, together with PssJ presence in soluble and membrane protein fractions indicated that the protein plays its role at the inner membrane interface and as a component of a larger complex.

EVI1 phosphorylation at S436 regulates interactions with CtBP1 and DNMT3A and promotes self-renewal

The transcriptional regulator EVI1 has an essential role in early development and haematopoiesis. However, acute myeloid leukaemia (AML) driven by aberrantly high EVI1 expression has very poor prognosis. To investigate the effects of post-translational modifications on EVI1 function, we carried out a mass spectrometry (MS) analysis of EVI1 in AML and detected dynamic phosphorylation at serine 436 (S436). Wild-type EVI1 (EVI1-WT) with S436 available for phosphorylation, but not non-phosphorylatable EVI1-S436A, conferred haematopoietic progenitor cell self-renewal and was associated with significantly higher organised transcriptional patterns. In silico modelling of EVI1-S436 phosphorylation showed reduced affinity to CtBP1, and CtBP1 showed reduced interaction with EVI1-WT compared with EVI1-S436A. The motif harbouring S436 is a target of CDK2 and CDK3 kinases, which interacted with EVI1-WT. The methyltransferase DNMT3A bound preferentially to EVI1-WT compared with EVI1-S436A, and a hypomethylated cell population associated by EVI1-WT expression in murine haematopoietic progenitors is not maintained with EVI1-S436A. These data point to EVI1-S436 phosphorylation directing functional protein interactions for haematopoietic self-renewal. Targeting EVI1-S436 phosphorylation may be of therapeutic benefit when treating EVI1-driven leukaemia.

Functional Distinction between Human and Mouse Sodium-Coupled Citrate Transporters and Its Biologic Significance: An Attempt for Structural Basis Using a Homology Modeling Approach

NaCT (SLC13A5; mINDY), a sodium-coupled citrate transporter, is the mammalian ortholog of Drosophila INDY. Loss-of-function mutations in human NaCT cause severe complications with neonatal epilepsy and encephalopathy (EIEE25). Surprisingly, mice lacking this transporter do not have this detrimental brain phenotype. The marked differences in transport kinetics between mouse and human NaCTs provide at least a partial explanation for this conundrum, but a structural basis for the differences is lacking. Neither human nor mouse NaCT has been crystallized, and any information known on their structures is based entirely on what was inferred from the structure of VcINDY, a related transporter in bacteria. Here, we highlight the functional features of human and mouse NaCTs and provide a plausible molecular basis for the differences based on a full-length homology modeling approach. The transport characteristics of human NaCT markedly differ from those of VcINDY. Therefore, the modeling with VcINDY as the template is flawed, but this is the best available option at this time. With the newly deduced model, we determined the likely locations of the disease-causing mutations and propose a new classification for the mutations based on their location and potential impact on transport function. This new information should pave the way for future design and development of novel therapeutics to restore the lost function of the mutant transporters as a treatment strategy for patients with EIEE25.

Functional insights from biophysical study of TREM2 interactions with apoE and Aβ1-42

INTRODUCTION

Triggering receptor expressed on myeloid cells-2 (TREM2) is an immune receptor expressed on microglia that also can become soluble (sTREM2). How TREM2 engages different ligands remains poorly understood.

METHODS

We used comprehensive biolayer interferometry (BLI) analysis to investigate TREM2 and sTREM2 interactions with apolipoprotein E (apoE) and monomeric amyloid beta (Aβ) (mAβ42).

RESULTS

TREM2 engagement of apoE was protein mediated with little effect of lipidation, showing slight affinity differences between isoforms (E4 > E3 > E2). Another family member, TREML2, did not bind apoE. Disease-linked TREM2 variants within a "basic patch" minimally impact apoE binding. Instead, TREM2 uses a unique hydrophobic surface to bind apoE, which requires the apoE hinge region. TREM2 and sTREM2 directly bind mAβ42 and potently inhibit Aβ42 polymerization, suggesting a potential role for soluble sTREM2 in preventing AD pathogenesis.

DISCUSSION

These findings demonstrate that TREM2 has at least two ligand-binding surfaces that might be therapeutic targets and uncovers a potential function for sTREM2 in directly inhibiting Aβ polymerization.

Identification of the nuclear and nucleolar localization signals of the Feline immunodeficiency virus Rev protein

Lentivirus genomes code for a regulatory protein essential for virus replication termed Rev. The Rev protein binds to partially spliced and unspliced viral RNAs and mediates their nuclear export. Therefore, Rev possesses functional domains that enable its shuttling between the cytoplasm and the nucleus. The Feline immunodeficiency virus (FIV), a lentivirus, can lead to an immunodeficiency syndrome after a long incubation period, similar to that associated with the human immunodeficiency virus type 1 (HIV-1). The FIV Rev functional domains have been predicted only by homology with those of HIV-1 Rev. In the present study, the nuclear and nucleolar localization signals (NLS and NoLS, respectively) of the FIV Rev were examined. A series of FIV Rev deletion mutants fused to the enhanced green fluorescent protein (EGFP) were used to localize the NLS in a region spanning amino acids (aa) 81-100. By using alanine substitution mutants, basic residues present between the amino acids (aa) 84-99 of the FIV Rev protein sequence were identified to form the NLS, whereas those between aa 82-95 were associated with the NoLS function. These results further enhance our understanding of how Rev exerts its role in the replication cycle of lentiviruses.

Statistical and molecular dynamics (MD) simulation approach to investigate the role of intrinsically disordered regions of shikimate dehydrogenase in microorganisms surviving at different temperatures

Hyperthermophiles, a subset of prokaryotes that thrive in adverse temperatures, potentially utilize the protein molecular biosystem for maintaining thermostability in a wide range of temperatures. Recent studies revealed that these organisms have smaller proportions of intrinsically disordered proteins. In this study, we performed sequence and structural analysis to investigate the maintenance of protein conformation and their stability at different temperatures. The sequence analysis reveals the higher proportion of charged amino acids are responsible for preventing the helix formation and, hence, become disordered regions. For structural analysis, we chose shikimate dehydrogenase from four species, namely Listeria monocytogenes, Escherichia coli, Thermus thermophilus, and Methanopyrus kandleri, and evaluated the protein adaptation at 283 K, 300 K, and 395 K temperatures. From this investigation, we found more residues of shikimate dehydrogenase prefer an order-to-disorder transition at 395 K only for hyperthermophilic species. The solvent-accessible surface area (SASA) and hydrogen-bond analysis revealed that the tertiary conformation and the number of hydrogen bonds for hyperthermophilic shikimate dehydrogenase are highly preserved at 395 K, compared to 300 K. Our simulation results conjointly provide shikimate dehydrogenase of hyperthermophile which resists high temperatures through stronger protein tertiary conformations.

The Splicing Factor hnRNP M Is a Critical Regulator of Innate Immune Gene Expression in Macrophages

While transcriptional control of innate immune gene expression is well characterized, almost nothing is known about how pre-mRNA splicing decisions influence, or are influenced by, macrophage activation. Here, we demonstrate that the splicing factor hnRNP M is a critical repressor of innate immune gene expression and that its function is regulated by pathogen sensing cascades. Loss of hnRNP M led to hyperinduction of a unique regulon of inflammatory and antimicrobial genes following diverse innate immune stimuli. While mutating specific serines on hnRNP M had little effect on its ability to control pre-mRNA splicing or transcript levels of housekeeping genes in resting macrophages, it greatly impacted the protein’s ability to dampen induction of specific innate immune transcripts following pathogen sensing. These data reveal a previously unappreciated role for pattern recognition receptor signaling in controlling splicing factor phosphorylation and establish pre-mRNA splicing as a critical regulatory node in defining innate immune outcomes.

West et al. report that hnRNP M represses expression of a cohort of innate immune transcripts in infected macrophages. IL6 splicing repression is relieved when hnRNP M is phosphorylated at specific residues, demonstrating that post-translational modification of splicing factors downstream of pathogen sensing can control maturation of innate immune mRNAs.

Atypical immune phenotype in severe combined immunodeficiency patients with novel mutations in IL2RG and JAK3

Mutations in the common gamma chain of the interleukin 2 receptor (IL2RG) or the associated downstream signaling enzyme Janus kinase 3 (JAK3) genes are typically characterized by a T cell-negative, B cell-positive, natural killer (NK) cell-negative (T^-B⁺NK^-) severe combined immunodeficiency (SCID) immune phenotype. We report clinical course, immunological, genetic and proteomic work-up of two patients with different novel mutations in the IL-2-JAK3 pathway with a rare atypical presentation of T^-B⁺NK^- SCID. Lymphocyte subpopulation revealed significant T cells lymphopenia, normal B cells, and NK cells counts (T-B+NK+SCID). Despite the presence of B cells, IgG levels were low and IgA and IgM levels were undetectable. T-cell proliferation in response to mitogens in patient 1 was very low and T-cell receptor V-beta chain repertoire in patient 2 was polyclonal. Whole-exome sequencing revealed novel mutations in both patients (patient 1-c.923delC frame-shift mutation in the IL2RG gene, patient 2-c.G172A a homozygous missense mutation in the JAK3 gene). Bioinformatic analysis of the JAK3 mutation indicated deleterious effect and 3D protein modeling located the mutation to a surface exposed alpha-helix structure. Our findings help to link between genotype and phenotype, which is a key factor for the diagnosis and treatment of SCID patients.

Phase Transitions by an Abundant Protein in the Anammox Extracellular Matrix Mediate Cell-to-Cell Aggregation and Biofilm Formation

By employing biophysical and liquid-liquid phase separation concepts, this study revealed how a highly abundant extracellular protein enhances the key environmental and industrial bioprocess of anaerobic ammonium oxidation (anammox). Extracellular proteins of environmental biofilms are understudied and poorly annotated in public databases. Understanding the function of extracellular proteins is also increasingly important for improving bioprocesses and resource recovery. Here, protein functions were assessed based on theoretical predictions of intrinsically disordered domains, known to promote adhesion and liquid-liquid phase separation, and available surface layer protein properties. A model is thus proposed to explain how the protein promotes aggregation and biofilm formation by extracellular matrix remodeling and phase transitions. This work provides a strong foundation for functional investigations of extracellular proteins involved in biofilm development.

This study describes the first direct functional assignment of a highly abundant extracellular protein from a key environmental and biotechnological biofilm performing an anaerobic ammonium oxidation (anammox) process. Expression levels of Brosi_A1236, belonging to a class of proteins previously suggested to be cell surface associated, were in the top one percentile of all genes in the “Candidatus Brocadia sinica”-enriched biofilm. The Brosi_A1236 structure was computationally predicted to consist of immunoglobulin-like anti-parallel β-strands, and circular dichroism conducted on the isolated surface protein indicated that β-strands are the dominant higher-order structure. The isolated protein was stained positively by the β-sheet-specific stain thioflavin T, along with cell surface- and matrix-associated regions of the biofilm. The surface protein has a large unstructured content, including two highly disordered domains at its C terminus. The disordered domains bound to the substratum and thereby facilitated the adhesion of negatively charged latex microspheres, which were used as a proxy for cells. The disordered domains and isolated whole surface protein also underwent liquid-liquid phase separation to form liquid droplets in suspension. Liquid droplets of disordered protein wet the surfaces of microspheres and bacterial cells and facilitated their coalescence. Furthermore, the surface layer protein formed gels as well as ordered crystalline structures. These observations suggest that biophysical remodeling through phase transitions promotes aggregation and biofilm formation.

Leiomodin creates a leaky cap at the pointed end of actin-thin filaments

Improper lengths of actin-thin filaments are associated with altered contractile activity and lethal myopathies. Leiomodin, a member of the tropomodulin family of proteins, is critical in thin filament assembly and maintenance; however, its role is under dispute. Using nuclear magnetic resonance data and molecular dynamics simulations, we generated the first atomic structural model of the binding interface between the tropomyosin-binding site of cardiac leiomodin and the N-terminus of striated muscle tropomyosin. Our structural data indicate that the leiomodin/tropomyosin complex only forms at the pointed end of thin filaments, where the tropomyosin N-terminus is not blocked by an adjacent tropomyosin protomer. This discovery provides evidence supporting the debated mechanism where leiomodin and tropomodulin regulate thin filament lengths by competing for thin filament binding. Data from experiments performed in cardiomyocytes provide additional support for the competition model; specifically, expression of a leiomodin mutant that is unable to interact with tropomyosin fails to displace tropomodulin at thin filament pointed ends and fails to elongate thin filaments. Together with previous structural and biochemical data, we now propose a molecular mechanism of actin polymerization at the pointed end in the presence of bound leiomodin. In the proposed model, the N-terminal actin-binding site of leiomodin can act as a "swinging gate" allowing limited actin polymerization, thus making leiomodin a leaky pointed-end cap. Results presented in this work answer long-standing questions about the role of leiomodin in thin filament length regulation and maintenance.

Slr4, a newly identified S-layer protein from marine Gammaproteobacteria, is a major biofilm matrix component

S‐layers are paracrystalline proteinaceous lattices that surround prokaryotic cells, forming a critical interface between the cells and their extracellular environment. Here, we report the discovery of a novel S‐layer protein present in the Gram‐negative marine organism, Pseudoalteromonas tunicata D2. An uncharacterized protein (EAR28894) was identified as the most abundant protein in planktonic cultures and biofilms. Bioinformatic methods predicted a beta‐helical structure for EAR28894 similar to the Caulobacter S‐layer protein, RsaA, despite sharing less than 20% sequence identity. Transmission electron microscopy revealed that purified EAR28894 protein assembled into paracrystalline sheets with a unique square lattice symmetry and a unit cell spacing of ~9.1 nm. An S‐layer was found surrounding the outer membrane in wild‐type cells and completely removed from cells in an EAR28894 deletion mutant. S‐layer material also appeared to be “shed” from wild‐type cells and was highly abundant in the extracellular matrix where it is associated with outer membrane vesicles and other matrix components. EAR28894 and its homologs form a new family of S‐layer proteins that are widely distributed in Gammaproteobacteria including species of Pseudoalteromonas and Vibrio, and found exclusively in marine metagenomes. We propose the name Slr4 for this novel protein family.

The ancillary N-terminal region of the yeast AP-1 transcription factor Yap8 contributes to its DNA binding specificity

Activator protein 1 (AP-1) is one of the largest families of basic leucine zipper (bZIP) transcription factors in eukaryotic cells. How AP-1 proteins achieve target DNA binding specificity remains elusive. In Saccharomyces cerevisiae, the AP-1-like protein (Yap) family comprises eight members (Yap1 to Yap8) that display distinct genomic target sites despite high sequence homology of their DNA binding bZIP domains. In contrast to the other members of the Yap family, which preferentially bind to short (7–8 bp) DNA motifs, Yap8 binds to an unusually long DNA motif (13 bp). It has been unclear what determines this unique specificity of Yap8. In this work, we use molecular and biochemical analyses combined with computer-based structural design and molecular dynamics simulations of Yap8–DNA interactions to better understand the structural basis of DNA binding specificity determinants. We identify specific residues in the N-terminal tail preceding the basic region, which define stable association of Yap8 with its target promoter. We propose that the N-terminal tail directly interacts with DNA and stabilizes Yap8 binding to the 13 bp motif. Thus, beside the core basic region, the adjacent N-terminal region contributes to alternative DNA binding selectivity within the AP-1 family.

FLAVi: An Enhanced Annotator for Viral Genomes of Flaviviridae

Responding to the ongoing and severe public health threat of viruses of the family Flaviviridae, including dengue, hepatitis C, West Nile, yellow fever, and Zika, demands a greater understanding of how these viruses emerge and spread. Updated phylogenies are central to this understanding. Most cladograms of Flaviviridae focus on specific lineages and ignore outgroups, hampering the efficacy of the analysis to test ingroup monophyly and relationships. This is due to the lack of annotated Flaviviridae genomes, which has gene content variation among genera. This variation makes analysis without partitioning difficult. Therefore, we developed an annotation pipeline for the genera of Flaviviridae (Flavirirus, Hepacivirus, Pegivirus, and Pestivirus, named "Fast Loci Annotation of Viruses" (FLAVi; http://flavi-web.com/), that combines ab initio and homology-based strategies. FLAVi recovered 100% of the genes in Flavivirus and Hepacivirus genomes. In Pegivirus and Pestivirus, annotation efficiency was 100% except for one partition each. There were no false positives. The combined phylogenetic analysis of multiple genes made possible by annotation has clear impacts over the tree topology compared to phylogenies that we inferred without outgroups or data partitioning. The final tree is largely congruent with previous hypotheses and adds evidence supporting the close phylogenetic relationship between dengue and Zika.

Mutation of CFAP57, a protein required for the asymmetric targeting of a subset of inner dynein arms in Chlamydomonas, causes primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is characterized by chronic airway disease, reduced fertility, and randomization of the left/right body axis. It is caused by defects of motile cilia and sperm flagella. We screened a cohort of affected individuals that lack an obvious axonemal defect for pathogenic variants using whole exome capture, next generation sequencing, and bioinformatic analysis assuming an autosomal recessive trait. We identified one subject with an apparently homozygous nonsense variant [(c.1762C>T), p.(Arg588*)] in the uncharacterized CFAP57 gene. Interestingly, the variant results in the skipping of exon 11 (58 amino acids), which may be due to disruption of an exonic splicing enhancer. In normal human nasal epithelial cells, CFAP57 localizes throughout the ciliary axoneme. Nasal cells from the PCD patient express a shorter, mutant version of CFAP57 and the protein is not incorporated into the axoneme. The missing 58 amino acids include portions of WD repeats that may be important for loading onto the intraflagellar transport (IFT) complexes for transport or docking onto the axoneme. A reduced beat frequency and an alteration in ciliary waveform was observed. Knockdown of CFAP57 in human tracheobronchial epithelial cells (hTECs) recapitulates these findings. Phylogenetic analysis showed that CFAP57 is highly conserved in organisms that assemble motile cilia. CFAP57 is allelic with the BOP2/IDA8/FAP57 gene identified previously in Chlamydomonas reinhardtii. Two independent, insertional fap57 Chlamydomonas mutant strains show reduced swimming velocity and altered waveforms. Tandem mass tag (TMT) mass spectroscopy shows that FAP57 is missing, and the "g" inner dyneins (DHC7 and DHC3) and the "d" inner dynein (DHC2) are reduced, but the FAP57 paralog FBB7 is increased. Together, our data identify a homozygous variant in CFAP57 that causes PCD that is likely due to a defect in the inner dynein arm assembly process.

Apcdd1 is a dual BMP/Wnt inhibitor in the developing nervous system and skin

Animal development and homeostasis depend on precise temporal and spatial intercellular signaling. Components shared between signaling pathways, generally thought to decrease specificity, paradoxically can also provide a solution to pathway coordination. Here we show that the Bone Morphogenetic Protein (BMP) and Wnt signaling pathways share Apcdd1 as a common inhibitor and that Apcdd1 is a taxon-restricted gene with novel domains and signaling functions. Previously, we showed that Apcdd1 inhibits Wnt signaling (Shimomura et al., 2010), here we find that Apcdd1 potently inhibits BMP signaling in body axis formation and neural differentiation in chicken, frog, zebrafish. Furthermore, we find that Apcdd1 has an evolutionarily novel protein domain. Our results from experiments and modeling suggest that Apcdd1 may coordinate the outputs of two signaling pathways that are central to animal development and human disease.

Genetic Variation of Multiple Serotypes of Enteroviruses Associated with Hand, Foot and Mouth Disease in Southern China

Enteroviruses (EVs) species A are a major public health issue in the Asia–Pacific region and cause frequent epidemics of hand, foot and mouth disease (HFMD) in China. Mild infections are common in children; however, HFMD can also cause severe illness that affects the central nervous system. To molecularly characterize EVs, a prospective HFMD virological surveillance program was performed in China between 2013 and 2016. Throat swabs, rectal swabs and stool samples were collected from suspected HFMD patients at participating hospitals. EVs were detected using generic real-time and nested reverse transcription-polymerase chain reactions (RT-PCRs). Then, the complete VP1 regions of enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16) and CVA6 were sequenced to analyze amino acid changes and construct a viral molecular phylogeny. Of the 2836 enrolled HFMD patients, 2,517 (89%) were EV positive. The most frequently detected EVs were CVA16 (32.5%, 819), CVA6 (31.2%, 785), and EV-A71 (20.4%, 514). The subgenogroups CVA16_B1b, CVA6_D3a and EV-A71_C4a were predominant in China and recombination was not observed in the VP1 region. Sequence analysis revealed amino acid variations at the 30, 29 and 44 positions in the VP1 region of EV-A71, CVA16 and CVA6 (compared to the respective prototype strains BrCr, G10 and Gdula), respectively. Furthermore, in 21 of 24 (87.5%) identified EV-A71 samples, a known amino acid substitution (D31N) that may enhance neurovirulence was detected. Our study provides insights about the genetic characteristics of common HFMD-associated EVs. However, the emergence and virulence of the described mutations require further investigation.

Structural basis of UCUU RNA motif recognition by splicing factor RBM20

The vertebrate splicing factor RBM20 (RNA binding motif protein 20) regulates protein isoforms important for heart development and function, with mutations in the gene linked to cardiomyopathy. Previous studies have identified the four nucleotide RNA motif UCUU as a common element in pre-mRNA targeted by RBM20. Here, we have determined the structure of the RNA Recognition Motif (RRM) domain from mouse RBM20 bound to RNA containing a UCUU sequence. The atomic details show that the RRM domain spans a larger region than initially proposed in order to interact with the complete UCUU motif, with a well-folded C-terminal helix encoded by exon 8 critical for high affinity binding. This helix only forms upon binding RNA with the final uracil, and removing the helix reduces affinity as well as specificity. We therefore find that RBM20 uses a coupled folding-binding mechanism by the C-terminal helix to specifically recognize the UCUU RNA motif.

Design of engineered surfaces for prospective detection of SARS-CoV-2 using quartz crystal microbalance-based techniques

INTRODUCTION

Rapid transmission of the severe acute respiratory syndrome coronavirus 2 has affected the whole world and forced it to a halt (lockdown). A fast and label-free detection method for the novel coronavirus needs to be developed along with the existing enzyme-linked immunosorbent assay (ELISA) and reverse transcription polymerase chain reaction (RT-PCR)-based methods.

AREAS COVERED

In this report, biophysical aspects of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein are outlined based on its recent reported electron microscopy structure. Protein binding sites are analyzed theoretically, which consisted of hydrophobic and positive charged amino acid residues. Different strategies to form mixed self-assembled monolayers (SAMs) of hydrophobic (CH₃) and negatively charged (COOH) groups are discussed to be used for the specific and strong interactions with spike protein. Bio-interfacial interactions between the spike protein and device (sensor) surface and its implications toward designing suitable engineered surfaces are summarized.

EXPERT OPINION

Implementation of the engineered surfaces in quartz crystal microbalance (QCM)-based detection techniques for the diagnosis of the novel coronavirus from oral swab samples is highlighted. The proposed strategy can be explored for the label-free and real-time detection with sensitivity up to ng level. These engineered surfaces can be reused after desorption.

In silico guided structural and functional analysis of genes with potential involvement in resistance to coffee leaf rust: A functional marker based approach

Physiology-based differentiation of S_H genes and Hemileia vastatrix races is the principal method employed for the characterization of coffee leaf rust resistance. Based on the gene-for-gene theory, nine major rust resistance genes (S_H1-9) have been proposed. However, these genes have not been characterized at the molecular level. Consequently, the lack of molecular data regarding rust resistance genes or candidates is a major bottleneck in coffee breeding. To address this issue, we screened a BAC library with resistance gene analogs (RGAs), identified RGAs, characterized and explored for any S_H related candidate genes. Herein, we report the identification and characterization of a gene (gene 11), which shares conserved sequences with other S_H genes and displays a characteristic polymorphic allele conferring different resistance phenotypes. Furthermore, comparative analysis of the two RGAs belonging to CC-NBS-LRR revealed more intense diversifying selection in tomato and grape genomes than in coffee. For the first time, the present study has unveiled novel insights into the molecular nature of the S_H genes, thereby opening new avenues for coffee rust resistance molecular breeding. The characterized candidate RGA is of particular importance for further biological function analysis in coffee.