YUP.SCX: coaxing atomic models into medium resolution electron density maps

The structures of large macromolecular complexes in different functional states can be determined by cryo-electron microscopy, which yields electron density maps of low to intermediate resolutions. The maps can be combined with high-resolution atomic structures of components of the complex, to produce a model for the complex that is more accurate than the formal resolution of the map. To this end, methods have been developed to dock atomic models into density maps rigidly or flexibly, and to refine a docked model so as to optimize the fit of the atomic model into the map. We have developed a new refinement method called YUP.SCX. The electron density map is converted into a component of the potential energy function to which terms for stereochemical restraints and volume exclusion are added. The potential energy function is then minimized (using simulated annealing) to yield a stereochemically-restrained atomic structure that fits into the electron density map optimally. We used this procedure to construct an atomic model of the 70S ribosome in the pre-accommodation state. Although some atoms are displaced by as much as 33 Å, they divide themselves into nearly rigid fragments along natural boundaries with smooth transitions between the fragments.     

Computational prediction of atomic structures of helical membrane proteins aided by EM maps

Integral membrane proteins pose a major challenge for protein-structure prediction because only approximately 100 high-resolution structures are available currently, thereby impeding the development of rules or empirical potentials to predict the packing of transmembrane alpha-helices. However, when an intermediate-resolution electron microscopy (EM) map is available, it can be used to provide restraints which, in combination with a suitable computational protocol, make structure prediction feasible. In this work we present such a protocol, which proceeds in three stages: 1), generation of an ensemble of alpha-helices by flexible fitting into each of the density rods in the low-resolution EM map, spanning a range of rotational angles around the main helical axes and translational shifts along the density rods; 2), fast optimization of side chains and scoring of the resulting conformations; and 3), refinement of the lowest-scoring conformations with internal coordinate mechanics, by optimizing the van der Waals, electrostatics, hydrogen bonding, torsional, and solvation energy contributions. In addition, our method implements a penalty term through a so-called tethering map, derived from the EM map, which restrains the positions of the alpha-helices. The protocol was validated on three test cases: GpA, KcsA, and MscL.     

Cryo‐EM model validation using independent map reconstructions

An increasing number of cryo‐electron microscopy (cryo‐EM) density maps are being generated with suitable resolution to trace the protein backbone and guide sidechain placement. Generating and evaluating atomic models based on such maps would be greatly facilitated by independent validation metrics for assessing the fit of the models to the data. We describe such a metric based on the fit of atomic models with independent test maps from single particle reconstructions not used in model refinement. The metric provides a means to determine the proper balance between the fit to the density and model energy and stereochemistry during refinement, and is likely to be useful in determining values of model building and refinement metaparameters quite generally.     

A disulfide conjugate between anti-tetanus antibodies and HIV (37-72)Tat neutralizes tetanus toxin inside chromaffin cells

The structure of the V1 ATPase from the tobacco hornworm Manduca sexta has been determined from electron micrographs of isolated, negatively stained specimens. The resulting images clearly show a pseudohexagonal arrangement of six equal-sized protein densities, presumably representing the three copies each of subunits A and B, which comprise the headpiece of the enzyme. A seventh density could be observed either centrally or asymmetrically to the hexamer. The maximum diameter of the V1 complex in the hexagonal projection is 13 nm with each of the six peripheral densities being 3-4 nm in diameter.     

Structure of the yeast vacuolar ATPase

The subunit architecture of the yeast vacuolar ATPase (V-ATPase) was analyzed by single particle transmission electron microscopy and electrospray ionization (ESI) tandem mass spectrometry. A three-dimensional model of the intact V-ATPase was calculated from two-dimensional projections of the complex at a resolution of 25 angstroms. Images of yeast V-ATPase decorated with monoclonal antibodies against subunits A, E, and G position subunit A within the pseudo-hexagonal arrangement in the V1, the N terminus of subunit G in the V1-V0 interface, and the C terminus of subunit E at the top of the V1 domain. ESI tandem mass spectrometry of yeast V1-ATPase showed that subunits E and G are most easily lost in collision-induced dissociation, consistent with a peripheral location of the subunits. An atomic model of the yeast V-ATPase was generated by fitting of the available x-ray crystal structures into the electron microscopy-derived electron density map. The resulting atomic model of the yeast vacuolar ATPase serves as a framework to help understand the role the peripheral stalk subunits are playing in the regulation of the ATP hydrolysis driven proton pumping activity of the vacuolar ATPase.     

