Abdul M Gbaj1*, Nisreen H Meiqal1, Inass A Sadawe1, Salah M Bensaber1, Abdulathim A A Alshoushan2, Massaud Salem Maamar3 and Anton Hermann4
Received: March 29, 2020; Published: June 16, 2020
*Corresponding author: Abdul M Gbaj, Associate Professor of Genetics and Biochemistry, Department of Medicinal Chemistry, Faculty of Pharmacy, University of Tripoli, LibyaK
A series of new ortho-phenylenediamine derivatives were designed. The crystal structure of the main protease monomer was used as a target protein for molecular docking of ortho- phenylenediamine derivatives and a protein-ligand interaction analysis was performed using Auto Dock 4.2 software. Based on the docking score and after additional three-dimensional similarity analysis, NHM7 [(10,10’-((1E,1’E)-(1,2-phenylenebis(azanylylidene))bis(methanylylidene))bis(anthracen-9(8aH)-one)] had the highest binding energy. The calculated binding energy of ortho- phenylenediamine indicated effective binding of the proposed inhibitors to COVID-19 proteinase.
Keywords: Coronavirus; severe acute respiratory syndrome; covid-19 protease inhibitors; lopinavir; n3 inhibitor; molecular docking
Coronaviruses (CoVs) comprise families of viruses pertained to the family of Coronaviridae. The viruses can circulate in humans and cause serious infections antinutrient the respiratory system [1,2]. Acute respiratory syndrome coronavirus (SARS-CoV) revealed that it can cause severe andoccasionallydeadly respiratory tract infections in humans [3-5]. The consecutive outbreaks, in addition, emphasize the threat of these viruses and caused a pandemic warning that has been declared a public health crisis of international anxiety [5-7].There are many potential targets against COVID-19 and among targets replication-related enzymes, such as protease3CL(pro) [also called SARS-CoV 3CL(pro)] are extremely conserved [8-10]. It has been reported that drugs that inhibit proteases are able preventing proliferation and replication of the virus by impeding with post-translational processing of vital viral polypeptides. In addition, they could also decrease the risk of drug-resistance produced by mutation. Following the protease inhibition approach two standard protease inhibitors were used as lead including Lopinavir and N3 (Figure 1) inhibitors identified of being capable to inhibit SARS-Covian protease . It has been reported that the SARS-CoV main protease has 96.1% homology with the COVID-19main protease andthereforemay be utilized as a homologous target for screening of ortho- phenylenediamine derivatives that could inhibit the proliferation and replication of COVID-19 [12-15]. The ortho- phenylenediamine derivatives are Schiff bases recognized for their therapeutic value as they were reported to have anti-inflammatory, analgesic, antiviral, antitumor, antifungal and antibacterial properties [16-21]. Molecular modeling is a recognized computational toolto aid early drug discovery and development. It is used to generate ideasof a compounds or macromolecules 3D conformation, protein–ligand interactions, and allows forecasts about biological activities. The integration of molecular modeling in drug or vaccine design can help in early drug or vaccine discoveries [22-24]. The main aim of this study is to further identify protease as a target, and by computational drug repurposing procedures to allocate appropriate inhibitory agents.
The starting geometry of the ortho-phenylenediamine derivatives was constructed using chem3D Ultra software (version 8.0, Cambridge soft Com., USA). The optimized geometry of orthophenylenediamine derivatives with the lowest energy was used for molecular dockings. The crystal structure of COVID-19 main protease in complex with the inhibitor N3 (6LU7) was downloaded from the Protein Data Bank http://www.rcsb.org/structure/6LU7 and COVID-19 and the main protease in apo form (6M03) was downloaded from the Protein Data Bank http://www.rcsb.org/ structure/6m03. Molecular dockings of ortho- phenylenediamine derivatives with 6LU7and 6M03 was accomplished by Auto Dock 4.2 software from the Scripps Research Institute (TSRI) ( http:// autodock.scripps.edu/). Firstly, polar hydrogen atoms were added into protein molecules. Then, partial atomic charges of the protease enzymes and ortho- phenylenediamine derivatives molecules were calculated using Kalman methods . In the process of molecular docking, the grid maps of dimensions: (60Å X 60Å X 60Å) and (36.8Å X 64.6Å X 60Å) for 6LU7 and 6M03, respectively, with a gridpoint spacing of 0.376Å and the grid boxes is centered. The number of genetic algorithms runs and the number of evaluations were set to 100. All other parameters were default settings. Cluster analysis was performed on the basis of docking results by using a root mean square (RMS) tolerance of 2.0Å, dependent on the binding free energy. Lastly, the dominating configuration of the binding complex of ortho- phenylenediamine derivatives and protease enzymes fragments with minimum energy of binding were determined which relied strongly on the information of 3D-structures of the protease binding site and ultimately generated a series of proteasebinding complexes.
Table 1 shows the binding energies of Lopinavir, N3 (Figure 1, as standards), ortho-phenylenediamine derivatives, and protease enzymes (6LU7 and 6M03) obtained by the molecular docking strategy. Molecular dockings of the ortho-phenylenediamine derivatives with protease enzymes (6LU7 and 6M03) were performed using Auto Dock 4.2 to obtain information about interaction forces between ortho phenylenediamine derivatives and protease enzymes (6LU7 and 6M03). ortho-phenylenediamine derivatives and protease enzymes (6LU7 and 6M03) were kept as flexible molecules and were docked into seven forms of rigid protease enzymes (6LU7 and 6M03) to obtain the preferential binding site to ortho-phenylenediamine derivatives on protease enzymes (6LU7 and 6M03). The molecular docking results are shown in Table 1. The modeling studies indicate an der Waals, hydrogen bonding (Table 1) and electrostatic interactions between ortho-phenylenediamine derivatives with protease enzymes (6LU7 and 6M03). The contribution of van der Waals and hydrogen bonding interaction is much greater than that of the electrostatic interaction because the sum of van der Waals energy, hydrogen bonding energy and desolation free energy is larger than the electrostatic energy, [26,27].The ortho-phenylenediamine derivatives, and protease enzymes (6LU7 and 6M03) interactions are shown in Figure 2. Ortho phenylenediamine derivatives provide higher binding energy (-8.1 to -11.0kcal/mol) compared to standard 6LU7 and 6M03 (-7.0 to -7.9 kcal/mol) (Table 1 & Figure 2) indicates four hydrogen bonds between NHM7 and 6LU7.In addition, NHM7 showed good docking interaction of -11.0 kcal/mol with the 6LU7 binding site (GLU166, VAL3, GLU166and LEU4) (Figure 2). Compound NHM7 has the highest binding energy of the series. This compound has an extra phenyl moiety attached to the naphthyl analogue of the phenylenediamine Schiff`s base derivative with a log P value of 7.49 indicating the importance of the lipophilicity for the interaction with the active site. The interaction of similar Schiff’s base ortho phenylenediamine derivatives with the proteases binding site of the enzyme is essential for effective inhibition as previously reported [28-31]. Therefore, ortho-phenylenediamine derivatives may be considered the most effective NHM7 and 6LU7 proteases inhibitors. The obtained results using computational drug repurposing is an efficient way to find novel applications for already known drugs . Molecular docking and binding free energy calculations for ortho- phenylenediamine derivatives can be used to forecast drugtarget interactions and binding affinity (Figure 3). The appearance of resistance to existing antiviral drugs or vaccines is a major challenge in antiviral drug development. The drug repurposing technique allows finding novel antiviral agents within a short period in order to overcome the challenges in antiviral therapy. Computational drug repurposing has previously been used to recognize drug candidates for viral infectious diseases like ZIKA, Ebola, influenza and dengue infections. These methods were also utilized to recognize possible drugs against MERS-CoV and SARSCoV [33,34] and following the COVID-19 outbreak, computational repurposing has been and are used for COVID-19.
Figure 2: Interaction model between NHM7 with COVID-19 main protease (6LU7) active site. NHM7 l green colour. Hydrogen bonds green broken line.
Table 1: Various energies in the binding process of ortho-phenylenediamine derivatives, N3 and Lopinavir with COVID-19 protease enzymes (6LU7, 6M03) obtained from molecular docking. The unit of all energies (ΔG) is kcal/mol.
In spite of the economic and societal shock of COVID-19 infections and the probability of future outbreaks of even more stern pathogenic COVID-19 in humans, there is still a lack of efficient antiviral strategies to treat COVID-19 and only few options are available to prevent COVID-19 infections. Rapid development and use of a broad-spectrum protease inhibitor alone or in combination with other potent inhibitors of proteases might fill the therapeutic gap spanning quarantine and hospital setting. Further elaborative work is necessary for better understanding the mechanisms of protease inhibition. According to modeling studies ortho-phenylenediamine derivatives may have the ability to inhibit COVID-19 proteases making them reasonable candidates for consideration of clinical trials and warrant further examination. Results presented in this study shall motivate future efforts in finding potent ortho-phenylenediamine derivatives that can be used for COVID-19 protease inhibition in vivo.
The authors gratefully acknowledge the technical support and valuable suggestions obtained from MS Amira Abdul Gbaj (Novelien Zone, Tripoli, Libya).
Bio chemistryUniversity of Texas Medical Branch, USA
Department of Criminal JusticeLiberty University, USA
Department of PsychiatryUniversity of Kentucky, USA
Department of MedicineGally International Biomedical Research & Consulting LLC, USA
Department of Urbanisation and AgriculturalMontreal university, USA
Oral & Maxillofacial PathologyNew York University, USA
Gastroenterology and HepatologyUniversity of Alabama, UK
Department of MedicineUniversities of Bradford, UK
OncologyCirculogene Theranostics, England
Radiation ChemistryNational University of Mexico, USA
Analytical ChemistryWentworth Institute of Technology, USA
Minimally Invasive SurgeryMercer University school of Medicine, USA
Pediatric DentistryUniversity of Athens , Greece
The annual scholar awards from Lupine Publishers honor a selected number Read More...