Favorable Binding Free Energy Research Articles

Microscopic Understanding of the Conformational Stability of the Aggregated Nonamyloid β Components of α-Synuclein.

Self-association of α-synuclein peptides into oligomeric species and ordered amyloid fibrils is associated with Parkinson's disease, a progressive neurodegenerative disorder. In particular, the peptide domain formed between the residues Glu-61 (or E61) and Val-95 (or V95) of α-synuclein, typically termed the "nonamyloid β component" (NAC), is known to play critical roles in forming aggregated structures. In this work, we have employed molecular dynamics simulations to explore the conformational properties and relative stabilities of aggregated protofilaments of different orders, namely, tetramer (P(4)), hexamer (P(6)), octamer (P(8)), decamer (P(10)), dodecamer (P(12)), and tetradecamer (P(14)), formed by the NAC domains of α-synuclein. Besides, center-of-mass pulling and umbrella sampling simulation methods have also been employed to characterize the mechanistic pathway of peptide association/dissociation and the corresponding free energy profiles. Structural analysis showed that the disordered C-terminal loop and the central core regions of the peptide units lead to more flexible and distorted structures of the lower order protofilaments (P(4) and P(6)) as compared to the higher order ones. Interestingly, our calculation shows the presence of multiple distinctly populated conformational states for the lower order protofilament P(4), which may drive the oligomerization process along multiple pathways to form different polymorphic α-synuclein fibrillar structures. It is further observed that the nonpolar interaction between the peptides and the corresponding nonpolar solvation free energy play a dominant role in stabilizing the aggregated protofilaments. Importantly, our result showed that reduced cooperativity during the binding of a peptide unit beyond a critical size of the protofilament (P(12)) leads to less favorable binding free energy of a peptide.

Potential of Onchocerca ochengi inosine-5′-monophosphate dehydrogenase (IMPDH) and guanosine-5′-monophosphate oxidoreductase (GMPR) as druggable and vaccine candidates: immunoinformatics screening

Onchocerciasis is a vector-borne disease caused by the filarial nematode Onchocerca volvulus, which is responsible for most of the visual impairments recorded in Africa, Asia and the Americas. It is known that O. volvulus has similar molecular and biological characteristics as Onchocerca ochengi in cattle. This study was designed to screen for immunogenic epitopes and binding pockets of O. ochengi IMPDH and GMPR ligands using immunoinformatic approaches. In this study, a total of 23 B cell epitopes for IMPDH and 7 B cell epitopes for GMPR were predicted using ABCpred tool, Bepipred 2.0 and Kolaskar and Tongaonkar methods. The CD4+ Th computational results showed 16 antigenic epitopes from IMPDH with strong binding affinity for DRB1_0301, DRB3_0101, DRB1_0103 and DRB1_1501 MHC II alleles while 8 antigenic epitopes from GMPR were predicted to bind DRB1_0101 and DRB1_0401 MHC II alleles, respectively. For the CD8+ CTLs analysis, 8 antigenic epitopes from IMPDH showed strong binding affinity to human leukocyte antigen HLA-A*26:01, HLA-A*03:01, HLA-A*24:02 and HLA-A*01:01 MHC I alleles while 2 antigenic epitopes from GMPR showed strong binding affinity to HLA-A*01:01 allele, respectively. The immunogenic B cell and T cell epitopes were further evaluated for antigenicity, non-alllergernicity, toxicity, IFN-gamma, IL4 and IL10. The docking score revealed favorable binding free energy with IMP and MYD scoring the highest binding affinity at −6.6 kcal/mol with IMPDH and −8.3 kcal/mol with GMPR. This study provides valuable insight on IMPDH and GMPR as potential drug targets and for the development of multiple epitope vaccine candidates. Communicated by Ramaswamy H. Sarma

A comparative study on inclusion complex formation between formononetin and β-cyclodextrin derivatives through multiscale classical and umbrella sampling simulations

