Identification and Modeling of a GT-A Fold in the α-Dystroglycan Glycosylating Enzyme LARGE1.

The acetylglucosaminyltransferase-like protein LARGE1 is an enzyme that is responsible for the final steps of the post-translational modifications of dystroglycan (DG), a membrane receptor that links the cytoskeleton with the extracellular matrix in the skeletal muscle and in a variety of other tissues. LARGE1 acts by adding the repeating disaccharide unit [-3Xyl-α1,3GlcAß1-] to the extracellular portion of the DG complex (α-DG); defects in the LARGE1 gene result in an aberrant glycosylation of α-DG and consequent impairment of its binding to laminin, eventually affecting the connection between the cell and the extracellular environment. In the skeletal muscle, this leads to degeneration of the muscular tissue and muscular dystrophy. So far, a few missense mutations have been identified within the LARGE1 protein and linked to congenital muscular dystrophy, and because no structural information is available on this enzyme, our understanding of the molecular mechanisms underlying these pathologies is still very limited. Here, we generated a 3D model structure of the two catalytic domains of LARGE1, combining different molecular modeling approaches. Furthermore, by using molecular dynamics simulations, we analyzed the effect on the structure and stability of the first catalytic domain of the pathological missense mutation S331F that gives rise to a severe form of muscle-eye-brain disease.

Identification and in silico characterization of a novel PKLR genotype in a Turkish newborn.

Pyruvate kinase deficiency (PKD) is the most common glycolytic defect leading to chronic nonspherocytic hemolytic anemia (CNSHA). Clinical manifestations of PKD reflect the symptoms and complications of the chronic hemolysis, including anemia, jaundice, bilirubin gallstones due to hyperbilirubinemia, splenomegaly and iron overload. In this study, we report the finding of a 5-months-old Turkish male newborn with moderate CNSHA and PKD. Mutation screening of Pyruvate Kinase Liver/Red (PKLR) gene revealed that the patient carried the known pathogenic variant (PV) c.1456C > T (p.Arg486Trp) and an unreported variant c.1067T > G (p.Met356Arg). Computational variant analysis (CVA) highlighted the deleterious structural effects on the mutant PK enzyme, suggesting its pathogenic role. In this patient, the molecular evaluation of PKD, that allowed the identification of the novel PKLR genotype, coupled with CVA led to the definitive and correct diagnosis of CNSHA.

The hemoglobin Gly16ß1Asp polymorphism in turbot (Scophthalmus maximus) is differentially distributed across European populations.

Turbot is an important flatfish widely distributed along the European coasts, whose fishery is centered in the North Sea. The commercial value of the species has boosted a successful aquaculture sector in Europe and China. Body growth is the main target of turbot breeding programs and is also a key trait related to local adaptation to temperature and salinity. Differences in growth rate and optimal growth temperature in turbot have been shown to be associated with a hemoglobin polymorphism reported more than 50 years ago. Here, we identified a Gly16Asp amino acid substitution in the ß1 globin subunit by searching for genetic variation in the five functional globin genes within the whole annotated turbot genome. We predicted increased stability of the turbot hemoglobin by the replacement of the conserved Gly with the negative charged Asp residue that is consistent with the higher rate of αß dimer assembly in the human J-Baltimore Gly16ß->Asp mutant than in normal HbA. The turbot Hbß1-Gly16 variant dominated in the northern populations examined, particularly in the Baltic Sea, while the Asp allele showed elevated frequencies in southern populations and was the prevalent variant in the Adriatic Sea. Body weight did not associate with the Hbß1 genotypes at farming conditions (i.e., high oxygen levels, feeding ad libitum) after analyzing 90 fish with high growth dispersal from nine turbot families. Nevertheless, all data at hand suggest that the turbot hemoglobin polymorphism has an adaptive significance in the variable wild conditions regarding temperature and oxygen availability.

Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod.

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4-Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC-specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy-Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6-2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long-distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co-occurrence of multiple adaptive alleles forming a 'supergene' in the NEAC population.

