Case report: A novel CACNA1S mutation associated with hypokalemic periodic paralysis

Background

Hypokalemic periodic paralysis (HypoKPP) is a rare neuromuscular genetic disorder causing recurrent episodes of flaccid paralysis. Most cases are associated with CACNA1S mutation, causing defect of calcium channel and subsequent impairment of muscle functions. Due to defined management approaches early diagnosis is crucial for promptly treatment and prevention new attacks.

Materials and methods

We report a case of HypoKPP associated with previously unreported mutation in CACNA1S gene (p.R900M). Molecular modeling of CaV1.1 was applied to evaluate its pathogenicity.

Results

As a patient referred between attacks neurological status, laboratory and neurophysiological examination were unremarkable. Molecular modeling predicted that the p.R900M mutation affects the process of calcium channels activation.

Conclusion

Novel CACNA1S mutation, associated with HypoKPP was identified. Monte-Carlo energy minimization of the CaV1.1 model supported the association of this mutation with this disease.

In Silico RNAseq and Biochemical Analyses of Glucose-6-Phosphate Dehydrogenase (G6PDH) from Sweet Pepper Fruits: Involvement of Nitric Oxide (NO) in Ripening and Modulation

Pepper (Capsicum annuum L.) fruit is a horticultural product consumed worldwide which has great nutritional and economic relevance. Besides the phenotypical changes that pepper fruit undergo during ripening, there are many associated modifications at transcriptomic, proteomic, biochemical, and metabolic levels. Nitric oxide (NO) is a recognized signal molecule that can exert regulatory functions in diverse plant processes including fruit ripening, but the relevance of NADPH as a fingerprinting of the crop physiology including ripening has also been proposed. Glucose-6-phosphate dehydrogenase (G6PDH) is the first and rate-limiting enzyme of the oxidative phase of the pentose phosphate pathway (oxiPPP) with the capacity to generate NADPH. Thus far, the available information on G6PDH and other NADPH-generating enzymatic systems in pepper plants, and their expression during the ripening of sweet pepper fruit, is very scarce. Therefore, an analysis at the transcriptomic, molecular and functional levels of the G6PDH system has been accomplished in this work for the first time. Based on a data-mining approach to the pepper genome and fruit transcriptome (RNA-seq), four G6PDH genes were identified in pepper plants and designated CaG6PDH1 to CaG6PDH4, with all of them also being expressed in fruits. While CaG6PDH1 encodes a cytosolic isozyme, the other genes code for plastid isozymes. The time-course expression analysis of these CaG6PDH genes during different fruit ripening stages, including green immature (G), breaking point (BP), and red ripe (R), showed that they were differentially modulated. Thus, while CaG6PDH2 and CaG6PDH4 were upregulated at ripening, CaG6PDH1 was downregulated, and CaG6PDH3 was slightly affected. Exogenous treatment of fruits with NO gas triggered the downregulation of CaG6PDH2, whereas the other genes were positively regulated. In-gel analysis using non-denaturing PAGE of a 50-75% ammonium-sulfate-enriched protein fraction from pepper fruits allowed for identifying two isozymes designated CaG6PDH I and CaG6PDH II, according to their electrophoretic mobility. In order to test the potential modulation of such pepper G6PDH isozymes, in vitro analyses of green pepper fruit samples in the presence of different compounds including NO donors (S-nitrosoglutathione and nitrosocysteine), peroxynitrite (ONOO-), a hydrogen sulfide (H2S) donor (NaHS, sodium hydrosulfide), and reducing agents such as reduced glutathione (GSH) and L-cysteine (L-Cys) were assayed. While peroxynitrite and the reducing compounds provoked a partial inhibition of one or both isoenzymes, NaHS exerted 100% inhibition of the two CaG6PDHs. Taken together these data provide the first data on the modulation of CaG6PDHs at gene and activity levels which occur in pepper fruit during ripening and after NO post-harvest treatment. As a consequence, this phenomenon may influence the NADPH availability for the redox homeostasis of the fruit and balance its active nitro-oxidative metabolism throughout the ripening process.

Synthesis and Biological Activity of a New Indenoisoquinoline Copper Derivative as a Topoisomerase I Inhibitor

Topoisomerases are interesting targets in cancer chemotherapy. Here, we describe the design and synthesis of a novel copper(II) indenoisoquinoline complex, WN198. The new organometallic compound exhibits a cytotoxic effect on five adenocarcinoma cell lines (MCF-7, MDA-MB-231, HeLa, HT-29, and DU-145) with the lowest IC50 (0.37 ± 0.04 μM) for the triple-negative MDA-MB-231 breast cancer cell line. Below 5 µM, WN198 was ineffective on non-tumorigenic epithelial breast MCF-10A cells and Xenopus oocyte G2/M transition or embryonic development. Moreover, cancer cell lines showed autophagy markers including Beclin-1 accumulation and LC3-II formation. The DNA interaction of this new compound was evaluated and the dose-dependent topoisomerase I activity starting at 1 μM was confirmed using in vitro tests and has intercalation properties into DNA shown by melting curves and fluorescence measurements. Molecular modeling showed that the main interaction occurs with the aromatic ring but copper stabilizes the molecule before binding and so can putatively increase the potency as well. In this way, copper-derived indenoisoquinoline topoisomerase I inhibitor WN198 is a promising antitumorigenic agent for the development of future DNA-damaging treatments.

Statistical Mechanics Metrics in Pairing and Parsing In Silico and Phenotypic Data of a Novel Genetic NFκB1 (c.T638A) Variant

(1) Background: Mutations in NFκB1, a transcriptional regulator of immunomodulating proteins, are a known cause of inborn errors of immunity. Our proband is a 22-year-old male with a diagnosis of common variable immunodeficiency (CVID), cytopenias with massive splenomegaly, and nodular regenerative hyperplasia of the liver. Genetic studies identified a novel, single-point mutation variant in NFκB1, c. T638A p. V213E. (2) Methods: Next-generation panel sequencing of the patient uncovered a novel single-point mutation in the NFκB1 gene that was modeled using the I-TASSER homology-modeling software, and molecular dynamics were assessed using the YASARA2 software (version 20.14.24). (3) Results: This variant replaces valine with glutamic acid at position 213 in the NFκB1 sequence. Molecular modeling and molecular dynamic studies showed altered dynamics in and around the rel homology domain, ankyrin regions, and death domain of the protein. We postulate that these changes alter overall protein function. (4) Conclusions: This case suggests the pathogenicity of a novel variant using protein-modeling techniques and molecular dynamic simulations.

Exploring the Interaction of New Pyridoquinazoline Derivatives with G-Quadruplex in the c-MYC Promoter Region

Novel amino-substituted pyridoquinazolinone derivatives have been designed and synthesized as potential c-MYC G-quadruplex (G4) ligands, employing an efficient methodology. All the new compounds exhibited moderate to good antiproliferative activity against the human osteosarcoma U2OS cell line. NMR and docking experiments revealed that the recently synthesized compounds interact with the Pu22 G-quadruplex in the c-MYC promoter region, establishing a 2:1 complex, with each molecule positioned over the tetrads at the 3'- and 5'-ends.

Curcumin-sulfobutyl-ether beta cyclodextrin inclusion complex: preparation, spectral characterization, molecular modeling, and antimicrobial activity

Urinary tract infections (UTIs) caused by Gram-negative bacteria E. coli is responsible for 80-90% of uncomplicated cases in women. The increased prevalence of antibiotic resistance has made the management of UTIs more challenging. Plant-derived compounds have long been used to treat various diseases, and constitute an alternative to antibiotic resistance. Curcumin (CUR), a naturally occurring polyphenolic phytoconstituent obtained from Curcuma longa is endowed with diverse medicinal properties. The present study aims to form a complex of CUR with Sulfobutyl ether-β-cyclodextrin (SBEβCD) to overcome the poor solubility and bioavailability of CUR and to evaluate the antimicrobial activity of CUR-SBEβCD. Phase solubility studies and spectral characterization showed the entrapment of CUR in the SBEβCD cavity. In silico docking studies performed to investigate the complexation process of CUR with SBEβCD, revealed that the methoxy group and OH group of CUR interacted with SBEβCD. The cytotoxicity and HET-CAM assays confirmed that CUR-SBEβCD was non-irritant. The prepared complex investigated with the disc diffusion method showed antimicrobial activity with a zone of inhibition (ZOI) of 13 mm against Escherichia coli (E. coli) and 11.5 mm against Staphylococcus aureus (S. aureus) whereas CUR alone did not show any ZOI. It can be concluded that prepared CUR-SBEβCD demonstrated superior antimicrobial activity and therefore can be a promising alternative for the treatment of UTIs.Communicated by Ramaswamy H. Sarma.

In-silico prediction, characterization, molecular docking and dynamic simulation studies for screening potential fungicides against leaf rust of Triticum aestivum

