Effect of Arginine-Rich Peptide Length on the Structure and Binding Strength of siRNA–Peptide Complexes

Neurobiochemical characteristics of arginine-rich peptides explain their potential therapeutic efficacy in neurodegenerative diseases

Neurodegenerative diseases, including Alzheimer̕ s disease (AD), Parkinson̕ s disease (PD), Huntington̕ s disease (HD), and Amyotrophic Lateral Sclerosis (ALS) require special attention to find new potential treatment methods. This review aims to summarize the current knowledge of the relationship between the biochemical properties of arginine-rich peptides (ARPs) and their neuroprotective effects to deal with the harmful effects of risk factors. It seems that ARPs have portrayed a promising and fantastic landscape for treating neurodegeneration-associated disorders. With multimodal mechanisms of action, ARPs play various unprecedented roles, including as the novel delivery platforms for entering the central nervous system (CNS), the potent antagonists for calcium influx, the invader molecules for targeting mitochondria, and the protein stabilizers. Interestingly, these peptides inhibit the proteolytic enzymes and block protein aggregation to induce pro-survival signaling pathways. ARPs also serve as the scavengers of toxic molecules and the reducers of oxidative stress agents. They also have anti-inflammatory, antimicrobial, and anti-cancer properties. Moreover, by providing an efficient nucleic acid delivery system, ARPs can play an essential role in developing various fields, including gene vaccines, gene therapy, gene editing, and imaging. ARP agents and ARP/cargo therapeutics can be raised as an emergent class of neurotherapeutics for neurodegeneration. Part of the aim of this review is to present recent advances in treating neurodegenerative diseases using ARPs as an emerging and powerful therapeutic tool. The applications and progress of ARPs-based nucleic acid delivery systems have also been discussed to highlight their usefulness as a broad-acting class of drugs.

In silico prediction of the in vitro behavior of polymeric gene delivery vectors

Non-viral gene delivery vectors have increasingly come under the spotlight, but their performaces are still far from being satisfactory. Therefore, there is an urgent need for forecasting tools and screening methods to enable the development of ever more effective transfectants. Here, coarse-grained (CG) models of gold standard transfectant poly(ethylene imine)s (PEIs) have been profitably used to investigate and highlight the effect of experimentally-relevant parameters, namely molecular weight (2 vs. 10 kDa) and topologies (linear vs. branched), protonation state, and ammine-to-phosphate ratios (N/Ps), on the complexation and the gene silencing efficiency of siRNA molecules. The results from the in vitro screening of cationic polymers and conditions were used to validate the in silico platform that we developed, such that the hits which came out of the CG models were of high practical relevance. We show that our in silico platform enables to foresee the most suitable conditions for the complexation of relevant siRNA-polycation assemblies, thereby providing a reliable predictive tool to test bench transfectants in silico, and foster the design and development of gene delivery vectors.

Coarse-grained Simulations of the Impact of Chain Length and Stiffness on the Formation and Aggregation of Polyelectrolyte Complexes

Polyelectrolyte complexes formed from nucleic acids and synthetic polycations have been studied because of their potential in gene delivery. Coarse-grained molecular dynamics simulations are performed to examine the impact of chain length and polyanion stiffness on polyplex formation and aggregation. Polyplexes containing single polyanion chain fall into three structural regimes depending on polyanion stiffness: flexible polyanions form collapsed complexes, semiflexible polyanions form various morphologies including toroids and hairpins, and stiff polyanions form rod-like structures. Polyplex size generally decreases as polycation length increases. Aggregation (i.e., formation of complexes containing multiple polyanions) is observed in some simulations containing multiple polyanions and an excess of short polycations. Aggregation is observed to only occur for semiflexible and stiff polyanions and is promoted by shorter polycation lengths. Simulations of short, stiff polyanions condensed by long polycations are used as a model for siRNA gene delivery complexes. These simulations show multiple polyanions are spaced out along the polycation with polyanion-polyanion interactions, usually limited to overlapping chain ends. These structures differ from aggregates of longer polyanions in which the polyanions are packed together in parallel, forming bundles.

Binding Dynamics of siRNA with Selected Lipopeptides: A Computer-Aided Study of the Effect of Lipopeptides' Functional Groups and Stereoisomerism

