Article Information
PubMed
Related Articles
- Can an Atomic Force Microscope Sequence DNA Using a Nanopore...Can an Atomic Force Microscope Sequence DNA Using a Nanopore?
Biophysical Journal, Volume 94, Issue 4, 15 February 2008, Pages 1233-1240
Shahid Qamar, Phil M. Williams and S.M. Lindsay
AbstractR. Bension has proposed that single molecules of DNA could be sequenced rapidly, in long sequential reads, by reading off the force required to pull a tightly fitting molecular ring over each base in turn using an atomic force microscope (AFM). We present molecular dynamics simulations that indicate that pulling DNA very rapidly (m/s) could generate large force peaks as each base is passed (∼1nN) with significant differences (∼0.5nN) between purine and pyrimidine. These speeds are six orders of magnitude faster than could be read out by a conventional AFM, and extending the calculations to accessible speeds using Kramers’ theory shows that thermal fluctuations dominate the process with the result that purine and pyrimidine cannot be distinguished with the pulling speeds attained by current AFM technology.
Abstract | Full Text | PDF (757 kb) - Cooperative hemoglobins: conserved fold, diverse quaternary ...Cooperative hemoglobins: conserved fold, diverse quaternary assemblies and allosteric mechanisms
Trends in Biochemical Sciences, Volume 26, Issue 5, 1 May 2001, Pages 297-304
William E Royer, James E Knapp, Kristen Strand and Holly A Heaslet
AbstractAssembly of hemoglobin subunits into cooperative complexes produces a remarkable variety of architectures, ranging in oligomeric state from dimers to complexes containing 144 hemoglobin subunits. Diverse stereochemical mechanisms for modulating ligand affinity through intersubunit interactions have been revealed from studies of three distinct hemoglobin assemblages. This mechanistic diversity, which occurs between assemblies of subunits that have the same fold, provides insight into the range of regulatory strategies that are available to protein molecules.
Abstract | Full Text | PDF (511 kb) - Structural and Dynamic Properties of the Homodimeric Hemoglo...Structural and Dynamic Properties of the Homodimeric Hemoglobin from Scapharca inaequivalvis Thr-72→Ile Mutant: Molecular Dynamics Simulation, Low Temperature Visible Absorption Spectroscopy, and Resonance Raman Spectroscopy Studies
Biophysical Journal, Volume 75, Issue 5, 1 November 1998, Pages 2489-2503
Mattia Falconi, Alessandro Desideri, Antonio Cupane, Maurizio Leone, Giovanni Ciccotti, Eric S. Peterson, Joel M. Friedman, Alessandra Gambacurta and Franca Ascoli
AbstractMolecular dynamics simulations, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy have been performed on a mutant of the homodimeric hemoglobin, where residue threonine 72, at the subunit interface, has been substituted by isoleucine. Molecular dynamics simulation indicates that in the Thr-72→Ile mutant several residues that have been shown to play a role in ligand binding fluctuate around orientations and distances similar to those observed in the x-ray structure of the CO derivative of the native hemoglobin, although the overall structure remains in the T state. Visible absorption spectroscopy data indicate that in the deoxy form the Soret band is less asymmetric in the mutant than in the native protein, suggesting a more planar heme structure; moreover, these data suggest a similar heme-solvent interaction in both the liganded and unliganded states of the mutant protein, at variance with that observed in the native protein. The “conformation sensitive” band III of the deoxy mutant protein is shifted to lower energy by >100cm with respect to the native one, about one-half of that observed in the low temperature photoproducts of both proteins, indicating a less polar or more hydrophobic heme environment. Resonance Raman spectroscopy data show a slight shift of the iron-proximal histidine stretching mode of the deoxy mutant toward lower frequency with respect to the native protein, which can be interpreted in terms of either a change in packing of the phenyl ring of Phe-97, as also observed from the simulation, or a loss of water in the heme pocket. In line with this latter interpretation, the number of water molecules that dynamically enters the intersubunit interface, as calculated by the molecular dynamics simulation, is lower in the mutant than in the native protein. The 10-ns photoproduct for the carbonmonoxy mutant derivative has a higher iron-proximal histidine stretching frequency than does the native protein. This suggests a subnanosecond relaxation that is slowed in the mutant, consistent with a stabilization of the R structure. Taken together, the molecular dynamics and the spectroscopic data indicate that the higher oxygen affinity displayed by the Thr-72→Ile mutant is mainly due to a local perturbation in the dimer interface that propagates to the heme region, perturbing the polarity of the heme environment and propionate interactions. These changes are consistent with a destabilization of the T state and a stabilization of the R state in the mutant relative to the native protein.
Abstract | Full Text | PDF (443 kb) - …more
Copyright © 2007 The Biophysical Society. All rights reserved.
Biophysical Journal, Volume 92, Issue 9, L85-L87, 1 May 2007
doi:10.1529/biophysj.106.101899
Biophysical Letters
A Milestoning Study of the Kinetics of an Allosteric Transition: Atomically Detailed Simulations of Deoxy Scapharca Hemoglobin
Ron Elber
, 
Department of Computer Science, Cornell University, Ithaca, New York
Address reprint requests and inquiries to Ron Elber, Tel.: 607-255-7416.Abstract
Atomically detailed simulations are used to compute the kinetics of the R-to-T transition in deoxy Scapharca hemoglobin. A computational approach called milestoning is utilized that combines 1), an efficient reaction path algorithm; and 2), a “fragment and glue” approach for classical trajectories. Milestoning computes the R-to-T transition kinetics on the microsecond timescale based on atomically detailed trajectories that rarely exceed a nanosecond. Eleven reference hypersurfaces (milestones) are constructed along the reaction coordinate, which is computed with a global path optimization algorithm. Two-hundred classical trajectories are calculated for each of the milestones to collect local distributions of first passage times. These local distributions are used in a non-Markovian theory to compute the overall timescale. Exponential enrichment of reactive trajectories, an important component of the milestoning approach, makes these calculations possible. The overall timescale of the reaction is estimated as 10±9μs, in accord with available experimental data. The barrier is not sharp and is spread over four milestones. Even after the most significant structural changes are completed (phenylalanine F4 ring flips), highly collective and activated motions continue. The calculations suggest an additional late free energy barrier.One of the important mechanisms of response to environmental changes in biology is that of allostery. Adjustments in protein packing and assembly affect basic protein activity in reaction to changes of ligand concentration, products, pH, and other effectors. Naturally these processes have attracted considerable experimental and theoretical interest. Most investigations focus on human (tetrameric) hemoglobin HbA for which we now have vast data, models, and some remaining questions due to the complexity of the system 1. Interestingly, a dimeric hemoglobin (Scapharca hemoglobin (HbI)) 2 shows a profoundly different allosteric mechanism than HbA while maintaining a tertiary fold typical of a globin. The allosteric transition in HbI is significantly simpler than the transition in HbA and more accessible to calculation and experiment. The most apparent structural changes observed in the allosteric transition of HbI are flips of two amino-acid rings and adjustments of water molecules at the interface of the homodimer. No large changes in subunit orientation are observed (in contrast to HbA), allowing for a recent study of the transition using time-dependent x-ray diffraction 3.
The time-resolved experiments provide a detailed picture of the progress of the reaction. However, the connection of the structural images to energetic, kinetic, and thermodynamic properties of the system is not obvious. Atomically detailed computer simulations hold the promise to bridge the gap. Unfortunately, computer simulations are difficult to perform in this case because: 1), the experimentally measured timescale is long (microseconds); and 2), there is no clear and narrow transition state or barrier domain. The first observation excludes the application of straightforward molecular dynamics simulations. Even if a few microsecond trajectories could be computed, they would not be sufficient to estimate the kinetics. The last observation makes it difficult to apply approaches appropriate for activated processes 4,5. These approaches require the actual transition time to be short (the short-time trajectories can be rare, leading to long timescale).
For the study of long time kinetics in complex systems with potentially broad spatial barriers and intermediates (experiment suggests a nanosecond intermediate in HbI 3), the method of milestoning 6,7 seems appropriate. There are three steps to a milestoning calculation:
Step 1. Sampling from local-equilibrium ensembles in cross-sections (milestones) between the reactant and product (Fig. 1);
Step 2. Computing trajectories between nearby milestones and estimating local first-passage time distributions (LFPTD) between them; and
Step 3. Computing the overall time dependence of the process using all the LFPTD-s in a stochastic model.
Figure 1
A schematic two-dimensional energy landscape illustrating the milestoning approach. A reaction coordinate, q, is computed between the two minima (R and P). A milestone i is a hyperplane perpendicular to the reaction coordinate at qi. Also shown is a sampled trajectory, which starts at milestone i and terminates at milestone i+1.Step 2. Computing trajectories between nearby milestones and estimating local first-passage time distributions (LFPTD) between them; and
Step 3. Computing the overall time dependence of the process using all the LFPTD-s in a stochastic model.
The two end structures (the structures of the R and T states of HbI) were prepared from 1HBI 8 and 4SDH 9 of the protein databank 10 by building missing hydrogens, and minimizing the energy for 50 steps. A discrete reaction coordinate with nine intermediates was computed with the self-penalty walk method 11 of the program MOIL 12. A milestone (hyperplane) is defined by its normal and a position in the plane, Xi. The calculation of the reaction coordinate includes only the protein and crystal waters. A movie of the computed reaction coordinate is included in the Supplementary Material .
Each configuration along the reaction coordinate was then solvated and equilibrated in a periodic box of water. The total number of atoms in the system was ∼17,100 (the number of water molecules varies slightly for different water boxes along the reaction coordinate). In the next step we sample configurations at 300K in each milestone by calculating a trajectory constrained to the hyperplane. Coordinates are saved every picosecond of a 200-ps trajectory, providing a total of 200 initial conditions for Step 2. Straightforward molecular dynamics trajectories were initiated at each milestone using the previously described sampling. Each trajectory, initiated at milestone i was examined during the calculation to determine the termination time at milestone i+1 or i−1. The termination times were binned to generate the LFPTD (or Kij(t)—the probability that a trajectory initiated at milestone i will terminate at milestone j after time t). In Fig. 2 we show K76(t) and K78(t).
Figure 2
The local first passage time distribution (LFPTD) of milestone 7. The times of trajectories initiated at milestone 7 and terminated at the previous or the next milestone are binned. The times from milestone 7 to 6 are black and from 7 to 8 are red. Moving forward (from 7 to 8) is significantly less likely, suggesting a free energy barrier (only nine trajectories made it to milestone 8).Once all the LFPTD (or Kij(t)) are computed, we calculate the time evolution of the system. In Faradjian and Elber 6 and West, Elber, and Shalloway 7, an integral equation was introduced for Pi(t), the probability of finding the system at milestone i at time t. We also introduced Qi(t), the probability of making a transition to milestone i exactly at time t. The outer integrals below are from zero to t and the inner integral from zero to t−t′,
 |
