Search PubMed⌕ Search

Biomedical subjects

Shoji Takada

Publications and source records attributed to Shoji Takada.

16 recordsLinked to original sources
PubMed2006

Secondary structure provides a template for the folding of nearby polypeptides.

Although protein structures are primarily encoded by their sequences, they are also critically dependent on environmental factors such as solvents and interactions with other molecules. Here we investigate how the folding-energy landscape of a short peptide is altered by interactions with another peptide, by performing atomistic replica-exchange molecular dynamics simulations of polyalanines in various environments. We analyzed the free-energy landscapes of Ala7 and Ala8 in isolation, near an alpha-helix template, and near a beta-strand template. The isolated Ala7 and Ala8 at 270 K were mainly in polyproline II helix conformations and in equilibrium between the alpha-helix and polyproline II helix, respectively, in harmony with the experiment. Interestingly, we found remarkably strong secondary-structure "templating"; namely, the alpha-helix template enhanced alpha-helix conformation and the beta-strand template induced beta-strand conformation in the simulated Ala8. The alpha-helix template lowered the nearby dielectric constant, which strengthened hydrogen bonds in the simulated Ala8, leading to alpha-helix stabilization. The beta-strand template provided hydrogen bond positions to the simulated Ala8, sharply inducing beta-strand structure. With or without templates, the energy landscape of Ala8 is always funnel-like and centered at the alpha-helix conformation, whereas entropic contribution disfavors the alpha-helix, leading to subtle competition. Secondary-structure templating may play a critical role in protein conformation dynamics in the cellular environment.

Computer Simulation↗
PubMed2006

Multiple-basin energy landscapes for large-amplitude conformational motions of proteins: Structure-based molecular dynamics simulations.

Biomolecules often undergo large-amplitude motions when they bind or release other molecules. Unlike macroscopic machines, these biomolecular machines can partially disassemble (unfold) and then reassemble (fold) during such transitions. Here we put forward a minimal structure-based model, the "multiple-basin model," that can directly be used for molecular dynamics simulation of even very large biomolecular systems so long as the endpoints of the conformational change are known. We investigate the model by simulating large-scale motions of four proteins: glutamine-binding protein, S100A6, dihydrofolate reductase, and HIV-1 protease. The mechanisms of conformational transition depend on the protein basin topologies and change with temperature near the folding transition. The conformational transition rate varies linearly with driving force over a fairly large range. This linearity appears to be a consequence of partial unfolding during the conformational transition.

Carrier Proteins↗
PubMed2006

Folding-based molecular simulations reveal mechanisms of the rotary motor F1-ATPase.

Biomolecular machines fulfill their function through large conformational changes that typically occur on the millisecond time scale or longer. Conventional atomistic simulations can only reach microseconds at the moment. Here, extending the minimalist model developed for protein folding, we propose the "switching Gō model" and use it to simulate the rotary motion of ATP-driven molecular motor F(1)-ATPase. The simulation recovers the unidirectional 120 degrees rotation of the gamma-subunit, the rotor. The rotation was induced solely by steric repulsion from the alpha(3)beta(3) subunits, the stator, which undergoes conformation changes during ATP hydrolysis. In silico alanine mutagenesis further elucidated which residues play specific roles in the rotation. Finally, regarding the mechanochemical coupling scheme, we found that the tri-site model does not lead to successful rotation but that the always-bi-site model produces approximately 30 degrees and approximately 90 degrees substeps, perfectly in accord with experiments. In the always-bi-site model, the number of sites occupied by nucleotides is always two during the hydrolysis cycle. This study opens up an avenue of simulating functional dynamics of huge biomolecules that occur on the millisecond time scales involving large-amplitude conformational change.

Catalysis↗
PubMed2006

Shaping up the protein folding funnel by local interaction: lesson from a structure prediction study.

Predicting protein tertiary structure by folding-like simulations is one of the most stringent tests of how much we understand the principle of protein folding. Currently, the most successful method for folding-based structure prediction is the fragment assembly (FA) method. Here, we address why the FA method is so successful and its lesson for the folding problem. To do so, using the FA method, we designed a structure prediction test of "chimera proteins." In the chimera proteins, local structural preference is specific to the target sequences, whereas nonlocal interactions are only sequence-independent compaction forces. We find that these chimera proteins can find the native folds of the intact sequences with high probability indicating dominant roles of the local interactions. We further explore roles of local structural preference by exact calculation of the HP lattice model of proteins. From these results, we suggest principles of protein folding: For small proteins, compact structures that are fully compatible with local structural preference are few, one of which is the native fold. These local biases shape up the funnel-like energy landscape.

Computer Simulation↗
PubMed2006

Protein grabs a ligand by extending anchor residues: molecular simulation for Ca2+ binding to calmodulin loop.