Formononetin, a naturally occurring isoflavone exhibits a wide range of therapeutic applications including antioxidant, anti-tumor, antiviral, anti-diabetic and neuroprotective activities. However, the low hydro-solubility of formononetin has limited its prospective use in cosmetic, neutraceutical and pharmaceutical industries. Cyclodextrins (CDs), especially β-CD and its derivatives have emerged as promising agents to improve the water solubility of poorly hydrosoluble compounds by the formation of inclusion complexes. We employed multiscale (1000 ns) explicit solvent and umbrella sampling molecular dynamics (MD) simulations to study the interactions and thermodynamic parameters of inclusion complex formation between formononetin and five most commonly used β-CD derivatives. Classical MD simulations revealed two possible binding conformations of formononetin inside the central cavity of hydroxypropyl-β-CD (HP-β-CD), randomly methylated-β-CD (ME-β-CD), and sulfobutylether-β-CD (SBE-β-CD). The binding conformation with the benzopyrone ring of formononetin inside the central cavity of β-CD derivatives was more frequent than the phenyl group occupying the hydrophobic cavity. These interactions were supported by a variety of non-bonded contacts including hydrogen bonds, pi-lone pair, pi-sigma, and pi-alkyl interactions. Formononetin showed favorable end-state MD-driven thermodynamic binding free energies with all the selected β-CD derivatives, except succinyl-β-CD (S-β-CD). Furthermore, umbrella sampling simulations were used to investigate the interactions and thermodynamic parameters of the host-guest inclusion complexes. The SBE-β-CD/formononetin inclusion complex showed the lowest binding energy signifying the highest affinity among all the selected host-guest inclusion complexes. Our study could be used as a standard for analyzing and comparing the ability of different β-CD derivatives to enhance the hydro-solubility of poorly soluble molecules.

An Attempt of Seeking Favorable Binding Free Energy Prediction Schemes Considering the Entropic Effect on Fis-DNA Binding.

Protein-DNA binding mechanisms in a complex manner are essential for understanding many biological processes. Over the past decades, numerous experiments and calculations have analyzed the specificity of protein-DNA binding. However, the accuracy of binding free energy prediction for multi-base DNA systems still needs to be improved. Fis is a DNA-binding protein that regulates various transcription and recombination reactions. In the present work, we tested several methods of predict binding free energy based on this system to find a favorable prediction scheme and explore the binding mechanism of Fis protein and DNA. Two solvent models (explicit and implicit solvent models) were chosen for the dynamics process, and the predicted binding free energy was more accurate under the explicit solvent model. When different Poisson-Boltzmann/Generalized Born (PB/GB) models were tested for DNA force fields (BSC1 and OL15), it was found that the binding free energy predicted by the selected OL15 force field performed better and the correlation between predicted and experimental values was improved with the increasing interior dielectric constant (Dk). Finally, using Dk = 8, the GBOBC1 model combined with interaction entropy (IE), which was calculated for entropic contribution (GBOBC1_IE_8), was screened out for the binding free energy prediction and analysis of the Fis-DNA system, and the validity of the method was further verified by testing the Cren7-DNA system. By performing conformational analysis of the minor groove, it was found that mutation of the DNA central sequence A/T to C/G and deletion of the guanine 2-amino group would change the minor groove width and thus affect the formation of the major groove, altering the interaction and atomic contact between the protein and the major groove, thus changing the binding affinity of Fis and DNA. Hopefully, the series of tests in this work can shed some light on the related studies of protein and DNA systems.

In silico investigation of falcipain-2 inhibition by hybrid benzimidazole-thiosemicarbazone antiplasmodial agents: A molecular docking, molecular dynamics simulation, and kinetics study.

