Allosteric effects of ligands and mutations on poliovirus RNA-dependent RNA polymerase

Potential Inhibition of COVID-19 RNA-dependent RNA Polymerase by Hepatitis C Virus Non-nucleoside Inhibitors: An In-silico Perspective

Background:

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a novel member of the genus betacoronavirus in the Coronaviridae family. It has been identified as the causative agent of coronavirus disease 2019 (COVID-19), spreading rapidly in Asia, America and Europe. Like some other RNA viruses, RNA replication and transcription of SARS-CoV-2 rely on its RNA-dependent RNA polymerase (RdRP), which is a therapeutic target of clinical importance. Crystal structure of SARS-CoV-2 was solved recently (PDB ID 6M71) with some missing residues.

Objective:

We used SARS-CoV-2 RdRP as a target protein to screen for possible chemical molecules with potential anti-viral effects.

Methods:

Here we modelled the missing residues 896-905 via homology modelling and then analysed the interactions of Hepatitis C virus allosteric non-nucleoside inhibitors (NNIs) in the reported NNIs binding sites in SARS-CoV-2 RdRP.

Results:

We found that MK-3281, filibuvir, setrobuvir and dasabuvir might be able to inhibit SARS-CoV-2 RdRP based on their binding affinities in the respective binding sites.

Conclusion:

Further in vitro and in vivo experimental research will be carried out to evaluate their effectiveness in COVID-19 treatment in the near future.

Allosteric inhibitors of Coxsackie virus A24 RNA polymerase

Coxsackie virus A24 (CVA24), a causative agent of acute hemorrhagic conjunctivitis, is a prototype of enterovirus (EV) species C. The RNA polymerase (3D(pol)) of CVA24 can uridylylate the viral peptide linked to the genome (VPg) from distantly related EV and is thus, a good model for studying this reaction. Once UMP is bound, VPgpU primes RNA elongation. Structural and mutation data have identified a conserved binding surface for VPg on the RNA polymerase (3D(pol)), located about 20Å from the active site. Here, computational docking of over 60,000 small compounds was used to select those with the lowest (best) specific binding energies (BE) for this allosteric site. Compounds with varying structures and low BE were assayed for their effect on formation of VPgU by CVA24-3D(pol). Two compounds with the lowest specific BE for the site inhibited both uridylylation and formation of VPgpolyU at 10-20μM. These small molecules can be used to probe the role of this allosteric site in polymerase function, and may be the basis for novel antiviral compounds.

Amiloride inhibits the initiation of Coxsackievirus and poliovirus RNA replication by inhibiting VPg uridylylation

The mechanism of amiloride inhibition of Coxsackievirus B3 (CVB3) and poliovirus type 1 (PV1) RNA replication was investigated using membrane-associated RNA replication complexes. Amiloride was shown to inhibit viral RNA replication and VPgpUpU synthesis. However, the drug had no effect on polymerase elongation activity during either (-) strand or (+) strand synthesis. These findings indicated that amiloride inhibited the initiation of RNA synthesis by inhibiting VPg uridylylation. In addition, in silico binding studies showed that amiloride docks in the VPg binding site on the back of the viral RNA polymerase, 3D(pol). Since VPg binding at this site on PV1 3D(pol) was previously shown to be required for VPg uridylylation, our results suggest that amiloride inhibits VPg binding to 3D(pol). In summary, our findings are consistent with a model in which amiloride inhibits VPgpUpU synthesis and viral RNA replication by competing with VPg for binding to 3D(pol).

Crystal structure of enterovirus 71 RNA-dependent RNA polymerase complexed with its protein primer VPg: implication for a trans mechanism of VPg uridylylation

Picornavirus RNA replication is initiated by VPg uridylylation, during which the hydroxyl group of the third tyrosine residue of the virally encoded protein VPg is covalently linked to two UMP molecules by RNA-dependent RNA polymerase (RdRp; also known as 3D(pol)). We previously identified site 311, located at the base of the palm domain of the enterovirus 71 (EV71) RdRp, to be the site for EV71 VPg binding and uridylylation. Here we report the crystal structure of EV71 3D(pol) complexed with VPg. VPg was anchored at the bottom of the palm domain of the 3D(pol) molecule and exhibited an extended V-shape conformation. The corresponding interface on 3D(pol) was mainly formed by residues within site 311 and other residues in the palm and finger domains. Mutations of the amino acids of 3D(pol) involved in the VPg interaction (3DL319A, 3DD320A, and 3DY335A) significantly disrupted VPg binding to 3D(pol), resulting in defective VPg uridylylation. In contrast, these mutations did not affect the RNA elongation activity of 3D(pol). In the context of viral genomic RNA, mutations that abolished VPg uridylylation activity were lethal for EV71 replication. Further in vitro analysis showed that the uridylylation activity was restored by mixing VPg-binding-defective and catalysis-defective mutants, indicating a trans mechanism for EV71 VPg uridylylation. Our results, together with previous results of other studies, demonstrate that different picornaviruses use distinct binding sites for VPg uridylylation.

