Combined approaches to flexible fitting and assessment in virus capsids undergoing conformational change

Fitting of atomic components into electron cryo-microscopy (cryoEM) density maps is routinely used to understand the structure and function of macromolecular machines. Many fitting methods have been developed, but a standard protocol for successful fitting and assessment of fitted models has yet to be agreed upon among the experts in the field. Here, we created and tested a protocol that highlights important issues related to homology modelling, density map segmentation, rigid and flexible fitting, as well as the assessment of fits. As part of it, we use two different flexible fitting methods (Flex-EM and iMODfit) and demonstrate how combining the analysis of multiple fits and model assessment could result in an improved model. The protocol is applied to the case of the mature and empty capsids of Coxsackievirus A7 (CAV7) by flexibly fitting homology models into the corresponding cryoEM density maps at 8.2 and 6.1 Å resolution. As a result, and due to the improved homology models (derived from recently solved crystal structures of a close homolog – EV71 capsid – in mature and empty forms), the final models present an improvement over previously published models. In close agreement with the capsid expansion observed in the EV71 structures, the new CAV7 models reveal that the expansion is accompanied by ∼5° counterclockwise rotation of the asymmetric unit, predominantly contributed by the capsid protein VP1. The protocol could be applied not only to viral capsids but also to many other complexes characterised by a combination of atomic structure modelling and cryoEM density fitting.     

Detection of secondary and supersecondary structures of proteins from cryo-electron microscopy

Recent advances in three-dimensional electron microscopy (3D EM) have enabled the quantitative visualization of the structural building blocks of proteins at improved resolutions. We provide algorithms to detect the secondary structures (α-helices and β-sheets) from proteins for which the volumetric maps are reconstructed at 6–10 Å resolution. Additionally, we show that when the resolution is coarser than 10 Å, some of the supersecondary structures can be detected from 3D EM maps. For both these algorithms, we employ tools from computational geometry and differential topology, specifically the computation of stable/unstable manifolds of certain critical points of the distance function induced by the molecular surface. Our results connect mathematically well-defined constructions with bio-chemically induced structures observed in proteins.     

Evolutionary bidirectional expansion for the tracing of alpha helices in cryo-electron microscopy reconstructions

Cryo-electron microscopy (cryo-EM) enables the imaging of macromolecular complexes in near-native environments at resolutions that often permit the visualization of secondary structure elements. For example, alpha helices frequently show consistent patterns in volumetric maps, exhibiting rod-like structures of high density. Here, we introduce VolTrac (Volume Tracer) – a novel technique for the annotation of alpha-helical density in cryo-EM data sets. VolTrac combines a genetic algorithm and a bidirectional expansion with a tabu search strategy to trace helical regions. Our method takes advantage of the stochastic search by using a genetic algorithm to identify optimal placements for a short cylindrical template, avoiding exploration of already characterized tabu regions. These placements are then utilized as starting positions for the adaptive bidirectional expansion that characterizes the curvature and length of the helical region. The method reliably predicted helices with seven or more residues in experimental and simulated maps at intermediate (4–10 Å) resolution. The observed success rates, ranging from 70.6% to 100%, depended on the map resolution and validation parameters. For successful predictions, the helical axes were located within 2 Å from known helical axes of atomic structures.     

BCL::EM-Fit: rigid body fitting of atomic structures into density maps using geometric hashing and real space refinement

Cryo-electron microscopy (cryoEM) can visualize large macromolecular assemblies at resolutions often below 10? and recently as good as 3.8-4.5 ?. These density maps provide important insights into the biological functioning of molecular machineries such as viruses or the ribosome, in particular if atomic-resolution crystal structures or models of individual components of the assembly can be placed into the density map. The present work introduces a novel algorithm termed BCL::EM-Fit that accurately fits atomic-detail structural models into medium resolution density maps. In an initial step, a "geometric hashing" algorithm provides a short list of likely placements. In a follow up Monte Carlo/Metropolis refinement step, the initial placements are optimized by their cross correlation coefficient. The resolution of density maps for a reliable fit was determined to be 10 ? or better using tests with simulated density maps. The algorithm was applied to fitting of capsid proteins into an experimental cryoEM density map of human adenovirus at a resolution of 6.8 and 9.0 ?, and fitting of the GroEL protein at 5.4 ?. In the process, the handedness of the cryoEM density map was unambiguously identified. The BCL::EM-Fit algorithm offers an alternative to the established Fourier/Real space fitting programs. BCL::EM-Fit is free for academic use and available from a web server or as downloadable binary file at http://www.meilerlab.org.     

