Molecular and Evolutionary Characteristics of Chicken Parvovirus (ChPV) Genomes Detected in Chickens with Runting-Stunting Syndrome.

Chicken Parvovirus (ChPV) belongs to the genus Aveparvovirus and is implicated in enteric diseases like runting-stunting syndrome (RSS) in poultry. In RSS, chicken health is affected by diarrhea, depression, and increased mortality, causing significant economic losses in the poultry industry. This study aimed to characterize the ChPV genomes detected in chickens with RSS through a metagenomic approach and compare the molecular and evolutionary characteristics within the Aveparvovirus galliform1 species. The intestinal content of broiler flocks affected with RSS was submitted to viral metagenomics. The assembled prevalent genomes were identified as ChPV after sequence and phylogenetic analysis, which consistently clustered separately from Turkey Parvovirus (TuPV). The strain USP-574-A presented signs of genomic recombination. The selective pressure analysis indicated that most of the coding genes in A. galliform1 are evolving under diversifying (negative) selection. Protein modeling of ChPV and TuPV viral capsids identified high conservancy over the VP2 region. The prediction of epitopes identified several co-localized antigenic peptides from ChPV and TuPV, especially for T-cell epitopes, highlighting the immunological significance of these sites. However, most of these peptides presented host-specific variability, obeying an adaptive scenario. The results of this study show the evolutionary path of ChPV and TuPV, which are influenced by diversifying events such as genomic recombination and selective pressure, as well as by adaptation processes, and their subsequent immunological impact.

Exploring the Structurally Conserved Regions and Functional Significance in Bacterial N-Terminal Nucleophile (Ntn) Amide-Hydrolases.

The initial adoption of penicillin as an antibiotic marked the start of exploring other compounds essential for pharmaceuticals, yet resistance to penicillins and their side effects has compromised their efficacy. The N-terminal nucleophile (Ntn) amide-hydrolases S45 family plays a key role in catalyzing amide bond hydrolysis in various compounds, including antibiotics like penicillin and cephalosporin. This study comprehensively analyzes the structural and functional traits of the bacterial N-terminal nucleophile (Ntn) amide-hydrolases S45 family, covering penicillin G acylases, cephalosporin acylases, and D-succinylase. Utilizing structural bioinformatics tools and sequence analysis, the investigation delineates structurally conserved regions (SCRs) and substrate binding site variations among these enzymes. Notably, sixteen SCRs crucial for substrate interaction are identified solely through sequence analysis, emphasizing the significance of sequence data in characterizing functionally relevant regions. These findings introduce a novel approach for identifying targets to enhance the biocatalytic properties of N-terminal nucleophile (Ntn) amide-hydrolases, while facilitating the development of more accurate three-dimensional models, particularly for enzymes lacking structural data. Overall, this research advances our understanding of structure-function relationships in bacterial N-terminal nucleophile (Ntn) amide-hydrolases, providing insights into strategies for optimizing their enzymatic capabilities.

Structural Insights into the Giardia lamblia Target of Rapamycin Homolog: A Bioinformatics Approach.

TOR proteins, also known as targets of rapamycin, are serine/threonine kinases involved in various signaling pathways that regulate cell growth. The protozoan parasite Giardia lamblia is the causative agent of giardiasis, a neglected infectious disease in humans. In this study, we used a bioinformatics approach to examine the structural features of GTOR, a G. lamblia TOR-like protein, and predict functional associations. Our findings confirmed that it shares significant similarities with functional TOR kinases, including a binding domain for the FKBP-rapamycin complex and a kinase domain resembling that of phosphatidylinositol 3-kinase-related kinases. In addition, it can form multiprotein complexes such as TORC1 and TORC2. These results provide valuable insights into the structure-function relationship of GTOR, highlighting its potential as a molecular target for controlling G. lamblia cell proliferation. Furthermore, our study represents a step toward rational drug design for specific anti-giardiasis therapeutic agents.

Structural Modeling of Nanobodies: A Benchmark of State-of-the-Art Artificial Intelligence Programs.

