Not just your average structures

Evaluating the impact of X-ray damage on conformational heterogeneity in room-temperature (277 K) and cryo-cooled protein crystals

Cryo-cooling has been nearly universally adopted to mitigate X-ray damage and facilitate crystal handling in protein X-ray crystallography. However, cryo X-ray crystallographic data provide an incomplete window into the ensemble of conformations that is at the heart of protein function and energetics. Room-temperature (RT) X-ray crystallography provides accurate ensemble information, and recent developments allow conformational heterogeneity (the experimental manifestation of ensembles) to be extracted from single-crystal data. Nevertheless, high sensitivity to X-ray damage at RT raises concerns about data reliability. To systematically address this critical issue, increasingly X-ray-damaged high-resolution data sets (1.02–1.52 Å resolution) were obtained from single proteinase K, thaumatin and lysozyme crystals at RT (277 K). In each case a modest increase in conformational heterogeneity with X-ray damage was observed. Merging data with different extents of damage (as is typically carried out) had negligible effects on conformational heterogeneity until the overall diffraction intensity decayed to ∼70% of its initial value. These effects were compared with X-ray damage effects in cryo-cooled crystals by carrying out an analogous analysis of increasingly damaged proteinase K cryo data sets (0.9–1.16 Å resolution). X-ray damage-associated heterogeneity changes were found that were not observed at RT. This property renders it difficult to distinguish real from artefactual conformations and to determine the conformational response to changes in temperature. The ability to acquire reliable heterogeneity information from single crystals at RT, together with recent advances in RT data collection at accessible synchrotron beamlines, provides a strong motivation for the widespread adoption of RT X-ray crystallography to obtain conformational ensemble information.

Instrumentation and experimental procedures for robust collection of X-ray diffraction data from protein crystals across physiological temperatures

Traditional X-ray diffraction data collected at cryo-temperatures have delivered invaluable insights into the three-dimensional structures of proteins, providing the backbone of structure-function studies. While cryo-cooling mitigates radiation damage, cryo-temperatures can alter protein conformational ensembles and solvent structure. Furthermore, conformational ensembles underlie protein function and energetics, and recent advances in room-temperature X-ray crystallography have delivered conformational heterogeneity information that can be directly related to biological function. Given this capability, the next challenge is to develop a robust and broadly applicable method to collect single-crystal X-ray diffraction data at and above room temperature. This challenge is addressed herein. The approach described provides complete diffraction data sets with total collection times as short as ∼5 s from single protein crystals, dramatically increasing the quantity of data that can be collected within allocated synchrotron beam time. Its applicability was demonstrated by collecting 1.09-1.54 Å resolution data over a temperature range of 293-363 K for proteinase K, thaumatin and lysozyme crystals at BL14-1 at the Stanford Synchrotron Radiation Lightsource. The analyses presented here indicate that the diffraction data are of high quality and do not suffer from excessive dehydration or radiation damage.

Validating Molecular Dynamics Simulations against Experimental Observables in Light of Underlying Conformational Ensembles

Far from the static, idealized conformations deposited into structural databases, proteins are highly dynamic molecules that undergo conformational changes on temporal and spatial scales that may span several orders of magnitude. These conformational changes, often intimately connected to the functional roles that proteins play, may be obscured by traditional biophysical techniques. Over the past 40 years, molecular dynamics (MD) simulations have complemented these techniques by providing the "hidden" atomistic details that underlie protein dynamics. However, there are limitations of the degree to which molecular simulations accurately and quantitatively describe protein motions. Here we show that although four molecular dynamics simulation packages (AMBER, GROMACS, NAMD, and ilmm) reproduced a variety of experimental observables for two different proteins (engrailed homeodomain and RNase H) equally well overall at room temperature, there were subtle differences in the underlying conformational distributions and the extent of conformational sampling obtained. This leads to ambiguity about which results are correct, as experiment cannot always provide the necessary detailed information to distinguish between the underlying conformational ensembles. However, the results with different packages diverged more when considering larger amplitude motion, for example, the thermal unfolding process and conformational states sampled, with some packages failing to allow the protein to unfold at high temperature or providing results at odds with experiment. While most differences between MD simulations performed with different packages are attributed to the force fields themselves, there are many other factors that influence the outcome, including the water model, algorithms that constrain motion, how atomic interactions are handled, and the simulation ensemble employed. Here four different MD packages were tested each using best practices as established by the developers, utilizing three different protein force fields and three different water models. Differences between the simulated protein behavior using two different packages but the same force field, as well as two different packages with different force fields but the same water models and approaches to restraining motion, show how other factors can influence the behavior, and it is incorrect to place all the blame for deviations and errors on force fields or to expect improvements in force fields alone to solve such problems.

