VirtuousPocketome: a computational tool for screening protein-ligand complexes to identify similar binding sites.

Protein residues within binding pockets play a critical role in determining the range of ligands that can interact with a protein, influencing its structure and function. Identifying structural similarities in proteins offers valuable insights into their function and activation mechanisms, aiding in predicting protein-ligand interactions, anticipating off-target effects, and facilitating the development of therapeutic agents. Numerous computational methods assessing global or local similarity in protein cavities have emerged, but their utilization is impeded by complexity, impractical automation for amino acid pattern searches, and an inability to evaluate the dynamics of scrutinized protein-ligand systems. Here, we present a general, automatic and unbiased computational pipeline, named VirtuousPocketome, aimed at screening huge databases of proteins for similar binding pockets starting from an interested protein-ligand complex. We demonstrate the pipeline's potential by exploring a recently-solved human bitter taste receptor, i.e. the TAS2R46, complexed with strychnine. We pinpointed 145 proteins sharing similar binding sites compared to the analysed bitter taste receptor and the enrichment analysis highlighted the related biological processes, molecular functions and cellular components. This work represents the foundation for future studies aimed at understanding the effective role of tastants outside the gustatory system: this could pave the way towards the rationalization of the diet as a supplement to standard pharmacological treatments and the design of novel tastants-inspired compounds to target other proteins involved in specific diseases or disorders. The proposed pipeline is publicly accessible, can be applied to any protein-ligand complex, and could be expanded to screen any database of protein structures.

Mechanical communication within the microtubule through network-based analysis of tubulin dynamics.

The identification of the mechanisms underlying the transfer of mechanical vibrations in protein complexes is crucial to understand how these super-assemblies are stabilized to perform specific functions within the cell. In this context, the study of the structural communication and the propagation of mechanical stimuli within the microtubule (MT) is important given the pivotal role of the latter in cell viability. In this study, we employed molecular modelling and the dynamical network analysis approaches to analyse the MT. The results highlight that ß -tubulin drives the transfer of mechanical information between protofilaments (PFs), which is altered at the seam due to a different interaction pattern. Moreover, while the key residues involved in the structural communication along the PF are generally conserved, a higher diversity was observed for amino acids mediating the lateral communication. Taken together, these results might explain why MTs with different PF numbers are formed in different organisms or with different ß -tubulin isotypes.

Toward a holographic brain paradigm: a lipid-centric model of brain functioning.

Due to the stimulation of neuronal membrane dipoles by action potentials, under suitable conditions coherent dipole oscillations can be formed. We argue that these dipole oscillations satisfy the weak Bose-Einstein condensate criteria of the Froehlich model of biological coherence. They can subsequently generate electromagnetic fields (EMFs) propagating in the inter-neuronal space. When neighboring neurons fire synchronously, EMFs can create interference patterns and hence form holographic images containing analog information about the sensory inputs that trigger neuronal activity. The mirror pattern projected by EMFs inside the neuron can encode information in the neuronal cytoskeleton. We outline an experimental verification of our hypothesis and its consequences for anesthesia, neurodegenerative diseases, and psychiatric states.

Toward a general and interpretable umami taste predictor using a multi-objective machine learning approach.

The umami taste is one of the five basic taste modalities normally linked to the protein content in food. The implementation of fast and cost-effective tools for the prediction of the umami taste of a molecule remains extremely interesting to understand the molecular basis of this taste and to effectively rationalise the production and consumption of specific foods and ingredients. However, the only examples of umami predictors available in the literature rely on the amino acid sequence of the analysed peptides, limiting the applicability of the models. In the present study, we developed a novel ML-based algorithm, named VirtuousUmami, able to predict the umami taste of a query compound starting from its SMILES representation, thus opening up the possibility of potentially using such a model on any database through a standard and more general molecular description. Herein, we have tested our model on five databases related to foods or natural compounds. The proposed tool will pave the way toward the rationalisation of the molecular features underlying the umami taste and toward the design of specific peptide-inspired compounds with specific taste properties.

Informed classification of sweeteners/bitterants compounds via explainable machine learning.

