Visualizing Biological Membrane Organization and Dynamics

Non-affine deformation analysis and 3D packing defects: A new way to probe membrane heterogeneity in molecular simulations

Here, we discuss a new framework developed over the last 5 years in our group to probe nanoscale membrane heterogeneity. The framework is based on the idea of characterizing lateral heterogeneity through non-affine deformation (NAD) measurements, transverse heterogeneity through three dimensional (3D) lipid packing defects, and using these approaches to formalize the seemingly trivial correlation between lateral organization and lipid packing in biological membranes. We find that measurements from NAD analysis, a prescription which is borrowed from Physics of glasses and granular material, can faithfully distinguish between liquid-ordered and disordered phases in membranes at molecular length scales and, can also be used to identify phase boundaries with high precision. Concomitantly, 3D-packing defects can not only distinguish between the two co-existing fluid phases based on their molecular scale packing (or membrane free volume), but also provide a route to connect the membrane domains to their functionality, such as exploring the molecular origins of inter-leaflet domain registration and peptide partitioning. The correlation between lateral membrane order and transverse packing presents novel molecular design-level features that can explain functions such as protein/peptide partitioning and small-molecule permeation dynamics in complex and heterogeneous membranes with high-fidelity. The framework allows us to explore the nature of lateral organization and molecular packing as a manifestation of intricate molecular interactions among a chemically rich variety of lipids and other molecules in a membrane with complex membrane composition and asymmetry across leaflets.

The Adsorption of P2X2 Receptors Interacting with IgG Antibodies Revealed by Combined AFM Imaging and Mechanical Simulation

The adsorption of proteins onto surfaces significantly impacts biomaterials, medical devices, and biological processes. This study aims to provide insights into the irreversible adsorption process of multiprotein complexes, particularly focusing on the interaction between anti-His6 IgG antibodies and the His6-tagged P2X2 receptor. Traditional approaches to understanding protein adsorption have centered around kinetic and thermodynamic models, often examining individual proteins and surface coverage, typically through Molecular Dynamics (MD) simulations. In this research, we introduce a computational approach employing Autodesk Maya 3D software for the investigation of multiprotein complexes' adsorption behavior. Utilizing Atomic Force Microscopy (AFM) imaging and Maya 3D-based mechanical simulations, our study yields real-time structural and kinetic observations. Our combined experimental and computational findings reveal that the P2X2 receptor-IgG antibody complex likely undergoes absorption in an 'extended' configuration. Whereas the P2X2 receptor is less adsorbed once is complexed to the IgG antibody compared to its individual state, the opposite is observed for the antibody. This insight enhances our understanding of the role of protein-protein interactions in the process of protein adsorption.

UNILIPID, a Methodology for Energetically Accurate Prediction of Protein Insertion into Implicit Membranes of Arbitrary Shape

The insertion of proteins into membranes is crucial for understanding their function in many biological processes. In this work, we present UNILIPID, a universal implicit lipid-protein description as a methodology for dealing with implicit membranes. UNILIPID is independent of the scale of representation and can be applied at the level of all atoms, coarse-grained particles down to the level of a single bead per amino acid. We provide example implementations for these scales and demonstrate the versatility of our approach by accurately reflecting the free energy of transfer for each amino acid. In addition to single membranes, we describe the analytical implementation of double membranes and show that UNILIPID is well suited for modeling at multiple scales. We generalize to membranes of arbitrary shape. With UNILIPID, we provide a methodological framework for a simple and general parameterization tuned to reproduce a selected reference hydrophobicity scale. The software we provide along with the methodological description is optimized for specific user features such as real-time response, live visual analysis, and virtual reality experiences.

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

Probing the Interaction Between Supercarrier RBC Membrane and Nanoparticles for Optimal Drug Delivery

Red blood cell (RBC) membrane-hitchhiking nanoparticles (NPs) have been an increasingly popular supercarrier for targeted drug delivery. However, the kinetic details of the shear-induced NP detachment process from RBC in blood flow remain unclear. Here, we perform detailed computational simulations of the traversal dynamics of an RBC-NP composite supercarrier with tunable properties. We show that the detachment of NPs from RBC occurs in a shear-dependent manner which is consistent with previous experiment results. We quantify the NP detachment rate in the microcapillary flow, and our simulation results suggest that there may be an optimal adhesion strength span of 25-40 μJ/m² for rigid spherical NPs to improve the supercarrier performance and targeting efficiency. In addition, we find that the stiffness and the shape of NPs alter the detachment efficiency by changing the RBC-NP contact areas. Together, these findings provide unique insights into the shear-dependent NP release from the RBC surface, facilitating the clinical utility of RBC-NP composite supercarriers in targeted and localized drug delivery with high precision and efficiency.

