Search

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.

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.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) emerge for their capability to better adapt to the human host aimed and enhance human-to-human transmission. Mutations in spike largely contributed to adaptation. Viral persistence is a prerequisite for intra-host virus evolution, and this likely occurred in immunocompromised patients who allow intra-host long-term viral replication. The underlying mechanism leading to the emergence of variants during viral persistence in the immunocompromised host is still unknown. Here, we show the existence of an ensemble of minor mutants in the early biological samples obtained from an immunocompromised patient and their dynamic interplay with the master mutant during a persistent and productive long-term infection. In particular, after 222 days of active viral replication, the original master mutant, named MB610, was replaced by a minor quasispecies (MB61222) expressing two critical mutations in spike, namely Q493K and N501T. Isolation of the two viruses allowed us to show that MB61222 entry into target cells occurred mainly by the fusion at the plasma membrane (PM), whereas endocytosis characterized the entry mechanism used by MB610. Interestingly, coinfection of two human cell lines of different origin with the SARS-CoV-2 isolates highlighted the early and dramatic predominance of MB61222 over MB610 replication. This finding may be explained by a faster replicative activity of MB61222 as compared to MB610 as well as by the capability of MB61222 to induce peculiar viral RNA-sensing mechanisms leading to an increased production of interferons (IFNs) and, in particular, of IFN-induced transmembrane protein 1 (IFITM1) and IFITM2. Indeed, it has been recently shown that IFITM2 is able to restrict SARS-CoV-2 entry occurring by endocytosis. In this regard, MB61222 may escape the antiviral activity of IFITMs by using the PM fusion pathway for entry into the target cell, whereas MB610 cannot escape this host antiviral response during MB61222 coinfection, since it has endocytosis as the main pathway of entry. Altogether, our data support the evidence of quasispecies fighting for host dominance by taking benefit from the cell machinery to restrict the productive infection of competitors in the viral ensemble. This finding may explain, at least in part, the extraordinary rapid worldwide turnover of VOCs that use the PM fusion pathway to enter into target cells over the original pandemic strain.

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.

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.

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.

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.

Significance Immune activation and inflammation are predictors of serious non-AIDS events even in virally suppressed HIV-1−infected individuals. This does not apply to HIV-2−infected patients, who experience a form of attenuated HIV-1 disease. Here, we show that the HIV-1 matrix protein 17 (p17) binds to and activates the common beta chain receptor (βCR). The βCR-activating epitope on p17 is expressed on the matrix protein of HIV-1 ancestors but not on that of HIV-2 and its ancestors. Our finding highlights this epitope as a signature tracing the HIV-1 evolutionary trajectory that may have represented a critical step to enhance the HIV-1 ancestors aggressiveness and early human-to-human transmission. Whether this functional epitope actually marks the pathogenic difference between HIV-1 and HIV-2 needs further investigation., The HIV-1 matrix protein p17 (p17) is a pleiotropic molecule impacting on different cell types. Its interaction with many cellular proteins underlines the importance of the viral protein as a major determinant of human specific adaptation. We previously showed the proangiogenic capability of p17. Here, by integrating functional analysis and receptor binding, we identify a functional epitope that displays molecular mimicry with human erythropoietin (EPO) and promotes angiogenesis through common beta chain receptor (βCR) activation. The functional EPO-like epitope was found to be present in the matrix protein of HIV-1 ancestors SIV originated in chimpanzees (SIVcpz) and gorillas (SIVgor) but not in that of HIV-2 and its ancestor SIVsmm from sooty mangabeys. According to biological data, evolution of the EPO-like epitope showed a clear differentiation between HIV-1/SIVcpz-gor and HIV-2/SIVsmm branches, thus highlighting this epitope on p17 as a divergent signature discriminating HIV-1 and HIV-2 ancestors. P17 is known to enhance HIV-1 replication. Similarly to other βCR ligands, p17 is capable of attracting and activating HIV-1 target cells and promoting a proinflammatory microenvironment. Thus, it is tempting to speculate that acquisition of an epitope on the matrix proteins of HIV-1 ancestors capable of triggering βCR may have represented a critical step to enhance viral aggressiveness and early human-to-human SIVcpz/gor dissemination. The hypothesis that the p17/βCR interaction and βCR abnormal stimulation may also play a role in sustaining chronic activation and inflammation, thus marking the difference between HIV-1 and HIV-2 in term of pathogenicity, needs further investigation.

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.

