Computational Modeling on Aquaporin-3 as Skin Cancer Target: A Virtual Screening Study

Methods for studying mammalian aquaporin biology

Aquaporins (AQPs), transmembrane water-conducting channels, have earned a great deal of scrutiny for their critical physiological roles in healthy and disease cell states, especially in the biomedical field. Numerous methods have been implemented to elucidate the involvement of AQP-mediated water transport and downstream signaling activation in eliciting whole cell, tissue, and organ functional responses. To modulate these responses, other methods have been employed to investigate AQP druggability. This review discusses standard in vitro, in vivo, and in silico methods for studying AQPs, especially for biomedical and mammalian cell biology applications. We also propose some new techniques and approaches for future AQP research to address current gaps in methodology.

Lipid-Based Nanoparticles for Drug/Gene Delivery: An Overview of the Production Techniques and Difficulties Encountered in Their Industrial Development

Over the past decade, the therapeutic potential of nanomaterials as novel drug delivery systems complementing conventional pharmacology has been widely acknowledged. Among these nanomaterials, lipid-based nanoparticles (LNPs) have shown remarkable pharmacological performance and promising therapeutic outcomes, thus gaining substantial interest in preclinical and clinical research. In this review, we introduce the main types of LNPs used in drug formulations such as liposomes, nanoemulsions, solid lipid nanoparticles, nanostructured lipid carriers, and lipid polymer hybrid nanoparticles, focusing on their main physicochemical properties and therapeutic potential. We discuss computational studies and modeling techniques to enhance the understanding of how LNPs interact with therapeutic cargo and to predict the potential effectiveness of such interactions in therapeutic applications. We also analyze the benefits and drawbacks of various LNP production techniques such as nanoprecipitation, emulsification, evaporation, thin film hydration, microfluidic-based methods, and an impingement jet mixer. Additionally, we discuss the major challenges associated with industrial development, including stability and sterilization, storage, regulatory compliance, reproducibility, and quality control. Overcoming these challenges and facilitating regulatory compliance represent the key steps toward LNP's successful commercialization and translation into clinical settings.

Aquaporin 3 inhibition suppresses the mitochondrial respiration rate and viability of multiple myeloma cells

Aquaporin 3 (AQP3) is a member of the aquaporin water channel family expressed by numerous cell types, including some cancer cells. Accumulating evidence suggests that AQP3 inhibition may impede cancer progression, but drugs targeting AQP3 are still in the early pre-clinical stage of development. Here, we examined the effect of AQP3 inhibition on multiple myeloma (MM), an incurable plasma cell malignancy. Four MM cell lines were cultured in the presence of an anti-AQP3 monoclonal antibody (mAb), the AQP3 inhibitor DFP00173, or corresponding controls, and the effects on cell viability, proliferation, apoptosis, and mitochondrial respiration capacity were compared. Both anti-AQP3 mAb and DFP00173 reduced cell growth, mitochondrial respiration rate, and electron transport chain complex I activity. Both agents also potentiated the antiproliferative efficacy of the anticancer drug venetoclax. Administration of the anti-AQP3 mAb to immunodeficient mice inoculated with RPMI8226 or KMS-11 MM cells significantly suppressed tumor growth. These data provide evidence that AQP3 blockade can suppress MM cell growth in vitro and tumor growth in mice. Thus, AQP3 inhibition may be an effective therapeutic strategy for MM.

Desmodin and isopongachromene as potential inhibitors of cyclin-dependent kinase 5: phytoconstituents targeting anticancer and neurological therapy

