Stabilization of the ADP/Metaphosphate Intermediate during ATP Hydrolysis in Pre-power Stroke Myosin: QUANTITATIVE ANATOMY OF AN ENZYME*

It has been proposed recently that ATP hydrolysis in ATPase enzymes proceeds via an initial intermediate in which the dissociated γ-phosphate of ATP is bound in the protein as a metaphosphate (P_γO₃⁻). A combined quantum/classical analysis of this dissociated nucleotide state inside myosin provides a quantitative understanding of how the enzyme stabilizes this unusual metaphosphate. Indeed, in vacuum, the energy of the ADP³⁻·P_γO₃⁻·Mg²⁺ complex is much higher than that of the undissociated ATP⁴⁻. The protein brings it to a surprisingly low value. Energy decomposition reveals how much each interaction in the protein stabilizes the metaphosphate state; backbone peptides of the P-loop contribute 50% of the stabilization energy, and the side chain of Lys-185⁺ contributes 25%. This can be explained by the fact that these groups make strong favorable interactions with the α- and β-phosphates, thus favoring the charge distribution of the metaphosphate state over that of the ATP state. Further stabilization (16%) is achieved by a hydrogen bond between the backbone C=O of Ser-237 (on loop Switch-1) and a water molecule perfectly positioned to attack the P_γO₃⁻ in the subsequent hydrolysis step. The planar and singly negative P_γO₃⁻ is a much better target for the subsequent nucleophilic attack by a negatively charged OH⁻ than the tetrahedral and doubly negative P_γO₄²⁻ group of ATP. Therefore, we argue that the present mechanism of metaphosphate stabilization is common to the large family of nucleotide-hydrolyzing enzymes. Methodologically, this work presents a computational approach that allows us to obtain a truly quantitative conception of enzymatic strategy.     

Impact of Recovery from Desensitization on Acid-sensing Ion Channel-1a (ASIC1a) Current and Response to High Frequency Stimulation

ASIC1a is a neuronal sodium channel activated by external H⁺ ions. To date, all the characterization of ASIC1a has been conducted applying long H⁺ stimuli lasting several seconds. Such experimental protocols weaken and even silence ASIC1a currents to repetitive stimulation. In this work, we examined ASIC1a currents by methods that use rapid application and removal of H⁺. We found that brief H⁺ stimuli, <100 ms, even if applied at high frequency, prevent desensitization thereby generate full and steady peak currents of human ASIC1a. Kinetic analysis of recovery from desensitization of hASIC1a revealed two desensitized states: short- and long-lasting with time constants of τ_Ds ≤0.5 and τ_Dl = 229 s, while in chicken ASIC1a the two desensitized states have similar values τ_D 4.5 s. It is the large difference in stability of the two desensitized states that makes hASIC1a desensitization more pronounced and complex than in cASIC1a. Furthermore, recovery from desensitization was unrelated to cytosolic variations in pH, ATP, PIP₂, or redox state but was dependent on the hydrophobicity of key residues in the first transmembrane segment (TM1). In conclusion, brief H⁺-stimuli maintain steady the magnitude of peak currents thereby the ASIC1a channel is well poised to partake in high frequency signals in the brain.     

A methyl <Superscript>1</Superscript>H double quantum CPMG experiment to study protein conformational exchange

Protein conformational changes play crucial roles in enabling function. The Carr–Purcell–Meiboom–Gill (CPMG) experiment forms the basis for studying such dynamics when they involve the interconversion between highly populated and sparsely formed states, the latter having lifetimes ranging from ~?0.5 to ~?5 ms. Among the suite of experiments that have been developed are those that exploit methyl group probes by recording methyl ¹H single quantum (Tugarinov and Kay in J Am Chem Soc 129:9514–9521, 2007) and triple quantum (Yuwen et al. in Angew Chem Int Ed Engl 55:11490–11494, 2016) relaxation dispersion profiles. Here we build upon these by developing a third experiment in which methyl ¹H double quantum coherences evolve during a CPMG relaxation element. By fitting single, double, and triple quantum datasets, akin to recording the single quantum dataset at static magnetic fields of B_o, 2B_o and 3B_o, we show that accurate exchange values can be obtained even in cases where exchange rates exceed 10,000 s^?1. The utility of the double quantum experiment is demonstrated with a pair of cavity mutants of T4 lysozyme (T4L) with ground and excited states interchanged and with exchange rates differing by fourfold (~?900 s^?1 and ~?3600 s^?1), as well as with a fast-folding domain where the unfolded state lifetime is ~?80 µs.     