Role of Protein Dynamics in Reaction Rate Enhancement by Enzymes

An integrated view of protein structure, dynamics, and function is emerging, where proteins are considered as dynamically active assemblies and internal motions are closely linked to function such as enzyme catalysis. Further, the motion of solvent bound to external regions of protein impacts internal motions and, therefore, protein function. Recently, we discovered a network of protein vibrations in enzyme cyclophilin A, coupled to its catalytic activity of peptidyl-prolyl cis-trans isomerization. Detailed studies suggest that this network, extending from surface regions to active site, is a conserved part of enzyme structure and has a role in promoting catalysis. In this report, theoretical investigations of concerted conformational fluctuations occurring on microsecond and longer time scales within the discovered network are presented. Using a new technique, kinetic energy was added to protein vibrational modes corresponding to conformational fluctuations in the network. The results reveal that protein dynamics promotes catalysis by altering transition state barrier crossing behavior of reaction trajectories. An increase in transmission coefficient and number of productive trajectories with increasing amounts of kinetic energy in vibrational modes is observed. Variations in active site enzyme-substrate interactions near transition state are found to be correlated with barrier recrossings. Simulations also showed that energy transferred from first solvation shell to surface residues impacts catalysis through network fluctuations. The detailed characterization of network presented here indicates that protein dynamics plays a role in rate enhancement by enzymes. Therefore, coupled networks in enzymes have wide implications in understanding allostericity and cooperative effects, as well as protein engineering and drug design.

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.

Triggering loops and enzyme function: identification of loops that trigger and modulate movements

Enzyme function often involves a conformational change. There is a general agreement that loops play a vital role in correctly positioning the catalytically important residues. Nevertheless, predicting the functional loops and most importantly their role in enzyme function remains a difficult task. A major reason for this difficulty is that loops that undergo conformational change are frequently not well conserved in their primary sequence. beta1,4-Galactosyltransferase is one such enzyme. There, the amino acid sequence of a long loop that undergoes a large conformational change upon substrate binding is not well conserved. Our molecular dynamics simulations show that the large conformational change in the long loop is brought about by a second, interacting loop. Interestingly, while the structural change of the second loop is much smaller than that of the long loop, its sequence (particularly glycine residues) is highly conserved. We further examine the generality of the proposition that there are loops that trigger movements but nevertheless show little or no structural changes in crystals. We focus on two other enzymes, enolase and lipase. We chose these enzymes, since they too undergo conformational change upon ligand binding, however, they have different folds and different functions. Through multiple sets of simulations we show that the conformational change of the functional loop(s) is brought about through communication of flexibility by triggering loops that have several glycine residues. We further propose that similar to the conservation of common favorable fold types and structural motifs, evolution has also conserved common "skillful" mechanisms. Mechanisms may be conserved across different folds, sequences and functions, with adaptation to specific enzymatic roles.

Structure determination of membrane-associated proteins from nuclear magnetic resonance data

This Review covers the delineation and optimization of protein-lipid systems for study using solution-state NMR spectroscopy. The first half presents the necessary background for a membrane protein biochemist to initiate collaboration with an NMR spectroscopist. The second half provides guidelines for the spectroscopist on data collection, analysis for obtaining conformational information, and structure generation and assessment. Although the emphasis is on the study of peptides in detergent micelles, methods are outlined for larger membrane-associated proteins and for use of other solubilizing agents.

Elbow flexibility and ligand-induced domain rearrangements in antibody Fab NC6.8: large effects of a small hapten

Four 700-ps molecular dynamics simulations were carried out to analyze the structural dynamics of the antigen-binding antibody fragment NC6.8, which is known to exhibit large structural changes upon complexation. The first simulation was started from the x-ray structure of the uncomplexed Fab and produced trajectory averages that closely match the crystallographic results. It allowed assessment of the flexibility of the Fab, revealing an elbow motion of the variable domains with respect to the constant domains. The second simulation was started from the uncomplexed x-ray structure after insertion of the ligand into the binding site. This perturbation resulted in a significantly altered trajectory, with quaternary structural changes corresponding in many aspects to the experimental differences between complexed and uncomplexed state. The observed trend toward a smaller elbow angle and a higher flexion of the H-chain could also be seen in the third simulation, which was started from the x-ray structure of the complex. The changes were revealed to be a clear consequence of the complexation with the ligand because in the fourth simulation (started from the experimental complex structure after removal of the hapten) the Fab remained close to its initial structure. Analyses of the quaternary structure and the binding site of Fab NC6.8 are presented for all four simulations, and possible interpretations are discussed.

Directed mutagenesis studies of the metal binding site at the subunit interface of Escherichia coli inorganic pyrophosphatase

