Detailed unfolding and folding of gaseous ubiquitin ions characterized by electron capture dissociation

The unfolding enthalpy of the native state of ubiquitin in solution is 5 to 8 times that of its gaseous ions, as determined by electron capture dissociation (ECD) mass spectrometry. Although two-state folding occurs in solution, the three-state gaseous process proposed for this by Clemmer and co-workers based on ion mobility data is supported in general by ECD mass spectra, including relative product yields, distinct Delta H(unfolding) values between states, site-specific melting temperatures, and folding kinetics indicating a cooperative process. ECD also confirms that the 13+ ions represent separate conformers, possibly with side-chain solvated alpha-helical structures. However, the ECD data on the noncovalent bonding in the 5+ to 13+ ions, determined overall in 69 of the 75 interresidue sites, shows that thermal unfolding proceeds via a diversity of intermediates whose conformational characteristics also depend on charge site locations. As occurs with increased acidity in solution, adding 6 protons to the 5+ ions completely destroys their tertiary noncovalent bonding. However, solvation of the newly protonated sites to the backbone instead increases the stability of the secondary structure (possibly an alpha-helix) of these gaseous ions, while in solution these new sites aid denaturation by solvation in the aqueous medium. Extensive ion equilibration can lead to even more compact and diverse conformers. The three-state unfolding of gaseous ubiquitin appears to involve ensembles of individual chain conformations in a "folding funnel" of parallel reaction paths. This also provides a further caution for characterizing solution conformers from their gas-phase behavior.     

Combined Chemical Modification and Collision Induced Unfolding Using Native Ion Mobility-Mass Spectrometry Provides Insights into Protein Gas-Phase Structure

Native mass spectrometry is now an important tool in structural biology. Thus, the nature of higher protein structure in the vacuum of the mass spectrometer is an area of significant interest. One of the major goals in the study of gas-phase protein structure is to elucidate the stabilising role of interactions at the level of individual amino acid residues. A strategy combining protein chemical modification together with collision induced unfolding (CIU) was developed and employed to probe the structure of compact protein ions produced by native electrospray ionisation. Tractable chemical modification was used to alter the properties of amino acid residues, and ion mobility-mass spectrometry (IM-MS) utilised to monitor the extent of unfolding as a function of modification. From these data the importance of specific intramolecular interactions for the stability of compact gas-phase protein structure can be inferred. Using this approach, and aided by molecular dynamics simulations, an important stabilising interaction between K6 and H68 in the protein ubiquitin was identified, as was a contact between the N-terminus and E22 in a ubiquitin binding protein UBA2.     

Molecular dynamics simulation of the structure, dynamics, and thermostability of the RNA hairpins uCACGg and cUUCGg

Classical replica-exchange molecular dynamics simulations are performed to study structure, dynamics and thermostability of the 14-mer RNA hairpins uCACGg and cUUCGg. Despite of the different sequence and closing base pair of the two systems, recent NMR studies have shown that the tetraloop CACG is strikingly similar in overall geometry and hydrogen bonding to the canonical UUCG tetraloop. On the other hand, the two systems differ significantly in their functionality and thermostability. The simulations confirm the structural similarities of the two RNA hairpins at room temperature but also reveal that the UUCG loop is more flexible than the CACG loop. Concerning the functionality, the CACG loop shows a stronger attitude to donate hydrogens than the UUCG loop, although their global solvent accessible surface is quite similar. The simulations qualitatively reproduce the experimentally found difference in melting temperatures (20 K). In the case of the uCACGg hairpin, the thermal unfolding occurs cooperatively in an all-or-none fashion, while the cUUCGg hairpin shows less cooperativity but exhibits intermediate states during the unfolding process.     

Coupling between hydration layer dynamics and unfolding kinetics of HP-36

