Fe as Hydrogen/Halogen Bond Acceptor in Square Pyramidal Fe(CO)(5)

Hydrogen bonded complexes formed between the square pyramidal Fe(CO)(5) with HX (X = F, Cl, Br), showing X-H center dot center dot center dot Fe interactions, have been investigated theoretically using density functional theory (DFT) including dispersion correction. Geometry, interaction energy, and large red shift of about 400 cm(-1) in the FIX stretching frequency confirm X-H center dot center dot center dot Fe hydrogen bond formation. In the (CO)(5)Fe center dot center dot center dot HBr complex, following the significant red shift, the HBr stretching mode is coupled with the carbonyl stretching modes. This clearly affects the correlation between frequency shift and binding energy, which is a hallmark of hydrogen bonds. Atoms in Molecule (AIM) theoretical analyses show the presence of a bond critical point between the iron and the hydrogen of FIX and significant mutual penetration. These X-H center dot center dot center dot Fe hydrogen bonds follow most but not all of the eight criteria proposed by Koch and Popelier (J. Phys. Chem. 1995, 99, 9747) based on their investigations on C-H center dot center dot center dot O hydrogen bonds. Natural bond orbital (NBO) analysis indicates charge transfer from the organometallic system to the hydrogen bond donor. However, there is no correlation between the extent of charge transfer and interaction,energy, contrary to what is proposed in the recent IUPAC recommendation (Pure Appl.. Chem. 2011, 83, 1637). The ``hydrogen bond radius'' for iron has been determined to be 1.60 +/- 0.02 angstrom, and not surprisingly it is between the covalent (127 angstrom) and van der Waals (2.0) radii of Fe. DFT and AIM theoretical studies reveal that Fe in square pyramidal Fe(CO)(5) can also form halogen bond with CIF and ClH as ``halogen bond donor''. Both these complexes show mutual penetration as well, though the Fe center dot center dot center dot Cl distance is closer to the sum of van der Waals radii of Fe and Cl in (CO)5Fe center dot center dot center dot ClH, and it is about 1 angstrom less in (CO)(5)Fe center dot center dot center dot ClF.

Quantum chemical study on influence of intermolecular hydrogen bonding on the geometry, the atomic charges and the vibrational dynamics of 2,6-dichlorobenzonitrile

FT-IR (4000-400 cm(-1)) and FT-Raman (4000-200 cm(-1)) spectral measurements on solid 2,6-dichlorobenzonitrile (2,6-DCBN) have been done. The molecular geometry, harmonic vibrational frequencies and bonding features in the ground state have been calculated by density functional theory at the B3LYP/6-311++G (d,p) level. A comparison between the calculated and the experimental results covering the molecular structure has been made. The assignments of the fundamental vibrational modes have been done on the basis of the potential energy distribution (PED). To investigate the influence of intermolecular hydrogen bonding on the geometry, the charge distribution and the vibrational spectrum of 2,6-DCBN; calculations have been done for the monomer as well as the tetramer. The intermolecular interaction energies corrected for basis set superposition error (BSSE) have been calculated using counterpoise method. Based on these results, the correlations between the vibrational modes and the structure of the tetramer have been discussed. Molecular electrostatic potential (MEP) contour map has been plotted in order to predict how different geometries could interact. The Natural Bond Orbital (NBO) analysis has been done for the chemical interpretation of hyperconjugative interactions and electron density transfer between occupied (bonding or lone pair) orbitals to unoccupied (antibonding or Rydberg) orbitals. UV spectrum was measured in methanol solution. The energies and oscillator strengths were calculated by Time Dependent Density Functional Theory (TD-DFT) and matched to the experimental findings. TD-DFT method has also been used for theoretically studying the hydrogen bonding dynamics by monitoring the spectral shifts of some characteristic vibrational modes involved in the formation of hydrogen bonds in the ground and the first excited state. The C-13 nuclear magnetic resonance (NMR) chemical shifts of the molecule were calculated by the Gauge independent atomic orbital (GIAO) method and compared with experimental results. Standard thermodynamic functions have been obtained and changes in thermodynamic properties on going from monomer to tetramer have been presented. (C) 2013 Elsevier B.V. All rights reserved.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

