Mechanics Poisson-Boltzmann Surface Area Research Articles

Identification of potential SARS-CoV-2 entry inhibitors by targeting the interface region between the spike RBD and human ACE2

Coronavirus disease 2019 (COVID-19) is a fatal infectious disease caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). The virus infection is initiated upon recognition and binding of the spike (S) protein receptor-binding domain (RBD) to the host cell surface receptor, angiotensin-converting enzyme 2 (ACE2). Blocking the interaction between S protein and ACE2 receptor is a novel approach to prevent the viral entry into the host cell. The present study is aimed at the identification of small molecules which can disrupt the interaction between SARS-CoV-2 S protein and human ACE2 receptor by binding to the interface region. A chemical library consisting of 1,36,191 molecules were screened for drug-like compounds based on Lipinski’s rule of five, Verber’s rule and in silico toxicity parameters. The filtered drug-like molecules were next subjected to molecular docking in the interface region of RBD. The best three hits viz; ZINC64023823, ZINC33039472 and ZINC00991597 were further taken for molecular dynamics (MD) simulation studies and binding free energy evaluations using Molecular mechanics-Poisson–Boltzmann surface area (MM-PBSA) and Molecular mechanics-Generalized Born surface area (MM-GBSA). The protein-ligand complexes showed stable trajectories throughout the simulation time. ZINC33039472 exhibited binding free energy value lower as compared to the control (emodin) with a higher contribution by gas-phase energy and van der Waals energy to the total binding free energy. Thus, ZINC33039472 is identified to be a promising interfacial binding molecule which can inhibit the interaction between the viral S protein and human ACE2 receptor which would consequently help in the management of the disease.

Detailed Characterization of the Cooperative Binding of Piperine with Heat Shock Protein 70 by Molecular Biophysical Approaches.

In this work, for the first time, details of the complex formed by heat shock protein 70 (HSP70) independent nucleotide binding domain (NBD) and piperine were characterized through experimental and computational molecular biophysical methods. Fluorescence spectroscopy results revealed positive cooperativity between the two binding sites. Circular dichroism identified secondary conformational changes. Molecular dynamics along with molecular mechanics Poisson Boltzmann surface area (MM/PBSA) reinforced the positive cooperativity, showing that the affinity of piperine for NBD increased when piperine occupied both binding sites instead of one. The spontaneity of the complexation was demonstrated through the Gibbs free energy (∆G < 0 kJ/mol) for different temperatures obtained experimentally by van’t Hoff analysis and computationally by umbrella sampling with the potential of mean force profile. Furthermore, the mean forces which drove the complexation were disclosed by van’t Hoff and MM/PBSA as being the non-specific interactions. In conclusion, the work revealed characteristics of NBD and piperine interaction, which may support further drug discover studies.

Caffeine and caffeine-containing pharmaceuticals as promising inhibitors for 3-chymotrypsin-like protease of SARS-CoV-2

In December 2019, a new coronavirus named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) led to the outbreak of a pulmonary disease called COVID-19, which killed thousands of people worldwide. Therefore, the necessity to find out the potential therapeutic pharmaceuticals is imperious. This study investigates the inhibitory effect of SARS-CoV-2 3-chymotrypsin-like protease (3CLpro) using caffeine and caffeine-containing pharmaceuticals (3CPs) based on molecular dynamics simulations and free energy calculations by means of molecular mechanics-Poisson–Boltzmann surface area (MMPBSA) and molecular mechanics-generalized-Born surface area (MMGBSA). Of these 3CPs, seven drugs approved by the US-Food and Drug Administration have shown a good binding affinity to the catalytic residues of 3CLpro of His41 and Cys145: caffeine, theophylline, dyphylline, pentoxifylline, linagliptin, bromotheophylline and istradefylline. Their binding affinity score ranged from –4.9 to –8.6 kcal/mol. The molecular dynamic simulation in an aqueous solution of docked complexes demonstrated that the 3CPs conformations bound to the active sites of 3CLpro during 200 ns molecular dynamics simulations. The free energy of binding also confirms the stability of the 3CPs–3CLpro complexes. To our knowledge, this in silico study shows for the first time very inexpensive drugs available in large quantities that can be potential inhibitors against 3CLpro. In particular, the repurposing of linagliptin, and caffeine are recommended for COVID-19 treatment after in vitro, in vivo and clinical trial validation. Communicated by Ramaswamy H. Sarma