We have performed atomistic molecular dynamics simulations of aqueous solutions of HP-36 at 300 K in its native state, as well as at high temperatures to explore the unfolding dynamics of the protein and its correlation with the motion of water around it. On increasing the temperature a partially unfolded molten globule state is formed where the smallest alpha helix (helix 2) unfolds into a coil. It is observed that the unfolding is initiated around the residue Phe-18 which shows a sharp displacement during unfolding. We have noticed that the unfolding of the protein affects the density of water near the protein surface. Besides, the dynamics of water in the protein hydration layer has been found to be strongly correlated with the time evolution of the unfolding process. We have introduced and calculated a displacement time correlation function to monitor the change in water motion relative to the protein backbone during unfolding. We find that the unfolding of helix 2 is associated with an increase in mobility of water around it as compared to water around the other two helices. We have also explored the microscopic aspects of secondary structure specific and site specific solvation dynamics of the protein. The calculations reveal that unfolding influences the solvation dynamics of the protein molecule in a heterogeneous manner depending on the location of the polar probe residues. This seems to be in agreement with recent experimental findings.     

Molecular dynamics simulation of thermal unfolding of Thermatoga maritima DHFR

Molecular dynamics simulations of the temperature-induced unfolding reaction of native dimeric dihydrofolate reductase from the hyperthermophile Thermatoga maritima (TmDHFR) and the experimentally inaccessible TmDHFR monomer were carried out at 400 K, 450 K and 500 K. The results revealed that the unfolding of TmDHFR subunits followed a similar path to that of the monomeric DHFR from the mesophile E. coli (EcDHFR). An initial collapse of the adenosine-binding domain (ABD) was followed by the loss of the N-terminal and loop domains (NDLD). Interestingly, the elements of the secondary structure of the isolated TmDHFR monomer were maintained for significantly longer periods of time for the hyperthermophilic enzyme, suggesting that subunit stability contributes to the enhanced resistance of TmDHFR to temperature-induced unfolding. The interactions between the subunits of the TmDHFR dimer led to a stabilisation of the NDLD. The hydrogen bonds between residues 140-143 in betaG of one subunit and residues 125-127 in betaF of the other subunit were retained for significant parts of the simulations at all temperatures. These intermolecular hydrogen bonds were lost after the unfolding of the individual subunits. The high stability of the dimer mediated by strong intersubunit contacts together with an intrinsically enhanced stability of the subunits compared to EcDHFR provides a molecular rational for the higher stability of the thermophilic enzyme. The computed unfolding pathways suggest that the partly folded dimer may be a genuine folding intermediate.     

Pulling direction as a reaction coordinate for the mechanical unfolding of single molecules

The folding and unfolding kinetics of single molecules, such as proteins or nucleic acids, can be explored by mechanical pulling experiments. Determining intrinsic kinetic information, at zero stretching force, usually requires an extrapolation by fitting a theoretical model. Here, we apply a recent theoretical approach describing molecular rupture in the presence of force to unfolding kinetic data obtained from coarse-grained simulations of ubiquitin. Unfolding rates calculated from simulations over a broad range of stretching forces, for different pulling directions, reveal a remarkable "turnover" from a force-independent process at low force to a force-dependent process at high force, akin to the "roll-over" in unfolding rates sometimes seen in studies using chemical denaturant. While such a turnover in rates is unexpected in one dimension, we demonstrate that it can occur for dynamics in just two dimensions. We relate the turnover to the quality of the pulling direction as a reaction coordinate for the intrinsic folding mechanism. A novel pulling direction, designed to be the most relevant to the intrinsic folding pathway, results in the smallest turnover. Our results are in accord with protein engineering experiments and simulations which indicate that the unfolding mechanism at high force can differ from the intrinsic mechanism. The apparent similarity between extrapolated and intrinsic rates in experiments, unexpected for different unfolding barriers, can be explained if the turnover occurs at low forces.     

Molecular dynamics simulation of NMR powder lineshapes of linear guests in structure I clathrate hydrates

