Energetics and dynamics in MbCN: CN--vibrational relaxation from molecular dynamics simulations

The dynamics of the cyanide anion bound to sperm-whale myoglobin is investigated using atomistic simulations. With density-functional theory, a 2D potential energy surface for the cyanide-heme complex is calculated. Two deep minima with a stabilization energy of approximately 50 kcal/mol corresponding to two different binding orientations (Fe-CN and Fe-NC) of the ligand are found. The Fe-CN conformation is favored over Fe-NC by several kcal/mol. Mixed quantum mechanics/molecular mechanics calculations show that the binding orientation affects the bond strength of the ligand, with a significantly different bond length and a 25 cm-1 shift in the fundamental CN-frequency. For the molecular dynamics (MD) simulations, a 3-center fluctuating charge model for the Fe-CN unit is developed that captures polarization and ligand-metal charge transfer. Stability arguments based on the energetics around the active site and the CN- frequency shifts suggest that the Fe-CN conformation with epsilon-protonation of His epsilon 64 are most likely, which is in agreement with experiment. Both equilibrium and nonequilibrium MD simulations are carried out to investigate the relaxation time scale and possible relaxation pathways in bound MbCN. The nonequilibrium MD simulations with a vibrationally excited ligand reveal that vibrational relaxation takes place on a time scale of hundreds of picoseconds within the active site. This finding supports the hypothesis that the experimentally observed relaxation rate (3.6 ps) reflects the repopulation of the electronic ground state.     

Active site structures and the redox properties of blue copper proteins: atomic resolution structure of azurin II and electronic structure calculations of azurin, plastocyanin and stellacyanin

Understanding how the active site structures of blue copper proteins determine their redox properties is the central structure-function relationship question of this important class of protein, also referred to as cupredoxins. We here describe both experimental and computational studies of azurin, plastocyanin and stellacyanin designed to define more accurately the geometric structures of the active site of the reduced and oxidized species, and thus to understand how these structures determine the redox potentials of these proteins. To this end the crystal structure of reduced azurin II has been determined at an atomic resolution of 1.13 Angstrom and is presented here. Co-ordinates and structure factors have been deposited in the RCSB Protein Data Bank with accession codes 2ccw and r2ccwsf respectively. The improved accuracy provided by the atomic resolution for the metal stereochemistry are utilised in conjunction with the EXAFS data for theoretical calculations. Multilevel calculations involving density functional theory and molecular mechanical potentials are used to predict both the geometric and electronic structure of the active sites of azurin, plastocyanin and stellacyanin and to estimate the relative redox potentials of these three proteins. We have also compared the relative energies of the structures obtained from experiment at varying resolutions, and from the isolated and embedded cluster calculations. We find significant energy differences between low and high (atomic) resolution structures arising primarily due to inaccuracies in the Cu-ligand distances in the lower resolution structures, emphasising the importance of accurate, very high resolution structural information. QM/MM structures are only approximately 1 kcal mol(-1) lower in energy than the 1.13 Angstrom structure while the optimized gas phase structure is 13.0 kcal mol(-1) lower in energy.     

AM1 calculations of bond dissociation energies. Allylic and benzylic C-H bonds

AM1 method and correlation dependence between electronic relaxation energy and valence change on the C atom of the breaking bond were used to calculate the bond dissociation energies in 50 compounds with allylic or benzylic C-H bonds. The average calculation error is 0.8 kcal/mol.     

Intramolecular hydrogen bond energy and cooperative interactions in α-, β-, and γ-cyclodextrin conformers

Accurate estimation of individual intramolecular hydrogen bond (H-bond) energies is an intricate task for multiply H-bonded systems. In such cases, the hydrogen bond strengths could be highly influenced by the cooperative interactions, for example, those between hydroxyl groups in sugars. In this work, we use the recently proposed molecular tailoring approach-based quantification (Deshmukh, Gadre, and Bartolotti, J Phys Chem A 2006, 110, 12519) to the extended systems of cyclodextrins (CDs). Further, the structure and stability of different conformers of α-, β-, and γ-CDs are explained based on the energetics and cooperative contribution to the strength of these H-bonds. The estimated O-H···O H-bond energies in the various CD conformers are found to vary widely from 1.1 to 8.3 kcal mol(-1). The calculated energy contributions to cooperativity toward the H-bond strengths fall in the range of 0.25-2.75 kcal mol(-1).     

Role of structural water molecule in HIV protease-inhibitor complexes: a QM/MM study

Structural water molecule 301 found at the interface of HIV protease-inhibitor complexes function as a hydrogen bond (H-bond) donor to carbonyl groups of the inhibitor as well as H-bond acceptor to amide/amine groups of the flap region of the protease. In this study, six systems of HIV protease-inhibitor complexes were analyzed, which have the presence of this "conserved" structural water molecule using a two-layer QM/MM ONIOM method. The combination of QM/MM and QM method enabled the calculation of strain energies of the bound ligands as well as the determination of their binding energies in the ligand-water and ligand-water-protease complexes. Although the ligand experiences considerable strain in the protein bound structure, the H-bond interactions through the structural water overcomes this strain effect to give a net stability in the range of 16-24 kcal/mol. For instance, in 1HIV system, the strain energy of the ligand was 12.2 kcal/mol, whereas the binding energy associated with the structural water molecule was 20.8 kcal/mol. In most of the cases, the calculated binding energy of structural water molecule showed the same trend as that of the experimental binding free energy values. Further, the classical MD simulations carried out on 1HVL system with and without structural water 301 showed that this conserved water molecule enhances the H-bond dynamics occurring at the Asp-bound active site region of the protease-inhibitor system, and therefore it will have a direct influence on the mechanism of drug action.     

Collision energy dependence of the O(1D) + HCl --> OH + Cl(2P) reaction studied by crossed beam scattering and quasiclassical trajectory calculations on ab initio potential energy surfaces

The dynamics of the O(1D) + HCl --> OH + Cl(2P) reaction are investigated by a crossed molecular beam ion-imaging method and quasiclassical trajectory calculations on the three ab initio potential energy surfaces, the ground 1(1)A' and two excited (1(1)A' and 2(1)A') states. The scattering experiment was carried out at collision energies of 4.2, 4.5, and 6.4 kcal/mol. The observed doubly differential cross sections (DCSs) for the Cl(2P) product exhibit almost no collision energy dependence over this inspected energy range. The nearly forward-backward symmetric DCS indicates that the reaction proceeds predominantly on the ground-state potential energy surface at these energies. Variation of the forward-backward asymmetry with collision energy is interpreted using an osculating complex model. Although the potential energy surfaces obtained by CASSCF-MRCI ab initio calculations exhibit relatively low potential barriers of 1.6 and 6.5 kcal/mol for 1(1)A' and 2(1)A', respectively, the dynamics calculations indicate that contributions of these excited states are small at the collision energies lower than 15.0 kcal/mol. Theoretical DCSs calculated for the ground-state reaction pathway agree well with the observed ones. These experimental and theoretical results suggest that the titled reaction at collision energies less than 6.5 kcal/mol is predominantly via the ground electronic state.     

1,3 Geminal interactions as the possible trend setting factors for C-H and C-C bond energies in alkanes. Support from a density functional theory based bond energy decomposition study

A bond energy decomposition analysis has been carried out to rationalize the well-established experimental fact that C-C and C-H bond energies decrease with increasing substitution on the carbon. It is shown that this trend is set by steric 1,3 repulsive interactions (geminal repulsion) that increase in the order 1,3 hydrogen-hydrogen < 1,3 hydrogen-carbon < 1,3 carbon-carbon. On the other hand, the radical stabilization energy has little influence on the observed trend for the C-H bond energy in H-CR(1)R(2)R(3) or the C-C bond energy in H(3)C-CR(1)R(2)R(3). Thus, it varies in H-CR(1)R(2)R(3) from -7.2 kcal/mol (H-CH3) to -6.5 kcal/mol (H-C(CH3)3) and in H(3)C-CR(1)R(2)R(3) from -19.0 kcal/mol (H(3)C-CH(3)) to -16.9 kcal/mol (H(3)C-C(CH(3))(3)). It was further found that the average intrinsic C-H bond energy in H-CR(1)R(2)R(3) of 129.2 kcal/mol is smaller than the average intrinsic C-C bond energy in H(3)C-CR(1)R(2)R(3) of 143.4 kcal/mol. However, after the inclusion of steric effects, the overall C-H bond becomes stronger than the C-C bond. The role of steric 1,3 repulsive interactions as the trend setting factor has most recently been suggested by Gronert (J. Org. Chem. 2006, 71, 1209) based on an empirical fit of alkane atomization energies.     

Diphtheria toxin catalyzed hydrolysis of NAD(+): molecular dynamics study of enzyme-bound substrate, transition state, and inhibitor.

The mechanism of the diphtheria toxin-catalyzed hydrolysis of NAD(+) was investigated by quantum chemical calculations and molecular dynamics simulations. Several effects that could explain the 6000-fold rate acceleration (Delta Delta G(++) approximately 5 kcal/mol) by the enzyme were considered. First, the carboxamide arm of the enzyme-bound NAD(+) adopts a trans conformation while the most stable conformation is cis. The most stable conformation for the nicotinamide product has the amide carbonyl trans. The activation energy for the cleavage of the ribosidic bond is reduced by 2 kcal/mol due to the relaxation of this ground state conformational stress in the transition state. Second, molecular dynamics simulations to the nanosecond time range revealed that the carboxylate of Glu148 forms a hydrogen bond to the substrate's 2' hydroxyl group in E.S (approximately 17% of the time) and E.TS (approximately 57% of the time) complexes. This interaction is not seen in crystal structures. The ApUp inhibitor is held more tightly by the enzyme than the transition state and the substrate. Analysis of correlated motions reveals differences in the pattern of anticorrelated motions for protein backbone atoms when the transition state occupies the active site as compared to the E.NAD(+) complex.     

Reliable structures and energetics for two new delocalized pi...pi prototypes: cyanogen dimer and diacetylene dimer

Two new prototype delocalized pi[dot dot dot]pi complexes are introduced: the dimers of cyanogen, (N[triple bond]C-C[triple bond]N)(2), and diacetylene, (HC[triple bond]C-C[triple bond]CH)(2). These dimers have properties similar to larger delocalized pi...pi systems such as benzene dimer but are small enough that they can be probed in far greater detail with high accuracy electronic structure methods. Parallel-slipped and T-shaped structures of both cyanogen dimer and diacetylene dimer have been optimized with 15 different procedures. The effects of basis set size, theoretical method, counterpoise correction, and the rigid monomer approximation on the structure and energetics of each dimer have been examined. MP2 and CCSD(T) optimized geometries for all four dimer structures are reported, as well as estimates of the CCSD(T) complete basis set (CBS) interaction energy for every optimized geometry. The data reported here suggest that future optimizations of delocalized pi[dot dot dot]pi clusters should be carried out with basis sets of triple-zeta quality. Larger basis sets and the expensive counterpoise correction to the molecular geometry are not necessary. The rigid monomer approximation has very little effect on structure and energetics of these dimers and may be used without consequence. Due to a consistent cancellation of errors, optimization with the MP2 method leads to CCSD(T)/CBS interaction energies that are within 0.2 kcal mol(-1) of those for structures optimized with the CCSD(T) method. Future studies that aim to resolve structures separated by a few tenths of a kcal mol(-1) should consider the effects of optimization with the CCSD(T) method.     

