Design, synthesis, and molecular dynamics simulation studies of quinoline derivatives as protease inhibitors against SARS-CoV-2

A new series of quinoline derivatives has been designed and synthesized as probable protease inhibitors (PIs) against severe acute respiratory syndrome coronavirus 2. In silico studies using DS v20.1.0.19295 software have shown that these compounds behaved as PIs while interacting at the allosteric site of target Mpro enzyme (6LU7). The designed compounds have shown promising docking results, which revealed that all compounds formed hydrogen bonds with His41, His164, Glu166, Tyr54, Asp187, and showed π-interaction with His41, the highly conserved amino acids in the target protein. Toxicity Prediction by Komputer Assisted Technology results confirmed that the compounds were found to be less toxic than the reference drug. Further, molecular dynamics simulations were performed on compound 5 and remdesivir with protease enzyme. Analysis of conformational stability, residue flexibility, compactness, hydrogen bonding, solvent accessible surface area (SASA), and binding free energy revealed comparable stability of protease:5 complex to the protease: remdesivir complex. The result of hydrogen bonding showed a large number of intermolecular hydrogen bonds formed between protein residues (Glu166 and Gln189) and ligand 5, indicating strong interaction, which validated the docking result. Further, compactness analysis, SASA and interactions like hydrogen-bonding demonstrated inhibitory properties of compound 5 similar to the existing reference drug. Thus, the designed compound 5 might act as a potential inhibitor against the protease enzyme.Communicated by Ramaswamy H. SarmaHighlightsQuinoline derivatives have been designed as protease inhibitors against SARS-CoV-2.The compounds were docked at the allosteric site of SARS-CoV-2-Mpro enzyme (PDB ID: 6LU7) to study the stability of protein-ligand complex.Docking studies indicated the stable ligand-protein complexes for all designed compounds.The Toxicity Prediction by Komputer Assisted Technology protocol in DS v20.1.0.19295 software was used to evaluate the toxicity of the designed quinoline derivatives.Molecular dynamics studies indicated the formation of stable ligand-Mpro complexes.

Density functional theory calculations and molecular docking of 2-phenylbenzimidazoles with estrogen receptor for quantitative structure-activity relationship studies

Benzimidazole derivatives, especially 2-phenylbenzimidazole with various substituents on the C-5, C-2 and C-6 positions, are so important in pharmaceutical chemistry. Multiple linear regression was applied to predict the activity of 27 novel 2-phenylbenzimidazole derivatives as anticancer agents. At first, we made an effort to create a QSAR model for a selected series of novel 2-phenylbenzimidazole with density functional theory and molecular docking descriptors. Then, we tried to investigate the nature of the interactions between 2-phenylbenzimidazole derivatives and the estrogen receptor using the molecular docking method. Six descriptors of MATS4e, GATS5e, R6v, R1v+, dipole moment, and torsional free energy were selected for modelling. Due to docking results, increase in the binding energy, and decrease in the dipole moment could increase inhibitor activity.

Alkylated benzimidazoles: Design, synthesis, docking, DFT analysis, ADMET property, molecular dynamics and activity against HIV and YFV

•

New benzimidazole analogs synthesized as antivirals against HIV-1 and yellow fever virus.

•

Molecular dynamics simulation studies indicated a stable ligand-protein complex of compound 3a within NNIBP of HIV-RT.

•

DFT analysis confirmed the stability of hydrogen bonding interaction between the TRP 229 residue of HIV-RT and compound 3a.

•

Molecules were tested for their anti-HIV and broad spectrum antiviral properties against different DNA and RNA viruses.

•

Antiviral properties and cytotoxicity determined using MTT assay.

•

Compound 3a showed anti-HIV activity and compound 2b showed excellent inhibition property against yellow fever virus.

A series of alkylated benzimidazole derivatives was synthesized and screened for their anti-HIV, anti-YFV, and broad-spectrum antiviral properties. The physicochemical parameters and drug-like properties of the compounds were assessed first, and then docking studies and MD simulations on HIV-RT allosteric sites were conducted to find the possible mode of their action. DFT analysis was also performed to confirm the nature of the hydrogen bonding interaction of active compounds. The in silico studies indicated that the molecules behaved like possible NNRTIs. The nature – polar or non-polar and position of the substituent present at fifth, sixth, and N-1 positions of the benzimidazole moiety played an important role in determining the antiviral properties of the compounds. Among the various compounds, 2-(5,6-dibromo-2-chloro-1H-benzimidazol-1-yl)ethan-1-ol (3a) showed anti-HIV activity with an appreciably low IC₅₀ value as 0.386 × 10⁻⁵μM. Similarly, compound 2b, 3-(2-chloro-5-nitro-1H-benzimidazol-1-yl) propan-1-ol, showed excellent inhibitory property against the yellow fever virus (YFV) with EC₅₀ value as 0.7824 × 10⁻²μM.

