Computational and theoretical exploration for clinical suitability of Remdesivir drug to SARS-CoV-2

Design single-stranded DNA aptamer of cluster of differentiation 47 protein by stochastic tunnelling-basin hopping-discrete molecular dynamics method

The formation of the Cluster of Differentiation 47 (CD47, PDB code: 2JJT)/signal regulatory protein α (SIRPα) complex is very important as it protects healthy cells from immune clearance while promoting macrophage phagocytosis for tumour elimination. Although several antibodies have been developed for cancer therapy, new function-blocking aptamers are still under development. This study aims to design the aptamer AptCD47, which can block the formation of the CD47/SIRPα complex. This study employs the MARTINI coarse-grained (CG) force field and the stochastic tunnelling-basin hopping-discrete molecular dynamics (STUN-BH-DMD) method to identify the most stable AptCD47/CD47 complexes. Coarse-grained molecular dynamics (CGMD) simulations were used to obtain root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) analyses. The results demonstrate that the formation of AptCD47/CD47 complexes renders the CD47 structure more stable than the single CD47 molecule in a water environment. The minimum energy pathway (MEP) obtained by the nudged elastic band (NEB) method indicates that the binding processes of 5'-ATTCAATTCC-3' and 5'-AGTGCAATCT-3' to CD47 are barrierless, which is much lower than the binding barrier of SIRPα to CD47 of about 14.23 kcal/mol. Therefore, these two AptCD47/CD47 complexes can create a high spatial binding barrier for SIRPα, preventing the formation of a stable CD47/SIRPα complex. The proposed numerical process with the MARTINI CG force field can be used to design CD47 aptamers that efficiently block SIRPα from binding to CD47.Communicated by Ramaswamy H. Sarma.

Prediction and validation of monoclonal antibodies separation in aqueous two-phase system using molecular dynamic simulation

In order to predict how mAbs partition in 20% ethylene oxide/80% propylene oxide (v/v) random copolymer (EO₂₀PO₈₀)/water aqueous two-phase system (ATPS), a molecular dynamic simulation model was developed using Gromacs and then validated by experiments. The ATPS was applied with seven kinds of salt, including buffer salt and strong dissociation salt that were commonly employed in the purification of protein. Na₂SO₄ was shown to have the best effects on lowering EO₂₀PO₈₀ content in the aqueous phase and enhancing recovery. The content of EO₂₀PO₈₀ in the sample solution was decreased to 0.62%±0.25% and the recovery of rituximab increased to 97.88%±0.95% by adding 300 mM Na₂SO₄ into back extraction ATPS. The viability determined by ELISA was 95.57% at the same time. A strategy for constructing a prediction model for the distribution of mAbs in ATPS was proposed in consideration of this finding. Partition of trastuzumab in ATPS was predicted by the model created using this method and the prediction result was further validated by experiments. The recovery of trastuzumab reached 95.63%±2.86% under the ideal extraction conditions suggested by the prediction model.

The role of S477N mutation in the molecular behavior of SARS-CoV-2 spike protein: An in-silico perspective

The attachment of SARA-CoV-2 happens between ACE2 and the receptor binding domain (RBD) on the spike protein. Mutations in this domain can affect the binding affinity of the spike protein for ACE2. S477N, one of the most common mutations reported in the recent variants, is located in the RBD. Today's computational approaches in biology, especially during the SARS-CoV-2 pandemic, assist researchers in predicting a protein's behavior in contact with other proteins in more detail. In this study, we investigated the interactions of the S477N-hACE2 in silico to find the impact of this mutation on its binding affinity for ACE2 and immunity responses using dynamics simulation, protein-protein docking, and immunoinformatics methods. Our computational analysis revealed an increased binding affinity of N477 for ACE2. Four new hydrogen and hydrophobic bonds in the mutant RBD-ACE2 were formed (with S19 and Q24 of ACE2), which do not exist in the wild type. Also, the protein spike structure in this mutation was associated with an increase in stabilization and a decrease in its fluctuations at the atomic level. N477 mutation can be considered as the cause of increased escape from the immune system through MHC-II.

Artificially expanded genetic information systems (AEGISs) as potent inhibitors of the RNA-dependent RNA polymerase of the SARS-CoV-2