Perception of taste is an emergent phenomenon arising from complex molecular interactions between chemical compounds and specific taste receptors. Among all the taste perceptions, the dichotomy of sweet and bitter tastes has been the subject of several machine learning studies for classification purposes. While previous studies have provided accurate sweeteners/bitterants classifiers, there is ample scope to enhance these models by enriching the understanding of the molecular basis of bitter-sweet tastes. Towards these goals, our study focuses on the development and testing of several machine learning strategies coupled with the novel SHapley Additive exPlanations (SHAP) for a rational sweetness/bitterness classification. This allows the identification of the chemical descriptors of interest by allowing a more informed approach toward the rational design and screening of sweeteners/bitterants. To support future research in this field, we make all datasets and machine learning models publicly available and present an easy-to-use code for bitter-sweet taste prediction.

Toward the design and development of peptidomimetic inhibitors of the Ataxin-1 aggregation pathway.

Spinocerebellar ataxia type 1 is a degenerative disorder caused by polyglutamine expansions and aggregation of Ataxin-1. The interaction between Capicua (CIC) and the AXH domain of Ataxin-1 protein has been suggested as a possible driver of aggregation for the expanded Ataxin-1 protein and the subsequent onset of spinocerebellar ataxia 1. Experimental studies have demonstrated that short constructs of CIC may prevent such aggregation and suggested this as a possible candidate to inspire the rational design of peptidomimetics. In this work, molecular modeling techniques, namely the alchemical mutation and force field-based molecular dynamics, have been employed to propose a pipeline for the rational design of a CIC-inspired inhibitor of the ataxin-1 aggregation pathway. In particular, this study has shown that the alchemical mutation can estimate the affinity between AXH and CIC with good correlation with experimental data, while molecular dynamics shed light on molecular mechanisms that occur for stabilization of the interaction between the CIC-inspired construct and the AXH domain of Ataxin-1. This work lays the foundation for a rational methodology for the in silico screening and design of peptidomimetics, which can expedite and streamline experimental studies to identify strategies for inhibiting the ataxin-1 aggregation pathway.

Reconstitution of microtubule into GTP-responsive nanocapsules.

Nanocapsules that collapse in response to guanosine triphosphate (GTP) have the potential as drug carriers for efficiently curing diseases caused by cancer and RNA viruses because GTP is present at high levels in such diseased cells and tissues. However, known GTP-responsive carriers also respond to adenosine triphosphate (ATP), which is abundant in normal cells as well. Here, we report the elaborate reconstitution of microtubule into a nanocapsule that selectively responds to GTP. When the tubulin monomer from microtubule is incubated at 37 °C with a mixture of GTP (17 mol%) and nonhydrolysable GTP* (83 mol%), a tubulin nanosheet forms. Upon addition of photoreactive molecular glue to the resulting dispersion, the nanosheet is transformed into a nanocapsule. Cell death results when a doxorubicin-containing nanocapsule, after photochemically crosslinked for properly stabilizing its shell, is taken up into cancer cells that overexpress GTP.

Alteration of lipid bilayer mechanics by volatile anesthetics: Insights from µs-long molecular dynamics simulations.

Very few drugs in clinical practice feature the chemical diversity, narrow therapeutic window, unique route of administration, and reversible cognitive effects of volatile anesthetics. The correlation between their hydrophobicity and their potency and the increasing amount of evidence suggesting that anesthetics exert their action on transmembrane proteins, justifies the investigation of their effects on phospholipid bilayers at the molecular level, given the strong functional and structural link between transmembrane proteins and the surrounding lipid matrix. Molecular dynamics simulations of a model lipid bilayer in the presence of ethylene, desflurane, methoxyflurane, and the nonimmobilizer 1,2-dichlorohexafluorocyclobutane (also called F6 or 2N) at different concentrations highlight the structural consequences of VA partitioning in the lipid phase, with a decrease of lipid order and bilayer thickness, an increase in overall lipid lateral mobility and area-per-lipid, and a marked reduction in the mechanical stiffness of the membrane, that strongly correlates with the compounds' hydrophobicity.

Cardiovascular risk prediction: from classical statistical methods to machine learning approaches.