An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population

Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.

Synthesis, docking, ADMET prediction, cytotoxicity and antimicrobial activity of oxathiadiazole derivatives

A series of molecules bearing oxathiadiazole, a five membered heterocyclic ring has been designed, synthesized and screened for antimicrobial activity. Molecules, 1a, 1b, 1d, 3a-b and 4a-b were found to be highly active (MIC value upto 1.5 μg/mL) against different human pathogens, namely S. aureus, B. cerus, P. aeruginosa and E. coli. Some of the compounds, 1a, 1b and 1d have also shown the antifungal activity (MIC value upto 6.2 μg/mL) against Candida albicans, Candida glubrate and Candida crusei. During in vitro cytotoxicity study, the oxathiadiazole derivatives showed less toxicity than the reference used against PBM, CEM and Vero (African green monkey kidney) cell lines. Docking studies suggested that all designed ligands interacted well within active site of PDF enzyme (PDB ID: 1G2A). Oxathiadiazole ring of all ligands formed H-bond with amino acid Leu91 at a distance ranging between 2.5-2.8 Å and also exhibited π - + and π - π interactions with amino acid residues Arg97 and His132, respectively. In silico ADMET evaluations of compounds showed more than 90% intestinal absorption for all compounds except 4b (87.45%), which too was greater than the reference drugs sulfamethoxazole (76.46%) and chloramphenicol (69.94%). TOPKAT results also supported the lower cytotoxicity of all compounds.

Design, Synthesis, and Antibacterial Activities of Novel Heterocyclic Arylsulphonamide Derivatives

Design, synthesis, and antibacterial activities of a series of arylsulphonamide derivatives as probable peptide deformylase (PDF) inhibitors have been discussed. Compounds have been designed following Lipinski's rule and after docking into the active site of PDF protein (PDB code: 1G2A) synthesized later on. Furthermore, to assess their antibacterial activity, screening of the compound was done in vitro conditions against Gram-positive and Gram-negative bacterial strains. In silico, studies revealed these compounds as potential antibacterial agents and this fact was also supported by their prominent scoring functions. Antibacterial results indicated that these molecules possessed a significant activity against Staphylococcus aureus, Bacillus cereus, Pseudomonas aeruginosa, and Escherichia coli with MIC values ranging from 0.06 to 0.29 μM. TOPKAT results showed that high LD₅₀ values and the compounds were assumed non-carcinogenic when various animal models were studied computationally.

Accessing gap-junction channel structure-function relationships through molecular modeling and simulations

Background

Gap junction channels (GJCs) are massive protein channels connecting the cytoplasm of adjacent cells. These channels allow intercellular transfer of molecules up to ~1 kDa, including water, ions and other metabolites. Unveiling structure-function relationships coded into the molecular architecture of these channels is necessary to gain insight on their vast biological function including electrical synapse, inflammation, development and tissular homeostasis. From early works, computational methods have been critical to analyze and interpret experimental observations. Upon the availability of crystallographic structures, molecular modeling and simulations have become a valuable tool to assess structure-function relationships in GJCs. Modeling different connexin isoforms, simulating the transport process, and exploring molecular variants, have provided new hypotheses and out-of-the-box approaches to the study of these important channels.

Methods

Here, we review foundational structural studies and recent developments on GJCs using molecular modeling and simulation techniques, highlighting the methods and the cross-talk with experimental evidence.

Results and discussion

By comparing results obtained by molecular modeling and simulations techniques with structural and functional information obtained from both recent literature and structural databases, we provide a critical assesment of structure-function relationships that can be obtained from the junction between theoretical and experimental evidence.

Docking and three-dimensional quantitative structure-activity relationship analyses of imidazole and thiazolidine derivatives as Aurora A kinase inhibitors

Aurora A kinase is involved in the inactivation of apoptosis leading to ovarian, breast, colon, and pancreatic cancers. Inhibitors of Aurora A kinase promote aberrant mitosis resulting in arrest at a pseudo G1 state to induce mitotic catastrophe, ultimately leading to apoptosis. In this study, ligand-based and docking-based three-dimensional quantitative structure-activity relationship (3D-QSAR) analyses of imidazole and thiazolidine derivatives as potential Aurora A kinase inhibitors were performed. The results provided highly reliable and predictive 3D-QSAR comparative molecular similarity index analysis (CoMSIA) models with a cross-validated q² value of 0.768, non-cross-validated r² value of 0.983, and predictive coefficient [Formula: see text] value of 0.978. CoMSIA contour maps suggested that the NH and benzyl hydroxy groups in R₉, and the CO group in the thiazolidine ring and pyridine ring were important components for biological activity. The maps also suggest that the introduction of hydroxy groups at C₂ of the imino-phenyl ring, C₅ in the pyridine ring, or the substitution of the imino-phenyl ring for the imino-2-pyridine ring could be applied to enhance biological activity.

