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.

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.

The resource consumption principle: attention and memory in volumes of neural tissue.

In the cerebral cortex, the small volume of the extracellular space in relation to the volume enclosed by synapses suggests an important functional role for this relationship. It is well known that there are atoms and molecules in the extracellular space that are absolutely necessary for synapses to function (e.g., calcium). I propose here the hypothesis that the rapid shift of these atoms and molecules from extracellular to intrasynaptic compartments represents the consumption of a shared, limited resource available to local volumes of neural tissue. Such consumption results in a dramatic competition among synapses for resources necessary for their function. In this paper, I explore a theory in which this resource consumption plays a critical role in the way local volumes of neural tissue operate. On short time scales, this principle of resource consumption permits a tissue volume to choose those synapses that function in a particular context and thereby helps to integrate the many neural signals that impinge on a tissue volume at any given moment. On longer time scales, the same principle aids in the stable storage and recall of information. The theory provides one framework for understanding how cerebral cortical tissue volumes integrate, attend to, store, and recall information. In this account, the capacity of neural tissue to attend to stimuli is intimately tied to the way tissue volumes are organized at fine spatial scales.

The invariant system of coordinates of antibody molecules: prediction of the "standard" C alpha framework of VL and VH domains.

A new approach of comparing protein structures that does not involve the procedure of superposition is suggested. An invariant system of coordinates for immunoglobulin molecules that is based on the geometrical symmetry inherent to the variable domain light-chain (VL)-heavy-chain (VH) complex is described. The coordinates of the Calpha atoms in 22 immunoglobulin structures are calculated in the invariant system of coordinates. We found that 76 identical positions in this Calpha framework are symmetrical about the twofold axis. Comparison of the identical positions in these molecules allows us to select 96 positions in the light chains and 87 positions in the heavy chains whose Calpha atom coordinates are approximately the same. To check whether the average coordinates of Calpha atoms in these positions complies with the stereochemical requirements, we calculated Calpha-Calpha distances. Seventy-three positions of the light chains and 72 positions of the heavy chains satisfy the Calpha-Calpha distance criterion. The Calpha atoms in these positions are used for constructing the "standard" Calpha framework of VL and VH complexes. The average coordinates of Calpha atoms are presented.

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.

Inhibitory sites in enzymes: Zinc removal and reactivation by thionein

Thionein (T) has not been isolated previously from biological material. However, it is generated transiently in situ by removal of zinc from metallothionein under oxidoreductive conditions, particularly in the presence of selenium compounds. T very rapidly activates a group of enzymes in which zinc is bound at an inhibitory site. The reaction is selective, as is apparent from the fact that T does not remove zinc from the catalytic sites of zinc metalloenzymes. T instantaneously reverses the zinc inhibition with a stoichiometry commensurate with its known capacity to bind seven zinc atoms in the form of clusters in metallothionein. The zinc inhibition is much more pronounced than was previously reported, with dissociation constants in the low nanomolar range. Thus, T is an effective, endogenous chelating agent, suggesting the existence of a hitherto unknown and unrecognized biological regulatory system. T removes the metal from an inhibitory zinc-specific enzymatic site with a resultant marked increase of activity. The potential significance of this system is supported by the demonstration of its operations in enzymes involved in glycolysis and signal transduction.

Dynamic expression of multiple scavenger receptor cysteine-rich genes in coelomocytes of the purple sea urchin

Coelomocytes, the heterogeneous population of sea urchin putative immune cells, were found to express a complex set of transcripts featuring scavenger receptor cysteine-rich (SRCR) repeats. SRCR domains define a metazoan superfamily of proteins, many of which are implicated in development and regulation of the immune system of vertebrates. Coelomocytes transcribe multiple SRCR genes from among a multigene family encoding an estimated number of 1,200 SRCR domains in specific patterns particular to each individual. Transcription levels for given SRCR genes may range from pronounced to undetectable, yet all tested animals harbor the genomic loci encoding these genes. Analysis of several SRCR genes revealed multiple loci corresponding to each type. In the case of one SRCR type, a cluster of at least three genes was detected within a 133-kb bacterial artificial chromosome insert, and conserved as well as unique regions were identified in sequences of three genomic clones derived from a single animal. Array hybridizations with repeated samples of coelomocyte messages revealed substantial alterations in levels of expression of many SRCR genes, with fluctuations of up to 10-fold in 1 week and up to 30-fold over a period of 3 months. This report is the first demonstration of genomic and transcriptional complexity in molecules expressed by invertebrate coelomocytes. The mechanisms controlling SRCR gene expression and the functional significance of this dynamic system await elucidation.