Cyclin-dependent kinase 5 (CDK5) is a proline-directed serine-threonine protein kinase vital for neuronal cell cycle arrest and differentiation. It activates by binding with p35 and p39 and is important for the functioning of the nervous system. A growing body of evidence suggests that CDK5 contributes to the onset and progression of neurodegeneration and tumorigenesis and represents itself as a potential therapeutic target. Our research illustrates virtual screening of phytochemicals from the IMPPAT (Indian Medicinal Plants, Phytochemistry and Therapeutics) library to search for potential inhibitors of CDK5. Initially, the compounds from the parent library were filtered out via their physicochemical properties following the Lipinski rule of five. Then sequentially, molecular docking-based virtual screening, PAINS filter, ADMET, PASS analysis, and molecular dynamics (MD) simulation were done using various computational tools to rule out adversities that can cause hindrances in the identification of potential inhibitors of CDK5. Finally, two compounds were selected via the extensive screening showing significant binding with CDK5 ATP-binding pocket and ultimately were selected as potent ATP-competitive inhibitors of CDK5. Finally, we propose that the elucidated compounds Desmodin and Isopongachromene can be used further in the drug discovery process and act as therapeutics in the medical industry to treat certain complex diseases, including cancer and neurodegeneration.Communicated by Ramaswamy H. Sarma.

Identification of Potential Druggable Targets and Structure-Based Virtual Screening for Drug-like Molecules against the Shrimp Pathogen Enterocytozoon hepatopenaei

Enterocytozoon hepatopenaei (EHP) causes slow growth syndrome in shrimp, resulting in huge economic losses for the global shrimp industry. Despite worldwide reports, there are no effective therapeutics for controlling EHP infections. In this study, five potential druggable targets of EHP, namely, aquaporin (AQP), cytidine triphosphate (CTP) synthase, thymidine kinase (TK), methionine aminopeptidase2 (MetAP2), and dihydrofolate reductase (DHFR), were identified via functional classification of the whole EHP proteome. The three-dimensional structures of the proteins were constructed using the artificial-intelligence-based program AlphaFold 2. Following the prediction of druggable sites, the ZINC15 and ChEMBL databases were screened against targets using docking-based virtual screening. Molecules with affinity scores ≥ 7.5 and numbers of interactions ≥ 9 were initially selected and subsequently enriched based on their ADMET properties and electrostatic complementarities. Five compounds were finally selected against each target based on their complex stabilities and binding energies. The compounds CHEMBL3703838, CHEMBL2132563, and CHEMBL133039 were selected against AQP; CHEMBL1091856, CHEMBL1162979, and CHEMBL525202 against CTP synthase; CHEMBL4078273, CHEMBL1683320, and CHEMBL3674540 against TK; CHEMBL340488, CHEMBL1966988, and ZINC000828645375 against DHFR; and CHEMBL3913373, ZINC000016682972, and CHEMBL3142997 against MetAP2.The compounds exhibited high stabilities and low binding free energies, indicating their abilities to suppress EHP infections; however, further validation is necessary for determining their efficacy.

Melanoma 2.0. Skin cancer as a paradigm for emerging diagnostic technologies, computational modelling and artificial intelligence

We live in an unprecedented time in oncology. We have accumulated samples and cases in cohorts larger and more complex than ever before. New technologies are available for quantifying solid or liquid samples at the molecular level. At the same time, we are now equipped with the computational power necessary to handle this enormous amount of quantitative data. Computational models are widely used helping us to substantiate and interpret data. Under the label of systems and precision medicine, we are putting all these developments together to improve and personalize the therapy of cancer. In this review, we use melanoma as a paradigm to present the successful application of these technologies but also to discuss possible future developments in patient care linked to them. Melanoma is a paradigmatic case for disruptive improvements in therapies, with a considerable number of metastatic melanoma patients benefiting from novel therapies. Nevertheless, a large proportion of patients does not respond to therapy or suffers from adverse events. Melanoma is an ideal case study to deploy advanced technologies not only due to the medical need but also to some intrinsic features of melanoma as a disease and the skin as an organ. From the perspective of data acquisition, the skin is the ideal organ due to its accessibility and suitability for many kinds of advanced imaging techniques. We put special emphasis on the necessity of computational strategies to integrate multiple sources of quantitative data describing the tumour at different scales and levels.

Inhibiting Cyclin-Dependent Kinase 6 by Taurine: Implications in Anticancer Therapeutics

