Conjugation with Phospholipids as a Modification Increasing Anticancer Activity of Phenolic Acids in Metastatic Melanoma—In Vitro and In Silico Studies

Phenolic acids possess many beneficial biological activities, including antioxidant and anti-inflammatory properties. Unfortunately, their low bioavailability restricts their potential medical uses, as it limits the concentration of phenolic acids achievable in the organism. The conjugation with phospholipids constitutes one of the most effective strategies to enhance compounds bioavailability in biological systems. In the present study, the conjugates of anisic (ANISA) and veratric acid (VA) with phosphatidylcholine (PC) were investigated. Since both ANISA and VA are inhibitors of tyrosinase, a melanocyte enzyme, the expression of which increases during tumorigenesis, anticancer potential of the conjugates was tested in several metastatic melanoma cell lines. The conjugates proved to be antiproliferative, apoptosis-inducing and cell-cycle-affecting agents, selective for cancerous cells and not affecting normal fibroblasts. The conjugates substituted by ANISA and VA, respectively, at both the sn-1 and sn-2 positions of PC, appeared the most promising, since they were effective against the vast majority of metastatic melanoma cell lines. Additionally, the conjugation of phenolic acids to PC increased their antioxidant activity. Molecular modeling was employed for the first time to estimate the features of the investigated conjugates relevant to their anticancer properties and membrane permeation. Again, the conjugates substituted by phenolic acid at both the sn-1 and sn-2 positions of PC seemed to be presumably most bioavailable.     

Small‐Molecule Mechanism of Action Studies in Caenorhabditis elegans

A general protocol for exogenous small‐molecule pull‐down experiments with Caenorhabditis elegans is described; it provides a link between small‐molecule screens in worms and existing mutant and RNAi technologies, thereby enabling organismal mechanism of action studies for the natural product clovanemagnolol. Forward chemical genetic screens followed by mechanism of action studies with C. elegans, when coupled with genetic validation of identified targets to reproduce the small molecule's phenotypic effects, provide a unique platform for discovering the biological targets of compounds that affect multicellular processes. First, the use of an immobilized FK506 derivative and soluble competition experiments with optimally prepared soluble C. elegans proteome successfully identified interactions with FK506 binding proteins 1 to 6. This approach was used to determine an unknown mechanism of action for clovanemagnolol, a small molecule that promotes axonal branching in both primary neuronal cultures and in vivo in C. elegans. Following the synthesis of an appropriately functionalized solid‐phase reagent bearing a clovanemagnolol analogue pull‐down experiments employing soluble competition identified kinesin light chain‐1 (KLC‐1), a protein involved in axonal cargo transport, as a putative target. This was corroborated through the use of mutant worms lacking klc‐1 and possessing GFP neuronal labeling, reproducing the axonal branching phenotype induced by the small molecule clovanemagnolol.     

Recent Advances on P-Glycoprotein (ABCB1) Transporter Modelling with In Silico Methods

ABC transporters play a critical role in both drug bioavailability and toxicity, and with the discovery of the P-glycoprotein (P-gp), this became even more evident, as it plays an important role in preventing intracellular accumulation of toxic compounds. Over the past 30 years, intensive studies have been conducted to find new therapeutic molecules to reverse the phenomenon of multidrug resistance (MDR) ), that research has found is often associated with overexpression of P-gp, the most extensively studied drug efflux transporter; in MDR, therapeutic drugs are prevented from reaching their targets due to active efflux from the cell. The development of P-gp inhibitors is recognized as a good way to reverse this type of MDR, which has been the subject of extensive studies over the past few decades. Despite the progress made, no effective P-gp inhibitors to reverse multidrug resistance are yet on the market, mainly because of their toxic effects. Computational studies can accelerate this process, and in silico models such as QSAR models that predict the activity of compounds associated with P-gp (or analogous transporters) are of great value in the early stages of drug development, along with molecular modelling methods, which provide a way to explain how these molecules interact with the ABC transporter. This review highlights recent advances in computational P-gp research, spanning the last five years to 2022. Particular attention is given to the use of machine-learning approaches, drug–transporter interactions, and recent discoveries of potential P-gp inhibitors that could act as modulators of multidrug resistance.     

