Upgrading nirmatrelvir to inhibit SARS-CoV-2 Mpro via DeepFrag and free energy calculations.

The first oral drug for the treatment of COVID-19, Paxlovid, has been authorized; however, nirmatrelvir, a major component of the drug, is reported to be associated with some side effects. Moreover, the appearance of many novel variants raises concerns about drug resistance, and designing new potent inhibitors to prevent viral replication is thus urgent. In this context, using a hybrid approach combining machine learning (ML) and free energy simulations, 6 compounds obtained by modifying nirmatrelvir were proposed to bind strongly to SARS-CoV-2 Mpro. The structural modification of nirmatrelvir significantly enhances the electrostatic interaction free energy between the protein and ligand and slightly decreases the vdW term. However, the vdW term is the most important factor in controlling the ligand-binding affinity. In addition, the modified nirmatrelvir might be less toxic to the human body than the original inhibitor.

Potential inhibitors for SARS-CoV-2 Mpro from marine compounds.

Preventing the biological activity of SARS-CoV-2 main protease using natural compounds is of great interest. In this context, using a combination of AutoDock Vina and fast pulling of ligand simulations, eleven marine fungi compounds were identified that probably play as highly potent inhibitors for preventing viral replication. In particular, four compounds including M15 (3-O-(6-O-α-l-arabinopyranosyl)-ß-d-glucopyranosyl-1,4-dimethoxyxanthone), M8 (wailupemycins H), M11 (cottoquinazolines B), and M9 (wailupemycins I) adopted the predicted ligand-binding free energy of -9.87, -9.82, -9.62, and -9.35 kcal mol-1, respectively, whereas the other adopted predicted ligand-binding free energies in the range from -8.54 to -8.94 kcal mol-1. The results were obtained using a combination of Vina and FPL simulations. Notably, although, AutoDock4 adopted higher accurate results in comparison with Vina, Vina is proven to be a more suitable technique for rapidly screening ligand-binding affinity with a large database of compounds since it requires much smaller computing resources. Furthermore, FPL is better than Vina to classify inhibitors upon ROC-AUC analysis.

Effective estimation of the inhibitor affinity of HIV-1 protease via a modified LIE approach.

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses. Computer-aided drug design (CADD) would enhance the discovery of new therapies, through which the estimation of ligand-binding affinity is critical to predict the most efficient inhibitor. A time-consuming binding free energy method would reduce the usefulness of CADD. The modified linear interaction energy (LIE) approach emerges as an appropriate protocol that performs this task. In particular, the polar interaction free energy, which is obtained via numerically resolving the linear Poisson-Boltzmann equation, plays as an important role in driving the binding mechanism of the HIV-1 PR + inhibitor complex. The electrostatic interaction energy contributes to the attraction between two molecules, but the vdW interaction acts as a repulsive factor between the ligand and the HIV-1 PR. Moreover, the ligands were found to adopt a very strong hydrophobic interaction with the HIV-1 PR. Furthermore, the results obtained corroborate the high accuracy and precision of computational studies with a large correlation coefficient value R = 0.83 and a small RMSE δ RMSE = 1.25 kcal mol-1. This method is less time-consuming than the other end-point methods, such as the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) and free energy perturbation (FEP) approaches. Overall, the modified LIE approach would provide ligand-binding affinity with HIV-1 PR accurately, precisely, and rapidly, resulting in a more efficient design of new inhibitors.

The influences of E22Q mutant on solvated 3Aß11-40 peptide: A REMD study.

The residue E22 plays a critical role in the aggregation process of Amyloid beta (Aß) peptides. The effect of E22Q mutant on the shapes of the solvated Aß11-40 trimer is clarified using a replica exchange molecular dynamics (REMD) simulation employing ∼20.6 µs of MD simulations with 48 disparate replicas. The increase of intramolecular polar contacts and salt bridge between the residue D23 to residues (24-29) was observed. The residual secondary structure of the mutated trimer is shifted in a similar way to the picture observed in previous investigations of F19W mutant. The free energy surface (FES) of the mutated E22Q system has a fewer number of minima in comparison with the wild-type trimer. The optimized shapes of the mutated E22Q form a significant increase in beta structure (47%) and serious decrease in coil content (46%) compared with the wild-type (of 36 and 56%, respectively). The binding affinity of constituting chains to the rest is of -43.7 ±â€¯6.5 kcal/mol, implying that the representative structure of E22Q is more stable than the wild-type one. Furthermore, the E22Q mutant increases the size of stable structures due to larger collision cross section (CCS) and solvent accessible area (SASA). The observed results may enhance the Aß inhibition throughout the contribution to the knowledge of the Aß oligomerization/aggregation.