Surprising similarities in structure comparison

Fitness landscape of substrate-adaptive mutations in evolved amino acid-polyamine-organocation transporters

The emergence of new protein functions is crucial for the evolution of organisms. This process has been extensively researched for soluble enzymes, but it is largely unexplored for membrane transporters, even though the ability to acquire new nutrients from a changing environment requires evolvability of transport functions. Here, we demonstrate the importance of environmental pressure in obtaining a new activity or altering a promiscuous activity in members of the amino acid-polyamine-organocation (APC)-type yeast amino acid transporters family. We identify APC members that have broader substrate spectra than previously described. Using in vivo experimental evolution, we evolve two of these transporter genes, AGP1 and PUT4, toward new substrate specificities. Single mutations on these transporters are found to be sufficient for expanding the substrate range of the proteins, while retaining the capacity to transport all original substrates. Nonetheless, each adaptive mutation comes with a distinct effect on the fitness for each of the original substrates, illustrating a trade-off between the ancestral and evolved functions. Collectively, our findings reveal how substrate-adaptive mutations in membrane transporters contribute to fitness and provide insights into how organisms can use transporter evolution to explore new ecological niches.

Comparative proteomics of a versatile, marine, iron-oxidizing chemolithoautotroph

This study conducted a comparative proteomic analysis to identify potential genetic markers for the biological function of chemolithoautotrophic iron oxidation in the marine bacterium Ghiorsea bivora. To date, this is the only characterized species in the class Zetaproteobacteria that is not an obligate iron-oxidizer, providing a unique opportunity to investigate differential protein expression to identify key genes involved in iron-oxidation at circumneutral pH. Over 1000 proteins were identified under both iron- and hydrogen-oxidizing conditions, with differentially expressed proteins found in both treatments. Notably, a gene cluster upregulated during iron oxidation was identified. This cluster contains genes encoding for cytochromes that share sequence similarity with the known iron-oxidase, Cyc2. Interestingly, these cytochromes, conserved in both Bacteria and Archaea, do not exhibit the typical β-barrel structure of Cyc2. This cluster potentially encodes a biological nanowire-like transmembrane complex containing multiple redox proteins spanning the inner membrane, periplasm, outer membrane, and extracellular space. The upregulation of key genes associated with this complex during iron-oxidizing conditions was confirmed by quantitative reverse transcription-PCR. These findings were further supported by electromicrobiological methods, which demonstrated negative current production by G. bivora in a three-electrode system poised at a cathodic potential. This research provides significant insights into the biological function of chemolithoautotrophic iron oxidation.

NMY-2, TOE-2 and PIG-1 regulate Caenorhabditis elegans asymmetric cell divisions

Asymmetric cell division is an important mechanism that generates cellular diversity during development. Not only do asymmetric cell divisions produce daughter cells of different fates, but many can also produce daughters of different sizes, which we refer to as Daughter Cell Size Asymmetry (DCSA). In Caenorhabditis elegans, apoptotic cells are frequently produced by asymmetric divisions that exhibit DCSA, where the smaller daughter dies. We focus here on the divisions of the Q.a and Q.p neuroblasts, which produce larger surviving cells and smaller apoptotic cells and divide with opposite polarity using both distinct and overlapping mechanisms. Several proteins regulate DCSA in these divisions. Previous studies showed that the PIG-1/MELK and TOE-2 proteins regulate DCSA in both the Q.a and Q.p divisions, and the non-muscle myosin NMY-2 regulates DCSA in the Q.a division but not the Q.p division. In this study, we examined endogenously tagged NMY-2, TOE-2, and PIG-1 reporters and characterized their distribution at the cortex during the Q.a and Q.p divisions. In both divisions, TOE-2 localized toward the side of the dividing cell that produced the smaller daughter, whereas PIG-1 localized toward the side that produced the larger daughter. As previously reported, NMY-2 localized to the side of Q.a that produced the smaller daughter and did not localize asymmetrically in Q.p. We used temperature-sensitive nmy-2 mutants to determine the role of nmy-2 in these divisions and were surprised to find that these mutants only displayed DCSA defects in the Q.p division. We generated double mutant combinations between the nmy-2 mutations and mutations in toe-2 and pig-1. Because previous studies indicate that DCSA defects result in the transformation of cells fated to die into their sister cells, the finding that the nmy-2 mutations did not significantly alter the Q.a and Q.p DCSA defects of toe-2 and pig-1 mutants but did alter the number of daughter cells produced by Q.a and Q.p suggests that nmy-2 plays a role in specifying the fates of the Q.a and Q.p that is independent of its role in DCSA.

Identification of domains in Plasmodium falciparum proteins of unknown function using DALI search on AlphaFold predictions

Plasmodium falciparum, the causative agent of malaria, poses a significant global health challenge, yet much of its biology remains elusive. A third of the genes in the P. falciparum genome lack annotations regarding their function, impeding our understanding of the parasite's biology. In this study, we employ structure predictions and the DALI search algorithm to analyse proteins encoded by uncharacterized genes in the reference strain 3D7 of P. falciparum. By comparing AlphaFold predictions to experimentally determined protein structures in the Protein Data Bank, we found similarities to known domains in 353 proteins of unknown function, shedding light on their potential functions. The lowest-scoring 5% of similarities were additionally validated using the size-independent TM-align algorithm, confirming the detected similarities in 88% of the cases. Notably, in over 70 P. falciparum proteins the presence of domains resembling heptatricopeptide repeats, which are typically involvement in RNA binding and processing, was detected. This suggests this family, which is important in transcription in mitochondria and apicoplasts, is much larger in Plasmodium parasites than previously thought. The results of this domain search provide a resource to the malaria research community that is expected to inform and enable experimental studies.

Bacillus subtilis YisK possesses oxaloacetate decarboxylase activity and exhibits Mbl-dependent localization

YisK is an uncharacterized protein in Bacillus subtilis previously shown to interact genetically with the elongasome protein Mbl. YisK overexpression leads to cell widening and lysis, phenotypes that are dependent on mbl and suppressed by mbl mutations. In the present work, we characterize YisK's localization, structure, and enzymatic activity. We show that YisK localizes as puncta that depend on Mbl. YisK belongs to the fumarylacetoacetate hydrolase (FAH) superfamily, and crystal structures revealed close structural similarity to two oxaloacetate (OAA) decarboxylases: human mitochondrial FAHD1 and Corynebacterium glutamicum Cg1458. We demonstrate that YisK can also catalyze the decarboxylation of OAA (K m = 134 µM, K cat = 31 min-1). A catalytic dead variant (YisK E148A, E150A) retains wild-type localization and still widens cells following overexpression, indicating these activities are not dependent on YisK catalysis. Conversely, a non-localizing variant (YisK E30A) retains wild-type enzymatic activity in vitro but localizes diffusely and no longer widens cells following overexpression. Together, these results suggest that YisK may be subject to spatial regulation that depends on the cell envelope synthesis machinery. IMPORTANCE The elongasome is a multiprotein complex that guides lengthwise growth in some bacteria. We previously showed that, in B. subtilis, overexpression of an uncharacterized putative enzyme (YisK) perturbed function of the actin-like elongasome protein Mbl. Here, we show that YisK exhibits Mbl-dependent localization. Through biochemical and structural characterization, we demonstrate that, like its mitochondrial homolog FAHD1, YisK can catalyze the decarboxylation of the oxaloacetate to pyruvate and CO2. YisK is the first example of an enzyme implicated in central carbon metabolism with subcellular localization that depends on Mbl.

From TgO/GABA-AT, GABA, and T-263 Mutant to Conception of Toxoplasma

