Electrostatic free energy calculations using the generalized solvent boundary potential method

Copyright 2002 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in The Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1507108.Free energyperturbation (FEP) calculations using all-atom molecular dynamics simulations with a large number of explicit solvent molecules are a powerful approach to study ligand–macromolecule association processes at the atomic level. One strategy to carry out FEP calculations efficiently and reduce computational time is to consider the explicit dynamics of only a small number of atoms in a localized region around the ligand. Such an approximation is motivated by the observation that the factors governing binding specificity are dominated by interactions in the vicinity of the ligand. However, a straightforward truncation of the system may yield inaccurate results as the influence exerted by the remote regions of the macromolecule and the surrounding solvent through long-range electrostatic effects may be significant. To obtain meaningful results, it is important to incorporate the influence of the remote regions of the ligand–macromolecule complex implicitly using some effective potential. The generalized solvent boundary potential (GSBP) that was developed recently [W. Im, S. Bernèche, and B. Roux, J. Chem. Phys. 114, 2924 (2001)] is an efficient computational method to represent the long-range electrostaticinteractions arising from remote (outer) regions in simulations of a localized (inner) region with a small number of explicit atoms. In the present work, FEP calculations combined with GSBP are used to illustrate the importance of these long-range electrostatic factors in estimation of the charging free energy of an aspartate ligand bound to the aspartyl-tRNA synthetase. Calculations with explicit spherical simulation inner regions of different radii are used to test the accuracy of the GSBP method and also illustrate the importance of explicit protein and solvent dynamics in the free energy estimation. The influence of the represented outer region is tested using separate simulations in which the reaction field and/or the protein static field are excluded. Both components are shown to be essential to obtain quantitatively meaningful results. The ability of implicitly treating the influence of protein fluctuations in the outer region using a protein dielectric constant is examined. It is shown that accurate charging free energy calculations can be performed for this system with a spherical region of 15 to 20 Å radius, which roughly corresponds to 1500–3500 moving atoms. The results indicate that GSBP in combination with FEP calculations is a precise and efficient approach to include long-range electrostatic effects in the study of ligand binding to large macromolecules

Starvation sensing by mycobacterial RelA/SpoT homologue through constitutive surveillance of translation

The stringent response, which leads to persistence of nutrient-starved mycobacteria, is induced by activation of the RelA/SpoT homolog (Rsh) upon entry of a deacylated-tRNA in a translating ribosome. However, the mechanism by which Rsh identifies such ribosomes in vivo remains unclear. Here, we show that conditions inducing ribosome hibernation result in loss of intracellular Rsh in a Clp protease-dependent manner. This loss is also observed in nonstarved cells using mutations in Rsh that block its interaction with the ribosome, indicating that Rsh association with the ribosome is important for Rsh stability. The cryo-EM structure of the Rsh-bound 70S ribosome in a translation initiation complex reveals unknown interactions between the ACT domain of Rsh and components of the ribosomal L7/L12 stalk base, suggesting that the aminoacylation status of A-site tRNA is surveilled during the first cycle of elongation. Altogether, we propose a surveillance model of Rsh activation that originates from its constitutive interaction with the ribosomes entering the translation cycle

Starvation Sensing by Mycobacterial RelA/SpoT Homologue Through Constitutive Surveillance of Translation

The stringent response, which leads to persistence of nutrient-starved mycobacteria, is induced by activation of the RelA/SpoT homolog (Rsh) upon entry of a deacylated-tRNA in a translating ribosome. However, the mechanism by which Rsh identifies such ribosomes in vivo remains unclear. Here, we show that conditions inducing ribosome hibernation result in loss of intracellular Rsh in a Clp protease-dependent manner. This loss is also observed in nonstarved cells using mutations in Rsh that block its interaction with the ribosome, indicating that Rsh association with the ribosome is important for Rsh stability. The cryo-EM structure of the Rsh-bound 70S ribosome in a translation initiation complex reveals unknown interactions between the ACT domain of Rsh and components of the ribosomal L7/L12 stalk base, suggesting that the aminoacylation status of A-site tRNA is surveilled during the first cycle of elongation. Altogether, we propose a surveillance model of Rsh activation that originates from its constitutive interaction with the ribosomes entering the translation cycle

Characterizing Structural Transitions Using Localized Free Energy Landscape Analysis

