MtsslWizard: In Silico Spin-Labeling and Generation of Distance Distributions in PyMOL

Connecting lipoxygenase function to structure by electron paramagnetic resonance

CONSPECTUS: Lipoxygenase enzymes insert oxygen in a polyunsaturated lipid, yielding a hydroperoxide product. When the acyl chain is arachidonate, with three cis-pentadiene units, 12 positionally and stereochemically different products might result. The plant lipids, linoleate and linolenate, have, respectively, four and eight potential oxygen insertion sites. The puzzle of how specificity is achieved in these reactions grows as more and more protein structures confirm the conservation of a lipoxygenase protein fold in plants, animals, and bacteria. Lipoxygenases are large enough (60-100 kDa) that they provide a protein shell completely surrounding an active site cavity that has the shape of a long acyl chain and contains a catalytic metal (usually iron). This Account summarizes electron paramagnetic resonance (EPR) spectroscopic, and other, experiments designed to bridge the gap between lipid-lipoxygenase interactions in solution and crystal structures. Experiments with spin-labeled lipids give a picture of bound lipids tethered to protein by an acyl chain, but with a polar end emerging from the cavity to solvent exposure, where the headgroup is highly flexible. The location of a spin on the polar end of a lysolecithin was determined by pulsed, dipolar EPR measurements, by representing the protein structure as a five-point grid of spin-labels with coordinates derived from 10 distance determinations between spin pairs. Distances from the lipid spin to each grid site completed a six-point representation of the enzyme with a bound lipid. Insight into the dynamics that allow substrate/product to enter/exit the cavity was obtained with a different set of spin-labeled protein mutants. Once substrate enters the cavity, the rate-limiting step of catalysis involves redox cycling at the metal center. Here, a mononuclear iron cycles between ferric and ferrous (high-spin) forms. Two helices provide pairs of side-chain ligands to the iron, resulting in characteristic EPR signals. Quantitative comparison of EPR spectra of plant and bacterial lipoxygenases has suggested conservation of a unique geometry of lipoxygenase iron centers. High frequency (94 GHz) EPR is consistent with a similar metal center in a manganese version of lipoxygenase. Overall, established and emerging EPR experiments have been developed and applied to the lipoxygenase family of enzymes to elucidate changes in the solution structures that are related to function.

BAX-induced apoptosis can be initiated through a conformational selection mechanism

BAX protein plays a key role in the mitochondria-mediated apoptosis. However, it remains unclear by what mechanism BAX is triggered to initiate apoptosis. Here, we reveal the mechanism using electron spin resonance (ESR) techniques. An inactive BAX monomer was found to exhibit conformational heterogeneity and exist at equilibrium in two conformations, one of which has never been reported. We show that upon apoptotic stimulus by BH3-only peptides, BAX can be induced to convert into either a ligand-bound monomer or an oligomer through a conformational selection mechanism. The kinetics of reaction is studied by means of time-resolved ESR, allowing a direct in situ observation for the transformation of BAX from the native to the bound states. In vitro mitochondrial assays provide further discrimination between the proposed BAX states, thereby revealing a population-shift allosteric mechanism in the process. BAX's apoptotic function is shown to critically depend on excursions between different structural conformations.

Closure of the cytoplasmic gate formed by TM5 and TM11 during transport in the oxalate/formate exchanger from Oxalobacter formigenes

OxlT, the oxalate/formate exchanger of Oxalobacter formigenes, is a member of the major facilitator superfamily of transporters. In the present work, substrate (oxalate) was found to enhance the reactivity of the cysteine mutant S336C on the cytoplasmic end of helix 11 to methanethiosulfonate ethyl carboxylate. In addition, S336C is found to spontaneously cross-link to S143C in TM5 in either native or reconstituted membranes under conditions that support transport. Continuous wave EPR measurements are consistent with this result and indicate that positions 143 and 336 are in close proximity in the presence of substrate. These two residues are localized within helix interacting GxxxG-like motifs (G₁₄₀LASG₁₄₄ and S₃₃₆DIFG₃₄₀) at the cytoplasmic poles of TM5 and TM11. Pulse EPR measurements were used to determine distances and distance distributions across the cytoplasmic or periplasmic ends of OxlT and were compared with the predictions of an inside-open homology model. The data indicate that a significant population of transporter is in an outside-open configuration in the presence of substrate; however, each end of the transporter exhibits significant conformational heterogeneity, where both inside-open and outside-open configurations are present. These data indicate that TM5 and TM11, which form part of the transport pathway, transiently close during transport and that there is a conformational equilibrium between inside-open and outside-open states of OxlT in the presence of substrate.

Modelling multi-protein complexes using PELDOR distance measurements for rigid body minimisation experiments using XPLOR-NIH

Crystallographic and NMR approaches have provided a wealth of structural information about protein domains. However, often these domains are found as components of larger multi domain polypeptides or complexes. Orienting domains within such contexts can provide powerful new insight into their function. The combination of site specific spin labelling and Pulsed Electron Double Resonance (PELDOR) provide a means of obtaining structural measurements that can be used to generate models describing how such domains are oriented. Here we describe a pipeline for modelling the location of thio-reactive nitroxyl spin locations to engineered sties on the histone chaperone Vps75. We then use a combination of experimentally determined measurements and symmetry constraints to model the orientation in which homodimers of Vps75 associate to form homotetramers using the XPLOR-NIH platform. This provides a working example of how PELDOR measurements can be used to generate a structural model.

Probing the structure of the mechanosensitive channel of small conductance in lipid bilayers with pulsed electron-electron double resonance

Mechanosensitive channel proteins are important safety valves against osmotic shock in bacteria, and are involved in sensing touch and sound waves in higher organisms. The mechanosensitive channel of small conductance (MscS) has been extensively studied. Pulsed electron-electron double resonance (PELDOR or DEER) of detergent-solubilized protein confirms that as seen in the crystal structure, the outer ring of transmembrane helices do not pack against the pore-forming helices, creating an apparent void. The relevance of this void to the functional form of MscS in the bilayer is the subject of debate. Here, we report PELDOR measurements of MscS reconstituted into two lipid bilayer systems: nanodiscs and bicelles. The distance measurements from multiple mutants derived from the PELDOR data are consistent with the detergent-solution arrangement of the protein. We conclude, therefore, that the relative positioning of the transmembrane helices is preserved in mimics of the cell bilayer, and that the apparent voids are not an artifact of detergent solution but a property of the protein that will have to be accounted for in any molecular mechanism of gating.

Identifying and quantitating conformational exchange in membrane proteins using site-directed spin labeling

Protein structures are not static but sample different conformations over a range of amplitudes and time scales. These fluctuations may involve relatively small changes in bond angles or quite large rearrangements in secondary structure and tertiary fold. The equilibrium between discrete structural substates on the microsecond to millisecond time scale is sometimes termed conformational exchange. Protein dynamics and conformational exchange are believed to provide the basis for many important activities, such as protein–protein and protein–ligand interactions, enzymatic activity and protein allostery; however, for many proteins, the dynamics and conformational exchange that lead to function are poorly defined.

Spectroscopic methods, such as NMR, are among the most important methods to explore protein dynamics and conformational exchange; however, they are difficult to implement in some systems and with some types of exchange events. Site-directed spin labeling (SDSL) is an EPR based approach that is particularly well-suited to high molecular-weight systems such as membrane proteins. Because of the relatively fast time scale for EPR spectroscopy, it is an excellent method to examine exchange. Conformations that are in exchange are captured as distinct populations in the EPR spectrum, and this feature when combined with the use of methods that can shift the free energy of conformational substates allows one to identify regions of proteins that are in dynamic exchange. In addition, modern pulse EPR methods have the ability to examine conformational heterogeneity, resolve discrete protein states, and identify the substates in exchange.

Protein crystallography has provided high-resolution models for a number of membrane proteins; but because of conformational exchange, these models do not always reflect the structures that are present when the protein is in a native bilayer environment. In the case of the Escherichia coli vitamin B₁₂ transporter, BtuB, the energy coupling segment of this protein undergoes a substrate-dependent unfolding, which acts to couple this outer-membrane protein to the inner-membrane protein TonB. EPR spectroscopy demonstrates that the energy coupling segment is in equilibrium between ordered and disordered states, and that substrate binding shifts this equilibrium to favor an unfolded state. However, in crystal structures of BtuB, this segment is resolved and folded within the protein, and neither the presence of this equilibrium nor the substrate-induced change is revealed. This is a result of the solute environment and the crystal lattice, both of which act to stabilize one conformational substate of the transporter.

Using SDSL, it can be shown that conformational exchange is present in other regions of BtuB and in other members of this transporter family. Conformational exchange has also been examined in systems such as the plasma membrane SNARE protein, syntaxin 1A, where dynamics are controlled by regulatory proteins such as munc18. Regulating the microsecond to millisecond time scale dynamics in the neuronal SNAREs is likely to be a key feature that regulates assembly of the SNAREs and neurotransmitter release.

Analysis of influenza A virus NS1 dimer interfaces in solution by pulse EPR distance measurements