Toxoplasma gondii causes morbidity, mortality, and disseminates widely via cat sexual stages. Here, we find T. gondii ornithine aminotransferase (OAT) is conserved across phyla. We solve TgO/GABA-AT structures with bound inactivators at 1.55 Å and identify an inactivator selective for TgO/GABA-AT over human OAT and GABA-AT. However, abrogating TgO/GABA-AT genetically does not diminish replication, virulence, cyst-formation, or eliminate cat's oocyst shedding. Increased sporozoite/merozoite TgO/GABA-AT expression led to our study of a mutagenized clone with oocyst formation blocked, arresting after forming male and female gametes, with "Rosetta stone"-like mutations in genes expressed in merozoites. Mutations are similar to those in organisms from plants to mammals, causing defects in conception and zygote formation, affecting merozoite capacitation, pH/ionicity/sodium-GABA concentrations, drawing attention to cyclic AMP/PKA, and genes enhancing energy or substrate formation in TgO/GABA-AT-related-pathways. These candidates potentially influence merozoite's capacity to make gametes that fuse to become zygotes, thereby contaminating environments and causing disease.

Novel pathogenic variant in MED12 causing non-syndromic dilated cardiomyopathy

BACKGROUND

Dilated cardiomyopathy (DCM) is a major cause of sudden cardiac death and heart failure. Up to 50% of all DCM cases have a genetic background, with variants in over 250 genes reported in association with DCM. Whole-exome sequencing (WES) is a powerful tool to identify variants underlying genetic cardiomyopathies. Via WES, we sought to identify DCM causes in a family with 2 affected patients.

METHODS

WES was performed on the affected members of an Iranian family to identify the genetic etiology of DCM. The candidate variant was segregated via polymerase chain reaction and Sanger sequencing. Computational modeling and protein-protein docking were performed to survey the impact of the variant on the structure and function of the protein.

RESULTS

A novel single-nucleotide substitution (G > A) in exon 9 of MED12, c.1249G > A: p.Val417Ile, NM_005120.3, was identified. The c.1249G > A variant was validated in the family. Bioinformatic analysis and computational modeling confirmed that c.1249G > A was the pathogenic variant responsible for the DCM phenotype.

CONCLUSION

We detected a novel DCM-causing variant in MED12 using WES. The variant in MED12 may decrease binding to cyclin-dependent kinase 8 (CDK8), affect its activation, and cause alterations in calcium-handling gene expression in the heart, leading to DCM.

The Kelch13 compartment contains highly divergent vesicle trafficking proteins in malaria parasites

Single amino acid changes in the parasite protein Kelch13 (K13) result in reduced susceptibility of P. falciparum parasites to artemisinin and its derivatives (ART). Recent work indicated that K13 and other proteins co-localising with K13 (K13 compartment proteins) are involved in the endocytic uptake of host cell cytosol (HCCU) and that a reduction in HCCU results in reduced susceptibility to ART. HCCU is critical for parasite survival but is poorly understood, with the K13 compartment proteins among the few proteins so far functionally linked to this process. Here we further defined the composition of the K13 compartment by analysing more hits from a previous BioID, showing that MyoF and MCA2 as well as Kelch13 interaction candidate (KIC) 11 and 12 are found at this site. Functional analyses, tests for ART susceptibility as well as comparisons of structural similarities using AlphaFold2 predictions of these and previously identified proteins showed that vesicle trafficking and endocytosis domains were frequent in proteins involved in resistance or endocytosis (or both), comprising one group of K13 compartment proteins. While this strengthened the link of the K13 compartment to endocytosis, many proteins of this group showed unusual domain combinations and large parasite-specific regions, indicating a high level of taxon-specific adaptation of this process. Another group of K13 compartment proteins did not influence endocytosis or ART susceptibility and lacked detectable vesicle trafficking domains. We here identified the first protein of this group that is important for asexual blood stage development and showed that it likely is involved in invasion. Overall, this work identified novel proteins functioning in endocytosis and at the K13 compartment. Together with comparisons of structural predictions it provides a repertoire of functional domains at the K13 compartment that indicate a high level of adaption of endocytosis in malaria parasites.

The RavA-ViaA chaperone complex modulates bacterial persistence through its association with the fumarate reductase enzyme

Regulatory ATPase variant A (RavA) is a MoxR AAA+ protein that functions together with a partner protein termed von Willebrand factor type A interacting with AAA+ ATPase (ViaA). RavA-ViaA are functionally associated with anaerobic respiration in Escherichia coli through interactions with the fumarate reductase (Frd) electron transport complex. Through this association, RavA and ViaA modulate the activity of the Frd complex and, hence, are proposed to have chaperone-like activity. However, the functional role of RavA-ViaA in the cell is not yet well established. We had demonstrated that RavA-ViaA can sensitize E. coli cells to sublethal concentrations of the aminoglycoside class of antibiotics. Since Frd has been associated with bacterial persistence against antibiotics, the relationship of RavA-ViaA and Frd was explored within this context. Experiments performed here reveal a function of RavA-ViaA in bacterial persistence upon treatment with antibiotics through the association of the chaperone complex with Frd. As part of this work, the NMR structure of the N-terminal domain of ViaA was solved. The structure reveals a novel alpha helical fold, which we name the VAN fold, that has not been observed before. We show that this domain is required for the function of the chaperone complex. We propose that modulating the levels of RavA-ViaA could enhance the susceptibility of Gram-negative bacteria to antibiotics.

The History and Science of the Major Birch Pollen Allergen Bet v 1

The term allergy was coined in 1906 by the Austrian scientist and pediatrician Clemens Freiherr von Pirquet. In 1976, Dietrich Kraft became the head of the Allergy and Immunology Research Group at the Department of General and Experimental Pathology of the University of Vienna. In 1983, Kraft proposed to replace natural extracts used in allergy diagnostic tests and vaccines with recombinant allergen molecules and persuaded Michael Breitenbach to contribute his expertise in molecular cloning as one of the mentors of this project. Thus, the foundation for the Vienna School of Molecular Allergology was laid. With the recruitment of Heimo Breiteneder as a young molecular biology researcher, the work began in earnest, resulting in the publication of the cloning of the first plant allergen Bet v 1 in 1989. Bet v 1 has become the subject of a very large number of basic scientific as well as clinical studies. Bet v 1 is also the founding member of the large Bet v 1-like superfamily of proteins with members-based on the ancient conserved Bet v 1 fold-being present in all three domains of life, i.e., archaea, bacteria and eukaryotes. This suggests that the Bet v 1 fold most likely already existed in the last universal common ancestor. The biological function of this protein was probably related to lipid binding. However, during evolution, a functional diversity within the Bet v 1-like superfamily was established. The superfamily comprises 25 families, one of which is the Bet v 1 family, which in turn is composed of 11 subfamilies. One of these, the PR-10-like subfamily of proteins, contains almost all of the Bet v 1 homologous allergens from pollen and plant foods. Structural and functional comparisons of Bet v 1 and its non-allergenic homologs of the superfamily will pave the way for a deeper understanding of the allergic sensitization process.

Feature architecture aware phylogenetic profiling indicates a functional diversification of type IVa pili in the nosocomial pathogen Acinetobacter baumannii

