help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Biophys J. BioFAST on February 16, 2007.
doi:10.1529/biophysj.106.098442
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.098442v1
92/10/3442    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Anselmi, M.
Right arrow Articles by Amadei, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Anselmi, M.
Right arrow Articles by Amadei, A.
Biophysical Journal 92:3442-3447 (2007)
© 2007 The Biophysical Society

Theoretical Characterization of Carbon Monoxide Vibrational Spectrum in Sperm Whale Myoglobin Distal Pocket

Massimiliano Anselmi *, Massimiliano Aschi {dagger}, Alfredo Di Nola * and Andrea Amadei {ddagger}

* Department of Chemistry, University of Rome "La Sapienza", Rome, Italy; {dagger} Department of Chemistry, Chemical Engineering and Materials, University of L'Aquila, L'Aquila, Italy; and {ddagger} Department of Chemistry, University of Rome "Tor Vergata", Rome, Italy

Correspondence: Address reprint requests to Dr. Andrea Amadei, Department of Chemistry, University of Rome "Tor Vergata", Via della Ricerca Scientifica 1, 00133 Rome, Italy. Tel.: 39-06-7259-4905; Fax: 39-06-7259-4328; E-mail: andrea.amadei{at}uniroma2.it.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 METHODS
 RESULTS
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
In this article we use the perturbed matrix method and an extended molecular dynamics sampling of the carbon monoxide (CO) in the myoglobin distal pocket to characterize the CO vibrational spectrum and hence to relate its spectroscopic features with the atomic-molecular behavior. Results show the accuracy of the method employed and confirm the assignment of the spectroscopic B1 and B2 states proposed by Lim et al.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
METHODS
 RESULTS
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Myoglobin (Mb), one of the most studied proteins so far, has been characterized by a great variety of techniques ranging from x-ray crystallography to molecular simulation (1Go–21Go). Myoglobin's small size, relative structural stability, and yet complex functional behavior (involving relevant conformational fluctuations as well as covalent binding of ligands) make this protein a virtually perfect model system to investigate, at the atomistic level, protein biochemical activity. Carbon monoxide (CO), one of the possible myoglobin ligands, has been used as a probe for the kinetics of the covalent binding process (2Go,22Go) (heme-CO binding) as well as for characterizing the Mb-ligand interaction during ligand diffusion in the protein matrix (23Go,24Go). In particular, CO (time-resolved) vibrational spectra have been extensively used in the last decade to determine CO behavior just after photolysis of its covalent bond to heme (14Go,19Go,25Go,26Go), characterizing the presence of two distinct spectroscopic states (B1, B2) probably corresponding to two different CO-heme relative orientations and still not uniquely assigned (27Go–33Go). Recent attempts to obtain vibrational (IR) spectra of CO in myoglobin by molecular dynamics (MD) simulations (34Go,35Go), providing a reasonable qualitative reproduction of experimental data, were based on a purely classical CO model, treating its intramolecular (quantum) vibrational mode as a classical degree of freedom and hence evaluating the IR spectrum via the CO dipole autocorrelation function. Such an approach is heavily dependent on the details of the semiempirical CO model used and may provide artifacts caused by the classical approximation of vibrational excitation.

Recently (36Go,37Go) we extended the perturbed matrix method (PMM), introduced previously (38Go,39Go), to model quantum-mechanically vibrational excitations of molecules in complex systems (i.e., including the effects of the atomic-molecular environment on IR spectra). In those articles we applied our method to CO in solution (water and chloroform), showing its efficiency and reliability. In this article we apply the same methodology to determine the CO IR spectrum in the Mb distal pocket, as obtained by PMM and MD simulations using the three-site "quadrupolar" CO model (40Go) employed in previous articles.


   METHODS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
RESULTS
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Initial coordinates for the simulation of Mb with photodissociated CO were taken from the 1.15 Å resolution crystal structure of CO-bound sperm whale Mb (PDB entry 1BZR) (41Go), in which we cut the CO-Fe bond. Thus, from the beginning of the simulation the system was modeled as an unliganded state of Mb with CO. According to recently published articles (33Go,42Go), His(E7)64 was modeled as the neutral tautomer, with hydrogen at the {varepsilon} position. The protein was solvated in a box with explicit SPC water molecules (43Go), large enough to contain the protein and 0.8 nm of solvent on all sides. The total number of atoms for the system was ~21,000.

