Tadaomi Furuta

Minimum free energy path of ligand-induced transition in adenylate kinase.

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.

The Tryptophan Residue at the Active Site Tunnel Entrance of Trichoderma reesei Cellobiohydrolase Cel7A Is Important for Initiation of Degradation of Crystalline Cellulose

Background: Mutation of Trp-40 in the Cel7A cellobiohydrolase from Trichoderma reesei (TrCel7A) causes a loss of crystalline cellulose-degrading ability. Results: Mutant W40A showed reduced specific activity for crystalline cellulose and diffused the cellulose chain from the entrance of the active site tunnel. Conclusion: Trp-40 is essential for chain end loading to initiate processive hydrolysis of TrCel7A. Significance: The mechanisms of crystalline polysaccharide degradation are clarified. The glycoside hydrolase family 7 cellobiohydrolase Cel7A from Trichoderma reesei is one of the best studied cellulases with the ability to degrade highly crystalline cellulose. The catalytic domain and the cellulose-binding domain (CBD) are both necessary for full activity on crystalline substrates. Our previous high-speed atomic force microscopy studies showed that mutation of Trp-40 at the entrance of the catalytic tunnel drastically decreases the ability to degrade crystalline cellulose. Here, we examined the activities of the WT enzyme and mutant W40A (with and without the CBD) for various substrates. Evaluation and comparison of the specific activities of the enzymes (WT, W40A, and the corresponding catalytic subunits (WTcat and W40Acat)) adsorbed on crystalline cellulose indicated that Trp-40 is involved in recruiting individual substrate chains into the active site tunnel to initiate processive hydrolysis. This was supported by molecular dynamics simulation study, i.e. the reducing end glucose unit was effectively loaded into the active site of WTcat, but not into that of W40Acat, when the simulation was started from subsite −7. However, when similar simulations were carried out starting from subsite −5, both enzymes held the substrate for 50 ns, indicating that the major difference between WTcat and W40Acat is the length of the free chain end of the substrate required to allow initiation of processive movements; this also reflects the difference between crystalline and amorphous celluloses. The CBD is important for enhancing the enzyme population on crystalline substrate, but it also decreases the specific activity of the adsorbed enzyme, possibly by attaching the enzyme to non-optimal places on the cellulose surface and/or hindering processive hydrolysis.

Analysis of the structural and functional roles of coupling helices in the ATP-binding cassette transporter MsbA through enzyme assays and molecular dynamics simulations.

ATP-binding cassette (ABC) transporters are constructed from some common structural units: the highly conserved nucleotide-binding domains (NBDs), which work as a nucleotide-dependent engine for driving substrate transport, the diverse transmembrane domains (TMDs), which create the translocation pathway, and the coupling helices (CHs), which are located at the NBD-TMD interface. Although the CHs are believed to be essential for NBD-TMD communication, their roles remain unclear. In this study, we performed enzyme assays and molecular dynamics (MD) simulations of the ABC transporter MsbA and two MsbA mutants in which the amino acid residues of one of the CHs were mutated to alanines: (i) wild type (Wt), (ii) CH1 mutant (Mt1), and (iii) CH2 mutant (Mt2). The experiments show that the CH2 mutation decreases the ATPase activity (kcat) compared with that of the Wt (a decrease of 32%), and a nearly equal degree of decrease in the ATP binding affinity (Km) was observed for both Mt1 and Mt2. The MD simulations successfully accounted for several structural and dynamical origins for these experimental observations. In addition, on the basis of collective motion and morphing analyses, we propose that the reverse-rotational motions and noddinglike motions between the NBDs and TMDs are indispensable for the conformational transition between the inward- and outward-facing conformations. In particular, CH2 is significantly important for the occurrence of the noddinglike motion. These findings provide important insights into the structure-function relationship of ABC transporters.

Different inhibitory potency of febuxostat towards mammalian and bacterial xanthine oxidoreductases: insight from molecular dynamics