The structural difference in proteins between unbound and bound forms directly suggests the importance of the conformational plasticity of proteins. However, pathways that connect two-end structures and how they are coupled to the binding reaction are not well understood at atomic resolution. Here, we analyzed the free-energy landscape, explicitly taking into account coupling between binding and conformational change by performing atomistic molecular dynamics simulations for Ca2+ binding to a calmodulin loop. Using the AMBER force field with explicit water solvent, we conducted umbrella sampling for the free-energy surface and steered molecular dynamics for the pathway search. We found that, at an early stage of binding, some key residue side chains extend their "arms" to catch Ca2+ and, after catching, they carry the Ca2+ to the center of the binding pocket. This grabbing motion resulted in smooth and stepwise exchange in coordination partners of Ca2+ from water oxygen to atoms in the calmodulin loop. The key residue that first caught the ion was one of the two acidic residues, which are highly conserved. In the pathway simulations, different pathways were observed between binding and dissociation reactions: The former was more diverse than the latter.

Amino Acids↗
PubMed2006

SimFold energy function for de novo protein structure prediction: consensus with Rosetta.

Predicting protein tertiary structures by in silico folding is still very difficult for proteins that have new folds. Here, we developed a coarse-grained energy function, SimFold, for de novo structure prediction, performed a benchmark test of prediction with fragment assembly simulations for 38 test proteins, and proposed consensus prediction with Rosetta. The SimFold energy consists of many terms that take into account solvent-induced effects on the basis of physicochemical consideration. In the benchmark test, SimFold succeeded in predicting native structures within 6.5 A for 12 of 38 proteins; this success rate was the same as that by the publicly available version of Rosetta (ab initio version 1.2) run with default parameters. We investigated which energy terms in SimFold contribute to structure prediction performance, finding that the hydrophobic interaction is the most crucial for the prediction, whereas other sequence-specific terms have weak but positive roles. In the benchmark, well-predicted proteins by SimFold and by Rosetta were not the same for 5 of 12 proteins, which led us to introduce consensus prediction. With combined decoys, we succeeded in prediction for 16 proteins, four more than SimFold or Rosetta separately. For each of 38 proteins, structural ensembles generated by SimFold and by Rosetta were qualitatively compared by mapping sampled structural space onto two dimensions. For proteins of which one of the two methods succeeded and the other failed in prediction, the former had a less scattered ensemble located around the native. For proteins of which both methods succeeded in prediction, often two ensembles were mixed up.

Alanine↗
PubMed2005

Applying a grid technology to protein structure predictor "ROKKY".

This paper describes a sub-project of BioGrid project called "HTC (High Throughput Computing) group." Generally, a protein structure prediction which requires large amount of computational resources is done by trial-and-error method. HTC group have been developing a high throughput computing system with a flexible workflow handling mechanism for a protein structure prediction. In this paper, we show how to apply our high throughput computing system to the protein structure predictor called "ROKKY."

Computational Biology↗
PubMed2004

Roles of physical interactions in determining protein-folding mechanisms: molecular simulation of protein G and alpha spectrin SH3.

Protein-folding mechanisms of two small globular proteins, IgG binding domain of protein G and alpha spectrin SH3 domain are investigated via Brownian dynamics simulations with a model made of coarse-grained physical energy functions responsible for sequence-specific interactions and weak Gō-like energies. The folding pathways of alpha spectrin SH3 are known to be mainly controlled by the native topology, while protein G folding is anticipated to be more sensitive to the sequence-specific effects than native topology. We found in the folding of protein G that the C terminal beta hairpin is formed earlier and is rigid, once ordered, in the presence of an intact C terminal turn. The alpha helix is found to exhibit repeated partial formations/deformations during folding and to be stabilized via the tertiary contact with preformed beta sheets. This predicted scenario is fully consistent with experimental phi value data. Moreover, we found that the folding route is critically affected when the hydrophobic interaction is excluded from physical energy terms, suggesting that the hydrophobicity critically contributes to the folding propensity of protein G. For the folding of alpha spectrin SH3, we found that the distal beta hairpin and diverging turn are parts formed early, fully in harmony with previous results of simple Gō-like and experimental analysis, supporting that the folding route of SH3 domain is robust and coded by the native topology. The hybrid method provides useful tools for analyzing roles of physical interactions in determining folding mechanisms.

Bacterial Proteins↗
PubMed2004

Optimizing physical energy functions for protein folding.

We optimize a physical energy function for proteins with the use of the available structural database and perform three benchmark tests of the performance: (1) recognition of native structures in the background of predefined decoy sets of Levitt, (2) de novo structure prediction using fragment assembly sampling, and (3) molecular dynamics simulations. The energy parameter optimization is based on the energy landscape theory and uses a Monte Carlo search to find a set of parameters that seeks the largest ratio deltaE(s)/DeltaE for all proteins in a training set simultaneously. Here, deltaE(s) is the stability gap between the native and the average in the denatured states and DeltaE is the energy fluctuation among these states. Some of the energy parameters optimized are found to show significant correlation with experimentally observed quantities: (1) In the recognition test, the optimized function assigns the lowest energy to either the native or a near-native structure among many decoy structures for all the proteins studied. (2) Structure prediction with the fragment assembly sampling gives structure models with root mean square deviation less than 6 A in one of the top five cluster centers for five of six proteins studied. (3) Structure prediction using molecular dynamics simulation gives poorer performance, implying the importance of having a more precise description of local structures. The physical energy function solely inferred from a structural database neither utilizes sequence information from the family of the target nor the outcome of the secondary structure prediction but can produce the correct native fold for many small proteins.

