1
|
Deshwal N, Singh MB, Bahadur I, Kaushik N, Kaushik NK, Singh P, Kumari K. A review on recent advancements on removal of harmful metal/metal ions using graphene oxide: Experimental and theoretical approaches. THE SCIENCE OF THE TOTAL ENVIRONMENT 2023; 858:159672. [PMID: 36306838 DOI: 10.1016/j.scitotenv.2022.159672] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2022] [Revised: 10/16/2022] [Accepted: 10/19/2022] [Indexed: 06/16/2023]
Abstract
Graphene oxide is a two-dimensional carbon nanomaterial and has gained huge popularity over the last decade. Because, the graphene oxide can be dispersed in water easily and it is one of the most researched two-dimensional materials in the current time. The extraordinary properties shown by graphene oxide (GO) are due to its unique chemical structure; includes various hydrophilic functional groups containing oxygen such as carboxyl, hydroxyl, carbonyl and tiny sp2 carbon domains surrounded by sp3 domains. These groups are very peculiar for various applications as they allow covalent functionalisation with a plethora of compounds. Large surface area, intrinsic fluorescence, excellent surface functionality, amphiphilicity, improved conductivity, high adsorption capacity and superior biocompatibility are some of the chemical properties have drawn research from various fields. Graphene oxide has various interactions such as coordination, chelation, hydrogen bonding, electrostatic interaction, hydrophobic effects, π-π interaction, acid base interaction etc., with various metal ions. This review is focused on the removal of metals and metal ions due to their interactions mentioned above. Further, potential of composites of graphene oxide in the removal of metal and metal ions is also discussed. Further, the current challenges in this field at industrial-scale are also discussed.
Collapse
Affiliation(s)
- Nidhi Deshwal
- Department of Chemistry, Atma Ram Sanatan Dharma College, University of Delhi, New Delhi, India
| | - Madhur Babu Singh
- Department of Chemistry, Atma Ram Sanatan Dharma College, University of Delhi, New Delhi, India
| | - Indra Bahadur
- Department of Chemistry, Faculty of Natural and Agricultural Sciences, North-West University, South Africa
| | - Neha Kaushik
- Department of Biotechnology, College of Engineering, The University of Suwon, Hwaseong 18323, South Korea
| | - Nagendra Kumar Kaushik
- Plasma Bioscience Research Center, Department of Electrical and Biological Physics, Kwangwoon University, Seoul 01897, South Korea.
| | - Prashant Singh
- Department of Chemistry, Atma Ram Sanatan Dharma College, University of Delhi, New Delhi, India.
| | - Kamlesh Kumari
- Department of Zoology, University of Delhi, Delhi, India.
| |
Collapse
|
2
|
Singh MB, Sharma R, Kumar D, Khanna P, Mansi, Khanna L, Kumar V, Kumari K, Gupta A, Chaudhary P, Kaushik N, Choi EH, Kaushik NK, Singh P. An understanding of coronavirus and exploring the molecular dynamics simulations to find promising candidates against the Mpro of nCoV to combat the COVID-19: A systematic review. J Infect Public Health 2022; 15:1326-1349. [PMID: 36288640 PMCID: PMC9579205 DOI: 10.1016/j.jiph.2022.10.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2022] [Revised: 10/08/2022] [Accepted: 10/12/2022] [Indexed: 11/10/2022] Open
Abstract
The first infection case of new coronavirus was reported at the end of 2019 and after then, the cases are reported in all nations across the world in a very short period. Further, the regular news of mutations in the virus has made life restricted with appropriate behavior. To date, a new strain (Omicron and its new subvariant Omicron XE) has brought fear amongst us due to a higher trajectory of increase in the number of cases. The researchers thus started giving attention to this viral infection and discovering drug-like candidates to cure the infections. Finding a drug for any viral infection is not an easy task and takes plenty of time. Therefore, computational chemistry/bioinformatics is followed to get promising molecules against viral infection. Molecular dynamics (MD) simulations are being explored to get drug candidates in a short period. The molecules are screened via molecular docking, which provides preliminary information which can be further verified by molecular dynamics (MD) simulations. To understand the change in structure, MD simulations generated several trajectories such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), hydrogen bonding, and radius of gyration for the main protease (Mpro) of the new coronavirus (nCoV) in the presence of small molecules. Additionally, change in free energy for the formation of complex of Mpro of nCoV with the small molecule can be determined by applying molecular mechanics with generalized born and surface area solvation (MM-GBSA). Thus, the promising molecules can be further explored for clinical trials to combat coronavirus disease-19 (COVID-19).
Collapse
Affiliation(s)
- Madhur Babu Singh
- Department of Chemistry, Atma Ram Sanatan Dharma College, University of Delhi, New Delhi, India
| | - Ritika Sharma
- Department of Biochemistry, University of Delhi, New Delhi, India
| | - Durgesh Kumar
- Department of Chemistry, Maitreyi College, University of Delhi, Delhi, India
| | - Pankaj Khanna
- Department of Chemistry, Acharya Narendra Dev College, University of Delhi, New Delhi, India
| | - Mansi
- University School of Basic and Applied Sciences, Guru Gobind Singh Indraprastha University, New Delhi, India
| | - Leena Khanna
- University School of Basic and Applied Sciences, Guru Gobind Singh Indraprastha University, New Delhi, India
| | - Vinod Kumar
- Special Centre for Nanoscience (SCNS), Jawaharlal Nehru University, New Delhi, India
| | - Kamlesh Kumari
- Department of Zoology, University of Delhi, New Delhi, India
| | - Akanksha Gupta
- Department of Chemistry, Sri Venkateswara College, University of Delhi, New Delhi, India
| | - Preeti Chaudhary
- Department of Chemistry, Atma Ram Sanatan Dharma College, University of Delhi, New Delhi, India
| | - Neha Kaushik
- Department of Biotechnology, College of Engineering, The University of Suwon, Hwaseong-si 18323, Republic of Korea.
| | - Eun Ha Choi
- Plasma Bioscience Research Center, Department of Electrical and Biological Physics, Kwangwoon University, Seoul 01897, Republic of Korea
| | - Nagendra Kumar Kaushik
- Plasma Bioscience Research Center, Department of Electrical and Biological Physics, Kwangwoon University, Seoul 01897, Republic of Korea.
| | - Prashant Singh
- Department of Chemistry, Atma Ram Sanatan Dharma College, University of Delhi, New Delhi, India.
| |
Collapse
|
3
|
Computational studies of anaplastic lymphoma kinase mutations reveal common mechanisms of oncogenic activation. Proc Natl Acad Sci U S A 2021; 118:2019132118. [PMID: 33674381 PMCID: PMC7958353 DOI: 10.1073/pnas.2019132118] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
High-risk tumors are genomically heterogeneous, harboring gene amplifications and mutations. The activation status of mutated proteins in cancer can profoundly impact disease progression, patient response, and drug sensitivity. Yet, outside of a few hotspot mutations, functional studies of clinically observed mutations are not commonly pursued. We report a combined experimental profiling and computational analysis of the effects of clinically observed and “test” mutations in the kinase domain of anaplastic lymphoma kinase (ALK), a known oncogenic driver in pediatric neuroblastoma. We find that the activation status of the mutated protein is a good indicator of the transforming ability in NIH 3T3 cells. We also report biophysical as well as data-driven models with predictive power to profile these mutant kinases in silico. Kinases play important roles in diverse cellular processes, including signaling, differentiation, proliferation, and metabolism. They are frequently mutated in cancer and are the targets of a large number of specific inhibitors. Surveys of cancer genome atlases reveal that kinase domains, which consist of 300 amino acids, can harbor numerous (150 to 200) single-point mutations across different patients in the same disease. This preponderance of mutations—some activating, some silent—in a known target protein make clinical decisions for enrolling patients in drug trials challenging since the relevance of the target and its drug sensitivity often depend on the mutational status in a given patient. We show through computational studies using molecular dynamics (MD) as well as enhanced sampling simulations that the experimentally determined activation status of a mutated kinase can be predicted effectively by identifying a hydrogen bonding fingerprint in the activation loop and the αC-helix regions, despite the fact that mutations in cancer patients occur throughout the kinase domain. In our study, we find that the predictive power of MD is superior to a purely data-driven machine learning model involving biochemical features that we implemented, even though MD utilized far fewer features (in fact, just one) in an unsupervised setting. Moreover, the MD results provide key insights into convergent mechanisms of activation, primarily involving differential stabilization of a hydrogen bond network that engages residues of the activation loop and αC-helix in the active-like conformation (in >70% of the mutations studied, regardless of the location of the mutation).
Collapse
|
4
|
Geronimo I, Vidossich P, Donati E, Vivo M. Computational investigations of polymerase enzymes: Structure, function, inhibition, and biotechnology. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1534] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Inacrist Geronimo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| | - Pietro Vidossich
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| | - Elisa Donati
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| | - Marco Vivo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| |
Collapse
|
5
|
Schlick T, Portillo-Ledesma S. Biomolecular modeling thrives in the age of technology. NATURE COMPUTATIONAL SCIENCE 2021; 1:321-331. [PMID: 34423314 PMCID: PMC8378674 DOI: 10.1038/s43588-021-00060-9] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 03/22/2021] [Indexed: 12/12/2022]
Abstract
The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, NY, USA
- Courant Institute of Mathematical Sciences, New York University, New York, NY, USA
- New York University–East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai, China
| | | |
Collapse
|
6
|
Bolhuis PG, Swenson DWH. Transition Path Sampling as Markov Chain Monte Carlo of Trajectories: Recent Algorithms, Software, Applications, and Future Outlook. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202000237] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Peter G. Bolhuis
- Amsterdam Center for Multiscale Modeling van 't Hoff Institute for Molecular Sciences University of Amsterdam PO Box 94157 1090 GD Amsterdam The Netherlands
| | - David W. H. Swenson
- Centre Blaise Pascal Ecole Normale Superieure 46, allée d'Italie 69364 Lyon Cedex 07 France
| |
Collapse
|
7
|
Levintov L, Paul S, Vashisth H. Reaction Coordinate and Thermodynamics of Base Flipping in RNA. J Chem Theory Comput 2021; 17:1914-1921. [PMID: 33594886 DOI: 10.1021/acs.jctc.0c01199] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.
Collapse
Affiliation(s)
- Lev Levintov
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, United States
| | - Sanjib Paul
- Department of Chemistry, New York University, New York 10003, New York, United States
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, United States
| |
Collapse
|
8
|
Hussain S, Haji-Akbari A. Studying rare events using forward-flux sampling: Recent breakthroughs and future outlook. J Chem Phys 2020; 152:060901. [DOI: 10.1063/1.5127780] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Sarwar Hussain
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
9
|
Paul S, Nair NN, Vashisth H. Phase space and collective variable based simulation methods for studies of rare events. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1634268] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| | - Nisanth N. Nair
- Department of Chemistry, Indian Institute of Technology, Kanpur, India
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| |
Collapse
|
10
|
Maffeo C, Chou HY, Aksimentiev A. Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations. ADVANCED THEORY AND SIMULATIONS 2019; 2:1800191. [PMID: 31728433 PMCID: PMC6855400 DOI: 10.1002/adts.201800191] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Indexed: 12/15/2022]
Abstract
Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.
Collapse
Affiliation(s)
- Christopher Maffeo
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Han-Yi Chou
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Aleksei Aksimentiev
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| |
Collapse
|
11
|
Paul S, Paul TK, Taraphder S. Orthogonal order parameters to model the reaction coordinate of an enzyme catalyzed reaction. J Mol Graph Model 2019; 90:18-32. [PMID: 30959266 DOI: 10.1016/j.jmgm.2019.03.023] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2018] [Revised: 03/01/2019] [Accepted: 03/26/2019] [Indexed: 01/03/2023]
Abstract
The choice of suitable collective variables in formulating an optimal reaction coordinate is a challenging task for activated transitions between a pair of stable states especially when dealing with biochemical changes such as enzyme catalyzed reactions. A detailed benchmarking study is carried out on the choice of collective variables that can distinguish between the stable states unambiguously. We specifically address the issue if these variables may be directly used to model the optimal reaction coordinate, or if it would be better to use their orthogonalized counterparts. The proposed computational scheme is applied to the rate determining intramolecular proton transfer step in the enzyme human carbonic anhydrase II. The optimum reaction coordinate is determined with and without orthogonalization of the collective variables pertinent to a key conformational fluctuation and the actual proton transfer step at the active site of the enzyme. Suitability of the predicted reaction coordinates in different processes is examined in terms of the free energy profile projected along the reaction coordinate, the rate constant of transition and the underlying molecular mechanism of barrier crossing. Our results indicate that a better agreement with earlier simulation and experimental data is obtained when the orthogonalized collective variables are used to model the reaction coordinate.
Collapse
Affiliation(s)
- Sanjib Paul
- Department of Chemistry, Indian Institute of Technology, Kharagpur, 721302, India
| | - Tanmoy Kumar Paul
- Department of Chemistry, Indian Institute of Technology, Kharagpur, 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology, Kharagpur, 721302, India.
| |
Collapse
|
12
|
Foley MC, Couto L, Rauf S, Boyke A. Insights into DNA polymerase δ’s mechanism for accurate DNA replication. J Mol Model 2019; 25:80. [DOI: 10.1007/s00894-019-3957-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2018] [Accepted: 02/05/2019] [Indexed: 11/28/2022]
|
13
|
Paul S, Paul TK, Taraphder S. Reaction Coordinate, Free Energy, and Rate of Intramolecular Proton Transfer in Human Carbonic Anhydrase II. J Phys Chem B 2018; 122:2851-2866. [DOI: 10.1021/acs.jpcb.7b10713] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Tanmoy Kumar Paul
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
14
|
|
15
|
Atis M, Johnson KA, Elber R. Pyrophosphate Release in the Protein HIV Reverse Transcriptase. J Phys Chem B 2017; 121:9557-9565. [PMID: 28926712 DOI: 10.1021/acs.jpcb.7b08320] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Enzymatic reactions usually occur in several steps: a step of substrate binding to the surface of the protein, a step of protein reorganization around the substrate and conduction of a chemical reaction, and a step of product release. The release of inorganic phosphate-PPi-from the matrix of the protein HIV reverse transcriptase is investigated computationally. Atomically detailed simulations with explicit solvent are analyzed to obtain the free energy profile, mean first passage time, and detailed molecular mechanisms of PPi escape. A challenge for the computations is of time scales. The experimental time scale of the process of interest is in milliseconds, and straightforward molecular dynamics simulations are in sub-microseconds. To overcome the time scale gap, we use the algorithm of Milestoning along a reaction coordinate to compute the overall free energy profile and rate. The methods of locally enhanced sampling and steered molecular dynamics determine plausible reaction coordinates. The observed molecular mechanism couples the transfer of the PPi to positively charged lysine side chains that are found on the exit pathway and to an exiting magnesium ion. In accord with experimental findings, the release rate is comparable to the chemical step, allowing for variations in substrate (DNA or RNA template) in which the release becomes rate determining.
Collapse
Affiliation(s)
- Murat Atis
- Institute for Computational Engineering and Sciences, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Kenneth A Johnson
- Department of Molecular Biosciences, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Ron Elber
- Institute for Computational Engineering and Sciences, The University of Texas at Austin , Austin, Texas 78712, United States.,Department of Chemistry, The University of Texas at Austin , Austin, Texas 78712, United States
| |
Collapse
|
16
|
Cao L, Cheng L, Tang T, Liu J, Zhang D, Xu P, Li G. Understanding enzyme reactions using enhanced sampling techniques. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2015.1132316] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
17
|
Kim T, Freudenthal BD, Beard WA, Wilson SH, Schlick T. Insertion of oxidized nucleotide triggers rapid DNA polymerase opening. Nucleic Acids Res 2016; 44:4409-24. [PMID: 27034465 PMCID: PMC4872097 DOI: 10.1093/nar/gkw174] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2016] [Accepted: 03/04/2016] [Indexed: 11/21/2022] Open
Abstract
A novel mechanism is unveiled to explain why a pro-mutagenic nucleotide lesion (oxidized guanine, 8-oxoG) causes the mammalian DNA repair polymerase-β (pol-β) to rapidly transition to an inactive open conformation. The mechanism involves unexpected features revealed recently in time-lapse crystallography. Specifically, a delicate water network associated with a lesion-stabilizing auxilliary product ion Mg(p) triggers a cascade of events that leads to poor active site geometry and the rupture of crucial molecular interactions between key residues in both the anti(8-oxoG:C) and syn(8-oxoG:A) systems. Once the base pairs in these lesioned systems are broken, dislocation of both Asp192 (a metal coordinating ligand) and the oxoG phosphate group (PO4) interfere with the hydrogen bonding between Asp192 and Arg258, whose rotation toward Asp192 is crucial to the closed-to-open enzyme transition. Energetically, the lesioned open states are similar in energy to those of the corresponding closed complexes after chemistry, in marked contrast to the unlesioned pol-β anti(G:C) system, whose open state is energetically higher than the closed state. The delicate surveillance system offers a fundamental protective mechanism in the cell that triggers DNA repair events which help deter insertion of oxidized lesions.
Collapse
Affiliation(s)
- Taejin Kim
- Department of Chemistry, New York University, 10th Floor Silver Center, 100 Washington Square East, New York, NY 10003, USA
| | - Bret D Freudenthal
- Genome Integrity and Structural Biology Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709, USA
| | - William A Beard
- Genome Integrity and Structural Biology Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709, USA
| | - Samuel H Wilson
- Genome Integrity and Structural Biology Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709, USA
| | - Tamar Schlick
- Department of Chemistry, New York University, 10th Floor Silver Center, 100 Washington Square East, New York, NY 10003, USA Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012, USA
| |
Collapse
|
18
|
Vuković L, Chipot C, Makino DL, Conti E, Schulten K. Molecular Mechanism of Processive 3' to 5' RNA Translocation in the Active Subunit of the RNA Exosome Complex. J Am Chem Soc 2016; 138:4069-78. [PMID: 26928279 PMCID: PMC4988868 DOI: 10.1021/jacs.5b12065] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Recent experimental studies revealed structural details of 3' to 5' degradation of RNA molecules, performed by the exosome complex. ssRNA is channeled through its multisubunit ring-like core into the active site tunnel of its key exonuclease subunit Rrp44, which acts both as an enzyme and a motor. Even in isolation, Rrp44 can pull and sequentially cleave RNA nucleotides, one at a time, without any external energy input and release a final 3-5 nucleotide long product. Using molecular dynamics simulations, we identify the main factors that control these processes. Our free energy calculations reveal that RNA transfer from solution into the active site of Rrp44 is highly favorable, but dependent on the length of the RNA strand. While RNA strands formed by 5 nucleotides or more correspond to a decreasing free energy along the translocation coordinate toward the cleavage site, a 4-nucleotide RNA experiences a free energy barrier along the same direction, potentially leading to incomplete cleavage of ssRNA and the release of short (3-5) nucleotide products. We provide new insight into how Rrp44 catalyzes a localized enzymatic reaction and performs an action distributed over several RNA nucleotides, leading eventually to the translocation of whole RNA segments into the position suitable for cleavage.
Collapse
Affiliation(s)
- Lela Vuković
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
- Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
- epartment of Chemistry, University of Texas at El Paso, El Paso, TX 79968, United States
| | - Christophe Chipot
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
- Laboratoire International Associé CNRS-University of Illinois, Université de Lorraine, Vandoeuvre-lès-Nancy, France
| | - Debora L. Makino
- Department of Structural Cell Biology, Max Planck Institute of Biochemistry, Martinsried, Germany
| | - Elena Conti
- Department of Structural Cell Biology, Max Planck Institute of Biochemistry, Martinsried, Germany
| | - Klaus Schulten
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
- Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States
| |
Collapse
|
19
|
Moscato B, Swain M, Loria JP. Induced Fit in the Selection of Correct versus Incorrect Nucleotides by DNA Polymerase β. Biochemistry 2015; 55:382-95. [PMID: 26678253 DOI: 10.1021/acs.biochem.5b01213] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
DNA polymerase β (Pol β) repairs single-nucleotide gapped DNA (sngDNA) by enzymatic incorporation of the Watson-Crick partner nucleotide at the gapped position opposite the templating nucleotide. The process by which the matching nucleotide is incorporated into a sngDNA sequence has been relatively well-characterized, but the process of discrimination from nucleotide misincorporation remains unclear. We report here NMR spectroscopic characterization of full-length, uniformly labeled Pol β in apo, sngDNA-bound binary, and ternary complexes containing matching and mismatching nucleotide. Our data indicate that, while binding of the correct nucleotide to the binary complex induces chemical shift changes consistent with the process of enzyme closure, the ternary Pol β complex containing a mismatching nucleotide exhibits no such changes and appears to remain in an open, unstable, binary-like conformation. Our findings support an induced-fit mechanism for polymerases in which a closed ternary complex can only be achieved in the presence of matching nucleotide.
Collapse
Affiliation(s)
- Beth Moscato
- Department of Chemistry, Yale University , 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Monalisa Swain
- Department of Chemistry, Yale University , 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - J Patrick Loria
- Department of Chemistry, Yale University , 225 Prospect Street, New Haven, Connecticut 06520, United States.,Department of Molecular Biophysics and Biochemistry, Yale University , 260 Whitney Avenue, New Haven, Connecticut 06520, United States
| |
Collapse
|
20
|
Perera L, Beard WA, Pedersen LG, Wilson SH. Applications of quantum mechanical/molecular mechanical methods to the chemical insertion step of DNA and RNA polymerization. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2014; 97:83-113. [PMID: 25458356 PMCID: PMC5573153 DOI: 10.1016/bs.apcsb.2014.10.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
We review theoretical attempts to model the chemical insertion reactions of nucleoside triphosphates catalyzed by the nucleic acid polymerases using combined quantum mechanical/molecular mechanical methodology. Due to an existing excellent database of high-resolution X-ray crystal structures, the DNA polymerase β system serves as a useful template for discussion and comparison. The convergence of structures of high-quality complexes and continued developments of theoretical techniques suggest a bright future for understanding the global features of nucleic acid polymerization.
Collapse
Affiliation(s)
- Lalith Perera
- Laboratory of Structural Biology, National Institution of Environmental Health Sciences, Research Triangle Park, North Carolina, USA.
| | - William A Beard
- Laboratory of Structural Biology, National Institution of Environmental Health Sciences, Research Triangle Park, North Carolina, USA
| | - Lee G Pedersen
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Samuel H Wilson
- Laboratory of Structural Biology, National Institution of Environmental Health Sciences, Research Triangle Park, North Carolina, USA
| |
Collapse
|
21
|
Yang Y, Qin Y, Ding Q, Bakhtina M, Wang L, Tsai MD, Zhong D. Ultrafast water dynamics at the interface of the polymerase-DNA binding complex. Biochemistry 2014; 53:5405-13. [PMID: 25105470 PMCID: PMC4148141 DOI: 10.1021/bi500810a] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
DNA polymerases slide on DNA during replication, and the interface must be mobile for various conformational changes. The role of lubricant interfacial water is not understood. In this report, we systematically characterize the water dynamics at the interface and in the active site of a tight binding polymerase (pol β) in its binary complex and ternary state using tryptophan as a local optical probe. Using femtosecond spectroscopy, we observed that upon DNA recognition the surface hydration water is significantly confined and becomes bound water at the interface, but the dynamics are still ultrafast and occur on the picosecond time scale. These interfacial water molecules are not trapped but are mobile in the heterogeneous binding nanospace. Combining our findings with our previous observation of ultrafast water motions at the interface of a loose binding polymerase (Dpo4), we conclude that the binding interface is dynamic and the water molecules in various binding clefts, channels, and caves are mobile and even fluid with different levels of mobility for loose or tight binding polymerases. Such a dynamic interface should be general to all DNA polymerase complexes to ensure the biological function of DNA synthesis.
Collapse
Affiliation(s)
- Yi Yang
- Department of Physics, ‡Department of Chemistry and Biochemistry, and §Programs of Biophysics, Chemical Physics, and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| | | | | | | | | | | | | |
Collapse
|
22
|
Ovchinnikov V, Karplus M. Investigations of α-helix↔β-sheet transition pathways in a miniprotein using the finite-temperature string method. J Chem Phys 2014; 140:175103. [PMID: 24811667 PMCID: PMC4032436 DOI: 10.1063/1.4871685] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2013] [Accepted: 04/07/2014] [Indexed: 12/17/2022] Open
Abstract
A parallel implementation of the finite-temperature string method is described, which takes into account the invariance of coordinates with respect to rigid-body motions. The method is applied to the complex α-helix↔β-sheet transition in a β-hairpin miniprotein in implicit solvent, which exhibits much of the complexity of conformational changes in proteins. Two transition paths are considered, one derived from a linear interpolant between the endpoint structures and the other derived from a targeted dynamics simulation. Two methods for computing the conformational free energy (FE) along the string are compared, a restrained method, and a tessellation method introduced by E. Vanden-Eijnden and M. Venturoli [J. Chem. Phys. 130, 194103 (2009)]. It is found that obtaining meaningful free energy profiles using the present atom-based coordinates requires restricting sampling to a vicinity of the converged path, where the hyperplanar approximation to the isocommittor surface is sufficiently accurate. This sampling restriction can be easily achieved using restraints or constraints. The endpoint FE differences computed from the FE profiles are validated by comparison with previous calculations using a path-independent confinement method. The FE profiles are decomposed into the enthalpic and entropic contributions, and it is shown that the entropy difference contribution can be as large as 10 kcal/mol for intermediate regions along the path, compared to 15-20 kcal/mol for the enthalpy contribution. This result demonstrates that enthalpic barriers for transitions are offset by entropic contributions arising from the existence of different paths across a barrier. The possibility of using systematically coarse-grained representations of amino acids, in the spirit of multiple interaction site residue models, is proposed as a means to avoid ad hoc sampling restrictions to narrow transition tubes.
Collapse
Affiliation(s)
- Victor Ovchinnikov
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Martin Karplus
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
23
|
Liu Y, Radhakrishnan R. Computational delineation of tyrosyl-substrate recognition and catalytic landscapes by the epidermal growth factor receptor tyrosine kinase domain. MOLECULAR BIOSYSTEMS 2014; 10:1890-904. [PMID: 24779031 DOI: 10.1039/c3mb70620f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase (RTK), which catalyzes protein phosphorylation reactions by transferring the γ-phosphoryl group from an ATP molecule to the hydroxyl group of tyrosine residues in protein substrates. EGFR is an important drug target in the treatment of cancers and a better understanding of the receptor function is critical to discern cancer mechanisms. We employ a suite of molecular simulation methods to explore the mechanism of substrate recognition and to delineate the catalytic landscape of the phosphoryl transfer reaction. Based on our results, we propose that a highly conserved region corresponding to Val852-Pro853-Ile854-Lys855-Trp856 in the EGFR tyrosine kinase domain (TKD) is essential for substrate binding. We also provide a possible explanation for the established experimental observation that protein tyrosine kinases (including EGFR) select substrates with a glutamic acid at the P - 1 position and a large hydrophobic amino acid at the P + 1 position. Furthermore, our mixed quantum mechanics/molecular mechanics (QM/MM) simulations show that the EGFR protein kinase favors the dissociative mechanism, although an alternative channel through the formation of an associative transition state is also possible. Our simulations establish some key molecular rules in the operation for substrate-recognition and for phosphoryl transfer in the EGFR TKD.
Collapse
Affiliation(s)
- Yingting Liu
- Department of Bioengineering, University of Pennsylvania, 240 Skirkanich, 210 S. 33rd Street, Philadelphia, PA 19104, USA.
| | | |
Collapse
|
24
|
Li Y, Freudenthal BD, Beard WA, Wilson SH, Schlick T. Optimal and variant metal-ion routes in DNA polymerase β's conformational pathways. J Am Chem Soc 2014; 136:3630-9. [PMID: 24511902 PMCID: PMC7032070 DOI: 10.1021/ja412701f] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To interpret recent structures of the R283K mutant of human DNA repair enzyme DNA polymerase β (pol β) differing in the number of Mg(2+) ions, we apply transition path sampling (TPS) to assess the effect of differing ion placement on the transition from the open one-metal to the closed two-metal state. We find that the closing pathway depends on the initial ion position, both in terms of the individual transition states and associated energies. The energy barrier of the conformational pathway varies from 25 to 58 kJ/mol, compared to the conformational energy barrier of 42 kJ/mol for the wild-type pol β reported previously. Moreover, we find a preferred ion route located in the center of the enzyme, parallel to the DNA. Within this route, the conformational pathway is similar to that of the overall open to closed transition of pol β, but outside it, especially when the ion starts near active site residues Arg258 and Asp190, the conformational pathway diverges significantly. Our findings should apply generally to pol β, since R283K is relatively far from the active site; further experimental and computational work are required to confirm this. Our studies also underscore the common feature that less active mutants have less stable closed states than their open states, in marked contrast to the wild-type enzyme, where the closed state is significantly more stable than the open form.
Collapse
Affiliation(s)
- Yunlang Li
- Department of Chemistry and Courant Institute of
Mathematical Sciences, New York University, 251 Mercer Street, New York, NY
10012
| | - Bret D. Freudenthal
- Laboratory of Structural Biology, National
Institute of Environmental Health Sciences, National Institutes of Health, Bethesda,
Research Triangle Park, NC 27709
| | - William A. Beard
- Laboratory of Structural Biology, National
Institute of Environmental Health Sciences, National Institutes of Health, Bethesda,
Research Triangle Park, NC 27709
| | - Samuel H. Wilson
- Laboratory of Structural Biology, National
Institute of Environmental Health Sciences, National Institutes of Health, Bethesda,
Research Triangle Park, NC 27709
| | - Tamar Schlick
- Department of Chemistry and Courant Institute of
Mathematical Sciences, New York University, 251 Mercer Street, New York, NY
10012
| |
Collapse
|
25
|
Moors SLC, Herdewijn P, Robben J, Ceulemans A. Cooperative dynamics of a DNA polymerase replicating complex. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:2554-63. [PMID: 24041502 DOI: 10.1016/j.bbapap.2013.09.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Received: 04/25/2013] [Revised: 08/22/2013] [Accepted: 09/06/2013] [Indexed: 11/15/2022]
Abstract
Engineered DNA polymerases continue to be the workhorses of many applications in biotechnology, medicine and nanotechnology. However, the dynamic interplay between the enzyme and the DNA remains unclear. In this study, we performed an extensive replica exchange with flexible tempering (REFT) molecular dynamics simulation of the ternary replicating complex of the archaeal family B DNA polymerase from the thermophile Thermococcus gorgonarius, right before the chemical step. The convoluted dynamics of the enzyme are reducible to rigid-body motions of six subdomains. Upon binding to the enzyme, the DNA double helix conformation changes from a twisted state to a partially untwisted state. The twisted state displays strong bending motion, whereby the DNA oscillates between a straight and a bent conformation. The dynamics of double-stranded DNA are strongly correlated with rotations of the thumb toward the palm, which suggests an assisting role of the enzyme during DNA translocation. In the complex, the primer-template duplex displays increased preference for the B-DNA conformation at the n-2 and n-3 dinucleotide steps. Interactions at the primer 3' end indicate that Thr541 and Asp540 are the acceptors of the first proton transfer in the chemical step, whereas in the translocation step both residues hold the primer 3' terminus in the vicinity of the priming site, which is crucial for high processivity.
Collapse
Affiliation(s)
- Samuel L C Moors
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium.
| | | | | | | |
Collapse
|
26
|
Maláč K, Barvík I. Substrate recognition by norovirus polymerase: microsecond molecular dynamics study. J Comput Aided Mol Des 2013; 27:373-88. [PMID: 23619980 DOI: 10.1007/s10822-013-9652-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2012] [Accepted: 04/18/2013] [Indexed: 01/02/2023]
Abstract
Molecular dynamics simulations of complexes between Norwalk virus RNA dependent RNA polymerase and its natural CTP and 2dCTP (both containing the O5'-C5'-C4'-O4' sequence of atoms bridging the triphosphate and sugar moiety) or modified coCTP (C5'-O5'-C4'-O4'), cocCTP (C5'-O5'-C4'-C4'') substrates were produced by means of CUDA programmable graphical processing units and the ACEMD software package. It enabled us to gain microsecond MD trajectories clearly showing that similar nucleoside triphosphates can bind surprisingly differently into the active site of the Norwalk virus RNA dependent RNA polymerase. It corresponds to their different modes of action (CTP-substrate, 2dCTP-poor substrate, coCTP-chain terminator, cocCTP-inhibitor). Moreover, extremely rare events-as repetitive pervasion of Arg182 into a potentially reaction promoting arrangement-were captured.
Collapse
Affiliation(s)
- Kamil Maláč
- Faculty of Mathematics and Physics, Institute of Physics, Charles University, Ke Karlovu 5, Prague 2, 121 16, Czech Republic
| | | |
Collapse
|
27
|
Ayyaswamy PS, Muzykantov V, Eckmann DM, Radhakrishnan R. Nanocarrier Hydrodynamics and Binding in Targeted Drug Delivery: Challenges in Numerical Modeling and Experimental Validation. J Nanotechnol Eng Med 2013; 4:101011-1010115. [PMID: 23917383 PMCID: PMC3708709 DOI: 10.1115/1.4024004] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2012] [Revised: 03/06/2013] [Indexed: 11/08/2022]
Abstract
This review discusses current progress and future challenges in the numerical modeling of targeted drug delivery using functionalized nanocarriers (NC). Antibody coated nanocarriers of various size and shapes, also called functionalized nanocarriers, are designed to be injected in the vasculature, whereby they undergo translational and rotational motion governed by hydrodynamic interaction with blood particulates as well as adhesive interactions mediated by the surface antibody binding to target antigens/receptors on cell surfaces. We review current multiscale modeling approaches rooted in computational fluid dynamics and nonequilibrium statistical mechanics to accurately resolve fluid, thermal, as well as adhesive interactions governing nanocarrier motion and their binding to endothelial cells lining the vasculature. We also outline current challenges and unresolved issues surrounding the modeling methods. Experimental approaches in pharmacology and bioengineering are discussed briefly from the perspective of model validation.
Collapse
Affiliation(s)
- Portonovo S. Ayyaswamy
- Department of Mechanical Engineering and Applied Mechanics,University of Pennsylvania,Philadelphia, PA 19104
| | - Vladimir Muzykantov
- Department of Pharmacology,and Center for Targeted Therapeutics and Translational Nanomedicine,University of Pennsylvania,Philadelphia, PA 19104
| | - David M. Eckmann
- Institute of Translational Medicine and Therapeutics,Department of Anesthesiology and Critical Care,and Department of Bioengineering,University of Pennsylvania,Philadelphia, PA 19104
| | - Ravi Radhakrishnan
- Institute of Translational Medicine and Therapeutics,Department of Bioengineering,Department of Chemical and Biomolecular Engineering,University of Pennsylvania,Philadelphia, PA 19104e-mail:
| |
Collapse
|
28
|
Abstract
Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.
Collapse
|
29
|
Arora K, Brooks CL. Multiple intermediates, diverse conformations, and cooperative conformational changes underlie the catalytic hydride transfer reaction of dihydrofolate reductase. Top Curr Chem (Cham) 2013; 337:165-87. [PMID: 23420416 PMCID: PMC4394636 DOI: 10.1007/128_2012_408] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
It has become increasingly clear that protein motions play an essential role in enzyme catalysis. However, exactly how these motions are related to an enzyme's chemical step is still intensely debated. This chapter examines the possible role of protein motions that display a hierarchy of timescales in enzyme catalysis. The linkage between protein motions and catalysis is investigated in the context of a model enzyme, E. coli dihydrofolate reductase (DHFR), that catalyzes the hydride transfer reaction in the conversion of dihydrofolate to tetrahydrofolate. The results of extensive computer simulations probing the protein motions that are manifest during different steps along the turnover cycle of DHFR are summarized. Evidence is presented that the protein motions modulate the catalytic efficacy of DHFR by generating a conformational ensemble conducive to the hydride transfer. The alteration of the equilibrium conformational ensemble rather than any protein dynamical effects is found to be sufficient to explain the rate-diminishing effects of mutation on the kinetics of the enzyme. These data support the view that the protein motions facilitate catalysis by establishing reaction competent conformations of the enzyme, but they do not directly couple to the chemical reaction itself. These findings have broad implications for our understanding of enzyme mechanisms and the design of novel protein catalysts.
Collapse
Affiliation(s)
- Karunesh Arora
- Department of Chemistry and Biophysics Program, University of Michigan, Ann Arbor, MI 48109
| | - Charles L. Brooks
- Department of Chemistry and Biophysics Program, University of Michigan, Ann Arbor, MI 48109
| |
Collapse
|
30
|
Schlick T, Arora K, Beard WA, Wilson SH. Perspective: pre-chemistry conformational changes in DNA polymerase mechanisms. Theor Chem Acc 2012; 131:1287. [PMID: 23459563 DOI: 10.1007/s00214-012-1287-7] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
In recent papers, there has been a lively exchange concerning theories for enzyme catalysis, especially the role of protein dynamics/pre-chemistry conformational changes in the catalytic cycle of enzymes. Of particular interest is the notion that substrate-induced conformational changes that assemble the polymerase active site prior to chemistry are required for DNA synthesis and impact fidelity (i.e., substrate specificity). High-resolution crystal structures of DNA polymerase β representing intermediates of substrate complexes prior to the chemical step are available. These structures indicate that conformational adjustments in both the protein and substrates must occur to achieve the requisite geometry of the reactive participants for catalysis. We discuss computational and kinetic methods to examine possible conformational change pathways that lead from the observed crystal structure intermediates to the final structures poised for chemistry. The results, as well as kinetic data from site-directed mutagenesis studies, are consistent with models requiring pre-chemistry conformational adjustments in order to achieve high fidelity DNA synthesis. Thus, substrate-induced conformational changes that assemble the polymerase active site prior to chemistry contribute to DNA synthesis even when they do not represent actual rate-determining steps for chemistry.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, 100 Washington Square East, Silver Building, New York, NY 10003, USA. Courant Institute of Mathematical Sciences, New York, University, 251 Mercer Street, New York, NY 10012, USA
| | | | | | | |
Collapse
|
31
|
Ram Prasad B, Kamerlin SCL, Florián J, Warshel A. Prechemistry barriers and checkpoints do not contribute to fidelity and catalysis as long as they are not rate limiting. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1288-6] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
32
|
Klvaňa M, Murphy DL, Jeřábek P, Goodman MF, Warshel A, Sweasy JB, Florián J. Catalytic effects of mutations of distant protein residues in human DNA polymerase β: theory and experiment. Biochemistry 2012; 51:8829-43. [PMID: 23013478 DOI: 10.1021/bi300783t] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We carried out free-energy calculations and transient kinetic experiments for the insertion of the right (dC) and wrong (dA) nucleotides by wild-type (WT) and six mutant variants of human DNA polymerase β (Pol β). Since the mutated residues in the point mutants, I174S, I260Q, M282L, H285D, E288K, and K289M, were not located in the Pol β catalytic site, we assumed that the WT and its point mutants share the same dianionic phosphorane transition-state structure of the triphosphate moiety of deoxyribonucleotide 5'-triphosphate (dNTP) substrate. On the basis of this assumption, we have formulated a thermodynamic cycle for calculating relative dNTP insertion efficiencies, Ω = (k(pol)/K(D))(mut)/(k(pol)/K(D))(WT) using free-energy perturbation (FEP) and linear interaction energy (LIE) methods. Kinetic studies on five of the mutants have been published previously using different experimental conditions, e.g., primer-template sequences. We have performed a presteady kinetic analysis for the six mutants for comparison with wild-type Pol β using the same conditions, including the same primer/template DNA sequence proximal to the dNTP insertion site used for X-ray crystallographic studies. This consistent set of kinetic and structural data allowed us to eliminate the DNA sequence from the list of factors that can adversely affect calculated Ω values. The calculations using the FEP free energies scaled by 0.5 yielded 0.9 and 1.1 standard deviations from the experimental log Ω values for the insertion of the right and wrong dNTP, respectively. We examined a hybrid FEP/LIE method in which the FEP van der Waals term for the interaction of the mutated amino acid residue with its surrounding environment was replaced by the corresponding van der Waals term calculated using the LIE method, resulting in improved 0.4 and 1.0 standard deviations from the experimental log Ω values. These scaled FEP and FEP/LIE methods were also used to predict log Ω for R283A and R283L Pol β mutants.
Collapse
Affiliation(s)
- Martin Klvaňa
- Department of Chemistry, Loyola University, Chicago, Illinois 60626, United States
| | | | | | | | | | | | | |
Collapse
|
33
|
Ovchinnikov V, Karplus M. Analysis and elimination of a bias in targeted molecular dynamics simulations of conformational transitions: application to calmodulin. J Phys Chem B 2012; 116:8584-603. [PMID: 22409258 PMCID: PMC3406239 DOI: 10.1021/jp212634z] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The popular targeted molecular dynamics (TMD) method for generating transition paths in complex biomolecular systems is revisited. In a typical TMD transition path, the large-scale changes occur early and the small-scale changes tend to occur later. As a result, the order of events in the computed paths depends on the direction in which the simulations are performed. To identify the origin of this bias, and to propose a method in which the bias is absent, variants of TMD in the restraint formulation are introduced and applied to the complex open ↔ closed transition in the protein calmodulin. Due to the global best-fit rotation that is typically part of the TMD method, the simulated system is guided implicitly along the lowest-frequency normal modes, until the large spatial scales associated with these modes are near the target conformation. The remaining portion of the transition is described progressively by higher-frequency modes, which correspond to smaller-scale rearrangements. A straightforward modification of TMD that avoids the global best-fit rotation is the locally restrained TMD (LRTMD) method, in which the biasing potential is constructed from a number of TMD potentials, each acting on a small connected portion of the protein sequence. With a uniform distribution of these elements, transition paths that lack the length-scale bias are obtained. Trajectories generated by steered MD in dihedral angle space (DSMD), a method that avoids best-fit rotations altogether, also lack the length-scale bias. To examine the importance of the paths generated by TMD, LRTMD, and DSMD in the actual transition, we use the finite-temperature string method to compute the free energy profile associated with a transition tube around a path generated by each algorithm. The free energy barriers associated with the paths are comparable, suggesting that transitions can occur along each route with similar probabilities. This result indicates that a broad ensemble of paths needs to be calculated to obtain a full description of conformational changes in biomolecules. The breadth of the contributing ensemble suggests that energetic barriers for conformational transitions in proteins are offset by entropic contributions that arise from a large number of possible paths.
Collapse
|
34
|
|
35
|
Li Y, Gridley CL, Jaeger J, Sweasy JB, Schlick T. Unfavorable electrostatic and steric interactions in DNA polymerase β E295K mutant interfere with the enzyme's pathway. J Am Chem Soc 2012; 134:9999-10010. [PMID: 22651551 PMCID: PMC3383778 DOI: 10.1021/ja300361r] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Mutations in DNA polymerase β (pol β) have been associated with approximately 30% of human tumors. The E295K mutation of pol β has been linked to gastric carcinoma via interference with base excision repair. To interpret the different behavior of E295K as compared to wild-type pol β in atomic and energetic detail, we resolve a binary crystal complex of E295K at 2.5 Å and apply transition path sampling (TPS) to delineate the closing pathway of the E295K pol β mutant. Conformational changes are important components in the enzymatic pathway that lead to and ready the enzyme for the chemical reaction. Our analyses show that the closing pathway of E295K mutant differs from the wild-type pol β in terms of the individual transition states along the pathway, associated energies, and the active site conformation in the final closed form of the mutant. In particular, the closed state of E295K has a more distorted active site than the active site in the wild-type pol β. In addition, the total energy barrier in the conformational closing pathway is 65 ± 11 kJ/mol, much higher than that estimated for both correct (e.g., G:C) and incorrect (e.g., G:A) wild-type pol β systems (42 ± 8 and 45 ± 7 kJ/mol, respectively). In particular, the rotation of Arg258 is the rate-limiting step in the conformational pathway of E295K due to unfavorable electrostatic and steric interactions. The distorted active site in the closed relative to open state and the high energy barrier in the conformational pathway may explain in part why the E295K mutant is observed to be inactive. Interestingly, however, following the closing of the thumb but prior to the rotation of Arg258, the E295K mutant complex has a similar energy level as compared to the wild-type pol β. This suggests that the E295K mutant may associate with DNA with similar affinity, but it may be hampered in continuing the process of chemistry. Supporting experimental data come from the observation that the catalytic activity of wild-type pol β is hampered when E295K is present: this may arise from the competition between E295K and wild-type enzyme for the DNA. These combined results suggest that the low insertion efficiency of E295K mutant as compared to wild-type pol β may be related to a closed form distorted by unfavorable electrostatic and steric interactions between Arg258 and other key residues. The active site is thus less competent for proceeding to the chemical reaction, which may also involve a higher reaction barrier than the wild-type or may not be possible in this mutant. Our analysis also suggests further experiments for other mutants to test the above hypothesis and dissect the roles of steric and electrostatic factors on enzyme behavior.
Collapse
Affiliation(s)
- Yunlang Li
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012
| | - Chelsea L. Gridley
- Department of Biomedical Sciences, School of Public Health, University at Albany, 1400 Washington Avenue, Albany, NY 12222, USA
| | - Joachim Jaeger
- Department of Biomedical Sciences, School of Public Health, University at Albany, 1400 Washington Avenue, Albany, NY 12222, USA
- Division of Genetics, Wadsworth Center NYS-DOH, New Scotland Avenue, Albany, NY 12208, USA
| | - Joann B. Sweasy
- Department of Therapeutic Radiology, Yale University School of Medicine, 333 Cedar Street, P.O. Box 208040, New Haven, CT 06520, USA
| | - Tamar Schlick
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012
| |
Collapse
|
36
|
Kirmizialtin S, Nguyen V, Johnson KA, Elber R. How conformational dynamics of DNA polymerase select correct substrates: experiments and simulations. Structure 2012; 20:618-27. [PMID: 22483109 PMCID: PMC3322391 DOI: 10.1016/j.str.2012.02.018] [Citation(s) in RCA: 100] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2012] [Revised: 02/21/2012] [Accepted: 02/23/2012] [Indexed: 11/26/2022]
Abstract
Nearly every enzyme undergoes a significant change in structure after binding it's substrate. Experimental and theoretical analyses of the role of changes in HIV reverse transcriptase structure in selecting a correct substrate are presented. Atomically detailed simulations using the Milestoning method predict a rate and free energy profile of the conformational change commensurate with experimental data. A large conformational change occurring on a millisecond timescale locks the correct nucleotide at the active site but promotes release of a mismatched nucleotide. The positions along the reaction coordinate that decide the yield of the reaction are not determined by the chemical step. Rather, the initial steps of weak substrate binding and protein conformational transition significantly enrich the yield of a reaction with a correct substrate, whereas the same steps diminish the reaction probability of an incorrect substrate.
Collapse
Affiliation(s)
- Serdal Kirmizialtin
- Department of Chemistry and Biochemistry and Institute for Computational Engineering and Sciences (ICES), The University of Texas at Austin, 1 University Station, Austin, Texas 78712
| | - Virginia Nguyen
- Institute for Cellular and Molecular Biology, The University of Texas at Austin, 2500 Speedway, MBB 3.122, Austin, Texas 78712, USA
| | - Kenneth A. Johnson
- Institute for Cellular and Molecular Biology, The University of Texas at Austin, 2500 Speedway, MBB 3.122, Austin, Texas 78712, USA
| | - Ron Elber
- Department of Chemistry and Biochemistry and Institute for Computational Engineering and Sciences (ICES), The University of Texas at Austin, 1 University Station, Austin, Texas 78712
| |
Collapse
|
37
|
Balbo PB, Wang ECW, Tsai MD. Kinetic mechanism of active site assembly and chemical catalysis of DNA polymerase β. Biochemistry 2011; 50:9865-75. [PMID: 22010960 DOI: 10.1021/bi200954r] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
It has been inferred from structural and computational studies that the mechanism of DNA polymerases involves subtle but important discrete steps that occur between binding and recognition of the correct dNTP and chemical catalysis. These steps potentially include local conformational changes involving active site residues, reorganization of Mg(2+)-coordinating ligands, and proton transfer. Here we address this broad issue by conducting extensive transient state kinetic analyses of DNA polymerase β (Pol β). We also performed kinetic simulations to evaluate alternative kinetic models. These studies provide some support for two-step subdomain closing and define constraints under which a kinetically significant prechemistry step can occur. To experimentally identify additional microscopic steps, we developed a stopped flow absorbance assay to measure proton formation that occurs during catalysis. These studies provide direct evidence that formation of the enzyme-bound 3'-O(-) nucleophile is rate determining for chemistry. We additionally show that at low pH the chemical step is rate limiting for catalysis, but at high pH, a postchemistry conformational step is rate limiting due to a pH-dependent increase in the rate of nucleotidyl transfer. Finally, we performed exhaustive analyses of [Mg(2+)] and pH effects. In contrast to published studies, the results suggest an irregular pH dependence of k(pol), which is consistent with general base catalysis involving cooperativity between two or more protonic residues. Overall, the results represent significant advancement in the kinetic mechanism of Pol β and also reconcile some computational and experimental findings.
Collapse
Affiliation(s)
- Paul B Balbo
- Institute of Biological Chemistry, Academia Sinica, Nangang, Taipei 115, Taiwan
| | | | | |
Collapse
|
38
|
Ram Prasad B, Warshel A. Prechemistry versus preorganization in DNA replication fidelity. Proteins 2011; 79:2900-19. [PMID: 21905114 DOI: 10.1002/prot.23128] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2011] [Revised: 06/30/2011] [Accepted: 07/05/2011] [Indexed: 01/30/2023]
Abstract
The molecular origin of nucleotide insertion catalysis and fidelity of DNA polymerases is explored by means of computational simulations. Special attention is paid to the examination of the validity of proposals that invoke prechemistry effects, checkpoints concepts, and dynamical effects. The simulations reproduce the observed fidelity in Pol β, starting with the relevant observed X-ray structures of the complex with the right (R) and wrong (W) nucleotides. The generation of free energy surfaces for the R and W systems also allowed us to analyze different proposals about the origin of the fidelity and to reach several important conclusions. It is found that the potential of mean force (PMF) obtained by proper sampling does not support QM/MM-based proposals of a large barrier before the prechemistry state. Furthermore, examination of dynamical proposals by the renormalization approach indicates that the motions from open to close configurations do not contribute to catalysis or fidelity. Finally we discuss and analyze the induced fit concept and show that, despite its importance, it does not explain fidelity. That is, the fidelity is apparently due to the change in the preorganization of the chemical site, as a result of the relaxation of the binding site upon binding of the incorrect nucleotide. Finally and importantly, since the issue is the barrier associated with the enzyme-substrate (ES)/DNA complex at the chemical transition state and not the path to this complex formation (unless this path involves rate determining steps), it is also not useful to invoke checkpoints while discussing fidelity.
Collapse
Affiliation(s)
- B Ram Prasad
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | |
Collapse
|
39
|
Shah M, Santiso EE, Trout BL. Computer Simulations of Homogeneous Nucleation of Benzene from the Melt. J Phys Chem B 2011; 115:10400-12. [DOI: 10.1021/jp203550t] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Manas Shah
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Erik E. Santiso
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Bernhardt L. Trout
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
40
|
Fuchigami S, Fujisaki H, Matsunaga Y, Kidera A. Protein Functional Motions: Basic Concepts and Computational Methodologies. ADVANCING THEORY FOR KINETICS AND DYNAMICS OF COMPLEX, MANY-DIMENSIONAL SYSTEMS: CLUSTERS AND PROTEINS 2011. [DOI: 10.1002/9781118087817.ch2] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
41
|
Schlick T, Collepardo-Guevara R, Halvorsen LA, Jung S, Xiao X. Biomolecularmodeling and simulation: a field coming of age. Q Rev Biophys 2011; 44:191-228. [PMID: 21226976 PMCID: PMC3700731 DOI: 10.1017/s0033583510000284] [Citation(s) in RCA: 112] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We assess the progress in biomolecular modeling and simulation, focusing on structure prediction and dynamics, by presenting the field’s history, metrics for its rise in popularity, early expressed expectations, and current significant applications. The increases in computational power combined with improvements in algorithms and force fields have led to considerable success, especially in protein folding, specificity of ligand/biomolecule interactions, and interpretation of complex experimental phenomena (e.g. NMR relaxation, protein-folding kinetics and multiple conformational states) through the generation of structural hypotheses and pathway mechanisms. Although far from a general automated tool, structure prediction is notable for proteins and RNA that preceded the experiment, especially by knowledge-based approaches. Thus, despite early unrealistic expectations and the realization that computer technology alone will not quickly bridge the gap between experimental and theoretical time frames, ongoing improvements to enhance the accuracy and scope of modeling and simulation are propelling the field onto a productive trajectory to become full partner with experiment and a field on its own right.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, 100 Washington Square East, Silver Building, New York, NY 10003, USA.
| | | | | | | | | |
Collapse
|
42
|
Xie P. A model for the dynamics of mammalian family X DNA polymerases. J Theor Biol 2011; 277:111-22. [PMID: 21377475 DOI: 10.1016/j.jtbi.2011.02.018] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2010] [Revised: 12/22/2010] [Accepted: 02/22/2011] [Indexed: 11/28/2022]
Abstract
Based on available structural studies, a model is presented for polymerization dynamics of mammalian family X DNA polymerases, including polymerases β, λ, μ, and terminal deoxynucleotidyl transferase (TdT). Using the model, distinct polymerization activities and processivities of the four polymerases acting on different forms of DNA substrate are analyzed and studied theoretically. A "gradient" of template dependence of polymerases β, λ, μ, and TdT is well explained. The much higher occurrence frequencies of the -1 frameshift DNA synthesis by pols λ and μ than that by pol β are well explained. The theoretical results on the polymerization processivities are also in agreement with the available experimental data.
Collapse
Affiliation(s)
- Ping Xie
- Key Laboratory of Soft Matter Physics and Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China.
| |
Collapse
|
43
|
Li Y, Schlick T. Modeling DNA polymerase μ motions: subtle transitions before chemistry. Biophys J 2011; 99:3463-72. [PMID: 21081096 DOI: 10.1016/j.bpj.2010.09.056] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2010] [Revised: 09/24/2010] [Accepted: 09/28/2010] [Indexed: 11/29/2022] Open
Abstract
To investigate whether an open-to-closed transition before the chemical step and induced-fit mechanism exist in DNA polymerase μ (pol μ), we analyze a series of molecular-dynamics simulations with and without the incoming nucleotide in various forms, including mutant systems, based on pol μ's crystal ternary structure. Our simulations capture no significant large-scale motion in either the DNA or the protein domains of pol μ. However, subtle residue motions can be distinguished, specifically of His(329) and Asp(330) to assemble in pol μ's active site, and of Gln(440) and Glu(443) to help accommodate the incoming nucleotide. Mutant simulations capture a DNA frameshift pairing and indicate the importance of Arg(444) and Arg(447) in stacking with the DNA template, and of Arg(448) and Gln(440) in helping to stabilize the position of both the DNA template and the incoming nucleotide. Although limited sampling in the molecular-dynamics simulations cannot be ruled out, our studies suggest an absence of a large-scale motion in pol μ. Together with the known crystallization difficulties of capturing the open form of pol μ, our studies also raise the possibility that a well-defined open form may not exist. Moreover, we suggest that residues Arg(448) and Gln(440) may be crucial for preventing insertion frameshift errors in pol μ.
Collapse
Affiliation(s)
- Yunlang Li
- Department of Chemistry, New York University, New York, NY, USA
| | | |
Collapse
|
44
|
Foley MC, Padow VA, Schlick T. DNA pol λ's extraordinary ability to stabilize misaligned DNA. J Am Chem Soc 2010; 132:13403-16. [PMID: 20822183 DOI: 10.1021/ja1049687] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
DNA polymerases have the venerable task of maintaining genome stability during DNA replication and repair. Errors, nonetheless, occur with error propensities that are polymerase specific. For example, DNA polymerase λ (pol λ) generates single-base deletions through template-strand slippage within short repetitive DNA regions much more readily than does the closely related polymerase β (pol β). Here we present in silico evidence to help interpret pol λ's greater tendency for deletion errors than pol β by its more favorable protein/DNA electrostatic interactions immediately around the extrahelical nucleotide on the template strand. Our molecular dynamics and free energy analyses suggest that pol λ provides greater stabilization to misaligned DNA than aligned DNA. Our study of several pol λ mutants of Lys544 (Ala, Phe, Glu) probes the interactions between the extrahelical nucleotide and the adjacent Lys544 to show that the charge of the 544 residue controls stabilization of the DNA misalignment. In addition, we identify other thumb residues (Arg538, Lys521, Arg517, and Arg514) that play coordinating roles in stabilizing pol λ's interactions with misaligned DNA. Interestingly, their aggregate stabilization effect is more important than that of any one component residue, in contrast to aligned DNA systems, as we determined from mutations of these key residues and energetic analyses. No such comparable network of stabilizing misaligned DNA exists in pol β. Evolutionary needs for DNA repair on substrates with minimal base-pairing, such as those encountered by pol λ in the non-homologous end-joining pathway, may have been solved by a greater tolerance to deletion errors. Other base-flipping proteins share similar binding properties and motions for extrahelical nucleotides.
Collapse
Affiliation(s)
- Meredith C Foley
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, New York 10012, USA
| | | | | |
Collapse
|
45
|
Betz K, Streckenbach F, Schnur A, Exner T, Welte W, Diederichs K, Marx A. Structures of DNA polymerases caught processing size-augmented nucleotide probes. Angew Chem Int Ed Engl 2010; 49:5181-4. [PMID: 20572212 DOI: 10.1002/anie.200905724] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Karin Betz
- Department of Chemistry, Konstanz Research School Chemical Biology, Universität Konstanz, Universitätsstrasse 10, 78457 Konstanz, Germany
| | | | | | | | | | | | | |
Collapse
|
46
|
Bhatt D, Zuckerman DM. Heterogeneous path ensembles for conformational transitions in semi-atomistic models of adenylate kinase. J Chem Theory Comput 2010; 6:3527-3539. [PMID: 21660120 PMCID: PMC3108504 DOI: 10.1021/ct100406t] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
We performed "weighted ensemble" path-sampling simulations of adenylate kinase, using several semi-atomistic protein models. The models have an all-atom backbone with various levels of residue interactions. The primary result is that full statistically rigorous path sampling required only a few weeks of single-processor computing time with these models, indicating the addition of further chemical detail should be readily feasible. Our semi-atomistic path ensembles are consistent with previous biophysical findings: the presence of two distinct pathways, identification of intermediates, and symmetry of forward and reverse pathways.
Collapse
Affiliation(s)
- Divesh Bhatt
- Department of Computational and Systems Biology, University of Pittsburgh
| | | |
Collapse
|
47
|
Betz K, Streckenbach F, Schnur A, Exner T, Welte W, Diederichs K, Marx A. Strukturen von DNA-Polymerasen mit 4′-alkylierten Nucleotiden im aktiven Zentrum. Angew Chem Int Ed Engl 2010. [DOI: 10.1002/ange.200905724] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
48
|
Foley MC, Schlick T. Relationship between conformational changes in pol lambda's active site upon binding incorrect nucleotides and mismatch incorporation rates. J Phys Chem B 2010; 113:13035-47. [PMID: 19572669 DOI: 10.1021/jp903172x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The correct replication and repair of DNA is critical for a cell's survival. Here, we investigate the fidelity of mammalian DNA polymerase lambda (pol lambda) utilizing dynamics simulation of the enzyme bound to incorrect incoming nucleotides including A:C, A:G, A(syn):G, A:A, A(syn):A, and T:G, all of which exhibit differing incorporation rates for pol lambda as compared to A:T bound to pol lambda. The wide range of DNA motion and protein residue side-chain motions observed in the mismatched systems demonstrates distinct differences when compared to the reference (correct base pair) system. Notably, Arg517's interactions with the DNA template strand bases in the active site are more limited, and Arg517 displays increased interactions with the incorrect dNTPs. This effect suggests that Arg517 helps provide a base-checking mechanism to discriminate correct from incorrect dNTPs. In addition, we find Tyr505 and Phe506 also play key roles in this base checking. A survey of the electrostatic potential landscape of the active sites and concomitant changes in electrostatic interaction energy between Arg517 and the dNTPs reveals that pol lambda binds incorrect dNTPs less tightly than the correct dNTP. These trends lead us to propose the following order for mismatch insertion by pol lambda: A:C > A:G > A(syn):G > T:G > A(syn):A > A:A. This sequence agrees with available kinetic data for incorrect nucleotide insertion opposite template adenine, with the exception of T:G, which may be more sensitive to the insertion context.
Collapse
Affiliation(s)
- Meredith C Foley
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, New York 10012, USA
| | | |
Collapse
|
49
|
Golosov AA, Warren JJ, Beese LS, Karplus M. The mechanism of the translocation step in DNA replication by DNA polymerase I: a computer simulation analysis. Structure 2010; 18:83-93. [PMID: 20152155 PMCID: PMC3325112 DOI: 10.1016/j.str.2009.10.014] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2009] [Revised: 10/20/2009] [Accepted: 10/22/2009] [Indexed: 11/17/2022]
Abstract
High-fidelity DNA polymerases copy DNA rapidly and accurately by adding correct deoxynucleotide triphosphates to a growing primer strand of DNA. Following nucleotide incorporation, a series of conformational changes translocate the DNA substrate by one base pair step, readying the polymerase for the next round of incorporation. Molecular dynamics simulations indicate that the translocation consists globally of a polymerase fingers-opening transition, followed by the DNA displacement and the insertion of the template base into the preinsertion site. They also show that the pyrophosphate release facilitates the opening transition and that the universally conserved Y714 plays a key role in coupling polymerase opening to DNA translocation. The transition involves several metastable intermediates in one of which the O helix is bent in the vicinity of G711. Completion of the translocation appears to require a gating motion of the O1 helix, perhaps facilitated by the presence of G715. These roles are consistent with the high level of conservation of Y714 and the two glycine residues at these positions. It is likely that a corresponding mechanism is applicable to other polymerases.
Collapse
Affiliation(s)
- Andrei A. Golosov
- Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford St., Cambridge, MA 02138
| | - Joshua J. Warren
- Department of Biochemistry, Duke University Medical Center, Box 3711, Durham, NC 27710
| | - Lorena S. Beese
- Department of Biochemistry, Duke University Medical Center, Box 3711, Durham, NC 27710
| | - Martin Karplus
- Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford St., Cambridge, MA 02138
- Institut de Science et d'Ingénierie Supramoléculaires Université Louis Pasteur, 67000 Strasbourg, France
| |
Collapse
|
50
|
Woo HJ, Jiang J, Lafer EM, Sousa R. ATP-induced conformational changes in Hsp70: molecular dynamics and experimental validation of an in silico predicted conformation. Biochemistry 2009; 48:11470-7. [PMID: 19883127 DOI: 10.1021/bi901256y] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The 70 kDa heat shock proteins (Hsp70s) play important roles in preventing the misfolding of proteins and repairing damage under stress by coupling ATP binding and hydrolysis to protein substrate release and binding, respectively. ATP binding is believed to induce closing of the Hsp70 nucleotide binding domain (NBD) around the nucleotide. We report here a combined computational-experimental study of this open-closed transition. All-atom molecular dynamics simulations were performed for isolated open state NBDs with and without bound ATP. The nucleotide-free NBD samples a wide range of open configurations exhibiting flexible rearrangements of its four subdomains (IA-IIB). In contrast, the ATP-bound Hsp70 NBD closes to a range of configurations that is substantially more closed than the conformation observed in crystals of ATP-complexed NBDs. The close approach of subdomains IB and IIB observed in the simulations results in a strong coordination of the fluorescence probe Trp90 of IB with Arg261 of IIB, a feature not seen in the crystal structures. To determine if this computationally observed conformation occurs in solution, we constructed an R261A mutant. The mutation was found to increase the K(m) and k(cat) for ATP and to significantly reduce the extent of the fluorescence quench observed upon ATP binding. Our results thus account for the previously unexplained ATP-driven change in Trp90 fluorescence seen in the isolated NBD.
Collapse
Affiliation(s)
- Hyung-June Woo
- Department of Chemistry, University of Nevada, Reno, Nevada 89557-0216, USA.
| | | | | | | |
Collapse
|