Structure of the Bacterial Cell Envelope and Interactions with Antimicrobials: Insights from Molecular Dynamics Simulations.

Bacteria have evolved over 3 billion years, shaping our intrinsic and symbiotic coexistence with these single-celled organisms. With rising populations of drug-resistant strains, the search for novel antimicrobials is an ongoing area of research. Advances in high-performance computing platforms have led to a variety of molecular dynamics simulation strategies to study the interactions of antimicrobial molecules with different compartments of the bacterial cell envelope of both Gram-positive and Gram-negative species. In this review, we begin with a detailed description of the structural aspects of the bacterial cell envelope. Simulations concerned with the transport and associated free energy of small molecules and ions through the outer membrane, peptidoglycan, inner membrane and outer membrane porins are discussed. Since surfactants are widely used as antimicrobials, a section is devoted to the interactions of surfactants with the cell wall and inner membranes. The review ends with a discussion on antimicrobial peptides and the insights gained from the molecular simulations on the free energy of translocation. Challenges involved in developing accurate molecular models and coarse-grained strategies that provide a trade-off between atomic details with a gain in sampling time are highlighted. The need for efficient sampling strategies to obtain accurate free energies of translocation is also discussed. Molecular dynamics simulations have evolved as a powerful tool that can potentially be used to design and develop novel antimicrobials and strategies to effectively treat bacterial infections.

Martini-3 Coarse-Grained Models for the Bacterial Lipopolysaccharide Outer Membrane of Escherichia coli.

The outer lipopolysaccharide (LPS) membrane of Gram-negative bacteria forms the main barrier for transport of antimicrobial molecules into the bacterial cell. In this study we develop coarse-grained models for the outer membrane of Escherichia coli in the Martini-3 framework. The coarse-grained model force field was parametrized and validated using all-atom simulations of symmetric membranes of lipid A and rough LPS as well as a complete asymmetric membrane of LPS with the O-antigen. The bonded parameters were obtained using an iterative refinement procedure with target bonded distributions obtained from all-atom simulations. The membrane thickness, area of the LPS, and density distributions for the different regions as well as the water and ion densities in Martini-3 simulations show excellent agreement with the all-atom data. Additionally the solvent accessible surface area for individual molecules in water was found to be in good agreement. The binding of calcium ions with phosphate and carboxylate moieties of LPS is accurately captured in the Martini-3 model, indicative of the integrity of the highly negatively charged LPS molecules in the outer membranes of Gram-negative bacteria. The melting transition of the coarse-grained lipid A membrane model was found to occur between 300 and 310 K, and the model captured variations in area per LPS, order parameter, and membrane thickness across the melting transition. Our study reveals that the proposed Martini-3 models for LPS are able to capture the physicochemical balance of the complex sugar architecture of the outer membrane of Escherichia coli. The coarse-grained models developed in this study would be useful for determining membrane protein interactions and permeation of potential antimicrobials through bacterial membranes at mesoscopic spatial and temporal scales.

Trajectory-Extending Kinetic Monte Carlo Simulations to Evaluate Pure and Gas Mixture Diffusivities through a Dense Polymeric Membrane.