The identification of multilocus imprinting disturbances (MLID) appears fundamental to uncover molecular pathways underlying imprinting disorders (IDs) and to complete clinical diagnosis of patients. However, MLID genetic associated mechanisms remain largely unknown. To characterize MLID in Beckwith-Wiedemann (BWS) and Silver-Russell (SRS) syndromes, we profiled by MassARRAY the methylation of 12 imprinted differentially methylated regions (iDMRs) in 21 BWS and 7 SRS patients with chromosome 11p15.5 epimutations. MLID was identified in 50% of BWS and 29% of SRS patients as a maternal hypomethylation syndrome. By next-generation sequencing, we searched for putative MLID-causative mutations in genes involved in methylation establishment/maintenance and found two novel missense mutations possibly causative of MLID: one in NLRP2, affecting ADP binding and protein activity, and one in ZFP42, likely leading to loss of DNA binding specificity. Both variants were paternally inherited. In silico protein modelling allowed to define the functional effect of these mutations. We found that MLID is very frequent in BWS/SRS. In addition, since MLID-BWS patients in our cohort show a peculiar pattern of BWS-associated clinical signs, MLID test could be important for a comprehensive clinical assessment. Finally, we highlighted the possible involvement of ZFP42 variants in MLID development and confirmed NLRP2 as causative locus in BWS-MLID.

Background Bacteria belonging to the Rhodococcus genus play an important role in the degradation of many contaminants, including methylbenzenes. These bacteria, widely distributed in the environment, are known to be a powerhouse of numerous degradation functions, due to their ability to metabolize a wide range of organic molecules including aliphatic, aromatic, polycyclic aromatic compounds (PAHs), phenols, and nitriles. In accordance with their immense catabolic diversity, Rhodococcus spp. possess large and complex genomes, which contain a multiplicity of catabolic genes, a high genetic redundancy of biosynthetic pathways and a sophisticated regulatory network. The present study aimed to identify genes involved in the o-xylene degradation in R. opacus strain R7 through a genome-based approach. Results Using genome-based analysis we identified all the sequences in the R7 genome annotated as dioxygenases or monooxygenases/hydroxylases and clustered them into two different trees. The akb, phe and prm sequences were selected as genes encoding respectively for dioxygenases, phenol hydroxylases and monooxygenases and their putative involvement in o-xylene oxidation was evaluated. The involvement of the akb genes in o-xylene oxidation was demonstrated by RT-PCR/qPCR experiments after growth on o-xylene and by the selection of the R7–50 leaky mutant. Although the akb genes are specifically activated for o-xylene degradation, metabolic intermediates of the pathway suggested potential alternative oxidation steps, possibly through monooxygenation. This led us to further investigate the role of the prm and the phe genes. Results showed that these genes were transcribed in a constitutive manner, and that the activity of the Prm monooxygenase was able to transform o-xylene slowly in intermediates as 3,4-dimethylphenol and 2-methylbenzylalcohol. Moreover, the expression level of phe genes, homologous to the phe genes of Rhodococcus spp. 1CP and UPV-1 with a 90% identity, could explain their role in the further oxidation of o-xylene and R7 growth on dimethylphenols. Conclusions These results suggest that R7 strain is able to degrade o-xylene by the Akb dioxygenase system leading to the production of the corresponding dihydrodiol. Likewise, the redundancy of sequences encoding for several monooxygenases/phenol hydroxylases, supports the involvement of other oxygenases converging in the o-xylene degradation pathway in R7 strain. Electronic supplementary material The online version of this article (10.1186/s12864-018-4965-6) contains supplementary material, which is available to authorized users.

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.