Biased coarse-grained molecular dynamics simulation approach for flexible fitting of X-ray structure into cryo electron microscopy maps

Several approaches have been introduced to interpret, in terms of high-resolution structure, low-resolution structural data as obtained from cryo-EM. As conformational changes are often observed in biological molecules, these techniques need to take into account the flexibility of proteins. Flexibility has been described in terms of movement between rigid domains and between rigid secondary structure elements, which present some limitations for studying dynamical properties. Normal mode analysis has also been used, but is limited to medium resolution data. All-atom molecular dynamics fitting techniques are more appropriate to fit structures into higher-resolution data as full protein flexibility is considered, but are cumbersome in terms of computational time. Here, we introduce a coarse-grained approach; a Go-model was used to represent biological molecules, combined with biased molecular dynamics to reproduce accurately conformational transitions. Illustrative examples on simulated data are shown. Accurate fittings can be obtained for resolution ranging from 5 to 20 Å. The approach was also tested on experimental data of Elongation Factor G and Escherichia coli RNA polymerase, where its validity is compared to previous models obtained from different techniques. This comparison demonstrates that quantitative flexible techniques, as opposed to manual docking, need to be considered to interpret low-resolution data.     

Scoring functions for cryoEM density fitting

In fitting atomic structures into cryoEM density maps of macromolecular assemblies, the cross-correlation function (CCF) is the most prevalent method of scoring the goodness-of-fit. However, there are still many possible, less studied ways of scoring fits. In this paper, we introduce four scores new to cryoEM fitting and compare their performance to three known scores. Our benchmark consists of (a) 4 protein assemblies with simulated maps at 5-20 ? resolution, including the heptameric ring of GroEL; and (b) 4 experimental maps of GroEL at ～6-23 ? resolution with corresponding fitted atomic models. We perturb each fit 1000 times and assess each new fit with each score. The correlation between a score and the Cα RMSD of each fit from the "correctly" fitted structure shows that the CCF is one of the best scores, but in certain situations could be augmented or even replaced by other scores. For instance, our implementation of a score based on mutual information outperforms or is comparable to the CCF in almost all test cases, and our new "envelope score" works as well as the CCF at sub-nanometer resolution but is an order of magnitude faster to calculate. The results also suggest that the width of the Gaussian function used to blur the atomic structure into a density map can significantly affect the fitting process. Finally, we show that our score-testing method, when combined with the Laplacian CCF or the mutual information scores, can be used as a statistical tool for improving cryoEM density fitting.     

Macromolecular docking restrained by a small angle X-ray scattering profile

While many structures of single protein components are becoming available, structural characterization of their complexes remains challenging. Methods for modeling assembly structures from individual components frequently suffer from large errors, due to protein flexibility and inaccurate scoring functions. However, when additional information is available, it may be possible to reduce the errors and compute near-native complex structures. One such type of information is a small angle X-ray scattering (SAXS) profile that can be collected in a high-throughput fashion from a small amount of sample in solution. Here, we present an efficient method for protein–protein docking with a SAXS profile (FoXSDock): generation of complex models by rigid global docking with PatchDock, filtering of the models based on the SAXS profile, clustering of the models, and refining the interface by flexible docking with FireDock. FoXSDock is benchmarked on 124 protein complexes with simulated SAXS profiles, as well as on 6 complexes with experimentally determined SAXS profiles. When induced fit is less than 1.5 Å interface C_α RMSD and the fraction residues of missing from the component structures is less than 3%, FoXSDock can find a model close to the native structure within the top 10 predictions in 77% of the cases; in comparison, docking alone succeeds in only 34% of the cases. Thus, the integrative approach significantly improves on molecular docking alone. The improvement arises from an increased resolution of rigid docking sampling and more accurate scoring.     

Rapid and simple preparation of mushroom DNA directly from colonies and fruiting bodies for PCR

We have optimized a simple and rapid preparation procedure for mushroom DNA extraction from colonies on media or from fruiting bodies for PCR amplification. The protocol combines microwaving twice for 1 min, cooling for 10 min, and centrifuging for 5 min. By using this procedure, more than 100 samples of mushroom DNA can be prepared within 1 h. The DNA obtained can be used for (1) identifying mushroom species by PCR and subsequent sequencing, (2) amplifying low copy number genes (at least 2,000 bp), and (3) screening genetic transformants. This technique will contribute to the mycology of mushroom species.     

Refinement of protein structures by iterative comparative modeling and CryoEM density fitting

