Evidence is emerging that the role of protein structure in disease needs to be rethought. Sequence mutations in proteins are often found to affect the rate at which a protein switches between structures. Modeling structural transitions in wildtype and variant proteins is central to understanding the molecular basis of disease. This paper investigates an efficient algorithmic realization of the stochastic roadmap simulation framework to model structural transitions in wildtype and variants of proteins implicated in human disorders. Our results indicate that the algorithm is able to extract useful information on the impact of mutations on protein structure and function.

The atomic force microscope (AFM) allows biomolecules to be observed and manipulated under native conditions. It produces images with an outstanding signal-to-noise ratio and addresses single molecules while the sample is in a buffer solution. Progress in sample preparation and instrumentation has led to topographs that reveal subnanometer details and the surface dynamics of biomolecules. Tethering single molecules between a support and a retracting AFM tip produces force–extension curves, giving information about the mechanical stability of secondary structural elements. For both imaging and force spectroscopy, the cantilever and its tip are critical:the mechanical properties of the cantilever dictate the force sensitivity and the scanning speed, whereas the tip shape determines the achievable lateral resolution.

Backbone dynamics of the basic/helix-loop-helix domain of Pho4 from Saccharomyces cerevisae have been probed by NMR techniques, in the absence of DNA, nonspecifically bound to DNA and bound to cognate DNA. Alpha proton chemical shift indices and nuclear Overhauser effect patterns were used to elucidate the secondary structure in these states. These secondary structures are compared to the co-crystal complex of Pho4 bound to a cognate DNA sequence (Shimizu T, Toumoto A, Ihara K, Shimizu M, Kyogou Y, Ogawa N, Oshima Y, Hakoshima T, 1997, EMBO J 15: 4689–4697). The dynamic information provides insight into the nature of this DNA binding domain as it progresses from free in solution to a specifically bound DNA complex. Relative to the unbound form, we show that formation of either the nonspecific and cognate DNA bound complexes involves a large change in conformation and backbone dynamics of the basic region. The nonspecific and cognate complexes, however, have nearly identical secondary structure and backbone dynamics. We also present evidence for conformational flexibility at a highly conserved glutamate basic region residue. These results are discussed in relation to the mechanism of sequence specific recognition and binding.

Conformational changes are essential for the activity of many proteins. If, or how fast, internal fluctuations are related to slow conformational changes that mediate protein function is not understood. In this study, we measure internal fluctuations of the transport protein lactose permease in the presence and absence of substrate by tryptophan fluorescence spectroscopy. We demonstrate that nanosecond fluctuations of α-helices are enhanced when the enzyme transports substrate. This correlates with previously published kinetic data from transport measurements showing that millisecond conformational transitions of the substrate-loaded carrier are faster than those in the absence of substrate. These findings corroborate the hypothesis of the hierarchical model of protein dynamics that predicts that slow conformational transitions are based on fast, thermally activated internal motions.

The contributions of backbone NH group dynamics to the conformational heat capacity of the B1 domain of Streptococcal protein G have been estimated from the temperature dependence of 15N NMR-derived order parameters. Longitudinal (R1) and transverse (R2) relaxation rates, transverse cross-relaxation rates (ηxy), and steady state {1H}–15N nuclear Overhauser effects were measured at temperatures of 0, 10, 20, 30, 40, and 50 °C for 89–100% of the backbone secondary amide nitrogen nuclei in the B1 domain. The ratio R2/ηxy was used to identify nuclei for which conformational exchange makes a significant contribution to R2. Relaxation data were fit to the extended model-free dynamics formalism, incorporating an axially symmetric molecular rotational diffusion tensor. The temperature dependence of the order parameter (S2) was used to calculate the contribution of each NH group to conformational heat capacity (Cp) and a characteristic temperature (T*), representing the density of conformational energy states accessible to each NH group. The heat capacities of the secondary structure regions of the B1 domain are significantly higher than those of comparable regions of other proteins, whereas the heat capacities of less structured regions are similar to those in other proteins. The higher local heat capacities are estimated to contribute up to ∼0.8 kJ/mol K to the total heat capacity of the B1 domain, without which the denaturation temperature would be ∼9 °C lower (78 °C rather than 87 °C). Thus, variation of backbone conformational heat capacity of native proteins may be a novel mechanism that contributes to high temperature stabilization of proteins.

The backbone dynamics of the four-helical bundle cytokine leukemia inhibitory factor (LIF) have been investigated using 15N NMR relaxation and amide proton exchange measurements on a murine–human chimera, MH35-LIF. For rapid backbone motions (on a time scale of 10 ps to 100 ns), as probed by 15N relaxation measurements, the dynamics parameters were calculated using the model-free formalism incorporating the model selection approach. The principal components of the inertia tensor of MH35-LIF, as calculated from its NMR structure, were 1:0.98:0.38. The global rotational motion of the molecule was, therefore, assumed to be axially symmetric in the analysis of its relaxation data. This yielded a diffusion anisotropy D∥/D⊥ of 1.31 and an effective correlation time (4D⊥ + 2D∥)−1 of 8.9 ns. The average values of the order parameters (S2) for the four helices, the long interhelical loops, and the N-terminus were 0.91, 0.84, and 0.65, respectively, indicating that LIF is fairly rigid in solution, except at the N-terminus. The S2 values for the long interhelical loops of MH35-LIF were higher than those of their counterparts in short-chain members of the four-helical bundle cytokine family. Residues involved in LIF receptor binding showed no consistent pattern of backbone mobilities, with S2 values ranging from 0.71 to 0.95, but residues contributing to receptor binding site III had relatively lower S2 values, implying higher amplitude motions than for the backbone of sites I and II. In the relatively slow motion regime, backbone amide exchange measurements showed that a number of amides from the helical bundle exchanged extremely slowly, persisting for several months in 2H2O at 37 °C. Evidence for local unfolding was considered, and correlations among various structure-related parameters and the backbone amide exchange rates were examined. Both sets of data concur in showing that LIF is one of the most rigid four-helical bundle cytokines.

A method is presented that allows the calculation of the lifetimes of tryptophan residues on the basis of spectral and structural data. It is applied to two different proteins. The calcium binding protein from the sarcoplasm of the muscles of the sand worm Nereis diversicolor (NSCP) changes its conformation upon binding of Ca2+ or Mg2+. NSCP contains three tryptophan residues at position 4, 57, and 170, respectively. The fluorescence lifetimes of W57 are investigated in a mutant in which W4 and W170 have been replaced. The time resolved fluorescence properties of W57 are linked to its different microconformations, which were determined by a molecular dynamics simulation map. Together with the determination of the radiative rate constant from the wavelength of maximum intensity of the decay associated spectra, it was possible to determine an exponential relation between the nonradiative rate constant and the distance between the indole CE3 atom and the carbonyl carbon of the peptide bond reflecting a mechanism of electron transfer as the main determinant of the value for the nonradiative rate constant. This result allows the calculation of the fluorescence lifetimes from the protein structure and the spectra. This method was further tested for the tryptophan of Ha-ras p21 (W32) and for W43 of the Tet repressor, which resulted in acceptable values for the predicted lifetimes.