With renewed interest in CO2 separations, carbon molecular sieving (CMS) membrane performance evaluation requires diffusion coefficients as inputs to have a reliable estimate of the permeability. An optimal material is desired to have both high selectivity and permeability. Gases diffusing through dense CMS and polymeric membranes experience extended subdiffusive regimes, which hinders reliable extraction of diffusion coefficients from mean squared displacement data. We improve the sampling of the diffusive landscape by implementing the trajectory-extending kinetic Monte Carlo (TEKMC) technique to efficiently extend molecular dynamics (MD) trajectories from ns to µs time scales. The obtained self-diffusion coefficient of pure CO2 in CMS membranes derived from a 6FDA/BPDA-DAM precursor polymer melt is found to linearly increase from 0.8-1.3 × 10-6 cm2 s-1 in the pressure range of 1-20 bar, which supports previous experimental findings. We also extended the TEKMC algorithm to evaluate the mixture diffusivities in binary mixtures to determine the permselectivity of CO2 in CH4 and N2 mixtures. The mixture diffusion coefficient of CO2 ranges from 1.3-7 × 10-6 cm2 s-1 in the binary mixture CO2/CH4, which is significantly higher than the pure gas diffusion coefficient. Robeson plot comparisons show that the permselectivity obtained from pure gas diffusion data is significantly lower than that predicted using mixture diffusivity data. Specifically, in the case of the CO2/N2 mixture, we find that using mixture diffusivities led to permselectivities lying above the Robeson limit highlighting the importance of using mixture diffusivity data for an accurate evaluation of the membrane performance. Combined with gas solubilities obtained from grand canonical Monte Carlo (GCMC) simulations, our work shows that simulations with the TEKMC method can be used to reliably evaluate the performance of materials for gas separations.

Selective Photonic Gasification of Strained Oxygen Clusters on Graphene for Tuning Pore Size in the Å Regime.

Controlling the size of single-digit pores, such as those in graphene, with an Å resolution has been challenging due to the limited understanding of pore evolution at the atomic scale. The controlled oxidation of graphene has led to Å-scale pores; however, obtaining a fine control over pore evolution from the pore precursor (i.e., the oxygen cluster) is very attractive. Herein, we introduce a novel "control knob" for gasifying clusters to form pores. We show that the cluster evolves into a core/shell structure composed of an epoxy group surrounding an ether core in a bid to reduce the lattice strain at the cluster core. We then selectively gasified the strained core by exposing it to 3.2 eV of light at room temperature. This allowed for pore formation with improved control compared to thermal gasification. This is because, for the latter, cluster-cluster coalescence via thermally promoted epoxy diffusion cannot be ruled out. Using the oxidation temperature as a control knob, we were able to systematically increase the pore density while maintaining a narrow size distribution. This allowed us to increase H2 permeance as well as H2 selectivity. We further show that these pores could differentiate CH4 from N2, which is considered to be a challenging separation. Dedicated molecular dynamics simulations and potential of mean force calculations revealed that the free energy barrier for CH4 translocation through the pores was lower than that for N2. Overall, this study will inspire research on the controlled manipulation of clusters for improved precision in incorporating Å-scale pores in graphene.

Cholesterol catalyzes unfolding in membrane-inserted motifs of the pore forming protein cytolysin A.

Plasma membrane-induced protein folding and conformational transitions play a central role in cellular homeostasis. Several transmembrane proteins are folded in the complex lipid milieu to acquire a specific structure and function. Bacterial pore forming toxins (PFTs) are proteins expressed by a large class of pathogenic bacteria that exploit the plasma membrane environment to efficiently undergo secondary structure changes, oligomerize, and form transmembrane pores. Unregulated pore formation causes ion imbalance, leading to cell death and infection. Determining the free energy landscape of these membrane-driven-driven transitions remains a challenging problem. Although cholesterol recognition is required for lytic activity of several proteins in the PFT family of toxins, the regulatory role of cholesterol for the α-PFT, cytolysin A expressed by Escherichia coli remains unexplained. In a recent free energy computation, we showed that the ß tongue, a critical membrane-inserted motif of the ClyA toxin, has an on-pathway partially unfolded intermediate that refolds into the helix-turn-helix motif of the pore state. To understand the molecular role played by cholesterol, we carry out string-method-based computations in membranes devoid of cholesterol, which reveals an increase of ∼30 times in the free energy barrier for the loss of ß sheet secondary structure when compared with membranes containing cholesterol. Specifically, the tyrosine-cholesterol interaction was found to be critical to creating the unfolded intermediate. Cholesterol also increases the packing and hydrophobicity of the bilayer, resulting in enhanced interactions of the bound protein before complete membrane insertion. Our study illustrates that cholesterol is critical to catalyzing and stabilizing the membrane-inserted unfolded state of the ß tongue motif of ClyA, opening up fresh insights into cholesterol-assisted unfolding of membrane proteins.