Molecular docking and molecular dynamics simulation of Bacillus thuringiensis dehalogenase against haloacids, haloacetates and chlorpyrifos

The high dependency and surplus use of agrochemical products have liberated enormous quantities of toxic halogenated pollutants into the environment and threaten the well-being of humankind. Herein, this study performed molecular docking, molecular dynamic (MD) simulations, molecular mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculations on the DehH2 from Bacillus thuringiensis, to identify the order of which the enzyme degrades different substrates, haloacids, haloacetate and chlorpyrifos. The study discovered that the DehH2 favored the degradation of haloacids and haloacetates (−3.3 − 4.6 kcal/mol) and formed three hydrogen bonds with Asp125, Arg201 and Lys202. Despite the inconclusive molecular docking result, chlorpyrifos was consistently shown to be the least favored substrate of the DehH2 in MD simulations and MM-PBSA calculations. Results of MD simulations revealed the DehH2-haloacid- (RMSD 0.15 − 0.25 nm) and DehH2-haloacetates (RMSF 0.05 − 0.25 nm) were more stable, with the DehH2-L-2CP complex being the most stable while the least was the DehH2-chlorpyrifos (RMSD 0.295 nm; RMSF 0.05 − 0.59 nm). The Molecular Mechanics Poisson-Boltzmann Surface Area calculations showed the DehH2-L-2CP complex (−24.27 kcal/mol) having the lowest binding energy followed by DehH2-MCA (−22.78 kcal/mol), DehH2-D-2CP (−21.82 kcal/mol), DehH2-3CP (−21.11 kcal/mol), DehH2-2,2-DCP (−18.34 kcal/mol), DehH2-2,3-DCP (−8.34 kcal/mol), DehH2-TCA (−7.62 kcal/mol), while chlorpyrifos was unable to spontaneously bind to DehH2 (+127.16 kcal/mol). In a nutshell, the findings of this study offer valuable insights into the rational tailoring of the DehH2 for expanding its substrate specificity and catalytic activity in the near future. Communicated by Ramaswamy H. Sarma

A Computational Study to Identify Potential Inhibitors of SARS-CoV-2 Main Protease (Mpro) from Eucalyptus Active Compounds

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was found to be a severe threat to global public health in late 2019. Nevertheless, no approved medicines have been found to inhibit the virus effectively. Anti-malarial and antiviral medicines have been reported to target the SARS-CoV-2 virus. This paper chose eight natural eucalyptus compounds to study their binding interactions with the SARS-CoV-2 main protease (Mpro) to assess their potential for becoming herbal drugs for the new SARS-CoV-2 infection virus. In-silico methods such as molecular docking, molecular dynamics (MD) simulations, and Molecular Mechanics Poisson Boltzmann Surface Area (MM/PBSA) analysis were used to examine interactions at the atomistic level. The results of molecular docking indicate that Mpro has good binding energy for all compounds studied. Three docked compounds, α-gurjunene, aromadendrene, and allo-aromadendrene, with highest binding energies of −7.34 kcal/mol (−30.75 kJ/mol), −7.23 kcal/mol (−30.25 kJ/mol), and −7.17 kcal/mol (−29.99 kJ/mol) respectively, were simulated with GROningen MAchine for Chemical Simulations (GROMACS) to measure the molecular interactions between Mpro and inhibitors in detail. Our MD simulation results show that α-gurjunene has the strongest binding energy of −20.37 kcal/mol (−85.21 kJ/mol), followed by aromadendrene with −18.99 kcal/mol (−79.45 kJ/mol), and finally allo-aromadendrene with −17.91 kcal/mol (−74.95 kJ/mol). The findings indicate that eucalyptus may be used to inhibit the Mpro enzyme as a drug candidate. This is the first computational analysis that gives an insight into the potential role of structural flexibility during interactions with eucalyptus compounds. It also sheds light on the structural design of new herbal medicinal products against Mpro.

In silico characterisation of olive phenolic compounds as potential cyclooxygenase modulators. Part 1

