Molecular Modeling of Signal Peptide Recognition by Eukaryotic Sec Complexes

Here, we review recent molecular modelling and simulation studies of the Sec translocon, the primary component/channel of protein translocation into the endoplasmic reticulum (ER) and bacterial periplasm, respectively. Our focus is placed on the eukaryotic Sec61, but we also mention modelling studies on prokaryotic SecY since both systems operate in related ways. Cryo-EM structures are now available for different conformational states of the Sec61 complex, ranging from the idle or closed state over an inhibited state with the inhibitor mycolactone bound near the lateral gate, up to a translocating state with bound substrate peptide in the translocation pore. For all these states, computational studies have addressed the conformational dynamics of the translocon with respect to the pore ring, the plug region, and the lateral gate. Also, molecular simulations are addressing mechanistic issues of insertion into the ER membrane vs. translocation into the ER, how signal-peptides are recognised at all in the translocation pore, and how accessory proteins affect the Sec61 conformation in the co- and post-translational pathways

How does Sec63 affect the conformation of Sec61 in yeast?

The Sec complex catalyzes the translocation of proteins of the secretory pathway into the endoplasmic reticulum and the integration of membrane proteins into the endoplasmic reticulum membrane. Some substrate peptides require the presence and involvement of accessory proteins such as Sec63. Recently, a structure of the Sec complex from Saccharomyces cerevisiae, consisting of the Sec61 channel and the Sec62, Sec63, Sec71 and Sec72 proteins was determined by cryo-electron microscopy (cryo-EM). Here, we show by co-precipitation that the Sec61 channel subunit Sbh1 is not required for formation of stable Sec63-Sec61 contacts. Molecular dynamics simulations started from the cryo-EM conformation of Sec61 bound to Sec63 and of unbound Sec61 revealed how Sec63 affects the conformation of Sec61 lateral gate, plug, pore region and pore ring diameter via three intermolecular contact regions. Molecular docking of SRP-dependent vs. SRP-independent signal peptide chains into the Sec61 channel showed that the pore regions affected by presence/absence of Sec63 play a crucial role in positioning the signal anchors of SRP-depen dent substrates nearby the lateral gate

TNF Receptor Type II as an Emerging Drug Target for the Treatment of Cancer, Autoimmune Diseases, and Graft-Versus-Host Disease: Current Perspectives and In Silico Search for Small Molecule Binders

There is now compelling evidence that TNF receptor type II (TNFR2) is predominantly expressed on CD4+Foxp3+ regulatory T cells (Tregs) and myeloid-derived suppressor cells (MDSCs), and plays a major role in the expansion and function of Tregs and MDSCs. Consequently, targeting of TNFR2 by either antagonists or agonists may represent a novel strategy in the treatment of cancer and autoimmune diseases, by downregulating or upregulating suppressor cell activity. The advance in the understanding of complex structure of TNFR2 and its binding with TNF at molecular levels offers opportunity for structure-guided drug discovery. This article reviews the current evidences regarding the decisive role of TNFR2 in immunosuppressive function of Tregs and MDSCs, and the current effort to develop novel TNFR2-targeting therapeutic agents in the treatment of cancer, autoimmune diseases, and graft-versus-host disease. To shed light on the potential TNFR2-targeting small molecules, we for the first time performed virtual screening of 400,000 natural compounds against the two TNF-binding sites, regions 3 and 4, of TNFR2. Our result showed that the top hits at region 4 had slightly higher docking energies than those at region 3. Nevertheless, free energy calculation from the TNF–TNFR2 molecular dynamics simulation revealed that the binding strength of TNF in region 3 is only one-tenth of that in region 4. This suggests that region 3 is a potentially more viable binding site to be targeted by small molecules than region 4. Therefore, the effectiveness in targeting region 3 of TNFR2 deserves further investigation

Quantitative Proteomics and Differential Protein Abundance Analysis after Depletion of Putative mRNA Receptors in the ER Membrane of Human Cells Identifies Novel Aspects of mRNA Targeting to the ER