The Gram-negative bacterial pathogen Acinetobacter baumannii is a major cause of hospital-acquired opportunistic infections. The increasing spread of pan-drug resistant strains makes A. baumannii top-ranking among the ESKAPE pathogens for which novel routes of treatment are urgently needed. Comparative genomics approaches have successfully identified genetic changes coinciding with the emergence of pathogenicity in Acinetobacter. Genes that are prevalent both in pathogenic and a-pathogenic Acinetobacter species were not considered ignoring that virulence factors may emerge by the modification of evolutionarily old and widespread proteins. Here, we increased the resolution of comparative genomics analyses to also include lineage-specific changes in protein feature architectures. Using type IVa pili (T4aP) as an example, we show that three pilus components, among them the pilus tip adhesin ComC, vary in their Pfam domain annotation within the genus Acinetobacter. In most pathogenic Acinetobacter isolates, ComC displays a von Willebrand Factor type A domain harboring a finger-like protrusion, and we provide experimental evidence that this finger conveys virulence-related functions in A. baumannii. All three genes are part of an evolutionary cassette, which has been replaced at least twice during A. baumannii diversification. The resulting strain-specific differences in T4aP layout suggests differences in the way how individual strains interact with their host. Our study underpins the hypothesis that A. baumannii uses T4aP for host infection as it was shown previously for other pathogens. It also indicates that many more functional complexes may exist whose precise functions have been adjusted by modifying individual components on the domain level.

In Vitro Susceptibility and Resistance of Mycoplasma genitalium to Nitroimidazoles

Mycoplasma genitalium is a sexually transmitted reproductive tract pathogen of men and women. M. genitalium infections are increasingly difficult to treat due to poor efficacy of doxycycline and acquired resistance to azithromycin and moxifloxacin. A recent clinical trial suggested that metronidazole may improve cure rates for women with pelvic inflammatory disease and reduced the detection of M. genitalium when included with standard doxycycline plus ceftriaxone treatment. As data regarding susceptibility of mycoplasmas to nitroimidazoles are lacking in the scientific literature, we determined the in vitro susceptibility of 10 M. genitalium strains to metronidazole, secnidazole, and tinidazole. MICs ranged from 1.6 to 12.5 μg/mL for metronidazole, 3.1 to 12.5 μg/mL for secnidazole, and 0.8 to 6.3 μg/mL for tinidazole. None of these agents was synergistic with doxycycline in checkerboard broth microdilution assays. Tinidazole was superior to metronidazole and secnidazole in terms of MIC and time-kill kinetics and was bactericidal (>99.9% killing) at concentrations below reported serum concentrations. Mutations associated with nitroimidazole resistance were identified by whole-genome sequencing of spontaneous resistant mutants, suggesting a mechanism for reductive activation of the nitroimidazole prodrug by a predicted NAD(P)H-dependent flavin mononucleotide (FMN) oxidoreductase. The presence of oxygen did not affect MICs of wild-type M. genitalium, but a nitroimidazole-resistant mutant was defective for growth under anaerobic conditions, suggesting that resistant mutants may have a fitness disadvantage in anaerobic genital sites. Clinical studies are needed to determine if nitroimidazoles, especially tinidazole, are effective for eradicating M. genitalium infections in men and women.

Glycation Interferes with the Activity of the Bi-Functional UDP-N-Acetylglucosamine 2-Epimerase/N-Acetyl-mannosamine Kinase (GNE)

Mutations in the gene coding for the bi-functional UDP-N-acetylglucosamine 2-epimerase/N-acetylmannosamine kinase (GNE), the key enzyme of the sialic acid biosynthesis, are responsible for autosomal-recessive GNE myopathy (GNEM). GNEM is an adult-onset disease with a yet unknown exact pathophysiology. Since the protein appears to work adequately for a certain period of time even though the mutation is already present, other effects appear to influence the onset and progression of the disease. In this study, we want to investigate whether the late onset of GNEM is based on an age-related effect, e.g., the accumulation of post-translational modifications (PTMs). Furthermore, we also want to investigate what effect on the enzyme activity such an accumulation would have. We will particularly focus on glycation, which is a PTM through non-enzymatic reactions between the carbonyl groups (e.g., of methylglyoxal (MGO) or glyoxal (GO)) with amino groups of proteins or other biomolecules. It is already known that the levels of both MGO and GO increase with age. For our investigations, we express each domain of the GNE separately, treat them with one of the glycation agents, and determine their activity. We demonstrate that the enzymatic activity of the N-acetylmannosamine kinase (GNE-kinase domain) decreases dramatically after glycation with MGO or GO-with a remaining activity of 13% ± 5% (5 mM MGO) and 22% ± 4% (5 mM GO). Whereas the activity of the UDP-N-acetylglucosamine 2-epimerase (GNE-epimerase domain) is only slightly reduced after glycation-with a remaining activity of 60% ± 8% (5 mM MGO) and 63% ± 5% (5 mM GO).

Structural and functional characterization of RNA dependent RNA polymerase of Macrobrachium rosenbergii nodavirus (MnRdRp)

Macrobrachium rosenbergii is a highly valued farmed freshwater species and its production has been affected globally by white tail disease caused by M. rosenbergii nodavirus (MrNV). MrNV is a single stranded positive sense RNA virus encoding RNA-dependent RNA polymerase (RdRp) for genome replication. Due to its essentiality for pathogenesis, it is an important drug target. The domain prediction of the complete sequence revealed the presence of two enzymatic regions namely methyl transferase and RdRp separated by transmembrane region. The predicted three-dimensional (3D) structure of MnRdRp using AlphaFold 2 shows that the structure is composed of three major sub-domains common for other polymerases namely fingers, palm and thumb. Structural similarity search revealed its similarity with other flaviviridea members especially with BVDV RdRp (BvdvRdRp). The structure of fingers and palm sub-domains is more conserved than the thumb sub-domain. A small α-helix named 'priming helix' having conserve Tyr was identified at position 829-833 with a potential role in de novo initiation. Analysis of electrostatic potential revealed that nucleotide and template channels are electropositive. Metal binding residues were identified as Asp599, Asp704 and Asp705. The α and β phosphates of incoming nucleotide interact with two Mn2+, Arg455 and Arg537. For recognition of 2'-OH of incoming rNTP, Asp604, Ser661 and Asn670 were identified which can form H-bond network with 2'-OH group. Docking study revealed that Dasabuvir can potentially inhibit MnRdRp. The study concluded that the overall structure and function of MnRdRp are similar to Flaviviridae polymerases and their inhibitors can work against this enzyme.Communicated by Ramaswamy H. Sarma.

Eight Unexpected Selenoprotein Families in Organometallic Biochemistry in Clostridium difficile, in ABC Transport, and in Methylmercury Biosynthesis

