Ab initio calculation of electronic coupling in the photosynthetic reaction center

We have carried out an ab initio electronic structure calculations of electron transfer couplings between chromophores in the bacterial photosynthetic reaction center. The couplings agree remarkably well with parameters obtained from recent quantum dynamical modeling of experimental data assuming an explicit intermediate mechanism. We also have computed couplings on the M-side of the reaction center and have found that the interaction of the primary donor to the M-side intermediate bacteriochlorophyll is quite small because of destructive interference of the two localized coupling matrix elements. This may explain the slow rate of electron transfer down the M-side of the reaction center.

From ab initio quantum mechanics to molecular neurobiology: A cation–π binding site in the nicotinic receptor

The nicotinic acetylcholine receptor is the prototype ligand-gated ion channel. A number of aromatic amino acids have been identified as contributing to the agonist binding site, suggesting that cation–π interactions may be involved in binding the quaternary ammonium group of the agonist, acetylcholine. Here we show a compelling correlation between: (i) ab initio quantum mechanical predictions of cation–π binding abilities and (ii) EC50 values for acetylcholine at the receptor for a series of tryptophan derivatives that were incorporated into the receptor by using the in vivo nonsense-suppression method for unnatural amino acid incorporation. Such a correlation is seen at one, and only one, of the aromatic residues—tryptophan-149 of the α subunit. This finding indicates that, on binding, the cationic, quaternary ammonium group of acetylcholine makes van der Waals contact with the indole side chain of α tryptophan-149, providing the most precise structural information to date on this receptor. Consistent with this model, a tethered quaternary ammonium group emanating from position α149 produces a constitutively active receptor.

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.

Theoretical investigation of the [1,2]-sigmatropic hydrogen migration in the mechanism of oxidation of 2-aminobenzoyl-CoA by 2-aminobenzoyl-CoA monooxygenase/reductase

The flavin hydroperoxide at the active site of the mixed-function oxidase 2-aminobenzoyl-CoA monooxygenase/reductase (Azoarcus evansii) transfers an oxygen to the 5-position of the 2-aminobenzoyl-CoA substrate to provide the alkoxide intermediate II−. Hydrogen migration from C5 to C6 follows this monooxygenation. The nature of the monooxygenation intermediate and plausible competing reactions leading to hydrogen migration have been considered. Ab initio molecular orbital theory has been used to calculate structures and electron distributions in intermediate and transition state structures. Electrostatic potential surface calculations establish that the transition state and product, associated with the C5 to C6 hydrogen transfer, are stabilized by electron distribution to the benzoyl-CoA thioester carbonyl oxygen. This is not so for the transition state and product associated with hydrogen transfer from C5 to C4. The activation energy for the 5,6-shift is 2.5 kcal/mol lower than that for the 5,4-shift. In addition, the product of the hydrogen 5,6-shift is more stable than is the product of the hydrogen 5,4-shift, by ≈6 kcal/mol. These results explain why only the shift of hydrogen from C5 to C6 is observed experimentally. Oxygen transfer and hydrogen migration almost coincide in the gas phase (activation energy of ≈0.6 kcal/mol, equivalent to a single bond vibration). Enzymatic formation of alkoxide II− requires its stabilization; thus, the rate constant for its breakdown would be slower than in the gas phase.

Ionophores and receptors using cation-π interactions: Collarenes

Cation-π interactions are important forces in molecular recognition by biological receptors, enzyme catalysis, and crystal engineering. We have harnessed these interactions in designing molecular systems with circular arrangement of benzene units that are capable of acting as ionophores and models for biological receptors. [n]Collarenes are promising candidates with high selectivity for a specific cation, depending on n, because of their structural rigidity and well-defined cavity size. The interaction energies of [n]collarenes with cations have been evaluated by using ab initio calculations. The selectivity of these [n]collarenes in aqueous solution was revealed by using statistical perturbation theory in conjunction with Monte Carlo and molecular dynamics simulations. It has been observed that in [n]collarenes the ratio of the interaction energies of a cation with it and the cation with the basic building unit (benzene) can be correlated to its ion selectivity. We find that collarenes are excellent and efficient ionophores that bind cations through cation-π interactions. [6]Collarene is found to be a selective host for Li+ and Mg2+, [8]collarene for K+ and Sr2+, and [10]collarene for Cs+ and Ba2+. This finding indicates that [10]collarene and [8]collarene could be used for effective separation of highly radioactive isotopes, 137Cs and 90Sr, which are major constituents of nuclear wastes. More interestingly, collarenes of larger cavity size can be useful in capturing organic cations. [12]Collarene exhibits a pronounced affinity for tetramethylammonium cation and acetylcholine, which implies that it could serve as a model for acetylcholinestrase. Thus, collarenes can prove to be novel and effective ionophores/model-receptors capable of heralding a new direction in molecular recognition and host-guest chemistry.