Chromophore Structure of Cyanobacterial Phytochrome Cph1 in the Pr State: Reconciling Structural and Spectroscopic Data by QM/MM Calculations

A quantum mechanics (QM)/molecular mechanics (MM) hybrid method was applied to the Pr state of the cyanobacterial phytochrome Cph1 to calculate the Raman spectra of the bound PCB cofactor. Two QM/MM models were derived from the atomic coordinates of the crystal structure. The models differed in the protonation site of His²⁶⁰ in the chromophore-binding pocket such that either the δ-nitrogen (M-HSD) or the ɛ-nitrogen (M-HSE) carried a hydrogen. The optimized structures of the two models display small differences specifically in the orientation of His²⁶⁰ with respect to the PCB cofactor and the hydrogen bond network at the cofactor-binding site. For both models, the calculated Raman spectra of the cofactor reveal a good overall agreement with the experimental resonance Raman (RR) spectra obtained from Cph1 in the crystalline state and in solution, including Cph1 adducts with isotopically labeled PCB. However, a distinctly better reproduction of important details in the experimental spectra is provided by the M-HSD model, which therefore may represent an improved structure of the cofactor site. Thus, QM/MM calculations of chromoproteins may allow for refining crystal structure models in the chromophore-binding pocket guided by the comparison with experimental RR spectra. Analysis of the calculated and experimental spectra also allowed us to identify and assign the modes that sensitively respond to chromophore-protein interactions. The most pronounced effect was noted for the stretching mode of the methine bridge A-B adjacent to the covalent attachment site of PCB. Due a distinct narrowing of the A-B methine bridge bond angle, this mode undergoes a large frequency upshift as compared with the spectrum obtained by QM calculations for the chromophore in vacuo. This protein-induced distortion of the PCB geometry is the main origin of a previous erroneous interpretation of the RR spectra based on QM calculations of the isolated cofactor.Abbreviations: Agp1, phytochrome from Agrobacterium tumefaciens; α-CPC, α-subunit of C-phycocyanin; BV, biliverdin IXα; B3LYP, three-parameter exchange functional according to Becke, Lee, Yang, and Parr; DFT, density functional theory; DrBphP, phytochrome from Deinococcus radiodurans; GAF, domain found in cGMP-specific phosphodiesterases; MM, molecular mechanics; MD, molecular dynamics; N-H ip, N-H in-plane bending; PCB, phycocyanobilin; PED, potential energy distribution; phyA, plant phytochrome; Pr, Pfr, red- and far-red absorbing parent states of phytochrome; PΦB, phytochromobilin; QM, quantum mechanics; RMSD, root mean-square deviation; RR, resonance Raman     

Elastic Torques about Membrane Edges: A Study of Pierced Egg Lecithin Vesicles

The shape of mechanically pierced giant vesicles is studied to obtain the elastic modulus of Gaussian curvature of egg lecithin bilayers. It is argued that such experiments are governed by an apparent modulus, ¯κ_app, not the true modulus of Gaussian curvature, ¯κ. A theory of ¯κ_app is proposed, regarding the pierced bilayer vesicle as a closed monolayer vesicle. The quantity measured, i.e. ¯κ_app/κ, where κ is the rigidity, agrees satisfactorily with the theory. We find ¯κ_app = -(1.9 ± 0.3) · 10^-12 erg (on the basis of κ = (2.3 ± 0.3) · 10^-12 erg). The result may have implications for bilayer fusion.     

Gating Kinetics of Single Large-Conductance Ca2+-Activated K+ Channels in High Ca2+ Suggest a Two-Tiered Allosteric Gating Mechanism?