Surface for catalysis by poliovirus RNA-dependent RNA polymerase

The poliovirus RNA-dependent RNA polymerase, 3Dpol, replicates the viral genomic RNA on the surface of virus-induced intracellular membranes. Macromolecular assemblies of 3Dpol form linear arrays of subunits that propagate along a strong protein-protein interaction called interface-I, as was observed in the crystal structure of wild-type poliovirus polymerase. These "filaments" recur with slight modifications in planar sheets and, with additional modifications that accommodate curvature, in helical tubes of the polymerase, by packing filaments together via a second set of interactions. Periodic variations of subunit orientations within 3Dpol tubes give rise to "ghost reflections" in diffraction patterns computed from electron cryomicrographs of helical arrays. The ghost reflections reveal that polymerase tubes are formed by bundles of four to five interface-I filaments, which are then connected to the next bundle of filaments with a perturbation of interface interactions between bundles. While enzymatically inactive polymerase is also capable of oligomerization, much thinner tubes that lack interface-I interactions between adjacent subunits are formed, suggesting that long-range allostery produces conformational changes that extend from the active site to the protein-protein interface. Macromolecular assemblies of poliovirus polymerase show repeated use of flexible interface interactions for polymerase lattice formation, suggesting that adaptability of polymerase-polymerase interactions facilitates RNA replication. In addition, the presence of a positively charged groove identified in polymerase arrays may help position and stabilize the RNA template during replication.

Interstitial contacts in an RNA-dependent RNA polymerase lattice

Catalytic activities can be facilitated by ordered enzymatic arrays that co-localize and orient enzymes and their substrates. The purified RNA-dependent RNA polymerase from poliovirus self-assembles to form two-dimensional lattices, possibly facilitating the assembly of viral RNA replication complexes on the cytoplasmic face of intracellular membranes. Creation of a two-dimensional lattice requires at least two different molecular contacts between polymerase molecules. One set of polymerase contacts, between the "thumb" domain of one polymerase and the back of the "palm" domain of another, has been previously defined. To identify the second interface needed for lattice formation and to test its function in viral RNA synthesis, we used a hybrid approach of electron microscopic and biochemical evaluation of both wild-type and mutant viral polymerases to evaluate computationally generated models of this second interface. A unique solution satisfied all constraints and predicted a two-dimensional structure formed from antiparallel arrays of polymerase fibers that use contacts from the flexible amino-terminal region of the protein. Enzymes that contained mutations in this newly defined interface did not form lattices and altered the structure of wild-type lattices. When reconstructed into virus, mutations that disrupt lattice assembly exhibited growth defects, synthetic lethality or both, supporting the function of the oligomeric lattice in infected cells. Understanding the structure of polymerase lattices within the multimeric RNA-dependent RNA polymerase complex should facilitate antiviral drug design and provide a precedent for other positive-strand RNA viruses.

Long-range interaction networks in the function and fidelity of poliovirus RNA-dependent RNA polymerase studied by nuclear magnetic resonance

The fidelity of the poliovirus RNA-dependent RNA polymerase (3D(pol)) plays a direct role in the genomic evolution and pathogenesis of the virus. A single site mutation (Gly64Ser) that is remote from the catalytic center results in a higher fidelity polymerase. NMR studies with [methyl-(13)C]methionine-labeled protein were used to compare the solution structure and dynamics of wild-type and Gly64Ser 3D(pol). The chemical shifts for the Met6 resonance were significantly different between wild-type and Gly64Ser 3D(pol) when bound in ternary complexes with RNA and incorrect, but not with correct, nucleotide, suggesting that the Gly64Ser mutation induces structural changes in the N-terminal β-strand when the enzyme is bound to incorrect but not correct nucleotide. We also observe changes in the transverse relaxation times for methionines near regions important for nucleotide and RNA binding and catalysis. Our strategy to assign the [methyl-(13)C]methionine resonances involved separately mutating each of the 17 methionines. Several substitutions produced additional resonances for both Met6 and Met187, a reporter for RNA binding, and conformational changes in the highly conserved motif B loop, even though these methionines are greater than 20 Å apart. The results for Gly64Ser and the other mutants are intriguing considering that they can result in structural and/or dynamic changes to methionines distant from the site of mutation. We propose that there is a long-distance network operating throughout 3D(pol) that coordinates ligand binding, conformational changes, and catalysis. Mutation of Gly64 results in structural and/or dynamic changes to the network that may affect polymerase fidelity.

