Molecular dynamics simulations of heterogeneous cell membranes in response to uniaxial membrane stretches at high loading rates

Effects of the structure of lipid-based agents in their complexation with a single stranded mRNA fragment: a computational study

In this work we employed fully atomistic molecular dynamics simulations, aiming towards a better understanding of the mechanisms associated with the formation and the stability of lipid-based RNA nanoassemblies, in an aqueous environment. We examined two groups of lipid-based complexation agents, differing in the degree of hydrophobicity and in the overall charge. The first group was comprised of cationic ionizable agents while the second included electrically neutral amphoteric phosphatidylcholine lipids. It was found that the overall charge of the complexation agents played the most decisive role in the energetics of the lipid/RNA association, while their degree of hydrophobicity affected their self-assembly and their complexation kinetics. The latter also affected the structural stability of the formed complexes since the water entrapped within the clusters of the less hydrophobic agents appeared to reduce the coherence of the lipid-RNA nanoassemblies. The combined effects of the aforementioned attributes dictated also the RNA conformation after complexation. The results from the present study provide thus new insight towards controlling the morphology, the energetic stability and the structural integrity of the formed complexes.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

Complexation of single stranded RNA with an ionizable lipid: an all-atom molecular dynamics simulation study

Complexation of a lipid-based ionizable cationic molecule (referred to as DML: see main text) with RNA in an aqueous medium was examined in detail by means of fully atomistic molecular dynamics simulations. The different stages of the DML-RNA association process were explored, while the structural characteristics of the final complex were described. The self-assembly process of the DML molecules was examined in the absence and in the presence of nucleotide sequences of different lengths. The formed DML clusters were described in detail in terms of their size and composition and were found to share common features in all the examined systems. Different timescales related to their self-assembly and their association with RNA were identified. It was found that beyond a time period of a few tens of ns, a conformationally stable DML-RNA complex was formed, characterized by DML clusters covering the entire contour of RNA. In a system with a 642-nucleotide sequence, the average size of the complex in the longest dimension was found to be close to 40 nm. The DML clusters were characterized by a rather low surface charge, while a propensity for the formation of larger size clusters close to RNA was noted. Apart from hydrophobic and electrostatic interactions, hydrogen bonding was found to play a key-role in the DML-DML and in the DML-RNA association. The information obtained regarding the structural features of the final complex, the timescales and the driving forces associated with the complexation and the self-assembly processes provide new insight towards a rational design of optimized lipid-based ionizable cationic gene delivery vectors.

Coarse-Grained Modeling of Pore Dynamics on the Red Blood Cell Membrane under Large Deformations

Transient pore formation on the membrane of red blood cells (RBCs) under high mechanical tensions is of great importance in many biomedical applications, such as RBC damage (hemolysis) and mechanoporation-based drug delivery. The dynamic process of pore formation, growth, and resealing is hard to visualize in experiments. We developed a mesoscale coarse-grained model to study the characteristics of transient pores on a patch of the lipid bilayer that is strengthened by an elastic meshwork representing the cytoskeleton. Unsteady molecular dynamics was used to study the pore formation and reseal at high strain rates close to the physiological ranges. The critical strain for pore formation, pore characteristics, and cytoskeleton effects were studied. Results show that the presence of the cytoskeleton increases the critical strain of pore formation and confines the pore growth. Moreover, the pore recovery process under negative strain rates (compression) is analyzed. Simulations show that pores can remain open for a long time during the high-speed tank-treading induced stretching and compression process that a patch of the RBC membrane usually experiences under high shear flow. Furthermore, complex loading conditions can affect the pore characteristics and result in denser pores. Finally, the effects of strain rate on pore formation are analyzed. Higher rate stretching of membrane patch can result in a significant increase in the critical areal strain and density of pores. Such a model reveals the dynamic molecular process of RBC damage in biomedical devices and mechanoporation that, to our knowledge, has not been reported before.

Localized Axolemma Deformations Suggest Mechanoporation as Axonal Injury Trigger

Traumatic brain injuries are a leading cause of morbidity and mortality worldwide. With almost 50% of traumatic brain injuries being related to axonal damage, understanding the nature of cellular level impairment is crucial. Experimental observations have so far led to the formulation of conflicting theories regarding the cellular primary injury mechanism. Disruption of the axolemma, or alternatively cytoskeletal damage has been suggested mainly as injury trigger. However, mechanoporation thresholds of generic membranes seem not to overlap with the axonal injury deformation range and microtubules appear too stiff and too weakly connected to undergo mechanical breaking. Here, we aim to shed a light on the mechanism of primary axonal injury, bridging finite element and molecular dynamics simulations. Despite the necessary level of approximation, our models can accurately describe the mechanical behavior of the unmyelinated axon and its membrane. More importantly, they give access to quantities that would be inaccessible with an experimental approach. We show that in a typical injury scenario, the axonal cortex sustains deformations large enough to entail pore formation in the adjoining lipid bilayer. The observed axonal deformation of 10–12% agree well with the thresholds proposed in the literature for axonal injury and, above all, allow us to provide quantitative evidences that do not exclude pore formation in the membrane as a result of trauma. Our findings bring to an increased knowledge of axonal injury mechanism that will have positive implications for the prevention and treatment of brain injuries.

Effect of pressure profile of shock waves on lipid membrane deformation

Use of shock waves to temporarily increase the permeability of the cell membrane is a promising approach in drug delivery and gene therapy to allow the translocation of macromolecules and small polar molecules into the cytoplasm. Our understanding of how the characteristics of the pressure profile of shock waves, such as peak pressure and pulse duration, influences membrane properties is limited. Here we study the response of lipid bilayer membranes to shock pulses with different pressure profiles using atomistic molecular dynamics simulations. From our simulation results, we find that the transient deformation/disordering of the membrane depends on both the magnitude and the pulse duration of the pressure profile of the shock pulse. For a low pressure impulse, peak pressure has a dominant effect on membrane structural changes, while for the high pressure impulse, we find that there exists an optimal pulse duration at which membrane deformation/disordering is maximized.

Cognition based bTBI mechanistic criteria; a tool for preventive and therapeutic innovations

Blast-induced traumatic brain injury has been associated with neurodegenerative and neuropsychiatric disorders. To date, although damage due to oxidative stress appears to be important, the specific mechanistic causes of such disorders remain elusive. Here, to determine the mechanical variables governing the tissue damage eventually cascading into cognitive deficits, we performed a study on the mechanics of rat brain under blast conditions. To this end, experiments were carried out to analyse and correlate post-injury oxidative stress distribution with cognitive deficits on a live rat exposed to blast. A computational model of the rat head was developed from imaging data and validated against in vivo brain displacement measurements. The blast event was reconstructed in silico to provide mechanistic thresholds that best correlate with cognitive damage at the regional neuronal tissue level, irrespectively of the shape or size of the brain tissue types. This approach was leveraged on a human head model where the prediction of cognitive deficits was shown to correlate with literature findings. The mechanistic insights from this work were finally used to propose a novel protective device design roadmap and potential avenues for therapeutic innovations against blast traumatic brain injury.