Crystallographic structure refinement of Chromatium high potential iron protein at two Angstroms resolution.

The structure of Chromatium high potential iron protein (HiPIP) has been refined by semiautomatic Fo-Fc (observed minus calculated structure amplitude Fourier methods to a convential R index, R=sum of the absolute value of Fo-Fc divided by the sum of Fo, of 24.7% for a model in which bond distances and angles are constrained to standard values. Bond length and angle constraints were applied only intermittenly during the computations. At a late stage of the refinement, atomic parameters for only the Fe4S4 cluster plus the 4 associated cystein S-gamma atoms were adjusted by least squares methods and kept fixed during the rest of the refinement. The refined model consists of 625 of the 632 nonhydrogen atoms in the protein plus 75 water molecules. Seven side chain atoms could not be located in the final electron density map. A computer program rather than visual inspection was used wherever possible in the refinement: for locating water molecules, for removing water molecules that too closely approach other atoms, for deleting atoms that lay in regions of low electron density, and for evaluating the progress of refinement. Fo-Fc Fourier refinement is sufficiently economical to be applied routinely in protein crystal structure determinations. The complete HiPIP refinement required approximately 12 hours of CDC 3600 computer time and cost less than $3000 starting from a "trial structure," based upon multipe isomorphoous replacement phases, which gave an R of 43%...     

Structure and function of chromophores in R-Phycoerythrin at 1.9 A resolution

The crystal structure of R-Phycoerythrin (R-PE) from Polysiphonia urceolata has been refined to a resolution of 1.9 A, based on the atomic coordinates of R-PE determined at 2.8 A resolution, through the use of difference Fourier method and steorochemistry parameters restrained refinement with model adjustment according to the electron density map. Crystallographic R-factor of the refined model is 0.195 (Rfree = 0.282) from 8-1.9 A. High resolution structure of R-PE showed precise interactions between the chromophores and protein residues, which explained the spectrum characteristic and function of chromophores. Four chiral atoms of phycourobilin (PUB) were identified as C(4)-S, C(16)-S, C(21)-S, and C(20)-R. In addition to the coupling distances of 19 A to 45 A between the chromophores which were observed and involved in the energy transfer pathway, high resolution structure of R-PE suggested other pathways of energy transfer, such as the ultrashort distance between alpha140a and beta155. It has been proposed that aromatic residues in linker proteins not only influence the conformation of chromophore, but may also bridge chromophores to improve the energy transfer efficiency.     

Sequencing a protein by x-ray crystallography. II. Refinement of yeast hexokinase B co-ordinates and sequence at 2.1 A resolution.

Although the amino acid sequence of yeast hexokinase B has not been determined by chemical means, crystallographic refinement of the hexokinase monomer was carried out at 2.1 Å resolution to improve both the atomic co-ordinates and the amino acid sequence, which had been obtained from a 2.5 Å electron density map. The atomic co-ordinates were adjusted by real-space refinement into a multiple isomorphous replacement map, followed by automated difference Fourier refinement, and restrained parameter structure factor least-squares refinement. The amino acid sequence was altered periodically after visual inspection of (F_o ? F_c) difference electron density maps. Evidence of the improvement in the amino acid sequence was provided by the better agreement between the X-ray and chemically derived amino acid compositions, and most importantly by the ability to locate two short peptides which had been chemically sequenced. While only 6 out of the 18 residues in these two peptides agree with the sequence of the original model, 12 residues agree with the sequence of the refined model and the others differ by only an atom or two. The refined model contains 3293 of of the 3596 non-hydrogen atoms expected from the amino acid composition and 152 bound water molecules. The crystallographic R factor at 2.1 Å is 0.25.We show that there are several advantages to refining the structure of even a protein of unknown sequence. (1) Improved phases can be obtained to the resolution limit of the diffraction pattern starting with a model derived from a 2.5 Å map. (2) The accuracy of the amino acid sequence derived by X-ray methods alone can be substantially improved. (3) Functionally important residues can be identified before chemical sequence information is available. (4) The improved X-ray sequence should greatly reduce the effort required to obtain a chemical sequence; since peptides as short as eight or nine residues can be located in the refined X-ray sequence, peptides do not need to be overlapped by chemical means.     

Three-dimensional structure of basal body triplet revealed by electron cryo-tomography

Basal bodies and centrioles play central roles in microtubule (MT)‐organizing centres within many eukaryotes. They share a barrel‐shaped cylindrical structure composed of nine MT triplet blades. Here, we report the structure of the basal body triplet at 33 Å resolution obtained by electron cryo‐tomography and 3D subtomogram averaging. By fitting the atomic structure of tubulin into the EM density, we built a pseudo‐atomic model of the tubulin protofilaments at the core of the triplet. The 3D density map reveals additional densities that represent non‐tubulin proteins attached to the triplet, including a large inner circular structure in the basal body lumen, which functions as a scaffold to stabilize the entire basal body barrel. We found clear longitudinal structural variations along the basal body, suggesting a sequential and coordinated assembly mechanism. We propose a model in which δ‐tubulin and other components participate in the assembly of the basal body.     

Conjoined Use of EM and NMR in RNA Structure Refinement

More than 40% of the RNA structures have been determined using nuclear magnetic resonance (NMR) technique. NMR mainly provides local structural information of protons and works most effectively on relatively small biomacromolecules. Hence structural characterization of large RNAs can be difficult for NMR alone. Electron microscopy (EM) provides global shape information of macromolecules at nanometer resolution, which should be complementary to NMR for RNA structure determination. Here we developed a new energy term in Xplor-NIH against the density map obtained by EM. We conjointly used NMR and map restraints for the structure refinement of three RNA systems — U2/U6 small-nuclear RNA, genome-packing motif (Ψ^CD)₂ from Moloney murine leukemia virus, and ribosome-binding element from turnip crinkle virus. In all three systems, we showed that the incorporation of a map restraint, either experimental or generated from known PDB structure, greatly improves structural precision and accuracy. Importantly, our method does not rely on an initial model assembled from RNA duplexes, and allows full torsional freedom for each nucleotide in the torsion angle simulated annealing refinement. As increasing number of macromolecules can be characterized by both NMR and EM, the marriage between the two techniques would enable better characterization of RNA three-dimensional structures.     

The atomic structure of crystalline porcine pancreatic elastase at 2.5 A resolution: comparisons with the structure of alpha-chymotrypsin

Three isomorphous heavy-atom derivatives have been used to calculate a 2.5 Å resolution electron density map of tosyl-elastase at pH 5.0, from which an accurate atomic model has been constructed. Atomic co-ordinates measured from this model have been refined using model building, real-space refinement and energy minimization programs. The three-dimensional conformation of the polypeptide chain is described in terms of conformational angles, hydrogen-bonding networks and the environment of different types of amino acid side-chain.Difference Fourier calculation of the high resolution structure of native elastase at pH 5.0 shows it to be virtually identical to that of the tosyl derivative, except near the tosyl group. The conformation of the catalytically important residues in native elastase is very similar to that of native α-chymotrypsin, except for the orientation of the active centre serine oxygen. The significance of important structural similarities and differences between these two enzymes is discussed.Elastase contains 25 internal water molecules which play an important role in stabilizing the active conformation of the enzyme. Many of these water molecules are in identical positions to those found in the interior of α-chymotrypsin     

Structure of integrin alpha5beta1 in complex with fibronectin

The membrane-distal headpiece of integrins has evolved to specifically bind large extracellular protein ligands, but the molecular architecture of the resulting complexes has not been determined. We used molecular electron microscopy to determine the three-dimensional structure of the ligand-binding headpiece of integrin alpha5beta1 complexed with fragments of its physiological ligand fibronectin. The density map for the unliganded alpha5beta1 headpiece shows a 'closed' conformation similar to that seen in the alphaVbeta3 crystal structure. By contrast, binding to fibronectin induces an 'open' conformation with a dramatic, approximately 80 degrees change in the angle of the hybrid domain of the beta subunit relative to its I-like domain. The fibronectin fragment binds to the interface between the beta-propeller and I-like domains in the integrin headpiece through the RGD-containing module 10, but direct contact of the synergy-region-containing module 9 to integrin is not evident. This finding is corroborated by kinetic analysis of real-time binding data, which shows that the synergy site greatly enhances k(on) but has little effect on the stability or k(off) of the complex.     

Polarizable Atomic Multipole X-Ray Refinement: Hydration Geometry and Application to Macromolecules

We recently developed a polarizable atomic multipole refinement method assisted by the AMOEBA force field for macromolecular crystallography. Compared to standard refinement procedures, the method uses a more rigorous treatment of x-ray scattering and electrostatics that can significantly improve the resultant information contained in an atomic model. We applied this method to high-resolution lysozyme and trypsin data sets, and validated its utility for precisely describing biomolecular electron density, as indicated by a 0.4-0.6% decrease in the R- and R_free-values, and a corresponding decrease in the relative energy of 0.4-0.8 Kcal/mol/residue. The re-refinements illustrate the ability of force-field electrostatics to orient water networks and catalytically relevant hydrogens, which can be used to make predictions regarding active site function, activity, and protein-ligand interaction energies. Re-refinement of a DNA crystal structure generates the zigzag spine pattern of hydrogen bonding in the minor groove without manual intervention. The polarizable atomic multipole electrostatics model implemented in the AMOEBA force field is applicable and informative for crystal structures solved at any resolution.     

Multiple subunit fitting into a low-resolution density map of a macromolecular complex using a gaussian mixture model

Recently, electron microscopy measurement of single particles has enabled us to reconstruct a low-resolution 3D density map of large biomolecular complexes. If structures of the complex subunits can be solved by x-ray crystallography at atomic resolution, fitting these models into the 3D density map can generate an atomic resolution model of the entire large complex. The fitting of multiple subunits, however, generally requires large computational costs; therefore, development of an efficient algorithm is required. We developed a fast fitting program, “gmfit”, which employs a Gaussian mixture model (GMM) to represent approximated shapes of the 3D density map and the atomic models. A GMM is a distribution function composed by adding together several 3D Gaussian density functions. Because our model analytically provides an integral of a product of two distribution functions, it enables us to quickly calculate the fitness of the density map and the atomic models. Using the integral, two types of potential energy function are introduced: the attraction potential energy between a 3D density map and each subunit, and the repulsion potential energy between subunits. The restraint energy for symmetry is also employed to build symmetrical origomeric complexes. To find the optimal configuration of subunits, we randomly generated initial configurations of subunit models, and performed a steepest-descent method using forces and torques of the three potential energies. Comparison between an original density map and its GMM showed that the required number of Gaussian distribution functions for a given accuracy depended on both resolution and molecular size. We then performed test fitting calculations for simulated low-resolution density maps of atomic models of homodimer, trimer, and hexamer, using different search parameters. The results indicated that our method was able to rebuild atomic models of a complex even for maps of 30 Å resolution if sufficient numbers (eight or more) of Gaussian distribution functions were employed for each subunit, and the symmetric restraints were assigned for complexes with more than three subunits. As a more realistic test, we tried to build an atomic model of the GroEL/ES complex by fitting 21-subunit atomic models into the 3D density map obtained by cryoelectron microscopy using the C7 symmetric restraints. A model with low root mean-square deviations (14.7 Å) was obtained as the lowest-energy model, showing that our fitting method was reasonably accurate. Inclusion of other restraints from biological and biochemical experiments could further enhance the accuracy.     

The three-dimensional structure of bovine rhodopsin determined by electron cryomicroscopy

G-protein-coupled receptors are integral membrane proteins that respond to environmental signals and initiate signal transduction pathways, which activate cellular processes. Rhodopsin, a well known member of the G-protein-coupled receptor family, is located in the disk membranes of the rod outer segment, where it is responsible for the visualization of dim light. Rhodopsin is the most extensively studied G-protein-coupled receptor, and knowledge about its structure serves as a template for other related receptors. We have gained detailed structural knowledge from the crystal structure (1), which was solved by x-ray crystallography in 2000 using three-dimensional crystals. Here we report a three-dimensional density map of bovine rhodopsin determined by electron cryomicroscopy of two-dimensional crystals with p22(1)2(1) symmetry. The usage of relatively small and disordered crystals made the process of structure determination challenging. Special attention was paid to the extraction of amplitudes and phases, since usable raw data were limited to a maximum tilt of 45 degrees. In the refinement process, an improved unbending procedure was applied. This led to a final resolution of 5.5 A in the membrane plane and approximately 13 A perpendicular to it, making our electron density map the most accurate map of a G-protein-coupled receptor currently available by electron microscopy. Most important is the information we gain about the center of the membrane plane and the orientation of the molecule relative to the bilayer. This information cannot be retrieved from the three-dimensional crystals. In our electron density map, all seven transmembrane helices were identified, and their arrangement is in agreement with the arrangement known from the crystal structure (1). In the retinal binding pocket, a density peak adjacent to helix 3 suggests the position of the beta-ionine ring of the chromophore, and in its vicinity several of the bigger amino acids can be identified.     

Refining the structure of the strongly bound actin-myosin complex by protein docking

The structure of the strongly bound complex of the globular myosin head and F-actin is a key for understanding some important details of the mechanism of the actin-myosin motor. Current knowledge about the structure is based on the docking of known atomic structures of actin and myosin heads into low-resolution EM electron density maps. To refine the structure, we suggested a new approach based on energy minimization using the ICM-Pro software. The minimization includes rigid-body movement of protein backbone and side chain optimization on the protein interface. Our best model structure is similar to that obtained from EM. It also provides the highest calculated interaction energy and agrees with a number of mutagenesis experiments. Using the structure, we suggest molecular explanations for actin activation of product release from myosin and actin-induced myosin dissociation.     

Structure of V-type ATPase from Clostridium fervidus by electron microscopy

F-type and V-type ATPases couple synthesis or hydrolysis of ATP to the translocation of H+ or Na+ across biological membranes and have similarities in structure and mechanism. In both types of enzymes three main parts can be distinguished: headpiece, membrane-bound piece and stalk region. We report on structural details of the membrane sector and stalk region, including the stator, of V-type ATPase from Clostridium fervidus, as determined by electron microscopy. Besides visualization of the stator structure, one of the main findings is that in certain projections the central stalk connecting V1 and V₀ makes an angle of about 70° with the membrane. Implications for the subunit arrangement in V-type and F-type ATPase are discussed.     

Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics

A novel method to flexibly fit atomic structures into electron microscopy (EM) maps using molecular dynamics simulations is presented. The simulations incorporate the EM data as an external potential added to the molecular dynamics force field, allowing all internal features present in the EM map to be used in the fitting process, while the model remains fully flexible and stereochemically correct. The molecular dynamics flexible fitting (MDFF) method is validated for available crystal structures of protein and RNA in different conformations; measures to assess and monitor the fitting process are introduced. The MDFF method is then used to obtain high-resolution structures of the E. coli ribosome in different functional states imaged by cryo-EM.     

X-ray structure of lipoamide dehydrogenase from Azotobacter vinelandii determined by a combination of molecular and isomorphous replacement techniques

The crystal structure of lipoamide dehydrogenase from Azotobacter vinelandii has been determined by a combination of molecular replacement and isomorphous replacement techniques yielding eventually a good-quality 2.8 A electron density map. Initially, the structure determination was attempted by molecular replacement procedures alone using a model of human glutathione reductase, which has 26% sequence identity with this bacterial dehydrogenase. The rotation function yielded the correct orientation of the model structure both when the glutathione reductase dimer and monomer were used as starting model. The translation function could not be solved, however. Consequently, data for two heavy-atom derivatives were collected using the Hamburg synchotron facilities. The derivatives had several sites in common, which was presumably a major reason why the electron density map obtained by isomorphous information alone was of poor quality. Application of solvent flattening procedures cleaned up the map considerably, however, showing clearly the outline of the lipoamide dehydrogenase dimer, which has a molecular weight of 100,000. Application of the "phased translation function", which combines the phase information of both isomorphous and molecular replacement, led to an unambiguous determination of the position of the model structure in the lipoamide dehydrogenase unit cell. The non-crystallographic 2-fold axis of the dimer was optimized by several cycles of constrained-restrained least-squares refinement and subsequently used for phase improvement by 2-fold density averaging. After ten cycles at 3.5 A, the resolution was gradually extended to 2.8 A in another 140 cycles. The 2.8 A electron density distribution obtained in this manner was of much improved quality and allowed building of an atomic model of A. vinelandii lipoamide dehydrogenase. It appears that in the orthorhombic crystals used each dimer is involved in contacts with eight surrounding dimers, leaving unexplained why the crystals are rather fragile. Contacts between subunits within one dimer, which are quite extensive, can be divided into two regions separated by a cavity. In one of the contact regions, the level of sequence identity with glutathione reductase is very low but it is quite high in the other. The folding of the polypeptide chain in each subunit is quite similar to that of glutathione reductase, as is the extended conformation of the co-enzyme FAD.(ABSTRACT TRUNCATED AT 400 WORDS)