Motif-VI loop acts as a nucleotide valve in the West Nile Virus NS3 Helicase.

The Orthoflavivirus NS3 helicase (NS3h) is crucial in virus replication, representing a potential drug target for pathogenesis. NS3h utilizes nucleotide triphosphate (ATP) for hydrolysis energy to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. Intermediate states along the ATP hydrolysis cycle and conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. Extensive molecular dynamics simulations of West Nile virus NS3h+ssRNA in the apo, ATP, ADP+Pi and ADP bound states were used to model the conformational ensembles along this cycle. Energetic and structural clustering analyses depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). Based on these results, MVIL mutants (D471L, D471N and D471E) were found to have a substantial reduction in ATPase activity and RNA replication compared to the wild-type. Simulations of the mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

Motif V regulates energy transduction between the flavivirus NS3 ATPase and RNA-binding cleft.

The unwinding of dsRNA intermediates is critical for the replication of flavivirus RNA genomes. This activity is provided by the C-terminal helicase domain of viral nonstructural protein 3 (NS3). As a member of the superfamily 2 (SF2) helicases, NS3 requires the binding and hydrolysis of ATP/NTP to translocate along and unwind double-stranded nucleic acids. However, the mechanism of energy transduction between the ATP- and RNA-binding pockets is not well-understood. Previous molecular dynamics simulations conducted by our group have identified Motif V as a potential "communication hub" for this energy transduction pathway. To investigate the role of Motif V in this process, here we combined molecular dynamics, biochemistry, and virology approaches. We tested Motif V mutations in both the replicon and recombinant protein systems to investigate viral genome replication, RNA-binding affinity, ATP hydrolysis activity, and helicase-mediated unwinding activity. We found that the T407A and S411A substitutions in NS3 reduce viral replication and increase the helicase-unwinding turnover rates by 1.7- and 3.5-fold, respectively, suggesting that flaviviruses may use suboptimal NS3 helicase activity for optimal genome replication. Additionally, we used simulations of each mutant to probe structural changes within NS3 caused by each mutation. These simulations indicate that Motif V controls communication between the ATP-binding pocket and the helical gate. These results help define the linkage between ATP hydrolysis and helicase activities within NS3 and provide insight into the biophysical mechanisms for ATPase-driven NS3 helicase function.

NMR and MD simulations reveal the impact of the V23D mutation on the function of yeast oligosaccharyltransferase subunit Ost4.

Asparagine-linked glycosylation, also known as N-linked glycosylation, is an essential and highly conserved co- and post-translational protein modification in eukaryotes and some prokaryotes. In the central step of this reaction, a carbohydrate moiety is transferred from a lipid-linked donor to the side-chain of a consensus asparagine in a nascent protein as it is synthesized at the ribosome. Complete loss of oligosaccharyltransferase (OST) function is lethal in eukaryotes. This reaction is carried out by a membrane-associated multisubunit enzyme, OST, localized in the endoplasmic reticulum. The smallest subunit, Ost4, contains a single membrane-spanning helix that is critical for maintaining the stability and activity of OST. Mutation of any residue from Met18 to Ile24 of Ost4 destabilizes the enzyme complex, affecting its activity. Here, we report solution nuclear magnetic resonance structures and molecular dynamics (MD) simulations of Ost4 and Ost4V23D in micelles. Our studies revealed that while the point mutation did not impact the structure of the protein, it affected its position and solvent exposure in the membrane mimetic environment. Furthermore, our MD simulations of the membrane-bound OST complex containing either WT or V23D mutant demonstrated disruption of most hydrophobic helix-helix interactions between Ost4V23D and transmembrane TM12 and TM13 of Stt3. This disengagement of Ost4V23D from the OST complex led to solvent exposure of the D23 residue in the hydrophobic pocket created by these interactions. Our study not only solves the structures of yeast Ost4 subunit and its mutant but also provides a basis for the destabilization of the OST complex and reduced OST activity.

A Hyperactive Kunjin Virus NS3 Helicase Mutant Demonstrates Increased Dissemination and Mortality in Mosquitoes.

