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

Blind testing of cross-linking/mass spectrometry hybrid methods in CASP11

Hybrid approaches combine computational methods with experimental data. The information contained in the experimental data can be leveraged to probe the structure of proteins otherwise elusive to computational methods. Compared with computational methods, the structures produced by hybrid methods exhibit some degree of experimental validation. In spite of these advantages, most hybrid methods have not yet been validated in blind tests, hampering their development. Here, we describe the first blind test of a specific cross-link based hybrid method in CASP. This blind test was coordinated by the CASP organizers and utilized a novel, high-density cross-linking/mass-spectrometry (CLMS) approach that is able to collect high-density CLMS data in a matter of days. This experimental protocol was developed in the Rappsilber laboratory. This approach exploits the chemistry of a highly reactive, photoactivatable cross-linker to produce an order of magnitude more cross-links than homobifunctional cross-linkers. The Rappsilber laboratory generated experimental CLMS data based on this protocol, submitted the data to the CASP organizers which then released this data to the CASP11 prediction groups in a separate, CLMS assisted modeling experiment. We did not observe a clear improvement of assisted models, presumably because the properties of the CLMS data-uncertainty in cross-link identification and residue-residue assignment, and uneven distribution over the protein-were largely unknown to the prediction groups and their approaches were not yet tailored to this kind of data. We also suggest modifications to the CLMS-CASP experiment and discuss the importance of rigorous blind testing in the development of hybrid methods. Proteins 2016; 84(Suppl 1):152-163. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Structural, dynamic, and functional aspects of helix association in membranes: a computational view

This review surveys recent achievements of molecular computer modeling in understanding spatial structure, dynamics, and mechanisms of functioning of transmembrane α-helical dimers in membranes. The factors driving self-association of hydrophobic helices in the membrane milieu are considered with examples of their applications to biologically relevant problems. The emphasis is made on the recent results, which help to understand important aspects of structure-function relations for these systems and their biological activity. Limitations and shortcomings of the methods, along with their perspectives in design of new membrane active agents, are discussed.

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.

Lamin B receptor: multi-tasking at the nuclear envelope

Lamin B receptor (LBR) is an integral membrane protein of the interphase nuclear envelope (NE). The N-terminal end resides in the nucleoplasm, binding to lamin B and heterochromatin, with the interactions disrupted during mitosis. The C-terminal end resides within the inner nuclear membrane, retreating with the ER away from condensing chromosomes during mitotic NE breakdown. Some of these properties are interpretable in terms of our current structural knowledge of LBR, but many of the structural features remain unknown. LBR apparently has an evolutionary history which brought together at least two ancient conserved structural domains (i.e., Tudor and sterol reductase). This convergence may have occurred with the emergence of the chordates and echinoderms. It is not clear what survival values have maintained LBR structure during evolution. But it seems likely that roles in post-mitotic nuclear reformation, interphase NE growth and compartmentalization of nuclear architecture might have provided some evolutionary advantage to preservation of the LBR gene.

EM-fold: De novo folding of alpha-helical proteins guided by intermediate-resolution electron microscopy density maps

In medium-resolution (7-10 A) cryo-electron microscopy (cryo-EM) density maps, alpha helices can be identified as density rods whereas beta-strand or loop regions are not as easily discerned. We are proposing a computational protein structure prediction algorithm "EM-Fold" that resolves the density rod connectivity ambiguity by placing predicted alpha helices into the density rods and adding missing backbone coordinates in loop regions. In a benchmark of 11 mainly alpha-helical proteins of known structure a native-like model is identified in eight cases (rmsd 3.9-7.9 A). The three failures can be attributed to inaccuracies in the secondary structure prediction step that precedes EM-Fold. EM-Fold has been applied to the approximately 6 A resolution cryo-EM density map of protein IIIa from human adenovirus. We report the first topological model for the alpha-helical 400 residue N-terminal region of protein IIIa. EM-Fold also has the potential to interpret medium-resolution density maps in X-ray crystallography.

