Dynamics of spontaneous flipping of a mismatched base in DNA duplex

The Thermodynamic and Kinetic Properties of the dA-rU DNA-RNA Hybrid Base Pair Investigated via Molecular Dynamics Simulations

DNA-RNA hybrid duplexes play essential roles during the reverse transcription of RNA viruses and DNA replication. The opening and conformation changes of individual base pairs are critical to their biological functions. However, the microscopic mechanisms governing base pair closing and opening at the atomic level remain poorly understood. In this study, we investigated the thermodynamic and kinetic parameters of the dA-rU base pair in a DNA-RNA hybrid duplex using 4 μs all-atom molecular dynamics (MD) simulations at different temperatures. Our results showed that the thermodynamic parameters of the dA-rU base pair aligned with the predictions of the nearest-neighbor model and were close to those of the AU base pair in RNA. The temperature dependence of the average lifetimes of both the open and the closed states, as well as the transition path times, were obtained. The free-energy barrier for a single base pair opening and closing arises from an increase in enthalpy due to the disruption of the base-stacking interactions and hydrogen bonding, along with an entropy loss attributed to the accompanying restrictions, such as torsional angle constraints and solvent viscosity.

Effect of molecular crowders on ligand binding kinetics with G-quadruplex DNA probed by fluorescence correlation spectroscopy

Guanine-rich single-stranded DNA folds into G-quadruplex DNA (GqDNA) structures, which play crucial roles in various biological processes. These structures are also promising targets for ligands, potentially inducing antitumor effects. While thermodynamic parameters of ligand/DNA interactions are well-studied, the kinetics of ligand interaction with GqDNA, particularly in cell-like crowded environments, remain less explored. In this study, we investigate the impact of molecular crowding agents (glucose, sucrose, and ficoll 70) at physiologically relevant concentrations (20% w/v) on the association and dissociation rates of the benzophenoxazine-core based ligand, cresyl violet (CV), with human telomeric antiparallel-GqDNA. We utilized fluorescence correlation spectroscopy (FCS) along with other techniques. Our findings reveal that crowding agents decrease the binding affinity of CV to GqDNA, with the most significant effect-a nearly three-fold decrease-observed with ficoll 70. FCS measurements indicate that this decrease is primarily due to a viscosity-induced slowdown of ligand association in the crowded environment. Interestingly, dissociation rates remain largely unaffected by smaller crowders, with only small effect observed in presence of ficoll 70 due to direct but weak interaction between the ligand and ficoll. These results along with previously reported data provide valuable insights into ligand/GqDNA interactions in cellular contexts, suggesting a conserved mechanism of saccharide crowder influence, regardless of variations in GqDNA structure and ligand binding mode. This underscores the importance of considering crowding effects in the design and development of GqDNA-targeted drugs for potential cancer treatment.

The role of nucleotide opening dynamics in facilitated target search by DNA-repair proteins

Preserving the genomic integrity stands a fundamental necessity, primarily achieved by the DNA repair proteins through their continuous patrolling on the DNA in search of lesions. However, comprehending how even a single base-pair lesion can be swiftly and specifically recognized amidst millions of base-pair sites remains a formidable challenge. In this study, we employ extensive molecular dynamics simulations using an appropriately tuned model of both protein and DNA to probe the underlying molecular principles. Our findings reveal that the dynamics of a non-canonical base generate an entropic signal that guides the one-dimensional search of a repair protein, thereby facilitating the recognition of the lesion site. The width of the funnel perfectly aligns with the one-dimensional diffusion length of DNA-binding proteins. The generic mechanism provides a physical basis for rapid recognition and specificity of DNA damage sensing and recognition.

Backbone Conformational Equilibrium in Mismatched DNA Correlates with Enzyme Activity

T:G mismatches in mammals arise primarily from the deamination of methylated CpG sites or the incorporation of improper nucleotides. The process by which repair enzymes such as thymine DNA glycosylase (TDG) identify a canonical DNA base in the incorrect pairing context remains a mystery. However, the abundant contacts of the repair enzymes with the DNA backbone suggest a role for protein-phosphate interaction in the recognition and repair processes, where conformational properties may facilitate the proper interactions. We have previously used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion and the effect of a mismatch or lesion compared to canonical DNA and found stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA. We have currently compared our results to substrate dependence for TDG, MBD4, M. HhaI, and CEBPβ, testing for correlations to sequence and base-pair dependence. We found strong correlations of our DNA phosphate backbone equilibrium (Keq) to different enzyme kinetics or binding parameters of these varied enzymes, suggesting that the backbone equilibrium may play an important role in mismatch recognition and/or conformational rearrangement and energetics during nucleotide flipping or other aspects of enzyme interrogation of the DNA substrate.

Register-Shifted Structures in Base-Flipped Uracil-Damaged DNA

We report the occurrence of register-shifted structures in simulations of uracil-containing dsDNA. These occur when the 3' base vicinal to uracil is thymine in U:A base-paired DNA. Upon base flipping of uracil, this 3' thymine hydrogen bonds with the adenine across the uracil instead of its complementary base. The register-shifted structure is persistent and sterically blocks re-entry of uracil into the helix stack. Register shifting might be important for DNA repair since the longer exposure of the lesion in register-shifted structures could facilitate enzymatic recognition and repair.

Dynamics of 8-Oxoguanine in DNA: Decisive Effects of Base Pairing and Nucleotide Context