Methods and Applications of In Silico Aptamer Design and Modeling

Aptamers are nucleic acid analogues of antibodies with high affinity to different targets, such as cells, viruses, proteins, inorganic materials, and coenzymes. Empirical approaches allow the design of in vitro aptamers that bind particularly to a target molecule with high affinity and selectivity. Theoretical methods allow significant expansion of the possibilities of aptamer design. In this study, we review theoretical and joint theoretical-experimental studies dedicated to aptamer design and modeling. We consider aptamers with different targets, such as proteins, antibiotics, organophosphates, nucleobases, amino acids, and drugs. During nucleic acid modeling and in silico design, a full set of in silico methods can be applied, such as docking, molecular dynamics (MD), and statistical analysis. The typical modeling workflow starts with structure prediction. Then, docking of target and aptamer is performed. Next, MD simulations are performed, which allows for an evaluation of the stability of aptamer/ligand complexes and determination of the binding energies with higher accuracy. Then, aptamer/ligand interactions are analyzed, and mutations of studied aptamers made. Subsequently, the whole procedure of molecular modeling can be reiterated. Thus, the interactions between aptamers and their ligands are complex and difficult to understand using only experimental approaches. Docking and MD are irreplaceable when aptamers are studied in silico.     

In Silico Discovery of Anticancer Peptides from Sanghuang

Anticancer peptide (ACP) is a short peptide with less than 50 amino acids that has been discovered in a variety of foods. It has been demonstrated that traditional Chinese medicine or food can help treat cancer in some cases, which suggests that ACP may be one of the therapeutic ingredients. Studies on the anti-cancer properties of Sanghuangporus sanghuang have concentrated on polysaccharides, flavonoids, triterpenoids, etc. The function of peptides has not received much attention. The purpose of this study is to use computer mining techniques to search for potential anticancer peptides from 62 proteins of Sanghuang. We used mACPpred to perform sequence scans after theoretical trypsin hydrolysis and discovered nine fragments with an anticancer probability of over 0.60. The study used AlphaFold 2 to perform structural modeling of the first three ACPs discovered, which had blast results from the Cancer PPD database. Using reverse docking technology, we found the target proteins and interacting residues of two ACPs with an unknown mechanism. Reverse docking results predicted the binding modes of the ACPs and their target protein. In addition, we determined the active part of ACPs by quantum chemical calculation. Our study provides a framework for the future discovery of functional peptides from foods. The ACPs discovered have the potential to be used as drugs in oncology clinical treatment after further research.     

In Silico Analysis of Enzyme Surface and Glycosylation Effect as a Tool for Efficient Covalent Immobilisation of CalB and PGA on Sepabeads®

This study presents a computational analysis of the structures of lipase B from Candida antarctica (CalB) and two penicillin G acylases (PGAs), from eukaryotic and prokaryotic sources, respectively. Molecular simulations were used to point out the regions of the enzymes that are prone to interact with immobilisation supports. In order to evaluate the accessibility of the active site, the location of the amino acid residues involved in the formation of covalent bonds with the polymers was visualised. The mapping of the distribution of hydrophobic and hydrophilic regions on the enzyme surface provided a view of the areas of the protein that can establish either hydrophobic or hydrophilic interactions with the carriers. Experimental data obtained from the immobilisation of the enzymes on supports bearing different chemical functionalities suggest the involvement of the glycan moiety in enzyme–polymer interactions. In the case of PGA the glycan moiety can constitute an extra site for the covalent linkage of the enzyme on the polymer.     

In Silico and Experimental ADAM17 Kinetic Modeling as Basis for Future Screening System for Modulators

