Identification of terpenoids as potential inhibitors of SARS-CoV-2 (main protease) and spike (RBD) via computer-aided drug design

The scientific community has been faced with a major challenge in the fight against the SARS-CoV-2 virus responsible for the COVID-19 pandemic, due to the lack of targeted antiviral drugs. To address this issue, we used an in silico approach to screen 23 natural compounds from the terpenoid class for their ability to target key SARS-CoV-2 therapeutic proteins. The results revealed that several compounds showed promising interactions with SARS-CoV-2 proteins, specifically the main protease and the spike receptor binding domain. The molecular docking analysis revealed the importance of certain residues, such as GLY143, SER144, CYS145 and GLU166, in the main protease of the SARS-CoV-2 protein, which play a crucial role in interactions with the ligand. In addition, our study highlighted the importance of interactions with residues GLY496, ARG403, SER494 and ARG393 of the spike receptor-binding domain within the SARS-CoV-2 protein. ADMET and drug similarity analyses were also performed, followed by molecular dynamics and MM-GBSA calculations, to identify potential drugs could be repurposed to combat COVID-19. Indeed, the results suggest that certain terpenoid compounds of plant origin have promising potential as therapeutic targets for SARS-CoV-2. However, additional experimental studies are required to confirm their efficacy as drugs against COVID-19.Communicated by Ramaswamy H. Sarma.

Molecular mechanism of inhibitor bindings to bromodomain-containing protein 9 explored based on molecular dynamics simulations and calculations of binding free energies

Recently, bromodomain-containing protein 9 (BRD9) has been a prospective therapeutic target for anticancer drug design. Molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were adopted to explore binding modes of three inhibitors (5SW, 5U2, and 5U6) to BRD9 and identify the hot spot of the inhibitor-BRD9 binding. The results indicate that the inhibitor 5SW has the strongest binding ability to BRD9 among the current three inhibitors. Furthermore, the rank of the binding free energies predicted by MM-GBSA approach agrees with that determined by the experimental values. In addition, inhibitor-residue interactions were computed by using residue-based free-energy decomposition method and the results suggest that residue His42 produces the CH-H interactions, residues Asn100, Ile53 and Val49 produce the CH-[Formula: see text] interactions with three inhibitors and Tyr106, Phe45 and Phe44 generate the π-π interactions with inhibitors. Notably, the residue Asn140 forms hydrogen bonding interactions with three inhibitors. This research is expected to provide useful molecular basis and dynamics information at atomic levels for the design of potent inhibitors inhibiting the activity of BRD9.

Inclusion of methoxy groups inverts the thermodynamic stabilities of DNA-RNA hybrid duplexes: A molecular dynamics simulation study

Modified nucleic acids have found profound applications in nucleic acid based technologies such as antisense and antiviral therapies. Previous studies on chemically modified nucleic acids have suggested that modifications incorporated in furanose sugar especially at 2'-position attribute special properties to nucleic acids when compared to other modifications. 2'-O-methyl modification to deoxyribose sugars of DNA-RNA hybrids is one such modification that increases nucleic acid stability and has become an attractive class of compounds for potential antisense applications. It has been reported that modification of DNA strands with 2'-O-methyl group reverses the thermodynamic stability of DNA-RNA hybrid duplexes. Molecular dynamics simulations have been performed on two hybrid duplexes (DR and RD) which differ from each other and 2'-O-methyl modified counterparts to investigate the effect of 2'-O-methyl modification on their duplex stability. The results obtained suggest that the modification drives the conformations of both the hybrid duplexes towards A-RNA like conformation. The modified hybrid duplexes exhibit significantly contrasting dynamics and hydration patterns compared to respective parent duplexes. In line with the experimental results, the relative binding free energies suggest that the introduced modifications stabilize the less stable DR hybrid, but destabilize the more stable RD duplex. Binding free energy calculations suggest that the increased hydrophobicity is primarily responsible for the reversal of thermodynamic stability of hybrid duplexes. Free energy component analysis further provides insights into the stability of modified duplexes.