The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Packing at the protein-water interface.

We have determined the packing efficiency at the protein-water interface by calculating the volumes of atoms on the protein surface and nearby water molecules in 22 crystal structures. We find that an atom on the protein surface occupies, on average, a volume approximately 7% larger than an atom of equivalent chemical type in the protein core. In these calculations, larger volumes result from voids between atoms and thus imply a looser or less efficient packing. We further find that the volumes of individual atoms are not related to their chemical type but rather to their structural location. More exposed atoms have larger volumes. Moreover, the packing around atoms in locally concave, grooved regions of protein surfaces is looser than that around atoms in locally convex, ridge regions. This as a direct manifestation of surface curvature-dependent hydration. The net volume increase for atoms on the protein surface is compensated by volume decreases in water molecules near the surface. These waters occupy volumes smaller than those in the bulk solvent by up to 20%; the precise amount of this decrease is directly related to the extent of contact with the protein.

Imaging of single molecule diffusion.

In recent years observations at the level of individual atoms and molecules became possible by microscopy and spectroscopy. Imaging of single fluorescence molecules has been achieved but has so far been restricted to molecules in the immobile state. Here we provide methodology for visualization of the motion of individual fluorescent molecules. It is applied to imaging of the diffusional path of single molecules in a phospholipid membrane by using phospholipids carrying one rhodamine dye molecule. For this methodology, fluorescence microscopy was carried to a sensitivity so that single fluorescent molecules illuminated for only 5 ms were resolvable at a signal/noise ratio of 28. Repeated illuminations permitted direct observation of the diffusional motion of individual molecules with a positional accuracy of 30 nm. Such capability has fascinating potentials in bioscience--for example, to correlate biological functions of cell membranes with movements, spatial organization, and stoichiometries of individual components.

Proteins on ribosome surface: Measurements of protein exposure by hot tritium bombardment technique

The hot tritium bombardment technique [Goldanskii, V. I., Kashirin, I. A., Shishkov, A. V., Baratova, L. A. & Grebenshchikov, N. I. (1988) J. Mol. Biol. 201, 567–574] has been applied to measure the exposure of proteins on the ribosomal surface. The technique is based on replacement of hydrogen by high energy tritium atoms in thin surface layer of macromolecules. Quantitation of tritium radioactivity of each protein has revealed that proteins S1, S4, S5, S7, S18, S20, and S21 of the small subunit, and proteins L7/L12, L9, L10, L11, L16, L17, L24, and L27 of the large subunit are well exposed on the surface of the Escherichia coli 70 S ribosome. Proteins S8, S10, S12, S16, S17, L14, L20, L29, L30, L31, L32, L33, and L34 have virtually no groups exposed on the ribosomal surface. The remaining proteins are found to be exposed to lesser degree than the well exposed ones. No additional ribosomal proteins was exposed upon dissociation of ribosomes into subunits, thus indicating the absence of proteins on intersubunit contacting surfaces.

Optical trapping and manipulation of neutral particles using lasers

The techniques of optical trapping and manipulation of neutral particles by lasers provide unique means to control the dynamics of small particles. These new experimental methods have played a revolutionary role in areas of the physical and biological sciences. This paper reviews the early developments in the field leading to the demonstration of cooling and trapping of neutral atoms in atomic physics and to the first use of optical tweezers traps in biology. Some further major achievements of these rapidly developing methods also are considered.

