Osmolyte perturbation reveals conformational equilibria in spin‐labeled proteins

Recent evidence suggests that proteins at equilibrium can exist in a manifold of conformational substates, and that these substates play important roles in protein function. Therefore, there is great interest in identifying regions in proteins that are in conformational exchange. Electron paramagnetic resonance spectra of spin‐labeled proteins containing the nitroxide side chain (R1) often consist of two (or more) components that may arise from slow exchange between conformational substates (lifetimes > 100 ns). However, crystal structures of proteins containing R1 have shown that multicomponent spectra can also arise from equilibria between rotamers of the side chain itself. In this report, it is shown that these scenarios can be distinguished by the response of the system to solvent perturbation with stabilizing osmolytes such as sucrose. Thus, site‐directed spin labeling (SDSL) emerges as a new tool to explore slow conformational exchange in proteins of arbitrary size, including membrane proteins in a native‐like environment. Moreover, equilibrium between substates with even modest differences in conformation is revealed, and the simplicity of the method makes it suitable for facile screening of multiple proteins. Together with previously developed strategies for monitoring picosecond to millisecond backbone dynamics, the results presented here expand the timescale over which SDSL can be used to explore protein flexibility.     

Hierarchical Conformational Analysis of Native Lysozyme Based on Sub-Millisecond Molecular Dynamics Simulations

Hierarchical organization of free energy landscape (FEL) for native globular proteins has been widely accepted by the biophysics community. However, FEL of native proteins is usually projected onto one or a few dimensions. Here we generated collectively 0.2 milli-second molecular dynamics simulation trajectories in explicit solvent for hen egg white lysozyme (HEWL), and carried out detailed conformational analysis based on backbone torsional degrees of freedom (DOF). Our results demonstrated that at micro-second and coarser temporal resolutions, FEL of HEWL exhibits hub-like topology with crystal structures occupying the dominant structural ensemble that serves as the hub of conformational transitions. However, at 100ns and finer temporal resolutions, conformational substates of HEWL exhibit network-like topology, crystal structures are associated with kinetic traps that are important but not dominant ensembles. Backbone torsional state transitions on time scales ranging from nanoseconds to beyond microseconds were found to be associated with various types of molecular interactions. Even at nanoseconds temporal resolution, the number of conformational substates that are of statistical significance is quite limited. These observations suggest that detailed analysis of conformational substates at multiple temporal resolutions is both important and feasible. Transition state ensembles among various conformational substates at microsecond temporal resolution were observed to be considerably disordered. Life times of these transition state ensembles are found to be nearly independent of the time scales of the participating torsional DOFs.     

Conformational substates of myoglobin detected by extrinsic dynamic fluorescence studies

The extent of conformational substates of two apomyoglobins, i.e., sperm whale and tuna apomyoglobin, was investigated by examining the fluorescence decay in the frequency domain of the extrinsic fluorophore TNS [6-(p-toluidino)-2-naphthalenesulfonic acid] bound to the heme binding site. Data analysis was performed in terms of a continuous, unimodal lifetime distribution having a Lorentzian shape. The results were compared with those for the free fluorophore in an isotropic nonviscous solvent. The incorporation of TNS into the protein matrix resulted in a broadening of the lifetime distribution due to the microenvironmental heterogeneity generated by structural fluctuations. The larger width of lifetime distribution observed for TNS bound to tuna apomyoglobin was related to a more extended conformational space accessible to the fluorophore in this protein compared to sperm whale myoglobin. A temperature increase from 15 to 40 degrees C produced a further broadening of the lifetime distributions of TNS bound to both proteins. This result can be explained by assuming the existence of conformational substates at high energy content or separated by high energy barriers, which are not populated at low temperature. The overall picture emerging from the reported data is that the lifetime distributions of TNS bound to apomyoglobins are determined largely by the number of conformational substates accessible to the protein matrix and, to a lesser extent, by the interconversion rates among these states.     

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.     

Osmolytes modulate conformational exchange in solvent‐exposed regions of membrane proteins

