Takeshi Kikuchi

Absence of binding between the human transferrin receptor and the transferrin complex of biological toxic trace element, aluminum, because of an incomplete open/closed form of the complex.

Human transferrin (Tf) very tightly binds two ferric ions to deliver iron to cells. Fe(III)2Tf (Fe2Tf) binds to the Tf receptor (TfR) at pH 7.4; however, iron-free Tf (apoTf) does not. Iron uptake is facilitated by endocytosis of the Fe2Tf–TfR complex. Tf can also bind aluminum ions, which cause toxic effects and are associated with many diseases. Since Al(III)2Tf (Al2Tf) does not bind to TfR, the uptake of aluminum by the cells does not occur through a TfR-mediated pathway. We have studied the absence of binding between Al2Tf and TfR by investigating the physicochemical characteristics of apoTf, Al2Tf, Fe2Tf, and TfR. The hydrodynamic radius of 38.8xa0Å for Al2Tf obtained by dynamic light scattering was between that of 42.6xa0Å for apoTf and 37.2xa0Å for Fe2Tf. The ζ potential of −11.3xa0mV for Al2Tf measured by capillary electrophoresis was close to −11.2xa0mV for apoTf as compared to −11.9xa0mV for Fe2Tf, indicating that the Al2Tf surface had a relatively scarce negative charge as the apoTf surface had. These results demonstrated that the structure of Al2Tf was a trade-off between the closed and open forms of Fe2Tf and apoTf, respectively. Consequently, it is suggested that Al2Tf cannot form specific ionic interresidual interactions, such as those formed by Fe2Tf, to bind to TfR, resulting in impossible complex formation between Al2Tf and TfR.

Prediction of folding pathway and kinetics among plant hemoglobins using an average distance map method

Computational methods, such as the ADM (average distance map) method, have been developed to predict folding of homologous proteins. In this work we used the ADM method to predict the folding pathway and kinetics among selected plant nonsymbiotic (nsHb), symbiotic (Lb), and truncated (tHb) hemoglobins (Hbs). Results predicted that (1) folding of plant Hbs occurs throughout the formation of compact folding modules mostly formed by helices A, B, and C, and E, F, G, and H (folding modules A/C and E/H, respectively), and (2) primitive (moss) nsHbs fold in the C→N direction, evolved (monocot and dicot) nsHbs fold either in the C→N or N→C direction, and Lbs and plant tHbs fold in the C→N direction. We also predicted relative folding rates of plant Hbs from qualitative analyses of the stability of subdomains and classified plant Hbs into fast and moderate folding. ADM analysis of nsHbs predicted that prehelix A plays a role during folding of the N‐terminal domain of Ceratodon nsHb, and that CD‐loop plays a role in folding of primitive (Physcomitrella and Ceratodon) but not evolved nsHbs. Modeling of the rice Hb1 A/C and E/H modules showed that module E/H overlaps to the Mycobacterium tuberculosis HbO two‐on‐two folding. This observation suggests that module E/H is an ancient tertiary structure in plant Hbs. Proteins 2005.

Engineering the substrate specificity of Alcaligenes D-aminoacylase useful for the production of D-amino acids by optical resolution.

D-Aminoacylase from Alcaligenes xylosoxydans subsp. xylosoxydans A-6 (AxD-NAase) offers a novel biotechnological application, the production of D-amino acid from the racemic mixture of N-acyl-DL-amino acids. However, its substrate specificity is biased toward certain N-acyl-D-amino acids. To construct mutant AxD-NAases with substrate specificities different from those of wild-type enzyme, the substrate recognition site of the AxD-NAase was rationally manipulated based on computational structural analysis and comparison of its primary structure with other D-aminoacylases with distinct substrate specificities. Mutations of amino acid residues, Phe191, Leu298, Tyr344, and Met346, which interact with the side chain of the substrate, induced marked changes in activities toward each substrate. For example, the catalytic efficiency (k(cat)/K(m)) of mutant F191W toward N-acetyl-D-Trp and N-acetyl-D-Ala was enhanced by 15.6- and 1.5-folds, respectively, compared with that of the wild-type enzyme, and the catalytic efficiency (k(cat)/K(m)) of mutant L298A toward N-acetyl-D-Trp was enhanced by 4.4-folds compared with that of the wild-type enzyme. Other enzymatic properties of both mutants, such as pH and temperature dependence, were the same as those of the wild-type enzyme. The F191W mutant in particular is considered to be useful for the enzymatic production of D-Trp which is an important building block of some therapeutic drugs.