Understanding the mechanisms of modulators’ action on enzymes is crucial for optimizing and designing pharmaceutical substances. The acute inflammatory response, in particular, is regulated mainly by a disintegrin and metalloproteinase (ADAM) 17. ADAM17 processes several disease mediators such as TNFα and APP, releasing their soluble ectodomains (shedding). A malfunction of this process leads to a disturbed inflammatory response. Chemical protease inhibitors such as TAPI-1 were used in the past to inhibit ADAM17 proteolytic activity. However, due to ADAM17′s broad expression and activity profile, the development of active-site-directed ADAM17 inhibitor was discontinued. New ‘exosite’ (secondary substrate binding site) inhibitors with substrate selectivity raised the hope of a substrate-selective modulation as a promising approach for inflammatory disease therapy. This work aimed to develop a high-throughput screen for potential ADAM17 modulators as therapeutic drugs. By combining experimental and in silico methods (structural modeling and docking), we modeled the kinetics of ADAM17 inhibitor. The results explain ADAM17 inhibition mechanisms and give a methodology for studying selective inhibition towards the design of pharmaceutical substances with higher selectivity.     

In Silico Analysis of the Molecular-Level Impact of SMPD1 Variants on Niemann-Pick Disease Severity

Sphingomyelin phosphodiesterase (SMPD1) is a key enzyme in the sphingolipid metabolism. Genetic SMPD1 variants have been related to the Niemann-Pick lysosomal storage disorder, which has different degrees of phenotypic severity ranging from severe symptomatology involving the central nervous system (type A) to milder ones (type B). They have also been linked to neurodegenerative disorders such as Parkinson and Alzheimer. In this paper, we leveraged structural, evolutionary and stability information on SMPD1 to predict and analyze the impact of variants at the molecular level. We developed the SMPD1-ZooM algorithm, which is able to predict with good accuracy whether variants cause Niemann-Pick disease and its phenotypic severity; the predictor is freely available for download. We performed a large-scale analysis of all possible SMPD1 variants, which led us to identify protein regions that are either robust or fragile with respect to amino acid variations, and show the importance of aromatic-involving interactions in SMPD1 function and stability. Our study also revealed a good correlation between SMPD1-ZooM scores and in vitro loss of SMPD1 activity. The understanding of the molecular effects of SMPD1 variants is of crucial importance to improve genetic screening of SMPD1-related disorders and to develop personalized treatments that restore SMPD1 functionality.     

Expression of Huntingtin and TDP-43 Derivatives in Fission Yeast Can Cause Both Beneficial and Toxic Effects

Many neurodegenerative disorders display protein aggregation as a hallmark, Huntingtin and TDP-43 aggregates being characteristic of Huntington disease and amyotrophic lateral sclerosis, respectively. However, whether these aggregates cause the diseases, are secondary by-products, or even have protective effects, is a matter of debate. Mutations in both human proteins can modulate the structure, number and type of aggregates, as well as their toxicity. To study the role of protein aggregates in cellular fitness, we have expressed in a highly tractable unicellular model different variants of Huntingtin and TDP-43. They each display specific patterns of aggregation and toxicity, even though in both cases proteins have to be very highly expressed to affect cell fitness. The aggregation properties of Huntingtin, but not of TDP-43, are affected by chaperones such as Hsp104 and the Hsp40 couple Mas5, suggesting that the TDP-43, but not Huntingtin, derivatives have intrinsic aggregation propensity. Importantly, expression of the aggregating form of Huntingtin causes a significant extension of fission yeast lifespan, probably as a consequence of kidnapping chaperones required for maintaining stress responses off. Our study demonstrates that in general these prion-like proteins do not cause toxicity under normal conditions, and in fact they can protect cells through indirect mechanisms which up-regulate cellular defense pathways.     

A Study on the Effect of the Substituent against PAK4 Inhibition Using In Silico Methods

