TAT-RasGAP317-326 kills cells by targeting inner-leaflet-enriched phospholipids.

TAT-RasGAP317-326 is a cell-penetrating peptide-based construct with anticancer and antimicrobial activities. This peptide kills a subset of cancer cells in a manner that does not involve known programmed cell death pathways. Here we have elucidated the mode of action allowing TAT-RasGAP317-326 to kill cells. This peptide binds and disrupts artificial membranes containing lipids typically enriched in the inner leaflet of the plasma membrane, such as phosphatidylinositol-bisphosphate (PIP2) and phosphatidylserine (PS). Decreasing the amounts of PIP2 in cells renders them more resistant to TAT-RasGAP317-326, while reducing the ability of cells to repair their plasma membrane makes them more sensitive to the peptide. The W317A TAT-RasGAP317-326 point mutant, known to have impaired killing activities, has reduced abilities to bind and permeabilize PIP2- and PS-containing membranes and to translocate through biomembranes, presumably because of a higher propensity to adopt an α-helical state. This work shows that TAT-RasGAP317-326 kills cells via a form of necrosis that relies on the physical disruption of the plasma membrane once the peptide targets specific phospholipids found on the cytosolic side of the plasma membrane.

A calibrated model with a single-generator simulating polysomnographically recorded periodic leg movements.

The aim of this study was to assess, with numerical simulations, if the complex mechanism of two (or more) interacting spinal/supraspinal structures generating periodic leg movements can be modelled with a single-generator approach. For this, we have developed the first phenomenological model to generate periodic leg movements in-silico. We defined the onset of a movement in one leg as the firing of a neuron integrating excitatory and inhibitory inputs from the central nervous system, while the duration of the movement was defined in accordance to statistical evidence. For this study, polysomnographic leg movement data from 32 subjects without periodic leg movements and 65 subjects with periodic leg movements were used. The proportion of single-leg and double-leg inputs, as well as their strength and frequency, were calibrated on the without periodic leg movements dataset. For periodic leg movements subjects, we added a periodic excitatory input common to both legs, and the distributions of the generator period and intensity were fitted to their dataset. Besides the many simplifying assumptions - the strongest being the stationarity of the generator processes during sleep - the model-simulated data did not differ significantly, to a large extent, from the real polysomnographic data. This represents convincing preliminary support for the validity of our single-generator model for periodic leg movements. Future model extensions will pursue the ambitious project of a supportive diagnostic and therapeutic tool, helping the specialist with realistic forecasting, and with cross-correlations and clustering with other patient meta-data.

Extended diffusion theory: Recovering dynamics from biased/accelerated molecular simulations.

Dynamical properties are of great importance in determining the behavior of synthetic and natural molecules, but capturing them by computational methods is a nontrivial task. Very often the time scales of the relevant phenomena are far beyond the typical time windows accessible by classical Molecular Dynamics (MD) simulations, currently limited to the order of microseconds on standard laboratory workstations. On the other hand, biased and accelerated simulations allow for fast and thorough exploration of the molecular conformational space, but they lose the dynamic information. The problem of recovering dynamics from biased/accelerated simulations is a very active field of research, but no totally robust/reliable solutions have been given yet. In this paper it is shown how the Smoluchowski equation, in the framework of Diffusion Theory (DT), can be used to bridge this gap, and dynamical properties, in the form of time correlation functions (TCFs), can be extracted also from such kind of simulations. DT is first extended (EDT) to express the mobility tensors entering the Smoluchowski operator in terms of a recently introduced unified and regularized Rotne-Prager-Yamakawa approximation, [P. J. Zuk, E. Wajnryb, K. A. Mizerski, P. Szymczak, J. Fluid. Mech. 2014, 741, R5, 1-13] also involving mixed rotation-translation contributions, and rotation-rotation terms beside the classical translation-translation ones, so far used in DT. Then, the method is applied to recover the dynamics of a nontrivial example of a peptide in explicit water from the first 200 ns of a Replica Exchange Molecular Dynamics simulation, which is a popular computational method that destroys the long time dynamics. EDT dynamics were found to favorably compare against those coming from a standard MD simulation of the same system, requiring a time window of 30 µs to converge. This result shows that EDT is a tool of practical value to recover the long time dynamics of systems in diffusive regimes from biased/accelerated simulations, to be exploited in those cases when direct evaluation by standard MD is unfeasible.

The Impact of Natural Compounds on S-Shaped Aß42 Fibril: From Molecular Docking to Biophysical Characterization.