The unwinding of double-stranded RNA intermediates is critical for the replication and packaging of flavivirus RNA genomes. This unwinding activity is achieved by the ATP-dependent nonstructural protein 3 (NS3) helicase. In previous studies, we investigated the mechanism of energy transduction between the ATP and RNA binding pockets using molecular dynamics simulations and enzymatic characterization. Our data corroborated the hypothesis that motif V is a communication hub for this energy transduction. More specifically, mutations T407A and S411A in motif V exhibit a hyperactive helicase phenotype, leading to the regulation of translocation and unwinding during replication. However, the effect of these mutations on viral infection in cell culture and in vivo is not well understood. Here, we investigated the role of motif V in viral replication using West Nile virus (Kunjin subtype) T407A and S411A mutants (T407A and S411A Kunjin, respectively) in cell culture and in vivo We were able to recover S411A Kunjin but unable to recover T407A Kunjin. Our results indicated that S411A Kunjin decreased viral infection and increased cytopathogenicity in cell culture compared to wild-type (WT) Kunjin. Similarly, decreased infection rates in surviving S411A Kunjin-infected Culex quinquefasciatus mosquitoes were observed, but S411A Kunjin infection resulted in increased mortality compared to WT Kunjin infection. Additionally, S411A Kunjin infection increased viral dissemination and saliva positivity rates in surviving mosquitoes compared to WT Kunjin infection. These data suggest that S411A Kunjin increases viral pathogenesis in mosquitoes. Overall, these data indicate that NS3 motif V may play a role in the pathogenesis, dissemination, and transmission efficiency of Kunjin virus.IMPORTANCE Kunjin and West Nile viruses belong to the arthropod-borne flaviviruses, which can result in severe symptoms, including encephalitis, meningitis, and death. Flaviviruses have expanded into new populations and emerged as novel pathogens repeatedly in recent years, demonstrating that they remain a global threat. Currently, there are no approved antiviral therapeutics against either Kunjin or West Nile viruses. Thus, there is a pressing need for understanding the pathogenesis of these viruses in humans. In this study, we investigated the role of the Kunjin virus helicase on infection in cell culture and in vivo This work provides new insight into how flaviviruses control pathogenesis and mosquito transmission through the nonstructural protein 3 helicase.

Allostery in the dengue virus NS3 helicase: Insights into the NTPase cycle from molecular simulations.

The C-terminus domain of non-structural 3 (NS3) protein of the Flaviviridae viruses (e.g. HCV, dengue, West Nile, Zika) is a nucleotide triphosphatase (NTPase) -dependent superfamily 2 (SF2) helicase that unwinds double-stranded RNA while translocating along the nucleic polymer. Due to these functions, NS3 is an important target for antiviral development yet the biophysics of this enzyme are poorly understood. Microsecond-long molecular dynamic simulations of the dengue NS3 helicase domain are reported from which allosteric effects of RNA and NTPase substrates are observed. The presence of a bound single-stranded RNA catalytically enhances the phosphate hydrolysis reaction by affecting the dynamics and positioning of waters within the hydrolysis active site. Coupled with results from the simulations, electronic structure calculations of the reaction are used to quantify this enhancement to be a 150-fold increase, in qualitative agreement with the experimental enhancement factor of 10-100. Additionally, protein-RNA interactions exhibit NTPase substrate-induced allostery, where the presence of a nucleotide (e.g. ATP or ADP) structurally perturbs residues in direct contact with the phosphodiester backbone of the RNA. Residue-residue network analyses highlight pathways of short ranged interactions that connect the two active sites. These analyses identify motif V as a highly connected region of protein structure through which energy released from either active site is hypothesized to move, thereby inducing the observed allosteric effects. These results lay the foundation for the design of novel allosteric inhibitors of NS3.

Self-Assembly of Perylenediimide-Single-Strand-DNA Conjugates: Employing Hydrophobic Interactions and DNA Base-Pairing To Create a Diverse Structural Space.

The self-assembly behavior of DNA conjugates possessing a perylenediimide (PDI) head group and an N-oligonucleotide tail has been investigated using a combination of optical spectroscopy and cryogenic transmission electron microscopy (cryo-TEM) imaging. To obtain insight into the interplay between PDI hydrophobic interactions and DNA base-pairing we employed systematic variation in the length and composition of the oligo tails. Conjugates with short (TA)n or (CG)n oligo tails (n≤3) form helical or nonhelical fibers constructed from π-stacked PDI head groups with pendent oligo tails in aqueous solution. Conjugates with longer (TA)n oligo tails also form stacks of PDI head groups, which are further aggregated by base-pairing between their oligo tails, leading to fiber bundling and formation of bilayers. The longer (CG)n conjugates form PDI end-capped duplexes, which further assemble into PDI-stacked arrays of duplexes leading to large scale ordered assemblies. Cryo-TEM imaging reveals that (CG)3 gives rise to both fibers and large assemblies, whereas (CG)5 assembles preferentially into large ordered structures.