8-Oxo-7,8-dihydroguanine (oxoG), an abundant DNA lesion, can mispair with adenine and induce mutations. To prevent this, cells possess DNA repair glycosylases that excise either oxoG from oxoG:C pairs (bacterial Fpg, human OGG1) or A from oxoG:A mispairs (bacterial MutY, human MUTYH). Early lesion recognition steps remain murky and may include enforced base pair opening or capture of a spontaneously opened pair. We adapted the CLEANEX-PM NMR protocol to detect DNA imino proton exchange and analyzed the dynamics of oxoG:C, oxoG:A, and their undamaged counterparts in nucleotide contexts with different stacking energy. Even in a poorly stacking context, the oxoG:C pair did not open easier than G:C, arguing against extrahelical base capture by Fpg/OGG1. On the contrary, oxoG opposite A significantly populated the extrahelical state, which may assist recognition by MutY/MUTYH.

Thermodynamics and kinetics of an A-U RNA base pair under force studied by molecular dynamics simulations

Mechanical force has been widely used to study RNA folding and unfolding. Understanding how the force affects the opening and closing of a single base pair, which is a basic step for RNA folding and unfolding and a fundamental behavior in some important biological activities, is crucial to understanding the mechanism of RNA folding and unfolding under mechanical force. In this work, we investigated the opening and closing process of an RNA base pair under mechanical force with constant-force stretching molecular dynamics simulations. It was found that high mechanical force results in overstretching, and the open state is a high-energy state. The enthalpy and entropy change of the base-pair opening-closing transition were obtained and the results at low forces were in good agreement with the nearest-neighbor model. The temperature and force dependence of the opening and closing rates were also obtained. The position of the transition state for the base-pair opening-closing transition under mechanical force was determined. The free energy barrier of opening a base pair without force is the enthalpy increase, and the work done by the force from the closed state to the transition state decreases the barrier and increases the opening rate. The free energy barrier of closing the base pair without force results from the entropy loss, and the work done by the force from the open state to the transition state increases the barrier and decreases the closing rate. The transition rates are strongly dependent on the temperature and force, while the transition path times are weakly dependent on force and temperature.

High-precision mapping reveals rare N6-deoxyadenosine methylation in the mammalian genome

N⁶-deoxyadenosine methylation (6mA) is the most widespread type of DNA modification in prokaryotes and is also abundantly distributed in some unicellular eukaryotes. However, 6mA levels are remarkably low in mammals. The lack of a precise and comprehensive mapping method has hindered more advanced investigations of 6mA. Here, we report a new method MM-seq (modification-induced mismatch sequencing) for genome-wide 6mA mapping based on a novel detection principle. We found that modified DNA bases are prone to form a local open region that allows capture by antibody, for example, via a DNA breathing or base-flipping mechanism. Specified endonuclease or exonuclease can recognize the antibody-stabilized mismatch-like structure and mark the exact modified sites for sequencing readout. Using this method, we examined the genomic positions of 6mA in bacteria (E. coli), green algae (C. reinhardtii), and mammalian cells (HEK239T, Huh7, and HeLa cells). In contrast to bacteria and green algae, human cells possess a very limited number of 6mA sites which are sporadically distributed across the genome of different cell types. After knocking out the RNA m⁶A methyltransferase METTL3 in mouse ES cells, 6mA becomes mostly diminished. Our results imply that rare 6mA in the mammalian genome is introduced by RNA m⁶A machinery via a non-targeted mechanism.

"Flexible hinge" dynamics in mismatched DNA revealed by fluorescence correlation spectroscopy

Altered unwinding/bending fluctuations at DNA lesion sites are implicated as plausible mechanisms for damage sensing by DNA-repair proteins. These dynamics are expected to occur on similar timescales as one-dimensional (1D) diffusion of proteins on DNA if effective in stalling these proteins as they scan DNA. We examined the flexibility and dynamics of DNA oligomers containing 3 base pair (bp) mismatched sites specifically recognized in vitro by nucleotide excision repair protein Rad4 (yeast ortholog of mammalian XPC). A previous Forster resonance energy transfer (FRET) study mapped DNA conformational distributions with cytosine analog FRET pair primarily sensitive to DNA twisting/unwinding deformations (Chakraborty et al. Nucleic Acids Res. 46: 1240-1255 (2018)). These studies revealed B-DNA conformations for nonspecific (matched) constructs but significant unwinding for mismatched constructs specifically recognized by Rad4, even in the absence of Rad4. The timescales of these unwinding fluctuations, however, remained elusive. Here, we labeled DNA with Atto550/Atto647N FRET dyes suitable for fluorescence correlation spectroscopy (FCS). With these probes, we detected higher FRET in specific, mismatched DNA compared with matched DNA, reaffirming unwinding/bending deformations in mismatched DNA. FCS unveiled the dynamics of these spontaneous deformations at ~ 300 µs with no fluctuations detected for matched DNA within the ~ 600 ns-10 ms FCS time window. These studies are the first to visualize anomalous unwinding/bending fluctuations in mismatched DNA on timescales that overlap with the < 500 µs "stepping" times of repair proteins on DNA. Such "flexible hinge" dynamics at lesion sites could arrest a diffusing protein to facilitate damage interrogation and recognition.

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol-1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol-1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson-Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