Nowadays, cardiovascular risk prediction scores are commonly used in primary prevention settings. Estimating the cardiovascular individual risk is of crucial importance for effective patient management and optimal therapy identification, with relevant consequences on secondary prevention settings. To reach this goal, a plethora of risk scores have been developed in the past, most of them assuming that each cardiovascular risk factor is linearly dependent on the outcome. However, the overall accuracy of these methods often remains insufficient to solve the problem at hand. In this scenario, machine learning techniques have repeatedly proved successful in improving cardiovascular risk predictions, being able to capture the non-linearity present in the data. In this concern, we present a detailed discussion concerning the application of classical versus machine learning-based cardiovascular risk scores in the clinical setting. This review aimed to give an overview of the current risk scores based on classical statistical approaches and machine learning techniques applied to predict the risk of several cardiovascular diseases, comparing them, discussing their similarities and differences, and highlighting their main drawbacks to aid the physician having a more critical understanding of these tools.

ALS2-Related Motor Neuron Diseases: From Symptoms to Molecules.

Infantile-onset Ascending Hereditary Spastic Paralysis, Juvenile Primary Lateral Sclerosis and Juvenile Amyotrophic Lateral Sclerosis are all motor neuron diseases related to mutations on the ALS2 gene, encoding for a 1657 amino acids protein named Alsin. This ~185 kDa multi-domain protein is ubiquitously expressed in various human tissues, mostly in the brain and the spinal cord. Several investigations have indicated how mutations within Alsin's structured domains may be responsible for the alteration of Alsin's native oligomerization state or Alsin's propensity to interact with protein partners. In this review paper, we propose a description of differences and similarities characterizing the above-mentioned ALS2-related rare neurodegenerative disorders, pointing attention to the effects of ALS2 mutation from molecule to organ and at the system level. Known cases were collected through a literature review and rationalized to deeply elucidate the neurodegenerative clinical outcomes as consequences of ALS2 mutations.

Fragmented blind docking: a novel protein-ligand binding prediction protocol.

In the present paper we propose a novel blind docking protocol based on Autodock-Vina. The developed docking protocol can provide binding site identification and binding pose prediction at the same time, by a systematical exploration of the protein volume performed with several preliminary docking calculations. In our opinion, this protocol can be successfully applied during the first steps of the virtual screening pipeline, because it provides binding site identification and binding pose prediction at the same time without visual evaluation of the binding site. After the binding pose prediction, MM/GBSA re-scoring rescoring procedures has been applied to improve the accuracy of the protein-ligand bound state. The FRAD protocol has been tested on 116 protein-ligand complexes of the Heat Shock Protein 90 - alpha, on 176 of Human Immunodeficiency virus protease 1, and on more than 100 protein-ligand system taken from the PDBbind dataset. Overall, the FRAD approach combined to MM/GBSA re-scoring can be considered as a powerful tool to increase the accuracy and efficiency with respect to other standard docking approaches when the ligand-binding site is unknown.Communicated by Ramaswamy H. Sarma.

Insights into the interaction dynamics between volatile anesthetics and tubulin through computational molecular modelling.

General anesthetics, able to reversibly suppress all conscious brain activity, have baffled medical science for decades, and little is known about their exact molecular mechanism of action. Given the recent scientific interest in the exploration of microtubules as putative functional targets of anesthetics, and the involvement thereof in neurodegenerative disorders, the present work focuses on the investigation of the interaction between human tubulin and four volatile anesthetics: ethylene, desflurane, halothane and methoxyflurane. Interaction sites on different tubulin isotypes are predicted through docking, along with an estimate of the binding affinity ranking. The analysis is expanded by Molecular Dynamics simulations, where the dimers are allowed to freely interact with anesthetics in the surrounding medium. This allowed for the determination of interaction hotspots on tubulin dimers, which could be linked to different functional consequences on the microtubule architecture, and confirmed the weak, Van der Waals-type interaction, occurring within hydrophobic pockets on the dimer. Both docking and MD simulations highlighted significantly weaker interactions of ethylene, consistent with its far lower potency as a general anesthetic. Overall, simulations suggest a transient interaction between anesthetics and microtubules in general anesthesia, and contact probability analysis shows interaction strengths consistent with the potencies of the four compounds.Communicated by Ramaswamy H. Sarma.

Genetic, cellular, and structural characterization of the membrane potential-dependent cell-penetrating peptide translocation pore.

Cell-penetrating peptides (CPPs) allow intracellular delivery of bioactive cargo molecules. The mechanisms allowing CPPs to enter cells are ill-defined. Using a CRISPR/Cas9-based screening, we discovered that KCNQ5, KCNN4, and KCNK5 potassium channels positively modulate cationic CPP direct translocation into cells by decreasing the transmembrane potential (Vm). These findings provide the first unbiased genetic validation of the role of Vm in CPP translocation in cells. In silico modeling and live cell experiments indicate that CPPs, by bringing positive charges on the outer surface of the plasma membrane, decrease the Vm to very low values (-150 mV or less), a situation we have coined megapolarization that then triggers formation of water pores used by CPPs to enter cells. Megapolarization lowers the free energy barrier associated with CPP membrane translocation. Using dyes of varying dimensions in CPP co-entry experiments, the diameter of the water pores in living cells was estimated to be 2 (-5) nm, in accordance with the structural characteristics of the pores predicted by in silico modeling. Pharmacological manipulation to lower transmembrane potential boosted CPP cellular internalization in zebrafish and mouse models. Besides identifying the first proteins that regulate CPP translocation, this work characterized key mechanistic steps used by CPPs to cross cellular membranes. This opens the ground for strategies aimed at improving the ability of cells to capture CPP-linked cargos in vitro and in vivo.

Before a drug can have its desired effect, it must reach its target tissue or organ, and enter its cells. This is not easy because cells are surrounded by the plasma membrane, a fat-based barrier that separates the cell from its external environment. The plasma membrane contains proteins that act as channels, shuttling specific molecules in and out of the cell, and it also holds charge, with its inside surface being more negatively charged than its outside surface. Cell-penetrating peptides are short sequences of amino acids (the building blocks that form proteins) that carry positive charges. These positive charges allow them to cross the membrane easily, but it is not well understood how. To find out how cell-penetrating peptides cross the membrane, Trofimenko et al. attached them to dyes of different sizes. This revealed that the cell-penetrating peptides enter the cell through temporary holes called water pores, which measure about two nanometres across. The water pores form when the membrane becomes 'megapolarized', this is, when the difference in charge between the inside and the outside of the membrane becomes greater than normal. This can happen when the negative charge on the inside surface or the positive charge on the outer surface of the membrane increase. Megapolarization depends on potassium channels, which transport positive potassium ions outside the cell, making the outside of the membrane positive. When cell-penetrating peptides arrive at the outer surface of the cell near potassium channels, they make it even more positive. This increases the charge difference between the inside and the outside of the cell, allowing water pores to form. Once the peptides pass through the pores, the charge difference between the inside and the outside of the cell membrane dissipates, and the pores collapse. Drug developers are experimenting with attaching cell-penetrating peptides to drugs to help them get inside their target cells. Currently there are several experimental medications of this kind in clinical trials. Understanding how these peptides gain entry, and what size of molecule they could carry with them, provides solid ground for further drug development.

Computational Study of Potential Galectin-3 Inhibitors in the Treatment of COVID-19.

Galectin-3 is a carbohydrate-binding protein and the most studied member of the galectin family. It regulates several functions throughout the body, among which are inflammation and post-injury remodelling. Recent studies have highlighted the similarity between Galectin-3's carbohydrate recognition domain and the so-called "galectin fold" present on the N-terminal domain of the S1 sub-unit of the SARS-CoV-2 spike protein. Sialic acids binding to the N-terminal domain of the Spike protein are known to be crucial for viral entry into humans, and the role of Galectin-3 as a mediator of lung fibrosis has long been the object of study since its levels have been found to be abnormally high in alveolar macrophages following lung injury. In this context, the discovery of a double inhibitor may both prevent viral entry and reduce post-infection pulmonary fibrosis. In this study, we use a database of 56 compounds, among which 37 have known experimental affinity with Galectin-3. We carry out virtual screening of this database with respect to Galectin-3 and Spike protein. Several ligands are found to exhibit promising binding affinity and interaction with the Spike protein's N-terminal domain as well as with Galectin-3. This finding strongly suggests that existing Galectin-3 inhibitors possess dual-binding capabilities to disrupt Spike-ACE2 interactions. Herein we identify the most promising inhibitors of Galectin-3 and Spike proteins, of which five emerge as potential dual effective inhibitors. Our preliminary results warrant further in vitro and in vivo testing of these putative inhibitors against SARS-CoV-2 with the hope of being able to halt the spread of the virus in the future.

Noncovalent Interactions with PAMAM and PPI Dendrimers Promote the Cellular Uptake and Photodynamic Activity of Rose Bengal: The Role of the Dendrimer Structure.

Rose bengal is an anionic dye considered as a potential photosensitizer for anticancer photodynamic therapy. The clinical utility of rose bengal is hampered by its short half-life, limited transmembrane transport, aggregation, and self-quenching; consequently, efficient drug carriers that overcome these obstacles are urgently required. In this study, we performed multilevel in vitro and in silico characterization of interactions between rose bengal and cationic poly(amidoamine) (PAMAM) and poly(propyleneimine) (PPI) dendrimers of the third and fourth generation and assessed the ability of the resultant complexes to modulate the photosensitizing properties of the drug. We focused on explaining the molecular basis of this phenomenon and proved that the generation- and structure-dependent binding of the dye by the dendrimers increases the cellular uptake and production of singlet oxygen and intracellular reactive oxygen species, leading to an increase in phototoxicity. We conclude that the application of dendrimer carriers could enable the design of efficient photodynamic therapies based on rose bengal.

Enhancing the activity of platinum-based drugs by improved inhibitors of ERCC1-XPF-mediated DNA repair.

PURPOSE: The ERCC1-XPF 5'-3' DNA endonuclease complex is involved in the nucleotide excision repair pathway and in the DNA inter-strand crosslink repair pathway, two key mechanisms modulating the activity of chemotherapeutic alkylating agents in cancer cells. Inhibitors of the interaction between ERCC1 and XPF can be used to sensitize cancer cells to such drugs. METHODS: We tested recently synthesized new generation inhibitors of this interaction and evaluated their capacity to sensitize cancer cells to the genotoxic activity of agents in synergy studies, as well as their capacity to inhibit the protein-protein interaction in cancer cells using proximity ligation assay. RESULTS: Compound B9 showed the best activity being synergistic with cisplatin and mitomycin C in both colon and lung cancer cells. Also, B9 abolished the interaction between ERCC1 and XPF in cancer cells as shown by proximity ligation assay. Results of different compounds correlated with values from our previously obtained in silico predictions. CONCLUSION: Our results confirm the feasibility of the approach of targeting the protein-protein interaction between ERCC1 and XPF to sensitize cancer cells to alkylating agents, thanks to the improved binding affinity of the newly synthesized compounds.

Prediction of Protein-Protein Interactions Between Alsin DH/PH and Rac1 and Resulting Protein Dynamics.

Alsin is a protein of 1,657 amino acids known for its crucial role in vesicular trafficking in neurons thanks to its ability to interact with two guanosine triphosphatases, Rac1 and Rab5. Evidence suggests that Rac1 can bind Alsin central region, composed by a Dbl Homology (DH) domain followed by a Pleckstrin Homology (PH) domain, leading to Alsin relocalization. However, Alsin three-dimensional structure and its relationship with known biological functions of this protein are still unknown. In this work, a homology model of the Alsin DH/PH domain was developed and studied through molecular dynamics both in the presence and in the absence of its binding partner, Rac1. Due to different conformations of DH domain, the presence of Rac1 seems to stabilize an open state of the protein, while the absence of its binding partner results in closed conformations. Furthermore, Rac1 interaction was able to reduce the fluctuations in the second conserved region of DH motif, which may be involved in the formation of a homodimer. Moreover, the dynamics of DH/PH was described through a Markov State Model to study the pathways linking the open and closed states. In conclusion, this work provided an all-atom model for the DH/PH domain of Alsin protein; moreover, molecular dynamics investigations suggested underlying molecular mechanisms in the signal transduction between Rac1 and Alsin, providing the basis for a deeper understanding of the whole structure-function relationship for Alsin protein.

Cardiovascular Risk Prediction in Ankylosing Spondylitis: From Traditional Scores to Machine Learning Assessment.