Algorithms↗
PubMed2003

How protein thermodynamics and folding mechanisms are altered by the chaperonin cage: molecular simulations.

How the Escherichia coli GroEL/ES chaperonin assists folding of a substrate protein remains to be uncovered. Recently, it was suggested that confinement into the chaperonin cage itself can significantly accelerate folding of a substrate. Performing comprehensive molecular simulations of eight proteins confined into various sizes L of chaperonin-like cage, we explore how and to what extent protein thermodynamics and folding mechanisms are altered by the cage. We show that a substrate protein is remarkably stabilized by confinement; the estimated increase in denaturation temperature DeltaTf is as large as approximately 60 degrees C. For a protein of size R0, the stabilization DeltaTf scales as (R0/L)nu, where nu approximately 3, which is consistent with a mean field theory of polymer. We also found significant free energy cost of confining a protein, which increases with R0/L, indicating that the confinement requires external work provided by the chaperonin system. In kinetic study, we show the folding is accelerated in a modestly well confined case, which is consistent with a recent experimental result on ribulose-1,5-bisphosphate carboxylase-oxygenase folding and simulation results of a beta hairpin. Interestingly, the acceleration of folding is likely to be larger for a protein with more complex topology, as quantified by the contact order. We also show how ensemble of folding pathways are altered by the chaperonin-like cage calculating a variant of value used in the study of spontaneous folding.

Chaperonins↗
PubMed2003

Competition between protein folding and aggregation with molecular chaperones in crowded solutions: insight from mesoscopic simulations.

The living cell is inherently crowded with proteins and macromolecules. To avoid aggregation of denatured proteins in the living cell, molecular chaperones play important roles. Here we introduce a simple model to describe crowded protein solutions with chaperone-like species based on a dynamic density functional theory. As predicted by others, our simulations show that macromolecular crowding enhances the association of proteins and chaperones. However, when the intrinsic folding rate of the protein is slow, it is possible that crowding also enhances aggregation of proteins. The results of simulation suggest that, when the concentration of the crowding agent is as high as that in the cell, the association of the protein and unbound chaperone becomes correlated with the aggregation process, and that the protein-bound chaperones efficiently destroy the potential nuclei of aggregates and thus prevent the aggregation.

Biophysics↗
PubMed2003

De novo design of foldable proteins with smooth folding funnel: automated negative design and experimental verification.

De novo sequence design of foldable proteins provides a way of investigating principles of protein architecture. We performed fully automated sequence design for a target structure having a three-helix bundle topology and synthesized the designed sequences. Our design principle is different from the conventional approach, in that instead of optimizing interactions within the target structure, we design the global shape of the protein folding funnel. This includes automated implementation of negative design by explicitly requiring higher free energy of the denatured state. The designed sequences do not have significant similarity to those of any natural proteins. The NMR and CD spectroscopic data indicated that one designed sequence has a well-defined three-dimensional structure as well as alpha-helical content consistent with the target.

Amino Acid Sequence↗
PubMed2002

Effects of macromolecular crowding on protein folding and aggregation studied by density functional theory: dynamics.

Inside the living cell is inherently crowded with proteins and other macromolecules. Thus, it is indispensable to take into account various interactions between the protein and other macromolecules for thorough understanding of protein functions in cellular contexts. Here we focus on the excluded volume interaction imposed on the protein by surrounding macromolecules or "crowding agents." We have presented a theoretical framework for describing equilibrium properties of proteins in crowded solutions [A. R. Kinjo and S. Takada, Phys. Rev. E (to be published)]. In the present paper, we extend the theory to describe nonequilibrium properties of proteins in crowded solutions. Dynamics simulations exhibit qualitatively different morphologies depending on the aggregating conditions, and it was found that macromolecular crowding accelerates the onset of aggregation while stabilizing the native protein in the quasiuniform phase before the onset of aggregation. It is also observed, however, that the aggregation may be kinetically inhibited in highly crowded conditions. The effects of crowding on folding and unfolding of proteins are also examined, and the results suggest that fast folding is an important factor in preventing aggregation of denatured proteins.

Biophysical Phenomena↗
PubMed2002

Effects of macromolecular crowding on protein folding and aggregation studied by density functional theory: statics.

Proteins are neither purified nor diluted inside the living cell. Thus it is indispensable to take into account various interactions between the protein of interest and other macromolecules for understanding the properties of proteins in physiological conditions. Here we focus on excluded volume interactions which are omnipresent in dense or crowded solutions of proteins and macromolecules or "crowding agents." A protein solution with macromolecular crowding agents is modeled by means of a density functional theory. Effects of macromolecular crowding on protein aggregation and stability are investigated in particular. Phase diagrams are obtained in various parameter spaces by solving the equation of state. Two generic features are found: the addition of the crowding agent (1) enhances the aggregation of the denatured proteins, and (2) stabilizes the native protein unless the aggregation occurs. The present theory is qualitatively in good agreement with experimental observations and unifies previous theories regarding the crowding effects on protein stability and aggregation.

Journal Article↗