Computational Structure Models of Apo and Diferric Transferrin–Transferrin Receptor Complexes

Complexation of transferrin (Tf) and its receptor (TfR) is an essential event for iron uptake by the cell. Much data has been accumulated regarding Tf-TfR complexation, such as results from mutagenesis. We created 3D structural models of apo-human Tf-TfR (apoTf-TfR) and Fe(III)2Tf-TfR (Fe2Tf-TfR) complexes by computational rigid body refinement. The models are consistent with published mutagenesis experiments. In our models, the C-lobes of apoTf and Fe2Tf bind to the helical domain of TfR, and the N-lobes are sandwiched between the ectodomain of TfR and the cell membrane as previously reported. Further, the molecules of apoTf and Fe2Tf are not forced to undergo large conformational changes upon complexation. The creation of the models led a new and important finding that a residue of TfR, R651, which is called a hot spot for Tf-TfR binding, interacts with Tf E385 when either apoTf or Fe2Tf bind to TfR. The models rationally interpret the iron release from Fe2Tf-TfR upon acidification, dissociation of apoTf from TfR at slightly alkaline pH, and metal specific recognition of TfR.

Analysis of 3D structural differences in the IgG-binding domains based on the interresidue average-distance statistics

It is well-known that the IgG-binding domain from staphylococcal protein A folds into a 3α helix bundle structure, while the IgG-binding domain of streptococcal protein G forms an (αxa0+xa0β) structure. Recently, He et al. (Biochemistry 44:14055–14061, 2005) made mutants of these proteins from the wild types of protein A and protein G strains. These mutants are referred to as protein A219 and protein G311, and it was showed that these two mutants have different 3D structures, i.e., the 3α helix bundle structure and the (αxa0+xa0β) structure, respectively, despite the high sequence identity (59%). The purpose of our study was to clarify how such 3D structural differences are coded in the sequences with high homology. To address this problem, we introduce a predicted contact map constructed based on the interresidue average-distance statistics for prediction of folding properties of a protein. We refer to this map as an average distance map (ADM). Furthermore, the statistics of interresidue distances can be converted to an effective interresidue potential. We calculated the contact frequency of each residue of a protein in random conformations with this effective interresidue potential, and then we obtained values similar to ϕ values. We refer to this contact frequency of each residue as a p(μ) value. The comparison of the p(μ) values to the ϕ values for a protein suggests that p(μ) values reveal the information on the folding initiation site. Using these techniques, we try to extract the information on the difference in the 3D structures of protein A219 and protein G311 coded in their amino acid sequences in the present work. The results show that the ADM analyses and the p(μ) value analyses predict the information of folding initiation sites, which can be used to detect the 3D difference in both proteins.

Sequence analysis on the information of folding initiation segments in ferredoxin-like fold proteins

BackgroundWhile some studies have shown that the 3D protein structures are more conservative than their amino acid sequences, other experimental studies have shown that even if two proteins share the same topology, they may have different folding pathways. There are many studies investigating this issue with molecular dynamics or Go-like model simulations, however, one should be able to obtain the same information by analyzing the proteins’ amino acid sequences, if the sequences contain all the information about the 3D structures. In this study, we use information about protein sequences to predict the location of their folding segments. We focus on proteins with a ferredoxin-like fold, which has a characteristic topology. Some of these proteins have different folding segments.ResultsDespite the simplicity of our methods, we are able to correctly determine the experimentally identified folding segments by predicting the location of the compact regions considered to play an important role in structural formation. We also apply our sequence analyses to some homologues of each protein and confirm that there are highly conserved folding segments despite the homologues’ sequence diversity. These homologues have similar folding segments even though the homology of two proteins’ sequences is not so high.ConclusionOur analyses have proven useful for investigating the common or different folding features of the proteins studied.

