Structural and dynamic characteristics of a partially folded state of ubiquitin revealed by hydrogen exchange mass spectrometry

Structural and dynamic properties of a partially folded conformation (A-state) of ubiquitin are studied using amide hydrogen exchange in solution (HDX) and mass spectrometric detection. A clear distinction between the native state of the protein and the A-state can be made when HDX is carried out in a semicorrelated regime. Convoluted exchange patterns are interpreted with the aid of HDX simulations in a three-state system (highly structured, partially unstructured, and fully unstructured states). The data clearly indicate a highly dynamic character of the non-native state. Furthermore, combination of HDX and protein ion fragmentation in the gas phase [by means of collision-induced dissociation (CAD)] is used to evaluate the conformational stability of various protein segments specifically in the molten globular state. Chain flexibility appears to be distributed very unevenly in this non-native conformation. The highest degree of structural disorder is displayed by the C-terminal segment (Gly(53)-Gly(76)), which was previously suggested to form a transient alpha-helix. The least dynamic segment of ubiquitin in the A-state is Thr(9)-Glu(18) (which was previously suggested to form a stable nativelike beta-strand), with the adjacent segments exhibiting somewhat diminished conformational stability. The study also demonstrates the power of mass spectrometry as a tool in providing conformer-specific information about the structure and dynamics of both native and non-native protein states coexisting in solution under equilibrium.     

Quantitative comparison of the hydrogen bond network of A-state and native ubiquitin by hydrogen bond scalar couplings

The backbone hydrogen bond (H-bond) network of the partially folded A-state of ubiquitin (60% methanol, 40% water, pH 2) has been characterized quantitatively by (h3)J(NC)(') H-bond scalar couplings between the (15)N nuclei of amino acid H-bond donors and the (13)C carbonyl nuclei of the acceptors. Results on (h3)J(NC)(') couplings and the amide proton ((1)H(N)) chemical shifts for the A-state are compared quantitatively to the native state. The (h3)J(NC)(') correlations of the A-state show intact, nativelike H-bonds of the first beta-hairpin beta1/beta2 and the alpha-helix, albeit at lower strength, whereas the H-bonds in the C-terminal part change from a pure beta-structure to an all alpha-helical H(N)(i)-->O(i-4) connectivity pattern. A residue-specific analysis reveals that the conformations within the conserved secondary structure segments are much more homogeneous in the A-state than in the native state. Thus, the strong asymmetry of (h3)J(NC)(') couplings and (1)H(N) chemical shifts between the interior and exterior sides of the native state alpha-helix vanishes in the A-state. This indicates that the bend of this helix around the native state hydrophobic core is released in the homogeneous solvent environment of the A-state. Similarly, an irregularity in the behavior of H-bond I3-->L15 in hairpin beta1/beta2, which results from strong contacts to strand beta5 in the native state, is absent in the A-state. These findings rationalize the behavior of the (1)H(N) chemical shifts in both states and indicate that the A-state is in many aspects similar to the onset of thermal denaturation of the native state.     

Hsc70-interacting HPD loop of the J domain of polyomavirus T antigens fluctuates in ps to ns and micros to ms

The backbone dynamics of the J domain from polyomavirus T antigens have been investigated using 15N NMR relaxation and molecular dynamics simulation. Model-free relaxation analysis revealed picosecond to nanosecond motions in the N terminus, the I-II loop, the C-terminal end of helix II through the HPD loop to the beginning of helix III, and the C-terminal end of helix III to the C terminus. The backbone dynamics of the HPD loop and termini are dominated by motions with moderately large amplitudes and correlation times of the order of a nanosecond or longer. Conformational exchange on the microsecond to millisecond timescale was identified in the HPD loop, the N and C termini, and the I-II loop. A 9.7ns MD trajectory manifested concerted swings of the HPD loop. Transitions between major and minor conformations of the HPD loop featured distinct patterns of change in backbone dihedral angles and hydrogen bonds. Fraying of the C-terminal end of helix II and the N-terminal end of helix III correlated with displacements of the HPD loop. Correlation of crankshaft motions of Gly46 and Gly47 with the collective motions of the HPD loop suggested an important role of the two glycine residues in the mobility of the loop. Fluctuations of the HPD loop correlated with relative reorientation of side-chains of Lys35 and Asp44 that interact with Hsc70.     