The intrinsic inductive properties of atoms or functional groups depend on the chemical properties of either electron-withdrawing groups (EWGs) or electron-donating groups (EDGs). This study aimed to evaluate in silico methods to determine whether changes in chemical properties of the compound by single atomic substitution affect the biological activity of target proteins and whether the results depend on the properties of the functional groups. We found an imidazo[4,5-b]pyridine-based PAK4 inhibitor, compound 1, as an initial hit compound with the well-defined binding mode for PAK4. In this study, we used both experimental and in silico methods to investigate the effect of atomic substitution on biological activity to optimize the initial hit compound. In biological assays, in the case of EWG, as the size of the halogen atom became smaller and the electronegativity increased, the biological activity IC₅₀ value ranged from 5150 nM to inactive; in the case of EDG, biological activity was inactive. Furthermore, we analyzed the interactions of PAK4 with compounds, focusing on the hinge region residues, L398 and E399, and gatekeeper residues, M395 and K350, of the PAK4 protein using molecular docking studies and fragment molecular orbital (FMO) methods to determine the differences between the effect of EWG and EDG on the activity of target proteins. These results of the docking score and binding energy did not explain the differences in biological activity. However, the pair-interaction energy obtained from the results of the FMO method indicated that there was a difference in the interaction energy between the EWG and EDG in the hinge region residues, L398 and E399, as well as in M395 and K350. The two groups with different properties exhibited opposite electrostatic energy and charge transfer energy between L398 and E399. Additionally, we investigated the electron distribution of the parts interacting with the hinge region by visualizing the molecular electrostatic potential (MEP) surface of the compounds. In conclusion, we described the properties of functional groups that affect biological activity using an in silico method, FMO.     

In Silico Study of Allosteric Communication Networks in GPCR Signaling Bias

Signaling bias is a promising characteristic of G protein-coupled receptors (GPCRs) as it provides the opportunity to develop more efficacious and safer drugs. This is because biased ligands can avoid the activation of pathways linked to side effects whilst still producing the desired therapeutic effect. In this respect, a deeper understanding of receptor dynamics and implicated allosteric communication networks in signaling bias can accelerate the research on novel biased drug candidates. In this review, we aim to provide an overview of computational methods and techniques for studying allosteric communication and signaling bias in GPCRs. This includes (i) the detection of allosteric communication networks and (ii) the application of network theory for extracting relevant information pipelines and highly communicated sites in GPCRs. We focus on the most recent research and highlight structural insights obtained based on the framework of allosteric communication networks and network theory for GPCR signaling bias.     

Involvement of Huntingtin in Development and Ciliary Beating Regulation of Larvae of the Sea Urchin,Hemicentrotus pulcherrimus

The multiple functions of the wild type Huntington’s disease protein of the sea urchin Hemicentrotus pulcherrimus (Hp-Htt) have been examined using the anti-Hp-Htt antibody (Ab) raised against synthetic oligopeptides. According to immunoblotting, Hp-Htt was detected as a single band at around the 350 kDa region at the swimming blastula stage to the prism larva stage. From the 2-arm pluteus stage (2aPL), however, an additional smaller band at the 165 kDa region appeared. Immunohistochemically, Hp-Htt was detected in the nuclei and the nearby cytoplasm of the ectodermal cells from the swimming blastula stage, and the blastocoelar cells from the mid-gastrula stage. The Ab-positive signal was converged to the ciliary band-associated strand (CBAS). There, it was accompanied by several CBAS-marker proteins in the cytoplasm, such as glutamate decarboxylase. Application of Hp-Htt morpholino (Hp-Htt-MO) has resulted in shortened larval arms, accompanied by decreased 5-bromo-2-deoxyuridin (BrdU) incorporation by the ectodermal cells of the larval arms. Hp-Htt-MO also resulted in lowered ciliary beating activity, accompanied by a disordered swirling pattern formation around the body. These Hp-Htt-MO-induced deficiencies took place after the onset of CBAS system formation at the larval arms. Thus, Hp-Htt is involved in cell proliferation and the ciliary beating pattern regulation signaling system in pluteus larvae.     

In Silico Prediction of the Metabolic Resistance of Vitamin D Analogs against CYP3A4 Metabolizing Enzyme