Enzymatic and nonenzymatic functions of viral RNA-dependent RNA polymerases within oligomeric arrays

Few antivirals are effective against positive-strand RNA viruses, primarily because the high error rate during replication of these viruses leads to the rapid development of drug resistance. One of the favored current targets for the development of antiviral compounds is the active site of viral RNA-dependent RNA polymerases. However, like many subcellular processes, replication of the genomes of all positive-strand RNA viruses occurs in highly oligomeric complexes on the cytosolic surfaces of the intracellular membranes of infected host cells. In this study, catalytically inactive polymerases were shown to participate productively in functional oligomer formation and catalysis, as assayed by RNA template elongation. Direct protein transduction to introduce either active or inactive polymerases into cells infected with mutant virus confirmed the structural role for polymerase molecules during infection. Therefore, we suggest that targeting the active sites of polymerase molecules is not likely to be the best antiviral strategy, as inactivated polymerases do not inhibit replication of other viruses in the same cell and can, in fact, be useful in RNA replication complexes. On the other hand, polymerases that could not participate in functional RNA replication complexes were those that contained mutations in the amino terminus, leading to altered contacts in the folded polymerase and mutations in a known polymerase-polymerase interaction in the two-dimensional protein lattice. Thus, the functional nature of multimeric arrays of RNA-dependent RNA polymerase supplies a novel target for antiviral compounds and provides a new appreciation for enzymatic catalysis on membranous surfaces within cells.

Picornavirus genome replication: roles of precursor proteins and rate-limiting steps in oriI-dependent VPg uridylylation

The 5' ends of all picornaviral RNAs are linked covalently to the genome-encoded peptide, VPg (or 3B). VPg linkage is thought to occur in two steps. First, VPg serves as a primer for production of diuridylylated VPg (VPg-pUpU) in a reaction catalyzed by the viral polymerase that is templated by an RNA element (oriI). It is currently thought that the viral 3AB protein is the source of VPg in vivo. Second, VPg-pUpU is transferred to the 3' end of plus- and/or minus-strand RNA and serves as primer for production of full-length RNA. Nothing is known about the mechanism of transfer. We present biochemical and biological evidence refuting the use of 3AB as the donor for VPg uridylylation. Our data are consistent with precursors 3BC and/or 3BCD being employed for uridylylation. This conclusion is supported by in vitro uridylylation of these proteins, the ability of a mutant replicon incapable of producing processed VPg to replicate in HeLa cells and cell-free extracts and corresponding precursor processing profiles, and the demonstration of 3BC-linked RNA in mutant replicon-transfected cells. These data permit elaboration of our model for VPg uridylylation to include the use of precursor proteins and invoke a possible mechanism for location of the diuridylylated, VPg-containing precursor at the 3' end of plus- or minus-strand RNA for production of full-length RNA. Finally, determinants of VPg uridylylation efficiency suggest formation and/or collapse or release of the uridylylated product as the rate-limiting step in vitro depending upon the VPg donor employed.

The crystal structure of coxsackievirus B3 RNA-dependent RNA polymerase in complex with its protein primer VPg confirms the existence of a second VPg binding site on Picornaviridae polymerases