Molecular dynamics simulations of the native and partially folded states of ubiquitin: influence of methanol cosolvent, pH, and temperature on the protein structure and dynamics

A series of explicit-solvent molecular dynamics simulations of the protein ubiquitin are reported, which investigate the effect of environmental factors (presence of methanol cosolvent in the aqueous solution, neutral or low pH value, room or elevated temperature) on the structure, stability, and dynamics of the protein. The simulations are initiated either from the native structure of the protein or from a model of a partially folded state (A-state) that is known to exist at low pH in methanol-water mixtures. The main results of the simulations are: (1) The ubiquitin native structure is remarkably stable at neutral pH in water; (2) the addition of the methanol cosolvent enhances the stability of the secondary structure but weakens tertiary interactions within the protein; (3) this influence of methanol on the protein structure is enhanced at low pH, while the effect of lowering the pH in pure water is limited; and (4) the A-state of ubiquitin can be described as a set of relatively rigid secondary structure elements (a native-like beta-sheet and native-like alpha-helix plus two nonnative alpha-helices) connected by flexible linkers.     

A comparison of the dynamic behavior of monomeric and dimeric insulin shows structural rearrangements in the active monomer

Molecular dynamics (MD) simulations (5-10ns in length) and normal mode analyses were performed for the monomer and dimer of native porcine insulin in aqueous solution; both starting structures were obtained from an insulin hexamer. Several simulations were done to confirm that the results obtained are meaningful. The insulin dimer is very stable during the simulation and remains very close to the starting X-ray structure; the RMS fluctuations calculated from the MD simulation agree with the experimental B-factors. Correlated motions were found within each of the two monomers; they can be explained by persistent non-bonded interactions and disulfide bridges. The correlated motions between residues B24 and B26 of the two monomers are due to non-bonded interactions between the side-chains and backbone atoms. For the isolated monomer in solution, the A chain and the helix of the B chain are found to be stable during 5ns and 10ns MD simulations. However, the N-terminal and the C-terminal parts of the B chain are very flexible. The C-terminal part of the B chain moves away from the X-ray conformation after 0.5-2.5ns and exposes the N-terminal residues of the A chain that are thought to be important for the binding of insulin to its receptor. Our results thus support the hypothesis that, when monomeric insulin is released from the hexamer (or the dimer in our study), the C-terminal end of the monomer (residues B25-B30) is rearranged to allow binding to the insulin receptor. The greater flexibility of the C-terminal part of the beta chain in the B24 (Phe-->Gly) mutant is in accord with the NMR results. The details of the backbone and side-chain motions are presented. The transition between the starting conformation and the more dynamic structure of the monomers is characterized by displacements of the backbone of Phe B25 and Tyr B26; of these, Phe B25 has been implicated in insulin activation.     

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin.

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.     

Investigation of the role of disulphide bond in modulating internal motions of BmK AGAP protein by molecular dynamics simulation

Elucidating structural determinants in the functional regions of toxins can provide useful knowledge for designing novel analgesic peptides. A series of 100 ns MD simulations were performed on the scorpion toxin BmK AGAP in native and disulphide bond broken states. The comparison of disulphide bond broken states with the native state showed the α-helix was found to be the key to the analgesic activity. Furthermore, our results revealed disulphide bonds have considerable influence on the functionally important essential modes of motions and the correlations between the motions of the Core domain and the C-terminal region which are involved in the analgesic activity. Therefore, we can conclude that disulphide bonds have a crucial role in modulating the function via adjusting the dynamics of scorpion toxin BmK AGAP molecule.     

Characterization of a partially denatured state of a protein by two-dimensional NMR: reduction of the hydrophobic interactions in ubiquitin