Graphene Nanopores Enhance Water Evaporation from Salt Solutions: Exploring the Effects of Ions and Concentration.

With increased water stress, the development of clean water technologies is an active area of research. Evaporation-based solutions offer the advantage of low energy consumption, and recently a 10-30 fold enhancement in water evaporation flux has been observed through Å-scale graphene nanopores (Lee, W.-C., et al., ACS Nano 2022, 16(9), 15382). Herein, using molecular dynamics simulations, we examine the suitability of Å-scale graphene nanopores in enhancing water evaporation from salt solutions (LiCl, NaCl, and KCl). Cation-π interactions between ions and the surface of nanoporous graphene are found to significantly influence ion populations in the nanopore vicinity, leading to varied water evaporation fluxes from different salt solutions. The highest water evaporation flux was observed for KCl solutions, followed by NaCl and LiCl solutions, with the differences reducing at lower concentrations. Relative to the bare liquid-vapor interface, 4.54 Å nanopores exhibit the highest evaporation flux enhancements ranging from 7 to 11, with an enhancement of 10.8 obtained for 0.6 M NaCl solution, which closely resembles seawater compositions. Functionalized nanopores induce short-lived water-water hydrogen bonds and reduce surface tension at the liquid-vapor interface, thereby lowering the free energy barrier for water evaporation with a negligible effect on the ion hydration dynamics. These findings can aid in developing green technologies for desalination and separation processes with low thermal energy input.

Conformational Flexibility Is a Key Determinant for the Lytic Activity of the Pore-Forming Protein, Cytolysin A.

Several bacterial infections are mediated by pore-forming toxins (PFTs), a subclass of proteins that oligomerize on mammalian cell membranes forming lytic nanopores. Cytolysin A (ClyA), an α-PFT, undergoes a dramatic conformational change restructuring its two membrane-binding motifs (the ß-tongue and the N-terminus helix), during pore formation. A complete molecular picture for this key transition and the driving force behind the secondary structure change upon membrane binding remain elusive. Using all-atom molecular dynamics (MD) simulations of the ClyA monomer and string method based free energy computations with path collective variables, we illustrate that an unfolded ß-tongue motif is an on-pathway intermediate during the transition to the helix-turn-helix motif of the protomer. An aggregate of 28 µs of all-atom thermal unfolding MD simulations of wild-type ClyA and its single point mutants reveal that the membrane-binding motifs of the ClyA protein display high structural flexibility in water. However, point mutations in these motifs lead to a distinct reduction in the flexibility, especially in the ß-tongue, thereby stabilizing the pretransition secondary structure. Resistance to unfolding was further corroborated by MD simulations of the ß-tongue mutant motif in the membrane. Combined with the thermal unfolding simulations, we posit that the ß-tongue as well as N-terminal mutants that lower the tendency to unfold and disorder the ß-tongue are detrimental to pore formation by ClyA and its lytic activity. Erythrocyte turbidity and vesicle leakage assays indeed reveal a loss of activity for the ß-tongue mutant, and delayed kinetics for the N-terminus mutants. On the other hand, a point mutation in the extracellular domain that did not abrogate lytic activity displayed similar unfolding characteristics as the wild type. Thus, attenuation of conformational flexibility in membrane-binding motifs correlates with reduced lytic and leakage activity. Combined with secondary structure changes observed in the membrane bound states, our study shows that the tendency to unfold in the ß-tongue region is a critical step in the conformational transition and bistability of the ClyA protein and mutants that disrupt this tendency reduced pore formation. Overall, our finding suggests that inherent flexibility in the protein could play a wider and hitherto unrecognized role in membrane-mediated conformational transitions of PFTs and other membrane protein transformations.

