Fusing simulation and experiment: The effect of mutations on the structure and activity of the influenza fusion peptide

During the infection process, the influenza fusion peptide (FP) inserts into the host membrane, playing a crucial role in the fusion process between the viral and host membranes. In this work we used a combination of simulation and experimental techniques to analyse the molecular details of this process, which are largely unknown. Although the FP structure has been obtained by NMR in detergent micelles, there is no atomic structure information in membranes. To answer this question, we performed bias-exchange metadynamics (BE-META) simulations, which showed that the lowest energy states of the membrane-inserted FP correspond to helical-hairpin conformations similar to that observed in micelles. BE-META simulations of the G1V, W14A, G12A/G13A and G4A/G8A/G16A/G20A mutants revealed that all the mutations affect the peptide's free energy landscape. A FRET-based analysis showed that all the mutants had a reduced fusogenic activity relative to the WT, in particular the mutants G12A/G13A and G4A/G8A/G16A/G20A. According to our results, one of the major causes of the lower activity of these mutants is their lower membrane affinity, which results in a lower concentration of peptide in the bilayer. These findings contribute to a better understanding of the influenza fusion process and open new routes for future studies

Probenecid Inhibits the Human Bitter Taste Receptor TAS2R16 and Suppresses Bitter Perception of Salicin

Bitter taste stimuli are detected by a diverse family of G protein-coupled receptors (GPCRs) expressed in gustatory cells. Each bitter taste receptor (TAS2R) responds to an array of compounds, many of which are toxic and can be found in nature. For example, human TAS2R16 (hTAS2R16) responds to β-glucosides such as salicin, and hTAS2R38 responds to thiourea-containing molecules such as glucosinolates and phenylthiocarbamide (PTC). While many substances are known to activate TAS2Rs, only one inhibitor that specifically blocks bitter receptor activation has been described. Here, we describe a new inhibitor of bitter taste receptors, p-(dipropylsulfamoyl)benzoic acid (probenecid), that acts on a subset of TAS2Rs and inhibits through a novel, allosteric mechanism of action. Probenecid is an FDA-approved inhibitor of the Multidrug Resistance Protein 1 (MRP1) transporter and is clinically used to treat gout in humans. Probenecid is also commonly used to enhance cellular signals in GPCR calcium mobilization assays. We show that probenecid specifically inhibits the cellular response mediated by the bitter taste receptor hTAS2R16 and provide molecular and pharmacological evidence for direct interaction with this GPCR using a non-competitive (allosteric) mechanism. Through a comprehensive analysis of hTAS2R16 point mutants, we define amino acid residues involved in the probenecid interaction that result in decreased sensitivity to probenecid while maintaining normal responses to salicin. Probenecid inhibits hTAS2R16, hTAS2R38, and hTAS2R43, but does not inhibit the bitter receptor hTAS2R31 or non-TAS2R GPCRs. Additionally, structurally unrelated MRP1 inhibitors, such as indomethacin, fail to inhibit hTAS2R16 function. Finally, we demonstrate that the inhibitory activity of probenecid in cellular experiments translates to inhibition of bitter taste perception of salicin in humans. This work identifies probenecid as a pharmacological tool for understanding the cell biology of bitter taste and as a lead for the development of broad specificity bitter blockers to improve nutrition and medical compliance

The Impact of Small Molecule Binding on the Energy Landscape of the Intrinsically Disordered Protein C-Myc

Intrinsically disordered proteins are attractive therapeutic targets owing to their prevalence in several diseases. Yet their lack of well-defined structure renders ligand discovery a challenging task. An intriguing example is provided by the oncoprotein c-Myc, a transcription factor that is over expressed in a broad range of cancers. Transcriptional activity of c-Myc is dependent on heterodimerization with partner protein Max. This protein-protein interaction is disrupted by the small molecule 10058-F4 (1), that binds to monomeric and disordered c-Myc. To rationalize the mechanism of inhibition, structural ensembles for the segment of the c-Myc domain that binds to 1 were computed in the absence and presence of the ligand using classical force fields and explicit solvent metadynamics molecular simulations. The accuracy of the computed structural ensembles was assessed by comparison of predicted and measured NMR chemical shifts. The small molecule 1 was found to perturb the composition of the apo equilibrium ensemble and to bind weakly to multiple distinct c-Myc conformations. Comparison of the apo and holo equilibrium ensembles reveals that the c-Myc conformations binding 1 are already partially formed in the apo ensemble, suggesting that 1 binds to c-Myc through an extended conformational selection mechanism. The present results have important implications for rational ligand design efforts targeting intrinsically disordered proteins

Association Between TAS2R38 Gene Polymorphisms and Colorectal Cancer Risk: A Case-Control Study in Two Independent Populations of Caucasian Origin