A stable, partially structured state of ubiquitin, the A-state, is formed at pH 2.0 in 60% methanol/40% water at 298 K. Detailed characterization of the structure of this state has been carried out by 2D NMR spectroscopy. Assignment of slowly exchanging amide resonances protected from the solvent in the native and A-state shows that gross structural reorganization of the protein has not occurred and that the A-state contains a subset of the interactions present in the native state (N-state). Vicinal coupling constants and NOESY data show the presence of the first two strands of the five-strand beta-sheet that is present in the native protein and part of the third beta-strand. The hydrophobic face of the beta-sheet in the A-state is covered by a partially structured alpha-helix, tentatively assigned to residues 24-34, that is considerably more flexible than the alpha-helix in the N-state. There is evidence for some fixed side-chain--side-chain interactions between these two units of structure. The turn-rich area of the protein, which contains seven reverse turns and a short piece of 3(10) helix, does not appear to be structured in the A-state and is approaching random coil.     

Dynamic and structural analysis of isotropically distributed molecular ensembles.

An efficient new method is presented for the characterization of motional correlations derived from a set of protein structures without requiring the separation of overall and internal motion. In this method, termed isotropically distributed ensemble (IDE) analysis, each structure is represented by an ensemble of isotropically distributed replicas corresponding to the situation found in an isotropic protein solution. This leads to a covariance matrix of the cartesian atomic positions with elements proportional to the ensemble average of scalar products of the position vectors with respect to the center of mass. Diagonalization of the covariance matrix yields eigenmodes and amplitudes that describe concerted motions of atoms, including overall rotational and intramolecular dynamics. It is demonstrated that this covariance matrix naturally distinguishes between "rigid" and "mobile" parts without necessitating a priori selection of a reference structure and an atom set for the orientational alignment process. The method was applied to the analysis of a 5-ns molecular dynamics trajectory of native ubiquitin and a 40-ns trajectory of a partially folded state of ubiquitin. The results were compared with essential dynamics analysis. By taking advantage of the spherical symmetry of the IDE covariance matrix, more than a 10-fold speed up is achieved for the computation of eigenmodes and mode amplitudes. IDE analysis is particularly suitable for studying the correlated dynamics of flexible and large molecules.     

Backbone dynamics of receptor binding and antigenic regions of a Pseudomonas aeruginosa pilin monomer

Pilin is the major structural protein that forms type IV pili of various pathogenic bacteria, including Pseudomonas aeruginosa. Pilin is involved in attachment of the bacterium to host cells during infection, in the initiation of immune response, and serves as a receptor for a variety of bacteriophage. We have used (15)N nuclear magnetic resonance relaxation measurements to probe the backbone dynamics of an N-terminally truncated monomeric pilin from P. aeruginosa strain K122-4. (15)N-T(1), -T(2), and [(1)H]-(15)N nuclear Overhauser enhancement measurements were carried out at three magnetic field strengths. The measurements were interpreted using the Lipari-Szabo model-free analysis, which reveals the amplitude of spatial restriction for backbone N-NH bond vectors with respect to nano- to picosecond time-scale motions. Regions of well-defined secondary structure exhibited consistently low-amplitude spatial fluctuations, while the terminal and loop regions showed larger amplitude motions in the subnano- to picosecond time-scale. Interestingly, the C-terminal disulfide loop region that contains the receptor binding domain was found to be relatively rigid on the pico- to nanosecond time-scale but exhibited motion in the micro- to millisecond time-scale. It is notable that this disulfide loop displays a conserved antigenic epitope and mediates binding to the asialo-GM(1) cell surface receptor. The present study suggests that a rigid backbone scaffold mediates attachment to the host cell receptor, and also maintains the conformation of the conserved antigenic epitope for antibody recognition. In addition, slower millisecond time-scale motions are likely to be crucial for conferring a range of specificity for these interactions. Characterization of pilin dynamics will aid in developing a detailed understanding of infection, and will facilitate the design of more efficient anti-adhesin synthetic vaccines and therapeutics against pathogenic bacteria containing type IV pili.     

An isolated helix persists in a sparsely populated form of KIX under native conditions