G-Tetrad-Selective Ligand Binding Kinetics in G-Quadruplex DNA Probed with Fluorescence Correlation Spectroscopy

Probing the kinetics of ligand binding to biomolecules is of paramount interest in biology and pharmacology. Measurements of such kinetic processes provide information on the rate-determining steps that control the binding affinity of ligands to biomolecules, thereby predicting the mechanism of the molecular interaction. In this context, ligand binding to G-quadruplex DNA (GqDNA) structures has attracted tremendous attention primarily because of their use in possible anticancer therapy. Although a large number of G-quadruplex-specific ligands have been proposed, probing the kinetics of G-tetrad-selective binding of (multiple) ligands within a G-quadruplex DNA (GqDNA) structure remains challenging. Most of the earlier studies focused on the thermodynamics of ligand binding; however, the kinetics of ligand association and dissociation with GqDNA, particularly binding of multiple ligands within a GqDNA structure, have not been explored. Here, we propose a simple fluorescence correlation spectroscopy-based method that measures the G-tetrad-selective association and dissociation rates of ligands within a GqDNA structure by correlating the fluorescence fluctuations of a site-specific (5' or 3' end-labeled) fluorophore (Cy3) in GqDNA due to quenching of Cy3 fluorescence, induced by the ligand binding to the G-tetrads. We show that well-known GqDNA ligands, BRACO19, TMPyP4, Hoechst 33258, and Hoechst 33342, have G-tetrad-selective association and dissociation rates, which suggest site-dependent variation of free energy barriers for binding/unbinding of the ligands with GqDNA. We also show that the measured kinetic rates depend not only on the G-tetrad site (5' vs 3' end) but also on the ligand and GqDNA structures.

A second DNA binding site on RFC facilitates clamp loading at gapped or nicked DNA

Clamp loaders place circular sliding clamp proteins onto DNA so that clamp-binding partner proteins can synthesize, scan, and repair the genome. DNA with nicks or small single-stranded gaps are common clamp-loading targets in DNA repair, yet these substrates would be sterically blocked given the known mechanism for binding of primer-template DNA. Here, we report the discovery of a second DNA binding site in the yeast clamp loader replication factor C (RFC) that aids in binding to nicked or gapped DNA. This DNA binding site is on the external surface and is only accessible in the open conformation of RFC. Initial DNA binding at this site thus provides access to the primary DNA binding site in the central chamber. Furthermore, we identify that this site can partially unwind DNA to create an extended single-stranded gap for DNA binding in RFC’s central chamber and subsequent ATPase activation. Finally, we show that deletion of the BRCT domain, a major component of the external DNA binding site, results in defective yeast growth in the presence of DNA damage where nicked or gapped DNA intermediates occur. We propose that RFC’s external DNA binding site acts to enhance DNA binding and clamp loading, particularly at DNA architectures typically found in DNA repair.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

Cryo-EM structures reveal high-resolution mechanism of a DNA polymerase sliding clamp loader

Sliding clamps are ring-shaped protein complexes that are integral to the DNA replication machinery of all life. Sliding clamps are opened and installed onto DNA by clamp loader AAA+ ATPase complexes. However, how a clamp loader opens and closes the sliding clamp around DNA is still unknown. Here, we describe structures of the Saccharomyces cerevisiae clamp loader Replication Factor C (RFC) bound to its cognate sliding clamp Proliferating Cell Nuclear Antigen (PCNA) en route to successful loading. RFC first binds to PCNA in a dynamic, closed conformation that blocks both ATPase activity and DNA binding. RFC then opens the PCNA ring through a large-scale 'crab-claw' expansion of both RFC and PCNA that explains how RFC prefers initial binding of PCNA over DNA. Next, the open RFC:PCNA complex binds DNA and interrogates the primer-template junction using a surprising base-flipping mechanism. Our structures indicate that initial PCNA opening and subsequent closure around DNA do not require ATP hydrolysis, but are driven by binding energy. ATP hydrolysis, which is necessary for RFC release, is triggered by interactions with both PCNA and DNA, explaining RFC's switch-like ATPase activity. Our work reveals how a AAA+ machine undergoes dramatic conformational changes for achieving binding preference and substrate remodeling.

Scanning Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy: Conformational Dynamics of DNA Holliday Junction from Microsecond to Subsecond

Single-molecule Förster resonance energy transfer (smFRET) is widely utilized to investigate the structural heterogeneity and dynamics of biomolecules. However, it has been difficult to simultaneously achieve a wide observation time window, a high structure resolution, and a high time resolution with the current smFRET methods. Herein, we introduce a new method utilizing two-dimensional fluorescence lifetime correlation spectroscopy (2D FLCS) and surface immobilization techniques. This method, scanning 2D FLCS, enables us to examine the structural heterogeneity and dynamics of immobilized biomolecules on a time scale from microsecond to subsecond by slowly scanning the sample stage at the rate of ∼1 μm/s. Application to the DNA Holliday junction (HJ) complex under various [Mg2+] conditions demonstrates that scanning 2D FLCS enables tracking reaction kinetics from 25 μs to 30 ms with a time resolution as high as 1 μs. Furthermore, the high structure resolution of scanning 2D FLCS allows us to unveil the ensemble nature of each isomer state and the heterogeneity of the dynamics of the HJ.

Tandem Substitutions in Somatic Hypermutation

