Modulation of Membrane Fluidity to Control Interfacial Water Structure and Dynamics in Saturated and Unsaturated Phospholipid Vesicles

Flexible RNA aptamers as inhibitors of Bacillus anthracis ribosomal protein S8: Insights from molecular dynamics simulations

Bacillus anthracis, the causative agent of anthrax, poses a substantial threat to public health and national security, and is recognized as a potential bioweapon due to its capacity to form resilient spores with enduring viability. Inhalation or ingestion of even minute quantities of aerosolized spores can lead to widespread illness and fatalities, underscoring the formidable lethality of the bacterium. With an untreated mortality rate of 100%, Bacillus anthracis is a disconcerting candidate for bioterrorism. In response to this critical scenario, we employed state-of-the-art computational tools to conceive and characterize flexible RNA aptamer therapeutics tailored for anthrax. The foundational structure of the flexible RNA aptamers was designed by removing the C2'-C3' in each nucleotide unit. Leveraging the crystal structure of Bacillus anthracis ribosomal protein S8 complexed with an RNA aptamer, we explored the structural, dynamic, and energetic aspects of the modified RNA aptamer - S8 protein complexes through extensive all-atom explicit-solvent molecular dynamics simulations (400 ns, 3 replicas each), followed by drawing comparisons to the control system. Our findings demonstrate the enhanced binding competencies of the flexible RNA aptamers to the S8 protein via better shape complementarity and improved H-bond network compared to the control RNA aptamer. This research offers valuable insights into the development of RNA aptamer therapeutics targeting Bacillus anthracis, paving the way for innovative strategies to mitigate the impact of this formidable pathogen.

Real Space and Time Imaging of Collective Headgroup Dipole Motions in Zwitterionic Lipid Bilayers

Lipid bilayers are supramolecular structures responsible for a range of processes, such as transmembrane transport of ions and solutes, and sorting and replication of genetic materials, to name just a few. Some of these processes are transient and currently, cannot be visualized in real space and time. Here, we developed an approach using 1D, 2D, and 3D Van Hove correlation functions to image collective headgroup dipole motions in zwitterionic phospholipid bilayers. We show that both 2D and 3D spatiotemporal images of headgroup dipoles are consistent with commonly understood dynamic features of fluids. However, analysis of the 1D Van Hove function reveals lateral transient and re-emergent collective dynamics of the headgroup dipoles-occurring at picosecond time scales-that transmit and dissipate heat at longer times, due to relaxation processes. At the same time, the headgroup dipoles also generate membrane surface undulations due a collective tilting of the headgroup dipoles. A continuous intensity band of headgroup dipole spatiotemporal correlations-at nanometer length and nanosecond time scales-indicates that dipoles undergo stretching and squeezing elastic deformations. Importantly, the above mentioned intrinsic headgroup dipole motions can be externally stimulated at GHz-frequency scale, enhancing their flexoelectric and piezoelectric capabilities (i.e., increased conversion efficiency of mechanical energy into electric energy). In conclusion, we discuss how lipid membranes can provide molecular-level insights about biological learning and memory, and as platforms for the development of the next generation of neuromorphic computers.

The effect of rhamnolipids on fungal membrane models as described by their interactions with phospholipids and sterols: An in silico study

Introduction: Rhamnolipids (RLs) are secondary metabolites naturally produced by bacteria of the genera Pseudomonas and Burkholderia with biosurfactant properties. A specific interest raised from their potential as biocontrol agents for crop culture protection in regard to direct antifungal and elicitor activities. As for other amphiphilic compounds, a direct interaction with membrane lipids has been suggested as the key feature for the perception and subsequent activity of RLs. Methods: Molecular Dynamics (MD) simulations are used in this work to provide an atomistic description of their interactions with different membranous lipids and focusing on their antifungal properties. Results and discussion: Our results suggest the insertion of RLs into the modelled bilayers just below the plane drawn by lipid phosphate groups, a placement that is effective in promoting significant membrane fluidification of the hydrophobic core. This localization is promoted by the formation of ionic bonds between the carboxylate group of RLs and the amino group of the phosphatidylethanolamine (PE) or phosphatidylserine (PS) headgroups. Moreover, RL acyl chains adhere to the ergosterol structure, forming a significantly higher number of van der Waals contact with respect to what is observed for phospholipid acyl chains. All these interactions might be essential for the membranotropic-driven biological actions of RLs.

Bicyclo-DNA mimics with enhanced protein binding affinities: insights from molecular dynamics simulations