We developed a method for structure characterization of assembly components by iterative comparative protein structure modeling and fitting into cryo-electron microscopy (cryoEM) density maps. Specifically, we calculate a comparative model of a given component by considering many alternative alignments between the target sequence and a related template structure while optimizing the fit of a model into the corresponding density map. The method relies on the previously developed Moulder protocol that iterates over alignment, model building, and model assessment. The protocol was benchmarked using 20 varied target-template pairs of known structures with less than 30% sequence identity and corresponding simulated density maps at resolutions from 5A to 25A. Relative to the models based on the best existing sequence profile alignment methods, the percentage of C(alpha) atoms that are within 5A of the corresponding C(alpha) atoms in the superposed native structure increases on average from 52% to 66%, which is half-way between the starting models and the models from the best possible alignments (82%). The test also reveals that despite the improvements in the accuracy of the fitness function, this function is still the bottleneck in reducing the remaining errors. To demonstrate the usefulness of the protocol, we applied it to the upper domain of the P8 capsid protein of rice dwarf virus that has been studied by cryoEM at 6.8A. The C(alpha) root-mean-square deviation of the model based on the remotely related template, bluetongue virus VP7, improved from 8.7A to 6.0A, while the best possible model has a C(alpha) RMSD value of 5.3A. Moreover, the resulting model fits better into the cryoEM density map than the initial template structure. The method is being implemented in our program MODELLER for protein structure modeling by satisfaction of spatial restraints and will be applicable to the rapidly increasing number of cryoEM density maps of macromolecular assemblies.     

Purification and crystallization of the ABC-type transport substrate-binding protein OppA from Thermoanaerobacter tengcongensis

Di- and oligopeptide- binding protein OppAs play important roles in solute and nutrient uptake, sporulation, biofilm formation, cell wall muropeptides recycling, peptide-dependent quorum-sensing responses, adherence to host cells, and a variety of other biological processes. Soluble OppA from Thermoanaerobacter tengcongensis was expressed in Escherichia coli. The protein was found to be >95% pure with SDS–PAGE after a series of purification steps and the purity was further verified by mass spectrometry. The protein was crystallized using the sitting-drop vapour-diffusion method with PEG 400 as the precipitant. Crystal diffraction extended to 2.25 Å. The crystal belonged to space group C222₁, with unit-cell parameters of a = 69.395, b = 199.572, c = 131.673 Å, and α = β = γ = 90°.     

Thermolabile antifreeze protein produced in Escherichia coli for structural analysis

Inclusion bodies of recombinant human growth hormone (r-hGH) were isolated from Escherichia coli, enriched and solubilized in 100 mM Tris buffer containing 6 M n-propanol and 2 M urea. Around 4 mg/ml of r-hGH from inclusion bodies were solubilized in 6 M n-propanol-based buffer containing 2 M urea. Existence of native-like secondary structure of r-hGH in 6 M n-propanol solution was confirmed by CD and fluorescence spectra. Solubilized r-hGH was subsequently refolded by pulsatile dilution, purified to homogeneity and found to be functionally active. Tris buffer containing 6 M n-propanol and 2 M urea also effectively solubilized a number of proteins expressed as inclusion bodies in E. coli. Mild solubilization of inclusion body proteins, chaotropic effect of n-propanol at high concentration and kosmotropic effect at lower concentration helped in improved refolding of the solubilized protein. Around 40% of the r-hGH in the form of inclusion body aggregates was refolded into bioactive form while using n-propanol as solubilization agent. Solubilization with 6 M n-propanol solution thus can be a viable alternative for achieving high throughput recovery of bioactive protein from inclusion bodies of E. coli.     

Flexible fitting in 3D-EM guided by the structural variability of protein superfamilies

A method for flexible fitting of molecular models into three-dimensional electron microscopy (3D-EM) reconstructions at a resolution range of 8-12 A is proposed. The approach uses the evolutionarily related structural variability existing among the protein domains of a given superfamily, according to structural databases such as CATH. A structural alignment of domains belonging to the superfamily, followed by a principal components analysis, is performed, and the first three principal components of the decomposition are explored. Using rigid body transformations for the secondary structure elements (SSEs) plus the cyclic coordinate descent algorithm to close the loops, stereochemically correct models are built for the structure to fit. All of the models are fitted into the 3D-EM map, and the best one is selected based on crosscorrelation measures. This work applies the method to both simulated and experimental data and shows that the flexible fitting was able to produce better results than rigid body fitting.     

Structural characterization of components of protein assemblies by comparative modeling and electron cryo-microscopy