We perform molecular dynamics simulations (up to 6 ns) for the structure I clathrate hydrates of linear molecules CS, CS(2), OCS, and C(2)H(2) in large cages at different temperatures in the stability range to determine the angular distribution and dynamics of the guests in the large cages. The long axes of linear guest molecules in the oblate large structure I clathrate hydrate cages are primarily confined near the equatorial plane of the cage rather than axial regions. This non-uniform spatial distribution leads to well-known anisotropic lineshapes in the solid-state NMR spectra of the guest species. We use the dynamic distribution of guest orientations in the cages during the MD simulations at different temperatures to predict the (13)C NMR powder lineshapes of the guests in the large cages. The length of the guests and intermolecular interactions of the guests in the water cages determine the angular distribution and the mobility of the guests in the sI large cages at different temperatures. At low temperatures the range of motion of the guests in the cages are limited and this is reflected in the skew of the predicted (13)C lineshapes. As the guest molecules reach the fast motion limit at higher temperatures, the lineshapes for CS, OCS, and C(2)H(2) are predicted to have the "standard" powder lineshapes of guest molecules.     

Computer simulations of the refolding of sperm whale apomyoglobin from high-temperature denaturated state

The refolding mechanism of apomyoglobin (apoMb) subsequent to high-temperature unfolding has been examined using computer simulations with atomic level detail. The folding of this protein has been extensively studied experimentally, providing a large database of folding parameters which can be probed using simulations. In the present study, 4-folding trajectories of apoMb were computed starting from coiled structures. A crystal structure of sperm whale myoglobin taken from the Protein Data Bank was used to construct the final native conformation by removal of the heme group followed by energy optimization. The initial unfolded conformations were obtained from high-temperature molecular dynamics simulations. Room-temperature refolding trajectories at neutral pH were obtained using the stochastic difference equation in length algorithm. The folding trajectories were compared with experimental results and two previous molecular dynamics studies at low pH. In contrast to the previous simulations, an extended intermediate with large helical content was not observed. In the present study, a structural collapse occurs without formation of helices or native contacts. Once the protein structure is more compact (radius of gyration<18 A) secondary and tertiary structures appear. These results suggest that apoMb follows a different folding pathway after high-temperature denaturation.     

Transfer of structural elements from compact to extended states in unsolvated ubiquitin

Multidimensional ion mobility spectrometry techniques (IMS-IMS and IMS-IMS-IMS) combined with mass spectrometry are used to study structural transitions of ubiquitin ions in the gas phase. It is possible to select and activate narrow distributions of compact and partially folded conformation types and examine new distributions of structures that are formed. Different compact conformations unfold, producing a range of new partially folded states and three resolvable peaks associated with elongated conformers. Under gentle activation conditions, the final populations of the three elongated forms depend on the initial structures of the selected ions. This requires that some memory of the compact state (most likely secondary structure) is preserved along the unfolding pathway. Activation of selected, partially folded intermediates (formed from specific compact states) leads to elongated state populations that are consistent with the initial selected compact form-evidence that intermediates not only retain elements of initial structure but also are capable of transmitting structure to final states.     

Quantum dynamics study of Li + HF reaction

We report rigorous quantum dynamics studies of the Li + HF reaction using the time-dependent wavepacket approach. The dynamics study is carried out on a recent ab initio potential energy surface, and state-selected reaction probabilities and cross sections are calculated up to 0.4 eV of collision energy. Many long-lived resonances (as long as 10 ps) at low collision energies (below 0.1 eV) are uncovered from the dynamics calculation. These long-lived resonances play a dominant role in the title reaction at low collision energies (below 0.1 eV). At higher energies, the direct reaction process becomes very important. The reaction probabilities from even rotational states exhibit a different energy dependence than those from odd rotational states. Our calculated integral cross section exhibits a broad maximum near the collision energy of 0.26 eV with small oscillations superimposed on the broad envelope which is reminiscent of the underlying resonance structures in reaction probabilities. The energy dependence of the present CS cross section is qualitatively different from the simple J-shifting approximation, in which a monotonic increase of cross section with collision energy was obtained. Received: 8 January 1997 / Accepted: 14 January 1997     

