Molecular dynamics (MD) simulations are powerful theoretical methods that can reveal biomolecular properties, such as structure, fluctuations, and ligand binding, at the atomic level. All-atom MD simulations elucidated a difference in the dynamic properties of RNA-dependent RNA polymerases (RdRps) in severe acute respiratory syndrom coronavirus 2 (SARS-CoV-2) and SARS-CoV, which may cause activity differences of these RdRps. RdRp is also a drug target for Coronavirus disease 2019. Nucleotide analogs, such as remdesivir and favipiravir, are considered to be taken up by RdRp and inhibit RNA replication. The recognition mechanism of RdRp for these drug molecules and adenosine triphosphate (ATP) was revealed by MD simulations at the atomic detail. In addition, various simulation studies on the complexes of SARS-CoV-2 RdRp with several nucleotide analogs are also presented.
To observe the tertiary-structure difference between nsp12s of SARS-CoV-2 and SARS-CoV, the average distances between Cα atoms in nsp12s were calculated as shown in Figure 3a and 3b. We can see that the two systems have the following in common: the NiRAN and palm domains are spatially close to each other, and the interface and fingers domains are close to each other. To clarify the differences in the average distances between nsp12s of SARS-CoV-2 and SARS-CoV, Ito et al. calculated the ratio of the difference (Figure 3c). The differences between nsp12s of SARS-CoV-2 and SARS-CoV are observed in the region indicated by the brown square. Blue lines (or blue meshes) are observed in residues around 430, 520, 560, 620, 690, 760, and 800. These results mean that the distances between all motifs of nsp12 in SARS-CoV are shorter than those of nsp12 in SARS-CoV-2. In particular, the distance between motifs F and G for SARS-CoV is up to 63% shorter than that for SARS-CoV-2.
Figure 3. The average distances between Cα atoms of nsp12 for (a) SARS-CoV-2 and (b) SARS-CoV. The borders between the domains in nsp12 are indicated by the green lines. (c) The ratios of the differences between the average distances for SARS-CoV nsp12 and those for SARS-CoV-2 nsp12. The brown square shows residues that have large differences. Reproduced with permission from Ref. [17].
In addition, dynamic cross-correlation (DCC) was calculated to investigate the correlation between domain motions. DCCs of SARS-CoV-2 nsp12 and SARS-CoV nsp12 are presented in Figure 4a and 4b. Here, red and blue indicate positive and negative correlations, respectively. The fact that there is a positive (negative) correlation between two residues indicates that the motions of these residues are in the same (opposite) direction. In both systems, positive correlations are found between most residues within the same domains. However, there are both positive and negative correlations in the interface domain of SARS-CoV nsp12. The boundary between these correlations is residue 330. Residues before and after residue 330 in the interface domain are positively correlated with the NiRAN domain and fingers domain, respectively. Figure 4c shows the differences in DCCs between SARS-CoV-2 and SARS-CoV nsp12s. As shown by the region surrounded by the brown lines, the differences are larger in the NiRAN and interface domains. These domains before residue 330 have a strong negative correlation with the fingers domain in SARS-CoV nsp12. That is, the regions before residue 330 move cooperatively with the fingers domain, moving closer and further away from each other.
Figure 4. DCCs of nsp12 for (a) SARS-CoV-2 and (b) SARS-CoV. The borders between the domains in nsp12 are indicated by the green lines. (c) Differences between DCCs for SARS-CoV nsp12 and those for SARS-CoV-2 nsp12. The region surrounded by the brown lines means residues with large differences. Reproduced with permission from Ref. [17].
As shown in Figure 3c, the distances between all motifs in SARS-CoV nsp12 are shorter compared to SARS-CoV-2 nsp12. This may enhance the RdRp activity of SARS-CoV. Furthermore, in SARS-CoV nsp12, the NiRAN and fingers domains move cooperatively toward and away from each other; because the removal of the NiRAN domain reduces the RdRp activity[26], the NiRAN domain is important for the RdRp activities. The cooperative movement of the NiRAN domain with the core (fingers) domain of RdRp may also enhance the activity of RdRp.
In this section, Tanimoto et al.'s MD simulations of RdRp with RemTP, FavTP, or ATP to clarify how RdRp recognizes the drugs and NTPs are presented[19].
As a result of the MD simulations, the ligand recognition process by RdRp was observed in all three systems of RemTP, FavTP, and ATP. First, the ligand recognition probability was calculated, as listed in Table 1. RemTP shows the highest probability, FavTP shows the second-highest probability, followed by ATP, although within the statistical errors. These results are in qualitative agreement with previous experimental studies[11][27]. In addition, MD simulations of the RdRp-RemTP complex using the free energy perturbation (FEP) method showed that RemTP is bound more strongly to RdRp than ATP[28], which is also consistent with the present results.
Table 1. The number of MD simulations in which RdRp recognized the ligands. Ligand recognition probability is also listed. Reproduced with permission from Ref. [19].
Ligand |
Ligand Recognition/Total |
Ligand Recognition Probability |
RemTP |
12/50 |
0.24 ± 0.07 |
FavTP |
9/50 |
0.18 ± 0.06 |
ATP |
7/50 |
0.14 ± 0.06 |
Next, to understand the mechanism of the ligand recognition by RdRp, the trajectories of the recognized ligands were examined. As a result, an interesting path was observed in which the lysine residues of RdRp carry ligands to the binding site like a “bucket brigade,” as shown in Figure 5. In this path, the phosphate groups of the ligands contacted LYS2 and LYS43 of nsp7 and LYS551, LYS621, and LYS798 of nsp12. Because nsp12 and nsp7 correspond to chain A and chain C, respectively, in the original cryo-EM structure, the residues are expressed here as “chain label + residue number + residue name”. These lysine residues have a positive charge, of which C2LYS, C43LYS, and A551LYS are in a line toward the binding site. In this process of ligand transportation, the phosphate groups of RemTP first interact with the side chain of C2LYS (state 1 (S1), Figure 5b). C2LYS passes RemTP to C43LYS, which is spatially close (state 2 (S2), Figure 5c). C43LYS then passes RemTP to A551LYS (state 3 (S3), Figure 5d). RemTP finally reaches the binding site (state 4 (S4), Figure 5e). The ligand also interacts electrically with A621LYS and A798LYS at the binding site. A similar process was also observed in the FavTP and ATP systems.
Figure 5. (a) “Bucket brigade” trajectory of RemTP recognized by RdRp. The black circles mean the positions at which RemTP has contact with RdRp residues. (b–e) Typical snapshot at each state (S1–S4). In (b–e), the lysine residues that contributed to the ligand recognition and RemTP are expressed as blue and red stick models, respectively. Reproduced with permission from Ref. [19].
These positively charged residues have been reported to be favorable for the NTP recognition[10][11]. Furthermore, the lysine residues, which contribute to the bucket-brigade ligand transportation, are highly conserved in RdRp of SARS-CoV[12]. Therefore, it is expected that for both SARS-CoV-2 and SARS-CoV RdRps, these linearly arranged lysine residues carry NTPs to the binding site, thereby enhancing the NTP recognition ability of RdRp.
It is important to understand the mechanism by which RemTP is bound to RdRp and inhibits the RNA replications. Zhang et al. examined how remdesivir integrated into the nascent RNA strand (N-Rem) inhibited RdRp from adding nucleotides to the strand[29]. They revealed that N-Rem led to a delayed chain termination, where the translocation of the nascent RNA strand is terminated once three nucleotides were added after the RemTP incorporation. It was clarified that the forward translocation of the nascent RNA strand was impeded by the electrostatic repulsion between ASP865 and the 1′-cyano group of N-Rem as well as the steric clash between SER861 and the 1′-cyano group in a position where N-Rem reaches after three nucleotide incorporations. They also found that N-Rem at this site greatly weakened the hydrogen bonds of base pairs with its template uracil due to the electrostatic attraction between LYS593 and the 1′-cyano group. Their simulation study showed that the 1′-cyano group on the ribose was essential for remdesivir to inhibit the RdRp function.
Using MD simulations and quantum mechanics/molecular mechanics simulations for SARS-CoV-2 RdRp with an RNA duplex, Aranda et al. reported the detailed mechanisms of the binding and incorporation of natural nucleotides and RemTP[30]. They found that RemTP was preferentially bound to RdRp over ATP, while it was incorporated into the nascent RNA strand with an efficiency only slightly lower than ATP. In addition, they reported that, unlike the simulation results obtained by Zhang et al.[29], no steric clash was detected between N-Rem and the residues of RdRp when the nascent RNA strand was translocated along the exit channel. Instead, they found that N-Rem was trapped at a position where the three nucleotides were incorporated after RemTP. Therefore, they suggested that either non-covalent or transient-covalent bonds between the 1′-cyano group of N-Rem at this position and hydroxyl group of SER861 could act as a trap for the nascent RNA strand and stall the translocation of the duplex.
Luo et al. performed MD simulations to elucidate the nascent-RNA-synthesis inhibition mechanism by remdesivir embedded in the template strand (T-Rem)[31]. Experimental observations have shown that T-Rem inhibits the synthesis of the nascent RNA strand[32]. They revealed that when T-Rem was at the binding site, the translocation of T-Rem was hampered by the hydrogen-bond formation between the 1′-cyano group of T-Rem and the backbone of GLY683 and the steric clash between the 1′-cyano group and the backbone of SER682.
Many simulation studies have been conducted on the inhibition mechanism of the RdRp function by nucleotide analogs other than remdesivir. Yuan et al. investigated the inhibitory effect of nucleotide analogs with various 2′ modifications against SARS-CoV-2 RdRp using MD simulations and FEP methods[33]. The nucleotide analogs included 2′-O-methyl uridine triphosphate (OMU-TP), sofosbuvir triphosphate (SFU-TP), 2′-C-methyl cytidine triphosphate (CMC-TP), Gemcitabine triphosphate (GMC-TP), and ara-uridine triphosphate (ARU-TP). Previous experimental studies reported that three of these, OMU-TP, SFU-TP, and CMC-TP, act as effective inhibitors, while GMC-TP and ARU-TP have no inhibitory effects[34][35][36]. They revealed that OMU decreased the binding probability of the subsequent NTP and consequently caused partial chain terminations due to the steric hindrance by its 2′-O-methyl modification. In addition, it was found that the bulky 2′-methyl substitutions in SFU and CMC largely disrupted the binding site, leading to the immediate chain termination. In contrast, GMC and ARU, which have smaller 2′ substitutions such as the fluorine atoms and ara-hydroxyl group, showed marginal effects on the polymerization process upon the incorporation. Their simulation results were consistent with previous experimental results[34][35][36] and elucidated the detailed inhibition mechanisms of 2′ substituted nucleotide analogs against SARS-CoV-2 RdRp.
Li et al. systematically investigated the inhibitory effects of ATP analogs possessing 2′ or 3′ ribose modifications against SARS-CoV-2 RdRp using MD simulations and FEP methods[37]. The analogs included clofarabine triphosphate, didanosine triphosphate, fludarabine triphosphate, vidarabine triphosphate, 2′-amino-2′-deoxyadenosine triphosphate, 2′,3′-didehydro-2′,3′-dideoxyadenosine triphosphate, and cordycepin triphosphate. They found that clofarabine and fludarabine could not form stable binding at the binding site and only had a minor effect on the next nucleotide incorporation into the nascent strand. It was also clarified that vidarabine and 2′-amino-2′-deoxyadenosine could not efficiently inhibit the incorporation of the next substrate, although they could be incorporated into the nascent strand as the substrate. Didanosine, 2′,3′-didehydro-2′,3′-dideoxyadenosine, and cordycepin could also be incorporated into the nascent strand and had the capability to terminate the next nucleotide addition while 2′,3′-didehydro-2′,3′-dideoxyadenosine triphosphate was less competitive than the other two analogs. Therefore, they concluded that substituting the 3′-hydroxyl group with one hydrogen atom would inherently inhibit the next nucleotide addition when it appears at the 3′ terminal of the nascent strand. They proposed that cordycepin and didanosine were promising nucleotide analogs as immediate terminators.
Nucleotide analogs inhibit the RNA replications by interfering with the addition of the next nucleotide (immediate chain termination)[34] or by interfering with the translocation of the nascent RNA strand after the incorporation of three nucleotides (delayed chain termination)[27]. It was shown that nucleotide analogs embedded in the template strand also inhibit the RNA replications (template-dependent inhibition)[32]. The simulation studies presented here elucidated the atomic level mechanisms of the delayed chain termination of remdesivir, immediate chain termination of 2′ and 3′ modified nucleotide analogs, and template-dependent inhibition of remdesivir. Furthermore, these simulation studies also suggested more promising nucleotide analogs for inhibiting the function of RdRp than remdesivir[33][37]. All simulation studies focused on the situation after NTPs or nucleotide analogs are incorporated into the binding site of RdRp. On the other hand, the simulation study described in Section 3[19] focused on the process by which ligands far from RdRp were incorporated into the binding site and revealed the bucket-brigade transport mechanism of NTPs and nucleotide analogs by lysine residues of RdRp. Overall, the simulation studies presented here help us in enhancing the understanding on how nucleotide analogs are recognized by RdRp and inhibit the RNA replication at the atomic level.
This entry is adapted from the peer-reviewed paper 10.3390/ijms231810358