We explore structural characterization of protein assemblies by a combination of electron cryo-microscopy (cryoEM) and comparative protein structure modeling. Specifically, our method finds an optimal atomic model of a given assembly subunit and its position within an assembly by fitting alternative comparative models into a cryoEM map. The alternative models are calculated by MODELLER [J. Mol. Biol. 234 (1993) 313] from different sequence alignments between the modeled protein and its template structures. The fitting of these models into a cryoEM density map is performed either by FOLDHUNTER [J. Mol. Biol. 308 (2001) 1033] or by a new density fitting module of MODELLER (Mod-EM). Identification of the most accurate model is based on the correlation between the model accuracy and the quality of fit into the cryoEM density map. To quantify this correlation, we created a benchmark consisting of eight proteins of different structural folds with corresponding density maps simulated at five resolutions from 5 to 15 angstroms, with three noise levels each. Each of the proteins in the set was modeled based on 300 different alignments to their remotely related templates (12-32% sequence identity), spanning the range from entirely inaccurate to essentially accurate alignments. The benchmark revealed that one of the most accurate models can usually be identified by the quality of its fit into the cryoEM density map, even for noisy maps at 15 angstroms resolution. Therefore, a cryoEM density map can be helpful in improving the accuracy of a comparative model. Moreover, a pseudo-atomic model of a component in an assembly may be built better with comparative models of the native subunit sequences than with experimentally determined structures of their homologs.     

Protein structure fitting and refinement guided by cryo-EM density

For many macromolecular assemblies, both a cryo-electron microscopy map and atomic structures of its component proteins are available. Here we describe a method for fitting and refining a component structure within its map at intermediate resolution (<15 A). The atomic positions are optimized with respect to a scoring function that includes the crosscorrelation coefficient between the structure and the map as well as stereochemical and nonbonded interaction terms. A heuristic optimization that relies on a Monte Carlo search, a conjugate-gradients minimization, and simulated annealing molecular dynamics is applied to a series of subdivisions of the structure into progressively smaller rigid bodies. The method was tested on 15 proteins of known structure with 13 simulated maps and 3 experimentally determined maps. At approximately 10 A resolution, Calpha rmsd between the initial and final structures was reduced on average by approximately 53%. The method is automated and can refine both experimental and predicted atomic structures.     

Redox properties of tyrosine and related molecules

Redox reactions of tyrosine play key roles in many biological processes, including water oxidation and DNA synthesis. We first review the redox properties of tyrosine (and other phenols) in small molecules and related polypeptides, then report work on (H20)/(Y48)-modified Pseudomonas aeruginosa azurin. The crystal structure of this protein (1.18 Å resolution) shows that H20 is strongly hydrogen bonded to Y48 (2.7–2.8 Å tyrosine-O to histidine-N distance). A firm conclusion is that proper tuning of the tyrosine potential by a proton-accepting base is critical for biological redox functions.     

Filaments from Ignicoccus hospitalis Show Diversity of Packing in Proteins Containing N-Terminal Type IV Pilin Helices

Bacterial motility is driven by the rotation of flagellar filaments that supercoil. The supercoiling involves the switching of coiled-coil protofilaments between two different states. In archaea, the flagellar filaments responsible for motility are formed by proteins with distinct homology in their N-terminal portion to bacterial Type IV pilins. The bacterial pilins have a single N-terminal hydrophobic α-helix, not the coiled coil found in flagellin. We have used electron cryo-microscopy to study the adhesion filaments from the archaeon Ignicoccus hospitalis. While I. hospitalis is non-motile, these filaments make transitions between rigid stretches and curved regions and appear morphologically similar to true archaeal flagellar filaments. A resolution of ~ 7.5 Å allows us to unambiguously build a model for the packing of these N-terminal α-helices, and this packing is different from several bacterial Type IV pili whose structure has been analyzed by electron microscopy and modeling. Our results show that the mechanism responsible for the supercoiling of bacterial flagellar filaments cannot apply to archaeal filaments.     

Crystal structure of human MTH1 and the 8-oxo-dGMP product complex

MTH1 hydrolyzes oxidized nucleotide triphosphates, thereby preventing them from being incorporated into DNA. We here present the structures of human MTH1 (1.9 Å) and its complex with the product 8-oxo-dGMP (1.8 Å). Unexpectedly MTH1 binds the nucleotide in the anti conformation with no direct interaction between the 8-oxo group and the protein. We suggest that the specificity depends on the stabilization of an enol tautomer of the 8-oxo form of dGTP. The binding of the product induces no major structural changes. The structures reveal the mode of nucleotide binding in MTH1 and provide the structural basis for inhibitor design.     