Unraveling the mystery of ATP hydrolysis in actin filaments.

Actin performs its myriad cellular functions by the growth and disassembly of its filamentous form. The hydrolysis of ATP in the actin filament has been shown to modulate properties of the filament, thus making it a pivotal regulator of the actin life cycle. Actin has evolved to selectively hydrolyze ATP in the filamentous form, F-actin, with an experimentally observed rate increase over the monomeric form, G-actin, of 4.3 × 10(4). The cause of this dramatic increase in rate is investigated in this paper using extensive QM/MM simulations of both G- and F-actin. To compute the free energy of hydrolysis in both systems, metadynamics is employed along two collective variables chosen to describe the reaction coordinates of hydrolysis. F-actin is modeled as a monomer with restraints applied to coarse-grained variables enforced to keep it in a filament-like conformation. The simulations reveal a barrier height reduction for ATP hydrolysis in F-actin as compared to G-actin of 8 ± 1 kcal/mol, in good agreement with the experimentally measured barrier height reduction of 7 ± 1 kcal/mol. The barrier height reduction is influenced by an enhanced rotational diffusion of water in F-actin as compared to G-actin and shorter water wires between Asp154 and the nucleophilic water in F-actin, leading to more rapid proton transport.

Thymine photodimer formation in DNA hairpins. Unusual conformations favor (6 - 4) vs. (2 + 2) adducts.

The photochemical reactions of eleven synthetic DNA hairpins possessing a single TT step either in a base-paired stem or in a hexanucleotide linker have been investigated. The major reaction products have been identified as the cis-syn (2 + 2) adduct and the (6 - 4) adduct on the basis of their spectroscopic properties including 1D and 2D NMR spectra, UV spectra and stability or instability to photochemical cleavage. Product quantum yields and ratios determined by HPLC analysis allow the behaviour of the eleven hairpins to be placed into three groups: Group I in which the (2 + 2) adduct is the major product, as is usually the case for DNA, Group II in which comparable amounts of (2 + 2) and (6 - 4) adducts are formed, and Group III in which the major product is the (6 - 4) adduct. The latter behaviour is without precedent in natural or synthetic DNA and appears to be related to the highly fluxional structures of the hairpin reactants. Molecular dynamics simulation of ground state conformations provides quantum yields and product ratios calculated using a single parameter model that are in reasonable agreement with most of the experimental results. Factors which may influence the observed product ratios are discussed.

Role of ATP Hydrolysis and Product Release in the Translocation Mechanism of SARS-CoV-2 NSP13.

In response to the emergence of COVID-19, caused by SARS-CoV-2, there has been a growing interest in understanding the functional mechanisms of the viral proteins to aid in the development of new therapeutics. Nonstructural protein 13 (nsp13) helicase is an attractive target for antivirals because it is essential for viral replication and has a low mutation rate, yet the structural mechanisms by which this enzyme binds and hydrolyzes ATP to cause unidirectional RNA translocation remain elusive. Using Gaussian accelerated molecular dynamics (GaMD), we generated comprehensive conformational ensembles of all substrate states along the ATP-dependent cycle. Shape-GMM clustering of the protein yields four protein conformations that describe an opening and closing of both the ATP pocket and the RNA cleft that is achieved through a combination of conformational selection and induction along the ATP hydrolysis cycle. Furthermore, three protein-RNA conformations are observed that implicate motifs Ia, IV, and V as playing a pivotal role in an ATP-dependent inchworm translocation mechanism. Finally, based on a linear discriminant analysis of protein conformations, we identify L405 as a pivotal residue for the opening and closing mechanism and propose a L405D mutation as a way to disrupt translocation. This research enhances our understanding of nsp13's role in viral replication and could contribute to the development of antiviral strategies.

Quantifying Unbiased Conformational Ensembles from Biased Simulations Using ShapeGMM.