Febuxostat, a drug recently approved in the US, European Union and Japan for treatment of gout, inhibits xanthine oxidoreductase (XOR)-mediated generation of uric acid during purine catabolism. It inhibits bovine milk XOR with a Ki in the picomolar-order, but we found that it is a much weaker inhibitor of Rhodobacter capsulatus XOR, even though the substrate-binding pockets of mammalian and bacterial XOR are well-conserved as regards to catalytically important residues and three-dimensional structure, and both permit the inhibitor to be accommodated in the active site, as indicated by computational docking studies. To clarify the reason for the difference of inhibitory potency towards the two XORs, we performed molecular dynamics simulations. The results indicate that differences in mobility of hydrophobic residues that do not directly interact with the substrate account for the difference in inhibitory potency.

ATP-induced conformational changes of nucleotide-binding domains in an ABC transporter. Importance of the water-mediated entropic force.

ATP binding cassette (ABC) proteins belong to a superfamily of active transporters. Recent experimental and computational studies have shown that binding of ATP to the nucleotide binding domains (NBDs) of ABC proteins drives the dimerization of NBDs, which, in turn, causes large conformational changes within the transmembrane domains (TMDs). To elucidate the active substrate transport mechanism of ABC proteins, it is first necessary to understand how the NBD dimerization is driven by ATP binding. In this study, we selected MalKs (NBDs of a maltose transporter) as a representative NBD and calculated the free-energy change upon dimerization using molecular mechanics calculations combined with a statistical thermodynamic theory of liquids, as well as a method to calculate the translational, rotational, and vibrational entropy change. This combined method is applied to a large number of snapshot structures obtained from molecular dynamics simulations containing explicit water molecules. The results suggest that the NBD dimerization proceeds with a large gain of water entropy when ATP molecules bind to the NBDs. The energetic gain arising from direct NBD-NBD interactions is canceled by the dehydration penalty and the configurational-entropy loss. ATP hydrolysis induces a loss of the shape complementarity between the NBDs, which leads to the dissociation of the dimer, due to a decrease in the water-entropy gain and an increase in the configurational-entropy loss. This interpretation of the NBD dimerization mechanism in concert with ATP, especially focused on the water-mediated entropy force, is potentially applicable to a wide variety of the ABC transporters.

Disturbed biopterin and folate metabolism in the Qdpr-deficient mouse

Quinonoid dihydropteridine reductase (QDPR) catalyzes the regeneration of tetrahydrobiopterin (BH4), a cofactor for monoamine synthesis, phenylalanine hydroxylation and nitric oxide production. Here, we produced and analyzed a transgenic Qdpr −/− mouse model. Unexpectedly, the BH4 contents in the Qdpr −/− mice were not decreased and even increased in some tissues, whereas those of the oxidized form dihydrobiopterin (BH2) were significantly increased. We demonstrated that unlike the wild‐type mice, dihydrofolate reductase regenerated BH4 from BH2 in the mutants. Furthermore, we revealed wide alterations in folate‐associated metabolism in the Qdpr −/− mice, which suggests an interconnection between folate and biopterin metabolism in the transgenic mouse model.

A Fluorescent Protein Scaffold for Presenting Structurally Constrained Peptides Provides an Effective Screening System to Identify High Affinity Target-Binding Peptides

