Structure Of Biomolecules Through Molecular Dynamics Simulations

In Silico Strategies for Designing of Peptide Inhibitors of Oncogenic K-Ras G12V Mutant: Inhibiting Cancer Growth and Proliferation

Ras plays a pivotal function in cell proliferation and is an important protein in signal transduction pathways. Mutations in genes encoding the Ras protein drive the signaling cascades essential for malignant transformation, tumour angiogenesis, and metastasis and are responsible for above 30% of all human cancers. There is evidence that N-Ras, K-Ras, and H-Ras play significant roles in human cancer. The mutated K-Ras protein is typically observed in malignant growths. Mutant K-Ras is the most common in lung, colon, and pancreatic cancers. The purpose of this research was to create peptides that inhibit K-Ras G12V. The crystal structure of the mutant K-Ras G12V-H-REV107 complex was obtained from a protein data bank. Further, we used a residue scan approach to create unique peptides from the reference peptide (H-REV107). AMBER molecular dynamics simulations were used to test the stability of the top four proposed peptides (based on binding free energies). Our findings showed that the top four selected peptides had stronger interactions with K-Ras than the reference peptide and have the ability to block the activation function of K-Ras. Our extensive analyses of binding affinities showed that our designed peptide possesses the potential to inhibit K-Ras and to reduce the progression of cancer.

Consensus docking and MM-PBSA computations identify putative furin protease inhibitors for developing potential therapeutics against COVID-19

The coronavirus disease 2019 (COVID-19) is a pandemic that has severely posed substantial health challenges and claimed millions of lives. Though vaccines have been produced to stem the spread of this disease, the death rate remains high since drugs used for treatment have therapeutic challenges. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes the disease, has a slew of potential therapeutic targets. Among them is the furin protease, which has a cleavage site on the virus’s spike protein. The cleavage site facilitates the entry of the virus into human cells via cell–cell fusion. This critical involvement of furin in the disease pathogenicity has made it a viable therapeutic strategy against the virus. This study employs the consensus docking approach using HYBRID and AutoDock Vina to virtually screen a pre-filtered library of 3942 natural product compounds of African origin against the human furin protease (PDB: 4RYD). Twenty of these compounds were selected as hits after meeting molecular docking cut-off of − 7 kcal.mol⁻¹, pose alignment inspection, and having favorable furin-ligand interactions. An area under the curve (AUC) value of 0.72 was computed from the receiver operator characteristic (ROC) curve, and Boltzmann-enhanced discrimination of the ROC curve (BEDROC) value of 0.65 showed that AutoDock Vina was a reasonable tool for selecting actives for this target. Seven of these hits were proposed as potential leads having had bonding interactions with catalytic triad residues Ser368, His194, and Asp153, and other essential residues in the active site with plausible binding free energies between − 189 and − 95 kJ/mol from the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations as well as favorable ADME/Tox properties. The molecules were also predicted as antiviral, anti-inflammatory, membrane permeability inhibitors, RNA synthesis inhibitors, cytoprotective, and hepatoprotective with probable activity (Pa) above 0.5 and probable inactivity values below 0.1. Some of them also have anti-influenza activity. Influenza virus has many similarities with SARS-CoV-2 in their mode of entry into human cells as both are facilitated by the furin protease. Pinobanksin 3-(E)-caffeate, one of the potential leads is a propolis compound. Propolis compounds have shown inhibitory effects against ACE2, TMPRSS2, and PAK1 signaling pathways of SARS-CoV-2 in previous studies. Likewise, quercitrin is structurally similar to isoquercetin, which is currently in clinical trials as possible medication for COVID-19.

Multi-spectroscopic, thermodynamic, and molecular docking/dynamic approaches for characterization of the binding interaction between calf thymus DNA and palbociclib

Studying the binding interaction between biological macromolecules and small molecules has formed the core of different research aspects. The interaction of palbociclib with calf thymus DNA at simulated physiological conditions (pH 7.4) was studied using different approaches, including spectrophotometry, spectrofluorimetry, FT-IR spectroscopy, viscosity measurements, ionic strength measurements, thermodynamic, molecular dynamic simulation, and docking studies. The obtained findings showed an apparent binding interaction between palbociclib and calf thymus DNA. Groove binding mode was confirmed from the findings of competitive binding studies with ethidium bromide or rhodamine B, UV–Vis spectrophotometry, and viscosity assessment. The binding constant (K_b) at 298 K calculated from the Benesi–Hildebrand equation was found to be 6.42 × 10³ M⁻¹. The enthalpy and entropy changes (∆H⁰ and ∆S⁰) were − 33.09 kJ mol⁻¹ and 61.78 J mol⁻¹ K⁻¹, respectively, showing that hydrophobic and hydrogen bonds constitute the primary binding forces. As indicated by the molecular docking results, palbociclib fits into the AT-rich region of the B-DNA minor groove with four base pairs long binding site. The dynamic performance and stability of the formed complex were also evaluated using molecular dynamic simulation studies. The in vitro study of the intermolecular binding interaction of palbociclib with calf thymus DNA could guide future clinical and pharmacological studies for the rational drug scheming with enhanced or more selective activity and greater efficacy.

Computational Analysis of Deleterious SNPs in NRAS to Assess Their Potential Correlation With Carcinogenesis

The NRAS gene is a well-known oncogene that acts as a major player in carcinogenesis. Mutations in the NRAS gene have been linked to multiple types of human tumors. Therefore, the identification of the most deleterious single nucleotide polymorphisms (SNPs) in the NRAS gene is necessary to understand the key factors of tumor pathogenesis and therapy. We aimed to retrieve NRAS missense SNPs and analyze them comprehensively using sequence and structure approaches to determine the most deleterious SNPs that could increase the risk of carcinogenesis. We also adopted structural biology methods and docking tools to investigate the behavior of the filtered SNPs. After retrieving missense SNPs and analyzing them using six in silico tools, 17 mutations were found to be the most deleterious mutations in NRAS. All SNPs except S145L were found to decrease NRAS stability, and all SNPs were found on highly conserved residues and important functional domains, except R164C. In addition, all mutations except G60E and S145L showed a higher binding affinity to GTP, implicating an increase in malignancy tendency. As a consequence, all other 14 mutations were expected to increase the risk of carcinogenesis, with 5 mutations (G13R, G13C, G13V, P34R, and V152F) expected to have the highest risk. Thermodynamic stability was ensured for these SNP models through molecular dynamics simulation based on trajectory analysis. Free binding affinity toward the natural substrate, GTP, was higher for these models as compared to the native NRAS protein. The Gly13 SNP proteins depict a differential conformational state that could favor nucleotide exchange and catalytic potentiality. A further application of experimental methods with all these 14 mutations could reveal new insights into the pathogenesis and management of different types of tumors.

