Germline NUP98 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.

Virtual Drug Repositioning as a Tool to Identify Natural Small Molecules That Synergize with Lumacaftor in F508del-CFTR Binding and Rescuing.

Cystic fibrosis is a hereditary disease mainly caused by the deletion of the Phe 508 (F508del) of the cystic fibrosis transmembrane conductance regulator (CFTR) protein that is thus withheld in the endoplasmic reticulum and rapidly degraded by the ubiquitin/proteasome system. Cystic fibrosis remains a potentially fatal disease, but it has become treatable as a chronic condition due to some CFTR-rescuing drugs that, when used in combination, increase in their therapeutic effect due to a synergic action. Also, dietary supplementation of natural compounds in combination with approved drugs could represent a promising strategy to further alleviate cystic fibrosis symptoms. On these bases, we screened by in silico drug repositioning 846 small synthetic or natural compounds from the AIFA database to evaluate their capacity to interact with the highly druggable lumacaftor binding site of F508del-CFTR. Among the identified hits, nicotinamide (NAM) was predicted to accommodate into the lumacaftor binding region of F508del-CFTR without competing against the drug but rather stabilizing its binding. The effective capacity of NAM to bind F508del-CFTR in a lumacaftor-uncompetitive manner was then validated experimentally by surface plasmon resonance analysis. Finally, the capacity of NAM to synergize with lumacaftor increasing its CFTR-rescuing activity was demonstrated in cell-based assays. This study suggests the possible identification of natural small molecules devoid of side effects and endowed with the capacity to synergize with drugs currently employed for the treatment of cystic fibrosis, which hopefully will increase the therapeutic efficacy with lower doses.

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 [...].

Cholenic acid derivative UniPR1331 impairs tumor angiogenesis via blockade of VEGF/VEGFR2 in addition to Eph/ephrin.

Angiogenesis, the formation of new blood vessels from preexisting ones, is crucial for tumor growth and metastatization, and is considered a promising therapeutic target. Unfortunately, drugs directed against a specific proangiogenic growth factor or receptor turned out to be of limited benefit for oncology patients, likely due to the high biochemical redundancy of the neovascularization process. In this scenario, multitarget compounds that are able to simultaneously tackle different proangiogenic pathways are eagerly awaited. UniPR1331 is a 3ß-hydroxy-Δ5-cholenic acid derivative, which is already known to inhibit Eph-ephrin interaction. Here, we employed an analysis pipeline consisting of molecular modeling and simulation, surface plasmon resonance spectrometry, biochemical assays, and endothelial cell models to demonstrate that UniPR1331 directly interacts with the vascular endothelial growth factor receptor 2 (VEGFR2) too. The binding of UniPR1331 to VEGFR2 prevents its interaction with the natural ligand vascular endothelial growth factor and subsequent autophosphorylation, signal transduction, and in vitro proangiogenic activation of endothelial cells. In vivo, UniPR1331 inhibits tumor cell-driven angiogenesis in zebrafish. Taken together, these data shed light on the pleiotropic pharmacological effect of UniPR1331, and point to Δ5-cholenic acid as a promising molecular scaffold for the development of multitarget antiangiogenic compounds.

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.

A Boron Delivery Antibody (BDA) with Boronated Specific Residues: New Perspectives in Boron Neutron Capture Therapy from an In Silico Investigation.

Boron Neutron Capture Therapy (BNCT) is a tumor cell-selective radiotherapy based on a nuclear reaction that occurs when the isotope boron-10 (10B) is radiated by low-energy thermal neutrons or epithermal neutrons, triggering a nuclear fission response and enabling a selective administration of irradiation to cells. Hence, we need to create novel delivery agents containing 10B with high tumor selectivity, but also exhibiting low intrinsic toxicity, fast clearance from normal tissue and blood, and no pharmaceutical effects. In the past, boronated monoclonal antibodies have been proposed using large boron-containing molecules or dendrimers, but with no investigations in relation to maintaining antibody specificity and structural and functional features. This work aims at improving the potential of monoclonal antibodies applied to BNCT therapy, identifying in silico the best native residues suitable to be substituted with a boronated one, carefully evaluating the effect of boronation on the 3D structure of the monoclonal antibody and on its binding affinity. A boronated monoclonal antibody was thus generated for specific 10B delivery. In this context, we have developed a case study of Boron Delivery Antibody Identification Pipeline, which has been tested on cetuximab. Cetuximab is an epidermal growth factor receptor (EGFR) inhibitor used in the treatment of metastatic colorectal cancer, metastatic non-small cell lung cancer, and head and neck cancer.

In silico drug repositioning on F508del-CFTR: A proof-of-concept study on the AIFA library.

