The use of remdesivir to treat COVID-19 will likely continue before clinical trials are completed. Due to the lengthening pandemic and evolving nature of the virus, predicting potential residues prone to mutation is crucial for the management of remdesivir resistance. Using a rational ligand-based interface design complemented with mutational mapping, we generated a total of 100,000 mutations and provided insight into the functional outcomes of mutations in the remdesivir-binding site in nsp12 subunit of RdRp. After designing 46 residues in the remdesivir-binding site of nsp12, the designs retained 97%–98% sequence identity, suggesting that very few mutations in nsp12 are required for SARS-CoV-2 to attain remdesivir resistance. Several mutants displayed decreased binding affinity to remdesivir, suggesting drug resistance. These hotspot residues had a higher probability of undergoing selective mutation and thus conferring remdesivir resistance. Identifying the potential residues prone to mutation improves our understanding of SARS-CoV-2 drug resistance and COVID-19 pathogenesis.

Several missense mutations in the coding region of angiogenin (ANG) gene have been identified in Amyotrophic Lateral Sclerosis (ALS) patients. These mutations lead to loss of either ribonucleolytic activity or nuclear translocation activity or both of ANG (protein encoded by ANG gene) causing ALS. We present here a cohesive and comprehensive picture of the molecular origins of functional loss of all ALS associated ANG mutants, emerging via extensive molecular dynamics simulations. Our method effectively predicts that conformational change of His114 results in loss of ribonucleolytic activity and that reduction of solvent accessible surface area of nuclear localization signal residues ³¹RRR³³ results in loss of nuclear translocation activity. These predictions hold true, without exception, for all ANG mutants studied and can be employed to infer whether a new ANG mutation is causative of ALS or benign ahead of experimental findings.

Mutations in the coding region of angiogenin (ANG) gene have been found in patients suffering from Amyotrophic Lateral Sclerosis (ALS). Neurodegeneration results from the loss of angiogenic ability of ANG (protein coded by ANG). In this work, we performed extensive molecular dynamics (MD) simulations of wild-type ANG and disease associated ANG variants to elucidate the mechanism behind the loss of ribonucleolytic activity and nuclear translocation activity, functions needed for angiogenesis. MD simulations were carried out to study the structural and dynamic differences in the catalytic site and nuclear localization signal residues between WT-ANG (Wild-type ANG) and six mutants. Variants K17I, S28N, P112L and V113I have confirmed association with ALS, while T195C and A238G single nucleotide polymorphisms (SNPs) encoding L35P and K60E mutants respectively, have not been associated with ALS. Our results show that loss of ribonucleolytic activity in K17I is caused by conformational switching of the catalytic residue His114 by 99°. The loss of nuclear translocation activity of S28N and P112L is caused by changes in the folding of the residues ³¹RRR³³ that result in the reduction in solvent accessible surface area (SASA). Consequently, we predict that V113I will exhibit loss of angiogenic properties by loss of nuclear translocation activity and L35P by loss of both ribonucleolytic activity and nuclear translocation activity. No functional loss was inferred for K60E. The MD simulation results were supported by hydrogen bond interaction analyses and molecular docking studies. Conformational switching of catalytic residue His114 seems to be the mechanism causing loss of ribonucleolytic activity and reduction in SASA of nuclear localization signal residues ³¹RRR³³ results in loss of nuclear translocation activity in ANG mutants. Therefore, we predict that L35P mutant, would exhibit loss of angiogenic functions, and hence would correlate with ALS while K60E would not show any loss.