The bioinformatics of a nine-gene locus, designated selenocysteine-assisted organometallic (SAO), was investigated after identifying six new selenoprotein families and constructing hidden Markov models (HMMs) that find and annotate members of those families. Four are selenoproteins in most SAO loci, including Clostridium difficile. They include two ABC transporter subunits, namely, permease SaoP, with selenocysteine (U) at the channel-gating position, and substrate-binding subunit SaoB. Cytosolic selenoproteins include SaoL, homologous to MerB organomercurial lyases from mercury resistance loci, and SaoT, related to thioredoxins. SaoL, SaoB, and surface protein SaoC (an occasional selenoprotein) share an unusual CU dipeptide motif, which is something rare in selenoproteins but found in selenoprotein variants of mercury resistance transporter subunit MerT. A nonselenoprotein, SaoE, shares homology with Cu/Zn efflux and arsenical efflux pumps. The organization of the SAO system suggests substrate interaction with surface-exposed selenoproteins, followed by import, metabolism that may cleave a carbon-to-heavy metal bond, and finally metal efflux. A novel type of mercury resistance is possible, but SAO instead may support fermentative metabolism, with selenocysteine-mediated formation of organometallic intermediates, followed by import, degradation, and metal efflux. Phylogenetic profiling shows SOA loci consistently co-occur with Stickland fermentation markers but even more consistently with 8Fe-9S cofactor-type double-cubane proteins. Hypothesizing that the SAO system forms organometallic intermediates, we investigated the known methylmercury formation protein families HgcA and HgcB. Both families contained overlooked selenoproteins. Most HgcAs have a CU motif N terminal to their previously accepted start sites. Seeking additional rare and overlooked selenoproteins may help reveal more cryptic aspects of microbial biochemistry. IMPORTANCE This work adds 8 novel prokaryotic selenoproteins to the 80 or so families previously known. It describes the SAO (selenocysteine-assisted organometallic) locus, with the most selenoproteins of any known system. The rare CU motif recurs throughout, suggesting the formation and degradation of organometallic compounds. That suggestion triggered a reexamination of HgcA and HcgB, which are methylmercury formation proteins that can adversely impact food safety. Both are selenoproteins, once corrected, with HgcA again showing a CU motif. The SAO system is plausibly a mercury resistance locus for selenium-dependent anaerobes. But instead, it may exploit heavy metals as cofactors in organometallic intermediate-forming pathways that circumvent high activation energies and facilitate the breakdown of otherwise poorly accessible nutrients. SAO could provide an edge that helps Clostridium difficile, an important pathogen, establish disease.

Database resources of the National Center for Biotechnology Information in 2023

The National Center for Biotechnology Information (NCBI) provides online information resources for biology, including the GenBank® nucleic acid sequence database and the PubMed® database of citations and abstracts published in life science journals. NCBI provides search and retrieval operations for most of these data from 35 distinct databases. The E-utilities serve as the programming interface for most of these databases. New resources include the Comparative Genome Resource (CGR) and the BLAST ClusteredNR database. Resources receiving significant updates in the past year include PubMed, PMC, Bookshelf, IgBLAST, GDV, RefSeq, NCBI Virus, GenBank type assemblies, iCn3D, ClinVar, GTR, dbGaP, ALFA, ClinicalTrials.gov, Pathogen Detection, antimicrobial resistance resources, and PubChem. These resources can be accessed through the NCBI home page at https://www.ncbi.nlm.nih.gov.

The genome wide analysis of Tryptophan Aminotransferase Related gene family, and their relationship with related agronomic traits in Brassica napus

Tryptophan Aminotransferase of Arabidopsis1/Tryptophan Aminotransferase-Related (TAA1/TAR) proteins are the enzymes that involved in auxin biosynthesis pathway. The TAA1/TAR gene family has been systematically characterized in several plants but has not been well reported in Brassica napus. In the present study, a total of 102 BnTAR genes with different number of introns were identified. It was revealed that these genes are distributed unevenly and occurred as clusters on different chromosomes except for A4, A5, A10 and C4 in B. napus. Most of the these BnTAR genes are conserved despite of existing of gene loss and gene gain. In addition, the segmental replication and whole-genome replication events were both play an important role in the BnTAR gene family formation. Expression profiles analysis indicated that the expression of BnTAR gene showed two patterns, part of them were mainly expressed in roots, stems and leaves of vegetative organs, and the others were mainly expressed in flowers and seeds of reproductive organs. Further analysis showed that many of BnTAR genes were located in QTL intervals of oil content or seed weight, for example BnAMI10 was located in cqOC-C5-4 and cqSW-A2-2, it indicated that some of the BnTAR genes might have relationship with these two characteristics. This study provides a multidimensional analysis of the TAA1/TAR gene family and a new insight into its biological function in B. napus.

Newly Discovered Peptides from the Coral Heliofungia actiniformis Show Structural and Functional Diversity

Scleractinian corals are crucially important to the health of some of the world's most biodiverse, productive, and economically important marine habitats. Despite this importance, analysis of coral peptidomes is still in its infancy. Here we show that the tentacle extract from the stony coral Heliofungia actiniformis is rich in peptides with diverse and novel structures. We have characterized the sequences and three-dimensional structures of four new peptides, three of which have no known homologues. We show that a 2 kDa peptide, Hact-2, promotes significant cell proliferation on human cells and speculate this peptide may be involved in the remarkable regenerative capacity of corals. We found a 3 kDa peptide, Hact-3, encoded within a fascin-like domain, and homologues of Hact-3 are present in the genomes of other coral species. Two additional peptides, Hact-4 and Hact-SCRiP1, with limited sequence similarity, both contain a beta-defensin-like fold and highlight a structural link with the small cysteine-rich proteins (SCRiP) family of proteins found predominantly in corals. Our results provide a first glimpse into the remarkable and unexplored structural diversity of coral peptides, providing insight into their diversity and putative functions and, given the ancient lineage of corals, potential insight into the evolution of structural motifs.

Soils and sediments host Thermoplasmata archaea encoding novel copper membrane monooxygenases (CuMMOs)

Copper membrane monooxygenases (CuMMOs) play critical roles in the global carbon and nitrogen cycles. Organisms harboring these enzymes perform the first, and rate limiting, step in aerobic oxidation of ammonia, methane, or other simple hydrocarbons. Within archaea, only organisms in the order Nitrososphaerales (Thaumarchaeota) encode CuMMOs, which function exclusively as ammonia monooxygenases. From grassland and hillslope soils and aquifer sediments, we identified 20 genomes from distinct archaeal species encoding divergent CuMMO sequences. These archaea are phylogenetically clustered in a previously unnamed Thermoplasmatota order, herein named the Ca. Angelarchaeales. The CuMMO proteins in Ca. Angelarchaeales are more similar in structure to those in Nitrososphaerales than those of bacteria, and contain all functional residues required for general monooxygenase activity. Ca. Angelarchaeales genomes are significantly enriched in blue copper proteins (BCPs) relative to sibling lineages, including plastocyanin-like electron carriers and divergent nitrite reductase-like (nirK) 2-domain cupredoxin proteins co-located with electron transport machinery. Ca. Angelarchaeales also encode significant capacity for peptide/amino acid uptake and degradation and share numerous electron transport mechanisms with the Nitrososphaerales. Ca. Angelarchaeales are detected at high relative abundance in some of the environments where their genomes originated from. While the exact substrate specificities of the novel CuMMOs identified here have yet to be determined, activity on ammonia is possible given their metabolic and ecological context. The identification of an archaeal CuMMO outside of the Nitrososphaerales significantly expands the known diversity of CuMMO enzymes in archaea and suggests previously unaccounted organisms contribute to critical global nitrogen and/or carbon cycling functions.

iCn3D: From Web-Based 3D Viewer to Structural Analysis Tool in Batch Mode

iCn3D was initially developed as a web-based 3D molecular viewer. It then evolved from visualization into a full-featured interactive structural analysis software. It became a collaborative research instrument through the sharing of permanent, shortened URLs that encapsulate not only annotated visual molecular scenes, but also all underlying data and analysis scripts in a FAIR manner. More recently, with the growth of structural databases, the need to analyze large structural datasets systematically led us to use Python scripts and convert the code to be used in Node. js scripts. We showed a few examples of Python scripts at https://github.com/ncbi/icn3d/tree/master/icn3dpython to export secondary structures or PNG images from iCn3D. Users just need to replace the URL in the Python scripts to export other annotations from iCn3D. Furthermore, any interactive iCn3D feature can be converted into a Node. js script to be run in batch mode, enabling an interactive analysis performed on one or a handful of protein complexes to be scaled up to analysis features of large ensembles of structures. Currently available Node. js analysis scripts examples are available at https://github.com/ncbi/icn3d/tree/master/icn3dnode. This development will enable ensemble analyses on growing structural databases such as AlphaFold or RoseTTAFold on one hand and Electron Microscopy on the other. In this paper, we also review new features such as DelPhi electrostatic potential, 3D view of mutations, alignment of multiple chains, assembly of multiple structures by realignment, dynamic symmetry calculation, 2D cartoons at different levels, interactive contact maps, and use of iCn3D in Jupyter Notebook as described at https://pypi.org/project/icn3dpy.

