Rand protease is a serine protease that shared common characteristics with members of the MEROPS S8 subtilisin family. It is thermostable, highly stable in organic solvent and broad in specificity. Many structures of homologous protein solved by X-ray crystallography and NMR have been deposited to Protein Data Bank (PDB) which allowed this study to rely on structure prediction by deep learning to build three-dimensional (3D) structure of full length and mature Rand protease (flRP and mRP). In silico cysteine mutation to 7 predicted high affinity Ca2+ coordinating residues were introduced, and the mutants were subjected to molecular dynamics simulation to study its effect on flRP and mRP. MD simulation showed a marked increase in flexibility of the pro-peptide segment indicating the impact of single cysteine substitution at high affinity Ca2+ coordinating residues to autolysis of flRP. MD simulation for mRP reaffirmed the role of Ca2+ coordinating sites in providing stability to Rand protease. In addition, these residues also affect the autolysis, folding and hydrophobicity of RP. Essential dynamics observed large contribution of the first few eigenvectors of flRP, mRP and their high affinity Ca2+ coordinating residues mutants to the TMSF values which indicates that these values account for a large portion of the overall atomic fluctuations. These results have given a more comprehensive understanding on the role of cysteine substituted Ca2+ coordinating surface loop to the structure of flRP and mRP which are important in contributing to the structural stability of subtilisin.Communicated by Ramaswamy H. Sarma.
To explore the new mode of action and reduce side effects, making conjugates of existing drugs is becoming an attractive tool in the realm of medicinal chemistry. In this work, we exploited this approach and synthesized new conjugates to assess their activities against the enzymes involved in different pathological conditions. Specifically, we design and synthesized conjugates involving acetylsalicylic acid and sulfa drugs, validating the newly crafted conjugates using techniques like IR, 1HNMR, 13CNMR, and elemental analysis. These conjugates underwent assessment for their ability to inhibit cyclooxygenase-2 (COX-2), urease enzymes, and their anti-inflammatory potential. A competitive mode of urease inhibition was observed for acetylsalicylic acid conjugated with sulfanilamide, sulfacetamide, and sulfadiazine with IC50 of 2.49 ± 0.35 µM, 6.21 ± 0.28 µM, and 6.57 ± 0.44 µM, respectively. Remarkably, the acetylsalicylic acid-sulfamethoxazole conjugate exhibited exceptional anti-inflammatory activity, effectively curtailing induced edema by 83.7%, a result akin to the reference anti-inflammatory drug indomethacin's performance (86.8%). Additionally, it demonstrated comparable COX-2 inhibition (75.8%) to the reference selective COX-2 inhibitor celecoxib that exhibited 77.1% inhibition at 10 µM concentration. To deepen our understanding, we employed molecular docking techniques to predict the binding interactions of competitive inhibitors with COX-2 and urease receptors. Additionally, MD simulations were carried out, confirming the stability of inhibitor-target complexes throughout the simulation period, devoid of significant conformational changes. Collectively, our research underscores the potential of coupling approved medicinal compounds to usher in novel categories of pharmacological agents, holding promise for addressing a wide spectrum of pathological disorders involving COX-2 and urease enzymes.Communicated by Ramaswamy H. Sarma.
The infectious Nipah virus (NiV) is categorized into NiV-M (Malaysia) and NiV-B (Bangladesh) groups based on its genome comparison, pathogenicity, and mortality rate. The development of therapeutic molecules has used NiV-M-derived data in multiple studies than NiV-B. In continuation with this, the protein level investigation is also less explored to understand the interaction with therapeutic neutralizing antibodies for NiV-B. So, this study focuses on understanding the impact of NiV-B-specific mutations on the interaction of therapeutic neutralizing antibodies with the G protein. The population-based comparative analysis of NiV-B G protein sequences with NiV-M sequence identified twenty-six mutations. These predominantly polar mutations were then used to model the mutant protein (G_MT). In a comparative study, the G protein G_MT and reference protein G_WT (Malaysian origin) were subjected to a protein docking with neutralizing human monoclonal antibody HENV26. The binding affinity and the free binding energy of the glycoprotein in complex with G-WT and G_MT were calculated using PRODIGY and MM/PBSA tools respectively. Based on the PRODIGY report, G-WT showed stronger binding (-13.8 kcal/mol) compared to that of the G_MT (-9.0 kcal/mol) with the HENV26 antibody. The stability of the complexes was evaluated using MM/PBSA which showed higher binding energy with HENV26 for G_WT (-75.11 kcal/mol) in contrast to G_MT (-41.66 kcal/mol). The results indicate that the mutant G protein has a reduced ability to bind to neutralizing antibodies, resulting in a decreased effectiveness against strains carrying these mutations.Communicated by Ramaswamy H. Sarma.
Extracellular signal-regulated kinase 2 (ERK-2) is a serine/threonine protein kinase in eukaryotic cells and belongs to the mitogen-activated protein kinase (MAPK) family. An activated form of ERK-2 phosphorylates substrates in the nucleus or cytoplasm and causes specific proteins to be expressed or activated, regulating cell proliferation, differentiation and other functions. Caffeic acid (3,4 - dihydroxy cinnamic acid), as previously reported, directly interacts with ERK-2 and reduces its effects in vitro. It is also reported to have a variety of pharmacological effects, including anti-inflammatory, immunomodulatory, antioxidant and anticancer activities. In the current study, a deep-learning protocol was employed to develop effective 100 compounds by modifying the chemical structure of DHC to improve its inhibitory performance against ERK-2. Calculations of physicochemical properties for those compounds revealed that 20 compounds had drug scores better than DHC (≥ 80%). Following that, molecular docking calculations were performed on the selected compounds and DHC. The obtained data revealed that five compounds had docking scores better than DHC (≥ -5.9 kcal/mol). Moreover, data from molecular mechanics and the Poisson - Boltzmann surface area (MM/PBSA) binding energy over 200 ns MD simulation confirmed that Cmd-1 and Cmd-4 exhibited higher stability with ΔGbinding of -40.8 and -49.1 kcal/mol, respectively, which is better than DHC (-35.1 kcal/mol). Finally, various energetic and structural studies showed the high stability of the two generated compounds within the active site of ERK-2. This study highlights the potential use of Cmd-1 and Cmd-4 as promising anti-ERK-2 drug candidates.Communicated by Ramaswamy H. Sarma.
During last decades, 3,5-disubstituted-tetrahydro-2H-thiadiazine-2-thione scaffold remains the center of interest due to their ease of preparation, diverse range substituents at N-3 and N-5 positions, and profound biological activities. In the current study, a series of 3,5-disubstituted-tetrahydro-2H-thiadiazine-2-thiones were synthesized in good to excellent yield, and the structure of the compounds were confirmed by various spectroscopic techniques such as FTIR, 1H-NMR, 13C-NMR and Mass spectrometry, and finally evaluated against Leishmania major. Whereas, all the evaluated compounds (1-33), demonstrate potential leishmanicidal activities with IC50 values in the range of (1.30- 149.98 uM). Among the evaluated compounds such as 3, 4, 6, and 10 exhibited excellent leishmanicidal activities with IC50 values of (2.17 μM), (2.39 μM), (2.00 μM), and (1.39 μM), respectively even better than the standard amphotericin B (IC50 = 0.50) and pentamidine (IC50 = 7.52). In order to investigate binding interaction of the most active compounds, molecular docking study was conducted with Leishmania major. Further molecular dynamic simulation study was also carried out to assess the stability and correct binding of the most active compound 10, within active site of the Leishamania major. Likewise, the physiochemical properties, drug likeness, and ADMET of the most active compounds were investigated, it was found that none of the compounds violate Lipiniski's rule of five, which show that this class of compounds had enough potential to be used as drug candidate in near future.Communicated by Ramaswamy H. Sarma.
Virtual screening aims to identify and rank compounds with drug/lead-like properties based on their affinity for the protein target. We developed a methodology that integrates structure- and ligand-based screening approaches to enhance hit rates against the TDP1 protein within a database of anthraquinone and chalcone derivatives, followed by evaluation of prioritized compounds through molecular simulations. This technique is particularly useful for training set imbalances. Four screening methods were used: QSAR, pharmacophore, shape similarity, and docking. Each method was individually trained to score compounds, and the scores were fused to create parallel Z-score fusion. The QSAR models exhibited satisfactory R2 values (0.84 to 0.75), whereas the pharmacophoric and shape similarity models demonstrated excellent performance (ROC:0.82-0.88). Docking enrichment analysis identified 6N0D as the optimal TDP1 crystal structure (ROC = 0.73). Remarkably, the consensus scoring method surpassed other screening methods, achieving the highest ROC value of 0.98. Docking screening prioritized compounds with binding modes resembling the co-crystallized ligands, whereas MMGBSA, consensus, and docking produced dynamic simulations that were as stable as the co-crystallized ligands. Additionally, the QSAR-selected compounds exhibited binding modes similar to those of commercially available TDP1 inhibitors. In this study, a strong correlation was found between the inhibitory concentrations and binding energy values of commercialized TDP1 inhibitors, indicating that the top-ranked compounds are expected to have potent inhibitory effects in the nano-/micromolar range. The results of this study establish that consensus scoring can be used as an adaptable mainstay virtual screening methodology, pending subsequent experimental validation for affirmation.Communicated by Ramaswamy H. Sarma.
Dengue fever is now one of the major global health concerns particularly for tropical and sub-tropical countries. However, there has been no FDA approved medication to treat dengue fever. Researchers are looking into DENV NS5 RdRp protease as a potential therapeutic target for discovering effective anti-dengue agents. The aim of this study to discover dengue virus inhibitor from a set of five compounds from Momordica charantia L. using a series of in-silico approaches. The compounds were docked into the active area of the DENV-2 NS5 RdRp protease to obtain the hit compounds. The successful compounds underwent additional testing for a study on drug-likeness similarity. Our study obtained Momordicoside-I as a lead compound which was further exposed to the Cytochrome P450 (CYP450) toxicity analysis to determine the toxicity based on docking scores and drug-likeness studies. Moreover, DFT studies were carried out to calculate the thermodynamic, molecular orbital and electrostatic potential properties for the lead compound. Moreover, the lead compound was next subjected to molecular dynamic simulation for 200 ns in order to confirm the stability of the docked complex and the binding posture discovered during docking experiment. Overall, the lead compound has demonstrated good medication like qualities, non-toxicity, and significant binding affinity towards the DENV-2 RdRp enzyme.Communicated by Ramaswamy H. Sarma.
Breast cancer continues to be a major health issue for women all over the world. Cancer medications like cisplatin, which are widely used, still have negative side effects. The novel complex was created as a potential anticancer medication candidate that is both effective and safe, with few side effects. The Cu(II) complex using the prolinedithiocarbamate ligands was synthesized in situ. The Cu(II) complexes Characterization by UV-Vis, FT-IR spectroscopy and melting point determination, conductivity, and HOMO-LUMO were studied. Computational NMR spectrum analysis was performed. The interaction of Cu(II)prolineditiocarbamate complex with cancer cell target protein (MCF-7) was confirmed by molecular docking and molecular dynamic. The pharmacokinetic/ADMET properties were also performed on the complex. Results of the cytotoxic complex test against cancer cells (MCF-7) undergoing apoptosis with an IC50 value of 13.64 µg/mL showed high anticancer activity in MCF-7 cancer cells. The in-vivo data for Cu(II)prolineditiocarbamate complex was predicted using the Protox online tool with an LD50 value of 2500 mg/kg and belonging to the GHS toxicity class 5, which means the compound has a low acute toxicity effect. The Cu(II) prolineitiocarbamate complex may pave the way for the development of essential metal-based chemotherapy for the treatment of breast cancer.Communicated by Ramaswamy H. Sarma.
Dengue hemorrhagic fever (DHF) is severe dengue with a hallmark of vascular leakage. β-tryptase has been found to promote vascular leakage in DHF patients, which could be a potential target for DHF treatment. This study aims to develop a theoretical background for designing and selecting human β-tryptase inhibitors through computational studies. Thirty-four α-keto-[1,2,3]-oxadiazoles scaffold-based compounds were used to generate 2D-QSAR models and for molecular docking studies with β-tryptase (PDB Code 4A6L). In addition, molecular dynamics (MD) simulation and molecular mechanics generalised born surface area (MM-GBSA) analysis on the binding of the reported most active compound, compound 11e, towards β-tryptase were performed. Finally, a structure-based pharmacophore model was generated. The selected 2D-QSAR models have statistically proven good models by internal and external validation as well as the y-randomization test. The docking results of compound 11e showed lower CDOCKER energy than the 4A6L co-crystallised ligand and a similar binding pattern as the 4A6L co-crystallised ligand. From molecular dynamics simulation, 4A6L in compound 11e bound state has RMSD below 2 Å throughout the 500 ns simulation, indicating the docked complex is stable. Besides, MM-GBSA analysis suggested the 4A6L-compound 11e docked complex (-66.04 Kcal/mol) is structurally as stable as the 4A6L-native ligand co-crystallized structure (-66.84 Kcal/mol). The best pharmacophore model identified features included hydrogen bond acceptor, ionic interaction, hydrophobic interaction, and aromatic ring, which contribute to the inhibitory potency of a compound. This study supplied insight and knowledge for developing novel chemical compounds with improved inhibition of β-tryptase.Communicated by Ramaswamy H. Sarma.
Interleukin 17 F is a member of IL-17 cytokine family with a 50% structural homology to IL-17A and plays a significant role either alone or in combination with IL-17A towards inflammation in Rheumatoid arthritis (RA). A growing number of drugs targeting IL-17 pathway are being tested against population specific disease markers. The major objective of this research was to investigate the anti-inflammatory effect of Anakinra (an IL-1 R1 inhibitor) and Ustekinumab (an IL-12 and IL-23 inhibitor) by targeting IL17F. The three dimensional structures of IL17F was taken from PDB while structures of drugs were taken from PubChem database. Docking was performed using MOE and Schrodinger ligand docking software and binding energies, including s-score using London-dG fitness function and glide score using glide internal energy function, between drug and targets were compared. Furthermore, Protein-Drug complex were subjected to 150 ns Molecular Dynamics (MD) Simulations using Schrodinger's Desmond Module. Docking and MD simulation results suggest anakinra as a more potent IL17F inhibitor and forming a more structurally stable complex.Communicated by Ramaswamy H. Sarma.
14-3-3σ protein is one of the seven isoforms from the highly conserved eukaryotic 14-3-3 protein family. Downregulation of 14-3-3σ expression has been observed in various tumors. TRIM25 is responsible for the proteolytic degradation of 14-3-3σ, in which abrogation of TRIM25 suppressed tumor growth through 14-3-3σ upregulation. However, to date, the exact 14-3-3σ interacting residues of TRIM25 have yet to be resolved. Thus, this study attempts to identify the peptide binding sequence of TRIM25 on 14-3-3σ via both bioinformatics and biophysical techniques. Multiple sequence alignment of the CC domain of TRIM25 revealed five potential peptide binding sequences (Peptide 1-5). Nuclear magnetic resonance (NMR) assay (1H CPMG) identified Peptide 1 as an important sequence for binding to 14-3-3σ. Competition NMR assay suggested that Peptide 1 binds to the amphipathic pocket of 14-3-3σ with an estimated KD of 116.4 µM by isothermal titration calorimetry. Further in silico docking and molecular dynamics simulations studies proposed that Peptide 1 is likely to interact with Lys49, Arg56, Arg129, and Tyr130 residues at the amphipathic pocket of 14-3-3σ. These results suggest that Peptide 1 may serve as a biological probe or a template to design inhibitors of TRIM25-14-3-3σ interaction as a potentially novel class of anticancer agents.Communicated by Ramaswamy H. Sarma.
Among the enzymes required for the efficient utilisation of pectin is polygalacturonase. Saccharobesus litoralis harbours two polygalacturonases belonging to glycoside hydrolase family 28 (GH28). One of them, PGQ1, cleaved polygalacturonate exolytically at the non-reducing end into monomeric units. It was most active at 60 °C and pH 8, with Km and kcat values of 2.3 mg/ml and 6.4 s-1 respectively. Its homology model of a right-handed parallel β-helix core consisted of Asp297 as the general acid and either Asp276 or Asp298 as the general base. By inferring the substrate binding modes at the -1 and +1 subsites from known crystal structures, a hexagalacturonate could be docked into the highly electropositive binding cleft. Interestingly, while no residues were present in the vicinity to make up the +2 and +4 subsites, Arg361 and Arg430 could readily bind to the carboxyl groups of the galacturonates at the +3 and +5 subsites respectively. Structural comparison suggested that this binding pattern with missing subsites might be unique to closely related exopolygalacturonases. As S. litoralis grew much more slowly on extracellular galacturonate due to the lack of a transporter for the monosaccharide, PGQ1 probably functioned in the periplasm to help degrade oligopectates completely.Communicated by Ramaswamy H. Sarma.
Meyerozyma guilliermondii is a rare yeast pathogen contributing to the deadly invasive candidiasis. M. guilliermondii strain SO, as a promising protein expression host, showed 99% proteome similarity with the clinically isolated ATCC 6260 (type strain) in a recent comparative genomic analysis. However, their in vitro virulence features and in vivo pathogenicity were uncharacterized. This study aimed to characterize the in vitro and in vivo pathogenicity of M. guilliermondii strain SO and analyze its Als proteins (MgAls) via comprehensive bioinformatics approaches. M. guilliermondii strain SO showed lower and higher sensitivity towards β-mercaptoethanol and lithium, respectively than the avirulent S. cerevisiae but exhibited the same tolerance towards cell wall-perturbing Congo Red with C. albicans. With 7.5× higher biofilm mass, M. guilliermondii strain SO also demonstrated 75% higher mortality rate in the zebrafish embryos with a thicker biofilm layer on the chorion compared to the avirulent S. cerevisiae. Being one of the most important Candida adhesins, sequence and structural analyses of four statistically identified MgAls showed that MgAls1056 was predicted to exhibit the most conserved amyloid-forming regions, tandem repeat domain and peptide binding cavity (PBC) compared to C. albicans Als3. Favoured from the predicted largest ligand binding site and druggable pockets, it showed the highest affinity towards hepta-threonine. Non-PBC druggable pockets in the most potent virulence contributing MgAls1056 provide new insights into developing antifungal drugs targeting non-albicans Candida spp. Virtual screening of available synthetic or natural bioactive compounds and MgAls1056 deletion from the fungal genome should be further performed and validated experimentally.Communicated by Ramaswamy H. Sarma.
Phosphodiesterase 4 (PDE4) has been established as a drug target for inflammatory diseases of respiratory tract like asthma and chronic obstructive pulmonary disease. The selective inhibitors of PDE4B, a subtype of PDE4, are devoid of adverse effects like nausea and vomiting commonly associated with non-selective PDE4B inhibitors. This makes the development of PDE4B subtype selective inhibitors a desirable research goal. Thus, in the present study, molecular docking, molecular dynamic simulations and binding free energy were performed to explore potential selective PDE4B inhibitors based on ginger phenolic compounds. The results of docking studies indicate that some of the ginger phenolic compounds demonstrate higher selective PDE4B inhibition than existing selective PDE4B inhibitors. Additionally, 6-gingerol showed the highest PDE4B inhibitory activity as well as selectivity. The comparison of binding mode of PDE4B/6-gingerol and PDE4D/6-gingerol complexes revealed that 6-gingerol formed additional hydrogen bond and hydrophobic interactions with active site and control region 3 (CR3) residues in PDE4B, which were primarily responsible for its PDE4B selectivity. The results of binding free energy demonstrated that electrostatic energy is the primary factor in elucidating the mechanism of PDE4B inhibition by 6-gingerol. Dynamic cross-correlation studies also supported the results of docking and molecular dynamics simulation. Finally, a small library of molecules were designed based on the identified structural features, majority of designed molecules showed higher PDE4B selectivity than 6-gingerol. These results provide important structural features for designing new selective PDE4B inhibitors as anti-inflammatory drugs and promising candidates for synthesis and pre-clinical pharmacological investigations.
The goal of this work is to use a variety of in-silico techniques to identify anti-diabetic agents against DPP-IV enzyme from five main curcumin analogues. To produce the successful molecules, five main curcumin analogues were docked into the active site of DPP-IV enzyme. In comparison to the control molecule (Saxagliptin, -6.9 kcal/mol), all the compounds have the highest binding affinity (-7.6 to -7.7 kcal/mol) for the DPP-IV enzyme. These compounds underwent further testing for studies on drug-likeness, pharmacokinetics, and acute toxicity to see the efficacy and safety of compounds. To assess the stability of the docking complex and the binding posture identified during the docking experiment, our study got THC as the lead compound, which was then exposed to 200 ns of molecular dynamic simulation and PCA analysis. Additionally, DFT calculations were conducted to determine the thermodynamic, molecular orbital, and electrostatic potential characteristics of lead compound. Overall, the lead chemical has shown strong drug-like properties, is non-toxic, and has a sizable affinity for the DPP-IV enzyme.Communicated by Ramaswamy H. Sarma.
Molecular docking, molecular dynamics (MD) simulation, atomic force microscopy (AFM) and multi-spectroscopic techniques were selected to unveil the molecular association between the hepatitis B virus (HBV) inhibitor, entecavir (ETR), and the major blood plasma transporter, human serum albumin (HSA). The entire docking and simulation analyses recognized ETR binding to subdomain IIA (Site I) of HSA through hydrogen bonds, hydrophobic and van der Waals forces while maintaining the complex's stability throughout the 100 ns. A gradual lessening in the Stern-Volmer quenching constant (Ksv) with rising temperatures registered ETR-induced quenching of HBV fluorescence as static quenching, thus advising complexation between ETR and HSA. The further advocation of this conclusion was seen from a larger value of the biomolecular quenching rate constant ((kq) > 1010 M-1s-1), changes in the spectra (UV-Vis absorption) of HSA following ETR inclusion and ETR-induced swelling of HSA in the AFM results. The ETR appeared to bind to HSA with moderate affinity (Ka=1.87-1.19×104 M-1) at 290, 300 and 310 K. Significant alterations in the protein's secondary and tertiary structures, including changes in the protein's Tyr/Trp microenvironment, were also detected by circular dichroism and three-dimensional fluorescence spectra when the protein was bound to ETR. The findings of the drug displacement study backed the docking results of Site I as ETR's preferred binding site in HSA.Communicated by Ramaswamy H. Sarma.
Cancer, an incurable global disease, demands urgent anti-cancer drug development. Marine alkaloids like Manzamine and its derivatives show promise as RSK inhibitors against cancer cell invasion. Replacing the hydrogen at the 8-position of Manzamine A with a hydroxyl group has been shown to significantly enhance its biological activity. In this article, we designed various functional group compounds (A1-A21) substituted 8-Manzamine A by docking, MM-GBSA, molecular dynamics (MD) simulation, and well-tempered metadynamics (WT-MetaD) simulations to evaluate their potential as RSK1 inhibitors. Ligands A1-A21 were docked in the RSK1 N-terminal kinase domain (PDB ID: 2Z7Q) using the Glide module. The calculation of binding energy was performed using Prime MM-GB/SA, while MD simulations were conducted with the Desmond module of Schrodinger suite 2023. Compound A5 exhibits the highest G-score (-7.01) compared to 8-Hydroxymanzamine A (-6.08). Additionally, compounds A6, A10, A12, A17, A11, A4, and A13 demonstrate increased activity against RSK1 when compared to both 8-Hydroxymanzamine A and Manzamine A. Residues LEU68, VAL76, LEU141, PHE143, LEU144, PHE150, ASP148, GLU191, and LEU194 of RSK1 protein play a key role in binding with ligands. An MD simulation of Compound A5 was carried out to explore the dynamic interactions within the protein-ligand complex. Furthermore, WT-MetaD simulations validated the docking study results and identified the most energetically favored conformations for the A5/RSK1 complex. Ligands A5, A6, A10, A12, A17, A11, A4, and A13, featuring diverse functional groups and good Glide scores, may have the potential for significant RSK1 activity and merit further development.Communicated by Ramaswamy H. Sarma.
The human serotonin transporters (hSERTs) are neurotransmitter sodium symporters of the aminergic G protein-coupled receptors, regulating the synaptic serotonin and neuropharmacological processes related to neuropsychiatric disorders, notably, depression. Selective serotonin reuptake inhibitors (SSRIs) such as fluoxetine and (S)-citalopram are competitive inhibitors of hSERTs and are commonly the first-line medications for major depressive disorder (MDD). However, treatment-resistance and unpleasant aftereffects constitute their clinical drawbacks. Interestingly, vilazodone emerged with polypharmacological (competitive and allosteric) inhibitions on hSERTs, amenable to improved efficacy. However, its application usually warrants adjuvant/combination therapy, another subject of critical adverse events. Thus, the discovery of alternatives with polypharmacological potentials (one-drug-multiple-target) and improved safety remains essential. In this study, carbazole analogues from chemical libraries were explored using docking and molecular dynamics (MD) simulation. Selectively, two IBScreen ligands, STOCK3S-30866 and STOCK1N-37454 predictively bound to the active pockets and expanded boundaries (extracellular vestibules) of the hSERTs more potently than vilazodone and (S)-citalopram. For instance, the two ligands showed docking scores of -9.52 and -9.59 kcal/mol and MM-GBSA scores of -92.96 and -65.66 kcal/mol respectively compared to vilazodone's respective scores of -7.828 and -59.27 against the central active site of the hSERT (PDB 7LWD). Similarly, the two ligands also docked to the allosteric pocket (PDB 5I73) with scores of -8.15 and -8.40 kcal/mol and MM-GBSA of -96.14 and -68.46 kcal/mol whereas (S)-citalopram has -6.90 and -69.39 kcal/mol respectively. The ligands also conferred conformational stability on the receptors during 100 ns MD simulations and displayed interesting ADMET profiles, representing promising hSERT modulators for MDD upon experimental validation.Communicated by Ramaswamy H. Sarma.
The human angiotensin-converting enzyme 2 (ACE-2) receptor is a metalloenzyme that plays an important role in regulating blood pressure by modulating angiotensin II. This receptor facilitates SARS-CoV-2 entry into human cells via receptor-mediated endocytosis, causing the global COVID-19 pandemic and a major health crisis. Kelulut honey (KH), one of Malaysian honey recently gained attention for its distinct flavour and taste while having many nutritional and medicinal properties. Recent study demonstrates the antiviral potential of KH against SARS-CoV-2 by inhibiting ACE-2 in vitro, but the bioactive compound pertaining to the ACE-2 inhibition is yet unknown. An ensemble docking-based virtual screening was employed to screen the phytochemical compounds from KH with high binding affinity against the 10 best representative structures of ACE-2 that mostly formed from MD simulation. From 110 phytochemicals previously identified in KH, 27 compounds passed the ADMET analysis and proceeded to docking. Among the docked compound, SDC and FMN consistently exhibited strong binding to ACE-2's active site (-9.719 and -9.473 kcal/mol) and allosteric site (-7.305 and -7.464 kcal/mol) as compared to potent ACE-2 inhibitor, MLN 4760. Detailed trajectory analysis of MD simulation showed stable binding interaction towards active and allosteric sites of ACE-2. KH's compounds show promise in inhibiting SARS-CoV-2 binding to ACE-2 receptors, indicating potential for preventive use or as a supplement to other COVID-19 treatments. Additional research is needed to confirm KH's antiviral effects and its role in SARS-CoV-2 therapy, including prophylaxis and adjuvant treatment with vaccination.Communicated by Ramaswamy H. Sarma.
Aquaporin is a water channel protein that facilitates the movement of water across the cell membrane. Aquaporin from the Antarctic region has been noted for its psychrophilic properties and its ability to perform at a lower temperature but there remains limited understanding of the water mechanism of Antarctic Pseudomonas sp. strain AMS3 However, studies regarding aquaporin isolated from psychrophilic Pseudomonas sp. are still scattered. Recently, the genome sequence of an Antarctic Pseudomonas sp. strain AMS3 revealed a gene sequence encoding for a putative aquaporin designated as AqpZ1 AMS3. In this study, structure analysis and a molecular dynamics (MD) simulation of a predicted model of a fully hydrated aquaporin tetramer embedded in a lipid bilayer was performed at different temperatures for structural flexibility and stability analysis. The MD simulation results revealed that the structures were able to remain stable at low to medium temperatures. The protein was observed to have high flexibility in the loop region as compared to the helices region throughout the simulated temperatures. The selectivity filter and NPA motifs play a major role in solute selectivity and the pore radius of the protein. The structural and functional characterization of this psychrophilic aquaporin provides new insights for the future applications of this protein.Communicated by Ramaswamy H. Sarma.