Molecular and Coarse-Grained Modeling to Characterize and Optimize Dendrimer-Based Nanocarriers for Short Interfering RNA Delivery

Dendrimer nanocarriers are unique hyper-branched polymers with biomolecule-like properties, representing a promising prospect as a nucleic acid delivery system. The design of effective dendrimer-based gene carriers requires considering several parameters, such as carrier morphology, size, molecular weight, surface chemistry, and flexibility/rigidity. In detail, the rational design of the dendrimer surface chemistry has been ascertained to play a crucial role on the efficiency of interaction with nucleic acids. Within this framework, advances in the field of organic chemistry have allowed us to design dendrimers with even small difference in the chemical structure of their surface terminals. In this study, we have selected two different cationic phosphorus dendrimers of generation 3 functionalized, respectively, with pyrrolidinium (DP) and morpholinium (DM) surface groups, which have demonstrated promising potential for short interfering RNA (siRNA) delivery. Despite DP and DM differing only for one atom in their chemical structure, in vitro and in vivo experiments have highlighted several differences between them in terms of siRNA complexation properties. In this context, we have employed coarse-grained molecular dynamics simulation techniques to shed light on the supramolecular characteristics of dendrimer-siRNA complexation, the so-called dendriplex formations. Our data provide important information on self-assembly dynamics driven by surface chemistry and competition mechanisms

Aminoacid substitutions in the glycine zipper affect the conformational stability of amyloid beta fibrils

The aggregation of amyloid-beta peptides is associated with the pathogenesis of Alzheimer’s disease. The hydrophobic core of the amyloid beta sequence contains a GxxxG repeated motif, called glycine zipper, which involves crucial residues for assuring stability and promoting the process of fibril formation. Mutations in this motif lead to a completely different oligomerization pathway and rate of fibril formation. In this work, we have tested G33L and G37L residue substitutions by molecular dynamics simulations. We found that both protein mutations may lead to remarkable changes in the fibril conformational stability. Results suggest the disruption of the glycine zipper as a possible strategy to reduce the aggregation propensity of amyloid beta peptides. On the basis of our data, further investigations may consider this key region as a binding site to design/discover novel effective inhibitors. Communicated by Ramaswamy H. Sarma

Stem cells and physical energies: can we really drive stem cell fate?

Adult stem cells are undifferentiated elements able to self-renew or differentiate to maintain tissue integrity. Within this context, stem cells are able to divide in a symmetric fashion, feature characterising all the somatic cells, or in an asymmetric way, which leads daughter cells to different fates. It is worth highlighting that cell polarity have a critical role in regulating stem cell asymmetric division and the proper control of cell division depends on different proteins involved in cell development, differentiation and maintenance of tissue homeostasis. Moreover, the interaction between cells and the extracellular matrix are crucial in influencing cell behavior, included in terms of mechanical properties as cytoskeleton plasticity and remodelling, and membrane tension. Finally, the activation of specific transcriptional program and epigenetic modifications contributes to cell fate determination, through modulation of cellular signalling cascades. It is well known that physical and mechanical stimuli are able to influence biological systems, and in this context, the effects of electromagnetic fields (EMFs) have already shown a considerable role, even though there is a lack of knowledge and much remains to be done around this topic. In this review, we summarize the historical background of EMFs applications and the main molecular mechanism involved in cellular remodelling, with particular attention to cytoskeleton elasticity and cell polarity, required for driving stem cell behavior

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

Free energy landscape of siRNA-polycation complexation: Elucidating the effect of molecular geometry, polymer flexibility, and charge neutralization

The success of medical threatments with DNA and silencing interference RNA is strongly related to the design of efficient delivery technologies. Cationic polymers represent an attractive strategy to serve as nucleic-acid carriers with the envisioned advantages of efficient complexation, low cost, ease of production, well-defined size, and low polydispersity index. However, the balance between efficacy and toxicity (safety) of these polymers is a challenge and in need of improvement. With the aim of designing more effective polycationic-based gene carriers, many parameters such as carrier morphology, size, molecular weight, surface chemistry, and flexibility/rigidity ratio need to be taken into consideration. In the present work, the binding mechanism of three cationic polymers (polyarginine, polylysine and polyethyleneimine) to a model siRNA target is computationally investigated at the atomistic level. In order to better understand the polycationic carrier-siRNA interactions, replica exchange molecular dynamic simulations were carried out to provide an exhaustive exploration of all the possible binding sites, taking fully into account the siRNA flexibility together with the presence of explicit solvent and ions. Moreover, well-tempered metadynamics simulations were employed to elucidate how molecular geometry, polycation flexibility, and charge neutralization affect the siRNA-polycations free energy landscape in term of low-energy binding modes and unbinding free energy barriers. Significant differences among polymer binding modes have been detected, revealing the advantageous binding properties of polyarginine and polylysine compared to polyethyleneimine

The role of structural polymorphism in driving the mechanical performance of the alzheimer's beta amyloid fibrils

Alzheimer's Disease (AD) is related with the abnormal aggregation of amyloid β-peptides Aβ1-40 and Aβ1-42, the latter having a polymorphic character which gives rise to U- or S-shaped fibrils. Elucidating the role played by the nanoscale-material architecture on the amyloid fibril stability is a crucial breakthrough to better understand the pathological nature of amyloid structures and to support the rational design of bio-inspired materials. The computational study here presented highlights the superior mechanical behavior of the S-architecture, characterized by a Young's modulus markedly higher than the U-shaped architecture. The S-architecture showed a higher mechanical resistance to the enforced deformation along the fibril axis, consequence of a better interchain hydrogen bonds' distribution. In conclusion, this study, focusing the attention on the pivotal multiscale relationship between molecular phenomena and material properties, suggests the S-shaped Aβ1-42 species as a target of election in computational screen/design/optimization of effective aggregation modulators

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

Thermodynamic and kinetic stability of the Josephin Domain closed arrangement: evidences from replica exchange molecular dynamics

Abstract Background Molecular phenomena driving pathological aggregation in neurodegenerative diseases are not completely understood yet. Peculiar is the case of Spinocerebellar Ataxia 3 (SCA3) where the conformational properties of the AT-3 N-terminal region, also called Josephin Domain (JD), play a key role in the first step of aggregation, having the JD an amyloidogenic propensity itself. For this reason, unraveling the intimate relationship between JD structural features and aggregation tendency may lead to a step forward in understanding the pathology and rationally design a cure. In this connection, computational modeling has demonstrated to be helpful in exploring the protein molecular dynamics and mechanism of action. Results Conformational dynamics of the JD is here finely investigated by replica exchange molecular dynamics simulations able to sample the microsecond time scale and to provide both a thermodynamic and kinetic description of the protein conformational changes. Accessible structural conformations of the JD have been identified in: open, intermediate and closed like arrangement. Data indicated the closed JD arrangement as the most likely protein arrangement. The protein transition from closed toward intermediate/open states was characterized by a rate constant higher than 700\ua0ns. This result also explains the inability of classical molecular dynamics to explore transitions from closed to open JD configuration on a time scale of hundreds of nanoseconds. Conclusion This work provides the first kinetic estimation of the JD transition pathway from open-like to closed-like arrangement and vice-versa, indicating the closed-like arrangement as the most likely configuration for a JD in water environment. More widely, the importance of our results is also underscored considering that the ability to provide a kinetic description of the protein conformational changes is a scientific challenge for both experimental and theoretical approaches to date. Reviewers This article was reviewed by Oliviero Carugo, Bojan Zagrovic

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