Computational drug repositioning is of growing interest to academia and industry, for its ability to rapidly screen a huge number of candidates in silico (exploiting comprehensive drug datasets) together with reduced development cost and time. The potential of drug repositioning has not been fully evaluated yet for cystic fibrosis (CF), a disease mainly caused by deletion of Phe 508 (F508del) of the cystic fibrosis transmembrane conductance regulator (CFTR) protein. F508del-CFTR is thus withheld in the endoplasmic reticulum and rapidly degraded by the ubiquitin/proteasome system. CF is still a fatal disease. Nowadays, it is treatable by some CFTR-rescuing drugs, but new-generation drugs with stronger therapeutic benefits and fewer side effects are still awaited. In this manuscript we report about the results of a pilot computational drug repositioning screening in search of F508del-CFTR-targeted drugs performed on AIFA library by means of a dedicated computational pipeline and surface plasmon resonance binding assay to experimentally validate the computational findings.

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.

Recent Strategic Advances in CFTR Drug Discovery: An Overview.

Cystic fibrosis transmembrane conductance regulator (CFTR)-rescuing drugs have already transformed cystic fibrosis (CF) from a fatal disease to a treatable chronic condition. However, new-generation drugs able to bind CFTR with higher specificity/affinity and to exert stronger therapeutic benefits and fewer side effects are still awaited. Computational methods and biosensors have become indispensable tools in the process of drug discovery for many important human pathologies. Instead, they have been used only piecemeal in CF so far, calling for their appropriate integration with well-tried CF biochemical and cell-based models to speed up the discovery of new CFTR-rescuing drugs. This review will give an overview of the available structures and computational models of CFTR and of the biosensors, biochemical and cell-based assays already used in CF-oriented studies. It will also give the reader some insights about how to integrate these tools as to improve the efficiency of the drug discovery process targeted to CFTR.

Heparin and heparan sulfate proteoglycans promote HIV-1 p17 matrix protein oligomerization: computational, biochemical and biological implications.

p17 matrix protein released by HIV+ cells interacts with leukocytes heparan sulfate proteoglycans (HSPGs), CXCR1 and CXCR2 exerting different cytokine-like activities that contribute to AIDS pathogenesis. Since the bioactive form of several cytokines is represented by dimers/oligomers and oligomerization is promoted by binding to heparin or HSPGs, here we evaluated if heparin/HSPGs also promote p17 oligomerization. Heparin favours p17 dimer, trimer and tetramer assembly, in a time- and biphasic dose-dependent way. Heparin-induced p17 oligomerization is of electrostatic nature, being it prevented by NaCl, by removing negative sulfated groups of heparin and by neutralizing positive lysine residues in the p17 N-terminus. A new computational protocol has been implemented to study heparin chains up to 24-mer accommodating a p17 dimer. Molecular dynamics show that, in the presence of heparin, two p17 molecules undergo conformational modifications creating a continuous "electropositive channel" in which heparin sulfated groups interact with p17 basic amino acids, promoting its dimerization. At the cell surface, HSPGs induce p17 oligomerization, as demonstrated by using B-lymphoblastoid Namalwa cells overexpressing the HSPG Syndecan-1. Also, HSPGs on the surface of BJAB and Raji human B-lymphoblastoid cells are required to p17 to induce ERK1/2 activation, suggesting that HS-induced oligomerization plays a role in p17-induced lymphoid dysregulation during AIDS.

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.

A single amino acid substitution confers B-cell clonogenic activity to the HIV-1 matrix protein p17.

Recent data highlight the presence, in HIV-1-seropositive patients with lymphoma, of p17 variants (vp17s) endowed with B-cell clonogenicity, suggesting a role of vp17s in lymphomagenesis. We investigated the mechanisms responsible for the functional disparity on B cells between a wild-type p17 (refp17) and a vp17 named S75X. Here, we show that a single Arginine (R) to Glycine (G) mutation at position 76 in the refp17 backbone (p17R76G), as in the S75X variant, is per se sufficient to confer a B-cell clonogenic potential to the viral protein and modulate, through activation of the PTEN/PI3K/Akt signaling pathway, different molecules involved in apoptosis inhibition (CASP-9, CASP-7, DFF-45, NPM, YWHAZ, Src, PAX2, MAPK8), cell cycle promotion and cancer progression (CDK1, CDK2, CDK8, CHEK1, CHEK2, GSK-3 beta, NPM, PAK1, PP2C-alpha). Moreover, the only R to G mutation at position 76 was found to strongly impact on protein folding and oligomerization by altering the hydrogen bond network. This generates a conformational shift in the p17 R76G mutant which enables a functional epitope(s), masked in refp17, to elicit B-cell growth-promoting signals after its interaction with a still unknown receptor(s). Our findings offer new opportunities to understand the molecular mechanisms accounting for the B-cell growth-promoting activity of vp17s.

New insights into selective PDE4D inhibitors: 3-(Cyclopentyloxy)-4-methoxybenzaldehyde O-(2-(2,6-dimethylmorpholino)-2-oxoethyl) oxime (GEBR-7b) structural development and promising activities to restore memory impairment.

