Symmetry in chemistry from the hydrogen atom to proteins

The last 2 decades have seen discoveries in highly excited states of atoms and molecules of phenomena that are qualitatively different from the “planetary” model of the atom, and the near-rigid model of molecules, characteristic of these systems in their low-energy states. A unified view is emerging in terms of approximate dynamical symmetry principles. Highly excited states of two-electron atoms display “molecular” behavior of a nonrigid linear structure undergoing collective rotation and vibration. Highly excited states of molecules described in the “standard molecular model” display normal mode couplings, which induce bifurcations on the route to molecular chaos. New approaches such as rigid–nonrigid correlation, vibrons, and quantum groups suggest a unified view of collective electronic motion in atoms and nuclear motion in molecules.

On the electronic nature of low-barrier hydrogen bonds in enzymatic reactions

The electronic nature of low-barrier hydrogen bonds (LBHBs) in enzymatic reactions is discussed based on combined low temperature neutron and x-ray diffraction experiments and on high level ab initio calculations by using the model substrate benzoylacetone. This molecule has a LBHB, as the intramolecular hydrogen bond is described by a double-well potential with a small barrier for hydrogen transfer. From an “atoms in molecules” analysis of the electron density, it is found that the hydrogen atom is stabilized by covalent bonds to both oxygens. Large atomic partial charges on the hydrogen-bonded atoms are found experimentally and theoretically. Therefore, the hydrogen bond gains stabilization from both covalency and from the normal electrostatic interactions found for long, weak hydrogen bonds. Based on comparisons with other systems having short-strong hydrogen bonds or LBHBs, it is proposed that all short-strong and LBHB systems possess similar electronic features of the hydrogen-bonded region, namely polar covalent bonds between the hydrogen atom and both heteroatoms in question.

Comparison of structures and energies of CH52+• with CH4+• and their possible role in superacidic methane activation

Contrary to previous theoretical studies at the UHF/6-31G* level, the methonium radical dication CH52+ is not a Cs symmetrical structure with a 2e—3c bond but a C2v symmetrical structure 1 with two 2e—3c bonds (at the UHF/6-31G**, UMP2/6-31G**, and UQCISD(T)/6-311G** levels). The Cs symmetrical structure is not even a minimum at the higher level of calculations. The four hydrogen atoms in 1 are bonded to the carbon atom by two 2e—3c bonds and the fifth hydrogen atom by a 2e—2c bond. The unpaired electron of 1 is located in a formal p-orbital (of the sp2-hybridized carbon atom) perpendicular to the plane of the molecule. Hydrogen scrambling in 1 is however extremely facile, as is in other C1 cations. It is found that the protonation of methane to CH5+ decreases the energy for subsequent homolytic cleavage resulting in the exothermic (24.1 kcal/mol) formation of CH4+•. Subsequent reaction with neutral methane while reforming CH5+ gives the methyl radical enabling reaction with excess methane to ethane and H2. The overall reaction is endothermic by 11.4 kcal/mol, but offers under conditions of oxidative removal of H2 an alternative to the more energetic carbocationic conversion of methane.

Designing safer chemicals: predicting the rates of metabolism of halogenated alkanes.

A computational model is presented that can be used as a tool in the design of safer chemicals. This model predicts the rate of hydrogen-atom abstraction by cytochrome P450 enzymes. Excellent correlations between biotransformation rates and the calculated activation energies (delta Hact) of the cytochrome P450-mediated hydrogen-atom abstractions were obtained for the in vitro biotransformation of six halogenated alkanes (1-fluoro-1,1,2,2-tetrachloroethane, 1,1-difluoro-1,2,2-trichloroethane, 1,1,1-trifluro-2,2-dichloroethane, 1,1,1,2-tetrafluoro-2-chloroethane, 1,1,1,2,2,-pentafluoroethane, and 2-bromo-2-chloro-1,1,1-trifluoroethane) with both rat and human enzyme preparations: In(rate, rat liver microsomes) = 44.99 - 1.79(delta Hact), r2 = 0.86; In(rate, human CYP2E1) = 46.99 - 1.77(delta Hact), r2 = 0.97 (rates are in nmol of product per min per nmol of cytochrome P450 and energies are in kcal/mol). Correlations were also obtained for five inhalation anesthetics (enflurane, sevoflurane, desflurane, methoxyflurane, and isoflurane) for both in vivo and in vitro metabolism by humans: In[F(-)]peak plasma = 42.87 - 1.57(delta Hact), r2 = 0.86. To our knowledge, these are the first in vivo human metabolic rates to be quantitatively predicted. Furthermore, this is one of the first examples where computational predictions and in vivo and in vitro data have been shown to agree in any species. The model presented herein provides an archetype for the methodology that may be used in the future design of safer chemicals, particularly hydrochlorofluorocarbons and inhalation anesthetics.

Structure at 2.7 Å resolution of the Paracoccus denitrificans two-subunit cytochrome c oxidase complexed with an antibody FV fragment