The emergence of artemisinin-resistant variants of Plasmodium falciparum necessitates the urgent search for novel antimalarial drugs. In this regard, an in silico study to screen antimalarial drug candidates from a series of benzimidazole-thiosemicarbazone hybrid molecules with interesting antiplasmodial properties and explore their falcipain-2 (FP2) inhibitory potentials has been undertaken herein. FP2 is a key cysteine protease that degrades hemoglobin in Plasmodium falciparum and is an important biomolecular target in the development of antimalarial drugs. Pharmacokinetic properties, ADMET profiles, MM/GBSA-based binding free energies, reaction mechanisms, and associated barrier heights have been investigated. DFT, molecular dynamics simulation, molecular docking, and ONIOM methods were used. From the results obtained, four 4N-substituted derivatives of the hybrid molecule (E)-2-(1-(5-chloro-1H-benzo[d]imidazol-2-yl)ethylidene)hydrazine-1-carbothioamide (1A) denoted 1B, 1C, 1D, and 1E are drug-like and promising inhibitors of FP2, exhibiting remarkably small inhibitory constants (5.94 × 10-14 - 2.59 × 10-04 [Formula: see text]M) and favorable binding free energies (-30.32 to -17.17kcal/mol). Moreover, the ONIOM results have revealed that 1B and possibly 1C and 1D may act as covalent inhibitors of FP2. The rate-determining step of the thermodynamically favorable covalent binding mechanism occurs across a surmountable barrier height of 24.18kcal/mol in water and 28.42kcal/mol in diethyl ether. Our findings are useful for further experimental investigations on the antimalarial activities of the hybrid molecules studied.

Screening and identification of potential MERS-CoV papain-like protease (PLpro) inhibitors; Steady-state kinetic and Molecular dynamic studies

MERS-CoV belongs to the coronavirus group. Recent years have seen a rash of coronavirus epidemics. In June 2012, MERS-CoV was discovered in the Kingdom of Saudi Arabia, with 2,591 MERSA cases confirmed by lab tests by the end of August 2022 and 894 deaths at a case-fatality ratio (CFR) of 34.5% documented worldwide. Saudi Arabia reported the majority of these cases, with 2,184 cases and 813 deaths (CFR: 37.2%), necessitating a thorough understanding of the molecular machinery of MERS-CoV. To develop antiviral medicines, illustrative investigation of the protein in coronavirus subunits are required to increase our understanding of the subject. In this study, recombinant expression and purification of MERS-CoV (PLpro), a primary goal for the development of 22 new inhibitors, were completed using a high throughput screening methodology that employed fragment-based libraries in conjunction with structure-based virtual screening. Compounds 2, 7, and 20, showed significant biological activity. Moreover, a docking analysis revealed that the three compounds had favorable binding mood and binding free energy. Molecular dynamic simulation demonstrated the stability of compound 2 (2-((Benzimidazol-2-yl) thio)-1-arylethan-1-ones) the strongest inhibitory activity against the PLpro enzyme. In addition, disubstitutions at the meta and para locations are the only substitutions that may boost the inhibitory action against PLpro. Compound 2 was chosen as a MERS-CoV PLpro inhibitor after passing absorption, distribution, metabolism, and excretion studies; however, further investigations are required.

Structural basis of Nrf2 activation by flavonolignans from silymarin

Several properties of silymarin (SM) extract have been attributed to their three major flavonolignans (silybin, silychristin, and silydianin) and their 2,3-dehydro derivatives (2,3-dehydrosilybin, 2,3-dehydrosilychristin, and 2,3-dehydrosilydianin). Experimental findings have suggested that the antioxidative and protective activities of these compounds could be due to their ability to activate nuclear factor erythroid 2-related factor 2 (Nrf2). The mechanism by which SM compounds exert their effect has been suggested to be by disrupting the complex between Nrf2 and Kelch-like ECH-associated protein 1 (Keap1). However, information about the structural and energetic basis of the inhibitory mechanism of SM compounds on the Nrf2-Keap1 pathway is lacking. We evaluated the binding properties of SM compounds because experimental findings have pointed to them as potential activators of Nrf2. Our study combined docking and molecular dynamics (MD) simulations with the Poisson–Boltzmann and generalized Born and surface area (MMPBSA and MMGBSA) methods and quantum mechanics-molecular mechanics (QMMM) calculations to investigate Keap1–ligand interactions. Our results predicted that silybinA and 2,3-dehydrosilybin bind to Keap1, forming interactions with the same pockets as those observed for the cocrystallized Keap1-Cpd16 complex but with more favorable binding free energies. These findings indicate that both natural compounds are potential activators of Nrf2.

