Elemental sulfur exists primarily as an S80 ring and serves as terminal electron acceptor for a variety of sulfur-fermenting bacteria. Hyperthermophilic archaea from black smoker vents are an exciting research tool to advance our knowledge of sulfur respiration under extreme conditions. Here, we use a hybrid method approach to demonstrate that the proteinaceous cavities of the S-layer nanotube of the hyperthermophilic archaeon Staphylothermus marinus act as a storage reservoir for cyclo-octasulfur S8. Fully atomistic molecular dynamics (MD) simulations were performed and the method of multiconfigurational thermodynamic integration was employed to compute the absolute free energy for transferring a ring of elemental sulfur S8 from an aqueous bath into the largest hydrophobic cavity of a fragment of archaeal tetrabrachion. Comparisons with earlier MD studies of the free energy of hydration as a function of water occupancy in the same cavity of archaeal tetrabrachion show that the sulfur ring is energetically favored over water.

Boron neutron capture therapy (BNCT) is a two-step therapeutic process that utilizes Boron-10 in combination with low energy neutrons to effectively eliminate targeted cells. This therapy is primarily used for difficult to treat head and neck carcinomas; recent advances have expanded this method to cover a broader range of carcinomas. However, it still remains an unconventional therapy where one of the barriers for widespread adoption is the adequate delivery of Boron-10 to target cells. In an effort to address this issue, we examined a unique nanoparticle drug delivery system based on a highly stable and modular proteinaceous nanotube. Initially, we confirmed and structurally analyzed ortho-carborane binding into the cavities of the nanotube. The high ratio of Boron to proteinaceous mass and excellent thermal stability suggest the nanotube system as a suitable candidate for drug delivery into cancer cells. The full physicochemical characterization of the nanotube then allowed for further mechanistic molecular dynamic studies of the ortho-carborane uptake and calculations of corresponding energy profiles. Visualization of the binding event highlighted the protein dynamics and the importance of the interhelical channel formation to allow movement of the boron cluster into the nanotube. Additionally, cell assays showed that the nanotube can penetrate outer membranes of cancer cells followed by localization around the cells’ nuclei. This work uses an integrative approach combining experimental data from structural, molecular dynamics simulations and biological experiments to thoroughly present an alternative drug delivery device for BNCT which offers additional benefits over current delivery methods.

Graphical abstract The identification of discrete subsets of residues which respond to a local injection of thermal energy into the catalytic intermediates or cofactors within the active sites of GSAM, using non-equilibrium MD simulations of thermal diffusion, shows that the architecture of GSAM potentially supports a number of specific pathways for the flow of vibrational energy within and between the two subunits of the enzyme. Assuming that the thermal excitation with a heat bath acts as a physical proxy for chemical activity, these pathways may form the basis of a communication network for coordinating and triggering allosteric conformational changes in the two subunits by shuttling information about chemical modifications of the catalytic intermediates within the two active sites back and forth between the two subunits.

Polycyclic aromatic hydrocarbons (PAHs) are toxic, mutagenic and among the most damaging chemical compounds with regard to living organisms. Because of their persistence and wide distribution removal from the environment is an important challenge. Here we report a new Nano container matrix based on the deep sea archaea-derived RHCC-Nanotube (RHCC-NT), which rapidly and preferentially binds low molecular weight PAHs. Under controlled-laboratory conditions and using fluorescence spectroscopy in combination with X-ray crystallography and MD simulations, we quantified the real-time binding of low molecular weight PAHs (2–4 rings) to our substrate. Binding coefficients ranged from 5.4 ± 1.6 (fluorene) to 32 ± 7.0 μM (acenaphthylene) and a binding capacity of 85 pmoles PAH per mg RHCC-NT, or 2.12 μmoles in a standard 25 mg sampler. The uptake rate of pyrene was calculated to be 1.59 nmol/hr∙mol RHCC-NT (at 10  C). Our results clearly show that RHCC-NT is uniquely suited as a monitoring matrix for low molecular weight PAHs.

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 8 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 = 400 K 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.