Crystal structure of a Staphylococcus aureus protein A domain complexed with the Fab fragment of a human IgM antibody: Structural basis for recognition of B-cell receptors and superantigen activity

Staphylococcus aureus produces a virulence factor, protein A (SpA), that contains five homologous Ig-binding domains. The interactions of SpA with the Fab region of membrane-anchored Igs can stimulate a large fraction of B cells, contributing to lymphocyte clonal selection. To understand the molecular basis for this activity, we have solved the crystal structure of the complex between domain D of SpA and the Fab fragment of a human IgM antibody to 2.7-Å resolution. In the complex, helices II and III of domain D interact with the variable region of the Fab heavy chain (VH) through framework residues, without the involvement of the hypervariable regions implicated in antigen recognition. The contact residues are highly conserved in human VH3 antibodies but not in other families. The contact residues from domain D also are conserved among all SpA Ig-binding domains, suggesting that each could bind in a similar manner. Features of this interaction parallel those reported for staphylococcal enterotoxins that are superantigens for many T cells. The structural homology between Ig VH regions and the T-cell receptor Vβ regions facilitates their comparison, and both types of interactions involve lymphocyte receptor surface remote from the antigen binding site. However, T-cell superantigens reportedly interact through hydrogen bonds with T-cell receptor Vβ backbone atoms in a primary sequence-independent manner, whereas SpA relies on a sequence-restricted conformational binding with residue side chains, suggesting that this common bacterial pathogen has adopted distinct molecular recognition strategies for affecting large sets of B and T lymphocytes.

An information theory model of hydrophobic interactions.

A molecular model of poorly understood hydrophobic effects is heuristically developed using the methods of information theory. Because primitive hydrophobic effects can be tied to the probability of observing a molecular-sized cavity in the solvent, the probability distribution of the number of solvent centers in a cavity volume is modeled on the basis of the two moments available from the density and radial distribution of oxygen atoms in liquid water. The modeled distribution then yields the probability that no solvent centers are found in the cavity volume. This model is shown to account quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes. The connection of information theory to statistical thermodynamics provides a basis for clarification of hydrophobic effects. The simplicity and flexibility of the approach suggest that it should permit applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.

X-ray absorption fine structure as a monitor of zinc coordination sites during oogenesis of Xenopus laevis.

The x-ray absorption fine structure (XAFS) zinc K-edge steps for intact stages I,II and V,VI Xenopus laevis oocytes demonstrate that the zinc concentration is about 3 and 1 mM, respectively. However, the chi(k) function for the early stage oocytes differs markedly from that for the late one. Analysis of the XAFS data for stage I,II oocytes indicates that zinc is bound to 2.0 +/- 0.5 sulfur atoms at an average coordination distance of 2.29 +/- 0.02 angstroms and 2.0 +/- 0.5 nitrogen or oxygen (N/O) atoms at 2.02 +/- 0.02 angstroms. In marked contrast, in stage V,VI oocytes, zinc is bound to 4.1 +/- 0.4 N/O atoms at an average distance of 1.98 +/- 0.01 angstroms. Our previous studies demonstrated that 90% of the zinc in stage VI oocytes is sequestered within yolk platelets, associated with a single molecule, lipovitellin, the proteolytically processed product of vitellogenin. XAFS analysis of yolk platelets, lipovitellin, and vitellogenin demonstrates that zinc is bound to 4.0 +/- 0.5 N/O ligands at an average distance of 1.98 +/- 0.01 angstroms in each case, identical to that of stage V,VI oocytes. The higher shell contributions in the Fourier transforms indicate that two of the N/O zinc ligands are His in both stage V,VI and I,II oocytes. The results show that in stage I,II oocytes, there is a high concentration of a zinc protein whose zinc coordination site likely is composed of (His)2(Cys)2, such as, e.g., TFIIIA. As the oocytes develop, the predominant zinc species becomes one that exhibits the (His)2(N/0)2 zinc site found in lipovitellin. Hence, the ligands to the zinc atoms in intact oocytes and the changes that take place as a function of oogenesis and after their fertilization, during embryogenesis, now can be examined and explored.