The number of applications for nanobodies is steadily expanding, positioning these molecules as fast-growing biologic products in the biotechnology market. Several of their applications require protein engineering, which in turn would greatly benefit from having a reliable structural model of the nanobody of interest. However, as with antibodies, the structural modeling of nanobodies is still a challenge. With the rise of artificial intelligence (AI), several methods have been developed in recent years that attempt to solve the problem of protein modeling. In this study, we have compared the performance in nanobody modeling of several state-of-the-art AI-based programs, either designed for general protein modeling, such as AlphaFold2, OmegaFold, ESMFold, and Yang-Server, or specifically designed for antibody modeling, such as IgFold, and Nanonet. While all these programs performed rather well in constructing the nanobody framework and CDRs 1 and 2, modeling CDR3 still represents a big challenge. Interestingly, tailoring an AI method for antibody modeling does not necessarily translate into better results for nanobodies.

Functional characterization of the p.(Gln195His) or Tainan and novel p.(Ser184Cys) or Toluca glucose-6-phosphate dehydrogenase (G6PD) gene natural variants identified through Mexican newborn screening for glucose-6-phosphate dehydrogenase deficiency.

BACKGROUND: Newborn screening for glucose-6-phosphate dehydrogenase deficiency (G6PDd) was implemented in Mexico beginning in 2017. In a Mexican population, genotyping analysis of G6PD as a second-tier method identified a previously unreported missense variant, p.(Ser184Cys), which we propose to call "Toluca", and the extremely rare p.(Gln195His) or "Tainan" variant, which was previously described in the Taiwanese population as a Class II allele through in silico evaluations. Here, we sought to perform in vitro biochemical characterizations of the Toluca and Tainan G6PD natural variants and describe their associated phenotypes. METHODS: The "Toluca" and "Tainan" variants were identified in three unrelated G6PDd newborn males, two of whom lacked evidence of acute hemolytic anemia (AHA) or neonatal hyperbilirubinemia (NHB). We constructed wild-type (WT), Tainan, and Toluca G6PD recombinant enzymes and performed in vitro assessments. RESULTS: Both variants had diminished G6PD expression, decreased affinities for glucose-6-phosphate and NADP+ substrates, significant decreases in catalytic efficiency (∼97 % with respect to WT-G6PD), and diminished thermostabilities that were partially rescued by NADP+. In silico protein modeling predicted that the variants would have destabilizing effects on the protein tertiary structure, potentially reducing the enzyme half-lives and/or catalytic efficiencies. CONCLUSION: Our data suggest that G6PD "Tainan" and "Toluca" are potential Class II natural variants, which agrees with the absence of chronic nonspherocytic hemolytic anemia (CNSHA) in our patients. It remains to be determined whether these variants represent high-risk genetic factors for developing CNSHA, AHA, and/or NHB.

Complete Genome Characterization of Reticuloendotheliosis Virus Detected in Chickens with Multiple Viral Coinfections.

Reticuloendotheliosis virus (REV) is a retroviral pathogen capable of infecting several avian hosts and is associated with immunosuppression, anemia, proventriculitis, neoplasia, and runting-stunting syndrome. Its genome contains the three major genes, gag, pol, and env, and two flanking long terminal repeat (LTR) regions. Complete genome sequences of REV are limited in terms of geographical origin. The aim of this study was to characterize the complete genome of REV detected in Brazilian chickens with multiple viral coinfections and analyze the polymorphisms in the deduced amino acids sequences corresponding to its encoded proteins. We tested the presence and completeness of REV as well as other viral pathogens in samples from Brazilian poultry farms by qPCR. The complete genomes of two REV strains were sequenced by overlapping fragments through the dideoxy method. Phylogenetic analysis, pairwise identity matrix, polymorphism identification and protein modeling were performed along the entire genome. We detected REV in 65% (26/40) of the tested samples. Concomitant viral infections were detected in 82.5% (33/40) of the samples and in 90% (9/10) of the farms. Multiple infections included up to seven viruses. Phylogenetic analysis classified both Brazilian strains into REV subtype 3, and the pairwise comparison indicated that strains from the USA and fowlpox virus (FWPV)-related strains were the most identical. The subdomain p18 in gag, the reverse transcriptase/ribonuclease H in pol, and the surface (SU) in the env protein were the most polymorphic in genomic comparisons. The relevant motifs for each protein were highly conserved, with fewer polymorphisms in the fusion peptide, immunosuppression domain, and disulfide bonds on the surface (SU) and transmembrane (TM) of env. This is the first study to include complete genomes of REV in Brazil and South America detected in farms with multiple viral coinfections. Our findings suggest an involvement of REV as an immunosuppressor and active agent in the emergence and progression of multiple infectious diseases. We also found a possible etiological relationship between Brazilian strains and the USA and FWPV recombinant strains. This information highlights the need for epidemiological vigilance regarding REV in association with another pathogens.

Structural Insights into the Substrate Transport Mechanisms in GTR Transporters through Ensemble Docking.

Glucosinolate transporters (GTRs) are part of the nitrate/peptide transporter (NPF) family, members of which also transport specialized secondary metabolites as substrates. Glucosinolates are defense compounds derived from amino acids. We selected 4-methylthiobutyl (4MTB) and indol-3-ylmethyl (I3M) glucosinolates to study how GTR1 from Arabidopsis thaliana transports these substrates in computational simulation approaches. The designed pipeline reported here includes massive docking of 4MTB and I3M in an ensemble of GTR1 conformations (in both inward and outward conformations) extracted from molecular dynamics simulations, followed by clustered and substrate-protein interactions profiling. The identified key residues were mutated, and their role in substrate transport was tested. We were able to identify key residues that integrate a major binding site of these substrates, which is critical for transport activity. In silico approaches employed here represent a breakthrough in the plant transportomics field, as the identification of key residues usually takes a long time if performed from a purely wet-lab experimental perspective. The inclusion of structural bioinformatics in the analyses of plant transporters significantly speeds up the knowledge-gaining process and optimizes valuable time and resources.

Bioinformatic Analysis of Two TOR (Target of Rapamycin)-Like Proteins Encoded by Entamoeba histolytica Revealed Structural Similarities with Functional Homologs.

The target of rapamycin (TOR), also known as FKBP-rapamycin associated protein (FRAP), is a protein kinase belonging to the PIKK (phosphatidylinositol 3-kinase (PI3K)-related kinases) family. TOR kinases are involved in several signaling pathways that control cell growth and proliferation. Entamoeba histolytica, the protozoan parasite that causes human amoebiasis, contains two genes encoding TOR-like proteins: EhFRAP and EhTOR2. To assess their potential as drug targets to control the cell proliferation of E. histolytica, we studied the structural features of EhFRAP and EhTOR2 using a biocomputational approach. The overall results confirmed that both TOR amoebic homologs share structural similarities with functional TOR kinases, and show inherent abilities to form TORC complexes and participate in protein-protein interaction networks. To our knowledge, this study represents the first in silico characterization of the structure-function relationships of EhFRAP and EhTOR2.

Further Evidence That Defects in Main Thyroid Dysgenesis-Related Genes Are an Uncommon Etiology for Primary Congenital Hypothyroidism in Mexican Patients: Report of Rare Variants in FOXE1, NKX2-5 and TSHR.

Mexico shows a high birth prevalence of congenital hypothyroidism (CH) due to thyroid dysgenesis (TD). PAX8 defects underlie only 1% of these cases and NKX2-1 does not seem to be involved. Here, we analyzed other TD-related genes in 128 non-related Mexican patients (females 77.3%; 6 months to 16.6 years) with non-syndromic CH-TD diagnosis established by clinical evaluation, thyroid hormone serum profiling, and scintigraphy (74%) or ultrasonography (26%). We performed Sanger sequencing of FOXE1, NKX2-5, and TSHR and evaluated copy number variations (CNVs) in TSHR, FOXE1, PAX8, and NKX2-1 by multiplex ligation-dependent probe amplification. Odds ratios for TD risk were explored for FOXE1 polyalanine stretches [polyAla-rs71369530] in cases and controls (N = 116). Five rare missense changes cataloged as benign (NKX2-5:p.(Ala119Ser)-rs137852684), of unknown significance (FOXE1:p.(Ala335Gly)-rs543372757; TSHR:p.(Asp118Asn)-rs1414102266), and likely pathogenic (FOXE1:p.(Gly124Arg)-rs774035532; TSHR:p.(Trp422Arg)-rs746029360) accounted for 1.5% (N = 2/128) of clinically relevant genotypes (supported in part by protein modeling) in CH-TD. No CNVs were identified, nor did polyAla > 14 alanines in FOXE1 significantly protect against TD. The present and previously published data collectively show that small clinically relevant germline variants in PAX8, FOXE1, and TSHR are found in only a very small proportion (2.5%) of isolated CH-TD Mexican patients.