Structural studies and SH3 domain binding properties of a human antiviral salivary proline-rich peptide.

Human saliva contains hundreds of small proline-rich peptides originated by the proteolytic cleavage of the salivary basic Proline-Rich Proteins. Nevertheless only for few of them a specific biological activity has been assigned to date. Among them, the 1932 Da peptide (p1932) has been patented as an anti-HIV agent. In order to shed light on the possible mechanism of action of this peptide, we assessed in this study, by means of molecular dynamics calculations, circular dichroism and FTIR spectroscopic techniques, that p1932 has an intrinsic propensity to adopt a polyproline-II helix arrangement. This structural feature combined with the presence of PxxP motifs in its primary structure, represents an essential property for the exploitation of several biological activities. Next to these findings, we recently demonstrated the ability of this peptide to be internalized within cells of the oral mucosa, thus we focused onto a possible intracellular target, represented by the SH3 domains family. Its ability to interact with selected SH3 domains was finally assayed by Surface Plasmon Resonance spectroscopy. As a result, only Fyn, Hck, and c-Src SH3 domains gave positive results in terms of interaction, showing dissociation constants ranging from nanomolar to micromolar values having the best performer a KD of 148 nM. It is noteworthy that all the interacting domains belong to the Src kinases family, suggesting a role for p1932 as a modulator of the signal transduction pathways mediated by these kinases. © 2016 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 714-725, 2016.

The conserved Phe GH5 of importance for hemoglobin intersubunit contact is mutated in gadoid fish.

BACKGROUND: Functionality of the tetrameric hemoglobin molecule seems to be determined by a few amino acids located in key positions. Oxygen binding encompasses structural changes at the interfaces between the α1ß2 and α2ß1 dimers, but also subunit interactions are important for the oxygen binding affinity and stability. The latter packing contacts include the conserved Arg B12 interacting with Phe GH5, which is replaced by Leu and Tyr in the αA and αD chains, respectively, of birds and reptiles. RESULTS: Searching all known hemoglobins from a variety of gnathostome species (jawed vertebrates) revealed the almost invariant Arg B12 coded by the AGG triplet positioned at an exon-intron boundary. Rare substitutions of Arg B12 in the gnathostome ß globins were found in pig, tree shrew and scaled reptiles. Phe GH5 is also highly conserved in the ß globins, except for the Leu replacement in the ß1 globin of five marine gadoid species, gilthead seabream and the Comoran coelacanth, while Cys and Ile were found in burbot and yellow croaker, respectively. Atlantic cod ß1 globin showed a Leu/Met polymorphism at position GH5 dominated by the Met variant in northwest-Atlantic populations that was rarely found in northeast-Atlantic cod. Site-specific analyses identified six consensus codons under positive selection, including 122ß(GH5), indicating that the amino acid changes identified at this position may offer an adaptive advantage. In fact, computational mutation analysis showed that the replacement of Phe GH5 with Leu or Cys decreased the number of van der Waals contacts essentially in the deoxy form that probably causes a slight increase in the oxygen binding affinity. CONCLUSIONS: The almost invariant Arg B12 and the AGG codon seem to be important for the packing contacts and pre-mRNA processing, respectively, but the rare mutations identified might be beneficial. The Leu122ß1(GH5)Met and Met55ß1(D6)Val polymorphisms in Atlantic cod hemoglobin modify the intradimer contacts B12-GH5 and H2-D6, while amino acid replacements at these positions in avian hemoglobin seem to be evolutionary adaptive in air-breathing vertebrates. The results support the theory that adaptive changes in hemoglobin functions are caused by a few substitutions at key positions.

Molecular interactions of hemoglobin with resveratrol: potential protective antioxidant role and metabolic adaptations of the erythrocyte.

