Refinement of atomic models in high resolution EM reconstructions using Flex-EM and local assessment.

As the resolutions of Three Dimensional Electron Microscopic reconstructions of biological macromolecules are being improved, there is a need for better fitting and refinement methods at high resolutions and robust approaches for model assessment. Flex-EM/MODELLER has been used for flexible fitting of atomic models in intermediate-to-low resolution density maps of different biological systems. Here, we demonstrate the suitability of the method to successfully refine structures at higher resolutions (2.5-4.5Å) using both simulated and experimental data, including a newly processed map of Apo-GroEL. A hierarchical refinement protocol was adopted where the rigid body definitions are relaxed and atom displacement steps are reduced progressively at successive stages of refinement. For the assessment of local fit, we used the SMOC (segment-based Manders' overlap coefficient) score, while the model quality was checked using the Qmean score. Comparison of SMOC profiles at different stages of refinement helped in detecting regions that are poorly fitted. We also show how initial model errors can have significant impact on the goodness-of-fit. Finally, we discuss the implementation of Flex-EM in the CCP-EM software suite.

An integrated approach for thermal stabilization of a mesophilic adenylate kinase.

Thermally stable proteins are desirable for research and industrial purposes, but redesigning proteins for higher thermal stability can be challenging. A number of different techniques have been used to improve the thermal stability of proteins, but the extents of stability enhancement were sometimes unpredictable and not significant. Here, we systematically tested the effects of multiple stabilization techniques including a bioinformatic method and structure-guided mutagenesis on a single protein, thereby providing an integrated approach to protein thermal stabilization. Using a mesophilic adenylate kinase (AK) as a model, we identified stabilizing mutations based on various stabilization techniques, and generated a series of AK variants by introducing mutations both individually and collectively. The redesigned proteins displayed a range of increased thermal stabilities, the most stable of which was comparable to a naturally evolved thermophilic homologue with more than a 25° increase in its thermal denaturation midpoint. We also solved crystal structures of three representative variants including the most stable variant, to confirm the structural basis for their increased stabilities. These results provide a unique opportunity for systematically analyzing the effectiveness and additivity of various stabilization mechanisms, and they represent a useful approach for improving protein stability by integrating the reduction of local structural entropy and the optimization of global noncovalent interactions such as hydrophobic contact and ion pairs.

Solution stability and variability in a simple model of globular proteins.

It is well known among molecular biologists that proteins with a common ancestor and that perform the same function in similar organisms, can have rather different amino-acid sequences. Mutations have altered the amino-acid sequences without affecting the function. A simple model of a protein in which the interactions are encoded by sequences of bits is introduced, and used to study how mutations can change these bits, and hence the interactions, while maintaining the stability of the protein solution. This stability is a simple minimal requirement on our model proteins which mimics part of the requirement on a real protein to be functional. The properties of our model protein, such as its second virial coefficient, are found to vary significantly from one model protein to another. It is suggested that this may also be the case for real proteins in vivo.

Refined structure of the complex between guanylate kinase and its substrate GMP at 2.0 A resolution.

The crystal structure of guanylate kinase from Saccharomyces cerevisiae complexed with its substrate GMP has been refined at a resolution of 2.0 A. The final crystallographic R-factor is 17.3% in the resolution range 7.0 A to 2.0 A for all reflections of the 100% complete data set. The final model has standard geometry with root-mean-square deviations of 0.016 A in bond lengths and 3.0 in bond angles. It consists of all 186 amino acid residues, the N-terminal acetyl group, the substrate GMP, one sulfate ion and 174 water molecules. Guanylate kinase is structurally related to adenylate kinases and G-proteins with respect to its central beta-sheet with connecting helices and the giant anion hole that binds nucleoside triphosphates. These nucleotides are ATP and GTP for the kinases and GTP for the G-proteins. The chain segment binding the substrate GMP of guanylate kinase differs grossly from the respective part of the adenylate kinases; it has no counterpart in the G-proteins. The binding mode of GMP is described in detail. Probably, the observed structure represents one of several structurally quite different intermediate states of the catalytic cycle.

The refined structure of the complex between adenylate kinase from beef heart mitochondrial matrix and its substrate AMP at 1.85 A resolution.

The crystal structure of the complex between adenylate kinase from bovine mitochondrial matrix and its substrate AMP has been refined at 1.85 A resolution (1 A = 0.1 nm). Based on 42,519 independent reflections of better than 10 A resolution, a final R-factor of 18.9% was obtained with a model obeying standard geometry within 0.016 A in bond lengths and 3.2 degrees in bond angles. There are two enzyme: substrate complexes in the asymmetric unit, each consisting of 226 amino acid residues, one AMP and one sulfate ion. A superposition of the two full-length polypeptides revealed deviations that can be described as small relative movements of three domains. Best superpositions of individual domains yielded a residual overall root-mean-square deviation of 0.3 A for the backbone atoms and 0.5 A for the sidechains. The final model contains 381 solvent molecules in the asymmetric unit, 2 x 72 = 144 of which occupy corresponding positions in both complexes.