Duchenne muscular dystrophy (DMD) is a debilitating fatal X‐linked muscle disorder. Recent findings indicate that IGFs play a central role in skeletal muscle regeneration and development. Among IGFs, insulinlike growth factor 2 (IGF2) is a key regulator of cell growth, survival, migration and differentiation. The type 2 IGF receptor (IGF2R) modulates circulating and tissue levels of IGF2 by targeting it to lysosomes for degradation. We found that IGF2R and the store‐operated Ca2+ channel CD20 share a common hydrophobic binding motif that stabilizes their association. Silencing CD20 decreased myoblast differentiation, whereas blockade of IGF2R increased proliferation and differentiation in myoblasts via the calmodulin/calcineurin/NFAT pathway. Remarkably, anti‐IGF2R induced CD20 phosphorylation, leading to the activation of sarcoplasmic/endoplasmic reticulum Ca2+‐ATPase (SERCA) and removal of intracellular Ca2+. Interestingly, we found that IGF2R expression was increased in dystrophic skeletal muscle of human DMD patients and mdx mice. Blockade of IGF2R by neutralizing antibodies stimulated muscle regeneration, induced force recovery and normalized capillary architecture in dystrophic mdx mice representing an encouraging starting point for the development of new biological therapies for DMD., IGF2R is over‐expressed in Duchenne Muscular Dystrophy (DMD) and mdx muscles. Blockade of IGF2R rescued the dystrophic muscle phenotype, ameliorated vascular architecture defects and improved muscle force.

Cystic fibrosis (CF) is mainly caused by the mutation F508del of the cystic fibrosis transmembrane conductance regulator (CFTR) that is thus retained in the endoplasmic reticulum and degraded. New drugs able to rescue F508del-CFTR trafficking and activity are eagerly awaited, a goal that requires the availability of computational and experimental models closely resembling the F508del-CFTR structure and environment in vivo. Here we describe the development of a biosensor based on F508del-CFTR in a lipid environment that proved to be endowed with a wider analytical potential in respect to the previous CFTR-based biosensors. Integrated with an appropriate computational model of the whole human F508del-CFTR in lipid environment and CFTR stability and functional assays, the new biosensor allowed the identification and characterization at the molecular level of the binding modes of some known F508del-CFTR-rescuing drugs and of a new aminoarylthiazole-Lumacaftor/Tezacaftor hybrid derivative endowed with promising F508del-CFTR-binding and rescuing activity.

Cystic fibrosis (CF) is mainly caused by the deletion of Phe 508 (DeltaF508) 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 DeltaF508-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 DeltaF508-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 DeltaF508-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.

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.

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.

Once released by HIV+ cells, p17 binds heparan sulfate proteoglycans (HSPGs) and CXCR1 on leukocytes causing their dysfunction. By exploiting an approach integrating computational modeling, site-directed mutagenesis of p17, chemical desulfation of heparin, and surface plasmon resonance, we characterized the interaction of p17 with heparin, a HSPG structural analog, and CXCR1. p17 binds to heparin with an affinity (Kd = 190 nm) that is similar to those of other heparin-binding viral proteins. Two stretches of basic amino acids (basic motifs) are present in p17 N and C termini. Neutralization (Arg→Ala substitution) of the N-terminal, but not of the C-terminal basic motif, causes the loss of p17 heparin-binding capacity. The N-terminal heparin-binding motif of p17 partially overlaps the CXCR1-binding domain. Accordingly, its neutralization prevents also p17 binding to the chemochine receptor. Competition experiments demonstrated that free heparin and heparan sulfate (HS), but not selectively 2-O-, 6-O-, and N-O desulfated heparins, prevent p17 binding to substrate-immobilized heparin, indicating that the sulfate groups of the glycosaminoglycan mediate p17 interaction. Evaluation of the p17 antagonist activity of a panel of biotechnological heparins derived by chemical sulfation of the Escherichia coli K5 polysaccharide revealed that the highly N,O-sulfated derivative prevents the binding of p17 to both heparin and CXCR1, thus inhibiting p17-driven chemotactic migration of human monocytes with an efficiency that is higher than those of heparin and HS. Here, we characterized at a molecular level the interaction of p17 with its cellular receptors, laying the basis for the development of heparin-mimicking p17 antagonists.

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.