MD simulations were performed with the Gromacs software package (44Go) using the GROMOS96 force field (45Go). The CO molecule was modeled with the three-site "quadrupolar" CO model (40Go). Simulations were carried out at constant temperature (300 K) using the isothermal temperature coupling (46Go) within a fixed-volume rectangular box and using periodic boundary conditions. The Lincs algorithm (47Go), to constrain bond lengths, and the rototranslational constraint algorithm (48Go), to stop protein rototranslational motions, were used. The initial velocities were taken randomly from a Maxwellian distribution at 300 K, and a time step of 2 fs was used in all simulations.

The particle mesh Ewald (PME) method (49Go) was used for the calculation of the long-range interactions with a grid spacing of 0.12 nm combined with a fourth-order B-spline interpolation to compute the potential and forces between grid points. A nonbond pair list cutoff of 9.0 Å was used for short-range interactions, and the pair list was updated every five time steps.

After thermalization and equilibration, we perform an initial 1-ns simulation to obtain a relaxed structure of the system with the ligand in the principal docking site of the distal pocket. Starting from this structure, we have performed five MD simulations using different initial velocities, given by 300 K Maxwellian distributions, each stopped at the first carbon monoxide escape from the distal pocket for a total of 21 ns of MD simulation of the photodissociated CO in the distal pocket.

Carbon monoxide vibrational states were obtained using the method described in detail in a previous article (36Go). In brief, the ground-state electronic energy along the internuclear distance of an isolated CO molecule was obtained by density functional theory (DFT) calculations. Note that, as described in a previous article (36Go), the internuclear distance range we utilized is about ± 0.15 Å around the equilibrium distance (~1.1 Å), which is, for a stiff vibrational mode as in CO, a proper fluctuation range for estimating the vibrational frequency within the harmonic approximation (36Go). Becke's three parameters exchange (50Go) and Lee, Yang, Parr correlation (51Go) function (B3LYP) were used for DFT calculations in conjunction with triple-{zeta} atomic basis set with polarization and diffuse functions, i.e., aug-cc-pv-tz (52Go). Configuration interactions (53Go) including single, double, and triple excitations (CISDT) calculations, using the above B3LYP/ aug-cc-pv-tz orbitals, were then carried out at each internuclear distance using an active space as large as 10 electrons in 35 orbitals for evaluating the unperturbed electronic states considered for PMM calculations. As shown in the previous article (36Go), such a computational procedure provides a very accurate description of vibrational and electronic states of the isolated CO molecule. All our quantum chemical calculations on isolated carbon monoxide were performed using the Gamess US package (54Go).

The essence of PMM is to use high-quality unperturbed electronic states as a basis set to express the Hamiltonian matrix of the quantum center (CO molecule) including the electric field perturbation resulting from the atomic environment (36Go–39Go), which is approximately equivalent to CI calculations including the perturbing electric field in the Hamiltonian operator.

Therefore, at each MD frame we obtained, by means of PMM, the corresponding perturbed electronic states providing the corresponding perturbed energy and dipole moment along the intramolecular coordinate (internuclear distance), hence allowing the evaluation of perturbed CO harmonic vibrational states ({phi}v) by solving (36Go)

Formula

In the previous equations, the vibrational Hamiltonian operator Hv is defined by the reduced mass µ', the intramolecular coordinate ß, and the harmonic force constant k obtained via quadratic fit of the perturbed electronic ground state energy in ß. Once the perturbed vibrational eigenstates and eigenvalues ({phi}v, Uv) were evaluated along the MD trajectory, we easily obtained the vibrational spectrum I({lambda}) (considering a unitary radiation density per unit frequency) via (36Go)

Formula
where B is the Einstein coefficient for the first perturbed vibrational excitation and {rho}({lambda}) is the corresponding probability density of excitation in the frequency-wavelength space. Note that I({lambda}), as expressed by the last equation, is not equivalent to the frequency probability density (in our case {rho}({lambda})) typically reported in other articles, as it involves the transition dipole effect.


   RESULTS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