Recent crystallographic studies on Escherichia coli inorganic pyrophosphatase (E-PPase) have identified three Mg2+ ions/enzyme hexamer in water-filled cavities formed by Asn24, Ala25, and Asp26 at the trimer-trimer interface (Kankare, J., Salminen, T., Lahti, R., Cooperman, B., Baykov, A. A., and Goldman, A. (1996) Biochemistry 35, 4670-4677). Here we show that D26S and D26N substitutions decrease the stoichiometry of tight Mg2+ binding to E-PPase by approximately 0.5 mol/mol monomer and increase hexamer stability in acidic medium. Mg2+ markedly decelerates the dissociation of enzyme hexamer into trimers at pH 5.0 and accelerates hexamer formation from trimers at pH 7.2 with wild type E-PPase and the N24D variant, in contrast to the D26S and D26N variants, when little or no effect is seen. The catalytic parameters describing the dependences of enzyme activity on substrate and Mg2+ concentrations are of the same magnitude for wild type E-PPase and the three variants. The affinity of the intertrimer site for Mg2+ at pH 7.2 is intermediate between those of two Mg2+ binding sites found in the E-PPase active site. It is concluded that the metal ion binding site found at the trimer-trimer interface of E-PPase is a high affinity site whose occupancy by Mg2+ greatly stabilizes the enzyme hexamer but has little effect on catalysis.

Ligand-induced domain movement in an antibody Fab: molecular dynamics studies confirm the unique domain movement observed experimentally for Fab NC6.8 upon complexation and reveal its segmental flexibility

Two molecular dynamics simulations were carried out for the antibody Fab NC6.8, both with and without the guanidinium sweetener ligand NC174, in order to assess the segmental flexibility as well as the conformational changes upon ligand binding. Trajectory analyses of the simulation of the uncomplexed Fab suggest low-amplitude motions of the Ig domains with respect to each other, most clearly reflected by a periodic alteration of the elbow angle within a range of 11 degrees. Upon insertion of the hapten into the binding site, the quaternary structure of the Fab exhibits considerable rearrangements: the elbow angle changes by almost 30 degrees, the light chain is elongated and the heavy chain becomes more flexed. Comparison with experiment reveals some interesting agreements with X-ray crystallographic results published previously.

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.

Efficient exchange of the primary quinone acceptor Q(A) in isolated reaction centers of Rhodopseudomonas viridis

A key step in the conversion of solar energy into chemical energy by photosynthetic reaction centers (RCs) occurs at the level of the two quinones, Q(A) and Q(B), where electron transfer couples to proton transfer. A great deal of our understanding of the mechanisms of these coupled reactions relies on the seminal work of Okamura et al. [Okamura, M. Y., Isaacson, R. A., & Feher, G. (1975) Proc. Natl. Acad. Sci. USA 88, 3491-3495], who were able to extract with detergents the firmly bound ubiquinone Q(A) from the RC of Rhodobacter sphaeroides and reconstitute the site with extraneous quinones. Up to now a comparable protocol was lacking for the RC of Rhodopseudomonas viridis despite the fact that its Q(A) site, which contains 2-methyl-3-nonaprenyl-1,4-naphthoquinone (menaquinone-9), has provided the best x-ray structure available. Fourier transform infrared difference spectroscopy, together with the use of isotopically labeled quinones, can probe the interaction of Q(A) with the RC protein. We establish that a simple incubation procedure of isolated RCs of Rp. viridis with an excess of extraneous quinone allows the menaquinone-9 in the Q(A) site to be almost quantitatively replaced either by vitamin K(1), a close analogue of menaquinone-9, or by ubiquinone. To our knowledge, this is the first report of quinone exchange in bacterial photosynthesis. The Fourier transform infrared data on the quinone and semiquinone vibrations show a close similarity in the bonding interactions of vitamin K(1) with the protein at the Q(A) site of Rp. viridis and Rb. sphaeroides, whereas for ubiquinone these interactions are significantly different. The results are interpreted in terms of slightly inequivalent quinone-protein interactions by comparison with the crystallographic data available for the Q(A) site of the two RCs.

Mutant RecA proteins which form hexamer-sized oligomers

We have analyzed the oligomeric properties of a number of mutant RecA proteins containing single amino acid substitutions within one region of the subunit interface. In contrast to wild-type RecA, which forms a heterogeneous population of different-sized oligomers, we find that many of these mutant proteins exist in a more homogeneous oligomeric form, which approximates to the size of a RecA hexamer. Some of these mutants have a significant level of activity in vivo for recombinational DNA repair and thus represent the first mutant RecA proteins identified which retain activity yet can exist in a discrete oligomeric state as free protein.