Investigation of the differences in activity between hydroxycycloalkyl N1 substituted pyrazole derivatives as inhibitors of B-Raf kinase by using docking, molecular dynamics, QM/MM, and fragment-based de novo design: study of binding mode of diastereomer compounds

N1 substituted pyrazole derivatives show diverse B-Raf kinase inhibitory activities when different hydroxy-substituted cycloalkyl groups are placed at this position. Docking, molecular dynamics (MD) simulations, and hybrid calculation methods (Quantum Mechanics/Molecular Mechanics (QM/MM)) were performed on the complexes, in order to explain these differences. Docking of the inhibitors showed the same orientation that X-ray crystal structure of the analogous (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime. MD simulations of the most active diastereomer compounds containing cis- and trans-3-hydroxycyclohexyl substituents showed stable interactions with residue Ile463 at the entrance of the B-Raf active site. On the other hand, the less active diastereomer compounds containing cis- and trans-2-hydroxycyclopentyl substituents showed interactions with inner residues Asn580 and Ser465. We found that the differences in activity can be explained by considering the dynamic interactions between the inhibitors and their surrounding residues within the B-Raf binding site. We also explained the activity trend by using a testing scoring function derived from more reliable QM/MM calculations. In addition, we search for new inhibitors from a virtual screening carried out by fragment-based de novo design. We generated a set of approximately 200 virtual compounds, which interact with Ile463 and fulfill druglikeness properties according to Lipinski, Veber, and Ghose rules.

Docking and quantitative structure-activity relationship studies for 3-fluoro-4-(pyrrolo[2,1-f][1,2,4]triazin-4-yloxy)aniline, 3-fluoro-4-(1H-pyrrolo[2,3-b]pyridin-4-yloxy)aniline, and 4-(4-amino-2-fluorophenoxy)-2-pyridinylamine derivatives as c-Met kinase inhibitors

We have performed docking of 3-fluoro-4-(pyrrolo[2,1-f][1,2,4]triazin-4-yloxy)aniline (FPTA), 3-fluoro-4-(1H-pyrrolo[2,3-b]pyridin-4-yloxy)aniline (FPPA), and 4-(4-amino-2-fluorophenoxy)-2-pyridinylamine (AFPP) derivatives complexed with c-Met kinase to study the orientations and preferred active conformations of these inhibitors. The study was conducted on a selected set of 103 compounds with variations both in structure and activity. Docking helped to analyze the molecular features which contribute to a high inhibitory activity for the studied compounds. In addition, the predicted biological activities of the c-Met kinase inhibitors, measured as IC(50) values were obtained by using quantitative structure-activity relationship (QSAR) methods: Comparative molecular similarity analysis (CoMSIA) and multiple linear regression (MLR) with topological vectors. The best CoMSIA model included steric, electrostatic, hydrophobic, and hydrogen bond-donor fields; furthermore, we found a predictive model containing 2D-autocorrelation descriptors, GETAWAY descriptors (GETAWAY: Geometry, Topology and Atom-Weight AssemblY), fragment-based polar surface area (PSA), and MlogP. The statistical parameters: cross-validate correlation coefficient and the fitted correlation coefficient, validated the quality of the obtained predictive models for 76 compounds. Additionally, these models predicted adequately 25 compounds that were not included in the training set.

Insights into the structure of urea-like compounds as inhibitors of the juvenile hormone epoxide hydrolase (JHEH) of the tobacco hornworm Manduca sexta: analysis of the binding modes and structure-activity relationships of the inhibitors by docking and CoMFA calculations

Substituted urea compounds are well-known as potent inhibitors of juvenile hormone epoxide hydrolase (JHEH) of the tobacco hornworm Manduca sexta. Docking simulations of 47 derivatives inside JHEH were performed to gain insight into the structural characteristics of these complexes. The obtained orientations show a strong similitude with the observed in the known X-ray crystal structures of human soluble epoxide hydrolase (sEH) complexed with dialkylurea inhibitors. In addition, the predicted inhibitor concentration (IC₅₀) of the above-mentioned compounds as JHEH inhibitors were obtained by a quantitative structure-activity relationship (QSAR) method by using comparative molecular field analysis (CoMFA) applied to aligned dataset. The best models included steric and electrostatic fields and had adequate predictive abilities. In addition, these models were used to predict the activity of an external test set of compounds that was not used for building the model. Furthermore, plots of the CoMFA fields allowed conclusions to be drawn for the choice of suitable inhibitors.