Diffusion network of CO in FeFe-Hydrogenase

Second and Outer Coordination Sphere Effects in Nitrogenase, Hydrogenase, Formate Dehydrogenase, and CO Dehydrogenase

Gases like H2, N2, CO2, and CO are increasingly recognized as critical feedstock in "green" energy conversion and as sources of nitrogen and carbon for the agricultural and chemical sectors. However, the industrial transformation of N2, CO2, and CO and the production of H2 require significant energy input, which renders processes like steam reforming and the Haber-Bosch reaction economically and environmentally unviable. Nature, on the other hand, performs similar tasks efficiently at ambient temperature and pressure, exploiting gas-processing metalloenzymes (GPMs) that bind low-valent metal cofactors based on iron, nickel, molybdenum, tungsten, and sulfur. Such systems are studied to understand the biocatalytic principles of gas conversion including N2 fixation by nitrogenase and H2 production by hydrogenase as well as CO2 and CO conversion by formate dehydrogenase, carbon monoxide dehydrogenase, and nitrogenase. In this review, we emphasize the importance of the cofactor/protein interface, discussing how second and outer coordination sphere effects determine, modulate, and optimize the catalytic activity of GPMs. These may comprise ionic interactions in the second coordination sphere that shape the electron density distribution across the cofactor, hydrogen bonding changes, and allosteric effects. In the outer coordination sphere, proton transfer and electron transfer are discussed, alongside the role of hydrophobic substrate channels and protein structural changes. Combining the information gained from structural biology, enzyme kinetics, and various spectroscopic techniques, we aim toward a comprehensive understanding of catalysis beyond the first coordination sphere.

Water Dynamics in a Peptide-appended Pillar[5]arene Artificial Channel in Lipid and Biomimetic Membranes

Peptide-appended Pillar[5]arene (PAP) is an artificial water channel that can be incorporated into lipid and polymeric membranes to achieve high permeability and enhanced selectivity for angstrom-scale separations [Shen et al. Nat. Commun.9:2294 (2018)]. In comparison to commonly studied rigid carbon nanotubes, PAP channels are conformationally flexible, yet these channels allow a high water permeability [Y. Liu and H. Vashisth Phys. Chem. Chem. Phys.21:22711 (2019)]. Using molecular dynamics (MD) simulations, we study water dynamics in PAP channels embedded in biological (lipid) and biomimetic (block-copolymer) membranes to probe the effect of the membrane environment on water transport characteristics of PAP channels. We have resolved the free energy surface and local minima for water diffusion within the channel in each type of membrane. We find that water follows single file transport with low free-energy barriers in regions surroundings the central ring of the PAP channel and the single file diffusivity of water correlates with the number of hydrogen bonding sites within the channel, as is known for other sub-nm pore-size synthetic and biological water channels [Horner et al. Sci. Adv.1:e1400083 (2015)].

Molecular dynamics simulations of ortho-carborane nano-diamond storage within the nonpolar channel cavities of a right-handed coiled-coil tetrabrachion nanotube

Molecular dynamics simulations have been performed on a complex in which clusters of boron in the form of molecules of the nanodiamond ortho-carborane ( C 2 B 10 H 12 ) have been inserted into the four large nonpolar cavities of a nanotube of the right-handed coiled-coil ( R H C C ) t e t r a b r a c h i o n . The techniques of multi-configurational thermodynamic integration, steered molecular dynamics and umbrella sampling have been combined to investigate the energetics of storage of ortho-carborane in the cavities and to map out the free energy landscape of the RHCC - t e t r a b r a c h i o n - o r t h o - c a r b o r a n e complex along the central channel and along directions transverse to the central channel. The purpose of the study was to explore potential pathways for the diffusion of ortho-carborane between the cavities and the solvent and to assess the stability of the complex as a possible drug delivery system for boron neutron capture therapy (BNCT). The investigation reveals a complex free energy landscape with a multitude of peaks and valleys, all of which can be related to specific architectural elements of the RHCC - n a n o t u b e , and the activation barriers for ortho-carborane capture and release support the requirements for rapid cargo uptake coupled with tight binding to the cavities.

Energetics of Storage and Diffusion of Water and Cyclo-Octasulfur for a Nonpolar Cavity of RHCC Tetrabrachion by Molecular Dynamics Simulations

Tetrabrachion forms the key component of the S-layer of Staphylothermus marinus. Molecular dynamics simulations have been used to study the energetics of occupancy of cavity 3 of the right-handed coiled-coil stalk of tetrabrachion by both water molecules and cyclooctasulfur S₈ crowns, as well as to determine possible pathways and free energy barriers for the diffusion of both water and cyclooctasulfur through the peptide walls of RHCC tetrabrachion between cavity 3 and bulk solvent. Calculations of the transfer free energy from solvent to cavity show that clusters of six, seven and eight water molecules are marginally stable in cavity 3, but that occupancy of the cavity by a cyclooctasulfur ring is favoured significantly over water clusters of all sizes. Thermal activation simulations at T = 400K revealed that water molecules diffusing through the wall pass through a sequence of metastable configurations where they are temporarily immobilized by forming networks of hydrogen bonds with specific wall residues. Calculations of the free energy of these metastable configurations using multi-configurational thermodynamic integration yielded a free energy profile with a principal free energy maximum ∆G~50 kJ/mol and a slight activation asymmetry in favour of the direction from cavity to solvent. Potential exit pathways for cyclooctasulfur were investigated with the methods of steered molecular dynamics and umbrella sampling. The cyclooctasulfur was steered through a gap in the tetrabrachion wall along a linear path from cavity 3 into the solvent and the resulting trajectory was subdivided into 22 sampling windows. The free energy profile created for the trajectory by umbrella sampling showed a sharp principal maximum as a function of the reaction coordinate with asymmetric free energy barriers ∆G_exit~220 kJ/mol and ∆G_entrance~100 kJ/mol for cavity exit and entrance, respectively.