Design, high-level expression, purification and characterization of soluble fragments of the hepatitis C virus NS3 RNA helicase suitable for NMR-based drug discovery methods and mechanistic studies

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.

Non-structural protein 3 (NS3) is a Flaviviridae (e.g. Hepatitis C, dengue, and Zika viruses) helicase that unwinds double stranded RNA while translocating along the nucleic polymer during viral genome replication. As a member of superfamily 2 (SF2) helicases, NS3 utilizes the free energy of nucleotide triphosphate (NTP) binding, hydrolysis, and product unbinding to perform its functions. While much is known about SF2 helicases, the pathways and mechanisms through which free energy is transduced between the NTP hydrolysis active site and RNA binding cleft remains elusive. Here we present a multiscale computational study to characterize the allosteric effects induced by the RNA and NTPase substrates (ATP, ADP, and P_i) as well as the pathways of short-range, residue-residue interactions that connect the two active sites. Results from this body of molecular dynamics simulations and electronic structure calculations are highlighted in context to the NTPase enzymatic cycle, allowing for development of testable hypotheses for validation of these simulations. Our insights, therefore, provide novel details about the biophysics of NS3 and guide the next generation of experimental studies.

A quantum mechanics/molecular mechanics study of the reaction mechanism of the hepatitis C virus NS3 protease with the NS5A/5B substrate

Combined quantum mechanics and molecular mechanics (QM/MM) calculations were carried out to characterize the reaction mechanism of the NS3 protease with its preferred substrate (NS5A/5B). The main purpose of this study was to locate the barrier states and intermediates along the distinguished coordinate path (DCP) involved in this process. These structures, and in particular the one corresponding to the first barrier state and intermediate (B1 and I1), could be a starting point for the synthesis of inhibitors of this protease, which could be used to treat hepatitis C. The two first steps of the reaction mechanism were studied, i.e., the acylation step and the breaking of the peptide bond. The first step takes place through a tetracoordinated intermediate, as suggested from previous works on other Serine proteases. The importance of the different amino acid residues was also considered (perturbation study where the MM charges of each residue were set to zero independently). The residues of the oxyanion hole were confirmed as the most important for the electrostatic stabilization of the tetracoordinate intermediate. Moreover, the role of other residues, e.g., Arg-155 and Asp-79, was also explained.

Double-stranded DNA-induced localized unfolding of HCV NS3 helicase subdomain 2

The NS3 helicase of the hepatitis C virus (HCV) unwinds double-stranded (ds) nucleic acid (NA) in an NTP-dependent fashion. Mechanistic details of this process are, however, largely unknown for the HCV helicase. We have studied the binding of dsDNA to an engineered version of subdomain 2 of the HCV helicase (d(2Delta)NS3h) by NMR and circular dichroism. Binding of dsDNA to d(2Delta)NS3h induces a local unfolding of helix (alpha(3)), which includes residues of conserved helicase motif VI (Q(460)RxxRxxR(467)), and strands (beta(1) and beta(8)) from the central beta-sheet. This also occurs upon lowering the pH (4.4) and introducing an R461A point mutation, which disrupt salt bridges with Asp 412 and Asp 427 in the protein structure. NMR studies on d(2Delta)NS3h in the partially unfolded state at low pH map the dsDNA binding site to residues previously shown to be involved in single-stranded DNA binding. Sequence alignment and structural comparison suggest that these Arg-Asp interactions are highly conserved in SF2 DEx(D/H) proteins. Thus, modulation of these interactions by dsNA may allow SF2 helicases to switch between conformations required for helicase function.

Two novel conserved motifs in the hepatitis C virus NS3 protein critical for helicase action