MolBridge: a program for identifying nonbonded interactions in small molecules and biomolecular structures

Identification and analysis of nonbonded interactions within a molecule and with the surrounding molecules are an essential part of structural studies, given the importance of these interactions in defining the structure and function of any supramolecular entity. MolBridge is an easy to use algorithm based purely on geometric criteria that can identify all possible nonbonded interactions, such as hydrogen bond, halogen bond, cation-pi, pi-pi and van der Waals, in small molecules as well as biomolecules. The user can either upload three-dimensional coordinate files or enter the molecular ID corresponding to the relevant database. The program is available in a standalone form and as an interactive web server with Jmol and JME incorporated into it. The program is freely downloadable and the web server version is also available at http://nucleix.mbu.iisc.ernet.in/molbridge/index.php.

Rotational spectra of propargyl alcohol dimer: A dimer bound with three different types of hydrogen bonds

Pure rotational spectra of the propargyl alcohol dimer and its three deuterium isotopologues have been observed in the 4 to 13 GHz range using a pulsed-nozzle Fourier transform microwave spectrometer. For the parent dimer, a total of 51 transitions could be observed and fitted within experimental uncertainty. For two mono-substituted and one bi-substituted deuterium isotopologues, a total of 14, 17, and 19 transitions were observed, respectively. The observed rotational constants for the parent dimer A = 2321.8335(4) MHz, B = 1150.4774(2) MHz, and C = 1124.8898(2) MHz] are close to those of the most stable structure predicted by ab initio calculations. Spectra of the three deuterated isotopologues and Kraitchman analysis positively confirm this structure. Geometrical parameters and ``Atoms in Molecules'' analysis on the observed structure reveal that the two propargyl alcohol units in the dimer are bound by three different types of hydrogen bonds: O-H center dot center dot center dot O, O-H center dot center dot center dot pi, and C-H center dot center dot center dot pi. To the best of our knowledge, propargyl alcohol seems to be the smallest molecule forming a homodimer with three different points of contact. (C) 2014 AIP Publishing LLC.

Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Dynamics of the chemical bond: inter- and intra-molecular hydrogen bond

In this discussion, we show that a static definition of a `bond' is not viable by looking at a few examples for both inter-and intra-molecular hydrogen bonding. This follows from our earlier work (Goswami and Arunan, Phys. Chem. Chem. Phys. 2009, 11, 8974) which showed a practical way to differentiate `hydrogen bonding' from `van der Waals interaction'. We report results from ab initio and atoms in molecules theoretical calculations for a series of Rg center dot center dot center dot HX complexes (Rg = He/Ne/Ar and X = F/Cl/Br) and ethane-1,2-diol. Results for the Rg center dot center dot center dot HX/DX complexes show that Rg center dot center dot center dot DX could have a `deuterium bond' even when Rg center dot center dot center dot HX is not `hydrogen bonded', according to the practical criterion given by Goswami and Arunan. Results for ethane-1,2-diol show that an `intra-molecular hydrogen bond' can appear during a normal mode vibration which is dominated by the O center dot center dot center dot O stretching, though a `bond' is not found in the equilibrium structure. This dynamical `bond' formation may nevertheless be important in ensuring the continuity of electron density across a molecule. In the former case, a vibration `breaks' an existing bond and in the later case, a vibration leads to `bond' formation. In both cases, the molecule/complex stays bound irrespective of what happens to this `hydrogen bond'. Both these cases push the borders on the recent IUPAC recommendation on hydrogen bonding (Arunan et al. Pure. Appl. Chem. 2011, 83 1637) and justify the inclusive nature of the definition.