To evaluate the equilibrium IR spectrum of CO within the distal pocket, we used the MD configurations as obtained by five independent MD trajectories all starting from the principal docking site (the most probable CO site just after photolysis) and interrupted at the first CO escape from distal pocket. Interestingly, no CO escaped outside the protein within the MD sampling achieved and the total 21 ns corresponding to CO within the distal pocket as obtained by the five MD trajectories, provide a distal pocket escape mean life of ~4–5 ns. Note that the configuration storing frequency (a configuration every 1 ps) guaranteed no time correlation and a good convergence for the excitation properties provided by PMM (55Go).

In Fig. 1 we show the IR spectrum and also report the noise (one standard deviation of the signal) for each bin used along the frequency axes (the bin size, defining our spectrum resolution, is 1 cm–1). Note that the calculated IR spectrum is ~60 cm–1 red-shifted with respect to the experimental one as a result of the ~60 cm–1 red shift provided by quantum chemical calculations used for the isolated CO (36Go) and in line with the accuracy limit of sophisticated quantum chemical calculations in determining the vibrational frequency. The figure clearly indicates the presence of two well-separated peaks, whose signal difference is far beyond the noise, corresponding to the experimentally observed B1 and B2 peaks, although with a lower frequency separation: ~2 cm–1 (our calculations) versus ~10 cm–1 (experimental) (24Go). Although it reproduces the experimental spectrum shape and width reasonably closely (within the noise) (32Go), the theoretical IR spectrum in Fig. 1 underestimates the peak shift as well as the absorption full range (~30 cm–1 versus ~60 cm–1). It is worth noting that for a molecule like CO, the IR spectrum broadening we compute can be ascribed only to the perturbing field fluctuations as provided by the environment atomic motions.


Figure 1
View larger version (9K):
[in this window]
[in a new window]

  		FIGURE 1  Carbon monoxide vibrational spectrum in myoglobin distal pocket as obtained by PMM and MD simulations. The error bars shown correspond to a standard deviation of the property, and time is expressed in atomic units.

 
Interestingly, evaluations of the vibrational frequency distribution in the myoglobin distal pocket, as obtained by completely classical models (10Go,34Go), provided an absorption full range between ~20 cm–1 and ~60 cm–1, indeed showing that for such classical methods the model details and/or the actual strategy employed to evaluate the vibrational frequencies (estimating the classical perturbed stretching constant or via the time autocorrelation function of the classically fluctuating dipole) may cause significant variations. The same methods also show discrepancies concerning the shape of the frequency distribution: one (10Go), as the authors correctly state, does not provide two clear peaks; the other (34Go) provides two main peaks (shifted by ~8 cm–1), but other relevant peaks are also present, thus raising serious doubts about the quantitative reliability of such evaluations. In fact, in both articles the frequency distribution shown is rather noisy, as expected by the limited sampling used, and no error bars are reported, making it impossible to judge the significance of the various peaks and hence very difficult to do a quantitative comparison with our results. Moreover, our quantum mechanically based calculations showed that the perturbed transition dipole for the vibrational excitation considered is not constant in the absorption frequency range, and the use of the frequency distribution alone provides a larger peak separation (4–5 cm–1), implying that disregarding the transition dipole is not really appropriate to describe spectroscopic features.

In recent literature on condensed phase spectroscopy (56Go) the relevance of the dynamic correlation of the excitation frequency has been pointed out. Analysis of time autocorrelation function of the excitation frequency as obtained from our MD simulations and shown in Fig. 2 provides a mean correlation time of ~200 fs, well matching a similar evaluation performed on the electronic excitation of solvated acetone (55Go). Such a short correlation time results from the fast relaxation of the perturbing electric field associated to the motions of the environment atoms. To quantitatively evaluate the B1 and B2 interconversion rates, we calculated the transition time distribution for crossing the ~2079 cm–1 border frequency (see Fig. 1) corresponding to the IR spectrum minimum between the B1 and B2 peaks. Fig. 3 shows, on a logarithmic scale, such distributions indicating that both B1->B2 and B2->B1 transitions may be well described by an exponential decay, with time constants of ~2 ps, well matching the experimental observations (32Go). These results demonstrate the PMM accuracy, the good quality of the GROMOS and quadrupolar three sites, CO force fields, and the importance of using a quantum mechanically based method for evaluating the excitation spectra.


