Biomedical subjects

Osamu Miyashita

Publications and source records attributed to Osamu Miyashita.

10 recordsLinked to original sources

PubMed2006

Conformational transitions of adenylate kinase: switching by cracking.

Conformational heterogeneity in proteins is known to often be the key to their function. We present a coarse grained model to explore the interplay between protein structure, folding and function which is applicable to allosteric or non-allosteric proteins. We employ the model to study the detailed mechanism of the reversible conformational transition of Adenylate Kinase (AKE) between the open to the closed conformation, a reaction that is crucial to the protein's catalytic function. We directly observe high strain energy which appears to be correlated with localized unfolding during the functional transition. This work also demonstrates that competing native interactions from the open and closed form can account for the large conformational transitions in AKE. We further characterize the conformational transitions with a new measure Phi(Func), and demonstrate that local unfolding may be due, in part, to competing intra-protein interactions.

Adenylate Kinase↗

PubMed2006

Exploring biomolecular machines: energy landscape control of biological reactions.

For almost 15 years, our Pathway model has been the most powerful model in terms of predicting the tunnelling mechanism for electron transfer (ET) in biological systems, particularly proteins. Going beyond the conventional Pathway models, we have generalized our method to understand how protein dynamics modulate not only the Franck-Condon factor, but also the tunnelling matrix element. We have demonstrated that when interference among pathways modulates the electron tunnelling interactions in proteins (particularly destructive interference), dynamical effects are of critical importance. Tunnelling can be controlled by protein conformations that lie far from equilibrium-those that minimize the effect of destructive interference during tunnelling, for example. In the opposite regime, electron tunnelling is mediated by one (or a few) constructively interfering pathway tubes and dynamical effects are modest. This new mechanism for dynamical modulation of the ET rate has been able to explain and/or predict several rates that were later confirmed by experiment. However, thermal fluctuations can also affect these molecular machines in many other ways. For example, we show how global transformations, which control protein functions such as allostery, may involve large-scale motion and possibly partial unfolding during the reaction event.

Azurin↗

PubMed2005

Coupled motions in the SH2 and kinase domains of Csk control Src phosphorylation.

The C-terminal Src kinase (Csk) phosphorylates and down-regulates Src family tyrosine kinases. The Csk-binding protein (Cbp) localizes Csk close to its substrates at the plasma membrane, and increases the specific activity of the kinase. To investigate this long-range catalytic effect, the phosphorylation of Src and the conformation of Csk were investigated in the presence of a high-affinity phosphopeptide derived from Cbp. This peptide binds tightly to the SH2 domain and enhances Src recognition (lowers K(m)) by increasing the apparent phosphoryl transfer rate in the Csk active site, a phenomenon detected in rapid quench flow experiments. Previous studies demonstrated that the regulation of Csk activity is linked to conformational changes in the enzyme that can be probed with hydrogen-deuterium exchange methods. We show that the Cbp peptide impacts deuterium incorporation into its binding partner (the SH2 domain), and into the SH2-kinase linker and several sequences in the kinase domain, including the glycine-rich loop in the active site. These findings, along with computational data from normal mode analyses, suggest that the SH2 domain moves in a cantilever fashion with respect to the small lobe of the kinase domain, ordering the active site for catalysis. The binding of a small Cbp-derived peptide to the SH2 domain of Csk modifies these motions, enhancing Src recognition.

Adaptor Proteins, Signal Transducing↗

PubMed2005

Interprotein electron transfer from cytochrome c2 to photosynthetic reaction center: tunneling across an aqueous interface.

