Lipid Bilayer Membrane Perturbation by Embedded Nanopores: A Simulation Study

Unlocking the specificity of antimicrobial peptide interactions for membrane-targeted therapies

Antimicrobial peptides (AMPs) are increasingly recognized as potent therapeutic agents, with their selective affinity for pathological membranes, low toxicity profile, and minimal resistance development making them particularly attractive in the pharmaceutical landscape. This study offers a comprehensive analysis of the interaction between specific AMPs, including magainin-2, pleurocidin, CM15, LL37, and clavanin, with lipid bilayer models of very different compositions that have been ordinarily used as biological membrane models of healthy mammal, cancerous, and bacterial cells. Employing unbiased molecular dynamics simulations and metadynamics techniques, we have deciphered the intricate mechanisms by which these peptides recognize pathogenic and pathologic lipid patterns and integrate into lipid assemblies. Our findings reveal that the transverse component of the peptide's hydrophobic dipole moment is critical for membrane interaction, decisively influencing the molecule's orientation and expected therapeutic efficacy. Our approach also provides insight on the kinetic and dynamic dependence on the peptide orientation in the axial and azimuthal angles when coming close to the membrane. The aim is to establish a robust framework for the rational design of peptide-based, membrane-targeted therapies, as well as effective quantitative descriptors that can facilitate the automated design of novel AMPs for these therapies using machine learning methods.

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.

Single-Molecule Electrical Profiling of Peptides and Proteins

In recent decades, there has been a significant increase in the application of single-molecule electrical analysis platforms in studying proteins and peptides. These advanced analysis methods have the potential for deep investigation of enzymatic working mechanisms and accurate monitoring of dynamic changes in protein configurations, which are often challenging to achieve in ensemble measurements. In this work, the prominent research progress in peptide and protein-related studies are surveyed using electronic devices with single-molecule/single-event sensitivity, including single-molecule junctions, single-molecule field-effect transistors, and nanopores. In particular, the successful commercial application of nanopores in DNA sequencing has made it one of the most promising techniques in protein sequencing at the single-molecule level. From single peptides to protein complexes, the correlation between their electrical characteristics, structures, and biological functions is gradually being established. This enables to distinguish different molecular configurations of these biomacromolecules through real-time electrical monitoring of their life activities, significantly improving the understanding of the mechanisms underlying various life processes.

Experimentally Probing the Effect of Confinement Geometry on Lipid Diffusion

The lateral mobility of molecules within the cell membrane is ultimately governed by the local environment of the membrane. Confined regions induced by membrane structures, such as protein aggregates or the actin meshwork, occur over a wide range of length scales and can impede or steer the diffusion of membrane components. However, a detailed picture of the origins and nature of these confinement effects remains elusive. Here, we prepare model lipid systems on substrates patterned with confined domains of varying geometries constructed with different materials to explore the influences of physical boundary conditions and specific molecular interactions on diffusion. We demonstrate a platform that is capable of significantly altering and steering the long-range diffusion of lipids by using simple oxide deposition approaches, enabling us to systematically explore how confinement size and shape impact diffusion over multiple length scales. While we find that a "boundary condition" description of the system captures underlying trends in some cases, we are also able to directly compare our systems to analytical models, revealing the unexpected breakdown of several approximate solutions. Our results highlight the importance of considering the length scale dependence when discussing properties such as diffusion.

Nanoporous Membranes of Densely Packed Carbon Nanotubes Formed by Lipid-Mediated Self-Assembly

Nanofiltration technology faces the competing challenges of achieving high fluid flux through uniformly narrow pores of a mechanically and chemically stable filter. Supported dense-packed 2D-crystals of single-walled carbon nanotube (CNT) porins with ∼1 nm wide pores could, in principle, meet these challenges. However, such CNT membranes cannot currently be synthesized at high pore density. Here, we use computer simulations to explore lipid-mediated self-assembly as a route toward densely packed CNT membranes, motivated by the analogy to membrane-protein 2D crystallization. In large-scale coarse-grained molecular dynamics (MD) simulations, we find that CNTs in lipid membranes readily self-assemble into large clusters. Lipids trapped between the CNTs lubricate CNT repacking upon collisions of diffusing clusters, thereby facilitating the formation of large ordered structures. Cluster diffusion follows the Saffman-Delbrück law and its generalization by Hughes, Pailthorpe, and White. On longer time scales, we expect the formation of close-packed CNT structures by depletion of the intervening shared annular lipid shell, depending on the relative strength of CNT-CNT and CNT-lipid interactions. Our simulations identify CNT length, diameter, and end functionalization as major factors for the self-assembly of CNT membranes.

