Conformational Dynamics and Molecular Characterization of Alsin MORN Monomer and Dimeric Assemblies.

Despite the ubiquity of membrane occupation recognition nexus (MORN) motifs across diverse species in both eukaryotic and prokaryotic organisms, these protein domains remain poorly characterized. Their significance is underscored in the context of the Alsin protein, implicated in the debilitating condition known as infantile-onset ascending hereditary spastic paralysis (IAHSP). Recent investigations have proposed that mutations within the Alsin MORN domain disrupt proper protein assembly, precluding the formation of the requisite tetrameric configuration essential for the protein's inherent biological activity. However, a comprehensive understanding of the relationship between the biological functions of Alsin and its three-dimensional molecular structure is hindered by the lack of available experimental structures. In this study, we employed and compared several protein structure prediction algorithms to identify a three-dimensional structure for the putative MORN of Alsin. Furthermore, inspired by experimental pieces of evidence from previous studies, we employed the developed models to predict and investigate two homo-dimeric assemblies, characterizing their stability. This study's insights into the three-dimensional structure of the Alsin MORN domain and the stability dynamics of its homo-dimeric assemblies suggest an antiparallel linear configuration stabilized by a noncovalent interaction network.

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.

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.

Fludarabine-Specific Molecular Interactions with Maltose-Modified Poly(propyleneimine) Dendrimer Enable Effective Cell Entry of the Active Drug Form: Comparison with Clofarabine.

Fludarabine is an anticancer antimetabolite essential for modern chemotherapy, but its efficacy is limited due to the complex pharmacokinetics. We demonstrated the potential use of maltose-modified poly(propyleneimine) dendrimer as drug delivery agent to improve the efficiency of therapy with fludarabine. In this study, we elaborated a novel synthesis technique for radioactively labeled fludarabine triphosphate to prove for the first time the direct ability of nucleotide-glycodendrimer complex to enter and kill leukemic cells, without the involvement of membrane nucleoside transporters and intracellular kinases. This will potentially allow to bypass the most common drug resistance mechanisms observed in the clinical setting. Further, we applied surface plasmon resonance and molecular modeling to elucidate the properties of the drug-dendrimer complexes. We showed that clofarabine, a more toxic nucleoside analogue drug, is characterized by significantly different molecular interactions with poly(propyleneimine) dendrimers than fludarabine, leading to different cellular outcomes (decreased rather than increased treatment efficiency). The most probable mechanistic explanation of uniquely dendrimer-enhanced fludarabine toxicity points to a crucial role of both an alternative cellular uptake pathway and the avoidance of intracellular phosphorylation of nucleoside drug form.

The synergistic effect of chlorotoxin-mApoE in boosting drug-loaded liposomes across the BBB.

We designed liposomes dually functionalized with ApoE-derived peptide (mApoE) and chlorotoxin (ClTx) to improve their blood-brain barrier (BBB) crossing. Our results demonstrated the synergistic activity of ClTx-mApoE in boosting doxorubicin-loaded liposomes across the BBB, keeping the anti-tumour activity of the drug loaded: mApoE acts promoting cellular uptake, while ClTx promotes exocytosis of liposomes.

Destabilizing the AXH Tetramer by Mutations: Mechanisms and Potential Antiaggregation Strategies.

The AXH domain of protein Ataxin 1 is thought to play a key role in the misfolding and aggregation pathway responsible for Spinocerebellar ataxia 1. For this reason, a molecular level understanding of AXH oligomerization pathway is crucial to elucidate the aggregation mechanism, which is thought to trigger the disease. This study employs classical and enhanced molecular dynamics to identify the structural and energetic basis of AXH tetramer stability. Results of this work elucidate molecular mechanisms behind the destabilizing effect of protein mutations, which consequently affect the AXH tetramer assembly. Moreover, results of the study draw attention for the first time, to our knowledge, to the R638 protein residue, which is shown to play a key role in AXH tetramer stability. Therefore, R638 might be also implicated in the AXH oligomerization pathway and stands out as a target for future experimental studies focused on self-association mechanisms and fibril formation of full-length ATX1.

Conformational Dynamics and Stability of U-Shaped and S-Shaped Amyloid ß Assemblies.

Alzheimer's disease is the most fatal neurodegenerative disorder characterized by the aggregation and deposition of Amyloid ß (Aß) oligomers in the brain of patients. Two principal variants of Aß exist in humans: Aß1-40 and Aß1-42. The former is the most abundant in the plaques, while the latter is the most toxic species and forms fibrils more rapidly. Interestingly, fibrils of Aß1-40 peptides can only assume U-shaped conformations while Aß1-42 can also arrange as S-shaped three-stranded chains, as recently discovered. As alterations in protein conformational arrangement correlate with cell toxicity and speed of disease progression, it is important to characterize, at molecular level, the conformational dynamics of amyloid fibrils. In this work, Replica Exchange Molecular Dynamics simulations were carried out to compare the conformational dynamics of U-shaped and S-shaped Aß17-42 small fibrils. Our computational results provide support for the stability of the recently proposed S-shaped model due to the maximized interactions involving the C-terminal residues. On the other hand, the U-shaped motif is characterized by significant distortions resulting in a more disordered assembly. Outcomes of our work suggest that the molecular architecture of the protein aggregates might play a pivotal role in formation and conformational stability of the resulting fibrils.