Cyclin-dependent kinase 6 (CDK6) is linked with a cyclin partner and plays a crucial role in the early stages of cancer development. It is currently a potential drug target for developing therapeutic molecules targeting cancer therapy. Here, we have identified taurine as an inhibitor of CDK6 using combined in silico and experimental studies. We performed various experiments to find the binding affinity of taurine with CDK6. Molecular docking analysis revealed critical residues of CDK6 that are involved in taurine binding. Fluorescence measurement studies showed that taurine binds to CDK6 with a significant binding affinity, with a binding constant of K = 0.7 × 10⁷ M^-1 for the CDK6-taurine complex. Enzyme inhibition assay suggested taurine as a good inhibitor of CDK6 possessing an IC₅₀ value of 4.44 μM. Isothermal titration calorimetry analysis further confirmed a spontaneous binding of taurine with CDK6 and delineated the thermodynamic parameters for the CDK6-taurine system. Altogether, this study established taurine as a CDK6 inhibitor, providing a base for using taurine and its derivatives in CDK6-associated cancer and other diseases.

Discovering Tuberosin and Villosol as Potent and Selective Inhibitors of AKT1 for Therapeutic Targeting of Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is a major cause of death in developing countries because of high tobacco consumption. RAC-alpha serine-threonine kinase (AKT1) is considered as an attractive drug target because its prolonged activation and overexpression are associated with cancer progression and metastasis. In addition, several AKT1 inhibitors are being developed to control OSCC and other associated forms of cancers. We performed a screening of the IMPPAT (Indian Medicinal Plants, Phytochemistry and Therapeutics) database to discover promising AKT1 inhibitors which pass through various important filters such as ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties, physicochemical properties, PAINS (pan-assay interference compounds) filters, PASS (prediction of activity spectra for substances) analysis, and specific interactions with AKT1. Molecules bearing admirable binding affinity and specificity towards AKT1 were selected for further analysis. Initially, we identified 30 natural compounds bearing appreciable affinity and specific interaction with AKT1. Finally, tuberosin and villosol were selected as potent and selective AKT1 inhibitors. To obtain deeper insights into binding mechanism and selectivity, we performed an all-atom molecular dynamics (MD) simulation and principal component analysis (PCA). We observed that both tuberosin and villosol strongly bind to AKT1, and their complexes were stable throughout the simulation trajectories. Our in-depth structure analysis suggested that tuberosin and villosol could be further exploited in the therapeutic targeting of OSCC and other cancers after further clinical validations.

Fast approximative methods for study of ligand transport and rational design of improved enzymes for biotechnologies

Acceleration of chemical reactions by the enzymes optimized using protein engineering represents one of the key pillars of the contribution of biotechnology towards sustainability. Tunnels and channels of enzymes with buried active sites enable the exchange of ligands, ions, and water molecules between the outer environment and active site pockets. The efficient exchange of ligands is a fundamental process of biocatalysis. Therefore, enzymes have evolved a wide range of mechanisms for repetitive conformational changes that enable periodic opening and closing. Protein-ligand interactions are traditionally studied by molecular docking, whereas molecular dynamics is the method of choice for studying conformational changes and ligand transport. However, computational demands make molecular dynamics impractical for screening purposes. Thus, several approximative methods have been recently developed to study interactions between a protein and ligand during the ligand transport process. Apart from identifying the best binding modes, these methods also provide information on the energetics of the transport and identify problematic regions limiting the ligand passage. These methods use approximations to simulate binding or unbinding events rapidly (calculation times from minutes to hours) and provide energy profiles that can be used to rank ligands or pathways. Here we provide a critical comparison of available methods, showcase their results on sample systems, discuss their practical applications in molecular biotechnologies and outline possible future developments.

Bioactive Phytoconstituents as Potent Inhibitors of Tyrosine-Protein Kinase Yes (YES1): Implications in Anticancer Therapeutics

Tyrosine-protein kinase Yes (YES1) belongs to the Tyrosine-protein kinase family and is involved in several biological activities, including cell survival, cell–cell adhesion, cell differentiation, and cytoskeleton remodeling. It is highly expressed in esophageal, lung, and bladder cancers, and thus considered as an attractive drug target for cancer therapy. In this study, we performed a virtual screening of phytoconstituents from the IMPPAT database to identify potential inhibitors of YES1. Initially, the molecules were retrieved on their physicochemical properties following the Lipinski rule of five. Then binding affinities calculation, PAINS filter, ADMET, and PASS analyses followed by an interaction analysis to select safe and clinically better hits. Finally, two compounds, Glabrene and Lupinisoflavone C (LIC), with appreciable affinities and a specific interaction towards the AlphaFold predicted structure of YES1, were identified. Their time-evolution analyses were carried out using an all-atom molecular dynamics (MD) simulation, principal component analysis, and free energy landscapes. Altogether, we propose that Glabrene and LIC can be further explored in clinical settings to develop anticancer therapeutics targeting YES1 kinase.

Exploring binding stability of hydroxy-3-(4-hydroxyphenyl)-5-(4-nitrophenyl)-5,5a,7,8,9,9a-hexahydrothiazolo[2,3-b] quinazolin-6-one with T790M/L858R EGFR-TKD

Cancer causes innumerable deaths every year globally. Breast cancer and non-small cell lung carcinoma are the most prevalent worldwide. EGFR-TKD is a neoplastic survival therapeutic target in a wide array of carcinoma cells. Various non-specific tyrosine kinase inhibitors lead to hyperphosphorylation and overexpression of EGFR-TKD and further mutations recognise deletion of exon 19. In this work, we study the binding affinity, binding stability, and strength of hydroxy-3-(4-hydroxyphenyl)-5-(4-nitrophenyl)-5,5a,7,8,9,9a-hexahydrothiazolo[2,3-b] quinazolin-6-one with TMLR mutated EGFR-TKD (T790M/L858R). The collective motions, residual mobility, and flexibility of TMLR mutated EGFR-TKD bound with reference and title molecule were calculated by principal component analysis. The meta-state conformations of both the simulated complexes were determined by Gibb's energy landscape analysis. The binding affinity exhibited by thiazolo-[2,3-b] quinazolinone and the reference molecule was found to be -7.95 ± 0.088 Kcal/mol and -9.13 ± 0.018 kcal/mol with TMLR mutated EGFR-TKD. The alignment of both the docked complexes was done by blosum40 matrix. Similar spatial orientations were exhibited by the synthesised ligand in the binding pocket of TMLR mutated EGFR-TKD, corresponding to the reference ligand. The ligand stability was computed for 100 ns. In addition, the radius of gyration, solvent accessible surface area, hydrogen bonds formed was calculated. The average ΔGbind of thiazolo-[2,3-b] quinazolinone was -41.212 ± 0.834 kJ/mol and for reference ligand -71.938 ± 0.367 kJ/mol, calculated by MM-PBSA. ADMET analysis concludes thiazolo-[2,3-b] quinazolinone derivative is safe. Further research work is encouraged to determine the efficacy of thiazolo-[2,3-b] quinazolinone against in vivo models.Communicated by Ramaswamy H. Sarma.

Aquaglyceroporin Modulators as Emergent Pharmacological Molecules for Human Diseases

Aquaglyceroporins, a sub-class of aquaporins that facilitate the diffusion of water, glycerol and other small uncharged solutes across cell membranes, have been recognized for their important role in human physiology and their involvement in multiple disorders, mostly related to disturbed energy homeostasis. Aquaglyceroporins dysfunction in a variety of pathological conditions highlighted their targeting as novel therapeutic strategies, boosting the search for potent and selective modulators with pharmacological properties. The identification of selective inhibitors with potential clinical applications has been challenging, relying on accurate assays to measure membrane glycerol permeability and validate effective functional blockers. Additionally, biologicals such as hormones and natural compounds have been revealed as alternative strategies to modulate aquaglyceroporins via their gene and protein expression. This review summarizes the current knowledge of aquaglyceroporins’ involvement in several pathologies and the experimental approaches used to evaluate glycerol permeability and aquaglyceroporin modulation. In addition, we provide an update on aquaglyceroporins modulators reported to impact disease, unveiling aquaglyceroporin pharmacological targeting as a promising approach for innovative therapeutics.

Mechanistic Insight into Binding of Huperzine A with Human Serum Albumin: Computational and Spectroscopic Approaches

Human serum albumin (HSA) is the most abundant protein in plasma synthesized by the liver and the main modulator of fluid distribution between body compartments. It has an amazing capacity to bind with multiple ligands, offering a store and transporter for various endogenous and exogenous compounds. Huperzine A (HpzA) is a natural sesquiterpene alkaloid found in Huperzia serrata and used in various neurological conditions, including Alzheimer's disease (AD). This study elucidated the binding of HpzA with HSA using advanced computational approaches such as molecular docking and molecular dynamic (MD) simulation followed by fluorescence-based binding assays. The molecular docking result showed plausible interaction between HpzA and HSA. The MD simulation and principal component analysis (PCA) results supported the stable interactions of the protein-ligand complex. The fluorescence assay further validated the in silico study, revealing significant binding affinity between HpzA and HSA. This study advocated that HpzA acts as a latent HSA binding partner, which may be investigated further in AD therapy in experimental settings.

Bioactive phytoconstituents as potent inhibitors of casein kinase-2: dual implications in cancer and COVID-19 therapeutics

Casein kinase 2 (CK2) is a conserved serine/threonine-protein kinase involved in hematopoietic cell survival, cell cycle control, DNA repair, and other cellular processes. It plays a significant role in cancer progression and viral infection. CK2 is considered a potential drug target in cancers and COVID-19 therapy. In this study, we have performed a virtual screening of phytoconstituents from the IMPPAT database to identify some potential inhibitors of CK2. The initial filter was the physicochemical properties of the molecules following the Lipinski rule of five. Then binding affinity calculation, PAINS filter, ADMET, and PASS analyses followed by interaction analysis were carried out to discover nontoxic and better hits. Finally, two compounds, stylopine and dehydroevodiamines with appreciable affinity and specific interaction towards CK2, were identified. Their time-evolution analyses were carried out using all-atom molecular dynamics simulation, principal component analysis and free energy landscape. Altogether, we propose that stylopine and dehydroevodiamines can be further explored in in vitro and in vivo settings to develop anticancer and antiviral therapeutics.

Showing protein–ligands interactions, electrostatic potential of CK2 bound to selected compounds, free energy landscapes of CK2-stylopine, and CK2-dehydroevodiamines complexes.

Identification of Phytoconstituents as Potent Inhibitors of Casein Kinase-1 Alpha Using Virtual Screening and Molecular Dynamics Simulations

Casein kinase-1 alpha (CK1α) is a multifunctional protein kinase that belongs to the serine/threonine kinases of the CK1α family. It is involved in various signaling pathways associated with chromosome segregation, cell metabolism, cell cycle progression, apoptosis, autophagy, etc. It has been known to involve in the progression of many diseases, including cancer, neurodegeneration, obesity, and behavioral disorders. The elevated expression of CK1α in diseased conditions facilitates its selective targeting for therapeutic management. Here, we have performed virtual screening of phytoconstituents from the IMPPAT database seeking potential inhibitors of CK1α. First, a cluster of compounds was retrieved based on physicochemical parameters following Lipinski's rules and PAINS filter. Further, high-affinity hits against CK1α were obtained based on their binding affinity score. Furthermore, the ADMET, PAINS, and PASS evaluation was carried out to select more potent hits. Finally, following the interaction analysis, we elucidated three phytoconstituents, Semiglabrinol, Curcusone_A, and Liriodenine, posturing considerable affinity and specificity towards the CK1α binding pocket. The result was further evaluated by molecular dynamics (MD) simulations, dynamical cross-correlation matrix (DCCM), and principal components analysis (PCA), which revealed that binding of the selected compounds, especially Semiglabrinol, stabilizes CK1α and leads to fewer conformational fluctuations. The MM-PBSA analysis suggested an appreciable binding affinity of all three compounds toward CK1α.