The recent outbreak of the SARS-CoV-2 infection has affected the lives and economy of more than 200 countries. The unavailability of virus-specific drugs has created an opportunity to identify potential therapeutic agents that can control the rapid transmission of this pandemic. Here, the mechanisms of the inhibition of the RNA-dependent RNA polymerase (RdRp), responsible for the replication of the virus in host cells, are examined by different ligands, such as Remdesivir (RDV), Remdesivir monophosphate (RMP), and several artificially expanded genetic information systems (AEGISs) including their different sequences by employing molecular docking, MD simulations, and MM/GBSA techniques. It is found that the binding of RDV to RdRp may block the RNA binding site. However, RMP would acquire a partially flipped conformation and may allow the viral RNA to enter into the binding site. The internal dynamics of RNA and RdRp may help RMP to regain its original position, where it may inhibit the RNA-chain elongation reaction. Remarkably, AEGISs are found to obstruct the binding site of RNA. It is shown that dPdZ, a two-nucleotide sequence containing P and Z would bind to RdRp very strongly and may occupy the positions of two nucleotides in the RNA strand, thereby denying access of the substrate-binding site to the viral RNA. Thus, it is proposed that the AEGISs may act as novel therapeutic candidates against the SARS-CoV-2. However, in vivo evaluations of their potencies and toxicities are needed before using them against COVID-19.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 spike evolutionary behaviors; simulation of N501Y mutation outcomes in terms of immunogenicity and structural characteristic

Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2), a large number of mutations in its genome have been reported. Some of the mutations occur in noncoding regions without affecting the pathobiology of the virus, while mutations in coding regions are significant. One of the regions where a mutation can occur, affecting the function of the virus is at the receptor‐binding domain (RBD) of the spike protein. RBD interacts with angiotensin‐converting enzyme 2 (ACE2) and facilitates the entry of the virus into the host cells. There is a lot of focus on RBD mutations, especially the displacement of N501Y which is observed in the UK/Kent, South Africa, and Brazilian lineages of SARS‐CoV‐2. Our group utilizes computational biology approaches such as immunoinformatics, protein–protein interaction analysis, molecular dynamics, free energy computation, and tertiary structure analysis to disclose the consequences of N501Y mutation at the molecular level. Surprisingly, we discovered that this mutation reduces the immunogenicity of the spike protein; also, displacement of Asn with Tyr reduces protein compactness and significantly increases the stability of the spike protein and its affinity to ACE2. Moreover, following the N501Y mutation secondary structure and folding of the spike protein changed dramatically.

Estrogen Receptor Modulators in Viral Infections Such as SARS-CoV-2: Therapeutic Consequences

COVID-19 is a pandemic respiratory disease caused by the SARS-CoV-2 coronavirus. The worldwide epidemiologic data showed higher mortality in males compared to females, suggesting a hypothesis about the protective effect of estrogens against severe disease progression with the ultimate end being patient's death. This article summarizes the current knowledge regarding the potential effect of estrogens and other modulators of estrogen receptors on COVID-19. While estrogen receptor activation shows complex effects on the patient's organism, such as an influence on the cardiovascular/pulmonary/immune system which includes lower production of cytokines responsible for the cytokine storm, the receptor-independent effects directly inhibits viral replication. Furthermore, it inhibits the interaction of IL-6 with its receptor complex. Interestingly, in addition to natural hormones, phytestrogens and even synthetic molecules are able to interact with the estrogen receptor and exhibit some anti-COVID-19 activity. From this point of view, estrogen receptor modulators have the potential to be included in the anti-COVID-19 therapeutic arsenal.

Target SARS-CoV-2: theoretical exploration on clinical suitability of certain drugs

We propose a unique theoretical methodology because of the global high priority rating to search for the repurposed drugs that outfit clinical suitability to SARS-CoV-2. The approach is based on the exploration of structural analysis, computation of biothermodynamics, interactions and the prediction of entropy sign successively via molecular dynamics. We tested this methodology for Favipiravir/Dolutegravir drugs on the apo form of SARS-CoV-2 main protease. This theoretical exploration not only suggested the presence of strong interactions between (SARS-CoV-2 + Favipiravir/Dolutegravir) but also emphasized the clinical suitability of Favipiravir over Dolutegravir to treat SARS-CoV-2 main protease. The supremacy of Favipiravir over Doultegravir is well supported by the results of global clinical trials on SARS-CoV-2 infection. Thus, this work will pave the way for incremental advancement towards future design and development of more specific inhibitors to treat SARS-CoV-2 infection in humans.Communicated by Ramaswamy H. Sarma.