In human cells, one-third of all polypeptides enter the secretory pathway at the endoplasmic reticulum (ER). The specificity and efficiency of this process are guaranteed by targeting of mRNAs and/or polypeptides to the ER membrane. Cytosolic SRP and its receptor in the ER membrane facilitate the cotranslational targeting of most ribosome-nascent precursor polypeptide chain (RNC) complexes together with the respective mRNAs to the Sec61 complex in the ER membrane. Alternatively, fully synthesized precursor polypeptides are targeted to the ER membrane post-translationally by either the TRC, SND, or PEX19/3 pathway. Furthermore, there is targeting of mRNAs to the ER membrane, which does not involve SRP but involves mRNA- or RNC-binding proteins on the ER surface, such as RRBP1 or KTN1. Traditionally, the targeting reactions were studied in cell-free or cellular assays, which focus on a single precursor polypeptide and allow the conclusion of whether a certain precursor can use a certain pathway. Recently, cellular approaches such as proximity-based ribosome profiling or quantitative proteomics were employed to address the question of which precursors use certain pathways under physiological conditions. Here, we combined siRNA-mediated depletion of putative mRNA receptors in HeLa cells with label-free quantitative proteomics and differential protein abundance analysis to characterize RRBP1- or KTN1- involving precursors and to identify possible genetic interactions between the various targeting pathways. Furthermore, we discuss the possible implications on the so-called TIGER domains and critically discuss the pros and cons of this experimental approach

Synthesis, Biological Evaluation and Docking Studies of Ring-Opened Analogues of Ipomoeassin F

The plant-derived macrocyclic resin glycoside ipomoeassin F (Ipom-F) binds to Sec61α and significantly disrupts multiple aspects of Sec61-mediated protein biogenesis at the endoplas mic reticulum, ultimately leading to cell death. However, extensive assessment of Ipom-F as a molecular tool and a therapeutic lead is hampered by its limited production scale, largely caused by intramolecular assembly of the macrocyclic ring. Here, using in vitro and/or in cellula biological assays to explore the first series of ring-opened analogues for the ipomoeassins, and indeed all resin glycosides, we provide clear evidence that macrocyclic integrity is not required for the cytotoxic inhibition of Sec61-dependent protein translocation by Ipom-F. Furthermore, our modeling suggests that open-chain analogues of Ipom-F can interact with multiple sites on the Sec61α subunit, most likely located at a previously identified binding site for mycolactone and/or the so-called lateral gate. Subsequent in silico-aided design led to the discovery of the stereochemically simplified analogue 3 as a potent, alternative lead compound that could be synthesized much more efficiently than Ipom-F and will accelerate future ipomoeassin research in chemical biology and drug discovery. Our work may also inspire further exploration of ring-opened analogues of other resin glycosides

Synthesis, Biological Evaluation and Docking Studies of Ring-Opened Analogues of Ipomoeassin F

The plant-derived macrocyclic resin glycoside ipomoeassin F (Ipom-F) binds to Sec61α and significantly disrupts multiple aspects of Sec61-mediated protein biogenesis at the endoplasmic reticulum, ultimately leading to cell death. However, extensive assessment of Ipom-F as a molecular tool and a therapeutic lead is hampered by its limited production scale, largely caused by intramolecular assembly of the macrocyclic ring. Here, using in vitro and/or in cellula biological assays to explore the first series of ring-opened analogues for the ipomoeassins, and indeed all resin glycosides, we provide clear evidence that macrocyclic integrity is not required for the cytotoxic inhibition of Sec61-dependent protein translocation by Ipom-F. Furthermore, our modeling suggests that open-chain analogues of Ipom-F can interact with multiple sites on the Sec61α subunit, most likely located at a previously identified binding site for mycolactone and/or the so-called lateral gate. Subsequent in silico-aided design led to the discovery of the stereochemically simplified analogue 3 as a potent, alternative lead compound that could be synthesized much more efficiently than Ipom-F and will accelerate future ipomoeassin research in chemical biology and drug discovery. Our work may also inspire further exploration of ring-opened analogues of other resin glycosides