Triticum aestivum is an important crop worldwide, which is a large source of food grain. T.aestivum demands on developed countries will grow every year, this increase in the demand is profoundly serious especially in the light climate change which would lead to a 29% reduction in final productivity. Rust fungus attacks the T.aestivum, specifically newly planted T.aestivum plants, which block the vascular system, stun, and finally damage grain and tillers. In present study we predict the 3D structure then find the binding pocket and conserved domains for MAPkinase-1 of Puccinia triticina. After that, screen the FungiPAD, PubChem, NPAtlas databases by physicochemical properties, docking, clustering, ADME (Absorption, distribution, metabolism, and excretion) and PAINS (pan assay interference compounds) filter analysis. Through this screening process screen the nine compounds, which are benzovindiflupyr, furametpyr, isopyrazam, fenaminstrobin, and flumorph from Fungicide database: zoxamide, vinclozolin, pentachloronitrobenzene, and dithianon from PubChem database, based on the binding energy, clustering, ADME and PAINS analysis. All these nine compounds bind in the same pocket and show the same pattern of interaction. Among these nine compounds, select the two compounds (PubChem:122087 (-6.96 kcal/mol) and FDBD02904 (-8.62 kcal/mol)) based on binding energy for 100 ns MD simulation and free energy calculation. MD simulation shows stability throughout the simulation, and it shows the sable interaction when compounds bind to the MAPKinase 1 protein which may help to protein kinase pathways in plant defense response. This result helps to design alternative fungicide against the wheat rust disease.Communicated by Ramaswamy H. Sarma.

Non-Covalent Interactions between dUTP C5-Substituents and DNA Polymerase Decrease PCR Efficiency

The approach based on molecular modeling was developed to study dNTP derivatives characterized by new polymerase-specific properties. For this purpose, the relative efficiency of PCR amplification with modified dUTPs was studied using Taq, Tth, Pfu, Vent, Deep Vent, Vent (exo-), and Deep Vent (exo-) DNA polymerases. The efficiency of PCR amplification with modified dUTPs was compared with the results of molecular modeling using the known 3D structures of KlenTaq polymerase-DNA-dNTP complexes. The dUTPs were C5-modified with bulky functional groups (the Cy5 dye analogs) or lighter aromatic groups. Comparing the experimental data and the results of molecular modeling revealed the decrease in PCR efficiency in the presence of modified dUTPs with an increase in the number of non-covalent bonds between the substituents and the DNA polymerase (about 15% decrease per one extra non-covalent bond). Generalization of the revealed patterns to all the studied polymerases of the A and B families is discussed herein. The number of non-covalent bonds between the substituents and polymerase amino acid residues is proposed to be a potentially variable parameter for regulating enzyme activity.

Characterization of erythroferrone oligomerization and its impact on BMP antagonism

Hepcidin, a peptide hormone that negatively regulates iron metabolism, is expressed by bone morphogenetic protein (BMP) signaling. Erythroferrone (ERFE) is an extracellular protein that binds and inhibits BMP ligands, thus positively regulating iron import by indirectly suppressing hepcidin. This allows for rapid erythrocyte regeneration after blood loss. ERFE belongs to the C1Q/TNF related protein (CTRP) family and is suggested to adopt multiple oligomeric forms: a trimer, a hexamer, and a high molecular weight species. The molecular basis for how ERFE binds BMP ligands and how the different oligomeric states impact BMP inhibition are poorly understood. In this study, we demonstrated that ERFE activity is dependent on the presence of stable dimeric or trimeric ERFE, and that larger species are dispensable for BMP inhibition. Additionally, we used an in-silico approach to identify a helix, termed the ligand binding domain (LBD), that was predicted to bind BMPs and occlude the type I receptor pocket. We provide evidence that the LBD is crucial for activity through luciferase assays and surface plasmon resonance (SPR) analysis. Our findings provide new insight into how ERFE oligomerization impacts BMP inhibition, while identifying critical molecular features of ERFE essential for binding BMP ligands.

Redesigning methionyl-tRNA synthetase for β-methionine activity with adaptive landscape flattening and experiments

Amino acids (AAs) with a noncanonical backbone would be a valuable tool for protein engineering, enabling new structural motifs and building blocks. To incorporate them into an expanded genetic code, the first, key step is to obtain an appropriate aminoacyl-tRNA synthetase. Currently, directed evolution is not available to optimize AAs with noncanonical backbones, since an appropriate selective pressure has not been discovered. Computational protein design (CPD) is an alternative. We used a new CPD method to redesign MetRS and increase its activity towards β-Met, which has an extra backbone methylene. The new method considered a few active site positions for design and used a Monte Carlo exploration of the corresponding sequence space. During the exploration, a bias energy was adaptively learned, such that the free energy landscape of the apo enzyme was flattened. Enzyme variants could then be sampled, in the presence of the ligand and the bias energy, according to their β-Met binding affinities. Eighteen predicted variants were chosen for experimental testing; 10 exhibited detectable activity for β-Met adenylation. Top predicted hits were characterized experimentally in detail. Dissociation constants, catalytic rates, and Michaelis constants for both α-Met and β-Met were measured. The best mutant retained a preference for α-Met over β-Met; however, the preference was reduced, compared to the wildtype, by a factor of 29. For this mutant, high resolution crystal structures were obtained in complex with both α-Met and β-Met, indicating that the predicted, active conformation of β-Met in the active site was retained.

Computational Prediction of 3,5-Diaryl-1H-Pyrazole and spiropyrazolines derivatives as potential acetylcholinesterase inhibitors for alzheimer disease treatment by 3D-QSAR, molecular docking, molecular dynamics simulation, and ADME-Tox

The efficacy of 40 synthesized variants of 3,5-diaryl-1H-pyrazole and spiropyrazoline' derivatives as acetylcholinesterase inhibitors is verified using a quantitative three-dimensional structure-activity relationship (3D-QSAR) by comparative molecular field analysis (CoMFA) and molecular similarity index analysis (CoMSIA) models. In this research, different field models proved that CoMSIA/SE model is the best model with high predictive power compared to several models (Qved2 = O.65; R2 = 0.980; R2test = 0.727). Also, contour maps produced by CoMSIA/SE model have been employed to prove the key structural needs of the activity. Consequently, six new compounds have been generated. Among these compounds, M4 and M5 were the most active but remained toxic and had poor absorption capacities. While the M1, M2, M3 and M6 remained highly active while respecting ADMET's characteristics. Molecular docking results showed compound M2 better with acetylcholinesterase than compound 22. The interactions are classical hydrogen bonding with residues TYR:124, TYR:72, and SER:293, which play a critical role in the biological activity as AChE inhibitors. MD results confirmed the docking results and showed that compound M2 had satisfactory stability with (ΔGbinding = -151.225 KJ/mol) in the active site of AChE receptor compared with compound 22 (ΔGbinding = -133.375 KJ/mol). In addition, both compounds had good stability regarding RMSD, Rg, and RMSF. The previous results show that the newly designed compound M2 is more active in the active site of AChE receptor than compound 22.Communicated by Ramaswamy H. Sarma.

Exploring the interactive mechanism of acarbose with the amylase SusG in the starch utilization system of the human gut symbiont Bacteroides thetaiotaomicron through molecular modeling

The α-amylase, SusG, is a principal component of the Bacteroides thetaiotaomicron (Bt) starch utilization system (Sus) used to metabolize complex starch molecules in the human gastrointestinal (GI) tract. We previously reported the non-microbicidal growth inhibition of Bt by the acarbose-mediated arrest of the Sus as a potential therapeutic strategy. Herein, we report a computational approach using density functional theory (DFT), molecular docking, and molecular dynamics (MD) simulation to explore the interactive mechanism between acarbose and SusG at the atomic level in an effort to understand how acarbose shuts down the Bt Sus. The docking analysis reveals that acarbose binds orthosterically to SusG with a binding affinity of -8.3 kcal/mol. The MD simulation provides evidence of conformational variability of acarbose at the active site of SusG and also suggests that acarbose interacts with the main catalytic residues via a general acid-base double-displacement catalytic mechanism. These results suggest that small molecule competitive inhibition against the SusG protein could impact the entire Bt Sus and eliminate or reduce the system's ability to metabolize starch. This computational strategy could serve as a potential avenue for structure-based drug design to discover other small molecules capable of inhibiting the Sus of Bt with high potency, thus providing a holistic approach for selective modulation of the GI microbiota.

Prediction of protein interactions is essential for studying biomolecular mechanisms

Structural characterization of protein interactions is essential for our ability to understand and modulate physiological processes. Computational approaches to modeling of protein complexes provide structural information that far exceeds capabilities of the existing experimental techniques. Protein structure prediction in general, and prediction of protein interactions in particular, has been revolutionized by the rapid progress in Deep Learning techniques. The work of Schweke et al. (Proteomics 2023, 23, 2200323) presents a community-wide study of an important problem of distinguishing physiological protein-protein complexes/interfaces (experimentally determined or modeled) from non-physiological ones. The authors designed and generated a large benchmark set of physiological and non-physiological homodimeric complexes, and evaluated a large set of scoring functions, as well as AlphaFold predictions, on their ability to discriminate the non-physiological interfaces. The problem of separating physiological interfaces from non-physiological ones is very difficult, largely due to the lack of a clear distinction between the two categories in a crowded environment inside a living cell. Still, the ability to identify key physiologically significant interfaces in the variety of possible configurations of a protein-protein complex is important. The study presents a major data resource and methodological development in this important direction for molecular and cellular biology.