Figure 2
View larger version (8K):
[in this window]
[in a new window]

  		FIGURE 2  Time autocorrelation function of the excitation frequency as provided by PMM and MD simulations.

 

Figure 3
View larger version (7K):
[in this window]
[in a new window]

  		FIGURE 3  Transition time distributions for the B1 and B2 interconversion (see text).

 
The B1 and B2 spectroscopic states have been extensively studied in recent years (23Go,24Go,27Go,29Go–31Go,33Go), leading to the widely accepted idea that such states correspond to the opposite orientations, in the principal docking site, of the CO molecule with respect to the iron atom. However, the assignment of the B1 and B2 states to CO orientations is still controversial. In fact, Lim et al. (28Go,29Go,31Go) proposed the B2 (low-frequency) state to be defined by the CO orientation with the carbon atom pointing toward the iron, whereas Nienhaus et al. (30Go,33Go) made the opposite assignment, and theoretical computational attempts, largely based on classical approximations in modeling vibrational excitation, did not provide conclusive quantitative results, although they showed qualitative agreement (10Go) with the Lim et al. assignment. To address this problem in detail, we first analyzed the CO behavior in the distal pocket in terms of its orientations with respect to the heme plane as defined by the corresponding CO polar angles (the out-of-plane angle {theta} and the in plane rotation angle {phi}) (34Go).

In Fig. 4 we show the distribution of the MD configurations on the {theta}-{phi} plane, clearly indicating the presence of two stable angular conformations both centered at {theta} {approx} 90° (i.e., CO parallel to the heme plane). These two angular states correspond, in the principal docking site, to the two opposite CO orientations toward the iron, with the state centered at {theta} {approx} 90°, {phi} {approx} –60° associated to the CO orientation with the carbon atom pointing toward the iron. Interestingly, within the 21 ns considered, the CO molecule mostly resided in the principal docking site as expected from previous computational and experimental data (5Go,6Go,9Go,10Go,13Go,57Go,58Go).


Figure 4
View larger version (40K):
[in this window]
[in a new window]

  		FIGURE 4  Projection of the MD sampling over polar angles plane defining the carbon monoxide orientation with respect to the heme plane. Note that the out-of-plane angle {theta} is zero when the CO molecule is perpendicular to the heme plane.

 
In Fig. 5 we show the IR mean excitation frequency as a function of the rotation angle {phi} along with the corresponding noise. It is evident that the largest frequency shift is for the two CO rotational orientations corresponding to {phi} {approx} –60° and {phi} {approx} 120°, indicating that the B1 and B2 peaks found in our spectrum are mainly caused by such rotational CO states as indeed clearly shown by Fig. 6, where we report the difference of the probability distributions in {phi} for the low- and high-frequency IR peaks (such distributions are obtained by the subpopulations belonging to the bins corresponding to the maxima, and the use of the distributions difference filters out the noise present in both subpopulations). These two CO rotational states defined by the {phi} angle do not exactly correspond to the CO orientations with the carbon or oxygen toward the iron but rather to the opposite CO dipole orientations in the heme plane, equivalent to the Lim et al. assignment when the CO molecule is located in the principal docking site. Such results clearly show that at physiological conditions the Lim et al. B1 and B2 states assignment is likely to be the correct one, although the B states should be properly described in terms of CO dipole orientations in the heme plane rather than orientations with either carbon or oxygen pointing toward the iron. Finally, in Fig. 7 we show the difference between the mean electric field component parallel to the carbon monoxide bond caused by each protein residue, the heme, and solvent, as provided by the two subpopulations corresponding to the B1 and B2 absorption maxima. It can be noted that the infrared absorption split between the B states, mainly a result of the CO dipole rotation with respect to the (perturbing) electric field, is largely determined by some key residues (in particular residues in the distal pocket), the heme group, and the solvent. Interestingly, the latter generates the largest electric field variation, hence showing the solvent relevance in the CO spectroscopic behavior in myoglobin.


Figure 5
View larger version (10K):
[in this window]
[in a new window]

  		FIGURE 5  Mean excitation frequency as a function of the rotational angle {phi}. The error bars shown correspond to a standard deviation of the property.

 

Figure 6
View larger version (15K):
[in this window]
[in a new window]

  		FIGURE 6  Difference of the probability distributions as a function of the polar angle {phi}, as obtained subtracting the high-frequency peak probability from that of the low-frequency peak (see text).

 