Non-steroidal anti-inflammatory drugs (NSAIDs) are widely used to reduce pain. These target cyclooxygenase (COX) enzymes which produce inflammatory mediators. Adverse effects associated with the use of traditional NSAIDs have led to a rise in the development of alternative therapies. Derived from Olea Europaea, olive oil is a main component of the Mediterranean diet, containing phenolic compounds that contribute to its antioxidant and anti-inflammatory properties. It has previously been found that oleocanthal, a phenolic compound derived from the olive, had similar effects to ibuprofen, a commonly used NSAID. There is an abundance of olive phenolic compounds that have yet to be investigated for their anti-inflammatory properties. In this study, it was sought to identify potential olive-derived compounds with the ability to inhibit COX enzymes, and study the mechanisms using in silico approaches. Molecular docking was employed to determine the COX inhibitory potential of an olive phenolic compound library. From docking, it was determined that 1-oleyltyrosol (1OL) and ligstroside derivative 2 (LG2) demonstrated the greatest binding affinity to both COX-1 and COX-2. Interactions with these compounds were further examined using molecular dynamics simulations. The residue contributions to binding free energy were computed using Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) methods, revealing that residues Leu93, Val116, Leu352, and Ala527 in COX-1 and COX-2 were key determinants of potential inhibition. Along with part 2 of this study, this work aims to identify and characterise novel phenolic compounds which may possess COX inhibitory properties.

Exploring the energetic basis of binding of currently used drugs against HIV-1 subtype CRF01_AE protease via molecular dynamics simulations

Non-B strains cause nearly 90% of the worldwide human immunodeficiency virus (HIV) infections. At the same time, the protease inhibitors (PIs) were designed for subtype B. Therefore, the use of PIs in the non-B subtype context requires further investigation. Herein, we have investigated the effectiveness of currently used four PIs, namely atazanavir, darunavir, lopinavir, and tipranavir against subtype CRF01_AE (PRCRF) by employing the MD/MMPBSA (molecular dynamics/molecular mechanics Poisson-Boltzmann surface area) scheme. Our investigation reveals that tipranavir is the most potent inhibitor against PRCRF while the other three PIs display a similar binding affinity. The energetic penalty arises due to the desolvation of polar groups always disfavor the association between PRCRF and PI, and this contribution is the least in the case of tipranavir/PRCRF compared to the other three PI-PRCRF complexes resulting in a better binding affinity for tipranavir. Further, it is revealed that the primary interaction controlling the binding of inhibitors with PRCRF is the van der Waals forces. The dynamic cross-correlation map and principal component analysis show that the anti-correlated motion at the flap region of PRCRF is diminished after the ligand binding. Further, our studies indicate that D25' forms a stable H-bond with darunavir, lopinavir, and tipranavir, while D25 forms a stable H-bond with atazanavir. The per-residue based decomposition of free energy reveals the actual residual origin of the binding free energy and identify the hotspot residues. Overall, the data presented in this study can guide the computer-aided rational design of more potent drugs targetting HIV-1 PRCRF. Communicated by Ramaswamy H. Sarma

Identification of Phosphatidylinositol 3-kinase δ (PI3Kδ) Inhibitor: Pharmacophore-based Virtual Screening and Molecular Dynamics Simulation

Phosphatidylinositol 3-kinase δ (PI3Kδ) is a validated drug target for the treatment of cancer. The present study aims to search for new inhibitors of PI3Kδ by employing pharmacophore modelling using LigandScout Advanced 4.3 software. The three hydrogen bond acceptors and two hydrophobic features were proposed as a pharmacophore model using LASW1976 structure. The model was then validated using the Area Under Curve (AUC) of Receiver Operating Characteristic (ROC) and GH score. It was used to screen new molecules in the ZINC database, which resulted in 599 hits. All 599 hits were then docked into PI3Kδ protein, and five best hits were submitted to 50 ns molecular dynamics simulations. Each hit complexed with PI3Kδ underwent minor conformational changes as indicated by the values of Root Mean Square Deviation (RMSD) and Root Mean Square Fluctuation (RMSF). Furthermore, prediction of the binding free energy using Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) method showed that five hits, i.e., Lig25/ZINC253496376, Lig682/ZINC98047241, Lig449/ZINC85878047, Lig554/ZINC253389510, and Lig199/ZINC12638303, had lower binding energy compared to LASW1976. This result indicated their potentials as new inhibitors of PI3Kδ.

Identification of novel and selective agonists for ABA receptor PYL3