Interactions of Surfactants with the Bacterial Cell Wall and Inner Membrane: Revealing the Link between Aggregation and Antimicrobial Activity.

Surfactants with their intrinsic ability to solubilize lipid membranes are widely used as antibacterial agents, and their interactions with the bacterial cell envelope are complicated by their differential aggregation tendencies. We present a combined experimental and molecular dynamics investigation to unravel the molecular basis for the superior antimicrobial activity and faster kill kinetics of shorter-chain fatty acid surfactant, laurate, when compared with the longer-chain surfactants studied in contact time assays with live Escherichia coli (E. coli). From all-atom molecular dynamics simulations, translocation events across peptidoglycan were the highest for laurate followed by sodium dodecyl sulfate, myristate, palmitate, oleate, and stearate. The translocation kinetics were positively correlated with the critical micellar concentration, which determined the free monomer surfactant concentration available for translocation across peptidoglycan. Interestingly, aggregates showed a lower propensity to translocate across the peptidoglycan layer and longer translocation times were observed for oleate, thereby revealing an intrinsic sieving property of the bacterial cell wall. Molecular dynamics simulations with surfactant-incorporated bacterial inner membranes revealed the greatest hydrophobic mismatch and membrane thinning in the presence of laurate when compared with the other surfactants. The enhanced antimicrobial efficacy of laurate over oleate was further verified by experiments with giant unilamellar vesicles, and electroporation molecular dynamics simulations revealed greater inner membrane poration tendency in the presence of laurate when compared with the longer-chain surfactants. Our study provides molecular insights into surfactant translocation across peptidoglycan and chain length-induced structural disruption of the inner membrane, which correlate with contact time kill efficacies observed as a function of chain length with E. coli. The insights gained from our study uncover unexplored barrier properties of the bacterial cell envelope to rationalize the development of antimicrobial formulations and therapeutics.

Finite temperature string method with umbrella sampling using path collective variables: application to secondary structure change in a protein.

The transition of an α-helix to a ß-sheet in proteins is among the most complex conformational changes seen in biomolecular systems. Due to long time scales involved in the transition, it is challenging to study such protein conformational changes using direct molecular dynamics simulations. This limitation is typically overcome using an indirect approach wherein one computes the free energy landscape associated with the transition. Computation of free energy landscapes, however, requires a suitable set of collective variables that describe the transition. In this work, we demonstrate the use of path collective variables [D. Branduardi, F. L. Gervasio and M. Parrinello, J. Chem. Phys., 2007, 126, 054103] and combine it with the finite temperature string (FTS) method [E. Weinan, W. Ren and E. Vanden-Eijnden, J. Phys. Chem. B, 2005, 109, 6688-6693] to determine the molecular mechanisms involved during the structural transition of the mini G-protein from an α-helix to a ß-hairpin. The transition from the α-helix proceeds via unfolding of the terminal residues, giving rise to a ß-turn unfolded intermediate to eventually form the ß-hairpin. Our proposed algorithm uses umbrella sampling simulations to simulate images along the string and the weighted histogram analysis to compute the free energy along the computed transition path. This work demonstrates that the string method in combination with path collective variables can be exploited to study complex protein conformational changes such as a complete change in the secondary structure.

Enhanced Water Evaporation from Å-Scale Graphene Nanopores.

Enhancing the kinetics of liquid-vapor transition from nanoscale confinements is an attractive strategy for developing evaporation and separation applications. The ultimate limit of confinement for evaporation is an atom thick interface hosting angstrom-scale nanopores. Herein, using a combined experimental/computational approach, we report highly enhanced water evaporation rates when angstrom sized oxygen-functionalized graphene nanopores are placed at the liquid-vapor interface. The evaporation flux increases for the smaller nanopores with an enhancement up to 35-fold with respect to the bare liquid-vapor interface. Molecular dynamics simulations reveal that oxygen-functionalized nanopores render rapid rotational and translational dynamics to the water molecules due to a reduced and short-lived water-water hydrogen bonding. The potential of mean force (PMF) reveals that the free energy barrier for water evaporation decreases in the presence of nanopores at the atomically thin interface, which further explains the enhancement in evaporation flux. These findings can enable the development of energy-efficient technologies relying on water evaporation.