Site‐directed spin labeling (SDSL) was used to investigate local structure and conformational exchange in two bacterial outer‐membrane TonB‐dependent transporters, BtuB and FecA. Protecting osmolytes, such as polyethylene glycols (PEGs) are known to modulate a substrate‐dependent conformational equilibrium in the energy coupling motif (Ton box) of BtuB. Here, we demonstrate that a segment that is N‐terminal to the Ton box in BtuB, is in conformational exchange between ordered and disordered states with or without substrate. Protecting osmolytes shift this equilibrium to favor the more ordered, folded state. However, a segment of BtuB that is C‐terminal to the Ton box that is not solvent exposed is insensitive to PEGs. Protecting osmolytes also modulate a conformational equilibrium in the Ton box of FecA, with larger molecular weight PEGs producing the largest shifts in the conformational free energy. These data indicate that solvent‐exposed regions of these transporters undergo conformational exchange and that regions of these transporters that are involved in protein–protein interactions sample multiple conformational substates. The sensitivity to solute provides an explanation for differences seen between two high‐resolution structures of BtuB, which each likely represent one conformation from a subset of states that are normally sampled by the protein. This work also illustrates how SDSL and osmolytes may be used to characterize and quantitate conformational equilibria in membrane proteins.     

Hydrogen exchange and the dynamic structure of proteins

Summary In native proteins, buried, labile protons undergo isotope exchange with solvent hydrogens, but the kinetics of exchange are markedly slower than in unfolded polypeptides. This indicates that, whereas buried protein atoms are shielded from solvent, the protein fluctuates around the time average structure and occasionally exposes buried sites to solvent. Generally, hydrogen exchange studies are designed to characterize the nature of the fluctuations between conformational substates, to monitor the shift in conformational equilibria among protein substates due to ligand binding or other factors, or to monitor the major cooperative denaturation transition. In this article, we review the recent reports of hydrogen exchange in proteins, focusing on recent advances in methodology, especially with regard to the implications of the results for the mechanism of hydrogen exchange in folded proteins.     

Heterogeneity and inaccuracy in protein structures solved by X-ray crystallography

Proteins are dynamic molecules, exhibiting structural heterogeneity in the form of anisotropic motion and discrete conformational substates, often of functional importance. In protein structure determination by X-ray crystallography, the observed diffraction pattern results from the scattering of X-rays by an ensemble of heterogeneous molecules, ordered and oriented by packing in a crystal lattice. The majority of proteins diffract to resolutions where heterogeneity is difficult to identify and model, and are therefore approximated by a single, average conformation with isotropic variance. Here we show that disregarding structural heterogeneity introduces degeneracy into the structure determination process, as many single, isotropic models exist that explain the diffraction data equally well. The large differences among these models imply that the accuracy of crystallographic structures has been widely overestimated. Further, it suggests that analyses that depend on small differences in the relative positions of atoms may be flawed.     

Conformational analysis of protein structures derived from NMR data

A study is presented of the conformational characteristics of NMR-derived protein structures in the Protein Data Bank compared to X-ray structures. Both ensemble and energy-minimized average structures are analyzed. We have addressed the problem using the methods developed for crystal structures by examining the distribution of ?, Ψ, and χ angles as indicators of global conformational irregularity. All these features in NMR structures occur to varying degrees in multiple conformational states. Some measures of local geometry are very tightly constrained by the methods used to generate the structure, e.g., proline ? angles, α-helix ?, Ψ angles, ω angles, and C_α chirality. The more lightly restrained torsion angles do show increasead clustering as the number of overall experimental observations increases. ?, Ψ, and χ₁ angle conformational heterogeneity is strongly correlated with accessibility but shows additional differences which reflect the differing number of observations possible in NMR for the various side chains (e.g., many for Trp, few for Ser). In general, we find that the core is defined to a notional resolution of 2.0 to 2.3 Å. Of real interest is the behavior of surface residues and in particular the side chains where multiple rotameric states in different structures can vary from 10% to 88%. Later generation structures show a much tighter definition which correlates with increasing use of J-coupling information, stereospecific assignments, and heteronumclear techniques. A suite of programs is being developed to address the special needs of NMR-derived structures which will take into account the existence of increased mobility in solution. © 1993 Wiley-Liss, Inc.     