The microsomal cytochrome P450 3A4 (CYP3A4) and mitochondrial cytochrome P450 24A1 (CYP24A1) hydroxylating enzymes both metabolize vitamin D and its analogs. The three-dimensional (3D) structure of the full-length native human CYP3A4 has been solved, but the respective structure of the main vitamin D hydroxylating CYP24A1 enzyme is unknown. The structures of recombinant CYP24A1 enzymes have been solved; however, from studies of the vitamin D receptor, the use of a truncated protein for docking studies of ligands led to incorrect results. As the structure of the native CYP3A4 protein is known, we performed rigid docking supported by molecular dynamic simulation using CYP3A4 to predict the metabolic conversion of analogs of 1,25-dihydroxyvitamin D2 (1,25D2). This is highly important to the design of novel vitamin D-based drug candidates of reasonable metabolic stability as CYP3A4 metabolizes ca. 50% of the drug substances. The use of the 3D structure data of human CYP3A4 has allowed us to explain the substantial differences in the metabolic conversion of the side-chain geometric analogs of 1,25D2. The calculated free enthalpy of the binding of an analog of 1,25D2 to CYP3A4 agreed with the experimentally observed conversion of the analog by CYP24A1. The metabolic conversion of an analog of 1,25D2 to the main vitamin D hydroxylating enzyme CYP24A1, of unknown 3D structure, can be explained by the binding strength of the analog to the known 3D structure of the CYP3A4 enzyme.     

Comprehensive In Silico Functional Prediction Analysis of CDKL5 by Single Amino Acid Substitution in the Catalytic Domain

Cyclin-dependent kinase-like 5 (CDKL5) is a serine/threonine protein kinase whose pathological mutations cause CDKL5 deficiency disorder. Most missense mutations are concentrated in the catalytic domain. Therefore, anticipating whether mutations in this region affect CDKL5 function is informative for clinical diagnosis. This study comprehensively predicted the pathogenicity of all 5700 missense substitutions in the catalytic domain of CDKL5 using in silico analysis and evaluating their accuracy. Each missense substitution was evaluated as “pathogenic” or “benign”. In silico tools PolyPhen-2 HumDiv mode/HumVar mode, PROVEAN, and SIFT were selected individually or in combination with one another to determine their performance using 36 previously reported mutations as a reference. Substitutions predicted as pathogenic were over 88.0% accurate using each of the three tools. The best performance score (accuracy, 97.2%; sensitivity, 100%; specificity, 66.7%; and Matthew’s correlation coefficient (MCC), 0.804) was achieved by combining PolyPhen-2 HumDiv, PolyPhen-2 HumVar, and PROVEAN. This provided comprehensive information that could accurately predict the pathogenicity of the disease, which might be used as an aid for clinical diagnosis.     

Screening of Novel Bioactive Peptides from Goat Casein: In Silico to In Vitro Validation

Food-derived bioactive peptides are of great interest to science and industry due to evolving drivers of food product innovation, including health and wellness. This study aims to draw attention through a critical study on how bioinformatics analysis is employed in the identification of bioactive peptides in the laboratory. An in silico analysis (PeptideRanker, BIOPEP, AHTpin, and mAHTPred) of a list of peptides from goat casein hydrolysate was performed to predict which sequences could potentially be bioactive. To validate the predictions, the in vitro antihypertensive potential of the five peptides with the highest potential was first measured. Then, for three of these, gastrointestinal digestion was simulated in vitro, followed by the analysis of the resulting ACE inhibitory activity as well as antioxidant capacity. We thus observed that the use of new computational biology technologies to predict peptide sequences is an important research tool, but they should not be used alone and complementarity with various in vitro and in vivo assays is essential.     

In Vitro and In Silico Characterization of an Antimalarial Compound with Antitumor Activity Targeting Human DNA Topoisomerase IB

Human DNA topoisomerase IB controls the topological state of supercoiled DNA through a complex catalytic cycle that consists of cleavage and religation reactions, allowing the progression of fundamental DNA metabolism. The catalytic steps of human DNA topoisomerase IB were analyzed in the presence of a drug, obtained by the open-access drug bank Medicines for Malaria Venture. The experiments indicate that the compound strongly and irreversibly inhibits the cleavage step of the enzyme reaction and reduces the cell viability of three different cancer cell lines. Molecular docking and molecular dynamics simulations suggest that the drug binds to the human DNA topoisomerase IB-DNA complex sitting inside the catalytic site of the enzyme, providing a molecular explanation for the cleavage-inhibition effect. For all these reasons, the aforementioned drug could be a possible lead compound for the development of an efficient anti-tumor molecule targeting human DNA topoisomerase IB.     