Hydrogen bond dynamics in heavy water studied with quantum dynamical simulations

The structure and dynamics of the hydrogen-bond network in heavy water (D(2)O) is studied as a function of the temperature using quantum dynamical simulations. Our approach combines an ab initio-based representation of the water interactions with an explicit quantum treatment of the molecular motion. A direct connection between the calculated linear and nonlinear vibrational spectra and the underlying molecular dynamics is made, which provides new insights into the rearrangement of the hydrogen-bond network in heavy water. A comparison with previous calculations on liquid H(2)O suggests that tunneling does not effectively contribute to the dynamics of the water hydrogen-bond network above the melting point. However, the effects of nuclear quantization are not negligible at all temperatures and become increasingly important near the melting point, which is in agreement with recent experimental analysis of the structural properties of liquid water as well as of the proton momentum distribution in supercooled water.     

Insights into the dynamics of HIV-1 protease: a kinetic network model constructed from atomistic simulations

The conformational dynamics in the flaps of HIV-1 protease plays a crucial role in the mechanism of substrate binding. We develop a kinetic network model, constructed from detailed atomistic simulations, to determine the kinetic mechanisms of the conformational transitions in HIV-1 PR. To overcome the time scale limitation of conventional molecular dynamics (MD) simulations, our method combines replica exchange MD with transition path theory (TPT) to study the diversity and temperature dependence of the pathways connecting functionally important states of the protease. At low temperatures the large-scale flap opening is dominated by a small number of paths; at elevated temperatures the transition occurs through many structurally heterogeneous routes. The expanded conformation in the crystal structure 1TW7 is found to closely mimic a key intermediate in the flap-opening pathways at low temperature. We investigated the different transition mechanisms between the semi-open and closed forms. The calculated relaxation times reveal fast semi-open ? closed transitions, and infrequently the flaps fully open. The ligand binding rate predicted from this kinetic model increases by 38-fold from 285 to 309 K, which is in general agreement with experiments. To our knowledge, this is the first application of a network model constructed from atomistic simulations together with TPT to analyze conformational changes between different functional states of a natively folded protein.     

Intramolecular vibrational redistribution in the vibrational predissociation dynamics of He-I2

