Mismatched base-pair simulations for ASFV Pol X/DNA complexes help interpret frequent G*G misincorporation

Genetic Variations of African Swine Fever Virus: Major Challenges and Prospects

African swine fever (ASF) is a contagious viral disease affecting pigs and wild boars. It typically presents as a hemorrhagic fever but can also manifest in various forms, ranging from acute to asymptomatic. ASF has spread extensively globally, significantly impacting the swine industry. The complex and highly variable character of the ASFV genome makes vaccine development and disease surveillance extremely difficult. The overall trend in ASFV evolution is towards decreased virulence and increased transmissibility. Factors such as gene mutation, viral recombination, and the strain-specificity of virulence-associated genes facilitate viral variations. This review deeply discusses the influence of these factors on viral immune evasion, pathogenicity, and the ensuing complexities encountered in vaccine development, disease detection, and surveillance. The ultimate goal of this review is to thoroughly explore the genetic evolution patterns and variation mechanisms of ASFV, providing a theoretical foundation for advancement in vaccine and diagnostic technologies.

Genome Plasticity of African Swine Fever Virus: Implications for Diagnostics and Live-Attenuated Vaccines

African swine fever (ASF) is a highly contagious transboundary viral hemorrhagic disease of domestic and wild pigs presenting a significant threat to the global swine industry. Following its introduction in Caucasus, Georgia, in 2007, the genome of the genotype II of African swine fever virus (ASFV) strain Georgia-07 and its derivatives accumulated significant mutations, resulting in the emergence of genetic variants within short epidemiological timescales as it spreads and infects different hosts in diverse ecosystems, causing outbreaks in Europe, South Asia, South East Asia and the Caribbean. This suggests that ASFV, with a comparatively large and complex DNA genome, is susceptible to genetic mutations including deletions and that although the virus is environmentally stable, it is genetically unstable. This has implications for the development of vaccines and diagnostic tests for disease detection and surveillance. Analysis of the ASFV genome revealed recombination hotspots, which in double-stranded DNA (dsDNA) viruses represent key drivers of genetic diversity. The ability of pox virus, a dsDNA virus with a genome complexity similar to ASFV, regaining virulence following the deletion of a virulence gene via gene amplification, coupled with the recent emergence and spread of live-attenuated ASFV vaccine strains causing disease and death in pigs in China, raise legitimate concerns around the use of live-attenuated ASFV vaccines in non-endemic regions to control the potential introduction. Further research into the risk of using live-attenuated ASFV in non-endemic regions is highly needed.

Comparative Analysis of Full Genome Sequences of African Swine Fever Virus Isolates Taken from Wild Boars in Russia in 2019

In this study, we report on the full genome phylogenetic analysis of four ASFV isolates obtained from wild boars in Russia. These samples originated from two eastern and two western regions of Russia in 2019. Phylogenetic analysis indicated that the isolates were assigned to genotype II and grouped according to their geographical origins. The two eastern isolates shared 99.99% sequence identity with isolates from China, Poland, Belgium, and Moldova, whereas the western isolates had 99.98% sequence identity with isolates from Lithuania and the original Georgia 2007 isolate. Based on the full genome phylogenies, we identified three single locus targets, MGF-360-10L, MGF-505-9R, and I267L, that yielded the same resolving power as the full genomes. The ease of alignment and a high level of variation make these targets a suitable selection as additional molecular markers in future ASFV phylogenetic practices.

Intra-epidemic genome variation in highly pathogenic African swine fever virus (ASFV) from the country of Georgia

BACKGROUND

African swine fever virus (ASFV) causes an acute hemorrhagic infection in suids with a mortality rate of up to 100%. No vaccine is available and the potential for catastrophic disease in Europe remains elevated due to the ongoing ASF epidemic in Russia and Baltic countries. To date, intra-epidemic whole-genome variation for ASFV has not been reported. To provide a more comprehensive baseline for genetic variation early in the ASF outbreak, we sequenced two Georgian ASFV samples, G-2008/1 and G-2008/2, derived from domestic porcine blood collected in 2008.

METHODS

