PubMed Health⌕ Search

Biomedical subjects

Yirong Mo

Publications and source records attributed to Yirong Mo.

9 recordsLinked to original sources
PubMed2005

Y-aromaticity: why is the trimethylenemethane dication more stable than the butadienyl dication?

[reactions: see text] Resonance energies of the trimethylenemethane dication (1) and the butadienyl dication (4) were evaluated using two independent computational methodologies to provide insight into the validity of Y-aromaticity. One methodology employed density functional theory calculations and examined the resonance contribution of the C=C double bond toward the double hydride abstraction enthalpies of methylpropene (6) and 2-butene (8), yielding 1 and 4, respectively. These resonance contributions by the double bond were determined by calculating the double hydride abstraction enthalpies of both the parallel and perpendicular conformations of vinylogues of 6 and 8, in which n = 1-4 vinyl units were inserted between the central carbon-carbon double bond and each of the reaction centers. Extrapolation of the resonance contribution in each vinylogue to n = 0 yielded the resonance contribution in the respective parent molecules. The second methodology employed an orbital deletion procedure (ODP), which effectively allowed us to examine the energies of individual resonance structures. The resonance energy of each dication is computed as the difference between the most stable resonance structure and that of the delocalized species. The two methodologies are in agreement, suggesting that the resonance energy of the trimethylenemethane dication is substantially greater than that of the butadienyl dication. The origin of this difference in resonance stabilization is discussed.

Journal Article↗
PubMed2005

Density functional study on dihydrogen activation at the H cluster in Fe-only hydrogenases.

Models simulating the catalytic diiron subcluster [FeFe](H) in Fe-only hydrogenases have often been designed for computational exploration of the catalytic mechanism of the formation and cleavage of dihydrogen. In this work, we extended the above models by explicitly considering the electron reservoir [4Fe-4S](H) which is linked to the diiron subcluster to form a whole H cluster ([6Fe-6S] = [4Fe-4S](H) + [FeFe](H)). Large-scale density functional theory (DFT) computations on the complete H cluster, together with simplified models in which the [4Fe-4S](H) subcluster is not directly involved in the reaction processes, have been performed to probe hydrogen activation on the Fe-only hydrogenases. A new intermediate state containing an Fe(p)...H...CN two-electron three-center bond is identified as a key player in the H2 formation/cleavage processes.

Binding Sites↗
PubMed2005

XMVB: a program for ab initio nonorthogonal valence bond computations.

An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available.

Journal Article↗
PubMed2004

Resonance effect in the allyl cation and anion: a revisit.

The interest over the magnitude of the conjugation effect in the allyl cation (1) and anion (2) has been revived recently by Barbour and Karty (J. Org. Chem. 2004, 69, 648-654), who derived the resonance energies of 20-22 and 17-18 kcal/mol for 1 and 2, respectively, using an empirical extrapolation approximation. This paper revisits the case by explicitly calculating the Pauling-Wheland resonance energy, which measures the stabilization from the most stable resonance structure to the delocalized energy-minimum state of a conjugated system, using our newly developed block-localized wave function (BLW) method. This BLW method has the geometrical optimization capability. The computations result in adiabatic resonance energies of 37 kcal/mol for 1 and 38 kcal/mol for 2. The significant disagreement between these values and Barbour and Karty's results originates from the neglect of structural and electronic variations in their derivation which are energy costing.

Journal Article↗
PubMed2004

Hyperconjugation effect in substituted methyl boranes: an orbital deletion procedure analysis.

The hyperconjugation effect in the substituted methyl boranes, XCH(2)BH(2) (X = H, CH(3), NH(2), PH(2), OH, SH, F, Cl, Br), has been quantitatively evaluated by using the orbital deletion procedure (ODP), where the p(pi) orbital on boron is deactivated. Except for the case of X = NH(2), which forms a three-membered ring, the magnitude of the hyperconjugative stabilization in all other substituted methylborane ranges from 6.8 to 3.4 kcal/mol. Significant structural changes are observed, particularly the shortening of the central B-C bond distance and the reducing of the corresponding XCB and HCB bond angles. In general, the strength of the hyperconjugative interaction between the occupied sigma(C-X) bond and the vacant p(pi) orbital on boron is correlated to the electronegativity of X, and the competition between the donation ability of the sigma(C-X) and the sigma(C-H) bonds determines the preference of the staggered or eclipsed structure as the energy minimum state. When the donation abilities of the C-X and C-H bonds are comparable, other factors such as electron correlation and steric effect may play elaborate roles in the geometrical propensity of the most stable structures.

