A Correspondence Between Solution-State Dynamics of an Individual Protein and the Sequence and Conformational Diversity of its Family

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.     

Protein-inhibitor flexible docking by a multicanonical sampling: native complex structure with the lowest free energy and a free-energy barrier distinguishing the native complex from the others

Flexible docking between a protein (lysozyme) and an inhibitor (tri-N-acetyl-D-glucosamine, tri-NAG) was carried out by an enhanced conformational sampling method, multicanonical molecular dynamics simulation. We used a flexible all-atom model to express lysozyme, tri-NAG, and water molecules surrounding the two bio-molecules. The advantages of this sampling method are as follows: the conformation of system is widely sampled without trapping at energy minima, a thermally equilibrated conformational ensemble at an arbitrary temperature can be reconstructed from the simulation trajectory, and the thermodynamic weight can be assigned to each sampled conformation. During the simulation, exchanges between the binding and free (i.e., unbinding) states of the protein and the inhibitor were repeatedly observed. The conformational ensemble reconstructed at 300 K involved various conformational clusters. The main outcome of the current study is that the most populated conformational cluster (i.e., the cluster of the lowest free energy) was assigned to the native complex structure (i.e., the X-ray complex structure). The simulation also produced non-native complex structures, where the protein and the inhibitor bound with different modes from that of the native complex structure, as well as the unbinding structures. A free-energy barrier (i.e., activation free energy) was clearly detected between the native complex structures and the other structures. The thermal fluctuations of tri-NAG in the lowest free-energy complex correlated well with the X-ray B-factors of tri-NAG in the X-ray complex structure. The existence of the free-energy barrier ensures that the lowest free-energy structure can be discriminated naturally from the other structures. In other words, the multicanonical molecular dynamics simulation can predict the native complex structure without any empirical objective function. The current study also manifested that the flexible all-atom model and the physico-chemically defined atomic-level force field can reproduce the native complex structure. A drawback of the current method is that it requires a time consuming computation due to the exhaustive conformational sampling. We discussed a possibility for combining the current method with conventional docking methods.     

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

Solution conformation of a peptide fragment representing a proposed RNA-binding site of a viral coat protein studied by two-dimensional NMR

The first 25 amino acids of the coat protein of cowpea chlorotic mottle virus are essential for binding the encapsidated RNA. Although an alpha-helical conformation has been predicted for this highly positively charged N-terminal region [Argos, P. (1981) Virology 110, 55-62; Vriend, G., Verduin, B. J. M., & Hemminga, M. A. (1986) J. Mol. Biol. 191, 453-460], no experimental evidence for this conformation has been presented so far. In this study, two-dimensional proton NMR experiments were performed on a chemically synthesized pentacosapeptide containing the first 25 amino acids of this coat protein [Ten Kortenaar, P. B. W., Krüse, J., Hemminga, M. A., & Tesser, G. I. (1986) Int. J. Pept. Protein Res. 27, 401-413]. All resonances could be assigned by a combined use of two-dimensional correlated spectroscopy and nuclear Overhauser enhancement spectroscopy carried out at four different temperatures. Various NMR parameters indicate the presence of a conformational ensemble consisting of helical structures rapidly converting into more extended states. Differences in chemical shifts and nuclear Overhauser effects indicate that lowering the temperature induces a shift of the dynamic equilibrium toward more helical structures. At 10 degrees C, a perceptible fraction of the conformational ensemble consists of structures with an alpha-helical conformation between residues 9 and 17, likely starting with a turnlike structure around Thr9 and Arg10. Both the conformation and the position of this helical region agree well with the secondary structure predictions mentioned above.     

OLDERADO: on-line database of ensemble representatives and domains. On Line Database of Ensemble Representatives And DOmains.