Genomic DNA was extracted directly from low-volume ASFV PCR-positive porcine blood samples and subjected to next generation sequencing on the Illumina Miseq platform. De novo and mapped sequence assemblies were performed using CLCBio software. Genomic illustrations, sequence alignments and assembly figures were generated using Geneious v10.2.4. Sequence repeat architecture was analyzed using DNASTAR GeneQuest 14.1.0.

RESULTS

The G-2008/1 and G-2008/2 genomes were distinguished from each other by coding changes in seven genes, including MGF 110-1 L, X69R, MGF 505-10R, EP364R, H233R, E199L, and MGF 360-21R in addition to eight homopolymer tract variations. The 2008/2 genome possessed a novel allele state at a previously undescribed intergenic repeat locus between genes C315R and C147L. The C315R/C147L locus represents the earliest observed variable repeat sequence polymorphism reported among isolates from this epidemic. No sequence variation was observed in conventional ASFV subtyping markers. The two genomes exhibited complete collinearity and identical gene content with the Georgia 2007/1 reference genome. Approximately 56 unique homopolymer A/T-tract variations were identified that were unique to the Georgia 2007/1 genome. In both 2008 genomes, within-sample sequence read heterogeneity was evident at six homopolymeric G/C-tracts confined to the known hypervariable ~ 7 kb region in the left terminal region of the genome.

CONCLUSIONS

This is the first intra-epidemic comparative genomic analysis reported for ASFV and provides insight into the intra-epidemic microevolution of ASFV. The genomes reported here, in addition to the G-2007/1 genome, provide an early baseline for future genome-level comparisons and epidemiological tracing efforts.

Conformational flexibility and base-pairing tendency of the tobacco carcinogen O6-[4-oxo-4-(3-pyridyl)butyl]guanine

The present work uses DFT calculations to characterize the conformational and hydrogen-bonding properties of O6-[4-oxo-4-(3-pyridyl)butyl]guanine (POB-G), a DNA adduct caused by tobacco. POB-G is found to adopt many isoenergetic conformations that allow for discrete interactions between the bulky moiety and the adducted G and/or pairing base. The calculated structure and stability of the hydrogen-bonded pairs between the Watson-Crick or Hoogsteen face of POB-G and the canonical DNA nucleobases fully rationalize the previously reported mutational spectra. Specifically, the stable, non-distorted pseudo-Watson-Crick POB-G:T pair explains the predominant G➔A mutations, while the stable, yet marginally distorted pairs between the Watson-Crick face of POB-G and A or C clarify the G➔T mutations and non-mutagenic replication. Finally, the stable, yet highly distorted Hoogsteen POB-G:G pair rationalizes the experimentally-observed insertion but lack of persistence of G opposite POB-G in DNA. Overall, these structural insights are critical for guiding future studies that strive to fully understand the adduct mutagenicity, including the accessible conformations and the replication of POB-G-adducted DNA.

Fidelity of Nucleotide Incorporation by the RNA-Dependent RNA Polymerase from Poliovirus

Using poliovirus (PV) and its RNA-dependent RNA polymerase (RdRp) as our primary model system, we have advanced knowledge fundamental to the chemistry and fidelity of nucleotide addition by nucleic acid polymerase. Two fidelity checkpoints exist prior to nucleotide addition. The first toggles the enzyme between a nucleotide binding-occluded state and a nucleotide binding-competent state. The second represents an ensemble of conformational states of conserved structural motifs that permits retention of the incoming nucleotide in a state competent for phosphoryl transfer long enough for chemistry to occur. Nucleophilic attack of the alpha-phosphorous atom of the incoming nucleotide produces a pentavalent transition state, collapse of which is facilitated by protonation of the pyrophosphate leaving group by a general acid. All of the relevant conformational states of the enzyme are controlled by a network of interacting residues that permits remote-site residues to control active-site function. The current state of the art for PV RdRp enzymology is such that mechanisms governing fidelity of this enzyme can now be targeted genetically and chemically for development of attenuated viruses and antiviral agents, respectively. Application of the knowledge obtained with the PV RdRp to the development of vaccines and antivirals for emerging RNA viruses represents an important goal for the future.

How DNA polymerase X preferentially accommodates incoming dATP opposite 8-oxoguanine on the template