Antiviral evaluation of 1,4-disubstituted-1,2,3-triazole derivatives against Chikungunya virus

Aim

This work aimed to investigate the antiviral activity of two 1,4-disubstituted-1,2,3-triazole derivatives (1 and 2) against Chikungunya virus (CHIKV) replication.

Materials & methods

Cytotoxicity was analyzed using colorimetric assays and the antiviral potential was evaluated using plaque assays and computational tools.

Results

Compound 2 showed antiviral activity against CHIKV 181-25 in BHK-21 and Vero cells. Also, this compound presented a higher activity against CHIKV BRA/RJ/18 in Vero cells, like compound 1. Compound 2 exhibited virucidal activity and inhibited virus entry while compound 1 inhibited virus release. Molecular docking suggested that these derivatives inhibit nsP1 protein while compound 1 may also target capsid protein.

Conclusion

Both compounds exhibit promising antiviral activity against CHIKV by blocking different steps of virus replication.

Synthesis of novel furan-based chalcone derivatives as anti-tuberculosis agents: in vitro, cytotoxicity assessment and in silico

Background: The aim of the study is to identify a novel furan-based chalcone derivative as potent inhibitor against the H37Rv strain. Materials & methods: The in silico pharmacokinetic characteristics, toxicity tests, molecular modeling, chemical synthesis and minimum inhibitory concentration (MIC; IC50) were carried out to evaluate the antitubercular potential of the synthesized furan-based chalcone analogues against H37Rv. Results & conclusion: Among the ten target compounds synthesized, DF02, DF05 and DF07 had MIC values of 1.6 μg/ml equivalent to isoniazid and DF10 showed MIC values of 3.25 μg/ml which is equipotent to pyrazinamide. All the other compounds had optimal concentrations 6.25-100 μg/ml against the H37Rv strain. Compounds DF02 and DF10 were further evaluated for cytotoxicity assay performed using HeLa cell lines.

Experimental and computational evaluation of dipeptidyl peptidase III inhibitors based on quinazolinone-Schiff's bases

Dipeptidyl peptidase III (DPP III) is a zinc-dependent enzyme that sequentially hydrolyzes biologically active peptides by cleaving dipeptides from their N-termini. Although its fundamental role is not been fully elucidated, human DPP III (hDPP III) has been recognized in several pathophysiological processes of interest for drug development. In this article 27 quinazolinone-Schiff's bases were studied for their inhibitory activity against hDPP III combining an in vitro experiment with a computational approach. The biochemical assay showed that most compounds exhibited inhibitory activity at the 100 μM concentration. The best QSAR model included descriptors from the following 2D descriptor groups: information content indices, 2D autocorrelations, and edge adjacency indices. Five compounds were found to be the most potent inhibitors with IC50 values below 10 µM, while molecular docking predicted that these compounds bind to the central enzyme cleft and interact with residues of the substrate binding subsites. Molecular dynamics simulations of the most potent inhibitor (IC50=0.96 µM) provided valuable information explaining the role of PHE109, ARG319, GLU327, GLU329, and ILE386 in the mechanism of the inhibitor binding and stabilization. This is the first study that gives insight into quinazolinone-Schiff's bases binding to this metalloenzyme.Communicated by Ramaswamy H. Sarma.

Changes in Hemoglobin Properties in Complex with Glutathione and after Glutathionylation

Hemoglobin is the main protein of red blood cells that provides oxygen transport to all cells of the human body. The ability of hemoglobin to bind the main low-molecular-weight thiol of the cell glutathione, both covalently and noncovalently, is not only an important part of the antioxidant protection of red blood cells, but also affects its affinity for oxygen in both cases. In this study, the properties of oxyhemoglobin in complex with reduced glutathione (GSH) and properties of glutathionylated hemoglobin bound to glutathione via an SS bond were characterized. For this purpose, the methods of circular dichroism, Raman spectroscopy, infrared spectroscopy, tryptophan fluorescence, differential scanning fluorimetry, and molecular modeling were used. It was found that the glutathionylation of oxyhemoglobin caused changes in the secondary structure of the protein, reducing the alpha helicity, but did not affect the heme environment, tryptophan fluorescence, and the thermostability of the protein. In the noncovalent complex of oxyhemoglobin with reduced glutathione, the secondary structure of hemoglobin remained almost unchanged; however, changes in the heme environment and the microenvironment of tryptophans, as well as a decrease in the protein's thermal stability, were observed. Thus, the formation of a noncovalent complex of hemoglobin with glutathione makes a more significant effect on the tertiary and quaternary structure of hemoglobin than glutathionylation, which mainly affects the secondary structure of the protein. The obtained data are important for understanding the functioning of glutathionylated hemoglobin, which is a marker of oxidative stress, and hemoglobin in complex with GSH, which appears to deposit GSH and release it during deoxygenation to increase the antioxidant protection of cells.

Molecular docking and molecular dynamics studies of natural products unravel potential inhibitors against OmpA of Acinetobacter baumannii

Emerging antimicrobial resistance has highlighted the need to design more effective antibiotics to treat deadly bacterial infections. Acinetobacter baumannii's outer membrane protein A (OmpA) is a critical virulence component involved in biofilm formation, immunomodulation, and antibiotic resistance, which characterizes it as a potential therapeutic target. The present study aimed to screen the natural product database (>1,00,000) to identify the potential inhibitor against OmpA. Molecular docking studies revealed that 10 compounds had good docking scores (≤ -7 kcal/mol) compared to the reported inhibitor epiestriol (-3.079). Further, these 10 compounds were subjected to ADME analysis and MMGBSA analysis. Based on MMGBSA results, we selected 5 compounds [NP-1 (MolPort-039-337-117), NP-5(MolPort-019-932-973), NP-6 (MolPort-005-948-336), NP-8(MolPort-042-673-978) and NP-9(MolPort-042-673-766)] with high binding affinity. Molecular dynamics simulation found that NP-5, NP-8, and NP-9 were stable after analysing their RMSD, RMSF, the radius of gyration, and hydrogen interactions of complexes. Our study revealed that NP-5, NP-8, and NP-9 bind perfectly with OmpA and can act as its potential inhibitors. The results of this study imply that the identified inhibitors have the potential for further investigation.Communicated by Ramaswamy H. Sarma.

Targeting FAK/PYK2 with SJP1602 for Anti-Tumor Activity in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) presents significant challenges due to its aggressive nature and limited treatment options. Focal adhesion kinase (FAK) has emerged as a critical factor promoting tumor growth and metastasis in TNBC. Despite encouraging results from preclinical and early clinical trials with various FAK inhibitors, none have yet achieved clinical success in TNBC treatment. This study investigates the therapeutic potential of a novel dual inhibitor of FAK and PYK2, named SJP1602, for TNBC. In vitro experiments demonstrate that SJP1602 effectively inhibits FAK and PYK2 activities, showing potent effects on both kinases. SJP1602 shows concentration-dependent inhibition of cell growth, migration, invasion, and 3D spheroid formation in TNBC cell lines, surpassing the efficacy of other FAK inhibitors. Pharmacokinetic studies in rats indicate favorable bioavailability and sustained plasma concentrations of SJP1602, supporting its potential as a therapeutic agent. Furthermore, in TNBC xenograft models, SJP1602 exhibits significant dose-dependent inhibition of tumor growth. These promising results emphasize the potential of SJP1602 as a potent dual inhibitor of FAK and PYK2, deserving further investigation in clinical trials for TNBC treatment.

Anionic Phospholipids Stimulate the Proton Pumping Activity of the Plant Plasma Membrane P-Type H+-ATPase

The activity of membrane proteins depends strongly on the surrounding lipid environment. Here, we characterize the lipid stimulation of the plant plasma membrane H+-ATPase Arabidopsis thaliana H+-ATPase isoform 2 (AHA2) upon purification and reconstitution into liposomes of defined lipid compositions. We show that the proton pumping activity of AHA2 is stimulated by anionic phospholipids, especially by phosphatidylserine. This activation was independent of the cytoplasmic C-terminal regulatory domain of the pump. Molecular dynamics simulations revealed several preferential contact sites for anionic phospholipids in the transmembrane domain of AHA2. These contact sites are partially conserved in functionally different P-type ATPases from different organisms, suggesting a general regulation mechanism by the membrane lipid environment. Our findings highlight the fact that anionic lipids play an important role in the control of H+-ATPase activity.

Computational assessment of lipid facilitated membrane permeation of vancomycin using force-probe molecular dynamic simulation