The pursuit for effective strategies inhibiting the amyloidogenic process in neurodegenerative disorders, such as Alzheimer's disease (AD), remains one of the main unsolved issues, and only a few drugs have demonstrated to delay the degeneration of the cognitive system. Moreover, most therapies induce severe side effects and are not effective at all stages of the illness. The need to find novel and reliable drugs appears therefore of primary importance. In this context, natural compounds have shown interesting beneficial effects on the onset and progression of neurodegenerative diseases, exhibiting a great inhibitory activity on the formation of amyloid aggregates and proving to be effective in many preclinical and clinical studies. However, their inhibitory mechanism is still unclear. In this work, ensemble docking and molecular dynamics simulations on S-shaped Aß42 fibrils have been carried out to evaluate the influence of several natural compounds on amyloid conformational behaviour. A deep understanding of the interaction mechanisms between natural compounds and Aß aggregates may play a key role to pave the way for design, discovery and optimization strategies toward an efficient destabilization of toxic amyloid assemblies.

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.

Leveraging Machine Learning-Guided Molecular Simulations Coupled with Experimental Data to Decipher Membrane Binding Mechanisms of Aminosterols.

Understanding the molecular mechanisms of the interactions between specific compounds and cellular membranes is essential for numerous biotechnological applications, including targeted drug delivery, elucidation of the drug mechanism of action, pathogen identification, and novel antibiotic development. However, estimation of the free energy landscape associated with solute binding to realistic biological systems is still a challenging task. In this work, we leverage the Time-lagged Independent Component Analysis (TICA) in combination with neural networks (NN) through the Deep-TICA approach for determining the free energy associated with the membrane insertion processes of two natural aminosterol compounds, trodusquemine (TRO), and squalamine (SQ). These compounds are particularly noteworthy because they interact with the outer layer of neuron membranes, protecting them from the toxic action of misfolded proteins involved in neurodegenerative disorders, in both their monomeric and oligomeric forms. We demonstrate how this strategy could be used to generate an effective collective variable for describing solute absorption in the membrane and for estimating free energy landscape of translocation via on-the-fly probability enhanced sampling (OPES) method. In this context, the computational protocol allowed an exhaustive characterization of the aminosterol entry pathway into a neuron-like lipid bilayer. Furthermore, it provided accurate prediction of membrane binding affinities, in close agreement with the experimental binding data obtained by using fluorescently labeled aminosterols and large unilamellar vesicles (LUVs). The findings contribute significantly to our understanding of aminosterol entry pathways and aminosterol-lipid membrane interactions. Finally, the computational methods deployed in this study further demonstrate considerable potential for investigating membrane binding processes.

Clickable polyethyleneimine incorporated into triblock copolymeric micelles as an efficient platform in the delivery of siRNA to NSCLC cells.

The efficacy of transfection vectors to cross the endosomal membrane into the cytosol is a central aspect in the development of nucleic acid-based therapeutics. The challenge remains the same: Delivery, Delivery, Delivery. Despite a rational and appropriate construct of triblock polymeric micelles, which could serve as an ideal platform for the co-delivery of siRNAs and hydrophobic anticancer drugs, we show here its inability to properly convey oligonucleotides to their final destination. In order to overcome biological barriers, a linear PEI comprising two orthogonal groups was synthesized, holding an appropriate balance between safety and efficacy. Micellar carriers were then formulated with this polymer to enhance endosomal siRNA release. This chemical technology also addresses the two major challenges to consider when developing novel micellar products for siRNA delivery, namely cytotoxicity of polycations and endosomal escape. Herein, we demonstrate successful release of siRNA using a polymer tailoring strategy combined with a relevant in vitro approach, considering STAT3 as a promising target in the treatment of non-small cell lung cancer (NSCLC).

Biophysical insights into cancer transformation and treatment.

Biological systems are hierarchically self-organized complex structures characterized by nonlinear interactions. Biochemical energy is transformed into work of physical forces required for various biological functions. We postulate that energy transduction depends on endogenous electrodynamic fields generated by microtubules. Microtubules and mitochondria colocalize in cells with microtubules providing tracks for mitochondrial movement. Besides energy transformation, mitochondria form a spatially distributed proton charge layer and a resultant strong static electric field, which causes water ordering in the surrounding cytosol. These effects create conditions for generation of coherent electrodynamic field. The metabolic energy transduction pathways are strongly affected in cancers. Mitochondrial dysfunction in cancer cells (Warburg effect) or in fibroblasts associated with cancer cells (reverse Warburg effect) results in decreased or increased power of the generated electromagnetic field, respectively, and shifted and rebuilt frequency spectra. Disturbed electrodynamic interaction forces between cancer and healthy cells may favor local invasion and metastasis. A therapeutic strategy of targeting dysfunctional mitochondria for restoration of their physiological functions makes it possible to switch on the natural apoptotic pathway blocked in cancer transformed cells. Experience with dichloroacetate in cancer treatment and reestablishment of the healthy state may help in the development of novel effective drugs aimed at the mitochondrial function.