In cases where the structure of a single protein is represented by an ensemble of conformations, there is often a need to determine the common features and to choose a "representative" conformation. This occurs, for example, with structures determined by NMR spectroscopy, analysis of the trajectory from a molecular dynamics simulation, or an ensemble of structures produced by comparative modeling. We reported previously automatic methods for (1) defining the atoms with low spatial variance across an ensemble (i.e., the "core" atoms) and the domains in which these atoms lie, and (2) clustering an ensemble into conformationally related subfamilies. To extend the utility of these methods, we have developed a freely available server on the World Wide Web at http:/(/)neon.chem.le.ac.uk/olderado/. This (1) contains an automatically generated database of representative structures, core atoms, and domains determined for 449 ensembles of NMR-derived protein structures in the Protein Data Bank (PDB) in May 1997, and (2) allows the user to upload a PDB-formatted file containing the coordinates of an ensemble of structures. The server returns in real time: (1) information on the residues constituting domains: (2) the structures that constitute each conformational subfamily; and (3) an interactive java-based three-dimensional viewer to visualise the domains and clusters. Such information is useful, for example, when selecting conformations to be used in comparative modeling and when choosing parts of structures to be used in molecular replacement. Here we describe the OLDERADO server.     

Ab initio simulation of a 57-residue protein in explicit solvent reproduces the native conformation in the lowest free-energy cluster

An enhanced conformational sampling method, multicanonical molecular dynamics (McMD), was applied to the ab intio folding of the 57-residue first repeat of human glutamyl- prolyl-tRNA synthetase (EPRS-R1) in explicit solvent. The simulation started from a fully extended structure of EPRS-R1 and did not utilize prior structural knowledge. A canonical ensemble, which is a conformational ensemble thermodynamically probable at an arbitrary temperature, was constructed by reweighting the sampled structures. Conformational clusters were obtained from the canonical ensemble at 300 K, and the largest cluster (i.e., the lowest free-energy cluster), which contained 34% of the structures in the ensemble, was characterized by the highest similarity to the NMR structure relative to all alternative clusters. This lowest free-energy cluster included native-like structures composed of two anti-parallel α-helices. The canonical ensemble at 300 K also showed that a short Gly-containing segment, which adopts an α-helix in the native structure, has a tendency to be structurally disordered. Atomic-level analyses demonstrated clearly that inter-residue hydrophobic interactions drive the helix formation of the Gly-containing segment, and that increasing the hydrophobic contacts accompanies exclusion of water molecules from the vicinity of this segment. This study has shown, for the first time, that the free-energy landscape of a structurally well-ordered protein of about 60 residues is obtainable with an all atom model in explicit water without prior structural knowledge.     

Structural Dynamics and Conformational Equilibria of SERCA Regulatory Proteins in Membranes by Solid-State NMR Restrained Simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.     

Structural Dynamics and Conformational Equilibria of SERCA Regulatory Proteins in Membranes by Solid-State NMR Restrained Simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.     

A switch I mutant of Cdc42 exhibits less conformational freedom

Cdc42 is a Ras-related small G-protein and functions as a molecular switch in signal transduction pathways linked with cell growth and differentiation. It is controlled by cycling between GTP-bound (active) and GDP-bound (inactive) forms. Nucleotide binding and hydrolysis are modulated by interactions with effectors and/or regulatory proteins. These interactions are centralized in two relatively flexible "Switch" regions as characterized by internal dynamics on multiple time scales [Loh, A. P., et al. (2001) Biochemistry 40, 4590-4600], and this flexibility may be essential for protein interactions. In the Switch I region, Thr(35) seems to be critical for function, as it is completely invariant in Ras-related proteins. To investigate the importance of conformational flexibility in Switch I of Cdc42, we mutated threonine to alanine, determined the solution structure, and characterized the backbone dynamics of the single-point mutant protein, Cdc42(T35A). Backbone dynamics data suggest that the mutation changes the time scale of the internal motions of several residues, with several resonances not being discernible in wild-type Cdc42 [Adams, P. D., and Oswald, R. E. (2007) Biomol. NMR Assignments 1, 225-227]. The mutation does not appear to affect the thermal stability of Cdc42, and chymotrypsin digestion data further suggest that changes in the conformational flexibility of Switch I slow proteolytic cleavage relative to that of the wild type. In vitro binding assays show less binding of Cdc42(T35A), relative to that of wild type, to a GTPase binding protein that inhibits GTP hydrolysis in Cdc42. These results suggest that the mutation of T(35) leads to the loss of conformational freedom in Switch I that could affect effector-regulatory protein interactions.     