In this study the efficacy of different edible lipids for drug permeation enhancement of vancomycin through biological membrane was investigated using molecular dynamic simulation. In this regard, at first the ability of the lipids for complex formation with the drug was evaluated for number of most common edible lipids including tripalmitin (TPA), trimyristin (TMY), labrafil (LAB), glycerol monostearate (GMS), glycerol monooleate (GMO), Distearoylphosphorylethanolamine (DSPE), dipalmitoylphosphatidylethanolamine (DPPE), Dipalmitoylphosphatidylcholine (DPPC), cholesterol (CL), stearic acid (SA), palmitic acid (PA) and oleic acid (OA). Then the complexes were pulled thorough a bilayer membrane while the changes in force were probed. The results showed that besides the SA, PA and OA the other examined lipids were able to perform a perfect molecular complex with the drug. Also the results of pulling simulation revealed that the least of force was needed for drug transmittance through the membrane when it was covered by LAB, TMY and DSPE. These results indicated that these lipids can be the excellent materials of choice as permeation enhancer for preparing a proper oral formulation of vancomycin.Communicated by Ramaswamy H. Sarma.

In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice

The aim of this study is to evaluate the anticonvulsant potential of schisandrin B, a main ingredient of Schisandra chinensis extracts. Schisandrin B showed anticonvulsant activity in the zebrafish larva pentylenetetrazole acute seizure assay but did not alter seizure thresholds in the intravenous pentylenetetrazole test in mice. Schisandrin B crosses the blood-brain barrier, which we confirmed in our in silico and in vivo analyses; however, the low level of its unbound fraction in the mouse brain tissue may explain the observed lack of anticonvulsant activity. Molecular docking revealed that the anticonvulsant activity of the compound in larval zebrafish might have been due to its binding to a benzodiazepine site within the GABAA receptor and/or the inhibition of the glutamate NMDA receptor. Although schisandrin B showed a beneficial anticonvulsant effect, toxicological studies revealed that it caused serious developmental impairment in zebrafish larvae, underscoring its teratogenic properties. Further detailed studies are needed to precisely identify the properties, pharmacological effects, and safety of schisandrin B.

The Potential Roles of Ficus carica Extract in the Management of COVID-19 Viral Infections: A Computer-Aided Drug Design Study

INTRODUCTION

The conventional processes of drug discovery are too expensive, time-consuming and the success rate is limited. Searching for alternatives that have evident safety and potential efficacy could save money, time and improve the current therapeutic regimen outcomes.

METHOD

Clinical phytotherapy implies the use of extracts of natural origin for prophylaxis, treatment, or management of human disorders. In this work, the potential role of common Fig (Ficus carica) in the management of COVID-19 infections has been explored. The antiviral effects of Cyanidin-3-rhamnoglucoside which is abundant in common Figs have been illustrated on COVID-19 targets. The immunomodulatory effect and the ability to ameliorate the cytokine storm associated with coronavirus infections have also been highlighted. This work involves various computational studies to investigate the potential roles of common figs in the management of COVID-19 viral infections.

RESULTS

Two molecular docking studies of all active ingredients in common Figs were conducted starting with MOE to provide initial insights, followed by Autodock Vina for further confirmation of the results of the top five compounds with the best docking score.

CONCLUSION

Finally, Molecular dynamic simulation alongside MMPBSA calculations were conducted using GROMACS to endorse and validate the entire work.

Deciphering Molecular Aspects of Potential α-Glucosidase Inhibitors within Aspergillus terreus: A Computational Odyssey of Molecular Docking-Coupled Dynamics Simulations and Pharmacokinetic Profiling