Temperature dependent 2nd derivative absorbance spectroscopy of aromatic amino acids as a probe of protein dynamics

Proteins display a broad peak in 250–300 nm region of their UV spectrum containing multiple overlapping bands arising from the aromatic rings of phenylalanine, tyrosine, and tryptophan residues. Employing high resolution 2^nd derivative absorbance spectroscopy, these overlapping absorption bands can be highly resolved and therefore provide a very sensitive measure of changes in the local microenvironment of the aromatic side chains. This has traditionally been used to detect both subtle and dramatic (i.e., unfolding) conformational alterations of proteins. Herein, we show that plots of the temperature dependent 2^nd derivative peak positions of aromatic residues have measurable slopes before protein unfolding and that these slopes are sensitive to the dielectric properties of the surrounding microenvironment. We further demonstrate that these slopes correlate with hydration of the buried aromatic residues in protein cores and can therefore be used as qualitative probes of protein dynamics.     

Heterogeneity of Protein Substates Visualized by Spin-label EPR

The energy landscape of proteins is characterized by a hierarchy of substates, which give rise to conformational heterogeneity at low temperatures. In multiply spin-labeled membranous Na,K-ATPase, this heterogeneous population of conformations is manifest by strong inhomogeneous broadening of the electron paramagnetic resonance (EPR) line shapes and nonexponential spin-echo decays, which undergo a transition to homogeneous broadening and exponential relaxation at higher temperatures (previous study). In this study, we apply these EPR methods to small water-soluble proteins, of the type for which the existence of conformational substates is well established. Both α-helical and β-sheet aqueous proteins that are spin-labeled on a single cysteine residue display spin-echo decays with a single phase-memory time T_2M and conventional EPR line shapes with predominantly homogeneous broadening, over a broad range of temperatures from 77 K to ∼ 250 K or higher. Above ∼ 200 K, the residual inhomogeneous broadening is reduced almost to zero. In contrast, both the proteins and the spin label alone, when in a glycerol-water mixture below the glass transition, display heterogeneity in spin-echo phase-memory time and a stronger inhomogeneous broadening of the conventional line shapes, similar to multiply spin-labeled membranous Na,K-ATPase below 200 K. Above 200 K (or the glass transition), a single phase-memory time and predominantly homogeneous broadening are found in both spin-label systems. The results are discussed in terms of solvent-mediated protein transitions, the ability of single spin-label sites to detect conformational heterogeneity, and the desirability of exploring multiple sites for proteins with the size and complexity of the Na,K-ATPase.     

Heterogeneity of Protein Substates Visualized by Spin-label EPR

The energy landscape of proteins is characterized by a hierarchy of substates, which give rise to conformational heterogeneity at low temperatures. In multiply spin-labeled membranous Na,K-ATPase, this heterogeneous population of conformations is manifest by strong inhomogeneous broadening of the electron paramagnetic resonance (EPR) line shapes and nonexponential spin-echo decays, which undergo a transition to homogeneous broadening and exponential relaxation at higher temperatures (previous study). In this study, we apply these EPR methods to small water-soluble proteins, of the type for which the existence of conformational substates is well established. Both α-helical and β-sheet aqueous proteins that are spin-labeled on a single cysteine residue display spin-echo decays with a single phase-memory time T_2M and conventional EPR line shapes with predominantly homogeneous broadening, over a broad range of temperatures from 77 K to ∼ 250 K or higher. Above ∼ 200 K, the residual inhomogeneous broadening is reduced almost to zero. In contrast, both the proteins and the spin label alone, when in a glycerol-water mixture below the glass transition, display heterogeneity in spin-echo phase-memory time and a stronger inhomogeneous broadening of the conventional line shapes, similar to multiply spin-labeled membranous Na,K-ATPase below 200 K. Above 200 K (or the glass transition), a single phase-memory time and predominantly homogeneous broadening are found in both spin-label systems. The results are discussed in terms of solvent-mediated protein transitions, the ability of single spin-label sites to detect conformational heterogeneity, and the desirability of exploring multiple sites for proteins with the size and complexity of the Na,K-ATPase.     

