Mechanochemical Model of the Power Stroke of the Single-Headed Motor Protein KIF1A

Binding patterns of inhibitors to different pockets of kinesin Eg5

The kinesin-5 family member, Eg5, plays very important role in the mitosis. As a mitotic protein, Eg5 is the target of various mitotic inhibitors. There are two targeting pockets in the motor domain of Eg5, which locates in the α2/L5/α3 region and the α4/α6 region respectively. We investigated the interactions between the different inhibitors and the two binding pockets of Eg5 by using all-atom molecular dynamics method. Combined the conformational analysis with the free-energy calculation, the binding patterns of inhibitors to the two binding pockets are shown. The α2/L5/α3 pocket can be divided into 4 regions. The structures and binding conformations of inhibitors in region 1 and 2 are highly conserved. The shape of α4/α6 pocket is alterable. The space of this pocket in ADP-binding state of Eg5 is larger than that in ADP·Pi-binding state due to the limitation of a hydrogen bond formed in the ADP·Pi-binding state. The results of this investigation provide the structural basis of the inhibitor-Eg5 interaction and offer a reference for the Eg5-targeted drug design.

Energy landscape quantifications of histone H3.3 recognition by chaperone DAXX reveal an uncoupled binding specificity and affinity

Histone variant H3.3 differs from the canonical histone H3.1 by only five amino acids, yet its chaperone death domain-associated protein (DAXX) can specifically recognize H3.3 over H3.1, despite having a large DAXX-interacting surface on the H3.3-H4 heterodimer common to that on the H3.1-H4 complex. This observation gives rise to the question of, from the binding energy point view, how high binding specificity may be achieved with small differences of the overall binding energy for protein-protein interactions in general. Here we investigate the mechanism of coupling of binding specificity and affinity in protein-protein interactions using the DAXX-H3.3-H4 complex as a model. Using a multi-scale method, we found that the hydrophobic interactions between DAXX and the H3.3-specific region contributed to their initial binding process. And the structural flexibility of the interacting partners contributed to the binding affinity after their encounter. By quantifying the free energy landscape, we revealed that the interaction between the specific residues of H3.3 and DAXX decreased the encounter barrier height while the folding of H3.3-H4 and DAXX increased the depth of the free energy basin of the final binding state. The encounter barrier height, which is not coupled to the thermodynamic stability of the final binding state, had a marked effect on the initial binding rate of flexible histones and chaperones. Based on the energy landscape theory, we found that the intrinsic binding energy funnel of this uncoupled recognition process was affected by the structural flexibility and the flexibility modulated the degree of coupling between binding specificity and affinity. Our work offers a biophysical explanation of the specific recognition between the histones and their chaperones, and also extends the use of energy landscape theory for understanding molecular recognitions in general.

ATP Acts as a Hydrotrope to Regulate the Phase Separation of NBDY Clusters

Nonannotated P-body dissociating polypeptide (NBDY) is a recently discovered human microprotein that has been found to be a novel component of the mRNA decapping complex. Previous studies have shown that the phosphorylation of NBDY promotes the liquid phase of the NBDY remixing in vitro. Typically, during the process of phosphorylation, a phosphate group is added to the protein through adenosine triphosphate (ATP) hydrolysis. It has been shown that ATP acts as a biological hydrotrope, affecting the phase separation of proteins in solution. In this study, we utilized simulation methods to investigate the dynamic properties of the NBDY clusters at various ATP concentrations. Our findings demonstrate that ATP can regulate the phase separation of NBDY clusters. Specifically, we identified a critical point in the concentration ratio between ATP and NBDY that exhibits a dual effect on the phase separation of NBDY. We observed that the nonsaturated ATP concentration can facilitate the formation of phase separation, while oversaturated ATP concentration promotes the diffusion of NBDY, and the oversaturated ATP-NBDY interaction impedes the phase separation of NBDY. Additionally, we found that ATPs can bind to the protein surface by aggregating into ATP clusters, which further hinders the diffusion of NBDY clusters. Our work provides general insight into the role of ATP in the phase separation of protein condensates.

Perspectives on the landscape and flux theory for describing emergent behaviors of the biological systems

We give a review on the landscape theory of the equilibrium biological systems and landscape-flux theory of the nonequilibrium biological systems as the global driving force. The emergences of the behaviors, the associated thermodynamics in terms of the entropy and free energy and dynamics in terms of the rate and paths have been quantitatively demonstrated. The hierarchical organization structures have been discussed. The biological applications ranging from protein folding, biomolecular recognition, specificity, biomolecular evolution and design for equilibrium systems as well as cell cycle, differentiation and development, cancer, neural networks and brain function, and evolution for nonequilibrium systems, cross-scale studies of genome structural dynamics and experimental quantifications/verifications of the landscape and flux are illustrated. Together, this gives an overall global physical and quantitative picture in terms of the landscape and flux for the behaviors, dynamics and functions of biological systems.