Abscisic acid (ABA) although complicated and expensive to produce, plays an important role in signalling responsible for regulation of developmental manifestations such as seed maturation and surviving through stress conditions. Hence, development of cost effective molecules with minimal side effects that mimic the functions of ABA is the need of the hour. In this agreement, we screened a series of 27 in-house synthesized 3-methyleneisoindolin-1-one molecules over three ABA receptors (PYR1, PYL1, and PYL3). The commercial ABA agonist Pyrabactin was taken as a standard ligand in this study. The top three molecules for each receptor were selected and further evaluated to estimate the dynamical contribution and complex stability via Molecular Mechanics-Poisson Boltzmann surface area calculations. Two molecules (Mol26 and Mol25) showed higher binding free energy and stable complex conformation for PYL3 in comparison to Pyrabactin. This study revealed the structural basis of the binding mechanism of 3-methyleneisoindolin-1-one molecules with ABA receptors. Mol26 and Mol25 were identified for the development of specific PYL3 agonists with a vast potential in agriculture to accentuate the ABA like action in plants.

Structural based study to identify new potential inhibitors for dual specificity tyrosine-phosphorylation- regulated kinase

Background and Objectives The Dual-specificity tyrosine-phosphorylation regulated kinase-1A (DYRK1A) a serine/threonine kinase that has freshly gained recognition as an essential drug target due to the discovery of its involvement in pathological diseases. The development of new potent inhibitors of DYRK1A would contribute to clarify the molecular mechanisms of associated diseases. It would administer a new lead compound for molecular-targeted protein, which was the primary focus of our study.Methods The library of in-house synthesized pyrrolone-fused benzosuberene (PBS) compounds was docked with DYRK1A receptor. Further, molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) estimations were conducted to confirm our docking outcomes and compared the stability of chosen complexes with the 2C3 (standard molecule) complex.Results This study reports Ligand15, Ligand14, and Ligand11 as potent inhibitors which showed higher ligand efficiency, binding affinity, lipophilic ligand efficiency, and favorable torsion values as compared to 2C3.Conclusion The stated methodologies revealed a unique mechanism of active site binding. The binding interactions within the active site showed that the chosen molecules had notable interactions than the standard molecule, which led to the generation of potential compounds to inhibit DYRK1A.

Molecular docking and molecular dynamics simulations of a mutant Acinetobacter haemolyticus alkaline-stable lipase against tributyrin

We previously reported on a mutant lipase KV1 (Mut-LipKV1) from Acinetobacter haemolyticus which optimal pH was raised from 8.0 to 11.0 after triple substitutions of surface aspartic acid (Asp) with lysine (Lys). Herein, this study further examined the Mut-LipKV1 by molecular docking, molecular dynamics (MD) simulations and molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) calculations to explore the structural requirements that participated in the effective binding of tributyrin and its catalytic triad (Ser165, Asp259 and His289) and identify detailed changes that occurred post mutation. Mut-LipKV1 bound favorably with tributyrin (−4.1 kcal/mol) and formed a single hydrogen bond with His289, at pH 9.0. Despite the incongruent docking analysis data, results of MD simulations showed configurations of both the tributyrin-Mut-LipKV1 (RMSD 0.3 nm; RMSF 0.05 − 0.3 nm) and the tributyrin-wildtype lipase KV1 (tributyrin-LipKV1) complexes (RMSD 0.35 nm; RMSF 0.05 − 0.4 nm) being comparably stable at pH 8.0. MM-PBSA analysis indicated that van der Waals interactions made the most contribution during the molecular binding process, with the Mut-LipKV1-tributyrin complex (−44.04 kcal/mol) showing relatively lower binding energy than LipKV1-tributyrin (−43.83 kcal/mol), at pH 12.0. All tributyrin-Mut-LipKV1 complexes displayed improved binding free energies over a broader pH range from 8.0 − 12.0, as compared to LipKV1-tributyrin. Future empirical works are thus, important to validate the improved alkaline-stability of Mut-LipKV1. In a nutshell, our research offered a considerable insight for further improving the alkaline tolerance of lipases. Communicated by Ramaswamy H. Sarma

Protein-Drug Binding Mode Prediction from the Apo-Protein Structure using a Molecular Dynamics-Based Pocket Generation Approach