Figure 7
View larger version (9K):
[in this window]
[in a new window]

  		FIGURE 7  Difference between the mean electric field component parallel to the carbon monoxide bond attributable to each protein residue, the Heme, and solvent, as obtained by the two subpopulations corresponding to the B1 and B2 absorption maxima.

 

   CONCLUSIONS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 CONCLUSIONS
ACKNOWLEDGEMENTS
 REFERENCES
 
The combined use of PMM with an extended MD sampling of the CO in the myoglobin distal pocket provided a clear assignment of the experimentally observed B states to the opposite CO rotational orientations in the heme plane, fully confirming the Lim et al. assignment based on spectroscopic data and in line with previous computational studies (10Go,34Go,35Go) that, although employing a purely classical approach, reproduced qualitatively the main experimental observations. Interestingly, such rotational states correspond to the opposite CO dipole orientations with respect to the perturbing electric field, hence suggesting the possibility of similar IR splitting in other myoglobin cavities.

Moreover, it emerged that such B-state splitting is largely determined by specific interactions including the CO-solvent one, which exerts the largest contribution.


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 CONCLUSIONS
 ACKNOWLEDGEMENTS
REFERENCES
 
We gratefully acknowledge "Centro Ricerche Studi Enrico Fermi" (Rome, Italy) and Consorzio interuniversitario per le Applicazioni di Supercalcolo per Università e Ricerca (Rome, Italy) for computational support.

Submitted on October 2, 2006; accepted for publication January 18, 2007.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Kendrew, J. C., and M. F. Perutz. 1957. X-ray studies of compounds of biological interest. Annu. Rev. Biochem. 26:327–372.[CrossRef][Medline]

2. Austin, R. H., K. W. Beeson, L. Eisenstein, H. Frauenfelder, and I. C. Gunsalus. 1975. Dynamics of ligand binding to myoglobin. Biochemistry. 14:5355–5373.[CrossRef][Medline]

3. Tilton, R. F., Jr., I. D. Kuntz, Jr., and G. A. Petsko. 1984. Cavities in proteins: structure of a metmyoglobin-xenon complex solved to 1.9 A. Biochemistry. 23:2849–2857.[CrossRef][Medline]

4. Elber, R., and M. Karplus. 1990. Enhanced sampling in molecular dynamics: use of the time-dependent Hartree approximation for a simulation of carbon monoxide diffusion through myoglobin. J. Am. Chem. Soc. 112:9161–9175.[CrossRef]

5. Schlichting, I., J. Berendzen, G. N. Phillips, Jr., and R. M. Sweet. 1994. Crystal structure of photolysed carbonmonoxy-myoglobin. Nature. 371:808–812.[CrossRef][Medline]

6. Teng, T. Y., V. Srajer, and K. Moffat. 1994. Photolysis-induced structural changes in single crystals of carbonmonoxy myoglobin at 40 K. Nat. Struct. Biol. 1:701–705.[CrossRef][Medline]

7. Srajer, V., T. Teng, T. Ursby, C. Pradervand, Z. Ren, S. Adachi, W. Schildkamp, D. Bourgeois, M. Wulff, and K. Moffat. 1996. Photolysis of the carbon monoxide complex of myoglobin: nanosecond time-resolved crystallography. Science. 274:1726–1729.[Abstract/Free Full Text]

8. Scott, E. E., and Q. H. Gibson. 1997. Ligand migration in sperm whale myoglobin. Biochemistry. 36:11909–11917.[CrossRef][Medline]

9. Vitkup, D., G. A. Petsko, and M. Karplus. 1997. A comparison between molecular dynamics and x-ray results for dissociated CO in myoglobin. Nat. Struct. Biol. 4:202–208.[CrossRef][Medline]

10. Meller, J., and R. Elber. 1998. Computer simulations of carbon monoxide photodissociation in myoglobin: structural interpretation of the B states. Biophys. J. 74:789–802.[Abstract/Free Full Text]

11. Schlichting, I., and K. Chu. 2000. Trapping intermediates in the crystal: ligand binding to myoglobin. Curr. Opin. Struct. Biol. 10:744–752.[CrossRef][Medline]