Upon antigen recognition, activation-induced cytosine deaminase initiates affinity maturation of the B-cell receptor by somatic hypermutation (SHM) through error-prone DNA repair pathways. SHM typically creates single nucleotide substitutions, but tandem substitutions may also occur. We investigated incidence and sequence context of tandem substitutions by massive parallel sequencing of V(D)J repertoires in healthy human donors. Mutation patterns were congruent with SHM-derived single nucleotide mutations, delineating initiation of the tandem substitution by AID. Tandem substitutions comprised 5,7% of AID-induced mutations. The majority of tandem substitutions represents single nucleotide juxtalocations of directly adjacent sequences. These observations were confirmed in an independent cohort of healthy donors. We propose a model where tandem substitutions are predominantly generated by translesion synthesis across an apyramidinic site that is typically created by UNG. During replication, apyrimidinic sites transiently adapt an extruded configuration, causing skipping of the extruded base. Consequent strand decontraction leads to the juxtalocation, after which exonucleases repair the apyramidinic site and any directly adjacent mismatched base pairs. The mismatch repair pathway appears to account for the remainder of tandem substitutions. Tandem substitutions may enhance affinity maturation and expedite the adaptive immune response by overcoming amino acid codon degeneracies or mutating two adjacent amino acid residues simultaneously.

Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond

Damaged or mismatched DNA bases are normally thought to be able to flip out of the helical stack, providing enzymes with access to the faulty genetic information otherwise hidden inside the helix. Thymine glycol (Tg) is one of the most common products of nucleic acid damage. However, the static and dynamic structures of DNA duplexes affected by 5R-Tg epimers are still not clearly understood, including the ability of these to undergo spontaneous base flipping. Structural effects of the 5R-Tg epimers on the duplex DNA are herein studied using molecular dynamics together with reliable DFT based calculations. In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA, and a barrier of 4.9 kcal mol^–1 is obtained by meta-eABF for cis-5R,6S-Tg base flipping out of the duplex DNA, comparable to the 5.4 kcal mol^–1 obtained for the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were identified, of which the most stable disrupts the Watson–Crick hydrogen-bonded G5/C20 base pair, leading to conformational distortion of the duplex. Interestingly, the relative barrier height of the 5R-Tg flipping is only 1.0 kcal mol^–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified to be essential for 5R-Tg flipping. The study clearly demonstrates the occurrence of partial trans-5R,6R-Tg epimer flipping in solution.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

Assessing the sequence dependence of pyrimidine-pyrimidone (6-4) photoproduct in a duplex double-stranded DNA: A pitfall for microsecond range simulation

Sequence dependence of the (6-4) photoproduct conformational landscape when embedded in six 25-bp duplexes is evaluated along extensive unbiased and enhanced (replica exchange with solute tempering, REST2) molecular dynamics simulations. The structural reorganization as the central pyrimidines become covalently tethered is traced back in terms of non-covalent interactions, DNA bending, and extrusion of adenines of the opposite strands. The close sequence pattern impacts the conformational landscape around the lesion, inducing different upstream and downstream flexibilities. Moreover, REST2 simulations allow us to probe structures possibly important for damaged DNA recognition.

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

Differentiating a Least-Stable Single Nucleotide Mismatch in DNA Via Metal Ion-Mediated Base Pairing and Using Thioflavin T as an Extrinsic Fluorophore

Monitoring the DNA dynamics in solution has great potential to develop new nucleic acid-based sensors and devices. With spectroscopic approaches, both at the ensemble average and single-molecule resolution, this study is directed to differentiate a single nucleotide mismatch (SNM) via a metal ion-stabilized mismatched base-pairing (C-Ag⁺-C/C-Cu²⁺-T) (C = cytosine, T = thymine) and site-selective extrinsic fluorophore, specifically, Thioflavin T (ThT). This is the first approach of its kind where dynamic quantities like molecular diffusion coefficients and diffusion times have been utilized to distinguish the least-stable SNM (CC & CT) formed by the most discriminating nucleobase, specifically, cytosine in a 20-mer duplex DNA. Additionally, this work also quantifies metal ions (Ag⁺ and Cu²⁺) at lower concentrations using fluorescence correlation spectroscopy. Our results can provide greater molecular-level insights into the mismatch-dependent metal-DNA interactions and also illuminate ThT as a new fluorophore to monitor the dynamics involved in DNA-metal composites.

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

Reaction Coordinate and Thermodynamics of Base Flipping in RNA

Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.

Molecular mechanism of networking among DegP, Skp and SurA in periplasm for biogenesis of outer membrane proteins

The biogenesis of outer membrane proteins (OMPs) is an extremely challenging process. In the periplasm of Escherichia coli, a group of quality control factors work together to exercise the safe-guard and quality control of OMPs. DegP, Skp and SurA are the three most prominent ones. Although extensive investigations have been carried out, the molecular mechanism regarding the networking among these proteins remains mostly mysterious. Our group has previously studied the molecular interactions of OMPs with SurA and Skp, using single-molecule detection (SMD). In this work, again using SMD, we studied how OmpC, a representative of OMPs, interacts with DegP, Skp and SurA collectively. Several important discoveries were made. The self-oligomerization of DegP to form hexamer occurs over hundred micromolars. When OmpC is in a monomer state at a low concentration, the OmpC·DegP6 and OmpC·DegP24 complexes form when the DegP concentration is around sub-micromolars and a hundred micromolars, respectively. High OmpC concentration promotes the binding affinity of DegP to OmpC by ∼100 folds. Skp and SurA behave differently when they interact synergistically with DegP in the presence of substrate. DegP can degrade SurA-protected OmpC, but Skp-protected OmpC forms the ternary complex OmpC·(Skp3)n·DegP6 (n = 1,2) to resist the DegP-mediated degradation. Combined with previous results, we were able to depict a comprehensive picture regarding the molecular mechanism of the networking among DegP, Skp and SurA in the periplasm for the OMPs biogenesis under physiological and stressed conditions.