Computational Mining and Characterization of Hypothetical Proteins of Mycobacterium bovis Toward the Identification of Probable Vaccine Candidates

A hypothetical protein (HP) is one that is known to exist only on the basis of a corresponding gene but without any function assigned to it. Many HPs have emerged as attractive vaccine candidates against prokaryotic and eukaryotic pathogens as well as against cancers. Mycobacterium bovis is a serious veterinary pathogen of tremendous zoonotic importance. This protocol describes a computational workflow for the identification of the HPs of M. bovis with vaccine potential and their subsequent structural and functional characterization.

Nonlethal presentations of CYP26B1-related skeletal anomalies and multiple synostoses syndrome

Retinoic acid exposures as well as defects in the retinoic acid-degrading enzyme CYP26B1 have teratogenic effects on both limb and craniofacial skeleton. An initial report of four individuals described a syndrome of fetal and infantile lethality with craniosynostosis and skeletal anomalies caused by homozygous pathogenic missense variants in CYP26B1. In contrast, a 22-year-old female was reported with a homozygous missense pathogenic variant in CYP26B1 with complex multisuture craniosynostosis and intellectual disability, suggesting that in some cases, biallelic pathogenic variants of CYP26B1 may be compatible with life. Here we describe four additional living individuals from two families with compound heterozygous pathogenic missense variants in CYP26B1. Structural assessment of these additional missense variants places them further from the catalytic site and supports a model consistent with milder nonlethal disease. In addition to previously reported findings of multisuture craniosynostosis, conductive hearing loss, joint contractures, long slender fingers, camptodactly, broad fingertips, and developmental delay/intellectual disability, skeletal imaging in our cases also revealed gracile long bones, gracile ribs, radioulnar synostosis, and carpal and/or tarsal fusions. These individuals broaden the phenotypic range of biallelic pathogenic variants in CYPB26B1 and most significantly clarify that mortality can range from perinatal lethality to survival into adulthood.

SWI/SNF subunit BAF155 N-terminus structure informs the impact of cancer-associated mutations and reveals a potential drug binding site

SWI/SNF (BAF) chromatin remodelling complexes are key regulators of gene expression programs, and attractive drug targets for cancer therapies. Here we show that the N-terminus of the BAF155/SMARCC1 subunit contains a putative DNA-binding MarR-like domain, a chromodomain and a BRCT domain that are interconnected to each other to form a distinct module. In this structure the chromodomain makes interdomain interactions and has lost its canonical function to bind to methylated lysines. The structure provides new insights into the missense mutations that target this module in cancer. This study also reveals two adjacent, highly-conserved pockets in a cleft between the domains that form a potential binding site, which can be targeted with small molecules, offering a new strategy to target SWI/SNF complexes.

Allen et al. determine crystal structure of the N-terminus of SWI/SNF chromatin remodelling factor subunit BAF155. They identify an interconnected structure of MarR-like, BRCT and chromodomains, visualise the potential effect of cancer-associated missense mutations and suggest a binding site and target for small molecule inhibitors.

NMR structure and dynamics of Q4DY78, a conserved kinetoplasid-specific protein from Trypanosoma cruzi

The 106-residue protein Q4DY78 (UniProt accession number) from Trypanosoma cruzi is highly conserved in the related kinetoplastid pathogens Trypanosoma brucei and Leishmania major. Given the essentiality of its orthologue in T. brucei, the high sequence conservation with other trypanosomatid proteins, and the low sequence similarity with mammalian proteins, Q4DY78 is an attractive protein for structural characterization. Here, we solved the structure of Q4DY78 by solution NMR and evaluated its backbone dynamics. Q4DY78 is composed of five α -helices and a small, two-stranded antiparallel β-sheet. The backbone RMSD is 0.22 ± 0.05 Å for the representative ensemble of the 20 lowest-energy structures. Q4DY78 is overall rigid, except for N-terminal residues (V⁸ to I¹⁰), residues at loop 4 (K⁵⁷ to G⁶⁵) and residues at the C-terminus (F⁸⁹ to F¹¹²). Q4DY78 has a short motif FPCAP that could potentially mediate interactions with the host cytoskeleton via interaction with EVH1 (Drosophila Enabled (Ena)/Vasodilator-stimulated phosphoprotein (VASP) homology 1) domains. Albeit Q4DY78 lacks calcium-binding motifs, its fold resembles that of eukaryotic calcium-binding proteins such as calcitracin, calmodulin, and polcacin Bet V4. We characterized this novel protein with a calcium binding fold without the capacity to bind calcium.

Demystifying Chronic Kidney Disease of Unknown Etiology (CKDu): Computational Interaction Analysis of Pesticides and Metabolites with Vital Renal Enzymes

Chronic kidney disease of unknown etiology (CKDu) has been recognized as a global non-communicable health issue. There are many proposed risk factors for CKDu and the exact reason is yet to be discovered. Understanding the inhibition or manipulation of vital renal enzymes by pesticides can play a key role in understanding the link between CKDu and pesticides. Even though it is very important to take metabolites into account when investigating the relationship between CKDu and pesticides, there is a lack of insight regarding the effects of pesticide metabolites towards CKDu. In this study, a computational approach was used to study the effects of pesticide metabolites on CKDu. Further, interactions of selected pesticides and their metabolites with renal enzymes were studied using molecular docking and molecular dynamics simulation studies. It was evident that some pesticides and metabolites have affinity to bind at the active site or at regulatory sites of considered renal enzymes. Another important discovery was the potential of some metabolites to have higher binding interactions with considered renal enzymes compared to the parent pesticides. These findings raise the question of whether pesticide metabolites may be a main risk factor towards CKDu.

Development of a TSR-Based Method for Protein 3-D Structural Comparison With Its Applications to Protein Classification and Motif Discovery

Development of protein 3-D structural comparison methods is important in understanding protein functions. At the same time, developing such a method is very challenging. In the last 40 years, ever since the development of the first automated structural method, ~200 papers were published using different representations of structures. The existing methods can be divided into five categories: sequence-, distance-, secondary structure-, geometry-based, and network-based structural comparisons. Each has its uniqueness, but also limitations. We have developed a novel method where the 3-D structure of a protein is modeled using the concept of Triangular Spatial Relationship (TSR), where triangles are constructed with the C_α atoms of a protein as vertices. Every triangle is represented using an integer, which we denote as “key,” A key is computed using the length, angle, and vertex labels based on a rule-based formula, which ensures assignment of the same key to identical TSRs across proteins. A structure is thereby represented by a vector of integers. Our method is able to accurately quantify similarity of structure or substructure by matching numbers of identical keys between two proteins. The uniqueness of our method includes: (i) a unique way to represent structures to avoid performing structural superimposition; (ii) use of triangles to represent substructures as it is the simplest primitive to capture shape; (iii) complex structure comparison is achieved by matching integers corresponding to multiple TSRs. Every substructure of one protein is compared to every other substructure in a different protein. The method is used in the studies of proteases and kinases because they play essential roles in cell signaling, and a majority of these constitute drug targets. The new motifs or substructures we identified specifically for proteases and kinases provide a deeper insight into their structural relations. Furthermore, the method provides a unique way to study protein conformational changes. In addition, the results from CATH and SCOP data sets clearly demonstrate that our method can distinguish alpha helices from beta pleated sheets and vice versa. Our method has the potential to be developed into a powerful tool for efficient structure-BLAST search and comparison, just as BLAST is for sequence search and alignment.