A Molecular Dynamics Study of Antimicrobial Peptide Interactions with the Lipopolysaccharides of the Outer Bacterial Membrane.

With rising bacterial resistance, antimicrobial peptides (AMPs) have been widely investigated as potential antibacterial molecules to replace conventional antibiotics. Our understanding of the molecular mechanisms for membrane disruption are largely based on AMP interactions with the inner phospholipid bilayers of both Gram-negative and Gram-positive bacteria. Mechanisms for AMP translocation across the outer membrane of Gram-negative bacteria composed of lipopolysaccharides and the asymmetric lipid bilayer are complicated by the secondary structure adopted by the peptide in the different membrane environments. We have employed atomistic molecular dynamics and umbrella-sampling simulations with an aggregate duration of [Formula: see text] 6 microseconds to obtain the free energy landscape of CM15 peptide translocating through the lipopolysaccharide region of Gram-negative bacteria, E. coli. The peptide has a favorable binding-free energy (- 130 kJ mol[Formula: see text]) in the O-antigen region with a large barrier (150 kJ mol[Formula: see text]) at the interface between the anionic core saccharides and upper bilayer leaflet made up of lipid-A molecules. Restraint-free molecular dynamics simulations show that the random coil structure is favored over the helix in both the extracellular aqueous region and the cation-rich core-saccharide regions of the outer membrane. The peptide and membrane properties are analyzed at each of the 100 ns duration of the umbrella-sampling windows to illustrate changes in peptide length, orientation, and hydration. Our study provides insights into the free energy landscape for the insertion of the AMP CM15 in the outer membrane of Gram-negative bacteria, and we discuss the implications of our findings with the broader question of how AMPs overcome this barrier during antimicrobial activity.

Influence of the extent of hydrophobicity on water organization and dynamics on 2D graphene oxide surfaces.

Graphene oxide (GO) nanomaterials are being extensively explored for a wide spectrum of applications, ranging from water desalination to fuel cell applications, due to their tunable mechanical, thermal, and electrical properties. In this paper, we have investigated the influence of the hydrophobic extent on the adsorption of water on 2D GO surfaces by performing a series of grand canonical Monte Carlo simulations at various relative pressures, P/P0, at 298 K and discuss the implications of our findings on proton transport characteristics. HR is defined as the ratio of the hydrophobic to hydrophilic areas on the GO surface. The structure of adsorbed water is studied by analyzing density distributions and hydrogen bonds. At moderate relative pressures of P/P0 < 0.6, a monolayer of adsorbed water, spanning the hydrophilic and hydrophobic regions of the GO surface, is observed for HR = 0, 0.5 and 1, and at higher pressures, a percolating hydrogen-bonded network is formed, which results in the formation of a thick water film. At intermediate water pressures, bridging water networks form across the hydrophobic regions. The GO surface of HR = 1 is seen to have a strong signature of a Janus surface, displaying increased fluctuations in adsorbed water molecules and hydrogen bonds. Our results suggest that if there is sufficient hydrophilicity on the GO surface, a relative humidity between 70 and 80% results in the formation of a fully formed contact water layer hydrogen-bonded with the surface functional groups along with a second layer of adsorbed water molecules. This coincides with hydration levels at which a maximum in the proton conductivity has been reported on 2D GO surfaces. Molecular dynamics simulations reveal a higher reorientational relaxation time at lower water hydration and the rotational entropy of interfacial water at lower hydration is higher than that of bulk water, indicating broader rotational phase space sampling.