The RNA-dependent RNA polymerase (RdRp) is a central piece in the replication machinery of RNA viruses. In picornaviruses this essential RdRp activity also uridylates the VPg peptide, which then serves as a primer for RNA synthesis. Previous genetic, binding, and biochemical data have identified a VPg binding site on poliovirus RdRp and have shown that is was implicated in VPg uridylation. More recent structural studies have identified a topologically distinct site on the closely related foot-and-mouth disease virus RdRp supposed to be the actual VPg-primer-binding site. Here, we report the crystal structure at 2.5-A resolution of active coxsackievirus B3 RdRp (also named 3D(pol)) in a complex with VPg and a pyrophosphate. The pyrophosphate is situated in the active-site cavity, occupying a putative binding site either for the coproduct of the reaction or an incoming NTP. VPg is bound at the base of the thumb subdomain, providing first structural evidence for the VPg binding site previously identified by genetic and biochemical methods. The binding mode of VPg to CVB3 3D(pol) at this site excludes its uridylation by the carrier 3D(pol). We suggest that VPg at this position is either uridylated by another 3D(pol) molecule or that it plays a stabilizing role within the uridylation complex. The CVB3 3D(pol)/VPg complex structure is expected to contribute to the understanding of the multicomponent VPg-uridylation complex essential for the initiation of genome replication of picornaviruses.

Structure-function relationships among RNA-dependent RNA polymerases

RNA-dependent RNA polymerases (RdRPs) play key roles in viral transcription and genome replication, as well as epigenetic and post-transcriptional control of cellular gene expression. In this article, we review the crystallographic, biochemical, and molecular genetic data available for viral RdRPs that have led to a detailed description of substrate and cofactor binding, fidelity of nucleotide selection and incorporation, and catalysis. It is likely that the cellular RdRPs will share some of the basic structural and mechanistic principles gleaned from studies of viral RdRPs. Therefore, studies of the viral RdRP establish a framework for the study of cellular RdRPs, an important yet understudied class of nucleic acid polymerases.

Membrane topography of the hydrophobic anchor sequence of poliovirus 3A and 3AB proteins and the functional effect of 3A/3AB membrane association upon RNA replication

Replication of poliovirus RNA takes place on the cytoplasmic surface of membranous vesicles that form after infection of the host cell. It is generally accepted that RNA polymerase 3D(pol) interacts with membranes in a complex with viral protein 3AB, which binds to membranes by means of a hydrophobic anchor sequence that is located near the C-terminus of the 3A domain. In this study, we used fluorescence and fluorescence quenching methods to define the topography of the anchor sequence in the context of 3A and 3AB proteins inserted in model membranes. Mutants with a single tryptophan near the center of the anchor sequence but lacking Trp elsewhere in 3A/3AB were constructed which, after the emergence of suppressor mutations, replicated well in HeLa cells. When a peptide containing the mutant anchor sequence was incorporated in model membrane vesicles, measurements of Trp depth within the lipid bilayer indicated formation of a transmembrane topography. However, rather than the 22-residue length predicted from hydrophobicity considerations, the transmembrane segment had an effective length of 16 residues, such that Gln64 likely formed the N-terminal boundary. Analogous experiments using full-length proteins bound to preformed model membrane vesicles showed that the anchor sequence formed a mixture of transmembrane and nontransmembrane topographies in the 3A protein but adopted only the nontransmembrane configuration in the context of 3AB protein. Studies of the function of 3A/3AB inserted into model membrane vesicles showed that membrane-bound 3AB is highly efficient in stimulating the activity of 3D(pol) in vitro while membrane-bound 3A totally lacks this activity. Moreover, in vitro uridylylation reactions showed that membrane-bound 3AB is not a substrate for 3D(pol), but free VPg released by cleavage of 3AB with proteinase 3CD(pro) could be uridylylated.

Characterization of protein-protein interactions critical for poliovirus replication: analysis of 3AB and VPg binding to the RNA-dependent RNA polymerase

Two critical interactions within the poliovirus RNA replication complex are those of the RNA-dependent RNA polymerase 3D with the viral proteins 3AB and VPg. 3AB is a membrane-binding protein responsible for the localization of the polymerase to the membranous vesicles at which replication occurs. VPg (a peptide comprising the 3B region of 3AB) is the 22-residue soluble product of 3AB cleavage and serves as the protein primer for RNA replication. The detailed interactions of these proteins with the RNA-dependent RNA polymerase 3D were analyzed to elucidate the precise roles of 3AB and VPg in the viral RNA replication complex. Using a membrane-based pull-down assay, we have identified a binding "hot-spot" spanning residues 100 to 104 in the 3B (VPg) region of 3AB which plays a critical role in mediating the interaction of 3AB with the polymerase. Isothermal titration calorimetry shows that the interaction of VPg with 3D is enthalpically driven, with a dissociation constant of 11 microM. Mutational analyses of VPg indicate that a subset of the residues important for 3AB-3D binding are also important for VPg-3D binding. Two residues in particular, P14 and R17, were shown to be absolutely critical for the binding interaction. This work provides the direct characterization of two binding interactions critical for the replication of this important class of viruses and identifies a conserved polymerase binding sequence responsible for targeting the polymerase.