The Ca²⁺-dependent gating mechanism of large-conductance calcium-activated K⁺ (BK) channels from cultured rat skeletal muscle was examined from low (4 μM) to high (1,024 μM) intracellular concentrations of calcium (Ca²⁺ _i) using single-channel recording. Open probability (P _o) increased with increasing Ca²⁺ _i (K _0.5 11.2 ± 0.3 μM at +30 mV, Hill coefficient of 3.5 ± 0.3), reaching a maximum of ∼0.97 for Ca²⁺ _i ∼ 100 μM. Increasing Ca²⁺ _i further to 1,024 μM had little additional effect on either P _o or the single-channel kinetics. The channels gated among at least three to four open and four to five closed states at high levels of Ca²⁺ _i (>100 μM), compared with three to four open and five to seven closed states at lower Ca²⁺ _i. The ability of kinetic schemes to account for the single-channel kinetics was examined with simultaneous maximum likelihood fitting of two-dimensional (2-D) dwell-time distributions obtained from low to high Ca²⁺ _i. Kinetic schemes drawn from the 10-state Monod-Wyman-Changeux model could not describe the dwell-time distributions from low to high Ca²⁺ _i. Kinetic schemes drawn from Eigen''s general model for a ligand-activated tetrameric protein could approximate the dwell-time distributions but not the dependency (correlations) between adjacent intervals at high Ca²⁺ _i. However, models drawn from a general 50 state two-tiered scheme, in which there were 25 closed states on the upper tier and 25 open states on the lower tier, could approximate both the dwell-time distributions and the dependency from low to high Ca²⁺ _i. In the two-tiered model, the BK channel can open directly from each closed state, and a minimum of five open and five closed states are available for gating at any given Ca²⁺ _i. A model that assumed that the apparent Ca²⁺-binding steps can reach a maximum rate at high Ca²⁺ _i could also approximate the gating from low to high Ca²⁺ _i. The considered models can serve as working hypotheses for the gating of BK channels.     

Magnetic field effects on radical pair intermediates in bacterial photosynthesis

We have investigated the effects of magnetic fields on the formation and decay of excited states in the photochemical reaction centers of Rhodopseudomonas sphaeroides. In chemically reduced reaction centers, a magnetic field decreases the fraction of the transient state P^F that decays by way of the bacteriochlorophyll triplet state P^R. At room temperature, a 2-kG field decreases the quantum yield of P^R by about 40%. In carotenoid-containing reaction centers, the yield of the carotenoid triplet state which forms via P^R is reduced similarly. The effect of the field depends monotonically on field-strength, saturating at about 1 kG. The effect decreases at lower temperatures, when the yield of P^R is higher. Magnetic fields do not significantly affect the formation of the triplet state of bacteriochlorophyll in vitro, the photooxidation of P-870 in reaction centers at moderate redox potential, or the decay kinetics of states P^F and P^R.The effects of magnetic fields support the view that state P^F is a radical pair which is born in a singlet state but undergoes a rapid transformation into a mixture of singlet and triplet states. A simple kinetic model can account for the effects of the field and relate them to the temperature dependence of the yield of P^R.     

Several transition states from C1 to skew conformations of β-d-glucopyranose

The interconversion pathways in the ring distortion of β-d-glucopyranose were investigated using density functional calculations. We examined the energies of several conformers of β-d-glucopyranose and tried to obtain the transition-state conformation and determine the pathway between a ⁴C₁ chair and some distorted ring conformers. The results showed that two E₃/²H₃ conformations and one E₃/⁴H₃ conformation were transition states in such ring puckering. The transition state with the lowest energy conformation is the E₃/²H₃ ring conformation with the side-chain conformation of r-ggG⁺. Intrinsic reaction coordinate calculations indicated that the E₃/²H₃ conformation with the lowest conformational energy is a transition state of the ring interconversion path between the conformations of ⁴C₁ and ²S_O/B_3,O. The energy barrier of this interconversion was 6.13 kcal/mol. As far as we know, this is the first example of finding pathways for an interconversion of glucopyranose ring puckering at the level of a quantum chemical calculation.     

Structural Interpretation of Paramagnetic Relaxation Enhancement-Derived Distances for Disordered Protein States