Structural bases that underline Trypanosoma cruzi calreticulin proinfective, antiangiogenic and antitumor properties.

Microbes have developed mechanisms to resist the host immune defenses and some elicit antitumor immune responses. About 6 million people are infected with Trypanosoma cruzi, the protozoan agent of Chagas' disease, the sixth neglected tropical disease worldwide. Eighty years ago, G. Roskin and N. Klyuyeva proposed that T. cruzi infection mediates an anti-cancer activity. This observation has been reproduced by several other laboratories, but no molecular basis has been proposed. We have shown that the highly pleiotropic chaperone calreticulin (TcCalr, formerly known as TcCRT), translocates from the parasite ER to the exterior, where it mediates infection. Similar to its human counterpart HuCALR (formerly known as HuCRT), TcCalr inhibits C1 in its capacity to initiate the classical pathway of complement activation. We have also proposed that TcCalr inhibits angiogenesis and it is a likely mediator of antitumor effects. We have generated several in silico structural TcCalr models to delimit a peptide (VC-TcCalr) at the TcCalr N-domain. Chemically synthesized VC-TcCalr did bind to C1q and was anti-angiogenic in Gallus gallus chorioallantoic membrane assays. These properties were associated with structural features, as determined in silico. VC-TcCalr, a strong dipole, interacts with charged proteins such as collagen-like tails and scavenger receptors. Comparatively, HuCALR has less polarity and spatial stability, probably due to at least substitutions of Gln for Gly, Arg for Lys, Arg for Asp and Ser for Arg that hinder protein-protein interactions. These differences can explain, at least in part, how TcCalr inhibits the complement activation pathway and has higher efficiency as an antiangiogenic and antitumor agent than HuCALR.

Structural determinant of functionality in acyl lipid desaturases.

Little is known about the structure-function relationship of membrane-bound lipid desaturases. Using a domain-swapping strategy, we found that the N terminus (comprising the two first transmembrane segments) region of Bacillus cereus DesA desaturase improves Bacillus subtilis Des activity. In addition, the replacement of the first two transmembrane domains from Bacillus licheniformis inactive open reading frame (ORF) BL02692 with the corresponding domain from DesA was sufficient to resurrect this enzyme. Unexpectedly, we were able to restore the activity of ORF BL02692 with a single substitution (Cys40Tyr) of a cysteine localized in the first transmembrane domain close to the lipid-water interface. Substitution of eight residues (Gly90, Trp104, Lys172, His228, Pro257, Leu275, Tyr282, and Leu284) by site-directed mutagenesis produced inactive variants of DesA. Homology modeling of DesA revealed that His228 is part of the metal binding center, together with the canonical His boxes. Trp104 shapes the hydrophobic tunnel, whereas Gly90 and Lys172 are probably involved in substrate binding/recognition. Pro257, Leu275, Tyr282, and Leu284 might be relevant for the structural arrangement of the active site or interaction with electron donors. This study reveals the role of the N-terminal region of Δ5 phospholipid desaturases and the individual residues necessary for the activity of this class of enzymes.

Compound heterozygous KCNQ1 mutations (A300T/P535T) in a child with sudden unexplained death: Insights into possible molecular mechanisms based on protein modeling.

