Treatment options for COVID‐19, caused by SARS‐CoV‐2, remain limited. Understanding viral pathogenesis at the molecular level is critical to develop effective therapy. Some recent studies have explored SARS‐CoV‐2–host interactomes and provided great resources for understanding viral replication. However, host proteins that functionally associate with SARS‐CoV‐2 are localized in the corresponding subnetwork within the comprehensive human interactome. Therefore, constructing a downstream network including all potential viral receptors, host cell proteases, and cofactors is necessary and should be used as an additional criterion for the validation of critical host machineries used for viral processing. This study applied both affinity purification mass spectrometry (AP‐MS) and the complementary proximity‐based labeling MS method (BioID‐MS) on 29 viral ORFs and 18 host proteins with potential roles in viral replication to map the interactions relevant to viral processing. The analysis yields a list of 693 hub proteins sharing interactions with both viral baits and host baits and revealed their biological significance for SARS‐CoV‐2. Those hub proteins then served as a rational resource for drug repurposing via a virtual screening approach. The overall process resulted in the suggested repurposing of 59 compounds for 15 protein targets. Furthermore, antiviral effects of some candidate drugs were observed in vitro validation using image‐based drug screen with infectious SARS‐CoV‐2. In addition, our results suggest that the antiviral activity of methotrexate could be associated with its inhibitory effect on specific protein–protein interactions.

Twelve members of the family 11 xylanases, including both mesophilic and thermophilic proteins, were studied using molecular dynamics (MD). Simulations of xylanases were carried out in an explicit water environment at four different temperatures, 300, 400, 500 and 600 K. A difference in thermotolerance between mesophilic and thermophilic xylanases became clear: thermophilic xylanases endured heat in higher simulation temperatures better than mesophilic ones. The unfolding pathways seemed to be similar for all simulations regardless of the protein. The unfolding initiates at the N-terminal region or alternatively from the α-helix region and proceeds to the ‘finger region’. Unfolding of these regions led to denaturated structures within the 4.5 ns simulation at 600 K. The results are in agreement with experimental mutant studies. The results show clearly that the stability of the protein is not evenly distributed over the whole structure. The MD analysis suggests regions in the protein structure which are more unstable and thus potential targets for mutation experiments to improve thermostability.

A mutated KRAS protein is frequently observed in human cancers. Traditionally, the oncogenic properties of KRAS missense mutants at position 12 (G12X) have been considered as equal. Here, by assessing the probabilities of occurrence of all KRAS G12X mutations and KRAS dynamics we show that this assumption does not hold true. Instead, our findings revealed an outstanding mutational bias. We conducted a thorough mutational analysis of KRAS G12X mutations and assessed to what extent the observed mutation frequencies follow a random distribution. Unique tissue-specific frequencies are displayed with specific mutations, especially with G12R, which cannot be explained by random probabilities. To clarify the underlying causes for the nonrandom probabilities, we conducted extensive atomistic molecular dynamics simulations (170 μs) to study the differences of G12X mutations on a molecular level. The simulations revealed an allosteric hydrophobic signaling network in KRAS, and that protein dynamics is altered among the G12X mutants and as such differs from the wild-type and is mutation-specific. The shift in long-timescale conformational dynamics was confirmed with Markov state modeling. A G12X mutation was found to modify KRAS dynamics in an allosteric way, which is especially manifested in the switch regions that are responsible for the effector protein binding. The findings provide a basis to understand better the oncogenic properties of KRAS G12X mutants and the consequences of the observed nonrandom frequencies of specific G12X mutations. The oncogene KRAS is frequently mutated in various cancers. When the amino acid glycine 12 is mutated, KRAS protein acquires oncogenic properties that result in tumor cell-growth and cancer progression. These mutations prevail especially in the pancreatic ductal adenocarcinoma, which is a cancer with an exceptionally dismal prognosis. To date, there is a limited understanding of the different mutations at the position 12, also regarding whether the different mutations would have different consequences. These discrepancies could have major implications for the future drug therapies targeting KRAS mutant harboring tumors. In this study, we made a critical assessment of the observed frequency of KRAS G12X mutations and the underlying causes for these frequencies. We also assessed KRAS G12X mutant discrepancies on an atomistic level by utilizing state-of-the-art molecular dynamics simulations. We found that the dynamics of the mutants does not only differ from the wild-type protein, but there is also a profound difference among the different mutants. These results emphasize that the different KRAS G12X mutations are not equal, and thereby they suggest that the future research related to mutant KRAS biology should account for these observations.