Quantifying the conformational ensembles of biomolecules is fundamental to describing mechanisms of processes such as protein folding, interconversion between folded states, ligand binding, and allosteric regulation. Accurate quantification of these ensembles remains a challenge for conventional molecular simulations of all but the simplest molecules due to insufficient sampling. Enhanced sampling approaches, such as metadynamics, were designed to overcome this challenge; however, the nonuniform frame weights that result from many of these approaches present an additional challenge to ensemble quantification techniques such as Markov State Modeling or structural clustering. Here, we present rigorous inclusion of nonuniform frame weights into a structural clustering method entitled shapeGMM. The result of frame-weighted shapeGMM is a high dimensional probability density and generative model for the unbiased system from which we can compute important thermodynamic properties such as relative free energies and configurational entropy. The accuracy of this approach is demonstrated by the quantitative agreement between GMMs computed by Hamiltonian reweighting and direct simulation of a coarse-grained helix model system. Furthermore, the relative free energy computed from a shapeGMM probability density of alanine dipeptide reweighted from a metadynamics simulation quantitatively reproduces the underlying free energy in the basins. Finally, the method identifies hidden structures along the actin globular to filamentous-like structural transition from a metadynamics simulation on a linear discriminant analysis coordinate trained on GMM states, illustrating how structural clustering of biased data can lead to biophysical insight. Combined, these results demonstrate that frame-weighted shapeGMM is a powerful approach to quantifying biomolecular ensembles from biased simulations.

Reaction Coordinates for Conformational Transitions Using Linear Discriminant Analysis on Positions.

In this work, we demonstrate that Linear Discriminant Analysis (LDA) applied to atomic positions in two different states of a biomolecule produces a good reaction coordinate between those two states. Atomic coordinates of a macromolecule are a direct representation of a macromolecular configuration, and yet, they are not used in enhanced sampling studies due to a lack of rotational and translational invariance. We resolve this issue using the technique of our prior work, whereby a molecular configuration is considered a member of an equivalence class in size-and-shape space, which is the set of all configurations that can be translated and rotated to a single point within a reference multivariate Gaussian distribution characterizing a single molecular state. The reaction coordinates produced by LDA applied to positions are shown to be good reaction coordinates both in terms of characterizing the transition between two states of a system within a long molecular dynamics (MD) simulation and also ones that allow us to readily produce free energy estimates along that reaction coordinate using enhanced sampling MD techniques.

The Role of ATP Hydrolysis and Product Release in the Translocation Mechanism of SARS-CoV-2 NSP13.

In response to the emergence of COVID-19, caused by SARS-CoV-2, there has been a growing interest in understanding the functional mechanisms of the viral proteins to aid in the development of new therapeutics. Non-structural protein 13 (Nsp13) helicase is an attractive target for antivirals because it is essential for viral replication and has a low mutation rate; yet, the structural mechanisms by which this enzyme binds and hydrolyzes ATP to cause unidirectional RNA translocation remain elusive. Using Gaussian accelerated molecular dynamics (GaMD), we generated a comprehensive conformational ensemble of all substrate states along the ATP-dependent cycle. ShapeGMM clustering of the protein yields four protein conformations that describe an opening and closing of both the ATP pocket and RNA cleft. This opening and closing is achieved through a combination of conformational selection and induction along the ATP cycle. Furthermore, three protein-RNA conformations are observed that implicate motifs Ia, IV, and V as playing a pivotal role in an ATP-dependent inchworm translocation mechanism. Finally, based on a linear discriminant analysis of protein conformations, we identify L405 as a pivotal residue for the opening and closing mechanism and propose a L405D mutation as a way of testing our proposed mechanism. This research enhances our understanding of nsp13's role in viral replication and could contribute to the development of antiviral strategies.

Motif-VI Loop Acts as a Nucleotide Valve in the West Nile Virus NS3 Helicase.