Peptides that have high affinity for target molecules on the surface of cancer cells are crucial for the development of targeted cancer therapies. However, unstructured peptides often fail to bind their target molecules with high affinity. To efficiently identify high-affinity target-binding peptides, we have constructed a fluorescent protein scaffold, designated gFPS, in which structurally constrained peptides are integrated at residues K131–L137 of superfolder green fluorescent protein. Molecular dynamics simulation supported the suitability of this site for presentation of exogenous peptides with a constrained structure. gFPS can present 4 to 12 exogenous amino acids without a loss of fluorescence. When gFPSs presenting human epidermal growth factor receptor type 2 (HER2)-targeting peptides were added to the culture medium of HER2-expressing cells, we could easily identify the peptides with high HER2-affinity and -specificity based on gFPS fluorescence. In addition, gFPS could be expressed on the yeast cell surface and applied for a high-throughput screening. These results demonstrate that gFPS has the potential to serve as a powerful tool to improve screening of structurally constrained peptides that have a high target affinity, and suggest that it could expedite the one-step identification of clinically applicable cancer cell-binding peptides.

Kirkwood–Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies

Cosolvents, such as urea, affect protein folding and binding, and the solubility of solutes. The modeling of cosolvents has been facilitated significantly by the rigorous Kirkwood-Buff (KB) theory of solutions, which can describe structural thermodynamics over the entire composition range of aqueous cosolvent mixtures based only on the solution density and the KB integrals (KBIs), i.e., the net excess radial distribution functions from the bulk. Using KBIs to describe solution thermodynamics has given rise to a clear guideline that an accurate prediction of KBIs is equivalent to accurate modeling of cosolvents. Taking urea as an example, here we demonstrate that an improvement in the prediction of KBIs comes from an improved reproduction of high-level quantum chemical (QC) electrostatic potential and molecular pairwise interaction energies. This rational approach to the improvement of the KBI prediction stems from a comparison of existing force fields, AMOEBA, and the generalized AMBER force field, as well as the further optimization of the former to enable better agreement with QC interaction energies. Such improvements would pave the way toward a rational and systematic determination of the transferable force field parameters for a number of important small molecule cosolvents.

In Silico Chaperonin-Like Cycle Helps Folding of Proteins for Structure Prediction

Currently, one of the most serious problems in protein-folding simulations for de novo structure prediction is conformational sampling of medium-to-large proteins. In vivo, folding of these proteins is mediated by molecular chaperones. Inspired by the functions of chaperonins, we designed a simple chaperonin-like simulation protocol within the framework of the standard fragment assembly method: in our protocol, the strength of the hydrophobic interaction is periodically modulated to help the protein escape from misfolded structures. We tested this protocol for 38 proteins and found that, using a certain defined criterion of success, our method could successfully predict the native structures of 14 targets, whereas only those of 10 targets were successfully predicted using the standard protocol. In particular, for non-alpha-helical proteins, our method yielded significantly better predictions than the standard approach. This chaperonin-inspired protocol that enhanced de novo structure prediction using folding simulations may, in turn, provide new insights into the working principles underlying the chaperonin system.

The mechanism of nucleotide-binding domain dimerization in the intact maltose transporter as studied by all-atom molecular dynamics simulations

The Escherichia coli maltose transporter MalFGK2‐E belongs to the protein superfamily of ATP‐binding cassette (ABC) transporters. This protein is composed of heterodimeric transmembrane domains (TMDs) MalF and MalG, and the homodimeric nucleotide‐binding domains (NBDs) MalK2. In addition to the TMDs and NBDs, the periplasmic maltose binding protein MalE captures maltose and shuttle it to the transporter. In this study, we performed all‐atom molecular dynamics (MD) simulations on the maltose transporter and found that both the binding of MalE to the periplasmic side of the TMDs and binding of ATP to the MalK2 are necessary to facilitate the conformational change from the inward‐facing state to the occluded state, in which MalK2 is completely dimerized. MalE binding suppressed the fluctuation of the TMDs and MalF periplasmic region (MalF‐P2), and thus prevented the incorrect arrangement of the MalF C‐terminal (TM8) helix. Without MalE binding, the MalF TM8 helix showed a tendency to intrude into the substrate translocation pathway, hindering the closure of the MalK2. This observation is consistent with previous mutagenesis experimental results on MalF and provides a new point of view regarding the understanding of the substrate translocation mechanism of the maltose transporter.