Ligand Lipophilicity Determines Molecular Mechanisms of Nanoparticle Adsorption to Lipid Bilayers.

The interactions of ligand-functionalized nanoparticles with the cell membrane affect cellular uptake, cytotoxicity, and related behaviors, but relating these interactions to ligand properties remains challenging. In this work, we perform coarse-grained molecular dynamics simulations to study how the adsorption of ligand-functionalized cationic gold nanoparticles (NPs) to a single-component lipid bilayer (as a model cell membrane) is influenced by ligand end group lipophilicity. A set of 2 nm diameter NPs, each coated with a monolayer of organic ligands that differ only in their end groups, was simulated to mimic NPs recently studied experimentally. Metadynamics calculations were performed to determine key features of the free energy landscape for adsorption as a function of the distance of the NP from the bilayer and the number of NP-lipid contacts. These simulations revealed that NP adsorption is thermodynamically favorable for all NPs due to the extraction of lipids from the bilayer and into the NP monolayer. To resolve ligand-dependent differences in adsorption behavior, string method calculations were performed to compute minimum free energy pathways for adsorption. These calculations revealed a surprising nonmonotonic dependence of the free energy barrier for adsorption on ligand end group lipophilicity. Large free energy barriers are predicted for the least lipophilic end groups because favorable NP-lipid contacts are initiated only through the unfavorable protrusion of lipid tail groups out of the bilayer. The smallest free energy barriers are predicted for end groups of intermediate lipophilicity which promote NP-lipid contacts by intercalating within the bilayer. Unexpectedly, large free energy barriers are also predicted for the most lipophilic end groups which remain sequestered within the ligand monolayer rather than intercalating within the bilayer. These trends are broadly in agreement with past experimental measurements and reveal how subtle variations in ligand lipophilicity dictate adsorption mechanisms and associated kinetics by influencing the interplay of lipid-ligand interactions.

Bacterial Quorum Sensing Signals Promote Large-Scale Remodeling of Lipid Membranes.

We report that N-acyl-l-homoserine lactones (AHLs), a class of nonionic amphiphiles that common bacteria use as signals to coordinate group behaviors, can promote large-scale remodeling in model lipid membranes. Characterization of supported lipid bilayers (SLBs) of the phospholipid 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) by fluorescence microscopy and quartz crystal microbalance with dissipation (QCM-D) reveals the well-studied AHL signal 3-oxo-C12-AHL and its anionic head group hydrolysis product (3-oxo-C12-HS) to promote the formation of long microtubules that can retract into hemispherical caps on the surface of the bilayer. These transformations are dynamic, reversible, and dependent upon the head group structure. Additional experiments demonstrate that 3-oxo-C12-AHL can promote remodeling to form microtubules in lipid vesicles and promote molecular transport across bilayers. Molecular dynamics (MD) simulations predict differences in thermodynamic barriers to translocation of these amphiphiles across a bilayer that are reflected in both the type and extent of reformation and associated dynamics. Our experimental observations can thus be interpreted in terms of accumulation and relief of asymmetric stresses in the inner and outer leaflets of a bilayer upon intercalation and translocation of these amphiphiles. Finally, experiments on Pseudomonas aeruginosa, a pathogen that uses 3-oxo-C12-AHL for cell-to-cell signaling, demonstrate that 3-oxo-C12-AHL and 3-oxo-C12-HS can promote membrane remodeling at biologically relevant concentrations and in the absence of other biosurfactants, such as rhamnolipids, that are produced at high population densities. Overall, these results have implications for the roles that 3-oxo-C12-AHL and its hydrolysis product may play in not only mediating intraspecies bacterial communication but also processes such as interspecies signaling and bacterial control of host-cell response. Our findings also provide guidance that could prove useful for the design of synthetic self-assembled materials that respond to bacteria in ways that are useful in the context of sensing, drug delivery, and in other fundamental and applied areas.

Analysis of Charged Peptide Loop-Flipping across a Lipid Bilayer Using the String Method with Swarms of Trajectories.

The hydrophobic core of the lipid bilayer is conventionally considered a barrier to the translocation of charged species such that the translocation of even single ions occurs on long time scales. In contrast, experiments have revealed that some materials, including peptides, proteins, and nanoparticles, can translocate multiple charged moieties across the bilayer on experimentally relevant time scales. Understanding the molecular mechanisms underlying this behavior is challenging because resolving corresponding free-energy landscapes with molecular simulation techniques is computationally expensive. To address this challenge, we use atomistic molecular dynamics simulations with the swarms-of-trajectories (SOT) string method to analyze charge translocation pathways across single-component lipid bilayers as a function of multiple collective variables. We first demonstrate that the SOT string method can reproduce the free-energy barrier for the translocation of a charged lysine amino acid analogue in good agreement with the literature. We then obtain minimum free-energy pathways for the translocation, or flipping, of charged peptide loops across the lipid bilayer by utilizing trajectories from prior temperature-accelerated molecular dynamics (TAMD) simulations as initial configurations. The corresponding potential of mean force calculations reveal that the protonation of a central membrane-exposed aspartate residue substantially reduces the free-energy barrier for flipping charged loops by modulating the water content of the bilayer. These results provide new insight into the thermodynamics underlying loop-flipping processes and highlight how the combination of TAMD and the SOT string method can be used to understand complex charge translocation mechanisms.