Rhodococcus opacus strain R7 (CIP107348) degrades several mono- and polycyclic aromatic hydrocarbons. Here, we present the high-quality draft genome sequence of strain R7, consisting of 10,118,052 bp, with a G+C content of 67.0%, 9,602 protein-coding genes, and 62 RNAs genes.

In the present study, the response of Rhodococcus aetherivorans BCP1 and Rhodococcus opacus R7 to various stress conditions and several antimicrobials was examined by PM in relation with genetic determinants, as revealed by annotation analysis of the two genomes. Comparison between metabolic activities and genetic features of BCP1 and R7 provided new insight into the environmental persistence of these two members of the genus Rhodococcus. (C) 2016 Institut Pasteur. Published by Elsevier Masson SAS. All rights reserved.

Alzheimer′s disease has recently emerged as a possible field of application for PDE4D inhibitors (PDE4DIs). The great structure similarity among the various PDE4 isoforms and, furthermore, the lack of the full length crystal structure of the enzyme, impaired the rational design of new selective PDE4DIs. In this paper, with the aim of exploring new insights into the PDE4D binding, we tackled the problem by performing a computational study based on docking simulations combined with molecular dynamics (D‐MD). Our work uniquely identified the binding mode and the key residues involved in the interaction with a number of in‐house catechol iminoether derivatives, acting as PDE4DIs. Moreover, the new binding mode was tested using a series of analogues previously reported by us and it was used to confirm their key structural features to allow PDE4D inhibition. The binding model disclosed within the current computational study may prove to be useful to further advance the design and synthesis of novel, more potent and selective, PDE4D inhibitors.

In this paper comparative genome and phenotype microarray analyses of Rhodococcus sp. BCP1 and Rhodococcus opacus R7 were performed. Rhodococcus sp. BCP1 was selected for its ability to grow on short-chain n-alkanes and R. opacus R7 was isolated for its ability to grow on naphthalene and on o-xylene. Results of genome comparison, includ- ing BCP1, R7, along with other Rhodococcus reference strains, showed that at least 30% of the genome of each strain presented unique sequences and only 50% of the predicted proteome was shared. To associate genomic features with metabolic capabilities of BCP1 and R7 strains, hundreds of different growth conditions were tested through Phenotype Microarray, by using Biolog plates and plates manually prepared with additional xenobiotic compounds. Around one-third of the surveyed carbon sources was utilized by both strains although R7 generally showed higher metabolic activity values compared to BCP1. More- over, R7 showed broader range of nitrogen and sulphur sources. Phenotype Microarray data were combined with genomic analysis to genetically support the metabolic features of the two strains. The genome analysis allowed to identify some gene clusters involved in the metabolism of the main tested xenobiotic compounds. Results show that R7 contains multi- ple genes for the degradation of a large set of aromatic and PAHs compounds, while a lower variability in terms of genes predicted to be involved in aromatic degradation was found in BCP1. This genetic feature can be related to the strong genetic pressure exerted by the two different environment from which the two strains were isolated. According to this, in the BCP1 genome the smo gene cluster involved in the short-chain n-alkanes degradation, is included in one of the unique regions and it is not conserved in the Rhodo- coccus strains compared in this work. Data obtained underline the great potential of these two Rhodococcus spp. strains for biodegradation and environmental decontamination processes.

Activated Protein C (APC) is a multifunctional serine protease, primarily known for its anticoagulant function in the coagulation system. Several studies have already elucidated its role in counteracting apoptosis and inflammation in cells, while significant effort is still ongoing for defining its involvement in sepsis. Earlier literature has shown that the antiseptic function of APC is mediated by its binding to leukocyte integrins, which is due to the presence of the integrin binding motif Arg-Gly-Asp at the N-terminus of the APC catalytic chain. Many natural mutants have been identified in patients with Protein C deficiency diagnosis including a variant of specificity pocket (Gly216Asp). In this work, we present a molecular model of the complex of APC with αVβ3 integrin obtained by protein-protein docking approach. A computational analysis of this variant is hereby presented, based on molecular dynamics and docking simulations, aiming at investigating the effects of the Gly216Asp mutation on the protein conformation and inferring its functional implications. Our study shows that such mutation is likely to impair the protease activity while preserving the overall protein fold. Moreover, superposition of the integrin binding motifs in wild-type and mutant forms suggests that the interaction with integrin can still occur and thus the mutant is likely to retain its antiseptic function related to the neutrophyl integrin binding. Therapeutic applications could result in this APC mutant which retains antiseptic function without anticoagulant side effects.

