Extracting Aggregation Free Energies of Mixed Clusters from Simulations of Small Systems: Application to Ionic Surfactant Micelles.

Micelle cluster distributions from molecular dynamics simulations of a solvent-free coarse-grained model of sodium octyl sulfate (SOS) were analyzed using an improved method to extract equilibrium association constants from small-system simulations containing one or two micelle clusters at equilibrium with free surfactants and counterions. The statistical-thermodynamic and mathematical foundations of this partition-enabled analysis of cluster histograms (PEACH) approach are presented. A dramatic reduction in computational time for analysis was achieved through a strategy similar to the selector variable method to circumvent the need for exhaustive enumeration of the possible partitions of surfactants and counterions into clusters. Using statistics from a set of small-system (up to 60 SOS molecules) simulations as input, equilibrium association constants for micelle clusters were obtained as a function of both number of surfactants and number of associated counterions through a global fitting procedure. The resulting free energies were able to accurately predict micelle size and charge distributions in a large (560 molecule) system. The evolution of micelle size and charge with SOS concentration as predicted by the PEACH-derived free energies and by a phenomenological four-parameter model fit, along with the sensitivity of these predictions to variations in cluster definitions, are analyzed and discussed.

X-ray spectroscopy of nitrile hydratase at pH 7 and 9.

The iron K-edge X-ray absorption spectrum of Rhodococcus sp. R312 (formerly Brevibacterium sp. R312) nitrile hydratase in frozen solutions at pH 7 and 9 has been analyzed to determine details of the iron coordination. EXAFS analysis implies two or three sulfur ligands per iron and overall six coordination; together with previous EPR and ENDOR results, this implies an N3S2O ligation sphere. The bond lengths from EXAFS analysis [rav(Fe-S) = 2.21 A at pH 7.3; rav(Fe-N/O) = 1.99 A] support cis coordination of two cysteine ligands and conclusively rule out nitric oxide coordination to the iron, a possibility proposed on the basis of an FTIR difference experiment [Noguchi, T., Honda, J., Nagamune, T., Sasabe, H., Inoue, Y., & Endo, I. (1995) FEBS Lett. 358, 9-12]. The higher-frequency EXAFS can be simulated well by inclusion of multiple scattering from two or three imidazole ligands; the fit to the data is improved if first-sphere multiple scattering pathways are also included. A slight shortening (by 0.02 +/- 0.01 A) of one or both Fe-S bonds when the pH is raised from 7.3 to 9.0 is consistent with shifts observed in the Raman spectrum [Brennan et al. (1996) Biochemistry 35, 10068-10077].

Covalent alteration of the prosthetic heme of human hemoglobin by BrCCl3. Cross-linking of heme to cysteine residue 93.

Recent studies have shown that a protein-bound heme adduct formed from the reaction of BrCCl3 with myoglobin was due to bonding of the proximal histidine residue through the ring I vinyl of a heme-CCl2 moiety. The present study reveals that BrCCl3 also reacts with the heme of reduced human hemoglobin to form two protein-bound heme adducts. Edman degradation and mass spectrometry provided evidence that these protein-bound heme adducts were addition products in which heme-CCL2 or heme-CCl3 were bound to cysteine residue 93 of the beta-chain of hemoglobin. It appeared that the cysteine residue was bonded regiospecifically to the ring I vinyl group of the altered heme moiety, because the nonprotein-bound products of the reaction included the beta-carboxyvinyl and alpha-hydroxy-beta-trichloromethylethyl derivatives of the ring I vinyl moiety of heme. The absorption spectra of the protein-bound adducts in both the oxidized and reduced states were highly similar to those described for hemichromes, which are thought to be involved in the formation of Heinz bodies and subsequent red cell lysis.