This article reports the role of resveratrol in the erythrocyte as a result of its interaction with hemoglobin and describes the effect of this interaction on the metabolism, the redox state, and the release of ATP. The drug crosses the erythrocyte membrane and binds to hemoglobin, altering its modulation and the release of ATP. Our study correlates the variation of the phosphorylation balance induced by resveratrol with the change in the intracellular concentration of ATP and with the decrease in ATP release from red blood cell and the consequent paracrine alteration on the vascular epithelium. Molecular docking calculations indicate larger specificity of binding for oxy-hemoglobin that correlates well with the stabilization of the R-quaternary structure and with the functional modulation of resveratrol on the protein. Finally, we locate a putative binding site at the central cavity of hemoglobin and characterize its key interacting residues with the drug. Computational results support the assumption that resveratrol may act as a protector agent against oxidative protein damage by interacting with hemoglobin.

Top-down HPLC-ESI-MS characterization of rat gliadoralin A, a new member of the family of rat submandibular gland glutamine-rich proteins and potential substrate of transglutaminase.

During HPLC-ESI-MS/MS analysis of rat submandibular saliva secreted under isoprenaline stimulation, a protein with an experimental [M+H](1+) = 10,544.24 m/z was detected (17.5 ± 0.7 min). The MS/MS fragmentation pattern, manually investigated, allowed establishing an internal sequence in agreement with a DNA-derived sequence of an unknown rat protein coded D3Z9M3 (Swiss-Prot). To match the experimental MS/MS fragmentation pattern and protein mass with theoretical data, the removal from the N terminus of the signal peptide and from the C terminus of three amino acid (a.a.) residues (Arg-Ala-Val) and the cyclization of the N-terminal glutamine in pyroglutamic had to be supposed, resulting in a mature protein of 90 a.a. HPLC-ESI-MS/MS of the trypsin digest ensured 100% sequence coverage. For the high glutamine content (34/90 = 37.8%) we propose to name this protein rat gliadoralin A 1-90. Low amounts of five different isoforms were sporadically detected, which did not significantly change their relative amounts after stimulation. Gliadoralin A is substrate for transglutaminase-2, having Lys 60 and different Gln residues as major determinants for enzyme recognition. In silico investigation of superior structures evidenced that a small part of the protein adopts an α-helical fold, whereas large segments are unfolded, suggesting an unordered conformation.

A Green Lipophilization Reaction of a Natural Antioxidant.

A natural antioxidant, widely spread in plants, chlorogenic acid (CGA), can be lipophilized through a heterogeneous, non-enzymatic, catalytic process. Thus, sulfonic resins under no solvent conditions allow to obtain a series of esters in up to 93% yield through reaction of CGA with fatty alcohols of different chain length. The reaction takes place in one single step under mild conditions with conversions up to 96% and selectivity up to 99%. Product recovery in high purity was very easy and the esters obtained were fully characterized with spectroscopic techniques and through the DPPH test to verify the preservation of antioxidant activity. According to this test, all of them showed increased activity with respect to the parent acid and anyway higher than butylated hydroxyanisole. An in-silico method also suggested their very low toxicity. The increased lipophilicity of the esters allows their formulation in cosmetic and nutraceutic lipid-based products.

Prediction of CD44 Structure by Deep Learning-Based Protein Modeling.

CD44 is a cell surface glycoprotein transmembrane receptor that is involved in cell-cell and cell-matrix interactions. It crucially associates with several molecules composing the extracellular matrix, the main one of which is hyaluronic acid. It is ubiquitously expressed in various types of cells and is involved in the regulation of important signaling pathways, thus playing a key role in several physiological and pathological processes. Structural information about CD44 is, therefore, fundamental for understanding the mechanism of action of this receptor and developing effective treatments against its aberrant expression and dysregulation frequently associated with pathological conditions. To date, only the structure of the hyaluronan-binding domain (HABD) of CD44 has been experimentally determined. To elucidate the nature of CD44s, the most frequently expressed isoform, we employed the recently developed deep-learning-based tools D-I-TASSER, AlphaFold2, and RoseTTAFold for an initial structural prediction of the full-length receptor, accompanied by molecular dynamics simulations on the most promising model. All three approaches correctly predicted the HABD, with AlphaFold2 outperforming D-I-TASSER and RoseTTAFold in the structural comparison with the crystallographic HABD structure and confidence in predicting the transmembrane helix. Low confidence regions were also predicted, which largely corresponded to the disordered regions of CD44s. These regions allow the receptor to perform its unconventional activity.