Protein structure search to support the development of protein structure prediction methods

Protein structure prediction is a long-standing unsolved problem in molecular biology that has seen renewed interest with the recent success of deep learning with AlphaFold at CASP13. While developing and evaluating protein structure prediction methods, researchers may want to identify the most similar known structures to their predicted structures. These predicted structures often have low sequence and structure similarity to known structures. We show how RUPEE, a purely geometric protein structure search, is able to identify the structures most similar to structure predictions, regardless of how they vary from known structures, something existing protein structure searches struggle with. RUPEE accomplishes this through the use of a novel linear encoding of protein structures as a sequence of residue descriptors. Using a fast Needleman-Wunsch algorithm, RUPEE is able to perform alignments on the sequences of residue descriptors for every available structure. This is followed by a series of increasingly accurate structure alignments from TM-align alignments initialized with the Needleman-Wunsch residue descriptor alignments to standard TM-align alignments of the final results. By using alignment normalization effectively at each stage, RUPEE also can execute containment searches in addition to full-length searches to identify structural motifs within proteins. We compare the results of RUPEE to the protein structure searches mTM-align, SSM, CATHEDRAL, and VAST using a benchmark derived from the protein structure predictions submitted to CASP13. RUPEE identifies better alignments on average with respect to TM-score as well as scores specific to SSM and CATHEDRAL, Q-score and SSAP-score, respectively.

The Structures of SctK and SctD from Pseudomonas aeruginosa Reveal the Interface of the Type III Secretion System Basal Body and Sorting Platform

Many Gram-negative bacterial pathogens use type III secretion systems (T3SS) to inject proteins into eukaryotic cells to subvert normal cellular functions. The T3SS apparatus (injectisome) shares a common architecture in all systems studied thus far, comprising three major components - the cytoplasmic sorting platform, envelope-spanning basal body and external needle with tip complex. The sorting platform consists of an ATPase (SctN) connected to "pods" (SctQ) having six-fold symmetry via radial spokes (SctL). These pods interface with the 24-fold symmetric SctD inner membrane ring (IR) via an adaptor protein (SctK). Here we report the first high-resolution structure of a SctK protein family member, PscK from Pseudomonas aeruginosa, as well as the structure of its interacting partner, the cytoplasmic domain of PscD (SctD). The cytoplasmic domain of PscD forms a forkhead-associated (FHA) fold, like that of its homologues from other T3SS. PscK, on the other hand, forms a helix-rich structure that does not resemble any known protein fold. Based on these structural findings, we present the first model for an interaction between proteins from the sorting platform and the IR. We also test the importance of the PscD residues predicted to mediate this electrostatic interaction using a two-hybrid analysis. The functional need for these residues in vivo was then confirmed by monitoring secretion of the effector ExoU. These structures will contribute to the development of atomic-resolution models of the entire sorting platform and to our understanding of the mechanistic interface between the sorting platform and the basal body of the injectisome.

pmTM-align: scalable pairwise and multiple structure alignment with Apache Spark and OpenMP

BACKGROUND

Structure comparison can provide useful information to identify functional and evolutionary relationship between proteins. With the dramatic increase of protein structure data in the Protein Data Bank, computation time quickly becomes the bottleneck for large scale structure comparisons. To more efficiently deal with informative multiple structure alignment tasks, we propose pmTM-align, a parallel protein structure alignment approach based on mTM-align/TM-align. pmTM-align contains two stages to handle pairwise structure alignments with Spark and the phylogenetic tree-based multiple structure alignment task on a single computer with OpenMP.

RESULTS

Experiments with the SABmark dataset showed that parallelization along with data structure optimization provided considerable speedup for mTM-align. The Spark-based structure alignments achieved near ideal scalability with large datasets, and the OpenMP-based construction of the phylogenetic tree accelerated the incremental alignment of multiple structures and metrics computation by a factor of about 2-5.

CONCLUSIONS

pmTM-align enables scalable pairwise and multiple structure alignment computing and offers more timely responses for medium to large-sized input data than existing alignment tools such as mTM-align.

Exploring Proteomic Drug Targets, Therapeutic Strategies and Protein - Protein Interactions in Cancer: Mechanistic View

Protein-Protein Interactions (PPIs) drive major signalling cascades and play critical role in cell proliferation, apoptosis, angiogenesis and trafficking. Deregulated PPIs are implicated in multiple malignancies and represent the critical targets for treating cancer. Herein, we discuss the key protein-protein interacting domains implicated in cancer notably PDZ, SH2, SH3, LIM, PTB, SAM and PH. These domains are present in numerous enzymes/kinases, growth factors, transcription factors, adaptor proteins, receptors and scaffolding proteins and thus represent essential sites for targeting cancer. This review explores the candidature of various proteins involved in cellular trafficking (small GTPases, molecular motors, matrix-degrading enzymes, integrin), transcription (p53, cMyc), signalling (membrane receptor proteins), angiogenesis (VEGFs) and apoptosis (BCL-2family), which could possibly serve as targets for developing effective anti-cancer regimen. Interactions between Ras/Raf; X-linked inhibitor of apoptosis protein (XIAP)/second mitochondria-derived activator of caspases (Smac/DIABLO); Frizzled (FRZ)/Dishevelled (DVL) protein; beta-catenin/T Cell Factor (TCF) have also been studied as prospective anticancer targets. Efficacy of diverse molecules/ drugs targeting such PPIs although evaluated in various animal models/cell lines, there is an essential need for human-based clinical trials. Therapeutic strategies like the use of biologicals, high throughput screening (HTS) and fragment-based technology could play an imperative role in designing cancer therapeutics. Moreover, bioinformatic/computational strategies based on genome sequence, protein sequence/structure and domain data could serve as competent tools for predicting PPIs. Exploring hot spots in proteomic networks represents another approach for developing targetspecific therapeutics. Overall, this review lays emphasis on a productive amalgamation of proteomics, genomics, biochemistry, and molecular dynamics for successful treatment of cancer.

Chemometric Models of Differential Amino Acids at the Navα and Navβ Interface of Mammalian Sodium Channel Isoforms

(1) Background: voltage-gated sodium channels (Na_vs) are integral membrane proteins that allow the sodium ion flux into the excitable cells and initiate the action potential. They comprise an α (Na_vα) subunit that forms the channel pore and are coupled to one or more auxiliary β (Na_vβ) subunits that modulate the gating to a variable extent. (2) Methods: after performing homology in silico modeling for all nine isoforms (Na_v1.1α to Na_v1.9α), the Na_vα and Na_vβ protein-protein interaction (PPI) was analyzed chemometrically based on the primary and secondary structures as well as topological or spatial mapping. (3) Results: our findings reveal a unique isoform-specific correspondence between certain segments of the extracellular loops of the Na_vα subunits. Precisely, loop S5 in domain I forms part of the PPI and assists Na_vβ1 or Na_vβ3 on all nine mammalian isoforms. The implied molecular movements resemble macroscopic springs, all of which explains published voltage sensor effects on sodium channel fast inactivation in gating. (4) Conclusions: currently, the specific functions exerted by the Na_vβ1 or Na_vβ3 subunits on the modulation of Na_vα gating remain unknown. Our work determined functional interaction in the extracellular domains on theoretical grounds and we propose a schematic model of the gating mechanism of fast channel sodium current inactivation by educated guessing.