Intramolecular hydrogen bonds involving organic fluorine in the derivatives of hydrazides: an NMR investigation substantiated by DFT based theoretical calculations

The rare examples of intramolecular hydrogen bonds (HB) of the type the N-H center dot center dot center dot F-C, detected in a low polarity solvent in the derivatives of hydrazides, by utilizing one and two-dimensional solution state multinuclear NMR techniques, are reported. The observation of through-space couplings, such as, (1h)J(FH), and (1h)J(FN), provides direct evidence for the existence of intra-molecular HB. Solvent induced perturbations and the variable temperature NMR experiments unambiguously establish the presence of intramolecular HB. The existence of multiple conformers in some of the investigated molecules is also revealed by two dimensional HOESY and N-15-H-1 HSQC experiments. The H-1 DOSY experimental results discard any possibility of self or cross dimerization of the molecules. The derived NMR experimental results are further substantiated by Density Function Theory (DFT) based Non Covalent Interaction (NCI), and Quantum Theory of Atom in Molecule (QTAIM) calculations. The NCI calculations served as a very sensitive tool for detection of non-covalent interactions and also confirm the presence of bifurcated HBs.

Isotopic fractionation in proteins as a measure of hydrogen bond length

If a deuterated molecule containing strong intramolecular hydrogen bonds is placed in a hydrogenated solvent, it may preferentially exchange deuterium for hydrogen. This preference is due to the difference between the vibrational zero-point energy for hydrogen and deuterium. It is found that the associated fractionation factor (I) is correlated with the strength of the intramolecular hydrogen bonds. This correlation has been used to determine the length of the H-bonds (donor-acceptor separation) in a diverse range of enzymes and has been argued to support the existence of short low-barrier H-bonds. Starting with a potential energy surface based on a simple diabatic state model for H-bonds, we calculate (I) as a function of the proton donor-acceptor distance R. For numerical results, we use a parameterization of the model for symmetric 0-H. ``.0 bonds R. H. McKenzie, Chem. Phys. Lett. 535, 196 (2012)]. We consider the relative contributions of the 0-H stretch vibration, O-H bend vibrations (both in plane and out of plane), tunneling splitting effects at finite temperature, and the secondary geometric isotope effect. We compare our total (I) as a function of R with NMR experimental results for enzymes, and in particular with an earlier model parametrization (D(R), used previously to determine bond lengths. (C) 2015 AIP Publishing LLC.

Studies of Organic Molecule Recognition by Synthetic Cyclophane Receptors in Aqueous Media

The behaviors of six new cyclophane receptors for organic guest molecules in aqueous media are reported. These new hosts are modifications of more basic parent structures, and the main goal of their examination has been to determine how the modifications affect host selectivity for cationic guests. In particular, we have been interested in determining how additional non-covalent binding interactions can complement the cation-π interactions active in the parent systems. Three types of modifications were made to these systems. Firstly, neutral methoxy and bromine substituents were added to produce four of the six new macrocycles. Secondly, two additional aromatic rings (relative to the parent host) capable of making cation-π interactions with charged guest species were appended. Thirdly, a negatively charged carboxyl group was attached to produce a cavity in which electrostatic interactions should enhance cationic guest binding. ^1H-NMR and circular dichroic techniques were employed to determine the binding affinities of a wide variety of organic guests for the parent and modified structures in aqueous media.

Bromination of the parent host greatly enhances its binding in a general fashion, primarily as the result of hydrophobic interactions. The addition of methoxy groups does not enhance binding, apparently as a result of a collapse of the hosts into a conformation that is not suitable for binding. The appendage of extra aromatic rings enhances the binding of positively charged guests, most likely in response to more complete encapsulation of guest species. The addition of a negatively charged carboxylate enhances the binding to only selective groups of cationic guests. AM1 calculations of the electrostatic potentials of several guests molecules suggests that the enhancements seen with the modified receptor compared to the parent are most likely the result of close contact between regions of highest potential on the guest and the appended carboxylate.

Hydrogen-Atom Transfer Photochemistry of Tetrakis(µ-pyrophosphito)diplatinate(II)

In many senses, the hydrogen-atom transfer reactions observed with the triplet excited state of pyrophosphito-bridged platinum(II) dimers resemble the reactions of organic ketone nπ* states. The first two chapters describe our attempts to understand the reactivity differences between these two chromophores. Reactivity of the metal dimers is strongly regulated by the detailed nature of the ligands that ring the axial site, the hydrogen-abstraction center. A hydrogen-bonded network linking the ligands facilitates H-atom transfer quenching with alcohols through the formation of a hydrogen-bonded complex between the alcohol and a dimer. For substrates of equal C-H bond strength that lack a hydroxyl group (e.g., benzyl hydrocarbons), the quenching rate is several orders of magnitude slower.

The shape and size of the axial site, as determined by the ligands, also discriminate among quenchers by their steric characteristics. Very small quenchers quench slowly because of high entropies of activation, while very large ones have large enthalpic barriers. The two effects find a balance with quenchers of "just the right size."

The third chapter discusses the design of a mass spectrometer that uses positron annihilation to ionize neutral molecules. The mass spectrometer creates positron-molecule adducts whose annihilation produces fragmentation products that may yield information on the bonding of positrons in such complexes.

On Ultrafast Time-Domain TeraHertz Spectroscopy in the Condensed Phase: Linear Spectroscopic Measurements of Hydrogen-Bond Dynamics of Astrochemical Ice Analogs and Nonlinear TeraHertz Kerr Effect Measurements of Vibrational Quantum Beats

Much of the chemistry that affects life on planet Earth occurs in the condensed phase. The TeraHertz (THz) or far-infrared (far-IR) region of the electromagnetic spectrum (from 0.1 THz to 10 THz, 3 cm^-1 to 300 cm^-1, or 3000 μm to 30 μm) has been shown to provide unique possibilities in the study of condensed-phase processes. The goal of this work is to expand the possibilities available in the THz region and undertake new investigations of fundamental interest to chemistry. Since we are fundamentally interested in condensed-phase processes, this thesis focuses on two areas where THz spectroscopy can provide new understanding: astrochemistry and solvation science. To advance these fields, we had to develop new instrumentation that would enable the experiments necessary to answer new questions in either astrochemistry or solvation science. We first developed a new experimental setup capable of studying astrochemical ice analogs in both the TeraHertz (THz), or far-Infrared (far-IR), region (0.3 - 7.5 THz; 10 - 250 cm^-1) and the mid-IR (400 - 4000 cm^-1). The importance of astrochemical ices lies in their key role in the formation of complex organic molecules, such as amino acids and sugars in space. Thus, the instruments are capable of performing variety of spectroscopic studies that can provide especially relevant laboratory data to support astronomical observations from telescopes such as the Herschel Space Telescope, the Stratospheric Observatory for Infrared Astronomy (SOFIA), and the Atacama Large Millimeter Array (ALMA). The experimental apparatus uses a THz time-domain spectrometer, with a 1750/875 nm plasma source and a GaP detector crystal, to cover the bandwidth mentioned above with ~10 GHz (~0.3 cm^-1) resolution.

Using the above instrumentation, experimental spectra of astrochemical ice analogs of water and carbon dioxide in pure, mixed, and layered ices were collected at different temperatures under high vacuum conditions with the goal of investigating the structure of the ice. We tentatively observe a new feature in both amorphous solid water and crystalline water at 33 cm^-1 (1 THz). In addition, our studies of mixed and layered ices show how it is possible to identify the location of carbon dioxide as it segregates within the ice by observing its effect on the THz spectrum of water ice. The THz spectra of mixed and layered ices are further analyzed by fitting their spectra features to those of pure amorphous solid water and crystalline water ice to quantify the effects of temperature changes on structure. From the results of this work, it appears that THz spectroscopy is potentially well suited to study thermal transformations within the ice.

To advance the study of liquids with THz spectroscopy, we developed a new ultrafast nonlinear THz spectroscopic technique: heterodyne-detected, ultrafast THz Kerr effect (TKE) spectroscopy. We implemented a heterodyne-detection scheme into a TKE spectrometer that uses a stilbazoiumbased THz emitter, 4-N,N-dimethylamino-4-N-methyl-stilbazolium 2,4,6-trimethylbenzenesulfonate (DSTMS), and high numerical aperture optics which generates THz electric field in excess of 300 kV/cm, in the sample. This allows us to report the first measurement of quantum beats at terahertz (THz) frequencies that result from vibrational coherences initiated by the nonlinear, dipolar interaction of a broadband, high-energy, (sub)picosecond THz pulse with the sample. Our instrument improves on both the frequency coverage, and sensitivity previously reported; it also ensures a backgroundless measurement of the THz Kerr effect in pure liquids. For liquid diiodomethane, we observe a quantum beat at 3.66 THz (122 cm^-1), in exact agreement with the fundamental transition frequency of the υ4 vibration of the molecule. This result provides new insight into dipolar vs. Raman selection rules at terahertz frequencies.

To conclude we discuss future directions for the nonlinear THz spectroscopy in the Blake lab. We report the first results from an experiment using a plasma-based THz source for nonlinear spectroscopy that has the potential to enable nonlinear THz spectra with a sub-100 fs temporal resolution, and how the optics involved in the plasma mechanism can enable THz pulse shaping. Finally, we discuss how a single-shot THz detection scheme could improve the acquisition of THz data and how such a scheme could be implemented in the Blake lab. The instruments developed herein will hopefully remain a part of the groups core competencies and serve as building blocks for the next generation of THz instrumentation that pushes the frontiers of both chemistry and the scientific enterprise as a whole.

Part I. Approximate Hartree-Fock wavefunctions one-electron properties and electronic structure of the water molecule. Part II. Perturbation-variational calculation of the nuclear spin-spin isotope coupling constant in HD

Part I

Several approximate Hartree-Fock SCF wavefunctions for the ground electronic state of the water molecule have been obtained using an increasing number of multicenter s, p, and d Slater-type atomic orbitals as basis sets. The predicted charge distribution has been extensively tested at each stage by calculating the electric dipole moment, molecular quadrupole moment, diamagnetic shielding, Hellmann-Feynman forces, and electric field gradients at both the hydrogen and the oxygen nuclei. It was found that a carefully optimized minimal basis set suffices to describe the electronic charge distribution adequately except in the vicinity of the oxygen nucleus. Our calculations indicate, for example, that the correct prediction of the field gradient at this nucleus requires a more flexible linear combination of p-orbitals centered on this nucleus than that in the minimal basis set. Theoretical values for the molecular octopole moment components are also reported.

Part II

The perturbation-variational theory of R. M. Pitzer for nuclear spin-spin coupling constants is applied to the HD molecule. The zero-order molecular orbital is described in terms of a single 1s Slater-type basis function centered on each nucleus. The first-order molecular orbital is expressed in terms of these two functions plus one singular basis function each of the types e^-r/r and e^-r ln r centered on one of the nuclei. The new kinds of molecular integrals were evaluated to high accuracy using numerical and analytical means. The value of the HD spin-spin coupling constant calculated with this near-minimal set of basis functions is J_HD = +96.6 cps. This represents an improvement over the previous calculated value of +120 cps obtained without using the logarithmic basis function but is still considerably off in magnitude compared with the experimental measurement of J_HD = +43 0 ± 0.5 cps.

The Design and Synthesis of Transition Metal Complexes Supported by Non-innocent Ligand Scaffolds for Small Molecule Activation

This dissertation focuses on the incorporation of non-innocent or multifunctional moieties into different ligand scaffolds to support one or multiple metal centers in close proximity. Chapter 2 focuses on the initial efforts to synthesize hetero- or homometallic tri- or dinuclear metal carbonyl complexes supported by para-terphenyl diphosphine ligands. A series of [M₂M’(CO)₄]-type clusters (M = Ni, Pd; M’ = Fe, Co) could be accessed and used to relate the metal composition to the properties of the complexes. During these studies it was also found that non-innocent behavior was observed in dinuclear Fe complexes that result from changes in oxidation state of the cluster. These studies led to efforts to rationally incorporate central arene moieties capable managing both protons and electrons during small molecule activation.

Chapter 3 discusses the synthesis of metal complexes supported by a novel para-terphenyl diphosphine ligand containing a non-innocent 1,4-hydroquinone moiety as the central arene. A Pd⁰-hydroquinone complex was found to mediate the activation of a variety of small molecules to form the corresponding Pd⁰-quinone complexes in a formal two proton ⁄ two electron transformation. Mechanistic investigations of dioxygen activation revealed a metal-first activation process followed by subsequent proton and electron transfer from the ligand. These studies revealed the capacity of the central arene substituent to serve as a reservoir for a formal equivalent of dihydrogen, although the stability of the M-quinone compounds prevented access to the PdII-quinone oxidation state, thus hindering of small molecule transformations requiring more than two electrons per equivalent of metal complex.

Chapter 4 discusses the synthesis of metal complexes supported by a ligand containing a 3,5-substituted pyridine moiety as the linker separating the phenylene phosphine donors. Nickel and palladium complexes supported by this ligand were found to tolerate a wide variety of pyridine nitrogen-coordinated electrophiles which were found to alter central pyridine electronics, and therefore metal-pyridine π-system interactions, substantially. Furthermore, nickel complexes supported by this ligand were found to activate H-B and H-Si bonds and formally hydroborate and hydrosilylate the central pyridine ring. These systems highlight the potential use of pyridine π-system-coordinated metal complexes to reversibly store reducing equivalents within the ligand framework in a manner akin to the previously discussed 1,4-hydroquinone diphosphine ligand scaffold.

Chapter 5 departs from the phosphine-based chemistry and instead focuses on the incorporation of hydrogen bonding networks into the secondary coordination sphere of [Fe₄(μ₄-O)]-type clusters supported by various pyrazolate ligands. The aim of this project is to stabilize reactive oxygenic species, such as oxos, to study their spectroscopy and reactivity in the context of complicated multimetallic clusters. Herein is reported this synthesis and electrochemical and Mössbauer characterization of a series of chloride clusters have been synthesized using parent pyrazolate and a 3-aminophenyl substituted pyrazolate ligand. Efforts to rationally access hydroxo and oxo clusters from these chloride precursors represents ongoing work that will continue in the group.

Appendix A discusses attempts to access [Fe₃Ni]-type clusters as models of the enzymatic active site of [NiFe] carbon monoxide dehydrogenase. Efforts to construct tetranuclear clusters with an interstitial sulfide proved unsuccessful, although a (μ₃-S) ligand could be installed through non-oxidative routes into triiron clusters. While [Fe₃Ni(μ₄-O)]-type clusters could be assembled, accessing an open heterobimetallic edge site proved challenging, thus prohibiting efforts to study chemical transformations, such as hydroxide attack onto carbon monoxide or carbon dioxide coordination, relevant to the native enzyme. Appendix B discusses the attempts to synthesize models of the full H-cluster of [FeFe]-hydrogenase using a bioinorganic approach. A synthetic peptide containing three cysteine donors was successfully synthesized and found to chelate a preformed synthetic [Fe₄S₄] cluster. However, efforts to incorporate the diiron subsite model complex proved challenging as the planned thioester exchange reaction was found to non-selectively acetylate the peptide backbone, thus preventing the construction of the full six-iron cluster.