Experimental and theoretical studies have showed that the native-state topology conceals a wealth of information about protein folding/unfolding. In this study, a method based on the Gaussian network model (GNM) is developed to study some properties of protein unfolding and explore the role of topology in protein unfolding process. The GNM has been successful in predicting atomic fluctuations around an energy minimum. However, in the GNM, the normal mode description is linear and cannot be accurate in studying protein folding/unfolding, which has many local minima in the energy landscape. To describe the nonlinearity of the conformational changes during protein unfolding, a method based on the iterative use of normal mode calculation is proposed. The protein unfolding process is mimicked through breaking the native contacts between the residues one by one according to the fluctuations of the distance between them. With this approach, the unfolding processes of two proteins, CI2 and barnase, are simulated. It is found that the sequence of protein unfolding events revealed by this method is consistent with that obtained from thermal unfolding by molecular dynamics and Monte Carlo simulations. The results indicate that this method is effective in studying protein unfolding. In this method, only the native contacts are considered, which implies that the native topology may play an important role in the protein unfolding process. The simulation results also show that the unfolding pathway is robust against the introduction of some noise, or stochastic characters. Furthermore, several conformations selected from the unfolding process are studied to show that the denatured state does not behave as a random coil, but seems to have highly cooperative motions, which may help and promote the polypeptide chain to fold into the native state correctly and speedily.

In this work, the mechanism of domain movements of glutamine-binding protein (GlnBP), especially the influence of the ligand on GlnBP dynamic behavior is investigated with the aid of a Gaussian network model (GNM) and an anisotropy elastic network model. The results show that the "open-closed" transition mainly appears as the large movement of the small domain, especially the top region including two alpha-helices and two beta-strands. The slowest mode of each three forms of GlnBP--ligand-free open, ligand-bound closed, and ligand-free closed GlnBP--shows that the open-closed motion of the two domains has a common hinge axis centered on Lys-87 and Gln-183. Accompanying the conformational transition, the residues within both large and small domains move in a highly coupled way. The peaks of the fast modes correspond to residues that were thought, in the GNM, to be important for the stability of the protein, and these residues may be involved in the interactions with the membrane-bound components. With the contacts between the large domain and the small domain increasing, the ability of the "open-closed" motion is decreased. All the results agree well with those of molecular dynamics simulations, and it is thought that the open-closed conformation transition is the nature of the topology structure of GlnBP. Also, the influence of the ligand on GlnBP is studied with a modified GNM method. The results obtained show that the ligand does not influence the closed-to-open transition tendency.

The interaction of HIV-1 trans-activator protein Tat with its cognate trans-activation response element (TAR) RNA is critical for viral transcription and replication. Therefore, it has long been considered as an attractive target for the development of antiviral compounds. Recently, the conformationally constrained cyclic peptide mimetics of Tat have been tested to be a promising family of lead peptides. Here, we focused on two representative cyclic peptides termed as L-22 and KP-Z-41, both of which exhibit excellent inhibitory potency against Tat and TAR interaction. By means of molecular dynamics simulations, we obtained a detailed picture of the interactions between them and HIV-1 TAR RNA. In results, it is found that the binding modes of the two cyclic peptides to TAR RNA are almost identical at or near the bulge regions, whereas the binding interfaces at the apical loop exhibit large conformational heterogeneity. In addition, it is revealed that electrostatic interaction energy contributes much more to KP-Z-41 complex formation than to L-22 complex, which is the main source of energy that results in a higher binding affinity of KP-Z-41 over-22 for TAR RNA. Furthermore, we identified a conserved motif RRK (Arg-Arg-Lys) that is shown to be essential for specific binding of this class of cyclic peptides to TAR RNA. This work can provide a useful insight into the design and modification of cyclic peptide inhibitors targeting the association of HIV-1 Tat and TAR RNA.

No abstract

AMPA receptor mediates the fast excitatory synaptic transmission in the central nervous system, and it is activated by the binding of glutamate that results in the opening of the transmembrane ion channel. In the present work, the thermodynamic method developed by our group was improved and then applied to identify the functionally key residues that regulate the glutamate-binding affinity of AMPA receptor. In our method, the key residues are identified as those whose perturbation largely changes the ligand binding free energy of the protein. It is found that besides the ligand binding sites, other residues distant from the binding cleft can also influence the glutamate binding affinity through a long-range allosteric regulation. These allosteric sites include the hinge region of the ligand binding cleft, the dimer interface of the ligand binding domain, the linkers between the ligand binding domain and the transmembrane domain, and the interface between the N-terminal domain and the ligand binding domain. Our calculation results are consistent with the available experimental data. The results are helpful for our understanding of the mechanism of long-range allosteric communication in the AMPA receptor and the mechanism of channel opening triggered by glutamate binding.