Prediction of disulfide bond engineering sites using a machine learning method

Disulfide bonds are covalently bonded sulfur atoms from cysteine pairs in protein structures. Due to the importance of disulfide bonds in protein folding and structural stability, artificial disulfide bonds are often engineered by cysteine mutation to enhance protein structural stability. To facilitate the experimental design, we implemented a method based on neural networks to predict amino acid pairs for cysteine mutations to form engineered disulfide bonds. The designed neural network was trained with high-resolution structures curated from the Protein Data Bank. The testing results reveal that the proposed method recognizes 99% of natural disulfide bonds. In the test with engineered disulfide bonds, the algorithm achieves similar accuracy levels with other state-of-the-art algorithms in published dataset and better performance for two comprehensively studied proteins with 70% accuracy, demonstrating potential applications in protein engineering. The neural network framework allows exploiting the full features in distance space, and therefore improves accuracy of the disulfide bond engineering site prediction. The source code and a web server are available at http://liulab.csrc.ac.cn/ssbondpre.

The Structural Basis of T4 Phage Lysis Control: DNA as the Signal for Lysis Inhibition

Optimal phage propagation depends on the regulation of the lysis of the infected host cell. In T4 phage infection, lysis occurs when the holin protein (T) forms lesions in the host membrane. However, the lethal function of T can be blocked by an antiholin (RI) during lysis inhibition (LIN). LIN sets if the infected cell undergoes superinfection, then the lysis is delayed until host/phage ratio becomes more favorable for the release of progeny. It has been thought that a signal derived from the superinfection is required to activate RI. Here we report structures that suggest a radically different model in which RI binds to T irrespective of superinfection, causing it to accumulate in a membrane as heterotetrameric 2RI-2T complex. Moreover, we show the complex binds non-specifically to DNA, suggesting that the gDNA from the superinfecting phage serves as the LIN signal and that stabilization of the complex by DNA binding is what defines LIN. Finally, we show that soluble domain of free RI crystallizes in a domain-swapped homotetramer, which likely works as a sink for RI molecules released from the RI-T complex to ensure efficient lysis. These results constitute the first structural basis and a new model not only for the historic LIN phenomenon but also for the temporal regulation of phage lysis in general.

2StrucCompare: a webserver for visualizing small but noteworthy differences between protein tertiary structures through interrogation of the secondary structure content

2StrucCompare is a webserver whose primary aim is to visualize subtle but functionally important differences between two related protein structures, either of the same protein or related homologues, with similar or functionally different tertiary structures. At the heart of the package is identifying and visualizing differences between conformations at the secondary structure and at the residue level, such as contact differences or side chain conformational differences found between two protein chains. The protein secondary structures are determined according to four established methods (DSSP, STRIDE, P-SEA and STICKS), and as each employs different assignment strategies, small conformational differences between the two structures can give rise to paired residues being denoted as having different secondary structure features with the different methods. 2StrucCompare captures both the large and more subtle differences found between structures, enabling visualization of these differences that could be key to an understanding of a proteins’ function.

2StrucCompare is freely accessible at http://2struccompare.cryst.bbk.ac.uk/index.php

NKX2-6 related congenital heart disease: Biallelic homeodomain-disrupting variants and truncus arteriosus

Congenital heart defects (CHD) are the most common birth defect and are both clinically and genetically heterogeneous. Truncus arteriosus (TA), characterized by a single arterial vessel arising from both ventricles giving rise to the coronary, pulmonary and systemic arteries, is rare and only responsible for 1% of all CHD. Two consanguineous families with TA were previously identified to have homozygous nonsense variants within the gene NKX2-6. NKX2-6 is a known downstream target of TBX1, an important transcriptional regulator implicated in the cardiac phenotype of 22q11.2 microdeletion syndrome. Herein, we report two siblings with TA presumably caused by compound heterozygous NKX2-6 variants without a history of consanguinity. Two in-house cohorts with conotruncal defects (CTD) were sequenced for variants in NKX2-6 and no additional cases of biallelic NKX2-6 variants were identified. The similar phenotype of these cases, and the clustering of variants that likely result in a truncated protein that disrupts the homeobox domain, suggest that biallelic loss of function for NKX2-6 is a rare genetic etiology for TA in particular, and possibly other types of CHD.

Neighborhood Preference of Amino Acids in Protein Structures and its Applications in Protein Structure Assessment

Amino acids form protein 3D structures in unique manners such that the folded structure is stable and functional under physiological conditions. Non-specific and non-covalent interactions between amino acids exhibit neighborhood preferences. Based on structural information from the protein data bank, a statistical energy function was derived to quantify amino acid neighborhood preferences. The neighborhood of one amino acid is defined by its contacting residues, and the energy function is determined by the neighboring residue types and relative positions. The neighborhood preference of amino acids was exploited to facilitate structural quality assessment, which was implemented in the neighborhood preference program NEPRE. The source codes are available via https://github.com/LiuLab-CSRC/NePre.

Biological applications of knowledge graph embedding models

Complex biological systems are traditionally modelled as graphs of interconnected biological entities. These graphs, i.e. biological knowledge graphs, are then processed using graph exploratory approaches to perform different types of analytical and predictive tasks. Despite the high predictive accuracy of these approaches, they have limited scalability due to their dependency on time-consuming path exploratory procedures. In recent years, owing to the rapid advances of computational technologies, new approaches for modelling graphs and mining them with high accuracy and scalability have emerged. These approaches, i.e. knowledge graph embedding (KGE) models, operate by learning low-rank vector representations of graph nodes and edges that preserve the graph's inherent structure. These approaches were used to analyse knowledge graphs from different domains where they showed superior performance and accuracy compared to previous graph exploratory approaches. In this work, we study this class of models in the context of biological knowledge graphs and their different applications. We then show how KGE models can be a natural fit for representing complex biological knowledge modelled as graphs. We also discuss their predictive and analytical capabilities in different biology applications. In this regard, we present two example case studies that demonstrate the capabilities of KGE models: prediction of drug-target interactions and polypharmacy side effects. Finally, we analyse different practical considerations for KGEs, and we discuss possible opportunities and challenges related to adopting them for modelling biological systems.

Biological Assembly Comparison with VAST

The VAST+ algorithm is an efficient, simple, and elegant solution to the problem of comparing the atomic structures of biological assemblies. Given two protein assemblies, it takes as input all the pairwise structural alignments of the component proteins. It then clusters the rotation matrices from the pairwise superpositions, with the clusters corresponding to subsets of the two assemblies that may be aligned and well superposed. It uses the Vector Alignment Search Tool (VAST) protein-protein comparison method for the input structural alignments, but other methods could be used, as well. From a chosen cluster, an "original" alignment for the assembly may be defined by simply combining the relevant input alignments. However, it is often useful to reduce/trim the original alignment, using a Monte Carlo refinement algorithm, which allows biologically relevant conformational differences to be more readily detected and observed. The method is easily extended to include RNA or DNA molecules. VAST+ results may be accessed via the URL https://www.ncbi.nlm.nih.gov/Structure , then entering a PDB accession or terms in the search box, and using the link [VAST+] in the upper right corner of the Structure Summary page.

DALI and the persistence of protein shape