In silico investigation of cytochrome bc1 molecular inhibition mechanism against Trypanosoma cruzi.

Chagas' disease is a neglected tropical disease caused by the kinetoplastid protozoan Trypanosoma cruzi. The only therapies are the nitroheterocyclic chemicals nifurtimox and benznidazole that cause various adverse effects. The need to create safe and effective medications to improve medical care remains critical. The lack of verified T. cruzi therapeutic targets hinders medication research for Chagas' disease. In this respect, cytochrome bc1 has been identified as a promising therapeutic target candidate for antibacterial medicines of medical and agricultural interest. Cytochrome bc1 belongs to the mitochondrial electron transport chain and transfers electrons from ubiquinol to cytochrome c1 by the action of two catalytic sites named Qi and Qo. The two binding sites are highly selective, and specific inhibitors exist for each site. Recent studies identified the Qi site of the cytochrome bc1 as a promising drug target against T. cruzi. However, a lack of knowledge of the drug mechanism of action unfortunately hinders the development of new therapies. In this context, knowing the cause of binding site selectivity and the mechanism of action of inhibitors and substrates is crucial for drug discovery and optimization processes. In this paper, we provide a detailed computational investigation of the Qi site of T. cruzi cytochrome b to shed light on the molecular mechanism of action of known inhibitors and substrates. Our study emphasizes the action of inhibitors at the Qi site on a highly unstructured portion of cytochrome b that could be related to the biological function of the electron transport chain complex.

Influence of dexamethasone on the interaction between glucocorticoid receptor and SOX9: A molecular dynamics study.

The glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. GR transcriptional activity is modulated by a series of ligands and coenzymes, where a ligand can act as an agonist or antagonist. GR agonists, such as the glucocorticoids dexamethasone (DEX) and prednisolone, are widely prescribed to patients with inflammatory and autoimmune diseases. DEX is also used to induce osteogenic differentiation in vitro. Recently, it has been highlighted that DEX induces changes in the osteogenic differentiation of human mesenchymal stromal cells by downregulating the transcription factor SRY-box transcription factor 9 (SOX9) and upregulating the peroxisome proliferator-activated receptor γ (PPARG). SOX9 is fundamental in the control of chondrogenesis, but also in osteogenesis by acting as a dominant-negative of RUNX2. Many processes remain to be clarified during cell fate determination, such as the interplay between the key transcription factors. The main objective pursued by this work is to shed light on the interaction between GR and SOX9 in the presence and absence of DEX at an atomic level of resolution using molecular dynamics simulations. The outcome of this research could help the understanding of possible molecular interactions between GR and SOX9 and their role in the determination of cell fate. The results highlight the key residues at the interface between GR and SOX9 involved in the complexation process and shed light on the mechanism through which DEX modulates GR-SOX9 binding and exerts its biological activity.

Reorganization of the outer layer of a model of the plasma membrane induced by a neuroprotective aminosterol.

Trodusquemine is an amphipathic aminosterol that has recently shown therapeutic benefit in neurodegenerative diseases altering the binding of misfolded proteins to the cell membrane. To unravel the underlying mechanism, we studied the interactions between Trodusquemine (TRO) and lipid monolayers simulating the outer layer of the plasma membrane. We selected two different compositions of dioleoylphosphatidylcholine (DOPC), sphingomyelin (SM), cholesterol (Chol) and monosialotetrahexosylganglioside (GM1) lipid mixture mimicking either a lipid-raft containing membrane (Ld+So phases) or a single-phase disordered membrane (Ld phase). Surface pressure-area isotherms and surface compressional modulus-area combined with Brewster Angle Microscopy (BAM) provided the thermodynamic and morphological information on the lipid monolayer in the presence of increasing amounts of TRO in the monolayer. Experiments revealed that TRO forms stable spreading monolayers at the buffer-air interface where it undergoes multiple reversible phase transitions to bi- and tri-layers at the interface. When TRO was spread at the interface with the lipid mixtures, we found that it distributes in the lipid monolayer for both the selected lipid compositions, but a maximum TRO uptake in the rafts-containing monolayer was observed for a Lipid/TRO molar ratio equal to 3:2. Statistical analysis of BAM images revealed that TRO induces a decrease in the size of the condensed domains, an increase in their number and in the thickness mismatch between the Ld and So phase. Experiments and MD simulations converge to indicate that TRO adsorbs preferentially at the border of the So domains. Removal of GM1 from the lipid Ld+So mixture resulted in an even greater TRO-mediated reduction of the size of the So domains suggesting that the presence of GM1 hinders the localization of TRO at the So domains boundaries. Taken together these observations suggest that Trodusquemine influences the organization of lipid rafts within the neuronal membrane in a dose-dependent manner whereas it evenly distributes in disordered expanded phases of the membrane model.