Structural changes in molecules are frequently observed during biological processes like replication, transcription and translation. These structural changes can usually be traced to specific distortions in the backbones of the macromolecules involved. Quantitative energetic characterization of such distortions can greatly advance the atomic-level understanding of the dynamic character of these biological processes.Molecular dynamics simulations combined with a variation of the Weighted Histogram Analysis Method for potential of mean force determination are applied to characterize localized structural changes for the test case of cytosine (underlined) base flipping in a GTCAGCGCATGG DNA duplex. Free energy landscapes for backbone torsion and sugar pucker degrees of freedom in the DNA are used to understand their behavior in response to the base flipping perturbation. By simplifying the base flipping structural change into a two-state model, a free energy difference of upto 14 kcal/mol can be attributed to the flipped state relative to the stacked Watson-Crick base paired state. This two-state classification allows precise evaluation of the effect of base flipping on local backbone degrees of freedom.The calculated free energy landscapes of individual backbone and sugar degrees of freedom expectedly show the greatest change in the vicinity of the flipping base itself, but specific delocalized effects can be discerned upto four nucleotide positions away in both 5' and 3' directions. Free energy landscape analysis thus provides a quantitative method to pinpoint the determinants of structural change on the atomic scale and also delineate the extent of propagation of the perturbation along the molecule. In addition to nucleic acids, this methodology is anticipated to be useful for studying conformational changes in all macromolecules, including carbohydrates, lipids, and proteins

Partial Base Flipping Is Sufficient for Strand Slippage near DNA Duplex Termini

Strand slippage is a structural mechanism by which insertion–deletion (indel) mutations are introduced during replication by polymerases. Three-dimensional atomic-resolution structural pathways are still not known for the decades-old template slippage description. The dynamic nature of the process and the higher energy intermediates involved increase the difficulty of studying these processes experimentally. In the present study, restrained and unrestrained molecular dynamics simulations, carried out using multiple nucleic acid force fields, are used to demonstrate that partial base-flipping can be sufficient for strand slippage at DNA duplex termini. Such strand slippage can occur in either strand, i.e. near either the 3′ or the 5′ terminus of a DNA strand, which suggests that similar structural flipping mechanisms can cause both primer and template slippage. In the repetitive mutation hot-spot sequence studied, non-canonical base-pairing with exposed DNA groove atoms of a neighboring G:C base-pair stabilizes a partially flipped state of the cytosine. For its base-pair partner guanine, a similar partially flipped metastable intermediate was not detected, and the propensity for sustained slippage was also found to be lower. This illustrates that a relatively small metastable DNA structural distortion in polymerase active sites could allow single base insertion or deletion mutations to occur, and stringent DNA groove molecular recognition may be required to maintain intrinsic DNA polymerase fidelity. The implications of a close relationship between base-pair dissociation, base unstacking, and strand slippage are discussed in the context of sequence dependence of indel mutations

Partial Base Flipping Is Sufficient for Strand Slippage near DNA Duplex Termini

Strand slippage is a structural mechanism by which insertion–deletion (indel) mutations are introduced during replication by polymerases. Three-dimensional atomic-resolution structural pathways are still not known for the decades-old template slippage description. The dynamic nature of the process and the higher energy intermediates involved increase the difficulty of studying these processes experimentally. In the present study, restrained and unrestrained molecular dynamics simulations, carried out using multiple nucleic acid force fields, are used to demonstrate that partial base-flipping can be sufficient for strand slippage at DNA duplex termini. Such strand slippage can occur in either strand, i.e. near either the 3′ or the 5′ terminus of a DNA strand, which suggests that similar structural flipping mechanisms can cause both primer and template slippage. In the repetitive mutation hot-spot sequence studied, non-canonical base-pairing with exposed DNA groove atoms of a neighboring G:C base-pair stabilizes a partially flipped state of the cytosine. For its base-pair partner guanine, a similar partially flipped metastable intermediate was not detected, and the propensity for sustained slippage was also found to be lower. This illustrates that a relatively small metastable DNA structural distortion in polymerase active sites could allow single base insertion or deletion mutations to occur, and stringent DNA groove molecular recognition may be required to maintain intrinsic DNA polymerase fidelity. The implications of a close relationship between base-pair dissociation, base unstacking, and strand slippage are discussed in the context of sequence dependence of indel mutations

Analyzing the Relationship between Single Base Flipping and Strand Slippage near DNA Duplex Termini

Insertion–deletion (indel) mutations are caused by strand slippage between pairing primer and template strands during nucleic acid strand extension. A possible causative factor for such strand slippage is base flipping in the primer strand or template strand, for insertion or deletion mutations, respectively. A simple mechanistic description is that the “hole” in the nucleic acid duplex left behind by a flipping base is occupied by a neighboring base on the same strand, resulting in slippage with respect to its paired strand. The extent of single base flipping required for occupation of its former place in the double helix by a neighboring base is not fully understood. The present study uses restrained molecular dynamics (MD) simulations along a pseudohedihedral base flipping parameter to construct two-dimensional free energy profiles along base flipping and strand slippage geometric parameters. These profiles, generated for both cytosine and guanine single base flipping in a short repetitive indel mutation hot-spot DNA sequence, illustrate the extent of single base flipping that can allow strand slippage by one base position. Relatively minor base flipping into both the major and minor grooves can result in strand slippage. Deconstruction of the collective variable strand slippage geometric parameter into its component distances illustrates the details of how strand slippage can accompany base flipping. The trans Watson–Crick:sugar edge interaction that stabilizes cytosine flipping in this hot-spot sequence is also characterized energetically. The impact of these results on understanding sequence dependence of indel errors in nucleic acid strand extension is discussed, along with a suggestion for future studies that can generalize the present findings to all nearest-neighbor sequence contexts

