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.

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.

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.

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.

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.

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.

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.

Structure-based virtual screening to identify novel carnitine acetyltransferase activators.

Carnitine acetyltransferase (CAT) is an attractive therapeutic target against fibrosis. We have identified few CAT activators through structure-based virtual screening followed by molecular dynamics simulations for assessment of the binding mode. A set of 10,000 drug-like molecules properly filtered from an initial chemical library of 13 M commercially available compounds were docked into the active site. Virtual hits were selected for in vitro experimental testing to validate the computational findings and the stability of the predicted complexes was evaluated by molecular dynamics simulations. Applied protocol led to the identification of three hit compounds showing promising activity, which can serve as potential scaffolds for further structural optimization. This is the first report of successful discovery of CAT activators through the use of structure-based virtual screening.

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.

Selective Inhibitors of T Cell Receptor Recognition of Antigen-MHC Complexes for Rheumatoid Arthritis.

Autoreactive T cells specific to human collagen type II have a crucial role in the development of rheumatoid arthritis (RA) in the context of MHC class II allele HLA-DRB1-*04. The protein-protein interactions between the T cell receptor (TCR) and the type II collagen bound to the allele MHC of class II may thus represent the target for the development of new drugs against RA. In this study, a structure-based pharmacophore model for potential small molecule inhibitors was developed from protein-protein interface structure. The 3D model obtained was used for a virtual screening workflow, which resulted in three hits for experimental follow up. Three compounds have been identified that interfere with the TCR/collagenII-MHCII (K i values below 10 µM) and open up new possibilities in the treatment of RA.

Identification and antitumor activity of a novel inhibitor of the NIMA-related kinase NEK6.

The NIMA (never in mitosis, gene A)-related kinase-6 (NEK6), which is implicated in cell cycle control and plays significant roles in tumorigenesis, is an attractive target for the development of novel anti-cancer drugs. Here we describe the discovery of a potent ATP site-directed inhibitor of NEK6 identified by virtual screening, adopting both structure- and ligand-based techniques. Using a homology-built model of NEK6 as well as the pharmacophoric features of known NEK6 inhibitors we identified novel binding scaffolds. Twenty-five compounds from the top ranking hits were subjected to in vitro kinase assays. The best compound, i.e. compound 8 ((5Z)-2-hydroxy-4-methyl-6-oxo-5-[(5-phenylfuran-2-yl)methylidene]-5,6-dihydropyridine-3-carbonitrile), was able to inhibit NEK6 with low micromolar IC50 value, also displaying antiproliferative activity against a panel of human cancer cell lines. Our results suggest that the identified inhibitor can be used as lead candidate for the development of novel anti-cancer agents, thus opening the possibility of new therapeutic strategies.

In silico investigation of the molecular effects caused by R123H variant in secretory phospholipase A2-IIA associated with ARDS.

Phospholipase A2-IIA catalyzes the hydrolysis of the sn-2 ester of glycerophospholipids. A rare c.428G > A (p.Arg143His) variant in PLA2G2A gene was found in two infants affected by acute respiratory distress syndrome (ARDS) by whole coding region and exon/intron boundaries sequencing. To obtain insights into the possible molecular effects of the rare R123H mutation in secretory PLA2-IIA (sPLA2-IIA), molecular modelling, molecular dynamics (MD) using principal component analysis (PCA) and continuum electrostatic calculations were conducted on the crystal structure of the wild type protein and on a generated model structure of the R123H mutant. Analysis of MD trajectories indicate that the overall stability of the protein is not affected by this mutation but nevertheless the catalytically crucial H-bond between Tyr51 and Asp91 as well as main electrostatic interactions in the region close to the mutation site are altered. PCA results indicate that the R123H replacement alter the internal molecular motions of the enzyme and that collective motions are increased. Electrostatic surface potential studies suggest that after mutation the interfacial binding to anionic phospholipid membranes and anionic proteins may be changed. The strengthening of electrostatic interactions may be propagated into the active site region thus potentially affecting the substrate recognition and enzymatic activity. Our findings provide the basis for further investigation and advances our understanding of the effects of mutations on sPLA2 structure and function.

Structural model of the full-length Ser/Thr protein kinase StkP from S. pneumoniae and its recognition of peptidoglycan fragments.

The unique eukaryotic-like Ser/Thr protein kinases of Streptococcus pneumoniae, StkP, plays a primary role in the cell division process. It is composed of an intracellular kinase domain, a transmembrane helix and four extracellular PASTA subunits. PASTA domains were shown to interact with cell wall fragments but the key questions related to the molecular mechanism governing ligand recognition remain unclear. To address this issue, the full-length structural model of StkP was generated by combining small-angle X-ray scattering data with the results of computer simulations. Docking and molecular dynamics studies on the generated three-dimensional model structure reveal the possibility of peptidoglycan fragment binding at the hinge regions between PASTA subunits with a preference for a bent hinge between PASTA3 and PASTA4.

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.

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.