Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Peroxo-iron and oxenoid-iron species as alternative oxygenating agents in cytochrome P450-catalyzed reactions: switching by threonine-302 to alanine mutagenesis of cytochrome P450 2B4.

Among biological catalysts, cytochrome P450 is unmatched in its multiplicity of isoforms, inducers, substrates, and types of chemical reactions catalyzed. In the present study, evidence is given that this versatility extends to the nature of the active oxidant. Although mechanistic evidence from several laboratories points to a hypervalent iron-oxenoid species in P450-catalyzed oxygenation reactions, Akhtar and colleagues [Akhtar, M., Calder, M. R., Corina, D. L. & Wright, J. N. (1982) Biochem. J. 201, 569-580] proposed that in steroid deformylation effected by P450 aromatase an iron-peroxo species is involved. We have shown more recently that purified liver microsomal P450 cytochromes, including phenobarbital-induced P450 2B4, catalyze the analogous deformylation of a series of xenobiotic aldehydes with olefin formation. The investigation presented here on the effect of site-directed mutagenesis of threonine-302 to alanine on the activities of recombinant P450 2B4 with N-terminal amino acids 2-27 deleted [2B4 (delta2-27)] makes use of evidence from other laboratories that the corresponding mutation in bacterial P450s interferes with the activation of dioxygen to the oxenoid species by blocking proton delivery to the active site. The rates of NADPH oxidation, hydrogen peroxide production, and product formation from four substrates, including formaldehyde from benzphetamine N-demethylation, acetophenone from 1-phenylethanol oxidation, cyclohexanol from cyclohexane hydroxylation, and cyclohexene from cyclohexane carboxaldehyde deformylation, were determined with P450s 2B4, 2B4 (delta2-27), and 2B4 (delta2-27) T302A. Replacement of the threonine residue in the truncated cytochrome gave a 1.6- to 2.5-fold increase in peroxide formation in the presence of a substrate, but resulted in decreased product formation from benzphetamine (9-fold), cyclohexane (4-fold), and 1-phenylethanol (2-fold). In sharp contrast, the deformylation of cyclohexane carboxaldehyde by the T302A mutant was increased about 10-fold. On the basis of these findings and our previous evidence that aldehyde deformylation is supported by added H202, but not by artificial oxidants, we conclude that the iron-peroxy species is the direct oxygen donor. It remains to be established which of the many other oxidative reactions involving P450 utilize this species and the extent to which peroxo-iron and oxenoid-iron function as alternative oxygenating agents with the numerous isoforms of this versatile catalyst.

BIND—The Biomolecular Interaction Network Database