Rhodococcus sp. strain BCP1 (DSM 44980) co-metabolizes chlorinated compounds and mineralizes a broad range of alkanes being highly tolerant to these toxic chemicals. Here, we present the high-quality draft genome sequence of strain BCP1 consisting of 6,231,823 bp, with a G+C content of 70.4%, 5,902 protein-coding genes, and 58 RNAs genes. Rhodococcus genus comprises Gram-positive, non-sporulating, aerobic bacteria that are widely distributed in the environment (1). Rhodococcus sp. strain BCP1 (formerly: Rhodococcus aetherovorans strain BCP1, DSM 44980) was selected from an aerobic butane-utilizing consortium as the prevailing isolate able to co-metabolize chloroform, vinyl chloride and trichloroethylene (2, 3). As BCP1 also catabolizes a wide range of aliphatic, alicyclic and carboxylated alkanes, it represents a strain of considerable environmental and industrial interest (4). The genome sequencing of Rhodococcus sp. BCP1 was performed using 454 sequencing technology (Roche GS FLX Titanium). The total number of sequence reads were 668,686 from one shotgun library and 353,744 from one paired-end library (8-kb inserts). All the reads were assembled using Newbler 2.6 into 123 contigs with an N50 length of 237,787 bp and an average genome coverage of 65x. Based on paired-end directional information, the contigs were further ordered into 3 scaffolds, giving a total genome size of 6.2 Mb with a G+C content of 70.4%.

Once released by HIV(+) cells, p17 binds heparan sulfate proteoglycans (HSPGs) and CXCR1 on leukocytes causing their dysfunction. By exploiting an approach integrating computational modeling, site-directed mutagenesis of p17, chemical desulfation of heparin, and surface plasmon resonance, we characterized the interaction of p17 with heparin, a HSPG structural analog, and CXCR1. p17 binds to heparin with an affinity (K(d) = 190 nm) that is similar to those of other heparin-binding viral proteins. Two stretches of basic amino acids (basic motifs) are present in p17 N and C termini. Neutralization (Arg→Ala substitution) of the N-terminal, but not of the C-terminal basic motif, causes the loss of p17 heparin-binding capacity. The N-terminal heparin-binding motif of p17 partially overlaps the CXCR1-binding domain. Accordingly, its neutralization prevents also p17 binding to the chemochine receptor. Competition experiments demonstrated that free heparin and heparan sulfate (HS), but not selectively 2-O-, 6-O-, and N-O desulfated heparins, prevent p17 binding to substrate-immobilized heparin, indicating that the sulfate groups of the glycosaminoglycan mediate p17 interaction. Evaluation of the p17 antagonist activity of a panel of biotechnological heparins derived by chemical sulfation of the Escherichia coli K5 polysaccharide revealed that the highly N,O-sulfated derivative prevents the binding of p17 to both heparin and CXCR1, thus inhibiting p17-driven chemotactic migration of human monocytes with an efficiency that is higher than those of heparin and HS. Here, we characterized at a molecular level the interaction of p17 with its cellular receptors, laying the basis for the development of heparin-mimicking p17 antagonists.