De novo inference of protein function from coarse-grained dynamics

Inference of molecular function of proteins is the fundamental task in the quest for understanding cellular processes. The task is getting increasingly difficult with thousands of new proteins discovered each day. The difficulty arises primarily due to lack of high-throughput experimental technique for assessing protein molecular function, a lacunae that computational approaches are trying hard to fill. The latter too faces a major bottleneck in absence of clear evidence based on evolutionary information. Here we propose a de novo approach to annotate protein molecular function through structural dynamics match for a pair of segments from two dissimilar proteins, which may share even <10% sequence identity. To screen these matches, corresponding 1 mu s coarse-grained (CG) molecular dynamics trajectories were used to compute normalized root-mean-square-fluctuation graphs and select mobile segments, which were, thereafter, matched for all pairs using unweighted three-dimensional autocorrelation vectors. Our in-house custom-built forcefield (FF), extensively validated against dynamics information obtained from experimental nuclear magnetic resonance data, was specifically used to generate the CG dynamics trajectories. The test for correspondence of dynamics-signature of protein segments and function revealed 87% true positive rate and 93.5% true negative rate, on a dataset of 60 experimentally validated proteins, including moonlighting proteins and those with novel functional motifs. A random test against 315 unique fold/function proteins for a negative test gave >99% true recall. A blind prediction on a novel protein appears consistent with additional evidences retrieved therein. This is the first proof-of-principle of generalized use of structural dynamics for inferring protein molecular function leveraging our custom-made CG FF, useful to all. (C) 2014 Wiley Periodicals, Inc

Pipeline for inferring protein function from dynamics using coarse-grained molecular mechanics forcefield

Dynamics is integral to the function of proteins, yet the use of molecular dynamics (MD) simulation as a technique remains under-explored for molecular function inference. This is more important in the context of genomics projects where novel proteins are determined with limited evolutionary information. Recently we developed a method to match the query protein's flexible segments to infer function using a novel approach combining analysis of residue fluctuation-graphs and auto-correlation vectors derived from coarse-grained (CG) MD trajectory. The method was validated on a diverse dataset with sequence identity between proteins as low as 3%, with high function-recall rates. Here we share its implementation as a publicly accessible web service, named DynFunc (Dynamics Match for Function) to query protein function from >= 1 mu s long CG dynamics trajectory information of protein subunits. Users are provided with the custom-developed coarse-grained molecular mechanics (CGMM) forcefield to generate the MD trajectories for their protein of interest. On upload of trajectory information, the DynFunc web server identifies specific flexible regions of the protein linked to putative molecular function. Our unique application does not use evolutionary information to infer molecular function from MD information and can, therefore, work for all proteins, including moonlighting and the novel ones, whenever structural information is available. Our pipeline is expected to be of utility to all structural biologists working with novel proteins and interested in moonlighting functions

Comparison of Biomolecular Force Fields for Alkanethiol Self-Assembled Monolayer Simulations

Reasonable all-atom or united-atom biomolecular force fields have been developed to represent the properties of proteins and lipid membranes in molecular dynamics simulations. However, since they have not been parametrized for self-assembled monolayers (SAMs), their utility in simulating SAMs and protein–SAM systems has not been confirmed. Here, we compared six popular biomolecular force fields, Lipid14, GAFF, L-OPLS, CHARMM36, Slipids, and GROMOS54a7, to simulate alkanethiol SAMs of short to long chains (C10–C18). Our results show that none of these force fields reproduce the chain length dependence of the tilt angle, and twist angle. Although the droplet contact angles on SAMs are well represented by all force fields, only GAFF and Lipid14 yield phase transition temperatures that are reasonably close to the experimental values. Overall, our comprehensive comparison suggests that GAFF and Lipid14 are better choices for SAM simulations; further improvements in the force field parameters for SAMs are required. (152 words