An Update on CFTR Drug Discovery: Opportunities and Challenges

The Biomolecules Special issue on “An Update on CFTR Drug Discovery: Opportunities and Challenges” includes three original research articles and a webinar session focusing on some recent findings concerning CFTR drug discovery [...

A Bittersweet Computational Journey among Glycosaminoglycans

Glycosaminoglycans (GAGs) are linear polysaccharides. In proteoglycans (PGs), they are attached to a core protein. GAGs and PGs can be found as free molecules, associated with the extracellular matrix or expressed on the cell membrane. They play a role in the regulation of a wide array of physiological and pathological processes by binding to different proteins, thus modulating their structure and function, and their concentration and availability in the microenvironment. Unfortunately, the enormous structural diversity of GAGs/PGs has hampered the development of dedicated analytical technologies and experimental models. Similarly, computational approaches (in particular, molecular modeling, docking and dynamics simulations) have not been fully exploited in glycobiology, despite their potential to demystify the complexity of GAGs/PGs at a structural and functional level. Here, we review the state-of-the art of computational approaches to studying GAGs/PGs with the aim of pointing out the “bitter” and “sweet” aspects of this field of research. Furthermore, we attempt to bridge the gap between bioinformatics and glycobiology, which have so far been kept apart by conceptual and technical differences. For this purpose, we provide computational scientists and glycobiologists with the fundamentals of these two fields of research, with the aim of creating opportunities for their combined exploitation, and thereby contributing to a substantial improvement in scientific knowledge

HIV-1 Tat and Heparan Sulfate Proteoglycans Orchestrate the Setup of in Cis and in Trans Cell-Surface Interactions Functional to Lymphocyte Trans-Endothelial Migration

HIV-1 transactivating factor Tat is released by infected cells. Extracellular Tat homodimerizes and engages several receptors, including integrins, vascular endothelial growth factor receptor 2 (VEGFR2) and heparan sulfate proteoglycan (HSPG) syndecan-1 expressed on various cells. By means of experimental cell models recapitulating the processes of lymphocyte trans-endothelial migration, here, we demonstrate that upon association with syndecan-1 expressed on lymphocytes, Tat triggers simultaneously the in cis activation of lymphocytes themselves and the in trans activation of endothelial cells (ECs). This “two-way” activation eventually induces lymphocyte adhesion and spreading onto the substrate and vascular endothelial (VE)-cadherin reorganization at the EC junctions, with consequent endothelial permeabilization, leading to an increased extravasation of Tat-presenting lymphocytes. By means of a panel of biochemical activation assays and specific synthetic inhibitors, we demonstrate that during the above-mentioned processes, syndecan-1, integrins, FAK, src and ERK1/2 engagement and activation are needed in the lymphocytes, while VEGFR2, integrin, src and ERK1/2 are needed in the endothelium. In conclusion, the Tat/syndecan-1 complex plays a central role in orchestrating the setup of the various in cis and in trans multimeric complexes at the EC/lymphocyte interface. Thus, by means of computational molecular modelling, docking and dynamics, we also provide a characterization at an atomic level of the binding modes of the Tat/heparin interaction, with heparin herein used as a structural analogue of the heparan sulfate chains of syndecan-1

Role of Q675H Mutation in Improving SARS-CoV-2 Spike Interaction with the Furin Binding Pocket

Genotype screening was implemented in Italy and showed a significant prevalence of new SARS-CoV-2 mutants carrying Q675H mutation, near the furin cleavage site of spike protein. Currently, this mutation, which is expressed on different SARS-CoV-2 lineages circulating worldwide, has not been thoughtfully investigated. Therefore, we performed phylogenetic and biocomputational analysis to better understand SARS-CoV-2 Q675H mutants’ evolutionary relationships with other circulating lineages and Q675H function in its molecular context. Our studies reveal that Q675H spike mutation is the result of parallel evolution because it arose independently in separate evolutionary clades. In silico data show that the Q675H mutation gives rise to a hydrogen-bonds network in the spike polar region. This results in an optimized directionality of arginine residues involved in interaction of spike with the furin binding pocket, thus improving proteolytic exposure of the viral protein. Furin was predicted to have a greater affinity for Q675H than Q675 substrate conformations. As a consequence, Q675H mutation could confer a fitness advantage to SARS-CoV-2 by promoting a more efficient viral entry. Interestingly, here we have shown that Q675H spike mutation is documented in all the VOCs. This finding highlights that VOCs are still evolving to enhance viral fitness and to adapt to the human host. At the same time, it may suggest Q675H spike mutation involvement in SARS-CoV-2 evolution

Germline <i>NUP98</i> Variants in Two Siblings with a Rothmund–Thomson-Like Spectrum: Protein Functional Changes Predicted by Molecular Modeling

Two adult siblings born to first-cousin parents presented a clinical phenotype reminiscent of Rothmund–Thomson syndrome (RTS), implying fragile hair, absent eyelashes/eyebrows, bilateral cataracts, mottled pigmentation, dental decay, hypogonadism, and osteoporosis. As the clinical suspicion was not supported by the sequencing of RECQL4, the RTS2-causative gene, whole exome sequencing was applied and disclosed the homozygous variants c.83G>A (p.Gly28Asp) and c.2624A>C (p.Glu875Ala) in the nucleoporin 98 (NUP98) gene. Though both variants affect highly conserved amino acids, the c.83G>A looked more intriguing due to its higher pathogenicity score and location of the replaced amino acid between phenylalanine-glycine (FG) repeats within the first NUP98 intrinsically disordered region. Molecular modeling studies of the mutated NUP98 FG domain evidenced a dispersion of the intramolecular cohesion elements and a more elongated conformational state compared to the wild type. This different dynamic behavior may affect the NUP98 functions as the minor plasticity of the mutated FG domain undermines its role as a multi-docking station for RNA and proteins, and the impaired folding can lead to the weakening or the loss of specific interactions. The clinical overlap of NUP98-mutated and RTS2/RTS1 patients, accounted by converging dysregulated gene networks, supports this first-described constitutional NUP98 disorder, expanding the well-known role of NUP98 in cancer

B-cell clonogenic activity of HIV-1 p17 variants is driven by PAR1-mediated EGF transactivation

Combined antiretroviral therapy (cART) for HIV-1 dramatically slows disease progression among HIV individuals. Currently, lymphoma represents the main cause of death among HIV-1-infected patients. Detection of p17 variants (vp17s) endowed with B-cell clonogenic activity in HIV-1-seropositive patients with lymphoma suggests their possible role in lymphomagenesis. Here, we demonstrate that the clonogenic activity of vp17s is mediated by their binding to PAR1 and to PAR1-mediated EGFR transactivation through Gq protein. The entire vp17s-triggered clonogenic process is MMPs dependent. Moreover, phosphoproteomic and bioinformatic analysis highlighted the crucial role of EGFR/PI3K/Akt pathway in modulating several molecules promoting cancer progression, including RAC1, ABL1, p53, CDK1, NPM, Rb, PTP-1B, and STAT1. Finally, we show that a peptide (F1) corresponding to the vp17s functional epitope is sufficient to trigger the PAR1/EGFR/PI3K/Akt pathway and bind PAR1. Our findings suggest novel potential therapeutic targets to counteract vp17-driven lymphomagenesis in HIV patients

Speeding Up the Identification of Cystic Fibrosis Transmembrane Conductance Regulator-Targeted Drugs: An Approach Based on Bioinformatics Strategies and Surface Plasmon Resonance

Cystic fibrosis (CF) is mainly caused by the deletion of Phe 508 (ΔF508) in the cystic fibrosis transmembrane conductance regulator (CFTR) protein that is thus withheld in the endoplasmic reticulum and rapidly degraded by the ubiquitin/proteasome system. New drugs able to rescue ΔF508-CFTR trafficking are eagerly awaited. An integrated bioinformatics and surface plasmon resonance (SPR) approach was here applied to investigate the rescue mechanism(s) of a series of CFTR-ligands including VX809, VX770 and some aminoarylthiazole derivatives (AAT). Computational studies tentatively identified a large binding pocket in the ΔF508-CFTR nucleotide binding domain-1 (NBD1) and predicted all the tested compounds to bind to three sub-regions of this main pocket. Noticeably, the known CFTR chaperone keratin-8 (K8) seems to interact with some residues located in one of these sub-pockets, potentially interfering with the binding of some ligands. SPR results corroborated all these computational findings. Moreover, for all the considered ligands, a statistically significant correlation was determined between their binding capability to ΔF508-NBD1 measured by SPR and the pockets availability measured by computational studies. Taken together, these results demonstrate a strong agreement between the in silico prediction and the SPR-generated binding data, suggesting a path to speed up the identification of new drugs for the treatment of cystic fibrosis

Phenotype Microarray PM with different Nitrogen sources.

<p>Metabolic differences among <i>R</i>. <i>opacus</i> R7 and <i>Rhodococcus</i> sp. BCP1 in presence of different Nitrogen sources (AI, AII, AIII). Based on activity values of phenotype microarray analysis, threshold values were established for every plates. Determined thresholds were high (green), upper middle (light green), lower middle (orange) and low (red) for high, upper middle, lower middle and low activity, respectively.</p

Phenotype Microarray PM with different carboxylic acids-alcohols, amines, amides, esters, fatty acids and polymers-amino acids as carbon sources.

<p>Metabolic differences among <i>R</i>. <i>opacus</i> R7 and <i>Rhodococcus</i> sp. BCP1 in presence of carboxylic acids (BI, BII), of alcohols, amines, amides, esters, fatty acids, polymers (C), and amino acids (D). Based on activity values of phenotype microarray analysis, threshold values were established for every plates. Determined thresholds were high (green), upper middle (light green), lower middle (orange) and low (red) for high, upper middle, lower middle and low activity, respectively.</p