Non-enzymatic and enzymatic hydrolysis of alkyl halides: A haloalkane dehalogenation enzyme evolved to stabilize the gas-phase transition state of an SN2 displacement reaction

The semiempirical PM3 method, calibrated against ab initio HF/6–31+G(d) theory, has been used to elucidate the reaction of 1,2-dichloroethane (DCE) with the carboxylate of Asp-124 at the active site of haloalkane dehalogenase of Xanthobacter autothropicus. Asp-124 and 13 other amino acid side chains that make up the active site cavity (Glu-56, Trp-125, Phe-128, Phe-172, Trp-175, Leu-179, Val-219, Phe-222, Pro-223, Val-226, Leu-262, Leu-263, and His-289) were included in the calculations. The three most significant observations of the present study are that: (i) the DCE substrate and Asp-124 carboxylate, in the reactive ES complex, are present as an ion-molecule complex with a structure similar to that seen in the gas-phase reaction of AcO− with DCE; (ii) the structures of the transition states in the gas-phase and enzymatic reaction are much the same where the structure formed at the active site is somewhat exploded; and (iii) the enthalpies in going from ground states to transition states in the enzymatic and gas-phase reactions differ by only a couple kcal/mol. The dehalogenase derives its catalytic power from: (i) bringing the electrophile and nucleophile together in a low-dielectric environment in an orientation that allows the reaction to occur without much structural reorganization; (ii) desolvation; and (iii) stabilizing the leaving chloride anion by Trp-125 and Trp-175 through hydrogen bonding.

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.

Femtosecond observation of benzyne intermediates in a molecular beam: Bergman rearrangement in the isolated molecule

In this communication, we report our femtosecond real-time observation of the dynamics for the three didehydrobenzene molecules (p-, m-, and o-benzyne) generated from 1,4-, 1,3-, and 1,2-dibromobenzene, respectively, in a molecular beam, by using femtosecond time-resolved mass spectrometry. The time required for the first and the second C-Br bond breakage is less than 100 fs; the benzyne molecules are produced within 100 fs and then decay with a lifetime of 400 ps or more. Density functional theory and high-level ab initio calculations are also reported herein to elucidate the energetics along the reaction path. We discuss the dynamics and possible reaction mechanisms for the disappearance of benzyne intermediates. Our effort focuses on the isolated molecule dynamics of the three isomers on the femtosecond time scale.

On the absorbance changes in the photocycle of the photoactive yellow protein: A quantum-chemical analysis

Spectral changes in the photocycle of the photoactive yellow protein (PYP) are investigated by using ab initio multiconfigurational second-order perturbation theory at the available structures experimentally determined. Using the dark ground-state crystal structure [Genick, U. K., Soltis, S. M., Kuhn, P., Canestrelli, I. L. & Getzoff, E. D. (1998) Nature (London) 392, 206–209], the ππ* transition to the lowest excited state is related to the typical blue-light absorption observed at 446 nm. The different nature of the second excited state (nπ*) is consistent with the alternative route detected at 395-nm excitation. The results suggest the low-temperature photoproduct PYPHL as the most plausible candidate for the assignment of the cryogenically trapped early intermediate (Genick et al.). We cannot establish, however, a successful correspondence between the theoretical spectrum for the nanosecond time-resolved x-ray structure [Perman, B., Šrajer, V., Ren, Z., Teng, T., Pradervand, C., et al. (1998) Science 279, 1946–1950] and any of the spectroscopic photoproducts known up to date. It is fully confirmed that the colorless light-activated intermediate recorded by millisecond time-resolved crystallography [Genick, U. K., Borgstahl, G. E. O., Ng, K., Ren, Z., Pradervand, C., et al. (1997) Science 275, 1471–1475] is protonated, nicely matching the spectroscopic features of the photoproduct PYPM. The overall contribution demonstrates that a combined analysis of high-level theoretical results and experimental data can be of great value to perform assignments of detected intermediates in a photocycle.

An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