12. Scott, E. E., Q. H. Gibson, and J. S. Olson. 2001. Mapping the pathways for O2 entry into and exit from myoglobin. J. Biol. Chem. 276:5177–5188.[Abstract/Free Full Text]

13. Srajer, V., Z. Ren, T. Y. Teng, M. Schmidt, T. Ursby, D. Bourgeois, C. Pradervand, W. Schildkamp, M. Wulff, and K. Moffat. 2001. Protein conformational relaxation and ligand migration in myoglobin: a nanosecond to millisecond molecular movie from time-resolved Laue x-ray diffraction. Biochemistry. 40:13802–13815.[CrossRef][Medline]

14. Lamb, D. C., K. Nienhaus, A. Arcovito, F. Draghi, A. E. Miele, M. Brunori, and G. U. Nienhaus. 2002. Structural dynamics of myoglobin: ligand migration among protein cavities studied by Fourier transform infrared/temperature derivative spectroscopy. J. Biol. Chem. 277:11636–11644.[Abstract/Free Full Text]

15. Bourgeois, D., B. Vallone, F. Schotte, A. Arcovito, A. E. Miele, G. Sciara, M. Wulff, P. Anfinrud, and M. Brunori. 2003. Complex landscape of protein structural dynamics unveiled by nanosecond Laue crystallography. Proc. Natl. Acad. Sci. USA. 100:8704–8709.[Abstract/Free Full Text]

16. Schotte, F., M. Lim, T. A. Jackson, A. V. Smirnov, J. Soman, J. S. Olson, G. N. Phillips, Jr., M. Wulff, and P. A. Anfinrud. 2003. Watching a protein as it functions with 150-ps time-resolved x-ray crystallography. Science. 300:1944–1947.[Abstract/Free Full Text]

17. Aschi, M., C. Zazza, R. Spezia, C. Bossa, A. Di Nola, M. Paci, and A. Amadei. 2004. Conformational fluctuations and electronic properties in myoglobin. J. Comput. Chem. 25:974–984.[CrossRef][Medline]

18. Bossa, C., M. Anselmi, D. Roccatano, A. Amadei, B. Vallone, M. Brunori, and A. Di Nola. 2004. Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin. Biophys. J. 86:3855–3862.[Abstract/Free Full Text]

19. Lamb, D. C., A. Arcovito, K. Nienhaus, O. Minkow, F. Draghi, M. Brunori, and G. U. Nienhaus. 2004. Structural dynamics of myoglobin: an infrared kinetic study of ligand migration in mutants YQR and YQRF. Biophys. Chem. 109:41–58.[CrossRef][Medline]

20. Bossa, C., A. Amadei, I. Daidone, M. Anselmi, B. Vallone, M. Brunori, and A. Di Nola. 2005. Molecular dynamics simulation of sperm whale myoglobin: effects of mutations and trapped CO on the structure and dynamics of cavities. Biophys. J. 89:465–474.[Abstract/Free Full Text]

21. Bourgeois, D., B. Vallone, A. Arcovito, G. Sciara, F. Schotte, P. A. Anfinrud, and M. Brunori. 2006. Extended subnanosecond structural dynamics of myoglobin revealed by Laue crystallography. Proc. Natl. Acad. Sci. USA. 103:4924–4929.[Abstract/Free Full Text]

22. Angeloni, L., and A. Feis. 2003. Protein relaxation in the photodissociation of myoglobin-CO complexes. Photochem. Photobiol. Sci. 2:730–740.[CrossRef][Medline]

23. Nienhaus, K., P. Deng, J. M. Kriegl, and G. U. Nienhaus. 2003. Structural dynamics of myoglobin: spectroscopic and structural characterization of ligand docking sites in myoglobin mutant L29W. Biochemistry. 42:9633–9646.[CrossRef][Medline]

24. Nienhaus, K., P. Deng, J. M. Kriegl, and G. U. Nienhaus. 2003. Structural dynamics of myoglobin: effect of internal cavities on ligand migration and binding. Biochemistry. 42:9647–9658.[CrossRef][Medline]

25. Li, T., M. L. Quillin, G. N. Phillips, Jr., and J. S. Olson. 1994. Structural determinants of the stretching frequency of CO bound to myoglobin. Biochemistry. 33:1433–1446.[CrossRef][Medline]