Incorporating into a Cα Go model the effects of geometrical restriction on Cα atoms caused by side chain orientations

Coarse‐grained Go models have been widely used for studying protein‐folding mechanisms. Despite the simplicity of the model, these can reproduce the essential features of the folding process of a protein. However, it is also known that side chains significantly contribute to the folding mechanism. Hence, it is desirable to incorporate the side chain effects into a coarse‐grained Go model. In this study, to distinguish the effects of side chain orientation and to understand how these effects contribute to folding mechanisms, we incorporate into a Cα Go model not only heterogeneous contact energies but also geometrical restraints around two Cα atoms in contact with each other. We confirm that the heterogeneity of contact energies governs the folding pathway of a protein and that the geometric constraints attributed to side chains reproduce cooperative transitions in folding. Proteins 2013; 81:1434–1445.

Analyses of the folding properties of ferredoxin-like fold proteins by means of a coarse-grained Gō model: Relationship between the free energy profiles and folding cores

The folding mechanisms of proteins with multi‐state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on‐going research issues. In this study, we investigate ferredoxin‐like fold proteins which have a simple topology and multi‐state transitions. We analyze the folding processes by means of a coarse‐grained Gō model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high‐free‐energy intermediate state, and ADA2h and S6, which fold into the native structure through two‐state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54‐55, are reproduced and compared with experimental observations. We show that the ferredoxin‐like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two‐state mechanism. Proteins 2014; 82:954–965.

In silico study on the substrate binding manner in human myo-inositol monophosphatase 2

The human IMPA2 gene encoding myo-inositol monophosphatase 2 is highly implicated with bipolar disorder but the substrates and the reaction mechanism of myo-inositol monophosphatase 2 have not been well elucidated.9 In the present study, we constructed 3D models of three- and two-Mg2+-ion bound myo-inositol monophosphatase 2, and studied substrate-binding manners using the docking program AutoDock3. The subsequent study showed that the three-metal-ion model could interact with myo-inositol monophosphates, as follows: The phosphate moiety coordinated three Mg2+ ions, and the inositol ring formed hydrogen bonds with the amino acids conserved in the family. Furthermore, the OH group vicinal to the phosphate group formed a hydrogen bond with a non-bridging oxygen atom of the phosphate. These interactions have been proposed as crucial for forming the transitional state, bipyramidal structure in the bovine myo-inositol monophosphatase. We therefore propose that the human myo-inositol monophosphatase 2 interacts with myo-inositol monophosphates in the three-metal-ion bound form, and proceeds the dephosphorylation through the three-metal-ion theory.

Topological and sequence information predict that foldons organize a partially overlapped and hierarchical structure

It has been suggested that proteins have substructures, called foldons, which can cooperatively fold into the native structure. However, several prior investigations define foldons in various ways, citing different foldon characteristics, thereby making the concept of a foldon ambiguous. In this study, we perform a Gō model simulation and analyze the characteristics of substructures that cooperatively fold into the native‐like structure. Although some results do not agree well with the experimental evidence due to the simplicity of our coarse‐grained model, our results strongly suggest that cooperatively folding units sometimes organize a partially overlapped and hierarchical structure. This view makes us easy to interpret some different proposal about the foldon as a difference of the hierarchical structure. On the basis of this finding, we present a new method to assign foldons and their hierarchy, using structural and sequence information. The results show that the foldons assigned by our method correspond to the intermediate structures identified by some experimental techniques. The new method makes it easy to predict whether a protein folds sequentially into the native structure or whether some foldons fold into the native structure in parallel. Proteins 2015; 83:1900–1913.