Sudden death in a child is a devastating event with important medical implications for surviving relatives. Because it may be the first manifestation of unknown inherited cardiac disease, molecular autopsy can be helpful to determine the cause of death and identify at risk family members. The aim of the study was to perform a molecular autopsy in a seven year-old girl with sudden unexplained death, to find evidence supporting the possible pathogenicity of mutations identified in inherited cardiac disease genes, and to clinically and genetically assess first-degree relatives. DNA from the index case was extracted from umbilical cord cells stored at birth, and DNA of first-degree relatives from blood samples. Targeted sequencing was performed using a Haloplex design including 81 cardiogenes. Possible functional consequences of the mutations were analyzed using protein modeling and structural mobility analyses. The child was compound heterozygous for KCNQ1 variants p.Ala300Thr and p.Pro535Thr. Ala300Thr is known to cause long QT syndrome in the homozygous state, while Pro535Thr is novel and of unknown clinical significance. The father and sibling were Ala300Thr heterozygous, and had normal QTc intervals at rest and during exercise. The asymptomatic mother was heterozygous for Pro535Thr, and showed borderline QTc at rest, but prolonged QTc during exercise. Protein modeling predicted that Ala300Thr alters the mobility profile of the Kv7.1 tetramer and Thr535 disrupts a calmodulin-binding site, probably causing co-assembly or trafficking defects of the mutant monomer. Altogether, the evidence strongly suggests that this child was affected with a recessive form of Romano Ward syndrome.

A novel pH-stable, endoglucanase (JqCel5A) isolated from a salt-lake microorganism, Jonesia quinghaiensis

Background: Endoglucanase, one of three type cellulases, can randomly cleave internal p-1,4-linkages in cellulose polymers. Thus, it could be applied in agricultural and industrial processes. Results: A novel endoglucanase gene (JqCel5A) was cloned from Jonesia quinghaiensis and functionally expressed in Escherichia coli Rosetta (DE3). It contained 1722 bp and encoded a 573-residue polypeptide consisting of a catalytic domain of glycoside hydrolase family 5 (GH5) and a type 2 carbohydrate-binding module (CBM2), together with a predicted molecular mass of 61.79 kD. The purified JqCel5A displayed maximum activity at 55°C and pH 7.0, with 21.7 U/mg, 26.19 U/mg and 4.81 U/mg towards the substrate carboxymethyl cellulose, barley glucan and filter paper, respectively. Interestingly, JqCel5A exhibited high pH stability over a broad pH range of pH (3-11), and had good tolerance to a wide variety of deleterious chemicals including heavy metals and detergent. The catalytic mechanism of JqCel5A was also investigated by site mutagenesis and homology-modeling in this study. Conclusions: It was believed that these properties might make JqCel5A to be potentially used in the suitable industrial catalytic condition, which has a broad pH fluctuation and/or chemical disturbance.

Comparative modeling of proteins: a method for engaging students' interest in bioinformatics tools.

The huge increase in data being produced in the genomic era has produced a need to incorporate computers into the research process. Sequence generation, its subsequent storage, interpretation, and analysis are now entirely computer-dependent tasks. Universities from all over the world have been challenged to seek a way of encouraging students to incorporate computational and bioinformatics skills since undergraduation in order to understand biological processes. The aim of this article is to report the experience of awakening students' interest in bioinformatics tools during a course focused on comparative modeling of proteins. The authors start by giving a full description of the course environmental context and students' backgrounds. Then they detail each class and present a general overview of the protein modeling protocol. The positive and negative aspects of the course are also reported, and some of the results generated in class and in projects outside the classroom are discussed. In the last section of the article, general perspectives about the course from students' point of view are given. This work can serve as a guide for professors who teach subjects for which bioinformatics tools are useful and for universities that plan to incorporate bioinformatics into the curriculum.

Activity of chalcones derived from 2,4,5-trimethoxybenzaldehyde against Meloidogyne exigua and in silico interaction of one chalcone with a putative caffeic acid 3-O-methyltransferase from Meloidogyne incognita.

Meloidogyne exigua is a parasitic nematode of plants that causes great losses to coffee farmers. In an effort to develop parasitic controls, 154 chalcones were synthesized and screened for activity against this nematode. The best results were obtained with (2E)-1-(4'-nitrophenyl)-3-(2,4,5-trimethoxyphenyl)prop-2-en-1-one (6) with a 50% lethal concentration (LC50) of 171 µg/ml against M. exigua second-stage juveniles, in comparison to the commercially-available nematicide carbofuran which had an LC50 of 260 µg/ml under the same conditions. When coffee plants were used, 6 reduced the nematode population to ~50% of that observed in control plants. To investigate the mechanism of action of 6, an in silico study was carried out, which indicated that 6 may act against M. exigua through inhibition of a putative caffeic acid 3-O-methyltransferase homodimer, the amino acid sequence of which was determined by examining the genome of Meloidogyne incognita.