Homology modeling and molecular dynamics simulations of MUC1-9/H-2Kb complex suggest novel binding interactions

International audienceHuman MUC1 is over-expressed in human adenocarcinomas and has been used as a target for immunotherapy studies. The 9-mer MUC1-9 peptide has been identified as one of the peptides which binds to murine MHC class I H-2K. The structure of MUC1-9 in complex with H-2K has been modeled and simulated with classical molecular dynamics, based on the x-ray structure of the SEV9 peptide/H-2K complex. Two independent trajectories with the solvated complex (10 ns in length) were produced. Approximately 12 hydrogen bonds were identified during both trajectories to contribute to peptide/MHC complex, as well as 1-2 water mediated hydrogen bonds. Stability of the complex was also confirmed by buried surface area analysis, although the corresponding values were about 20% lower than those of the original x-ray structure. Interestingly, a bulged conformation of the peptide's central region, partially characterized as a -turn, was found exposed form the binding groove. In addition, P1 and P9 residues remained bound in the A and F binding pockets, even though there was a suggestion that P9 was more flexible. The complex lacked numerous water mediated hydrogen bonds that were present in the reference peptide x-ray structure. Moreover, local displacements of residues Asp4, Thr5 and Pro9 resulted in loss of some key interactions with the MHC molecule. This might explain the reduced affinity of the MUC1-9 peptide, relatively to SEV9, for the MHC class I H-2K

Evolution of protein-coupled RNA dynamics during hierarchical assembly of ribosomal complexes

Assembly of 30S ribosomes involves the hierarchical addition of ribosomal proteins that progressively stabilize the folded 16S rRNA. Here, we use three-color single molecule FRET to show how combinations of ribosomal proteins uS4, uS17 and bS20 in the 16S 5' domain enable the recruitment of protein bS16, the next protein to join the complex. Analysis of real-time bS16 binding events shows that bS16 binds both native and non-native forms of the rRNA. The native rRNA conformation is increasingly favored after bS16 binds, explaining how bS16 drives later steps of 30S assembly. Chemical footprinting and molecular dynamics simulations show that each ribosomal protein switches the 16S conformation and dampens fluctuations at the interface between rRNA subdomains where bS16 binds. The results suggest that specific protein-induced changes in the rRNA dynamics underlie the hierarchy of 30S assembly and simplify the search for the native ribosome structure

Mechanism of eIF6 release from the nascent 60S ribosomal subunit.

SBDS protein (deficient in the inherited leukemia-predisposition disorder Shwachman-Diamond syndrome) and the GTPase EFL1 (an EF-G homolog) activate nascent 60S ribosomal subunits for translation by catalyzing eviction of the antiassociation factor eIF6 from nascent 60S ribosomal subunits. However, the mechanism is completely unknown. Here, we present cryo-EM structures of human SBDS and SBDS-EFL1 bound to Dictyostelium discoideum 60S ribosomal subunits with and without endogenous eIF6. SBDS assesses the integrity of the peptidyl (P) site, bridging uL16 (mutated in T-cell acute lymphoblastic leukemia) with uL11 at the P-stalk base and the sarcin-ricin loop. Upon EFL1 binding, SBDS is repositioned around helix 69, thus facilitating a conformational switch in EFL1 that displaces eIF6 by competing for an overlapping binding site on the 60S ribosomal subunit. Our data reveal the conserved mechanism of eIF6 release, which is corrupted in both inherited and sporadic leukemias.Supported by a Federation of European Biochemical Societies Long term Fellowship (to FW), Specialist Programme from Bloodwise [12048] (AJW), the Medical Research Council [MC_U105161083] (AJW) and [U105115237] (RRK), Wellcome Trust strategic award to the Cambridge Institute for Medal Research [100140], Tesni Parry Trust (AJW), Ted’s Gang (AJW) and the Cambridge NIHR Biomedical Research Centre.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/nsmb.311

Substrate binding and translocation of the serotonin transporter studied by docking and molecular dynamics simulations