Thermal movements of proteins: small-scale fluctuation and conformation substates

A comparison is made of the published data on different physical methods concerning small-scale (small-amplitude) thermal motions, conformers and conformational substates of myoglobin and some other proteins in crystal and solution. Because of the underestimation of the statical and dynamical rotational disorder of the crystal lattice, the X-ray diffraction method may somewhat overestimate the amplitude values of the intramolecular thermal motions, which leads to the underestimation of the extent of conformational rigidity of the protein molecule. In solution the rigid-body rotation of the macromolecules makes a significant contribution to the apparent amplitude of atomic displacement. The damping of these motions upon the increase of viscosity or freezing of solution, crystallization or intramolecular crosslinking can partially simulate the depression of intramolecular thermal motions. According to infrared spectroscopy data the crystallization of carboxymyoglobin as well as temperature decrease down to 5K are not accompanied by preferential stabilization of one of the several discrete conformers of the active site or by diminished distribution breadth of conformational substates within individual conformers. The structural differences between these discrete conformers and the more so between conformational substates are not sizable and are confined within a limit of the amplitude of thermal motion in the polymeric monocrystal including their paracrystalline regions. The dehydration of the hemoglobin and myoglobin films leads to the loss of conformational rigidity of the active site and to the distortion of its native conformation, i.e. to the inactivation of protein.     

Probing Substates in Sperm Whale Myoglobin Using High-Pressure Crystallography

Pressures in the 100 MPa range are known to have an enormous number of effects on the action of proteins, but straightforward means for determining the structural basis of these effects have been lacking. Here, crystallography has been used to probe effects of pressure on sperm whale myoglobin structure. A comparison of pressure effects with those seen at low pH suggests that structural changes under pressure are interpretable as a shift in the populations of conformational substates. Furthermore, a novel high-pressure protein crystal-cooling method has been used to show low-temperature metastability, providing an alternative to room temperature, beryllium pressure cell-based techniques. The change in protein structure due to pressure is not purely compressive and involves conformational changes important to protein activity. Correlation with low-pH structures suggests observed structural changes are associated with global conformational substates. Methods developed here open up a direct avenue for exploration of the effects of pressure on proteins.     

Beyond rotamers: a generative,probabilistic model of side chains in proteins

Background  

Accurately covering the conformational space of amino acid side chains is essential for important applications such as protein design, docking and high resolution structure prediction. Today, the most common way to capture this conformational space is through rotamer libraries - discrete collections of side chain conformations derived from experimentally determined protein structures. The discretization can be exploited to efficiently search the conformational space. However, discretizing this naturally continuous space comes at the cost of losing detailed information that is crucial for certain applications. For example, rigorously combining rotamers with physical force fields is associated with numerous problems.     

Conformational substates in enzyme mechanism: the 120 K structure of alpha-lytic protease at 1.5 A resolution.

Insight into the dynamic properties of alpha-lytic protease (alpha LP) has been obtained through the use of low-temperature X-ray crystallography and multiple-conformation refinement. Previous studies of alpha LP have shown that the residues around the active site are able to move significantly to accommodate substrates of different sizes. Here we show a link between the ability to accommodate ligands and the dynamics of the binding pocket. Although the structure of alpha LP at 120 K has B-factors with a uniformly low value of 4.8 A2 for the main chain, four regions stand out as having significantly higher B-factors. Because thermal motion should be suppressed at cryogenic temperatures, the high B-factors are interpreted as the result of trapped conformational substates. The active site residues that are perturbed during accommodation of different substrates are precisely those showing conformational substates, implying that substrate binding selects a subset of conformations from the ensemble of accessible states. To better characterize the precise nature of these substates, a protein model consisting of 16 structures has been refined and evaluated. The model reveals a number of features that could not be well-described by conventional B-factors: for example, 40% of the main-chain residue conformations are distributed asymmetrically or in discrete clusters. Furthermore, these data demonstrate an unexpected correlation between motions on either side of the binding pocket that we suggest is a consequence of "dynamic close packing." These results provide strong evidence for the role of protein dynamics in substrate binding and are consistent with the results of dynamic studies of ligand binding in myoglobin and ribonuclease A.     