Self-Assembled Ligands Targeting TLR7: A Molecular Level Investigation.

Toll-like receptors (TLRs) are pattern recognition transmembrane proteins that play an important role in innate immunity. In particular, TLR7 plays a role in detecting nucleic acids derived from viruses and bacteria. The huge number of pathologies in which TLR7 is involved has led to an increasing interest in developing new compounds targeting this protein. Several conjugation strategies were proposed for TLR7 agonists to increase the potency while maintaining a low toxicity. In this work, we focus the attention on two promising classes of TLR7 compounds derived from the same pharmacophore conjugated with phospholipid and polyethylene glycol (PEG). A multidisciplinary investigation has been carried out by molecular dynamics (MD), dynamic light scattering (DLS), electron paramagnetic resonance (EPR), and cytotoxicity assessment. DLS and MD indicated how only the phospholipid conjugation provides the compound abilities to self-assemble in an orderly fashion with a maximal pharmacophore exposition to the solvent. Further EPR and cytotoxicity experiments highlighted that phospholipid compounds organize in stable aggregates and well interact with TLR7, whereas PEG conjugation was characterized by poorly stable aggregates at the cells surface. The methodological framework proposed in this study may be used to investigate, at a molecular level, the interactions generally occurring between aggregated ligands, to be used as drugs, and protein receptors.

Josephin Domain Structural Conformations Explored by Metadynamics in Essential Coordinates.

The Josephin Domain (JD), i.e. the N-terminal domain of Ataxin 3 (At3) protein, is an interesting example of competition between physiological function and aggregation risk. In fact, the fibrillogenesis of Ataxin 3, responsible for the spinocerebbellar ataxia 3, is strictly related to the JD thermodynamic stability. Whereas recent NMR studies have demonstrated that different JD conformations exist, the likelihood of JD achievable conformational states in solution is still an open issue. Marked differences in the available NMR models are located in the hairpin region, supporting the idea that JD has a flexible hairpin in dynamic equilibrium between open and closed states. In this work we have carried out an investigation on the JD conformational arrangement by means of both classical molecular dynamics (MD) and Metadynamics employing essential coordinates as collective variables. We provide a representation of the free energy landscape characterizing the transition pathway from a JD open-like structure to a closed-like conformation. Findings of our in silico study strongly point to the closed-like conformation as the most likely for a Josephin Domain in water.

Explaining the Microtubule Energy Balance: Contributions Due to Dipole Moments, Charges, van der Waals and Solvation Energy.

Microtubules are the main components of mitotic spindles, and are the pillars of the cellular cytoskeleton. They perform most of their cellular functions by virtue of their unique dynamic instability processes which alternate between polymerization and depolymerization phases. This in turn is driven by a precise balance between attraction and repulsion forces between the constituents of microtubules (MTs)-tubulin dimers. Therefore, it is critically important to know what contributions result in a balance of the interaction energy among tubulin dimers that make up microtubules and what interactions may tip this balance toward or away from a stable polymerized state of tubulin. In this paper, we calculate the dipole-dipole interaction energy between tubulin dimers in a microtubule as part of the various contributions to the energy balance. We also compare the remaining contributions to the interaction energies between tubulin dimers and establish a balance between stabilizing and destabilizing components, including the van der Waals, electrostatic, and solvent-accessible surface area energies. The energy balance shows that the GTP-capped tip of the seam at the plus end of microtubules is stabilized only by - 9 kcal/mol, which can be completely reversed by the hydrolysis of a single GTP molecule, which releases + 14 kcal/mol and destabilizes the seam by an excess of + 5 kcal/mol. This triggers the breakdown of microtubules and initiates a disassembly phase which is aptly called a catastrophe.

Conformational fluctuations of the AXH monomer of Ataxin-1.

In this paper, we report the results of molecular dynamics simulations of AXH monomer of Ataxin-1. The AXH domain plays a crucial role in Ataxin-1 aggregation, which accompanies the initiation and progression of Spinocerebellar ataxia type 1. Our simulations involving both classical and replica exchange molecular dynamics, followed by principal component analysis of the trajectories obtained, reveal substantial conformational fluctuations of the protein structure, especially in the N-terminal region. We show that these fluctuations can be generated by thermal noise since the free energy barriers between conformations are small enough for thermally stimulated transitions. In agreement with the previous experimental findings, our results can be considered as a basis for a future design of ataxin aggregation inhibitors that will require several key conformations identified in the present study as molecular targets for ligand binding.

Characterization of the AXH domain of Ataxin-1 using enhanced sampling and functional mode analysis.