Paramagnetic relaxation enhancement (PRE) is a powerful technique for studying transient tertiary organizations of unfolded and partially folded proteins. The heterogeneous and dynamic nature of disordered protein states, together with the r⁻⁶ dependence of PRE, presents significant challenges for reliable structural interpretation of PRE-derived distances. Without additional knowledge of accessible conformational substates, ensemble-simulation-based protocols have been used to calculate structure ensembles that appear to be consistent with the PRE distance restraints imposed on the ensemble level with the proper r⁻⁶ weighting. However, rigorous assessment of the reliability of such protocols has been difficult without intimate knowledge of the true nature of disordered protein states. Here we utilize sets of theoretical PRE distances derived from simulated structure ensembles that represent the folded, partially folded and unfolded states of a small protein to investigate the efficacy of ensemble-simulation-based structural interpretation of PRE distances. The results confirm a critical limitation that, due to r⁻⁶ weighting, only one or a few members need to satisfy the distance restraints and the rest of the ensemble are essentially unrestrained. Consequently, calculated structure ensembles will appear artificially heterogeneous no matter whether the PRE distances are derived from the folded, partially unfolded or unfolded state. Furthermore, the nature of the heterogeneous ensembles is largely determined by the protein model employed in structure calculation and reflects little on the true nature of the underlying disordered state. These findings suggest that PRE measurements on disordered protein states alone generally do not contain enough information for a reliable structural interpretation and that the latter will require additional knowledge of accessible conformational substates. Interestingly, when a very large number of PRE measurements is available, faithful structural interpretation might be possible with intermediate ensemble sizes under ideal conditions.     

Steric Crowding of the Turn Region Alters the Tertiary Fold of Amyloid-β18–35 and Makes It Soluble

Aβ self-assembles into parallel cross-β fibrillar aggregates, which is associated with Alzheimer''s disease pathology. A central hairpin turn around residues 23–29 is a defining characteristic of Aβ in its aggregated state. Major biophysical properties of Aβ, including this turn, remain unaltered in the central fragment Aβ_18–35. Here, we synthesize a single deletion mutant, ΔG25, with the aim of sterically hindering the hairpin turn in Aβ_18–35. We find that the solubility of the peptide goes up by more than 20-fold. Although some oligomeric structures do form, solution state NMR spectroscopy shows that they have mostly random coil conformations. Fibrils ultimately form at a much higher concentration but have widths approximately twice that of Aβ_18–35, suggesting an opening of the hairpin bend. Surprisingly, two-dimensional solid state NMR shows that the contact between Phe¹⁹ and Leu³⁴ residues, observed in full-length Aβ and Aβ_18–35, is still intact in these fibrils. This is possible if the monomers in the fibril are arranged in an antiparallel β-sheet conformation. Indeed, IR measurements, supported by tyrosine cross-linking experiments, provide a characteristic signature of the antiparallel β-sheet. We conclude that the self-assembly of Aβ is critically dependent on the hairpin turn and on the contact between the Phe¹⁹ and Leu³⁴ regions, making them potentially sensitive targets for Alzheimer''s therapeutics. Our results show the importance of specific conformations in an aggregation process thought to be primarily driven by nonspecific hydrophobic interactions.     

Why the Lysogenic State of Phage λ Is So Stable: A Mathematical Modeling Approach

We develop a mathematical model of the phage λ lysis/lysogeny switch, taking into account recent experimental evidence demonstrating enhanced cooperativity between the left and right operator regions. Model parameters are estimated from available experimental data. The model is shown to have a single stable steady state for these estimated parameter values, and this steady state corresponds to the lysogenic state. When the CI degradation rate (γ_cI) is slightly increased from its normal value (γ_cI 0.0 min⁻¹), two additional steady states appear (through a saddle-node bifurcation) in addition to the lysogenic state. One of these new steady states is stable and corresponds to the lytic state. The other steady state is an (unstable) saddle node. The coexistence these two globally stable steady states (the lytic and lysogenic states) is maintained with further increases of γ_cI until γ_cI 0.35 min⁻¹, when the lysogenic steady state and the saddle node collide and vanish (through a reverse saddle node bifurcation) leaving only the lytic state surviving. These results allow us to understand the high degree of stability of the lysogenic state because, normally, it is the only steady state. Further implications of these results for the stability of the phage λ switch are discussed, as well as possible experimental tests of the model.     

NMR characterization of copper and lipid interactions of the C2B domain of synaptotagmin I—relevance to the non-classical secretion of the human acidic fibroblast growth factor (hFGF-1)

