The Kubas interaction in M(II) (M = Ti, V, Cr) hydrazine-based hydrogen storage materials: a DFT study

The Cr(II) binding sites of an experimentally realised hydrazine linked hydrogen storage material have been studied computationally using density functional theory. Both the experimentally determined rise in H(2) binding enthalpy upon alteration of the ancillary ligand from bis[(trimethylsilyl)methyl] to hydride, and the number of H(2) molecules per Cr centre, are reproduced reasonably well. Comparison with analogous Ti(II), V(II) and Mn(II) systems suggests that future experiments should focus on the earliest 3d metals, and also suggests that 5 and 7 wt% H(2) storage may be possible for V(II) and Ti(II) respectively. Alteration of the metal does not have a large effect on the M-H(2) interaction energy, while alteration of the ancillary ligand bound to the metal centre, from bis[(trimethylsilyl)methyl] or hydride to two hydride ligands, THF and only hydrazine based ligands, indicates that ancillary ligands that are poor π-acceptors give stronger M-H(2) interactions. Good evidence is found that the M-H(2) interaction is Kubas type. Orbitals showing σ-donation from H(2) to the metal and π-back-donation from the metal to the dihydrogen are identified, and atoms-in-molecules analysis indicates that the electron density at the bond critical points of the bound H(2) is similar to that of classical Kubas systems. The Kubas interaction is dominated by σ-donation from the H(2) to the metal for Cr(II), but is more balanced between σ-donation and π-back-donation for the Ti(II) and V(II) analogues. This difference in behaviour is traced to a lowering in energy of the metal 3d orbitals across the transition series.     

Hydride-induced amplification of performance and binding enthalpies in chromium hydrazide gels for Kubas-type hydrogen storage

Hydrogen is the ideal fuel because it contains the most energy per gram of any chemical substance and forms water as the only byproduct of consumption. However, storage still remains a formidable challenge because of the thermodynamic and kinetic issues encountered when binding hydrogen to a carrier. In this study, we demonstrate how the principal binding sites in a new class of hydrogen storage materials based on the Kubas interaction can be tuned by variation of the coordination sphere about the metal to dramatically increase the binding enthalpies and performance, while also avoiding the shortcomings of hydrides and physisorpion materials, which have dominated most research to date. This was accomplished through hydrogenation of chromium alkyl hydrazide gels, synthesized from bis(trimethylsilylmethyl) chromium and hydrazine, to form materials with low-coordinate Cr hydride centers as the principal H(2) binding sites, thus exploiting the fact that metal hydrides form stronger Kubas interactions than the corresponding metal alkyls. This led to up to a 6-fold increase in storage capacity at room temperature. The material with the highest capacity has an excess reversible storage of 3.23 wt % at 298 K and 170 bar without saturation, corresponding to 40.8 kg H(2)/m(3), comparable to the 2015 DOE system goal for volumetric density (40 kg/m(3)) at a safe operating pressure. These materials possess linear isotherms and enthalpies that rise on coverage, retain up to 100% of their adsorption capacities on warming from 77 to 298 K, and have no kinetic barrier to adsorption or desorption. In a practical system, these materials would use pressure instead of temperature as a toggle and can thus be used in compressed gas tanks, currently employed in the majority of hydrogen test vehicles, to dramatically increase the amount of hydrogen stored, and therefore range of any vehicle.     

Mononuclear homoleptic allyl complexes of the first row transition metals: species with unusual metal electronic configurations