In Silico Food-Drug Interaction: A Case Study of Eluxadoline and Fatty Meal

Food-drug interaction is an infrequently considered aspect in clinical practice. Usually, drugs are taken together with meals and what follows may adversely affect pharmacokinetic and pharmacodynamic properties, and hence, the therapeutic effects. In this study, a computational protocol was proposed to explain the different assimilations of two µ-receptors agonists, eluxadoline and loperamide, with a peculiar pharmacokinetic profile. Compared to loperamide, eluxadoline is absorbed less after the intake of a fatty meal, and the LogP values do not explain this event. Firstly, keeping in mind the different pH in the intestinal tract, the protonation states of both compounds were calculated. Then, all structures were subjected to a conformational search by using MonteCarlo and Molecular Dynamics methods, with solvation terms mimicking the water and weak polar solvent (octanol). Both computational results showed that eluxadoline has less conformational freedom in octanol, unlike loperamide, which exhibits constant behavior in both solvents. Therefore, we hypothesize that fatty meal causes the “closure” of the eluxadoline molecule to prevent the exposure of the polar groups and their interaction with water, necessary for the drug absorption. Based on our results, this work could be a reasonable “case study”, useful for future investigation of the drug pharmacokinetic profile.     

Pharmacological and In Silico Analysis of Oat Avenanthramides as EGFR Inhibitors: Effects on EGF-Induced Lung Cancer Cell Growth and Migration

Avena sativa L. is a wholegrain cereal and an important edible crop. Oats possesses high nutritional and health promoting values and contains high levels of bioactive compounds, including a group of phenolic amides, named avenanthramides (Avns), exerting antioxidant, anti-inflammatory, and anticancer activities. Epidermal growth factor receptor (EGFR) represents one of the most known oncogenes and it is frequently up-regulated or mutated in human cancers. The oncogenic effects of EGFR include enhanced cell growth, angiogenesis, and metastasis, and down-regulation or inhibition of EGFR signaling has therapeutic benefit. Front-line EGFR tyrosine kinase inhibitor therapy is the standard therapy for patients with EGFR-mutated lung cancer. However, the clinical effects of EGFR inhibition may be lost after a few months of treatment due to the onset of resistance. Here, we showed the anticancer activity of Avns, focusing on EGFR activation and signaling pathway. Lung cancer cellular models have been used to evaluate the activity of Avns on tumor growth, migration, EMT, and anoikis induced by EGF. In addition, docking and molecular dynamics simulations showed that the Avns bind with high affinity to a region in the vicinity of αC-helix and the DGF motif of EGFR, jeopardizing the target biological function. Altogether, our results reveal a new pharmacological activity of Avns as EGFR tyrosine kinase inhibitors.     

Antiviral Activity of Two Marine Carotenoids against SARS-CoV-2 Virus Entry In Silico and In Vitro

The marine carotenoids fucoxanthin and siphonaxanthin are powerful antioxidants that are attracting focused attention to identify a variety of health benefits and industry applications. In this study, the binding energy of these carotenoids with the SARS-CoV-2 Spike-glycoprotein was predicted by molecular docking simulation, and their inhibitory activity was confirmed with SARS-CoV-2 pseudovirus on HEK293 cells overexpressing angiotensin-converting enzyme 2 (ACE2). Siphonaxanthin from Codium fragile showed significant antiviral activity with an IC50 of 87.4 μM against SARS-CoV-2 pseudovirus entry, while fucoxanthin from Undaria pinnatifida sporophyll did not. The acute toxicities were predicted to be relatively low, and pharmacokinetic predictions indicate GI absorption. Although further studies are needed to elucidate the inhibition of viral infection by siphonaxanthin, these results provide useful information in the application of these marine carotenoids for the treatment and prevention of COVID-19.