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.

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.

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.