A number of evanescent unsubstituted homoleptic allyl derivatives M(C(3)H(5))(n) of the first row transition metals have been reported in the literature. In addition, the much more thermally stable silylated derivatives M[C(3)H(3)(SiMe(3))(2)](2) (M = Cr, Fe, Co, Ni) are reported to survive vacuum sublimation without significant decomposition. In this connection, the complete series of homoleptic allyl derivatives M(C(3)H(5))(n) (n = 2, 3; M = Sc, Ti, V, Cr, Mn, Fe, Co, Ni) have been studied theoretically using density functional theory. In most of the lowest energy predicted M(C(3)H(5))(n) structures all of the allyl groups are bonded as trihapto η(3)-C(3)H(5) ligands and the metals have considerably less than the normally favored 18-electron configuration. Such ligands can be considered formally as bidentate ligands with the metal atom connected to the centers of the two C-C bonds of the η(3)-C(3)H(5) group. The later transition metal diallyls M(C(3)H(5))(2) (M = Cr, Mn, Fe, Co, Ni) form two stereoisomers of similar relative energies, namely the C(2h) staggered isomer and the C(2v) eclipsed isomer with the orientation of the η(3)-C(3)H(5) groups corresponding to square planar metal coordination of the bidentate η(3)-C(3)H(5) ligands. The staggered and eclipsed Ni(C(3)H(5))(2) isomers have been observed experimentally by NMR. Less symmetrical M(C(3)H(5))(2) structures are found for the earlier transition metals Sc, Ti, and V in which the orientation of the allyl groups corresponds to tetrahedral metal coordination. The triallylmetal derivatives M(C(3)H(5))(3) are predicted to be thermodynamically viable with respect to allyl loss to give the corresponding diallylmetal derivatives, except for triallylnickel. The lowest energy Ni(C(3)H(5))(3) structure has two trihaptoallyl ligands and one monohaptoallyl ligand, whereas the lowest energy Mn(C(3)H(5))(3) structures have only one trihaptoallyl ligand and two monohaptoallyl ligands. Otherwise, the M(C(3)H(5))(3) complexes have structures with three trihaptoallyl ligands corresponding formally to octahedral metal coordination. The M(C(3)H(5))(3) complexes (M = Cr, Co) thus correspond to a well-known series of "classical" octahedral coordination complexes, namely, those of the d(3) Cr(III) and the d(6) Co(III), respectively.     

Enhanced H2 adsorption in isostructural metal-organic frameworks with open metal sites: strong dependence of the binding strength on metal ions

Metal-organic frameworks (MOFs) with open metal sites exhibit a much stronger H2 binding strength than classical MOFs, due to the direct interaction between H2 and the coordinately unsaturated metal ions. Here we report a systematic study of the H2 adsorption on a series of isostructural MOFs, M2(dhtp) (M = Mg, Mn, Co, Ni, Zn). The experimental, initial isosteric heats of adsorption for H2 (Qst) of these MOFs range from 8.5 to 12.9 kJ/mol, with increasing Qst in the following order: Zn, Mn, Mg, Co, and Ni. The H2 binding energies derived from first-principles calculation follow the same trend as the experimental observation on Qst, confirming the electrostatic Coulomb attraction between the H2 and the open metals being the major interaction. We also found a strong correlation between the metal ion radius, the M-H2 distance, and the H2 binding strength, which provides a viable, empirical method to predict the relative H2 binding strength of different open metals.     

Interaction of molecular hydrogen with open transition metal centers for enhanced binding in metal-organic frameworks: a computational study

