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

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.

Molecular dynamics simulation studies on binding of activator and inhibitor to Munc13-1 C1 in the presence of membrane

Munc13-1 is a presynaptic active zone protein that plays a critical role in priming the synaptic vesicle and releasing neurotransmitters in the brain. Munc13-1 acts as a scaffold and is activated when diacylglycerol (DAG)/phorbol ester binds to its C1 domain in the plasma membrane. Our previous studies showed that bryostatin 1 activated the Munc13-1, but resveratrol inhibited the phorbol ester-induced Munc13-1 activity. To gain structural insights into the binding of the ligand into Munc13-1 C1 in the membrane, we conducted 1.0 μs molecular dynamics (MD) simulation on Munc13-1 C1-ligand-lipid ternary system using phorbol 13-acetate, bryostatin 1 and resveratrol as ligands. Munc13-1 C1 shows higher conformational stability and less mobility along membrane with phorbol 13-acetate and bryostatin 1 than with resveratrol. Bryostatin 1 and phorbol ester remained in the protein active site, but resveratrol moved out of Munc13-1 C1 during the MD simulation. While bryostatin 1-bound Munc13-1 C1 showed two different positioning in the membrane, phorbol 13-acetate and resveratrol-bound Munc13-1 C1 only showed one positioning. Phorbol 13-acetate formed hydrogen bond with Ala-574 and Gly-589. Bryostatin 1 had more hydrogen bonds with Trp-588 and Arg-592 than with other residues. Resveratrol formed hydrogen bond with Ile-590. This study suggests that different ligands control Munc13-1 C1's mobility and positioning in the membrane differently. Ligand also has a critical role in the interaction between Munc13-1 C1 and lipid membrane. Our results provide structural basis of the pharmacological activity of the ligands and highlight the importance of membrane in Munc13-1 activity.Communicated by Ramaswamy H. Sarma.

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.

Electric Field-Mediated Fibronectin-Hydroxyapatite Interaction: A Molecular Insight

In experimental research-driven biomaterials science, the influence of different material properties (elastic stiffness, surface energy, etc.) and, to a relatively lesser extent, biophysical stimulation (electric/magnetic) on cell-material interactions has been extensively investigated. Despite the central importance of protein adsorption on cell-material interactions, the quantitative analysis to probe into the role of physicochemical factors in protein adsorption remains largely unexplored in biomaterials science. In recent studies, the critical role of electric field stimulation toward the modulation of cell functionality in implantable biomaterials has been experimentally demonstrated. Given this background, we investigated the influence of external electric field stimulation (upto 1.00 V/nm) on fibronectin (FN) adsorption on a hydroxyapatite (HA) (001) surface at 300 K using the all-atom molecular dynamics (MD) simulation method. FN adsorption was found to be governed by attractive electrostatic interactions, which changed with the electric field strength. Nonmonotonous changes in the structural integrity of FN were recorded with the change in the field strength and direction. This can be attributed to the spatial rearrangement of the positions of local charges and the global structural changes of proteins. The dipole moment vectors of FN, water, and HA quantitatively exhibited a similar pattern of orienting themselves parallel to the field direction, with field strength-dependent increase in their magnitudes. No significant change has been recorded for the radial distribution function of water surrounding FN. Field-dependent variation in the salt bridge nets and the number of hydrogen bonds between FN and HA were also examined. One of the important results in the context of cell-material interaction is that the RGD (Arg-Gly-Asp) sequence of FN was exposed to the solvent side when the field was applied along an outward direction perpendicular to the HA (001) surface. In summary, the present study provides molecular insights into the influence of electric field stimulation on phenomenological interactions involved in FN adsorption on the HA surface.