A brief history of visualizing membrane systems in molecular dynamics simulations

Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.

Quantitative characterization of the path of glucose diffusion facilitated by human glucose transporter 1

Glucose transporter GLUT1 is ubiquitously expressed in the human body from the red cells to the blood-brain barrier to the skeletal muscles. It is physiologically relevant to understand how GLUT1 facilitates diffusion of glucose across the cell membrane. It is also pathologically relevant because GLUT1 deficiency causes neurological disorders and anemia and because GLUT1 overexpression fuels the abnormal growth of cancer cells. This article presents a quantitative investigation of GLUT1 based on all-atom molecular-dynamics (MD) simulations of the transporter embedded in lipid bilayers of asymmetric inner-and-outer-leaflet lipid compositions, subject to asymmetric intra-and-extra-cellular environments. This is in contrast with the current literature of MD studies that have not considered both of the aforementioned asymmetries of the cell membrane. The equilibrium (unbiased) dynamics of GLUT1 shows that it can facilitate glucose diffusion across the cell membrane without undergoing large-scale conformational motions. The Gibbs free-energy profile, which is still lacking in the current literature of GLUT1, quantitatively characterizes the diffusion path of glucose from the periplasm, through an extracellular gate of GLUT1, on to the binding site, and off to the cytoplasm. This transport mechanism is validated by the experimental data that GLUT1 has low water-permeability, uptake-efflux symmetry, and 10 kcal/mol Arrhenius activation barrier around 37 °C.

Design - a new way to look at old molecules

We discuss how design enriches molecular science, particularly structural biology and bioinformatics. We present two use cases, one in academic practice and the other to design for outreach. The first case targets the representation of ion channels and their dynamic properties. In the second, we document a transition process from a research environment to general-purpose designs. Several testimonials from practitioners are given. By describing the design process of abstracted shapes, exploded views of molecular structures, motion-averaged slices, 360-degree panoramic projections, and experiments with lit sphere shading, we document how designers help make scientific data accessible without betraying its meaning, and how a creative mind adds value over purely data-driven visualizations. A similar conclusion was drawn for public outreach, as we found that comic-book-style drawings are better suited for communicating science to a broad audience.

Skeletal animation for visualizing dynamic shapes of macromolecules

We apply a skeletal animation technique developed for general computer graphics animation to display the dynamic shape of protein molecules. Polygon-based models for macromolecules such as atomic representations, surface models, and protein ribbon models are deformed by the motion of skeletal bones that provide coarse-grained descriptions of detailed computer graphics models. Using the animation software Blender, we developed methods to generate the skeletal bones for molecules. Our example of the superposition of normal modes demonstrates the thermal fluctuating motion obtained from normal mode analysis. The method is also applied to display the motions of protein molecules using trajectory coordinates of a molecular dynamics simulation. We found that a standard motion capture file was practical and useful for describing the motion of the molecule using available computer graphics tools.

The CELLmicrocosmos Tools: A Small History of Java-Based Cell and Membrane Modelling Open Source Software Development

For more than one decade, CELLmicrocosmos tools are being developed. Here, we discus some of the technical and administrative hurdles to keep a software suite running so many years. The tools were being developed during a number of student projects and theses, whereas main developers refactored and maintained the code over the years. The focus of this publication is laid on two Java-based Open Source Software frameworks. Firstly, the CellExplorer with the PathwayIntegration combines the mesoscopic and the functional level by mapping biological networks onto cell components using database integration. Secondly, the MembraneEditor enables users to generate membranes of different lipid and protein compositions using the PDB format. Technicalities will be discussed as well as the historical development of these tools with a special focus on group-based development. In this way, university-associated developers of Integrative Bioinformatics applications should be inspired to go similar ways. All tools discussed in this publication can be downloaded and installed from https://www.CELLmicrocosmos.org.