Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Instability of surface-temperature filaments in strain and shear

The effects of uniform straining and shearing on the stability of a surface quasi-geostrophic temperature filament are investigated. Straining is shown to stabilize perturbations for wide filaments but only for a finite time until the filament thins to a critical width, after which some perturbations can grow. No filament can be stabilized in practice, since there are perturbations that can grow large for any strain rate. The optimally growing perturbations, defined as solutions that reach a certain threshold amplitude first, are found numerically for a wide range of parameter values. The radii of the vortices formed through nonlinear roll-up are found to be proportional to θ/s, where θ is the temperature anomaly of the filament and s the strain rate, and are not dependent on the initial size of the filament. Shearing is shown to reduce the normal-mode growth rates, but it cannot stabilize them completely when there are temperature discontinuities in the basic state; smooth filaments can be stabilized completely by shearing and a simple scaling argument provides the shear rate required. Copyright © 2010 Royal Meteorological Society