Hydrogen bond dynamics in membrane protein function

Changes in inter-helical hydrogen bonding are associated with the conformational dynamics of membrane proteins. The function of the protein depends on the surrounding lipid membrane. Here we review through specific examples how dynamical hydrogen bonds can ensure an elegant and efficient mechanism of long-distance intra-protein and protein-lipid coupling, contributing to the stability of discrete protein conformational substates and to rapid propagation of structural perturbations. This article is part of a Special Issue entitled: Protein Folding in Membranes.     

Energy Landscape and Dynamics of Biomolecules Extended Abstract

Proteins are not isolated homogeneous systems. Each protein can exist in a very large number of conformations (conformational substates) that are characterized by an energy landscape. The main conformational motions, similar to the α and β fluctuations in glasses, are linked to fluctuations in the bulk solvent and the hydration shell.     

Frustration‐guided motion planning reveals conformational transitions in proteins

Proteins exist as conformational ensembles, exchanging between substates to perform their function. Advances in experimental techniques yield unprecedented access to structural snapshots of their conformational landscape. However, computationally modeling how proteins use collective motions to transition between substates is challenging owing to a rugged landscape and large energy barriers. Here, we present a new, robotics‐inspired motion planning procedure called dCC‐RRT that navigates the rugged landscape between substates by introducing dynamic, interatomic constraints to modulate frustration. The constraints balance non‐native contacts and flexibility, and instantaneously redirect the motion towards sterically favorable conformations. On a test set of eight proteins determined in two conformations separated by, on average, 7.5 Å root mean square deviation (RMSD), our pathways reduced the Cα atom RMSD to the goal conformation by 78%, outperforming peer methods. We then applied dCC‐RRT to examine how collective, small‐scale motions of four side‐chains in the active site of cyclophilin A propagate through the protein. dCC‐RRT uncovered a spatially contiguous network of residues linked by steric interactions and collective motion connecting the active site to a recently proposed, non‐canonical capsid binding site 25 Å away, rationalizing NMR and multi‐temperature crystallography experiments. In all, dCC‐RRT can reveal detailed, all‐atom molecular mechanisms for small and large amplitude motions. Source code and binaries are freely available at https://github.com/ExcitedStates/KGS/ .     

Conformational substates of calmodulin revealed by single-pair fluorescence resonance energy transfer: influence of solution conditions and oxidative modification

A calmodulin (CaM) mutant (T34,110C-CaM) doubly labeled with fluorescence probes AlexaFluor 488 and Texas Red in opposing domains (CaM-DA) has been used to examine conformational heterogeneity in CaM by single-pair fluorescence resonance energy transfer (spFRET). Burst-integrated FRET efficiencies of freely diffusing CaM-DA single molecules yielded distributions of distance between domains of CaM-DA. We recently reported distinct conformational substates of Ca(2+)-CaM-DA and apoCaM-DA, with peaks in the distance distributions centered at approximately 28 A, 34-38 A, and 55 A [Slaughter et al. (2004) J. Phys. Chem. B 108, 10388-10397]. In the present study, shifts in the amplitudes and center distances of the conformational substates were detected with variation in solution conditions. The amplitude of an extended conformation was observed to change as a function of Ca(2+) over a free Ca(2+) range that is consistent with binding to the high affinity, C-terminal Ca(2+) binding sites, suggesting the existence of communication between lobes of CaM. Lowering pH shifted the relative amplitudes of the conformations, with a marked increase in the presence of the compact conformations and an almost complete absence of the extended conformation. In addition, the single-molecule distance distribution of apoCaM-DA at reduced ionic strength was shifted to longer distance and showed evidence of an increase in conformational heterogeneity relative to apoCaM-DA at physiological ionic strength. Oxidation of methionine residues in CaM-DA produced a substantial increase in the amplitude of the extended conformation relative to the more compact conformation. The results are considered in light of a hypothesis that suggests that electrostatic interactions between charged amino acid side chains play an important role in determining the most stable CaM conformation under varying solution conditions.