Virtual screening and molecular dynamics simulations provide insight into repurposing drugs against SARS-CoV-2 variants Spike protein/ACE2 interface.

After over two years of living with Covid-19 and hundreds of million cases worldwide there is still an unmet need to find proper treatments for the novel coronavirus, due also to the rapid mutation of its genome. In this context, a drug repositioning study has been performed, using in silico tools targeting Delta Spike protein/ACE2 interface. To this aim, it has been virtually screened a library composed by 4388 approved drugs through a deep learning-based QSAR model to identify protein-protein interactions modulators for molecular docking against Spike receptor binding domain (RBD). Binding energies of predicted complexes were calculated by Molecular Mechanics/Generalized Born Surface Area from docking and molecular dynamics simulations. Four out of the top twenty ranking compounds showed stable binding modes on Delta Spike RBD and were evaluated also for their effectiveness against Omicron. Among them an antihistaminic drug, fexofenadine, revealed very low binding energy, stable complex, and interesting interactions with Delta Spike RBD. Several antihistaminic drugs were found to exhibit direct antiviral activity against SARS-CoV-2 in vitro, and their mechanisms of action is still debated. This study not only highlights the potential of our computational methodology for a rapid screening of variant-specific drugs, but also represents a further tool for investigating properties and mechanisms of selected drugs.

Non-synonymous variation and protein structure of candidate genes associated with selection in farm and wild populations of turbot (Scophthalmus maximus).

Non-synonymous variation (NSV) of protein coding genes represents raw material for selection to improve adaptation to the diverse environmental scenarios in wild and livestock populations. Many aquatic species face variations in temperature, salinity and biological factors throughout their distribution range that is reflected by the presence of allelic clines or local adaptation. The turbot (Scophthalmus maximus) is a flatfish of great commercial value with a flourishing aquaculture which has promoted the development of genomic resources. In this study, we developed the first atlas of NSVs in the turbot genome by resequencing 10 individuals from Northeast Atlantic Ocean. More than 50,000 NSVs where detected in the ~ 21,500 coding genes of the turbot genome, and we selected 18 NSVs to be genotyped using a single Mass ARRAY multiplex on 13 wild populations and three turbot farms. We detected signals of divergent selection on several genes related to growth, circadian rhythms, osmoregulation and oxygen binding in the different scenarios evaluated. Furthermore, we explored the impact of NSVs identified on the 3D structure and functional relationship of the correspondent proteins. In summary, our study provides a strategy to identify NSVs in species with consistently annotated and assembled genomes to ascertain their role in adaptation.

A Modified Basophil Activation Test for the Clinical Management of Immediate Hypersensitivity Reactions to Paclitaxel: A Proof-of-Concept Study.