DNA-protein interactions occur at all levels of DNA expression and replication and are crucial determinants for the survival of a cell. Several modified nucleotides have been utilized to manipulate these interactions and have implications in drug discovery. In the present article, we evaluated the binding of bicyclo-nucleotides (generated by forming a methylene bridge between C1' and C5' in sugar, leading to a bicyclo system with C2' axis of symmetry at the nucleotide level) to proteins. We utilized four ssDNA-protein complexes with experimentally known binding free energies and investigated the binding of modified nucleotides to proteins via all-atom explicit solvent molecular dynamics (MD) simulations (200 ns), and compared the binding with control ssDNA-protein systems. The modified ssDNA displayed enhanced binding to proteins as compared to the control ssDNA, as seen by means of MD simulations followed by MM-PBSA calculations. Further, the Delphi-based electrostatic estimation revealed that the high binding of modified ssDNA to protein might be related to the enhanced electrostatic complementarity displayed by the modified ssDNA molecules in all the four systems considered for the study. The improved binding achieved with modified nucleotides can be utilized to design and develop anticancer/antisense molecules capable of targeting proteins or ssRNAs.Communicated by Ramaswamy H. Sarma.

How Far Is "Bulk Water" from Interfaces? Depends on the Nature of the Surface and What We Measure

Using systematic molecular dynamics (MD) simulations, we revisit the question: At what distance from an interface do the properties of "bulk water" get recovered? We have considered three different kinds of interfaces: nonpolar (hydrophobic; isooctane-water interface), charged (negative; AOT bilayer), and polar (zwitterionic; POPC bilayer). In order to interrogate the extent of perturbation of the interfacial water molecules as a function of the distance from the interface, we utilize a diverse range of structural and dynamical parameters. To capture the structural perturbations, we look into local density (translational order), local tetrahedral order parameter, and dipolar orientation of the water molecules. We also explore the anisotropic diffusion of the water molecules in the direction perpendicular to the interface as well as the planar diffusion parallel to the interface in a distance dependent manner. In addition, the orientational time correlation functions have been computed to understand the extent of slowdown in the rotational dynamics. As expected, the electrostatic field emanating from the charged AOT interface seems to have the highest long-range effect on the orientational order and dynamics of the water molecules, whereas specific interactions like hydrogen bonding and electrostatic interaction lead to significant trapping and kinetic slowdown for both AOT and POPC (zwitterionic) very close to the interface. Our analysis highlights that not only the length-scale of perturbation depends on the nature of the interfaces and specific interactions but also the type of water property that we measure/calculate. Different water properties seem to have widely different length-scale of perturbation. Orientational order parameters seem to be perturbed to a much longer length-scale as compared to translational order parameters. The global orientational order of water can be perturbed even up to ∼4-5 nm near the negatively charged AOT surface in the absence of any extra electrolyte. This observation has significant implication toward the interpretation of experimental measurements as well since different spectroscopic techniques would probe different parameters or water properties with possible mutual disagreement and inconsistency between different types of measurements. Thus, our study provides a broader and unifying perspective toward the aspect of "context dependent" structural and dynamical perturbation of "interfacial water".

Dehydration induced dynamical heterogeneity and ordering mechanism of lipid bilayers

Understanding the influence of dehydration on the membrane structure is crucial to control membrane functionality related to domain formation and cell fusion under anhydrobiosis conditions. To this end, we perform all-atom molecular dynamic simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dimyristoylphosphatidylcholine lipid membranes at different hydration levels at 308 K. As dehydration increases, the lipid area per head group decreases with an increase in bilayer thickness and lipid order parameters indicating bilayer ordering. Concurrently, translational and rotational dynamics of interfacial water (IW) molecules near membranes slow down. On the onset of bilayer ordering, the IW molecules exhibit prominent features of dynamical heterogeneity evident from non-Gaussian parameters and one-dimensional van Hove correlation functions. At a fully hydrated state, diffusion constants (D) of the IW follow a scaling relation, D∼τ_α ^-1, where the α relaxation time (τ_α) is obtained from self-intermediate scattering functions. However, upon dehydration, the relation breaks and the D of the IW follows a power law behavior as D∼τ_α ^-0.57, showing the signature of glass dynamics. τ_α and hydrogen bond lifetime calculated from intermittent hydrogen bond auto-correlation functions undergo a similar crossover in association with bilayer ordering on dehydration. The bilayer ordering is accompanied with an increase in fraction of caged lipids spanned over the bilayer surface and a decrease in fraction of mobile lipids due to the non-diffusive dynamics. Our analyses reveal that the microscopic mechanism of lipid ordering by dehydration is governed by dynamical heterogeneity. The fundamental understanding from this study can be applied to complex bio-membranes to trap functionally relevant gel-like domains at room temperature.

A Comparative Study on DMSO-Induced Modulation of the Structural and Dynamical Properties of Model Bilayer Membranes

Modulating the structures and properties of biomembranes via permeation of small amphiphilic molecules is immensely important, having diverse applications in cell biology, biotechnology, and pharmaceuticals, because their physiochemical and biological interactions lead to new pathways for transdermal drug delivery and administration. In this work, we have elucidated the role of dimethyl sulfoxide (DMSO), broadly used as a penetration-enhancing agent and cryoprotective agent on model lipid membranes, using a combination of fluorescence microscopy and time-resolved fluorescence spectroscopy. Spatially resolved fluorescence lifetime imaging microscopy (FLIM) has been employed to unravel how the fluidity of the DMSO-induced bilayer regulates the structural alteration of the vesicles. Moreover, we have also shown that the dehydration effect of DMSO leads to weakening of the hydrogen bond between lipid headgroups and water molecules and results in faster solvation dynamics as demonstrated by femtosecond time-resolved fluorescence spectroscopy. It has been gleaned that the water dynamics becomes faster because bilayer rigidity decreases in the presence of DMSO, which is also supported by time-resolved rotational anisotropy measurements. The enhanced diffusivity and increased membrane fluidity in the presence of DMSO are further ratified at the single-molecule level through fluorescence correlation spectroscopy (FCS) measurements. Our results indicate that while the presence of DMSO significantly affects the 1,2-dimyristoyl-rac-glycero-3-phosphocholine (DMPC) and 1,2-dipalmitoyl-rac-glycero-3-phosphatidylcholine (DPPC) bilayers, it has a weak effect on 1,2-dimyristoyl-sn-glycero-3-phospho-rac-glycerol (DMPG) vesicles, which might explain the preferential interaction of DMSO with the positively charged choline group present in DMPC and DPPC vesicles. The experimental findings have also been further verified with molecular dynamics simulation studies. Moreover, it has been observed that DMSO is likely to have a differential effect on heterogeneous bilayer membranes depending on the structure and composition of their headgroups. Our results illuminate the importance of probing the lipid structure and composition of cellular membranes in determining the effects of cryoprotective agents.