PubMed Health⌕ Search

Biomedical subjects

Mitsunori Ikeguchi

Publications and source records attributed to Mitsunori Ikeguchi.

14 recordsLinked to original sources
PubMed2006

Structural changes in the cytoplasmic domain of phospholamban by phosphorylation at Ser16: a molecular dynamics study.

Phospholamban is a 52-residue integral membrane protein that regulates the activity of the sarcoplasmic reticulum calcium pump in cardiac muscle. Its inhibitory action is relieved when phospholamban is phosphorylated at Ser16 by cAMP-dependent protein kinase. To computationally explore all possible conformations of the phosphorylated form, and thereby to understand the structural effects of phosphorylation, replica-exchange molecular dynamics (REMD) was applied to the cytoplasmic domain that includes Ser16. The simulations showed that (i) without phosphorylation, the region from Lys3 to Ser16 takes all alpha-helical conformations; (ii) when phosphorylated, the alpha-helix is partially unwound in the C-terminal part (from Ser10 to Ala15) resulting in less extended conformations; (iii) the phosphate at Ser16 forms salt bridges with Arg9, Arg13, and/or Arg14; and (iv) the salt bridges with Arg13 and Arg14 distort the alpha-helix and induce unwinding of the C-terminal part. We then applied conventional all-atom molecular dynamics simulations to the full-length phospholamban in the phospholipid bilayer. The results were consistent with those obtained with REMD simulations, suggesting that the transmembrane part of phospholamban and the lipid bilayer itself have only minor effects on the conformational changes in the cytoplasmic domain. The distortions caused by the salt bridges involving the phosphate at Ser16 readily explain the relief of the inhibitory effect of phospholamban by phosphorylation, as they will substantially reduce the population of all helical conformations, which are presumably required for the binding to the calcium pump. This will also be the mechanism for releasing the phosphorylated phospholamban from kinase.

Calcium-Binding Proteins↗
PubMed2006

A novel magnetic resonance-based method to measure gene expression in living cells.

In unicellular and multicellular eukaryotes, elaborate gene regulatory mechanisms facilitate a broad range of biological processes from cell division to morphological differentiation. In order to fully understand the gene regulatory networks involved in these biological processes, the spatial and temporal patterns of expression of many thousands of genes will need to be determined in real time in living organisms. Currently available techniques are not sufficient to achieve this goal; however, novel methods based on magnetic resonance (MR) imaging may be particularly useful for sensitive detection of gene expression in opaque tissues. This report describes a novel reporter gene system that monitors gene expression dynamically and quantitatively, in yeast cells, by measuring the accumulation of inorganic polyphosphate (polyP) using MR spectroscopy (MRS) or MR spectroscopic imaging (MRI). Because this system is completely non-invasive and does not require exogenous substrates, it is a powerful tool for studying gene expression in multicellular organisms, as well.

Blotting, Northern↗
PubMed2005

Comparative simulations of aquaporin family: AQP1, AQPZ, AQP0 and GlpF.

Molecular dynamics simulations were performed for four members of the aquaporin family (AQP1, AQPZ, AQP0, and GlpF) in the explicit membrane environment. The single-channel water permeability, pf, was evaluated to be GlpF approximately AQPZ > AQP1 >> AQP0, while their relative pore sizes were GlpF >> AQP1 > AQPZ >> AQP0. This relation between pf and pore size indicates that water permeability was determined not only by the channel radius, but also another competing factor. Analysis of water dynamics revealed that this factor was the single-file nature of water transport.

Amino Acid Sequence↗
PubMed2005

Atomically detailed description of the unfolding of alpha-lactalbumin by the combined use of experiments and simulations.

The recombinant form of goat alpha-lactalbumin has a significantly faster unfolding rate compared to the authentic form, although the two molecules differ only in an extra methionine at the N terminus of the recombinant. The mechanism of the destabilization caused by this residue was investigated through the combined use of kinetic experiments and molecular dynamics simulations. Unfolding simulations for the authentic and recombinant forms at 398 K (ten trajectories of 5 ns for each form, 100 ns total) precisely reproduced the experimentally observed differences in unfolding behavior. In addition, experiments reproduced the destabilization of a mutant protein, T38A, faithfully as predicted by the simulations. This bidirectional verification between experiments and simulations enabled the atomically detailed description of the role of the extra methionine residue in the unfolding process.