Pulsed electron–electron double resonance (PELDOR) is an electron paramagnetic resonance (EPR) spectroscopy technique for nanometer distance measurements between paramagnetic centers such as radicals. PELDOR has been recognized as a valuable tool to approach structural questions in biological systems. In this manuscript, we demonstrate the value of distance measurements for differentiating competing structural models on the dimerization of the effector domain (ED) of the non-structural protein 1 (NS1) of the influenza A virus. Our results show NS1 to be well amenable to nanometer distance measurements by EPR, yielding high quality data. In combination with mutants perturbing protein dimerization and in silico prediction based on crystal structures, we can exclude one of two potential dimerization interfaces. Furthermore, our results lead to a viable hypothesis of a NS1 ED:ED interface which is flexible through rotation around the vector interconnecting the two native cysteines. These results prove the high value of pulse EPR as a complementary method for structural biology.

High-resolution crystal structure of spin labelled (T21R1) azurin from Pseudomonas aeruginosa: a challenging structural benchmark for in silico spin labelling algorithms

Background

EPR-based distance measurements between spin labels in proteins have become a valuable tool in structural biology. The direct translation of the experimental distances into structural information is however often impaired by the intrinsic flexibility of the spin labelled side chains. Different algorithms exist that predict the approximate conformation of the spin label either by using pre-computed rotamer libraries of the labelled side chain (rotamer approach) or by simply determining its accessible volume (accessible volume approach). Surprisingly, comparisons with many experimental distances have shown that both approaches deliver the same distance prediction accuracy of about 3 Å.

Results

Here, instead of comparing predicted and experimental distances, we test the ability of both approaches to predict the actual conformations of spin labels found in a new high-resolution crystal structure of spin labelled azurin (T21R1). Inside the crystal, the label is found in two very different environments which serve as a challenging test for the in silico approaches.

Conclusions

Our results illustrate why simple and more sophisticated programs lead to the same prediciton error. Thus, a more precise treatment of the complete environment of the label and also its interactions with the environment will be needed to increase the accuracy of in silico spin labelling algorithms.

The histone chaperones Vps75 and Nap1 form ring-like, tetrameric structures in solution

NAP-1 fold histone chaperones play an important role in escorting histones to and from sites of nucleosome assembly and disassembly. The two NAP-1 fold histone chaperones in budding yeast, Vps75 and Nap1, have previously been crystalized in a characteristic homodimeric conformation. In this study, a combination of small angle X-ray scattering, multi angle light scattering and pulsed electron–electron double resonance approaches were used to show that both Vps75 and Nap1 adopt ring-shaped tetrameric conformations in solution. This suggests that the formation of homotetramers is a common feature of NAP-1 fold histone chaperones. The tetramerisation of NAP-1 fold histone chaperones may act to shield acidic surfaces in the absence of histone cargo thus providing a ‘self-chaperoning’ type mechanism.

Perturbations of the straight transmembrane α-helical structure of the amyloid precursor protein affect its processing by γ-secretase

The amyloid precursor protein (APP) is a widely expressed type I transmembrane (TM) glycoprotein present at the neuronal synapse. The proteolytic cleavage by γ-secretase of its C-terminal fragment produces amyloid-β (Aβ) peptides of different lengths, the deposition of which is an early indicator of Alzheimer disease. At present, there is no consensus on the conformation of the APP-TM domain at the biological membrane. Although structures have been determined by NMR in detergent micelles, their conformation is markedly different. Here we show by using molecular simulations that the APP-TM region systematically prefers a straight α-helical conformation once embedded in a membrane bilayer. However, APP-TM is highly flexible, and its secondary structure is strongly influenced by the surrounding lipid environment, as when enclosed in detergent micelles. This behavior is confirmed when analyzing in silico the atomistic APP-TM population observed by residual dipolar couplings and double electron-electron resonance spectroscopy. These structural and dynamic features are critical in the proteolytic processing of APP by the γ-secretase enzyme, as suggested by a series of Gly(700) mutants. Affecting the hydration and flexibility of APP-TM, these mutants invariantly show an increase in the production of Aβ38 compared with Aβ40 peptides, which is reminiscent of the effect of γ-secretase modulators inhibitors.

Paramagnetic metal ions in pulsed ESR distance distribution measurements