The serotonin (5-HT) transporter (SERT) plays an important role in the termination of 5-HT-mediated neurotransmission by transporting 5-HT away from the synaptic cleft and into the presynaptic neuron. In addition, SERT is the main target for antidepressant drugs, including the selective serotonin reuptake inhibitors (SSRIs). The three-dimensional (3D) structure of SERT has not yet been determined, and little is known about the molecular mechanisms of substrate binding and transport, though such information is very important for the development of new antidepressant drugs. In this study, a homology model of SERT was constructed based on the 3D structure of a prokaryotic homologous leucine transporter (LeuT) (PDB id: 2A65). Eleven tryptamine derivates (including 5-HT) and the SSRI (S)-citalopram were docked into the putative substrate binding site, and two possible binding modes of the ligands were found. To study the conformational effect that ligand binding may have on SERT, two SERT–5-HT and two SERT–(S)-citalopram complexes, as well as the SERT apo structure, were embedded in POPC lipid bilayers and comparative molecular dynamics (MD) simulations were performed. Our results show that 5-HT in the SERT–5-HTB complex induced larger conformational changes in the cytoplasmic parts of the transmembrane helices of SERT than any of the other ligands. Based on these results, we suggest that the formation and breakage of ionic interactions with amino acids in transmembrane helices 6 and 8 and intracellular loop 1 may be of importance for substrate translocation

Characterizing Structural Transitions Using Localized Free Energy Landscape Analysis

Structural changes in molecules are frequently observed during biological processes like replication, transcription and translation. These structural changes can usually be traced to specific distortions in the backbones of the macromolecules involved. Quantitative energetic characterization of such distortions can greatly advance the atomic-level understanding of the dynamic character of these biological processes.Molecular dynamics simulations combined with a variation of the Weighted Histogram Analysis Method for potential of mean force determination are applied to characterize localized structural changes for the test case of cytosine (underlined) base flipping in a GTCAGCGCATGG DNA duplex. Free energy landscapes for backbone torsion and sugar pucker degrees of freedom in the DNA are used to understand their behavior in response to the base flipping perturbation. By simplifying the base flipping structural change into a two-state model, a free energy difference of upto 14 kcal/mol can be attributed to the flipped state relative to the stacked Watson-Crick base paired state. This two-state classification allows precise evaluation of the effect of base flipping on local backbone degrees of freedom.The calculated free energy landscapes of individual backbone and sugar degrees of freedom expectedly show the greatest change in the vicinity of the flipping base itself, but specific delocalized effects can be discerned upto four nucleotide positions away in both 5' and 3' directions. Free energy landscape analysis thus provides a quantitative method to pinpoint the determinants of structural change on the atomic scale and also delineate the extent of propagation of the perturbation along the molecule. In addition to nucleic acids, this methodology is anticipated to be useful for studying conformational changes in all macromolecules, including carbohydrates, lipids, and proteins

Structural Analysis of a Peptide Fragment of Transmembrane Transporter Protein Bilitranslocase

Using a combination of genomic and post-genomic approaches is rapidly altering the number of identified human influx carriers. A transmembrane protein bilitranslocase (TCDB 2.A.65) has long attracted attention because of its function as an organic anion carrier. It has also been identified as a potential membrane transporter for cellular uptake of several drugs and due to its implication in drug uptake, it is extremely important to advance the knowledge about its structure. However, at present, only the primary structure of bilitranslocase is known. In our work, transmembrane subunits of bilitranslocase were predicted by a previously developed chemometrics model and the stability of these polypeptide chains were studied by molecular dynamics (MD) simulation. Furthermore, sodium dodecyl sulfate (SDS) micelles were used as a model of cell membrane and herein we present a high-resolution 3D structure of an 18 amino acid residues long peptide corresponding to the third transmembrane part of bilitranslocase obtained by use of multidimensional NMR spectroscopy. It has been experimentally confirmed that one of the transmembrane segments of bilitranslocase has alpha helical structure with hydrophilic amino acid residues oriented towards one side, thus capable of forming a channel in the membrane

Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations

Molecular dynamics simulations of membrane proteins have provided deeper insights into their functions and interactions with surrounding environments at the atomic level. However, compared to solvation of globular proteins, building a realistic protein/membrane complex is still challenging and requires considerable experience with simulation software. Membrane Builder in the CHARMM-GUI website (http://www.charmm-gui.org) helps users to build such a complex system using a web browser with a graphical user interface. Through a generalized and automated building process including system size determination as well as generation of lipid bilayer, pore water, bulk water, and ions, a realistic membrane system with virtually any kinds and shapes of membrane proteins can be generated in 5 minutes to 2 hours depending on the system size. Default values that were elaborated and tested extensively are given in each step to provide reasonable options and starting points for both non-expert and expert users. The efficacy of Membrane Builder is illustrated by its applications to 12 transmembrane and 3 interfacial membrane proteins, whose fully equilibrated systems with three different types of lipid molecules (DMPC, DPPC, and POPC) and two types of system shapes (rectangular and hexagonal) are freely available on the CHARMM-GUI website. One of the most significant advantages of using the web environment is that, if a problem is found, users can go back and re-generate the whole system again before quitting the browser. Therefore, Membrane Builder provides the intuitive and easy way to build and simulate the biologically important membrane system

Computational Fragment-Based Binding Site Identification by Ligand Competitive Saturation

Fragment-based drug discovery using NMR and x-ray crystallographic methods has proven utility but also non-trivial time, materials, and labor costs. Current computational fragment-based approaches circumvent these issues but suffer from limited representations of protein flexibility and solvation effects, leading to difficulties with rigorous ranking of fragment affinities. To overcome these limitations we describe an explicit solvent all-atom molecular dynamics methodology (SILCS: Site Identification by Ligand Competitive Saturation) that uses small aliphatic and aromatic molecules plus water molecules to map the affinity pattern of a protein for hydrophobic groups, aromatic groups, hydrogen bond donors, and hydrogen bond acceptors. By simultaneously incorporating ligands representative of all these functionalities, the method is an in silico free energy-based competition assay that generates three-dimensional probability maps of fragment binding (FragMaps) indicating favorable fragment∶protein interactions. Applied to the two-fold symmetric oncoprotein BCL-6, the SILCS method yields two-fold symmetric FragMaps that recapitulate the crystallographic binding modes of the SMRT and BCOR peptides. These FragMaps account both for important sequence and structure differences in the C-terminal halves of the two peptides and also the high mobility of the BCL-6 His116 sidechain in the peptide-binding groove. Such SILCS FragMaps can be used to qualitatively inform the design of small-molecule inhibitors or as scoring grids for high-throughput in silico docking that incorporate both an atomic-level description of solvation and protein flexibility

Proton-Binding Sites of Acid-Sensing Ion Channel 1

Acid-sensing ion channels (ASICs) are proton-gated cation channels that exist throughout the mammalian central and peripheral nervous systems. ASIC1 is the most abundant of all the ASICs and is likely to modulate synaptic transmission. Identifying the proton-binding sites of ASCI1 is required to elucidate its pH-sensing mechanism. By using the crystal structure of ASIC1, the protonation states of each titratable site of ASIC1 were calculated by solving the Poisson-Boltzmann equation under conditions wherein the protonation states of all these sites are simultaneously in equilibrium. Four acidic-acidic residue pairs—Asp238-Asp350, Glu220-Asp408, Glu239-Asp346, and Glu80-Glu417—were found to be highly protonated. In particular, the Glu80-Glu417 pair in the inner pore was completely protonated and possessed 2 H+, implying its possible importance as a proton-binding site. The pKa of Glu239, which forms a pair with a possible pH-sensing site Asp346, differs among each homo-trimer subunit due to the different H-bond pattern of Thr237 in the different protein conformations of the subunits. His74 possessed a pKa of ≈6–7. Conservation of His74 in the proton-sensitive ASIC3 that lacks a residue corresponding to Asp346 may suggest its possible pH-sensing role in proton-sensitive ASICs