The intramolecular vibrational redistribution (IVR) process is investigated in wave packet simulations of the vibrational predissociation dynamics of He-I(2)(B,upsilon') in the region of high upsilon' levels, upsilon' = 35-65. The simulations indicate that for upsilon' < or = 45 the dynamics is dominated by direct predissociation, whereas for higher upsilon' levels the onset of IVR appears and becomes increasingly important. The IVR process occurs via coupling of the initial state in the upsilon' manifold to intermediate long-lived resonances belonging to the lower upsilon < upsilon' vibrational manifolds. The IVR dynamics manifests itself in multiexponential behavior and oscillations in the time-dependent population curves associated with the He-I(2)(B,upsilon') initial state, the He-I(2)(B,upsilon < upsilon') intermediate complexes, and the final product states. The population curves corresponding to the upsilon'- 1 intermediate resonances located below the He + I(2)(B,upsilon'-1,j=0) dissociation limit are analyzed. It is found that initial population is transferred to all the intermediate resonance states considered, including those more separated in energy from the initial one. The results obtained for population transfer between the initial and the intermediate states can be explained by the intensity of the matrix elements coupling the initial and the intermediate resonances, in combination with the Rabi's formula for population exchange between two coupled states.     

Semiclassical nonadiabatic dynamics based on quantum trajectories for the O(3P,1D) + H2 system

The O(3P,1D) + H2 --> OH + H reaction is studied using trajectory dynamics within the approximate quantum potential approach. Calculations of the wave-packet reaction probabilities are performed for four coupled electronic states for total angular momentum J = 0 using a mixed coordinate/polar representation of the wave function. Semiclassical dynamics is based on a single set of trajectories evolving on an effective potential-energy surface and in the presence of the approximate quantum potential. Population functions associated with each trajectory are computed for each electronic state. The effective surface is a linear combination of the electronic states with the contributions of individual components defined by their time-dependent average populations. The wave-packet reaction probabilities are in good agreement with the quantum-mechanical results. Intersystem crossing is found to have negligible effect on reaction probabilities summed over final electronic states.     

Charge-state dependent compaction and dissociation of protein complexes: insights from ion mobility and molecular dynamics

Collapse to compact states in the gas phase, with smaller collision cross sections than calculated for their native-like structure, has been reported previously for some protein complexes although not rationalized. Here we combine experimental and theoretical studies to investigate the gas-phase structures of four multimeric protein complexes during collisional activation. Importantly, using ion mobility-mass spectrometry (IM-MS), we find that all four macromolecular complexes retain their native-like topologies at low energy. Upon increasing the collision energy, two of the four complexes adopt a more compact state. This collapse was most noticeable for pentameric serum amyloid P (SAP) which contains a large central cavity. The extent of collapse was found to be highly correlated with charge state, with the surprising observation that the lowest charge states were those which experience the greatest degree of compaction. We compared these experimental results with in vacuo molecular dynamics (MD) simulations of SAP, during which the temperature was increased. Simulations showed that low charge states of SAP exhibited compact states, corresponding to collapse of the ring, while intermediate and high charge states unfolded to more extended structures, maintaining their ring-like topology, as observed experimentally. To simulate the collision-induced dissociation (CID) of different charge states of SAP, we used MS to measure the charge state of the ejected monomer and assigned this charge to one subunit, distributing the residual charges evenly among the remaining four subunits. Under these conditions, MD simulations captured the unfolding and ejection of a single subunit for intermediate charge states of SAP. The highest charge states recapitulated the ejection of compact monomers and dimers, which we observed in CID experiments of high charge states of SAP, accessed by supercharging. This strong correlation between theory and experiment has implications for further studies as well as for understanding the process of CID and for applications to gas-phase structural biology more generally.     

Proteins in a shear flow

The conformational dynamics of a single protein molecule in a shear flow is investigated using Brownian dynamics simulations. A structure-based coarse grained model of a protein is used. We consider two proteins, ubiquitin and integrin, and find that at moderate shear rates they unfold through a sequence of metastable states-a pattern which is distinct from a smooth unraveling found in homopolymers. Full unfolding occurs only at very large shear rates. Furthermore, the hydrodynamic interactions between the amino acids are shown to hinder the shear flow unfolding. The characteristics of the unfolding process depend on whether a protein is anchored or not, and if it is, on the choice of an anchoring point.     

Temperature-dependent H/D exchange of compact and elongated cytochrome c ions in the gas phase

Isotopic exchange reactions of compact and elongated conformations of gaseous cytochrome c ions (+5 and +9 states) with D2O have been measured as a function of temperature (from 300 to approximately 440 K) using ion mobility techniques. Rate constants for those sites that exchange at high temperatures (>400 K) are about an order of magnitude smaller than rate constants for sites that exchange at 300 K. Although the exchange rates decrease, the maximum exchange levels for rapidly exchanging sites increase with temperature. At 300 K, exchange levels of 53 +/- 3 and 63 +/- 3 are measured for the compact and elongated states, respectively. From 300 to 335 K, the exchange levels increase slightly to approximately 60 to 70 hydrogens. Above 335 K, the levels increase to a value of approximately 200 for the +5 state and approximately 190 for the +9 state, near the maximum possible levels, 200 and 204 for these respective charge states. Molecular dynamics simulations have been carried out on structures having calculated cross sections that are near the experimental values in order to explore the exchange process. Overall, it appears that charge site and exchange site proximities are important factors in the exchange profiles for the elongated +9 ion and the compact +5 ion.     

The energy landscape of unsolvated peptides: helix formation and cold denaturation in Ac-A4G7A4 + H+

Ion mobility measurements and molecular dynamics simulations were performed for unsolvated A4G7A4 + H+ and Ac-A4G7A4 + H+ (Ac = acetyl, A = alanine, G = glycine) peptides. As expected, A4G7A4 + H+ adopts a globular conformation (a compact, random-looking, three-dimensional structure) over the entire temperature range examined (100-410 K). Ac-A4G7A4 + H+ on the other hand is designed to have a flat energy landscape with a marginally stable helical state. This peptide shows at least four different conformations at low temperatures (<230 K). The two conformations with the largest cross sections are attributed to - and partial -helices, while the one with the smallest cross section is globular. The other main conformation may be partially helical. Ac-A4G7A4 + H+ becomes predominantly globular at intermediate temperatures and then becomes more helical as the temperature is raised further. This unexpected behavior may be due to the helix having a higher vibrational entropy than the globular state, as predicted by some recent calculations (Ma, B.; Tsai, C.-J.; Nussinov, R. Biophys. J. 2000, 79, 2739-2753).     

Guanidinium chloride molecular diffusion in aqueous and mixed water-ethanol solutions

Solutions containing guanidinium chloride (GdmCl), or equivalently guanidine hydrochloride (GdnHCl), are commonly used to denature macromolecules such as proteins and DNA in, for example, microfluidics studies of protein unfolding. To design and study such applications, it is necessary to know the diffusion coefficients for GdmCl in the solution. To this end, we use molecular dynamics simulations to calculate the diffusion coefficients of GdmCl in water and in water-ethanol solutions, for which no direct experimental measurements exist. The fully atomistic simulations show that the guandinium cation Gdm (+) diffusion decreases as the concentration of both Gdm (+) and ethanol in the solution increases. The simulations are validated against available literature data, both transformed measured viscosity values and computed diffusion coefficients, and we show that a prudent choice of water model, namely TIP4P-Ew, gives calculated diffusion coefficients in good agreement with the transformed measured viscosity values. The calculated Gdm (+) diffusion behavior is explained as a dynamic mixture of free cation, stacked cation, and ion-paired species in solution, with weighted contributions to Gdm (+) diffusion from the stacked and paired states helping explain measured viscosity data in terms of atom-scale dynamics.     

Acetonitrile-induced unfolding of the photosystem II manganese-stabilizing protein studied by electrospray mass spectrometry

In this paper an acetonitrile-induced unfolding of the manganese-stabilizing protein (MSP) of photosystem II was discovered. More distinct unfolding states of MSP were identified than previously by using mainly electrospray ionization mass spectrometry (ESI-MS), together with fluorescence spectra and far-UV circular dichroism (CD) at pH 2.0, 6.2 or 11.6, and with acetonitrile concentrations from 0 to 50%. At pH 6.2 with acetonitrile concentration changing from 0 to 10%, relatively broad charge-state distributions and poor intensity were observed in ESI-MS, indicating the presence of coexisting conformers. It was concluded that the structure of the MSP protein is unlikely to be a tightly folded form. When the concentration of acetonitrile was 20-40%, simulating the state in the biological membrane, changes in the state of unfolding of MSP were observed to a certain extent using ESI-MS, fluorescence and CD spectroscopy. The charge-state distribution in ESI-MS was found to move toward high states (from 13+ to 27+ to 15+ to 31+) with increasing acetonitrile concentration. At pH 2.0, the MSP structure is rearranged into an unfolded state, and at pH 11.6 the MSP structure is induced to assume another unordered state by deprotonation of appropriate residues. An interesting observation was that a second peak envelope emerged with 20-50% acetonitrile in the medium at pH 11.6.