Interprotein electron transfer (ET) reactions play an important role in biological energy conversion processes. One of these reactions, the ET between cytochrome c(2) (cyt) and reaction center from photosynthetic bacteria, is the focus of this theoretical study. The changes in the ET rate constant at fixed distances during the association process were calculated as the cyt moved from the electrostatically stabilized encounter complex to the bound state having short range van der Waals contacts in the tunneling region. Multiple conformations of the protein were generated by molecular dynamics simulations including explicit water molecules. For each of these conformations, the ET rate was calculated by using the Pathways model. The ET rate increased smoothly as the cyt approached from the encounter complex to the bound state, with a tunneling decay factor beta = 1.1 A(-1). This relatively efficient coupling between redox centers is due to the ability of interfacial water molecules to form multiple strong hydrogen bonding pathways connecting tunneling pathways on the surfaces of the two proteins. The ET rate determined for the encounter complex ensemble of states is only about a factor of 100 slower than that of the bound state (tau = 100 micros, compared with 1 micros), because of fluctuations of the cyt within the encounter complex ensemble through configurations having strong tunneling pathways. The ET rate for the encounter complex is in agreement with rates observed in mutant reaction centers modified to remove shortrange hydrophobic interactions, suggesting that in this case, ET occurs within the solvent-separated, electrostatically stabilized encounter complex.

Cytochromes c2↗

PubMed2005

Simple energy landscape model for the kinetics of functional transitions in proteins.

It is evident that protein conformational transitions play important roles in biological machinery; however, detailed pictures of these transition processes capable of making kinetic prediction are not yet available. For a full description of these transitions, we first need to describe kinematically movements between stable states. Then, more importantly, a free energy profile associated with the conformational change needs to be obtained. Recently, a new model to describe the energy landscape of protein conformational changes was applied to the conformational transition of adenylate kinase. In this model, the conformational change coupled to the ligand binding is described as a switching between two energy surfaces that correspond to ligand bound and unbound states. The nonlinearity of the protein conformational changes is described through an iterative usage of normal mode calculations. In addition, another kind of nonlinearity enters the dynamics of the conformational transitions due to cracking, or partial unfolding, which may occur during the conformational transitions. The consequences of this theoretical model are explored in greater detail. An improved model for the cracking that includes the cooperativity of the partial unfolding in analogy to nucleation is introduced.

Allosteric Site↗

PubMed2004

Transition state and encounter complex for fast association of cytochrome c2 with bacterial reaction center.

Electrostatic interactions strongly enhance the electron transfer reaction between cytochrome (Cyt) c(2) and reaction center (RC) from photosynthetic bacteria, yielding a second-order rate constant, k(2) approximately 10(9) s(-1).M(-1), close to the diffusion limit. The proposed mechanism involves an encounter complex (EC) stabilized by electrostatic interactions, followed by a transition state (TS), leading to the bound complex active in electron transfer. The effect of electrostatic interactions was previously studied by Tetreault et al. [Tetreault, M., Cusanovich, M., Meyer, T., Axelrod, H. & Okamura, M. Y. (2002) Biochemistry 41, 5807-5815] by measuring k(2) for RC and Cyt molecules with modified charged residues at the binding interface. The present work is a computational analysis of this kinetic study to determine the ensemble of configurations of the TS and EC. Changes in the TS energies due to different mutations were compared with differences in the calculated electrostatic energies for a wide range of Cyt/RC configurations. The TS ensemble, obtained from structures having the highest correlation coefficients in the comparison with experimental data, has the Cyt displaced by approximately 10 A from its position in x-ray crystal structure, close to the average position of the EC ensemble, with strong electrostatic interactions between Cyt on the M subunit side of the RC surface. The heme of the Cyt is oriented toward Tyr L162 on the RC, the tunneling contact in the bound final state on the RC. The similarity between the structures of the EC, TS, and bound state can account for the rapid rate of association responsible for fast diffusion-controlled electron transfer.

Cytochromes c2↗

PubMed2004

Dynamic coupling between the SH2 domain and active site of the COOH terminal Src kinase, Csk.

The SH2 domain is required for high catalytic activity in the COOH-terminal Src kinase (Csk). Previous solution studies suggest that a short peptide sequence, the SH2-kinase linker, provides a functional connection between the active site and the distal SH2 domain that could underlie this catalytic phenomenon. Substitutions in Phe183 (tyrosine, alanine, and glycine), a critical hydrophobic residue in the linker, result in large decreases in substrate turnover and large increases in the K(m) for ATP. Indeed, F183G possesses kinetic parameters that are similar to that for a truncated form of Csk lacking the SH2 domain, suggesting that a single mutation disrupts communication between this domain and the active site. Based on equilibrium and stopped-flow fluorescence experiments, the elevated K(m) values for the mutants are due to changes in the rates of phosphoryl transfer and not to reduced ATP-binding affinities. Based on hydrogen-deuterium exchange experiments, glycine substitution reduces flexibility in several polypeptide regions in Csk, tyrosine substitution increases flexibility, and alanine substitution leads to mixed effects compared to wild-type. Normal mode analysis indicates that Phe183 and its environment are under strain, a theoretical finding that supports the results of mutations. Overall, the data indicate that domain-domain interactions, controlled through the SH2-kinase linker, provide a dynamic balance within the Csk framework that is ideal for efficient phosphoryl transfer in the active site.