Human fibroblast growth factor (hFGF-1) is a ∼ 17 kDa heparin binding cytokine. It lacks the conventional hydrophobic N-terminal signal sequence and is secreted through non-classical secretion routes. Under stress, hFGF-1 is released as a multiprotein complex consisting of hFGF-1, S100A13 (a calcium binding protein), and p40 synaptotagmin (Syt1). Copper (Cu²⁺) is shown to be required for the formation of the multiprotein hFGF-1 release complex (Landriscina et al. ,2001; Di Serio et al., 2008). Syt1, containing the lipid binding C2B domain, is believed to play an important role in the eventual export of the hFGF-1 across the lipid bilayer. In this study, we characterize Cu²⁺ and lipid interactions of the C2B domain of Syt1 using multidimensional NMR spectroscopy. The results highlight how Cu²⁺ appears to stabilize the protein bound to pS vesicles. Cu²⁺ and lipid binding interface mapped using 2D ¹H-¹⁵N heteronuclear single quantum coherence experiments reveal that residues in β-strand I contributes to the unique Cu²⁺ binding site in the C2B domain. In the absence of metal ions, residues located in Loop II and β-strand IV contribute to binding to unilamelar pS vesicles. In the presence of Cu²⁺, additional residues located in Loops I and III appear to stabilize the protein-lipid interactions. The results of this study provide valuable information towards understanding the molecular mechanism of the Cu²⁺-induced non-classical secretion of hFGF-1.     

Contrasting Effects of Two Lipid Cofactors of Prion Replication on the Conformation of the Prion Protein

Recent studies introduced two experimental protocols for converting full-length recombinant prion protein (rPrP) purified from E.coli into the infectious prion state (PrP^Sc) with high infectivity titers. Both protocols employed protein misfolding cyclic amplification (PMCA) for generating PrP^Sc de novo, but used two different lipids, 1-palmitoyl-2-oleolyl-sn-glycero-3-phospho(1’-rac-glycerol) (POPG) or phosphatidylethanolamine (PE), as conversion cofactors. The current study compares the effect of POPG and PE on the physical properties of native, α-helical full-length mouse rPrP under the solvent conditions used for converting rPrP into PrP^Sc. Surprisingly, the effects of POPG and PE on rPrP physical properties, including its conformation, thermodynamic stability, aggregation state and interaction with a lipid, were found to be remarkably different. PE was shown to have minimal, if any, effects on rPrP thermodynamic stability, cooperativity of unfolding, immediate solvent environment or aggregation state. In fact, little evidence indicates that PE interacts with rPrP directly. In contrast, POPG was found to bind to and induce dramatic changes in rPrP structure, including a loss of α-helical conformation and formation of large lipid-protein aggregates that were resistant to partially denaturing conditions. These results suggest that the mechanisms by which lipids assist conversion of rPrP into PrP^Sc might be fundamentally different for POPG and PE.     

High Resolution 3D Imaging of Ex-Vivo Biological Samples by Micro CT

Non-destructive volume visualization can be achieved only by tomographic techniques, of which the most efficient is the x-ray micro computerized tomography (μCT).High resolution μCT is a very versatile yet accurate (1-2 microns of resolution) technique for 3D examination of ex-vivo biological samples^{1, 2}. As opposed to electron tomography, the μCT allows the examination of up to 4 cm thick samples. This technique requires only few hours of measurement as compared to weeks in histology. In addition, μCT does not rely on 2D stereologic models, thus it may complement and in some cases can even replace histological methods^{3, 4}, which are both time consuming and destructive. Sample conditioning and positioning in μCT is straightforward and does not require high vacuum or low temperatures, which may adversely affect the structure. The sample is positioned and rotated 180° or 360°between a microfocused x-ray source and a detector, which includes a scintillator and an accurate CCD camera, For each angle a 2D image is taken, and then the entire volume is reconstructed using one of the different available algorithms^5-7. The 3D resolution increases with the decrease of the rotation step. The present video protocol shows the main steps in preparation, immobilization and positioning of the sample followed by imaging at high resolution.     

Logical identification of all steady states: The concept of feedback loop characteristic states

Biological regulatory systems can be described in terms of non-linear differential equations or in logical terms (using an “infinitely non-linear” approximation). Until recently, only part of the steady states of a system could be identified on logical grounds. The reason was that steady states frequently have one or more variable located on a threshold (see below); those steady states were not detected because so far no logical status was assigned to threshold values. This is why we introduced logical scales with values 0,¹θ, 1²θ, 2, ..., in which¹θ,²θ, ... are the logical values assigned to the successive thresholds of the scale. We thus have, in addition to the regular logical states,singular states in which one or more variables is located on a threshold. This permits identifyingall the steady states on logical grounds. It was noticed that each feedback loop (or reunion of disjointed loops) can be characterized by a logical state located at the thresholds at which the variables of the loop operate. This led to the concept ofloop-characteristic state, which, as we will see, enormously simplifies the analysis.The core of this paper is a formal demonstration that among the singular states of a system, only loop-characteristic states can be steady. Reciprocally, given a loop-characteristic state, there are parameter values for which this state is steady; in this case, the loop is effective (i.e. it generates multistationarity if it is a positive loop, homeostasis if it is a negative loop). This not only results in the above-mentioned radical simplification of the identification of the steady states, but in an entirely new view of the relation between feedback loops and steady states.     