In the drug development process, accurate prediction of the binding configuration between a target protein and a drug candidate is essential for optimizing the chemical structure of promising lead compounds. While we have successfully predicted the binding modes for several kinds of protein-drug pairs, prediction from the apo-protein structure remains to be challenging. This is because drug binding often induces significant protein conformational changes. Here, we report a computational workflow that incorporates a novel pocket generation method. First, the closed protein pocket is expanded by repeatedly filling virtual atoms during molecular dynamics (MD) simulations. Second, after ligand docking toward the sampled pocket structures, binding mode candidates are ranked by MD/Molecular Mechanics Poisson-Boltzmann Surface Area. We validated our workflow using the CDK2 kinase, which has an especially-closed ATP-binding pocket in the apo-form, and several inhibitors. The crystallographic pose coincided with the top-ranked docking pose for 59% (34/58) of all assessed compounds and was within the top five-ranked ones for 88% (51/58), while those estimated by a conventional prediction protocol were 9% (5/58) and 50% (29/58), respectively. Our study demonstrates that the prediction accuracy is significantly improved by preceding pocket expansion, leading to generation of conformationally-diverse binding mode candidates.

Identification of a novel binding mechanism of Quinoline based molecules with lactate dehydrogenase of Plasmodium falciparum

The deadliest disease caused by the Plasmodium species is malaria. Among other species, the infection caused by Plasmodium falciparum (Pf) is life-threatening. The biological function and three-dimensional structure of PfLDH and human LDH are very similar. Any treatment aiming to inhibit the PfLDH can also affect the activity of human LDH. This study aims to identify molecules that show high selectivity for PfLDH without having a profound effect on the activity of human LDH. In this study, 30 in-house synthesized Quinolines based molecules were docked with both PfLDH and human LDH. Based on molecular docking results, molecules 3j, 4b, 4h, 4m were showed selectivity towards PfLDH. All the four molecules had better binding affinity, ligand efficiency, lipophilic ligand efficiency, and torsion than Chloroquine (referenced inhibitor). Further, molecular dynamics (MD) simulations and molecular mechanics Poisson Boltzmann surface area (MM-PBSA) calculations were carried out to validate the docking results and to compare the stability of selected complexes against the complex with Chloroquine. MD simulations showed stable dynamic behavior of all the selected molecules in comparison to Chloroquine. Finally, on the basis of MM-PBSA analysis, molecule 3j was revealed novel binding mechanism and was selected as the potent and selective inhibitor of PfLDH.Communicated by Ramaswamy H. Sarma

Computational Modeling Explains the Multi Sterol Ligand Specificity of the N-Terminal Domain of Niemann-Pick C1-Like 1 Protein.

Niemann–Pick C1 like 1 (NPC1L1) is a sterol transporter expressed in the apical membrane of enterocytes and hepatocytes. NPC1L1 resembles the lysosomal NPC1 protein including an N-terminal domain (NTD), which binds a variety of sterols. The molecular mechanisms underlying this multiligand specificity of the NTD of NPC1L1 (NPC1L1–NTD) are not known. On the basis of the crystal structure of NPC1L1–NTD, we have investigated the structural details of protein–sterol interactions using molecular mechanics Poisson Boltzmann surface area calculations here. We found a good agreement between experimental and calculated binding affinities with similar ranking of various sterol ligands. We defined hydrogen bonding of sterol ligands via the 3′-β-hydroxy group inside the binding pose as instrumental in stabilizing the interaction. A leucine residue (LEU213) at the mouth of the binding pocket transiently opens to allow for the access of sterol into the binding pose. Our calculations also predict that NPC1L1–NTD binds polyene sterols, such as dehydroergosterol or cholestatrienol with high affinity, which validates their use in future experiments as close intrinsically fluorescent cholesterol analogs. A free energy decomposition and computational mutation analysis revealed that the binding of various sterols to NPC1L1–NTD depends critically on specific amino acid residues within the binding pocket. Some of these residues were previously detected as being relevant for intestinal cholesterol absorption. We show that clinically known mutations in the NPC1L1–NTD associated with lowered risk of coronary heart disease result in strongly reduced binding energies, providing a molecular explanation for the clinical phenotype.

Modeling antimalarial and antihuman African trypanosomiasis compounds: a ligand- and structure-based approaches.

This study examines the interaction of 137 antimalarial and antihuman African trypanosomiasis compounds [bis(2-aminoimidazolines), bisguanidinediphenyls and polyamines] on three different in vitro assays (Trypanosoma brucei rhodesiense (T.b.r.), Plasmodium falciparum (P.f.) and cytotoxicity-L6 cells). ΔTm values, wherever available, were also examined for the considered ligands. Eight DNA-ligand complexes and one DNA structure without ligand were selected from protein data bank (PDB) based on the structural similarity. Geometry optimization of all the considered ligands was carried out at the B3LYP/6-31G(d) level of theory. The AutoDock4 tool was utilized for the docking of these molecules at the minor groove of nine selected DNA crystal structures. We observed DT20, DA6, DT8 and DT19 residues generally interact with most of the considered ligands. Molecular dynamics simulations, molecular mechanics-generalized born surface area and molecular mechanics-Poisson Boltzmann surface area calculations indicate that the docked poses are generally stable and docked ligands do not show much deviation in the minor groove of DNA until 10ns simulation. Efficient and statistically significant quantitative structure-activity relationship models for T.b.r., P.f., C-L6 and ΔTm values were developed. All the generated models are internally and externally validated. We predicted a few ligands with significant IC50 values against P.f. based on the developed models. These results may help to design new and potent antimalarial and antihuman African trypanosomal compounds.

Prediction of novel inhibitors for Crotalus adamanteus l-amino acid oxidase by repurposing FDA-approved drugs: a virtual screening and molecular dynamics simulation investigation

One of the deadliest enzymes in the snake venom is l-amino acid oxidase (LAAO) which plays an important role in the pathophysiological effects during snake envenomation. Some effects of this enzyme on the human body are apoptosis, platelet aggregation, edema, hemorrhage, and cytotoxicity. Hence, inhibiting the enzyme activity to reduce its degradation effects is of great medical and pharmacological importance. On the other hand, drug repurposing is a process to find the new existing drug for a new medical indication. Since Crotalus adamanteus LAAO has no crystal structure in the protein data bank, first, its 3D structure was constructed by homology modeling using 1REO as the template and then modeled structure was evaluated by several algorithms. We screened the FDA-approved drugs by structure-based virtual screening, molecular dynamics (MD) simulation, and Molecular Mechanics Poisson Boltzmann Surface Area (MM/PBSA) to identify new inhibitors for the snake venom LAAO. Interestingly, docking results revealed that half of the hits belong to the propionic acid derivatives drugs. In addition, MD simulation was performed to assess the interaction profile of the docked protein–hits complexes. Meanwhile, Arg88, Gln112, Lys345, Trp356 form consistent hydrogen bond interactions with Dexketoprofen, Flurbiprofen, Ketoprofen, Morphine, and Citric acid during simulation. According to the results, each of the four compounds can be an appropriate inhibitor of LAAO and since our study was based on drug repurposing could be evaluated in phase II clinical trials.

Design, synthesis, biological evaluation and molecular dynamics studies of 4-thiazolinone derivatives as protein tyrosine phosphatase 1B (PTP1B) inhibitors

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin signaling pathway, and more and more studies have shown that it is a potential target for the treatment of type 2 diabetes mellitus (T2DM). In this study, 17 new 4-thiazolinone derivatives were designed and synthesized as novel PTP1B inhibitors, and ADMET prediction confirmed that these compounds were to be drug-like. In vitro enzyme activity experiments were performed on these compounds, and it was found that a plurality of compounds had good inhibitory activity and high selectivity against PTP1B protein. Among them, compound 7p exhibited the best inhibitory activity with an IC50 of 0.92 μM. The binding mode of compound 7p and PTP1B protein was explored, revealing the reason for its high efficiency. In addition, molecular dynamics simulations for the PTP1BWT and PTP1Bcomp#7p systems revealed the effects of compound 7p on PTP1B protein at the molecular level. In summary, the study reported for the first time that 4-thiazolinone derivatives as a novel PTP1B inhibitor had good inhibitory activity and selectivity for the treatment of T2DM, providing more options for the development of PTP1B inhibitors.AbbreviationsBBBblood-brain barrierCDC25Bcell division cycle 25 homolog BCYP2D6Cytochrome P450 2D6 bindingDCCMdynamic cross-correlation mapDSDiscovery StudioH bondhydrogen bondHIAhuman intestinal absorptionLARleukocyte antigen-related phosphataseMDmolecular dynamicsMEG-2maternal-effect germ-cell defective 2MM-PBSAmolecular mechanics Poisson Boltzmann surface area)PCAprincipal component analysisPDBProtein Data BankpNPPp–nitrophenyl phosphatePPBplasma protein bindingPTP1Bprotein tyrosine phosphotase 1BRMSDroot mean square deviationRMSFroot mean square fluctuationSHP-1src homologous phosphatase-1SHP-2src homologous phosphatase-2SPCsingle-point chargeTCPTPT cell protein tyrosine phosphataseT2DMType 2 diabetes mellitusVDWvan der WaalsCommunicated by Ramaswamy H. Sarma