The modified base 8-oxo-7,8-dihydro-2'-deoxyguanosine (oxoG) is a common DNA adduct produced by the oxidation of DNA by reactive oxygen species. Kinetic data reveal that DNA polymerase X (pol X) from the African swine fever virus incorporates adenine (dATP) opposite to oxoG with higher efficiency than the non-damaged G:C basepair. To help interpret the kinetic data, we perform molecular dynamics simulations of pol X/DNA complexes, in which the template base opposite to the incoming dNTP (dCTP, dATP, dGTP) is oxoG. Our results suggest that pol X accommodates the oxoGsyn:A mispair by sampling closed active conformations that mirror those observed in traditional Watson-Crick complexes. Moreover, for both the oxoGsyn:A and oxoG:C ternary complexes, conformational sampling of the polymerase follows previously described large subdomain movements, local residue motions, and active site reorganization. Interestingly, the oxoGsyn:A system exhibits superior active site geometry in comparison to the oxoG:C system. Simulations for the other mismatch basepair complexes reveal large protein subdomain movement for all systems, except for oxoG:G, which samples conformations close to the open state. In addition, active site geometry and basepairing of the template base with the incoming nucleotide, reveal distortions and misalignments that range from moderate (i.e., oxoG:Asyn) to extreme (i.e., oxoGanti/syn:G). These results agree with the available kinetic data for pol X and provide structural insights regarding the mechanism by which this polymerase can accommodate incoming nucleotides opposite oxoG. Our simulations also support the notion that α-helix E is involved both in DNA binding and active site stabilization. Our proposed mechanism by which pol X can preferentially accommodate dATP opposite template oxoG further underscores the role that enzyme dynamics and conformational sampling operate in polymerase fidelity and function.

Atomistic nature of the DPT tautomerisation of the biologically important C·C* DNA base mispair containing amino and imino tautomers of cytosine: a QM and QTAIM approach

A theoretical study of tautomerisation of the biologically important cytosine·cytosine* (C·C*) DNA mismatch with a propeller-like structure (|C4N3N3C4| = 32.4°; C1 symmetry) and cis-oriented N1H glycosidic bonds, formed by the amino and imino tautomers of the C nucleobase, via the asynchronous concerted double proton transfer (DPT) along two H-bonds through the transition state (TSC·C*↔C*·C) (|C4N3N3C4| = 48.5°; C1 symmetry) into the C*·C mispair was carried out for the first time. It was established that the C·C*/C*·C DNA base mispair is associated by the antiparallel N4H···N4 (6.66 kcal mol(-1)), N3H···N3 (6.47 kcal mol(-1)) H-bonds and the O2···O2 van der Waals (vdW) contact (0.33 kcal mol(-1)), while the zwitterionic TSC·C*↔C*·C is stabilized by the parallel N4(+)H···N4(-) (13.55 kcal mol(-1)), N3(+)H···N3(-) (13.20 kcal mol(-1)) H-bonds and the O2(+)···O2(-) vdW contact (0.60 kcal mol(-1)). It was shown that the C·C* ↔ C*·C tautomerisation via the DPT is assisted by the O2···O2 vdW contact, that in contrast to the two others N4H···N4 and N3H···N3 H-bonds exists along the entire intrinsic reaction coordinate (IRC) range. The positive values of the Grunenberg's compliance constants (30.919 and 21.384 Å mdyn(-1) for C·C*/C*·C and TSC·C*↔C*·C, respectively) indicate that the O2···O2 vdW contact is a stabilizing closed-shell interaction. It was found that the middle N3H···N3 H-bond is anti-cooperative with the upper N4H···N4 H-bond and cooperative with the lower O2···O2 vdW contact. The 9 key points, which can be considered as electron-topological "fingerprints" of the asynchronous concerted C·C* ↔ C*·C tautomerisation process via the DPT were revealed along the IRC and examined in detail. It was shown that the C·C*/C*·C base mispair is a thermodynamically and dynamically stable structure. Its lifetime is equal to 1.53 × 10(-7) s at the MP2/cc-pVQZ//B3LYP/6-311++G(d,p) level of theory in vacuum. All 6 low-frequency intermolecular vibrations are able to develop during this time span.

How a low-fidelity DNA polymerase chooses non-Watson-Crick from Watson-Crick incorporation

