Arylamines QSAR-Based Design and Molecular Dynamics of New Phenylthiophene and Benzimidazole Derivatives with Affinity for the C111, Y268, and H73 Sites of SARS-CoV-2 PLpro Enzyme

Gianfranco Sabadini,Marco Mellado,César Morales,Jaime Mella
DOI: https://doi.org/10.3390/ph17050606
IF: 4.6
2024-05-10
Pharmaceuticals
Abstract:A non-structural SARS-CoV-2 protein, PLpro, is involved in post-translational modifications in cells, allowing the evasion of antiviral immune response mechanisms. In this study, potential PLpro inhibitory drugs were designed using QSAR, molecular docking, and molecular dynamics. A combined QSAR equation with physicochemical and Free-Wilson descriptors was formulated. The r2, q2, and r2test values were 0.833, 0.770, and 0.721, respectively. From the equation, it was found that the presence of an aromatic ring and a basic nitrogen atom is crucial for obtaining good antiviral activity. Then, a series of structures for the binding sites of C111, Y268, and H73 of PLpro were created. The best compounds were found to exhibit pIC50 values of 9.124 and docking scoring values of −14 kcal/mol. The stability of the compounds in the cavities was confirmed by molecular dynamics studies. A high number of stable contacts and good interactions over time were exhibited by the aryl-thiophenes Pred14 and Pred15, making them potential antiviral candidates.
pharmacology & pharmacy,chemistry, medicinal
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to design new inhibitors of the PLpro enzyme of the SARS - CoV - 2 virus. Specifically, the researchers used quantitative structure - activity relationship (QSAR), molecular docking and molecular dynamics simulation methods to design a series of new benzothiophene and benzimidazole derivatives. These compounds can bind to the C111, Y268 and H73 sites of the PLpro enzyme, thereby inhibiting the activity of this enzyme. The main goal of the study was to discover potential drug candidates with high affinity and good antiviral activity to combat SARS - CoV - 2 virus infection. ### Research Background The non - structural protein PLpro of the SARS - CoV - 2 virus participates in post - translational modification in cells and helps the virus evade the host's antiviral immune response mechanism. Therefore, the development of effective PLpro inhibitors is considered one of the important strategies for the treatment of SARS - CoV - 2 infection. ### Research Methods 1. **QSAR Analysis**: The researchers constructed a combined QSAR equation containing physicochemical properties and Free - Wilson descriptors to predict the antiviral activity of compounds. Through cross - validation, a five - component partial least squares (PLS) model was obtained, and its internal and external validation parameters all showed good statistical performance. 2. **Molecular Docking**: Molecular docking technology was used to evaluate the binding mode and affinity of the newly designed compounds with the PLpro enzyme. 3. **Molecular Dynamics Simulation**: Molecular dynamics simulation was further used to verify the stability and interaction of the compounds at the binding sites of the PLpro enzyme. ### Main Results - **QSAR Model**: The final QSAR equation is: \[ \text{pIC}_{50} = 3.84264 - 0.23411nF - 0.79971nCl - 0.21961nO + 0.16803\text{ALogP} + 0.04589\text{ALogP}^2 + 0.00796\text{AMR} + 0.00119\text{DPSA - 3} + 0.2659n\text{Ar} + 0.17442n\text{Basic} + 0.22252n\text{Acid} + 0.07801n\text{HBDon} + 0.00225\text{TopoPSA} + 0.00145\text{VABC} \] where the meanings of each descriptor are as follows: - ALogP: Hydrophobicity - ALogP2: Square of ALogP - AMR: Molar refractivity - DPSA - 3: Difference between PPSA - 3 and PNSA - 3 - nAr: Number of aromatic rings - nBasic: Number of basic nitrogen atoms - nAcid: Number of carboxylic acid groups - nF: Number of fluorine atoms - nCl: Number of chlorine atoms - nHBDon: Number of hydrogen - bond - donor atoms - nO: Number of oxygen atoms - TopoPSA: Topological polar surface area - VABC: Van der Waals volume considering atomic and bond contributions - **Molecular Docking and Dynamics Simulation**: The binding affinities of the best compounds Pred14 and Pred15 at the C111, Y268 and H73 sites were pIC50 values of 9.124 and docking score values of - 14 kcal/mol, respectively. Molecular dynamics simulation showed that these compounds had high stability and good interactions in the binding pocket. ### Conclusion Through QSAR, molecular docking and molecular dynamics simulation, the researchers successfully designed a series of new benzothiophene and benzimidazole derivatives. These compounds have high affinity and good antiviral activity for the PLpro enzyme of the SARS - CoV - 2 virus. In particular, Pred14 and Pred15 were identified as potential antiviral drug candidates and are worthy of further experimental verification.