NMR relaxation dispersion techniques were used to investigate conformational exchange of the three-helix bundle protein KIX under native conditions. These experiments provide site-resolved kinetic information about microsecond-to-millisecond time scale motions along with structural (chemical shift) information without requiring a perturbation of the equilibrium. All kinetic data are consistent with an apparent two-state transition between natively folded KIX and a partially unfolded high-energy state that is populated to 3.0 +/- 0.2% at 27 degrees C. By combining (13)C- and (15)N-based experiments that probe specific structural aspects, we show that the sparsely populated high-energy state displays a strong conformational preference. An isolated secondary structural element, C-terminal helix alpha3, is highly populated, while the hydrophobic core of the domain and the remainder of the protein backbone, including helices alpha1 and alpha2, are disordered and devoid of specific interactions. This high-energy state presumably represents the equilibrium analogue of a folding intermediate that is transiently populated in stopped-flow kinetic experiments [Horng, J. C., Tracz, S. M., Lumb, K. J., and Raleigh, D. P. (2002) Biochemistry 44, 627-634].     

Reorientational contact-weighted elastic network model for the prediction of protein dynamics: comparison with NMR relaxation

A new model for the prediction of protein backbone motions is presented. The model, termed reorientational contact-weighted elastic network model, is based on a multidimensional reorientational harmonic potential of the backbone amide bond vector orientations and it is applied to the interpretation of dynamics parameters obtained from NMR relaxation data. The individual energy terms are weighted as a function of the intervector distances and by the contact strengths of each bond vector with respect to its local environment. Correlated reorientational motional properties of the bond vectors are obtained by means of normal mode analysis. Application to a set of proteins with known three-dimensional structures yields good to excellent agreement between predicted and experimental NMR order parameters presenting an improvement over the local contact model. The reorientational eigenmodes of the reorientational contact-weighted elastic network model method provide direct information on the collective nature of protein backbone motions. The dominant eigenmodes have a notably low collectivity, which is consistent with the behavior found for reorientational eigenmodes from molecular dynamics simulations.     

Thermodynamics of staphylococcal nuclease denaturation. II. The A-state.

Staphylococcal nuclease, at low pH and in the presence of high salt concentrations, has previously been proposed to exist in a partially folded or molten globule form called the "A-state" (Fink et al., 1993, Protein Sci 2:1155-1160). We have found that the A-state of nuclease at pH 2.1 in the presence of moderate to high salt concentrations and at low temperature exists in a substantially folded form structurally more similar to a native state. The A-state has the far-UV circular dichroism spectra characteristic of the native protein, which indicates that it has a large degree of secondary structure. Upon heating, the A-state denatures with a sigmoidal change in far-UV ellipticity and an observable peak in a differential scanning calorimeter trace, indicating that it is thermodynamically distinct from the denatured state. Three different mutations in a residue normally buried in the protein's core stabilize or destabilize the A-state in the same way as they affect the denaturation of the native state. The A-state must, therefore, contain at least some tertiary packing of side chains. Unlike the native state, which shows cold denaturation at low temperatures, the A-state is most stable at temperatures below 0 degrees C.     

In-cell protein dynamics

The native intracellular environment of proteins is crowded with metabolites and macromolecules. However, most biophysical information concerning proteins is acquired in dilute solution. To determine whether there are differences in dynamics, nuclear magnetic resonance spectroscopy can be used to measure 15N relaxation in uniformly 15N-enriched apocytochrome b5 inside living Escherichia coli and in dilute solution. Such data can then be used to compare the fast backbone dynamics of the partially folded protein in cells to its dynamics in dilute solution by using Lipari-Szabo analysis. It appears that the intracellular environment does not alter the protein's structure, or significantly change its fast dynamics. Specifically, the cytosol does not change the amplitude of fast backbone motions, but does increase the average timescale of these motions, most likely due to the increase in viscosity of the cytosol.     

Secondary and tertiary structure of the A-state of cytochrome c from resonance Raman spectroscopy.