Hybrid approaches: applying computational methods in cryo-electron microscopy

Recent advances in cryo-electron microscopy have led to an increasing number of high (3-5A) to medium (5-10A) resolution cryoEM density maps. These density maps contain valuable information about the protein structure but frequently require computational algorithms to aid their structural interpretation. It is these hybrid approaches between cryoEM and computational protein structure prediction algorithms that will shape protein structure elucidation from density maps.

Structural refinement of membrane proteins by restrained molecular dynamics and solvent accessibility data

We present an approach for incorporating solvent accessibility data from electron paramagnetic resonance experiments in the structural refinement of membrane proteins through restrained molecular dynamics simulations. The restraints have been parameterized from oxygen (PiO(2)) and nickel-ethylenediaminediacetic acid (PiNiEdda) collision frequencies, as indicators of lipid or aqueous exposed spin-label sites. These are enforced through interactions between a pseudoatom representation of the covalently attached Nitroxide spin-label and virtual "solvent" particles corresponding to O(2) and NiEdda in the surrounding environment. Interactions were computed using an empirical potential function, where the parameters have been optimized to account for the different accessibilities of the spin-label pseudoatoms to the surrounding environment. This approach, "pseudoatom-driven solvent accessibility refinement", was validated by refolding distorted conformations of the Streptomyces lividans potassium channel (KcsA), corresponding to a range of 2-30 A root mean-square deviations away from the native structure. Molecular dynamics simulations based on up to 58 electron paramagnetic resonance restraints derived from spin-label mutants were able to converge toward the native structure within 1-3 A root mean-square deviations with minimal computational cost. The use of energy-based ranking and structure similarity clustering as selection criteria helped in the convergence and identification of correctly folded structures from a large number of simulations. This approach can be applied to a variety of integral membrane protein systems, regardless of oligomeric state, and should be particularly useful in calculating conformational changes from a known reference crystal structure.

Electron microscopy of integrins

Integrins are a family of heterodimeric, cell-surface receptors that mediate interactions between the cytoskeleton and the extracellular matrix. We have used electron microscopy and single-particle image analysis combined with molecular modeling to investigate the structures of the full-length integrin alpha(IIb)beta(3) and the ectodomain of alpha(V)beta(3) in a complex with fibronectin. The full-length integrin alpha(IIb)beta(3) is purified from human platelets by ion exchange and gel filtration chromatography in buffers containing the detergent octyl-beta-D-glucopyranoside, whereas the recombinant ectodomain of alpha(V)beta(3) is soluble in aqueous buffer. Transmission electron microscopy is performed either in negative stain, where the protein is embedded in a heavy metal such as uranyl acetate, or in the frozen-hydrated state, where the sample is flash-frozen such that the buffer is vitrified and native conditions are preserved. Individual integrin particles are selected from low-dose micrographs, either by manual identification or an automated method using a cross-correlation search of the micrograph against a set of reference images. Due to the small size of integrin heterodimers (approximately 250 kDa) and the low electron dose required to minimize beam damage, the signal-to-noise level of individual particles is quite low, both by negative-stain electron microscopy and electron cryomicroscopy. Consequently, it is necessary to average many particle images with equivalent views. The particle images are subjected to reference-free alignment and classification, in which the particles are aligned to a common view and further grouped by statistical methods into classes with common orientations. Assessment of the structure from a set of two-dimensional averaged projections is often difficult, and a further three-dimensional (3D) reconstruction analysis is performed to classify each particle as belonging to a specific projection from a single 3D model. The 3D reconstruction algorithm is an iterative projection-matching routine in which the classified particles are used to construct a new, 3D map for the next iteration. Docking of known high-resolution structures of individual subdomains within the molecular envelope of the 3D EM map is used to derive a pseudoatomic model of the integrin complex. This approach of 3D EM image analysis and pseudoatomic modeling is a powerful strategy for exploring the structural biology of transmembrane signaling by integrins because it is likely that multiple conformational states will be difficult to crystallize, whereas the different states should be amenable to electron cryomicroscopy.