CB2 receptor belongs to the large family of G-protein coupled receptors (GPCRs) controlling a wide variety of signal transduction. The recent crystallographic determination of human β2 adrenoreceptor and its high sequence similarity with human CB2 receptor (hCB2) prompted us to compute a theoretical model of hCB2 based also on β2 adrenoreceptor coordinates. This model has been employed to perform docking and molecular dynamic simulations on WIN-55,212-2 (CB2 agonist commonly used in binding experiments), in order to identify the putative CB2 receptor agonist binding site, followed by molecular docking studies on a series of indol-3-yl-tetramethylcyclopropyl ketone derivatives, a novel class of potent CB2 agonists. Successively, docking-based Comparative Molecular Fields Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) studies were also performed. The CoMSIA model resulted to be the more predictive, showing r(ncv)(2) = 0.96, r(cv)(2) = 0.713, SEE = 0.193, F = 125.223, and r(2)(pred) = 0.78. The obtained 3D-QSAR models allowed us to derive more complete guidelines for the design of new analogues with improved potency so as to synthesize new indoles showing high CB2 affinity.

Background The design of mutants in protein functional regions, such as the ligand binding sites, is a powerful approach to recognize the determinants of specific protein activities in cellular pathways. For an exhaustive analysis of selected positions of protein structure large scale mutagenesis techniques are often employed, with laborious and time consuming experimental set-up. 'In silico' mutagenesis and screening simulation represents a valid alternative to laboratory methods to drive the 'in vivo' testing toward more focused objectives. Results We present here a high performance computational procedure for large-scale mutant modelling and subsequent evaluation of the effect on ligand binding affinity. The mutagenesis was performed with a 'saturation' approach, where all 20 natural amino acids were tested in positions involved in ligand binding sites. Each modelled mutant was subjected to molecular docking simulation and stability evaluation. The simulated protein-ligand complexes were screened for their impairment of binding ability based on change of calculated Ki compared to the wild-type. An example of application to the Endothelial Protein C Receptor residues involved in lipid binding is reported. Conclusion The computational pipeline presented in this work is a useful tool for the design of structurally stable mutants with altered affinity for ligand binding, considerably reducing the number of mutants to be experimentally tested. The saturation mutagenesis procedure does not require previous knowledge of functional role of the residues involved and allows extensive exploration of all possible substitutions and their pairwise combinations. Mutants are screened by docking simulation and stability evaluation followed by a rationally driven selection of those presenting the required characteristics. The method can be employed in molecular recognition studies and as a preliminary approach to select models for experimental testing.

The H5N1 virus neuraminidase structure was solved in two different conformations depending on the inhibitor concentration. In the absence of oseltamivir or at a low concentration, the neuraminidase structure assumes an open form that closes at a high oseltamivir concentration due to the shift of the so-called 150-loop near the active site. Although the close conformation is similar to all the other structurally known neuraminidase types, it doesn't appear to be the most likely physiological condition for N1. To investigate the specific ligand binding properties of the open form, we screened by docking simulation, a large dataset of ligands and compared the results with closed form. The virtual screening procedure was implemented in a docking pipeline that also performs a step-by-step, target specific, filtering approach for data reduction. The selected ligands display binding ability involving multiple sites of interaction including the active site and an adjacent cavity made available by the 150-loop shift. Two ligands are especially interesting and are proposed as substituents to design oseltamivir derivatives specifically suited for the open conformation.

Protein C (PC) is a key regulator of blood clotting and inflammation. Its inherited deficiency is associated with venous thromboembolism, and recombinant activated PC is currently used to increase survival in severe sepsis. The molecular basis of inherited PC deficiency is heterogeneous. Due to its multiple physiologic interactions and functions, and its modular structure, natural variants aid in the understanding of the relationship between critical residues and discrete functions. This knowledge has important therapeutic implications in the planning of a recombinant activated PC with a specific therapeutic target and devoid of major collateral effects. A way of predicting important functional consequences of residue variation is the use of molecular modeling and structural interpretation of amino acidic substitutions. A study of 21 out of 32 identified PC gene (PROC) variants is presented. For three of them, localized in the active site, electrostatic potential variation was calculated. For more than half of the studied variants, an explanation for the functional impairment could be derived from computational analysis, allowing a focused choice of which variants it is worthwhile pursuing. Hum Mutat 28(4), 345–355, 2007. © 2006 Wiley-Liss, Inc.