Ferricytochrome c can be converted to the partially folded A-state at pH 2.2 in the presence of 1.5 M NaCl. The structure of the A-state has been studied in comparison with the native and unfolded states, using resonance Raman spectroscopy with visible and ultraviolet excitation wavelengths. Spectra obtained with 200 nm excitation show a decrease in amide II intensity consistent with loss of structure for the 50s and 70s helices. The 230-nm spectra contain information on vibrational modes of the single Trp 59 side chain and the four tyrosine side chains (Tyr 48, 67, 74, and 97). The Trp 59 modes indicate that the side chain remains in a hydrophobic environment but loses its tertiary hydrogen bond and is rotationally disordered. The tyrosine modes Y8b and Y9a show disruption of tertiary hydrogen bonding for the Tyr 48, 67, and 74 side chains. The high-wavenumber region of the 406.7-nm resonance Raman spectrum reveals a mixed spin heme iron atom, which arises from axial coordination to His 18 and a water molecule. The low-frequency spectral region reports on heme distortions and indicates a reduced degree of interaction between the heme and the polypeptide chain. A structural model for the A-state is proposed in which a folded protein subdomain, consisting of the heme and the N-terminal, C-terminal, and 60s helices, is stabilized through nonbonding interactions between helices and with the heme.     

Efficient RMSD measures for the comparison of two molecular ensembles. Root-mean-square deviation

Quantitative measures are presented for comparing the conformations of two molecular ensembles. The measures are based on Kabsch's formula for the root-mean-square deviation (RMSD) and the covariance matrix of atomic positions of isotropically distributed ensembles (IDE). By using a Taylor series expansion, it is shown that the RMSD can be expressed solely in terms of the IDE matrices. A fast approximate method is introduced for the pairwise RMSD determination whose computational cost scales linearly with the number of structures. A similarity measure for two structural ensembles that is based on the trace metric of the differences of powers of the IDE matrices is presented. The measures are illustrated for conformational ensembles generated by a molecular dynamics computer simulation of a partially folded A-state analog of ubiquitin.     

Accurate and efficient description of protein vibrational dynamics: comparing molecular dynamics and Gaussian models

Current all-atom potential based molecular dynamics (MD) allows the identification of a protein's functional motions on a wide-range of timescales, up to few tens of nanoseconds. However, functional, large-scale motions of proteins may occur on a timescale currently not accessible by all-atom potential based MD. To avoid the massive computational effort required by this approach, several simplified schemes have been introduced. One of the most satisfactory is the Gaussian network approach based on the energy expansion in terms of the deviation of the protein backbone from its native configuration. Here, we consider an extension of this model that captures in a more realistic way the distribution of native interactions due to the introduction of effective side-chain centroids. Since their location is entirely determined by the protein backbone, the model is amenable to the same exact and computationally efficient treatment as previous simpler models. The ability of the model to describe the correlated motion of protein residues in thermodynamic equilibrium is established through a series of successful comparisons with an extensive (14 ns) MD simulation based on the AMBER potential of HIV-1 protease in complex with a peptide substrate. Thus, the model presented here emerges as a powerful tool to provide preliminary, fast yet accurate characterizations of protein near-native motion.     

Pressure effects on the ensemble dynamics of ubiquitin inspected with molecular dynamics simulations and isotropic reorientational eigenmode dynamics

According to NMR chemical shift data, the ensemble of ubiquitin is a mixture of “open” and “closed” conformations at rapid equilibrium. Pressure perturbations provide the means to study the transition between the two conformers by imposing an additional constraint on the system's partial molar volume. Here we use nanosecond-timescale molecular dynamics simulations to characterize the network of correlated motions accessible to the conformers at low- and high-pressure conditions. Using the isotropic reorientational eigenmode dynamics formalism to analyze our simulation trajectories, we reproduce NMR relaxation data without fitting any parameters of our model. Comparative analysis of our results suggests that the two conformations behave very differently. The dynamics of the “closed” conformation are almost unaffected by pressure and are dominated by large-amplitude correlated motions of residues 23-34 in the extended α-helix. The “open” conformation under conditions of normal pressure displays increased mobility, focused on the loop residues 17-20, 46-55, and 58-59 at the bottom of the core of the structure, as well as the C-terminal residues 69-76, that directly participate in key protein-protein interactions. For the same conformation, a pressure increase induces a loss of separability between molecular tumbling and internal dynamics, while motions between different backbone sites become uncorrelated.