Theoretical investigation of the reactivity of copper atoms with OCS: comparison with CS2 and CO2

The reaction mechanism of the Cu atom with OCS and CO2 has been studied by means of density functional method (B3LYP). The overall energetics has been refined at the CCSD(T) level. In the case of the Cu + OCS reaction, the CS insertion route is found much more favorable than the CO insertion one. This later reaction is direct and involves an activation energy of 83.3 kcal/mol and is endothermic by 50.0 kcal/mol at the CCSD(T) level. The insertion into the CS bond proceeds through the eta1s and eta2cs coordination species as intermediates and is found exothermic by about 20 kcal/mol. The highest transition structure along this route is only 11.5 kcal/mol higher in energy than the reactant's ground states. In the case of the Cu + CO2 reaction, the insertion route into the CO bond is also found direct but with a lower endothermicity (30.6 kcal/mol) and smaller activation energy (61.1 kcal/mol) than that into the CO bond of OCS. In all cases, the insertion mechanism proceeds simultaneously with electron transfer from the Cu atom to OCS (or CO2) molecule.     

Estimates of the energy of intramolecular hydrogen bonds

A method for the estimation of the energy of intramolecular hydrogen bonds in conjugated systems existing in a variety of conformations is presented. The method is applied to determine the intramolecular hydrogen bond energy in 3-aminopropenal and 3-aminopropenthial. According to the proposed estimation scheme, the intramolecular H-bond energies are found to be of the order of 5-7 kcal/mol. These results are compared with those obtained by using other estimation schemes as well as with the recent results by other authors. Also, the H-bond energies in dimers and trimers of the two molecules are calculated and compared with the corresponding data for internally hydrogen-bonded monomers. This comparison shows that the bond equalization effect is primarily due to proton donor-proton acceptor proximity. In comparison with intermolecular hydrogen bonds, the rigidity of the chelate skeleton enhances this proximity effect. The same effect can be seen in systems with intermolecular hydrogen bonds, although its magnitude is diminished because of the absence of additional forces which pull the proton donor and proton acceptor groups toward each other. No specific resonance-assisted origin of the intramolecular hydrogen bond energy seems to be needed to elucidate the energetics of these bonds.     

Base-pair interactions in the gas-phase proton-bonded complexes of C+G and C+GC

Interactions involved in the formation of gas-phase proton-bonded molecular complexes of cytosine (C) and guanine (G) were theoretically investigated for the case of C(+)G and C(+)GC using B3LYP density functional theory. In this study, particular focus was on the dimeric interaction of proton-bonded C(+)G, where a proton bond and a hydrogen bond are cooperatively involved. The dimer interaction energy in terms of dissociation energy (D(e)) was predicted to be 41.8 kcal/mol. The lowest (frozen) energy structure for the C(+)G dimeric complex was found to be CH(+)...G rather than C...H(+)G in spite of the lower proton affinity of the cytosine moiety, which was more stable by 3.3 kcal/mol. The predicted harmonic vibrational frequencies and bond lengths suggest that the combined contributions of proton and hydrogen bonding may determine the resultant stability of each complex structure. In contrast to the dimer case, in the case of the isolated C(+)GC triplet, the two minimum energy structures of CH(+)...GC and C...H(+)GC were predicted to be almost equivalent in total energy. The dissociation energy (D(e)) for the C(+)G pairing in the C(+)GC triplet was 43.7 kcal/mol. Other energetics are also reported. As for the proton-transfer reaction in the proton-bond axis, the forward proton-transfer barriers for the dimer and trimer complexes were also predicted to be very low, 3.6 and 1.5 kcal/mol (DeltaE(e)(PT)), respectively.     

Determination of the C4-H bond dissociation energies of NADH models and their radical cations in acetonitrile