Identification of Novel AXL Kinase Inhibitors Using Ligand-Based Pharmacophore Screening and Molecular Dynamics Simulations

AXL kinase is a promising target in novel drug discovery for cancer. A ligand-based pharmacophore model was generated with the Pharmit web server. Its inbuilt PubChem molecule database was screened and led to 408 candidate molecules. Docking of the AXL kinase active sites with the identified list of candidate molecules was carried out with Autodock Vina docking software. This resulted in four compounds selected for further investigation. Molecular dynamics simulation of two ligands (PubChem-122421875 and PubChem-78160848) showed considerable binding with AXL kinase. From the MM-PBSA binding free energies investigation, the PubChem-122421875 (G = −179.3 kJ/mol) and PubChem-78160848 (G = −208.3 kJ/mol) ligands had favorable protein-ligand complex stability and binding free energy. Hence, PubChem-122421875 and PubChem-78160848 molecules identified in this work could be a potent starting point for developing novel AXL kinase inhibitor molecules.

A Molecular Dynamics Study of Antimicrobial Peptide Interactions with the Lipopolysaccharides of theOuterBacterial Membrane.

With rising bacterial resistance, antimicrobial peptides (AMPs) have been widely investigated as potential antibacterial molecules to replace conventional antibiotics. Our understanding of the molecular mechanisms for membrane disruption are largely based on AMP interactions with the inner phospholipid bilayers of both Gram-negative and Gram-positive bacteria. Mechanisms for AMP translocation across the outer membrane of Gram-negative bacteria composed of lipopolysaccharides and the asymmetric lipid bilayer are complicated by the secondary structure adopted by the peptide in the different membrane environments. We have employed atomistic molecular dynamics and umbrella-sampling simulations with an aggregate duration of [Formula: see text] 6 microseconds to obtain the free energy landscape of CM15 peptide translocating through the lipopolysaccharide region of Gram-negative bacteria, E. coli. The peptide has a favorable binding-free energy (-130 kJ mol[Formula: see text]) in the O-antigen region with a large barrier (150 kJ mol[Formula: see text]) at the interface between the anionic core saccharides and upper bilayer leaflet made up of lipid-A molecules. Restraint-free molecular dynamics simulations show that the random coil structure is favored over the helix in both the extracellular aqueous region and the cation-rich core-saccharide regions of the outer membrane. The peptide and membrane properties are analyzed at each of the 100 ns duration of the umbrella-sampling windows toillustrate changes in peptide length, orientation, and hydration. Our study provides insights into the free energy landscape for the insertion of the AMP CM15 in the outer membrane of Gram-negative bacteria, and we discuss the implications of our findings with the broader question of how AMPs overcome this barrier during antimicrobial activity.

Antiviral activity of ouabain against a Brazilian Zika virus strain