Automated AFM analysis of DNA bending reveals initial lesion sensing strategies of DNA glycosylases

Base excision repair is the dominant DNA repair pathway of chemical modifications such as deamination, oxidation, or alkylation of DNA bases, which endanger genome integrity due to their high mutagenic potential. Detection and excision of these base lesions is achieved by DNA glycosylases. To investigate the remarkably high efficiency in target site search and recognition by these enzymes, we applied single molecule atomic force microscopy (AFM) imaging to a range of glycosylases with structurally different target lesions. Using a novel, automated, unbiased, high-throughput analysis approach, we were able to resolve subtly different conformational states of these glycosylases during DNA lesion search. Our results lend support to a model of enhanced lesion search efficiency through initial lesion detection based on altered mechanical properties at lesions. Furthermore, its enhanced sensitivity and easy applicability also to other systems recommend our novel analysis tool for investigations of diverse, fundamental biological interactions.

Ligand Recognition in Viral RNA Necessitates Rare Conformational Transitions

Ribonucleic acids (RNAs) are conformationally flexible molecules that fold into three-dimensional structures and play an important role in different cellular processes as well as in the development of many diseases. RNA has therefore become an important target for developing novel therapeutic approaches. The biophysical processes underlying RNA function are often associated with rare structural transitions that play a key role in ligand recognition. In this work, we probe these rarely occurring transitions using nonequilibrium simulations by characterizing the dissociation of a ligand molecule from an HIV-1 viral RNA element. Specifically, we observed base-flipping rare events that are coupled with ligand binding/unbinding and also provided mechanistic details underlying these transitions.

Single-Molecule Electrical Detection: A Promising Route toward the Fundamental Limits of Chemistry and Life Science

The ultimate limit of analytical chemistry is single-molecule detection, which allows one to visualize the dynamic processes of chemical/biological interactions with single-molecule or single-event sensitivity and hence enables the study of stochastic fluctuations under equilibrium conditions and the observation of time trajectories and reaction pathways of individual species in nonequilibrated systems. In addition, such studies may also allow the direct observation of novel microscopic quantum effects and fundamental discoveries of underlying molecular mechanisms in organic reactions and biological processes that are not accessible in ensemble experiments, thus providing unique opportunities to solve the key problems of physical, chemical, and life sciences. Consequently, the field of single-molecule detection has received considerable attention and has witnessed tremendous advances in different directions in combination with other disciplines. This Account describes our ongoing work on the development of groundbreaking methods (termed "single-molecule electrical approaches") of translating the detailed processes of chemical reactions or biological functions into detectable electrical signals at the single-event level on the platform of single-molecule electronic devices, with a particular focus on graphene-molecule-graphene single-molecule junctions (GMG-SMJs) and silicon-nanowire-based single-molecule electrical nanocircuits. These nanocircuit-based architectures are complementary to conventional optical or mechanical techniques but exhibit obvious advantages such as the absence of problems associated with bleaching and fluorescent labeling. Dash-line lithography (DLL) is an efficient lithographic method of cutting graphene and forming carboxylic-acid-functionalized nanogapped graphene point contact arrays developed to address the formidable challenges of molecular device fabrication difficulty and poor stability. Molecules of interest terminated by amines on both ends can be covalently sandwiched between graphene point contacts to create high-throughput robust GMG-SMJs containing only one molecule as the conductive element. In conjunction with the ease of device fabrication and device stability, this feature distinguishes GMG-SMJs as a new testbed platform for single-molecule analysis characterized by high temporal resolution and superior signal-to-noise ratios. By exploiting the DLL method, we have fabricated molecular devices that are sensitive to external stimuli and are capable of transducing chemical/biochemical events into electrical signals at the single-molecule level, with notable examples including host-guest interaction, hydrogen bond dynamics, DNA intercalation, photoinduced conformational transition, carbocation formation, nucleophilic addition, and stereoelectronic effect. In addition to GMG-SMJs and considering compatibility with the silicon-based industry, we have also developed a reliable method of point-functionalizing silicon-nanowire-based nanotransistors to afford single-molecule electrical nanocircuits. This approach proved to be a robust platform for single-molecule electrical analysis capable of probing fast dynamic processes such as single-protein detection, DNA hybridization/polymorphism, and motor rotation dynamics. The above systematic investigations emphasize the importance and unique advantages of universal single-molecule electrical approaches for realizing direct, label-free, real-time electrical measurements of reaction dynamics with single-event sensitivity. These approaches promise a fascinating mainstream platform to explore the dynamics of stochastic processes in chemical/biological systems as well as gain information in fields ranging from reaction chemistry for elucidating the intrinsic mechanisms to genomics or proteomics for accurate molecular and even point-of-care clinical diagnoses.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

Molecular Basis of Substrate Recognition of Endonuclease Q from the Euryarchaeon Pyrococcus furiosus