Molecular sensing in the lingual mucosa and in the gastro-intestinal tract play a role in the detection of ingested harmful drugs and toxins. Therefore, genetic polymorphisms affecting the capability of initiating these responses may be critical for the subsequent efficiency of avoiding and/or eliminating possible threats to the organism. By using a tagging approach in the region of Taste Receptor 2R38 (TAS2R38) gene, we investigated all the common genetic variation of this gene region in relation to colorectal cancer risk with a case-control study in a German population (709 controls and 602 cases) and in a Czech population (623 controls and 601 cases). We found that there were no significant associations between individual SNPs of the TAS2R38 gene and colorectal cancer in the Czech or in the German population, nor in the joint analysis. However, when we analyzed the diplotypes and the phenotypes we found that the non-taster group had an increased risk of colorectal cancer in comparison to the taster group. This association was borderline significant in the Czech population, (OR = 1.28, 95% CI 0.99–1.67; Pvalue = 0.058) and statistically significant in the German population (OR = 1.36, 95% CI 1.06–1.75; Pvalue = 0.016) and in the joint analysis (OR = 1.34, 95% CI 1.12–1.61; Pvalue = 0.001). In conclusion, we found a suggestive association between the human bitter tasting phenotype and the risk of CRC in two different populations of Caucasian origin

Visualizing the Reaction Coordinate of an O-GlcNAc Hydrolase

N-Acetylglucosamine β-O-linked to serine and threonine residues of nucleocytoplasmic proteins (O-GlcNAc) has been linked to neurodegeneration, cellular stress response, and transcriptional regulation. Removal of O-GlcNAc is catalyzed by O-GlcNAcase (OGA) using a substrate-assisted catalytic mechanism. Here we define the reaction coordinate using chemical approaches and directly observe both a Michaelis complex and the oxazoline intermediate

Bioinformatics and molecular modeling in glycobiology

The field of glycobiology is concerned with the study of the structure, properties, and biological functions of the family of biomolecules called carbohydrates. Bioinformatics for glycobiology is a particularly challenging field, because carbohydrates exhibit a high structural diversity and their chains are often branched. Significant improvements in experimental analytical methods over recent years have led to a tremendous increase in the amount of carbohydrate structure data generated. Consequently, the availability of databases and tools to store, retrieve and analyze these data in an efficient way is of fundamental importance to progress in glycobiology. In this review, the various graphical representations and sequence formats of carbohydrates are introduced, and an overview of newly developed databases, the latest developments in sequence alignment and data mining, and tools to support experimental glycan analysis are presented. Finally, the field of structural glycoinformatics and molecular modeling of carbohydrates, glycoproteins, and protein–carbohydrate interaction are reviewed

METAGUI. A VMD interface for analyzing metadynamics and molecular dynamics simulations

We present a new computational tool, METAGUI, which extends the VMD program with a graphical user interface that allows constructing a thermodynamic and kinetic model of a given process simulated by large-scale molecular dynamics. The tool is specially designed for analyzing metadynamics based simulations. The huge amount of diverse structures generated during such a simulation is partitioned into a set of microstates (i.e. structures with similar values of the collective variables). Their relative free energies are then computed by a weighted-histogram procedure and the most relevant free energy wells are identified by diagonalization of the rate matrix followed by a commitor analysis. All this procedure leads to a convenient representation of the metastable states and long-time kinetics of the system which can be compared with experimental data. The tool allows to seamlessly switch between a collective variables space representation of microstates and their atomic structure representation, which greatly facilitates the set-up and analysis of molecular dynamics simulations. METAGUI is based on the output format of the PLUMED plugin, making it compatible with a number of different molecular dynamics packages like AMBER, NAMD, GROMACS and several others. The METAGUI source files can be downloaded from the PLUMED web site (http://www.plumed-code.org). RI Pietrucci, Fabio/C-2326-200

The conformational free energy landscape of beta-D-glucopyranose. implications for substrate preactivation in beta-glucoside hydrolases

Using ab initio metadynamics we have computed the conformational free energy landscape of beta-D-glucopyranose as a function of the puckering coordinates. We show that the correspondence between the free energy and the Stoddard's pseudorotational itinerary for the system is rather poor. The number of free energy minima (9) is smaller than the number of ideal structures (13). Moreover, only six minima correspond to a canonical conformation. The structural features, the electronic properties, and the relative stability of the predicted conformers permit the rationalization of the occurrence of distorted sugar conformations in all the available X-ray structures of beta-glucoside hydrolase Michaelis complexes. We show that these enzymes recognize the most stable distorted conformers of the isolated substrate and at the same time the ones better prepared for catalysis in terms of bond elongation/shrinking and charge distribution. This suggests that the factors governing the distortions present in these complexes are largely dictated by the intrinsic properties of a single glucose unit