Molecular simulations of lipid membrane partitioning and translocation by bacterial quorum sensing modulators.

Quorum sensing (QS) is a bacterial communication process mediated by both native and non-native small-molecule quorum sensing modulators (QSMs), many of which have been synthesized to disrupt QS pathways. While structure-activity relationships have been developed to relate QSM structure to the activation or inhibition of QS receptors, less is known about the transport mechanisms that enable QSMs to cross the lipid membrane and access intracellular receptors. In this study, we used atomistic MD simulations and an implicit solvent model, called COSMOmic, to analyze the partitioning and translocation of QSMs across lipid bilayers. We performed umbrella sampling at atomistic resolution to calculate partitioning and translocation free energies for a set of naturally occurring QSMs, then used COSMOmic to screen the water-membrane partition and translocation free energies for 50 native and non-native QSMs that target LasR, one of the LuxR family of quorum-sensing receptors. This screening procedure revealed the influence of systematic changes to head and tail group structures on membrane partitioning and translocation free energies at a significantly reduced computational cost compared to atomistic MD simulations. Comparisons with previously determined QSM activities suggest that QSMs that are least likely to partition into the bilayer are also less active. This work thus demonstrates the ability of the computational protocol to interrogate QSM-bilayer interactions which may help guide the design of new QSMs with engineered membrane interactions.

Bacterial Quorum Sensing Signals Self-Assemble in Aqueous Media to Form Micelles and Vesicles: An Integrated Experimental and Molecular Dynamics Study.

Many species of common bacteria communicate and coordinate group behaviors, including toxin production and surface fouling, through a process known as quorum sensing (QS). In Gram-negative bacteria, QS is regulated by N-acyl-l-homoserine lactones (AHLs) that possess a polar homoserine lactone headgroup and a nonpolar aliphatic tail. Past studies demonstrate that AHLs can aggregate in water or adsorb at interfaces, suggesting that molecular self-assembly could play a role in processes that govern bacterial communication. We used a combination of biophysical characterization and atomistic molecular dynamics (MD) simulations to characterize the self-assembly behaviors of 12 structurally related AHLs. We used static light scattering and measurements of surface tension to characterize the assembly of four naturally occurring AHLs (3-oxo-C8-AHL, 3-oxo-C12-AHL, C12-AHL, and C16-AHL) in aqueous media and determine their critical aggregation concentrations (CACs). MD simulations and alchemical free energy calculations were used to predict thermodynamically preferred aggregate structures for each AHL. Those calculations predicted that AHLs with 10 or 12 tail carbon atoms should form spherical micelles and that AHLs with 14 or 16 tail carbon atoms should form vesicles in solution. Characterization of solutions of AHLs using negative stain transmission electron microscopy (TEM) and dynamic light scattering (DLS) revealed aggregates with sizes consistent with spherical micelles or small unilamellar vesicles for 3-oxo-C12-AHL and C12-AHL and the formation of large vesicles (∼250 nm) in solutions of C16-AHL. These experimental findings are in general agreement with our simulation predictions. Overall, our results provide insight into processes of self-assembly that can occur in this class of bacterial amphiphiles and, more broadly, provide a potential basis for understanding how AHL structure could influence processes that bacteria use to drive important group behaviors.

Characterizing the Molecular Mechanisms for Flipping Charged Peptide Flanking Loops across a Lipid Bilayer.

The cell membrane largely prevents the passive diffusion of charged molecules due to the large free energy barrier associated with translocating charged groups across the hydrophobic lipid bilayer core. Despite this barrier, some peptides can interconvert between transmembrane and surface-adsorbed states by "flipping" charged flanking loops across the bilayer on a surprisingly rapid second-minute time scale. The transmembrane helices of some multispanning membrane proteins undergo similar reorientation processes, suggesting that loop-flipping may be a mechanism for regulating membrane protein topology; however, the molecular mechanisms underlying this behavior remain unknown. In this work, we study the loop-flipping behavior exhibited by a peptide with a hydrophobic transmembrane helix, charged flanking loops, and a central, membrane-exposed aspartate residue of varying protonation state. We utilize all-atom temperature accelerated molecular dynamics simulations to predict the likelihood of loop-flipping without predefining specific loop-flipping pathways. We demonstrate that this approach can identify multiple possible flipping pathways, with the prevalence of each pathway depending on the protonation state of the central residue. In particular, we find that a charged central residue facilitates loop-flipping by stabilizing membrane water defects, enabling the "self-catalysis" of charge translocation. These findings provide detailed molecular-level insights into charged loop-flipping pathways that may generalize to other charge translocation processes, such as lipid flip-flop or the large-scale conformational rearrangements of multispanning membrane proteins.