Endonuclease Q (EndoQ), a DNA repair endonuclease, was originally identified in the hyperthermophilic euryarchaeon Pyrococcus furiosus in 2015. EndoQ initiates DNA repair by generating a nick on DNA strands containing deaminated bases and an abasic site. Although EndoQ is thought to be important for maintaining genome integrity in certain bacteria and archaea, the underlying mechanism catalyzed by EndoQ remains unclear. Here, we provide insights into the molecular basis of substrate recognition by EndoQ from P. furiosus (PfuEndoQ) using biochemical approaches. Our results of the substrate specificity range and the kinetic properties of PfuEndoQ demonstrate that PfuEndoQ prefers the imide structure in nucleobases along with the discovery of its cleavage activity toward 5,6-dihydrouracil, 5-hydroxyuracil, 5-hydroxycytosine, and uridine in DNA. The combined results for EndoQ substrate binding and cleavage activity analyses indicated that PfuEndoQ flips the target base from the DNA duplex, and the cleavage activity is highly dependent on spontaneous base flipping of the target base. Furthermore, we find that PfuEndoQ has a relatively relaxed substrate specificity; therefore, the role of EndoQ in restriction modification systems was explored. The activity of the EndoQ homolog from Bacillus subtilis was found not to be inhibited by the uracil glycosylase inhibitor from B. subtilis bacteriophage PBS1, whose genome is completely replaced by uracil instead of thymine. Our findings suggest that EndoQ not only has additional functions in DNA repair but also could act as an antiviral enzyme in organisms with EndoQ.IMPORTANCE Endonuclease Q (EndoQ) is a lesion-specific DNA repair enzyme present in certain bacteria and archaea. To date, it remains unclear how EndoQ recognizes damaged bases. Understanding the mechanism of substrate recognition by EndoQ is important to grasp genome maintenance systems in organisms with EndoQ. Here, we find that EndoQ from the euryarchaeon Pyrococcus furiosus recognizes the imide structure in nucleobases by base flipping, and the cleavage activity is enhanced by the base pair instability of the target base, along with the discovery of its cleavage activity toward 5,6-dihydrouracil, 5-hydroxyuracil, 5-hydroxycytosine, and uridine in DNA. Furthermore, a potential role of EndoQ in Bacillus subtilis as an antiviral enzyme by digesting viral genome is demonstrated.

Polymorphic G:G mismatches act as hotspots for inducing right-handed Z DNA by DNA intercalation

DNA mismatches are highly polymorphic and dynamic in nature, albeit poorly characterized structurally. We utilized the antitumour antibiotic Co^II(Chro)₂ (Chro = chromomycin A3) to stabilize the palindromic duplex d(TTGGCGAA) DNA with two G:G mismatches, allowing X-ray crystallography-based monitoring of mismatch polymorphism. For the first time, the unusual geometry of several G:G mismatches including syn–syn, water mediated anti–syn and syn–syn-like conformations can be simultaneously observed in the crystal structure. The G:G mismatch sites of the d(TTGGCGAA) duplex can also act as a hotspot for the formation of alternative DNA structures with a GC/GA-5′ intercalation site for binding by the GC-selective intercalator actinomycin D (ActiD). Direct intercalation of two ActiD molecules to G:G mismatch sites causes DNA rearrangements, resulting in backbone distortion to form right-handed Z-DNA structures with a single-step sharp kink. Our study provides insights on intercalators-mismatch DNA interactions and a rationale for mismatch interrogation and detection via DNA intercalation.

Enhanced sampling in molecular dynamics

Although molecular dynamics simulations have become a useful tool in essentially all fields of chemistry, condensed matter physics, materials science, and biology, there is still a large gap between the time scale which can be reached in molecular dynamics simulations and that observed in experiments. To address the problem, many enhanced sampling methods were introduced, which effectively extend the time scale being approached in simulations. In this perspective, we review a variety of enhanced sampling methods. We first discuss collective-variables-based methods including metadynamics and variationally enhanced sampling. Then, collective variable free methods such as parallel tempering and integrated tempering methods are presented. At last, we conclude with a brief introduction of some newly developed combinatory methods. We summarize in this perspective not only the theoretical background and numerical implementation of these methods but also the new challenges and prospects in the field of the enhanced sampling.

Base-flipping dynamics from an intrahelical to an extrahelical state exerted by thymine DNA glycosylase during DNA repair process

Thymine DNA glycosylase (TDG) is a DNA repair enzyme that excises a variety of mismatched or damaged nucleotides (nts), e.g. dU, dT, 5fC and 5caC. TDG is shown to play essential roles in maintaining genome integrity and correctly programming epigenetic modifications through DNA demethylation. After locating the lesions, TDG employs a base-flipping strategy to recognize the damaged nucleobases, whereby the interrogated nt is extruded from the DNA helical stack and binds into the TDG active site. The dynamic mechanism of the base-flipping process at an atomistic resolution, however, remains elusive. Here, we employ the Markov State Model (MSM) constructed from extensive all-atom molecular dynamics (MD) simulations to reveal the complete base-flipping process for a G.T mispair at a tens of microsecond timescale. Our studies identify critical intermediates of the mispaired dT during its extrusion process and reveal the key TDG residues involved in the inter-state transitions. Notably, we find an active role of TDG in promoting the intrahelical nt eversion, sculpturing the DNA backbone, and penetrating into the DNA minor groove. Three additional TDG substrates, namely dU, 5fC, and 5caC, are further tested to evaluate the substituent effects of various chemical modifications of the pyrimidine ring on base-flipping dynamics.

