Enhancement of protein mechanical stability: Correlated deformations are handcuffed by ligand binding.

As revealed by previous experiments, protein mechanical stability can be effectively regulated by ligand binding with the binding site distant from the force-bearing region. However, the mechanism for such long-range allosteric control of protein mechanics is still largely unknown. In this work, we use protein topology-based elastic network model (ENM) and all-atomic steered molecular dynamics (SMD) simulations to study the impact of ligand binding on protein mechanical stability in two systems, i.e., GB1 and CheY-binding P2-domain of CheA (CBDCheA). Both ENM and SMD results show that the ligand binding has considerable and negligible effects on the mechanical stability of these two proteins, respectively. These results are consistent with the experimental observations. A physical mechanism for the enhancement of protein mechanical stability was then proposed: the correlated deformations of the force-bearing region and the binding site are handcuffed by the binding of ligand. The handcuff effect suppresses the propagation of internal force in the force-bearing region, thus improving the resistance to the loading force. Our study indicates that ENM method can effectively identify the structure motifs allosterically related to the deformation in the force bearing region, as well as the force propagation pathway within the structure of the studied proteins. Hence, it should be helpful to understand the molecular origin of the different mechanical properties in response to ligand binding for GB1 and CBDCheA.

Energy transport pathway in proteins: Insights from non-equilibrium molecular dynamics with elastic network model.

Intra-molecular energy transport between distant functional sites plays important roles in allosterically regulating the biochemical activity of proteins. How to identify the specific intra-molecular signaling pathway from protein tertiary structure remains a challenging problem. In the present work, a non-equilibrium dynamics method based on the elastic network model (ENM) was proposed to simulate the energy propagation process and identify the specific signaling pathways within proteins. In this method, a given residue was perturbed and the propagation of energy was simulated by non-equilibrium dynamics in the normal modes space of ENM. After that, the simulation results were transformed from the normal modes space to the Cartesian coordinate space to identify the intra-protein energy transduction pathways. The proposed method was applied to myosin and the third PDZ domain (PDZ3) of PSD-95 as case studies. For myosin, two signaling pathways were identified, which mediate the energy transductions form the nucleotide binding site to the 50 kDa cleft and the converter subdomain, respectively. For PDZ3, one specific signaling pathway was identified, through which the intra-protein energy was transduced from ligand binding site to the distant opposite side of the protein. It is also found that comparing with the commonly used cross-correlation analysis method, the proposed method can identify the anisotropic energy transduction pathways more effectively.

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.