PubMed Health⌕ Search

Biomedical subjects

Akio Kitao

Publications and source records attributed to Akio Kitao.

10 recordsLinked to original sources

Gap compression/extension mechanism of bacterial flagellar hook as the molecular universal joint.

Bacterial flagellar hook acts as a molecular universal joint, transmitting torque produced by the flagellar basal body, a rotary motor, to the flagellar filament. The hook forms polymorphic supercoil structures and can be considered as an assembly of 11 circularly arranged protofilaments. We investigated the molecular mechanism of the universal joint function of the hook by a approximately two-million-atom molecular dynamics simulation. On the inner side of the supercoil, protein subunits are highly packed along the protofilament and no gaps remain for further compression, whereas subunits are slightly separated and are hydrogen bonded through one layer of water molecules on the outer side. As for the intersubunit interactions between protofilaments, subunits are packed along the 6-start helix in a left-handed supercoil whereas they are highly packed along the 5-start helix in a right-handed supercoil. We conclude that the supercoiled structures of the hook in the left- and right-handed forms make maximal use of the gaps between subunits, which we call "gap compression/extension mechanism". Mutual sliding of subunits at the subunit interface accompanying rearrangements of intersubunit hydrogen bonds is interpreted as a mechanism to allow continuous structural change of the hook during flagellar rotation at low energy cost.

Amino Acid Sequence↗

Molecular dynamics simulations of NAD+-induced domain closure in horse liver alcohol dehydrogenase.

Horse liver alcohol dehydrogenase is a homodimer, the protomer having a coenzyme-binding domain and a catalytic domain. Using all available x-ray structures and 50 ns of molecular dynamics simulations, we investigated the mechanism of NAD+-induced domain closure. When the well-known loop at the domain interface was modeled to its conformation in the closed structure, the NAD+-induced domain closure from the open structure could be simulated with remarkable accuracy. Native interactions in the closed structure between Arg369, Arg47, His51, Ala317, Phe319, and NAD+ were seen to form at different stages during domain closure. Removal of the Arg369 side-chain charge resulted in the loss of the tendency to close, verifying that specific interactions do help drive the domains closed. Further simulations and a careful analysis of x-ray structures suggest that the loop prevents domain closure in the absence of NAD+, and a cooperative mechanism operates between the subunits for domain closure. This cooperative mechanism explains the role of the loop as a block to closure because in the absence of NAD+ it would prevent the occurrence of an unliganded closed subunit when the other subunit closes on NAD+. Simulations that started with one subunit open and one closed supported this.

Alcohol Dehydrogenase↗

Switch interactions control energy frustration and multiple flagellar filament structures.

Bacterial flagellar filament is a macromolecular assembly consisting of a single protein, flagellin. Bacterial swimming is controlled by the conformational transitions of this filament between left- and right-handed supercoils induced by the flagellar motor torque. We present a massive molecular dynamics simulation that was successful in constructing the atomic-level supercoil structures consistent with various experimental data and further in elucidating the detailed underlying molecular mechanisms of the polymorphic supercoiling. We have found that the following three types of interactions are keys to understanding the supercoiling mechanism. "Permanent" interactions are always maintained between subunits in the various supercoil structures. "Sliding" interactions are formed between variable hydrophilic or hydrophobic residue pairs, allowing intersubunit shear without large change in energy. The formation and breakage of "switch" interactions stabilize inter- and intrasubunit interactions, respectively. We conclude that polymorphic supercoiling is due to the energy frustration between them. The transition between supercoils is achieved by a "transform and relax" mechanism: the filament structure is geometrically transformed rapidly and then slowly relaxes to energetically metastable states by rearranging interactions.

Computer Simulation↗

Amplitudes and directions of internal protein motions from a JAM analysis of 15N relaxation data.

