The Solvent-Exposed C-Terminus of the Cytolysin A Pore-Forming Toxin Directs Pore Formation and Channel Function in Membranes

Polymerized porin as a novel delivery platform for coronavirus vaccine

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), seriously threatens human life and health. The correct folding and polymerization of the receptor-binding domain (RBD) protein of coronavirus in Escherichia coli may reduce the cost of SARS-CoV-2 vaccines. In this study, we constructed this nanopore by using the principle of ClyA porin polymerization triggered by the cell membrane. We used surfactants to "pick" the ClyA-RBD nanopore from the bacterial outer membrane. More importantly, the polymerized RBD displayed on the ClyA-RBD polymerized porin (RBD-PP) already displays some correct spatial conformational epitopes that can induce neutralizing antibodies. The nanostructures of RBD-PP can target lymph nodes and promote antigen uptake and processing by dendritic cells, thereby effectively eliciting the production of anti-SARS-CoV-2 neutralizing antibodies, systemic cellular immune responses, and memory T cells. We applied this PP-based vaccine platform to fabricate an RBD-based subunit vaccine against SARS-CoV-2, which will provide a foundation for the development of inexpensive coronavirus vaccines. The development of a novel vaccine delivery system is an important part of innovative drug research. This novel PP-based vaccine platform is likely to have additional applications, including other viral vaccines, bacterial vaccines, tumor vaccines, drug delivery, and disease diagnosis.

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.

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.

Mechanistic Insights into Pore Formation by an α-Pore Forming Toxin: Protein and Lipid Bilayer Interactions of Cytolysin A

Pore forming toxins (PFTs) are the largest class of bacterial toxins playing a central role in bacterial pathogenesis. They are proteins specifically designed to form nanochannels in the membranes of target cells, ultimately resulting in cell death and establishing infection. PFTs are broadly classified as α- and β-PFTs, depending on secondary structures that form the transmembrane channel. A unique feature about this class of proteins is the drastic conformational changes and complex oligomerization pathways that occur upon exposure to the plasma membrane. A molecular understanding of pore formation has implications in designing novel intervention strategies to combat rising antimicrobial resistance, targeted-cancer therapy, as well as designing nanopores for specialized technologies. Central to unraveling the pore formation pathway is the availability of high resolution crystal structures. In this regard, β-toxins are better understood, when compared with α-toxins whose pore forming mechanisms are complicated by an incomplete knowledge of the driving forces for amphiphatic membrane-inserted helices to organize into functional pores. With the publication of the first crystal structure for an α-toxin, cytolysin A (ClyA), in 2009 we embarked on an extensive multiscale study to unravel its pore forming mechanism. This Account represents the collective mechanistic knowledge gained in our laboratories using a variety of experimental and theoretical techniques which include large scale molecular dynamics (MD) simulations, kinetic modeling studies, single-molecule fluorescence imaging, and super-resolution spectroscopy. We reported MD simulations of the ClyA protomer, oligomeric intermediates, and full pore complex in a lipid bilayer and mapped the conformational transitions that accompany membrane binding. Using single-molecule fluorescence imaging, the conformational transition was experimentally verified by analysis of various diffusion states of membrane bound ClyA. Importantly, we have uncovered a hitherto unknown putative cholesterol binding motif in the membrane-inserted helix of ClyA. Distinct binding pockets for cholesterol formed by adjacent membrane-inserted helices are revealed in MD simulations. Cholesterol appears to play a dual role by stabilizing both the membrane-inserted protomer as well as oligomeric intermediates. Molecular dynamics simulations and kinetic modeling studies suggest that the membrane-inserted arcs oligomerize reversibly to form the predominant transmembrane oligomeric intermediates during pore formation. We posit that this mechanistic understanding of the complex action of α-PFTs has implications in unraveling pore assembly across the wider family of bacterial toxins. With emerging antimicrobial resistance, alternate therapies may rely on disrupting pore functionality or oligomerization of these pathogenic determinants utilized by bacteria, and our study includes assessing the potential for dendrimers as pore blockers.

Assessing the extent of the structural and dynamic modulation of membrane lipids due to pore forming toxins: insights from molecular dynamics simulations