Immediate hypersensitivity reactions (iHSRs) to taxanes are observed in 6% and 4% of gynecologic and breast cancer patients, respectively. Drug desensitization is the only option, as no comparable alternative therapy is available. Surfactants in the taxane formulation have been implicated in the immunopathogenesis of iHSRs, although sporadic skin test (ST) positivity and iHSRs to nab-paclitaxel have suggested the involvement of the taxane moiety and/or IgE-mediated pathomechanisms. In vitro diagnostic tests might offer insights into mechanisms underlying iHSRs to taxanes. The aim of the present study was to address this unmet need by developing a novel basophil activation test (BAT). The study included patients (n = 31) undergoing paclitaxel/carboplatin therapy. Seventeen patients presented with iHSRs to paclitaxel (iHSR-Taxpos), and eleven were tolerant (iHSR-Taxneg). Fourteen patients presented with iHSRs to carboplatin (iHSR-Plpos), and fourteen were tolerant (iHSR-Plneg). The BAT median stimulation index (SI) values were 1.563 (range, 0.02-4.11; n = 11) and -0.28 (range -4.88-0.07, n = 11) in iHSR-Taxpos and iHSR-Taxneg, respectively. The BAT median SI values were 4.45 (range, 0.1-26.7; n = 14) and 0 (range, -0.51-1.65; n = 12) in iHSR-Plpos and iHSR-Plneg, respectively. SI levels were not associated with iHSR severity grading. Comparing BAT results in iHSR-Taxpos and iHSR-Taxneg showed the area under the receiver operator characteristic (ROC) curve to be 0.9752 (p = 0.0002). The cutoff calculated by the maximized likelihood ratio identified 90.91% of iHSR-Taxpos patients and 90.91% of iHSR-Taxneg patients. Comparing BAT results for iHSR-Plpos and iHSR-Plneg showed the area under the ROC curve to be 0.9286 (p = 0.0002). The cutoff calculated by the maximized likelihood ratio identified 78.57% of iHSR-Plpos patients and 91.67% of iHSR-Plneg patients. Most iHSR-Taxpos patients for which ST was available (10/11) scored ST-negative and BAT-positive, whereas most iHSR-Plpos patients for which ST was available (14/14) scored both BAT- and ST-positive. This suggested the intervention of non-IgE-mediated mechanisms in iHSR-Taxpos patients. Consistent with this view, an in silico molecular docking analysis predicted the high affinity of paclitaxel to the degranulation-competent MRGPRX2 receptor. This hypothesis warrants further in vitro investigations. In conclusion, the present study provides preliminary proof-of-concept evidence that this novel BAT has potential utility in understanding mechanisms underlying iHSRs to taxanes.

A TLR/CD44 axis regulates T cell trafficking in experimental and human multiple sclerosis.

In the pathogenesis of autoimmune disorders, the modulation of leukocytes' trafficking plays a central role, still poorly understood. Here, we focused on the effect of TLR2 ligands in trafficking of T helper cells through reshuffling of CD44 isoforms repertoire. Concurrently, strain background and TLR2 haplotype affected Wnt/ß-catenin signaling pathway and expression of splicing factors. During EAE, mCD44 v9- v 10 was specifically enriched in the forebrain and showed an increased ability to bind stably to osteopontin. Similarly, we observed that hCD44 v7 was highly enriched in cells of cerebrospinal fluid from MS patients with active lesions. Moreover, TLRs engagement modulated the composition of CD44 variants also in human T helper cells, supporting the hypothesis that pathogens or commensals, through TLRs, in turn modulate the repertoire of CD44 isoforms, thereby controlling the distribution of lesions in the CNS. The interference with this mechanism(s) represents a potential tool for prevention and treatment of autoimmune relapses and exacerbations.

Restricted T-Cell Repertoire in the Epicardial Adipose Tissue of Non-ST Segment Elevation Myocardial Infarction Patients.