Amino Acid Substitution↗
PubMed2005

Protonation of the acidic residues in the transmembrane cation-binding sites of the ca(2+) pump.

The ionization states of the acidic residues around the Ca2+-binding sites of sarcoplasmic reticulum Ca2+ ATPase are studied by continuum electrostatic calculations and all-atom molecular dynamics simulations with explicit solvent and phospholipids. The two methods consistently indicate that Glu58 and Glu908 are protonated at neutral pH. The Ca2+ coordination and the H-bonds formed by the protonation of Glu58 and Glu908 are stable in an MD simulation, whereas the H-bonds are disrupted and the Ca2+ coordination geometry is severely altered in another simulation treating these residues unprotonated. The results clearly indicate that the H-bonds formed by protonation of Glu58 and Glu908 provide extra stability for the Ca2+-binding sites of Ca2+ ATPase.

Binding Sites↗
PubMed2005

Protein structural change upon ligand binding: linear response theory.

A simple formula based on linear response theory is proposed to explain and predict the structural change of proteins upon ligand binding. By regarding ligand binding as an external perturbation, the structural change as a response is described by atomic fluctuations in the ligand-free form and the protein-ligand interactions. The results for three protein systems of various sizes are consistent with the observations in the crystal structures, confirming the validity of the linear relationship between the equilibrium fluctuations and the structural change upon ligand binding.

Acetyl Coenzyme A↗
PubMed2004

Phylogeny of protein-folding trajectories reveals a unique pathway to native structure.

To scrutinize how a protein folds at atomic resolution, we performed 200 molecular dynamics simulations (each of 50 ns) of the miniprotein Trp-cage on the computational grid. Within the trajectories, 58 folding and 31 unfolding events were identified and subjected to extensive comparison and classification. Based on an analogy with biological sequences, the folding and unfolding trajectories (arrays of sequential snapshots of structures) were aligned by dynamic programming allowing gaps. A phylogenetic tree derived from the alignments revealed four distinct groups of the trajectories, characterized by the Trp side-chain motions and the main-chain motions. It was found that only one group attained the native structure and that the other three led to pseudonative structures having the correct main-chain trace but different nonnative Trp side-chain rotamers, indicating that those four folded structures were each attained through a unique folding pathway.

Computational Biology↗
PubMed2004

Partial rigid-body dynamics in NPT, NPAT and NPgammaT ensembles for proteins and membranes.

A partial rigid-body method of molecular dynamics simulations for proteins and membranes is presented. In this method, the symplectic integrator for rigid bodies is combined with the equations of motion for the NPT ensemble. The standard NPT ensemble is extended to the membrane-specific ensembles, the NPAT (constant normal pressure and lateral surface area of membranes and constant temperature) and NPgammaT (constant normal pressure and lateral surface tension of membranes and constant temperature) ensembles. By more than 30-ns simulations of aqueous proteins and hydrated lipid bilayers, the results of the partial rigid-body method demonstrated excellent conservation of total energy and consistent behavior with the traditional constraint method in terms of structural distribution and fluctuation of proteins and lipids. The efficient implementation of the partial rigid-body method in parallel computation is presented, which is shown to work well in large-scale molecular dynamics simulations.

Algorithms↗
PubMed2004

The crystal structure of the tryptophan synthase beta subunit from the hyperthermophile Pyrococcus furiosus. Investigation of stabilization factors.

The structure of the tryptophan synthase beta2 subunit (Pfbeta2) from the hyperthermophile, Pyrococcus furiosus, was determined by X-ray crystallographic analysis at 2.2 A resolution, and its stability was examined by DSC. This is the first report of the X-ray structure of the tryptophan synthase beta2 subunit alone, although the structure of the tryptophan synthase alpha2beta2 complex from Salmonella typhimurium has already been reported. The structure of Pfbeta2 was essentially similar to that of the beta2 subunit (Stbeta2) in the alpha2beta2 complex from S. typhimurium. The sequence alignment with secondary structures of Pfbeta and Stbeta in monomeric form showed that six residues in the N-terminal region and three residues in the C-terminal region were deleted in Pfbeta, and one residue at Pro366 of Stbeta and at Ile63 of Pfbeta was inserted. The denaturation temperature of Pfbeta2 was higher by 35 degrees C than the reported values from mesophiles at approximately pH 8. On the basis of structural information on both proteins, the analyses of the contributions of each stabilization factor indicate that: (a) the higher stability of Pfbeta2 is not caused by either a hydrophobic interaction or an increase in ion pairs; (b) the number of hydrogen bonds involved in the main chains of Pfbeta is greater by about 10% than that of Stbeta, indicating that the secondary structures of Pfbeta are more stabilized than those of Stbeta and (c) the sequence of Pfbeta seems to be better fitted to an ideally stable structure than that of Stbeta, as assessed from X-ray structure data.