Three-dimensional structure of the detergent-solubilized Vibrio cholerae cytolysin (VCC) heptamer by electron cryomicroscopy

Vibrio cholerae cytolysin (VCC) is a pore-forming toxin that inserts a lytic water-filled channel into susceptible host membranes. Assembly of the toxin on cell surfaces may be enhanced by two tandem lectin domains, in addition to direct interactions with lipids and cholesterol within the membrane itself. We used single-particle electron cryomicroscopy (cryoEM) to generate a low-resolution molecular structure of the detergent-solubilized VCC oligomer to 20 Å resolution. After confirming a heptameric arrangement of individual protomers, sevenfold averaging around the central pore was utilized to improve the structure. Docking of the previously determined VCC protoxin crystal structure was possible with rigid-body rearrangements between the cytolytic and lectin domains. A second cryoEM reconstruction of a truncated VCC mutant supported the topology of our model in which the carboxyl-terminal lectin domain forms “spikes” around the toxin core with the putative carbohydrate receptor-binding site accessible on the surface of the oligomer. This finding points to an assembly mechanism in which lectin domains may remain bound to receptors on the cell surface throughout assembly of the cytolytic toxin core and explains the hemagglutinating activity of purified toxin. Our model provides an insight into the structural rearrangements that accompany VCC-mediated cytolysis and may aid in the engineering of novel pore-forming toxins to attack specific cells towards therapeutic ends.     

The utility of geometrical and chemical restraint information extracted from predicted ligand-binding sites in protein structure refinement

Exhaustive exploration of molecular interactions at the level of complete proteomes requires efficient and reliable computational approaches to protein function inference. Ligand docking and ranking techniques show considerable promise in their ability to quantify the interactions between proteins and small molecules. Despite the advances in the development of docking approaches and scoring functions, the genome-wide application of many ligand docking/screening algorithms is limited by the quality of the binding sites in theoretical receptor models constructed by protein structure prediction. In this study, we describe a new template-based method for the local refinement of ligand-binding regions in protein models using remotely related templates identified by threading. We designed a Support Vector Regression (SVR) model that selects correct binding site geometries in a large ensemble of multiple receptor conformations. The SVR model employs several scoring functions that impose geometrical restraints on the Cα positions, account for the specific chemical environment within a binding site and optimize the interactions with putative ligands. The SVR score is well correlated with the RMSD from the native structure; in 47% (70%) of the cases, the Pearson’s correlation coefficient is >0.5 (>0.3). When applied to weakly homologous models, the average heavy atom, local RMSD from the native structure of the top-ranked (best of top five) binding site geometries is 3.1 Å (2.9 Å) for roughly half of the targets; this represents a 0.1 (0.3) Å average improvement over the original predicted structure. Focusing on the subset of strongly conserved residues, the average heavy atom RMSD is 2.6 Å (2.3 Å). Furthermore, we estimate the upper bound of template-based binding site refinement using only weakly related proteins to be ～2.6 Å RMSD. This value also corresponds to the plasticity of the ligand-binding regions in distant homologues. The Binding Site Refinement (BSR) approach is available to the scientific community as a web server that can be accessed at http://cssb.biology.gatech.edu/bsr/.     

Online immobilized enzyme microreactor for the glucose oxidase enzymolysis and enzyme inhibition assay

An online immobilized glucose oxidase (GOx) capillary microreactor was developed based on an enzymatic redox reaction with 1,4-benzoquinone as an acceptor of electrons, replacing the molecular oxygen typically used in a GOx reaction to achieve direct ultraviolet detection without derivation. A high efficiency of enzymolysis was obtained at 1 mg ml^?1 1,4-benzoquinone for 5 min of incubation at 25 °C, and baseline separation of the substrate and product could be achieved with a resolution of 3.85 by employing 20 mM phosphate buffer (pH 8.0) containing 40 mg ml^?1 sulfated β-cyclodextrin as an additive, a constant voltage of 15 kV, and a detection wavelength of 220 nm. In addition, an online enzyme inhibition study was performed on the immobilized GOx microreactor with metal ions Ag⁺ and Cu²⁺ used as model inhibitors. The results indicate that Ag⁺ (IC₅₀ = 69.16 μM) has a markedly higher inhibitory effect than Cu²⁺ (IC₅₀ = 1.33 mM). The protocol described can be applied in high-throughput screening of enzyme reactions and inhibitors.