Virtual and solution conformations of oligosaccharides

The possibility that observed nuclear Overhauser enhancements and bulk longitudinal relaxation times, parameters measured by 1H NMR and often employed in determining the preferred solution conformation of biologically important molecules, are the result of averaging over many conformational states is quantitatively evaluated. Of particular interest was to ascertain whether certain 1H NMR determined conformations are "virtual" in nature; i.e., the fraction of the population of molecules actually found at any time within the subset of conformational space defined as the "solution conformation" is vanishingly small. A statistical mechanics approach was utilized to calculate an ensemble average relaxation matrix from which (NOE)'s and (T1)'s are calculated. Model glycosidic linkages in four oligosaccharides were studied. The solution conformation at any glycosidic linkage is properly represented by a normalized, Boltzmann distribution of conformers generated from an appropriate potential energy surface. The nature of the resultant population distributions is such that 50% of the molecular population is found within 1% of available microstates, while 99% of the molecular population occupies about 10% of the ensemble microstates, a number roughly equal to that sterically allowed. From this analysis we conclude that in many cases quantitative interpretation of NMR relaxation data, which attempts to define a single set of allowable torsion angle values consistent with the observed data, will lead to solution conformations that are either virtual or reflect torsion angle values possessed by a minority of the molecular population. On the other hand, calculation of ensemble average NMR relaxation data yields values in agreement with experimental results. Observed values of NMR relaxation data are the result of the complex interdependence of the population distribution and NOE (or T1) surfaces in conformational space. In conformational analyses, NMR data can therefore be used to test different population distributions calculated from empirical potential energy functions.     

The N(0)-binding region of the vesicular stomatitis virus phosphoprotein is globally disordered but contains transient α-helices

The phosphoprotein (P) of vesicular stomatitis virus (VSV) interacts with nascent nucleoprotein (N), forming the N(0)-P complex that is indispensable for the correct encapsidation of newly synthesized viral RNA genome. In this complex, the N-terminal region (P(NTR)) of P prevents N from binding to cellular RNA and keeps it available for encapsidating viral RNA genomes. Here, using nuclear magnetic resonance (NMR) spectroscopy and small-angle X-ray scattering (SAXS), we show that an isolated peptide corresponding to the 60 first N-terminal residues of VSV P (P(60)) and encompassing P(NTR) has overall molecular dimensions and a dynamic behavior characteristic of a disordered protein but transiently populates conformers containing α-helices. The modeling of P(60) as a conformational ensemble by the ensemble optimization method using SAXS data correctly reproduces the α-helical content detected by NMR spectroscopy and suggests the coexistence of subensembles of different compactness. The populations and overall dimensions of these subensembles are affected by the addition of stabilizing (1M trimethylamine-N-oxide) or destabilizing (6M guanidinium chloride) cosolvents. Our results are interpreted in the context of a scenario whereby VSV P(NTR) constitutes a molecular recognition element undergoing a disorder-to-order transition upon binding to its partner when forming the N(0)-P complex.     

Slow conformational dynamics in the hamster prion protein

