Impacts of the charged residues mutation S48E/N62H on the thermostability and unfolding behavior of cold shock protein: insights from molecular dynamics simulation with Go model.

The cold shock protein from the hyperthermophile Thermotoga maritima (Tm-Csp) exhibits significantly higher thermostability than its homologue from the thermophile Bacillus caldolyticus (Bc-Csp). Experimental studies have shown that the electrostatic interactions unique to Tm-Csp are responsible for improving its thermostability. In the present work, the favorable charged residues in Tm-Csp were grafted into Bc-Csp by a double point mutation of S48E/N62H, and the impacts of the mutation on the thermostability and unfolding/folding behavior of Bc-Csp were then investigated by using a modified Go model, in which the electrostatic interactions between charged residues were considered in the model. Our simulation results show that this Tm-Csp-like charged residue mutation can effectively improve the thermostability of Bc-Csp without changing its two-state folding mechanism. Besides that, we also studied the unfolding kinetics and unfolding/folding pathway of the wild-type Bc-Csp and its mutant. It is found that this charged residue mutation obviously enhanced the stability of the C-terminal region of Bc-Csp, which decreases the unfolding rate and changes the unfolding/folding pathway of the protein. Our studies indicate that the thermostability, unfolding kinetics and unfolding/folding pathway of Bc-Csp can be artificially changed by introducing Tm-Csp-like favorable electrostatic interactions into Bc-Csp.

Analysis of conformational motions and related key residue interactions responsible for a specific function of proteins with elastic network model.

Protein collective motions play a critical role in many biochemical processes. How to predict the functional motions and the related key residue interactions in proteins is important for our understanding in the mechanism of the biochemical processes. Normal mode analysis (NMA) of the elastic network model (ENM) is one of the effective approaches to investigate the structure-encoded motions in proteins. However, the motion modes revealed by the conventional NMA approach do not necessarily correspond to a specific function of protein. In the present work, a new analysis method was proposed to identify the motion modes responsible for a specific function of proteins and then predict the key residue interactions involved in the functional motions by using a perturbation approach. In our method, an internal coordinate that accounts for the specific function was introduced, and the Cartesian coordinate space was transformed into the internal/Cartesian space by using linear approximation, where the introduced internal coordinate serves as one of the axes of the coordinate space. NMA of ENM in this internal/Cartesian space was performed and the function-relevant motion modes were identified according to their contributions to the specific function of proteins. Then the key residue interactions important for the functional motions of the protein were predicted as the interactions whose perturbation largely influences the fluctuation along the internal coordinate. Using our proposed methods, the maltose transporter (MalFGK2) from E. Coli was studied. The functional motions and the key residue interactions that are related to the channel-gating function of this protein were successfully identified.

Conformational Motions and Functionally Key Residues for Vitamin B12 Transporter BtuCD-BtuF Revealed by Elastic Network Model with a Function-Related Internal Coordinate.

BtuCD-BtuF from Escherichia coli is a binding protein-dependent adenosine triphosphate (ATP)-binding cassette (ABC) transporter system that uses the energy of ATP hydrolysis to transmit vitamin B12 across cellular membranes. Experimental studies have showed that during the transport cycle, the transporter undergoes conformational transitions between the "inward-facing" and "outward-facing" states, which results in the open-closed motions of the cytoplasmic gate of the transport channel. The opening-closing of the channel gate play critical roles for the function of the transporter, which enables the substrate vitamin B12 to be translocated into the cell. In the present work, the extent of opening of the cytoplasmic gate was chosen as a function-related internal coordinate. Then the mean-square fluctuation of the internal coordinate, as well as the cross-correlation between the displacement of the internal coordinate and the movement of each residue in the protein, were calculated based on the normal mode analysis of the elastic network model to analyze the function-related motions encoded in the structure of the system. In addition, the key residues important for the functional motions of the transporter were predicted by using a perturbation method. In order to facilitate the calculations, the internal coordinate was introduced as one of the axes of the coordinate space and the conventional Cartesian coordinate space was transformed into the internal/Cartesian space with linear approximation. All the calculations were carried out in this internal/Cartesian space. Our method can successfully identify the functional motions and key residues for the transporter BtuCD-BtuF, which are well consistent with the experimental observations.

Prediction of allosteric sites on protein surfaces with an elastic-network-model-based thermodynamic method.

Allostery is a rapid and efficient way in many biological processes to regulate protein functions, where binding of an effector at the allosteric site alters the activity and function at a distant active site. Allosteric regulation of protein biological functions provides a promising strategy for novel drug design. However, how to effectively identify the allosteric sites remains one of the major challenges for allosteric drug design. In the present work, a thermodynamic method based on the elastic network model was proposed to predict the allosteric sites on the protein surface. In our method, the thermodynamic coupling between the allosteric and active sites was considered, and then the allosteric sites were identified as those where the binding of an effector molecule induces a large change in the binding free energy of the protein with its ligand. Using the proposed method, two proteins, i.e., the 70 kD heat shock protein (Hsp70) and GluA2 alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) receptor, were studied and the allosteric sites on the protein surface were successfully identified. The predicted results are consistent with the available experimental data, which indicates that our method is a simple yet effective approach for the identification of allosteric sites on proteins.