1
|
Qi Z, Lu P, Long X, Cao X, Wu M, Xin K, Xue T, Gao X, Huang Y, Wang Q, Jiang C, Xu JR, Liu H. Adaptive advantages of restorative RNA editing in fungi for resolving survival-reproduction trade-offs. SCIENCE ADVANCES 2024; 10:eadk6130. [PMID: 38181075 PMCID: PMC10776026 DOI: 10.1126/sciadv.adk6130] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 12/05/2023] [Indexed: 01/07/2024]
Abstract
RNA editing in various organisms commonly restores RNA sequences to their ancestral state, but its adaptive advantages are debated. In fungi, restorative editing corrects premature stop codons in pseudogenes specifically during sexual reproduction. We characterized 71 pseudogenes and their restorative editing in Fusarium graminearum, demonstrating that restorative editing of 16 pseudogenes is crucial for germ tissue development in fruiting bodies. Our results also revealed that the emergence of premature stop codons is facilitated by restorative editing and that premature stop codons corrected by restorative editing are selectively favored over ancestral amino acid codons. Furthermore, we found that ancestral versions of pseudogenes have antagonistic effects on reproduction and survival. Restorative editing eliminates the survival costs of reproduction caused by antagonistic pleiotropy and provides a selective advantage in fungi. Our findings highlight the importance of restorative editing in the evolution of fungal complex multicellularity and provide empirical evidence that restorative editing serves as an adaptive mechanism enabling the resolution of genetic trade-offs.
Collapse
Affiliation(s)
- Zhaomei Qi
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Ping Lu
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Xinyuan Long
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Xinyu Cao
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Mengchun Wu
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Kaiyun Xin
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Tuan Xue
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Xinlong Gao
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Yi Huang
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Qinhu Wang
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Cong Jiang
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| | - Jin-Rong Xu
- Department of Botany and Plant Pathology, Purdue University, West Lafayette, IN 47907, USA
| | - Huiquan Liu
- State Key Laboratory for Crop Stress Resistance and High-Efficiency Production, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China
| |
Collapse
|
2
|
Imai H, Miyoshi T, Murakami R, Ito K, Ishino Y, Uchiumi T. Functional role of the C-terminal tail of the archaeal ribosomal stalk in recruitment of two elongation factors to the sarcin/ricin loop of 23S rRNA. Genes Cells 2015; 20:613-24. [DOI: 10.1111/gtc.12256] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2015] [Accepted: 04/27/2015] [Indexed: 11/28/2022]
Affiliation(s)
- Hirotatsu Imai
- Department of Biology; Faculty of Science; Niigata University; 8050 Ikarashi 2-no-cho Nishi-ku Niigata 950-2181 Japan
| | - Tomohiro Miyoshi
- Department of Biology; Faculty of Science; Niigata University; 8050 Ikarashi 2-no-cho Nishi-ku Niigata 950-2181 Japan
| | - Ryo Murakami
- Department of Biology; Faculty of Science; Niigata University; 8050 Ikarashi 2-no-cho Nishi-ku Niigata 950-2181 Japan
| | - Kosuke Ito
- Department of Biology; Faculty of Science; Niigata University; 8050 Ikarashi 2-no-cho Nishi-ku Niigata 950-2181 Japan
| | - Yoshizumi Ishino
- Department of Bioscience and Biotechnology; Kyushu University; 6-10-1 Hakozaki Higashi-ku Fukuoka 812-8581 Japan
| | - Toshio Uchiumi
- Department of Biology; Faculty of Science; Niigata University; 8050 Ikarashi 2-no-cho Nishi-ku Niigata 950-2181 Japan
| |
Collapse
|
3
|
Ito K, Honda T, Suzuki T, Miyoshi T, Murakami R, Yao M, Uchiumi T. Molecular insights into the interaction of the ribosomal stalk protein with elongation factor 1α. Nucleic Acids Res 2014; 42:14042-52. [PMID: 25428348 PMCID: PMC4267659 DOI: 10.1093/nar/gku1248] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
In all organisms, the large ribosomal subunit contains multiple copies of a flexible protein, the so-called ‘stalk’. The C-terminal domain (CTD) of the stalk interacts directly with the translational GTPase factors, and this interaction is required for factor-dependent activity on the ribosome. Here we have determined the structure of a complex of the CTD of the archaeal stalk protein aP1 and the GDP-bound archaeal elongation factor aEF1α at 2.3 Å resolution. The structure showed that the CTD of aP1 formed a long extended α-helix, which bound to a cleft between domains 1 and 3 of aEF1α, and bridged these domains. This binding between the CTD of aP1 and the aEF1α•GDP complex was formed mainly by hydrophobic interactions. The docking analysis showed that the CTD of aP1 can bind to aEF1α•GDP located on the ribosome. An additional biochemical assay demonstrated that the CTD of aP1 also bound to the aEF1α•GTP•aminoacyl-tRNA complex. These results suggest that the CTD of aP1 interacts with aEF1α at various stages in translation. Furthermore, phylogenetic perspectives and functional analyses suggested that the eukaryotic stalk protein also interacts directly with domains 1 and 3 of eEF1α, in a manner similar to the interaction of archaeal aP1 with aEF1α.
Collapse
Affiliation(s)
- Kosuke Ito
- Department of Biology, Faculty of Science, Niigata University, 8050 Ikarashi 2-no-cho, Nishi-ku, Niigata 950-2181, Japan
| | - Takayoshi Honda
- Department of Biology, Faculty of Science, Niigata University, 8050 Ikarashi 2-no-cho, Nishi-ku, Niigata 950-2181, Japan
| | - Takahiro Suzuki
- Department of Biology, Faculty of Science, Niigata University, 8050 Ikarashi 2-no-cho, Nishi-ku, Niigata 950-2181, Japan
| | - Tomohiro Miyoshi
- Department of Biology, Faculty of Science, Niigata University, 8050 Ikarashi 2-no-cho, Nishi-ku, Niigata 950-2181, Japan
| | - Ryo Murakami
- Department of Biology, Faculty of Science, Niigata University, 8050 Ikarashi 2-no-cho, Nishi-ku, Niigata 950-2181, Japan
| | - Min Yao
- Faculty of Advanced Life Science, Hokkaido University, Kita-ku, Kita-10, Nishi-8, Sapporo 060-0810, Japan
| | - Toshio Uchiumi
- Department of Biology, Faculty of Science, Niigata University, 8050 Ikarashi 2-no-cho, Nishi-ku, Niigata 950-2181, Japan
| |
Collapse
|
4
|
Simonson T, Aleksandrov A, Satpati P. Electrostatic free energies in translational GTPases: Classic allostery and the rest. Biochim Biophys Acta Gen Subj 2014; 1850:1006-1016. [PMID: 25047891 DOI: 10.1016/j.bbagen.2014.07.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2014] [Revised: 07/03/2014] [Accepted: 07/05/2014] [Indexed: 12/31/2022]
Abstract
GTPases typically switch between an inactive, OFF conformation and an active, ON conformation when a GDP ligand is replaced by GTP. Their ON/OFF populations and activity thus depend on the stabilities of four protein complexes, two apo-protein forms, and GTP/GDP in solution. A complete characterization is usually not possible experimentally and poses major challenges for simulations. We review the most important methodological challenges and we review thermodynamic data for two GTPases involved in translation of the genetic code: archaeal Initiation Factors 2 and 5B (aIF2, aIF5B). One main challenge is the multiplicity of states and conformations, including those of GTP/GDP in solution. Another is force field accuracy, especially for interactions of GTP/GDP with co-bound divalent Mg(2+) ions. The calculation of electrostatic free energies also poses specific challenges, and requires careful protocols. For aIF2, experiments and earlier simulations showed that it is a "classic" GTPase, with distinct ON/OFF conformations that prefer to bind GTP and GDP, respectively. For aIF5B, we recently proposed a non-classic mechanism, where the ON/OFF states differ only in the protonation state of Glu81 in the nucleotide binding pocket. This model is characterized here using free energy simulations. The methodological analysis should help future studies, while the aIF2, aIF5B examples illustrate the diversity of ATPase/GTPase mechanisms. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Thomas Simonson
- Laboratoire de Biochimie (CNRS unit 7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France.
| | - Alexey Aleksandrov
- Laboratoire de Biochimie (CNRS unit 7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France
| | - Priyadarshi Satpati
- Laboratoire de Biochimie (CNRS unit 7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France; Department of Cell and Molecular Biology, University of Uppsala, Sweden
| |
Collapse
|
5
|
Simonson T, Satpati P. Nucleotide recognition by the initiation factor aIF5B: free energy simulations of a neoclassical GTPase. Proteins 2012; 80:2742-57. [PMID: 22887821 DOI: 10.1002/prot.24158] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2012] [Revised: 07/16/2012] [Accepted: 07/28/2012] [Indexed: 12/13/2022]
Abstract
The GTPase aIF5B is a universally conserved initiation factor that assists ribosome assembly. Crystal structures of its nucleotide complexes, X-ray(GTP) and X-ray(GDP), are similar in the nucleotide vicinity, but differ in the orientation of a distant domain IV. This has led to two, contradictory, mechanistic models. One postulates that X-ray(GTP) and X-ray(GDP) are, respectively, the active, "ON" and the inactive, "OFF" states; the other postulates that both structures are OFF, whereas the ON state is still uncharacterized. We study GTP/GDP binding using molecular dynamics and a continuum electrostatic free energy method. We predict that X-ray(GTP) has a ≈ 3 kcal/mol preference to bind GDP, apparently contradicting its assignment as ON. However, the preference arises mainly from a single, nearby residue from the switch 2 motif: Glu81, which becomes protonated upon GTP binding, with a free energy cost of about 4 kcal/mol. We then propose a different model, where Glu81 protonation/deprotonation defines the ON/OFF states. With this model, the X-ray(GTP):GTP complex, with its protonated Glu81, is ON, whereas X-ray(GTP):GDP is OFF. The model postulates that distant conformational changes such as domain IV rotation are "uncoupled" from GTP/GDP exchange and do not affect the relative GTP/GDP binding affinities. We analyze the model using a general thermodynamic framework for GTPases. It yields rather precise predictions for the nucleotide specificities of each state, and the state specificities of each nucleotide, which are roughly comparable to the homologues IF2 and aIF2, despite the lack of any conformational switching in the model.
Collapse
Affiliation(s)
- Thomas Simonson
- Department of Biology, Laboratoire de Biochimie (CNRS UMR7654), Ecole Polytechnique, 91128 Palaiseau, France.
| | | |
Collapse
|
6
|
Blaby-Haas CE, Flood JA, Crécy-Lagard VD, Zamble DB. YeiR: a metal-binding GTPase from Escherichia coli involved in metal homeostasis. Metallomics 2012; 4:488-97. [PMID: 22511334 DOI: 10.1039/c2mt20012k] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
A comparative genomic analysis predicted that many members of the under-characterized COG0523 subfamily of putative P-loop GTPases function in metal metabolism. In this work we focused on the uncharacterized Escherichia coli protein YeiR by studying both the physiology of a yeiR mutant and the in vitro biochemical properties of YeiR expressed as a fusion with the maltose-binding protein (YeiR-MBP). Our results demonstrate that deletion of yeiR increases the sensitivity of E. coli to EDTA or cadmium, and this phenotype is linked to zinc depletion. In vitro, the tagged protein binds several Zn(2+) ions with nanomolar affinity and oligomerizes in the presence of metal. The GTPase activity of YeiR is similar to that measured for other members of the group, but GTP hydrolysis is enhanced by Zn(2+) binding. These results support the predicted connection between the COG0523 P-loop GTPases and roles in metal homeostasis.
Collapse
Affiliation(s)
- Crysten E Blaby-Haas
- Department of Microbiology & Cell Science, University of Florida, PO Box 110700, Gainesville, FL 32611-0700, USA.
| | | | | | | |
Collapse
|
7
|
Satpati P, Simonson T. Conformational selection through electrostatics: Free energy simulations of GTP and GDP binding to archaeal initiation factor 2. Proteins 2012; 80:1264-82. [PMID: 22275120 DOI: 10.1002/prot.24023] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2011] [Revised: 11/29/2011] [Accepted: 12/13/2011] [Indexed: 11/05/2022]
Abstract
Archaeal Initiation Factor 2 is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and their dependence on the conformational state, molecular dynamics free energy simulations were performed. The ON state specificity was predicted to be weak, with a GTP/GDP binding free energy difference of -1 kcal/mol, favoring GTP. The OFF state specificity is larger, 4 kcal/mol, favoring GDP. The overall effects result from a competition among many interactions in several complexes. To interpret them, we use a simpler, dielectric continuum model. Several effects are robust with respect to the model details. Both nucleotides have a net negative charge, so that removing them from solvent into the binding pocket carries a desolvation penalty, which is large for the ON state, and strongly disfavors GTP binding compared to GDP. Short-range interactions between the additional GTP phosphate group and ionized sidechains in the binding pocket offset most, but not all of the desolvation penalty; more distant groups also contribute significantly, and the switch 1 loop only slightly. The desolvation penalty is lower for the more open, wetter OFF state, and the GTP/GDP difference much smaller. Short-range interactions in the binding pocket and with more distant groups again make a significant contribution. Overall, the simulations help explain how conformational selection is achieved with a single phosphate group.
Collapse
Affiliation(s)
- Priyadarshi Satpati
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, CNRS, 91128 Palaiseau Cedex, France
| | | |
Collapse
|
8
|
Satpati P, Simonson T. Conformational selection by the aIF2 GTPase: a molecular dynamics study of functional pathways. Biochemistry 2011; 51:353-61. [PMID: 22165972 DOI: 10.1021/bi201675n] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Archaeal initiation factor 2 (aIF2) is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To improve our understanding of the role of each conformational state in the aIF2 "life cycle", we start from the state immediately after GTP hydrolysis, ON:GDP:P(i) (where P(i) is inorganic phosphate), and consider the possible next steps on the pathway to the OFF:GDP product. The first possibility is P(i) dissociation, leading to ON:GDP, which could then relax into OFF:GDP. We use molecular dynamics simulations to compute the P(i) dissociation free energy and show that dissociation is highly favorable. The second possibility is conformational relaxation into the OFF state before P(i) dissociation, to form OFF:GDP:P(i). We estimate the corresponding free energy approximately, 2 ± 3.5 kcal/mol, so that this is an uphill or weakly downhill process. A third possibility is relaxation into another conformation, neither ON nor OFF. Indeed, a third, "MIXED" conformation was seen recently in a crystal structure of the aIF2:GDP:P(i) complex. For this conformational state, P(i) dissociation is weakly unfavorable, in contrast to the ON and OFF states. From this, we will deduce that if the MIXED:GDP complex is not too unstable, the ON:GDP:P(i) → MIXED:GDP:P(i) transformation is a downhill process, which can occur spontaneously. This suggests that the MIXED state could be a functional intermediate.
Collapse
Affiliation(s)
- Priyadarshi Satpati
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, 91128 Palaiseau, France
| | | |
Collapse
|
9
|
Dmitriev SE, Stolboushkina EA, Terenin IM, Andreev DE, Garber MB, Shatsky IN. Archaeal translation initiation factor aIF2 can substitute for eukaryotic eIF2 in ribosomal scanning during mammalian 48S complex formation. J Mol Biol 2011; 413:106-14. [PMID: 21884705 DOI: 10.1016/j.jmb.2011.08.026] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2011] [Revised: 08/10/2011] [Accepted: 08/11/2011] [Indexed: 02/05/2023]
Abstract
Heterotrimeric translation initiation factor (IF) a/eIF2 (archaeal/eukaryotic IF 2) is present in both Eukarya and Archaea. Despite strong structural similarity between a/eIF2 orthologs from the two domains of life, their functional relationship is obscure. Here, we show that aIF2 from Sulfolobus solfataricus can substitute for its mammalian counterpart in the reconstitution of eukaryotic 48S initiation complexes from purified components. aIF2 is able to correctly place the initiator Met-tRNA(i) into the P-site of the 40S ribosomal subunit and accompany the entire set of eukaryotic translation IFs in the process of cap-dependent scanning and AUG codon selection. However, it seems to be unable to participate in the following step of ribosomal subunit joining. In accordance with this, aIF2 inhibits rather than stimulates protein synthesis in mammalian cell-free system. The ability of recombinant aIF2 protein to direct ribosomal scanning suggests that some archaeal mRNAs may utilize this mechanism during translation initiation.
Collapse
Affiliation(s)
- Sergey E Dmitriev
- Belozersky Institute of Physico-Chemical Biology, Moscow State University, Building A, Moscow 119992, Russia.
| | | | | | | | | | | |
Collapse
|
10
|
English BP, Hauryliuk V, Sanamrad A, Tankov S, Dekker NH, Elf J. Single-molecule investigations of the stringent response machinery in living bacterial cells. Proc Natl Acad Sci U S A 2011; 108:E365-73. [PMID: 21730169 PMCID: PMC3150888 DOI: 10.1073/pnas.1102255108] [Citation(s) in RCA: 203] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
The RelA-mediated stringent response is at the heart of bacterial adaptation to starvation and stress, playing a major role in the bacterial cell cycle and virulence. RelA integrates several environmental cues and synthesizes the alarmone ppGpp, which globally reprograms transcription, translation, and replication. We have developed and implemented novel single-molecule tracking methodology to characterize the intracellular catalytic cycle of RelA. Our single-molecule experiments show that RelA is on the ribosome under nonstarved conditions and that the individual enzyme molecule stays off the ribosome for an extended period of time after activation. This suggests that the catalytically active part of the RelA cycle is performed off, rather than on, the ribosome, and that rebinding to the ribosome is not necessary to trigger each ppGpp synthesis event. Furthermore, we find fast activation of RelA in response to heat stress followed by RelA rapidly being reset to its inactive state, which makes the system sensitive to new environmental cues and hints at an underlying excitable response mechanism.
Collapse
Affiliation(s)
- Brian P. English
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Vasili Hauryliuk
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
- University of Tartu, Institute of Technology, Tartu, Estonia
| | - Arash Sanamrad
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Stoyan Tankov
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
- University of Tartu, Institute of Technology, Tartu, Estonia
| | - Nynke H. Dekker
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
- Department of Bionanoscience, Kavli Institute of Nanoscience, Delft University of Technology, Delft, The Netherlands; and
| | - Johan Elf
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
- Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| |
Collapse
|
11
|
Satpati P, Clavaguéra C, Ohanessian G, Simonson T. Free energy simulations of a GTPase: GTP and GDP binding to archaeal initiation factor 2. J Phys Chem B 2011; 115:6749-63. [PMID: 21534562 PMCID: PMC3097523 DOI: 10.1021/jp201934p] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
![]()
Archaeal initiation factor 2 (aIF2) is a protein involved in the initiation of protein biosynthesis. In its GTP-bound, “ON” conformation, aIF2 binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, “OFF” conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and its dependence on the ON or OFF conformational state of aIF2, molecular dynamics free energy simulations (MDFE) are a tool of choice. However, the validity of the computed free energies depends on the simulation model, including the force field and the boundary conditions, and on the extent of conformational sampling in the simulations. aIF2 and other GTPases present specific difficulties; in particular, the nucleotide ligand coordinates a divalent Mg2+ ion, which can polarize the electronic distribution of its environment. Thus, a force field with an explicit treatment of electronic polarizability could be necessary, rather than a simpler, fixed charge force field. Here, we begin by comparing a fixed charge force field to quantum chemical calculations and experiment for Mg2+:phosphate binding in solution, with the force field giving large errors. Next, we consider GTP and GDP bound to aIF2 and we compare two fixed charge force fields to the recent, polarizable, AMOEBA force field, extended here in a simple, approximate manner to include GTP. We focus on a quantity that approximates the free energy to change GTP into GDP. Despite the errors seen for Mg2+:phosphate binding in solution, we observe a substantial cancellation of errors when we compare the free energy change in the protein to that in solution, or when we compare the protein ON and OFF states. Finally, we have used the fixed charge force field to perform MDFE simulations and alchemically transform GTP into GDP in the protein and in solution. With a total of about 200 ns of molecular dynamics, we obtain good convergence and a reasonable statistical uncertainty, comparable to the force field uncertainty, and somewhat lower than the predicted GTP/GDP binding free energy differences. The sign and magnitudes of the differences can thus be interpreted at a semiquantitative level, and are found to be consistent with the experimental binding preferences of ON- and OFF-aIF2.
Collapse
Affiliation(s)
- Priyadarshi Satpati
- Laboratoire de Biochimie (CNRS UMR7654), Department of Biology, Ecole Polytechnique, CNRS, Palaiseau, France
| | | | | | | |
Collapse
|