The Biomolecular Interaction Network Database (BIND; http://binddb.org) is a database designed to store full descriptions of interactions, molecular complexes and pathways. Development of the BIND 2.0 data model has led to the incorporation of virtually all components of molecular mechanisms including interactions between any two molecules composed of proteins, nucleic acids and small molecules. Chemical reactions, photochemical activation and conformational changes can also be described. Everything from small molecule biochemistry to signal transduction is abstracted in such a way that graph theory methods may be applied for data mining. The database can be used to study networks of interactions, to map pathways across taxonomic branches and to generate information for kinetic simulations. BIND anticipates the coming large influx of interaction information from high-throughput proteomics efforts including detailed information about post-translational modifications from mass spectrometry. Version 2.0 of the BIND data model is discussed as well as implementation, content and the open nature of the BIND project. The BIND data specification is available as ASN.1 and XML DTD.

A role for intermolecular disulfide bonds in prion diseases?

The key event in prion diseases seems to be the conversion of the prion protein PrP from its normal cellular isoform (PrPC) to an aberrant “scrapie” isoform (PrPSc). Earlier studies have detected no covalent modification in the scrapie isoform and have concluded that the PrPC → PrPSc conversion is a purely conformational transition involving no chemical reactions. However, a reexamination of the available biochemical data suggests that the PrPC → PrPSc conversion also involves a covalent reaction of the (sole) intramolecular disulfide bond of PrPC. Specifically, the data are consistent with the hypothesis that infectious prions are composed of PrPSc polymers linked by intermolecular disulfide bonds. Thus, the PrPC → PrPSc conversion may involve not only a conformational transition but also a thiol/disulfide exchange reaction between the terminal thiolate of such a PrPSc polymer and the disulfide bond of a PrPC monomer. This hypothesis seems to account for several unusual features of prion diseases.

Femtochemistry

The topic of femtochemistry is surveyed from both theoretical and experimental points of view. A time-dependent wave packet description of the photodissociation of the O—C—S molecule reveals vibrational motion in the transition-state region and suggests targets for direct experimental observation. Theoretical approaches for treating femtosecond chemical phenomena in condensed phases are featured along with prospects for laser-controlled chemical reactions by using tailored ultrashort chirped pulses. An experimental study of the photoisomerization of retinal in the protein bacteriorhodopsin is discussed with an aim to gain insight into the potential energy surfaces on which this remarkably efficient and selective reactions proceeds. Finally, a prospective view of new frontiers in femtochemistry is given.

Manganese oxide minerals: Crystal structures and economic and environmental significance

Manganese oxide minerals have been used for thousands of years—by the ancients for pigments and to clarify glass, and today as ores of Mn metal, catalysts, and battery material. More than 30 Mn oxide minerals occur in a wide variety of geological settings. They are major components of Mn nodules that pave huge areas of the ocean floor and bottoms of many fresh-water lakes. Mn oxide minerals are ubiquitous in soils and sediments and participate in a variety of chemical reactions that affect groundwater and bulk soil composition. Their typical occurrence as fine-grained mixtures makes it difficult to study their atomic structures and crystal chemistries. In recent years, however, investigations using transmission electron microscopy and powder x-ray and neutron diffraction methods have provided important new insights into the structures and properties of these materials. The crystal structures for todorokite and birnessite, two of the more common Mn oxide minerals in terrestrial deposits and ocean nodules, were determined by using powder x-ray diffraction data and the Rietveld refinement method. Because of the large tunnels in todorokite and related structures there is considerable interest in the use of these materials and synthetic analogues as catalysts and cation exchange agents. Birnessite-group minerals have layer structures and readily undergo oxidation reduction and cation-exchange reactions and play a major role in controlling groundwater chemistry.

Putting ion channels to work: Mechanoelectrical transduction, adaptation, and amplification by hair cells

As in other excitable cells, the ion channels of sensory receptors produce electrical signals that constitute the cellular response to stimulation. In photoreceptors, olfactory neurons, and some gustatory receptors, these channels essentially report the results of antecedent events in a cascade of chemical reactions. The mechanoelectrical transduction channels of hair cells, by contrast, are coupled directly to the stimulus. As a consequence, the mechanical properties of these channels shape our hearing process from the outset of transduction. Channel gating introduces nonlinearities prominent enough to be measured and even heard. Channels provide a feedback signal that controls the transducer's adaptation to large stimuli. Finally, transduction channels participate in an amplificatory process that sensitizes and sharpens hearing.

Analysis of metabolic pathways by the growth of cells in the presence of organic solvents.

A new approach to the analysis of metabolic pathways involving poorly water-soluble intermediates is proposed. It relies upon the ability of the hydrophobic intermediates formed by a sequence of intracellular reactions to cross the membrane(s) and partition between aqueous and organic phases, when cells are incubated in the presence of a nonpolar and nontoxic organic solvent. As a result of this thermodynamically driven efflux of the formed intermediates from the cell, they accumulate in the organic medium in sufficient quantities for GC-MS analysis and identification. This enables direct determination of the sequence of chemical reactions involved with no requirement for the isolation of each individual metabolite from a cell-free extract. The feasibility of the proposed methodology has been demonstrated by the elucidation of the biosynthesis of (R)-gamma-decalactone from (R)-ricinoleic acid catalyzed by the yeast Sporidiobolus ruinenii grown in the presence of decane. The corresponding 4-hydroxy-acid intermediates, formed in the course of beta-oxidation of (R)-ricinoleic acid, were simultaneously observed in a single experiment on the same chromatogram. Potential applications of this proposed methodology are briefly discussed.

Negative activation enthalpies in the kinetics of protein folding.

Although the rates of chemical reactions become faster with increasing temperature, the converse may be observed with protein-folding reactions. The rate constant for folding initially increases with temperature, goes through a maximum, and then decreases. The activation enthalpy is thus highly temperature dependent because of a large change in specific heat (delta Cp). Such a delta Cp term is usually presumed to be a consequence of a large decrease in exposure of hydrophobic surfaces to water as the reaction proceeds from the denatured state to the transition state for folding: the hydrophobic side chains are surrounded by "icebergs" of water that melt with increasing temperature, thus making a large contribution to the Cp of the denatured state and a smaller one to the more compact transition state. The rate could also be affected by temperature-induced changes in the conformational population of the ground state: the heat required for the progressive melting of residual structure in the denatured state will contribute to delta Cp. By examining two proteins with different refolding mechanisms, we are able to find both of these two processes; barley chymotrypsin inhibitor 2, which refolds from a highly unfolded state, fits well to a hydrophobic interaction model with a constant delta Cp of activation, whereas barnase, which refolds from a more structured denatured state, deviates from this ideal behavior.

Catalytically critical nucleotide in domain 5 of a group II intron.

Domain 5 (D5) is a small hairpin structure within group II introns. A bimolecular assay system depends on binding by D5 to an intron substrate for self-splicing activity. In this study, mutations in D5 identify two among six nearly invariant nucleotides as being critical for 5' splice junction hydrolysis but unimportant for binding. A mutation at another site in D5 blocks binding. Thus, mutations can distinguish two D5 functions: substrate binding and catalysis. The secondary structure of D5 may resemble helix I formed by the U2 and U6 small nuclear RNAs in the eukaryotic spliceosome. Our results support a revision of the previously proposed correspondence between D5 and helix I on the basis of the critical trinucleotide 5'-AGC-3' present in both. We suggest that this trinucleotide plays a similar role in promoting the chemical reactions for both splicing systems.

Mechanical and chemical unfolding of a single protein: A comparison

Is the mechanical unraveling of protein domains by atomic force microscopy (AFM) just a technological feat or a true measurement of their unfolding? By engineering a protein made of tandem repeats of identical Ig modules, we were able to get explicit AFM data on the unfolding rate of a single protein domain that can be accurately extrapolated to zero force. We compare this with chemical unfolding rates for untethered modules extrapolated to 0 M denaturant. The unfolding rates obtained by the two methods are the same. Furthermore, the transition state for unfolding appears at the same position on the folding pathway when assessed by either method. These results indicate that mechanical unfolding of a single protein by AFM does indeed reflect the same event that is observed in traditional unfolding experiments. The way is now open for the extensive use of AFM to measure folding reactions at the single-molecule level. Single-molecule AFM recordings have the added advantage that they define the reaction coordinate and expose rare unfolding events that cannot be observed in the absence of chemical denaturants.

Chemical camouflage of antigenic determinants: Stealth erythrocytes

In a number of clinical circumstances it would be desirable to artificially conceal cellular antigenic determinants to permit survival of heterologous donor cells. A case in point is the problem encountered in transfusions of patients with rare blood types or chronically transfused patients who become allosensitized to minor blood group determinants. We have tested the possibility that chemical modification of the red blood cell (RBC) membrane might serve to occlude antigenic determinants, thereby minimizing transfusion reactions. To this end, we have covalently bound methoxy(polyethylene glycol) (mPEG) to the surface of mammalian RBC via cyanuric chloride coupling. Human RBC treated with this technique lose ABO blood group reactivity as assessed by solution–phase antisera agglutination. In accord with this, we also find a profound decrease in anti-blood group antibody binding. Furthermore, whereas human monocytes avidly phagocytose untreated sheep RBC, mPEG-derivatized sheep RBC are ineffectively phagocytosed. Surprisingly, human and mouse RBC appear unaffected by this covalent modification of the cell membrane. Thus, mPEG-treated RBC are morphologically normal, have normal osmotic fragility, and mPEG-derivatized murine RBC have normal in vivo survival, even following repeated infusions. Finally, in preliminary experiments, mPEG-modified sheep RBC intraperitoneally transfused into mice show significantly improved (up to 360-fold) survival when compared with untreated sheep RBC. We speculate that similar chemical camouflage of intact cells may have significant clinical applications in both transfusion (e.g., allosensitization and autoimmune hemolytic disease) and transplantation (e.g., endothelial cells and pancreatic β cells) medicine.

Chemical synthesis of the precursor molecule of the Aequorea green fluorescent protein, subsequent folding, and development of fluorescence

The present paper describes the total chemical synthesis of the precursor molecule of the Aequorea green fluorescent protein (GFP). The molecule is made up of 238 amino acid residues in a single polypeptide chain and is nonfluorescent. To carry out the synthesis, a procedure, first described in 1981 for the synthesis of complex peptides, was used. The procedure is based on performing segment condensation reactions in solution while providing maximum protection to the segment. The effectiveness of the procedure has been demonstrated by the synthesis of various biologically active peptides and small proteins, such as human angiogenin, a 123-residue protein analogue of ribonuclease A, human midkine, a 121-residue protein, and pleiotrophin, a 136-residue protein analogue of midkine. The GFP precursor molecule was synthesized from 26 fully protected segments in solution, and the final 238-residue peptide was treated with anhydrous hydrogen fluoride to obtain the precursor molecule of GFP containing two Cys(acetamidomethyl) residues. After removal of the acetamidomethyl groups, the product was dissolved in 0.1 M Tris⋅HCl buffer (pH 8.0) in the presence of DTT. After several hours at room temperature, the solution began to emit a green fluorescence (λmax = 509 nm) under near-UV light. Both fluorescence excitation and fluorescence emission spectra were measured and were found to have the same shape and maxima as those reported for native GFP. The present results demonstrate the utility of the segment condensation procedure in synthesizing large protein molecules such as GFP. The result also provides evidence that the formation of the chromophore in GFP is not dependent on any external cofactor.

Electrogenic light reactions in photosystem I: resolution of electron-transfer rates between the iron-sulfur centers.

Flash-induced voltage changes (electrogenic events) in photosystem I particles from spinach, oriented in a phospholipid layer, have been studied at room temperature on a time scale ranging from 1 micros to several seconds. A phospholipid layer containing photosystem I particles was adsorbed to a Teflon film separating two aqueous compartments. Voltage changes were measured across electrodes immersed in the compartments. In the absence of added electron donors and acceptors, a multiphasic voltage increase, associated with charge separation, was followed by a decrease, associated with charge recombination. Several kinetic phases were resolved: a rapid (<1 micros) increase, ascribed to electron transfer from the primary electron donor P700 to the iron-sulfur electron acceptor FB, was followed by a slower, biphasic increase with time constants of 30 and 200 micros. The 30-micros phase is assigned to electron transfer from FB to the iron-sulfur center FA. The voltage decrease had a time constant of 90 ms, ascribed to charge recombination from FA to P700. Upon chemical prereduction of FA and FB the 30- and 200-micros phases disappeared and the decay time constant was accelerated to 330 micros, assigned to charge recombination from the phylloquinone electron acceptor (A1) or the iron-sulfur center FX to P700.