1
|
Gade M, Gardner JM, Jain P, Laurino P. Nucleoside-Driven Specificity of DNA Methyltransferase. Chembiochem 2023; 24:e202300094. [PMID: 37548117 DOI: 10.1002/cbic.202300094] [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: 02/05/2023] [Revised: 08/03/2023] [Accepted: 08/05/2023] [Indexed: 08/08/2023]
Abstract
We have studied the adenosine binding specificities of two bacterial DNA methyltransferases, Taq methyltransferase (M.TaqI), and HhaI methyltransferase (M.HhaI). While they have similar cofactor binding pocket interactions, experimental data showed different specificity for novel S-nucleobase-l-methionine cofactors (SNMs; N=guanosyl, cytidyl, uridyl). Protein dynamics corroborate the experimental data on the cofactor specificities. For M.TaqI the specificity for S-adenosyl-l-methionine (SAM) is governed by the tight binding on the nucleoside part of the cofactor, while for M.HhaI the degree of freedom of the nucleoside chain allows the acceptance of other bases. The experimental data prove catalytically productive methylation by the M.HhaI binding pocket for all the SNMs. Our results suggest a new route for successful design of unnatural SNM analogues for methyltransferases as a tool for cofactor engineering.
Collapse
Affiliation(s)
- Madhuri Gade
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, 1919-1 Tancha, Onna, Okinawa, 904-0495, Japan
| | - Jasmine M Gardner
- Department of Chemistry - BMC, Uppsala University, Box 576, 751 23, Uppsala, Sweden
| | - Prashant Jain
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, 1919-1 Tancha, Onna, Okinawa, 904-0495, Japan
| | - Paola Laurino
- Protein Engineering and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, 1919-1 Tancha, Onna, Okinawa, 904-0495, Japan
| |
Collapse
|
2
|
Dasgupta S, Dev A, Chongdar N, Basak P, Dastidar SG, Basu G. Signatures of tRNA Glx -specificity in proteobacterial glutamyl-tRNA synthetases. Proteins 2023. [PMID: 37953434 DOI: 10.1002/prot.26634] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Revised: 10/16/2023] [Accepted: 10/25/2023] [Indexed: 11/14/2023]
Abstract
The canonical function of glutamyl-tRNA synthetase (GluRS) is to glutamylate tRNAGlu . Yet not all bacterial GluRSs glutamylate tRNAGlu ; many glutamylate both tRNAGlu and tRNAGln , while some glutamylate only tRNAGln and not the cognate substrate tRNAGlu . Understanding the basis of the unique specificity of tRNAGlx is important. Mutational studies have hinted at hotspot residues, both on tRNAGlx and GluRS, which play crucial roles in tRNAGlx -specificity. However, its underlying structural basis remains unexplored. The majority of biochemical studies related to tRNAGlx -specificity have been performed on GluRS from Escherichia coli and other proteobacterial species. However, since the early crystal structures of GluRS and tRNAGlu -bound GluRS were from non-proteobacterial species (Thermus thermophilus), proteobacterial biochemical data have often been interpreted in the context of non-proteobacterial GluRS structures. Marked differences between proteobacterial and non-proteobacterial GluRSs have been demonstrated; therefore, it is important to understand tRNAGlx -specificity vis-a-vis proteobacterial GluRS structures. To this end, we solved the crystal structure of a double mutant GluRS from E. coli. Using the solved structure and several other currently available proteo- and non-proteobacterial GluRS crystal structures, we probed the structural basis of the tRNAGlx -specificity of bacterial GluRSs. Specifically, our analyses suggest a unique role played by the tRNAGlx D-helix contacting loop of GluRS in the modulation of tRNAGln -specificity. While earlier studies have identified functional hotspots on tRNAGlx that control the tRNAGlx -specificity of GluRS, this is the first report of complementary signatures of tRNAGlx -specificity in GluRS.
Collapse
Affiliation(s)
- Saumya Dasgupta
- Department of Biophysics, Bose Institute, Kolkata, India
- Department of Chemistry, Amity Institute of Applied Sciences, Amity University Kolkata, Kolkata, India
| | - Aditya Dev
- Department of Biophysics, Bose Institute, Kolkata, India
| | - Nipa Chongdar
- Department of Biophysics, Bose Institute, Kolkata, India
- Interdisciplinary School of Life Sciences, Indian Institute of Technology, Ponda, India
| | - Premananda Basak
- Department of Biological Sciences, Bose Institute, Kolkata, India
| | | | - Gautam Basu
- Department of Biophysics, Bose Institute, Kolkata, India
- Department of Biological Sciences, Indian Institute of Science Education and Research Kolkata, Nadia, India
| |
Collapse
|
3
|
Issa NT, Byers SW, Dakshanamurthy S. ES-Screen: A Novel Electrostatics-Driven Method for Drug Discovery Virtual Screening. Int J Mol Sci 2022; 23:ijms232314830. [PMID: 36499162 PMCID: PMC9736079 DOI: 10.3390/ijms232314830] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 11/23/2022] [Accepted: 11/24/2022] [Indexed: 12/02/2022] Open
Abstract
Electrostatic interactions drive biomolecular interactions and associations. Computational modeling of electrostatics in biomolecular systems, such as protein-ligand, protein-protein, and protein-DNA, has provided atomistic insights into the binding process. In drug discovery, finding biologically plausible ligand-protein target interactions is challenging as current virtual screening and adjuvant techniques such as docking methods do not provide optimal treatment of electrostatic interactions. This study describes a novel electrostatics-driven virtual screening method called 'ES-Screen' that performs well across diverse protein target systems. ES-Screen provides a unique treatment of electrostatic interaction energies independent of total electrostatic free energy, typically employed by current software. Importantly, ES-Screen uses initial ligand pose input obtained from a receptor-based pharmacophore, thus independent of molecular docking. ES-Screen integrates individual polar and nonpolar replacement energies, which are the energy costs of replacing the cognate ligand for a target with a query ligand from the screening. This uniquely optimizes thermodynamic stability in electrostatic and nonpolar interactions relative to an experimentally determined stable binding state. ES-Screen also integrates chemometrics through shape and other physicochemical properties to prioritize query ligands with the greatest physicochemical similarities to the cognate ligand. The applicability of ES-Screen is demonstrated with in vitro experiments by identifying novel targets for many drugs. The present version includes a combination of many other descriptor components that, in a future version, will be purely based on electrostatics. Therefore, ES-Screen is a first-in-class unique electrostatics-driven virtual screening method with a unique implementation of replacement electrostatic interaction energies with broad applicability in drug discovery.
Collapse
|
4
|
Yamaguchi S, Nakashima H, Moriwaki Y, Terada T, Shimizu K. Prediction of protein mononucleotide binding sites using AlphaFold2 and machine learning. Comput Biol Chem 2022; 100:107744. [DOI: 10.1016/j.compbiolchem.2022.107744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Revised: 07/12/2022] [Accepted: 07/22/2022] [Indexed: 11/26/2022]
|
5
|
Ganguly S, Murugan NA, Ghosh D, Narayanaswamy N, Govindaraju T, Basu G. DNA Minor Groove-Induced cis- trans Isomerization of a Near-Infrared Fluorescent Probe. Biochemistry 2021; 60:2084-2097. [PMID: 34142803 DOI: 10.1021/acs.biochem.1c00281] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The discovery of small molecules that exhibit turn-on far-red or near-infrared (NIR) fluorescence upon DNA binding and understanding how they bind DNA are important for imaging and bioanalytical applications. Here we report the DNA-bound structure and the DNA binding mechanism of quinone cyanine dithiazole (QCy-DT), a recently reported AT-specific turn-on NIR fluorescent probe for double-stranded DNA. The nuclear magnetic resonance (NMR)-derived structure showed minor groove binding but no specific ligand-DNA interactions, consistent with an endothermic and entropy-driven binding mechanism deduced from isothermal titration calorimetry. Minor groove binding is typically fast because it minimally perturbs the DNA structure. However, QCy-DT exhibited unusually slow DNA binding. The cyanine-based probe is capable of cis-trans isomerization due to overlapping methine bridges, with 16 possible slowly interconverting cis/trans isomers. Using NMR, density functional theory, and free energy calculations, we show that the DNA-free and DNA-bound environments of QCy-DT prefer distinctly different isomers, indicating that the origin of the slow kinetics is a cis-trans isomerization and that the minor groove preferentially selects an otherwise unstable cis/trans isomer of QCy-DT. Flux analysis showed the conformational selection pathway to be the dominating DNA binding mechanism at low DNA concentrations, which switches to the induced fit pathway at high DNA concentrations. This report of cis/trans isomerization of a ligand, upon binding the DNA minor groove, expands the prevailing understanding of unique discriminatory powers of the minor groove and has an important bearing on using polymethine cyanine dyes to probe the kinetics of molecular interactions.
Collapse
Affiliation(s)
- Sudakshina Ganguly
- Department of Biophysics, Bose Institute, P-1/12 CIT Scheme VIIM, Kolkata 700054, West Bengal, India
| | - N Arul Murugan
- Division of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology, S-106 91 Stockholm, Sweden
| | - Debasis Ghosh
- Bioorganic Chemistry Laboratory, New Chemistry Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bengaluru 560064, Karnataka, India
| | - Nagarjun Narayanaswamy
- Bioorganic Chemistry Laboratory, New Chemistry Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bengaluru 560064, Karnataka, India
| | - Thimmaiah Govindaraju
- Bioorganic Chemistry Laboratory, New Chemistry Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur P.O., Bengaluru 560064, Karnataka, India
| | - Gautam Basu
- Department of Biophysics, Bose Institute, P-1/12 CIT Scheme VIIM, Kolkata 700054, West Bengal, India
| |
Collapse
|
6
|
Vashisht K, Singh P, Verma S, Dixit R, Mishra N, Pandey KC. The nucleotide specificity of succinyl-CoA synthetase of Plasmodium falciparum is not determined by charged gatekeeper residues alone. FEBS Open Bio 2021; 11:578-587. [PMID: 33174373 PMCID: PMC7931218 DOI: 10.1002/2211-5463.13034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Revised: 10/09/2020] [Accepted: 11/09/2020] [Indexed: 02/05/2023] Open
Abstract
Substrate specificity of an enzyme is an important characteristic of its mechanism of action. Investigation of the nucleotide specificity of Plasmodium falciparum succinyl‐CoA synthetase (SCS; PfSCS) would provide crucial insights of its substrate recognition. Charged gatekeeper residues have been shown to alter the substrate specificity via electrostatic interactions with approaching substrates. The enzyme kinetics of recombinant PfSCS (wild‐type), generated by refolding of the individual P. falciparum SCSβ and Blastocystis SCSα subunits, demonstrated ADP‐forming activity (KmATP = 48 µm). Further, the introduction of charged gatekeeper residues, either positive (Lys and Lys) or negative (Glu and Asp), resulted in significant reductions in the ATP affinity of PfSCS. It is interesting to note that the recombinant PfSCSβ subunit can be refolded to a functional enzyme conformation using Blastocystis SCSα, indicating the possibility of subunits swapping among different organisms. These results concluded that electrostatic interactions at the gatekeeper region alone are insufficient to alter the substrate specificity of PfSCS, and further structural analysis with a particular focus on binding site architecture is required.
Collapse
Affiliation(s)
- Kapil Vashisht
- Protein Biochemistry and Engineering Lab, Parasite-Host Biology Group, ICMR-National Institute of Malaria Research, New Delhi, India
| | - Pallavi Singh
- Protein Biochemistry and Engineering Lab, Parasite-Host Biology Group, ICMR-National Institute of Malaria Research, New Delhi, India
| | - Sonia Verma
- Protein Biochemistry and Engineering Lab, Parasite-Host Biology Group, ICMR-National Institute of Malaria Research, New Delhi, India
| | - Rajnikant Dixit
- Protein Biochemistry and Engineering Lab, Parasite-Host Biology Group, ICMR-National Institute of Malaria Research, New Delhi, India
| | - Neelima Mishra
- Protein Biochemistry and Engineering Lab, Parasite-Host Biology Group, ICMR-National Institute of Malaria Research, New Delhi, India
| | - Kailash C Pandey
- Protein Biochemistry and Engineering Lab, Parasite-Host Biology Group, ICMR-National Institute of Malaria Research, New Delhi, India
| |
Collapse
|
7
|
Vashisht K, Verma S, Gupta S, Lynn AM, Dixit R, Mishra N, Valecha N, Hamblin KA, Maytum R, Pandey KC, van der Giezen M. Engineering Nucleotide Specificity of Succinyl-CoA Synthetase in Blastocystis: The Emerging Role of Gatekeeper Residues. Biochemistry 2017; 56:534-542. [PMID: 27478903 PMCID: PMC5404824 DOI: 10.1021/acs.biochem.6b00098] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
![]()
Charged,
solvent-exposed residues at the entrance to the substrate
binding site (gatekeeper residues) produce electrostatic dipole interactions
with approaching substrates, and control their access by a novel mechanism
called “electrostatic gatekeeper effect”. This proof-of-concept
study demonstrates that the nucleotide specificity can be engineered
by altering the electrostatic properties of the gatekeeper residues
outside the binding site. Using Blastocystis succinyl-CoA
synthetase (SCS, EC 6.2.1.5), we demonstrated that the gatekeeper
mutant (ED) resulted in ATP-specific SCS to show high GTP specificity.
Moreover, nucleotide binding site mutant (LF) had no effect on GTP
specificity and remained ATP-specific. However, via combination of
the gatekeeper mutant with the nucleotide binding site mutant (ED+LF),
a complete reversal of nucleotide specificity was obtained with GTP,
but no detectable activity was obtained with ATP. This striking result
of the combined mutant (ED+LF) was due to two changes; negatively
charged gatekeeper residues (ED) favored GTP access, and nucleotide
binding site residues (LF) altered ATP binding, which was consistent
with the hypothesis of the “electrostatic gatekeeper effect”.
These results were further supported by molecular modeling and simulation
studies. Hence, it is imperative to extend the strategy of the gatekeeper
effect in a different range of crucial enzymes (synthetases, kinases,
and transferases) to engineer substrate specificity for various industrial
applications and substrate-based drug design.
Collapse
Affiliation(s)
- Kapil Vashisht
- Host-parasite Interaction Biology Group, National Institute of Malaria Research, ICMR , New Delhi 110077, India
| | - Sonia Verma
- Host-parasite Interaction Biology Group, National Institute of Malaria Research, ICMR , New Delhi 110077, India
| | - Sunita Gupta
- School of Computational and Integrative Sciences, Jawaharlal Nehru University , New Delhi 110067, India
| | - Andrew M Lynn
- School of Computational and Integrative Sciences, Jawaharlal Nehru University , New Delhi 110067, India
| | - Rajnikant Dixit
- Host-parasite Interaction Biology Group, National Institute of Malaria Research, ICMR , New Delhi 110077, India
| | - Neelima Mishra
- Host-parasite Interaction Biology Group, National Institute of Malaria Research, ICMR , New Delhi 110077, India
| | - Neena Valecha
- Host-parasite Interaction Biology Group, National Institute of Malaria Research, ICMR , New Delhi 110077, India
| | | | - Robin Maytum
- School of Life Sciences, University of Bedfordshire , University Square, Luton LU1 3JU, U.K
| | - Kailash C Pandey
- Host-parasite Interaction Biology Group, National Institute of Malaria Research, ICMR , New Delhi 110077, India
| | | |
Collapse
|
8
|
Sinha A, Ray A, Ganguly S, Ghosh Dastidar S, Sarkar S. Variation in the ribosome interacting loop of the Sec61α from Giardia lamblia. Biol Direct 2015; 10:56. [PMID: 26424409 PMCID: PMC4588681 DOI: 10.1186/s13062-015-0087-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2015] [Accepted: 09/24/2015] [Indexed: 11/17/2022] Open
Abstract
The interaction between the ribosome and the endoplasmic reticulum-located Sec61 protein translocon is mediated through an arginine residue of Sec61α, which is conserved in all prokaryotic and eukaryotic orthologues characterized to date. Using in silico approaches we report that instead of arginine, this ribosome-interaction function is most likely discharged by a lysine residue in the protist Giardia lamblia. This functional substitution of the R with a K in GlSec61α may have taken place to accommodate a G-rich rRNA.
Collapse
Affiliation(s)
- Abhishek Sinha
- Department of Biochemistry, Bose Institute, P-1/12 CIT Road, Scheme VII M, Kolkata, 700054,, West Bengal, India.
| | - Atrayee Ray
- Department of Biochemistry, Bose Institute, P-1/12 CIT Road, Scheme VII M, Kolkata, 700054,, West Bengal, India.
| | - Sandipan Ganguly
- Molecular Parasitology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T Road, Scheme XM, Kolkata, 700010,, West Bengal, India.
| | - Shubhra Ghosh Dastidar
- Bioinformatics Center, Bose Institute, P-1/12 CIT Scheme VII M, Kolkata, 700054,, West Bengal, India.
| | - Srimonti Sarkar
- Department of Biochemistry, Bose Institute, P-1/12 CIT Road, Scheme VII M, Kolkata, 700054,, West Bengal, India.
| |
Collapse
|
9
|
Das M, Basu G. Protein-protein association rates captured in a single geometric parameter. Proteins 2015; 83:1557-62. [DOI: 10.1002/prot.24860] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2015] [Revised: 06/04/2015] [Accepted: 07/02/2015] [Indexed: 02/05/2023]
Affiliation(s)
- Madhurima Das
- Department of Biophysics; Bose Institute; P-1/12 CIT Scheme VIIM Kolkata 700054 India
| | - Gautam Basu
- Department of Biophysics; Bose Institute; P-1/12 CIT Scheme VIIM Kolkata 700054 India
| |
Collapse
|
10
|
Usha S, Selvaraj S. Structure-wise discrimination of adenine and guanine by proteins on the basis of their nonbonded interactions. J Biomol Struct Dyn 2014; 33:1474-92. [DOI: 10.1080/07391102.2014.958759] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
11
|
Wicky BIM, Stenta M, Dal Peraro M. TAL effectors specificity stems from negative discrimination. PLoS One 2013; 8:e80261. [PMID: 24282528 PMCID: PMC3840011 DOI: 10.1371/journal.pone.0080261] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2013] [Accepted: 10/01/2013] [Indexed: 01/24/2023] Open
Abstract
Transcription Activator-Like (TAL) effectors are DNA-binding proteins secreted by phytopathogenic bacteria that interfere with native cellular functions by binding to plant DNA promoters. The key element of their architecture is a domain of tandem-repeats with almost identical sequences. Most of the polymorphism is located at two consecutive amino acids termed Repeat Variable Diresidue (RVD). The discovery of a direct link between the RVD composition and the targeted nucleotide allowed the design of TAL-derived DNA-binding tools with programmable specificities that revolutionized the field of genome engineering. Despite structural data, the molecular origins of this specificity as well as the recognition mechanism have remained unclear. Molecular simulations of the recent crystal structures suggest that most of the protein-DNA binding energy originates from non-specific interactions between the DNA backbone and non-variable residues, while RVDs contributions are negligible. Based on dynamical and energetic considerations we postulate that, while the first RVD residue promotes helix breaks – allowing folding of TAL as a DNA-wrapping super-helix – the second provides specificity through a negative discrimination of matches. Furthermore, we propose a simple pharmacophore-like model for the rationalization of RVD-DNA interactions and the interpretation of experimental findings concerning shared affinities and binding efficiencies. The explanatory paradigm presented herein provides a better comprehension of this elegant architecture and we hope will allow for improved designs of TAL-derived biotechnological tools.
Collapse
Affiliation(s)
- Basile I. M. Wicky
- Laboratory for Biomolecular Modeling, Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Marco Stenta
- Laboratory for Biomolecular Modeling, Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Matteo Dal Peraro
- Laboratory for Biomolecular Modeling, Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Swiss Institute of Bioinformatics (SIB), Lausanne, Switzerland
- * E-mail:
| |
Collapse
|
12
|
Chakraborty S, Rao BJ, Baker N, Asgeirsson B. Structural phylogeny by profile extraction and multiple superimposition using electrostatic congruence as a discriminator. INTRINSICALLY DISORDERED PROTEINS 2013; 1. [PMID: 25364645 PMCID: PMC4212511 DOI: 10.4161/idp.25463] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Phylogenetic analysis of proteins using multiple sequence alignment (MSA) assumes an underlying evolutionary relationship in these proteins which occasionally remains undetected due to considerable sequence divergence. Structural alignment programs have been developed to unravel such fuzzy relationships. However, none of these structure based methods have used electrostatic properties to discriminate between spatially equivalent residues. We present a methodology for MSA of a set of related proteins with known structures using electrostatic properties as an additional discriminator (STEEP). STEEP first extracts a profile, then generates a multiple structural superimposition providing a consolidated spatial framework for comparing residues and finally emits the MSA. Residues that are aligned differently by including or excluding electrostatic properties can be targeted by directed evolution experiments to transform the enzymatic properties of one protein into another. We have compared STEEP results to those obtained from a MSA program (ClustalW) and a structural alignment method (MUSTANG) for chymotrypsin serine proteases. Subsequently, we used PhyML to generate phylogenetic trees for the serine and metallo-β-lactamase superfamilies from the STEEP generated MSA, and corroborated the accepted relationships in these superfamilies. We have observed that STEEP acts as a functional classifier when electrostatic congruence is used as a discriminator, and thus identifies potential targets for directed evolution experiments. In summary, STEEP is unique among phylogenetic methods for its ability to use electrostatic congruence to specify mutations that might be the source of the functional divergence in a protein family. Based on our results, we also hypothesize that the active site and its close vicinity contains enough information to infer the correct phylogeny for related proteins.
Collapse
Affiliation(s)
- Sandeep Chakraborty
- Department of Biological Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400 005, India
| | - Basuthkar J Rao
- Department of Biological Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400 005, India
| | - Nathan Baker
- Pacific Northwest National Laboratory, 902 Battelle Blvd, Richland, Washington 99354, United States
| | - Bjarni Asgeirsson
- Science Institute, Department of Biochemistry, University of Iceland, Dunhaga 3, IS-107 Reykjavik, Iceland
| |
Collapse
|
13
|
Kondabagil K, Draper B, Rao VB. Adenine recognition is a key checkpoint in the energy release mechanism of phage T4 DNA packaging motor. J Mol Biol 2011; 415:329-42. [PMID: 22100308 DOI: 10.1016/j.jmb.2011.11.013] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2011] [Revised: 11/03/2011] [Accepted: 11/07/2011] [Indexed: 11/25/2022]
Abstract
ATP is the source of energy for numerous biochemical reactions in all organisms. Tailed bacteriophages use ATP to drive powerful packaging machines that translocate viral DNA into a procapsid and compact it to near-crystalline density. Here we report that a complex network of interactions dictates adenine recognition and ATP hydrolysis in the pentameric phage T4 large "terminase" (gp17) motor. The network includes residues that form hydrogen bonds at the edges of the adenine ring (Q138 and Q143), base-stacking interactions at the plane of the ring (I127 and R140), and cross-talking bonds between adenine, triphosphate, and Walker A P-loop (Y142, Q143, and R140). These interactions are conserved in other translocases such as type I/type III restriction enzymes and SF1/SF2 helicases. Perturbation of any of these interactions, even the loss of a single hydrogen bond, leads to multiple defects in motor functions. Adenine recognition is therefore a key checkpoint that ensures efficient ATP firing only when the fuel molecule is precisely engaged with the motor. This may be a common feature in the energy release mechanism of ATP-driven molecular machines that carry out numerous biomolecular reactions in biological systems.
Collapse
Affiliation(s)
- Kiran Kondabagil
- Department of Biology, The Catholic University of America, 620 Michigan Avenue NE, Washington, DC 20064, USA
| | | | | |
Collapse
|
14
|
Chauhan JS, Mishra NK, Raghava GPS. Prediction of GTP interacting residues, dipeptides and tripeptides in a protein from its evolutionary information. BMC Bioinformatics 2010; 11:301. [PMID: 20525281 PMCID: PMC3098072 DOI: 10.1186/1471-2105-11-301] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2010] [Accepted: 06/03/2010] [Indexed: 11/17/2022] Open
Abstract
Background Guanosine triphosphate (GTP)-binding proteins play an important role in regulation of G-protein. Thus prediction of GTP interacting residues in a protein is one of the major challenges in the field of the computational biology. In this study, an attempt has been made to develop a computational method for predicting GTP interacting residues in a protein with high accuracy (Acc), precision (Prec) and recall (Rc). Result All the models developed in this study have been trained and tested on a non-redundant (40% similarity) dataset using five-fold cross-validation. Firstly, we have developed neural network based models using single sequence and PSSM profile and achieved maximum Matthews Correlation Coefficient (MCC) 0.24 (Acc 61.30%) and 0.39 (Acc 68.88%) respectively. Secondly, we have developed a support vector machine (SVM) based models using single sequence and PSSM profile and achieved maximum MCC 0.37 (Prec 0.73, Rc 0.57, Acc 67.98%) and 0.55 (Prec 0.80, Rc 0.73, Acc 77.17%) respectively. In this work, we have introduced a new concept of predicting GTP interacting dipeptide (two consecutive GTP interacting residues) and tripeptide (three consecutive GTP interacting residues) for the first time. We have developed SVM based model for predicting GTP interacting dipeptides using PSSM profile and achieved MCC 0.64 with precision 0.87, recall 0.74 and accuracy 81.37%. Similarly, SVM based model have been developed for predicting GTP interacting tripeptides using PSSM profile and achieved MCC 0.70 with precision 0.93, recall 0.73 and accuracy 83.98%. Conclusion These results show that PSSM based method performs better than single sequence based method. The prediction models based on dipeptides or tripeptides are more accurate than the traditional model based on single residue. A web server "GTPBinder" http://www.imtech.res.in/raghava/gtpbinder/ based on above models has been developed for predicting GTP interacting residues in a protein.
Collapse
Affiliation(s)
- Jagat S Chauhan
- Bioinformatics Centre, Institute of Microbial Technology (IMTECH), Chandigarh, India
| | | | | |
Collapse
|
15
|
Kahraman A, Morris RJ, Laskowski RA, Favia AD, Thornton JM. On the diversity of physicochemical environments experienced by identical ligands in binding pockets of unrelated proteins. Proteins 2010; 78:1120-36. [PMID: 19927322 DOI: 10.1002/prot.22633] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Most function prediction methods that identify cognate ligands from binding site analyses work on the assumption of molecular complementarity. These approaches build on the conjectured complementarity of geometrical and physicochemical properties between ligands and binding sites so that similar binding sites will bind similar ligands. We found that this assumption does not generally hold for protein-ligand interactions and observed that it is not the chemical composition of ligand molecules that dictates the complementarity between protein and ligand molecules, but that the ligand's share within the functional mechanism of a protein determines the degree of complementarity. Here, we present for a set of cognate ligands a descriptive analysis and comparison of the physicochemical properties that each ligand experiences in various nonhomologous binding pockets. The comparisons in each ligand set reveal large variations in their experienced physicochemical properties, suggesting that the same ligand can bind to distinct physicochemical environments. In some protein ligand complexes, the variation was found to correlate with the electrochemical characteristic of ligand molecules, whereas in others it was disclosed as a prerequisite for the biochemical function of the protein. To achieve binding, proteins were observed to engage in subtle balancing acts between electrostatic and hydrophobic interactions to generate stabilizing free energies of binding. For the presented analysis, a new method for scoring hydrophobicity from molecular environments was developed showing high correlations with experimental determined desolvation energies. The presented results highlight the complexities of molecular recognition and underline the challenges of computational structural biology in developing methods to detect these important subtleties.
Collapse
Affiliation(s)
- Abdullah Kahraman
- European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridgeshire CB10 1SD, United Kingdom.
| | | | | | | | | |
Collapse
|
16
|
Das M, Basu G. Coulomb energies of protein-protein complexes with monopole-free charge distributions. J Mol Graph Model 2008; 27:846-51. [PMID: 19167253 DOI: 10.1016/j.jmgm.2008.12.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2008] [Revised: 12/11/2008] [Accepted: 12/12/2008] [Indexed: 01/06/2023]
Abstract
From an analysis of Coulomb energy distributions of a large set of protein-protein complexes we show that the positive tail in the energy distribution disappears when the monopole-monopole term, the only energy term independent of inter-subunit orientations, is removed. This indicates that unfavorable Coulomb energies associated with subunit orientations are excluded in protein-protein complexes. The overall result remained unchanged when solvent effects were included. Our results have important bearing on the restriction of subunit orientations in protein-protein complexes and complement a recent work [K. Brock, K. Talley, K. Coley, P. Kundrotas, E. Alexov, Optimization of electrostatic interactions in protein-protein complexes, Biophys. J. 93 (2007) 3340-3352.] which showed that Coulomb energy of interaction in protein-protein complexes is sequence-optimized.
Collapse
Affiliation(s)
- Madhurima Das
- Bioinformatics Centre, Bose Institute, P-1/12 CIT Scheme VIIM, Kolkata 700054, India
| | | |
Collapse
|
17
|
Crystal structure of a translation termination complex formed with release factor RF2. Proc Natl Acad Sci U S A 2008; 105:19684-9. [PMID: 19064930 DOI: 10.1073/pnas.0810953105] [Citation(s) in RCA: 187] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We report the crystal structure of a translation termination complex formed by the Thermus thermophilus 70S ribosome bound with release factor RF2, in response to a UAA stop codon, solved at 3 A resolution. The backbone of helix alpha5 and the side chain of serine of the conserved SPF motif of RF2 recognize U1 and A2 of the stop codon, respectively. A3 is unstacked from the first 2 bases, contacting Thr-216 and Val-203 of RF2 and stacking on G530 of 16S rRNA. The structure of the RF2 complex supports our previous proposal that conformational changes in the ribosome in response to recognition of the stop codon stabilize rearrangement of the switch loop of the release factor, resulting in docking of the universally conserved GGQ motif in the PTC of the 50S subunit. As seen for the RF1 complex, the main-chain amide nitrogen of glutamine in the GGQ motif is positioned to contribute directly to catalysis of peptidyl-tRNA hydrolysis, consistent with mutational studies, which show that most side-chain substitutions of the conserved glutamine have little effect. We show that when the H-bonding capability of the main-chain N-H of the conserved glutamine is eliminated by substitution with proline, peptidyl-tRNA esterase activity is abolished, consistent with its proposed role in catalysis.
Collapse
|
18
|
Hamblin K, Standley DM, Rogers MB, Stechmann A, Roger AJ, Maytum R, van der Giezen M. Localization and nucleotide specificity of Blastocystis succinyl-CoA synthetase. Mol Microbiol 2008; 68:1395-405. [PMID: 18452512 PMCID: PMC2440562 DOI: 10.1111/j.1365-2958.2008.06228.x] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
The anaerobic lifestyle of the intestinal parasite Blastocystis raises questions about the biochemistry and function of its mitochondria-like organelles. We have characterized the Blastocystis succinyl-CoA synthetase (SCS), a tricarboxylic acid cycle enzyme that conserves energy by substrate-level phosphorylation. We show that SCS localizes to the enigmatic Blastocystis organelles, indicating that these organelles might play a similar role in energy metabolism as classic mitochondria. Although analysis of residues inside the nucleotide-binding site suggests that Blastocystis SCS is GTP-specific, we demonstrate that it is ATP-specific. Homology modelling, followed by flexible docking and molecular dynamics simulations, indicates that while both ATP and GTP fit into the Blastocystis SCS active site, GTP is destabilized by electrostatic dipole interactions with Lys 42 and Lys 110, the side-chains of which lie outside the nucleotide-binding cavity. It has been proposed that residues in direct contact with the substrate determine nucleotide specificity in SCS. However, our results indicate that, in Blastocystis, an electrostatic gatekeeper controls which ligands can enter the binding site.
Collapse
Affiliation(s)
- Karleigh Hamblin
- School of Biological and Chemical Sciences, Queen Mary, University of London, Mile End Road, London E1 4NS, UK
| | | | | | | | | | | | | |
Collapse
|
19
|
Lee JE, Bae E, Bingman CA, Phillips GN, Raines RT. Structural basis for catalysis by onconase. J Mol Biol 2007; 375:165-77. [PMID: 18001769 DOI: 10.1016/j.jmb.2007.09.089] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2007] [Accepted: 09/20/2007] [Indexed: 11/18/2022]
Abstract
Onconase (ONC) is a homolog of bovine pancreatic ribonuclease (RNase A) from the frog Rana pipiens. ONC displays antitumoral activity and is in advanced clinical trials for the treatment of cancer. Here, we report the first atomic structures of ONC-nucleic acid complexes: a T89N/E91A ONC-5'-AMP complex at 1.65 A resolution and a wild-type ONC-d(AUGA) complex at 1.90 A resolution. The latter structure and site-directed mutagenesis were used to reveal the atomic basis for substrate recognition and turnover by ONC. The residues in ONC that are proximal to the scissile phosphodiester bond (His10, Lys31, and His97) and uracil nucleobase (Thr35, Asp67, and Phe98) are conserved from RNase A and serve to generate a similar bell-shaped pH versus k(cat)/K(M) profile for RNA cleavage. Glu91 of ONC forms two hydrogen bonds with the guanine nucleobase in d(AUGA), and Thr89 is in close proximity to that nucleobase. Installing a neutral or cationic residue at position 91 or an asparagine residue at position 89 virtually eliminated the 10(2)-fold guanine:adenine preference of ONC. A variant that combined such substitutions, T89N/E91A ONC, actually preferred adenine over guanine. In contrast, installing an arginine residue at position 91 increased the guanine preference and afforded an ONC variant with the highest known k(cat)/K(M) value. These data indicate that ONC discriminates between guanine and adenine by using Coulombic interactions and a network of hydrogen bonds. The structure of the ONC-d(AUGA) complex was also used to probe other aspects of catalysis. For example, the T5R substitution, designed to create a favorable Coulombic interaction between ONC and a phosphoryl group in RNA, increased ribonucleolytic activity by twofold. No variant, however, was more toxic to human cancer cells than wild-type ONC. Together, these findings provide a cynosure for understanding catalysis of RNA cleavage in a system of high medicinal relevance.
Collapse
Affiliation(s)
- J Eugene Lee
- Department of Biochemistry, University of Wisconsin-Madison, Madison, WI 53706-1544, USA
| | | | | | | | | |
Collapse
|
20
|
Lee AM, Singleton SF. Intersubunit electrostatic complementarity in the RecA nucleoprotein filament regulates nucleotide substrate specificity and conformational activation. Biochemistry 2006; 45:4514-29. [PMID: 16584187 DOI: 10.1021/bi052296x] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Escherichia coli RecA protein is the prototypical member of a family of molecular motors that transduces ATP binding and hydrolysis for mechanical function. While many general mechanistic features of RecA action are known, specific structural and functional insights into the molecular basis of RecA activation remain elusive. Toward a more complete understanding of the interdependence between ATP and DNA binding by RecA, we report the characterization of a mutant RecA protein wherein the aspartate residue at position 100 within the ATP binding site is replaced by arginine. Physiologically, D100R RecA was characterized by an inducible, albeit reduced, activation of the SOS response and a diminished ability to promote cellular survival after UV irradiation. Biochemically, the D100R substitution caused a surprisingly modest perturbation of RecA-ATP interactions and an unexpected and significant decrease in the affinity of RecA for ssDNA. Moreover, in vitro assays of RecA activities requiring the coordinated processing of ATP and DNA revealed (1) a 2-5-fold decrease in steady-state turnover of ATP; (2) no formation of mixed nucleoprotein filaments when wild-type and D100R RecA compete for limiting ssDNA; and (3) no formation of strand exchange reaction products. Taken together, these results suggest that the D100R mutational effects on isolated RecA activities combine synergistically to perturb its higher-order functions. We conclude that the replacement of Asp100 resulted in a change in the electrostatic complementarity between RecA monomers during active filament assembly that prevents the protein from fully accessing the active multimeric state.
Collapse
Affiliation(s)
- Andrew M Lee
- School of Pharmacy, Division of Medicinal Chemistry & Natural Products, University of North Carolina at Chapel Hill, CB #7360, Chapel Hill, North Carolina 27599-7360, USA
| | | |
Collapse
|
21
|
Oparina NJ, Kalinina OV, Gelfand MS, Kisselev LL. Common and specific amino acid residues in the prokaryotic polypeptide release factors RF1 and RF2: possible functional implications. Nucleic Acids Res 2005; 33:5226-34. [PMID: 16162810 PMCID: PMC1214553 DOI: 10.1093/nar/gki841] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Termination of protein synthesis is promoted in ribosomes by proper stop codon discrimination by class 1 polypeptide release factors (RFs). A large set of prokaryotic RFs differing in stop codon specificity, RF1 for UAG and UAA, and RF2 for UGA and UAA, was analyzed by means of a recently developed computational method allowing identification of the specificity-determining positions (SDPs) in families composed of proteins with similar but not identical function. Fifteen SDPs were identified within the RF1/2 superdomain II/IV known to be implicated in stop codon decoding. Three of these SDPs had particularly high scores. Five residues invariant for RF1 and RF2 [invariant amino acid residues (IRs)] were spatially clustered with the highest-scoring SDPs that in turn were located in two zones within the SDP/IR area. Zone 1 (domain II) included PxT and SPF motifs identified earlier by others as ‘discriminator tripeptides’. We suggest that IRs in this zone take part in the recognition of U, the first base of all stop codons. Zone 2 (domain IV) possessed two SDPs with the highest scores not identified earlier. Presumably, they also take part in stop codon binding and discrimination. Elucidation of potential functional role(s) of the newly identified SDP/IR zones requires further experiments.
Collapse
Affiliation(s)
- Nina J. Oparina
- To whom correspondence should be addressed. Tel: +7 095 1351419; Fax: +7 095 1351405;
| | - Olga V. Kalinina
- Department of Bioengineering and Bioinformatics, Moscow State UniversityVorob'evy Gory, 1-73, Moscow 119992, Russia
| | - Mikhail S. Gelfand
- Institute for Information Transmission Problems, Russian Academy of SciencesBolshoi Karetnyi per., 19, Moscow 127994, Russia
- State Scientific Centre GosNIIGenetika1st Dorozhny pr. 1, Moscow, 113545, Russia
| | | |
Collapse
|