The use of pulsed electron spin resonance (ESR) to measure interspin distance distributions has advanced biophysical research. The three major techniques that use pulsed ESR are relaxation rate based distance measurements, double quantum coherence (DQC), and double electron electron resonance (DEER). Among these methods, the DEER technique has become particularly popular largely because it is easy to implement on commercial instruments and because programs are available to analyze experimental data. Researchers have widely used DEER to measure the structure and conformational dynamics of molecules labeled with the methanethiosulfonate spin label (MTSSL). Recently, researchers have exploited endogenously bound paramagnetic metal ions as spin probes as a way to determine structural constraints in metalloproteins. In this context Cu(2+) has served as a useful paramagnetic metal probe at X-band for DEER based distance measurements. Sample preparation is simple, and a coordinated-Cu(2+) ion offers limited spatial flexibility, making it an attractive probe for DEER experiments. On the other hand, Cu(2+) has a broad absorption ESR spectrum at low temperature, which leads to two potential complications. First, the Cu(2+)-based DEER time domain data has lower signal to noise ratio compared with MTSSL. Second, accurate distance distribution analysis often requires high-quality experimental data at different external magnetic fields or with different frequency offsets. In this Account, we summarize characteristics of Cu(2+)-based DEER distance distribution measurements and data analysis methods. We highlight a novel application of such measurements in a protein-DNA complex to identify the metal ion binding site and to elucidate its chemical mechanism of function. We also survey the progress of research on other metal ions in high frequency DEER experiments.

PELDOR in rotationally symmetric homo-oligomers

Nanometre distance measurements by pulsed electron-electron double resonance (PELDOR) spectroscopy have become an increasingly important tool in structural biology. The theoretical underpinning of the experiment is well defined for systems containing two nitroxide spin-labels (spin pairs); however, recently experiments have been reported on homo-oligomeric membrane proteins consisting of up to eight spin-labelled monomers. We have explored the theory behind these systems by examining model systems based on multiple spins arranged in rotationally symmetric polygons. The results demonstrate that with a rising number of spins within the test molecule, increasingly strong distortions appear in distance distributions obtained from an analysis based on the simple spin pair approach. These distortions are significant over a range of system sizes and remain so even when random errors are introduced into the symmetry of the model. We present an alternative approach to the extraction of distances on such systems based on a minimisation that properly treats multi-spin correlations. We demonstrate the utility of this approach on a spin-labelled mutant of the heptameric Mechanosensitive Channel of Small Conductance of E. coli.

RosettaEPR: rotamer library for spin label structure and dynamics

An increasingly used parameter in structural biology is the measurement of distances between spin labels bound to a protein. One limitation to these measurements is the unknown position of the spin label relative to the protein backbone. To overcome this drawback, we introduce a rotamer library of the methanethiosulfonate spin label (MTSSL) into the protein modeling program Rosetta. Spin label rotamers were derived from conformations observed in crystal structures of spin labeled T4 lysozyme and previously published molecular dynamics simulations. Rosetta’s ability to accurately recover spin label conformations and EPR measured distance distributions was evaluated against 19 experimentally determined MTSSL labeled structures of T4 lysozyme and the membrane protein LeuT and 73 distance distributions from T4 lysozyme and the membrane protein MsbA. For a site in the core of T4 lysozyme, the correct spin label conformation (Χ₁ and Χ₂) is recovered in 99.8% of trials. In surface positions 53% of the trajectories agree with crystallized conformations in Χ₁ and Χ₂. This level of recovery is on par with Rosetta performance for the 20 natural amino acids. In addition, Rosetta predicts the distance between two spin labels with a mean error of 4.4 Å. The width of the experimental distance distribution, which reflects the flexibility of the two spin labels, is predicted with a mean error of 1.3 Å. RosettaEPR makes full-atom spin label modeling available to a wide scientific community in conjunction with the powerful suite of modeling methods within Rosetta.

Technological advances in site-directed spin labeling of proteins

Molecular flexibility over a wide time range is of central importance to the function of many proteins, both soluble and membrane. Revealing the modes of flexibility, their amplitudes, and time scales under physiological conditions is the challenge for spectroscopic methods, one of which is site-directed spin labeling EPR (SDSL-EPR). Here we provide an overview of some recent technological advances in SDSL-EPR related to investigation of structure, structural heterogeneity, and dynamics of proteins. These include new classes of spin labels, advances in measurement of long range distances and distance distributions, methods for identifying backbone and conformational fluctuations, and new strategies for determining the kinetics of protein motion.

Conformational dynamics and distribution of nitroxide spin labels

Long-range distance measurements based on paramagnetic relaxation enhancement (PRE) in NMR, quantification of surface water dynamics near biomacromolecules by Overhauser dynamic nuclear polarization (DNP) and sensitivity enhancement by solid-state DNP all depend on introducing paramagnetic species into an otherwise diamagnetic NMR sample. The species can be introduced by site-directed spin labeling, which offers precise control for positioning the label in the sequence of a biopolymer. However, internal flexibility of the spin label gives rise to dynamic processes that potentially influence PRE and DNP behavior and leads to a spatial distribution of the electron spin even in solid samples. Internal dynamics of spin labels and their static conformational distributions have been studied mainly by electron paramagnetic resonance spectroscopy and molecular dynamics simulations, with a large body of results for the most widely applied methanethiosulfonate spin label MTSL. These results are critically discussed in a unifying picture based on rotameric states of the group that carries the spin label. Deficiencies in our current understanding of dynamics and conformations of spin labeled groups and of their influence on NMR observables are highlighted and directions for further research suggested.

mtsslSuite: In silico spin labelling, trilateration and distance-constrained rigid body docking in PyMOL

Nanometer distance measurements based on electron paramagnetic resonance methods in combination with site-directed spin labelling are powerful tools for the structural analysis of macromolecules. The software package mtsslSuite provides scientists with a set of tools for the translation of experimental distance distributions into structural information. The package is based on the previously published mtsslWizard software for in silico spin labelling. The mtsslSuite includes a new version of MtsslWizard that has improved performance and now includes additional types of spin labels. Moreover, it contains applications for the trilateration of paramagnetic centres in biomolecules and for rigid-body docking of subdomains of macromolecular complexes. The mtsslSuite is tested on a number of challenging test cases and its strengths and weaknesses are evaluated.

The architecture of EssB, an integral membrane component of the type VII secretion system

The membrane-bound EssB is an integral and essential component of the bacterial type VII secretion system that can contribute to pathogenicity. The architecture of Geobacillus thermodenitrificans EssB has been investigated by combining crystallographic and EPR spectroscopic methods. The protein forms a dimer that straddles the cytoplasmic membrane. A helical fold is observed for the C-terminal segment, which is positioned on the exterior of the membrane. This segment contributes most to dimer formation. The N-terminal segment displays a structure related to the pseudokinase fold and may contribute to function by recognizing substrates or secretion system partners. The remaining part of EssB may serve as an anchor point for the secretion apparatus, which is embedded in the cytoplasmic membrane with the C-terminal domain protruding out to interact with partner proteins or components of peptidoglycan.

Conformational state of the MscS mechanosensitive channel in solution revealed by pulsed electron-electron double resonance (PELDOR) spectroscopy

The heptameric mechanosensitive channel of small conductance (MscS) provides a critical function in Escherichia coli where it opens in response to increased bilayer tension. Three approaches have defined different closed and open structures of the channel, resulting in mutually incompatible models of gating. We have attached spin labels to cysteine mutants on key secondary structural elements specifically chosen to discriminate between the competing models. The resulting pulsed electron-electron double resonance (PELDOR) spectra matched predicted distance distributions for the open crystal structure of MscS. The fit for the predictions by structural models of MscS derived by other techniques was not convincing. The assignment of MscS as open in detergent by PELDOR was unexpected but is supported by two crystal structures of spin-labeled MscS. PELDOR is therefore shown to be a powerful experimental tool to interrogate the conformation of transmembrane regions of integral membrane proteins.

xPyder: a PyMOL plugin to analyze coupled residues and their networks in protein structures

A versatile method to directly identify and analyze short- or long-range coupled or communicating residues in a protein conformational ensemble is of extreme relevance to achieve a complete understanding of protein dynamics and structural communication routes. Here, we present xPyder, an interface between one of the most employed molecular graphics systems, PyMOL, and the analysis of dynamical cross-correlation matrices (DCCM). The approach can also be extended, in principle, to matrices including other indexes of communication propensity or intensity between protein residues, as well as the persistence of intra- or intermolecular interactions, such as those underlying protein dynamics. The xPyder plugin for PyMOL 1.4 and 1.5 is offered as Open Source software via the GPL v2 license, and it can be found, along with the installation package, the user guide, and examples, at http://linux.btbs.unimib.it/xpyder/.

Characterization of Protein Conformational Changes with Sparse Spin-Label Distance Constraints

The combination of site-directed spin labeling with pulse EPR distance measurements can provide a moderate number of distance constraints on the nanometer length scale for proteins in different states. By adapting an existing algorithm (Zheng, W.; Brooks, B. R. Biophys. J. 2006, 90, 4327) to the problem, we address the question to what extent conformational change can be characterized when the protein structure is known for one of the states, whereas only a sparse set of distance constraints between spin labels is available for the other state. We find that the type and general direction of the conformational change can be recognized, while the amplitude may be uncertain.