Heterolytic and homolytic bond dissociation energies of the C4-H bonds in ten NADH models (seven 1,4-dihydronicotinamide derivatives, two Hantzsch 1,4-dihydropyridine derivatives, and 9,10-dihydroacridine) and their radical cations in acetonitrile were evaluated by titration calorimetry and electrochemistry, according to the four thermodynamic cycles constructed from the reactions of the NADH models with N,N,N',N'-tetramethyl-p-phenylenediamine radical cation perchlorate in acetonitrile (note: C9-H bond rather than C4-H bond for 9,10-dihydroacridine; however, unless specified, the C9-H bond will be described as a C4-H bond for convenience). The results show that the energetic scales of the heterolytic and homolytic bond dissociation energies of the C4-H bonds cover ranges of 64.2-81.1 and 67.9-73.7 kcal mol(-1) for the neutral NADH models, respectively, and the energetic scales of the heterolytic and homolytic bond dissociation energies of the (C4-H)(.+) bonds cover ranges of 4.1-9.7 and 31.4-43.5 kcal mol(-1) for the radical cations of the NADH models, respectively. Detailed comparison of the two sets of C4-H bond dissociation energies in 1-benzyl-1,4-dihydronicotinamide (BNAH), Hantzsch 1,4-dihydropyridine (HEH), and 9,10-dihydroacridine (AcrH(2)) (as the three most typical NADH models) shows that for BNAH and AcrH(2), the heterolytic C4-H bond dissociation energies are smaller (by 3.62 kcal mol(-1)) and larger (by 7.4 kcal mol(-1)), respectively, than the corresponding homolytic C4-H bond dissociation energy. However, for HEH, the heterolytic C4-H bond dissociation energy (69.3 kcal mol(-1)) is very close to the corresponding homolytic C4-H bond dissociation energy (69.4 kcal mol(-1)). These results suggests that the hydride is released more easily than the corresponding hydrogen atom from BNAH and vice versa for AcrH(2), and that there are two almost equal possibilities for the hydride and the hydrogen atom transfers from HEH. Examination of the two sets of the (C4-H)(.+) bond dissociation energies shows that the homolytic (C4-H)(.+) bond dissociation energies are much larger than the corresponding heterolytic (C4-H)(.+) bond dissociation energies for the ten NADH models by 23.3-34.4 kcal mol(-1); this suggests that if the hydride transfer from the NADH models is initiated by a one-electron transfer, the proton transfer should be more likely to take place than the corresponding hydrogen atom transfer in the second step. In addition, some elusive structural information about the reaction intermediates of the NADH models was obtained by using Hammett-type linear free-energy analysis.     

Why are the Homoleptic Diyl Complexes M(InR)4 with M = Ni and Pt Stable Compounds while only Ni(CO)4 but not Pt(CO)4 can become Isolated? A Theoretical Study of M(EMe)44 and M(CO)4 (M = Ni,Pd, Pt; E = B,Al, Ga,In, Tl)

The equilibrium geometries and first bond dissociation energies of the homoleptic complexes M(EMe)₄ and M(CO)₄ with M = Ni, Pd, Pt and E = B, Al, Ga, In, Tl have been calculated at the gradient corrected DFT level using the BP86 functionals. The electronic structure of the metal‐ligand bonds has been examined with the topologial analysis of the electron density distribution. The nature of the bonding is revealed by partitioning the metal‐ligand interaction energies into contributions by electrostatic attraction, covalent bonding and Pauli repulsion. The calculated data show that the M‐CO and M‐EMe bonding is very similar. However, the M‐EMe bonds of the lighter elements E are much stronger than the M‐CO bonds. The bond energies of the latter are as low or even lower than the M‐TlMe bonds. The main reason why Pd(CO)₄ and Pt(CO)₄ are unstable at room temperature in a condensed phase can be traced back to the already rather weak bond energy of the Ni‐CO bond. The Pd‐L bond energies of the complexes with L = CO and L = EMe are always 10 — 20 kcal/mol lower than the Ni‐L bond energies. The calculated bond energy of Ni(CO)₄ is only D_o = 27 kcal/mol. Thus, the bond energy of Pd(CO)₄ is only D_o = 12 kcal/mol. The first bond dissociation energy of Pt(CO)₄ is low because the relaxation energy of the Pt(CO)₃ fragment is rather high. The low bond energies of the M‐CO bonds are mainly caused by the relatively weak electrostatic attraction and by the comparatively large Pauli repulsion. The σ and π contributions to the covalent M‐CO interactions have about the same strength. The π bonding in the M‐EMe bonds is less than in the M‐CO bonds but it remains an important part of the bond energy. The trends of the electrostatic and covalent contributions to the bond energies and the σ and π bonding in the metal‐ligand bonds are discussed.     