INTRODUCTION: The performance of seven cardiovascular (CV) risk algorithms is evaluated in a multicentric cohort of ankylosing spondylitis (AS) patients. Performance and calibration of traditional CV predictors have been compared with the novel paradigm of machine learning (ML). METHODS: A retrospective analysis of prospectively collected data from an AS cohort has been performed. The primary outcome was the first CV event. The discriminatory ability of the algorithms was evaluated using the area under the receiver operating characteristic (ROC) curve (AUC), which is like the concordance-statistic (c-statistic). Three ML techniques were considered to calculate the CV risk: support vector machine (SVM), random forest (RF), and k-nearest neighbor (KNN). RESULTS: Of 133 AS patients enrolled, 18 had a CV event. c-statistic scores of 0.71, 0.61, 0.66, 0.68, 0.66, 0.72, and 0.67 were found, respectively, for SCORE, CUORE, FRS, QRISK2, QRISK3, RRS, and ASSIGN. AUC values for the ML algorithms were: 0.70 for SVM, 0.73 for RF, and 0.64 for KNN. Feature analysis showed that C-reactive protein (CRP) has the highest importance, while SBP and hypertension treatment have lower importance. CONCLUSIONS: All of the evaluated CV risk algorithms exhibit a poor discriminative ability, except for RRS and SCORE, which showed a fair performance. For the first time, we demonstrated that AS patients do not show the traditional ones used by CV scores and that the most important variable is CRP. The present study contributes to a deeper understanding of CV risk in AS, allowing the development of innovative CV risk patient-specific models.

In silico Investigations of the Mode of Action of Novel Colchicine Derivatives Targeting ß-Tubulin Isotypes: A Search for a Selective and Specific ß-III Tubulin Ligand.

The cardinal role of microtubules in cell mitosis makes them interesting drug targets for many pharmacological treatments, including those against cancer. Moreover, different expression patterns between cell types for several tubulin isotypes represent a great opportunity to improve the selectivity and specificity of the employed drugs and to design novel compounds with higher activity only on cells of interest. In this context, tubulin isotype ßIII represents an excellent target for anti-tumoral therapies since it is overexpressed in most cancer cells and correlated with drug resistance. Colchicine is a well-known antimitotic agent, which is able to bind the tubulin dimer and to halt the mitotic process. However, it shows high toxicity also on normal cells and it is not specific for isotype ßIII. In this context, the search for colchicine derivatives is a matter of great importance in cancer research. In this study, homology modeling techniques, molecular docking, and molecular dynamics simulations have been employed to characterize the interaction between 55 new promising colchicine derivatives and tubulin isotype ßIII. These compounds were screened and ranked based on their binding affinity and conformational stability in the colchicine binding site of tubulin ßIII. Results from this study point the attention on an amide of 4-chlorine thiocolchicine. This colchicine-derivative is characterized by a unique mode of interaction with tubulin, compared to all other compounds considered, which is primarily characterized by the involvement of the α-T5 loop, a key player in the colchicine binding site. Information provided by the present study may be particularly important in the rational design of colchicine-derivatives targeting drug resistant cancer phenotypes.

Using the Gibbs Function as a Measure of Human Brain Development Trends from Fetal Stage to Advanced Age.

We propose to use a Gibbs free energy function as a measure of the human brain development. We adopt this approach to the development of the human brain over the human lifespan: from a prenatal stage to advanced age. We used proteomic expression data with the Gibbs free energy to quantify human brain's protein-protein interaction networks. The data, obtained from BioGRID, comprised tissue samples from the 16 main brain areas, at different ages, of 57 post-mortem human brains. We found a consistent functional dependence of the Gibbs free energies on age for most of the areas and both sexes. A significant upward trend in the Gibbs function was found during the fetal stages, which is followed by a sharp drop at birth with a subsequent period of relative stability and a final upward trend toward advanced age. We interpret these data in terms of structure formation followed by its stabilization and eventual deterioration. Furthermore, gender data analysis has uncovered the existence of functional differences, showing male Gibbs function values lower than female at prenatal and neonatal ages, which become higher at ages 8 to 40 and finally converging at late adulthood with the corresponding female Gibbs functions.