26. Nienhaus, G. U., and K. Nienhaus. 2002. Infrared study of carbon monoxide migration among internal cavities of myoglobin mutant L29W. J. Biol. Phys. 28:163–172.[CrossRef]

27. Lim, M., T. A. Jackson, and P. A. Anfinrud. 1995. Binding of CO to myoglobin from a heme pocket docking site to form nearly linear Fe-C-O. Science. 269:962–966.[Abstract/Free Full Text]

28. Lim, M., T. A. Jackson, and P. A. Anfinrud. 1995. Mid-infrared vibrational spectrum of CO after photodissociation from heme: evidence for a ligand docking site in the heme pocket of hemoglobin and myoglobin. J. Chem. Phys. 102:4355–4366.[CrossRef]

29. Lim, M., T. A. Jackson, and P. A. Anfinrud. 1997. Ultrafast rotation and trapping of carbon monoxide dissociated from myoglobin. Nat. Struct. Biol. 4:209–214.[CrossRef][Medline]

30. Nienhaus, K., P. Deng, J. S. Olson, J. J. Warren, and G. U. Nienhaus. 2003. Structural dynamics of myoglobin: ligand migration and binding in valine 68 mutants. J. Biol. Chem. 278:42532–42544.[Abstract/Free Full Text]

31. Lim, M., T. A. Jackson, and P. A. Anfinrud. 2004. Orientational distribution of CO before and after photolysis of MbCO and HbCO: a determination using time-resolved polarized Mid-IR spectroscopy. J. Am. Chem. Soc. 126:7946–7957.[CrossRef][Medline]

32. Kim, S., and M. Lim. 2005. Picosecond dynamics of ligand interconversion in the primary docking site of heme proteins. J. Am. Chem. Soc. 127:5786–5787.[CrossRef][Medline]

33. Nienhaus, K., J. S. Olson, S. Franzen, and G. U. Nienhaus. 2005. The origin of stark splitting in the initial photoproduct state of MbCO. J. Am. Chem. Soc. 127:40–41.[CrossRef][Medline]

34. Nutt, D. R., and M. Meuwly. 2003. Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin. Biophys. J. 85:3612–3623.[Abstract/Free Full Text]

35. Nutt, D. R., and M. Meuwly. 2004. CO migration in native and mutant myoglobin: atomistic simulations for the understanding of protein function. Proc. Natl. Acad. Sci. USA. 101:5998–6002.[Abstract/Free Full Text]

36. Amadei, A., F. Marinelli, M. D'Abramo, M. D'Alessandro, M. Anselmi, A. Di Nola, and M. Aschi. 2005. Theoretical modeling of vibroelectronic quantum states in complex molecular systems: solvated carbon monoxide, a test case. J. Chem. Phys. 122: article 124506.

37. D'Alessandro, M., F. Marinelli, M. D'Abramo, M. Aschi, A. Di Nola, and A. Amadei. 2005. Ground and excited electronic state thermodynamics of aqueous carbon monoxide: a theoretical study. J. Chem. Phys. 122: article 124507.

38. Aschi, M., R. Spezia, A. Di Nola, and A. Amadei. 2001. A first-principles method to model perturbed electronic wavefunctions: the effect of an external homogeneous electric field. Chem. Phys. Lett. 344:374–380.[CrossRef]

39. Spezia, R., M. Aschi, A. Di Nola, and A. Amadei. 2002. Extension of the perturbed matrix method: application to a water molecule. Chem. Phys. Lett. 365:450–456.[CrossRef]

40. Straub, J. E., and M. Karplus. 1991. Molecular dynamics study of the photodissociation of carbon monoxide from myoglobin: ligand dynamics in the first 10 ps. Chem. Phys. 158:221–248.[CrossRef]

41. Kachalova, G. S., A. N. Popov, and H. D. Bartunik. 1999. A steric mechanism for inhibition of CO binding to heme proteins. Science. 284:473–476.[Abstract/Free Full Text]

42. De Angelis, F., A. A. Jarzecki, R. Car, and T. G. Spiro. 2005. Quantum chemical evaluation of protein control over heme ligation: CO/O2 discrimination in myoglobin. J. Phys. Chem. B. 109:3065–3070.[Medline]