Thermodynamic stability, in-vitro permeability, and in-silico molecular modeling of the optimal Elaeis guineensis leaves extract water-in-oil nanoemulsion

Nanoemulsion is a delivery system used to enhance bioavailability of plant-based compounds across the stratum corneum. Elaeis guineensis leaves are rich source of polyphenolic antioxidants, viz. gallic acid and catechin. The optimal E. guineensis leaves extract water-in-oil nanoemulsion was stable against coalescence, but it was under significant influence of Ostwald ripening over 90 days at 25 °C. The in-vitro permeability revealed a controlled and sustained release of the total phenolic compounds (TPC) of EgLE with a cumulative amount of 1935.0 ± 45.7 µgcm-2 after 8 h. The steady-state flux and permeation coefficient values were 241.9 ± 5.7 µgcm-2 h-1 and 1.15 ± 0.03 cm.h-1, respectively. The kinetic release mechanism for TPC of EgLE was best described by the Korsmeyer-Peppas model due to the highest linearity of R2 = 0.9961, indicating super case II transport mechanism. The in-silico molecular modelling predicted that the aquaporin-3 protein in the stratum corneum bonded preferably to catechin over gallic acid through hydrogen bonds due to the lowest binding energies of - 57.514 kcal/mol and - 8.553 kcal/mol, respectively. Thus, the in-silico study further verified that catechin could improve skin hydration. Therefore, the optimal nanoemulsion could be used topically as moisturizer to enhance skin hydration based on the in-silico prediction.

MAP/Microtubule Affinity Regulating Kinase 4 Inhibitory Potential of Irisin: A New Therapeutic Strategy to Combat Cancer and Alzheimer's Disease

Irisin is a clinically significant protein playing a valuable role in regulating various diseases. Irisin attenuates synaptic and memory dysfunction, highlighting its importance in Alzheimer's disease. On the other hand, Microtubule Affinity Regulating Kinase 4 (MARK4) is associated with various cancer types, uncontrolled neuronal migrations, and disrupted microtubule dynamics. In addition, MARK4 has been explored as a potential drug target for cancer and Alzheimer's disease therapy. Here, we studied the binding and subsequent inhibition of MARK4 by irisin. Irisin binds to MARK4 with an admirable affinity (K = 0.8 × 107 M-1), subsequently inhibiting its activity (IC50 = 2.71 µm). In vitro studies were further validated by docking and simulations. Molecular docking revealed several hydrogen bonds between irisin and MARK4, including critical residues, Lys38, Val40, and Ser134. Furthermore, the molecular dynamic simulation showed that the binding of irisin resulted in enhanced stability of MARK4. This study provides a rationale to use irisin as a therapeutic agent to treat MARK4-associated diseases.

Structure-Based Virtual Screening and Molecular Dynamics Simulation to Identify Potential SARS-CoV-2 Spike Receptor Inhibitors from Natural Compound Database

The outbreak of respiratory disease, COVID-19 caused by SARS-CoV-2 has now been spread globally and the number of new infections is rising every moment. There are no specific medications that are currently available to combat the disease. The spike receptor of SARS-CoV-2 facilitates the viral entry into a host cell and initiation of infection. Targeting the viral entry at the initial stage has a better advantage than inhibiting it in later stages of the viral life cycle. This study deals with identification of the potential natural molecule or its derivatives from MolPort Databank as SARS-CoV-2 spike receptor inhibitors using structure-based virtual screening followed by molecular dynamics simulation. On the basis of ADME properties, docking score, MMGBSAbinding energy, 150 ns molecular docking studies, and final molecular dynamics analysis, two natural compounds - 3 (MolPort-002-535-004) docking score -9.10 kcal mol-1 and 4 (MolPort-005-910-183) docking score -8.5 kcal mol-1, are selected as potential in-silico spike receptor inhibitors. Both hits are commercially available and can be further used for in-vitro and in-vivo studies. Findings of this study can facilitate rational drug design against SARS-CoV-2 spike receptor.