Aims: Human epicardial adipose tissue, a dynamic source of multiple bioactive factors, holds a close functional and anatomic relationship with the epicardial coronary arteries and communicates with the coronary artery wall through paracrine and vasocrine secretions. We explored the hypothesis that T-cell recruitment into epicardial adipose tissue (EAT) in patients with non-ST segment elevation myocardial infarction (NSTEMI) could be part of a specific antigen-driven response implicated in acute coronary syndrome onset and progression. Methods and Results: We enrolled 32 NSTEMI patients and 34 chronic coronary syndrome (CCS) patients undergoing coronary artery bypass grafting (CABG) and 12 mitral valve disease (MVD) patients undergoing surgery. We performed EAT proteome profiling on pooled specimens from three NSTEMI and three CCS patients. We performed T-cell receptor (TCR) spectratyping and CDR3 sequencing in EAT and peripheral blood mononuclear cells of 29 NSTEMI, 31 CCS, and 12 MVD patients. We then used computational modeling studies to predict interactions of the TCR beta chain variable region (TRBV) and explore sequence alignments. The EAT proteome profiling displayed a higher content of pro-inflammatory molecules (CD31, CHI3L1, CRP, EMPRINN, ENG, IL-17, IL-33, MMP-9, MPO, NGAL, RBP-4, RETN, VDB) in NSTEMI as compared to CCS (P < 0.0001). CDR3-beta spectratyping showed a TRBV21 enrichment in EAT of NSTEMI (12/29 patients; 41%) as compared with CCS (1/31 patients; 3%) and MVD (none) (ANOVA for trend P < 0.001). Of note, 11/12 (92%) NSTEMI patients with TRBV21 perturbation were at their first manifestation of ACS. Four patients with the first event shared a distinctive TRBV21-CDR3 sequence of 178 bp length and 2/4 were carriers of the human leukocyte antigen (HLA)-A*03:01 allele. A 3D analysis predicted the most likely epitope able to bind HLA-A3*01 and interact with the TRBV21-CDR3 sequence of 178 bp length, while the alignment results were consistent with microbial DNA sequences. Conclusions: Our study revealed a unique immune signature of the epicardial adipose tissue, which led to a 3D modeling of the TCRBV/peptide/HLA-A3 complex, in acute coronary syndrome patients at their first event, paving the way for epitope-driven therapeutic strategies.

A gain-of-glycosylation mutation associated with myoclonus-dystonia syndrome affects trafficking and processing of mouse ε-sarcoglycan in the late secretory pathway.

Missense mutations in the SGCE gene encoding ε-sarcoglycan account for approximately 15% of SGCE-positive cases of myoclonus-dystonia syndrome (MDS) in humans. In this study, we show that while the majority of MDS-associated missense mutants modeled with a murine ε-sarcoglycan cDNA are substrates for endoplasmic reticulum-associated degradation, one mutant, M68T (analogous to human c.275T>C, p.M92T), located in the Ig-like domain of ε-sarcoglycan, results in a gain-of-glycosylation mutation producing a protein that is targeted to the plasma membrane, albeit at reduced levels compared to wild-type ε-sarcoglycan. Removal of the ectopic N-linked glycan failed to restore efficient plasma membrane targeting of M68T demonstrating that the substitution rather than the glycan was responsible for the trafficking defect of this mutant. M68T also colocalized with CD63-positive vesicles in the endosomal-lysosomal system and was found to be more susceptible to lysosomal proteolysis than wild-type ε-sarcoglycan. Finally, we demonstrate impaired ectodomain shedding of M68T, a process that occurs physiologically for ε-sarcoglycan resulting in the lysosomal trafficking of the intracellular C-terminal domain of the protein. Our findings show that functional analysis of rare missense mutations can provide a mechanistic insight into the pathogenesis of MDS and the physiological role of ε-sarcoglycan.

A new beta-chain haemoglobin variant with increased oxygen affinity: Hb Roma [beta115(g17)Ala-->Val].