The aa3 type cytochrome c oxidase consisting of the core subunits I and II only was isolated from the soil bacterium Paracoccus denitrificans and crystallized as complex with a monoclonal antibody Fv fragment. Crystals could be grown in the presence of a number of different nonionic detergents. However, only undecyl-β-d-maltoside and cyclohexyl-hexyl-β-d-maltoside yielded well-ordered crystals suitable for high resolution x-ray crystallographic studies. The crystals belong to space group P212121 and diffract x-rays to at least 2.5 Å (1 Å = 0.1 nm) resolution using synchrotron radiation. The structure was determined to a resolution of 2.7 Å using molecular replacement and refined to a crystallographic R-factor of 20.5% (Rfree = 25.9%). The refined model includes subunits I and II and the 2 chains of the Fv fragment, 2 heme A molecules, 3 copper atoms, and 1 Mg/Mn atom, a new metal (Ca) binding site, 52 tentatively identified water molecules, and 9 detergent molecules. Only four of the water molecules are located in the cytoplasmic half of cytochrome c oxidase. Most of them are near the interface of subunits I and II. Several waters form a hydrogen-bonded cluster, including the heme propionates and the Mg/Mn binding site. The Fv fragment binds to the periplasmic polar domain of subunit II and is critically involved in the formation of the crystal lattice. The crystallization procedure is well reproducible and will allow for the analysis of the structures of mechanistically interesting mutant cytochrome c oxidases.

The change in hydrogen bond strength accompanying charge rearrangement: Implications for enzymatic catalysis

The equilibrium for formation of the intramolecular hydrogen bond (KHB) in a series of substituted salicylate monoanions was investigated as a function of ΔpKa, the difference between the pKa values of the hydrogen bond donor and acceptor, in both water and dimethyl sulfoxide. The dependence of log KHB upon ΔpKa is linear in both solvents, but is steeper in dimethyl sulfoxide (slope = 0.73) than in water (slope = 0.05). Thus, hydrogen bond strength can undergo substantially larger increases in nonaqueous media than aqueous solutions as the charge density on the donor or acceptor atom increases. These results support a general mechanism for enzymatic catalysis, in which hydrogen bonding to a substrate is strengthened as charge rearranges in going from the ground state to the transition state; the strengthening of the hydrogen bond would be greater in a nonaqueous enzymatic active site than in water, thus providing a rate enhancement for an enzymatic reaction relative to the solution reaction. We suggest that binding energy of an enzyme is used to fix the substrate in the low-dielectric active site, where the strengthening of the hydrogen bond in the course of a reaction is increased.

Physical reasons for the unusual α-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH2 in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

Polar residues drive association of polyleucine transmembrane helices

Although many polar residues are directly involved in transmembrane protein functions, the extent to which they contribute to more general structural features is still unclear. Previous studies have demonstrated that asparagine residues can drive transmembrane helix association through interhelical hydrogen bonding [Choma, C., Gratkowski, H., Lear, J. D. & DeGrado, W. F. (2000) Nat. Struct. Biol. 7, 161–166; and Zhou, F. X., Cocco, M. J., Russ, W. P., Brunger, A. T. & Engelman, D. M. (2000) Nat. Struct. Biol. 7, 154–160]. We have studied the ability of other polar residues to promote helix association in detergent micelles and in biological membranes. Our results show that polyleucine sequences with Asn, Asp, Gln, Glu, and His, residues capable of being simultaneously hydrogen bond donors and acceptors, form homo- or heterooligomers. In contrast, polyleucine sequences with Ser, Thr, and Tyr do not associate more than the polyleucine sequence alone. The results therefore provide experimental evidence that interactions between polar residues in the helices of transmembrane proteins may serve to provide structural stability and oligomerization specificity. Furthermore, such interactions can allow structural flexibility required for the function of some membrane proteins.

Barriers to protein folding: formation of buried polar interactions is a slow step in acquisition of structure.

In the MYL mutant of the Arc repressor dimer, sets of partially buried salt-bridge and hydrogen-bond interactions mediated by Arg-31, Glu-36, and Arg-40 in each subunit are replaced by hydrophobic interactions between Met-31, Tyr-36, and Leu-40. The MYL refolding/dimerization reaction differs from that of wild type in being 10- to 1250-fold faster, having an earlier transition state, and depending upon viscosity but not ionic strength. Formation of the wild-type salt bridges in a hydrophobic environment clearly imposes a kinetic barrier to folding, which can be lowered by high salt concentrations. The changes in the position of the transition state and viscosity dependence can be explained if denatured monomers interact to form a partially folded dimeric intermediate, which then continues folding to form the native dimer. The second step is postulated to be rate limiting for wild type. Replacing the salt bridge with hydrophobic interactions lowers this barrier for MYL. This makes the first kinetic barrier rate limiting for MYL refolding and creates a downhill free-energy landscape in which most molecules which reach the intermediate state continue to form native dimers.