The Assembly Switch Mechanism of FtsZ Filament Revealed by All-Atom Molecular Dynamics Simulations and Coarse-Grained Models.

Bacterial cytoskeletal protein FtsZ binds and hydrolyzes GTP, and assembles into dynamic filaments that are essential for cell division. Here, we used a multi-scale computational strategy that combined all-atom molecular dynamics (MD) simulations and coarse-grained models to reveal the conformational dynamics of assembled FtsZ. We found that the top end of a filament is highly dynamic and can undergo T-to-R transitions in both GTP- and GDP-bound states. We observed several subcategories of nucleation related dimer species, which leading to a feasible nucleation pathway. In addition, we observed that FtsZ filament exhibits noticeable amounts of twisting, indicating a substantial helicity of the FtsZ filament. These results agree with the previously models and experimental data. Anisotropy network model (ANM) analysis revealed a polymerization enhanced assembly cooperativity, and indicated that the cooperative motions in FtsZ are encoded in the structure. Taken together, our study provides a molecular-level understanding of the diversity of the structural states of FtsZ and the relationships among polymerization, hydrolysis, and cooperative assembly, which should shed new light on the molecular basis of FtsZ's cooperativity.

Molecular basis for ligand activation of the human KCNQ2 channel.

The voltage-gated potassium channel KCNQ2 is responsible for M-current in neurons and is an important drug target to treat epilepsy, pain and several other diseases related to neuronal hyper-excitability. A list of synthetic compounds have been developed to directly activate KCNQ2, yet our knowledge of their activation mechanism is limited, due to lack of high-resolution structures. Here, we report cryo-electron microscopy (cryo-EM) structures of the human KCNQ2 determined in apo state and in complex with two activators, ztz240 or retigabine, which activate KCNQ2 through different mechanisms. The activator-bound structures, along with electrophysiology analysis, reveal that ztz240 binds at the voltage-sensing domain and directly stabilizes it at the activated state, whereas retigabine binds at the pore domain and activates the channel by an allosteric modulation. By accurately defining ligand-binding sites, these KCNQ2 structures not only reveal different ligand recognition and activation mechanisms, but also provide a structural basis for drug optimization and design.

An efficient perturbation method to predict the functionally key sites of glutamine binding protein.

Glutamine-Binding Protein (GlnBP) of Escherichia coli, an important member of the periplasmic binding protein family, is responsible for the first step in the active transport of glutamine across the cytoplasmic membrane. In this work, the functionally key regulation sites of GlnBP were identified by utilizing a perturbation method proposed by our group, in which the residues whose perturbations markedly change the binding free energy between GlnBP and glutamine are considered to be functionally key residues. The results show that besides the substrate binding sites, some other residues distant from the binding pocket, including the ones in the hinge regions between the two domains, the front- and back- door channels and the exposed region, are important for the function of glutamine binding and transport. The predicted results are well consistent with the theoretical and experimental data, which indicates that our method is an effective approach to identify the key residues important for both ligand binding and long-range allosteric signal transmission. This work can provide some insights into the function performance of GlnBP and the physical mechanism of its allosteric regulation.

Approach to the unfolding and folding dynamics of add A-riboswitch upon adenine dissociation using a coarse-grained elastic network model.

Riboswitches are noncoding mRNA segments that can regulate the gene expression via altering their structures in response to specific metabolite binding. We proposed a coarse-grained Gaussian network model (GNM) to examine the unfolding and folding dynamics of adenosine deaminase (add) A-riboswitch upon the adenine dissociation, in which the RNA is modeled by a nucleotide chain with interaction networks formed by connecting adjoining atomic contacts. It was shown that the adenine binding is critical to the folding of the add A-riboswitch while the removal of the ligand can result in drastic increase of the thermodynamic fluctuations especially in the junction regions between helix domains. Under the assumption that the native contacts with the highest thermodynamic fluctuations break first, the iterative GNM simulations showed that the unfolding process of the adenine-free add A-riboswitch starts with the denature of the terminal helix stem, followed by the loops and junctions involving ligand binding pocket, and then the central helix domains. Despite the simplified coarse-grained modeling, the unfolding dynamics and pathways are shown in close agreement with the results from atomic-level MD simulations and the NMR and single-molecule force spectroscopy experiments. Overall, the study demonstrates a new avenue to investigate the binding and folding dynamics of add A-riboswitch molecule which can be readily extended for other RNA molecules.