Molecular hydrogen is known to form stable, "nonclassical" sigma complexes with transition metal centers that are stabilized by donor-acceptor interactions and electrostatics. In this computational study, we establish that strong H2 sorption sites can be obtained in metal-organic frameworks by incorporating open transition metal sites on the organic linkers. Using density functional theory and energy decomposition analysis, we investigate the nature and characteristics of the H2 interaction with models of exposed open metal binding sites {half-sandwich piano-stool shaped complexes of the form (Arene)ML(3- n)(H2)n [M=Cr, Mo, V(-), Mn(+); Arene = C6H5X (X=H, F, Cl, OCH3, NH2, CH3, CF3) or C6H3Y2X (Y=COOH, X=CF3, Cl; L=CO; n=1-3]}. The metal-H2 bond dissociation energy of the studied complexes is calculated to be between 48 and 84 kJ/mol, based on the introduction of arene substituents, changes to the metal core, and of charge-balancing ligands. Thus, design of the binding site controls the H2 binding affinity and could be potentially used to control the magnitude of the H2 interaction energy to achieve reversible sorption characteristics at ambient conditions. Energy decomposition analysis illuminates both the possibilities and present challenges associated with rational materials design.     

Understanding the Chloride Affinity of Barbiturates for Anion Receptor Design

Due to their potential binding sites, barbituric acid (BA) and its derivatives have been used in metal coordination chemistry. Yet their abilities to recognize anions remain unexplored. In this work, we were able to identify four structural features of barbiturates that are responsible for a certain anion affinity. The set of coordination interactions can be finely tuned with covalent decorations at the methylene group. DFT-D computations at the BLYP-D3(BJ)/aug-cc-pVDZ level of theory show that the C−H bond is as effective as the N−H bond to coordinate chloride. An analysis of the electron charge density at the C−H⋅⋅⋅Cl⁻ and N−H⋅⋅⋅Cl⁻ bond critical points elucidates their similarities in covalent character. Our results reveal that the special acidity of the C−H bond shows up when the methylene group moves out of the ring plane and it is mainly governed by the orbital interaction energy. The amide and carboxyl groups are the best choices to coordinate the ion when they act together with the C−H bond. We finally show how can we use this information to rationally improve the recognition capability of a small cage-like complex that is able to coordinate NaCl.     

Origins of dihydrogen binding to metal-inserted porphyrins: electric polarization and Kubas interaction

Using density functional theory calculations, we have investigated the interactions between hydrogen molecules and metalloporphyrins. A metal atom, such as Ca or Ti, is introduced for incorporation in the central N(4) cavity. Within local density approximation (generalized gradient approximation), we find that the average binding energy of H(2) to the Ca atom is about 0.25 (0.1) eV/H(2) up to four H(2) molecules, whereas that to the Ti atom is about 0.6 (0.3) eV per H(2) up to two H(2) molecules. Our analysis of orbital hybridization between the inserted metal atom and molecular hydrogen shows that H(2) binds weakly to Ca-porphyrin through a weak electric polarization in dihydrogen, but is strongly hybridized with Ti-porphyrin through the Kubas interaction. The presence of d orbitals in Ti may explain the difference in the interaction types.     

Electronic structure of [(CpCr)[(CO)3M]]mu-Cot (M = Cr, Fe): a theoretical study

The electronic structure of two cyclooctatetraene-bridged dinuclear first-row transition metal complexes of the type [(CpM)[(CO)3M']]mu-Cot (M = Cr; M' = Fe (1), Cr (2)) was investigated by complete active space self-consistent field (CASSCF) calculations. In this context the differences in the binding capabilities of the complex fragments CpM and (CO)3M are discussed on the basis of extended Huckel molecular orbital (MO) calculations. The geometries used for the CASSCF calculations for complex 1 were obtained from the crystal structure. For 2 a model structure was established by geometry optimization using density functional methods. The CASSCF results agree well with the experimental findings and provide insight into the binding situation of the two compounds. Complex 1 can be regarded as being composed of a chromocene-like subunit CpCr(eta5-C5H5) and the fragment (CO)3Fe(eta3-C3H3). A direct metal-metal bond is found, involving one initially singly occupied orbital of each fragment, leading to a doublet ground state for 1 with the remaining unpaired electron localized at the chromium center. For 2 no such direct metal-metal bond can be recognized. A very weak direct metal-metal interaction is induced by electron donation from the Cot2- ligand into a formally unoccupied metal-metal binding orbital combination. In the quartet ground state all three unpaired electrons are localized at the chromium center of the formally doubly positive charged CpCr unit, on which complex fragment [(CO)3Cr(eta5-Cot)]2- acts like a cyclopentadienyl ligand. The coordination sphere of the chromium center of the CpCr unit resembles that of a metallocene metal center and its metal 3d occupation scheme corresponds to that of vanadocene.     

First-principles study of hydrogen adsorption in metal-doped COF-10

Covalent organic frameworks (COFs), due to their low-density, high-porosity, and high-stability, have promising applications in gas storage. In this study we have explored the potential of COFs doped with Li and Ca metal atoms for storing hydrogen under ambient thermodynamic conditions. Using density functional theory we have performed detailed calculations of the sites Li and Ca atoms occupy in COF-10 and their interaction with hydrogen molecules. The binding energy of Li atom on COF-10 substrate is found to be about 1.0 eV and each Li atom can adsorb up to three H(2) molecules. However, at high concentration, Li atoms cluster and, consequently, their hydrogen storage capacity is reduced due to steric hindrance between H(2) molecules. On the other hand, due to charge transfer from Li to the substrate, O sites provide additional enhancement for hydrogen adsorption. With increasing concentration of doped metal atoms, the COF-10 substrate provides an additional platform for storing hydrogen. Similar conclusions are reached for Ca doped COF-10.     

Investigating bidentate and tridentate carbamoylmethylphosphine oxide ligand interactions with rare-Earth elements using electrospray ionization quadrupole ion trap mass spectrometry

Electrospray ionization (ESI) quadrupole ion trap mass spectrometry (QIT-MS) and collisionally activated dissociation (CAD) were used to evaluate the rare-earth binding properties of two hydrophobic carbamoylmethylphosphine oxide (CMPO) ligands, the normal bidentate variety, (t-BuC6H4)2P(O)CH2C(O)N(i-Bu)2 (A), a new potentially tridentate extractant, (t-BuC6H4)2P(O)CH[CH2C(O)N(i-Bu)2]C(O)N(i-Bu)2 (B), and tributyl phosphate. The mass spectral results obtained from analysis of 1% HNO3/methanol solution containing the ligands and dissolved lanthanide salts reveal that the favorable stoichiometries of the ligand/metal/nitrate complexes are 2:1:2 for the bidentate ligand A, 1:1:2 for the tridentate ligand B, and 3:1:2 for the monodentate tributyl phosphate. These observed stoichiometries correlate with the number of available binding sites on each ligand as well as with potential steric effects. Energy-variable collisionally activated dissociation experiments showed that for the 2:1:2 complexes involving ligand A or B, as the ionic radius of the bound metal decreased, the removal of nitric acid required less energy and resulted in less extensive spontaneous solvent coordination. This experimental trend suggests that, as the ionic radius of the lanthanide ion decreases, a pair of the carbamoylmethylphosphine ligands is able to more completely solvate the bound metal ion thereby weakening the nitrate-metal interaction.     

Conformational/configurational analysis of all the binding geometries of cobalt(III) bleomycin

No crystal structure of metallobleomycin (BLM) exists, and the exact coordination mode of the ligand is unknown. To date, spectroscopic investigations of BLM complexes and crystal structures of BLM models have been used to propose its metal coordination sites. This has led to contradictory interpretations of the metal coordination sphere in BLM. Inorganic molecular mechanics and configurational/conformational searches were used to analyze HOO-CoBLM A2, H2O-CoBLM A2, and HOO-CoPEP with commonly proposed binding geometries. The lowest energy binding geometry found was one with the mannose carbamoyl bound to the cobalt ion. The Monte Carlo dihedral and translational variational searches were able to find most of the configurations available to cobalt(III) bleomycin in the three binding geometries examined.     

H2 adsorption on 3d transition metal clusters: a combined infrared spectroscopy and density functional study

The adsorption of H2 on a series of gas-phase transition metal (scandium, vanadium, iron, cobalt, and nickel) clusters containing up to 20 metal atoms is studied using IR-multiple photon dissociation spectroscopy complemented with density functional theory based calculations. Comparison of the experimental and calculated spectra gives information on hydrogen-bonding geometries. The adsorption of H2 is found to be exclusively dissociative on Sc(n)O+, V(n)+, Fe(n)+, and Co(n)+, and both atomic and molecularly chemisorbed hydrogen is present in Ni(n)H(m)+ complexes. It is shown that hydrogen adsorption geometries depend on the elemental composition as well as on the cluster size and that the adsorption sites are different for clusters and extended surfaces. In contrast to what is observed for extended metal surfaces, where hydrogen has a preference for high coordination sites, hydrogen can be both 2- or 3-fold coordinated to cationic metal clusters.     

Titanium-capped carbon chains as promising new hydrogen storage media

The capacity of Ti-capped sp carbon atomic chains for use as hydrogen storage media is studied using first-principles density functional theory. The Ti atom is strongly attached at one end of the carbon chains via d-p hybridization, forming stable TiC(n) complexes. We demonstrate that the number of adsorbed H(2) molecules on Ti through Kubas interactions depends upon the chain types. For polyyne (n even) or cumulene (n odd) structures, each Ti atom can hold up to five or six H(2) molecules, respectively. Furthermore, the TiC(5) chain effectively terminated on a C(20) fullerene can store hydrogen with an optimal binding energy of 0.52 eV per H(2) molecule. Our results reveal a possible way to explore high-capacity hydrogen storage materials in truly one-dimensional carbon structures.     

Homo‐Helical Rod Packing as a Path Toward the Highest Density of Guest‐Binding Metal Sites in Metal–Organic Frameworks

In porous materials, metal sites with coordinate solvents offer opportunities for many applications, especially those promoted by host–guest chemistry, but such sites are especially hard to create for Li‐based materials, because unlike transition metals, lithium does not usually possess a high‐enough coordination number for both framework construction and guest binding. This challenge is addressed by mimicking the functional group ratio and metal‐to‐ligand charge ratio in MOF‐74. A family of rod‐packing lithium–organic frameworks (CPM‐47, CPM‐48, and CPM‐49) were obtained. These materials exhibit an extremely high density of guest‐binding lithium sites. Also unusual is the homo‐helical rod‐packing in the CPM series, as compared to the hetero‐helical rod packing by helices of opposite handedness in MOF‐74. This work demonstrates new chemical and structural possibilities in developing a record‐setting high density of guest‐binding metal sites in inorganic–organic porous materials.     

Density Functional Theory Study of Hydrogen Adsorption in a Ti‐Decorated Mg‐Based Metal–Organic Framework‐74

The Ti‐binding energy and hydrogen adsorption energy of a Ti‐decorated Mg‐based metal–organic framework‐74 (Mg‐MOF‐74) were evaluated by using first‐principles calculations. Our results revealed that only three Ti adsorption sites were found to be stable. The adsorption site near the metal oxide unit is the most stable. To investigate the hydrogen‐adsorption properties of Ti‐functionalized Mg‐MOF‐74, the hydrogen‐binding energy was determined. For the most stable Ti adsorption site, we found that the hydrogen adsorption energy ranged from 0.26 to 0.48 eV H₂^?1. This is within the desirable range for practical hydrogen‐storage applications. Moreover, the hydrogen capacity was determined by using ab initio molecular dynamics simulations. Our results revealed that the hydrogen uptake by Ti‐decorated Mg‐MOF‐74 at temperatures of 77, 150, and 298 K and ambient pressure were 1.81, 1.74, and 1.29 H₂ wt %, respectively.     

Ligand-assisted enhancement of CO2 capture in metal-organic frameworks

Using density functional theory with a van der Waals-corrected functional, we elucidate how CO(2) binds to a novel "BTT-type" metal-organic framework (MOF) featuring open metal centers. We show that CO(2) binds most favorably to open metal cation sites, but with an adsorption energy that can be three times more sensitive to the choice of the bridging ligand than to metal cation choice. A strong, three-site interaction between CO(2) and the open-metal site is predicted, with the binding energy enhanced by up to a factor of 2, depending on the ligand. The CO(2)-MOF binding can be attributed to a combination of electrostatics and vdW dispersive interactions, both of which are critically sensitive to the local environment, and both of which contribute nearly equally to the overall binding strength. We show that a judicious choice of the organic linker and the metal center allows the binding energy to be tuned from 34.8 kJ/mol (for CaBTTri) to a maximum of 64.5 kJ/mol (MgBTT).     

Kubas-type hydrogen storage in V(III) polymers using tri- and tetradentate bridging ligands

Oxalic acid, oxamide, glycolic acid, and glycolamide were employed as 2-carbon linkers to synthesize a series of one-dimensional V(III) polymers from trismesityl vanadium(III)·THF containing a high concentration of low-valent metal sites that can be exploited for Kubas binding in hydrogen storage. Synthesized materials were characterized by powder X-ray diffraction (PXRD), nitrogen adsorption (BET), X-ray photoelectron spectroscopy (XPS), infrared spectroscopy (IR), Raman spectroscopy, thermogravimetric analysis, and elemental analysis. Because each of these organic linkers possesses a different number of protons and coordinating atoms, the products in each case were expected to have different stoichiometries with respect to the number of mesityl groups eliminated and also a different geometry about the V(III) centers. For example, the oxalate and glycolate polymers contained residual mesityl groups; however, these could be exchanged with hydride via hydrogenolysis. The highest adsorption capacity was recorded on the product of trismesityl vanadium(III)·THF with oxamide (3.49 wt % at 77 K and 85 bar). As suggested by the high enthalpy of adsorption (17.9 kJ/mol H(2)), a substantial degree of performance of the vanadium metal centers was retained at room temperature (25%), corresponding to a gravimetric adsorption of 0.87 wt % at 85 bar, close to the performance of MOF-177 at this temperature and pressure. This is remarkable given the BET surface area of this material is only 9 m(2)/g. A calculation on the basis of thermogravimetric results provides 0.88 hydrogen molecule per vanadium center under these conditions. Raman studies with H(2) and D(2) showed the first unequivocal evidence for Kubas binding on a framework metal in an extended solid, and IR studies demonstrated H(D) exchange of the vanadium hydride with coordinated D(2). These spectroscopic observations are sufficient to assign the rising trends in isosteric heats of hydrogen adsorption observed previously by our group in several classes of materials containing low-valent transition metals to the Kubas interaction.     

金属串配合物(n,m)[Cr3(PhPyF)4Cl2](n=2, 3, 4; m=2, 1, 0)的配位结构及其与电场的关系

应用密度泛函理论BP86 方法研究具有分子导线潜在应用的金属串配合物(n, m)[Cr₃(PhPyF)₄Cl₂](HPhPyF=N, N'-苯基吡啶基甲脒; n=2, 3, 4; m=2, 1, 0)的配位结构及其受电场作用的影响, n、m分别表示PhPyF-的苯环在左侧和在右侧的配体个数. 结果表明: (1) 零电场下, 四个PhPyF-的(2, 2)、(3, 1)和(4, 0)三种配位方式能量差别很小, 为竞争态, (2, 2)最稳定. (4, 0)结构中两端轴向配体Cl 均可与Cr 配位, 且Cl₄―Cr1 键比Cl5―Cr3键更强, 若作为分子器件可与电极结合, 这与(4, 0)[CuCuM(npa)₄Cl][PF₆](M=Pd, Pt; Hnpa=2-萘啶苯胺)靠近苯环一端的轴向配体无法与M配位不同. (2) 在(2, 2)、(3, 1)和(4, 0)中, Cr₃⁶⁺链均具有三中心三电子离域σ键, 但离域性逐渐减弱. 随四个PhPyF-配位方式趋于一致, 分子极性逐渐增大, 由Cl₄指向Cl5(Z)方向, Cr1的α自旋密度增大, Cr2 的β和Cr3 的α自旋密度减小. (3) 分子的几何结构和电子结构在电场下发生规律性变化, 在-Z方向电场作用下, (3, 1)、(4, 0)电子移动方向与极性方向相同, 使分子的键长、自旋密度、电荷和能隙变化显著性均大于Z方向电场, 且极性越大变化越显著, 有利于提高分子导电性.     

Thermodynamics of pyridine coordination in 1,4-phenylene bridged bimetallic (Pd, Pt) complexes containing two N,C,N' motifs, 1,4-M2-[C6(CH2NR2)4-2,3,5,6

The thermodynamics of pyridine coordination in 1,4-phenylene-bridged binuclear palladium and platinum organometallic complexes [1,4-(MOTf)2-&{C6(CH2NR2)4-2,3,5,6}] (11, M =Pd, Pt; R =CH3, C2H5, R2 = -(CH2)5-) are measured by 1H NMR in DMSO-d6. The coordination of substituted pyridines by bimetallic complexes 11 or 12 in DMSO is found to proceed via two effectively independent metalligand binding events, and the association constants for pyridine coordination and rate constants for pyridine exchange are nearly identical to those measured previously on monometallic analogs. A linear free energy relationship between the association constant for pyridine coordination and the inductive Hammett constant of the pyridine substituent is observed, and the sensitivity (rho = -1.7 to -2.1) in DMSO depends only slightly on metal (Pd vs Pt) and spectator ligand (pincer dialkylamine vs triarylphosphine). The association constant for a particular pyridine ligand, however, varies by roughly 3 orders of magnitude across the series of metal complexes. The effective independence of the two coordination sites and the range of available thermodynamic and kinetic behaviors of the coordination guide the use of these versatile building blocks in metallosupramolecular applications.     

First-second shell interactions in metal binding sites in proteins: a PDB survey and DFT/CDM calculations

The role of the second shell in the process of metal binding and selectivity in metalloproteins has been elucidated by combining Protein Data Bank (PDB) surveys of Mg, Mn, Ca, and Zn binding sites with density functional theory/continuum dielectric methods (DFT/CDM). Peptide backbone groups were found to be the most common second-shell ligand in Mg, Mn, Ca, and Zn binding sites, followed (in decreasing order) by Asp/Glu, Lys/Arg, Asn/Gln, and Ser/Thr side chains. Aromatic oxygen- or nitrogen-containing side chains (Tyr, His, and Trp) and sulfur-containing side chains (Cys and Met) are seldom found in the second coordination layer. The backbone and Asn/Gln side chain are ubiquitous in the metal second coordination layer as their carbonyl oxygen and amide hydrogen can act as a hydrogen-bond acceptor and donor, respectively, and can therefore partner practically every first-shell ligand. The second most common outer-shell ligand, Asp/Glu, predominantly hydrogen bonds to a metal-bound water or Zn-bound histidine and polarizes the H-O or H-N bond. In certain cases, a second-shell Asp/Glu could affect the protonation state of the metal ligand. It could also energetically stabilize a positively charged metal complex more than a neutral ligand such as the backbone and Asn/Gln side chain. As for the first shell, the second shell is predicted to contribute to the metal selectivity of the binding site by discriminating between metal cations of different ionic radii and coordination geometries. The first-shell-second-shell interaction energies decay rapidly with increasing solvent exposure of the metal binding site. They are less favorable but are of the same order of magnitude as compared to the respective metal-first-shell interaction energies. Altogether, the results indicate that the structure and properties of the second shell are dictated by those of the first layer. The outer shell is apparently designed to stabilize/protect the inner-shell and complement/enhance its properties.