Crystal structure of poliovirus 3CD protein: virally encoded protease and precursor to the RNA-dependent RNA polymerase

Poliovirus 3CD is a multifunctional protein that serves as a precursor to the protease 3C(pro) and the viral polymerase 3D(pol) and also plays a role in the control of viral replication. Although 3CD is a fully functional protease, it lacks polymerase activity. We have solved the crystal structures of 3CD at a 3.4-A resolution and the G64S fidelity mutant of 3D(pol) at a 3.0-A resolution. In the 3CD structure, the 3C and 3D domains are joined by a poorly ordered polypeptide linker, possibly to facilitate its cleavage, in an arrangement that precludes intramolecular proteolysis. The polymerase active site is intact in both the 3CD and the 3D(pol) G64S structures, despite the disruption of a network proposed to position key residues in the active site. Therefore, changes in molecular flexibility may be responsible for the differences in fidelity and polymerase activities. Extensive packing contacts between symmetry-related 3CD molecules and the approach of the 3C domain's N terminus to the VPg binding site suggest how 3D(pol) makes biologically relevant interactions with the 3C, 3CD, and 3BCD proteins that control the uridylylation of VPg during the initiation of viral replication. Indeed, mutations designed to disrupt these interfaces have pronounced effects on the uridylylation reaction in vitro.

Picornavirus genome replication: assembly and organization of the VPg uridylylation ribonucleoprotein (initiation) complex

All picornaviruses have a protein, VPg, covalently linked to the 5'-ends of their genomes. Uridylylated VPg (VPg-pUpU) is thought to serve as the protein primer for RNA synthesis. VPg-pUpU can be produced in vitro by the viral polymerase, 3Dpol, in a reaction in which a single adenylate residue of a stem-loop structure, termed oriI, templates processive incorporation of UMP into VPg by using a "slide-back" mechanism. This reaction is greatly stimulated by viral precursor protein 3CD or its processed derivative, 3C; both contain RNA-binding and protease activities. We show that the 3C domain encodes specificity for oriI, and the 3D domain enhances the overall affinity for oriI. Thus, 3C(D) stimulation exhibits an RNA length dependence. By using a minimal system to evaluate the mechanism of VPg uridylylation, we show that the active complex contains polymerase, oriI, and 3C(D) at stoichiometry of 1:1:2. Dimerization of 3C(D) is supported by physical and structural data. Polymerase recruitment to and retention in this complex require a protein-protein interaction between the polymerase and 3C(D). Physical and functional data for this interaction are provided for three picornaviruses. VPg association with this complex is weak, suggesting that formation of a complex containing all necessary components of the reaction is rate-limiting for the reaction. We suggest that assembly of this complex in vivo would be facilitated by use of precursor proteins instead of processed proteins. These data provide a glimpse into the organization of the ribonucleoprotein complex that catalyzes this key step in picornavirus genome replication.

Role of RNA structure and RNA binding activity of foot-and-mouth disease virus 3C protein in VPg uridylylation and virus replication

The uridylylation of the VPg peptide primer is the first stage in the replication of picornavirus RNA. This process can be achieved in vitro using purified components, including 3B (VPg) with the RNA dependent RNA polymerase (3Dpol), the precursor 3CD, and an RNA template containing the cre/bus. We show that certain RNA sequences within the foot-and-mouth disease virus (FMDV) 5' untranslated region but outside of the cre/bus can enhance VPg uridylylation activity. Furthermore, we have shown that the FMDV 3C protein alone can substitute for 3CD, albeit less efficiently. In addition, the VPg precursors, 3B(3)3C and 3B(123)3C, can function as substrates for uridylylation in the absence of added 3C or 3CD. Residues within the FMDV 3C protein involved in interaction with the cre/bus RNA have been identified and are located on the face of the protein opposite from the catalytic site. These residues within 3C are also essential for VPg uridylylation activity and efficient virus replication.

Novel, structure-based mechanism for uridylylation of the genome-linked peptide (VPg) of picornaviruses