Zika virus (ZIKV) is an emerging arbovirus associated with neurological disorders. Currently, no specific vaccines or antivirals are available to treat the ZIKV infection. Ouabain, a cardiotonic steroid known as Na+/K+-ATPase inhibitor, has been previously described as an immunomodulatory substance by our group. Here, we evaluated for the first time the antiviral activity of this promising substance against a Brazilian ZIKV strain. Vero cells were treated with different concentrations of ouabain before and after the infection with ZIKV. The antiviral effect was evaluated by the TCID50 method and RT-qPCR. Ouabain presented a dose-dependent inhibitory effect against ZIKV, mainly when added post infection. The reduction of infectious virus was accompanied by a decrease in ZIKV RNA levels, suggesting that the mechanism of ZIKV inhibition by ouabain occurred at the replication step. In addition, our in silico data demonstrated a conformational stability and favorable binding free energy of ouabain in the biding sites of the NS5-RdRp and NS3-helicase proteins, which could be related to its mechanism of action. Taken together, these data demonstrate the antiviral activity of ouabain against a Brazilian ZIKV strain and evidence the potential of cardiotonic steroids as promising antiviral agents.

Exploring the Potential Role of Theaflavin-3,3′-Digallate in Inhibiting Various Stages of SARS-CoV-2 Life Cycle: An In-Silico Approach

IntroductionTheaflavins belong to the class of polyphenols that are predominantly found in black tea. The major derivatives of theaflavins found in black tea are theaflavin (TF1), theaflavin-3-gallate (TF2A), theaflavin-3′-gallate (TF2B), and theaflavin-3,3′-digallate (TF3). Theaflavin-3,3′-digallate (TF3) is a natural compound present in black tea and known to possess antiviral activity. This study had attempted to explore the potential role of TF3 in inhibiting various stages of the SARS-CoV-2 life cycle.MethodsMolecular docking studies of TF3 along with positive controls was performed on eight different targets of SARS-CoV-2 followed by binding free energy (MM-GBSA) calculations. The docked complexes with favourable docking and binding free energy results were subjected to molecular dynamics (MD) simulation studies to assess the stability of the dock complex. Finally, TF3 and all the positive controls were taken for ADMET analysis.ResultsThe docking and binding free energy results of TF3 was compared against the positive controls. TF3 showed the highest binding energy against all the targets and formed more stable interactions for a longer duration on MD simulations with CLpro, RdRp, helicase and spike protein. Also, the promising in-silico ADMET profile further warrants the exploration of this compound through in-vitro and in-vivo methods.ConclusionThrough this study, we tried to evaluate the role of theaflavin-3,3’-digallate on multiple targets of SARS-CoV-2, and the positive in-silico results which were obtained on various pharmacodynamic and pharmacokinetic parameters, give a ray of hope as a potential therapeutic drug to this rapidly spreading disease. The search for a curative therapy for SARS-CoV-2 is still ongoing. The favourable preliminary results of TF3 through in-silico analysis offers a ray of hope in ending this devasting pandemic.Supplementary InformationThe online version contains supplementary material available at 10.1007/s42250-022-00376-7.

Synthesis and characterization of single‐walled carbon nanotube: Cyto‐genotoxicity in Allium cepa root tips and molecular docking studies

Herein, single-walled carbon nanotubes (SWCNTs) were synthesized by the thermal chemical vapor deposition (CVD) method, and characterized by scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDX), transmission electron microscopy (TEM), Raman spectroscopy, dynamic light scattering (DLS), and thermo-gravimetric analysis (TGA). The results indicated that obtained nanotubes were SWCNTs with high crystallinity and their average diameter was 10.15 ± 3nm. Allium cepa ana-telophase and comet assays on the root meristem were employed to evaluate the cytotoxic and genotoxic effects of SWCNTs by examining mitotic phases, mitotic index (MI), chromosomal aberrations (CAs), and DNA damage. A. cepa root tip cells were exposed to SWCNTs at concentrations of 12.5, 25, 50, and 100 μg/ml for 4h. Distilled water and methyl methanesulfonate (MMS, 10μg/ml) were used as the negative and positive control groups, respectively. It was observed that MIs decreased statistically significantly for all applied doses. Besides, CAs such as chromosome laggards, disturbed anaphase-telophase, stickiness and bridges and also DNA damage increased in the presence of SWCNTs in a concentration-dependent manner. In the molecular docking study, the SWCNT were found to be a strong DNA major groove binder showing an energetically very favorable binding free energy of -21.27 kcal/mol. Furthermore, the SWCNT interacted effectively with the nucleotides on both strands of DNA primarily via hydrophobic π and electrostatic interactions. As a result, cytotoxic and genotoxic effects of SWCNTs in A. cepa root meristematic cells which is a reliable system for assessment of nanoparticle toxicology were demonstrated in this study.