Bio-Guided Isolation of SARS-CoV-2 Main Protease Inhibitors from Medicinal Plants: In Vitro Assay and Molecular Dynamics

Since the emergence of the pandemic of the coronavirus disease (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the discovery of antiviral phytoconstituents from medicinal plants against SARS-CoV-2 has been comprehensively researched. In this study, thirty-three plants belonging to seventeen different families used traditionally in Saudi Arabia were tested in vitro for their ability to inhibit the SARS-CoV-2 main protease (M^PRO). Major constituents of the bio-active extracts were isolated and tested for their inhibition potential against this enzyme; in addition, their antiviral activity against the SARS-CoV-2 Egyptian strain was assessed. Further, the thermodynamic stability of the best active compounds was studied through focused comparative insights for the active metabolites regarding ligand–target binding characteristics at the molecular level. Additionally, the obtained computational findings provided useful directions for future drug optimization and development. The results revealed that Psiadia punctulata, Aframomum melegueta, and Nigella sativa extracts showed a high percentage of inhibition of 66.4, 58.7, and 31.5%, against SARS-CoV-2 M^PRO, respectively. The major isolated constituents of these plants were identified as gardenins A and B (from P. punctulata), 6-gingerol and 6-paradol (from A. melegueta), and thymoquinone (from N. sativa). These compounds are the first to be tested invitro against SARS-CoV-2 M^PRO. Among the isolated compounds, only thymoquinone (THY), gardenin A (GDA), 6-gingerol (GNG), and 6-paradol (PAD) inhibited the SARS-CoV-2 M^PRO enzyme with inhibition percentages of 63.21, 73.80, 65.2, and 71.8%, respectively. In vitro assessment of SARS-CoV-2 (hCoV-19/Egypt/NRC-03/2020 (accession number on GSAID: EPI_ISL_430820) revealed a strong-to-low antiviral activity of the isolated compounds. THY showed relatively high cytotoxicity and was anti-SARS-CoV-2, while PAD demonstrated a cytotoxic effect on the tested VERO cells with a selectivity index of CC₅₀/IC₅₀ = 1.33 and CC₅₀/IC₅₀ = 0.6, respectively. Moreover, GNG had moderate activity at non-cytotoxic concentrations in vitro with a selectivity index of CC₅₀/IC₅₀ = 101.3/43.45 = 2.3. Meanwhile, GDA showed weak activity with a selectivity index of CC₅₀/IC₅₀ = 246.5/83.77 = 2.9. The thermodynamic stability of top-active compounds revealed preferential stability and SARS-CoV-2 M^PRO binding affinity for PAD through molecular-docking-coupled molecular dynamics simulation. The obtained results suggest the treating potential of these plants and/or their active metabolites for COVID-19. However, further in-vivo and clinical investigations are required to establish the potential preventive and treatment effectiveness of these plants and/or their bio-active compounds in COVID-19.

Mining of Marburg Virus Proteome for Designing an Epitope-Based Vaccine

Marburg virus (MARV) is one of the most harmful zoonotic viruses with deadly effects on both humans and nonhuman primates. Because of its severe outbreaks with a high rate of fatality, the world health organization put it as a risk group 4 pathogen and focused on the urgent need for the development of effective solutions against that virus. However, up to date, there is no effective vaccine against MARV in the market. In the current study, the complete proteome of MARV (seven proteins) was analyzed for the antigenicity score and the virulence or physiological role of each protein where we nominated envelope glycoprotein (Gp), Transcriptional activator (VP30), and membrane-associated protein (VP24) as the candidates for epitope prediction. Following that, a vaccine construct was designed based on CTL, HTL, and BCL epitopes of the selected protein candidates and to finalize the vaccine construct, several amino acid linkers, β-defensin adjuvant, and PADRE peptides were incorporated. The generated potential vaccine was assessed computationally for several properties such as antigenicity, allergenicity, stability, and other structural features where the outcomes of these assessments nominated this potential vaccine to be validated for its binding affinity with two molecular targets TLR-8 and TLR-4. The binding score and the stability of the vaccine-receptor complex, which was deeply studied through molecular docking-coupled dynamics simulation, supported the selection of our designed vaccine as a putative solution for MARV that should be validated through future wet-lab experiments. Here, we describe the computational approach for designing and analysis of this potential vaccine.

Rational Discovery of Antimicrobial Peptides by Means of Artificial Intelligence

Antibiotic resistance is a worldwide public health problem due to the costs and mortality rates it generates. However, the large pharmaceutical industries have stopped searching for new antibiotics because of their low profitability, given the rapid replacement rates imposed by the increasingly observed resistance acquired by microorganisms. Alternatively, antimicrobial peptides (AMPs) have emerged as potent molecules with a much lower rate of resistance generation. The discovery of these peptides is carried out through extensive in vitro screenings of either rational or non-rational libraries. These processes are tedious and expensive and generate only a few AMP candidates, most of which fail to show the required activity and physicochemical properties for practical applications. This work proposes implementing an artificial intelligence algorithm to reduce the required experimentation and increase the efficiency of high-activity AMP discovery. Our deep learning (DL) model, called AMPs-Net, outperforms the state-of-the-art method by 8.8% in average precision. Furthermore, it is highly accurate to predict the antibacterial and antiviral capacity of a large number of AMPs. Our search led to identifying two unreported antimicrobial motifs and two novel antimicrobial peptides related to them. Moreover, by coupling DL with molecular dynamics (MD) simulations, we were able to find a multifunctional peptide with promising therapeutic effects. Our work validates our previously proposed pipeline for a more efficient rational discovery of novel AMPs.

Engineered Design of the E-Helix Structure on Ferritin Nanoparticles

Insertion of an immunogenic epitope at the C-terminus of ferritin has shown the potential to produce a stable and efficacious vaccine. There is however limited understanding of how C-terminus insertion affects ferritin protein stability. The E-helix at the C-terminus has attracted interest because there are contradictory reports as to whether it has a role in protein stabilization. Here, we report, for the first time, combining molecular dynamics simulation (MDS) with experiment to engineer the design of the E-helix at the C-terminus of engineered human ferritin heavy chain (F₁) inserted with Epstein-Barr nuclear antigen 1 (EBNA1, E₁) and flexible linker (L₃) residues (to afford F₁L₃E₁). Hot spots on the E-helix of the C-terminus were predicted by MDS at aa 167 (Glu) and aa 171 (Asp). Five (5) variants of F₁L₃E₁ were constructed by considering hot spots and alteration of electrostatic or hydrophobic interfaces, namely, (1) C1, hot spots substituted with noncharged residue Gln; (2) C2, hot spots substituted with positively charged residue Arg; (3) C3, hydrophobic residues substituted with the most hydrophobic residues Val and Ile; (4) C4, hydrophobic residues substituted with the most hydrophilic residues Gln and Asn; and (5) C5, a heptad repeat structure in the E-helix disrupted by substituting "a" and "d" heptad residues with noncharged polar residue Gln. It was found that the E-helix is essential to maintain integrated protein stability and that changing the hydrophobic interface (C3 and C4) had more significant effects on protein folding and stability than changing the electrostatic interface (C1 and C2). It was confirmed by both MDS and experiment that variants C1, C2, and C5 were able to fold to form stable conformational structures with protein surface hydrophobicity similar to that of F₁L₃E₁. However, they are less thermally stable than F₁L₃E₁. Significant changes in hydrophobicity drove significant protein aggregation for variants C3 and C4. It is concluded that the molecular design of the C-terminus in engineered ferritin, especially the E-helix, is important to ensure the epitope-based chimeric vaccine is safe (aggregate free) and efficacious.

Translocating Peptides of Biomedical Interest Obtained from the Spike (S) Glycoprotein of the SARS-CoV-2

At the beginning of 2020, the pandemic caused by the SARS-CoV-2 virus led to the fast sequencing of its genome to facilitate molecular engineering strategies to control the pathogen’s spread. The spike (S) glycoprotein has been identified as the leading therapeutic agent due to its role in localizing the ACE2 receptor in the host’s pulmonary cell membrane, binding, and eventually infecting the cells. Due to the difficulty of delivering bioactive molecules to the intracellular space, we hypothesized that the S protein could serve as a source of membrane translocating peptides. AHB-1, AHB-2, and AHB-3 peptides were identified and analyzed on a membrane model of DPPC (dipalmitoylphosphatidylcholine) using molecular dynamics (MD) simulations. An umbrella sampling approach was used to quantify the energy barrier necessary to cross the boundary (13.2 to 34.9 kcal/mol), and a flat-bottom pulling helped to gain a deeper understanding of the membrane’s permeation dynamics. Our studies revealed that the novel peptide AHB-1 exhibited comparable penetration potential of already known potent cell-penetrating peptides (CPPs) such as TP2, Buforin II, and Frenatin 2.3s. Results were confirmed by in vitro analysis of the peptides conjugated to chitosan nanoparticles, demonstrating its ability to reach the cytosol and escape endosomes, while maintaining high biocompatibility levels according to standardized assays.

Structure-based identification of potential SARS-CoV-2 main protease inhibitors

To address coronavirus disease (COVID-19), currently, no effective drug or vaccine is available. In this regard, molecular modeling approaches are highly useful to discover potential inhibitors of the main protease (Mpro) enzyme of SARS-CoV-2. Since, the Mpro enzyme plays key roles in mediating viral replication and transcription; therefore, it is considered as an attractive drug target to control SARS-CoV-2 infection. By using structure-based drug design, pharmacophore modeling, and virtual high throughput drug screening combined with docking and all-atom molecular dynamics simulation approach, we have identified five potential inhibitors of SARS-CoV-2 Mpro. MD simulation studies revealed that compound 54035018 binds to the Mpro with high affinity (ΔGbind -37.40 kcal/mol), and the complex is more stable in comparison with other protein-ligand complexes. We have identified promising leads to fight COVID-19 infection as these compounds fulfill all drug-likeness properties. However, experimental and clinical validations are required for COVID-19 therapy.Communicated by Ramaswamy H. Sarma.

Soyasapogenol-B as a Potential Multitarget Therapeutic Agent for Neurodegenerative Disorders: Molecular Docking and Dynamics Study

Neurodegenerative disorders involve various pathophysiological pathways, and finding a solution for these issues is still an uphill task for the scientific community. In the present study, a combination of molecular docking and dynamics approaches was applied to target different pathways leading to neurodegenerative disorders such as Alzheimer’s disease. Initially, abrineurin natural inducers were screened using physicochemical properties and toxicity assessment. Out of five screened compounds, a pentacyclic triterpenoid, i.e., Soyasapogenol B appeared to be the most promising after molecular docking and simulation analysis. Soyasapogenol B showed low TPSA (60.69), high absorption (82.6%), no Lipinski rule violation, and no toxicity. Docking interaction analysis revealed that Soyasapogenol B bound effectively to all of the targeted proteins (AChE, BuChE MAO-A, MAO-B, GSK3β, and NMDA), in contrast to other screened abrineurin natural inducers and inhibitors. Importantly, Soyasapogenol B bound to active site residues of the targeted proteins in a similar pattern to the native ligand inhibitor. Further, 100 ns molecular dynamics simulations analysis showed that Soyasapogenol B formed stable complexes against all of the targeted proteins. RMSD analysis showed that the Soyasapogenol B–protein complex exhibited average RMSD values of 1.94 Å, 2.11 Å, 5.07 Å, 2.56 Å, 3.83 Å and 4.07 Å. Furthermore, the RMSF analysis and secondary structure analysis also indicated the stability of the Soyasapogenol B–protein complexes.

Identification of the most damaging nsSNPs in the human CFL1 gene and their functional and structural impacts on cofilin-1 protein

The cofilin-1 protein, encoded by CFL1, is an actin-binding protein that regulates F-actin depolymerization and nucleation activity through phosphorylation and dephosphorylation. CFL1 has been implicated in the development of neurodegenerative diseases (Alzheimer's disease and Huntington's disease), neuronal migration disorders (lissencephaly, epilepsy, and schizophrenia), and neural tube closure defects. Mutations in CFL1 have been associated with impaired neural crest cell migration and neural tube closure defects. In our study, various computational approaches were utilized to explore single-nucleotide polymorphisms (SNPs) in CFL1. The Variation Viewer and gnomAD databases were used to retrieve CFL1 SNPs, including 46 nonsynonymous SNPs (nsSNPs). The functional and structural annotation of SNPs was performed using 12 sequence-based web applications, which identified 20 nsSNPs as being the most likely to be deleterious or disease-causing. The conservation of cofilin-1 protein structures was illustrated using the ConSurf and PROSITE web servers, which projected the 12 most deleterious nsSNPs onto conserved domains, with the potential to disrupt the protein's functionality. These 12 nsSNPs were selected for protein structure construction, and the DynaMut/DUET servers predicted that the protein variants V7G, L84P, and L99A were the most likely to be damaging to the cofilin-1 protein structure or function. The evaluation of molecular docking studies demonstrated that the L99A and L84P cofilin-1 variants reduce the binding affinity for actin compared with the native cofilin-1 structure, and molecular dynamic simulation studies confirmed that these variants might destabilize the protein structure. The consequences of putative mutations on protein-protein interactions and post-translational modification sites in the cofilin-1 protein structure were analyzed. This study represents the first complete approach to understanding the effects of nsSNPs within the actin-depolymerizing factor/cofilin family, which suggested that SNPs resulting in L84P (rs199716082) and L99A (rs267603119) variants represent significant CFL1 mutations associated with disease development.

Identification of acridinedione scaffolds as potential inhibitor of DENV-2 C protein: An in silico strategy to combat dengue

Dengue is a prominent viral disease transmitted by mosquitoes to humans that affects mainly tropical and subtropical countries worldwide. The global spread of dengue virus (DENV) is mainly occurred by Aedes aegypti and Aedes albopictus mosquitoes. The dengue virus serotypes-2 (DENV-2) is a widely prevalent serotype of DENV, that causes the hemorrhagic fever and bleeding in the mucosa, which can be fatal. In the life cycle of DENV-2, a structural capsid (DENV-2 C) protein forms the nucleocapsid assembly and bind to the viral progeny RNA. For DENV-2 maturation, the nucleocapsid is a vital component. We used virtual ligand screening to filter out the best in-house synthesized acridinedione analogs (DSPD molecules) that could efficiently bind to DENV-2 C protein. The molecular docking and dynamics simulations studies were performed to analyze the effect of DSPD molecules on DENV-2 C protein after binding. Our findings showed that DSPD molecules strongly interacted with DENV-2 C protein, as evident from molecular interactions and several time-dependent molecular dynamics-driven analyses. Moreover, this study was also supported by the thermodynamic binding free energy and steered molecular dynamics simulations. Therefore, we intend to suggest that the DSPD3 molecule could be used as a potential therapeutic molecule against dengue complications as compared to the cocrystallized inhibitor ST-148. However, further studies are required to demonstrate the ability of DSPD3 to induce DENV-2 C tetramer formation.

In silico Designing of an Epitope-Based Vaccine Against Common E. coli Pathotypes

Escherichia coli (E. coli) is a Gram-negative bacterium that belongs to the family Enterobacteriaceae. While E. coli can stay as an innocuous resident in the digestive tract, it can cause a group of symptoms ranging from diarrhea to live threatening complications. Due to the increased rate of antibiotic resistance worldwide, the development of an effective vaccine against E. coli pathotypes is a major health priority. In this study, a reverse vaccinology approach along with immunoinformatics has been applied for the detection of potential antigens to develop an effective vaccine. Based on our screening of 5,155 proteins, we identified lipopolysaccharide assembly protein (LptD) and outer membrane protein assembly factor (BamA) as vaccine candidates for the current study. The conservancy of these proteins in the main E. coli pathotypes was assessed through BLASTp to make sure that the designed vaccine will be protective against major E. coli pathotypes. The multitope vaccine was constructed using cytotoxic T lymphocyte (CTL), helper T lymphocyte (HTL), and B cell lymphocyte (BCL) epitopes with suitable linkers and adjuvant. Following that, it was analyzed computationally where it was found to be antigenic, soluble, stable, and non-allergen. Additionally, the adopted docking study, as well as all-atom molecular dynamics simulation, illustrated the promising predicted affinity and free binding energy of this constructed vaccine against the human Toll-like receptor-4 (hTLR-4) dimeric state. In this regard, wet lab studies are required to prove the efficacy of the potential vaccine construct that demonstrated promising results through computational validation.

Dissociation Pathways of the p53 DNA Binding Domain from DNA and Critical Roles of Key Residues Elucidated by dPaCS-MD/MSM

p53 is a transcriptional factor that regulates cell response to a variety of stresses. About a half of all human tumors contain p53 mutations, and the accumulation of mutations in the DNA binding domain of p53 (p53-DBD) can cause destabilization of p53 and its complex with DNA. To identify the key residues of the p53-DBD/DNA binding and to understand the dissociation mechanisms of the p53-DBD/DNA complex, the dissociation process of p53-DBD from a DNA duplex that contains the consensus sequence (the specific target of p53-DBD) was investigated by a combination of dissociation parallel cascade selection molecular dynamics (dPaCS-MD) and the Markov state model (MSM). This combination (dPaCS-MD/MSM) enabled us to simulate dissociation of the two large molecules based on an all-atom model with a short simulation time (11.2 ± 2.2 ns per trial) and to analyze dissociation pathways, free energy landscape (FEL), and binding free energy. Among 75 trials of dPaCS-MD, p53-DBD dissociated first from the major groove and then detached from the minor groove in 93% of the cases, while 7% of the cases unbinding from the minor groove occurred first. Minor groove binding is mainly stabilized by R248, identified as the most important residue that tightly binds deep inside the minor groove. The standard binding free energy calculated from the FEL was -10.9 ± 0.4 kcal/mol, which agrees with an experimental value of -11.1 kcal/mol. These results indicate that the dPaCS-MD/MSM combination can be a powerful tool to investigate dissociation mechanisms of two large molecules. Analysis of the p53 key residues for DNA binding indicates high correlations with cancer-related mutations, confirming that impairment of the interactions between p53-DBD and DNA can be frequently related to cancer.

Repurposing α-Adrenoreceptor Blockers as Promising Anti-Virulence Agents in Gram-Negative Bacteria

Antimicrobial resistance is among the world’s most urgent public health problems. Diminishing of the virulence of bacteria is a promising approach to decrease the development of bacterial resistance. Quorum sensing (QS) systems orchestrate the bacterial virulence in inducer–receptors manner. Bacteria can spy on the cells of the host by sensing adrenergic hormones and other neurotransmitters, and in turn, these neurotransmitters can induce bacterial pathogenesis. In this direction, α-adrenergic blockers were proposed as an anti-virulence agents through inhibiting the bacterial espionage. The current study aimed to explore the α-blockers’ anti-QS activities. Within comprehensive in silico investigation, the binding affinities of seven α-adrenoreceptor blockers were evaluated towards structurally different QS receptors. From the best docked α-blockers into QS receptors, terazosin was nominated to be subjected for further in vivo and in vitro anti-QS and anti-virulence activities against Chromobacterium violaceum and Pseudomonas aeruginosa. Terazosin showed a significant ability to diminish the QS-controlled pigment production in C. violaceum. Moreover, Terazosin decreased the P. aeruginosa biofilm formation and down-regulated its QS-encoding genes. Terazosin protected mice from the P. aeruginosa pathogenesis. In conclusion, α-adrenergic blockers are proposed as promising anti-virulence agents as they hinder QS receptors and inhibit bacterial espionage.

Computational and Biological Evaluation of β-Adrenoreceptor Blockers as Promising Bacterial Anti-Virulence Agents

Bacterial resistance to antibiotics is an increasing public health threat as it has the potential to affect people at any stage of life, as well as veterinary. Various approaches have been proposed to counteract the bacterial resistance development. Tackling bacterial virulence is one of the most promising approaches that confer several merits. The bacterial virulence is mainly regulated by a communication system known as quorum sensing (QS) system. Meanwhile, bacteria can sense the adrenergic hormones and eavesdrops on the host cells to establish their infection, adrenergic hormones were shown to enhance the bacterial virulence. In this study, β-adrenoreceptor blockers were proposed not only to stop bacterial espionage on our cells but also as inhibitors to the bacterial QS systems. In this context, a detailed in silico study has been conducted to evaluate the affinities of twenty-two β-blockers to compete on different structural QS receptors. Among the best docked and thermodynamically stable β-blockers; atenolol, esmolol, and metoprolol were subjected to further in vitro and in vivo investigation to evaluate their anti-QS activities against Chromobacterium violaceum, Pseudomonas aeruginosa and Salmonella typhimurium. The three tested β-blockers decreased the production of QS-controlled C. violaceum, and the formation of biofilm by P. aeruginosa and S. typhimurium. Additionally, the tested β-blockers down-regulated the P. aeruginosa QS-encoding genes and S. typhimurium sensor kinase encoding genes. Furthermore, metoprolol protected mice against P. aeruginosa and S. typhimurium. Conclusively, these investigated β-blockers are promising anti-virulence agents antagonizing adrenergic hormones induced virulence, preventing bacterial espionage, and blocking bacterial QS systems.

In Silico and In Vitro Studies for Benzimidazole Anthelmintics Repurposing as VEGFR-2 Antagonists: Novel Mebendazole-Loaded Mixed Micelles with Enhanced Dissolution and Anticancer Activity

Cancer is a leading cause of death worldwide and its incidence is unfortunately anticipated to rise in the next years. On the other hand, vascular endothelial growth factor receptor 2 (VEGFR-2) is highly expressed in tumor-associated endothelial cells, where it affects tumor-promoting angiogenesis. Therefore, VEGFR-2 is considered one of the most promising therapeutic targets for cancer treatment. Furthermore, some FDA-approved benzimidazole anthelmintics have already shown potential anticancer activities. Therefore, repurposing them against VEGFR-2 can provide a rapid and effective alternative that can be implicated safely for cancer treatment. Hence, 13 benzimidazole anthelmintic drugs were subjected to molecular docking against the VEGFR-2 receptor. Among the tested compounds, fenbendazole (FBZ, 1), mebendazole (MBZ, 2), and albendazole (ABZ, 3) were proposed as potential VEGFR-2 antagonists. Furthermore, molecular dynamics simulations were carried out at 200 ns, giving more information on their thermodynamic and dynamic properties. Besides, the anticancer activity of the aforementioned drugs was tested in vitro against three different cancer cell lines, including liver cancer (HUH7), lung cancer (A549), and breast cancer (MCF7) cell lines. The results depicted potential cytotoxic activity especially against both HUH7 and A549 cell lines. Furthermore, to improve the aqueous solubility of MBZ, it was formulated in the form of mixed micelles (MMs) which showed an enhanced drug release with better promising cytotoxicity results compared to the crude MBZ. Finally, an in vitro quantification for VEGFR-2 concentration in treated HUH7 cells has been conducted based on the enzyme-linked immunosorbent assay. The results disclosed that FBZ, MBZ, and ABZ significantly (p < 0.001) reduced the concentration of VEGFR-2, while the lowest inhibition was achieved in MBZ-loaded MMs, which was even much better than the reference drug sorafenib. Collectively, the investigated benzimidazole anthelmintics could be encountered as lead compounds for further structural modifications and thus better anticancer activity, and that was accomplished through studying their structure-activity relationships.

Computer-aided identification of potential inhibitors against Necator americanus glutathione S-transferase 3

Hookworm infection is caused by the blood-feeding hookworm gastrointestinal nematodes. Its harmful effects include anemia and retarded growth and are common in the tropics. A current control method involves the mass drug administration of synthetic drugs, mainly albendazole and mebendazole. There are however concerns of low efficacy and drug resistance due to their repeated and excessive use. Although, Necator americanus glutathione S-transferase 3 (Na-GST-3) is a notable target, using natural product libraries for computational elucidation of promising leads is underexploited. This study sought to use pharmacoinformatics techniques to identify compounds of natural origins with the potential to be further optimized as promising inhibitors. A compendium of 3182 African natural products together with five known helminth GST inhibitors including Cibacron blue was screened against the active sites of the Na-GST-3 structure (PDB ID: 3W8S). The hit compounds were profiled to ascertain the mechanisms of binding, anthelmintic bioactivity, physicochemical and pharmacokinetic properties. The AutoDock Vina docking protocol was validated by obtaining 0.731 as the area under the curve calculated via the receiver operating characteristics curve. Four compounds comprising ZINC85999636, ZINC35418176, ZINC14825190, and Dammarane Triterpene13 were identified as potential lead compounds with binding energies less than -9.0 kcal/mol. Furthermore, the selected compounds formed key intermolecular interactions with critical residues Tyr95, Gly13 and Ala14. Notably, ZINC85999636, ZINC14825190, and dammarane triterpene13 were predicted as anthelmintics, whilst all the four molecules shared structural similarities with known inhibitors. Molecular modelling showed that the compounds had reasonably good binding free energies. More so, they had high binding affinities when screened against other variants of the Na-GST, namely Na-GST-1 and Na-GST-2. Ligand quality assessment using ligand efficiency dependent lipophilicity, ligand efficiency, ligand efficiency scale and fit quality scale showed the molecules are worthy candidates for further optimization. The inhibitory potentials of the molecules warrant in vitro studies to evaluate their effect on the heme regulation mechanisms.

Calendulaglycoside A showing potential activity against SARS-CoV-2 main protease: Molecular docking, molecular dynamics, and SAR studies

BACKGROUND AND AIM

The discovery of drugs capable of inhibiting SARS-CoV-2 is a priority for human beings due to the severity of the global health pandemic caused by COVID-19. To this end, natural products can provide therapeutic alternatives that could be employed as an effective safe treatment for COVID-19.

EXPERIMENTAL PROCEDURE

Twelve compounds were isolated from the aerial parts of C. officinalis L. and investigated for their inhibitory activities against SARS-CoV-2 Mpro compared to its co-crystallized N3 inhibitor using molecular docking studies. Furthermore, a 100 ns MD simulation was performed for the most active two promising compounds, Calendulaglycoside A (SAP5) and Osteosaponin-I (SAP8).

RESULTS AND CONCLUSION

At first, molecular docking studies showed interesting binding scores as compared to the N3 inhibitor. Calendulaglycoside A (SAP5) achieved a superior binding than the co-crystallized inhibitor indicating promising affinity and intrinsic activity towards the Mpro of SARS-CoV-2 as well. Moreover, findings illustrated preferential stability for SAP5 within the Mpro pocket over that of N3 beyond the 40 ns MD simulation course. Structural preferentiality for triterpene-Mpro binding highlights the significant role of 17β-glucosyl and carboxylic 3α-galactosyl I moieties through high electrostatic interactions across the MD simulation trajectories. Furthermore, this study clarified a promising SAR responsible for the antiviral activity against the SARS-CoV-2 Mpro and the design of new drug candidates targeting it as well. The above findings could be promising for fast examining the previously isolated triterpenes both pre-clinically and clinically for the treatment of COVID-19.

β-Blockers bearing hydroxyethylamine and hydroxyethylene as potential SARS-CoV-2 Mpro inhibitors: rational based design, in silico, in vitro, and SAR studies for lead optimization

The global COVID-19 pandemic became more threatening especially after the introduction of the second and third waves with the current large expectations for a fourth one as well. This urged scientists to rapidly develop a new effective therapy to combat SARS-CoV-2. Based on the structures of β-adrenergic blockers having the same hydroxyethylamine and hydroxyethylene moieties present in the HIV-1 protease inhibitors which were found previously to inhibit the replication of SARS-CoV, we suggested that they may decrease the SARS-CoV-2 entry into the host cell through their ability to decrease the activity of RAAS and ACE2 as well. Herein, molecular docking of twenty FDA-approved β-blockers was performed targeting SARS-CoV-2 Mpro. Results showed promising inhibitory activities especially for Carvedilol (CAR) and Nebivolol (NEB) members. Moreover, these two drugs together with Bisoprolol (BIS) as an example from the lower active ones were subjected to molecular dynamics simulations at 100 ns. Great stability across the whole 100 ns timeframe was observed for the top docked ligands, CAR and NEB, over BIS. Conformational analysis of the examined drugs and hydrogen bond investigation with the pocket's crucial residues confirm the great affinity and confinement of CAR and NEB within the Mpro binding site. Moreover, the binding-free energy analysis and residue-wise contribution analysis highlight the nature of ligand-protein interaction and provide guidance for lead development and optimization. Furthermore, the examined three drugs were tested for their in vitro inhibitory activities towards SARS-CoV-2. It is worth mentioning that NEB achieved the most potential anti-SARS-CoV-2 activity with an IC50 value of 0.030 μg ml-1. Besides, CAR was found to have a promising inhibitory activity with an IC50 of 0.350 μg ml-1. Also, the IC50 value of BIS was found to be as low as 15.917 μg ml-1. Finally, the SARS-CoV-2 Mpro assay was performed to evaluate and confirm the inhibitory effects of the tested compounds (BIS, CAR, and NEB) towards the SARS-CoV-2 Mpro enzyme. The obtained results showed very promising SARS-CoV-2 Mpro inhibitory activities of BIS, CAR, and NEB (IC50 = 118.50, 204.60, and 60.20 μg ml-1, respectively) compared to lopinavir (IC50 = 73.68 μg ml-1) as a reference standard.

Investigating the structure-activity relationship of marine natural polyketides as promising SARS-CoV-2 main protease inhibitors

Since its first report in December 2019, the novel coronavirus virus, SARS-CoV-2, has caused an unprecedented global health crisis and economic loss imposing a tremendous burden on the worldwide finance, healthcare system, and even daily life. Even with the introduction of different preventive vaccines, there is still a dire need for effective antiviral therapeutics. Nature has been considered as the historical trove of drug discovery and development, particularly in cases of worldwide crises. Herein, a comprehensive in silico investigation of a highly focused chemical library of 34 pederin-structurally related marine compounds, belonging to four polyketides families, was initiated against the SARS-CoV-2 main protease, Mpro, being the key replicating element of the virus and main target in many drugs development programs. Two of the most potent SARS-CoV-2 Mpro co-crystallized inhibitors, O6K and N3, were added to the tested database as reference standards. Through molecular docking simulation, promising compounds including Pederin (1), Dihydro-onnamide A (11), Onnamide C (14), Pseudo-onnamide A (17), and Theopederin G (29) have been identified from different families based on their superior ligand-protein energies and relevant binding profiles with the key Mpro pocket residues. Thermodynamic behaviors of the identified compounds were investigated through 200 ns all-atom molecular dynamics simulation illustrating their significant stability and pocket accommodation. Furthermore, structural activity preferentiality was identified for the pederin-based marine compounds highlighting the importance of the terminal guanidine and cyclic hemiacetal linker, and the length of the sidechain. Our findings highlight the challenges of targeting SARS-CoV-2 Mpro as well as recommending further in vitro and in vivo studies regarding the examined marine products either alone or in combination paving the way for promising lead molecules.

Development of Effective Therapeutic Molecule from Natural Sources against Coronavirus Protease

The SARS-CoV-2 main protease (M^pro) is one of the molecular targets for drug design. Effective vaccines have been identified as a long-term solution but the rate at which they are being administered is slow in several countries, and mutations of SARS-CoV-2 could render them less effective. Moreover, remdesivir seems to work only with some types of COVID-19 patients. Hence, the continuous investigation of new treatments for this disease is pivotal. This study investigated the inhibitory role of natural products against SARS-CoV-2 M^pro as repurposable agents in the treatment of coronavirus disease 2019 (COVID-19). Through in silico approach, selected flavonoids were docked into the active site of M^pro. The free energies of the ligands complexed with M^pro were computationally estimated using the molecular mechanics-generalized Born surface area (MM/GBSA) method. In addition, the inhibition process of SARS-CoV-2 Mpro with these ligands was simulated at 100 ns in order to uncover the dynamic behavior and complex stability. The docking results showed that the selected flavonoids exhibited good poses in the binding domain of M^pro. The amino acid residues involved in the binding of the selected ligands correlated well with the residues involved with the mechanism-based inhibitor (N3) and the docking score of Quercetin-3-O-Neohesperidoside (−16.8 Kcal/mol) ranked efficiently with this inhibitor (−16.5 Kcal/mol). In addition, single-structure MM/GBSA rescoring method showed that Quercetin-3-O-Neohesperidoside (−87.60 Kcal/mol) is more energetically favored than N3 (−80.88 Kcal/mol) and other ligands (Myricetin 3-Rutinoside (−87.50 Kcal/mol), Quercetin 3-Rhamnoside (−80.17 Kcal/mol), Rutin (−58.98 Kcal/mol), and Myricitrin (−49.22 Kcal/mol). The molecular dynamics simulation (MDs) pinpointed the stability of these complexes over the course of 100 ns with reduced RMSD and RMSF. Based on the docking results and energy calculation, together with the RMSD of 1.98 ± 0.19 Å and RMSF of 1.00 ± 0.51 Å, Quercetin-3-O-Neohesperidoside is a better inhibitor of M^pro compared to N3 and other selected ligands and can be repurposed as a drug candidate for the treatment of COVID-19. In addition, this study demonstrated that in silico docking, free energy calculations, and MDs, respectively, are applicable to estimating the interaction, energetics, and dynamic behavior of molecular targets by natural products and can be used to direct the development of novel target function modulators.

Structure and Thermal Stability of wtRop and RM6 Proteins through All-Atom Molecular Dynamics Simulations and Experiments

In the current work we study, via molecular simulations and experiments, the folding and stability of proteins from the tertiary motif of 4-α-helical bundles, a recurrent motif consisting of four amphipathic α-helices packed in a parallel or antiparallel fashion. The focus is on the role of the loop region in the structure and the properties of the wild-type Rop (wtRop) and RM6 proteins, exploring the key factors which can affect them, through all-atom molecular dynamics (MD) simulations and supporting by experimental findings. A detailed investigation of structural and conformational properties of wtRop and its RM6 loopless mutation is presented, which display different physical characteristics even in their native states. Then, the thermal stability of both proteins is explored showing RM6 as more thermostable than wtRop through all studied measures. Deviations from native structures are detected mostly in tails and loop regions and most flexible residues are indicated. Decrease of hydrogen bonds with the increase of temperature is observed, as well as reduction of hydrophobic contacts in both proteins. Experimental data from circular dichroism spectroscopy (CD), are also presented, highlighting the effect of temperature on the structural integrity of wtRop and RM6. The central goal of this study is to explore on the atomic level how a protein mutation can cause major changes in its physical properties, like its structural stability.

Ligand-Dependent Conformational Transitions in Molecular Dynamics Trajectories of GPCRs Revealed by a New Machine Learning Rare Event Detection Protocol

Central among the tools and approaches used for ligand discovery and design are Molecular Dynamics (MD) simulations, which follow the dynamic changes in molecular structure in response to the environmental condition, interactions with other proteins, and the effects of ligand binding. The need for, and successes of, MD simulations in providing this type of essential information are well documented, but so are the challenges presented by the size of the resulting datasets encoding the desired information. The difficulty of extracting information on mechanistically important state-to-state transitions in response to ligand binding and other interactions is compounded by these being rare events in the MD trajectories of complex molecular machines, such as G-protein-coupled receptors (GPCRs). To address this problem, we have developed a protocol for the efficient detection of such events. We show that the novel Rare Event Detection (RED) protocol reveals functionally relevant and pharmacologically discriminating responses to the binding of different ligands to the 5-HT_2AR orthosteric site in terms of clearly defined, structurally coherent, and temporally ordered conformational transitions. This information from the RED protocol offers new insights into specific ligand-determined functional mechanisms encoded in the MD trajectories, which opens a new and rigorously reproducible path to understanding drug activity with application in drug discovery.

Molecular Docking and Dynamics Simulation Study of Hyrtios erectus Isolated Scalarane Sesterterpenes as Potential SARS-CoV-2 Dual Target Inhibitors

Presently, the world is under the toll of pandemic coronavirus disease-2019 (COVID-19) outbreak caused by SARS-CoV-2. Lack of effective and safe therapeutics has stressed the scientific community for developing novel therapeutics capable of alleviating and stopping this pandemic. Within the presented study, molecular docking, ADME properties and all-atom molecular dynamic (MD) simulation, along with two standard antiviral agents (lopinavir and benzopurpurin-4B), were applied to investigate 15 scalaranes sesterterpenes natural compounds, purified from the Red Sea marine sponge Hyrtios erectus, as potential COVID-19 dual-target inhibitors. Following multi-step docking within COVID-19 main protease and Nsp15 endoribonuclease cavities, nine promising drug-like compounds exhibited higher docking scores as well as better interactions with the target's crucial residues than those of reference ligands. Compounds 2, 6, 11, and 15, were predicted to simultaneously subdue the activity of the two COVID-19 targets. Dynamics behavior of the best-docked molecules, compounds 15 and 6, within COVID-19 target pockets showed substantial stability of ligand-protein complexes as presented via several MD simulation parameters. Furthermore, calculated free-binding energies from MD simulation illustrated significant ligand's binding affinity towards respective target pockets. All provided findings supported the utility of scalarane-based sesterterpenes, particularly compounds 15 and 6, as promising lead candidates guiding the development of effective therapeutics against SARS-CoV-2.

ABBV-744 as a potential inhibitor of SARS-CoV-2 main protease enzyme against COVID-19

A new pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide and become pandemic with thousands new deaths and infected cases globally. To address coronavirus disease (COVID-19), currently no effective drug or vaccine is available. This necessity motivated us to explore potential lead compounds by considering drug repurposing approach targeting main protease (Mpro) enzyme of SARS-CoV-2. This enzyme considered to be an attractive drug target as it contributes significantly in mediating viral replication and transcription. Herein, comprehensive computational investigations were performed to identify potential inhibitors of SARS-CoV-2 Mpro enzyme. The structure-based pharmacophore modeling was developed based on the co-crystallized structure of the enzyme with its biological active inhibitor. The generated hypotheses were applied for virtual screening based PhaseScore. Docking based virtual screening workflow was used to generate hit compounds using HTVS, SP and XP based Glide GScore. The pharmacological and physicochemical properties of the selected lead compounds were characterized using ADMET. Molecular dynamics simulations were performed to explore the binding affinities of the considered lead compounds. Binding energies revealed that compound ABBV-744 binds to the Mpro with strong affinity (ΔGbind -45.43 kcal/mol), and the complex is more stable in comparison with other protein-ligand complexes. Our study classified three best compounds which could be considered as promising inhibitors against main protease SARS-CoV-2 virus.

In silico approach: biological prediction of nordentatin derivatives as anticancer agent inhibitors in the cAMP pathway

A combination of computational techniques has been carried out to predict the binding of nordentatin derivatives based on pyranocoumarin semi-synthesis with the target protein from the expression of the PDE4B gene. The inhibition of the cAMP pathway is the main target of anti-cancer drugs, which is responsible for uncontrolled cell division in cancer. Modeling was done using a combination of semi-empirical methods and the density functional theory (PM3-DFT/6-31G*/B3LYP) to obtain the optimal structure of a small ligand that could be modeled. Studies on the interaction of the ligands and amino acid residues on protein targets were carried out using a combination of molecular docking and molecular dynamic simulation. Molecular docking based on functional grid scores showed a very good native ligand pose with an RMSD of 0.93 Å in determining the initial coordinates of the ligand-receptor interactions. Furthermore, the amino acid residues responsible for interaction through H-bonds were Tyr103, His104, His177, Met217, and Gln313. The binding free energy (kcal mol-1) results of the candidates were PS-1 (-36.84 ± 0.31), PS-2 (-35.34 ± 0.28), PS-3 (-26.65 ± 0.30), PS-5 (-42.66 ± 0.26), PS-7 (-35.33 ± 0.23), and PS-9 (-32.57 ± 0.20), which are smaller than that of the native ligand Z72 (-24.20 ± 0.19), and thus these have good potential as drugs that can inhibit the cAMP pathway. These results provide theoretical information for the efficient inhibition of the cAMP pathway in the future.

Investigation of dielectric constants of water in a nano-confined pore

We report the dielectric properties of pure water confined in a silica hydrophilic nanopore determined from the computation of the density of liquid in the confined phase by the grand canonical Monte Carlo (GCMC) simulations. The silica cylindrical nanopore is divided into n concentric radial shells to get a better insight into the dielectric properties of nanoconfined water. We find that the average values of the dielectric constants are very close and almost independent of the number of concentric radial shells. The decrease in the dielectric constant of confined pure water is due to the strong orientation of water dipoles in the vicinity of the surface while water dipoles do not exhibit any preferential orientation in bulk phase.

Investigation of dielectric properties in a silica hydrophilic pore with two water reservoirs.