Computational modeling reveals a hydrophobic force-induced pore-forming mechanism for cholesterol-dependent cytolysins

During the pore-forming process, cholesterol-dependent cytolysins (CDCs) bind to cholesterol-rich membranes and subsequently undergo a series of conformational changes, predominantly involving in the collapse of the protein and the transformation from helices to {beta}-hairpins to form a large hydrophobic pore. In the current study, we reconstructed a structural model for both the prepore and pore-forming complexes of PFO based on the cryo-EM data of pneumolysin and performed molecular dynamics simulations and free energy calculations to study the conformational changes in the PFO prepore-to-pore conversion. Our simulations indicate that D2 cannot collapse spontaneously due to the hydrogen bonding and pi-pi interactions between domains D2 and D3, which are partially weakened by binding to cell membranes and oligomerization. The free energy landscape for the prepore-to-pore conversion reveals that an additional force is required for the collapse of D2 to overcome an energy barrier of [~] 24 kcal/mol. Based on these computational results, we proposed a hydrophobic force-induced pore-forming mechanism to explain the pore-forming process of CDCs. In this mechanism, the hydrophobic interactions between the TMHs and membranes are essential for the prepore-to-pore conversion. The hydrophobic force generated by the TMHs-membrane interactions drives the conformational changes in domains D2 and D3. These findings well explain how the conformational changes within two distant domains synergistically occur, and fits well for the previous biophysical and biochemical data.