Ataxin-1 is the protein responsible for the Spinocerebellar ataxia type 1, an incurable neurodegenerative disease caused by polyglutamine expansion. The AXH domain plays a pivotal role in physiological functions of Ataxin-1. In Spinocerebellar ataxia 1, the AXH domain is involved in the misfolding and aggregation pathway. Here molecular modeling is applied to investigate the protein-protein interactions contributing to the AXH dimer stability. Particular attention is focused on: (i) the characterization of AXH monomer-monomer interface; (ii) the molecular description of the AXH monomer-monomer interaction dynamics. Technically, an approach based on functional mode analysis, here applied to replica exchange molecular dynamics trajectories, was employed. The findings of this study are consistent with previous experimental results and elucidate the pivotal role of the I580 residue in mediating the AXH monomer-monomer interaction dynamics.

Structure Based Modeling of Small Molecules Binding to the TLR7 by Atomistic Level Simulations.

Toll-Like Receptors (TLR) are a large family of proteins involved in the immune system response. Both the activation and the inhibition of these receptors can have positive effects on several diseases, including viral pathologies and cancer, therefore prompting the development of new compounds. In order to provide new indications for the design of Toll-Like Receptor 7 (TLR7)-targeting drugs, the mechanism of interaction between the TLR7 and two important classes of agonists (imidazoquinoline and adenine derivatives) was investigated through docking and Molecular Dynamics simulations. To perform the computational analysis, a new model for the dimeric form of the receptors was necessary and therefore created. Qualitative and quantitative differences between agonists and inactive compounds were determined. The in silico results were compared with previous experimental observations and employed to define the ligand binding mechanism of TLR7.

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.

Multiscale Computational Analysis of the Effect of Taxol on Microtubule Mechanics.

Microtubules (MTs) are widely recognized as targets for cancer therapies. They are directly related to unique mechanical properties, closely dependent on MT architecture and tubulin molecular features. Taxol is known to affect tubulin interactions resulting in the stabilization of the MT lattice, and thus the hierarchical organization stability, mechanics, and function. A deeper understanding of the molecular mechanisms through which taxol modulates intertubulin interactions in the MT lattice, and consequently, its stability and mechanical response is crucial to characterize how MT properties are regulated by environmental factors, such as interacting ligands. In this study, a computational analysis of the effect of taxol on the MT was performed at different scales, combining molecular dynamics simulation, dynamical network analysis, and elastic network modeling. The results show that the taxol-induced conformational differences at the M-loop region increase the stability of the lateral interactions and the amount of surface in contact between laterally coupled tubulins. Moreover, the conformational rearrangements in the taxane binding site result in a different structural communication pattern. Finally, the different conformation of the tubulin heterodimers and the stabilized lateral interactions resulted in a tendency toward higher deformation of the vibrating MT in the presence of taxol. Overall, this work provides additional insights into taxol-induced stabilization and relates the conformational changes at the tubulin level to the MT mechanics. Besides providing useful insights into taxol effect on MT mechanics, a methodological framework that could be used to characterize the effects of other MT stabilizing agents is presented.

Predicting multiple taste sensations with a multiobjective machine learning method.

Taste perception plays a pivotal role in guiding nutrient intake and aiding in the avoidance of potentially harmful substances through five basic tastes - sweet, bitter, umami, salty, and sour. Taste perception originates from molecular interactions in the oral cavity between taste receptors and chemical tastants. Hence, the recognition of taste receptors and the subsequent perception of taste heavily rely on the physicochemical properties of food ingredients. In recent years, several advances have been made towards the development of machine learning-based algorithms to classify chemical compounds' tastes using their molecular structures. Despite the great efforts, there remains significant room for improvement in developing multi-class models to predict the entire spectrum of basic tastes. Here, we present a multi-class predictor aimed at distinguishing bitter, sweet, and umami, from other taste sensations. The development of a multi-class taste predictor paves the way for a comprehensive understanding of the chemical attributes associated with each fundamental taste. It also opens the potential for integration into the evolving realm of multi-sensory perception, which encompasses visual, tactile, and olfactory sensations to holistically characterize flavour perception. This concept holds promise for introducing innovative methodologies in the rational design of foods, including pre-determining specific tastes and engineering complementary diets to augment traditional pharmacological treatments.

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.

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.

Multiscale modeling of cellular actin filaments: from atomistic molecular to coarse-grained dynamics.

In this article, we present a computational multiscale model for the characterization of subcellular proteins. The model is encoded inside a simulation tool that builds coarse-grained (CG) force fields from atomistic simulations. Equilibrium molecular dynamics simulations on an all-atom model of the actin filament are performed. Then, using the statistical distribution of the distances between pairs of selected groups of atoms at the output of the MD simulations, the force field is parameterized using the Boltzmann inversion approach. This CG force field is further used to characterize the dynamics of the protein via Brownian dynamics simulations. This combination of methods into a single computational tool flow enables the simulation of actin filaments with length up to 400 nm, extending the time and length scales compared to state-of-the-art approaches. Moreover, the proposed multiscale modeling approach allows to investigate the relationship between atomistic structure and changes on the overall dynamics and mechanics of the filament and can be easily (i) extended to the characterization of other subcellular structures and (ii) used to investigate the cellular effects of molecular alterations due to pathological conditions.