Adenosine Triphosphate↗

PubMed2004

Flexible multi-scale fitting of atomic structures into low-resolution electron density maps with elastic network normal mode analysis.

A novel method is presented for the quantitative flexible docking of high-resolution structure into low-resolution maps of macromolecular complexes from electron microscopy. This method uses a linear combination of low-frequency normal modes from elastic network description of the molecular framework in an iterative manner to deform the structure optimally to conform to the low-resolution electron density map. The methodology utilizes gradient following techniques in collective normal modes to locally optimize the overall correlation coefficient between computed and measured electron density. To evaluate the performance of our approach, several proteins, which undergo large conformational changes, have been studied. We demonstrate that refinement based on normal mode analysis provides an accurate and fast alternative for the flexible fitting of high-resolution structure into a low-resolution density map. Additionally, we show that lower resolution (multi-scale) structural models can be used for the normal mode searching in lieu of fully atomic models with little loss of overall accuracy.

Algorithms↗

PubMed2004

Normal mode based flexible fitting of high-resolution structure into low-resolution experimental data from cryo-EM.

A new method for the flexible fitting of high-resolution structures into low-resolution maps of macromolecular complexes from electron microscopy has been recently described in applications to simulated electron density maps. This method uses a linear combination of low-frequency normal modes in an iterative manner to deform the structure optimally to conform to the low-resolution electron density map. Gradient-following techniques in the coordinate space of collective normal modes are used to optimize the overall correlation coefficient between computed and measured electron densities. With this approach, multi-scale flexible fitting can be performed using all-atoms or Calpha atoms. In this paper, illustrative studies of normal mode based flexible fitting to experimental cryo-EM maps are presented for three different systems. Large, functionally relevant conformational changes for elongation factor G bound to the ribosome, Escherichia coli RNA polymerase and cowpea chlorotic mottle virus are elucidated as the result of the application of NMFF from high-resolution structures to cryo-electron microscopy maps.

Cryoelectron Microscopy↗

PubMed2003

Continuum electrostatic model for the binding of cytochrome c2 to the photosynthetic reaction center from Rhodobacter sphaeroides.

Electrostatic interactions are important for protein-protein association. In this study, we examined the electrostatic interactions between two proteins, cytochrome c(2) (cyt c(2)) and the reaction center (RC) from the photosynthetic bacterium Rhodobacter sphaeroides, that function in intermolecular electron transfer in photosynthesis. Electrostatic contributions to the binding energy for the cyt c(2)-RC complex were calculated using continuum electrostatic methods based on the recent cocrystal structure [Axelrod, H. L., et al. (2002) J. Mol. Biol. 319, 501-515]. Calculated changes in binding energy due to mutations of charged interface residues agreed with experimental results for a protein dielectric constant epsilon(in) of 10. However, the electrostatic contribution to the binding energy for the complex was close to zero due to unfavorable desolvation energies that compensate for the favorable Coulomb attraction. The electrostatic energy calculated as a function of displacement of the cyt c(2) from the bound position showed a shallow minimum at a position near but displaced from the cocrystal configuration. These results show that although electrostatic steering is present, other short-range interactions must be present to contribute to the binding energy and to determine the structure of the complex. Calculations made to model the experimental data on association rates indicate a solvent-separated transition state for binding in which the cyt c(2) is displaced approximately 8 A above its position in the bound complex. These results are consistent with a two-step model for protein association: electrostatic docking of the cyt c(2) followed by desolvation to form short-range van der Waals contacts for rapid electron transfer.

Binding Sites↗