The VPg peptide, which is found in poliovirus infected cells either covalently bound to the 5'-end of both plus and minus strand viral RNA, or in a uridylylated free form, is essential for picornavirus replication. Combining experimental structure and mutation results with molecular modeling suggests a new mechanism for VPg uridylylation, which assigns an additional function, that of scaffold, to the polymerase. The polarity of the NMR structure of VPg is complementary to the binding site on the surface of poliovirus polymerase determined previously by mutagenesis. Docking VPg at this position places the reactive tyrosinate close to the 5'-end of Poly(A)7 RNA when this is bound with its 3'-end in the active site of the polymerase. The triphosphate tail of a UTP moiety, base paired with the 5'-end of the RNA, projects back over the Tyr3-OH and is held in position by conserved positively charged side-chains of VPg. Other conserved residues mediate binding to the polymerase surface and serve as ligands for metal ion catalyzed transphosphorylation. Additional viral proteins or a second polymerase molecule may aid in stabilizing the components of the reaction. In the model complex, VPg can direct its own uridylylation before entering the polymerase active site.

NMR structure of the viral peptide linked to the genome (VPg) of poliovirus

VPgs are essential for replication of picornaviruses, which cause diseases such as poliomyelitis, foot and mouth disease, and the common cold. VPg in infected cells is covalently linked to the 5' end of the viral RNA, or, in a uridylylated form, free in the cytoplasm. We show here the first solution structure for a picornaviral VPg, that of the 22-residue peptide from poliovirus serotype 1. VPg in buffer is inherently flexible, but a single conformer was obtained by adding trimethylamine N-oxide (TMAO). TMAO had only minor effects on the TOCSY spectrum. However, it increased the amount of structured peptide, as indicated by more peaks in the NOESY spectrum and an up to 300% increase in the ratio of normalized NOE cross peak intensities to that in buffer. The data for VPg in TMAO yielded a well defined structure bundle with 0.6 A RMSD (versus 6.6 A in buffer alone), with 10-30 unambiguous constraints per residue. The structure consists of a large loop region from residues 1 to 14, from which the reactive tyrosinate projects outward, and a C-terminal helix from residues 18 to 21 that aligns the sidechains of conserved residues on one face. The structure has a stable docking position at an area on the poliovirus polymerase crystal structure identified as a VPg binding site by mutagenesis studies. Further, UTP and ATP dock in a base-specific manner to the reactive face of VPg, held in place by residues conserved in all picornavirus VPgs.

The structure of a protein primer-polymerase complex in the initiation of genome replication

Picornavirus RNA replication is initiated by the covalent attachment of a UMP molecule to the hydroxyl group of a tyrosine in the terminal protein VPg. This reaction is carried out by the viral RNA-dependent RNA polymerase (3D). Here, we report the X-ray structure of two complexes between foot-and-mouth disease virus 3D, VPg1, the substrate UTP and divalent cations, in the absence and in the presence of an oligoadenylate of 10 residues. In both complexes, VPg fits the RNA binding cleft of the polymerase and projects the key residue Tyr3 into the active site of 3D. This is achieved by multiple interactions with residues of motif F and helix alpha8 of the fingers domain and helix alpha13 of the thumb domain of the polymerase. The complex obtained in the presence of the oligoadenylate showed the product of the VPg uridylylation (VPg-UMP). Two metal ions and the catalytic aspartic acids of the polymerase active site, together with the basic residues of motif F, have been identified as participating in the priming reaction.

Nucleotide channel of RNA-dependent RNA polymerase used for intermolecular uridylylation of protein primer

Poliovirus VPg is a 22 amino acid residue peptide that serves as the protein primer for replication of the viral RNA genome. VPg is known to bind directly to the viral RNA-dependent RNA polymerase, 3D, for covalent uridylylation, yielding mono and di-uridylylated products, VPg-pU and VPg-pUpU, which are subsequently elongated. To model the docking of the VPg substrate to a putative VPg-binding site on the 3D polymerase molecule, we performed a variety of structure-based computations followed by experimental verification. First, potential VPg folded structures were identified, yielding a suite of predicted beta-hairpin structures. These putative VPg structures were then docked to the region of the polymerase implicated by genetic experiments to bind VPg, using grid-based and fragment-based methods. Residues in VPg predicted to affect binding were identified through molecular dynamics simulations, and their effects on the 3D-VPg interaction were tested computationally and biochemically. Experiments with mutant VPg and mutant polymerase molecules confirmed the predicted binding site for VPg on the back side of the polymerase molecule during the uridylylation reaction, opposite to that predicted to bind elongating RNA primers.