43. Berendsen, H. J. C., J. P. M. Postma, W. F. van Gunsteren, and J. Hermans. 1981. Interaction models for water in relation to protein hydration. In Intermolecular Forces. Pullman, B., editor. D. Reidel Publishing Company, Dordrecht, The Netherlands. 331–342.

44. Berendsen, H. J. C., D. van der Spoel, and R. van Drunen. 1995. GROMACS: A message-passing parallel molecular dynamics implementation. Comp. Phys. Comm. 91:43–56.[CrossRef]

45. van Gunsteren, W. F., S. Billeter, A. Eising, P. Hunenberger, P. Kruger, A. E. Mark, W. Scott, and I. Tironi. 1996. Biomolecular Simulations: The GROMOS96 Manual and User Guide. BIOMOS bv, Zurich, Groningen.

46. Evans, D. J., and G. P. Morriss. 1990. Statistical Mechanics of Nonequilibrium Liquids. Academic Press, London.

47. Hess, B., H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije. 1997. LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 18:1463–1472.[CrossRef]

48. Amadei, A., G. Chillemi, M. A. Ceruso, A. Grottesi, and A. Di Nola. 2000. Molecular dynamics simulations with constrained roto-translational motions: Theoretical basis and statistical mechanical consistency. J. Chem. Phys. 112:9–23.[CrossRef]

49. Essmann, U., L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen. 1995. A smooth particle mesh Ewald method. J. Chem. Phys. 103:8577–8593.[CrossRef]

50. Becke, A. D. 1993. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98:5648–5652.[CrossRef]

51. Lee, C., W. Yang, and R. G. Parr. 1988. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B. 37:785–789.[CrossRef]

52. Dunning, T. H., Jr. 1989. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 90:1007–1023.[CrossRef]

53. Brooks, B. R., and H. F. Schaefer III. 1979. The graphical unitary group approach to the electron correlation problem. Methods and preliminary applications. J. Chem. Phys. 70:5092–5106.[CrossRef]

54. Schmidt, M. W., K. K. Baldridge, J. A. Boatz, S. T. Elbert, M. S. Gordon, J. H. Jensen, S. Koseki, N. Matsunaga, K. A. Nguyen, S. Su, T. L. Windus, M. Dupuis, and J. A. Montgomery, Jr. 1993. General atomic and molecular electronic structure system. J. Comput. Chem. 14:1347–1363.[CrossRef]

55. D'Abramo, M., M. Aschi, A. Di Nola, and A. Amadei. 2006. On the importance of configurational sampling in theoretical calculation of electronic properties of complex molecular systems: Acetone in water. Chem. Phys. Lett. 424:289–294.[CrossRef]

56. Lawrence, C. P., and J. L. Skinner. 2005. Quantum corrections in vibrational and electronic condensed phase spectroscopy: line shapes and echoes. Proc. Natl. Acad. Sci. USA. 102:6720–6725.[Abstract/Free Full Text]

57. Hartmann, H., S. Zinser, P. Komninos, R. T. Schneider, G. U. Nienhaus, and F. Parak. 1996. X-ray structure determination of a metastable state of carbonmonoxy myoglobin after photodissociation. Proc. Natl. Acad. Sci. USA. 93:7013–7016.[Abstract/Free Full Text]

58. Ma, J., S. Huo, and J. E. Straub. 1997. Molecolar dynamics simulation study of the B-states of solvated carbon monoxymyoglobin. J. Am. Chem. Soc. 119:2541–2551.[CrossRef]




This article has been cited by other articles:


Home page
 Biophys. JHome page
M. Anselmi, A. Di Nola, and A. Amadei
The Kinetics of Ligand Migration in Crystallized Myoglobin as Revealed by Molecular Dynamics Simulations
Biophys. J., June 1, 2008; 94(11): 4277 - 4281.
[Abstract] [Full Text] [PDF]

Home page
 Biophys. JHome page
N. Plattner and M. Meuwly
The Role of Higher CO-Multipole Moments in Understanding the Dynamics of Photodissociated Carbonmonoxide in Myoglobin
Biophys. J., April 1, 2008; 94(7): 2505 - 2515.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.098442v1
92/10/3442    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Anselmi, M.
Right arrow Articles by Amadei, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Anselmi, M.
Right arrow Articles by Amadei, A.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2007 by the Biophysical Society.