A dogma for DNA polymerase catalysis is that the enzyme binds DNA first, followed by MgdNTP. This mechanism contributes to the selection of correct dNTP by Watson-Crick base pairing, but it cannot explain how low-fidelity DNA polymerases overcome Watson-Crick base pairing to catalyze non-Watson-Crick dNTP incorporation. DNA polymerase X from the deadly African swine fever virus (Pol X) is a half-sized repair polymerase that catalyzes efficient dG:dGTP incorporation in addition to correct repair. Here we report the use of solution structures of Pol X in the free, binary (Pol X:MgdGTP), and ternary (Pol X:DNA:MgdGTP with dG:dGTP non-Watson-Crick pairing) forms, along with functional analyses, to show that Pol X uses multiple unprecedented strategies to achieve the mutagenic dG:dGTP incorporation. Unlike high fidelity polymerases, Pol X can prebind purine MgdNTP tightly and undergo a specific conformational change in the absence of DNA. The prebound MgdGTP assumes an unusual syn conformation stabilized by partial ring stacking with His115. Upon binding of a gapped DNA, also with a unique mechanism involving primarily helix αE, the prebound syn-dGTP forms a Hoogsteen base pair with the template anti-dG. Interestingly, while Pol X prebinds MgdCTP weakly, the correct dG:dCTP ternary complex is readily formed in the presence of DNA. H115A mutation disrupted MgdGTP binding and dG:dGTP ternary complex formation but not dG:dCTP ternary complex formation. The results demonstrate the first solution structural view of DNA polymerase catalysis, a unique DNA binding mode, and a novel mechanism for non-Watson-Crick incorporation by a low-fidelity DNA polymerase.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.

Atomistic understanding of the C·T mismatched DNA base pair tautomerization via the DPT: QM and QTAIM computational approaches

It was established that the cytosine·thymine (C·T) mismatched DNA base pair with cis-oriented N1H glycosidic bonds has propeller-like structure (|N3C4C4N3| = 38.4°), which is stabilized by three specific intermolecular interactions-two antiparallel N4H…O4 (5.19 kcal mol(-1)) and N3H…N3 (6.33 kcal mol(-1)) H-bonds and a van der Waals (vdW) contact O2…O2 (0.32 kcal mol(-1)). The C·T base mispair is thermodynamically stable structure (ΔG(int) = -1.54 kcal mol(-1) ) and even slightly more stable than the A·T Watson-Crick DNA base pair (ΔG(int) = -1.43 kcal mol(-1)) at the room temperature. It was shown that the C·T ↔ C*·T* tautomerization via the double proton transfer (DPT) is assisted by the O2…O2 vdW contact along the entire range of the intrinsic reaction coordinate (IRC). The positive value of the Grunenberg's compliance constants (31.186, 30.265, and 22.166 Å/mdyn for the C·T, C*·T*, and TS(C·T ↔ C*·T*), respectively) proves that the O2…O2 vdW contact is a stabilizing interaction. Based on the sweeps of the H-bond energies, it was found that the N4H…O4/O4H…N4, and N3H…N3 H-bonds in the C·T and C*·T* base pairs are anticooperative and weaken each other, whereas the middle N3H…N3 H-bond and the O2…O2 vdW contact are cooperative and mutually reinforce each other. It was found that the tautomerization of the C·T base mispair through the DPT is concerted and asynchronous reaction that proceeds via the TS(C·T ↔ C*·T*) stabilized by the loosened N4-H-O4 covalent bridge, N3H…N3 H-bond (9.67 kcal mol(-1) ) and O2…O2 vdW contact (0.41 kcal mol(-1) ). The nine key points, describing the evolution of the C·T ↔ C*·T* tautomerization via the DPT, were detected and completely investigated along the IRC. The C*·T* mispair was revealed to be the dynamically unstable structure with a lifetime 2.13·× 10(-13) s. In this case, as for the A·T Watson-Crick DNA base pair, activates the mechanism of the quantum protection of the C·T DNA base mispair from its spontaneous mutagenic tautomerization through the DPT.

"Gate-keeper" residues and active-site rearrangements in DNA polymerase μ help discriminate non-cognate nucleotides