Synthesis, crystal, and Hirschfeld surface, DFT and molecular docking studies of 6-(3‑chloro-4-fluorophenyl)-4-ethoxy-2-(4-methoxyphenyl)quinazoline derivative

The 4‑chloro-6-iodo-2-(4-methoxyphenyl)quinazoline 2 undergoes nucleophilic substitution with sodium ethoxide followed by the Suzuki–Miyaura reactions to afford novel 6-(3‑chloro-4-fluorophenyl)-4-ethoxy-2-(4-methoxyphenyl)quinazoline 4. Structural elucidation was carried out using a combination of spectroscopic (NMR, IR and MS) and single crystal X-ray diffraction (XRD) techniques. Compound 4 crystallizes in monoclinic space group C2/c with crystal parameters a = 28.116(4) Å, b = 7.2267(9) Å, c = 20.975(3) Å, β = 116.727(4)°, V = 3806.6(8) Å3 and Z = 8. The crystal packing is stabilized by intermolecular hydrogen bond, C–H⋯F, C–H⋯Cl, C–H⋯π, and π⋯π interactions. The halogen Cl⋯Cl contacts were found to be of the type I having the angles ϴ1 and ϴ2 = 154ᵒ. DFT analysis of the optimized geometry is in agreement with the solid phase results. The Hirshfeld surface was used to analyze the properties of the intermolecular interactions of 4. The energy-framework analysis reveals that π···π and C–H···π interactions energies are mainly dispersive and are the most important forces in the crystal structure. The in-silico evaluation of compounds 4 and vandetanib with the vascular endothelial growth factor receptor-2 (VEGFR-2) was carried out. Molecular docking results showed that 4 has the most favourable binding free energy (−9.39 kcal/mol) compared with that of vandetanib (−8.31 kcal/mol). The MD simulation results reveal there is a presence of hydrogen bonds between compound 4 and the VEGFR-2 protein.

Coupling of acceptor-substituted diazo compounds and tertiary thioamides: synthesis of enamino carbonyl compounds and their pharmacological evaluation.

This paper describes the synthesis of enamino carbonyl compounds by the copper(i)-catalyzed coupling of acceptor-substituted diazo compounds and tertiary thioamides. We plan to use this method to synthesize indolizidine (−)-237D analogs to find α6-selective antismoking agents. Therefore, we also performed in silico α6-nAchRs binding studies of selected products. Compounds with low root-mean-square deviation values showed more favorable binding free energies. We also report preliminary pharmacokinetic data on indolizidine (−)-237D and found it to have weak activity at CYP3A4. In addition, as enamino carbonyl compounds are also known for antimicrobial properties, we screened previously reported and new enamino carbonyl compounds for antibacterial, antimicrobial, and antifungal properties. Eleven compounds showed significant antimicrobial activities.

Molecular mechanism of anionic stabilizer for telomere G-quadruplex.

