Sotaro Fuchigami

Slow dynamics in protein fluctuations revealed by time-structure based independent component analysis: The case of domain motions

Protein dynamics on a long time scale was investigated using all-atom molecular dynamics (MD) simulation and time-structure based independent component analysis (tICA). We selected the lysine-, arginine-, ornithine-binding protein (LAO) as a target protein and focused on its domain motions in the open state. A MD simulation of the LAO in explicit water was performed for 600 ns, in which slow and large-amplitude domain motions of the LAO were observed. After extracting domain motions by rigid-body domain analysis, the tICA was applied to the obtained rigid-body trajectory, yielding slow modes of the LAOs domain motions in order of decreasing time scale. The slowest mode detected by the tICA represented not a closure motion described by a largest-amplitude mode determined by the principal component analysis but a twist motion with a time scale of tens of nanoseconds. The slow dynamics of the LAO were well described by only the slowest mode and were characterized by transitions between two basins. The results show that tICA is promising for describing and analyzing slow dynamics of proteins.

Classification and Annotation of the Relationship between Protein Structural Change and Ligand Binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank-one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.

Slow dynamics of a protein backbone in molecular dynamics simulation revealed by time-structure based independent component analysis.

We recently proposed the method of time-structure based independent component analysis (tICA) to examine the slow dynamics involved in conformational fluctuations of a protein as estimated by molecular dynamics (MD) simulation [Y. Naritomi and S. Fuchigami, J. Chem. Phys. 134, 065101 (2011)]. Our previous study focused on domain motions of the protein and examined its dynamics by using rigid-body domain analysis and tICA. However, the protein changes its conformation not only through domain motions but also by various types of motions involving its backbone and side chains. Some of these motions might occur on a slow time scale: we hypothesize that if so, we could effectively detect and characterize them using tICA. In the present study, we investigated slow dynamics of the protein backbone using MD simulation and tICA. The selected target protein was lysine-, arginine-, ornithine-binding protein (LAO), which comprises two domains and undergoes large domain motions. MD simulation of LAO in explicit water was performed for 1 μs, and the obtained trajectory of C(α) atoms in the backbone was analyzed by tICA. This analysis successfully provided us with slow modes for LAO that represented either domain motions or local movements of the backbone. Further analysis elucidated the atomic details of the suggested local motions and confirmed that these motions truly occurred on the expected slow time scale.

Gas-Phase Structure of the Histone Multimers Characterized by Ion Mobility Mass Spectrometry and Molecular Dynamics Simulation

The minimum structural unit of chromatin is the nucleosome core particle (NCP), consisting of 146 bp of DNA wrapped around a histone octamer, which itself contains two H2A/H2B dimers and one (H3/H4)2 tetramer. These multimers possess functionally important tail regions that are intrinsically disordered. In order to elucidate the mechanisms behind NCP assembly and disassembly processes, which are highly related to gene expression, structural characterization of the H2A/H2B dimer and (H3/H4)2 tetramer will be of importance. In the present study, human histone multimers with disordered tail regions were characterized by electrospray ionization (ESI) ion mobility-mass spectrometry (IM-MS) and molecular dynamics (MD) simulation. Experimentally obtained arrival times of these histone multimer ions showed rather wide distributions, implying that multiple conformers exist for each histone multimer in the gas phase. To examine their structures, MD simulations of the histone multimers were performed first in solution and then in vacuo at four temperatures, resulting in a variety of histone multimer structures. Theoretical collision cross-section (CCS) values calculated for the simulated structures revealed that structural models with smaller CCS values had more compact tail regions than those with larger CCS values. This implied that variation of the CCS values of the histone multimers were primarily due to the random behaviors of the tail regions in the gas phase. The combination of IM-MS and MD simulation enabled clear and comprehensive characterization of the gas-phase structures of histone multimers containing disordered tails.

Linear response theory in dihedral angle space for protein structural change upon ligand binding

Coupling between proteins motion and ligand binding can be well explained by the linear response theory (Ikeguchi, M.; Ueno, J.; Sato, M.; Kidera, A. Phys Rev Lett 2005, 94, 078102.), in which the structural change is treated as a response to ligand binding. The prediction accuracy of structural change upon ligand binding has been improved by replacing the variables in the linear response theory from Cartesian coordinates to dihedral angles. The dihedral angle theory can more accurately describe the rotational motions of protein domains compared with the Cartesian theory, which tends to shift the coordinate to the tangential direction of the domain rotation. In this study, the ligand‐bound form of Ferric‐binding protein was predicted from its ligand‐free form using the dihedral linear response theory. When the variance–covariance matrix, the key component in the linear response theory, was derived by linear conversion from Cartesian coordinates to dihedral angles, the dihedral linear response theory gave an improvement in the prediction. Therefore, the description of the rotational motion by dihedral angles is crucial for accurate prediction of protein structural change.

Nickel binding to NikA: An additional binding site reconciles spectroscopy, calorimetry and crystallography

