1
|
Pan X, Van R, Pu J, Nam K, Mao Y, Shao Y. Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions. J Chem Theory Comput 2023; 19:8234-8244. [PMID: 37943896 PMCID: PMC10835707 DOI: 10.1021/acs.jctc.3c00973] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2023]
Abstract
In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.
Collapse
Affiliation(s)
- Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Richard Van
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
- Laboratory of Computational Biology, National, Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| |
Collapse
|
2
|
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
|
3
|
Furlan V, Bren U. Insight into Inhibitory Mechanism of PDE4D by Dietary Polyphenols Using Molecular Dynamics Simulations and Free Energy Calculations. Biomolecules 2021; 11:479. [PMID: 33806914 PMCID: PMC8004924 DOI: 10.3390/biom11030479] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2021] [Revised: 03/19/2021] [Accepted: 03/21/2021] [Indexed: 12/11/2022] Open
Abstract
Phosphodiesterase 4 (PDE4), mainly present in immune, epithelial, and brain cells, represents a family of key enzymes for the degradation of cyclic adenosine monophosphate (cAMP), which modulates inflammatory response. In recent years, the inhibition of PDE4 has been proven to be an effective therapeutic strategy for the treatment of neurological disorders. PDE4D constitutes a high-interest therapeutic target primarily for the treatment of Alzheimer's disease, as it is highly involved in neuroinflammation, learning ability, and memory dysfunctions. In the present study, a thorough computational investigation consisting of molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations based on the linear response approximation (LRA) method was performed to study dietary polyphenols as potential PDE4D inhibitors. The obtained results revealed that curcumin, 6-gingerol, capsaicin, and resveratrol represent potential PDE4D inhibitors; however, the predicted binding free energies of 6-gingerol, capsaicin, and resveratrol were less negative than in the case of curcumin, which exhibited the highest inhibitory potency in comparison with a positive control rolipram. Our results also revealed that the electrostatic component through hydrogen bonding represents the main driving force for the binding and inhibitory activity of curcumin, 6-gingerol, and resveratrol, while the van der Waals component through shape complementarity plays the most important role in capsaicin's inhibitory activity. All investigated compounds form hydrophobic interactions with residues Gln376 and Asn602 as well as hydrogen bonds with nearby residues Asp438, Met439, and Ser440. The binding mode of the studied natural compounds is consequently very similar; however, it significantly differs from the binding of known PDE4 inhibitors. The uncovered molecular inhibitory mechanisms of four investigated natural polyphenols, curcumin, 6-gingerol, capsaicin, and resveratrol, form the basis for the design of novel PDE4D inhibitors for the treatment of Alzheimer's disease with a potentially wider therapeutic window and fewer adverse side effects.
Collapse
Affiliation(s)
- Veronika Furlan
- Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova 17, SI-2000 Maribor, Slovenia;
| | - Urban Bren
- Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova 17, SI-2000 Maribor, Slovenia;
- Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, SI-6000 Koper, Slovenia
| |
Collapse
|
4
|
|
5
|
Walker AR, Cisneros GA. Computational Simulations of DNA Polymerases: Detailed Insights on Structure/Function/Mechanism from Native Proteins to Cancer Variants. Chem Res Toxicol 2017; 30:1922-1935. [PMID: 28877429 PMCID: PMC5696005 DOI: 10.1021/acs.chemrestox.7b00161] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
![]()
Genetic information is vital in the
cell cycle of DNA-based organisms.
DNA polymerases (DNA Pols) are crucial players in transactions dealing
with these processes. Therefore, the detailed understanding of the
structure, function, and mechanism of these proteins has been the
focus of significant effort. Computational simulations have been applied
to investigate various facets of DNA polymerase structure and function.
These simulations have provided significant insights over the years.
This perspective presents the results of various computational studies
that have been employed to research different aspects of DNA polymerases
including detailed reaction mechanism investigation, mutagenicity
of different metal cations, possible factors for fidelity synthesis,
and discovery/functional characterization of cancer-related mutations
on DNA polymerases.
Collapse
Affiliation(s)
- Alice R Walker
- Department of Chemistry, University of North Texas , 1155 Union Circle, Denton, Texas 76203, United States
| | - G Andrés Cisneros
- Department of Chemistry, University of North Texas , 1155 Union Circle, Denton, Texas 76203, United States
| |
Collapse
|
6
|
Maximoff SN, Kamerlin SCL, Florián J. DNA Polymerase λ Active Site Favors a Mutagenic Mispair between the Enol Form of Deoxyguanosine Triphosphate Substrate and the Keto Form of Thymidine Template: A Free Energy Perturbation Study. J Phys Chem B 2017; 121:7813-7822. [PMID: 28732447 DOI: 10.1021/acs.jpcb.7b04874] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
Human DNA polymerase λ is an intermediate fidelity member of the X family, which plays a role in DNA repair. Recent X-ray diffraction structures of a ternary complex of a loop-deletion mutant of polymerase λ, a deoxyguanosine triphosphate analogue, and a gapped DNA show that guanine and thymine form a mutagenic mispair with an unexpected Watson-Crick-like geometry rather than a wobble geometry. Hence, there is an intriguing possibility that either thymine in the DNA or guanine in the deoxyguanosine triphosphate analogue may spend a substantial fraction of time in a deprotonated or enol form (both are minor species in aqueous solution) in the active site of the polymerase λ mutant. The experiments do not determine particular forms of the nucleobases that contribute to this mutagenic mispair. Thus, we investigate the thermodynamics of formation of various mispairs between guanine and thymine in the ternary complex at a neutral pH using classical molecular dynamics simulations and the free energy perturbation method. Our free energy calculations, as well as a comparison of the experimental and computed structures of mispairs, indicate that the Watson-Crick-like mispair between the enol tautomer of guanine and the keto tautomer of thymine is dominant. The wobble mispair between the keto forms of guanine and thymine and the Watson-Crick-like mispair between the keto tautomer of guanine and the enol tautomer of thymine are less prevalent, and mispairs that involve deprotonated guanine or thymine are thermodynamically unlikely. These findings are consistent with the experiment and relevant for understanding mechanisms of spontaneous mutagenesis.
Collapse
Affiliation(s)
- Sergey N Maximoff
- Department of Chemistry and Biochemistry, Loyola University Chicago , 1032 W. Sheridan Road, Chicago, Illinois 60660, United States
| | | | - Jan Florián
- Department of Chemistry and Biochemistry, Loyola University Chicago , 1032 W. Sheridan Road, Chicago, Illinois 60660, United States
| |
Collapse
|
7
|
Klvaňa M, Bren U, Florián J. Uniform Free-Energy Profiles of the P-O Bond Formation and Cleavage Reactions Catalyzed by DNA Polymerases β and λ. J Phys Chem B 2016; 120:13017-13030. [PMID: 27992186 PMCID: PMC5217713 DOI: 10.1021/acs.jpcb.6b08581] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
![]()
Human
X-family DNA polymerases β (Polβ) and λ
(Polλ) catalyze the nucleotidyl-transfer reaction in the base
excision repair pathway of the cellular DNA damage response. Using
empirical valence bond and free-energy perturbation simulations, we
explore the feasibility of various mechanisms for the deprotonation
of the 3′-OH group of the primer DNA strand, and the subsequent
formation and cleavage of P–O bonds in four Polβ, two
truncated Polλ (tPolλ), and two tPolλ Loop1 mutant
(tPolλΔL1) systems differing in the initial X-ray crystal
structure and nascent base pair. The average calculated activation
free energies of 14, 18, and 22 kcal mol–1 for Polβ,
tPolλ, and tPolλΔL1, respectively, reproduce the
trend in the observed catalytic rate constants. The most feasible
reaction pathway consists of two successive steps: specific base (SB)
proton transfer followed by rate-limiting concerted formation and
cleavage of the P–O bonds. We identify linear free-energy relationships
(LFERs) which show that the differences in the overall activation
and reaction free energies among the eight studied systems are determined
by the reaction free energy of the SB proton transfer. We discuss
the implications of the LFERs and suggest pKa of the 3′-OH group as a predictor of the catalytic
rate of X-family DNA polymerases.
Collapse
Affiliation(s)
- Martin Klvaňa
- Laboratory of Physical Chemistry and Chemical Thermodynamics, Faculty of Chemistry and Chemical Technology, University of Maribor , Smetanova ulica 17, 2000 Maribor, Slovenia.,Department of Chemistry and Biochemistry, Loyola University Chicago , 1032 W. Sheridan Road, Chicago, Illinois 60660, United States
| | - Urban Bren
- Laboratory of Physical Chemistry and Chemical Thermodynamics, Faculty of Chemistry and Chemical Technology, University of Maribor , Smetanova ulica 17, 2000 Maribor, Slovenia.,Laboratory for Molecular Modeling, National Institute of Chemistry , Hajdrihova ulica 19, 1001 Ljubljana, Slovenia
| | - Jan Florián
- Department of Chemistry and Biochemistry, Loyola University Chicago , 1032 W. Sheridan Road, Chicago, Illinois 60660, United States
| |
Collapse
|
8
|
|
9
|
Lameira J, Ram Prasad B, Chu ZT, Warshel A. Methyltransferases do not work by compression, cratic, or desolvation effects, but by electrostatic preorganization. Proteins 2015; 83:318-30. [PMID: 25388538 PMCID: PMC4300294 DOI: 10.1002/prot.24717] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2014] [Revised: 10/30/2014] [Accepted: 11/03/2014] [Indexed: 11/05/2022]
Abstract
The enzyme catechol O-methyltransferase (COMT) catalyzes the transfer of a methyl group from S-adenosylmethionine to dopamine and related catechols. The search for the origin of COMT catalysis has led to different proposals and hypothesis, including the entropic, the NAC, and the compression proposals as well as the more reasonable electrostatic idea. Thus, it is important to understand the catalytic power of this enzyme and to examine the validity of different proposals and in particular the repeated recent implication of the compression idea. The corresponding analysis should be done by well-defined physically-based considerations that involve computations rather than circular interpretations of experimental results. Thus, we explore here the origin of the catalytic efficiency of COMT by using the empirical valence bond and the linear response approximation approaches. The results demonstrate that the catalytic effect of COMT is mainly due to electrostatic preorganization effects. It is also shown that the compression, NAC and entropic proposals do not account for the catalytic effect.
Collapse
Affiliation(s)
- Jeronimo Lameira
- University of Southern California, Department of Chemistry, SGM 418, 3620 McClintosk Avenue, Los Angeles, California 90089, United States
- Faculdade de Biotecnologia e Laboratório de Planejamento e Desenvolvimento de Fármacos; Universidade Federal do Pará, 66075-110, Belém, PA, Brazil
| | - B Ram Prasad
- University of Southern California, Department of Chemistry, SGM 418, 3620 McClintosk Avenue, Los Angeles, California 90089, United States
| | - Zhen T. Chu
- University of Southern California, Department of Chemistry, SGM 418, 3620 McClintosk Avenue, Los Angeles, California 90089, United States
| | - Arieh Warshel
- University of Southern California, Department of Chemistry, SGM 418, 3620 McClintosk Avenue, Los Angeles, California 90089, United States
| |
Collapse
|
10
|
Sharir-Ivry A, Varatharaj R, Shurki A. Challenges within the Linear Response Approximation When Studying Enzyme Catalysis and Effects of Mutations. J Chem Theory Comput 2014; 11:293-302. [DOI: 10.1021/ct500751f] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Avital Sharir-Ivry
- Institute for Drug Design,
School of Pharmacy, The Lise Meitner-Minerva Center for Computational
Quantum Chemistry, The Hebrew University of Jerusalem, Jerusalem 91120, Israel
| | - Rajapandian Varatharaj
- Institute for Drug Design,
School of Pharmacy, The Lise Meitner-Minerva Center for Computational
Quantum Chemistry, The Hebrew University of Jerusalem, Jerusalem 91120, Israel
| | - Avital Shurki
- Institute for Drug Design,
School of Pharmacy, The Lise Meitner-Minerva Center for Computational
Quantum Chemistry, The Hebrew University of Jerusalem, Jerusalem 91120, Israel
| |
Collapse
|
11
|
Bren U, Fuchs JE, Oostenbrink C. Cooperative binding of aflatoxin B1 by cytochrome P450 3A4: a computational study. Chem Res Toxicol 2014; 27:2136-47. [PMID: 25398138 DOI: 10.1021/tx5004062] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Aflatoxin B1 (AFB1)-the most potent natural carcinogen known to men-is metabolized by cytochrome P450 3A4 (CYP3A4), either to the genotoxic AFB1 exo-8,9-epoxide or to the detoxified 3α-hydroxy AFB1. The activation of the procarcinogen proceeds in a highly cooperative fashion, which differs from common allosteric regulation in the sense that it can be attributed to simultaneous occupancy of a single large and malleable active site by multiple ligand molecules. Unfortunately, unlike in the case of ketoconazole, there is currently no experimental structure available for the doubly ligated CYP3A4-AFB1 complex. Therefore, we employed a sequential molecular docking protocol to create various possible doubly ligated complexes and subsequently performed molecular dynamics simulations and free-energy calculations to check for their consistency with the available experimental data on regio- and stereoselectivity of both AFB1 oxidations as well as with available kinetic data. Only the system in which the first AFB1 molecule was bound in a face-on C8-C9 epoxidation mode and the second AFB1 molecule was bound in a side-on 3α-hydroxylation mode-a result of an unconstrained molecular docking protocol-has successfully fulfilled all the imposed criteria and is therefore proposed as the most likely structure of the doubly ligated complex of CYP3A4 with AFB1. The empirical Linear Interaction Energy method revealed that shape complementarity through nonpolar dispersion interactions between the two bound AFB1 molecules is the main source of the experimentally observed positive homotropic cooperativity. The reported study represents a nice example of how state-of-the-art molecular modeling techniques can be used to study complicated macromolecular complexes, whose structures have not yet been experimentally determined, and to validate these against the available experimental data. The proposed structure will facilitate future studies on the rational design of successful AFB1 modulators or on human subpopulations characterized by specific CYP3A4 polymorphisms that are especially sensitive to AFB1.
Collapse
Affiliation(s)
- Urban Bren
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Muthgasse 18, AT-1190 Vienna, Austria
| | | | | |
Collapse
|
12
|
Zuo ZL, Guo L, Mancera RL. Free energy of binding of coiled-coil complexes with different electrostatic environments: the influence of force field polarisation and capping. NATURAL PRODUCTS AND BIOPROSPECTING 2014; 4:285-295. [PMID: 25159896 PMCID: PMC4199946 DOI: 10.1007/s13659-014-0036-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/29/2014] [Accepted: 08/08/2014] [Indexed: 06/03/2023]
Abstract
Coiled-coils are well known protein-protein interaction motifs, with the leucine zipper region of activator protein-1 (AP-1) consisting of the c-Jun and c-Fos proteins being a typical example. Molecular dynamics (MD) simulations using the MM/GBSA method have been used to predict the free energy of interaction of these proteins. The influence of force field polarisation and capping on the predicted free energy of binding of complexes with different electrostatic environments (net charge) were investigated. Although both force field polarisation and peptide capping are important for the prediction of the absolute free energy of binding, peptide capping has the largest influence on the predicted free energy of binding. Polarisable simulations appear better suited to determine structural properties of the complexes of these proteins while non-polarisable simulations seem to give better predictions of the associated free energies of binding.
Collapse
Affiliation(s)
- Zhi-Li Zuo
- State Key Laboratory of Phytochemistry and Plant Resources in West China, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, 650201 China
| | - Ling Guo
- College of Animal Husbandry & Veterinary, Liaoning Medical University, Jinzhou, 121001 China
| | - Ricardo L. Mancera
- School of Biomedical Sciences, CHIRI Biosciences, Curtin University, GPO Box U1987, Perth, WA 6845 Australia
| |
Collapse
|
13
|
Elias AA, Cisneros GA. Computational study of putative residues involved in DNA synthesis fidelity checking in Thermus aquaticus DNA polymerase I. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2014; 96:39-75. [PMID: 25443954 DOI: 10.1016/bs.apcsb.2014.06.003] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
A fidelity-checking site for DNA polymerase I has been proposed based on recent single-molecule Förster resonance energy transfer studies. The checking site is believed to ensure proper base pairing of the newly inserted nucleotide. Computational studies have been utilized to predict residues involved in this putative checking site on the Klenow and Bacillus fragments. Here, we employ energy decomposition analysis, electrostatic free energy response, and noncovalent interaction plots to identify the residues involved in the hypothesized checking site in the homologous Klenow fragment from Thermus aquaticus (Klentaq). Our results indicate multiple protein residues that show altered interactions for three mispairs compared to the correctly paired DNA dimer. Many of these residues are also conserved along A family polymerases.
Collapse
Affiliation(s)
- Angela A Elias
- Department of Chemistry, Wayne State University, Detroit, Michigan, USA
| | - G Andrés Cisneros
- Department of Chemistry, Wayne State University, Detroit, Michigan, USA.
| |
Collapse
|
14
|
Sund J, Lind C, Åqvist J. Binding site preorganization and ligand discrimination in the purine riboswitch. J Phys Chem B 2014; 119:773-82. [PMID: 25014157 DOI: 10.1021/jp5052358] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The progress of RNA research has suggested a wide variety of RNA molecules as possible targets for pharmaceutical drug molecules. Structure-based computational methods for predicting binding modes and affinities are now important tools in drug discovery, but these methods have mainly been focused on protein targets. Here we employ molecular dynamics free-energy perturbation calculations and the linear interaction energy method to analyze the energetics of ligand binding to purine riboswitches. Calculations are carried out for 14 different purine complexes with the guanine and adenine riboswitches in order to examine their ligand recognition principles. The simulations yield binding affinities in good agreement with experimental data and rationalize the selectivity of the riboswitches for different ligands. In particular, it is found that these receptors have an unusually high degree of electrostatic preorganization for their cognate ligands, and this effect is further quantified by explicit free-energy calculations, which show that the standard electrostatic linear interaction energy parametrization is suboptimal in this case. The adenine riboswitch specifically uses the electrostatic preorganization to discriminate against guanine by preventing the formation of a G-U wobble base pair.
Collapse
Affiliation(s)
- Johan Sund
- Department of Cell and Molecular Biology, Uppsala University , Biomedical Center, Box 596, SE-751 24 Uppsala, Sweden
| | | | | |
Collapse
|
15
|
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
|
16
|
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
|
17
|
Mishra SK, Sund J, Åqvist J, Koča J. Computational prediction of monosaccharide binding free energies to lectins with linear interaction energy models. J Comput Chem 2012; 33:2340-50. [DOI: 10.1002/jcc.23081] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2012] [Revised: 06/25/2012] [Accepted: 07/09/2012] [Indexed: 01/28/2023]
|
18
|
Graham SE, Syeda F, Cisneros GA. Computational prediction of residues involved in fidelity checking for DNA synthesis in DNA polymerase I. Biochemistry 2012; 51:2569-78. [PMID: 22397306 DOI: 10.1021/bi201856m] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Recent single-molecule Förster resonance energy transfer studies of DNA polymerase I have led to the proposal of a postinsertion fidelity-checking site. This site is hypothesized to ensure proper base pairing of the newly inserted nucleotide. To help test this hypothesis, we have used energy decomposition, electrostatic free energy response, and noncovalent interaction analysis analyses to identify residues involved in this putative checking site. We have used structures of DNA polymerase I from two different organisms, the Klenow fragment from Escherichia coli and the Bacillus fragment from Bacillus stearothermophilus. Our results point to several residues that show altered interactions for three mispairs compared to the correctly paired DNA dimer. Furthermore, many of these residues are conserved among A family polymerases. The identified residues provide potential targets for mutagenesis studies for investigation of the fidelity-checking site hypothesis.
Collapse
Affiliation(s)
- Sarah E Graham
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | | | | |
Collapse
|
19
|
Frushicheva MP, Warshel A. Towards quantitative computer-aided studies of enzymatic enantioselectivity: the case of Candida antarctica lipase A. Chembiochem 2011; 13:215-23. [PMID: 22190449 DOI: 10.1002/cbic.201100600] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2011] [Indexed: 11/05/2022]
Abstract
The prospect for consistent computer-aided refinement of stereoselective enzymes is explored by simulating the hydrolysis of enantiomers of an α-substituted ester by wild-type and mutant Candida antarctica lipase A, using several strategies. In particular, we focused on the use of the empirical valence bond (EVB) method in a quantitative screening for enantioselectivity, and evaluate both k(cat) and k(cat)/K(M) of the R and S stereoisomers. We found that an extensive sampling is essential for obtaining converging results. This requirement points towards possible problems with approaches that use a limited conformational sampling. However, performing the proper sampling appears to give encouraging results and to offer a powerful tool for the computer-aided design of enantioselective enzymes. We also explore faster strategies for identifying mutations that will help in augmenting directed-evolution experiments, but these approaches require further refinement.
Collapse
Affiliation(s)
- Maria P Frushicheva
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089-1062, USA
| | | |
Collapse
|
20
|
Warshel A, Dryga A. Simulating electrostatic energies in proteins: perspectives and some recent studies of pKas, redox, and other crucial functional properties. Proteins 2011; 79:3469-84. [PMID: 21910139 DOI: 10.1002/prot.23125] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2011] [Revised: 05/09/2011] [Accepted: 06/09/2011] [Indexed: 01/30/2023]
Abstract
Electrostatic energies provide what is arguably the most effective tool for structure-function correlation of biological molecules. Here, we provide an overview of the current state-of-the-art simulations of electrostatic energies in macromolecules, emphasizing the microscopic perspective but also relating it to macroscopic approaches. We comment on the convergence issue and other problems of the microscopic models and the ways of keeping the microscopic physics while moving to semi-macroscopic directions. We discuss the nature of the protein dielectric "constants" reiterating our long-standing point that the dielectric "constants" in semi-macroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering different functional properties including pK(a)'s, redox potentials, ion and proton channels, enzyme catalysis, ligand binding, and protein stability. We emphasize the microscopic overcharging approach for studying pK(a) 's of internal groups in proteins and give a demonstration of power of this approach. We also emphasize recent advances in coarse grained models with a physically based electrostatic treatment and provide some examples including further directions in treating voltage activated ion channels.
Collapse
Affiliation(s)
- Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA.
| | | |
Collapse
|
21
|
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
|
22
|
Klvaňa M, Jeřábek P, Goodman MF, Florián J. An abridged transition state model to derive structure, dynamics, and energy components of DNA polymerase β fidelity. Biochemistry 2011; 50:7023-32. [PMID: 21739967 DOI: 10.1021/bi200790s] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We show how a restricted reaction surface can be used to facilitate the calculation of biologically important contributions of active site geometries and dynamics to DNA polymerase fidelity. Our analysis, using human DNA polymerase beta (pol β), is performed within the framework of an electrostatic linear free energy response (EFER) model. The structure, dynamics, and energetics of pol β-DNA-dNTP interactions are computed between two points on the multidimensional reaction free energy surface. "Point 1" represents a ground state activation intermediate (GSA), which is obtained by deprotonating the terminal 3'OH group of the primer DNA strand. "Point 2" is the transition state (PTS) for the attack of the 3'O(-) (O(nuc)) on the P(α) atom of dNTP substrate, having the electron density of a dianionic phosphorane intermediate. Classical molecular dynamics simulations are used to compute the geometric and dynamic contributions to the formation of right and wrong O(nuc)-P chemical bonds. Matched dCTP·G and mismatched dATP·G base pairs are used to illustrate the analysis. Compared to the dCTP·G base pair, the dATP·G mismatch has fewer GSA configurations with short distances between O(nuc) and P(α) atoms and between the oxygen in the scissile P-O bond (O(lg)) and the nearest structural water. The thumb subdomain conformation of the GSA complex is more open for the mismatch, and the H-bonds in the mispair become more extended during the nucleophilic attack than in the correct pair. The electrostatic contributions of pol β and DNA residues to catalysis of the right and wrong P-O(nuc) bond formation are 5.3 and 3.1 kcal/mol, respectively, resulting in an 80-fold contribution to fidelity. The EFER calculations illustrate the considerable importance of Arg183 and an O(lg)-proximal water molecule to pol β fidelity.
Collapse
Affiliation(s)
- Martin Klvaňa
- Department of Chemistry, Loyola University Chicago, Chicago, Illinois 60626, United States
| | | | | | | |
Collapse
|
23
|
Cerón-Carrasco JP, Zúñiga J, Requena A, Perpète EA, Michaux C, Jacquemin D. Combined effect of stacking and solvation on the spontaneous mutation in DNA. Phys Chem Chem Phys 2011; 13:14584-9. [DOI: 10.1039/c1cp20946a] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
24
|
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
|
25
|
Singh N, Warshel A. Absolute binding free energy calculations: on the accuracy of computational scoring of protein-ligand interactions. Proteins 2010; 78:1705-23. [PMID: 20186976 DOI: 10.1002/prot.22687] [Citation(s) in RCA: 100] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Calculating the absolute binding free energies is a challenging task. Reliable estimates of binding free energies should provide a guide for rational drug design. It should also provide us with deeper understanding of the correlation between protein structure and its function. Further applications may include identifying novel molecular scaffolds and optimizing lead compounds in computer-aided drug design. Available options to evaluate the absolute binding free energies range from the rigorous but expensive free energy perturbation to the microscopic linear response approximation (LRA/beta version) and related approaches including the linear interaction energy (LIE) to the more approximated and considerably faster scaled protein dipoles Langevin dipoles (PDLD/S-LRA version) as well as the less rigorous molecular mechanics Poisson-Boltzmann/surface area (MM/PBSA) and generalized born/surface area (MM/GBSA) to the less accurate scoring functions. There is a need for an assessment of the performance of different approaches in terms of computer time and reliability. We present a comparative study of the LRA/beta, the LIE, the PDLD/S-LRA/beta, and the more widely used MM/PBSA and assess their abilities to estimate the absolute binding energies. The LRA and LIE methods perform reasonably well but require specialized parameterization for the nonelectrostatic term. The PDLD/S-LRA/beta performs effectively without the need of reparameterization. Our assessment of the MM/PBSA is less optimistic. This approach appears to provide erroneous estimates of the absolute binding energies because of its incorrect entropies and the problematic treatment of electrostatic energies. Overall, the PDLD/S-LRA/beta appears to offer an appealing option for the final stages of massive screening approaches.
Collapse
Affiliation(s)
- Nidhi Singh
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-1062, USA
| | | |
Collapse
|
26
|
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]
|
27
|
Rucker R, Oelschlaeger P, Warshel A. A binding free energy decomposition approach for accurate calculations of the fidelity of DNA polymerases. Proteins 2010; 78:671-80. [PMID: 19842163 DOI: 10.1002/prot.22596] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
DNA polymerase beta (pol beta) is a small eukaryotic enzyme with the ability to repair short single-stranded DNA gaps that has found use as a model system for larger replicative DNA polymerases. For all DNA polymerases, the factors determining their catalytic power and fidelity are the interactions between the bases of the base pair, amino acids near the active site, and the two magnesium ions. In this report, we study effects of all three aspects on human pol beta transition state (TS) binding free energies by reproducing a consistent set of experimentally determined data for different structures. Our calculations comprise the combination of four different base pairs (incoming pyrimidine nucleotides incorporated opposite both matched and mismatched purines) with four different pol beta structures (wild type and three mutants). We generate three fragments of the incoming deoxynucleoside 5'-triphosphate-TS and run separate calculations for the neutral base part and the highly charged triphosphate part, using different dielectric constants in order to account for the specific dielectric response. This new approach improves our ability to predict the effect of matched and mismatched base pairing and of mutations in DNA polymerases on fidelity and may be a useful tool in studying the potential of DNA polymerase mutations in the development of cancer. It also supports our point of view with regards to the origin of the structural control of fidelity, allowing for a quantified description of the fidelity of DNA polymerases.
Collapse
Affiliation(s)
- Robert Rucker
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | | |
Collapse
|
28
|
Bren U, Lah J, Bren M, Martínek V, Florián J. DNA duplex stability: the role of preorganized electrostatics. J Phys Chem B 2010; 114:2876-85. [PMID: 20131770 PMCID: PMC2841231 DOI: 10.1021/jp9064246] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The insertion of a DNA base moiety at the end of a DNA duplex to form a Watson-Crick or wobble pair during DNA annealing or replication is a step of fundamental biological importance. Therefore, we investigated the energetics of a formation of the terminal G x C, G x T, and G x A base pairs in DNA containing a 5'-dangling G adjacent to the base insertion point using differential scanning calorimetry and computer simulations. The energies calculated along classical molecular dynamics trajectories in aqueous solution were analyzed in the framework of linear-response approximation (LRA) to obtain relative free energies for the base insertion and their electrostatic, van der Waals, and preorganization components. Using the generic set of LRA parameters, the calculated free energies disfavored the mispair formation by 2.5 (G x C --> G x T) and 1.7 (G x C --> G x A) kcal/mol, in reasonable agreement with the experimental free energy differences of 1.8 and 1.4 kcal/mol, respectively. The calculated preorganization components of these free energies of 0.6 (G x C --> G x T) and -0.1 (G x C --> G x A) kcal/mol show that electrostatic preorganization, which is an important source of DNA replication fidelity, plays a lesser role in the mispair destabilization in the absence of DNA polymerase.
Collapse
Affiliation(s)
- Urban Bren
- National Institute of Chemistry, Hajdrihova 19, 1001 Ljubljana, Slovenia
| | - Jurij Lah
- University of Ljubljana, Faculty of Chemistry and Chemical Technology, Aškerčeva 5, 1000 Ljubljana, Slovenia
| | - Matevž Bren
- Institute of Mathematics, Physics and Mechanics, Jadranska 19, 1000 Ljubljana, Slovenia
| | - Václav Martínek
- Department of Chemistry, Loyola University Chicago, Chicago, IL 60626, USA
| | - Jan Florián
- Department of Chemistry, Loyola University Chicago, Chicago, IL 60626, USA
| |
Collapse
|
29
|
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
|
30
|
Cerón-Carrasco JP, Requena A, Zúñiga J, Michaux C, Perpète EA, Jacquemin D. Intermolecular Proton Transfer in Microhydrated Guanine−Cytosine Base Pairs: a New Mechanism for Spontaneous Mutation in DNA. J Phys Chem A 2009; 113:10549-56. [DOI: 10.1021/jp906551f] [Citation(s) in RCA: 69] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
| | | | | | - C. Michaux
- Laboratorie de Chimie Biologique Structurale, Facultés Universitaires Notre-Dame de la Paix, Rue de Bruxelles, 61, 5000 Namur, Belgium
| | | | | |
Collapse
|
31
|
Exner TE. Insights into the high fidelity of a DNA polymerase I mutant. J Mol Model 2009; 15:1271-80. [PMID: 19333629 DOI: 10.1007/s00894-009-0491-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2008] [Accepted: 03/12/2009] [Indexed: 10/21/2022]
Abstract
Mutants of DNA polymerase I from Thermus aquaticus (Taq) with higher fidelity compared to the wild type enzyme were identified in an earlier study by Summerer et al. (Angew Chem Int Ed 44:4712-4715, 2005). Here, one of these mutants, PLQ (consensus residues 879-881), was analysed using molecular dynamics simulations. This was done by calculating the structures of the ternary complex comprising the enzyme, the DNA primer and template as well as the incoming nucleotide before the chemical reaction for the Watson-Crick and different mismatched base pairings. The results show that the high fidelity of the mutant can be explained partly by different specific interactions between the amino acids of the enzyme and the DNA primer end as well as, in some mismatches, a displacement of the primer relative to the incoming deoxyribonucleoside triphosphate and the catalytic magnesium ion. This displacement is facilitated by reduced steric interactions between the enzyme and the DNA.
Collapse
Affiliation(s)
- Thomas E Exner
- Department of Chemistry and Zukunftskolleg, University of Konstanz, Germany.
| |
Collapse
|
32
|
Perdih A, Bren U, Solmajer T. Binding free energy calculations of N-sulphonyl-glutamic acid inhibitors of MurD ligase. J Mol Model 2009; 15:983-96. [PMID: 19198900 DOI: 10.1007/s00894-009-0455-8] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2008] [Accepted: 01/06/2009] [Indexed: 10/21/2022]
Abstract
The increasing incidence of bacterial resistance to most available antibiotics has underlined the urgent need for the discovery of novel efficacious antibacterial agents. The biosynthesis of bacterial peptidoglycan, where the MurD enzyme is involved in the intracellular phase of UDP-MurNAc-pentapeptide formation, represents a collection of highly selective targets for novel antibacterial drug design. Structural studies of N-sulfonyl-glutamic acid inhibitors of MurD have made possible the examination of binding modes of this class of compounds, providing valuable information for the lead optimization phase of the drug discovery cycle. Binding free energies were calculated for a series of MurD N-sulphonyl-Glu inhibitors using the linear interaction energy (LIE) method. Analysis of interaction energy during the 20-ns MD trajectories revealed non-polar van der Waals interactions as the main driving force for the binding of these inhibitors, and excellent agreement with the experimental free energies was obtained. Calculations of binding free energies for selected moieties of compounds in this structural class substantiated even deeper insight into the source of inhibitory activity. These results constitute new valuable information to further assist the lead optimization process.
Collapse
Affiliation(s)
- Andrej Perdih
- Laboratory for Molecular Modeling and NMR Spectroscopy, National Institute of Chemistry, Hajdrihova 19, 1001, Ljubljana, Slovenia
| | | | | |
Collapse
|
33
|
Guo Q, Jureller JE, Warren JT, Solomaha E, Florián J, Tang WJ. Protein-protein docking and analysis reveal that two homologous bacterial adenylyl cyclase toxins interact with calmodulin differently. J Biol Chem 2008; 283:23836-45. [PMID: 18583346 DOI: 10.1074/jbc.m802168200] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
Calmodulin (CaM), a eukaryotic calcium sensor that regulates diverse biological activities, consists of N- and C-terminal globular domains (N-CaM and C-CaM, respectively). CaM serves as the activator of CyaA, a 188-kDa adenylyl cyclase toxin secreted by Bordetella pertussis, which is the etiologic agent for whooping cough. Upon insertion of the N-terminal adenylyl cyclase domain (ACD) of CyaA to its targeted eukaryotic cells, CaM binds to this domain tightly ( approximately 200 pm affinity). This interaction activates the adenylyl cyclase activity of CyaA, leading to a rise in intracellular cAMP levels to disrupt normal cellular signaling. We recently solved the structure of CyaA-ACD in complex with C-CaM to elucidate the mechanism of catalytic activation. However, the structure of the interface between N-CaM and CyaA, the formation of which contributes a 400-fold increase of binding affinity between CyaA and CaM, remains elusive. Here, we used site-directed mutations and molecular dynamic simulations to generate several working models of CaM-bound CyaA-ACD. The validity of these models was evaluated by disulfide bond cross-linking, point mutations, and fluorescence resonance energy transfer experiments. Our study reveals that a beta-hairpin region (amino acids 259-273) of CyaA-ACD likely makes contacts with the second calcium binding motif of the extended CaM. This mode of interaction differs from the interaction of N-CaM with anthrax edema factor, which binds N-CaM via its helical domain. Thus, two structurally conserved, bacterial adenylyl cyclase toxins have evolved to utilize distinct binding surfaces and modes of activation in their interaction with CaM, a highly conserved eukaryotic signaling protein.
Collapse
Affiliation(s)
- Qing Guo
- Ben May Department for Cancer Biology, The University of Chicago, 929 East 57th Street, Chicago, IL 60637, USA
| | | | | | | | | | | |
Collapse
|
34
|
Calculation of protein-ligand binding free energy by using a polarizable potential. Proc Natl Acad Sci U S A 2008; 105:6290-5. [PMID: 18427113 DOI: 10.1073/pnas.0711686105] [Citation(s) in RCA: 183] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The binding of charged ligands benzamidine and diazamidine to trypsin was investigated by using a polarizable potential energy function and explicit-water molecular dynamics simulations. The binding free energies were computed from the difference between the free energies of decoupling the ligand from water and protein environments. Both the absolute and the relative free energies from the perturbation simulations agree with experimental measurements to within 0.5 kcal.mol(-1). Comparison of free-energy components sampled from different thermodynamic paths indicates that electrostatics is the main driving force behind benzamidine recognition of trypsin. The contribution of electronic polarization to binding appears to be crucial. By computing the free-energy contribution caused by the polarization between the ligand and its surroundings, we found that polarization has the opposite effect in dissimilar environments. Although polarization favors ligand solvation in water, it weakens the protein-ligand attraction by screening the electrostatic interaction between trypsin and benzamidine. We also examined the relative binding free energies of a benzamidine analog diazamidine to trypsin. The changes in free energy on benzamidine-diazamidine substitution were tens of kilocalories in both water and trypsin environments; however, the change in the total binding free energy is <2 kcal.mol(-1) because of cancellation, consistent with the experimental results. Overall, our results suggest that the use of a polarizable force field, given adequate sampling, is capable of achieving chemical accuracy in molecular simulations of protein-ligand recognition.
Collapse
|
35
|
Xiang Y, Goodman MF, Beard WA, Wilson SH, Warshel A. Exploring the role of large conformational changes in the fidelity of DNA polymerase beta. Proteins 2008; 70:231-47. [PMID: 17671961 PMCID: PMC2365506 DOI: 10.1002/prot.21668] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
The relationships between the conformational landscape, nucleotide insertion catalysis and fidelity of DNA polymerase beta are explored by means of computational simulations. The simulations indicate that the transition states for incorporation of right (R) and wrong (W) nucleotides reside in substantially different protein conformations. The protein conformational changes that reproduce the experimentally observed fidelity are significantly larger than the small rearrangements that usually accompany motions from the reactant state to the transition state in common enzymatic reactions. Once substrate binding has occurred, different constraints imposed on the transition states for insertion of R and W nucleotides render it highly unlikely that both transition states can occur in the same closed structure, because the predicted fidelity would then be many orders of magnitude too large. Since the conformational changes reduce the transition state energy of W incorporation drastically they decrease fidelity rather than increase it. Overall, a better agreement with experimental data is attained when the R is incorporated through a transition state in a closed conformation and W is incorporated through a transition state in one or perhaps several partially open conformations. The generation of free energy surfaces for R and W also allow us to analyze proposals about the relationship between induced fit and fidelity.
Collapse
Affiliation(s)
- Yun Xiang
- Department of Chemistry, University of Southern California, Los Angeles, California 90089
- *Correspondence to: Yun Xiang, Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, CA 90089. E-mail: or Arieh Warshel, Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, CA 90089. E-mail:
| | - Myron F. Goodman
- Department of Chemistry, University of Southern California, Los Angeles, California 90089
- Department of Biological Sciences, University of Southern California, Los Angeles, California 90089
| | - William A. Beard
- Laboratory of Structural Biology, NIEHS, National Institutes of Health, DHHS, Research Triangle Park, North California 27709
| | - Samuel H. Wilson
- Laboratory of Structural Biology, NIEHS, National Institutes of Health, DHHS, Research Triangle Park, North California 27709
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, California 90089
- *Correspondence to: Yun Xiang, Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, CA 90089. E-mail: or Arieh Warshel, Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, CA 90089. E-mail:
| |
Collapse
|
36
|
Warshel A, Kato M, Pisliakov AV. Polarizable Force Fields: History, Test Cases, and Prospects. J Chem Theory Comput 2007; 3:2034-45. [DOI: 10.1021/ct700127w] [Citation(s) in RCA: 278] [Impact Index Per Article: 16.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Arieh Warshel
- University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, California 90089-1062
| | - Mitsunori Kato
- University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, California 90089-1062
| | - Andrei V. Pisliakov
- University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, California 90089-1062
| |
Collapse
|
37
|
Warshel A, Sharma PK, Kato M, Xiang Y, Liu H, Olsson MHM. Electrostatic basis for enzyme catalysis. Chem Rev 2007; 106:3210-35. [PMID: 16895325 DOI: 10.1021/cr0503106] [Citation(s) in RCA: 991] [Impact Index Per Article: 58.3] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Arieh Warshel
- Department of Chemistry, University of Southern California, SGM Building 418, 3620 McClintock Avenue, Los Angeles, California 90089-1062, USA.
| | | | | | | | | | | |
Collapse
|
38
|
Bren U, Martínek V, Florian J. Free energy simulations of uncatalyzed DNA replication fidelity: structure and stability of T.G and dTTP.G terminal DNA mismatches flanked by a single dangling nucleotide. J Phys Chem B 2007; 110:10557-66. [PMID: 16722767 DOI: 10.1021/jp060292b] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A reference system for DNA replication fidelity was studied by free energy perturbation (FEP) and linear interaction energy (LIE) methods. The studied system included a hydrated duplex DNA with the 5'-CG dangling end of the templating strand, and dCTP4-.Mg2+ or dTTP4-.Mg2+ inserted opposite the dangling G to form a correct (i.e., Watson-Crick) or incorrect (i.e., wobble) base pair, respectively. The average distance between the 3'-terminal oxygen of the primer strand and the alpha-phosphorus of dNTP was found to be 0.2 A shorter for the correct base pair than for the incorrect base pair. Binding of the incorrect dNTP was found to be disfavored by 0.4 kcal/mol relative to the correct dNTP. We estimated that improved binding and more near-attack configurations sampled by the correct base pair should translate in aqueous solution and in the absence of DNA polymerase into a six times faster rate for the incorporation of the correct dNTP into DNA. The accuracy of the calculated binding free energy difference was verified by examining the relative free energy for melting duplex DNA containing GC and GT terminal base pairs flanked by a 5' dangling C. The calculated LIE and FEP free energies of 1.7 and 1.1 kcal/mol, respectively, compared favorably with the experimental estimate of 1.4 kcal/mol obtained using the nearest neighbor parameters. To decompose the calculated free energies into additive electrostatic and van der Waals contributions and to provide a set of rigorous theoretical data for the parametrization of the LIE method, we suggested a variant of the FEP approach, for which we coined a binding-relevant free energy (BRFE) acronym. BRFE approach is characterized by its unique perturbation pathway and by its exclusion of the intramolecular energy of a rigid part of the ligand from the total potential energy.
Collapse
Affiliation(s)
- Urban Bren
- Department of Chemistry, Loyola University Chicago, Chicago, Illinois 60660, USA
| | | | | |
Collapse
|
39
|
Bren U, Martínek V, Florián J. Decomposition of the solvation free energies of deoxyribonucleoside triphosphates using the free energy perturbation method. J Phys Chem B 2007; 110:12782-8. [PMID: 16800613 DOI: 10.1021/jp056623m] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Free energy perturbation (FEP) calculations using the Amber 95 force field and the TIP3P water model were carried out to evaluate the solvation free energy of deoxyribonucleoside triphosphates in aqueous solution. Solvation free energies of -307.5, -311.5, -314.1, and -317.0 kcal/mol were calculated for the (Mg x dTTP)2-, (Mg x dATP)2-, (Mg x dCTP)2-, and (Mg x dGTP)2- complexes, respectively. Structural origins of the relative solvation free energies of deoxyribonucleoside phosphates were examined by calculating the contribution of the interaction of the base moiety with its surroundings. We showed that for each nucleobase the magnitude of this contribution is unaffected by substituting the 5'-OH group of the corresponding nucleoside with the charged mono- or triphosphate groups. This free energy contribution was further decomposed into the sum of free energies originating from the interactions of the base with itself, its substituent, water, and Na+ ions. Although the sum of these components was nearly constant over a wide range of solutes the individual free energy constituents varied significantly. Furthermore, this decomposition showed a high degree of additivity. Computational conditions necessary for obtaining additive free energy decomposition for the systems studied here within the framework of the FEP method included the use of a single mutation pathway and a subdivision of the FEP protocol into 51 or more windows.
Collapse
Affiliation(s)
- Urban Bren
- Department of Chemistry, Loyola University Chicago, Chicago, Illinois 60626, USA
| | | | | |
Collapse
|
40
|
Ngo JCK, Gullingsrud J, Giang K, Yeh MJ, Fu XD, Adams JA, McCammon JA, Ghosh G. SR protein kinase 1 is resilient to inactivation. Structure 2007; 15:123-33. [PMID: 17223538 DOI: 10.1016/j.str.2006.11.011] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2006] [Revised: 11/27/2006] [Accepted: 11/29/2006] [Indexed: 11/25/2022]
Abstract
SR protein kinase 1 (SRPK1) is a constitutively active kinase, which processively phosphorylates multiple serines within its substrates, ASF/SF2. We describe crystallographic, molecular dynamics, and biochemical results that shed light on how SRPK1 preserves its constitutive active conformation. Our structure reveals that unlike other known active kinase structures, the activation loop remains in an active state without any specific intraprotein interactions. Moreover, SRPK1 remains active despite extensive mutation to the activation segment. Molecular dynamics simulations reveal that SRPK1 partially absorbs the effect of mutations by forming compensatory interactions that maintain a catalytically competent chemical environment. Furthermore, SRPK1 is similarly resistant to deletion of its spacer loop region. Based upon a model of SRPK1 bound to a segment encompassing the docking motif and active-site peptide of ASF/SF2, we suggest a mechanism for processive phosphorylation and propose that the atypical resiliency we observed is critical for SRPK1's processive activity.
Collapse
Affiliation(s)
- Jacky Chi Ki Ngo
- Department of Chemistry & Biochemistry, University of California, San Diego, La Jolla, CA 92093, USA
| | | | | | | | | | | | | | | |
Collapse
|
41
|
|
42
|
Martínek V, Bren U, Goodman MF, Warshel A, Florián J. DNA polymerase beta catalytic efficiency mirrors the Asn279-dCTP H-bonding strength. FEBS Lett 2007; 581:775-80. [PMID: 17286973 PMCID: PMC2001272 DOI: 10.1016/j.febslet.2007.01.042] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2006] [Revised: 01/17/2007] [Accepted: 01/18/2007] [Indexed: 10/23/2022]
Abstract
Ternary complexes of wild type or mutant form of human DNA polymerase beta (pol beta) bound to DNA and dCTP substrates were studied by molecular dynamics (MD) simulations. The occurrences of contact configurations (CC) of structurally important atom pairs were sampled along the MD trajectories, and converted into free-energy differences, DeltaG(CC). DeltaG(CC) values were correlated with the experimental binding and catalytic free energies for the wild type pol beta and its Arg183Ala, Tyr271Ala, Asp276Val, Lys280Gly, Arg283Ala, and Glu295Ala mutants. The correlation coefficients show that the strength of the H-bond between dCTP and Asn279 is a strong predictor of the mutation-induced changes in the catalytic efficiency of pol beta. This finding is consistent with the view that enzyme preorganization plays a major role in controlling DNA polymerase specific activity.
Collapse
Affiliation(s)
- Václav Martínek
- Department of Chemistry, Loyola University Chicago, Chicago, IL 60626
- Department of Biochemistry, Faculty of Science, Charles University, Albertov 2030, 12840 Prague, Czech Republic
| | - Urban Bren
- Department of Chemistry, Loyola University Chicago, Chicago, IL 60626
- Institute of Physics, Faculty of Mathematics and Physics, Charles University, 12116 Prague, Czech Republic
| | - Myron F. Goodman
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, CA 90089
| | - Jan Florián
- Department of Chemistry, Loyola University Chicago, Chicago, IL 60626
- Institute of Physics, Faculty of Mathematics and Physics, Charles University, 12116 Prague, Czech Republic
- Send correspondence to Jan Florián, Department of Chemistry, Loyola University Chicago, Chicago, IL 60626.
| |
Collapse
|
43
|
Lin P, Pedersen LC, Batra VK, Beard WA, Wilson SH, Pedersen LG. Energy analysis of chemistry for correct insertion by DNA polymerase beta. Proc Natl Acad Sci U S A 2006; 103:13294-9. [PMID: 16938895 PMCID: PMC1569157 DOI: 10.1073/pnas.0606006103] [Citation(s) in RCA: 81] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
X-ray crystallographic structures of human DNA polymerase beta with nonhydrolyzable analogs containing all atoms in the active site required for catalysis provide a secure starting point for a theoretical analysis (quantum mechanics/molecular mechanics) of the mechanism of chemistry without biasing of modeling assumptions as required in previous studies. These structures provide the basis for a detailed quantum mechanics/molecular mechanics study of the path for the complete transfer of a monophosphate nucleoside donor to the sugar acceptor in the active site. The reaction is largely associative with the main energetic step preceded by proton transfer from the terminal primer deoxyribose O3' to Asp-256. The key residues that provide electrostatic stabilization of the transition state are identified and compared with those identified by mutational studies.
Collapse
Affiliation(s)
- Ping Lin
- Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599; and
| | - Lars C. Pedersen
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709-2233
| | - Vinod K. Batra
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709-2233
| | - William A. Beard
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709-2233
| | - Samuel H. Wilson
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709-2233
| | - Lee G. Pedersen
- Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599; and
- Laboratory of Structural Biology, National Institute of Environmental Health Sciences, National Institutes of Health, P.O. Box 12233, Research Triangle Park, NC 27709-2233
- To whom correspondence should be addressed. E-mail:
| |
Collapse
|
44
|
Warshel A, Sharma PK, Kato M, Parson WW. Modeling electrostatic effects in proteins. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2006; 1764:1647-76. [PMID: 17049320 DOI: 10.1016/j.bbapap.2006.08.007] [Citation(s) in RCA: 424] [Impact Index Per Article: 23.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2006] [Revised: 08/17/2006] [Accepted: 08/18/2006] [Indexed: 10/24/2022]
Abstract
Electrostatic energies provide what is perhaps the most effective tool for structure-function correlation of biological molecules. This review considers the current state of simulations of electrostatic energies in macromolecules as well as the early developments of this field. We focus on the relationship between microscopic and macroscopic models, considering the convergence problems of the microscopic models and the fact that the dielectric 'constants' in semimacroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering a wide range of functional properties including pK(a)'s, redox potentials, ion and proton channels, enzyme catalysis, ligand binding and protein stability. We conclude by pointing out that, despite the current problems and the significant misunderstandings in the field, there is an overall progress that should lead eventually to quantitative descriptions of electrostatic effects in proteins and thus to quantitative descriptions of the function of proteins.
Collapse
Affiliation(s)
- Arieh Warshel
- University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA 90089-1062, USA.
| | | | | | | |
Collapse
|
45
|
Griep M, Kotera C, Nelson R, Viljoen H. Kinetics of the DNA polymerase pyrococcus kodakaraensis. Chem Eng Sci 2006. [DOI: 10.1016/j.ces.2005.12.032] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
46
|
Archontis G, Watson KA, Xie Q, Andreou G, Chrysina ED, Zographos SE, Oikonomakos NG, Karplus M. Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues. Proteins 2006; 61:984-98. [PMID: 16245298 DOI: 10.1002/prot.20641] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.
Collapse
Affiliation(s)
- G Archontis
- Department of Physics, University of Cyprus, Cyprus.
| | | | | | | | | | | | | | | |
Collapse
|
47
|
Gorb L, Shishkin O, Leszczynski J. Charges of phosphate groups. A role in stabilization of 2'-deoxyribonucleotides. A DFT investigation. J Biomol Struct Dyn 2005; 22:441-54. [PMID: 15588107 DOI: 10.1080/07391102.2005.10507015] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
We have analyzed the relative stabilities and Gibbs tautomeric free energy for tautomeric transitions of neutral 2'-deoxyribonucleotides and its mono- and di-protonated forms. Geometry optimizations of these nucleic acid constituents have been performed at the DFT/B3LYP level using the standard 6-31G(d) basis set. The prediction of relative stabilities, Gibbs tautomeric free energy has been made at the B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) level of theory. For each nucleoside four major conformers, i.e., north/anti, north/syn, south/anti, and south/syn have been taken into consideration. We have found the substantial effect of the uncompensated charge on the relative stability of 2'-deoxyribonucleotides. In particular, when the charge of 2'-deoxyribonucleotide anions is completely compensated by protons, the syn conformations have been found to be the global minima due to stabilization provided by intramolecular hydrogen bonds. However, the negative charge that appears due to the successive removal of the protons from the phosphate group destabilizes these syn conformations and stabilizes preferably the south/anti conformations (except of 2'-deoxyguanosine phosphate). Only 2'-deoxyribonucleotides, possessing south/anti and north/anti orientations, containing guanine and cytosine can contribute significantly to the rate of spontaneous point mutations due to the formation of biologically relevant amounts of 'rare' tautomers. However, we found strong influence of uncompensated negative charge for 2'-deoxyribonucleotides which possess syn conformations. Finally we have found that the proton transfer could result in the spontaneous change of 2'-deoxyribonucleotides conformations. We conclude that this phenomenon could be considered as a new way for the stabilization of 'rare' isomers for such DNA bases as cytosine and thymine.
Collapse
Affiliation(s)
- Leonid Gorb
- Computational Center for Molecular Structure and Interactions, Department of Chemistry, Jackson State University, P.O. Box 17910, 1325 Lynch Street, Jackson, MS 39217, USA.
| | | | | |
Collapse
|
48
|
Florián J, Goodman MF, Warshel A. Computer simulations of protein functions: searching for the molecular origin of the replication fidelity of DNA polymerases. Proc Natl Acad Sci U S A 2005; 102:6819-24. [PMID: 15863620 PMCID: PMC1100748 DOI: 10.1073/pnas.0408173102] [Citation(s) in RCA: 96] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2004] [Indexed: 11/18/2022] Open
Abstract
The use of computers to simulate the functions of complex biological macromolecules is essential to achieve a microscopic description of biological processes and to model and interpret experimental data. Here we apply theoretical computational approaches to investigate the fidelity of T7 DNA polymerase, divided into discrete steps that include contributions from substrate binding, pK(a) shifts, and rate constants for the PO bond-breaking and bond-making processes. We begin by defining the discrimination between right and wrong nucleotides in terms of the free energy landscape for the dNMP incorporation reaction. We then use the linear response approximation and the empirical valence bond methods to obtain converging results for the contribution of the binding and chemical steps to the overall fidelity. These approaches are successful in reproducing general trends in the observed polymerase incorporation fidelity. The calculations demonstrate the potential for further integration of theoretical and experimental studies to analyze high- and low-fidelity DNA polymerases.
Collapse
Affiliation(s)
- Jan Florián
- Department of Chemistry, Loyola University, Chicago, IL 60626, USA.
| | | | | |
Collapse
|
49
|
Archontis G, Simonson T. Proton binding to proteins: a free-energy component analysis using a dielectric continuum model. Biophys J 2005; 88:3888-904. [PMID: 15821163 PMCID: PMC1305621 DOI: 10.1529/biophysj.104.055996] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Proton binding plays a critical role in protein structure and function. We report pK(a) calculations for three aspartates in two proteins, using a linear response approach, as well as a "standard" Poisson-Boltzmann approach. Averaging over conformations from the two endpoints of the proton-binding reaction, the protein's atomic degrees of freedom are explicitly modeled. Treating macroscopically the protein's electronic polarizability and the solvent, a meaningful model is obtained, without adjustable parameters. It reproduces qualitatively the electrostatic potentials, proton-binding free energies, Marcus reorganization free energies, and pK(a) shifts from explicit solvent molecular dynamics simulations, and the pK(a) shifts from experiment. For thioredoxin Asp-26, which has a large pK(a) upshift, we correctly capture the balance between unfavorable carboxylate desolvation and favorable interactions with a nearby lysine; similarly for RNase A Asp-14, which has a large pK(a) downshift. For the unshifted thioredoxin Asp-20, desolvation by the protein cavity is overestimated by 2.9 pK(a) units; several effects could explain this. "Standard" Poisson-Boltzmann methods sidestep this problem by using a large, ad hoc protein dielectric; but protein charge-charge interactions are then incorrectly downscaled, giving an unbalanced description of the reaction and a large error for the shifted pK(a) values of Asp-26 and Asp-14.
Collapse
|
50
|
Griep M, Whitney S, Nelson M, Viljoen H. DNA polymerase chain reaction: A model of error frequencies and extension rates. AIChE J 2005. [DOI: 10.1002/aic.10604] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|