1
|
Natarajan L, De Sciscio ML, Nardi AN, Sekhar A, Del Giudice A, D’Abramo M, Naganathan AN. A finely balanced order-disorder equilibrium sculpts the folding-binding landscape of an antibiotic sequestering protein. Proc Natl Acad Sci U S A 2024; 121:e2318855121. [PMID: 38709926 PMCID: PMC11098121 DOI: 10.1073/pnas.2318855121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Accepted: 03/24/2024] [Indexed: 05/08/2024] Open
Abstract
TipA, a MerR family transcription factor from Streptomyces lividans, promotes antibiotic resistance by sequestering broad-spectrum thiopeptide-based antibiotics, thus counteracting their inhibitory effect on ribosomes. TipAS, a minimal binding motif which is expressed as an isoform of TipA, harbors a partially disordered N-terminal subdomain that folds upon binding multiple antibiotics. The extent and nature of the underlying molecular heterogeneity in TipAS that shapes its promiscuous folding-function landscape is an open question and is critical for understanding antibiotic-sequestration mechanisms. Here, combining equilibrium and time-resolved experiments, statistical modeling, and simulations, we show that the TipAS native ensemble exhibits a pre-equilibrium between binding-incompetent and binding-competent substates, with the fully folded state appearing only as an excited state under physiological conditions. The binding-competent state characterized by a partially structured N-terminal subdomain loses structure progressively in the physiological range of temperatures, swells on temperature increase, and displays slow conformational exchange across multiple conformations. Binding to the bactericidal antibiotic thiostrepton follows a combination of induced-fit and conformational-selection-like mechanisms, via partial binding and concomitant stabilization of the binding-competent substate. These ensemble features are evolutionarily conserved across orthologs from select bacteria that infect humans, underscoring the functional role of partial disorder in the native ensemble of antibiotic-sequestering proteins belonging to the MerR family.
Collapse
Affiliation(s)
- Lawanya Natarajan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai600036, India
| | | | | | - Ashok Sekhar
- Molecular Biophysics Unit, Indian Institute of Science Bangalore, Bengaluru560 012, India
| | | | - Marco D’Abramo
- Department of Chemistry, Sapienza University of Rome, Rome00185, Italy
| | - Athi N. Naganathan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai600036, India
| |
Collapse
|
2
|
Lukose B, Maruno T, Faidh M, Uchiyama S, Naganathan A. Molecular and thermodynamic determinants of self-assembly and hetero-oligomerization in the enterobacterial thermo-osmo-regulatory protein H-NS. Nucleic Acids Res 2024; 52:2157-2173. [PMID: 38340344 PMCID: PMC10954469 DOI: 10.1093/nar/gkae090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 01/03/2024] [Accepted: 01/29/2024] [Indexed: 02/12/2024] Open
Abstract
Environmentally regulated gene expression is critical for bacterial survival under stress conditions, including extremes in temperature, osmolarity and nutrient availability. Here, we dissect the thermo- and osmo-responsory behavior of the transcriptional repressor H-NS, an archetypal nucleoid-condensing sensory protein, ubiquitous in enterobacteria that infect the mammalian gut. Through experiments and thermodynamic modeling, we show that H-NS exhibits osmolarity, temperature and concentration dependent self-association, with a highly polydisperse native ensemble dominated by monomers, dimers, tetramers and octamers. The relative population of these oligomeric states is determined by an interplay between dimerization and higher-order oligomerization, which in turn drives a competition between weak homo- versus hetero-oligomerization of protein-protein and protein-DNA complexes. A phosphomimetic mutation, Y61E, fully eliminates higher-order self-assembly and preserves only dimerization while weakening DNA binding, highlighting that oligomerization is a prerequisite for strong DNA binding. We further demonstrate the presence of long-distance thermodynamic connectivity between dimerization and oligomerization sites on H-NS which influences the binding of the co-repressor Cnu, and switches the DNA binding mode of the hetero-oligomeric H-NS:Cnu complex. Our work thus uncovers important organizational principles in H-NS including a multi-layered thermodynamic control, and provides a molecular framework broadly applicable to other thermo-osmo sensory proteins that employ similar mechanisms to regulate gene expression.
Collapse
Affiliation(s)
- Bincy Lukose
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | | | - Mohammed A Faidh
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | | | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
3
|
Kannan A, Chaurasiya DK, Naganathan AN. Conflicting Interfacial Electrostatic Interactions as a Design Principle to Modulate Long-Range Interdomain Communication. ACS BIO & MED CHEM AU 2024; 4:53-67. [PMID: 38404745 PMCID: PMC10885104 DOI: 10.1021/acsbiomedchemau.3c00047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 10/19/2023] [Accepted: 10/23/2023] [Indexed: 02/27/2024]
Abstract
The extent and molecular basis of interdomain communication in multidomain proteins, central to understanding allostery and function, is an open question. One simple evolutionary strategy could involve the selection of either conflicting or favorable electrostatic interactions across the interface of two closely spaced domains to tune the magnitude of interdomain connectivity. Here, we study a bilobed domain FF34 from the eukaryotic p190A RhoGAP protein to explore one such design principle that mediates interdomain communication. We find that while the individual structural units in wild-type FF34 are marginally coupled, they exhibit distinct intrinsic stabilities and low cooperativity, manifesting as slow folding. The FF3-FF4 interface harbors a frustrated network of highly conserved electrostatic interactions-a charge troika-that promotes the population of multiple, decoupled, and non-native structural modes on a rugged native landscape. Perturbing this network via a charge-reversal mutation not only enhances stability and cooperativity but also dampens the fluctuations globally and speeds up the folding rate by at least an order of magnitude. Our work highlights how a conserved but nonoptimal network of interfacial electrostatic interactions shapes the native ensemble of a bilobed protein, a feature that could be exploited in designing molecular systems with long-range connectivity and enhanced cooperativity.
Collapse
Affiliation(s)
- Adithi Kannan
- Department of Biotechnology,
Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Dhruv Kumar Chaurasiya
- Department of Biotechnology,
Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Athi N. Naganathan
- Department of Biotechnology,
Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
4
|
Dani R, Pawloski W, Chaurasiya DK, Srilatha NS, Agarwal S, Fushman D, Naganathan AN. Conformational Tuning Shapes the Balance between Functional Promiscuity and Specialization in Paralogous Plasmodium Acyl-CoA Binding Proteins. Biochemistry 2023; 62:2982-2996. [PMID: 37788430 PMCID: PMC10774088 DOI: 10.1021/acs.biochem.3c00449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
Paralogous proteins confer enhanced fitness to organisms via complex sequence-conformation codes that shape functional divergence, specialization, or promiscuity. Here, we dissect the underlying mechanism of promiscuous binding versus partial subfunctionalization in paralogues by studying structurally identical acyl-CoA binding proteins (ACBPs) from Plasmodium falciparum that serve as promising drug targets due to their high expression during the protozoan proliferative phase. Combining spectroscopic measurements, solution NMR, SPR, and simulations on two of the paralogues, A16 and A749, we show that minor sequence differences shape nearly every local and global conformational feature. A749 displays a broader and heterogeneous native ensemble, weaker thermodynamic coupling and cooperativity, enhanced fluctuations, and a larger binding pocket volume compared to A16. Site-specific tryptophan probes signal a graded reduction in the sampling of substates in the holo form, which is particularly apparent in A749. The paralogues exhibit a spectrum of binding affinities to different acyl-CoAs with A749, the more promiscuous and hence the likely ancestor, binding 1000-fold stronger to lauroyl-CoA under physiological conditions. We thus demonstrate how minor sequence changes modulate the extent of long-range interactions and dynamics, effectively contributing to the molecular evolution of contrasting functional repertoires in paralogues.
Collapse
Affiliation(s)
- Rahul Dani
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Westley Pawloski
- Center for Biomolecular Structure & Organization, Department of Chemistry & Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Dhruv Kumar Chaurasiya
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | | | - Sonal Agarwal
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - David Fushman
- Center for Biomolecular Structure & Organization, Department of Chemistry & Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
5
|
Thermodynamic architecture and conformational plasticity of GPCRs. Nat Commun 2023; 14:128. [PMID: 36624096 PMCID: PMC9829892 DOI: 10.1038/s41467-023-35790-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 12/29/2022] [Indexed: 01/11/2023] Open
Abstract
G-protein-coupled receptors (GPCRs) are ubiquitous integral membrane proteins involved in diverse cellular signaling processes. Here, we carry out a large-scale ensemble thermodynamic study of 45 ligand-free GPCRs employing a structure-based statistical mechanical framework. We find that multiple partially structured states co-exist in the GPCR native ensemble, with the TM helices 1, 6 and 7 displaying varied folding status, and shaping the conformational landscape. Strongly coupled residues are anisotropically distributed, accounting for only 13% of the residues, illustrating that a large number of residues are inherently dynamic. Active-state GPCRs are characterized by reduced conformational heterogeneity with altered coupling-patterns distributed throughout the structural scaffold. In silico alanine-scanning mutagenesis reveals that extra- and intra-cellular faces of GPCRs are coupled thermodynamically, highlighting an exquisite structural specialization and the fluid nature of the intramolecular interaction network. The ensemble-based perturbation methodology presented here lays the foundation for understanding allosteric mechanisms and the effects of disease-causing mutations in GCPRs.
Collapse
|
6
|
Loss of stability and unfolding cooperativity in hPGK1 upon gradual structural perturbation of its N-terminal domain hydrophobic core. Sci Rep 2022; 12:17200. [PMID: 36229482 PMCID: PMC9561527 DOI: 10.1038/s41598-022-22088-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Accepted: 10/10/2022] [Indexed: 01/06/2023] Open
Abstract
Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype-phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.
Collapse
|
7
|
Ooka K, Liu R, Arai M. The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics. Molecules 2022; 27:molecules27144460. [PMID: 35889332 PMCID: PMC9319528 DOI: 10.3390/molecules27144460] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Revised: 07/07/2022] [Accepted: 07/08/2022] [Indexed: 11/16/2022] Open
Abstract
Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.
Collapse
Affiliation(s)
- Koji Ooka
- Department of Physics, Graduate School of Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
- Komaba Organization for Educational Excellence, College of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
| | - Runjing Liu
- Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
| | - Munehito Arai
- Department of Physics, Graduate School of Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
- Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan;
- Correspondence:
| |
Collapse
|
8
|
Pacheco-Garcia JL, Loginov DS, Anoz-Carbonell E, Vankova P, Palomino-Morales R, Salido E, Man P, Medina M, Naganathan AN, Pey AL. Allosteric Communication in the Multifunctional and Redox NQO1 Protein Studied by Cavity-Making Mutations. Antioxidants (Basel) 2022; 11:antiox11061110. [PMID: 35740007 PMCID: PMC9219786 DOI: 10.3390/antiox11061110] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2022] [Revised: 05/28/2022] [Accepted: 05/30/2022] [Indexed: 02/01/2023] Open
Abstract
Allosterism is a common phenomenon in protein biochemistry that allows rapid regulation of protein stability; dynamics and function. However, the mechanisms by which allosterism occurs (by mutations or post-translational modifications (PTMs)) may be complex, particularly due to long-range propagation of the perturbation across protein structures. In this work, we have investigated allosteric communication in the multifunctional, cancer-related and antioxidant protein NQO1 by mutating several fully buried leucine residues (L7, L10 and L30) to smaller residues (V, A and G) at sites in the N-terminal domain. In almost all cases, mutated residues were not close to the FAD or the active site. Mutations L→G strongly compromised conformational stability and solubility, and L30A and L30V also notably decreased solubility. The mutation L10A, closer to the FAD binding site, severely decreased FAD binding affinity (≈20 fold vs. WT) through long-range and context-dependent effects. Using a combination of experimental and computational analyses, we show that most of the effects are found in the apo state of the protein, in contrast to other common polymorphisms and PTMs previously characterized in NQO1. The integrated study presented here is a first step towards a detailed structural–functional mapping of the mutational landscape of NQO1, a multifunctional and redox signaling protein of high biomedical relevance.
Collapse
Affiliation(s)
- Juan Luis Pacheco-Garcia
- Departamento de Química Física, Universidad de Granada, Av. Fuentenueva s/n, 18071 Granada, Spain
- Correspondence: (J.L.P.-G.); (A.L.P.); Tel.: +34-958243173 (A.L.P.)
| | - Dmitry S. Loginov
- Institute of Microbiology—BioCeV, Academy of Sciences of the Czech Republic, Prumyslova 595, 252 50 Vestec, Czech Republic; (D.S.L.); (P.M.)
| | - Ernesto Anoz-Carbonell
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI) (GBsC-CSIC Joint Unit), Universidad de Zaragoza, 50009 Zaragoza, Spain; (E.A.-C.); (M.M.)
| | - Pavla Vankova
- Institute of Biotechnology—BioCeV, Academy of Sciences of the Czech Republic, Prumyslova 595, 252 50 Vestec, Czech Republic;
- Department of Biochemistry, Faculty of Science, Charles University, Hlavova 2030/8, 128 43 Prague, Czech Republic
| | - Rogelio Palomino-Morales
- Departamento de Bioquímica y Biología Molecular I, Facultad de Ciencias y Centro de Investigaciones Biomédicas (CIBM), Universidad de Granada, 18016 Granada, Spain;
| | - Eduardo Salido
- Center for Rare Diseases (CIBERER), Hospital Universitario de Canarias, Universidad de la Laguna, 38320 Tenerife, Spain;
| | - Petr Man
- Institute of Microbiology—BioCeV, Academy of Sciences of the Czech Republic, Prumyslova 595, 252 50 Vestec, Czech Republic; (D.S.L.); (P.M.)
| | - Milagros Medina
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI) (GBsC-CSIC Joint Unit), Universidad de Zaragoza, 50009 Zaragoza, Spain; (E.A.-C.); (M.M.)
| | - Athi N. Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras (IITM), Chennai 600036, India;
| | - Angel L. Pey
- Departamento de Química Física, Unidad de Excelencia en Química Aplicada a Biomedicina y Medioambiente e Instituto de Biotecnología, Universidad de Granada, Av. Fuentenueva s/n, 18071 Granada, Spain
- Correspondence: (J.L.P.-G.); (A.L.P.); Tel.: +34-958243173 (A.L.P.)
| |
Collapse
|
9
|
Naganathan AN, Dani R, Gopi S, Aranganathan A, Narayan A. Folding Intermediates, Heterogeneous Native Ensembles and Protein Function. J Mol Biol 2021; 433:167325. [PMID: 34695380 DOI: 10.1016/j.jmb.2021.167325] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Revised: 10/17/2021] [Accepted: 10/18/2021] [Indexed: 01/01/2023]
Abstract
Single domain proteins fold via diverse mechanisms emphasizing the intricate relationship between energetics and structure, which is a direct consequence of functional constraints and demands imposed at the level of sequence. On the other hand, elucidating the interplay between folding mechanisms and function is challenging in large proteins, given the inherent shortcomings in identifying metastable states experimentally and the sampling limitations associated with computational methods. Here, we show that free energy profiles and surfaces of large systems (>150 residues), as predicted by a statistical mechanical model, display a wide array of folding mechanisms with ubiquitous folding intermediates and heterogeneous native ensembles. Importantly, residues around the ligand binding or enzyme active site display a larger tendency to partially unfold and this manifests as intermediates or excited states along the folding coordinate in ligand binding domains, transcription repressors, and representative enzymes from all the six classes, including the SARS-CoV-2 receptor binding domain (RBD) of the spike protein and the protease Mpro. It thus appears that it is relatively easier to distill the imprints of function on the folding landscape of larger proteins as opposed to smaller systems. We discuss how an understanding of energetic-entropic features in ordered proteins can pinpoint specific avenues through which folding mechanisms, populations of partially structured states and function can be engineered.
Collapse
Affiliation(s)
- Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India.
| | - Rahul Dani
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India. https://twitter.com/Soundha
| | - Akashnathan Aranganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Abhishek Narayan
- Department of Physics, University of Alberta, Edmonton, AB T6G 2E1, Canada
| |
Collapse
|
10
|
Gamiz-Arco G, Risso VA, Gaucher EA, Gavira JA, Naganathan AN, Ibarra-Molero B, Sanchez-Ruiz JM. Combining Ancestral Reconstruction with Folding-Landscape Simulations to Engineer Heterologous Protein Expression. J Mol Biol 2021; 433:167321. [PMID: 34687715 DOI: 10.1016/j.jmb.2021.167321] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 10/01/2021] [Accepted: 10/17/2021] [Indexed: 11/30/2022]
Abstract
Obligate symbionts typically exhibit high evolutionary rates. Consequently, their proteins may differ considerably from their modern and ancestral homologs in terms of both sequence and properties, thus providing excellent models to study protein evolution. Also, obligate symbionts are challenging to culture in the lab and proteins from uncultured organisms must be produced in heterologous hosts using recombinant DNA technology. Obligate symbionts thus replicate a fundamental scenario of metagenomics studies aimed at the functional characterization and biotechnological exploitation of proteins from the bacteria in soil. Here, we use the thioredoxin from Candidatus Photodesmus katoptron, an uncultured symbiont of flashlight fish, to explore evolutionary and engineering aspects of protein folding in heterologous hosts. The symbiont protein is a standard thioredoxin in terms of 3D-structure, stability and redox activity. However, its folding outside the original host is severely impaired, as shown by a very slow refolding in vitro and an inefficient expression in E. coli that leads mostly to insoluble protein. By contrast, resurrected Precambrian thioredoxins express efficiently in E. coli, plausibly reflecting an ancient adaptation to unassisted folding. We have used a statistical-mechanical model of the folding landscape to guide back-to-ancestor engineering of the symbiont protein. Remarkably, we find that the efficiency of heterologous expression correlates with the in vitro (i.e., unassisted) folding rate and that the ancestral expression efficiency can be achieved with only 1-2 back-to-ancestor replacements. These results demonstrate a minimal-perturbation, sequence-engineering approach to rescue inefficient heterologous expression which may potentially be useful in metagenomics efforts targeting recent adaptations.
Collapse
Affiliation(s)
- Gloria Gamiz-Arco
- Departamento de Quimica Fisica, Facultad de Ciencias, Unidad de Excelencia de Quimica Aplicada a Biomedicina y Medioambiente (UEQ), Universidad de Granada, 18071 Granada, Spain
| | - Valeria A Risso
- Departamento de Quimica Fisica, Facultad de Ciencias, Unidad de Excelencia de Quimica Aplicada a Biomedicina y Medioambiente (UEQ), Universidad de Granada, 18071 Granada, Spain
| | - Eric A Gaucher
- Department of Biology, Georgia State University, Atlanta, GA 30303, USA
| | - Jose A Gavira
- Laboratorio de Estudios Cristalograficos, Instituto Andaluz de Ciencias de la Tierra, CSIC, Unidad de Excelencia de Quimica Aplicada a Biomedicina y Medioambiente (UEQ), Universidad de Granada, Avenida de las Palmeras 4, Armilla, Granada 18100, Spain. https://twitter.com/Gavirius
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India.
| | - Beatriz Ibarra-Molero
- Departamento de Quimica Fisica, Facultad de Ciencias, Unidad de Excelencia de Quimica Aplicada a Biomedicina y Medioambiente (UEQ), Universidad de Granada, 18071 Granada, Spain.
| | - Jose M Sanchez-Ruiz
- Departamento de Quimica Fisica, Facultad de Ciencias, Unidad de Excelencia de Quimica Aplicada a Biomedicina y Medioambiente (UEQ), Universidad de Granada, 18071 Granada, Spain.
| |
Collapse
|
11
|
Naganathan AN, Kannan A. A hierarchy of coupling free energies underlie the thermodynamic and functional architecture of protein structures. Curr Res Struct Biol 2021; 3:257-267. [PMID: 34704074 PMCID: PMC8526763 DOI: 10.1016/j.crstbi.2021.09.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Revised: 09/08/2021] [Accepted: 09/30/2021] [Indexed: 12/22/2022] Open
Abstract
Protein sequences and structures evolve by satisfying varied physical and biochemical constraints. This multi-level selection is enabled not just by the patterning of amino acids on the sequence, but also via coupling between residues in the native structure. Here, we employ an energetically detailed statistical mechanical model with millions of microstates to extract such long-range structural correlations, i.e. thermodynamic coupling free energies, from a diverse family of protein structures. We find that despite the intricate and anisotropic distribution of coupling patterns, the majority of residues (>70%) are only marginally coupled contributing to functional motions and catalysis. Physical origins of ‘sectors’, determinants of native ensemble heterogeneity in extant, ancient and designed proteins, and the basis for allostery emerge naturally from coupling free energies. The statistical framework highlights how evolutionary selection and optimization occur at the level of global interaction network for a given protein fold impacting folding, function, and allosteric outputs. Evolution of protein structures occurs at the level of global interaction network. More than 70% of the protein residues are weakly or marginally coupled. Functional ‘sector’ regions are a manifestation of marginal coupling. Coupling indices vary across the entire proteins in extant-ancient and natural-designed pairs. The proposed methodology can be used to understand allostery and epistasis.
Collapse
Affiliation(s)
- Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, 600036, India
| | - Adithi Kannan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, 600036, India
| |
Collapse
|
12
|
Pacheco-Garcia JL, Anoz-Carbonell E, Vankova P, Kannan A, Palomino-Morales R, Mesa-Torres N, Salido E, Man P, Medina M, Naganathan AN, Pey AL. Structural basis of the pleiotropic and specific phenotypic consequences of missense mutations in the multifunctional NAD(P)H:quinone oxidoreductase 1 and their pharmacological rescue. Redox Biol 2021; 46:102112. [PMID: 34537677 PMCID: PMC8455868 DOI: 10.1016/j.redox.2021.102112] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Revised: 07/21/2021] [Accepted: 08/17/2021] [Indexed: 10/31/2022] Open
Abstract
The multifunctional nature of human flavoproteins is critically linked to their ability to populate multiple conformational states. Ligand binding, post-translational modifications and disease-associated mutations can reshape this functional landscape, although the structure-function relationships of these effects are not well understood. Herein, we characterized the structural and functional consequences of two mutations (the cancer-associated P187S and the phosphomimetic S82D) on different ligation states which are relevant to flavin binding, intracellular stability and catalysis of the disease-associated NQO1 flavoprotein. We found that these mutations affected the stability locally and their effects propagated differently through the protein structure depending both on the nature of the mutation and the ligand bound, showing directional preference from the mutated site and leading to specific phenotypic manifestations in different functional traits (FAD binding, catalysis and inhibition, intracellular stability and pharmacological response to ligands). Our study thus supports that pleitropic effects of disease-causing mutations and phosphorylation events on human flavoproteins may be caused by long-range structural propagation of stability effects to different functional sites that depend on the ligation-state and site-specific perturbations. Our approach can be of general application to investigate these pleiotropic effects at the flavoproteome scale in the absence of high-resolution structural models.
Collapse
Affiliation(s)
- Juan Luis Pacheco-Garcia
- Departamento de Química Física, Universidad de Granada, Av. Fuentenueva s/n, 18071, Granada, Spain
| | - Ernesto Anoz-Carbonell
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Instituto de Biocomputación y Física de Sistemas Complejos (GBsC-CSIC and Joint Unit), Universidad de Zaragoza, 50009, Zaragoza, Spain
| | - Pavla Vankova
- Institute of Microbiology, Academy of Sciences of the Czech Republic, Videnska 1083, Prague 4, 142 20, Czech Republic; Department of Biochemistry, Faculty of Science, Charles University, Hlavova 2030/8, Prague 2, 128 43, Czech Republic
| | - Adithi Kannan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras (IITM), Chennai, 600036, India
| | - Rogelio Palomino-Morales
- Departmento de Bioquímica y Biología Molecular I, Facultad de Ciencias y Centro de Investigaciones Biomédicas (CIBM), Universidad de Granada, Granada, Spain
| | - Noel Mesa-Torres
- Departamento de Química Física, Universidad de Granada, Av. Fuentenueva s/n, 18071, Granada, Spain
| | - Eduardo Salido
- Center for Rare Diseases (CIBERER), Hospital Universitario de Canarias, Universidad de la Laguna, 38320, Tenerife, Spain
| | - Petr Man
- Institute of Microbiology, Academy of Sciences of the Czech Republic, Videnska 1083, Prague 4, 142 20, Czech Republic
| | - Milagros Medina
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Instituto de Biocomputación y Física de Sistemas Complejos (GBsC-CSIC and Joint Unit), Universidad de Zaragoza, 50009, Zaragoza, Spain
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras (IITM), Chennai, 600036, India
| | - Angel L Pey
- Departamento de Química Física, Unidad de Excelencia en Química Aplicada a Biomedicina y Medioambiente e Instituto de Biotecnología, Universidad de Granada, Av. Fuentenueva s/n, 18071, Granada, Spain.
| |
Collapse
|
13
|
Sannigrahi A, Chowdhury S, Das B, Banerjee A, Halder A, Kumar A, Saleem M, Naganathan AN, Karmakar S, Chattopadhyay K. The metal cofactor zinc and interacting membranes modulate SOD1 conformation-aggregation landscape in an in vitro ALS model. eLife 2021; 10:e61453. [PMID: 33825682 PMCID: PMC8087447 DOI: 10.7554/elife.61453] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2020] [Accepted: 04/01/2021] [Indexed: 12/22/2022] Open
Abstract
Aggregation of Cu-Zn superoxide dismutase (SOD1) is implicated in the motor neuron disease, amyotrophic lateral sclerosis (ALS). Although more than 140 disease mutations of SOD1 are available, their stability or aggregation behaviors in membrane environment are not correlated with disease pathophysiology. Here, we use multiple mutational variants of SOD1 to show that the absence of Zn, and not Cu, significantly impacts membrane attachment of SOD1 through two loop regions facilitating aggregation driven by lipid-induced conformational changes. These loop regions influence both the primary (through Cu intake) and the gain of function (through aggregation) of SOD1 presumably through a shared conformational landscape. Combining experimental and theoretical frameworks using representative ALS disease mutants, we develop a 'co-factor derived membrane association model' wherein mutational stress closer to the Zn (but not to the Cu) pocket is responsible for membrane association-mediated toxic aggregation and survival time scale after ALS diagnosis.
Collapse
Affiliation(s)
- Achinta Sannigrahi
- Structural Biology & Bio-Informatics Division, CSIR-Indian Institute of Chemical BiologyKolkataIndia
| | - Sourav Chowdhury
- Chemistry and Chemical Biology, Harvard UniversityCambridgeUnited States
| | - Bidisha Das
- Structural Biology & Bio-Informatics Division, CSIR-Indian Institute of Chemical BiologyKolkataIndia
- Academy of Scientific and Innovative Research (AcSIR), CSIR- Human Resource development Centre CampusGhaziabadIndia
| | | | | | - Amaresh Kumar
- School of Biological Sciences, National Institute of Science Education and Research (NISER)BhubaneswarIndia
| | - Mohammed Saleem
- School of Biological Sciences, National Institute of Science Education and Research (NISER)BhubaneswarIndia
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology MadrasChennaiIndia
| | | | - Krishnananda Chattopadhyay
- Structural Biology & Bio-Informatics Division, CSIR-Indian Institute of Chemical BiologyKolkataIndia
- Academy of Scientific and Innovative Research (AcSIR), CSIR- Human Resource development Centre CampusGhaziabadIndia
| |
Collapse
|
14
|
Gopi S, Lukose B, Naganathan AN. Diverse Native Ensembles Dictate the Differential Functional Responses of Nuclear Receptor Ligand-Binding Domains. J Phys Chem B 2021; 125:3546-3555. [PMID: 33818099 DOI: 10.1021/acs.jpcb.1c00972] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Native states of folded proteins are characterized by a large ensemble of conformations whose relative populations and interconversion dynamics determine the functional output. This is more apparent in transcription factors that have evolved to be inherently sensitive to small perturbations, thus fine-tuning gene expression. To explore the extent to which such functional features are imprinted on the folding landscape of transcription factor ligand-binding domains (LBDs), we characterize paralogous LBDs of the nuclear receptor (NR) family employing an energetically detailed and ensemble-based Ising-like statistical mechanical model. We find that the native ensembles of the LBDs from glucocorticoid receptor, PPAγ, and thyroid hormone receptor display a remarkable diversity in the width of the native wells, the number and nature of partially structured states, and hence the degree of conformational order. Monte Carlo simulations employing the full state representation of the ensemble highlight that many of the functional conformations coexist in equilibrium, whose relative populations are sensitive to both temperature and the strength of ligand binding. Allosteric modulation of the degree of structure at a coregulator binding site on ligand binding is shown to arise via a redistribution of populations in the native ensembles of glucocorticoid and PPAγ LBDs. Our results illustrate how functional requirements can drive the evolution of conformationally diverse native ensembles in paralogs.
Collapse
Affiliation(s)
- Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Bincy Lukose
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
15
|
Garg S, Sagar A, Singaraju GS, Dani R, Bari NK, Naganathan AN, Rakshit S. Weakening of interaction networks with aging in tip-link protein induces hearing loss. Biochem J 2021; 478:121-134. [PMID: 33270084 PMCID: PMC7813477 DOI: 10.1042/bcj20200799] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Revised: 12/02/2020] [Accepted: 12/03/2020] [Indexed: 02/08/2023]
Abstract
Age-related hearing loss (ARHL) is a common condition in humans marking the gradual decrease in hearing with age. Perturbations in the tip-link protein cadherin-23 that absorbs the mechanical tension from sound and maintains the integrity of hearing is associated with ARHL. Here, in search of molecular origins for ARHL, we dissect the conformational behavior of cadherin-23 along with the mutant S47P that progresses the hearing loss drastically. Using an array of experimental and computational approaches, we highlight a lower thermodynamic stability, significant weakening in the hydrogen-bond network and inter-residue correlations among β-strands, due to the S47P mutation. The loss in correlated motions translates to not only a remarkable two orders of magnitude slower folding in the mutant but also to a proportionately complex unfolding mechanism. We thus propose that loss in correlated motions within cadherin-23 with aging may trigger ARHL, a molecular feature that likely holds true for other disease-mutations in β-strand-rich proteins.
Collapse
Affiliation(s)
- Surbhi Garg
- Department of Chemical Sciences, Indian Institute of Science Education and Research Mohali, Punjab, India
| | - Amin Sagar
- Centre de Biochimie Structurale INSERM, CNRS, Université de Montpellier, Montpellier, France
| | - Gayathri S. Singaraju
- Department of Chemical Sciences, Indian Institute of Science Education and Research Mohali, Punjab, India
| | - Rahul Dani
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Naimat K. Bari
- Institute of Nano Science and Technology (INST), Phase-10, Sector-64, Mohali, Punjab 160062, India
| | - Athi N. Naganathan
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sabyasachi Rakshit
- Department of Chemical Sciences, Indian Institute of Science Education and Research Mohali, Punjab, India
- Centre for Protein Science Design and Engineering, Indian Institute of Science Education and Research Mohali, Punjab, India
| |
Collapse
|
16
|
Narayan A, Gopi S, Lukose B, Naganathan AN. Electrostatic Frustration Shapes Folding Mechanistic Differences in Paralogous Bacterial Stress Response Proteins. J Mol Biol 2020; 432:4830-4839. [DOI: 10.1016/j.jmb.2020.06.026] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 06/30/2020] [Accepted: 06/30/2020] [Indexed: 01/06/2023]
|