Intracellular nickel is required by Escherichia coli as a cofactor for a number of enzymes and is necessary for anaerobic respiration. However, high concentrations of nickel are toxic, so both import and export systems have evolved to control the cellular level of the metal. The nik operon in E. coli encodes a nickel-uptake system that includes the periplasmic nickel-binding protein NikA. The crystal structures of wild-type NikA both bound to nickel and in the apo form have been solved previously. The liganded structure appeared to show an unusual interaction between the nickel and the protein in which no direct bonds are formed. The highly unusual nickel coordination suggested by the crystal structure contrasted strongly with earlier X-ray spectroscopic studies. The known nickel-binding site has been probed by extensive mutagenesis and isothermal titration calorimetry and it has been found that even large numbers of disruptive mutations appear to have little effect on the nickel affinity. The crystal structure of a binding-site mutant with nickel bound has been solved and it is found that nickel is bound to two histidine residues at a position distant from the previously characterized binding site. This novel site immediately resolves the conflict between the crystal structures and other biophysical analyses. The physiological relevance of the two binding sites is discussed.

Latent dynamics of a protein molecule observed in dihedral angle space

Dihedral angles are alternative set of variables to Cartesian coordinates for representing protein dynamics. The two sets of variables exhibit extremely different behavior. Motions in dihedral angle space are characterized by latent dynamics, in which motion induced in each dihedral angle is always compensated for by motions of many other dihedral angles, in order to maintain a rigid globular shape. Using molecular dynamics simulations, we propose a molecular mechanism for the latent dynamics in dihedral angle space. It was found that, due to the unique structure of dihedral principal components originating in the globular shape of the protein, the dihedral principal components with large (small) amplitudes are highly correlated with the eigenvectors of the metric matrix with small (large) eigenvalues. Such an anticorrelation in the eigenmode structures minimizes the mean square displacement of Cartesian coordinates upon rotation of dihedral angles. In contrast, a short peptide, deca-alanine in this study, does not show such behavior of the latent dynamics in the dihedral principal components, but shows similar behaviors to those of the Cartesian principal components, due to the absence of constraints to maintain a rigid globular shape.

Multivariate frequency domain analysis of protein dynamics.

Multivariate frequency domain analysis (MFDA) is proposed to characterize collective vibrational dynamics of protein obtained by a molecular dynamics (MD) simulation. MFDA performs principal component analysis (PCA) for a bandpass filtered multivariate time series using the multitaper method of spectral estimation. By applying MFDA to MD trajectories of bovine pancreatic trypsin inhibitor, we determined the collective vibrational modes in the frequency domain, which were identified by their vibrational frequencies and eigenvectors. At near zero temperature, the vibrational modes determined by MFDA agreed well with those calculated by normal mode analysis. At 300 K, the vibrational modes exhibited characteristic features that were considerably different from the principal modes of the static distribution given by the standard PCA. The influences of aqueous environments were discussed based on two different sets of vibrational modes, one derived from a MD simulation in water and the other from a simulation in vacuum. Using the varimax rotation, an algorithm of the multivariate statistical analysis, the representative orthogonal set of eigenmodes was determined at each vibrational frequency.

CLASP2 Has Two Distinct TOG Domains That Contribute Differently to Microtubule Dynamics.

CLIP-associated proteins CLASPs are mammalian microtubule (MT) plus-end tracking proteins (+TIPs) that promote MT rescue in vivo. Their plus-end localization is dependent on other +TIPs, EB1 and CLIP-170, but in the leading edge of the cell, CLASPs display lattice-binding activity. MT association of CLASPs is suggested to be regulated by multiple TOG (tumor overexpressed gene) domains and by the serine-arginine (SR)-rich region, which contains binding sites for EB1. Here, we report the crystal structures of the two TOG domains of CLASP2. Both domains consist of six HEAT repeats, which are similar to the canonical paddle-like tubulin-binding TOG domains, but have arched conformations. The degrees and directions of curvature are different between the two TOG domains, implying that they have distinct roles in MT binding. Using biochemical, molecular modeling and cell biological analyses, we have investigated the interactions between the TOG domains and αβ-tubulin and found that each domain associates differently with αβ-tubulin. Our findings suggest that, by varying the degrees of domain curvature, the TOG domains may distinguish the structural conformation of the tubulin dimer, discriminate between different states of MT dynamic instability and thereby function differentially as stabilizers of MTs.

Molecular modeling of the HAMP domain of sensory rhodopsin II transducer from Natronomonas pharaonis

The halobacterial transducer of sensory rhodopsin II (HtrII) is a photosignal transducer associated with phototaxis in extreme halophiles. The HAMP domain, a linker domain in HtrII, is considered to play an important role in transferring the signal from the membrane to the cytoplasmic region, although its structure in the complex remains undetermined. To establish the structural basis for understanding the mechanism of signal transduction, we present an atomic model of the structure of the N-terminal HAMP domain from Natronomonas pharaonis (HtrII: 84–136), based on molecular dynamics (MD) simulations. The model was built by homology modeling using the NMR structure of Af1503 from Archaeoglobus fulgidus as a template. The HAMP domains of Af1503 and HtrII were stable during MD simulations over 100 ns. Quantitative analyses of inter-helical packing indicated that the Af1503 HAMP domain stably maintained unusual knobs-to-knobs packing, as observed in the NMR structure, while the bulky side-chains of HtrII shifted the packing state to canonical knobs-into-holes. The role of the connector loop in maintaining structural stability was also discussed using MD simulations of loop deletion mutants.