Telomere DNA assumes a high-order G-quadruplex (G4) structure, stabilization of which prevents telomere lengthening by telomerase in cancer. Through applying combined molecular simulation methods, an investigation on the selective binding mechanism of anionic phthalocyanine 3,4',4'',4'''-tetrasulfonic acid (APC) and human hybrid (3 + 1) G4s was firstly performed at the atomic level. Compared to the groove binding mode of APC and the hybrid type I (hybrid-I) telomere G4, APC preferred to bind to the hybrid type II (hybrid-II) telomere G4 via end-stacking interactions, which showed much more favorable binding free energies. Analyses of the non-covalent interaction and binding free energy decomposition revealed a decisive role of van der Waals interaction in the binding of APC and telomere hybrid G4s. And the binding of APC and hybrid-II G4 that showed the highest binding affinity adopted the end-stacking binding mode to form the most extensive van der Waals interactions. These findings add new knowledge to the design of selective stabilizers targeting telomere G4 in cancer.

Pharmacophore mapping of the crucial mediators of dual inhibitor activity of PanK and PyrG in tuberculosis disease

ABSTRACT Tuberculosis (TB) is an infection caused by Mycobacterium tuberculosis, responsible for 1.66 million documented deaths in 2020. The increase in reported cases of M. tuberculosis resistance to the main drugs shows the need to develop new and efficient drugs for better TB control. A strategy to overcome the resistance problem could be the polypharmacology approach, whereby researchers develop a single molecule that acts on different targets. Polypharmacology could have features that make it more effective than classical polypharmacy, in which investigators pull various drugs with high affinity for one target together. Based on the effectiveness of multitargeting, this work aimed to use combined in silico techniques to identify potential dual inhibitors for Pantothenate Kinases (PanK) and CTP synthetase PyrG. In silico drug design allows for screening of potential leads, thus decreasing time and resource consumption. This study demonstrates an optimised and proven screening technique to discover a potential small-molecule inhibitor of Mycobacterium tuberculosis PanK and PyrG. The favourable binding free energy values and interaction profiles show ZINC53995076 as a potent dual inhibitor of the Mycobacterium tuberculosis PanK and PyrG.

Substrate Access Mechanism in a Novel Membrane-Bound Phospholipase A of Pseudomonas aeruginosa Concordant with Specificity and Regioselectivity.

PlaF is a cytoplasmic membrane-bound phospholipase A1 from Pseudomonas aeruginosa that alters the membrane glycerophospholipid (GPL) composition and fosters the virulence of this human pathogen. PlaF activity is regulated by a dimer-to-monomer transition followed by tilting of the monomer in the membrane. However, how substrates reach the active site and how the characteristics of the active site tunnels determine the activity, specificity, and regioselectivity of PlaF for natural GPL substrates have remained elusive. Here, we combined unbiased and biased all-atom molecular dynamics (MD) simulations and configurational free-energy computations to identify access pathways of GPL substrates to the catalytic center of PlaF. Our results map out a distinct tunnel through which substrates access the catalytic center. PlaF variants with bulky tryptophan residues in this tunnel revealed decreased catalysis rates due to tunnel blockage. The MD simulations suggest that GPLs preferably enter the active site with the sn-1 acyl chain first, which agrees with the experimentally demonstrated PLA1 activity of PlaF. We propose that the acyl chain-length specificity of PlaF is determined by the structural features of the access tunnel, which results in favorable free energy of binding of medium-chain GPLs. The suggested egress route conveys fatty acid (FA) products to the dimerization interface and, thus, contributes to understanding the product feedback regulation of PlaF by FA-triggered dimerization. These findings open up opportunities for developing potential PlaF inhibitors, which may act as antibiotics against P. aeruginosa.

Molecular docking and simulation studies of synthetic protease inhibitors against COVID-19: a computational study