A method has been developed for characterizing dynamic structures of proteins in solution by using nuclear magnetic resonance (NMR) restraints and 15N relaxation data. This method is based on the concept of the jumping-among-minima (JAM) model. In this model we assume that protein dynamics can be described on the basis of conformational substates, and involves intra- and inter-substate motion. A set of substates is created by picking energy-minimized conformations from the conformational space consistent with the geometric NMR restraints. Intra-substate motions, which occur on the timescale of approximately 10 ps, are simulated with molecular dynamics (MD) calculations with force-field energy terms. Statistical weights of the conformational substates are determined to reproduce the NMR relaxation parameters. The refinement procedure consists of four stages: (i) determination of the ensemble of structures that satisfy NMR restraints, (ii) determination of intra-substate fluctuation, (iii) determination of statistical weights of conformational substates to reproduce model-free relaxation parameters, and (iv) analysis of the resulting dynamic structure to determine amplitudes and directions of internal protein motions. This method was employed to investigate structure and dynamics of the adhesion domain of human CD2 (hCD2) in solution. Two major collective modes, whose contributions to atomic mean-square fluctuations are 77.1% in total, are identified by the refinement. The first mode is interpreted as a rigid-body motion of a protein segment consisting of a part of the B--C loop, a part of the F strand, and the F--G loop. Another type of smaller-amplitude mode is indicated for the C'--C'' loop. The motions affect primarily the curvature of the slightly concave counterreceptor-binding site and represent transitions between a concave (closed) and flat (open) binding face. By comparing the ensemble of structures in solution to the complex structure with counterreceptor CD58, we found that these two types of motions resemble the change upon counterreceptor binding.

Models, Molecular↗

Protein boson peak originated from hydration-related multiple minima energy landscape.

The boson peak is a broad peak found in the low-frequency region of inelastic neutron and Raman scattering spectra in many glassy materials, including biopolymers below approximately 200 K. Here, we give a novel insight into the origins of the protein boson peak, which may also be valid for materials other than proteins. Molecular simulation reveals that the structured water molecules around a protein molecule increase the number of local minima in the protein energy landscape, which plays a key role in the origin of the boson peak. The peak appears when the protein dynamics are trapped within a local energy minimum at cryogenic temperatures. This trapping causes very low frequency collective motions to shift to higher frequencies. We demonstrate that the characteristic frequency of such systems shifts higher as the temperature decreases also in model one-dimensional energy surfaces with multiple minima.

Animals↗

Structure of the bacterial flagellar hook and implication for the molecular universal joint mechanism.

The bacterial flagellum is a motile organelle, and the flagellar hook is a short, highly curved tubular structure that connects the flagellar motor to the long filament acting as a helical propeller. The hook is made of about 120 copies of a single protein, FlgE, and its function as a nano-sized universal joint is essential for dynamic and efficient bacterial motility and taxis. It transmits the motor torque to the helical propeller over a wide range of its orientation for swimming and tumbling. Here we report a partial atomic model of the hook obtained by X-ray crystallography of FlgE31, a major proteolytic fragment of FlgE lacking unfolded terminal regions, and by electron cryomicroscopy and three-dimensional helical image reconstruction of the hook. The model reveals the intricate molecular interactions and a plausible switching mechanism for the hook to be flexible in bending but rigid against twisting for its universal joint function.

Bacterial Proteins↗

How can free energy component analysis explain the difference in protein stability caused by amino acid substitutions? Effect of three hydrophobic mutations at the 56th residue on the stability of human lysozyme.

To elucidate the molecular mechanism of thermal stability, it is essential to determine what are the major free energy components that contribute significantly to the total free energy difference caused by amino acid mutations. In this work, we carried out free energy calculations based on all-atom molecular dynamics simulations to investigate the effect of three hydrophobic mutations at the same position, I56A, I56V and I56F of human lysozyme. The calculated free energy differences are in good agreement with the experimental values in all cases. From free energy component analysis, we found that small changes in stability in the I56A and I56V mutants originate from the short-range Lennard-Jones interactions, whereas the I56F mutant is largely destabilized owing to the changes in the long-range electrostatic interactions. The calculated results are also compared with the free energy components determined by an empirical relationship based on the native-state structure and thermodynamic data. Although this relationship has been shown to be very successful in reproducing the stability changes caused by various amino acid substitutions in several proteins, the changes of stability in I56V and I56F mutants are not reproduced very well. By comparing the free energy components calculated by these two approaches, we showed that the effect of the long-range interaction on the stability changes may be underestimated in the empirical relationships when the structural change caused by mutation is relatively small, as in I56F. It is also suggested that estimation of the change in accessible surface area, deltadeltaASA, may be overestimated if the structure around the mutation site in the denatured state is native-like, which would cause overestimation of the free energy change as in the case of I56V. Our results clearly show that the combined approach of the free energy calculation based on the all-atom molecular dynamics simulation and the empirical relationships is very useful for understanding the detailed mechanism of protein stability.

Amino Acid Substitution↗