Amino Acid Sequence↗
PubMed2003

Two different Swi5-containing protein complexes are involved in mating-type switching and recombination repair in fission yeast.

Homologous recombination is an important biological process that occurs in all organisms and facilitates genome rearrangements and repair of DNA double-strand breaks. Eukaryotic Rad51 proteins (Rad51sp or Rhp51 in fission yeast) are functional and structural homologs of bacterial RecA protein, an evolutionarily conserved protein that plays a key role in homologous pairing and strand exchange between homologous DNA molecules in vitro. Here we show that the fission yeast swi5+ gene, which was originally identified as a gene required for normal mating-type switching, encodes a protein conserved among eukaryotes and is involved in a previously uncharacterized Rhp51 (Rad51sp)-dependent recombination repair pathway that does not require the Rhp55/57 (Rad55/57sp) function. Protein interactions with both Swi5 and Rhp51 were found to be mediated by a domain common to Swi2 and Sfr1 (Swi five-dependent recombination repair protein 1, a previously uncharacterized protein with sequence similarity to the C-terminal part of Swi2). Genetic epistasis analyses suggest that the Swi5-Sfr1-Rhp51 interactions function specifically in DNA recombination repair, whereas the Swi5-Swi2-Rhp51 interactions may function, together with chromodomain protein Swi6 (HP1 homolog), in mating-type switching.

Amino Acid Sequence↗
PubMed2002

New method for parallel computation of Hessian matrix of conformational energy function in internal coordinates.

A new algorithm for parallel calculation of the second derivatives (Hessian) of the conformational energy function of biomolecules in internal coordinates is proposed. The basic scheme of this algorithm is the division of the entire calculation of the Hessian matrix (called "task") into subtasks and the optimization of the assignment of processors to each subtask by considering both the load balancing and reduction of the communication cost. A genetic algorithm is used for this optimization considering the dependencies between subtasks. We applied this method to a glutaminyl transfer RNA (Gln-tRNA) molecule for which the scalability of our previously developed parallel algorithm was significantly decreased when the large number of processors was used. The speedup for the calculation was 32.6 times with 60 processors, which is considerably better than the speedup for our previously reported parallel algorithm. The elapsed time for the calculation of subtasks, data sending, and data receiving was analyzed, and the effect of the optimization using the genetic algorithm is discussed.

Algorithms↗
PubMed2001

Solution structure of the fibronectin type III domain from Bacillus circulans WL-12 chitinase A1.

Growing evidence suggests that horizontal gene transfer plays an integral role in the evolution of bacterial genomes. One of the debated examples of horizontal gene transfer from animal to prokaryote is the fibronectin type III domain (FnIIID). Certain extracellular proteins of soil bacteria contain an unusual cluster of FnIIIDs, which show sequence similarity to those of animals and are likely to have been acquired horizontally from animals. Here we report the solution structure of the FnIIID of chitinase A1 from Bacillus circulans WL-12. To the best of our knowledge, this is the first tertiary structure to be reported for an FnIIID from a bacterial protein. The structure of the domain shows significant similarity to FnIIIDs from animal proteins. Sequence comparisons with FnIIIDs from other soil bacteria proteins show that the core-forming residues are highly conserved and, thus, are under strong evolutionary pressure. Striking similarities in the tertiary structures of bacterial FnIIIDs and their mammalian counterparts may support the hypothesis that the evolution of the FnIIID in bacterial carbohydrases occurred horizontally. The total lack of surface-exposed aromatic residues also suggests that the role of this FnIIID is different from those of other bacterial beta-sandwich domains, which function as carbohydrate-binding modules.

Amino Acid Sequence↗