The first equation provides the probability of transition to i exactly at time t. The right-hand side is a sum of the probabilities of past transitions to milestones i±1 and a transition at the right time into i. The second equation is for the probability of being at i at time t. The right-hand side is the probability of making a transition to i some time in the past and remaining there until t. These integral equations can be solved directly with the initial condition Pi(0) and setting the last milestone to be absorbing. They can also be solved in terms of moments as has been shown in Rousseau et al. 13. Of special interest is the first moment, the so-called mean first passage time 〈〈τ〉〉 that is widely used to assess timescales of stochastic processes: 〈〈τ〉〉=ʃτ · dPf/dτ · dτ, where f denotes the last absorbing boundary. The first passage time was estimated as 10±9μs. Although the error bars are large, the results are within experimental values 14,15.
With an agreement with experiment on the overall timescale at hand, the atomically detailed simulations allow us to investigate mechanistic questions. For example, are the phenylalanine flips the cause of the free energy barrier, or are other processes contributing significantly to the free energy profile? In Fig. 3 we show the distribution of the distances between the two phenylalanine rings that are useful indicators of their R (or T) states.
Figure 3
The distributions of the Cz–Cz distances of the two Phenylalanines Phe-97 for milestones 3 (black), 5 (red), 7 (green), and 9 (blue). This distance is a useful measure for the R (or T) transition. Most of the side-chain transitions are performed in the neighborhood of milestone 5. The distribution of milestone 7 is similar to milestone 9, suggesting that the side-chain transitions are complete at milestone 7.Fig. 3 demonstrates the interplay between different steps of the calculations. The calculation of the reaction coordinate (or the determination of the positions of the milestones) suggests that the final rearrangement of the protein structure are of a global nature and not associated with the side-chain transition since the side-chain transition is essentially complete in milestone 7 (see also Fig. 4). This observation was confirmed by the thermal sampling in the milestone plane. But how relevant is the last phase, which does not include the phenylalanines? The first passage calculations at milestone 7 indicate that the free energy barrier for the last phase is significant and impacts the kinetics (in addition to an earlier barrier associated with the rotation of the side chains). Structural differences between milestone 7 and the deoxy state are shown in Fig. 4.
Figure 4
The difference in the packing of Phe-97 (white) against the heme (red) in milestone 7 (top) and the deoxy state (bottom). While the phenylalanine packing is essentially complete, there remain numerous adjustments. The lysine and the water molecules should come higher and Phe-97 requires a slight rotation. The carboxyl groups of the heme also change orientation.What are the factors that allowed milestoning to complete a calculation not accessible to straightforward molecular dynamics? There are three factors:
1. The existence of milestones allows bootstrapping of rare trajectories. If we place a milestone halfway to the top of the barrier, then climbing to the top will be easier when initiating trajectories half the way.
2. Since the trajectories are diffusive and keep moving forward and backward, placing milestones reduces the number of backward and forward cycles we need to compute explicitly with the trajectories (the cycles are considered in the integral equations).
3. The large number of short trajectories is easy to distribute on a cluster of computers, a common facility in computationally oriented laboratories.
2. Since the trajectories are diffusive and keep moving forward and backward, placing milestones reduces the number of backward and forward cycles we need to compute explicitly with the trajectories (the cycles are considered in the integral equations).
3. The large number of short trajectories is easy to distribute on a cluster of computers, a common facility in computationally oriented laboratories.
In conclusion, we illustrated the capacity of milestoning to study the kinetics in a complex system and point out a late activated phase in the allosteric transition of HbI, a step of nonlocal structural rearrangements.
Acknowledgments
I thank Anthony West for his help in setting up the calculations and to Quentin H. Gibson for many useful discussions.
This research was supported by National Institutes of Health grant No. GM05979605.
References and foodnotes
1. (1998). The stereochemical mechanism of the cooperative effects in hemoglobin revisited. Annu. Rev. Biophys. Biomol. Struct. 27, 1–34. CrossRef | PubMed
2. (1984). Kinetics of ligand binding and quaternary conformational changes in the homodimeric hemoglobin from Scapharca inaequivalvis. J. Biol. Chem. 259, 6730–6738. PubMed
3. (2006). Allosteric transition in real time: Time-resolved crystallographic studies of cooperative dimeric hemoglobin. Proc. Natl. Acad. Sci. USA 103, 7649–7654. CrossRef | PubMed
4. (1935). The activated complex in chemical reactions. J. Chem. Phys. 3, 107–115. CrossRef | PubMed
5. (2002). Transition path sampling. Adv. Chem. Phys. 123, 1–78. PubMed
6. (2004). Computing timescales from reaction coordinates by milestoning. J. Chem. Phys. 120, 10880–10889. CrossRef | PubMed
7. West, A. M. A., R. Elber, and D. Shalloway. Extending molecular dynamics timescales with milestoning: example of complex kinetics in a solvated peptide. J. Chem. Phys. In press..
8. (1994). Crystal structure of oxygenated Scapharca dimeric hemoglobin at 1.7Ångström resolution. J. Biol. Chem. 269, 25259–25267. PubMed
9. (1994). High-resolution crystallographic analysis of co-operative dimeric hemoglobin. J. Mol. Biol. 235, 657–681. CrossRef | PubMed
10. (2000). The Protein Data Bank. Nucleic Acids Res. 28, 235–242. CrossRef | PubMed
11. (1990). Self avoiding walk between two fixed points as a tool to calculate reaction paths in large molecular system. Int. J. Quantum Chem. 24, 167–186. PubMed
12. (1995). MOIL: a program for simulations of macromolecules. Comput. Phys. Comm. 91, 159–189. PubMed
13. (1993). Heme-Hem interactions in homodimeric cooperative hemoglobin—evidence from transient Raman scattering. J. Biol. Chem. 268, 5719–5723. PubMed
14. Nichols, J. C., W. E. Royer, Jr., and Q. H. Gibson. An optical signal correlated with the allosteric transition in Scapharca inaequivalvis HbI. Biochemistry. 45:15748–15755..
15. (2006). Efficient computation of the first passage time distribution of the generalized master equation by steady-state relaxation. J. Chem. Phys. 124, , Art. No. 054112. PubMed
Publication Information
Received: November 25, 2006
Accepted: February 20, 2007
