David Butcher

Structures of the kinetically trapped i-motif DNA intermediates

In the present work, the conformational dynamics and folding pathways of i-motif DNA were studied in solution and in the gas-phase as a function of the solution pH conditions using circular dichroism (CD), photoacoustic calorimetry analysis (PAC), trapped ion mobility spectrometry-mass spectrometry (TIMS-MS), and molecular dynamics (MD). Solution studies showed at thermodynamic equilibrium the existence of a two-state folding mechanism, whereas during the pH = 7.0 → 4.5 transition a fast and slow phase (ΔHfast + ΔHslow = 43 ± 7 kcal mol-1) with a volume change associated with the formation of hemiprotonated cytosine base pairs and concomitant collapse of the i-motif oligonucleotide into a compact conformation were observed. TIMS-MS experiments showed that gas-phase, kinetically trapped i-motif DNA intermediates produced by nanoESI are preserved, with relative abundances depending on the solution pH conditions. In particular, a folded i-motif DNA structure was observed in nanoESI-TIMS-MS for low charge states in both positive and negative ion mode (e.g., z = ±3 to ±5) at low pH conditions. As solution pH increases, the cytosine neutralization leads to the loss of cytosine-cytosine+ (C·CH+) base pairing in the CCC strands and in those conditions we observe partially unfolded i-motif DNA conformations in nanoESI-TIMS-MS for higher charge states (e.g., z = -6 to -9). Collisional induced activation prior to TIMS-MS showed the existence of multiple local free energy minima, associated with the i-motif DNA unfolding at z = -6 charge state. For the first time, candidate gas-phase structures are proposed based on mobility measurements of the i-motif DNA unfolding pathway. Moreover, the inspection of partially unfolded i-motif DNA structures (z = -7 and z = -8 charge states) showed that the presence of inner cations may or may not induce conformational changes in the gas-phase. For example, incorporation of ammonium adducts does not lead to major conformational changes while sodium adducts may lead to the formation of sodium mediated bonds between two negatively charged sides inducing the stabilization towards more compact structures in new local, free energy minima in the gas-phase.

Role of Ionic Strength and pH in Modulating Thermodynamic Profiles Associated with CO Escape from Rice Nonsymbiotic Hemoglobin 1

Type 1 nonsymbiotic hemoglobins are found in a wide variety of land plants and exhibit very high affinities for exogenous gaseous ligands. These proteins are presumed to have a role in protecting plant cells from oxidative stress under etiolated/hypoxic conditions through NO dioxygenase activity. In this study we have employed photoacoustic calorimetry, time-resolved absorption spectroscopy, and classical molecular dynamics simulations in order to elucidate thermodynamics, kinetics, and ligand migration pathways upon CO photodissociation from WT and a H73L mutant of type 1 nonsymbiotic hemoglobin from Oryza sativa (rice). We observe a temperature dependence of the resolved thermodynamic parameters for CO photodissociation from CO-rHb1 which we attribute to temperature dependent formation of a network of electrostatic interactions in the vicinity of the heme propionate groups. We also observe slower ligand escape from the protein matrix under mildly acidic conditions in both the WT and H73L mutant (τ = 134 ± 19 and 90 ± 15 ns). Visualization of transient hydrophobic channels within our classical molecular dynamics trajectories allows us to attribute this phenomenon to a change in the ligand migration pathway which occurs upon protonation of the distal His73, His117, and His152. Protonation of these residues may be relevant to the functioning of the protein in vivo given that etiolation/hypoxia can cause a decrease in intracellular pH in plant cells.

Solution and gas-phase modifiers effect on heme proteins environment and conformational space

The molecular environment is known to impact the secondary and tertiary structure of biomolecules, shifting the equilibrium between different conformational and oligomerization states. In the present study, the effect of solution additives and gas-phase modifiers on the molecular environment of two common heme proteins, bovine cytochrome c and equine myoglobin, is investigated as a function of the time after desolvation (e.g., 100 - 500 ms) using trapped ion mobility spectrometry – mass spectrometry. Changes in the mobility profiles are observed depending on the starting solution composition (i.e., in aqueous solution at neutral pH or in the presence of organic content: methanol, acetone, or acetonitrile) depending on the protein. In the presence of gas-phase modifiers (i.e., N2 containing methanol, acetone, or acetonitrile), a shift in the mobility profiles driven by the gas-modifier mass and size and changes in the relative abundances and number of IMS bands are observed. We attribute these changes in the mobility profiles in the presence of gas-phase modifiers to a clustering/declustering mechanism by which organic molecules adsorb to the protein ion surface and lower energetic barriers for interconversion between conformational states, thus redefining the free energy landscape and equilibria between conformers. These structural biology experiments open new avenues for manipulation and interrogation of biomolecules in the gas-phase with the potential to emulate a large suite of solution conditions, ultimately including conditions that more accurately reflect a variety of intracellular environments.

Non-symbiotic Hemoglobin Conformational Space Dependence on the HEME Coordination using NESI-TIMS-TOF MS

In this study, for the first time, the conformational space of the rice non-symbiotic hemoglobin type 1 (rHb1) was studied as a function of the starting solution pH using trapped ion mobility spectrometry coupled to mass spectrometry (TIMS-MS) and molecular dynamics. Comparison of the charge state distribution, apo to holo form ratio, and the collision cross section ( ) profiles as a function of the solution pH showed higher stability of the rHb1 wild-type (WT) when compared with the H73L mutant at mildly acidic conditions. Comparison of the profiles of the rHb1 WT and H73L holo and apo form showed that only the initial unfolding pathways involved the heme cavity, with and without a heme loss, followed by unfolding pathways not necessarily involving the environment of the heme prosthetic group. Candidate structures for the nine transitions observed in the profiles were proposed using molecular nfolding dynamic simulations based on the profiles, UV absorption spectroscopy and circular dichroism data as way to describe a potential unfolding pathway. The described unfolding pathway suggests that the rHb1 unfolding is driven by initial distancing of the A, B, and H helices, while the heme cavity and heme group remains intact, followed by the distancing of the E, F, and G helices and subsequent loss of the -helical structure leading to a final random coil conformation.