Although the mechanism of the conformational conversion from the cellular (PrP(C)) to the scrapie (PrP(Sc)) form of animal prion proteins has yet to be elucidated, evidence is accumulating that may provide insight into the conversion process at atomic resolution. Here we show critical aspects of the slow fluctuation dynamics of the recombinant hamster prion protein, rPrP(90-231), based on NMR relaxation analysis using Carr-Purcell-Meiboom-Gill (CPMG) experiments, and compare them in detail with results from high-pressure NMR. Residues exhibiting slow fluctuations on the time scale of microseconds to milliseconds are mainly localized on helices B and C (172-193 and 200-227), which include locally disordered regions in an intermediate conformer, PrP*, identified previously by high-pressure NMR [Kuwata, K., et al., (2002) Biochemistry 41, 12277-12283]. Moreover, chemical shift differences between two putative exchanging conformers obtained by the CPMG relaxation analysis and the linear component of the pressure-induced chemical shift changes are reasonably correlated at individual residue sites. These observations suggest that both the CMPG relaxation and the pressure shifts reflect slow conformational fluctuations and that these slow motions in PrP(C) are related to the trajectories leading to the transition to PrP*.     

Analysis of 31P MAS NMR spectra and transversal relaxation of bacteriophage M13 and tobacco mosaic virus.

Phosphorus magic angle spinning nuclear magnetic resonance (NMR) spectra and transversal relaxation of M13 and TMV are analyzed by use of a model, which includes both local backbone motions of the encapsulated nucleic acid molecules and overall rotational diffusion of the rod-shaped virions about their length axis. Backbone motions influence the sideband intensities by causing a fast restricted reorientation of the phosphodiesters. To evaluate their influence on the observed sideband patterns, we extend the model that we used previously to analyze nonspinning 31P NMR lineshapes (Magusin, P.C.M.M., and M. A. Hemminga. 1993a. Biophys. J. 64:1861-1868) to magic angle spinning NMR experiments. Backbone motions also influence the conformation of the phosphodiesters, causing conformational averaging of the isotropic chemical shift, which offers a possible explanation for the various linewidths of the centerband and the sidebands observed for M13 gels under various conditions. The change of the experimental lineshape of M13 as a function of temperature and hydration is interpreted in terms of fast restricted fluctuation of the dihedral angles between the POC and the OCH planes on both sides of the 31P nucleus in the nucleic acid backbone. Backbone motions also seem to be the main cause of transversal relaxation measured at spinning rates of 4 kHz or higher. At spinning rates less than 2 kHz, transversal relaxation is significantly faster. This effect is assigned to slow, overall rotation of the rod-shaped M13 phage about its length axis. Equations are derived to simulate the observed dependence of T2e on the spinning rate.     

Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free‐energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography‐based restraints into the MD protocol. Because these restraints are aimed at the ensemble‐average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native‐like fashion as dictated by the original force field. To validate this approach, we have used the data from solid‐state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well‐established model protein, ubiquitin. The ensemble‐restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid‐state chemical shifts, and backbone order parameters. The predictions for ¹⁵N relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high‐resolution crystallographic data can be viewed as protein‐specific empirical corrections to the standard force fields.     

Picometer-scale conformational heterogeneity separates functional from nonfunctional states of a photoreceptor protein

Protein structural fluctuations occur over a wide spatial scale, ranging from minute, picometer-scale displacements, to large, interdomain motions and partial unfolding. While large-scale protein structural changes and their effects on protein function have been the focus of much recent attention, small-scale fluctuations have been less well studied, and are generally assumed to have proportionally smaller effects. Here we use the bacterial photoreceptor photoactive yellow protein (PYP) to test if subtle structural changes do, indeed, imply equally subtle functional effects. We flash froze crystals of PYP to trap the protein's conformational ensemble, and probed the molecules in this ensemble for their ability to facilitate PYP's biological function (i.e., light-driven isomerization of its chromophore). Our results indicate that the apparently homogeneous structural state observed in a 0.82 A crystal structure in fact comprises an ensemble of conformational states, in which subpopulations with nearly identical structures display dramatically different functional properties.     

Molecular dynamics simulation of dimeric and monomeric forms of human prion protein: insight into dynamics and properties

