1
|
The Role of Active-Site Plasticity in Damaged-Nucleotide Recognition by Human Apurinic/Apyrimidinic Endonuclease APE1. Molecules 2020; 25:molecules25173940. [PMID: 32872297 PMCID: PMC7504742 DOI: 10.3390/molecules25173940] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2020] [Revised: 08/25/2020] [Accepted: 08/26/2020] [Indexed: 12/19/2022] Open
Abstract
Human apurinic/apyrimidinic (AP) endonuclease APE1 hydrolyzes phosphodiester bonds on the 5′ side of an AP-site, and some damaged nucleotides such as 1,N6-ethenoadenosine (εA), α-adenosine (αA), and 5,6-dihydrouridine (DHU). To investigate the mechanism behind the broad substrate specificity of APE1, we analyzed pre-steady-state kinetics of conformational changes in DNA and the enzyme during DNA binding and damage recognition. Molecular dynamics simulations of APE1 complexes with one of damaged DNA duplexes containing εA, αA, DHU, or an F-site (a stable analog of an AP-site) revealed the involvement of residues Asn229, Thr233, and Glu236 in the mechanism of DNA lesion recognition. The results suggested that processing of an AP-site proceeds faster in comparison with nucleotide incision repair substrates because eversion of a small abasic site and its insertion into the active site do not include any unfavorable interactions, whereas the insertion of any target nucleotide containing a damaged base into the APE1 active site is sterically hindered. Destabilization of the α-helix containing Thr233 and Glu236 via a loss of the interaction between these residues increased the plasticity of the damaged-nucleotide binding pocket and the ability to accommodate structurally different damaged nucleotides. Nonetheless, the optimal location of εA or αA in the binding pocket does not correspond to the optimal conformation of catalytic amino acid residues, thereby significantly decreasing the cleavage efficacy for these substrates.
Collapse
|
2
|
Substrate specificity of human apurinic/apyrimidinic endonuclease APE1 in the nucleotide incision repair pathway. Nucleic Acids Res 2019; 46:11454-11465. [PMID: 30329131 PMCID: PMC6265485 DOI: 10.1093/nar/gky912] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2018] [Accepted: 10/10/2018] [Indexed: 12/14/2022] Open
Abstract
Human apurinic/apyrimidinic (AP) endonuclease APE1 catalyses the hydrolysis of phosphodiester bonds on the 5′ side of an AP-site (in the base excision repair pathway) and of some damaged nucleotides (in the nucleotide incision repair pathway). The range of substrate specificity includes structurally unrelated damaged nucleotides. Here, to examine the mechanism of broad substrate specificity of APE1, we performed pulsed electron–electron double resonance (PELDOR) spectroscopy and pre-steady-state kinetic analysis with Förster resonance energy transfer (FRET) detection of DNA conformational changes during DNA binding and lesion recognition. Equilibrium PELDOR and kinetic FRET data revealed that DNA binding by APE1 leads to noticeable damage-dependent bending of a DNA duplex. Molecular dynamics simulations showed that the damaged nucleotide is everted from the DNA helix and placed into the enzyme’s binding pocket, which is formed by Asn-174, Asn-212, Asn-229, Ala-230, Phe-266 and Trp-280. Nevertheless, no damage-specific contacts were detected between these amino acid residues in the active site of the enzyme and model damaged substrates containing 1,N6-ethenoadenosine, α-adenosine, 5,6-dihydrouridine or F-site. These data suggest that the substrate specificity of APE1 is controlled by the ability of a damaged nucleotide to flip out from the DNA duplex in response to an enzyme-induced DNA distortion.
Collapse
|
3
|
Molecular dynamics simulation of the opposite-base preference and interactions in the active site of formamidopyrimidine-DNA glycosylase. BMC STRUCTURAL BIOLOGY 2017; 17:5. [PMID: 28482831 PMCID: PMC5422863 DOI: 10.1186/s12900-017-0075-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/12/2016] [Accepted: 04/20/2017] [Indexed: 01/20/2023]
Abstract
Background Formamidopyrimidine-DNA glycosylase (Fpg) removes abundant pre-mutagenic 8-oxoguanine (oxoG) bases from DNA through nucleophilic attack of its N-terminal proline at C1′ of the damaged nucleotide. Since oxoG efficiently pairs with both C and A, Fpg must excise oxoG from pairs with C but not with A, otherwise a mutation occurs. The crystal structures of several Fpg–DNA complexes have been solved, yet no structure with A opposite the lesion is available. Results Here we use molecular dynamic simulation to model interactions in the pre-catalytic complex of Lactococcus lactis Fpg with DNA containing oxoG opposite C or A, the latter in either syn or anti conformation. The catalytic dyad, Pro1–Glu2, was modeled in all four possible protonation states. Only one transition was observed in the experimental reaction rate pH dependence plots, and Glu2 kept the same set of interactions regardless of its protonation state, suggesting that it does not limit the reaction rate. The adenine base opposite oxoG was highly distorting for the adjacent nucleotides: in the more stable syn models it formed non-canonical bonds with out-of-register nucleotides in both the damaged and the complementary strand, whereas in the anti models the adenine either formed non-canonical bonds or was expelled into the major groove. The side chains of Arg109 and Phe111 that Fpg inserts into DNA to maintain its kinked conformation tended to withdraw from their positions if A was opposite to the lesion. The region showing the largest differences in the dynamics between oxoG:C and oxoG:A substrates was unexpectedly remote from the active site, located near the linker joining the two domains of Fpg. This region was also highly conserved among 124 analyzed Fpg sequences. Three sites trapping water molecules through multiple bonds were identified on the protein–DNA interface, apparently helping to maintain enzyme-induced DNA distortion and participating in oxoG recognition. Conclusion Overall, the discrimination against A opposite to the lesion seems to be due to incorrect DNA distortion around the lesion-containing base pair and, possibly, to gross movement of protein domains connected by the linker. Electronic supplementary material The online version of this article (doi:10.1186/s12900-017-0075-y) contains supplementary material, which is available to authorized users.
Collapse
|
4
|
Effects of mono- and divalent metal ions on DNA binding and catalysis of human apurinic/apyrimidinic endonuclease 1. MOLECULAR BIOSYSTEMS 2016; 12:1527-39. [PMID: 27063150 DOI: 10.1039/c6mb00128a] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Here, we used stopped-flow fluorescence techniques to conduct a comparative kinetic analysis of the conformational transitions in human apurinic/apyrimidinic endonuclease 1 (APE1) and in DNA containing an abasic site in the course of their interaction. Effects of monovalent (K(+)) and divalent (Mg(2+), Mn(2+), Ca(2+), Zn(2+), Cu(2+), and Ni(2+)) metal ions on DNA binding and catalytic stages were studied. It was shown that the first step of substrate binding (corresponding to formation of a primary enzyme-substrate complex) does not depend on the concentration (0.05-5.0 mM) or the nature of divalent metal ions. In contrast, the initial DNA binding efficiency significantly decreased at a high concentration (5-250 mM) of monovalent K(+) ions, indicating the involvement of electrostatic interactions in this stage. It was also shown that Cu(2+) ions abrogated the DNA binding ability of APE1, possibly, due to a strong interaction with DNA bases and the sugar-phosphate backbone. In the case of Ca(2+) ions, the catalytic activity of APE1 was lost completely with retention of binding potential. Thus, the enzymatic activity of APE1 is increased in the order Zn(2+) < Ni(2+) < Mn(2+) < Mg(2+). Circular dichroism spectra and calculation of the contact area between APE1 and DNA reveal that Mg(2+) ions stabilize the protein structure and the enzyme-substrate complex.
Collapse
|
5
|
55 MDTRA: unified approach to molecular dynamics data preparation and analysis. J Biomol Struct Dyn 2015. [DOI: 10.1080/07391102.2015.1032671] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
6
|
Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species. BMC STRUCTURAL BIOLOGY 2014; 14:23. [PMID: 25315147 PMCID: PMC4209456 DOI: 10.1186/s12900-014-0023-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/06/2014] [Accepted: 10/03/2014] [Indexed: 11/13/2022]
Abstract
Background The identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures. Results The current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified. Conclusions The results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.
Collapse
|
7
|
Thermodynamics of the DNA damage repair steps of human 8-oxoguanine DNA glycosylase. PLoS One 2014; 9:e98495. [PMID: 24911585 PMCID: PMC4049573 DOI: 10.1371/journal.pone.0098495] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2014] [Accepted: 05/03/2014] [Indexed: 11/19/2022] Open
Abstract
Human 8-oxoguanine DNA glycosylase (hOGG1) is a key enzyme responsible for initiating the base excision repair of 7,8-dihydro-8-oxoguanosine (oxoG). In this study a thermodynamic analysis of the interaction of hOGG1 with specific and non-specific DNA-substrates is performed based on stopped-flow kinetic data. The standard Gibbs energies, enthalpies and entropies of specific stages of the repair process were determined via kinetic measurements over a temperature range using the van’t Hoff approach. The three steps which are accompanied with changes in the DNA conformations were detected via 2-aminopurine fluorescence in the process of binding and recognition of damaged oxoG base by hOGG1. The thermodynamic analysis has demonstrated that the initial step of the DNA substrates binding is mainly governed by energy due to favorable interactions in the process of formation of the recognition contacts, which results in negative enthalpy change, as well as due to partial desolvation of the surface between the DNA and enzyme, which results in positive entropy change. Discrimination of non-specific G base versus specific oxoG base is occurring in the second step of the oxoG-substrate binding. This step requires energy consumption which is compensated by the positive entropy contribution. The third binding step is the final adjustment of the enzyme/substrate complex to achieve the catalytically competent state which is characterized by large endothermicity compensated by a significant increase of entropy originated from the dehydration of the DNA grooves.
Collapse
|
8
|
DNA damage processing by human 8-oxoguanine-DNA glycosylase mutants with the occluded active site. J Biol Chem 2013; 288:28936-47. [PMID: 23955443 DOI: 10.1074/jbc.m113.487322] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
8-Oxoguanine-DNA glycosylase (OGG1) removes premutagenic lesion 8-oxoguanine (8-oxo-G) from DNA and then nicks the nascent abasic (apurinic/apyrimidinic) site by β-elimination. Although the structure of OGG1 bound to damaged DNA is known, the dynamic aspects of 8-oxo-G recognition are not well understood. To comprehend the mechanisms of substrate recognition and processing, we have constructed OGG1 mutants with the active site occluded by replacement of Cys-253, which forms a wall of the base-binding pocket, with bulky leucine or isoleucine. The conformational dynamics of OGG1 mutants were characterized by single-turnover kinetics and stopped-flow kinetics with fluorescent detection. Additionally, the conformational mobility of wild type and the mutant OGG1 substrate complex was assessed using molecular dynamics simulations. Although pocket occlusion distorted the active site and greatly decreased the catalytic activity of OGG1, it did not fully prevent processing of 8-oxo-G and apurinic/apyrimidinic sites. Both mutants were notably stimulated in the presence of free 8-bromoguanine, indicating that this base can bind to the distorted OGG1 and facilitate β-elimination. The results agree with the concept of enzyme plasticity, suggesting that the active site of OGG1 is flexible enough to compensate partially for distortions caused by mutation.
Collapse
|
9
|
109 Molecular dynamics done quick: efficient trajectory analysis software. J Biomol Struct Dyn 2013. [DOI: 10.1080/07391102.2013.786351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
10
|
Thermodynamics of the multi-stage DNA lesion recognition and repair by formamidopyrimidine-DNA glycosylase using pyrrolocytosine fluorescence--stopped-flow pre-steady-state kinetics. Nucleic Acids Res 2012; 40:7384-92. [PMID: 22584623 PMCID: PMC3424566 DOI: 10.1093/nar/gks423] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Formamidopyrimidine-DNA glycosylase, Fpg protein from Escherichia coli, initiates base excision repair in DNA by removing a wide variety of oxidized lesions. In this study, we perform thermodynamic analysis of the multi-stage interaction of Fpg with specific DNA-substrates containing 7,8-dihydro-8-oxoguanosine (oxoG), or tetrahydrofuran (THF, an uncleavable abasic site analog) and non-specific (G) DNA-ligand based on stopped-flow kinetic data. Pyrrolocytosine, highly fluorescent analog of the natural nucleobase cytosine, is used to record multi-stage DNA lesion recognition and repair kinetics over a temperature range (10–30°C). The kinetic data were used to obtain the standard Gibbs energy, enthalpy and entropy of the specific stages using van’t Hoff approach. The data suggest that not only enthalpy-driven exothermic oxoG recognition, but also the desolvation-accompanied entropy-driven enzyme-substrate complex adjustment into the catalytically active state play equally important roles in the overall process.
Collapse
|
11
|
Adenine and guanine recognition of stop codon is mediated by different N domain conformations of translation termination factor eRF1. Nucleic Acids Res 2011; 39:7134-46. [PMID: 21602268 PMCID: PMC3167606 DOI: 10.1093/nar/gkr376] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Positioning of release factor eRF1 toward adenines and the ribose-phosphate backbone of the UAAA stop signal in the ribosomal decoding site was studied using messenger RNA (mRNA) analogs containing stop signal UAA/UAAA and a photoactivatable cross-linker at definite locations. The human eRF1 peptides cross-linked to these analogs were identified. Cross-linkers on the adenines at the 2nd, 3rd or 4th position modified eRF1 near the conserved YxCxxxF loop (positions 125–131 in the N domain), but cross-linker at the 4th position mainly modified the tripeptide 26-AAR-28. This tripeptide cross-linked also with derivatized 3′-phosphate of UAA, while the same cross-linker at the 3′-phosphate of UAAA modified both the 26–28 and 67–73 fragments. A comparison of the results with those obtained earlier with mRNA analogs bearing a similar cross-linker at the guanines indicates that positioning of eRF1 toward adenines and guanines of stop signals in the 80S termination complex is different. Molecular modeling of eRF1 in the 80S termination complex showed that eRF1 fragments neighboring guanines and adenines of stop signals are compatible with different N domain conformations of eRF1. These conformations vary by positioning of stop signal purines toward the universally conserved dipeptide 31-GT-32, which neighbors guanines but is oriented more distantly from adenines.
Collapse
|
12
|
Three distinct peptides from the N domain of translation termination factor eRF1 surround stop codon in the ribosome. RNA (NEW YORK, N.Y.) 2010; 16:1902-14. [PMID: 20688868 PMCID: PMC2941099 DOI: 10.1261/rna.2066910] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2009] [Accepted: 06/27/2010] [Indexed: 05/07/2023]
Abstract
To study positioning of the polypeptide release factor eRF1 toward a stop signal in the ribosomal decoding site, we applied photoactivatable mRNA analogs, derivatives of oligoribonucleotides. The human eRF1 peptides cross-linked to these short mRNAs were identified. Cross-linkers on the guanines at the second, third, and fourth stop signal positions modified fragment 31-33, and to lesser extent amino acids within region 121-131 (the "YxCxxxF loop") in the N domain. Hence, both regions are involved in the recognition of the purines. A cross-linker at the first uridine of the stop codon modifies Val66 near the NIKS loop (positions 61-64), and this region is important for recognition of the first uridine of stop codons. Since the N domain distinct regions of eRF1 are involved in a stop-codon decoding, the eRF1 decoding site is discontinuous and is not of "protein anticodon" type. By molecular modeling, the eRF1 molecule can be fitted to the A site proximal to the P-site-bound tRNA and to a stop codon in mRNA via a large conformational change to one of its three domains. In the simulated eRF1 conformation, the YxCxxxF motif and positions 31-33 are very close to a stop codon, which becomes also proximal to several parts of the C domain. Thus, in the A-site-bound state, the eRF1 conformation significantly differs from those in crystals and solution. The model suggested for eRF1 conformation in the ribosomal A site and cross-linking data are compatible.
Collapse
MESH Headings
- Base Sequence
- Codon, Terminator/genetics
- Codon, Terminator/metabolism
- Cross-Linking Reagents
- Humans
- In Vitro Techniques
- Models, Molecular
- Mutagenesis, Site-Directed
- Peptide Chain Termination, Translational
- Peptide Fragments/chemistry
- Peptide Fragments/genetics
- Peptide Fragments/metabolism
- Peptide Mapping
- Peptide Termination Factors/chemistry
- Peptide Termination Factors/genetics
- Peptide Termination Factors/metabolism
- Protein Conformation
- Protein Structure, Tertiary
- RNA, Messenger/chemistry
- RNA, Messenger/genetics
- RNA, Messenger/metabolism
- RNA, Transfer/chemistry
- RNA, Transfer/genetics
- RNA, Transfer/metabolism
- Recombinant Proteins/chemistry
- Recombinant Proteins/genetics
- Recombinant Proteins/metabolism
- Ribosomes/genetics
- Ribosomes/metabolism
Collapse
|
13
|
Pre-steady-state kinetic study of substrate specificity of Escherichia coli formamidopyrimidine--DNA glycosylase. Biochemistry 2007; 46:424-35. [PMID: 17209553 DOI: 10.1021/bi060787r] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Formamidopyrimidine-DNA glycosylase (Fpg) is responsible for removal of 8-oxoguanine (8-oxoG) and other oxidized purine lesions from DNA and can also excise some oxidatively modified pyrimidines [such as dihydrouracil (DHU)]. Fpg is also specific for a base opposite the lesion, efficiently excising 8-oxoG paired with C but not with A. We have applied stopped-flow kinetics using intrinsic tryptophan fluorescence of the enzyme and fluorescence of 2-aminopurine-labeled DNA to analyze the conformational dynamics of Escherichia coli Fpg during processing of good substrates (8-oxoG.C), poor substrates (8-oxoG.A), and substrates of unclear specificity (such as DHU and 8-oxoG opposite T or G). The analysis of fluorescence traces allows us to conclude that when the enzyme encounters its true substrate, 8-oxoG.C, the complex enters the productive catalytic reaction after approximately 50 ms, partitioning the substrate away from the competing dissociation process, while poor substrates linger in the initial encounter complex for longer. Several intermediate ES complexes were attributed to different structures that exist along the reaction pathway. A likely sequence of events is that the damaged base is first destabilized by the enzyme binding and then everted from DNA, followed by insertion of several amino acid residues into DNA and isomerization of the enzyme into a pre-excision complex. We conclude that rejection of the incorrect substrates occurs mostly at the early stage of formation of the pre-eversion recognition complex, supporting the role of indirect readout in damage recognition.
Collapse
|
14
|
Free energies of protein decoys provide insight into determinants of protein stability. Protein Sci 2001; 10:2498-506. [PMID: 11714917 PMCID: PMC2374037 DOI: 10.1110/ps.15501] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2001] [Revised: 08/29/2001] [Accepted: 09/06/2001] [Indexed: 10/21/2022]
Abstract
We have calculated the stability of decoy structures of several proteins (from the CASP3 models and the Park and Levitt decoy set) relative to the native structures. The calculations were performed with the force field-consistent ES/IS method, in which an implicit solvent (IS) model is used to calculate the average solvation free energy for snapshots from explicit simulations (ESs). The conformational free energy is obtained by adding the internal energy of the solute from the ESs and an entropic term estimated from the covariance positional fluctuation matrix. The set of atomic Born radii and the cavity-surface free energy coefficient used in the implicit model has been optimized to be consistent with the all-atom force field used in the ESs (cedar/gromos with simple point charge (SPC) water model). The decoys are found to have a consistently higher free energy than that of the native structure; the gap between the native structure and the best decoy varies between 10 and 15 kcal/mole, on the order of the free energy difference that typically separates the native state of a protein from the unfolded state. The correlation between the free energy and the extent to which the decoy structures differ from the native (as root mean square deviation) is very weak; hence, the free energy is not an accurate measure for ranking the structurally most native-like structures from among a set of models. Analysis of the energy components shows that stability is attained as a result of three major driving forces: (1) minimum size of the protein-water surface interface; (2) minimum total electrostatic energy, which includes solvent polarization; and (3) minimum protein packing energy. The detailed fit required to optimize the last term may underlie difficulties encountered in recovering the native fold from an approximate decoy or model structure.
Collapse
|
15
|
ES/IS: estimation of conformational free energy by combining dynamics simulations with explicit solvent with an implicit solvent continuum model. Biophys Chem 1999; 78:195-205. [PMID: 10343388 DOI: 10.1016/s0301-4622(98)00230-0] [Citation(s) in RCA: 71] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
This paper reviews a recently developed method for calculating the total conformational free energy of a solute macromolecule in water solvent. The method consists of a relatively short simulation by molecular dynamics with explicit solvent molecules (ES) to produce a set of microstates of the macroscopic conformation. Conformational internal solute energy and entropy are obtained from the simulation, the latter in the quasi-harmonic approximation by analysis of the covariance matrix. The implicit solvent (IS) surface energy-dielectric continuum model is used to calculate the average solvation free energy as the sum of the free energies of creating the solute-size hydrophobic cavity, of the van der Waals solute-solvent interactions and of the polarization of water solvent by the solute's charges. We have earlier applied this method to calculate the conformational free energy of native and intentionally misfolded globular conformations of proteins (the EMBL set of deliberately misfolded proteins), and have obtained good discrimination in favor of the native conformations in all instances. These results are summarized and further analyzed to show that, on average, three major component terms of the free energy all contribute in favor of discrimination. We discuss possible improvements of the ES/IS method. It is shown how the force field can be made self-consistent by adapting the parameters for calculation of surface and polarization free energies closely to the molecular mechanics force field used in the dynamics simulation, using established simulation methods to compute free energies for cavity formation and a charging process with the molecular mechanics force field to provide a set of (quasi-experimental) reference data that can be used to refine the parameters of the continuum models. The molecular surface area together with a microscopic surface free energy near 70 cal/(mol A2) is found to be a consistent descriptor of the cavity free energy. Preliminary results indicate that a linear-response approximation for the polarization of water solvent reaction near typical polar and charged protein groups is accurate to within approximately 90%.
Collapse
|
16
|
Abstract
A theoretical study to identify the conformational preferences of lysine-based oligopeptides has been carried out. The solvation free energy and free energy of ionization of the oligopeptides have been calculated by using a fast multigrid boundary element method that considers the coupling between the conformation of the molecule and the ionization equilibria explicitly, at a given pH value. It has been found experimentally that isolated alanine and lysine residues have somewhat small intrinsic helix-forming tendencies; however, results from these simulations indicate that conformations containing right-handed alpha-helical turns are energetically favorable at low values of pH for lysine-based oligopeptides. Also, unusual patterns of interactions among lysine side chains with large hydrophobic contacts and close proximity (5-6 A) between charged NH3+ groups are observed. Similar arrangements of charged groups have been seen for lysine and arginine residues in experimentally determined structures of proteins available from the Protein Data Bank. The lowest-free-energy conformation of the sequence Ac-(LYS)6-NMe from these simulations showed large pKalpha shifts for some of the NH3+ groups of the lysine residues. Such large effects are not observed in the lowest-energy conformations of oligopeptide sequences with two, three, or four lysine residues. Calculations on the sequence Ac-LYS-(ALA)4-LYS-NMe also reveal low-energy alpha-helical conformations with interactions of one of the LYS side chains with the helix backbone in an arrangement quite similar to the one described recently by (Proc. Natl. Acad. Sci. U.S.A. 93:4025-4029). The results of this study provide a sound basis with which to discuss the nature of the interactions, such as hydrophobicity, charge-charge interaction, and solvent polarization effects, that stabilize right-handed alpha-helical conformations.
Collapse
|
17
|
Discrimination between native and intentionally misfolded conformations of proteins: ES/IS, a new method for calculating conformational free energy that uses both dynamics simulations with an explicit solvent and an implicit solvent continuum model. Proteins 1998; 32:399-413. [PMID: 9726412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
A new method for calculating the total conformational free energy of proteins in water solvent is presented. The method consists of a relatively brief simulation by molecular dynamics with explicit solvent (ES) molecules to produce a set of microstates of the macroscopic conformation. Conformational energy and entropy are obtained from the simulation, the latter in the quasi-harmonic approximation by analysis of the covariance matrix. The implicit solvent (IS) dielectric continuum model is used to calculate the average solvation free energy as the sum of the free energies of creating the solute-size hydrophobic cavity, of the van der Waals solute-solvent interactions, and of the polarization of water solvent by the solute's charges. The reliability of the solvation free energy depends on a number of factors: the details of arrangement of the protein's charges, especially those near the surface; the definition of the molecular surface; and the method chosen for solving the Poisson equation. Molecular dynamics simulation in explicit solvent relaxes the protein's conformation and allows polar surface groups to assume conformations compatible with interaction with solvent, while averaging of internal energy and solvation free energy tend to enhance the precision. Two recently developed methods--SIMS, for calculation of a smooth invariant molecular surface, and FAMBE, for solution of the Poisson equation via a fast adaptive multigrid boundary element--have been employed. The SIMS and FAMBE programs scale linearly with the number of atoms. SIMS is superior to Connolly's MS (molecular surface) program: it is faster, more accurate, and more stable, and it smooths singularities of the molecular surface. Solvation free energies calculated with these two programs do not depend on molecular position or orientation and are stable along a molecular dynamics trajectory. We have applied this method to calculate the conformational free energy of native and intentionally misfolded globular conformations of proteins (the EMBL set of deliberately misfolded proteins) and have obtained good discrimination in favor of the native conformations in all instances.
Collapse
|
18
|
Abstract
SIMS, a new method of calculating a smooth invariant molecular dot surface, is presented. The SIMS method generates the smooth molecular surface by rolling two probe spheres. A solvent probe sphere is rolled over the molecule and produces a Richards-Connolly molecular surface (MS), which envelops the solvent-excluded volume of the molecule. In deep crevices, Connolly's method of calculating the MS has two deficiencies. First, it produces self-intersecting parts of the molecular surface, which must be removed to obtain the correct MS. Second, the correct MS is not smooth, i.e., the direction of the normal vector of the MS is not continuous, and some points of the MS are singular. We present an exact method for removing self-intersecting parts and smoothing the singular regions of the MS. The singular MS is smoothed by rolling a smoothing probe sphere over the inward side of the singular MS. The MS in the vicinity of singularities is replaced with the reentrant surface of the smoothing probe sphere. The smoothing method does not disturb the topology of a singular MS, and the smooth MS is a better approximation of the dielectric border between high dielectric solvent and the low dielectric molecular interior. The SIMS method generates a smooth molecular dot surface, which has a quasi-uniform dot distribution in two orthogonal directions on the molecular surface, which is invariant with molecular rotation and stable under changes in the molecular conformation, and which can be used in a variety of implicit methods of modeling solvent effects. The SIMS program is faster than the Connolly MS program, and in a matter of seconds generates a smooth dot MS of a 200-residue protein. The program is available from the authors on request (see http:@femto.med.unc.edu/SIMS).
Collapse
|
19
|
Coupling between folding and ionization equilibria: effects of pH on the conformational preferences of polypeptides. J Mol Biol 1996; 264:770-83. [PMID: 8980685 DOI: 10.1006/jmbi.1996.0676] [Citation(s) in RCA: 73] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
A new approach to the conformational study of polypeptides is presented. It considers explicitly the coupling between the conformation of the molecule and the ionization equilibria at a given pH value. Calculations of the solvation free energy and free energy of ionization of a 17-residue polypeptide are carried out using a fast multigrid boundary element method (MBE). The MBE method uses an adaptive tessellation of the molecular surface by boundary elements with non-regular size to solve the Poisson equation rapidly, and with a high degree of accuracy. The MBE method is integrated into the ECEPP (Empirical Conformational Energy Program for Peptides) algorithm to compute the coupling between the ionization state and the conformation of the molecule. This approach has been applied to study the conformational preference of a short polypeptide for which the available NMR and CD experimental data indicate that conformations containing a right-handed alpha-helical segment are energetically more favorable at low values of pH. The results of calculations using the present method agree quite well with experiments, in contrast to previous applications with standard techniques (using pre-assigned charges at each pH) that were not able to reproduce the experimental findings. Also, it is shown how the coupling to the conformation leads to different degrees of ionization of a given type of residue, for example glutamic acid, at different positions in the amino acid sequence, at any given pH. The results of this study provide a sound basis to discuss the origin of the stability of polypeptide conformations, and its dependence on the environmental conditions.
Collapse
|
20
|
Block-units method for conformational calculations of large nucleic acid chains. II. The two-hierarchical approach and its application to conformational arrangement of the unusual T psi C loop of rabbit tRNA(Val). Biopolymers 1990; 29:1519-29. [PMID: 2386804 DOI: 10.1002/bip.360291203] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The two-level hierarchical methodology is suggested for conformational calculations of large fragments of nucleic acids. The method of the first level is intended for performing a fast screening of the conformational phase space. The high-level method may be used to refine structurally important conformations. The method of the first level is the block-units method, which has been developed specially for these purposes (see part I). It has been shown that the block-units method allows the satisfactory calculation of the structure parameters of the optimal conformations of polynucleotides. The results of the conformational rearrangement calculations of the T psi C loop of the tRNA(Phe) after modification of its sequence are represented.
Collapse
|