Coupling of mitochondrial population evolution to microtubule dynamics in fission yeast cells: a kinetic Monte Carlo study.

Mitochondrial populations in cells are maintained by cycles of fission and fusion events. Perturbation of this balance has been observed in several diseases such as cancer and neurodegeneration. In fission yeast cells, the association of mitochondria with microtubules inhibits mitochondrial fission [Mehta et al., J. Biol. Chem., 2019, 294, 3385], illustrating the intricate coupling between mitochondria and the dynamic population of microtubules within the cell. In order to understand this coupling, we carried out kinetic Monte Carlo (KMC) simulations to predict the evolution of mitochondrial size distributions for different cases; wild-type cells, cells with short and long microtubules, and cells without microtubules. Comparisons are made with mitochondrial distributions reported in experiments with fission yeast cells. Using experimentally determined mitochondrial fission and fusion frequencies, simulations implemented without the coupling of microtubule dynamics predicted an increase in the mean number of mitochondria, equilibrating within 50 s. The mitochondrial length distribution in these models also showed a higher occurrence of shorter mitochondria, implying a greater tendency for fission, similar to the scenario observed in the absence of microtubules and cells with short microtubules. Interestingly, this resulted in overestimating the mean number of mitochondria and underestimating mitochondrial lengths in cells with wild-type and long microtubules. However, coupling mitochondria's fission and fusion events to the microtubule dynamics effectively captured the mitochondrial number and size distributions in wild-type and cells with long microtubules. Thus, the model provides greater physical insight into the temporal evolution of mitochondrial populations in different microtubule environments, allowing one to study both the short-time evolution as observed in the experiments (<5 minutes) as well as their transition towards a steady-state (>15 minutes). Our study illustrates the critical role of microtubules in mitochondrial dynamics and coupling microtubule growth and shrinkage dynamics is critical to predicting the evolution of mitochondrial populations within the cell.

A generic force field for simulating native protein structures using dissipative particle dynamics.

A coarse-grained force field for molecular dynamics simulations of the native structures of proteins in a dissipative particle dynamics (DPD) framework is developed. The parameters for bonded interactions are derived by mapping the bonds and angles for 20 amino acids onto target distributions obtained from fully atomistic simulations in explicit solvent. A dual-basin potential is introduced for stabilizing backbone angles, to cover a wide spectrum of protein secondary structures. The backbone dihedral potential enables folding of the protein from an unfolded initial state to the folded native structure. The proposed force field is validated by evaluating the structural properties of several model peptides and proteins including the SARS-CoV-2 fusion peptide, consisting of α-helices, ß-sheets, loops and turns. Detailed comparisons with fully atomistic simulations are carried out to assess the ability of the proposed force field to stabilize the different secondary structures present in proteins. The compact conformations of the native states were evident from the radius of gyration and the high intensity peaks of the root mean square deviation histograms, which were found to be within 0.4 nm. The Ramachandran-like energy landscape on the phase space of backbone angles (θ) and dihedrals (ϕ) effectively captured the conformational phase space of α-helices at ∼(ϕ = 50°,θ = 90°) and ß-strands at ∼(ϕ = ±180°,θ = 90-120°). Furthermore, the residue-residue native contacts were also well reproduced by the proposed DPD model. The applicability of the model to multidomain complexes was assessed using lysozyme and a large α-helical bacterial pore-forming toxin, cytolysin A. Our study illustrates that the proposed force field is generic, and can potentially be extended for efficient in silico investigations of membrane bound polypeptides and proteins using DPD simulations.

Pore Forming Protein Induced Biomembrane Reorganization and Dynamics: A Focused Review.