A central theme in prion protein research is the detection of the process that underlies the conformational transition from the normal cellular prion form (PrP(C)) to its pathogenic isoform (PrP(Sc)). Although the three-dimensional structures of monomeric and dimeric human prion protein (HuPrP) have been revealed by NMR spectroscopy and x-ray crystallography, the process underlying the conformational change from PrP(C) to PrP(Sc) and the dynamics and functions of PrP(C) remain unknown. The dimeric form is thought to play an important role in the conformational transition. In this study, we performed molecular dynamics (MD) simulations on monomeric and dimeric HuPrP at 300 K and 500 K for 10 ns to investigate the differences in the properties of the monomer and the dimer from the perspective of dynamic and structural behaviors. Simulations were also undertaken with Asp178Asn and acidic pH, which is known as a disease-associated factor. Our results indicate that the dynamics of the dimer and monomer were similar (e.g., denaturation of helices and elongation of the beta-sheet). However, additional secondary structure elements formed in the dimer might result in showing the differences in dynamics and properties between the monomer and dimer (e.g., the greater retention of dimeric than monomeric tertiary structure).     

Determination of the secondary structure and molecular topology of interleukin-1 beta by use of two- and three-dimensional heteronuclear 15N-1H NMR spectroscopy

A study of the regular secondary structure elements of recombinant human interleukin-1 beta has been carried out using NMR spectroscopy. Using a randomly 15N labeled sample, a number of heteronuclear three- and two-dimensional NMR experiments have been performed, which have enabled a complete analysis of short-, medium-, and long-range NOEs between protons of the polypeptide backbone, based on the sequence-specific resonance assignments that have been reported previously [Driscoll, P. C., Clore, G. M., Marion, D., Wingfield, P. T., & Gronenborn, A. M. (1990) Biochemistry 29, 3542-3556]. In addition, accurate measurements of a large number of 3JHN alpha coupling constants have been carried out by two-dimensional heteronuclear multiple-quantum-coherence-J spectroscopy. Amide NH solvent exchange rates have been measured by following the time dependence of the 15N-1H correlation spectrum of interleukin-1 beta on dissolving the protein in D2O solution. Analysis of these data indicate that the structure of interleukin-1 beta consists of 12 extended beta-strands aligned in a single extended network of antiparallel beta-sheet structure that in part folds into a skewed six-stranded beta-barrel. In the overall structure the beta-strands are connected by tight turns, short loops, and long loops in a manner that displays approximate pseudo-three-fold symmetry. The secondary structure analysis is discussed in the light of the unrefined X-ray structure of interleukin-1 beta at 3-A resolution [Priestle, J. P., Sch?r, H.-P., & Grütter, M. G. (1988) EMBO J. 7, 339-343], as well as biological activity data. Discernible differences between the two studies are highlighted. Finally, we have discovered conformational heterogeneity in the structure of interleukin-1 beta, which is characterized by an exchange rate that is slow on the NMR chemical shift time scale.     

Transient interaction with nanoparticles "freezes" a protein in an ensemble of metastable near-native conformations

It is well-known that adsorption of proteins on interfaces often induces substantial alterations of the protein structure. However, very little is known about whether these conformational changes have any consequence for the protein conformation after desorption from the interface. To investigate this matter, we have selected a protein-particle system in which the enzyme human carbonic anhydrase I (HCAI) alternates between the adsorbed and free state upon interaction with the silica nanoparticles. High-resolution NMR analysis of the protein with the particles present in the sample shows a spectrum that indicates a molten globular-like structure. Removal of particles results in refolding of virtually all HCAI molecules to a fully active form. However, the two-dimensional NMR analysis shows that refolding does not result in a single well-defined protein structure but rather provides an ensemble of protein molecules with near-native conformations. A detailed comparative chemical shift analysis of 108 amide signals in (1)H-(15)N HSQC spectra of native and desorbed HCAI reveals that the most profound effects are located at beta-strands in the center of the molecule. The observation of very slow H-D exchange in the central beta-strands of HCAI [Kjellsson, A., Sethson, I., and Jonsson, B. H. (2003) Biochemistry 42, 363-374] in conjunction with our results indicates that the kinetic barriers for conformational rearrangements in the central core of the protein are low in the presence of nanoparticles but are very high under native conditions.