Search

Identification of the internal dynamics of multidomain proteins is crucial for clarifying the mechanism of their functions. The neutron spin echo (NSE) technique is well suited for studying internal dynamics. However, the requirement for relatively high protein concentrations and the lack of appropriate analytical methods have impeded the investigation of the internal dynamics with NSE. To overcome these difficulties, we employed a unique approach to study the internal dynamics of a multidomain protein, EcoO109I, whose dynamics was anticipated to be pertinent to DNA degradation. We anticipated a synergetic effect between the NSE measurement at interference-free protein concentration and all-atom molecular dynamics simulation. Through this approach, the internal dynamics of EcoO109I was successfully observed within temporal and spatial scales. Additionally, principal component analysis (PCA) was applied to the internal dynamics trajectory to identify the dominant motion of the internal dynamics. The first PCA mode, which was the most cooperative among all PCA modes, mainly explained the internal dynamics. This dominant mode of EcoO109I exhibited the motion which facilitated both the access of DNA to the recognition site and the cleavage of DNA. Therefore, our approach can identify the functionally relevant internal dynamics of multidomain proteins.

Aggregates cause a fatal problem in the structural analysis of a biomacromolecule in solution using small-angle X-ray or neutron scattering (SAS): they deteriorate the scattering profile of the target molecule and lead to an incorrect structure. Recently, an integrated method of analytical ultracentrifugation (AUC) and SAS, abbreviated AUC–SAS, was developed as a new approach to overcome this problem. However, the original version of AUC–SAS does not offer a correct scattering profile of the target molecule when the weight fraction of aggregates is higher than ca 10%. In this study, the obstacle point in the original AUC–SAS approach is identified. The improved AUC–SAS method is then applicable to a solution with a relatively larger weight fraction of aggregates (≤20%).

Three domain fragments of a multi-domain protein, ER-60, were ligated in two short linker regions using asparaginyl endopeptidase not involving denaturation. To identify appropriate ligation sites, by selecting several potential ligation sites with fewer mutations around two short linker regions, their ligation efficiencies and the functions of the ligated ER-60s were examined experimentally. To evaluate the dependence of ligation efficiencies on the ligation sites computationally, steric hinderances around the sites for the ligation were calculated through molecular dynamics simulations. Utilizing the steric hindrance, a site-dependent ligation potential index was introduced as reproducing the experimental ligation efficiency. Referring to this index, the reconstruction of ER-60 was succeeded by the ligation of the three domains for the first time. In addition, the new ligation potential index well-worked for application to other domain ligations. Therefore, the index may serve as a more time-effective tool for multi-site ligations.

Solving structural ensembles of flexible biomolecules is a challenging research area. Here, we propose a method to obtain possible structural ensembles of a biomolecule based on small-angle X-ray scattering (SAXS) and molecular dynamics simulations. Our idea is to clip a time series that matches a SAXS profile from a simulation trajectory. To examine its practicability, we applied our idea to a multi-domain protein ER-60 and successfully extracted time series longer than 1 micro second from trajectories of coarse-grained molecular dynamics simulations. In the extracted time series, the domain conformation was distributed continuously and smoothly in a conformational space. Preferred domain conformations were also observed. Diversity among scattering curves calculated from each ER-60 structure was interpreted to reflect an open-close motion of the protein. Although our approach did not provide a unique solution for the structural ensemble of the biomolecule, each extracted time series can be an element of the real behavior of ER-60. Considering its low computational cost, our approach will play a key role to identify biomolecular dynamics by integrating SAXS, simulations, and other experiments.

In the cyanobacterial circadian clock system, KaiA, KaiB and KaiC periodically assemble into a large complex. Here we determined the overall structure of their fully assembled complex by integrating experimental and computational approaches. Small-angle X-ray and inverse contrast matching small-angle neutron scatterings coupled with size-exclusion chromatography provided constraints to highlight the spatial arrangements of the N-terminal domains of KaiA, which were not resolved in the previous structural analyses. Computationally built 20 million structural models of the complex were screened out utilizing the constrains and then subjected to molecular dynamics simulations to examine their stabilities. The final model suggests that, despite large fluctuation of the KaiA N-terminal domains, their preferential positionings mask the hydrophobic surface of the KaiA C-terminal domains, hindering additional KaiA-KaiC interactions. Thus, our integrative approach provides a useful tool to resolve large complex structures harboring dynamically fluctuating domains.