Partial Base Flipping Is Sufficient for Strand Slippage near DNA Duplex Termini

Strand slippage is a structural mechanism by which insertion–deletion (indel) mutations are introduced during replication by polymerases. Three-dimensional atomic-resolution structural pathways are still not known for the decades-old template slippage description. The dynamic nature of the process and the higher energy intermediates involved increase the difficulty of studying these processes experimentally. In the present study, restrained and unrestrained molecular dynamics simulations, carried out using multiple nucleic acid force fields, are used to demonstrate that partial base-flipping can be sufficient for strand slippage at DNA duplex termini. Such strand slippage can occur in either strand, i.e. near either the 3′ or the 5′ terminus of a DNA strand, which suggests that similar structural flipping mechanisms can cause both primer and template slippage. In the repetitive mutation hot-spot sequence studied, non-canonical base-pairing with exposed DNA groove atoms of a neighboring G:C base-pair stabilizes a partially flipped state of the cytosine. For its base-pair partner guanine, a similar partially flipped metastable intermediate was not detected, and the propensity for sustained slippage was also found to be lower. This illustrates that a relatively small metastable DNA structural distortion in polymerase active sites could allow single base insertion or deletion mutations to occur, and stringent DNA groove molecular recognition may be required to maintain intrinsic DNA polymerase fidelity. The implications of a close relationship between base-pair dissociation, base unstacking, and strand slippage are discussed in the context of sequence dependence of indel mutations

Partial Base Flipping Is Sufficient for Strand Slippage near DNA Duplex Termini

Strand slippage is a structural mechanism by which insertion–deletion (indel) mutations are introduced during replication by polymerases. Three-dimensional atomic-resolution structural pathways are still not known for the decades-old template slippage description. The dynamic nature of the process and the higher energy intermediates involved increase the difficulty of studying these processes experimentally. In the present study, restrained and unrestrained molecular dynamics simulations, carried out using multiple nucleic acid force fields, are used to demonstrate that partial base-flipping can be sufficient for strand slippage at DNA duplex termini. Such strand slippage can occur in either strand, i.e. near either the 3′ or the 5′ terminus of a DNA strand, which suggests that similar structural flipping mechanisms can cause both primer and template slippage. In the repetitive mutation hot-spot sequence studied, non-canonical base-pairing with exposed DNA groove atoms of a neighboring G:C base-pair stabilizes a partially flipped state of the cytosine. For its base-pair partner guanine, a similar partially flipped metastable intermediate was not detected, and the propensity for sustained slippage was also found to be lower. This illustrates that a relatively small metastable DNA structural distortion in polymerase active sites could allow single base insertion or deletion mutations to occur, and stringent DNA groove molecular recognition may be required to maintain intrinsic DNA polymerase fidelity. The implications of a close relationship between base-pair dissociation, base unstacking, and strand slippage are discussed in the context of sequence dependence of indel mutations

Analyzing the Relationship between Single Base Flipping and Strand Slippage near DNA Duplex Termini

Insertion–deletion (indel) mutations are caused by strand slippage between pairing primer and template strands during nucleic acid strand extension. A possible causative factor for such strand slippage is base flipping in the primer strand or template strand, for insertion or deletion mutations, respectively. A simple mechanistic description is that the “hole” in the nucleic acid duplex left behind by a flipping base is occupied by a neighboring base on the same strand, resulting in slippage with respect to its paired strand. The extent of single base flipping required for occupation of its former place in the double helix by a neighboring base is not fully understood. The present study uses restrained molecular dynamics (MD) simulations along a pseudohedihedral base flipping parameter to construct two-dimensional free energy profiles along base flipping and strand slippage geometric parameters. These profiles, generated for both cytosine and guanine single base flipping in a short repetitive indel mutation hot-spot DNA sequence, illustrate the extent of single base flipping that can allow strand slippage by one base position. Relatively minor base flipping into both the major and minor grooves can result in strand slippage. Deconstruction of the collective variable strand slippage geometric parameter into its component distances illustrates the details of how strand slippage can accompany base flipping. The trans Watson–Crick:sugar edge interaction that stabilizes cytosine flipping in this hot-spot sequence is also characterized energetically. The impact of these results on understanding sequence dependence of indel errors in nucleic acid strand extension is discussed, along with a suggestion for future studies that can generalize the present findings to all nearest-neighbor sequence contexts