Determination of Base-Flipping Free-Energy Landscapes from Nonequilibrium Stratification

Correct calculation of the variation of free energy upon base flipping is crucial in understanding the dynamics of DNA systems. The free-energy landscape along the flipping pathway gives the thermodynamic stability and the flexibility of base-paired states. Although numerous free-energy simulations are performed in the base flipping cases, no theoretically rigorous nonequilibrium techniques are devised and employed to investigate the thermodynamics of base flipping. In the current work, we report a general nonequilibrium stratification scheme for the efficient calculation of the free-energy landscape of base flipping in DNA duplex. We carefully monitor the convergence behavior of the equilibrium sampling based free-energy simulation and the nonequilibrium stratification and determine the empirical length of time blocks required for converged sampling. Comparison between the performances of the equilibrium umbrella sampling and the nonequilibrium stratification is given. The results show that nonequilibrium free-energy simulation achieves similar accuracy and efficiency compared with the equilibrium enhanced sampling technique in the base flipping cases. We further test a convergence criterion we previously proposed and it comes out that the convergence determined by this criterion agrees with those given by the time-invariant behavior of PMF and the nonlinear dependence of standard deviation on the sample size.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Assessing AMBER force fields for protein folding in an implicit solvent

Molecular dynamics (MD) simulation implemented with a state-of-the-art protein force field and implicit solvent model is an attractive approach to investigate protein folding, one of the most perplexing problems in molecular biology. But how well can force fields developed independently of implicit solvent models work together in reproducing diverse protein native structures and measuring the corresponding folding thermodynamics is not always clear. In this work, we performed enhanced sampling MD simulations to assess the ability of six AMBER force fields (FF99SBildn, FF99SBnmr, FF12SB, FF14ipq, FF14SB, and FF14SBonlysc) as coupled with a recently improved pair-wise GB-Neck2 model in modeling the folding of two helical and two β-sheet peptides. Whilst most of the tested force fields can yield roughly similar features for equilibrium conformational ensembles and detailed folding free-energy profiles for short α-helical TC10b in an implicit solvent, the measured counterparts are significantly discrepant in the cases of larger or β-structured peptides (HP35, 1E0Q, and GTT). Additionally, the calculated folding/unfolding thermodynamic quantities can only partially match the experimental data. Although a combination of the force fields and GB-Neck2 implicit model able to describe all aspects of the folding transitions towards the native structures of all the considered peptides was not identified, we found that FF14SBonlysc coupled with the GB-Neck2 model seems to be a reasonably balanced combination to predict peptide folding preferences.

Role of Conformational Fluctuations of Protein toward Methylation in DNA by Cytosine-5-methyltransferase

Methylation of cytosine is the common epigenetic modification in genomes ranging from bacteria to mammals, and aberrant methylation leads to human diseases including cancer. Recognition of a cognate DNA sequence by DNA methyltransferases and flipping of a target base into the enzyme active site pocket are the key steps in DNA methylation. Using molecular dynamics simulations and enhanced sampling techniques here we elucidate the role of conformational fluctuations of protein and active or passive involvement of protein elements that mediate base flipping and formation of the closed catalytic complex. The free energy profiles for the flipping of target cytosine into the enzyme active site support the major groove base eversion pathway; and the results show that the closed state of enzyme increases the free energy barrier, whereas the open state reduces it. We found that the interactions of the key loop residues of protein with cognate DNA altered the protein motions, and modulation of protein fluctuations relates to the closed catalytic complex formation. Methylation of cytosine in the active site of the closed complex destabilizes the interactions of catalytic loop residues with cognate DNA and reduces the stability of the closed state. Our study provides microscopic insights on the base flipping mechanism coupled with enzyme's loop motions and provides evidence for the role of conformational fluctuations of protein in the enzyme-catalyzed DNA processing mechanism.

Interactions of the DNA Repair Enzyme Human Thymine DNA Glycosylase with Cognate and Noncognate DNA

Glycosylases specifically recognize and flip their target base out of the DNA helix into the enzyme's active site. Our simulations show that a partially flipped state, already present in free DNA carrying a T:G mispair, becomes the more probable state compared to the closed state after binding of thymine DNA glycosylase (TDG). Paired thymine (T:A) or methyl-cytosine (mC:G) does not exhibit a partially flipped state in free or complexed DNA. Important enzyme-DNA interactions exhibit significant strength in the intrahelical and extrahelical TDG-DNA complexes. The computed binding free energy differences suggest these interactions account for the stabilization of the partially flipped state, thereby driving the T:G mispair toward base flip. In the fully flipped state, the cognate base thymine is significantly better accommodated in the enzyme's active site than noncognate bases are, suggesting the hydrolysis step as the last of several stages at which base recognition can be achieved.

Energetics of base flipping at a DNA mismatch site confined at the latch constriction of α-hemolysin