DALI is a popular resource for comparing protein structures. The software is based on distance-matrix alignment. The associated web server provides tools to navigate, integrate and organize some data pushed out by genomics and structural genomics. The server has been running continuously for the past 25 years. Structural biologists routinely use DALI to compare a new structure against previously known protein structures. If significant similarities are discovered, it may indicate a distant homology, that is, that the structures are of shared origin. This may be significant in determining the molecular mechanisms, as these may remain very similar from a distant predecessor to the present day, for example, from the last common ancestor of humans and bacteria. Meta-analysis of independent reference-based evaluations of alignment accuracy and fold discrimination shows DALI at top rank in six out of 12 studies. The web server and standalone software are available from http://ekhidna2.biocenter.helsinki.fi/dali.

The cell cycle-regulated DNA adenine methyltransferase CcrM opens a bubble at its DNA recognition site

The Caulobacter crescentus cell cycle-regulated DNA methyltransferase (CcrM) methylates the adenine of hemimethylated GANTC after replication. Here we present the structure of CcrM in complex with double-stranded DNA containing the recognition sequence. CcrM contains an N-terminal methyltransferase domain and a C-terminal nonspecific DNA-binding domain. CcrM is a dimer, with each monomer contacting primarily one DNA strand: the methyltransferase domain of one molecule binds the target strand, recognizes the target sequence, and catalyzes methyl transfer, while the C-terminal domain of the second molecule binds the non-target strand. The DNA contacts at the 5-base pair recognition site results in dramatic DNA distortions including bending, unwinding and base flipping. The two DNA strands are pulled apart, creating a bubble comprising four recognized base pairs. The five bases of the target strand are recognized meticulously by stacking contacts, van der Waals interactions and specific Watson–Crick polar hydrogen bonds to ensure high enzymatic specificity.

CcrM is a cell cycle-regulated DNA methyltransferase that methylates an adenine within a specific sequence following replication in the gram negative bacterium Caulobacter crescentus. Here the authors present a crystal structure of DNA-bound CcrM that reveals the molecular mechanism leading to sequence-specific methylation.

Collection of antimicrobial peptides database and its derivatives: Applications and beyond

Collection of antimicrobial peptides (CAMP), CAMPSign, and ClassAMP are open-access resources that have been developed to enhance research on antimicrobial peptides (AMPs). Comprehensive information on AMPs and machine learning-based predictive models are made available for users through these resources. As of date, CAMP_R3 has 10,247 sequences, 757 structures, and 114 family-specific signatures of AMPs along with associated tools for AMP sequence and structure analysis. CAMPSign uses family-specific sequence conservation, in the form of patterns and hidden Markov models for identification of AMPs. ClassAMP can be used to classify AMPs as antibacterial, antifungal, or antiviral based on sequence information. Here we describe CAMP and its derivatives and illustrate, with a few examples, the contribution of these online resources to the advancement of our current understanding of AMPs.

A DNA-Binding Protein Tunes Septum Placement during Bacillus subtilis Sporulation

Bacillus subtilis is a bacterium capable of differentiating into a spore form more resistant to environmental stress. Early in sporulation, each cell possesses two copies of a circular chromosome. A polar FtsZ ring (Z ring) directs septation over one of the chromosomes, generating two cell compartments. The smaller "forespore" compartment initially contains only 25 to 30% of one chromosome, and this transient genetic asymmetry is required for differentiation. Timely assembly of polar Z rings and precise capture of the chromosome in the forespore both require the DNA-binding protein RefZ. To mediate its role in chromosome capture, RefZ must bind to specific DNA motifs (RBMs) that localize near the poles at the time of septation. Cells artificially induced to express RefZ during vegetative growth cannot assemble Z rings, an effect that also requires DNA binding. We hypothesized that RefZ-RBM complexes mediate precise chromosome capture by modulating FtsZ function. To investigate, we isolated 10 RefZ loss-of-function (rLOF) variants unable to inhibit cell division yet still capable of binding RBMs. Sporulating cells expressing the rLOF variants in place of wild-type RefZ phenocopied a ΔrefZ mutant, suggesting that RefZ acts through an FtsZ-dependent mechanism. The crystal structure of RefZ was solved, and wild-type RefZ and the rLOF variants were further characterized. Our data suggest that RefZ's oligomerization state and specificity for the RBMs are critical determinants influencing RefZ's ability to affect FtsZ dynamics. We propose that RBM-bound RefZ complexes function as a developmentally regulated nucleoid occlusion system for fine-tuning the position of the septum relative to the chromosome during sporulation.IMPORTANCE The bacterial nucleoid forms a large, highly organized structure. Thus, in addition to storing the genetic code, the nucleoid harbors positional information that can be leveraged by DNA-binding proteins to spatially constrain cellular activities. During B. subtilis sporulation, the nucleoid undergoes reorganization, and the cell division protein FtsZ assembles polarly to direct septation over one chromosome. The TetR family protein RefZ binds DNA motifs (RBMs) localized near the poles at the time of division and is required for both timely FtsZ assembly and precise capture of DNA in the future spore compartment. Our data suggest that RefZ exploits nucleoid organization by associating with polarly localized RBMs to modulate the positioning of FtsZ relative to the chromosome during sporulation.

Predicting Protein-Protein Interfaces that Bind Intrinsically Disordered Protein Regions

A long-standing goal in biology is the complete annotation of function and structure on all protein-protein interactions, a large fraction of which is mediated by intrinsically disordered protein regions (IDRs). However, knowledge derived from experimental structures of such protein complexes is disproportionately small due, in part, to challenges in studying interactions of IDRs. Here, we introduce IDRBind, a computational method that by combining gradient boosted trees and conditional random field models predicts binding sites of IDRs with performance approaching state-of-the-art globular interface predictions, making it suitable for proteome-wide applications. Although designed and trained with a focus on molecular recognition features, which are long interaction-mediating-elements in IDRs, IDRBind also predicts the binding sites of short peptides more accurately than existing specialized predictors. Consistent with IDRBind's specificity, a comparison of protein interface categories uncovered uniform trends in multiple physicochemical properties, positioning molecular recognition feature interfaces between peptide and globular interfaces.

Structural and immunological characterization of a new nucleotidyltransferase-like antigen from Paracoccidioides brasiliensis

Pb27 antigen is an interesting alternative to immunological diagnosis of Paracoccidioidomycosis (PCM) and has demonstrated to be protective in experimental PCM. Its tertiary structure and possible function remained unknown till now. To study Pb27 at the atomic level, the recombinant protein was expressed in Escherichia coli BL21(DE3), purified, and its three-dimensional structure was solved by X-ray crystallography. Based on this structure, we performed a residue correlation analysis and in silico ligand search assays to address a possible biological function to Pb27. We identified Pb27 as a member of the extensive nucleotidyltransferase superfamily. The protein has an αβαβαβ topology with two domains (N- and C-terminal domains) and adopts a monomeric form as its biological unit in solution. Structural comparisons with similar members of the superfamily clearly indicate Pb27 C-terminal domain is singular and may play an important role in its biological function. Bioinformatics analysis suggested that Pb27 might bind to ATP and CTP. This suggestion is corroborated by the fact that a magnesium cation is coordinated by two aspartic acid residues present at the active site (between N- and C-terminal domains), as evidenced by X-ray diffraction data. Besides, NMR assays (¹H-¹⁵N HSQC spectra) confirmed the binding of CTP to Pb27, demonstrating for the first time an interaction between a nucleotide and this protein. Moreover, we evaluated the reactivity of sera from patients with Paracoccidioides brasiliensis infection against the recombinant form of Pb27 and showed that it was recognized by sera from infected and treated patients. Predicted B and T cell epitopes were synthesized and further evaluated against sera of PCM patients, providing information of the most reactive peptides in Pb27 primary structure which interact with specific Pb27 antibodies.