The engineering issues pertaining to nanoparticle systems toward targeted gene therapies have not been fully probed. Recent experiments have identified specific structural characteristics of a novel class of lipopeptides (LP) that may lead to potent nanocarriers intended as RNAi therapeutics, albeit the molecular mechanism that underlies their performance remains unexplored. We conducted molecular dynamics simulations in atomistic detail coupled with free energy computations to study the dynamics and thermodynamics of an acrylate- and an epoxide-derived LP, members of the aforesaid class, upon their binding to siRNA in aqueous solution aiming at examining structure-potency relations. We found that the entropic part of the free energy of binding predominates; moreover, the first LP class tends to disrupt the Watson-Crick base pairing of siRNA, whereas the latter leaves the double helix intact. Moreover, the identified tug-of-war effect between LP-water and LP-siRNA hydrogen bonding in the supramolecular complex can underpin synthesis routes toward tuning the association dynamics. Our simulations on two diastereomers of the epoxide-derived LP showed significant structural and energetics differences upon binding, as a result of steric effects imposed by the different absolute configurations at their chiral centers. These findings may serve as crucial design parameters toward modulating the interplay between complex stability and ease of releasing the nucleic acid drug into the cell.

Driving Force for the Complexation of Charged Polypeptides

The phase separation of oppositely charged polyelectrolytes in solution is of current interest. In this work, we study the driving force for polyelectrolyte complexation using molecular dynamics simulations. We calculate the potential of mean force between poly(lysine) and poly(glutamate) oligomers using three different force fields, an atomistic force field and two coarse-grained force fields. There is qualitative agreement between all force fields, i.e., the sign and magnitude of the free energy and the nature of the driving force are similar, which suggests that the molecular nature of water does not play a significant role. For fully charged peptides, we find that the driving force for association is entropic in all cases when small ions either neutralize the poly ions, or are in excess. The removal of all counterions switches the driving force, making complexation energetic. This suggests that the entropy of complexation is dominated by the counterions. When only 6 residues of a 11-mer are charged, however, the driving force is energetic in the abscence of excess salt. The simulations shed insight into the mechanism of complex coacervation and the importance of realistic models for the polyions.

Polymers for siRNA Delivery: A Critical Assessment of Current Technology Prospects for Clinical Application

The number of polymer-based vectors for siRNA delivery in clinical trials lags behind other delivery strategies; however, the molecular architectures and chemical compositions available to polymers make them attractive candidates for further exploration. Polymer vectors are extensively investigated in academic laboratories worldwide with fundamental progress having recently been made in the areas of high-throughput screening, synthetic methods, cellular internalization, endosomal escape and computational prediction and analysis. This review assesses recent advances within the field and highlights relevant developments from within the complementary fields of nanotechnology and protein chemistry with the intent to propose future work that addresses key gaps within the current body of knowledge, potentially advancing the development of the next generation of polymeric vectors.

Martini Coarse-Grained Force Field: Extension to DNA

We systematically parameterized a coarse-grained (CG) model for DNA that is compatible with the Martini force field. The model maps each nucleotide into six to seven CG beads and is parameterized following the Martini philosophy. The CG nonbonded interactions are based on partitioning of the nucleobases between polar and nonpolar solvents as well as base-base potential of mean force calculations. The bonded interactions are fit to single-stranded DNA (ssDNA) atomistic simulations and an elastic network is used to retain double-stranded DNA (dsDNA) and other specific DNA conformations. We present the implementation of the Martini DNA model and demonstrate the properties of individual bases, ssDNA as well as dsDNA, and DNA-protein complexes. The model opens up large-scale simulations of DNA interacting with a wide range of other (bio)molecules that are available within the Martini framework.

The binding and insertion of imidazolium-based ionic surfactants into lipid bilayers: the effects of the surfactant size and salt concentration

Imidazolium-based ionic surfactants with hydrocarbon tails of different sizes were simulated with lipid bilayers at different salt concentrations.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease

BACKGROUND

Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end.

METHODS

We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total).

RESULTS

The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique.

CONCLUSIONS

We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data.

GENERAL SIGNIFICANCE

Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.