Effect of Binding-Affinity and ATPase Activity on the Velocities of Kinesins Using Ratchet Models

We use two-state ratchet models containing single and coupled Brownian motors to understand the role of motor-microtubule binding, ATPase reaction rate and dimerisation on the translational velocities of Kinesin motors. We use model parameters derived from the experimental measurements on KIF1A, KIF13A, KIF13B, and KIF16B motors to compute velocities in μm/s. We observe that both the models show the same trend in velocities (KIF1A > KIF13A > KIF13B > KIF16B) as the experimental results. However, the models significantly underpredict the velocities when compared with the experiments. The predictions of the coupled-motor model are closer to the experiments than those of the single-motor model. Our results indicate that the variation of ATPase reaction rate governs the trend in velocities for the above four motors. The variation of motor-microtubule binding affinity and the coupling strength between the motor domains may only have a secondary effect. More rigorous models that incorporate the power-stroke mechanism are necessary for better quantitative compliance with the experiments.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

A kinetic dissection of the fast and superprocessive kinesin-3 KIF1A reveals a predominant one-head-bound state during its chemomechanical cycle

The kinesin-3 family contains the fastest and most processive motors of the three neuronal transport kinesin families, yet the sequence of states and rates of kinetic transitions that comprise the chemomechanical cycle and give rise to their unique properties are poorly understood. We used stopped-flow fluorescence spectroscopy and single-molecule motility assays to delineate the chemomechanical cycle of the kinesin-3, KIF1A. Our bacterially expressed KIF1A construct, dimerized via a kinesin-1 coiled-coil, exhibits fast velocity and superprocessivity behavior similar to WT KIF1A. We established that the KIF1A forward step is triggered by hydrolysis of ATP and not by ATP binding, meaning that KIF1A follows the same chemomechanical cycle as established for kinesin-1 and -2. The ATP-triggered half-site release rate of KIF1A was similar to the stepping rate, indicating that during stepping, rear-head detachment is an order of magnitude faster than in kinesin-1 and kinesin-2. Thus, KIF1A spends the majority of its hydrolysis cycle in a one-head-bound state. Both the ADP off-rate and the ATP on-rate at physiological ATP concentration were fast, eliminating these steps as possible rate-limiting transitions. Based on the measured run length and the relatively slow off-rate in ADP, we conclude that attachment of the tethered head is the rate-limiting transition in the KIF1A stepping cycle. Thus, KIF1A's activity can be explained by a fast rear-head detachment rate, a rate-limiting step of tethered-head attachment that follows ATP hydrolysis, and a relatively strong electrostatic interaction with the microtubule in the weakly bound post-hydrolysis state.

How Kinesin-1 Utilize the Energy of Nucleotide: The Conformational Changes and Mechanochemical Coupling in the Unidirectional Motion of Kinesin-1

Kinesin-1 is a typical motile molecular motor and the founding member of the kinesin family. The most significant feature in the unidirectional motion of kinesin-1 is its processivity. To realize the fast and processive movement on the microtubule lattice, kinesin-1 efficiently transforms the chemical energy of nucleotide binding and hydrolysis to the energy of mechanical movement. The chemical and mechanical cycle of kinesin-1 are coupled to avoid futile nucleotide hydrolysis. In this paper, the research on the mechanical pathway of energy transition and the regulating mechanism of the mechanochemical cycle of kinesin-1 is reviewed.

Multiscale landscape of molecular mechanism of SIRT1 activation by STACs

Sirtuins are a family of highly conserved NAD-dependent deacetylase that are involved in multiple biological processes in both prokaryotes and eukaryotes. Many sirtuin-activating compounds (STACs) have been reported for SIRT1, which is the best-characterized sirtuin. However, the molecular mechanism of SIRT1 activation by STACs remains controversial. Here, we developed a multiscale simulation model to explore this mechanism. By quantifying the free energy landscape for the closed conformation of a SIRT1-FdL peptide-resveratrol complex, we found a positive correlation between the barrier height of the active free energy basin and the experimentally determined fluctuations in the rate of SIRT1 deacetylation by resveratrol. In addition, by monitoring dynamics, we found that the open conformation of a SIRT1-p53-STAC-1 complex had a faster rate of conformational change than the closed structure. We also determined the structural properties of each thermodynamic or dynamic state and found that two potential activating factors, the stability of FdL peptide (the p53 peptide substrate including an AMC fluorophore group) binding and the stability of the SIRT1 conformation, were weakly correlated under certain conditions. These results address the controversial question of whether the AMC fluorophore group and native hydrophobic residues have similar roles in the SIRT1 activation process. Finally, we captured the global landscape of the transition, including less stable and more stable states, and proposed a global physical landscape for the mechanism of SIRT1 activation by STACs.