The flavivirus NS3 helicase (NS3h), a highly conserved protein, plays a pivotal role in virus replication and thus represents a potential drug target for flavivirus pathogenesis. NS3h utilizes nucleotide triphosphate, such as ATP, for hydrolysis energy (ATPase) to translocate on single-stranded nucleic acids, which is an important step in the unwinding of double-stranded nucleic acids. The intermediate states along the ATP binding and hydrolysis cycle, as well as the conformational changes between these states, represent important yet difficult-to-identify targets for potential inhibitors. We use extensive molecular dynamics simulations of apo, ATP, ADP+Pi, and ADP bound to WNV NS3h+ssRNA to model the conformational ensembles along this cycle. Energetic and structural clustering analyses on these trajectories depict a clear trend of differential enthalpic affinity of NS3h with ADP, demonstrating a probable mechanism of hydrolysis turnover regulated by the motif-VI loop (MVIL). These findings were experimentally corroborated using viral replicons encoding three mutations at the D471 position. Replication assays using these mutants demonstrated a substantial reduction in viral replication compared to the wild-type. Molecular simulations of the D471 mutants in the apo state indicate a shift in MVIL populations favoring either a closed or open 'valve' conformation, affecting ATP entry or stabilization, respectively. Combining our molecular modeling with experimental evidence highlights a conformation-dependent role for MVIL as a 'valve' for the ATP-pocket, presenting a promising target for antiviral development.

Size-and-Shape Space Gaussian Mixture Models for Structural Clustering of Molecular Dynamics Trajectories.

Determining the optimal number and identity of structural clusters from an ensemble of molecular configurations continues to be a challenge. Recent structural clustering methods have focused on the use of internal coordinates due to the innate rotational and translational invariance of these features. The vast number of possible internal coordinates necessitates a feature space supervision step to make clustering tractable but yields a protocol that can be system type-specific. Particle positions offer an appealing alternative to internal coordinates but suffer from a lack of rotational and translational invariance, as well as a perceived insensitivity to regions of structural dissimilarity. Here, we present a method, denoted shape-GMM, that overcomes the shortcomings of particle positions using a weighted maximum likelihood alignment procedure. This alignment strategy is then built into an expectation maximization Gaussian mixture model (GMM) procedure to capture metastable states in the free-energy landscape. The resulting algorithm distinguishes between a variety of different structures, including those indistinguishable by root-mean-square displacement and pairwise distances, as demonstrated on several model systems. Shape-GMM results on an extensive simulation of the fast-folding HP35 Nle/Nle mutant protein support a four-state folding/unfolding mechanism, which is consistent with previous experimental results and provides kinetic details comparable to previous state-of-the art clustering approaches, as measured by the VAMP-2 score. Currently, training of shape-GMMs is recommended for systems (or subsystems) that can be represented by ≲200 particles and ≲100k configurations to estimate high-dimensional covariance matrices and balance computational expense. Once a shape-GMM is trained, it can be used to predict the cluster identities of millions of configurations.

Conserved motifs in the flavivirus NS3 RNA helicase enzyme.

Flaviviruses are a major health concern because over half of the world population is at risk of infection and there are very few antiviral therapeutics to treat diseases resulting from infection. Replication is an essential part of the flavivirus survival. One of the viral proteins, NS3 helicase, is critical for unwinding the double stranded RNA intermediate during flaviviral replication. The helicase performs the unwinding of the viral RNA intermediate structure in an ATP-dependent manner. NS3 helicase is a member of the Viral/DEAH-like subfamily of the superfamily 2 helicase containing eight highly conserved structural motifs (I, Ia, II, III, IV, IVa, V, and VI) localized between the ATP-binding and RNA-binding pockets. Of these structural motifs only three are well characterized for function in flaviviruses (I, II, and VI). The roles of the other structural motifs are not well understood for NS3 helicase function, but comparison of NS3 with other superfamily 2 helicases within the viral/DEAH-like, DEAH/RHA, and DEAD-box subfamilies can be used to elucidate the roles of these structural motifs in the flavivirus NS3 helicase. This review aims to summarize the role of each conserved structural motif within flavivirus NS3 in RNA helicase function. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA in Disease and Development > RNA in Disease.

Modeling Catalysis in Allosteric Enzymes: Capturing Conformational Consequences.