Pore forming proteins are a broad class of pathogenic proteins secreted by organisms as virulence factors due to their ability to form pores on the target cell membrane. Bacterial pore forming toxins (PFTs) belong to a subclass of pore forming proteins widely implicated in bacterial infections. Although the action of PFTs on target cells have been widely investigated, the underlying membrane response of lipids during membrane binding and pore formation has received less attention. With the advent of superresolution microscopy as well as the ability to carry out molecular dynamics (MD) simulations of the large protein membrane assemblies, novel microscopic insights on the pore forming mechanism have emerged over the last decade. In this review, we focus primarily on results collated in our laboratory which probe dynamic lipid reorganization induced in the plasma membrane during various stages of pore formation by two archetypal bacterial PFTs, cytolysin A (ClyA), an α-toxin and listeriolysin O (LLO), a ß-toxin. The extent of lipid perturbation is dependent on both the secondary structure of the membrane inserted motifs of pore complex as well as the topological variations of the pore complex. Using confocal and superresolution stimulated emission depletion (STED) fluorescence correlation spectroscopy (FCS) and MD simulations, lipid diffusion, cholesterol reorganization and deviations from Brownian diffusion are correlated with the oligomeric state of the membrane bound protein as well as the underlying membrane composition. Deviations from free diffusion are typically observed at length scales below ∼130 nm to reveal the presence of local dynamical heterogeneities that emerge at the nanoscale-driven in part by preferential protein binding to cholesterol and domains present in the lipid membrane. Interrogating the lipid dynamics at the nanoscale allows us further differentiate between binding and pore formation of ß- and α-PFTs to specific domains in the membrane. The molecular insights gained from the intricate coupling that occurs between proteins and membrane lipids and receptors during pore formation are expected to improve our understanding of the virulent action of PFTs.

Bacterial protein listeriolysin O induces nonmonotonic dynamics because of lipid ejection and crowding.

Membrane-bound protein complexes involving pore forming toxins (PFTs) released by virulent bacteria are known to form transmembrane pores leading to host cell lysis. Developing alternative strategies against PFT mediated bacterial virulence factors requires an understanding of the cellular membrane response. However, membrane disruption and related lipid reorganization events during attack by PFTs remain largely unexplored. We report counterintuitive and nonmonotonic variations in lipid diffusion, measured using confocal fluorescence correlation spectroscopy, due to interplay of lipid ejection and crowding by membrane-bound oligomers of a prototypical cholesterol-dependent cytolysin, listeriolysin O (LLO). The observed dynamical crossover is correlated with concentration dependent transitions of LLO oligomeric state populations from rings to arc-like pore complexes, predicted using a proposed two-state free area-based diffusion model. At low PFT concentrations, a hitherto unexplored regime of increased lipid diffusivity is attributed to lipid ejection events because of a preponderance of ring-like pore states. At higher protein concentrations in which membrane-inserted arc-like pores dominate, lipid ejection is less efficient and the ensuing crowding results in a lowering of lipid diffusion. These variations in lipid dynamics are corroborated by macroscopic rheological response measurements of PFT bound vesicles. Our study correlates PFT oligomeric state transitions, membrane remodeling, and mechanical property variations, providing unique insights into the pore forming mechanisms of cholesterol-dependent cytolysins.

Evaluating Coarse-Grained MARTINI Force-Fields for Capturing the Ripple Phase of Lipid Membranes.

Phospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (Pß') membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained (CG) resolutions. Here, with an aggregate ∼2.5 µs all-atom and ∼122 µs CGMD simulations, we systematically assess the ability of six CG MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid and water force-field (FF) variants, parametrized to capture the DPPC gel and fluid phases, for their ability to capture the Pß' phase, and compared observations with those from an all-atom FF. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (>2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated Pß' and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the Pß' phase, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically trapped structures caused by finite size effects rather than being representative of true Pß' phases. We showed that a MARTINI FF variant that could capture the tilted Lß' gel phase, a prerequisite for stabilizing the Pß' phase, was unable to capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs (including MARTINI3) may require specific reparametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase.

Hydrophobic Gating and 1/f Noise of the Anthrax Toxin Channel.

"Pink" or 1/f noise is a natural phenomenon omnipresent in physics, economics, astrophysics, biology, and even music and languages. In electrophysiology, the stochastic activity of a number of biological ion channels and artificial nanopores could be characterized by current noise with a 1/f power spectral density. In the anthrax toxin channel (PA63), it appears as fast voltage-independent current interruptions between conducting and nonconducting states. This behavior hampers potential development of PA63 as an ion-channel biosensor. On the bright side, the PA63 flickering represents a mesmerizing phenomenon to investigate. Notably, similar 1/f fluctuations are observed in the channel-forming components of clostridial binary C2 and iota toxins, which share functional and structural similarities with the anthrax toxin channel. Similar to PA63, they are evolved to translocate the enzymatic components of the toxins into the cytosol. Here, using high-resolution single-channel lipid bilayer experiments and all-atom molecular dynamic simulations, we suggest that the 1/f noise in PA63 occurs as a result of "hydrophobic gating" at the ϕ-clamp region, the phenomenon earlier observed in several water-filled channels "fastened" inside by the hydrophobic belts. The ϕ-clamp is a narrow "hydrophobic ring" in the PA63 lumen formed by seven or eight phenylalanine residues at position 427, conserved in the C2 and iota toxin channels, which catalyzes protein translocation. Notably, the 1/f noise remains undetected in the F427A PA63 mutant. This finding can elucidate the functional purpose of 1/f noise and its possible role in the transport of the enzymatic components of binary toxins.

Using multiscale molecular dynamics simulations to obtain insights into pore forming toxin mechanisms.

Pore forming toxins (PFTs) are virulent proteins released by several species, including many strains of bacteria, to attack and kill host cells. In this article, we focus on the utility of molecular dynamics (MD) simulations and the molecular insights gleaned from these techniques on the pore forming pathways of PFTs. In addition to all-atom simulations which are widely used, coarse-grained MARTINI models and structure-based models have also been used to study PFTs. Here, the emphasis is on methods and techniques involved while setting up, monitoring, and evaluating properties from MD simulations of PFTs in a membrane environment. We draw from several case studies to illustrate how MD simulations have provided molecular insights into protein-protein and protein-lipid interactions, lipid dynamics, conformational transitions and structures of both the oligomeric intermediates and assembled pore structures.

Predicting interfacial hot-spot residues that stabilize protein-protein interfaces in oligomeric membrane-toxin pores through hydrogen bonds and salt bridges.

Pore forming toxins (PFTs) are proteins which form unregulated oligomeric pores on target plasma membranes to cause ion leakage and cell death and represent the largest class of bacterial virulence factors. With increasing antibiotic-resistant bacterial strains, alternate therapies are being developed to target toxin pore formation rather than the bacteria themselves. One strategy is to undermine the stability of these multimeric pores, whose subunits are held together by complex amino acid interaction networks, by identifying key residues in these networks which could be plausible drug or mutagenesis targets. However, this requires a quantitative assessment of per residue contributions towards pore stability, which cannot be reliably gleaned from static crystal/cryo-EM pore structures. In this study, we overcome this limitation by developing a computational screening algorithm that employs fully atomistic molecular dynamics simulations coupled with energy-based screening that can predict 'hot-spot' residues which engage in persistent and stabilizing hydrogen bonds or salt bridges across protein-protein interfaces. Application of this algorithm to prototypical α-PFT (cytolysin A) and ß-PFT (α-hemolysin) membrane-inserted pores yielded a small predicted set of highly interacting residues, blocking of which could destabilize pore complexes. Previous mutagenesis studies validate some of our in silico predictions. The algorithm could be applied to all pores with known structures to generate a database of destabilizing mutations, which could then serve as a basis for experimental validation and rational structure-based inhibitor design.Communicated by Ramaswamy H. Sarma.