1
|
Gasparini P, Philot EA, Pantaleão SQ, Torres-Bonfim NESM, Kliousoff A, Quiroz RCN, Perahia D, Simões RP, Magro AJ, Scott AL. Unveiling mutation effects on the structural dynamics of the main protease from SARS-CoV-2 with hybrid simulation methods. J Mol Graph Model 2023; 121:108443. [PMID: 36870228 PMCID: PMC9945984 DOI: 10.1016/j.jmgm.2023.108443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Revised: 02/15/2023] [Accepted: 02/21/2023] [Indexed: 02/25/2023]
Abstract
The main protease of SARS-CoV-2 (called Mpro or 3CLpro) is essential for processing polyproteins encoded by viral RNA. Several Mpro mutations were found in SARS-CoV-2 variants, which are related to higher transmissibility, pathogenicity, and resistance to neutralization antibodies. Macromolecules adopt several favored conformations in solution depending on their structure and shape, determining their dynamics and function. In this study, we used a hybrid simulation method to generate intermediate structures along the six lowest frequency normal modes and sample the conformational space and characterize the structural dynamics and global motions of WT SARS-CoV-2 Mpro and 48 mutations, including mutations found in P.1, B.1.1.7, B.1.351, B.1.525 and B.1.429+B.1.427 variants. We tried to contribute to the elucidation of the effects of mutation in the structural dynamics of SARS-CoV-2 Mpro. A machine learning analysis was performed following the investigation regarding the influence of the K90R, P99L, P108S, and N151D mutations on the dimeric interface assembling of the SARS-CoV-2 Mpro. The parameters allowed the selection of potential structurally stable dimers, which demonstrated that some single surface aa substitutions not located at the dimeric interface (K90R, P99L, P108S, and N151D) are able to induce significant quaternary changes. Furthermore, our results demonstrated, by a Quantum Mechanics method, the influence of SARS-CoV-2 Mpro mutations on the catalytic mechanism, confirming that only one of the chains of the WT and mutant SARS-CoV-2 Mpros are prone to cleave substrates. Finally, it was also possible to identify the aa residue F140 as an important factor related to the increasing enzymatic reactivity of a significant number of SARS-CoV-2 Mpro conformations generated by the normal modes-based simulations.
Collapse
Affiliation(s)
- P Gasparini
- Computational Biology and Biophysics Laboratory, Federal University of ABC - UFABC, Santo André, São Paulo, Brazil
| | - E A Philot
- Computational Biology and Biophysics Laboratory, Federal University of ABC - UFABC, Santo André, São Paulo, Brazil
| | - S Q Pantaleão
- Computational Biology and Biophysics Laboratory, Federal University of ABC - UFABC, Santo André, São Paulo, Brazil
| | - N E S M Torres-Bonfim
- Computational Biology and Biophysics Laboratory, Federal University of ABC - UFABC, Santo André, São Paulo, Brazil
| | - A Kliousoff
- Computational Biology and Biophysics Laboratory, Federal University of ABC - UFABC, Santo André, São Paulo, Brazil
| | - R C N Quiroz
- Computational Biology and Biophysics Laboratory, Federal University of ABC - UFABC, Santo André, São Paulo, Brazil
| | - D Perahia
- École Normale Supérieure Paris-Saclay, LBPA, Scaly, France
| | - R P Simões
- Department of Bioprocesses and Biotechnology, School of Agriculture (FCA), Unesp, Botucatu, São Paulo, Brazil
| | - A J Magro
- Department of Bioprocesses and Biotechnology, School of Agriculture (FCA), Unesp, Botucatu, São Paulo, Brazil; Institute of Biotechnology (IBTEC), Unesp, Botucatu, São Paulo, Brazil
| | - A L Scott
- Computational Biology and Biophysics Laboratory, Federal University of ABC - UFABC, Santo André, São Paulo, Brazil.
| |
Collapse
|
2
|
Quiroz RCN, Philot EA, General IJ, Perahia D, Scott AL. Effect of phosphorylation on the structural dynamics, thermal stability of human dopamine transporter: A simulation study using normal modes, molecular dynamics and Markov State Model. J Mol Graph Model 2023; 118:108359. [PMID: 36279761 DOI: 10.1016/j.jmgm.2022.108359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 10/04/2022] [Accepted: 10/06/2022] [Indexed: 11/06/2022]
Abstract
The Human Dopamine Transporter (hDAT) plays an essential role in modulating the Influx/Efflux of dopamine, and it is involved in the mechanism of certain neurodegenerative diseases such as Parkinson's disease. Several studies have reported important states for Dopamine transport: outward-facing open state (OFo), the outward-facing closed state (OFc), the holo-occluded state closed (holo), and the inward-facing open state (IFo). Furthermore, experimental assays have shown that different phosphorylation conditions in hDAT can affect the rate of dopamine absorption. We present a protocol using hybrid simulation methods to study the conformational dynamics and stability of states of hDAT under different phosphorylation sites. With this protocol, we explored the conformational space of hDAT, identified the states, and evaluated the free energy differences and the transition probabilities between them in each of the phosphorylation cases. We also presented the conformational changes and correlated them with those described in the literature. There is a thesis/hypothesis that the phosphorylation condition corresponding to NP-333 system (where all sites Ser/Thr from residue 2 to 62 and 254 to 613 are phosphorylated, except residue 333) would decrease the rate of dopamine transport from the extracellular medium to the intracellular medium by hDAT as previously described in the literature by Lin et al., 2003. Our results corroborated this thesis/hypothesis and the data reported. It is probably due to the affectation/changes/alteration of the conformational dynamics of this system that makes the intermediate states more likely and makes it difficult to initial states associated with the uptake of dopamine in the extracellular medium, corroborating the experimental results. Furthermore, our results showed that just single phosphorylation/dephosphorylation could alter intrinsic protein motions affecting the sampling of one or more states necessary for dopamine transport. In this sense, the modification of phosphorylation influences protein movements and conformational preferences, affecting the stability of states and the transition between them and, therefore, the transport.
Collapse
Affiliation(s)
- R C N Quiroz
- Biossistemas, Universidade Federal do ABC, CCNH, Santo André, Brazil; Centro de Matemática, Computação e Cognição. Laboratório de Biofísica e Biologia Computacional. Universidade Federal do ABC, Santo André, São Paulo, Brazil
| | - E A Philot
- Centro de Matemática, Computação e Cognição. Laboratório de Biofísica e Biologia Computacional. Universidade Federal do ABC, Santo André, São Paulo, Brazil
| | - I J General
- School of Science and Technology, Universidad Nacional de San Martin, ICIFI and CONICET, 25 de Mayo y Francia, San Martín, 1650, Buenos Aires, Argentina
| | - D Perahia
- Laboratoire de Biologie et Pharmacologie Appliquée, Ecole Normale Supérieure Paris-Saclay, UMR 8113, CNRS, 4 avenue des Sciences, 91190 Gif-sur-Yvette, France
| | - A L Scott
- UFABC - Universidade Federal Do ABC, Centro de Matemática, Computação e Cognição, Laboratório de Biofísica e Biologia Computacional, Brazil.
| |
Collapse
|
3
|
Wang J, Miao Y. Mechanistic Insights into Specific G Protein Interactions with Adenosine Receptors. J Phys Chem B 2019; 123:6462-6473. [PMID: 31283874 DOI: 10.1021/acs.jpcb.9b04867] [Citation(s) in RCA: 75] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Coupling between G-protein-coupled receptors (GPCRs) and the G proteins is a key step in cellular signaling. Despite extensive experimental and computational studies, the mechanism of specific GPCR-G protein coupling remains poorly understood. This has greatly hindered effective drug design of GPCRs that are primary targets of ∼1/3 of currently marketed drugs. Here, we have employed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method to decipher the mechanism of the GPCR-G protein interactions. Adenosine receptors (ARs) were used as model systems based on very recently determined cryo-EM structures of the A1AR and A2AAR coupled with the Gi and Gs proteins, respectively. Changing the Gi protein to the Gs led to increased fluctuations in the A1AR and agonist adenosine (ADO), while agonist 5'-N-ethylcarboxamidoadenosine (NECA) binding in the A2AAR could be still stabilized upon changing the Gs protein to the Gi. Free energy calculations identified one stable low-energy conformation for each of the A1AR-Gi and A2AAR-Gs complexes as in the cryo-EM structures, similarly for the A2AAR-Gi complex. In contrast, the ADO agonist and Gs protein sampled multiple conformations in the A1AR-Gs system. GaMD simulations thus indicated that the A1AR preferred to couple with the Gi protein to the Gs, while the A2AAR could couple with both the Gs and Gi proteins, being highly consistent with experimental findings of the ARs. More importantly, detailed analysis of the atomic simulations showed that the specific AR-G protein coupling resulted from remarkably complementary residue interactions at the protein interface, involving mainly the receptor transmembrane 6 helix and the Gα α5 helix and α4-β6 loop. In summary, the GaMD simulations have provided unprecedented insights into the dynamic mechanism of specific GPCR-G protein interactions at an atomistic level.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences , University of Kansas , Lawrence , Kansas 66047 , United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences , University of Kansas , Lawrence , Kansas 66047 , United States
| |
Collapse
|
4
|
Renault P, Louet M, Marie J, Labesse G, Floquet N. Molecular Dynamics Simulations of the Allosteric Modulation of the Adenosine A2A Receptor by a Mini-G Protein. Sci Rep 2019; 9:5495. [PMID: 30940903 PMCID: PMC6445292 DOI: 10.1038/s41598-019-41980-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2018] [Accepted: 03/19/2019] [Indexed: 01/14/2023] Open
Abstract
Through their coupling to G proteins, G Protein-Coupled Receptors (GPCRs) trigger cellular responses to various signals. Some recent experiments have interestingly demonstrated that the G protein can also act on the receptor by favoring a closed conformation of its orthosteric site, even in the absence of a bound agonist. In this work, we explored such an allosteric modulation by performing extensive molecular dynamics simulations on the adenosine A2 receptor (A2aR) coupled to the Mini-Gs protein. In the presence of the Mini-Gs, we confirmed a restriction of the receptor’s agonist binding site that can be explained by a modulation of the intrinsic network of contacts of the receptor. Of interest, we observed similar effects with the C-terminal helix of the Mini-Gs, showing that the observed effect on the binding pocket results from direct local contacts with the bound protein partner that cause a rewiring of the whole receptor’s interaction network.
Collapse
Affiliation(s)
- Pedro Renault
- Institut des Biomolécules Max Mousseron (IBMM), CNRS UMR5247, Université de Montpellier, ENSCM, 34090, Montpellier, France.,Centre de Biochimie Structurale, Université de Montpellier, CNRS, INSERM, 34090, Montpellier, France
| | - Maxime Louet
- Institut des Biomolécules Max Mousseron (IBMM), CNRS UMR5247, Université de Montpellier, ENSCM, 34090, Montpellier, France
| | - Jacky Marie
- Institut des Biomolécules Max Mousseron (IBMM), CNRS UMR5247, Université de Montpellier, ENSCM, 34090, Montpellier, France
| | - Gilles Labesse
- Centre de Biochimie Structurale, Université de Montpellier, CNRS, INSERM, 34090, Montpellier, France
| | - Nicolas Floquet
- Institut des Biomolécules Max Mousseron (IBMM), CNRS UMR5247, Université de Montpellier, ENSCM, 34090, Montpellier, France.
| |
Collapse
|
5
|
Girola N, Resende-Lara PT, Figueiredo CR, Massaoka MH, Azevedo RA, Cunha RLOR, Polonelli L, Travassos LR. Molecular, Biological and Structural Features of V L CDR-1 Rb44 Peptide, Which Targets the Microtubule Network in Melanoma Cells. Front Oncol 2019; 9:25. [PMID: 30740361 PMCID: PMC6355703 DOI: 10.3389/fonc.2019.00025] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2018] [Accepted: 01/08/2019] [Indexed: 12/21/2022] Open
Abstract
Microtubules are important drug targets in tumor cells, owing to their role in supporting and determining the cell shape, organelle movement and cell division. The complementarity-determining regions (CDRs) of immunoglobulins have been reported to be a source of anti-tumor peptide sequences, independently of the original antibody specificity for a given antigen. We found that, the anti-Lewis B mAb light-chain CDR1 synthetic peptide Rb44, interacted with microtubules and induced depolymerization, with subsequent degradation of actin filaments, leading to depolarization of mitochondrial membrane-potential, increase of ROS, cell cycle arrest at G2/M, cleavage of caspase-9, caspase-3 and PARP, upregulation of Bax and downregulation of Bcl-2, altogether resulting in intrinsic apoptosis of melanoma cells. The in vitro inhibition of angiogenesis was also an Rb44 effect. Peritumoral injection of Rb44L1 delayed growth of subcutaneously grafted melanoma cells in a syngeneic mouse model. L1-CDRs from immunoglobulins and their interactions with tubulin-dimers were explored to interpret effects on microtubule stability. The opening motion of tubulin monomers allowed for efficient L1-CDR docking, impairment of dimer formation and microtubule dissociation. We conclude that Rb44 VL-CDR1 is a novel peptide that acts on melanoma microtubule network causing cell apoptosis in vitro and melanoma growth inhibition in vivo.
Collapse
Affiliation(s)
- Natalia Girola
- Department of Microbiology, Immunology and Parasitology, Experimental Oncology Unit, Federal University of São Paulo, São Paulo, Brazil
| | - Pedro T Resende-Lara
- Computational Biology and Bioinformatics Laboratory, Federal University of ABC, Santo André, Brazil
| | - Carlos R Figueiredo
- Department of Microbiology, Immunology and Parasitology, Experimental Oncology Unit, Federal University of São Paulo, São Paulo, Brazil.,Department of Molecular and Clinical Cancer Medicine, University of Liverpool, Liverpool, United Kingdom
| | | | - Ricardo A Azevedo
- Department of Microbiology, Immunology and Parasitology, Experimental Oncology Unit, Federal University of São Paulo, São Paulo, Brazil
| | - Rodrigo L O R Cunha
- Chemical Biology Laboratory, Natural and Human Sciences Center, Federal University of ABC, Santo André, Brazil
| | - Luciano Polonelli
- Unit of Biomedical, Biotechnological and Translational Sciences, Department of Medicine and Surgery, Universitá degli Studi di Parma, Parma, Italy
| | - Luiz R Travassos
- Department of Microbiology, Immunology and Parasitology, Experimental Oncology Unit, Federal University of São Paulo, São Paulo, Brazil.,Recepta Biopharma, São Paulo, Brazil
| |
Collapse
|
6
|
Hong L, Vani BP, Thiede EH, Rust MJ, Dinner AR. Molecular dynamics simulations of nucleotide release from the circadian clock protein KaiC reveal atomic-resolution functional insights. Proc Natl Acad Sci U S A 2018; 115:E11475-E11484. [PMID: 30442665 PMCID: PMC6298084 DOI: 10.1073/pnas.1812555115] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
The cyanobacterial clock proteins KaiA, KaiB, and KaiC form a powerful system to study the biophysical basis of circadian rhythms, because an in vitro mixture of the three proteins is sufficient to generate a robust ∼24-h rhythm in the phosphorylation of KaiC. The nucleotide-bound states of KaiC critically affect both KaiB binding to the N-terminal domain (CI) and the phosphotransfer reactions that (de)phosphorylate the KaiC C-terminal domain (CII). However, the nucleotide exchange pathways associated with transitions among these states are poorly understood. In this study, we integrate recent advances in molecular dynamics methods to elucidate the structure and energetics of the pathway for Mg·ADP release from the CII domain. We find that nucleotide release is coupled to large-scale conformational changes in the KaiC hexamer. Solvating the nucleotide requires widening the subunit interface leading to the active site, which is linked to extension of the A-loop, a structure implicated in KaiA binding. These results provide a molecular hypothesis for how KaiA acts as a nucleotide exchange factor. In turn, structural parallels between the CI and CII domains suggest a mechanism for allosteric coupling between the domains. We relate our results to structures observed for other hexameric ATPases, which perform diverse functions.
Collapse
Affiliation(s)
- Lu Hong
- Graduate Program in Biophysical Sciences, The University of Chicago, Chicago, IL 60637
| | - Bodhi P Vani
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
| | - Erik H Thiede
- Department of Chemistry, The University of Chicago, Chicago, IL 60637
| | - Michael J Rust
- Department of Molecular Genetics and Cell Biology, The University of Chicago, Chicago, IL 60637;
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- Institute for Genomics and Systems Biology, The University of Chicago, Chicago, IL 60637
| | - Aaron R Dinner
- Department of Chemistry, The University of Chicago, Chicago, IL 60637;
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL 60637
- James Franck Institute, The University of Chicago, Chicago, IL 60637
| |
Collapse
|
7
|
Figueiredo CR, Azevedo RA, Mousdell S, Resende-Lara PT, Ireland L, Santos A, Girola N, Cunha RLOR, Schmid MC, Polonelli L, Travassos LR, Mielgo A. Blockade of MIF-CD74 Signalling on Macrophages and Dendritic Cells Restores the Antitumour Immune Response Against Metastatic Melanoma. Front Immunol 2018; 9:1132. [PMID: 29875777 PMCID: PMC5974174 DOI: 10.3389/fimmu.2018.01132] [Citation(s) in RCA: 96] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2018] [Accepted: 05/04/2018] [Indexed: 12/15/2022] Open
Abstract
Mounting an effective immune response against cancer requires the activation of innate and adaptive immune cells. Metastatic melanoma is the most aggressive form of skin cancer. While immunotherapies have shown a remarkable success in melanoma treatment, patients develop resistance by mechanisms that include the establishment of an immune suppressive tumor microenvironment. Thus, understanding how metastatic melanoma cells suppress the immune system is vital to develop effective immunotherapies against this disease. In this study, we find that macrophages (MOs) and dendritic cells (DCs) are suppressed in metastatic melanoma and that the Ig-CDR-based peptide C36L1 is able to restore MOs and DCs' antitumorigenic and immunogenic functions and to inhibit metastatic growth in lungs. Specifically, C36L1 treatment is able to repolarize M2-like immunosuppressive MOs into M1-like antitumorigenic MOs, and increase the number of immunogenic DCs, and activated cytotoxic T cells, while reducing the number of regulatory T cells and monocytic myeloid-derived suppressor cells in metastatic lungs. Mechanistically, we find that C36L1 directly binds to the MIF receptor CD74 which is expressed on MOs and DCs, disturbing CD74 structural dynamics and inhibiting MIF signaling on these cells. Interfering with MIF-CD74 signaling on MOs and DCs leads to a decrease in the expression of immunosuppressive factors from MOs and an increase in the capacity of DCs to activate cytotoxic T cells. Our findings suggest that interfering with MIF-CD74 immunosuppressive signaling in MOs and DCs, using peptide-based immunotherapy can restore the antitumor immune response in metastatic melanoma. Our study provides the rationale for further development of peptide-based therapies to restore the antitumor immune response in metastatic melanoma.
Collapse
Affiliation(s)
- Carlos R Figueiredo
- Department of Molecular and Clinical Cancer Medicine, University of Liverpool, Liverpool, United Kingdom.,Experimental Oncology Unit (UNONEX), Department of Microbiology, Immunology and Parasitology, Federal University of São Paulo (UNIFESP), São Paulo, Brazil
| | - Ricardo A Azevedo
- Experimental Oncology Unit (UNONEX), Department of Microbiology, Immunology and Parasitology, Federal University of São Paulo (UNIFESP), São Paulo, Brazil
| | - Sasha Mousdell
- Department of Molecular and Clinical Cancer Medicine, University of Liverpool, Liverpool, United Kingdom
| | - Pedro T Resende-Lara
- Laboratory of Computational Biology and Bioinformatics, Federal University of ABC, Santo André, Brazil.,Laboratoire de Biologie et Pharmacologie Appliquées (LBPA), UMR 8113, Ecole Normale Supérieure, Cachan, France
| | - Lucy Ireland
- Department of Molecular and Clinical Cancer Medicine, University of Liverpool, Liverpool, United Kingdom
| | - Almudena Santos
- Department of Molecular and Clinical Cancer Medicine, University of Liverpool, Liverpool, United Kingdom
| | - Natalia Girola
- Experimental Oncology Unit (UNONEX), Department of Microbiology, Immunology and Parasitology, Federal University of São Paulo (UNIFESP), São Paulo, Brazil
| | - Rodrigo L O R Cunha
- Chemical Biology Laboratory, Natural and Human Sciences Center, Federal University of ABC, Santo André, Brazil
| | - Michael C Schmid
- Department of Molecular and Clinical Cancer Medicine, University of Liverpool, Liverpool, United Kingdom
| | - Luciano Polonelli
- Unit of Biomedical, Biotechnological and Translational Sciences, Department of Medicine and Surgery, Universitá degli Studi di Parma, Parma, Italy
| | - Luiz R Travassos
- Experimental Oncology Unit (UNONEX), Department of Microbiology, Immunology and Parasitology, Federal University of São Paulo (UNIFESP), São Paulo, Brazil
| | - Ainhoa Mielgo
- Department of Molecular and Clinical Cancer Medicine, University of Liverpool, Liverpool, United Kingdom
| |
Collapse
|
8
|
Iglesias J, Saen‐oon S, Soliva R, Guallar V. Computational structure‐based drug design: Predicting target flexibility. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1367] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
| | | | | | - Victor Guallar
- Life Science DepartmentBarcelonaSpain
- ICREA, Passeig Lluís Companys 23BarcelonaSpain
| |
Collapse
|
9
|
Meissner GO, de Resende Lara PT, Scott LPB, Braz ASK, Chaves-Moreira D, Matsubara FH, Soares EM, Trevisan-Silva D, Gremski LH, Veiga SS, Chaim OM. Molecular cloning and in silico characterization of knottin peptide, U2-SCRTX-Lit2, from brown spider (Loxosceles intermedia) venom glands. J Mol Model 2016; 22:196. [PMID: 27488102 DOI: 10.1007/s00894-016-3067-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2016] [Accepted: 07/10/2016] [Indexed: 01/16/2023]
Abstract
Inhibitor cystine knots (ICKs) are a family of structural peptides with a large number of cysteine residues that form intramolecular disulfide bonds, resulting in a knot. These peptides are involved in a variety of biological functions including predation and defense, and are found in various species, such as spiders, scorpions, sea anemones, and plants. The Loxosceles intermedia venom gland transcriptome identified five groups of ICK peptides that represent more than 50 % of toxin-coding transcripts. Here, we describe the molecular cloning of U2-Sicaritoxin-Lit2 (U2-SCRTX-Lit2), bioinformatic characterization, structure prediction, and molecular dynamic analysis. The sequence of U2-SCRTX-Lit2 obtained from the transcriptome is similar to that of μ-Hexatoxin-Mg2, a peptide that inhibits the insect Nav channel. Bioinformatic analysis of sequences classified as ICK family members also showed a conservation of cysteine residues among ICKs from different spiders, with the three dimensional molecular model of U2-SCRTX-Lit2 similar in structure to the hexatoxin from μ-hexatoxin-Mg2a. Molecular docking experiments showed the interaction of U2-SCRTX-Lit2 to its predictable target-the Spodoptera litura voltage-gated sodium channel (SlNaVSC). After 200 ns of molecular dynamic simulation, the final structure of the complex showed stability in agreement with the experimental data. The above analysis corroborates the existence of a peptide toxin with insecticidal activity from a novel ICK family in L. intermedia venom and demonstrates that this peptide targets Nav channels.
Collapse
Affiliation(s)
- Gabriel Otto Meissner
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil
| | - Pedro Túlio de Resende Lara
- Laboratory of Computational Biology and Bioinformatics, Federal University of ABC, Santo André, São Paulo, Brazil
| | - Luis Paulo Barbour Scott
- Laboratory of Computational Biology and Bioinformatics, Federal University of ABC, Santo André, São Paulo, Brazil
| | - Antônio Sérgio Kimus Braz
- Laboratory of Computational Biology and Bioinformatics, Federal University of ABC, Santo André, São Paulo, Brazil
| | - Daniele Chaves-Moreira
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil
| | - Fernando Hitomi Matsubara
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil
| | - Eduardo Mendonça Soares
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil
| | - Dilza Trevisan-Silva
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil
| | - Luiza Helena Gremski
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil.,Laboratory of Molecular Immunopathology, Department of Clinical Pathology, Clinical Hospital of Federal University of Paraná, Curitiba, Paraná, Brazil
| | - Silvio Sanches Veiga
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil
| | - Olga Meiri Chaim
- Department of Cell Biology, Federal University of Paraná, Jardim das Américas, 81531-990, Curitiba, Paraná, Brazil.
| |
Collapse
|
10
|
Computational Simulation of the Activation Cycle of Gα Subunit in the G Protein Cycle Using an Elastic Network Model. PLoS One 2016; 11:e0159528. [PMID: 27483005 PMCID: PMC4970668 DOI: 10.1371/journal.pone.0159528] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2016] [Accepted: 07/04/2016] [Indexed: 01/13/2023] Open
Abstract
Agonist-activated G protein-coupled receptors (GPCRs) interact with GDP-bound G protein heterotrimers (Gαβγ) promoting GDP/GTP exchange, which results in dissociation of Gα from the receptor and Gβγ. The GTPase activity of Gα hydrolyzes GTP to GDP, and the GDP-bound Gα interacts with Gβγ, forming a GDP-bound G protein heterotrimer. The G protein cycle is allosterically modulated by conformational changes of the Gα subunit. Although biochemical and biophysical methods have elucidated the structure and dynamics of Gα, the precise conformational mechanisms underlying the G protein cycle are not fully understood yet. Simulation methods could help to provide additional details to gain further insight into G protein signal transduction mechanisms. In this study, using the available X-ray crystal structures of Gα, we simulated the entire G protein cycle and described not only the steric features of the Gα structure, but also conformational changes at each step. Each reference structure in the G protein cycle was modeled as an elastic network model and subjected to normal mode analysis. Our simulation data suggests that activated receptors trigger conformational changes of the Gα subunit that are thermodynamically favorable for opening of the nucleotide-binding pocket and GDP release. Furthermore, the effects of GTP binding and hydrolysis on mobility changes of the C and N termini and switch regions are elucidated. In summary, our simulation results enabled us to provide detailed descriptions of the structural and dynamic features of the G protein cycle.
Collapse
|
11
|
Abstract
Molecular dynamics (MD) simulations have become a powerful and popular method for the study of protein allostery, the widespread phenomenon in which a stimulus at one site on a protein influences the properties of another site on the protein. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric binding sites and the determination of the mechanistic basis for allostery. These results can provide a foundation for applications including rational drug design and protein engineering. Here, we provide an introduction to the investigation of protein allostery using molecular dynamics simulation. We emphasize the importance of designing simulations that include appropriate perturbations to the molecular system, such as the addition or removal of ligands or the application of mechanical force. We also demonstrate how the bidirectional nature of allostery-the fact that the two sites involved influence one another in a symmetrical manner-can facilitate such investigations. Through a series of case studies, we illustrate how these concepts have been used to reveal the structural basis for allostery in several proteins and protein complexes of biological and pharmaceutical interest.
Collapse
|
12
|
Floquet N, Costa MGS, Batista PR, Renault P, Bisch PM, Raussin F, Martinez J, Morris MC, Perahia D. Conformational Equilibrium of CDK/Cyclin Complexes by Molecular Dynamics with Excited Normal Modes. Biophys J 2015; 109:1179-89. [PMID: 26255588 PMCID: PMC4576171 DOI: 10.1016/j.bpj.2015.07.003] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2015] [Revised: 06/16/2015] [Accepted: 07/01/2015] [Indexed: 01/21/2023] Open
Abstract
Cyclin-dependent kinases (CDKs) and their associated regulatory cyclins are central for timely regulation of cell-cycle progression. They constitute attractive pharmacological targets for development of anticancer therapeutics, since they are frequently deregulated in human cancers and contribute to sustained, uncontrolled tumor proliferation. Characterization of their structural/dynamic features is essential to gain in-depth insight into structure-activity relationships. In addition, the identification of druggable pockets or key intermediate conformations yields potential targets for the development of novel classes of inhibitors. Structural studies of CDK2/cyclin A have provided a wealth of information concerning monomeric/heterodimeric forms of this kinase. There is, however, much less structural information for other CDK/cyclin complexes, including CDK4/cyclin D1, which displays an alternative (open) position of the cyclin partner relative to CDK, contrasting with the closed CDK2/cyclin A conformation. In this study, we carried out normal-mode analysis and enhanced sampling simulations with our recently developed method, molecular dynamics with excited normal modes, to understand the conformational equilibrium on these complexes. Interestingly, the lowest-frequency normal mode computed for each complex described the transition between the open and closed conformations. Exploration of these motions with an explicit-solvent representation using molecular dynamics with excited normal modes confirmed that the closed conformation is the most stable for the CDK2/cyclin A complex, in agreement with their experimentally available structures. On the other hand, we clearly show that an open↔closed equilibrium may exist in CDK4/cyclin D1, with closed conformations resembling that captured for CDK2/cyclin A. Such conformational preferences may result from the distinct distributions of frustrated contacts in each complex. Using the same approach, the putative roles of the Thr(160) phosphoryl group and the T-loop conformation were investigated. These results provide a dynamic view of CDKs revealing intermediate conformations not yet characterized for CDK members other than CDK2, which will be useful for the design of inhibitors targeting critical conformational transitions.
Collapse
Affiliation(s)
- Nicolas Floquet
- Institut des Biomolécules Max Mousseron (IBMM), Centre National de la Recherche Scientifique UMR 5247, Université de Montpellier, Ecole Normale Supérieure de Chimie de Montpellier, Faculté de Pharmacie, Montpellier, France.
| | - Mauricio G S Costa
- Programa de Computação Científica, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Paulo R Batista
- Programa de Computação Científica, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Pedro Renault
- Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | - Paulo M Bisch
- Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | - Florent Raussin
- Institut des Biomolécules Max Mousseron (IBMM), Centre National de la Recherche Scientifique UMR 5247, Université de Montpellier, Ecole Normale Supérieure de Chimie de Montpellier, Faculté de Pharmacie, Montpellier, France
| | - Jean Martinez
- Institut des Biomolécules Max Mousseron (IBMM), Centre National de la Recherche Scientifique UMR 5247, Université de Montpellier, Ecole Normale Supérieure de Chimie de Montpellier, Faculté de Pharmacie, Montpellier, France
| | - May C Morris
- Institut des Biomolécules Max Mousseron (IBMM), Centre National de la Recherche Scientifique UMR 5247, Université de Montpellier, Ecole Normale Supérieure de Chimie de Montpellier, Faculté de Pharmacie, Montpellier, France
| | - David Perahia
- Laboratoire de Biologie et de Pharmacologie Appliquée, Ecole Normale Supérieure de Cachan, Centre National de la Recherche Scientifique, Cachan, France.
| |
Collapse
|
13
|
Louet M, Seifert C, Hensen U, Gräter F. Dynamic Allostery of the Catabolite Activator Protein Revealed by Interatomic Forces. PLoS Comput Biol 2015; 11:e1004358. [PMID: 26244893 PMCID: PMC4526232 DOI: 10.1371/journal.pcbi.1004358] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2014] [Accepted: 05/28/2015] [Indexed: 11/23/2022] Open
Abstract
The Catabolite Activator Protein (CAP) is a showcase example for entropic allostery. For full activation and DNA binding, the homodimeric protein requires the binding of two cyclic AMP (cAMP) molecules in an anti-cooperative manner, the source of which appears to be largely of entropic nature according to previous experimental studies. We here study at atomic detail the allosteric regulation of CAP with Molecular dynamics (MD) simulations. We recover the experimentally observed entropic penalty for the second cAMP binding event with our recently developed force covariance entropy estimator and reveal allosteric communication pathways with Force Distribution Analyses (FDA). Our observations show that CAP binding results in characteristic changes in the interaction pathways connecting the two cAMP allosteric binding sites with each other, as well as with the DNA binding domains. We identified crucial relays in the mostly symmetric allosteric activation network, and suggest point mutants to test this mechanism. Our study suggests inter-residue forces, as opposed to coordinates, as a highly sensitive measure for structural adaptations that, even though minute, can very effectively propagate allosteric signals. The Catabolite Activator Protein (CAP) is a well-studied example for how cellular catabolite levels are integrated into the gene regulation. Its affinity for a specific stretch of DNA can be switched on by the binding of two nucleotide molecules termed cAMP to its two protomers. Even though the nucleotides occupy structurally identical binding pockets, the second cAMP binding occurs at an affinity orders of magnitude lower than the first cAMP binding. The question arises how, in the absence of structural changes, the first binding can affect the second. An answer from experiments has been that the communication is largely of entropic nature, i.e. the second cAMP binding would lead to a pronounced reduction in atomic fluctuations of the protein without affecting the atomic mean positions. We here revisited this question by performing Molecular Dynamics simulations. By measuring correlations of forces, a newly derived method outperforming the more common coordinate-based approach, we could recover the previously determined entropic penalty. In addition, however, we observed unobtrusive structural changes of side-chain interactions leading to the occlusion of the second binding pocket that add a critical ‘enthalpic’ component hitherto overlooked. Our study provides a mechanistic view onto the intriguing anti-cooperativity of CAP.
Collapse
Affiliation(s)
- Maxime Louet
- Heidelberg Institutes for Theoretical Studies gGmbH, Heidelberg, Germany
| | - Christian Seifert
- Heidelberg Institutes for Theoretical Studies gGmbH, Heidelberg, Germany
| | - Ulf Hensen
- Eidgenössische Technische Hochschule Zürich, Department of Biosystem Science and Engineering, Basel, Switzerland
| | - Frauke Gräter
- Heidelberg Institutes for Theoretical Studies gGmbH, Heidelberg, Germany; CAS-MPG Partner Institute and Key Laboratory for Computational Biology, Shanghai, China
| |
Collapse
|
14
|
Dror RO, Mildorf TJ, Hilger D, Manglik A, Borhani DW, Arlow DH, Philippsen A, Villanueva N, Yang Z, Lerch MT, Hubbell WL, Kobilka BK, Sunahara RK, Shaw DE. SIGNAL TRANSDUCTION. Structural basis for nucleotide exchange in heterotrimeric G proteins. Science 2015; 348:1361-5. [PMID: 26089515 DOI: 10.1126/science.aaa5264] [Citation(s) in RCA: 208] [Impact Index Per Article: 23.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
G protein-coupled receptors (GPCRs) relay diverse extracellular signals into cells by catalyzing nucleotide release from heterotrimeric G proteins, but the mechanism underlying this quintessential molecular signaling event has remained unclear. Here we use atomic-level simulations to elucidate the nucleotide-release mechanism. We find that the G protein α subunit Ras and helical domains-previously observed to separate widely upon receptor binding to expose the nucleotide-binding site-separate spontaneously and frequently even in the absence of a receptor. Domain separation is necessary but not sufficient for rapid nucleotide release. Rather, receptors catalyze nucleotide release by favoring an internal structural rearrangement of the Ras domain that weakens its nucleotide affinity. We use double electron-electron resonance spectroscopy and protein engineering to confirm predictions of our computationally determined mechanism.
Collapse
Affiliation(s)
- Ron O Dror
- D. E. Shaw Research, New York, NY 10036, USA.
| | | | - Daniel Hilger
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Aashish Manglik
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305, USA
| | | | | | | | - Nicolas Villanueva
- Department of Pharmacology, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - Zhongyu Yang
- Jules Stein Eye Institute and Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095, USA
| | - Michael T Lerch
- Jules Stein Eye Institute and Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095, USA
| | - Wayne L Hubbell
- Jules Stein Eye Institute and Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095, USA
| | - Brian K Kobilka
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Roger K Sunahara
- Department of Pharmacology, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - David E Shaw
- D. E. Shaw Research, New York, NY 10036, USA. Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA.
| |
Collapse
|
15
|
Pachov DV, van den Bedem H. Nullspace Sampling with Holonomic Constraints Reveals Molecular Mechanisms of Protein Gαs. PLoS Comput Biol 2015; 11:e1004361. [PMID: 26218073 PMCID: PMC4517867 DOI: 10.1371/journal.pcbi.1004361] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2014] [Accepted: 05/22/2015] [Indexed: 11/19/2022] Open
Abstract
Proteins perform their function or interact with partners by exchanging between conformational substates on a wide range of spatiotemporal scales. Structurally characterizing these exchanges is challenging, both experimentally and computationally. Large, diffusional motions are often on timescales that are difficult to access with molecular dynamics simulations, especially for large proteins and their complexes. The low frequency modes of normal mode analysis (NMA) report on molecular fluctuations associated with biological activity. However, NMA is limited to a second order expansion about a minimum of the potential energy function, which limits opportunities to observe diffusional motions. By contrast, kino-geometric conformational sampling (KGS) permits large perturbations while maintaining the exact geometry of explicit conformational constraints, such as hydrogen bonds. Here, we extend KGS and show that a conformational ensemble of the α subunit Gαs of heterotrimeric stimulatory protein Gs exhibits structural features implicated in its activation pathway. Activation of protein Gs by G protein-coupled receptors (GPCRs) is associated with GDP release and large conformational changes of its α-helical domain. Our method reveals a coupled α-helical domain opening motion while, simultaneously, Gαs helix α5 samples an activated conformation. These motions are moderated in the activated state. The motion centers on a dynamic hub near the nucleotide-binding site of Gαs, and radiates to helix α4. We find that comparative NMA-based ensembles underestimate the amplitudes of the motion. Additionally, the ensembles fall short in predicting the accepted direction of the full activation pathway. Taken together, our findings suggest that nullspace sampling with explicit, holonomic constraints yields ensembles that illuminate molecular mechanisms involved in GDP release and protein Gs activation, and further establish conformational coupling between key structural elements of Gαs.
Collapse
Affiliation(s)
- Dimitar V. Pachov
- Department of Chemistry, Stanford University, Stanford, California, United States of America
| | - Henry van den Bedem
- Joint Center for Structural Genomics, Stanford Synchrotron Radiation Lightsource, Stanford University, Stanford, California, United States of America
| |
Collapse
|
16
|
Duc NM, Kim HR, Chung KY. Structural mechanism of G protein activation by G protein-coupled receptor. Eur J Pharmacol 2015; 763:214-22. [PMID: 25981300 DOI: 10.1016/j.ejphar.2015.05.016] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2015] [Revised: 03/03/2015] [Accepted: 05/11/2015] [Indexed: 12/17/2022]
Abstract
G protein-coupled receptors (GPCRs) are a family of membrane receptors that regulate physiology and pathology of various organs. Consequently, about 40% of drugs in the market targets GPCRs. Heterotrimeric G proteins are composed of α, β, and γ subunits, and act as the key downstream signaling molecules of GPCRs. The structural mechanism of G protein activation by GPCRs has been of a great interest, and a number of biochemical and biophysical studies have been performed since the late 80's. These studies investigated the interface between GPCR and G proteins and the structural mechanism of GPCR-induced G protein activation. Recently, arrestins are also reported to be important molecular switches in GPCR-mediated signal transduction, and the physiological output of arrestin-mediated signal transduction is different from that of G protein-mediated signal transduction. Understanding the structural mechanism of the activation of G proteins and arrestins would provide fundamental information for the downstream signaling-selective GPCR-targeting drug development. This review will discuss the structural mechanism of GPCR-induced G protein activation by comparing previous biochemical and biophysical studies.
Collapse
Affiliation(s)
- Nguyen Minh Duc
- School of Pharmacy, Sungkyunkwan University, 2066 Seobu-ro, Jangan-gu, Suwon 440-746, Republic of Korea
| | - Hee Ryung Kim
- School of Pharmacy, Sungkyunkwan University, 2066 Seobu-ro, Jangan-gu, Suwon 440-746, Republic of Korea
| | - Ka Young Chung
- School of Pharmacy, Sungkyunkwan University, 2066 Seobu-ro, Jangan-gu, Suwon 440-746, Republic of Korea.
| |
Collapse
|
17
|
Yao XQ, Grant BJ. Domain-opening and dynamic coupling in the α-subunit of heterotrimeric G proteins. Biophys J 2014; 105:L08-10. [PMID: 23870276 DOI: 10.1016/j.bpj.2013.06.006] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2013] [Revised: 06/03/2013] [Accepted: 06/05/2013] [Indexed: 11/30/2022] Open
Abstract
Heterotrimeric G proteins are conformational switches that turn on intracellular signaling cascades in response to the activation of G-protein-coupled receptors. Receptor activation by extracellular stimuli promotes a cycle of GTP binding and hydrolysis on the G protein α-subunit (Gα). Important conformational transitions occurring during this cycle have been characterized from extensive crystallographic studies of Gα. However, the link between the observed conformations and the mechanisms involved in G-protein activation and effector interaction remain unclear. Here we describe a comprehensive principal component analysis of available Gα crystallographic structures supplemented with extensive unbiased conventional and accelerated molecular dynamics simulations that together characterize the response of Gα to GTP binding and hydrolysis. Our studies reveal details of activating conformational changes as well as the intrinsic flexibility of the α-helical domain that includes a large-scale 60° domain opening under nucleotide-free conditions. This result is consistent with the recently reported open crystal structure of Gs, the stimulatory G protein for adenylyl cyclase, in complex with the α2 adrenergic receptor. Sets of unique interactions potentially important for the conformational transition are also identified. Moreover simulations reveal nucleotide-dependent dynamical couplings of distal regions and residues potentially important for the allosteric link between functional sites.
Collapse
Affiliation(s)
- Xin-Qiu Yao
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA
| | | |
Collapse
|
18
|
Alexander NS, Preininger AM, Kaya AI, Stein RA, Hamm HE, Meiler J. Energetic analysis of the rhodopsin-G-protein complex links the α5 helix to GDP release. Nat Struct Mol Biol 2014; 21:56-63. [PMID: 24292645 PMCID: PMC3947367 DOI: 10.1038/nsmb.2705] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2013] [Accepted: 10/02/2013] [Indexed: 01/19/2023]
Abstract
We present a model of interaction of Gi protein with the activated receptor (R*) rhodopsin, which pinpoints energetic contributions to activation and reconciles the β2 adrenergic receptor-Gs crystal structure with new and previously published experimental data. In silico analysis demonstrated energetic changes when the Gα C-terminal helix (α5) interacts with the R* cytoplasmic pocket, thus leading to displacement of the helical domain and GDP release. The model features a less dramatic domain opening compared with the crystal structure. The α5 helix undergoes a 63° rotation, accompanied by a 5.7-Å translation, that reorganizes interfaces between α5 and α1 helices and between α5 and β6-α5. Changes in the β6-α5 loop displace αG. All of these movements lead to opening of the GDP-binding pocket. The model creates a roadmap for experimental studies of receptor-mediated G-protein activation.
Collapse
Affiliation(s)
- Nathan S Alexander
- 1] Department of Chemistry, Vanderbilt University, Nashville, Tennessee, USA. [2]
| | - Anita M Preininger
- 1] Department of Pharmacology, Vanderbilt University, Nashville, Tennessee, USA. [2]
| | - Ali I Kaya
- Department of Pharmacology, Vanderbilt University, Nashville, Tennessee, USA
| | - Richard A Stein
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, Tennessee, USA
| | - Heidi E Hamm
- Department of Pharmacology, Vanderbilt University, Nashville, Tennessee, USA
| | - Jens Meiler
- 1] Department of Chemistry, Vanderbilt University, Nashville, Tennessee, USA. [2] Department of Pharmacology, Vanderbilt University, Nashville, Tennessee, USA
| |
Collapse
|
19
|
Decipher the mechanisms of protein conformational changes induced by nucleotide binding through free-energy landscape analysis: ATP binding to Hsp70. PLoS Comput Biol 2013; 9:e1003379. [PMID: 24348227 PMCID: PMC3861046 DOI: 10.1371/journal.pcbi.1003379] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2013] [Accepted: 10/15/2013] [Indexed: 11/30/2022] Open
Abstract
ATP regulates the function of many proteins in the cell by transducing its binding and hydrolysis energies into protein conformational changes by mechanisms which are challenging to identify at the atomic scale. Based on molecular dynamics (MD) simulations, a method is proposed to analyze the structural changes induced by ATP binding to a protein by computing the effective free-energy landscape (FEL) of a subset of its coordinates along its amino-acid sequence. The method is applied to characterize the mechanism by which the binding of ATP to the nucleotide-binding domain (NBD) of Hsp70 propagates a signal to its substrate-binding domain (SBD). Unbiased MD simulations were performed for Hsp70-DnaK chaperone in nucleotide-free, ADP-bound and ATP-bound states. The simulations revealed that the SBD does not interact with the NBD for DnaK in its nucleotide-free and ADP-bound states whereas the docking of the SBD was found in the ATP-bound state. The docked state induced by ATP binding found in MD is an intermediate state between the initial nucleotide-free and final ATP-bound states of Hsp70. The analysis of the FEL projected along the amino-acid sequence permitted to identify a subset of 27 protein internal coordinates corresponding to a network of 91 key residues involved in the conformational change induced by ATP binding. Among the 91 residues, 26 are identified for the first time, whereas the others were shown relevant for the allosteric communication of Hsp70 s in several experiments and bioinformatics analysis. The FEL analysis revealed also the origin of the ATP-induced structural modifications of the SBD recently measured by Electron Paramagnetic Resonance. The pathway between the nucleotide-free and the intermediate state of DnaK was extracted by applying principal component analysis to the subset of internal coordinates describing the transition. The methodology proposed is general and could be applied to analyze allosteric communication in other proteins. The precise biophysical characterization of the mechanisms of the protein conformational changes controlled by a nucleotide remains a challenge in biology. Molecular dynamics simulations of proteins in different nucleotide-binding states contain information on the nucleotide-dependent conformational dynamics. However, it is difficult to extract relevant information about the conformation-induced mechanism from the raw molecular dynamics data. Herein, we addressed this issue for the major ATP-dependent molecular chaperones Hsp70 s, which contribute to crucial cellular processes and are involved in several neurodegenerative diseases and in cancer. To function, Hsp70 undergoes several conformational changes controlled by the state of its nucleotide-binding domain. We demonstrated that the analysis of the effective free-energy landscape of the protein projected along the amino-acid sequence and computed from the molecular dynamics simulations of Hsp70 in different nucleotide-binding states, holds the key to identify the key residues of the conformational induced pathway. Identification of the key residues involved in the propagation of the structural changes induced by ATP binding offer alternative druggable specific sites other than the ligand binding clefts. The methodology developed for Hsp70 is general and can be adapted to any ligand induced conformational change in proteins.
Collapse
|
20
|
Louet M, Karakas E, Perret A, Perahia D, Martinez J, Floquet N. Conformational restriction of G-proteins Coupled Receptors (GPCRs) upon complexation to G-proteins: a putative activation mode of GPCRs? FEBS Lett 2013; 587:2656-61. [PMID: 23851072 DOI: 10.1016/j.febslet.2013.06.052] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2013] [Revised: 06/13/2013] [Accepted: 06/29/2013] [Indexed: 10/26/2022]
Abstract
GPCRs undergo large conformational changes during their activation. Starting from existing X-ray structures, we used Normal Modes Analyses to study the collective motions of the agonist-bound β2-adrenergic receptor both in its isolated "uncoupled" and G-protein "coupled" conformations. We interestingly observed that the receptor was able to adopt only one major motion in the protein:protein complex. This motion corresponded to an anti-symmetric rotation of both its extra- and intra-cellular parts, with a key role of previously identified highly conserved proline residues. Because this motion was also retrieved when performing NMA on 7 other GPCRs which structures were available, it is strongly suspected to possess a significant biological role, possibly being the "activation mode" of a GPCR when coupled to G-proteins.
Collapse
Affiliation(s)
- Maxime Louet
- Institut des Biomolécules Max Mousseron (IBMM, CNRS UMR5247), Faculté de Pharmacie, 15 avenue Charles Flahault, BP 14491, 34093 Montpellier Cedex 05, France
| | | | | | | | | | | |
Collapse
|
21
|
Preininger AM, Meiler J, Hamm HE. Conformational flexibility and structural dynamics in GPCR-mediated G protein activation: a perspective. J Mol Biol 2013; 425:2288-98. [PMID: 23602809 DOI: 10.1016/j.jmb.2013.04.011] [Citation(s) in RCA: 80] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2013] [Revised: 04/08/2013] [Accepted: 04/10/2013] [Indexed: 12/16/2022]
Abstract
Structure and dynamics of G proteins and their cognate receptors, both alone and in complex, are becoming increasingly accessible to experimental techniques. Understanding the conformational changes and timelines that govern these changes can lead to new insights into the processes of ligand binding and associated G protein activation. Experimental systems may involve the use of, or otherwise stabilize, non-native environments. This can complicate our understanding of structural and dynamic features of processes such as the ionic lock, tryptophan toggle, and G protein flexibility. While elements in the receptor's transmembrane helices and the C-terminal α5 helix of Gα undergo well-defined structural changes, regions subject to conformational flexibility may be important in fine-tuning the interactions between activated receptors and G proteins. The pairing of computational and experimental approaches will continue to provide powerful tools to probe the conformation and dynamics of receptor-mediated G protein activation.
Collapse
Affiliation(s)
- Anita M Preininger
- Department of Pharmacology, Vanderbilt University Medical Center, Nashville, TN 37232-6600, USA.
| | | | | |
Collapse
|
22
|
Louet M, Charlier L, Martinez J, Floquet N. Dissociation of membrane-anchored heterotrimeric G-protein induced by G(α) subunit binding to GTP. J Chem Inf Model 2012; 52:3022-7. [PMID: 23094679 DOI: 10.1021/ci3003717] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Heterotrimeric G-proteins' activation on the intracellular side of the cell membrane is initiated by stimulation of the G-Protein Coupled Receptors (GPCRs) extra-cellular part. This two-step activation mechanism includes (1) an exchange between GDP and GTP molecules in the G(α) subunit and (2) a dissociation of the whole G(αβγ) complex into two membrane-anchored blocks, namely the isolated G(α) and G(βγ) subunits. Although X-ray data are available for both inactive G(αβγ):GDP and active G(α):GTP complexes, intermediate steps involved in the molecular mechanism of the dissociation have not yet been addressed at the molecular level. In this study, we first built a membrane-anchored intermediate G(iαβγ):GTP complex. This model was then equilibrated by molecular dynamics simulations before the Targeted Molecular Dynamics (TMD) technique was used to force the G(α) subunit to evolve from its inactive (GDP-bound) to its active (GTP-bound) conformations, as described by available X-ray data. The TMD constraint was applied only to the G(α) subunit so that the resulting global rearrangements acting on the whole G(αβγ):GTP heterotrimer could be analyzed. We showed how these mainly local conformational changes of G(α) could initiate large domain:domain motions of the whole complex, the G(βγ) behaving as an almost quasi-rigid block. This separation of the two G(α):GTP and G(βγ) subunits required the loss of several interactions at the G(α):G(βγ) interface that were reported. This study provided an atomistic view of the crucial intermediate step of the G-proteins activation, e.g., the dissociation, that could hardly be elucidated by the experiment.
Collapse
Affiliation(s)
- Maxime Louet
- Institut des Biomolécules Max Mousseron (IBMM), CNRS UMR5247, Université Montpellier 1, Université Montpellier 2, Faculté de Pharmacie, 15 av. Charles Flahault, BP 14491, 34093 Montpellier Cedex 05, France
| | | | | | | |
Collapse
|
23
|
Louet M, Martinez J, Floquet N. GDP release preferentially occurs on the phosphate side in heterotrimeric G-proteins. PLoS Comput Biol 2012; 8:e1002595. [PMID: 22829757 PMCID: PMC3400569 DOI: 10.1371/journal.pcbi.1002595] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2012] [Accepted: 05/21/2012] [Indexed: 11/18/2022] Open
Abstract
After extra-cellular stimulation of G-Protein Coupled Receptors (GPCRs), GDP/GTP exchange appears as the key, rate limiting step of the intracellular activation cycle of heterotrimeric G-proteins. Despite the availability of a large number of X-ray structures, the mechanism of GDP release out of heterotrimeric G-proteins still remains unknown at the molecular level. Starting from the available X-ray structure, extensive unconstrained/constrained molecular dynamics simulations were performed on the complete membrane-anchored Gi heterotrimer complexed to GDP, for a total simulation time overcoming 500 ns. By combining Targeted Molecular Dynamics (TMD) and free energy profiles reconstruction by umbrella sampling, our data suggest that the release of GDP was much more favored on its phosphate side. Interestingly, upon the forced extraction of GDP on this side, the whole protein encountered large, collective motions in perfect agreement with those we described previously including a domain to domain motion between the two ras-like and helical sub-domains of Gα. Despite the availability of many structural and biochemical data, the activation of G-proteins remains to be understood at the molecular level. We used a computation tool to decipher the first limiting step of this activation: GDP release. Combining different methods of analysis, we propose that the GDP exit occurs on its phosphate side. This study helped to rationalize some experimental observations from the literature and opens many perspectives concerning the study of G-proteins activation and their putative inhibition.
Collapse
Affiliation(s)
| | | | - Nicolas Floquet
- Institut des Biomolécules Max Mousseron (IBMM), CNRS UMR5247, Université Montpellier 1, Université Montpellier 2, Faculté de Pharmacie, Montpellier, France
- * E-mail:
| |
Collapse
|
24
|
Ho BK, Perahia D, Buckle AM. Hybrid approaches to molecular simulation. Curr Opin Struct Biol 2012; 22:386-93. [PMID: 22633678 DOI: 10.1016/j.sbi.2012.05.005] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2012] [Revised: 05/07/2012] [Accepted: 05/08/2012] [Indexed: 10/28/2022]
Abstract
Molecular dynamics (MD) simulation is an established method for studying the conformational changes that are important for protein function. Recent advances in hardware and software have allowed MD simulations over the same timescales as experiment, improving the agreement between theory and experiment to a large extent. However, running such simulations are costly, in terms of resources, storage, and trajectory analysis. There is still a place for techniques that involve short MD simulations. In order to overcome the sampling paucity of short time-scales, hybrid methods that include some form of MD simulation can exploit certain features of the system of interest, often combining experimental information in surprising ways. Here, we review some recent hybrid approaches to the simulation of proteins.
Collapse
Affiliation(s)
- Bosco K Ho
- Department of Biochemistry and Molecular Biology, Monash University, Clayton, Victoria 3800, Australia
| | | | | |
Collapse
|