Uncovering the mechanisms of cyclic peptide self-assembly in membranes with the chirality-aware MA(R/S)TINI forcefield

Cyclic peptides (CPs) formed by alternation of D- and L-amino acids (D,L-CPs) can self-assemble into nanotubes (SCPNs) by parallel or/and antiparallel stacking. Different applications have been attributed to these nanotubes, including the disruption of lipid bilayers of specific compositions and the selective transport of ions throughout membranes. Molecular dynamics (MD) simulations have significantly contributed to understand the interaction between CPs, including the structural, dynamic and transport properties of their supramolecular aggregates. The high computational cost of atomic resolution forcefields makes them impractical for simulating the self-assembly of macromolecules, so coarse-grained (CG) models might represent a more feasible solution for this purpose. However, general CG models used for the simulation of biomolecules such as the MARTINI forcefield do not explicitly consider the non-covalent interactions leading to the formation of secondary structure patterns in proteins. This becomes particularly important in the case of CPs due to the D- and L-chirality alternation in their sequence, leading to opposite orientations of the backbone polar groups on both sides of the cyclic ring plane. In order to overcome this limitation, we have extended the MARTINI forcefield to introduce chirality in each residue of the CPs. The new parametrization, which we have called MA(R/S)TINI, reproduces the expected self-assembly patterns for several CP sequences in the presence of different membrane models, explicitly considering the chirality of the CPs and with no significant extra computational cost. Our simulations provide new mechanistic information of how these systems self-assemble in presence of different lipid scenarios, showing that the CP-CP and CP-membrane interactions are sensitive to the peptide sequence chirality. This opens the door to design new bioactive CPs based on CG-MD simulations. A web-based tool for the automatic parameterization of new CP sequences using MA(R/S)TINI, among other functionalities, is under construction (see http://cyclopep.com).

Ultrasensitive Detection of MCF-7 Cells with a Carbon Nanotube-Based Optoelectronic-Pulse Sensor Framework

Biosensors are of vital significance for healthcare by supporting the management of infectious diseases for preventing pandemics and the diagnosis of life-threatening conditions such as cancer. However, the advancement of the field can be limited by low sensing accuracy. Here, we altered the bioelectrical signatures of the cells using carbon nanotubes (CNTs) via structural loosening effects. Using an alternating current (AC) pulse under light irradiation, we developed a photo-assisted AC pulse sensor based on CNTs to differentiate between healthy breast epithelial cells (MCF-10A) and luminal breast cancer cells (MCF-7) within a heterogeneous cell population. We observed a previously undemonstrated increase in current contrast for MCF-7 cells with CNTs compared to MCF-10A cells with CNTs under light exposure. Moreover, we obtained a detection limit of ∼1.5 × 103 cells below a baseline of ∼1 × 104 cells for existing electrical-based sensors for an adherent, heterogeneous cell population. All-atom molecular dynamics (MD) simulations reveal that interactions between the embedded CNT and cancer cell membranes result in a less rigid lipid bilayer structure, which can facilitate CNT translocation for enhancing current. This as-yet unconsidered cancer cell-specific method based on the unique optoelectrical properties of CNTs represents a strategy for unlocking the detection of a small population of cancer cells and provides a promising route for the early diagnosis, monitoring, and staging of cancer.

Role of polyplex charge density in lipopolyplexes

Lipopolyplexes have received extensive attention lately in gene therapy delivery. However, the interactions between the polyplex and the liposome and their underlying molecular mechanisms remain to be elucidated. Here, we adopted a simple model, mainly to illustrate the impact of polyplex charge density on the self-assembly of liposomes (containing DOPE and CHEMS lipids) using coarse-grained molecular dynamics simulations. Our simulation results show that when the charge density increases in the polyplex, more lipids, especially CHEMS (a negatively charged helper lipid) lipids, are attracted to the polyplex (positively charged) surface, and meanwhile nearby water molecules are driven away from the polyplex, resulting in a less spherical liposome. Energy decomposition analyses further reveal that, at higher charge densities, the polyplex exhibits much stronger interactions with CHEMS lipids than with water molecules, with the majority contribution from electrostatic interactions. In addition, the mobility of lipids, especially CHEMS, is reduced as the polyplex charge density increases, indicating a more rigid liposome. Overall, our molecular dynamics simulations elucidate the influence of polyplex charge density on the liposome self-assembly process at the atomic level, which provides a complementary approach to experiments for a better understanding of this promising gene therapy delivery system.

Unravelling hierarchical levels of structure in lipid membranes

•

A set of hierarchical levels of structure can be proposed for lipid bilayers.

•

The composition of a lipid bilayer is identified as a fluid version of primary structure.

•

The interaction between leaflets is taken as the secondary structure of lipid bilayers.

•

A method to identify membrane domains and their interaction is proposed as a third level of structure.

•

The highly specific lipid perturbation around embedded macromolecules is taken as the quaternary level of structure.

In analogy with the hierarchical levels typically used to describe the structure of nucleic acids or proteins and keeping in mind that lipid bilayers are not just mere envelopers for biological material but directly responsible for many important functions of life, it is discussed here how membrane models can also be interpreted in terms of different hierarchies in their structure. Namely, lipid composition, interaction between leaflets, existence and interaction of domains arising from the coordinate behavior of lipids and their properties, plus the manifest and specific perturbation of the lipid organization around macromolecules embedded in a membrane are hereby used to define the primary, secondary, tertiary and quaternary structures, respectively. Molecular Dynamics simulations are used to illustrate this proposal. Alternative levels of organization and methods to define domains can be proposed but the final aim is to highlight the paradigm arising from this description which is expected to have significant consequences on deciphering the underlying factors governing membranes and their interactions with other molecules.

Triggered Assembly of a DNA-Based Membrane Channel

Chemistry is in a powerful position to synthetically replicate biomolecular structures. Adding functional complexity is key to increase the biomimetics' value for science and technology yet is difficult to achieve with poorly controlled building materials. Here, we use defined DNA blocks to rationally design a triggerable synthetic nanopore that integrates multiple functions of biological membrane proteins. Soluble triggers bind via molecular recognition to the nanopore components changing their structure and membrane position, which controls the assembly into a defined channel for efficient transmembrane cargo transport. Using ensemble, single-molecule, and simulation analysis, our activatable pore provides insight into the kinetics and structural dynamics of DNA assembly at the membrane interface. The triggered channel advances functional DNA nanotechnology and synthetic biology and will guide the design of controlled nanodevices for sensing, cell biological research, and drug delivery.

Molecular Dynamics Simulations of Transmembrane Cyclic Peptide Nanotubes Using Classical Force Fields, Hydrogen Mass Repartitioning, and Hydrogen Isotope Exchange Methods: A Critical Comparison

Self-assembled cyclic peptide nanotubes with alternating D- and L-amino acid residues in the sequence of each subunit have attracted a great deal of attention due to their potential for new nanotechnology and biomedical applications, mainly in the field of antimicrobial peptides. Molecular dynamics simulations can be used to characterize these systems with atomic resolution at different time scales, providing information that is difficult to obtain via wet lab experiments. However, the performance of classical force fields typically employed in the simulation of biomolecules has not yet been extensively tested with this kind of highly constrained peptide. Four different classical force fields (AMBER, CHARMM, OPLS, and GROMOS), using a nanotube formed by eight D,L-α-cyclic peptides inserted into a lipid bilayer as a model system, were employed here to fill this gap. Significant differences in the pseudo-cylindrical cavities formed by the nanotubes were observed, the most important being the diameter of the nanopores, the number and location of confined water molecules, and the density distribution of the solvent molecules. Furthermore, several modifications were performed on GROMOS54a7, aiming to explore acceleration strategies of the MD simulations. The hydrogen mass repartitioning (HMR) and hydrogen isotope exchange (HIE) methods were tested to slow down the fastest degrees of freedom. These approaches allowed a significant increase in the time step employed in the equation of the motion integration algorithm, from 2 fs up to 5-7 fs, with no serious changes in the structural and dynamical properties of the nanopores. Subtle differences with respect to the simulations with the unmodified force fields were observed in the concerted movements of the cyclic peptides, as well as in the lifetime of several H-bonds. All together, these results are expected to contribute to better understanding of the behavior of self-assembled cyclic peptide nanotubes, as well as to support the methods tested to speed up general MD simulations; additionally, they do provide a number of quantitative descriptors that are expected to be used as a reference to design new experiments intended to validate and complement computational studies of antimicrobial cyclic peptides.

Membrane Nanopores Induced by Nanotoroids via an Insertion and Pore-Forming Pathway

The formation of membrane nanopores is one of the crucial activities of cells and has attracted considerable attention. However, the understanding of their types and mechanisms is still limited. Herein, we report a novel nanopore formation phenomenon achieved through the insertion of polymeric nanotoroids into the cellular membrane. As revealed by theoretical simulations, the nanotoroid can embed in the membrane, leaving a nanopore on the cell. The through-the-cavity wrapping of lipids is critical for the retention of the nanotoroid in the membrane, which is attributed to both a relatively large inner cavity of the nanotoroid and a moderate attraction between the nanotoroid and membrane lipids. Under the guidance of the simulation predictions, experiments using polypeptide toroids as pore-forming agents were performed, confirming the unique biophysical phenomenon. This work demonstrates a distinctive pore-forming pathway, deepens the understanding of the membrane nanopore phenomenon, and assists in the design of advanced pore-forming materials.

Transmembrane Self-Assembled Cyclic Peptide Nanotubes Based on α-Residues and Cyclic δ-Amino Acids: A Computational Study

Self-assembling cyclic peptide nanotubes have been shown to function as synthetic, integral transmembrane channels. The combination of natural and nonnatural aminoacids in the sequence of cyclic peptides enables the control not only of their outer surface but also of the inner cavity behavior and properties, affecting, for instance, their permeability to different molecules including water and ions. Here, a thorough computational study on a new class of self-assembling peptide motifs, in which δ-aminocycloalkanecarboxylic acids are alternated with natural α-amino acids, is presented. The presence of synthetic δ-residues creates hydrophobic regions in these α,δ-SCPNs, which makes them especially attractive for their potential implementation in the design of new drug or diagnostic agent carrier systems. Using molecular dynamics simulations, the behavior of water molecules, different ions (Li+, Na+, K+, Cs+, and Ca2+), and their correspondent counter Cl- anions is extensively investigated in the nanoconfined environment. The structure and dynamics are mutually combined in a diving immersion inside these transmembrane channels to discover a fascinating submarine nanoworld where star-shaped water channels guide the passage of cations and anions therethrough.

Effect of Water Models on Transmembrane Self-Assembled Cyclic Peptide Nanotubes

Self-assembling cyclic peptide nanotubes can form nanopores when they are inserted in lipid bilayers, acting as ion and/or water permeable channels. In order to improve the versatility of these systems, it is possible to specifically design cyclic peptides with a combination of natural and non-natural amino acids, enabling the control of the nature of the inner cavity of the channels. Here, the behavior of two types of self-assembling peptide motifs, alternating α-amino acids with γ- or δ-aminocycloalkanecarboxylic acids, is studied via molecular dynamics (MD) simulations. The behavior of water molecules in nanopores is expected to affect the properties of these channels and therefore merits detailed examination. A number of water models commonly used in MD simulations have been validated by how well they reproduce bulk water properties. However, it is less clear how these water models behave in the nanoconfined condition inside a channel. The behavior of four different water models-TIP3P, TIP4P, TIP4P/2005, and OPC-are evaluated in MD simulations of self-assembled cyclic peptide nanotubes of distinct composition and diameter. The dynamic behavior of the water molecules and ions in these designed artificial channels depends subtly on the water model used. TIP3P water molecules move faster than those of TIP4P, TIP4P/2005, and OPC. This demeanor is clearly observed in the filling of the nanotube, in water diffusion within the pore, and in the number and stability of hydrogen bonds of the peptides with water. It was also shown that the water model influences the simulated ion flux through the nanotubes, with TIP3P producing the greatest ion flux. Additionally, the two more recent models, TIP4P/2005 and OPC, which are known to reproduce the experimental self-diffusion coefficient of bulk water quite well, exhibit very similar results under the nanoconfined conditions studied here. Because none of these models have been parametrized specifically for waters confined in peptide nanotubes, this study provides a point of reference for further validation.

Molecular mechanisms underlying the role of the puckered surface in the biocompatibility of black phosphorus

As a newly emerging two-dimensional material, black phosphorus (BP) has received broad attention in the field of biomedical applications. Prior to its clinical application, its cytotoxicity to cells should be carefully evaluated; however, this field is still in its infancy. Motivated by this, we performed molecular dynamics (MD) simulations to systematically investigate the potential mechanisms of the cytotoxicity of BP to the lipid membrane, including lipid extraction, penetration into the membrane, and the impacts of BP on the physical properties of the membrane. Surprisingly, we observed that BP could not extract lipid molecules from the membrane. The thermodynamic analyses suggested that the puckered surface structure could weaken the interactions between BP and lipid molecules, thus inhibiting the lipid extraction. Additionally, through simulating the spontaneous interaction modes between BP and the lipid membrane, we found that the "passivated" edges of BP prohibited it from penetrating into the membrane. As a result, BP could only spontaneously lie parallel on the surface of the membrane, in which manner BP exerted little influence on the properties of the lipid membrane. To comprehensively appraise the cytotoxicity, we even artificially inserted BP into the membrane and compared the effects of BP and graphene on the properties of the membrane. Simulation results showed that the influences of the inserted BP on the lipid properties were much milder than those of graphene. Overall, the present work suggests that BP possesses distinctive biocompatibility benefiting from its puckered surface structure. This work provides a better understanding of the interactions between BP and the membrane, which may offer some useful suggestions for exploring strategies to improve the biocompatibility of nanomaterials.

Tim4 recognizes carbon nanotubes and mediates phagocytosis leading to granuloma formation

Macrophage recognition and phagocytosis of crystals is critical for the associated fibrosis and cancer. Of note, multi-walled carbon nanotubes (MWCNTs), the highly representative products of nanotechnology, induce macrophage NLRP3 inflammasome activation and cause asbestosis-like pathogenesis. However, it remains largely unknown how macrophages efficiently recognize MWCNTs on their cell surfaces. Here, we identify by a targeted screening of phagocyte receptors the phosphatidylserine receptors T cell immunoglobulin mucin 4 (Tim4) and Tim1 as the pattern-recognition receptors for carbon crystals. Docking simulation studies reveal spatiotemporally stable interfaces between aromatic residues in the extracellular IgV domain of Tim4 and one-dimensional carbon crystals. Further, CRISPR-Cas9-mediated deletion of Tim4 and Tim1 reveals that Tim4, but not Tim1, critically contributes to the recognition of MWCNTs by peritoneal macrophages and to granuloma development in a mouse model of direct mesothelium exposure to MWCNTs. These results suggest that Tim4 recognizes MWCNTs through aromatic interactions and mediates phagocytosis leading to granulomas.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Membrane interactions in drug delivery: Model cell membranes and orthogonal techniques

To generate the desired effect in the human body, the active pharmaceutical ingredient usually needs to interact with a receptor located on the cell membrane or inside the cell. Thus, understanding membrane interactions is of great importance when it comes to the development and testing of new drug molecules or new drug delivery systems. Nowadays, there is a tremendous selection of both model cell membranes and of techniques that can be used to characterize interactions between selected model cell membranes and a drug molecule, an excipient, or a drug delivery system. Having such a wide selection of model cell membranes and techniques available makes it sometimes challenging to select the optimal combination for a specific study. Furthermore, it is difficult to compare results obtained using different model cell membranes and techniques, and not all in vitro studies translate as well to an estimation of the in vivo biological activity or understanding of mode of action. This review provides an overview of the available lipid bilayer-based model cell membranes and of the most widely employed techniques for studying membrane interactions. Finally, the need for employing complimentary characterization techniques in order to acquire more reliable and in-depth information is highlighted.

Parallel Versus Antiparallel β-Sheet Structure in Cyclic Peptide Hybrids Containing γ- or δ-Cyclic Amino Acids

Cyclic peptides with disc-shaped structures have emerged as potent building blocks for the preparation of new biomaterials in fields ranging from biological to material science. In this work, we analyze in depth the self-assembling properties of a new type of cyclic peptides based on the alternation of α-residues and cyclic δ-amino acids (α,δ-CPs). To examine the preferred stacking properties adopted by cyclic peptides bearing this type of amino acids, we carried out a synergistic in vitro/in silico approximation by using simple dimeric models and then extended to nanotubes. Although these new cyclic peptides (α,δ-CPs) can interact either in a parallel or antiparallel fashion, our results confirm that although the parallel β-sheet is more stable, it can be switched to the antiparallel stacking by choosing residues that can establish favorable cross-strand interactions. Moreover, the subsequent comparison by using the same methodology but applied to α,γ-CPs models, up to the moment assumed as antiparallel-like d,l-α-CPs, led to unforeseen conclusions that put into question preliminary conjectures about these systems. Surprisingly, they tend to adopt a parallel β-sheet directed by the skeleton interactions. These results imply a change of paradigm with respect to cyclic peptide designs that should be considered for dimers and nanotubes.

Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif

Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.

Rings, Hexagons, Hetals, and Dipolar Moment Sink-Sources: The Fanciful Behavior of Water around Cyclodextrin Complexes

The basket-like geometry of cyclodextrins (CDs), with a cavity able to host hydrophobic groups, makes these molecules well suited for a large number of fundamental and industrial applications. Most of the established CD-based applications rely on trial and error studies, often ignoring key information at the atomic level that could be employed to design new products and to optimize their use. Computational simulations are well suited to fill this gap, especially in the case of CD systems due to their low number of degrees of freedom compared with typical macromolecular systems. Thus, the design and validation of solid and efficient methods to simulate and analyze CD-based systems is key to contribute to this field. The behavior of supramolecular complexes critically depends on the media where they are embedded, so the detailed characterization of the solvent is required to fully understand these systems. In the present work, we use the inclusion complex formed by two α-CDs and one sodium dodecyl sulfate molecule to test eight different parameterizations of the GROMOS and AMBER force fields, including several methods aimed to increase the conformational sampling in computational molecular dynamics simulation trajectories. The system proved to be extremely sensitive to the employed force field, as well as to the presence of a water/air interface. In agreement with previous experiments and in contrast to the results obtained with AMBER, the analysis of the simulations using GROMOS showed a quick adsorption of the complex to the interface as well as an extremely exotic behavior of the water molecules surrounding the structure both in the bulk aqueous solution and at the water surface. The chirality of the CD molecule seems to play an important role in this behavior. All together, these results are expected to be useful to better understand the behavior of CD-based supramolecular complexes such as adsorption or aggregation driving forces, as well as to introduce new methods able to speed up general MD simulations.

Unexpected Cholesterol-Induced Destabilization of Lipid Membranes near Transmembrane Carbon Nanotubes

Cholesterol is a crucial component of mammalian cell membranes that takes part in many vital processes. It is generally accepted that cholesterol stabilizes the membrane and induces transitions into ordered states. In contrast to expectations, we demonstrate that cholesterol can destabilize the membrane by creating a nanodomain around a perpendicularly embedded ultrashort carbon nanotube (CNT), and we show that cholesterol triggers the translocation of an ultrashort CNT through the cell membrane. Using atomistic simulations, we report the existence of a nanoscale domain around an ultrashort carbon nanotube within a crossover distance of 0.9 nm from the surface of the nanotube, where the properties of the bilayer are different from the bulk: the domain is characterized by increased fluctuations, increased thickness, and increased order of the lipids with respect to the bulk. Cholesterol decreases the thickness and order of lipids and increases the fluctuations with respect to a pure lipid bilayer. Experimentally, we confirm that cholesterol nanodomains provoke spontaneous translocation of nanotubes through a lipid bilayer even for low membrane tensions. A specially designed microfluidic device allows us to trace the kinetic pathway of the translocation process and establish the threshold cholesterol concentration of 20% for translocation. The reported nanoscale cholesterol-induced membrane restructuring near the ultrashort CNT in lipid membranes enables precise control and specific targeting of a membrane using cholesterol. As an example, it may allow for specific targeting between cholesterol-rich mammalian cells and cholesterol-poor bacterial cells.

Mild lipid extraction and anisotropic cell membrane penetration of α-phase phosphorene carbide nanoribbons by molecular dynamics simulation studies

α-PC penetrates the interior of membrane efficiently only along its zigzag direction rather than its armchair direction.

DNA fragment translocation through the lipid membrane assisted by carbon nanotube

DNA delivery through cell membrane is a fundamental step for efficiency gene therapy. As a potential DNA carrier, carbon nanotubes (CNTs) have been extensively studied due to its unique properties. However, the mechanism of DNA translocation with CNTs through cell membrane is still not well understood. In this study, the DNA translocation process through POPC (1-palmitoyl-2-oleoylphosphatidylcholine) membrane with the assistance of CNTs was explored by molecular dynamics (MD) simulation. Our simulation results demonstrated that the CNTs could insert steadily into the POPC membrane, and DNA molecules tends to insert into the inner space of CNTs. With the assistance of CNTs, the free energy of nucleotides passing through the POPC membrane decreases. Moreover, the free energy of nucleotides (DA (deoxyadenosine), DT (deoxythymidine), DC (deoxycytidine), and DG (deoxyguanosine)) passing through POPC membrane follows the order: DA (deoxyadenosine) > DG (deoxyguanosine) > DC (deoxycytidine) > DT (deoxythymidine). These results may promote the design and application of CNT-based gene delivery system.

Molecular dynamics study of membrane permeabilization by wild-type and mutant lytic peptides from the non-enveloped Flock House virus

Flock House virus (FHV) serves as a model system for understanding infection mechanisms utilized by non-enveloped viruses to transport across cellular membranes. During the infection cycle of FHV, a fundamental stage involves disruption of the endosomal membrane by membrane active peptides, following externalization of the peptides from the capsid interior. The FHV lytic agents are the 44 C-terminal amino acids residues of the capsid protein, which are auto-catalytically cleaved during the capsid maturation process. The cleaved peptides are termed γ peptides. In this study, we perform multi-scale molecular dynamics simulations including 40 μs all-atom molecular dynamics simulations to study the behavior of pre-inserted transmembrane lytic peptides at a high concentration in a neutral membrane. We study the dynamical organization among peptides to form oligomeric bundles in four systems including the wild-type γ peptide and three mutant forms; namely, a truncation mutant in which the 23 C-terminal residues are deleted (γ₁), a construct where the 8 C-terminal residues of γ are fused to γ₁ (Δ385-399 γ) and a single-point mutant (F402A γ), all of which have been experimentally shown to drastically affect infectivity and lytic activity compared to the wild-type γ. Our results shed light on the actions of varied forms of the FHV lytic peptide including membrane insertion, trans-membrane stability, peptide oligomerization, water permeation activity and dynamic pore formation. Findings from this study provide detailed structural information and rationale for the differences in lytic activity among variants of FHV γ.

Conformational dynamics and interfacial interactions of peptide-appended pillar[5]arene water channels in biomimetic membranes

Peptide appended pillar[5]arene (PAP) is an artificial water channel resembling biological water channel proteins, which has shown a significant potential for designing bioinspired water purification systems. Given that PAP channels need to be incorporated at a high density in membrane matrices, it is critical to examine the role of channel-channel and channel-membrane interactions in governing the structural and functional characteristics of channels. To resolve the atomic-scale details of these interactions, we have carried out atomistic molecular dynamics (MD) simulations of multiple PAP channels inserted in a lipid or a block-copolymer (BCP) membrane matrix. Classical MD simulations on a sub-microsecond timescale showed clustering of channels only in the lipid membrane, but enhanced sampling MD simulations showed thermodynamically-favorable dimerized states of channels in both lipid and BCP membranes. The dimerized configurations of channels, with an extensive buried surface area, were stabilized via interactions between the aromatic groups in the peptide arms of neighboring channels. The conformational metrics characterizing the orientational and structural changes in channels revealed a higher flexibility in the lipid membrane as opposed to the BCP membrane although hydrogen bonds between the channel and the membrane molecules were not a major contributor to the stability of channels in the BCP membrane. We also found that the channels undergo wetting/dewetting transitions in both lipid and BCP membranes with a marginally higher probability of undergoing a dewetting transition in the BCP membrane. Collectively, these results highlight the role of channel dynamics in governing channel-channel and channel-membrane interfacial interactions, and provide atomic-scale insights needed to design stable and functional biomimetic membranes for efficient separations.

Helical Inclusions in Phospholipid Membranes: Lipid Adaptation and Chiral Order

The lipid bilayer is a flexible matrix that is able to adapt in response to the perturbation induced by inclusions, such as peptides and proteins. Here we use molecular dynamics simulations with a coarse-grained model to investigate the effect of a helical inclusion on a lipid bilayer in the liquid disordered phase. We show that the helical inclusion induces a collective tilt of acyl chains, with a small, yet unambiguous difference between a right- and a left-handed inclusion. This behavior is rationalized using the elastic continuum theory: The magnitude of the chiral (twist) deformation of the bilayer is determined by the interaction at the lipid/inclusion interface, and the decay length is controlled by the elastic properties of the bilayer. The lipid reorganization can thus be identified as a generic mechanism that, together with specific interactions, contributes to chiral recognition in phospholipid bilayers. An enhanced response is expected in highly ordered environments, such as rafts in biomembranes, with a potential impact on membrane-mediated interactions between inclusions.

Addressable Peptide Self-Assembly on the Cancer Cell Membrane for Sensitizing Chemotherapy of Renal Cell Carcinoma

Chemotherapy has been validated unavailable for treatment of renal cell carcinoma (RCC) in clinic due to its intrinsic drug resistance. Sensitization of chemo-drug response plays a crucial role in RCC treatment and increase of patient survival. Herein, a recognition-reaction-aggregation (RRA) cascaded strategy is utilized to in situ construct peptide-based superstructures on the renal cancer cell membrane, enabling specifically perturbing the permeability of cell membranes and enhancing chemo-drug sensitivity in vitro and in vivo. First, P1-DBCO can specifically recognize renal cancer cells by targeting carbonic anhydrase IX. Subsequently, P2-N₃ is introduced and efficiently reacts with P1-DBCO to form a peptide P3, which exhibits enhanced hydrophobicity and simultaneously aggregates into a superstructure. Interestingly, the superstructure retains on the cell membrane and perturbs its integrity/permeability, allowing more doxorubicin (DOX) uptaken by renal cancer cells. Owing to this increased influx, the IC₅₀ is significantly reduced by nearly 3.5-fold compared with that treated with free DOX. Finally, RRA strategy significantly inhibits the tumor growth of xenografted mice with a 3.2-fold enhanced inhibition rate compared with that treated with free DOX. In summary, this newly developed RRA strategy will open a new avenue for chemically engineering cell membranes with diverse biomedical applications.

Colorimetric in situ assay of membrane-bound enzyme based on lipid bilayer inhibition of ion transport

Membrane-bound enzymes (MBEs), which make up a very high proportion of intracellular enzymes, catalyze a variety of activities that are currently analyzed by various techniques after purification. However, due to their amphipathic character, the purification of MBEs is difficult. Therefore, the most productive approach represents in situ analysis of MBEs in the cellular membrane. Methods: In this study, using membrane-bound α-glucosidase (α-Glu) as an example, we have developed a colorimetric in situ assay for MBEs based on the inhibitory effect of lipid bilayer on ion transport. The enzyme substrate could mediate the self-assembly of phospholipid PEG derivative around magnetic nanospheres that were modified with boronic acid. The formation of lipid bilayer could inhibit the leaking of iron ions under acidic conditions. However, the product of the catalytic reaction had no capability for self-assembly of the lipid bilayer, leading to the release of iron ions from the magnetic nanospheres under acidic pH. Results: The colorimetric in situ assay for MBEs could not only analyze the activity of membrane-bound α-Glu located on Caco-2 cells but could also evaluate the α-Glu inhibitors in cell medium. Conclusions: The simple, economic, and efficient method proposed here offers a potential application for high-throughput testing of α-Glu and its inhibitors. Our study also outlines a strategy for exploring the colorimetric method to detect the activities of MBEs in situ.

Real-time dynamics of carbon nanotube porins in supported lipid membranes visualized by high-speed atomic force microscopy

In-plane mobility of proteins in lipid membranes is one of the fundamental mechanisms supporting biological functionality. Here we use high-speed atomic force microscopy (HS-AFM) to show that a novel type of biomimetic channel-carbon nanotube porins (CNTPs)-is also laterally mobile in supported lipid membranes, mimicking biological protein behaviour. HS-AFM can capture real-time dynamics of CNTP motion in the supported lipid bilayer membrane, build long-term trajectories of the CNTP motion and determine the diffusion coefficients associated with this motion. Our analysis shows that diffusion coefficients of CNTPs fall into the same range as those of proteins in supported lipid membranes. CNTPs in HS-AFM experiments often exhibit 'directed' diffusion behaviour, which is common for proteins in live cell membranes.This article is part of the themed issue 'Membrane pores: from structure and assembly, to medicine and technology'.

Communication between the leaflets of asymmetric membranes revealed from coarse-grain molecular dynamics simulations

We use coarse-grain molecular simulations to investigate the structural and dynamics differences between an asymmetric and a symmetrical membrane, both containing beta barrel transmembrane proteins. We find in where the dynamics of the two leaflets differ greatly, the slowest leaflet dominates the structural effects and importance of protein-lipid interactions.

Molecular dynamics simulations of carbon nanotube porins in lipid bilayers

Carbon nanotube porins embedded in lipid membranes are studied by molecular dynamics simulations.

Sub-nanometer pore formation in single-molecule-thick polyurea molecular-sieving membrane: a computational study

The surprisingly narrow sub-nm-pore-size distribution and urea-versus-glucose selectivity of a single-molecule-thick polyurea membrane are explained by Monte Carlo simulations on simple molecular models.

Ordering of lipid membranes altered by boron nitride nanosheets

Boron nitride nanosheets are novel promising nanomaterials with a lower cytotoxicity than graphene making them a better candidate for biomedical applications.

Insight of Transmembrane Processes of Self-Assembling Nanotubes Based on a Cyclic Peptide Using Coarse Grained Molecular Dynamics Simulation

Transmembrane self-assembling cyclic peptide (SCP) nanotubes are promising candidates for delivering specific molecules through cell membranes. The detailed mechanisms behind the transmembrane processes, as well as stabilization factors of transmembrane structures, are difficult to elucidate through experiments. In this study, the effects of peptide sequence and oligomeric state on the transmembrane capabilities of SCP nanotubes and the perturbation of embedded SCP nanotubes acting on the membrane were investigated based on coarse grained molecular dynamics simulation. The simulation results reveal that hydrophilic SCP oligomers result in the elevation of the energy barrier while the oligomerization of hydrophobic SCPs causes the reduction of the energy barrier, further leading to membrane insertion. Once SCP nanotubes are embedded, membrane properties such as density, thickness, ordering state and lateral mobility are adjusted along the radial direction. This study provides insight into the transmembrane strategy of SCP nanotubes and sheds light on designing novel transport systems.