Underlying mechanistic insights into the structural properties of melamine and uric acid complexes with compositional variation under ambient conditions

The structural properties of melamine-uric acid complexes (which are responsible for kidney stones) with compositional variations are examined using a series of classical molecular dynamics simulations. The preferential interaction parameters imply that melamine interacts more strongly with uric acid than with other melamine molecules present in the system, whereas uric acid preferentially interacts with other uric acid molecules rather than with melamine. The stronger interactions among uric acid molecules produce higher-order uric acid clusters, which “drag” neighboring melamine molecules to be added to a cluster. Determination of orientational preferences between aromatic planes reveals that π–π stacking is responsible for uric acid self-association but less significant for melamine-melamine and melamine-uric acid accumulation. Cluster structure analyses suggest that higher concentrations of melamine, uric acid, or both result in a large insoluble melamine-uric acid complex cluster. Molecular mechanics-Poisson Boltzmann surface area calculations give a negative binding energy, indicating favorable complexation between melamine and uric acid molecules. Moreover, the overall complexation energy [ΔG0(mel-mel)+ ΔG0(uri-uri)+ ΔG0(mel-uri)] is more negative than ΔG0bind(mel-uri). The lifetime of melamine dimers is quite low compared with those of uric acid-uric acid and melamine-uric acid dimers, resulting in a low percentage of larger clusters for melamine-melamine interaction and a significant percentage of higher-order melamine-uric acid and uric acid-uric acid clusters with longer lifetimes. Furthermore, melamine and uric acid form strong hydrogen bonds, and melamine-melamine interactions are dominated by hydrogen bonding, whereas uric acid forms only a small number of hydrogen bonds with other uric acid molecules.

A profound computational study to prioritize the natural compound inhibitors against the P. falciparum orotidine-5-monophosphate decarboxylase enzyme

In the current contribution, a multicomplex-based pharmacophore modeling approach was employed on the structural proteome of Plasmodium falciparum orotidine-5-monophosphate decarboxylase enzyme (PfOMPDC). Among the constructed pharmacophore models, the representative hypotheses were selected as the primary filter to screen the molecules with the complementary features responsible for showing inhibition. Thereafter, auxiliary evaluations were performed on the screened candidates via drug-likeness and molecular docking studies. Subsequently, the stability of the docked protein-ligand complexes was scrutinized by employing molecular dynamics simulations and molecular mechanics-Poisson Boltzmann surface area based free binding energy calculations. The stability the docked candidates was compared with the highly active crystallized inhibitor (3S9Y-FNU) to seek more potential candidates. All the docked molecules displayed stable dynamic behavior and high binding free energy in comparison to 3S9Y-FNU. The employed workflow resulted in the retrieval of five drug-like candidates with diverse scaffolds that may show inhibitory activity against PfOMPDC and could be further used as the novel scaffold to develop novel antimalarials.Communicated by Ramaswamy H. Sarma

Homology modeling, docking and structure-based virtual screening for new inhibitor identification of Klebsiella pneumoniae heptosyltransferase-III

Klebsiella pneumoniae (K. pneumoniae) is a Gram-negative opportunistic pathogen commonly associated with hospital-acquired infections that are often resistant even to antibiotics. Heptosyltransferase (HEP) belongs to the family of glycosyltransferase-B (GT-B) and plays an important in the synthesis of lipopolysaccharides (LPS) essential for the formation of bacterial cell membrane. HEP-III participates in the transfer of heptose sugar to the outer surface of bacteria to synthesize LPS. LPS truncation increases the bacterial sensitivity to hydrophobic antibiotics and detergents, making the HEP as a novel drug target. In the present study, we report the 3D homology model of K. pneumoniae HEP-III and its structure validation. Active site was identified based on similarities with known structures using Dali server, and structure-based pharmacophore model was developed for the active site substrate ADP. The generated pharmacophore model was used as a 3D search query for virtual screening of the ASINEX database. The hit compounds were further filtered based on fit value, molecular docking, docking scores, molecular dynamics (MD) simulations of HEP-III complexed with hit molecules, followed by binding free energy calculations using Molecular Mechanics-Poisson–Boltzmann Surface Area (MM-PBSA). The insights obtained in this work provide the rationale for design of novel inhibitors targeting K. pneumoniae HEP-III and the mechanistic aspects of their binding. Communicated by Ramaswamy H. Sarma