Greater understanding of enzymatic mechanisms aids the discovery of new targets for biologics, the development of biocatalytic transformations, and de novo enzyme design. Methods using quantum mechanical (QM) potentials, such as Density Functional Theory (DFT), have enabled complex multistep enzymatic mechanisms to be studied, often in quantitative detail. Nevertheless, the dynamic interconversion of enzyme conformations between active and inactive catalytic forms, involving length- and timescales inaccessible to QM treatments, presents a formidable challenge for the development of computational models for allosterically modulated enzymes. We present an overview of the key concepts underlying multistate models of enzyme catalysis, enzyme allostery, and the challenge that large-scale conformational changes pose for methods using QM, QM/MM, and MM potentials. Structural clustering is highlighted as a valuable approach to bridge molecular dynamics conformational sampling of MM potentials and quantum chemical cluster models of catalysis. Particularly relevant to this discussion is structural allostery, which serves as the exemplar of conformational consequences. Here, a well-characterized allosteric enzyme, Imidazole Glycerol Phosphate Synthase (IGPS), is used to showcase the importance of multiple conformations and guide a new direction for qualitative understanding and quantitative modeling in enzyme catalysis.

DNA-based optomechanical molecular motor.

An azobenzene-capped DNA hairpin coupled to an AFM is presented as an optically triggered single-molecule motor. The photoinduced trans to cis isomerization of azobenzene affects both the overall length of the molecule and the ability of the DNA bases to hybridize. Using a combination of molecular dynamics simulations and free energy calculations the unfolding of both isomers along the O5'-O3' extension coordinate is monitored. The potentials of mean force (PMFs) along this coordinate indicate that there are two major differences induced by photoisomerization. The first is that the interbase hydrogen bond and stacking interactions are stable for a greater range of extensions in the trans system than in the cis system. The second difference is due to a decreased chain length of the cis isomer with respect to the trans isomer. These differences are exploited to extract work in optomechanical cycles. The disruption of the hairpin structure gives a maximum of 3.4 kcal mol(-1) of extractable work per cycle with an estimated maximum efficiency of 2.4%. Structure-function insights into the operation of this motor are provided, and the effect of the cantilever stiffness on the extractable work is characterized.

Electron donor-acceptor interactions with flanking purines influence the efficiency of thymine photodimerization.

Quantum yields for thymine photodimerization (Φ(TT)) have been determined for a series of short DNA single-strand and base-paired hairpin structures possessing a single thymine-thymine step with flanking purines. Values of Φ(TT) are strongly dependent upon the oxidation potential of the flanking purine, decreasing in the order: inosine > adenine > guanine > deazaguanine. The dependence of Φ(TT) on the ionization potential of the flanking purine is more pronounced when the purine of lower oxidation potential is located at the 5'- versus 3'-position in either a single strand or a hairpin. Molecular dynamics simulations for hairpin structures indicate that the TT step is π-stacked with both the 5' and 3' purine, but that there is little π-stacking with either purine in single-strand structures. The observation of moderately intense long-wavelength UV absorption features for hairpins having 5'-Z or G flanking purines suggests that excitation of ground state donor-acceptor complexes may account for more extensive quenching of dimerization by 5'- versus 3'-purines. The "purine effect" on Φ(TT) is attributed to a combination of ground state conformation, ground state electron donor-acceptor interactions, and excited state exciplex formation.

Role of ATP in the RNA Translocation Mechanism of SARS-CoV-2 NSP13 Helicase.

The COVID-19 pandemic has demonstrated the need to develop potent and transferable therapeutics to treat coronavirus infections. Numerous antiviral targets are being investigated, but nonstructural protein 13 (nsp13) stands out as a highly conserved and yet understudied target. Nsp13 is a superfamily 1 (SF1) helicase that translocates along and unwinds viral RNA in an ATP-dependent manner. Currently, there are no available structures of nsp13 from SARS-CoV-1 or SARS-CoV-2 with either ATP or RNA bound, which presents a significant hurdle to the rational design of therapeutics. To address this knowledge gap, we have built models of SARS-CoV-2 nsp13 in Apo, ATP, ssRNA and ssRNA+ATP substrate states. Using 30 µs of a Gaussian-accelerated molecular dynamics simulation (at least 6 µs per substrate state), these models were confirmed to maintain substrate binding poses that are similar to other SF1 helicases. A Gaussian mixture model and linear discriminant analysis structural clustering protocol was used to identify key structural states of the ATP-dependent RNA translocation mechanism. Namely, four RNA-nsp13 structures are identified that exhibit ATP-dependent populations and support the inchworm mechanism for translocation. These four states are characterized by different RNA-binding poses for motifs Ia, IV, and V and suggest a power stroke-like motion of domain 2A relative to domain 1A. This structural and mechanistic insight of nsp13 RNA translocation presents novel targets for the further development of antivirals.