1
|
Španinger E, Potočnik U, Bren U. Molecular Dynamics Simulations Predict That rSNP Located in the HNF‑1α Gene Promotor Region Linked with MODY3 and Hepatocellular Carcinoma Promotes Stronger Binding of the HNF‑4α Transcription Factor. Biomolecules 2020; 10:biom10121700. [PMID: 33371430 PMCID: PMC7767403 DOI: 10.3390/biom10121700] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Revised: 12/06/2020] [Accepted: 12/18/2020] [Indexed: 12/24/2022] Open
Abstract
Our study aims to investigate the impact of the Maturity-onset diabetes of the young 3 disease-linked rSNP rs35126805 located in the HNF-1α gene promotor on the binding of the transcription factor HNF-4α and consequently on the regulation of HNF-1α gene expression. Our focus is to calculate the change in the binding affinity of the transcription factor HNF-4α to the DNA, caused by the regulatory single nucleotide polymorphism (rSNP) through molecular dynamics simulations and thermodynamic analysis of acquired results. Both root-mean-square difference (RMSD) and the relative binding free energy ΔΔGbind reveal that the HNF-4α binds slightly more strongly to the DNA containing the mutation (rSNP) making the complex more stable/rigid, and thereby influencing the expression of the HNF-1α gene. The resulting disruption of the HNF-4α/HNF-1α pathway is also linked to hepatocellular carcinoma metastasis and enhanced apoptosis in pancreatic cancer cells. To the best of our knowledge, this represents the first study where thermodynamic analysis of the results obtained from molecular dynamics simulations is performed to uncover the influence of rSNP on the protein binding to DNA. Therefore, our approach can be generally applied for studying the impact of regulatory single nucleotide polymorphisms on the binding of transcription factors to the DNA.
Collapse
Affiliation(s)
- Eva Španinger
- Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova ulica 17, SI-2000 Maribor, Slovenia; (E.Š.); (U.P.)
| | - Uroš Potočnik
- Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova ulica 17, SI-2000 Maribor, Slovenia; (E.Š.); (U.P.)
- Faculty of Medicine, University of Maribor, Taborska 8, SI-2000 Maribor, Slovenia
| | - Urban Bren
- Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova ulica 17, SI-2000 Maribor, Slovenia; (E.Š.); (U.P.)
- Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, SI-6000 Koper, Slovenia
- Correspondence: ; Tel.: +386-2-2294-421
| |
Collapse
|
2
|
Kathuria P, Singh P, Sharma P, Wetmore SD. Replication of the Aristolochic Acid I Adenine Adduct (ALI-N6-A) by a Model Translesion Synthesis DNA Polymerase: Structural Insights on the Induction of Transversion Mutations from Molecular Dynamics Simulations. Chem Res Toxicol 2020; 33:2573-2583. [DOI: 10.1021/acs.chemrestox.0c00183] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Preetleen Kathuria
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Prebhleen Singh
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| |
Collapse
|
3
|
NMR and computational methods for molecular resolution of allosteric pathways in enzyme complexes. Biophys Rev 2019; 12:155-174. [PMID: 31838649 DOI: 10.1007/s12551-019-00609-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2019] [Accepted: 12/05/2019] [Indexed: 12/30/2022] Open
Abstract
Allostery is a ubiquitous biological mechanism in which a distant binding site is coupled to and drastically alters the function of a catalytic site in a protein. Allostery provides a high level of spatial and temporal control of the integrity and activity of biomolecular assembles composed of proteins, nucleic acids, or small molecules. Understanding the physical forces that drive allosteric coupling is critical to harnessing this process for use in bioengineering, de novo protein design, and drug discovery. Current microscopic models of allostery highlight the importance of energetics, structural rearrangements, and conformational fluctuations, and in this review, we discuss the synergistic use of solution NMR spectroscopy and computational methods to probe these phenomena in allosteric systems, particularly protein-nucleic acid complexes. This combination of experimental and theoretical techniques facilitates an unparalleled detection of subtle changes to structural and dynamic equilibria in biomolecules with atomic resolution, and we provide a detailed discussion of specialized NMR experiments as well as the complementary methods that provide valuable insight into allosteric pathways in silico. Lastly, we highlight two case studies to demonstrate the adaptability of this approach to enzymes of varying size and mechanistic complexity.
Collapse
|
4
|
Gladovic M, Spaninger E, Bren U. Nucleic Bases Alkylation with Acrylonitrile and Cyanoethylene Oxide: A Computational Study. Chem Res Toxicol 2018; 31:97-104. [DOI: 10.1021/acs.chemrestox.7b00268] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Martin Gladovic
- Faculty
of Chemistry and Chemical Technology, University of Maribor, Smetanova
17, SI-2000 Maribor, Slovenia
- Faculty
of Chemistry and Chemical Technology, University of Ljubljana, Vecna pot
113, SI-1000 Ljubljana, Slovenia
| | - Eva Spaninger
- Faculty
of Chemistry and Chemical Technology, University of Maribor, Smetanova
17, SI-2000 Maribor, Slovenia
| | - Urban Bren
- Faculty
of Chemistry and Chemical Technology, University of Maribor, Smetanova
17, SI-2000 Maribor, Slovenia
- National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
5
|
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
|
6
|
Ogrizek M, Konc J, Bren U, Hodošček M, Janežič D. Role of magnesium ions in the reaction mechanism at the interface between Tm1631 protein and its DNA ligand. Chem Cent J 2016; 10:41. [PMID: 27398092 PMCID: PMC4939058 DOI: 10.1186/s13065-016-0188-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2016] [Accepted: 06/27/2016] [Indexed: 12/24/2022] Open
Abstract
A protein, Tm1631 from the hyperthermophilic organism Thermotoga maritima belongs to a domain of unknown function protein family. It was predicted that Tm1631 binds with the DNA and that the Tm1631–DNA complex is an endonuclease repair system with a DNA repair function (Konc et al. PLoS Comput Biol 9(11): e1003341, 2013). We observed that the severely bent, strained DNA binds to the protein for the entire 90 ns of classical molecular dynamics (MD) performed; we could observe no significant changes in the most distorted region of the DNA, where the cleavage of phosphodiester bond occurs. In this article, we modeled the reaction mechanism at the interface between Tm1631 and its proposed ligand, the DNA molecule, focusing on cleavage of the phosphodiester bond. After addition of two Mg2+ ions to the reaction center and extension of classical MD by 50 ns (totaling 140 ns), the DNA ligand stayed bolted to the protein. Results from density functional theory quantum mechanics/molecular mechanics (QM/MM) calculations suggest that the reaction is analogous to known endonuclease mechanisms: an enzyme reaction mechanism with two Mg2+ ions in the reaction center and a pentacovalent intermediate. The minimum energy pathway profile shows that the phosphodiester bond cleavage step of the reaction is kinetically controlled and not thermodynamically because of a lack of any energy barrier above the accuracy of the energy profile calculation. The role of ions is shown by comparing the results with the reaction mechanisms in the absence of the Mg2+ ions where there is a significantly higher reaction barrier than in the presence of the Mg2+ ions.A protein, Tm1631 from the hyperthermophilic organism Thermotoga maritima belongs to a domain of unknown function protein family. We modeled the reaction mechanism at the interface between Tm1631 and its proposed ligand, the DNA molecule, focusing on cleavage of the phosphodiester bond ![]()
Collapse
Affiliation(s)
- Mitja Ogrizek
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
| | - Janez Konc
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia ; Laboratory for Physical Chemistry and Thermodynamics, Faculty of Chemistry and Chemical Technology, University of Maribor, Smetanova ulica 17, 2000 Maribor, Slovenia ; Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, 6000 Koper, Slovenia
| | - Urban Bren
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia ; Laboratory for Physical Chemistry and Thermodynamics, Faculty of Chemistry and Chemical Technology, University of Maribor, Smetanova ulica 17, 2000 Maribor, Slovenia ; Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, 6000 Koper, Slovenia
| | - Milan Hodošček
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
| | - Dušanka Janežič
- Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, 6000 Koper, Slovenia
| |
Collapse
|
7
|
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
|
8
|
Jeřábek P, Florián J, Stiborová M, Martínek V. Flexible docking-based molecular dynamics/steered molecular dynamics calculations of protein-protein contacts in a complex of cytochrome P450 1A2 with cytochrome b5. Biochemistry 2014; 53:6695-705. [PMID: 25313797 DOI: 10.1021/bi500814t] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Formation of transient complexes of cytochrome P450 (P450) with another protein of the endoplasmic reticulum membrane, cytochrome b5 (cyt b5), dictates the catalytic activities of several P450s. Therefore, we examined formation and binding modes of the complex of human P450 1A2 with cyt b5. Docking of soluble domains of these proteins was performed using an information-driven flexible docking approach implemented in HADDOCK. Stabilities of the five unique binding modes of the P450 1A2-cyt b5 complex yielded by HADDOCK were evaluated using explicit 10 ns molecular dynamics (MD) simulations in aqueous solution. Further, steered MD was used to compare the stability of the individual P450 1A2-cyt b5 binding modes. The best binding mode was characterized by a T-shaped mutual orientation of the porphyrin rings and a 10.7 Å distance between the two redox centers, thus satisfying the condition for a fast electron transfer. Mutagenesis studies and chemical cross-linking, which, in the absence of crystal structures, were previously used to deduce specific P450-cyt b5 interactions, indicated that the negatively charged convex surface of cyt b5 binds to the positively charged concave surface of P450. Our simulations further elaborate structural details of this interface, including nine ion pairs between R95, R100, R138, R362, K442, K455, and K465 side chains of P450 1A2 and E42, E43, E49, D65, D71, and heme propionates of cyt b5. The universal heme-centric system of internal coordinates was proposed to facilitate consistent classification of the orientation of the two porphyrins in any protein complex.
Collapse
Affiliation(s)
- Petr Jeřábek
- Department of Biochemistry, Faculty of Science, Charles University in Prague , Albertov 2030, 128 43 Prague 2, Czech Republic
| | | | | | | |
Collapse
|
9
|
C(α) torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation. J Mol Model 2014; 20:2196. [PMID: 24728650 DOI: 10.1007/s00894-014-2196-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2014] [Accepted: 03/02/2014] [Indexed: 02/02/2023]
Abstract
Given the increasing complexity of simulated molecular systems, and the fact that simulation times have now reached milliseconds to seconds, immense amounts of data (in the gigabyte to terabyte range) are produced in current molecular dynamics simulations. Manual analysis of these data is a very time-consuming task, and important events that lead from one intermediate structure to another can become occluded in the noise resulting from random thermal fluctuations. To overcome these problems and facilitate a semi-automated data analysis, we introduce in this work a measure based on C(α) torsion angles: torsion angles formed by four consecutive C(α) atoms. This measure describes changes in the backbones of large systems on a residual length scale (i.e., a small number of residues at a time). Cluster analysis of individual C(α) torsion angles and its fuzzification led to continuous time patches representing (meta)stable conformations and to the identification of events acting as transitions between these conformations. The importance of a change in torsion angle to structural integrity is assessed by comparing this change to the average fluctuations in the same torsion angle over the complete simulation. Using this novel measure in combination with other measures such as the root mean square deviation (RMSD) and time series of distance measures, we performed an in-depth analysis of a simulation of the open form of DNA polymerase I. The times at which major conformational changes occur and the most important parts of the molecule and their interrelations were pinpointed in this analysis. The simultaneous determination of the time points and localizations of major events is a significant advantage of the new bottom-up approach presented here, as compared to many other (top-down) approaches in which only the similarity of the complete structure is analyzed.
Collapse
|
10
|
Dubey KD, Chaubey AK, Ojha RP. Stability of trimeric DENV envelope protein at low and neutral pH: An insight from MD study. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:53-64. [DOI: 10.1016/j.bbapap.2012.08.014] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/27/2012] [Revised: 08/13/2012] [Accepted: 08/14/2012] [Indexed: 01/29/2023]
|
11
|
Bren U, Janežič D. Individual degrees of freedom and the solvation properties of water. J Chem Phys 2012; 137:024108. [PMID: 22803529 DOI: 10.1063/1.4732514] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using molecular dynamics simulations in conjunction with home-developed Split Integration Symplectic Method we effectively decouple individual degrees of freedom of water molecules and connect them to corresponding thermostats. In this way, we facilitate elucidation of structural, dynamical, spectral, and hydration properties of bulk water at any given combination of rotational, translational, and vibrational temperatures. Elevated rotational temperature of the water medium is found to severely hinder hydration of polar molecules, to affect hydration of ionic species in a nonmonotonous way and to somewhat improve hydration of nonpolar species. As proteins consist of charged, polar, and nonpolar amino-acid residues, the developed methodology is also applied to critically evaluate the hypothesis that the overall decrease in protein hydration and the change in the subtle balance between hydration of various types of amino-acid residues provide a plausible physical mechanism through which microwaves enhance aberrant protein folding and aggregation.
Collapse
Affiliation(s)
- Urban Bren
- Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1001 Ljubljana, Slovenia
| | | |
Collapse
|
12
|
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
|
13
|
Soares ROS, Caliri A. Stereochemical features of the envelope protein Domain III of dengue virus reveals putative antigenic site in the five-fold symmetry axis. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2012; 1834:221-30. [PMID: 23009809 DOI: 10.1016/j.bbapap.2012.09.007] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2012] [Revised: 09/07/2012] [Accepted: 09/11/2012] [Indexed: 11/18/2022]
Abstract
We bring to attention a characteristic parasitic pattern present in the dengue virus: it undergoes several intensive thermodynamic variations due to host environmental changes, from a vector's digestive tract, through the human bloodstream and intracellular medium. Comparatively, among the known dengue serotypes, we evaluate the effects that these medium variations may induce to the overall structural characteristics of the Domain III of the envelope (E) protein, checking for stereochemical congruences that could lead to the identification of immunologic relevant regions. We used molecular dynamics and principal component analysis to study the protein in solution, for all four dengue serotypes, under distinct pH and temperature. We stated that, while the core of Domain III is remarkably rigid and effectively unaffected by most of the mentioned intensive variations, the loops account for major and distinguishable flexibilities. Therefore, the rigidity of the Domain III core provides a foothold that projects specifically two of these high flexible loop regions towards the inner face of the envelope pores, which are found at every five-fold symmetry axis of the icosahedron-shaped mature virus. These loops bear a remarkable low identity though with high occurrence of ionizable residues, including histidines. Such stereochemical properties can provide very particular serotype-specific electrostatic surface patterns, suggesting a viral fingerprint region, on which other specific molecules and ions can establish chemical interactions in an induced fit mechanism. We assert that the proposed regions share enough relevant features to qualify for further immunologic and pharmacologic essays, such as target peptide synthesis and phage display using dengue patients' sera.
Collapse
Affiliation(s)
- R O S Soares
- Departamento de Física e Química, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. do Café, S/N. 14040-903, Ribeirão Preto, São Paulo, Brazil.
| | | |
Collapse
|
14
|
Figueroa H, Peddi D, Osborne JM, Wilson BM, Pesaru RR, Kurva B, Ramaraju S, Milletti MC, Heyl DL. Modeling the Interface between Islet Amyloid Polypeptide and Insulin-Based Aggregation Inhibitors: Correlation to Aggregation Kinetics and Membrane Damage. J Chem Inf Model 2012; 52:1298-307. [DOI: 10.1021/ci300119c] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Hector Figueroa
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Durgaprasad Peddi
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Joshua M. Osborne
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Brenan M. Wilson
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Ranadheer Reddy Pesaru
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Balakrishna Kurva
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Swathi Ramaraju
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Maria C. Milletti
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| | - Deborah L. Heyl
- Chemistry Department, Eastern Michigan University, Ypsilanti, Michigan 48197, United States
| |
Collapse
|
15
|
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
|
16
|
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
|
17
|
|
18
|
Fulle S, Christ NA, Kestner E, Gohlke H. HIV-1 TAR RNA spontaneously undergoes relevant apo-to-holo conformational transitions in molecular dynamics and constrained geometrical simulations. J Chem Inf Model 2010; 50:1489-501. [PMID: 20726603 DOI: 10.1021/ci100101w] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We report all-atom molecular dynamics and replica exchange molecular dynamics simulations on the unbound human immunodeficiency virus type-1 (HIV-1) transactivation responsive region (TAR) RNA structure and three TAR RNA structures in bound conformations of, in total, approximately 250 ns length. We compare the extent of observed conformational sampling with that of the conceptually simpler and computationally much cheaper constrained geometrical simulation approach framework rigidity optimized dynamic algorithm (FRODA). Atomic fluctuations obtained by replica-exchange molecular dynamics (REMD) simulations agree quantitatively with those obtained by molecular dynamics (MD) and FRODA simulations for the unbound TAR structure. Regarding the stereochemical quality of the generated conformations, backbone torsion angles and puckering modes of the sugar-phosphate backbone were reproduced equally well by MD and REMD simulations, but further improvement is needed in the case of FRODA simulations. Essential dynamics analysis reveals that all three simulation approaches show a tendency to sample bound conformations when starting from the unbound TAR structure, with MD and REMD simulations being superior with respect to FRODA. These results are consistent with the experimental view that bound TAR RNA conformations are transiently sampled in the free ensemble, following a conformation selection model. The simulation-generated TAR RNA conformations have been successfully used as receptor structures for docking. This finding has important implications for RNA-ligand docking in that docking into an ensemble of simulation-generated RNA structures is shown to be a valuable means to cope with large apo-to-holo conformational transitions of the receptor structure.
Collapse
Affiliation(s)
- Simone Fulle
- Department of Biological Sciences, Goethe-University, Frankfurt, Germany
| | | | | | | |
Collapse
|
19
|
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
|
20
|
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
|
21
|
Udommaneethanakit T, Rungrotmongkol T, Bren U, Frecer V, Stanislav M. Dynamic behavior of avian influenza A virus neuraminidase subtype H5N1 in complex with oseltamivir, zanamivir, peramivir, and their phosphonate analogues. J Chem Inf Model 2009; 49:2323-32. [PMID: 19780597 DOI: 10.1021/ci900277r] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The outbreak of avian influenza A subtype H5N1 virus has raised a global concern for both animal as well as human health. Recently, drug resistance in H5N1 infections has been widely reported due to neuraminidase mutations. Consequently, the understanding of inhibitor-neuraminidase interactions at the molecular level represents the main goal of our study. Molecular dynamics simulations were carried out for the neuraminidase N1 in complex with six inhibitors--oseltamivir, zanamivir, peramivir, and their phosphonate analogues. Molecular dynamics trajectories were extensively analyzed in terms of important interactions between inhibitors and the enzyme target. Results show that open and closed forms (defined by the relative position of the flexible 150-loop) of neuraminidase N1 interchange during the course of 20 ns molecular dynamics simulation of the protein-inhibitor complexes. Reported free energies of closing indicate that the carboxylate inhibitors prefer the closed form more than their phosphonate analogues. This can be understood in view of the negative total charge (-1 e0) of the phosphonate inhibitors, which repels the Asp151 residue of the loop away from the inhibitor and drives the complex into the open form. Obtained results constitute new valuable information to assist further drug development of inhibitors against the H5N1 avian influenza A virus and could also inspire similar studies for other systems of the influenza family such as the 2009 influenza A (H1N1) virus.
Collapse
Affiliation(s)
- Thanyarat Udommaneethanakit
- International Centre for Science and High Technology, UNIDO, AREA Science Park, Padriciano 99, I-34012 Trieste, Italy
| | | | | | | | | |
Collapse
|
22
|
Yamtich J, Sweasy JB. DNA polymerase family X: function, structure, and cellular roles. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2009; 1804:1136-50. [PMID: 19631767 DOI: 10.1016/j.bbapap.2009.07.008] [Citation(s) in RCA: 109] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/05/2009] [Revised: 07/02/2009] [Accepted: 07/03/2009] [Indexed: 10/20/2022]
Abstract
The X family of DNA polymerases in eukaryotic cells consists of terminal transferase and DNA polymerases beta, lambda, and mu. These enzymes have similar structural portraits, yet different biochemical properties, especially in their interactions with DNA. None of these enzymes possesses a proofreading subdomain, and their intrinsic fidelity of DNA synthesis is much lower than that of a polymerase that functions in cellular DNA replication. In this review, we discuss the similarities and differences of three members of Family X: polymerases beta, lambda, and mu. We focus on biochemical mechanisms, structural variation, fidelity and lesion bypass mechanisms, and cellular roles. Remarkably, although these enzymes have similar three-dimensional structures, their biochemical properties and cellular functions differ in important ways that impact cellular function.
Collapse
Affiliation(s)
- Jennifer Yamtich
- Department of Therapeutic Radiology, Yale University School of Medicine, 333 Cedar Street, New Haven, CT 06520, USA
| | | |
Collapse
|
23
|
Vardi-Kilshtain A, Roca M, Warshel A. The empirical valence bond as an effective strategy for computer-aided enzyme design. Biotechnol J 2009; 4:495-500. [PMID: 19229886 DOI: 10.1002/biot.200800299] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The ability of the empirical valence bond (EVB) to be used in screening active site residues in enzyme design is explored in a preliminary way. This validation is done by comparing the ability of this approach to evaluate the catalytic contributions of various residues in chorismate mutase. It is demonstrated that the EVB model can serve as an accurate tool in the final stages of computer-aided enzyme design (CAED). The ability of the model to predict quantitatively the catalytic power of enzymes should augment the capacity of current approaches for enzyme design.
Collapse
|
24
|
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
|
25
|
Sucato CA, Upton TG, Kashemirov BA, Osuna J, Oertell K, Beard WA, Wilson SH, Florián J, Warshel A, McKenna CE, Goodman MF. DNA polymerase beta fidelity: halomethylene-modified leaving groups in pre-steady-state kinetic analysis reveal differences at the chemical transition state. Biochemistry 2007; 47:870-9. [PMID: 18161950 DOI: 10.1021/bi7014162] [Citation(s) in RCA: 73] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The mechanism of DNA polymerase beta-catalyzed nucleotidyl transfer consists of chemical steps involving primer 3' OH deprotonation, nucleophilic attack, and pyrophosphate leaving-group elimination, preceded by dNTP binding which induces a large-amplitude conformational change for Watson-Crick nascent base pairs. Ambiguity in the nature of the rate-limiting step and active-site structural differences between correct and incorrect base-paired transition states remain obstacles to understanding DNA replication fidelity. Analogues of dGTP where the beta-gamma bridging oxygen is replaced with fluorine-substituted methylene groups have been shown to probe the contribution of leaving-group elimination to the overall catalytic rate (Biochemistry 46, 461-471). Here, the analysis is expanded substantially to include a broad range of halogen substituents with disparate steric and electronic properties. Evaluation of linear free energy relationships for incorporation of dGTP analogues opposite either template base C or T reveals a strong correlation of log(kpol) to leaving group pKa. Significantly different kpol behavior is observed with a subset of the analogues, with magnitude dependent on the identity of the nascent base pair. This observation, and the absence of an analogous effect on ground state analogue binding (Kd values), points to active-site structural differences at the chemical transition state. Reduced catalysis with bulky halo-containing substrates is manifested in the fidelity of T-G incorporation, where the CCl2-bridging analogue shows a 27-fold increase in fidelity over the natural dGTP. Solvent pH and deuterium isotope-effect data are also used to evaluate mechanistic differences between correct and mispaired incorporation.
Collapse
Affiliation(s)
- Christopher A Sucato
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|