Journal Article↗
PubMed2004

Charge transfer in the electron donor-acceptor complex BH3NH3.

As a simple yet strongly binding electron donor-acceptor (EDA) complex, BH(3)NH(3) serves as a good example to study the electron pair donor-acceptor complexes. We employed both the ab initio valence bond (VB) and block-localized wave function (BLW) methods to explore the electron transfer from NH(3) to BH(3). Conventionally, EDA complexes have been described by two diabatic states: one neutral state and one ionic charge-transferred state. Ab initio VB self-consistent field (VBSCF) computations generate the energy profiles of the two diabatic states together with the adiabatic (ground) state. Our calculations evidently demonstrated that the electron transfer between NH(3) and BH(3) falls in the abnormal regime where the reorganization energy is less than the exoergicity of the reaction. The nature of the NH(3)-BH(3) interaction is probed by an energy decomposition scheme based on the BLW method. We found that the variation of the charge-transfer energy with the donor-acceptor distance is insensitive to the computation levels and basis sets, but the estimation of the amount of electron transferred heavily depends on the population analysis procedures. The recent resurgence of interest in the nature of the rotation barrier in ethane prompted us to analyze the conformational change of BH(3)NH(3), which is an isoelectronic system with ethane. We found that the preference of the staggered structure over the eclipsed structure of BH(3)NH(3) is dominated by the Pauli exchange repulsion.

Journal Article↗
PubMed2004

Enzymatic mechanism of Fe-only hydrogenase: density functional study on H-H making/breaking at the diiron cluster with concerted proton and electron transfers.

The mechanism of the enzymatic hydrogen bond forming/breaking (2H(+) + 2e<==>H(2)) and the plausible charge and spin states of the catalytic diiron subcluster [FeFe](H) of the H cluster in Fe-only hydrogenases are probed computationally by the density functional theory. It is found that the active center [FeFe](H) can be rationally simulated as [[H](CH(3)S)(CO)(CN(-))Fe(p)(CO(b))(mu-SRS)Fe(d)(CO)(CN(-))L], where the monovalence [H] stands for the [4Fe4S](H)(2+) subcluster bridged to the [FeFe](H) moiety, (CH(3)S) represents a Cys-S, and (CO(b)) represents a bridging CO. L could be a CO, H(2)O, H(-), H(2), or a vacant coordination site on Fe(d). Model structures of possible redox states are optimized and compared with the X-ray crystallographic structures and FTIR experimental data. On the basis of the optimal structures, we study the most favorable path of concerted proton transfer and electron transfer in H(2)-forming/breaking reactions at [FeFe](H). Previous mechanisms derived from quantum chemical computations of Fe-only hydrogenases (Cao, Z.; Hall, M. B. J. Am. Chem. Soc. 2001, 123, 3734; Fan, H.; Hall, M. B. J. Am. Chem. Soc. 2001, 123, 3828) involved an unidentified bridging residue (mu-SRS), which is either a propanedithiolate or dithiomethylamine. Our proposed mechanism, however, does not require such a ligand but makes use of a shuttle of oxidation states of the iron atoms and a reaction site between the two iron atoms. Therefore, the hydride H(b)(-) (bridged to Fe(p) and Fe(d)) and eta(2)-H(2) at Fe(p) or Fe(d) most possibly play key roles in the dihydrogen reversible oxidation at the [FeFe](H) active center. This suggested way of H(2) formation/splitting is reminiscent of the mechanism of [NiFe] hydrogenases and therefore would unify the mechanisms of the two related enzymes.

Catalysis↗
PubMed2002

Cation-pi interactions: an energy decomposition analysis and its implication in delta-opioid receptor-ligand binding.

The nature and strength of the cation-pi interaction in protein-ligand binding are modeled by considering a series of nonbonded complexes involving N-substituted piperidines and substituted monocylic aromatics that mimic the delta-opioid receptor-ligand binding. High-level ab initio quantum mechanical calculations confirm the importance of such cation-pi interactions, whose intermolecular interaction energy ranges from -6 to -12 kcal/mol. A better understanding of the electrostatics, polarization, and other intermolecular interactions is obtained by appropriately decomposing the total interaction energy into their individual components. The energy decomposition analysis is also useful for parametrizing existing molecular mechanics force fields that could then account for energetic contributions arising out of cation-pi interactions in biomolecules. The present results further provide a framework for interpreting experimental results from point mutation reported for the delta-opioid receptor.

Cations↗