Structural analysis on mutation residues and interfacial water molecules for human TIM disease understanding

10.1186/1471-2105-14-S16-S11BMC Bioinformatics14SUPPL16-BBMI

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

Order through Disorder: Hyper-Mobile C-Terminal Residues Stabilize the Folded State of a Helical Peptide. A Molecular Dynamics Study

Conventional wisdom has it that the presence of disordered regions in the three-dimensional structures of polypeptides not only does not contribute significantly to the thermodynamic stability of their folded state, but, on the contrary, that the presence of disorder leads to a decrease of the corresponding proteins' stability. We have performed extensive 3.4 µs long folding simulations (in explicit solvent and with full electrostatics) of an undecamer peptide of experimentally known helical structure, both with and without its disordered (four residue long) C-terminal tail. Our simulations clearly indicate that the presence of the apparently disordered (in structural terms) C-terminal tail, increases the thermodynamic stability of the peptide's folded (helical) state. These results show that at least for the case of relatively short peptides, the interplay between thermodynamic stability and the apparent structural stability can be rather subtle, with even disordered regions contributing significantly to the stability of the folded state. Our results have clear implications for the understanding of peptide energetics and the design of foldable peptides

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

Mesodynamics in the SARS nucleocapsid measured by NMR field cycling

Protein motions on all timescales faster than molecular tumbling are encoded in the spectral density. The dissection of complex protein dynamics is typically performed using relaxation rates determined at high and ultra-high field. Here we expand this range of the spectral density to low fields through field cycling using the nucleocapsid protein of the SARS coronavirus as a model system. The field-cycling approach enables site-specific measurements of R1 at low fields with the sensitivity and resolution of a high-field magnet. These data, together with high-field relaxation and heteronuclear NOE, provide evidence for correlated rigid-body motions of the entire β-hairpin, and corresponding motions of adjacent loops with a time constant of 0.8 ns (mesodynamics). MD simulations substantiate these findings and provide direct verification of the time scale and collective nature of these motions

Molecular Dynamics Analysis of Apolipoprotein-D - Lipid Hydroperoxide Interactions: Mechanism for Selective Oxidation of Met-93

Background: Recent studies suggest reduction of radical-propagating fatty acid hydroperoxides to inert hydroxides by interaction with apolipoprotein-D (apoD) Met93 may represent an antioxidant function for apoD. The nature and structural consequences of this selective interaction are unknown. Methodology/Principal Findings: Herein we used molecular dynamics (MD) analysis to address these issues. Longtimescale simulations of apoD suggest lipid molecules are bound flexibly, with the molecules free to explore multiple conformations in a binding site at the entrance to the classical lipocalin ligand-binding pocket. Models of 5s- 12s- and 15s hydroperoxyeicosatetraenoic acids were created and the lipids found to wrap around Met93 thus providing a plausible mechanism by which eicosatetraenoic acids bearing hydroperoxides on different carbon atoms can interact with Met93 to yield Met93 sulfoxide (Met93SO). Simulations of glycosylated apoD indicated that a second solvent exposed Met at position 49 was shielded by a triantennerary N-glycan attached to Asn45 thereby precluding lipid interactions. MD simulations of apoD showed B-factors of the loop containing Met93SO were higher in the oxidized protein, indicating increased flexibility that is predicted to destabilize the protein and promote self-association. Conclusions/Significance: These studies provide novel insights into the mechanisms that may contribute to the antioxidant function of apoD and the structural consequences that result if Met93SO is not redox-cycled back to its native state

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