COVID-19 is the most recent threat to global health. Many people preferred treatment in case of infection instead of vaccination. The inhibition of viral replication is a good strategy for the treatment of COVID-19 infection. 3CLpro and PLpro are two important viral proteases responsible for proteolysis, infection, and replication of the virus. Therefore, targeting of these two enzymes is an attractive way to deal with COVID-19. The aim of this study was to screen some synthetic protease inhibitors to determine an appropriate hit molecule against COVID-19 using molecular docking and molecular dynamic simulations. The strategy depends on docking existing synthetic compounds mostly HIV protease inhibitors against two COVID-19 proteases to identify promising drugs for the treatment of COVID-19. We used protein data bank to obtain the X-ray crystal structure of the most important COVID-19 proteases 3CL pro (PDB ID: 6M2N) and PL pro (PDB ID: 6WX4). In this conceptual context, an attempt has been made to suggest an in silico computational relationship between 50 synthetic protease inhibitors and COVID-19 proteases. Out of 50 screened compounds, the best docking scores were found for these five protease inhibitors BDBM7021, BDBM698, BDBM694, BDBM93239, BDBM700. A 100-ns MD simulation was carried out to assess the stability of COVID-19 proteases and inhibitors, revealing an average RMSD value of 0.7 and favorable binding free energy (MM-GBSA) for all complexes confirming their potency as powerful binders in the COVID-19 proteases’ binding pocket. Furthermore, the current results must be confirmed using in-vitro and in-vivo antiviral methods. Communicated by Ramaswamy H. Sarma

Potential Novel Thioether-Amide or Guanidine-Linker Class of SARS-CoV-2 Virus RNA-Dependent RNA Polymerase Inhibitors Identified by High-Throughput Virtual Screening Coupled to Free-Energy Calculations.

SARS-CoV-2, or severe acute respiratory syndrome coronavirus 2, represents a new pathogen from the family of Coronaviridae that caused a global pandemic of COVID-19 disease. In the absence of effective antiviral drugs, research of novel therapeutic targets such as SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) becomes essential. This viral protein is without a human counterpart and thus represents a unique prospective drug target. However, in vitro biological evaluation testing on RdRp remains difficult and is not widely available. Therefore, we prepared a database of commercial small-molecule compounds and performed an in silico high-throughput virtual screening on the active site of the SARS-CoV-2 RdRp using ensemble docking. We identified a novel thioether-amide or guanidine-linker class of potential RdRp inhibitors and calculated favorable binding free energies of representative hits by molecular dynamics simulations coupled with Linear Interaction Energy calculations. This innovative procedure maximized the respective phase-space sampling and yielded non-covalent inhibitors representing small optimizable molecules that are synthetically readily accessible, commercially available as well as suitable for further biological evaluation and mode of action studies.

Studies of Interaction Mechanism between Pyrido [3,4-d] Pyrimidine Inhibitors and Mps1.

Monopolar spindle 1 (Mps1), a dual-specific kinase, is related to the proper execution of chromosome biorientation and mitotic checkpoint signaling. The overexpression of Mps1 promotes the occurrence of cancer or the survival of aneuploid cancer cells, in other words, the reduction of Mps1 will severely reduce the viability of human cancer cells. Therefore, Mps1 is a potential target for cancer treatment. Recently, a series of novel pyrido [3,4-d] pyrimidine derivatives targeting Mps1 with high biological activity were synthesized. The crystal structure of Mps1 in complex with pyrido [3,4-d] pyrimidine derivatives was also reported, but there were no specific mechanism studies for this series of small molecule inhibitors. In this study, complexes binding modes were probed by molecular docking and further validated by molecular dynamics simulations and the molecular mechanics/generalized Born surface area (MM/GBSA) method. The results indicated that the van der Waals interactions and the nonpolar solvation energies were responsible to the basis for favorable binding free energies, all inhibitors interacted with residues I531, V539, M602, C604, N606, I607, L654, I663, and P673 of Mps1. By analyzing the hydrogen bonds, we found the residues G605 and K529 in Mps1 formed stable hydrogen bonds with compounds, it was more conducive to activities of Mps1 inhibitors. According to the above analysis, we further designed five new compounds. We found that compounds IV and V were better potential Mps1 inhibitors through docking and ADMET prediction. The obtained new insights not only were helpful in understanding the binding mode of inhibitors in Mps1, but also provided important references for further rational design of Mps1 inhibitors.