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 wildtype 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 Nterminal 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.

The transition of an α-helix to a β-sheet in proteins is among the most complex conformational changes seen in bio-molecular systems. Currently, it is difficult to study such protein conformational changes in a direct molecular dynamics simulation. 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 [ J. Chem. Phys. 126, 054103 (2007)] and combine it with the finite temperature string (FTS) method [J. Phys. Chem. B 109, 6688-6693 (2005)] 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.

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 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 E. coli is less understood. In a recent free energy computation, we have shown 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 have carried 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 stabilizing the unfolded intermediate. In the absence of cholesterol the membrane was found to undergo large curvature deformations in both leaflets of the membrane accompanied by bilayer thinning. Our study with the α-toxin, ClyA illustrates that cholesterol is critical to catalyzing and stabilizing the unfolded state of the β-tongue in the membrane, opening up fresh insights into cholesterol assisted unfolding of membrane proteins.

Bacterial pore-forming toxins (PFTs) bind and oligomerize on mammalian cell membranes forming nanopores, that cause cell lysis to promote a wide range of bacterial infections. Cytolysin A (ClyA), an alpha(α)-PFT, is known to undergo one of the largest conformational changes during its transition from a water soluble monomeric form to the membrane embedded dodecameric nanopore assembly. Despite extensive work on the structure and assembly of ClyA, a complete molecular picture of the interplay between the protein segments and membrane lipids driving this transformation remains elusive. In this study, we combine experiments and all-atom molecular dynamics (MD) simulations of ClyA and its mutants to unravel the role of the two key membrane interacting motifs, namely, the β-tongue and N-terminus helix, in facilitating this critical transition. Erythrocyte turbidity and vesicle leakage assays reveal a loss of activity for β-tongue mutant (Y178F), and delayed kinetics for the N-terminus mutants (Y27A and Y27F). All atom, thermal unfolding molecular dynamics simulations of the monomer carried out at 310, 350 and 400 K reveal a distinct reduction in the flexibility in both the β-tongue and N-terminal regions of the mutants when compared with the wild type. This decreased loss of conformational flexbility correlates positively with the reduced lytic and leakage activity observed in experiments, indicating that the tendency to lose secondary structure in the β-tongue region is an important step in the conformational transition bistability of the ClyA protein. Simulations with the membrane inserted oligomeric arcs representing the pore state reveal a greater destabilization tendency among the β-tongue mutant as inferred from secondary structure and N-terminal positioning. Our combined experimental and simulation study, reveals that conformational flexibility is indispensable for the outward movement of the β-tongue and the tendency to induce disorder in the β-tongue is an important step in the transition to the membrane mediated helix-turn-helix motif integral to ClyA pore formation. This observed loss of secondary structure is akin to the structural transitions observed in intrinsically disordered proteins (IDPs) to support protein function. Our finding suggest that inherent flexibility in the protein could play a wider and hitherto unrecognized role in the membrane mediated conformational transitions of PFTs in general.

No abstract