Non-Gaussian noise-optimized intracellular cytosolic calcium oscillations

We have numerically studied the effect of a particular kind of non-Gaussian colored noise (NGN), characterized by the deviation q from Gaussian noise (q = 1), on intracellular cytosolic calcium (Ca²⁺) oscillations. It is found that, as q is increased, the Ca²⁺ oscillation regularity increases and reaches a best performance at an optimal q, and then decreases with further increasing q, which represents the occurrence of coherence resonance, i.e., the most regular Ca²⁺ oscillations. Similar phenomena occur for different values of noise intensity and correlation time of the NGN. This phenomenon of deviation-optimized Ca²⁺ oscillations show that, external non-Gaussian noises of different types can enhance and even optimize the intracellular Ca²⁺ oscillations. This result provides new insights into the constructive roles and potential applications of non-Gaussian noises in intracellular cytosolic Ca²⁺ oscillations.     

Structural Basis for Conformational Dynamics of GTP-bound Ras Protein

Ras family small GTPases assume two interconverting conformations, “inactive” state 1 and “active” state 2, in their GTP-bound forms. Here, to clarify the mechanism of state transition, we have carried out x-ray crystal structure analyses of a series of mutant H-Ras and M-Ras in complex with guanosine 5′-(β,γ-imido)triphosphate (GppNHp), representing various intermediate states of the transition. Crystallization of H-RasT35S-GppNHp enables us to solve the first complete tertiary structure of H-Ras state 1 possessing two surface pockets unseen in the state 2 or H-Ras-GDP structure. Moreover, determination of the two distinct crystal structures of H-RasT35S-GppNHp, showing prominent polysterism in the switch I and switch II regions, reveals a pivotal role of the guanine nucleotide-mediated interaction between the two switch regions and its rearrangement by a nucleotide positional change in the state 2 to state 1 transition. Furthermore, the ³¹P NMR spectra and crystal structures of the GppNHp-bound forms of M-Ras mutants, carrying various H-Ras-type amino acid substitutions, also reveal the existence of a surface pocket in state 1 and support a similar mechanism based on the nucleotide-mediated interaction and its rearrangement in the state 1 to state 2 transition. Intriguingly, the conformational changes accompanying the state transition mimic those that occurred upon GDP/GTP exchange, indicating a common mechanistic basis inherent in the high flexibility of the switch regions. Collectively, these results clarify the structural features distinguishing the two states and provide new insights into the molecular basis for the state transition of Ras protein.     

The spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β

To provide molecular-level insights into the spontaneous replication error and the mismatch discrimination mechanisms of human DNA polymerase β (polβ), we report four crystal structures of polβ complexed with dG•dTTP and dA•dCTP mismatches in the presence of Mg²⁺ or Mn²⁺. The Mg²⁺-bound ground-state structures show that the dA•dCTP-Mg²⁺ complex adopts an ‘intermediate’ protein conformation while the dG•dTTP-Mg²⁺ complex adopts an open protein conformation. The Mn²⁺-bound ‘pre-chemistry-state’ structures show that the dA•dCTP-Mn²⁺ complex is structurally very similar to the dA•dCTP-Mg²⁺ complex, whereas the dG•dTTP-Mn²⁺ complex undergoes a large-scale conformational change to adopt a Watson–Crick-like dG•dTTP base pair and a closed protein conformation. These structural differences, together with our molecular dynamics simulation studies, suggest that polβ increases replication fidelity via a two-stage mismatch discrimination mechanism, where one is in the ground state and the other in the closed conformation state. In the closed conformation state, polβ appears to allow only a Watson–Crick-like conformation for purine•pyrimidine base pairs, thereby discriminating the mismatched base pairs based on their ability to form the Watson–Crick-like conformation. Overall, the present studies provide new insights into the spontaneous replication error and the replication fidelity mechanisms of polβ.