Unique, two-state modulating current signatures are observed when a cytosine-cytosine mismatch pair is confined at the 2.4 nm latch constriction of the α-hemolysin (αHL) nanopore. We have previously speculated that the modulation is due to base flipping at the mismatch site. Base flipping is a biologically significant mechanism in which a single base is rotated out of the DNA helical stack by 180°. It is the mechanism by which enzymes are able to access bases for repair operations without disturbing the global structure of the helix. Here, temperature dependent ion channel recordings of individual double-stranded DNA duplexes inside αHL are used to derive thermodynamic (ΔH, ΔS) and kinetic (E_A) parameters for base flipping of a cytosine at an unstable cytosine-cytosine mismatch site. The measured activation energy for flipping a cytosine located at the latch of αHL out of the helix (18 ± 1 kcal mol^-1) is comparable to that previously reported for base flipping at mismatch sites from NMR measurements and potential mean force calculations. We propose that the αHL nanopore is a useful tool for measuring conformational changes in dsDNA at the single molecule level.

The effects of implicit modeling of nonpolar solvation on protein folding simulations

Implicit solvent models, in which the polar and nonpolar solvation free-energies of solute molecules are treated separately, have been widely adopted for molecular dynamics simulation of protein folding. While the development of the implicit models is mainly focused on the methodological improvement and key parameter optimization for polar solvation, nonpolar solvation has been either ignored or described by a simplistic surface area (SA) model. In this work, we performed the folding simulations of multiple β-hairpin and α-helical proteins with varied surface tension coefficients embedded in the SA model to clearly demonstrate the effects of nonpolar solvation treated by a popular SA model on protein folding. The results indicate that the change in the surface tension coefficient does not alter the ability of implicit solvent simulations to reproduce a protein native structure but indeed controls the components of the equilibrium conformational ensemble and modifies the energetic characterization of the folding transition pathway. The suitably set surface tension coefficient can yield explicit solvent simulations and/or experimentally suggested folding mechanism of protein. In addition, the implicit treatment of both polar and nonpolar components of solvation free-energy contributes to the overestimation of the secondary structure in implicit solvent simulations.

Surface Transient Binding-Based Fluorescence Correlation Spectroscopy (STB-FCS), a Simple and Easy-to-Implement Method to Extend the Upper Limit of the Time Window to Seconds

Fluorescence correlation spectroscopy is a powerful single-molecule tool that is able to capture kinetic processes occurring at the nanosecond time scale. However, the upper limit of its time window is restricted by the dwell time of the molecule of interest in the confocal detection volume, which is usually around submilliseconds for a freely diffusing biomolecule. Here, we present a simple and easy-to-implement method, named surface transient binding-based fluorescence correlation spectroscopy (STB-FCS), which extends the upper limit of the time window to seconds. We further demonstrated that STB-FCS enables capture of both intramolecular and intermolecular kinetic processes whose time scales cross several orders of magnitude.

Single-Molecule Titration in a Protein Nanoreactor Reveals the Protonation/Deprotonation Mechanism of a C:C Mismatch in DNA

Measurement of single-molecule reactions can elucidate microscopic mechanisms that are often hidden from ensemble analysis. Herein, we report the acid-base titration of a single DNA duplex confined within the wild-type α-hemolysin (α-HL) nanopore for up to 3 h, while monitoring the ionic current through the nanopore. Modulation between two states in the current-time trace for duplexes containing the C:C mismatch in proximity to the latch constriction of α-HL is attributed to the base flipping of the C:C mismatch. As the pH is lowered, the rate for the C:C mismatch to flip from the intra-helical state to the extra-helical state ( k_intra-extra) decreases, while the rate for base flipping from the extra-helical state to the intra-helical state ( k_extra-intra) remains unchanged. Both k_intra-extra and k_extra-intra are on the order of 1 × 10^-2 s^-1 to 1 × 10^-1 s^-1 and remain stable over the time scale of the measurement (several hours). Analysis of the pH-dependent kinetics of base flipping using a hidden Markov kinetic model demonstrates that protonation/deprotonation occurs while the base pair is in the intra-helical state. We also demonstrate that the rate of protonation is limited by transport of H⁺ into the α-HL nanopore. Single-molecule kinetic isotope experiments exhibit a large kinetic isotope effect (KIE) for k_intra-extra ( k_H/ k_D ≈ 5) but a limited KIE for k_extra-intra ( k_H/ k_D ≈ 1.3), supporting our model. Our experiments correspond to the longest single-molecule measurements performed using a nanopore, and demonstrate its application in interrogating mechanisms of single-molecule reactions in confined geometries.

Direct single-molecule dynamic detection of chemical reactions

Single-molecule detection can reveal time trajectories and reaction pathways of individual intermediates/transition states in chemical reactions and biological processes, which is of fundamental importance to elucidate their intrinsic mechanisms. We present a reliable, label-free single-molecule approach that allows us to directly explore the dynamic process of basic chemical reactions at the single-event level by using stable graphene-molecule single-molecule junctions. These junctions are constructed by covalently connecting a single molecule with a 9-fluorenone center to nanogapped graphene electrodes. For the first time, real-time single-molecule electrical measurements unambiguously show reproducible large-amplitude two-level fluctuations that are highly dependent on solvent environments in a nucleophilic addition reaction of hydroxylamine to a carbonyl group. Both theoretical simulations and ensemble experiments prove that this observation originates from the reversible transition between the reactant and a new intermediate state within a time scale of a few microseconds. These investigations open up a new route that is able to be immediately applied to probe fast single-molecule physics or biophysics with high time resolution, making an important contribution to broad fields beyond reaction chemistry.