BACKGROUND: Haemoglobin Roma [beta115(G17)Ala-->Val] is a new adult haemoglobin variant found in a patient presenting a mild hypochromia and microcytosis. We studied this previously uncharacterised variant in order to evaluate the effect on the structural and funcional properties of the Ala-->Val substitution at the alpha1beta1 interface. METHODS AND RESULTS: The variant chain was identified by direct DNA sequencing of the beta-globin gene, which revealed a GCC-->GTC mutation in codon 115. This mutation was confirmed by mass spectrometric analysis of the tetramers and peptides. The oxygen-binding properties of the haemoglobin A/haemoglobin Roma mixture, in which the variant makes up 25% of the haemoglobins, showed a significant increase in oxygen affinity with respect to normal haemoglobin A, both in the absence and presence of 2,3-bisphosphoglycerate. The role of the betaG17 position, situated at the alpha(1)beta(1) interface, has been examined using computational models of haemoglobin Roma and other known betaG17 variants, in comparison with normal haemoglobin A. CONCLUSIONS: This study suggests that the beta115(G17)Ala-->Val substitution at the alpha1beta1 interface is responsible for increased oxygen affinity and mild destabilisation of the haemoglobin Roma. GENERAL SIGNIFICANCE: An amino acid substitution at the G17 position of the alpha1beta1 interface may result in stabilisation of the high affinity R-state of the haemoglobin molecule.

Polymorphism, selection and tandem duplication of transferrin genes in Atlantic cod (Gadus morhua)--conserved synteny between fish monolobal and tetrapod bilobal transferrin loci.

BACKGROUND: The two homologous iron-binding lobes of transferrins are thought to have evolved by gene duplication of an ancestral monolobal form, but any conserved synteny between bilobal and monolobal transferrin loci remains unexplored. The important role played by transferrin in the resistance to invading pathogens makes this polymorphic gene a highly valuable candidate for studying adaptive divergence among local populations. RESULTS: The Atlantic cod genome was shown to harbour two tandem duplicated serum transferrin genes (Tf1, Tf2), a melanotransferrin gene (MTf), and a monolobal transferrin gene (Omp). Whereas Tf1 and Tf2 were differentially expressed in liver and brain, the Omp transcript was restricted to the otoliths. Fish, chicken and mammals showed highly conserved syntenic regions in which monolobal and bilobal transferrins reside, but contrasting with tetrapods, the fish transferrin genes are positioned on three different linkage groups. Sequence alignment of cod Tf1 cDNAs from Northeast (NE) and Northwest (NW) Atlantic populations revealed 22 single nucleotide polymorphisms (SNP) causing the replacement of 16 amino acids, including eight surface residues revealed by the modelled 3D-structures, that might influence the binding of pathogens for removal of iron. SNP analysis of a total of 375 individuals from 14 trans-Atlantic populations showed that the Tf1-NE variant was almost fixed in the Baltic cod and predominated in the other NE Atlantic populations, whereas the NW Atlantic populations were more heterozygous and showed high frequencies of the Tf-NW SNP alleles. CONCLUSIONS: The highly conserved synteny between fish and tetrapod transferrin loci infers that the fusion of tandem duplicated Omp-like genes gave rise to the modern transferrins. The multiple nonsynonymous substitutions in cod Tf1 with putative structural effects, together with highly divergent allele frequencies among different cod populations, strongly suggest evidence for positive selection and local adaptation in trans-Atlantic cod populations.

Targeting SARS-CoV-2 Spike Protein/ACE2 Protein-Protein Interactions: a Computational Study.

The spike glycoprotein (S) of the SARS-CoV-2 virus surface plays a key role in receptor binding and virus entry. The S protein uses the angiotensin converting enzyme (ACE2) for entry into the host cell and binding to ACE2 occurs at the receptor binding domain (RBD) of the S protein. Therefore, the protein-protein interactions (PPIs) between the SARS-CoV-2 RBD and human ACE2, could be attractive therapeutic targets for drug discovery approaches designed to inhibit the entry of SARS-CoV-2 into the host cells. Herein, with the support of machine learning approaches, we report structure-based virtual screening as an effective strategy to discover PPIs inhibitors from ZINC database. The proposed computational protocol led to the identification of a promising scaffold which was selected for subsequent binding mode analysis and that can represent a useful starting point for the development of new treatments of the SARS-CoV-2 infection.