Hyperglycemia, as a hallmark of the metabolic malady diabetes mellitus, has been an overwhelming healthcare burden owing to its high rates of comorbidity and mortality, as well as prospective complications affecting different body organs. Available therapeutic agents, with α-glucosidase inhibitors as one of their cornerstone arsenal, control stages of broad glycemia while showing definitive characteristics related to their low clinical efficiency and off-target complications. This has propelled the academia and industrial section into discovering novel and safer candidates. Herein, we provided a thorough computational exploration of identifying candidates from the marine-derived Aspergillus terreus isolates. Combined structural- and ligand-based approaches using a chemical library of 275 metabolites were adopted for pinpointing promising α-glucosidase inhibitors, as well as providing guiding insights for further lead optimization and development. Structure-based virtual screening through escalating precision molecular docking protocol at the α-glucosidase canonical pocket identified 11 promising top-docked hits, with several being superior to the market drug reference, acarbose. Comprehensive ligand-based investigations of these hits' pharmacokinetics ADME profiles, physiochemical characterizations, and obedience to the gold standard Lipinski's rule of five, as well as toxicity and mutagenicity profiling, proceeded. Under explicit conditions, a molecular dynamics simulation identified the top-stable metabolites: butyrolactone VI (SK-44), aspulvinone E (SK-55), butyrolactone I 4''''-sulfate (SK-72), and terrelumamide B (SK-173). They depicted the highest free binding energies and steadiest thermodynamic behavior. Moreover, great structural insights have been revealed, including the advent of an aromatic scaffold-based interaction for ligand-target complex stability. The significance of introducing balanced hydrophobic/polar moieties, like triazole and other bioisosteres of carboxylic acid, has been highlighted across docking, ADME/Tox profiling, and molecular dynamics studies for maximizing binding interactions while assuring safety and optimal pharmacokinetics for targeting the intestinal-localized α-glucosidase enzyme. Overall, this study provided valuable starting points for developing new α-glucosidase inhibitors based on nature-derived unique scaffolds, as well as guidance for prospective lead optimization and development within future pre-clinical and clinical investigations.

Design, Synthesis, and Biological Evaluation of Novel Indanone Derivatives as Cholinesterase Inhibitors for Potential Use in Alzheimer's Disease

Indanone derivatives containing meta/para-substituted aminopropoxy benzyl/benzylidene moieties were designed based on the structures of donepezil and ebselen analogs as the cholinesterase inhibitors. The designed compounds were synthesized and their acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) inhibitory activities were measured. Inhibitory potencies (IC50 values) for the synthesized compounds ranged from 0.12 to 11.92 μM and 0.04 to 24.36 μM against AChE and BChE, respectively. Compound 5 c showed the highest AChE inhibitory potency with IC50 value of 0.12 μM, whereas the highest BChE inhibition was achieved by structure 7 b (IC50 =0.04 μM). Structure-activity relationship (SAR) analysis revealed that there is no significant difference between meta and para-substituted derivatives in AChE and BChE inhibition. However, the most potent AChE inhibitor 5 c belongs to meta-substituted compounds, while the most active BChE inhibitor is para-substituted derivative 7 b. The order of enzyme inhibition potency based on the substituted amine group is dimethyl amine>piperidine>morpholine. Compounds containing C=C linkage are more potent AChE inhibitors than the corresponding saturated structures. Molecular docking studies indicated that 5 c interacts with AChE in a very similar way to that observed experimentally for donepezil. The introduced indanone-aminopropoxy benzylidenes could be used in drug-discovery against Alzheimer's disease.

Toward α-1,3/4 fucosyltransferases targeted drug discovery: In silico uncovering of promising natural inhibitors of fucosyltransferase 6

Sialyl Lewis X (sLex ) antigen is a fucosylated cell-surface glycan that is normally involved in cell-cell interactions. The enhanced expression of sLex on cell surface glycans, which is attributed to the upregulation of fucosyltransferase 6 (FUT6), has been implicated in facilitating metastasis in human colorectal, lung, prostate, and oral cancers. The role that the upregulated FUT6 plays in the progression of tumor to malignancy, with reduced survival rates, makes it a potential target for anticancer drugs. Unfortunately, the lack of experimental structures for FUT6 has hampered the design and development of its inhibitors. In this study, we used in silico techniques to identify potential FUT6 inhibitors. We first modeled the three-dimensional structure of human FUT6 using AlphaFold. Then, we screened the natural compound libraries from the COCONUT database to sort out potential natural products (NPs) with best affinity toward the FUT6 model. As a result of these simulations, we identified three NPs for which we predicted binding affinities and interaction patterns quite similar to those we calculated for two experimentally tested FUT6 inhibitors, that is, fucose mimetic-1 and a GDP-triazole derived compound. We also performed molecular dynamics (MD) simulations for the FUT6 complexes with identified NPs, to investigate their stability. Analysis of the MD simulations showed that the identified NPs establish stable contacts with FUT6 under dynamics conditions. On these grounds, the three screened compounds appear as promising natural alternatives to experimentally tested FUT6 synthetic inhibitors, with expected comparable binding affinity. This envisages good prospects for future experimental validation toward FUT6 inhibition.

Antiviral activity of myricetin glycosylated compounds isolated from Marcetia taxifolia against chikungunya virus

The chikungunya virus (CHIKV) has produced epidemic outbreaks of significant public health impact. The clinical symptoms of this disease are fever, polyarthralgia, and skin rash, generally self-limiting, although patients may develop a chronic disabling condition or suffer lethal complications. Unfortunately, there is no specific treatment or vaccine available. Thus, the search for effective therapies to control CHIKV infection is an urgent need. This study evaluated the antiviral activity of flavonoids isolated from Marcetia taxifolia by in vitro and in silico analysis. Cytotoxicity of compounds was determined by MTT assay and viral load was assessed in cell substrates supernatants by plaque-forming and RT-qPCR assays. Selected molecules were analyzed by molecular docking assays. Myricetin 3-rhamnoside (MR) and myricetin 3-(6-rhamnosylgalactoside) (MRG) were tested for antiviral assays and analyzed by the TCID50 method and RT-qPCR. MR exhibited dose-dependent antiviral activity, reducing viral titer at concentrations of 150-18.8 μg/mL by at least 1-log. Similarly, MRG showed a significant decrease in viral titer at concentrations of 37.5, 9.4, and 2.3 μg/mL. RT-qPCR analysis also displayed a substantial reduction of CHIKV RNA for both flavonoids. Furthermore, molecular docking of the selected flavonoids proposed the nsP3 macrodomain as a possible target of action. Our study reveals that MR and MRG could be considered promising anti-CHIKV therapeutic agents. Molecular modeling studies showed MR and MRG ligands with a high affinity for the N-terminal region of the nsP3 macrodomain, postulating them as a potential target of action for the CHIKV control.

Synthesis and Structure Determination of Substituted Thiazole Derivatives as EGFR/BRAFV600E Dual Inhibitors Endowed with Antiproliferative Activity

2,3,4-trisubstituted thiazoles 3a-i, having a methyl group in position four, were synthesized by the reaction of 1,4-disubstituted thiosemicarbazides with chloroacetone in ethyl acetate/Et₃N at room temperature or in ethanol under reflux. The structures of new compounds were determined using NMR spectroscopy, mass spectrometry, and elemental analyses. Moreover, the structure of compound 3a was unambiguously confirmed with X-ray analysis. The cell viability assay of 3a-i at 50 µM was greater than 87%, and none of the tested substances were cytotoxic. Compounds 3a-i demonstrated good antiproliferative activity, with GI₅₀ values ranging from 37 to 86 nM against the four tested human cancer cell lines, compared to the reference erlotinib, which had a GI₅₀ value of 33 nM. The most potent derivatives were found to be compounds 3a, 3c, 3d, and 3f, with GI₅₀ values ranging from 37 nM to 54 nM. The EGFR-TK and BRAF^V600E inhibitory assays' results matched the antiproliferative assay's results, with the most potent derivatives, as antiproliferative agents, also being the most potent EGFR and BRAF^V600E inhibitors. The docking computations were employed to investigate the docking modes and scores of compounds 3a, 3c, 3d, and 3f toward BRAF^V600E and EGFR. Docking computations demonstrated the good affinity of compound 3f against BRAF^V600E and EGFR, with values of -8.7 and -8.5 kcal/mol, respectively.

Computational Prediction of the Interaction of Ivermectin with Fibrinogen

Hypercoagulability and formation of extensive and difficult-to-lyse microclots are a hallmark of both acute COVID-19 and long COVID. Fibrinogen, when converted to fibrin, is responsible for clot formation, but abnormal structural and mechanical clot properties can lead to pathologic thrombosis. Recent experimental evidence suggests that the spike protein (SP) from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) may directly bind to the blood coagulation factor fibrinogen and induce structurally abnormal blood clots with heightened proinflammatory activity. Accordingly, in this study, we used molecular docking and molecular dynamics simulations to explore the potential activity of the antiparasitic drug ivermectin (IVM) to prevent the binding of the SARS-CoV-2 SP to fibrinogen and reduce the occurrence of microclots. Our computational results indicate that IVM may bind with high affinity to multiple sites on the fibrinogen peptide, with binding more likely in the central, E region, and in the coiled-coil region, as opposed to the globular D region. Taken together, our in silico results suggest that IVM may interfere with SP-fibrinogen binding and, potentially, decrease the formation of fibrin clots resistant to degradation. Additional in vitro studies are warranted to validate whether IVM binding to fibrinogen is sufficiently stable to prevent interaction with the SP, and potentially reduce its thrombo-inflammatory effect in vivo.

Computational investigation of imidazo[2,1-b]oxazole derivatives as potential mutant BRAF kinase inhibitors: 3D-QSAR, molecular docking, molecular dynamics simulation, and ADMETox studies

BRAF inhibitors are known to be an effective therapeutic target for treating melanoma and other types of cancer. Using 3D-QSAR, molecular docking, and MD simulations, this study evaluated various imidazo[2,1-b]oxazole derivatives that function as mutant BRAF kinase inhibitors. Comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) were used to create the 3D-QSAR models. CoMSIA/SEHA model has solid predictive power across several models (Q² = 0.578; R² = 0.828; R²_pred = 0.74) and is the best model according to the numerous field models generated. The created model's predictive power was evaluated through external validation using a test set. CoMSIA/SEHA contour maps collect information that can be used to identify critical regions with solid anticancer activity. We developed four inhibitors with high predicted activity due to these observations. ADMET prediction was used to assess the toxicity of the proposed imidazo[2,1-b]oxazole compounds. The predictive molecules (T₁-T₄) demonstrated good ADMET properties, excluding the toxic active compounds 11r from the database. Molecular docking was also used to determine the patterns and modes of interactions between imidazo[2,1-b]oxazole ligands and receptors, which revealed that the proposed imidazo[2,1-b]oxazole scaffold was stable in the receptor's active site (PDB code: 4G9C). The suggested compounds (T1-T4) were subjected to molecular dynamics simulations lasting 100 ns to determine their binding free energies. The results showed that T₂ had a more favorable binding free energy (-149.552 kJ/mol) than T₁ (-112.556 kJ/mol), T₃ (-115.503 kJ/mol), and T₄ (-102.553 kJ/mol). The results suggest that the imidazo[2,1-b]oxazole compounds investigated in this study have potential as inhibitors of BRAF kinase and could be further developed as anticancer drugs. Highlights22 imidazo[2,1-b]oxazole compounds were subjected to research on three-dimensional quantitative conformational relationships.Using contour maps from 3D-QSAR models as a guide was used to figure out the areas and strategies for structural optimization.Combined molecular docking, molecular dynamics simulations, and binding free energy calculations to verify the inhibitor activity of the proposed 22 imidazo[2,1-b]oxazole compounds.Four potential B-RAF Kinase inhibitors were discovered, providing theoretical clues for developing a highly anticancer agent.Communicated by Ramaswamy H. Sarma.

Effect of Linker Elongation on the VGSC Affinity and Anticonvulsant Activity among 4-Alkyl-5-aryl-1,2,4-triazole-3-thione Derivatives

The main aim of the current project was to investigate the effect of the linker size in 4-alkyl-5-aryl-1,2,4-triazole-3-thione derivatives, known as a group of antiepileptic drug candidates, on their affinity towards voltage-gated sodium channels (VGSCs). The rationale of the study was based both on the SAR observations and docking simulations of the interactions between the designed ligands and the binding site of human VGSC. HYDE docking scores, which describe hydrogen bonding, desolvation, and hydrophobic effects, obtained for 5-[(3-chlorophenyl)ethyl]-4-butyl/hexyl-1,2,4-triazole-3-thiones, justified their beneficial sodium channel blocking activity. The results of docking simulations were verified using a radioligand binding assay with [³H]batrachotoxin. Unexpectedly, although the investigated triazole-based compounds acted as VGSC ligands, their affinities were lower than those of the respective analogs containing shorter alkyl linkers. Since numerous sodium channel blockers are recognized as antiepileptic agents, the obtained 1,2,4-triazole derivatives were examined for antiepileptic potential using an experimental model of tonic-clonic seizures in mice. Median effective doses (ED₅₀) of the compounds examined in MES test reached 96.6 ± 14.8 mg/kg, while their median toxic doses (TD₅₀), obtained in the rotarod test, were even as high as 710.5 ± 47.4 mg/kg.

Discovery of Novel Coumarin-Schiff Base Hybrids as Potential Acetylcholinesterase Inhibitors: Design, Synthesis, Enzyme Inhibition, and Computational Studies

To discover anti-acetylcholinesterase agents for the treatment of Alzheimer's disease (AD), a series of novel Schiff base-coumarin hybrids was rationally designed, synthesized successfully, and structurally characterized using Fourier transform infrared (FTIR), Nuclear magnetic resonance (NMR), and High-Resolution Mass Spectrometry (HRMS) analyses. These hybrids were evaluated for their potential inhibitory effect on acetylcholinesterase (AChE). All of them exhibited excellent inhibitory activity against AChE. The IC₅₀ values ranged from 87.84 to 515.59 μg/mL; hybrids 13c and 13d with IC₅₀ values of 0.232 ± 0.011 and 0.190 ± 0.004 µM, respectively, showed the most potent activity as acetylcholinesterase inhibitors (AChEIs). The reference drug, Galantamine, yielded an IC₅₀ of 1.142 ± 0.027 µM. Reactivity descriptors, including chemical potential (μ), chemical hardness (η), electrophilicity (ω), condensed Fukui function, and dual descriptors are calculated at wB97XD/6-311++ G (d,p) to identify reactivity changes of the designed compounds. An in-depth investigation of the natural charge pattern of the studied compounds led to a deep understanding of the important interaction centers between these compounds and the biological receptors of AChE. The molecular electrostatic surface potential (MESP) of the most active site in these derivatives was determined using high-quality information and visualization. Molecular docking analysis was performed to predict binding sites and binding energies. The structure-activity-property relationship studies indicated that the proposed compounds exhibit good oral bioavailability properties. To explore the stability and dynamic behavior of the ligand-receptor complexes, molecular dynamics simulations (MDS) were performed for 100 ns on the two best docked derivatives, 13c and 13d, with the AChE (4EY7) receptor. A popular method for determining the free binding energies (MM/GBSA) is performed using snapshots taken from the systems' trajectories at 100 ns. These results revealed that the complex system of compound 13d acquired a relatively more stable conformation and exhibited better descriptors than the complex system of compound 13c and the Galantamine drug, suggesting its potential as an effective inhibiting drug. The binding free energy analysis revealed that the 13d-4EY7 complex exhibited greater stability with AChE receptors compared to other complexes.

Reconstructing the transport cycle in the sugar porter superfamily using coevolution-powered machine learning

Sugar porters (SPs) represent the largest group of secondary-active transporters. Some members, such as the glucose transporters (GLUTs), are well known for their role in maintaining blood glucose homeostasis in mammals, with their expression upregulated in many types of cancers. Because only a few sugar porter structures have been determined, mechanistic models have been constructed by piecing together structural states of distantly related proteins. Current GLUT transport models are predominantly descriptive and oversimplified. Here, we have combined coevolution analysis and comparative modeling, to predict structures of the entire sugar porter superfamily in each state of the transport cycle. We have analyzed the state-specific contacts inferred from coevolving residue pairs and shown how this information can be used to rapidly generate free-energy landscapes consistent with experimental estimates, as illustrated here for the mammalian fructose transporter GLUT5. By comparing many different sugar porter models and scrutinizing their sequence, we have been able to define the molecular determinants of the transport cycle, which are conserved throughout the sugar porter superfamily. We have also been able to highlight differences leading to the emergence of proton-coupling, validating, and extending the previously proposed latch mechanism. Our computational approach is transferable to any transporter, and to other protein families in general.

Targeted intracellular delivery of antitubercular bioactive(s) to Mtb infected macrophages via transferrin functionalized nanoliposomes

The packaging of antimicrobials/chemotherapeutics into nanoliposomes can enhance their activity while minimizing toxicity. However, their use is still limited owing to inefficient/inadequate loading strategies. Several bioactive(s) which are non ionizable, and poorly aqueous soluble cannot be easily encapsulated into aqueous core of liposomes by using conventional means. Such bioactive(s) however could be encapsulated in the liposomes by forming their water soluble molecular inclusion complex with cyclodextrins. In this study, we developed Rifampicin (RIF) - 2-hydroxylpropyl-β-cyclodextrin (HP-β-CD) molecular inclusion complex. The HP-β-CD-RIF complex interaction was assessed by using computational analysis (molecular modeling). The HP-β-CD-RIF complex and Isoniazid were co-loaded in the small unilamellar vesicles (SUVs). Further, the developed system was functionalized with transferrin, a targeting moiety. Transferrin functionalized SUVs (Tf-SUVs) could preferentially deliver their payload intracellularly in the endosomal compartment of macrophages. In in vitro study on infected Raw 264.7 macrophage cells revealed that the encapsulated bioactive(s) could eradicate the pathogen more efficiently than free bioactive(s). In vivo studies further revealed that the Tf-SUVs could accumulate and maintain intracellular bioactive(s) concentrations in macrophages. The study suggests Tf-SUVs as a promising module for targeted delivery of a drug combination with improved/optimal therapeutic index and effective clinical outcomes.

pyCHARMM: Embedding CHARMM Functionality in a Python Framework

CHARMM is rich in methodology and functionality as one of the first programs addressing problems of molecular dynamics and modeling of biological macromolecules and their partners, e.g., small molecule ligands. When combined with the highly developed CHARMM parameters for proteins, nucleic acids, small molecules, lipids, sugars, and other biologically relevant building blocks, and the versatile CHARMM scripting language, CHARMM has been a trendsetting platform for modeling studies of biological macromolecules. To further enhance the utility of accessing and using CHARMM functionality in increasingly complex workflows associated with modeling biological systems, we introduce pyCHARMM, Python bindings, functions, and modules to complement and extend the extensive set of modeling tools and methods already available in CHARMM. These include access to CHARMM function-generated variables associated with the system (psf), coordinates, velocities and forces, atom selection variables, and force field related parameters. The ability to augment CHARMM forces and energies with energy terms or methods derived from machine learning or other sources, written in Python, CUDA, or OpenCL and expressed as Python callable routines is introduced together with analogous functions callable during dynamics calculations. Integration of Python-based graphical engines for visualization of simulation models and results is also accessible. Loosely coupled parallelism is available for workflows such as free energy calculations, using MBAR/TI approaches or high-throughput multisite λ-dynamics (MSλD) free energy methods, string path optimization calculations, replica exchange, and molecular docking with a new Python-based CDOCKER module. CHARMM accelerated platform kernels through the CHARMM/OpenMM API, CHARMM/DOMDEC, and CHARMM/BLaDE API are also readily integrated into this Python framework. We anticipate that pyCHARMM will be a robust platform for the development of comprehensive and complex workflows utilizing Python and its extensive functionality as well as an optimal platform for users to learn molecular modeling methods and practices within a Python-friendly environment such as Jupyter Notebooks.

Novel Isosteviol-Based FXa Inhibitors: Molecular Modeling, In Silico Design and Docking Simulation

Direct oral anticoagulants are an important and relatively new class of synthetic anticoagulant drugs commonly used for the pharmacotherapy of thromboembolic disorders. However, they still have some limitations and serious side effects, which continuously encourage medicinal chemists to search for new active compounds acting as human-activated coagulation factor X (FXa) inhibitors. Isosteviol is a nontoxic hydrolysis product of naturally occurring stevioside and possesses a wide range of therapeutic properties, including anticoagulant activity. The present contribution describes the in silico design of novel oxime ether isosteviol derivatives as well as a molecular modeling approach based on QSAR analysis and a docking simulation for searching for novel isosteviol-based compounds as potential FXa inhibitors. The elaborated ANN model, encompassing topological and geometrical information, exhibited a significant correlation with FXa-inhibitory activity. Moreover, the docking simulation indicated six of the most promising isosteviol-like compounds for further investigation. Analysis showed that the most promising derivatives contain heterocyclic, aromatic, five-membered moieties, with substituents containing chlorine or fluorine atoms. It is anticipated that the findings reported in the present work may provide useful information for designing effective FXa inhibitors as anticoagulant agents.

Network pharmacology and multitarget analysis of Nigella sativa in the management of diabetes and obesity: a computational study

Obesity and diabetes are commonly associated with one another and represent a significant global health issue, with a recent surge in disease incidence. Nigella sativa, also known as black cumin, is believed to possess several health benefits, including anti-diabetic, anticancer, antioxidant, antimicrobial, and anti-obesity properties. In this study, we aimed to identify the active compounds derived from N. sativa, which can potentially inhibit key protein targets and signaling pathways associated with diabesity treatment. We employed an exhaustive in silico search, which led to the identification of 22 potential compounds. Out of these, only five hits were found to be non-toxic, including Arabic and ascorbic acids, dihydrocodeine, catechin, and kaempferol. Our analysis revealed that these hits were associated with genes such as AKT1, IL6, SRC, and EGFR. Finally, we conducted molecular docking and molecular dynamics simulations, which identified kaempferol as the best binder for AKT1 in comparison to the reference molecule. Overall, our in silico integrated pipeline provides a useful approach to identify non-toxic phytocompounds as promising drug candidates to treat diabetes and obesity.Communicated by Ramaswamy H. Sarma.

Emerging Trends of Computational Chemistry and Molecular Modeling in Froth Flotation: A Review

Froth flotation is the most versatile process in mineral beneficiation, extensively used to concentrate a wide range of minerals. This process comprises mixtures of more or less liberated minerals, water, air, and various chemical reagents, involving a series of intermingled multiphase physical and chemical phenomena in the aqueous environment. Today's main challenge facing the froth flotation process is to gain atomic-level insights into the properties of its inherent phenomena governing the process performance. While it is often challenging to determine these phenomena via trial-and-error experimentations, molecular modeling approaches not only elicit a deeper understanding of froth flotation but can also assist experimental studies in saving time and budget. Thanks to the rapid development of computer science and advances in high-performance computing (HPC) infrastructures, theoretical/computational chemistry has now matured enough to successfully and gainfully apply to tackle the challenges of complex systems. In mineral processing, however, advanced applications of computational chemistry are increasingly gaining ground and demonstrating merit in addressing these challenges. Accordingly, this contribution aims to encourage mineral scientists, especially those interested in rational reagent design, to become familiarized with the necessary concepts of molecular modeling and to apply similar strategies when studying and tailoring properties at the molecular level. This review also strives to deliver the state-of-the-art integration and application of molecular modeling in froth flotation studies to assist either active researchers in this field to disclose new directions for future research or newcomers to the field to initiate innovative works.

Identification of potential biogenic chalcones against antibiotic resistant efflux pump (AcrB) via computational study

The cases of bacterial multidrug resistance are increasing every year and becoming a serious concern for human health. Multidrug efflux pumps are key players in the formation of antibiotic resistance, which transfer out a broad spectrum of drugs from the cell and convey resistance to the host. Efflux pumps have significantly reduced the efficacy of the previously available antibiotic armory, thereby increasing the frequency of therapeutic failures. In gram-negative bacteria, the AcrAB-TolC efflux pump is the principal transporter of the substrate and plays a major role in the formation of antibiotic resistance. In the current work, advanced computer-aided drug discovery approaches were utilized to find hit molecules from the library of biogenic chalcones against the bacterial AcrB efflux pump. The results of the performed computational studies via molecular docking, drug-likeness prediction, pharmacokinetic profiling, pharmacophore mapping, density functional theory, and molecular dynamics simulation study provided ZINC000004695648, ZINC000014762506, ZINC000014762510, ZINC000095099506, and ZINC000085510993 as stable hit molecules against the AcrB efflux pumps. Identified hits could successfully act against AcrB efflux pumps after optimization as lead molecules.Communicated by Ramaswamy H. Sarma.

Solvent-Free Loading of Vitamin A Palmitate into β-Cyclodextrin Metal-Organic Frameworks for Stability Enhancement

Cyclodextrin metal-organic frameworks (CD-MOFs) exhibit a high structural diversity, which contributes to their functional properties. In this study, we have successfully synthesized a novel type of β-cyclodextrin metal-organic framework (β-CD-POF(I)) that exhibits excellent drug adsorption capacity and enhances stability. Single-crystal X-ray diffraction analysis revealed that β-CD-POF(I) possessed the dicyclodextrin channel moieties and long-parallel tubular cavities. Compared with the reported β-CD-MOFs, the β-CD-POF(I) has a more promising drug encapsulation capability. Here, the stability of vitamin A palmitate (VAP) was effectively improved by the solvent-free method. Molecular modeling and other characterization techniques like synchrotron radiation Fourier transform infrared spectroscopy (SR-FTIR), differential scanning calorimetry (DSC), powder X-ray diffraction (PXRD), thermogravimetric analysis (TGA), and nitrogen adsorption isotherm were applied to confirm that the VAP was successfully encapsulated into the channel formed by the dicyclodextrin pairs. Furthermore, the mechanism of stability enhancement for VAP was determined to be due to the constraint and separation effects of β-CD pairs on VAP. Therefore, β-CD-POF(I) is capable of trapping and stabilizing certain unstable drug molecules, offering benefits and application possibilities. One kind of cyclodextrin particle with characteristic shapes of dicyclodextrin channel moieties and parallel tubular cavities, which was synthesized by a facile process. Subsequently, the spatial structure and characteristics of the β-CD-POF(I) were primarily confirmed. The structure of β-CD-POF(I) was then compared to that of KOH-β-CD-MOF, and a better material for vitamin A palmitate (VAP) encapsulation was determined. VAP was successfully loaded into the particles by solvent-free method. The arrangement of spatial structure made cyclodextrin molecular cavity encapsulation in β-CD-POF(I) more stable for VAP capture than that of KOH-β-CD-MOF.

Harmonic Vibrational Frequency Simulation of Pharmaceutical Molecules via a Novel Multi-Molecular Fragment Interception Method

By means of a computational method based on Density Functional Theory (DFT), using commercially available software, a novel method for simulating equilibrium geometry harmonic vibrational frequencies is proposed. Finasteride, Lamivudine, and Repaglinide were selected as model molecules to study the adaptability of the new method. Three molecular models, namely the single-molecular, central-molecular, and multi-molecular fragment models, were constructed and calculated by Generalized Gradient Approximations (GGAs) with the PBE functional via the Material Studio 8.0 program. Theoretical vibrational frequencies were assigned and compared to the corresponding experimental data. The results indicated that the traditional single-molecular calculation and scaled spectra with scale factor exhibited the worst similarity for all three pharmaceutical molecules among the three models. Furthermore, the central-molecular model with a configuration closer to the empirical structure resulted in a reduction of mean absolute error (MAE) and root mean squared error (RMSE) in all three pharmaceutics, including the hydrogen-bonded functional groups. However, the improvement in computational accuracy for different drug molecules using the central-molecular model for vibrational frequency calculation was unstable. Whereas, the new multi-molecular fragment interception method showed the best agreement with experimental results, exhibiting MAE and RMSE values of 8.21 cm^-1 and 18.35 cm^-1 for Finasteride, 15.95 cm^-1 and 26.46 cm^-1 for Lamivudine, and 12.10 cm^-1 and 25.82 cm^-1 for Repaglinide. Additionally, this work provides comprehensive vibrational frequency calculations and assignments for Finasteride, Lamivudine, and Repaglinide, which have never been thoroughly investigated in previous research.

In vivo and in silico comparison analyses of Cry toxin activities toward the sugarcane giant borer

The sugarcane giant borer, Telchin licus licus, is an insect pest that causes significant losses in sugarcane crops and in the sugar-alcohol sector. Chemical and manual control methods are not effective. As an alternative, in the current study, we have screened Bacillus thuringiensis (Bt) Cry toxins with high toxicity against this insect. Bioassays were conducted to determine the activity of four Cry toxins (Cry1A (a, b, and c) and Cry2Aa) against neonate T. licus licus larvae. Notably, the Cry1A family toxins had the lowest LC₅₀ values, in which Cry1Ac presented 2.1-fold higher activity than Cry1Aa, 1.7-fold larger than Cry1Ab, and 9.7-fold larger than Cry2Aa toxins. In silico analyses were performed as a perspective to understand putative interactions between T. licus licus receptors and Cry1A toxins. The molecular dynamics and docking analyses for three putative aminopeptidase N (APN) receptors (TlAPN1, TlAPN3, and TlAPN4) revealed evidence for the amino acids that may be involved in the toxin-receptor interactions. Notably, the properties of Cry1Ac point to an interaction site that increases the toxin's affinity for the receptor and likely potentiate toxicity. The interacting amino acid residues predicted for Cry1Ac in this work are probably those shared by the other Cry1A toxins for the same region of APNs. Thus, the presented data extend the existing knowledge of the effects of Cry toxins on T. licus licus and should be considered in further development of transgenic sugarcane plants resistant to this major occurring insect pest in sugarcane fields.

In silico design and computational evaluation of novel 2-arylaminopyrimidine-based compounds as potential multi-targeted protein kinase inhibitors: application for the native and mutant (T315I) Bcr-Abl tyrosine kinase

An integrated computational approach to drug discovery was used to identify novel potential inhibitors of the native and mutant (T315I) Bcr-Abl tyrosine kinase, the enzyme playing a key role in the pathogenesis of chronic myeloid leukemia (CML). This approach included i) design of chimeric molecules based on the 2-arylaminopyrimidine fragment, the main pharmacophore of the Abl kinase inhibitors imatinib and nilotinib used in the clinic for the CML treatment, ii) molecular docking of these compounds with the ATP-binding site of the native and mutant Abl kinase, iii) refinement of the ligand-binding poses by the quantum chemical method PM7, iv) molecular dynamics simulations of the ligand/Abl complexes, and v) prediction of the ligand/Abl binding affinity in terms of scoring functions of molecular docking, machine learning, quantum chemistry, and molecular dynamics. As a result, five top-ranking compounds able to effectively block the enzyme catalytic site were identified. According to the data obtained, these compounds exhibit close modes of binding to the Abl kinase active site that are mainly provided by hydrogen bonds and multiple van der Waals contacts. The identified compounds show high binding affinity to the native and mutant Abl kinase comparable with the one calculated for the FDA-approved kinase-targeted inhibitors imatinib, nilotinib, and ponatinib used in the calculations as a positive control. The results obtained testify to the predicted drug candidates against CML may serve as good scaffolds for the design of novel anticancer agents able to target the ATP-binding pocket of the native and mutant Abl kinase.Communicated by Ramaswamy H. Sarma.

Cesium activates the neurotransmitter receptor for glycine

The monovalent cations sodium and potassium are crucial for the proper functioning of excitable cells, but, in addition, other monovalent alkali metal ions such as cesium and lithium can also affect neuronal physiology. For instance, there have been recent reports of adverse effects resulting from self-administered high concentrations of cesium in disease conditions, prompting the Food and Drug Administration (FDA) to issue an alert concerning cesium chloride. As we recently found that the monovalent cation NH4+ activates glycine receptors (GlyRs), we investigated the effects of alkali metal ions on the function of the GlyR, which belongs to one of the most widely distributed neurotransmitter receptors in the peripheral and central nervous systems. Whole-cell voltage clamp electrophysiology was performed with HEK293T cells transiently expressing different splice and RNA-edited variants of GlyR α2 and α3 homopentameric channels. By examining the influence of various milli- and sub-millimolar concentrations of lithium, sodium, potassium, and cesium on these GlyRs in comparison to its natural ligand glycine (0.1 mM), we could show that cesium activates GlyRs in a concentration- and post-transcriptional-dependent way. Additionally, we conducted atomistic molecular dynamic simulations on GlyR α3 embedded in a membrane bilayer with potassium and cesium, respectively. The simulations revealed slightly different GlyR-ion binding profiles for potassium and cesium, identifying interactions near the glycine binding pocket (potassium and cesium) and close to the RNA-edited site (cesium) in the extracellular GlyR domain. Together, these findings show that cesium acts as an agonist of GlyRs.

Mitigating the risk of antimalarial resistance via covalent dual-subunit inhibition of the Plasmodium proteasome

The Plasmodium falciparum proteasome constitutes a promising antimalarial target, with multiple chemotypes potently and selectively inhibiting parasite proliferation and synergizing with the first-line artemisinin drugs, including against artemisinin-resistant parasites. We compared resistance profiles of vinyl sulfone, epoxyketone, macrocyclic peptide, and asparagine ethylenediamine inhibitors and report that the vinyl sulfones were potent even against mutant parasites resistant to other proteasome inhibitors and did not readily select for resistance, particularly WLL that displays covalent and irreversible binding to the catalytic β2 and β5 proteasome subunits. We also observed instances of collateral hypersensitivity, whereby resistance to one inhibitor could sensitize parasites to distinct chemotypes. Proteasome selectivity was confirmed using CRISPR/Cas9-edited mutant and conditional knockdown parasites. Molecular modeling of proteasome mutations suggested spatial contraction of the β5 P1 binding pocket, compromising compound binding. Dual targeting of P. falciparum proteasome subunits using covalent inhibitors provides a potential strategy for restoring artemisinin activity and combating the spread of drug-resistant malaria.

Conformationally Selective 2-Aminotetralin Ligands Targeting the alpha2A- and alpha2C-Adrenergic Receptors

Many important physiological processes are mediated by alpha2A- and alpha2C-adrenergic receptors (α2Rs), a subtype of class A G protein-coupled receptors (GPCRs). However, α2R signaling is poorly understood, and there are few approved medications targeting these receptors. Drug discovery aimed at α2Rs is complicated by the high degree of binding pocket homology between α2AR and α2CR, which confounds ligand-mediated selective activation or inactivation of signaling associated with a particular subtype. Meanwhile, α2R signaling is complex and it is reported that activating α2AR is beneficial in many clinical contexts, while activating α2CR signaling may be detrimental to these positive effects. Here, we report on a novel 5-substituted-2-aminotetralin (5-SAT) chemotype that, depending on substitution, has diverse pharmacological activities at α2Rs. Certain lead 5-SAT analogues act as partial agonists at α2ARs, while functioning as inverse agonists at α2CRs, a novel pharmacological profile. Leads demonstrate high potency (e.g., EC50 < 2 nM) at the α2AR and α2CRs regarding Gαi-mediated inhibition of adenylyl cyclase and production of cyclic adenosine monophosphate (cAMP). To help understand the molecular basis of 5-SAT α2R multifaceted functional activity, α2AR and α2CR molecular models were built from the crystal structures and 1 μs molecular dynamics (MD) simulations and molecular docking experiments were performed for a lead 5-SAT with α2AR agonist and α2CR inverse agonist activity, i.e., (2S)-5-(2'-fluorophenyl)-N,N-dimethyl-1,2,3,4-tetrahydronaphthalen-2-amine (FPT), in comparison to the FDA-approved (for opioid withdrawal symptoms) α2AR/α2CR agonist lofexidine. Results reveal several interactions between FPT and α2AR and α2CR amino acids that may impact the functional activity. The computational data in conjunction with experimental in vitro affinity and function results provide information to understand ligand stabilization of functionally distinct GPCR conformations regarding α2AR and α2CRs.

Design and Identification of Inhibitors for the Spike-ACE2 Target of SARS-CoV-2

When an epidemic started in the Chinese city of Wuhan in December 2019, coronavirus was identified as the cause. Infection by the virus occurs through the interaction of viral S protein with the hosts' angiotensin-converting enzyme 2. By leveraging resources such as the DrugBank database and bioinformatics techniques, ligands with potential activity against the SARS-CoV-2 spike protein were designed and identified in this investigation. The FTMap server and the Molegro software were used to determine the active site of the Spike-ACE2 protein's crystal structure. Virtual screening was performed using a pharmacophore model obtained from antiparasitic drugs, obtaining 2000 molecules from molport^®. The ADME/Tox profiles were used to identify the most promising compounds with desirable drug characteristics. The binding affinity investigation was then conducted with selected candidates. A molecular docking study showed five structures with better binding affinity than hydroxychloroquine. Ligand_003 showed a binding affinity of -8.645 kcal·mol^-1, which was considered an optimal value for the study. The values presented by ligand_033, ligand_013, ligand_044, and ligand_080 meet the profile of novel drugs. To choose compounds with favorable potential for synthesis, synthetic accessibility studies and similarity analyses were carried out. Molecular dynamics and theoretical IC₅₀ values (ranging from 0.459 to 2.371 µM) demonstrate that these candidates are promising for further tests. Chemical descriptors showed that the candidates had strong molecule stability. Theoretical analyses here show that these molecules have potential as SARS-CoV-2 antivirals and therefore warrant further investigation.

Blood pH Analysis in Combination with Molecular Medical Tools in Relation to COVID-19 Symptoms

The global outbreak of SARS-CoV-2/COVID-19 provided the stage to accumulate an enormous biomedical data set and an opportunity as well as a challenge to test new concepts and strategies to combat the pandemic. New research and molecular medical protocols may be deployed in different scientific fields, e.g., glycobiology, nanopharmacology, or nanomedicine. We correlated clinical biomedical data derived from patients in intensive care units with structural biology and biophysical data from NMR and/or CAMM (computer-aided molecular modeling). Consequently, new diagnostic and therapeutic approaches against SARS-CoV-2 were evaluated. Specifically, we tested the suitability of incretin mimetics with one or two pH-sensitive amino acid residues as potential drugs to prevent or cure long-COVID symptoms. Blood pH values in correlation with temperature alterations in patient bodies were of clinical importance. The effects of biophysical parameters such as temperature and pH value variation in relation to physical-chemical membrane properties (e.g., glycosylation state, affinity of certain amino acid sequences to sialic acids as well as other carbohydrate residues and lipid structures) provided helpful hints in identifying a potential Achilles heel against long COVID. In silico CAMM methods and in vitro NMR experiments (including ³¹P NMR measurements) were applied to analyze the structural behavior of incretin mimetics and SARS-CoV fusion peptides interacting with dodecylphosphocholine (DPC) micelles. These supramolecular complexes were analyzed under physiological conditions by ¹H and ³¹P NMR techniques. We were able to observe characteristic interaction states of incretin mimetics, SARS-CoV fusion peptides and DPC membranes. Novel interaction profiles (indicated, e.g., by ³¹P NMR signal splitting) were detected. Furthermore, we evaluated GM1 gangliosides and sialic acid-coated silica nanoparticles in complex with DPC micelles in order to create a simple virus host cell membrane model. This is a first step in exploring the structure-function relationship between the SARS-CoV-2 spike protein and incretin mimetics with conserved pH-sensitive histidine residues in their carbohydrate recognition domains as found in galectins. The applied methods were effective in identifying peptide sequences as well as certain carbohydrate moieties with the potential to protect the blood-brain barrier (BBB). These clinically relevant observations on low blood pH values in fatal COVID-19 cases open routes for new therapeutic approaches, especially against long-COVID symptoms.

Synthesis, Computational Insights, and Evaluation of Novel Sigma Receptors Ligands

The development of diazabicyclo[4.3.0]nonane and 2,7-diazaspiro[3.5]nonane derivatives as sigma receptors (SRs) ligands is reported. The compounds were evaluated in S1R and S2R binding assays, and modeling studies were carried out to analyze the binding mode. The most notable compounds, 4b (AD186, K_iS1R = 2.7 nM, K_iS2R = 27 nM), 5b (AB21, K_iS1R = 13 nM, K_iS2R = 102 nM), and 8f (AB10, K_iS1R = 10 nM, K_iS2R = 165 nM), have been screened for analgesic effects in vivo, and their functional profile was determined through in vivo and in vitro models. Compounds 5b and 8f reached the maximum antiallodynic effect at 20 mg/kg. The selective S1R agonist PRE-084 completely reversed their action, indicating that the effects are entirely dependent on the S1R antagonism. Conversely, compound 4b sharing the 2,7-diazaspiro[3.5]nonane core as 5b was completely devoid of antiallodynic effect. Interestingly, compound 4b fully reversed the antiallodynic effect of BD-1063, indicating that 4b induces an S1R agonistic in vivo effect. The functional profiles were confirmed by the phenytoin assay. Our study might establish the importance of 2,7-diazaspiro[3.5]nonane core for the development of S1R compounds with specific agonist or antagonist profile and the role of the diazabicyclo[4.3.0]nonane in the development of novel SR ligands.