Multiscale simulation unravel the kinetic mechanisms of inflammasome assembly

Protein assembly and crowding simulations

Proteins encounter frequent molecular interactions in biological environments. Computer simulations have become an increasingly important tool in providing mechanistic insights into how such interactions in vivo relate to their biological function. The review here focuses on simulations describing protein assembly and molecular crowding effects as two important aspects that are distinguished mainly by how specific and long-lived protein contacts are. On the topic of crowding, recent simulations have provided new insights into how crowding affects protein folding and stability, modulates enzyme activity, and affects diffusive properties. Recent studies of assembly processes focus on assembly pathways, especially for virus capsids, amyloid aggregation pathways, and the role of multivalent interactions leading to phase separation. Also, discussed are technical challenges in achieving increasingly realistic simulations of interactions in cellular environments.

Computational Simulation of Holin S105 in Membrane Bilayer and Its Dimerization Through a Helix-Turn-Helix Motif

During the final step of the bacteriophage infection cycle, the cytoplasmic membrane of host cells is disrupted by small membrane proteins called holins. The function of holins in cell lysis is carried out by forming a highly ordered structure called lethal lesion, in which the accumulation of holins in the cytoplasmic membrane leads to the sudden opening of a hole in the middle of this oligomer. Previous studies showed that dimerization of holins is a necessary step to induce their higher order assembly. However, the molecular mechanism underlying the holin-mediated lesion formation is not well understood. In order to elucidate the functions of holin, we first computationally constructed a structural model for our testing system: the holin S105 from bacteriophage lambda. All atom molecular dynamic simulations were further applied to refine its structure and study its dynamics as well as interaction in lipid bilayer. Additional simulations on association between two holins provide supportive evidence to the argument that the C-terminal region of holin plays a critical role in regulating the dimerization. In detail, we found that the adhesion of specific nonpolar residues in transmembrane domain 3 (TMD3) in a polar environment serves as the driven force of dimerization. Our study therefore brings insights to the design of binding interfaces between holins, which can be potentially used to modulate the dynamics of lesion formation.

Understanding the impacts of cellular environments on ligand binding of membrane receptors by computational simulations

Binding of cell surface receptors with their extracellular ligands initiates various intracellular signaling pathways. However, our understanding of the cellular functions of these receptors is very limited due to the fact that in vivo binding between ligands and receptors has only been successfully measured in a very small number of cases. In living cells, receptors are anchored on surfaces of the plasma membrane, which undergoes thermal undulations. Moreover, it has been observed in various systems that receptors can be organized into oligomers prior to ligand binding. It is not well understood how these cellular factors play roles in regulating the dynamics of ligand-receptor interactions. Here, we tackled these problems by using a coarse-grained kinetic Monte Carlo simulation method. Using this method, we demonstrated that the membrane undulations cause a negative effect on ligand-receptor interactions. We further found that the preassembly of membrane receptors on the cell surface can not only accelerate the kinetics of ligand binding but also reduce the noises during the process. In general, our study highlights the importance of membrane environments in regulating the function of membrane receptors in cells. The simulation method can be potentially applied to specific receptor systems involved in cell signaling.

Understand the Functions of Scaffold Proteins in Cell Signaling by a Mesoscopic Simulation Method

Scaffold proteins are central players in regulating the spatial-temporal organization of many important signaling pathways in cells. They offer physical platforms to downstream signaling proteins so that their transient interactions in a crowded and heterogeneous environment of cytosol can be greatly facilitated. However, most scaffold proteins tend to simultaneously bind more than one signaling molecule, which leads to the spatial assembly of multimeric protein complexes. The kinetics of these protein oligomerizations are difficult to quantify by traditional experimental approaches. To understand the functions of scaffold proteins in cell signaling, we developed a, to our knowledge, new hybrid simulation algorithm in which both spatial organization and binding kinetics of proteins were implemented. We applied this new technique to a simple network system that contains three molecules. One molecule in the network is a scaffold protein, whereas the other two are its binding targets in the downstream signaling pathway. Each of the three molecules in the system contains two binding motifs that can interact with each other and are connected by a flexible linker. By applying the new simulation method to the model, we show that the scaffold proteins will promote not only thermodynamics but also kinetics of cell signaling given the premise that the interaction between the two signaling molecules is transient. Moreover, by changing the flexibility of the linker between two binding motifs, our results suggest that the conformational fluctuations in a scaffold protein play a positive role in recruiting downstream signaling molecules. In summary, this study showcases the capability of computational simulation in understanding the general principles of scaffold protein functions.