Infections caused by many virulent bacterial strains are triggered by the release of pore forming toxins (PFTs), which form oligomeric transmembrane pore complexes on the target plasma membrane. The spatial extent of the perturbation to the surrounding lipids during pore formation is relatively unexplored. Using all-atom molecular dynamics simulations, we investigate the changes in the structure and dynamics of lipids in a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer in the presence of contrasting PFTs. Cytolysin A (ClyA), an α toxin with its inserted wedge shaped bundle of inserted α helices, induces significant asymmetry across the membrane leaflets in comparison with α hemolysin (AHL), a β toxin. Despite the differences in hydrophobic mismatch and uniquely different topologies of the two oligomers, perturbations to lipid order as reflected in the tilt angle and order parameters and membrane thinning are short ranged, lying within ∼2.5 nm from the periphery of either pore complex, and commensurate with distances typically associated with van der Waals forces. In contrast, the spatial extent of perturbations to the lipid dynamics extends outward to at least 4 nm for both proteins, and the continuous survival probabilities reveal the presence of a tightly bound shell of lipids in this region. Displacement probability distributions show long tails and the distinctly non-Gaussian features reflect the induced dynamic heterogeneity. A detailed profiling of the protein-lipid contacts with tyrosine, tryptophan, lysine and arginine residues shows increased non-polar contacts in the cytoplasmic leaflet for both PFTs, with a higher number of atomic contacts in the case of AHL in the extracellular leaflet due to the mushroom-like topology of the pore complex. The short ranged nature of the perturbations observed in this simple one component membrane suggests inherent plasticity of membrane lipids enabling the recovery of the structure and membrane fluidity even in the presence of these large oligomeric transmembrane protein assemblies. This observation has implications in membrane repair processes such as budding or vesicle fusion events used to mitigate PFT virulence, where the underlying lipid dynamics and fluidity in the vicinity of the pore complex are expected to play an important role.

Correlated protein conformational states and membrane dynamics during attack by pore-forming toxins

Pore-forming toxins (PFTs) are a class of proteins implicated in a wide range of virulent bacterial infections and diseases. These toxins bind to target membranes and subsequently oligomerize to form functional pores that eventually lead to cell lysis. While the protein undergoes large conformational changes on the bilayer, the connection between intermediate oligomeric states and lipid reorganization during pore formation is largely unexplored. Cholesterol-dependent cytolysins (CDCs) are a subclass of PFTs widely implicated in food poisoning and other related infections. Using a prototypical CDC, listeriolysin O (LLO), we provide a microscopic connection between pore formation, lipid dynamics, and leakage kinetics by using a combination of Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) measurements on single giant unilamellar vesicles (GUVs). Upon exposure to LLO, two distinct populations of GUVs with widely different leakage kinetics emerge. We attribute these differences to the existence of oligomeric intermediates, sampling various membrane-bound conformational states of the protein, and their intimate coupling to lipid rearrangement and dynamics. Molecular dynamics simulations capture the influence of various membrane-bound conformational states on the lipid and cholesterol dynamics, providing molecular interpretations to the FRET and FCS experiments. Our study establishes a microscopic connection between membrane binding and conformational changes and their influence on lipid reorganization during PFT-mediated cell lysis. Additionally, our study provides insights into membrane-mediated protein interactions widely implicated in cell signaling, fusion, folding, and other biomolecular processes.

Cholesterol promotes Cytolysin A activity by stabilizing the intermediates during pore formation

Pore-forming toxins (PFTs) form nanoscale pores across target membranes causing cell death. Cytolysin A (ClyA) from Escherichia coli is a prototypical α-helical toxin that contributes to cytolytic phenotype of several pathogenic strains. It is produced as a monomer and, upon membrane exposure, undergoes conformational changes and finally oligomerizes to form a dodecameric pore, thereby causing ion imbalance and finally cell death. However, our current understanding of this assembly process is limited to studies in detergents, which do not capture the physicochemical properties of biological membranes. Here, using single-molecule imaging and molecular dynamics simulations, we study the ClyA assembly pathway on phospholipid bilayers. We report that cholesterol stimulates pore formation, not by enhancing initial ClyA binding to the membrane but by selectively stabilizing a protomer-like conformation. This was mediated by specific interactions by cholesterol-interacting residues in the N-terminal helix. Additionally, cholesterol stabilized the oligomeric structure using bridging interactions in the protomer-protomer interfaces, thereby resulting in enhanced ClyA oligomerization. This dual stabilization of distinct intermediates by cholesterol suggests a possible molecular mechanism by which ClyA achieves selective membrane rupture of eukaryotic cell membranes. Topological similarity to eukaryotic membrane proteins suggests evolution of a bacterial α-toxin to adopt eukaryotic motifs for its activation. Broad mechanistic correspondence between pore-forming toxins hints at a wider prevalence of similar protein membrane insertion mechanisms.

Nanoscale Heterogeneities Drive Enhanced Binding and Anomalous Diffusion of Nanoparticles in Model Biomembranes

Interaction of functional nanoparticles with cells and model biomembranes has been widely studied to evaluate the effectiveness of the particles as potential drug delivery vehicles and bioimaging labels as well as in understanding nanoparticle cytotoxicity effects. Charged nanoparticles, in particular, with tunable surface charge have been found to be effective in targeting cellular membranes as well as the subcellular matrix. However, a microscopic understanding of the underlying physical principles that govern nanoparticle binding, uptake, or diffusion on cells is lacking. Here, we report the first experimental studies of nanoparticle diffusion on model biomembranes and correlate this to the existence of nanoscale dynamics and structural heterogeneities using super-resolution stimulated emission depletion (STED) microscopy. Using confocal and STED microscopy coupled with fluorescence correlation spectroscopy (FCS), we provide novel insight on why these nanoparticles show enhanced binding on two-component lipid bilayers as compared to single-component membranes and how binding and diffusion is correlated to subdiffraction nanoscale dynamics and structure. The enhanced binding is also dictated, in part, by the presence of structural and dynamic heterogeneity, as revealed by STED-FCS studies, which could potentially be used to understand enhanced nanoparticle binding in raft-like domains in cell membranes. In addition, we also observe a clear correlation between the enhanced nanoparticle diffusion on membranes and the extent of membrane penetration by the nanoparticles. Our results not only have a significant impact on our understanding of nanoparticle binding and uptake as well as diffusion in cell and biomembranes, but have very strong implications for uptake mechanisms and diffusion of other biomolecules, like proteins on cell membranes and their connections to functional membrane nanoscale platform.

Pathways for creation and annihilation of nanoscale biomembrane domains reveal alpha and beta-toxin nanopore formation processes

Raft-like functional domains with putative sizes of 20–200 nm and which are evolving dynamically are believed to be the most crucial regions in cellular membranes which determine cell signaling and various functions of cells.

Assessing the Structure and Stability of Transmembrane Oligomeric Intermediates of an α-Helical Toxin

Protein membrane interactions play an important role in our understanding of diverse phenomena ranging from membrane-assisted protein aggregation to oligomerization and folding. Pore-forming toxins (PFTs) are the primary vehicle for infection by several strains of bacteria. These proteins which are expressed in a water-soluble form (monomers) bind to the target membrane and conformationally transform (protomers) and self-assemble to form a multimer transmembrane pore complex through a process of oligomerization. On the basis of the structure of the transmembrane domains, PFTs are broadly classified into β or α toxins. In contrast to β-PFTs, the paucity of available crystal structures coupled with the amphipathic nature of the transmembrane domains has hindered our understanding of α-PFT pore formation. In this article, we use molecular dynamics (MD) simulations to examine the process of pore formation of the bacterial α-PFT, cytolysin A from Escherichia coli (ClyA) in lipid bilayer membranes. Using atomistic MD simulations ranging from 50 to 500 ns, we show that transmembrane oligomeric intermediates or "arcs" form stable proteolipidic complexes consisting of protein arcs with toroidal lipids lining the free edges. By creating initial conditions where the lipids are contained within the arcs, we study the dynamics of spontaneous lipid evacuation and toroidal edge formation. This process occurs on the time scale of tens of nanoseconds, suggesting that once protomers oligomerize, transmembrane arcs are rapidly stabilized to form functional water channels capable of leakage. Using umbrella sampling with a coarse-grained molecular model, we obtain the free energy of insertion of a single protomer into the membrane. A single inserted protomer has a stabilization free energy of -52.9 ± 1.2 kJ/mol and forms a stable transmembrane water channel capable of leakage. Our simulations reveal that arcs are stable and viable intermediates that can occur during the pore-formation pathway for ClyA.