Pi and sigma-phenylethynyl radicals and their isomers o-, m-, and p-ethynylphenyl: structures, energetics, and electron affinities

Molecular structures, energetics, vibrational frequencies, and electron affinities are predicted for the phenylethynyl radical and its isomers. Electron affinities are computed using density functional theory, -namely, the BHLYP, BLYP, B3LYP, BP86, BPW91, and B3PW91 functionals-, employing the double-zeta plus polarization DZP++ basis set; this level of theory is known to perform well for the computation of electron affinities. Furthermore, ab initio computations employing perturbation theory, coupled cluster with single and double excitations [CCSD], and the inclusion of perturbative triples [CCSD(T)] are performed to determine the relative energies of the isomers. These higher level computations are performed with the correlation consistent family of basis sets cc-pVXZ (X = D, T, Q, 5). Three electronic states are probed for the phenylethynyl radical. In C2v symmetry, the out-of-plane (2B1) radical is predicted to lie about 10 kcal/mol below the in-plane (2B2) radical by DFT methods, which becomes 9.4 kcal/mol with the consideration of the CCSD(T) method. The energy difference between the lowest pi and sigma electronic states of the phenylethynyl radical is also about 10 kcal/mol according to DFT; however, CCSD(T) with the cc-pVQZ basis set shows this energy separation to be just 1.8 kcal/mol. The theoretical electron affinities of the phenylethynyl radical are predicted to be 3.00 eV (B3LYP/DZP++) and 3.03 eV (CCSD(T)/DZP++//MP2/DZP++). The adiabatic electron affinities (EAad) of the three isomers of phenylethynyl, that is, the ortho-, meta-, and para-ethynylphenyl, are predicted to be 1.45, 1.40, and 1.43 eV, respectively. Hence, the phenylethynyl radical binds an electron far more effectively than the three other radicals studied. Thermochemical predictions, such as the bond dissociation energies of the aromatic and ethynyl C-H bonds and the proton affinities of the phenylethynyl and ethynylphenyl anions, are also reported.     

A theoretical study of surface reduction mechanisms of CeO(2)(111) and (110) by H(2).

Reaction mechanisms for the interactions between CeO(2)(111) and (110) surfaces are investigated using periodic density functional theory (DFT) calculations. Both standard DFT and DFT+U calculations to examine the effect of the localization of Ce 4f states on the redox chemistry of H(2)-CeO(2) interactions are described. For mechanistic studies, molecular and dissociative local minima are initially located by placing an H(2) molecule at various active sites of the CeO(2) surfaces. The binding energies of physisorbed species optimized using the DFT and DFT+U methods are very weak. The dissociative adsorption reactions producing hydroxylated surfaces are all exothermic; exothermicities at the DFT level range from 4.1 kcal mol(-1) for the (111) to 26.5 kcal mol(-1) for the (110) surface, while those at the DFT+U level are between 65.0 kcal mol(-1) for the (111) and 81.8 kcal mol(-1) for the (110) surface. Predicted vibrational frequencies of adsorbed OH and H(2)O species on the surfaces are in line with available experimental and theoretical results. Potential energy profiles are constructed by connecting molecularly adsorbed and dissociatively adsorbed intermediates on each CeO(2) surface with tight transition states using the nudged elastic band (NEB) method. It is found that the U correction method plays a significant role in energetics, especially for the intermediates of the exit channels and products that are partially reduced. The surface reduction reaction on CeO(2)(110) is energetically much more favorable. Accordingly, oxygen vacancies are more easily formed on the (110) surface than on the (111) surface.     

Role of protein structure and the role of individual fingers in zinc finger protein–DNA recognition: a molecular dynamics simulation study and free energy calculations

Molecular dynamics and MM_GBSA energy calculations on various zinc finger proteins containing three and four fingers bound to their target DNA gave insights into the role of each finger in the DNA binding process as part of the protein structure. The wild type Zif 268 (PDB code: 1AAY) gave a ΔG value of ??76.1 (14) kcal/mol. Zinc fingers ZF1, ZF2 and ZF3 were mutated in one experiment and in another experiment one finger was cut and the rest of the protein was studied for binding. The ΔΔG values for the Zinc Finger protein with both ZF1 and ZF2 mutated was +?80 kcal/mol, while mutating only ZF1 the ΔΔG value was +?52 kcal/mol (relative to the wild type). Cutting ZF3 and studying the protein consisting only of ZF1 linked to ZF2 gave a ΔΔG value of +?68 kcal/mol. Upon cutting ZF1, the resulting ZF2 linked to ZF3 protein gave a ΔΔG value of +?41 kcal/mol. The above results shed light on the importance of each finger in the binding process, especially the role of ZF1 as the anchoring finger followed in importance by ZF2 and ZF3. The energy difference between the binding of the wild type protein Zif268 (1AAY) and that for individual finger binding to DNA according to the formula: ΔΔG_{linkers, otherstructuralfactors}?=?ΔG_zif268???(ΔG_F1+F2+F3) gave a value?=???44.5 kcal/mol. This stabilization can be attributed to the contribution of linkers and other structural factors in the intact protein in the DNA binding process. DNA binding energies of variant proteins of the wild type Zif268 which differ in their ZF1 amino acid sequence gave evidence of a good relationship between binding energy and recognition and specificity, this finding confirms the reported vital role of ZF1 in the ZF protein scanning and anchoring to the target DNA sequence. The role of hydrogen bonds in both specific and nonspecific amino acid-DNA contacts is discussed in relation to mutations. The binding energies of variant Zinc Finger proteins confirmed the role of ZF1 in the recognition, specificity and anchoring of the zinc finger protein to DNA.     

Homonuclear transition-metal trimers

Density-functional theory has been used to determine the ground-state geometries and electronic states for homonuclear transition-metal trimers constrained to equilateral triangle geometries. This represents the first application of consistent theoretical methods to all of the ten 3d block transition-metal trimers, from scandium to zinc. A search of the potential surfaces yields the following electronic ground states and bond lengths: Sc3(2A1',2.83 A), Ti3(7E',2.32 A), V3(2E",2.06 A), Cr3(17E',2.92 A), Mn3(16A2',2.73 A), Fe3(11E",2.24 A), Co3(6E",2.18 A), Ni3(3A2",2.23 A), Cu3(2E',2.37 A), and Zn3(1A1',2.93 A). Vibrational frequencies, several low-lying electronic states, and trends in bond lengths and atomization energies are discussed. The predicted dissociation energies DeltaE(M3-->M2+M) are 49.4 kcal mol(-1)(Sc3), 64.3 kcal mol(-1)(Ti3), 60.7 kcal mol(-1)(V3), 11.5 kcal mol(-1)(Cr3), 32.4 kcal mol(-1)(Mn3), 61.5 kcal mol(-1)(Fe3), 78.0 kcal mol(-1)(Co3), 86.1 kcal mol(-1)(Ni3), 26.8 kcal mol(-1)(Cu3), and 4.5 kcal mol(-1)(Zn3).