Indium-cadmium-oxide films having exceptional electrical conductivity and optical transparency: Clues for optimizing transparent conductors

Materials with high electrical conductivity and optical transparency are needed for future flat panel display, solar energy, and other opto-electronic technologies. InxCd1-xO films having a simple cubic microstructure have been grown on amorphous glass substrates by a straightforward chemical vapor deposition process. The x = 0.05 film conductivity of 17,000 S/cm, carrier mobility of 70 cm2/Vs, and visible region optical transparency window considerably exceed the corresponding parameters for commercial indium-tin oxide. Ab initio electronic structure calculations reveal small conduction electron effective masses, a dramatic shift of the CdO band gap with doping, and a conduction band hybridization gap caused by extensive Cd 5s + In 5s mixing.

Conformation of coenzyme pyrroloquinoline quinone and role of Ca2+ in the catalytic mechanism of quinoprotein methanol dehydrogenase

The ab initio structures of 2,7,9-tricarboxypyrroloquinoline quinone (PQQ), semiquinone (PQQH), and dihydroquinone (PQQH2) have been determined and compared with ab initio structures of the (PQQ)Ca2+, (PQQH)Ca2+, and (PQQH2)Ca2+ complexes as well as the x-ray structure of (PQQ)Ca2+ bound at the active site of the methanol dehydrogenase (MDH) of methyltropic bacteria. Plausible mechanisms for the MDH oxidation of methanol involving the (PQQ)Ca2+ complex are explored via ab initio computations and discussed. Considering the reaction of methanol with PQQ in the absence of Ca2+, nucleophilic addition of methanol to the PQQ C-5 carbonyl followed by a retro-ene elimination is deemed unlikely due to large energy barrier. A much more favorable disposition of the methanol C-5 adduct to provide formaldehyde involves proton ionization of the intermediate followed by elimination of methoxide concerted with hydride transfer to the oxygen of the C-4 carbonyl. Much the same transition state is reached if one searches for the transition state beginning with Asp-303–CO2−general-base removal of the methanol proton of the (PQQ)Ca2+O(H)CH3 complex concerted with hydride transfer to the oxygen at C-4. For such a mechanism the role of the Ca2+ moiety would be to (i) contribute to the formation of the ES complex (ii) provide a modest decrease in the pKa of methanol substrate,; and (iii) polarize the oxygen at C-5.

Site-specific conformational determination in thermal unfolding studies of helical peptides using vibrational circular dichroism with isotopic substitution

Understanding the detailed mechanism of protein folding requires dynamic, site-specific stereochemical information. The short time response of vibrational spectroscopies allows evaluation of the distribution of populations in rapid equilibrium as the peptide unfolds. Spectral shifts associated with isotopic labels along with local stereochemical sensitivity of vibrational circular dichroism (VCD) allow determination of the segment sequence of unfolding. For a series of alanine-rich peptides that form α-helices in aqueous solution, we used isotopic labeling and VCD to demonstrate that the α-helix noncooperatively unwinds from the ends with increasing temperature. For these blocked peptides, the C-terminal is frayed at 5°C. Ab initio level theoretical simulations of the IR and VCD band shapes are used to analyze the spectra and to confirm the conformation of the labeled components. The VCD signals associated with the labeled residues are amplified by coupling to the nonlabeled parts of the molecule. Thus small labeled segments are detectable and stereochemically defined in moderately large peptides in this report of site-specific peptide VCD conformational analysis.

RefSeq and LocusLink: NCBI gene-centered resources

Thousands of genes have been painstakingly identified and characterized a few genes at a time. Many thousands more are being predicted by large scale cDNA and genomic sequencing projects, with levels of evidence ranging from supporting mRNA sequence and comparative genomics to computing ab initio models. This, coupled with the burgeoning scientific literature, makes it critical to have a comprehensive directory for genes and reference sequences for key genomes. The NCBI provides two resources, LocusLink and RefSeq, to meet these needs. LocusLink organizes information around genes to generate a central hub for accessing gene-specific information for fruit fly, human, mouse, rat and zebrafish. RefSeq provides reference sequence standards for genomes, transcripts and proteins; human, mouse and rat mRNA RefSeqs, and their corresponding proteins, are discussed here. Together, RefSeq and LocusLink provide a non-redundant view of genes and other loci to support research on genes and gene families, variation, gene expression and genome annotation. Additional information about LocusLink and RefSeq is available at http://www.ncbi.nlm.nih.gov/LocusLink/.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.