The hepatitis C virus (HCV) NS3 helicase shares several conserved motifs with other superfamily 2 (SF2) helicases. Besides these sequences, several additional helicase motifs are conserved among the various HCV genotypes and quasispecies. The roles of two such motifs are examined here. The first motif (YRGXDV) forms a loop that connects SF2 helicase motifs 4 and 5, at the tip of which is Arg393. When Arg393 is changed to Ala, the resulting protein (R393A) retains a nucleic acid stimulated ATPase but cannot unwind RNA. R393A also unwinds DNA more slowly than wild type and translocates poorly on single-stranded DNA (ssDNA). DNA and RNA stimulate ATP hydrolysis catalyzed by R393A like the wild type, but the mutant protein binds ssDNA more weakly both in the presence and absence of the non-hydrolyzable ATP analog ADP(BeF3). The second motif (DFSLDPTF) forms a loop that connects two anti-parallel sheets between SF2 motifs 5 and 6. When Phe444 in this Phe-loop is changed to Ala, the resulting protein (F444A) is devoid of all activities. When Phe438 is changed to Ala, the protein (F438A) retains nucleic acid-stimulated ATPase, but does not unwind RNA. F438A unwinds DNA and translocates on ssDNA at about half the rate of the wild type. Equilibrium binding data reveal that this uncoupling of ATP hydrolysis and unwinding is due to the fact that the F438A mutant does not release ssDNA upon ATP binding like the wild type. A model is presented explaining the roles of the Arg-clamp and the Phe-loop in the unwinding reaction.

Solubility engineering of the HhaI methyltransferase

DNA methylation is involved in epigenetic control of numerous cellular processes in eukaryotes, however, many mechanistic aspects of this phenomenon are not yet understood. A bacterial prototype cytosine-C5 methyltransferase, M.HhaI, serves as a paradigm system for structural and mechanistic studies of biological DNA methylation, but further analysis of the 37 kDa protein is hampered by its insufficient solubility (0.15 mM). To overcome this problem, three hydrophobic patches on the surface of M.HhaI that are not involved in substrate interactions were subjected to site-specific mutagenesis. Residues M51 or V213 were substituted by polar amino acids of a similar size, and/or the C-terminal tetrapeptide FKPY was replaced by a single glycine residue (Delta324G). Two out of six mutants, delta324G and V213S/delta324G, showed improved solubility in initial analyses and were purified to homogeneity using a newly developed procedure. Biochemical studies of the engineered methyltransferases showed that the deletion mutant delta324G retained identical DNA binding, base flipping and catalytic properties as the wild-type enzyme. In contrast, the engineered enzyme showed (i) a significantly increased solubility (>0.35 mM), (ii) high-quality 2D-[(15)N,(1)H] TROSY NMR spectra, and (iii) (15)N spin relaxation times evidencing the presence of a monomeric well-folded protein in solution.

Solution structure and backbone dynamics of an engineered arginine-rich subdomain 2 of the hepatitis C virus NS3 RNA helicase

The NS3 protein of the hepatitis C virus (HCV) is a 631 amino acid residue bifunctional enzyme with a serine protease localized to the N-terminal 181 residues and an RNA helicase located in the C-terminal 450 residues. The HCV NS3 RNA helicase consists of three well-defined subdomains which all contribute to its helicase activity. The second subdomain of the HCV helicase is flexibly linked to the remainder of the NS3 protein and could undergo rigid-body movements during the unwinding of double-stranded RNA. It also contains several motifs that are implicated in RNA binding and in coupling NTP hydrolysis to nucleic acid unwinding and translocation. As part of our efforts to use NMR techniques to assist in deciphering the enzyme's structure-function relationships and developing specific small molecule inhibitors, we have determined the solution structure of an engineered subdomain 2 of the NS3 RNA helicase of HCV, d(2Delta)-HCVh, and studied the backbone dynamics of this protein by (15)N-relaxation experiments using a model-free approach. The NMR studies on this 142-residue construct reveal that overall subdomain 2 of the HCV helicase is globular and well structured in solution even in the absence of the remaining parts of the NS3 protein. Its solution structure is very similar to the corresponding parts in the X-ray structures of the HCV NS3 helicase domain and intact bifunctional HCV NS3 protein. Slow hydrogen-deuterium exchange rates map to a well-structured, stable hydrophobic core region away from the subdomain interfaces. In contrast, the regions facing the subdomain interfaces in the HCV NS3 helicase domain are less well structured in d(2Delta)-HCVh, show fast hydrogen-deuterium exchange rates, and the analysis of the dynamic properties of d(2Delta)-HCVh reveals that these regions of the protein show distinct dynamical features. In particular, residues in motif V, which may be involved in transducing allosteric effects of nucleotide binding and hydrolysis on RNA binding, exhibit slow conformational exchange on the milli- to microsecond time-scale. The intrinsic conformational flexibility of this loop region may facilitate conformational changes required for helicase function.