Incorporating the cognate instead of non-cognate substrates is crucial for DNA polymerase function. Here we analyze molecular dynamics simulations of DNA polymerase μ (pol μ) bound to different non-cognate incoming nucleotides including A:dCTP, A:dGTP, A(syn):dGTP, A:dATP, A(syn):dATP, T:dCTP, and T:dGTP to study the structure-function relationships involved with aberrant base pairs in the conformational pathway; while a pol μ complex with the A:dTTP base pair is available, no solved non-cognate structures are available. We observe distinct differences of the non-cognate systems compared to the cognate system. Specifically, the motions of active-site residue His329 and Asp330 distort the active site, and Trp436, Gln440, Glu443 and Arg444 tend to tighten the nucleotide-binding pocket when non-cognate nucleotides are bound; the latter effect may further lead to an altered electrostatic potential within the active site. That most of these “gate-keeper” residues are located farther apart from the upstream primer in pol μ, compared to other X family members, also suggests an interesting relation to pol μ's ability to incorporate nucleotides when the upstream primer is not paired. By examining the correlated motions within pol μ complexes, we also observe different patterns of correlations between non-cognate systems and the cognate system, especially decreased interactions between the incoming nucleotides and the nucleotide-binding pocket. Altered correlated motions in non-cognate systems agree with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies propose the following order for difficulty of non-cognate system insertions by pol μ: T:dGTP<A(syn):dATP<T:dCTP<A:dGTP<A(syn):dGTP<A:dCTP<A:dATP. This sequence agrees with available kinetic data for non-cognate nucleotide insertions, with the exception of A:dGTP, which may be more sensitive to the template sequence. The structures and conformational aspects predicted here are experimentally testable.

DNA polymerase μ (pol μ) is an enzyme that participates in DNA repair and thus has a central role in maintaining the integrity of genetic information. To efficiently repair the DNA, discriminating the cognate instead of non-cognate nucleotides (“fidelity-checking”) is required. Here we analyze molecular dynamics simulations of pol μ bound to different non-cognate nucleotides to study the structure-function relationships involved in the fidelity-checking mechanism of pol μ on the atomic level. Our results suggest that His329, Asp330, Trp436, Gln440, Glu443, and Arg444 are of great importance for pol μ's fidelity-checking mechanism. We also observe altered patterns of correlated motions within pol μ complex when non-cognate instead of cognate nucleotides are bound, which agrees with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies help interpret the available kinetic data of various non-cognate nucleotide insertions by pol μ. We also suggest experimentally testable predictions; for example, a point mutation like E443M may reduce the ability of pol μ to insert the cognate more than of non-cognate nucleotides. Our studies suggest an interesting relation to pol μ's unique ability to incorporate nucleotides when the upstream primer is not paired.

Relationship between conformational changes in pol lambda's active site upon binding incorrect nucleotides and mismatch incorporation rates

The correct replication and repair of DNA is critical for a cell's survival. Here, we investigate the fidelity of mammalian DNA polymerase lambda (pol lambda) utilizing dynamics simulation of the enzyme bound to incorrect incoming nucleotides including A:C, A:G, A(syn):G, A:A, A(syn):A, and T:G, all of which exhibit differing incorporation rates for pol lambda as compared to A:T bound to pol lambda. The wide range of DNA motion and protein residue side-chain motions observed in the mismatched systems demonstrates distinct differences when compared to the reference (correct base pair) system. Notably, Arg517's interactions with the DNA template strand bases in the active site are more limited, and Arg517 displays increased interactions with the incorrect dNTPs. This effect suggests that Arg517 helps provide a base-checking mechanism to discriminate correct from incorrect dNTPs. In addition, we find Tyr505 and Phe506 also play key roles in this base checking. A survey of the electrostatic potential landscape of the active sites and concomitant changes in electrostatic interaction energy between Arg517 and the dNTPs reveals that pol lambda binds incorrect dNTPs less tightly than the correct dNTP. These trends lead us to propose the following order for mismatch insertion by pol lambda: A:C > A:G > A(syn):G > T:G > A(syn):A > A:A. This sequence agrees with available kinetic data for incorrect nucleotide insertion opposite template adenine, with the exception of T:G, which may be more sensitive to the insertion context.