Phosphodiesterase type 4D (PDE4D) has been indicated as a promising target for treating neurodegenerative pathologies such as Alzheimer's Disease (AD). By preventing cAMP hydrolysis, PDE4 inhibitors (PDE4Is) increase the cAMP response element-binding protein (CREB) phosphorylation, synaptic plasticity and long-term memory formation. Pharmacological and behavioral studies on our hit GEBR-7b demonstrated that selective PDE4DIs could improve memory without causing emesis and sedation. The hit development led to new molecule series, herein reported, characterized by a catechol structure bonded to five member heterocycles. Molecular modeling studies highlighted the pivotal role of a polar alkyl chain in conferring selective enzyme interaction. Compound 8a showed PDE4D3 selective inhibition and was able to increase intracellular cAMP levels in neuronal cells, as well as in the hippocampus of freely moving rats. Furthermore, 8a was able to readily cross the blood-brain barrier and enhanced memory performance in mice without causing any emetic-like behavior. These data support the view that PDE4D is an adequate molecular target to restore memory deficits in different neuropathologies, including AD, and also indicate compound 8a as a promising candidate for further preclinical development.

Exome sequencing identifies variants in two genes encoding the LIM-proteins NRAP and FHL1 in an Italian patient with BAG3 myofibrillar myopathy.

Myofibrillar myopathies (MFMs) are genetically heterogeneous dystrophies characterized by the disintegration of Z-disks and myofibrils and are associated with mutations in genes encoding Z-disk or Z-disk-related proteins. The c.626 C > T (p.P209L) mutation in the BAG3 gene has been described as causative of a subtype of MFM. We report a sporadic case of a 26-year-old Italian woman, affected by MFM with axonal neuropathy, cardiomyopathy, rigid spine, who carries the c.626 C > T mutation in the BAG3 gene. The patient and her non-consanguineous healthy parents and brother were studied with whole exome sequencing (WES) to further investigate the genetic basis of this complex phenotype. In the patient, we found that the BAG3 mutation is associated with variants in the NRAP and FHL1 genes that encode muscle-specific, LIM domain containing proteins. Quantitative real time PCR, immunohistochemistry and Western blot analysis of the patient's muscular biopsy showed the absence of NRAP expression and FHL1 accumulation in aggregates in the affected skeletal muscle tissue. Molecular dynamic analysis of the mutated FHL1 domain showed a modification in its surface charge, which could affect its capability to bind its target proteins. To our knowledge this is the first study reporting, in a BAG3 MFM, the simultaneous presence of genetic variants in the BAG3 and FHL1 genes (previously described as independently associated with MFMs) and linking the NRAP gene to MFM for the first time.

Homology modeling, docking studies and molecular dynamic simulations using graphical processing unit architecture to probe the type-11 phosphodiesterase catalytic site: a computational approach for the rational design of selective inhibitors.

Phosphodiesterase 11 (PDE11) is the latest isoform of the PDEs family to be identified, acting on both cyclic adenosine monophosphate and cyclic guanosine monophosphate. The initial reports of PDE11 found evidence for PDE11 expression in skeletal muscle, prostate, testis, and salivary glands; however, the tissue distribution of PDE11 still remains a topic of active study and some controversy. Given the sequence similarity between PDE11 and PDE5, several PDE5 inhibitors have been shown to cross-react with PDE11. Accordingly, many non-selective inhibitors, such as IBMX, zaprinast, sildenafil, and dipyridamole, have been documented to inhibit PDE11. Only recently, a series of dihydrothieno[3,2-d]pyrimidin-4(3H)-one derivatives proved to be selective toward the PDE11 isoform. In the absence of experimental data about PDE11 X-ray structures, we found interesting to gain a better understanding of the enzyme-inhibitor interactions using in silico simulations. In this work, we describe a computational approach based on homology modeling, docking, and molecular dynamics simulation to derive a predictive 3D model of PDE11. Using a Graphical Processing Unit architecture, it is possible to perform long simulations, find stable interactions involved in the complex, and finally to suggest guideline for the identification and synthesis of potent and selective inhibitors.

Modelling the interaction of steroid receptors with endocrine disrupting chemicals.

BACKGROUND: The organic polychlorinated compounds like dichlorodiphenyltrichloroethane with its metabolites and polychlorinated biphenyls are a class of highly persistent environmental contaminants. They have been recognized to have detrimental health effects both on wildlife and humans acting as endocrine disrupters due to their ability of mimicking the action of the steroid hormones, and thus interfering with hormone response. There are several experimental evidences that they bind and activate human steroid receptors. However, despite the growing concern about the toxicological activity of endocrine disrupters, molecular data of the interaction of these compounds with biological targets are still lacking. RESULTS: We have used a flexible docking approach to characterize the molecular interaction of seven endocrine disrupting chemicals with estrogen, progesterone and androgen receptors in the ligand-binding domain. All ligands docked in the buried hydrophobic cavity corresponding to the hormone steroid pocket. The interaction was characterized by multiple hydrophobic contacts involving a different number of residues facing the binding pocket, depending on ligands orientation. The EDC ligands did not display a unique binding mode, probably due to their lipophilicity and flexibility, which conferred them a great adaptability into the hydrophobic and large binding pocket of steroid receptors. CONCLUSION: Our results are in agreement with toxicological data on binding and allow to describe a pattern of interactions for a group of ECD to steroid receptors suggesting the requirement of a hydrophobic cavity to accommodate these chlorine carrying compounds. Although the affinity is lower than for hormones, their action can be brought about by a possible synergistic effect.