Dataset Information

Ferrocene Orientation Determined Intramolecular Interactions Using Energy Decomposition Analysis.

ABSTRACT: Two very different quantum mechanically based energy decomposition analyses (EDA) schemes are employed to study the dominant energy differences between the eclipsed and staggered ferrocene conformers. One is the extended transition state (ETS) based on the Amsterdam Density Functional (ADF) package and the other is natural EDA (NEDA) based in the General Atomic and Molecular Electronic Structure System (GAMESS) package. It reveals that in addition to the model (theory and basis set), the fragmentation channels more significantly affect the interaction energy terms (ΔE) between the conformers. It is discovered that such an interaction energy can be absorbed into the pre-partitioned fragment channels so that to affect the interaction energies in a particular conformer of Fc. To avoid this, the present study employs a complete fragment channel-the fragments of ferrocene are individual neutral atoms. It therefore discovers that the major difference between the ferrocene conformers is due to the quantum mechanical Pauli repulsive energy and orbital attractive energy, leading to the eclipsed ferrocene the energy preferred structure. The NEDA scheme further indicates that the sum of attractive (negative) polarization (POL) and charge transfer (CL) energies prefers the eclipsed ferrocene. The repulsive (positive) deformation (DEF) energy, which is dominated by the cyclopentadienyle (Cp) rings, prefers the staggered ferrocene. Again, the cancellation results in a small energy residue in favour of the eclipsed ferrocene, in agreement with the ETS scheme. Further Natural Bond Orbital (NBO) analysis indicates that all NBO energies, total Lewis (no Fe) and lone pair (LP) deletion all prefer the eclipsed Fc conformer. The most significant energy preferring the eclipsed ferrocene without cancellation is the interactions between the donor lone pairs (LP) of the Fe atom and the acceptor antibond (BD*) NBOs of all C-C and C-H bonds in the ligand, LP(Fe)-BD*(C-C & C-H), which strongly stabilizes the eclipsed (D5h) conformation by -457.6 kcal·mol-1.

SUBMITTER: Wang F

PROVIDER: S-EPMC5458881 | biostudies-other | 2015 Nov

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:Energy decomposition analysis (EDA) is a well-established approach to dissect the interaction energy into chemically sound components. Despite the inherent requirement of reference states has been a long-standing object of debate, the direct relation with the molecular orbital analysis helps in building up predictive models. The alternative molecular energy decomposition schemes that decompose the total energy into atomic and diatomic contributions, such as the interacting quantum atoms (IQA), has no external reference requirements and also the intra- and intermolecular interactions are treated on equal footing. However, a connection with heuristic chemical models are limited, bringing about a somewhat narrower predictive power. While efforts to reconcile the bonding picture obtained by both methodologies have been discussed in the past, a synergic combination of them has not been tackled yet. Herein, we present the use of IQA decomposition of the individual terms arising from the EDA in the context of intermolecular interactions, henceforth EDA-IQA. The method is applied to a molecular set covering a wide range of interaction types, including hydrogen bonding, charge-dipole, π-π and halogen interactions. We find that the electrostatic energy from EDA, entirely seen as intermolecular, leads to meaningful and non-negligible intra-fragment contributions upon IQA decomposition, originated from charge penetration. EDA-IQA also affords the decomposition of the Pauli repulsion term into intra- and inter-fragment contributions. The intra-fragment term is destabilizing, particularly for the moieties that are net acceptors of charge, while the inter-fragment Pauli term is actually stabilizing. In the case of the orbital interaction term, the sign and magnitude of the intra-fragment contribution at equilibrium geometries is largely driven by the amount of charge transfer, while the inter-fragment contribution is clearly stabilizing. EDA-IQA terms show a smooth behavior along the intermolecular dissociation path of selected systems. The new EDA-IQA methodology provides a richer energy decomposition scheme that aims at bridging the gap between the two main distinct real-space and Hilbert-space methodologies. Via this approach, the partitioning can be used directionally on all the EDA terms aiding in identifying the causal effects on geometries and/or reactivity.

Project description:A detailed analysis of the composition and properties of hydrophobic nuclei and microclusters in pancreatic ribonuclease A (RNase A) has been carried out. Distance calculations for all noncovalently bonded atoms revealed that the average number of nonpolar contacts between a side chain of an amino acid and its neighbors is substantially larger if it involves hydrophobic residues rather than nonhydrophobic ones. However, the difference decreased when the number of contacts per nonpolar group and/or atom were calculated. Three main nuclei and five microclusters were identified, and their quantitative parameters were calculated. These nuclei include hydrophobic residues with a substantial number of nonpolar contacts with the environment (Phe 8, Phe 120, Phe 46, Tyr 25, Tyr 97, Ile 107, Leu 35, Ile 81, Val 54, Val 108, Met 29, Met 30). Hydrophobic nuclei of RNase A differ in shape and in composition, in the number of intranuclear contacts and of associated residues, as well as in their internal mobility. All eight cysteine residues are involved in nonpolar interactions with amino acid residues of hydrophobic nuclei. Active site amino acid residues of RNase A form a noncovalent contact network comprised of themselves, as well as of many conserved residues from hydrophobic nuclei. Sequence alignment with some other members of the RNase A family of proteins shows remarkable similarity in positions and in conservation of the main nonpolar residues, comprising cores of two (out of three) hydrophobic nuclei. A correlation was shown to exist between the average density of contacts for side-chain atoms and the number of amino acids to be found in the appropriate positions in the sequences of related mammalian ribonucleases. However, there are certain amino acid positions in the third, smaller nucleus, which are highly variable within the family. Taking into account that this nucleus is composed of residues belonging to different elements of the secondary structure, it is likely that the mutual orientation of these elements can be somehow different for these proteins.

Dataset Information

Ferrocene Orientation Determined Intramolecular Interactions Using Energy Decomposition Analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets