Solution structure of the loops of bacteriorhodopsin closely resembles the crystal structure

Bacteriorhodopsin is one of very few transmembrane proteins for which high resolution structures have been solved. The structure shows a bundle of seven helices connected by six turns. Some turns in proteins are stabilized by short range interactions and can behave as small domains. These observations suggest that peptides containing the sequence of the turns in a membrane protein such as bacteriorhodopsin may form stable turn structures in solution. To test this hypothesis, we determined the solution structure of three peptides each containing the sequence of one of the turns in bacteriorhodopsin. The solution structures of the peptides closely resemble the structures of the corresponding turns in the high resolution structures of the intact protein.     

Structures of the transmembrane helices of the G-protein coupled receptor, rhodopsin.

An hypothesis is tested that individual peptides corresponding to the transmembrane helices of the membrane protein, rhodopsin, would form helices in solution similar to those in the native protein. Peptides containing the sequences of helices 1, 4 and 5 of rhodopsin were synthesized. Two peptides, with overlapping sequences at their termini, were synthesized to cover each of the helices. The peptides from helix 1 and helix 4 were helical throughout most of their length. The N- and C-termini of all the peptides were disordered and proline caused opening of the helical structure in both helix 1 and helix 4. The peptides from helix 5 were helical in the middle segment of each peptide, with larger disordered regions in the N- and C-termini than for helices 1 and 4. These observations show that there is a strong helical propensity in the amino acid sequences corresponding to the transmembrane domain of this G-protein coupled receptor. In the case of the peptides from helix 4, it was possible to superimpose the structures of the overlapping sequences to produce a construct covering the whole of the sequence of helix 4 of rhodopsin. As similar superposition for the peptides from helix 1 also produced a construct, but somewhat less successfully because of the disordering in the region of sequence overlap. This latter problem was more severe for helix 5 and therefore a single peptide was synthesized for the entire sequence of this helix, and its structure determined. It proved to be helical throughout. Comparison of all these structures with the recent crystal structure of rhodopsin revealed that the peptide structures mimicked the structures seen in the whole protein. Thus similar studies of peptides may provide useful information on the secondary structure of other transmembrane proteins built around helical bundles.     

Detergent-free membrane protein crystallization.

A comprehensive understanding of structure-function relationships of proteins requires their structures to be elucidated to high resolution. With most membrane proteins this has not been accomplished so far, mainly because of their notoriously poor crystallizability. Here we present a completely detergent-free procedure for the incorporation of a native purple membrane into a monoolein-based lipidic cubic phase, and subsequent crystallization of three-dimensional bacteriorhodopsin crystals therein. These crystals exhibit comparable X-ray diffraction quality and mosaicity, and identical crystal habit and space group to those of bacteriorhodopsin crystals that are grown from detergent-solubilized protein in cubic phase.     

Optimal bundling of transmembrane helices using sparse distance constraints

We present a two-step approach to modeling the transmembrane spanning helical bundles of integral membrane proteins using only sparse distance constraints, such as those derived from chemical cross-linking, dipolar EPR and FRET experiments. In Step 1, using an algorithm, we developed, the conformational space of membrane protein folds matching a set of distance constraints is explored to provide initial structures for local conformational searches. In Step 2, these structures refined against a custom penalty function that incorporates both measures derived from statistical analysis of solved membrane protein structures and distance constraints obtained from experiments. We begin by describing the statistical analysis of the solved membrane protein structures from which the theoretical portion of the penalty function was derived. We then describe the penalty function, and, using a set of six test cases, demonstrate that it is capable of distinguishing helical bundles that are close to the native bundle from those that are far from the native bundle. Finally, using a set of only 27 distance constraints extracted from the literature, we show that our method successfully recovers the structure of dark-adapted rhodopsin to within 3.2 A of the crystal structure.     

Solid-state NMR structures of integral membrane proteins

Abstract

Solid-state NMR is unique for its ability to obtain three-dimensional structures and to measure atomic-resolution structural and dynamic information for membrane proteins in native lipid bilayers. An increasing number and complexity of integral membrane protein structures have been determined by solid-state NMR using two main methods. Oriented sample solid-state NMR uses macroscopically aligned lipid bilayers to obtain orientational restraints that define secondary structure and global fold of embedded peptides and proteins and their orientation and topology in lipid bilayers. Magic angle spinning (MAS) solid-state NMR uses unoriented rapidly spinning samples to obtain distance and torsion angle restraints that define tertiary structure and helix packing arrangements. Details of all current protein structures are described, highlighting developments in experimental strategy and other technological advancements. Some structures originate from combining solid- and solution-state NMR information and some have used solid-state NMR to refine X-ray crystal structures. Solid-state NMR has also validated the structures of proteins determined in different membrane mimetics by solution-state NMR and X-ray crystallography and is therefore complementary to other structural biology techniques. By continuing efforts in identifying membrane protein targets and developing expression, isotope labelling and sample preparation strategies, probe technology, NMR experiments, calculation and modelling methods and combination with other techniques, it should be feasible to determine the structures of many more membrane proteins of biological and biomedical importance using solid-state NMR. This will provide three-dimensional structures and atomic-resolution structural information for characterising ligand and drug interactions, dynamics and molecular mechanisms of membrane proteins under physiological lipid bilayer conditions.     

Characterization and crystal packing of three-dimensional bacteriorhodopsin crystals

The three-dimensional crystals of the integral membrane protein bacteriorhodopsin have been characterized by X-ray diffraction and freeze-fracture electron microscopy: the needle-like form A crystals belong to space group P 1 (pseudohexagonal) with seven molecules per crystallographic unit cell forming one turn of a non-crystallographic helix. The probable arrangement of the bacteriorhodopsin molecules is derived from freeze-fracture electron micrographs and chromophore orientation. Membrane-like structures are not present. The same helices of bacteriorhodopsin molecules found in crystal form A also make up the cube-like crystal form B. They are now arranged in all three mutually perpendicular directions. These cubes are always highly disordered, since the unit cell length corresponds to 6.7 molecules of the 7-fold helix. Very often, conversion of bacteriorhodopsin from the three-dimensional crystals into filamentous material occurs.     

Structure and dynamics of dark-state bovine rhodopsin revealed by chemical cross-linking and high-resolution mass spectrometry

Recent work using chemical cross-linking to define interresidue distance constraints in proteins has shown that these constraints are useful for testing tertiary structural models. We applied this approach to the G-protein-coupled receptor bovine rhodopsin in its native membrane using lysine- and cysteine-targeted bifunctional cross-linking reagents. Cross-linked proteolytic peptides of rhodopsin were identified by combined liquid chromatography and FT-ICR mass spectrometry with automated data-reduction and assignment software. Tandem mass spectrometry was used to verify cross-link assignments and locate the exact sites of cross-link attachment. Cross-links were observed to form between 10 pairs of residues in dark-state rhodopsin. For each pair, cross-linkers with a range of linker lengths were tested to determine an experimental distance-of-closest-approach (DCA) between reactive side-chain atoms. In all, 28 cross-links were identified using seven different cross-linking reagents. Molecular mechanics procedures were applied to published crystal structure data to calculate energetically achievable theoretical DCAs between reactive atoms without altering the position of the protein backbone. Experimentally measured DCAs are generally in good agreement with the theoretical DCAs. However, a cross-link between C316 and K325 in the C-terminal region cannot be rationalized by DCA simulations and suggests that backbone reorientation relative to the crystal coordinates occurs on the timescale of cross-linking reactions. Biochemical and spectroscopic data from other studies have found that the C-terminal region is highly mobile in solution and not fully represented by X-ray crystallography data. Our results show that chemical cross-linking can provide reliable three-dimensional structural information and insight into local conformational dynamics in a membrane protein.     

Calculation of protein conformations by proton-proton distance constraints. A new efficient algorithm

We have developed a method to determine the three-dimensional structure of a protein molecule from such a set of distance constraints as can be determined by nuclear magnetic resonance studies. The currently popular methods for distance geometry based on the use of the metric matrix are applicable only to small systems. The method developed here is applicable to large molecules, such as proteins, with all atoms treated explicitly. This method works in the space of variable dihedral angles and determines a three-dimensional structure by minimization of a target function. We avoid difficulties hitherto inherent in this type of approach by two new devices: the use of variable target functions; and a method of rapid calculation of the gradient of the target functions. The method is applied to the determination of the structures of a small globular protein, bovine pancreatic trypsin inhibitor, from several artificial sets of distance constraints extracted from the X-ray crystal structure of this molecule. When a good set of constraints was available for both short- and long-range distances, the crystal structure was regenerated nearly exactly. When some ambiguities, such as those expected in experimental information, are allowed, the protein conformation can be determined up to a few local deformations. These ambiguities are mainly associated with the low resolving power of the short-range information.     

Three-dimensional solution structure of Ca(2+)-loaded porcine calbindin D9k determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of native, intact porcine calbindin D9k has been determined by distance geometry and restrained molecular dynamics calculations using distance and dihedral angle constraints obtained from 1H NMR spectroscopy. The protein has a well-defined global fold consisting of four helices oriented in a pairwise antiparallel manner such that two pairs of helix-loop-helix motifs (EF-hands) are joined by a linker segment. The two EF-hands are further coupled through a short beta-type interaction between the two Ca(2+)-binding loops. Overall, the structure is very similar to that of the highly homologous native, minor A form of bovine calbindin D9k determined by X-ray crystallography [Szebenyi, D. M. E., & Moffat, K. (1986) J. Biol. Chem. 261, 8761-8776]. A model structure built from the bovine calbindin D9k crystal structure shows several deviations larger than 2 A from the experimental distance constraints for the porcine protein. These structural differences are efficiently removed by subjecting the model structure to the experimental distance and dihedral angle constraints in a restrained molecular dynamics protocol, thereby generating a model that is very similar to the refined distance geometry derived structures. The N-terminal residues of the intact protein that are absent in the minor A form appear to be highly flexible and do not influence the structure of other regions of the protein. This result is important because it validates the conclusions drawn from the wide range of studies that have been carried out on minor A forms rather than the intact calbindin D9k.     

Membrane Protein Structural Validation by Oriented Sample Solid-State NMR: Diacylglycerol Kinase

The validation of protein structures through functional assays has been the norm for many years. Functional assays perform this validation for water-soluble proteins very well, but they need to be performed in the same environment as that used for the structural analysis. This is difficult for membrane proteins that are often structurally characterized in detergent environments, although functional assays for these proteins are most frequently performed in lipid bilayers. Because the structure of membrane proteins is known to be sensitive to the membrane mimetic environment, such functional assays are appropriate for validating the protein construct, but not the membrane protein structure. Here, we compare oriented sample solid-state NMR spectral data of diacylglycerol kinase previously published with predictions of such data from recent structures of this protein. A solution NMR structure of diacylglycerol kinase has been obtained in detergent micelles and three crystal structures have been obtained in a monoolein cubic phase. All of the structures are trimeric with each monomer having three transmembrane and one amphipathic helices. However, the solution NMR structure shows typical perturbations induced by a micelle environment that is reflected in the predicted solid-state NMR resonances from the structural coordinates. The crystal structures show few such perturbations, especially for the wild-type structure and especially for the monomers that do not have significant crystal contacts. For these monomers the predicted and observed data are nearly identical. The thermostabilized constructs do show more perturbations, especially the A41C mutation that introduces a hydrophilic residue into what would be the middle of the lipid bilayer inducing additional hydrogen bonding between trimers. These results demonstrate a general technique for validating membrane protein structures with minimal data obtained from membrane proteins in liquid crystalline lipid bilayers by oriented sample solid-state NMR.     

Analysis of secondary structure in proteins by chemical cross-linking coupled to MS

Chemical cross-linking is an attractive technique for the study of the structure of protein complexes due to its low sample consumption and short analysis time. Furthermore, distance constraints obtained from the identification of cross-linked peptides by MS can be used to construct and validate protein models. If a sufficient number of distance constraints are obtained, then determining the secondary structure of a protein can allow inference of the protein's fold. In this work, we show how the distance constraints obtained from cross-linking experiments can identify secondary structures within the protein sequence. Molecular modeling of alpha helices and beta sheets reveals that each secondary structure presents different cross-linking possibilities due to the topological distances between reactive residues. Cross-linking experiments performed with amine reactive cross-linkers with model alpha helix containing proteins corroborated the molecular modeling predictions. The cross-linking patterns established here can be extended to other cross-linkers with known lengths for the determination of secondary structures in proteins.     

High-resolution solution structure of the 18 kDa substrate-binding domain of the mammalian chaperone protein Hsc70.

The three-dimensional structure for the substrate-binding domain of the mammalian chaperone protein Hsc70 of the 70 kDa heat shock class (HSP70) is presented. This domain includes residues 383-540 (18 kDa) and is necessary for the binding of the chaperone with substrate proteins and peptides. The high-resolution NMR solution structure is based on 4150 experimental distance constraints leading to an average root-mean-square precision of 0.38 A for the backbone atoms and 0.76 A for all atoms in the beta-sandwich sub-domain. The protein is observed to bind residue Leu539 in its hydrophobic substrate-binding groove by intramolecular interaction. The position of a helical latch differs dramatically from what is observed in the crystal and solution structures of the homologous prokaryotic chaperone DnaK. In the Hsc70 structure, the helix lies in a hydrophobic groove and is anchored by a buried salt-bridge. Residues involved in this salt-bridge appear to be important for the allosteric functioning of the protein. A mechanism for interdomain allosteric modulation of substrate-binding is proposed. It involves large-scale movements of the helical domain, redefining the location of the hinge area that enables such motions.     

Time-resolved microspectroscopy on a single crystal of bacteriorhodopsin reveals lattice-induced differences in the photocycle kinetics

The determination of the intermediate state structures of the bacteriorhodopsin photocycle has lead to an unprecedented level of understanding of the catalytic process exerted by a membrane protein. However, the crystallographic structures of the intermediate states are only relevant if the working cycle is not impaired by the crystal lattice. Therefore, we applied visible and Fourier transform infrared spectroscopy (FTIR) microspectroscopy with microsecond time resolution to compare the photoreaction of a single bacteriorhodopsin crystal to that of bacteriorhodopsin residing in the native purple membrane. The analysis of the FTIR difference spectra of the resolved intermediate states reveals great similarity in structural changes taking place in the crystal and in PM. However, the kinetics of the photocycle are significantly altered in the three-dimensional crystal as compared to PM. Strikingly, the L state decay is accelerated in the crystal, whereas the M decay is delayed. The physical origin of this deviation and the implications for trapping of intermediate states are discussed. As a methodological advance, time-resolved step-scan FTIR spectroscopy on a single protein crystal is demonstrated for the first time which may be used in the future to gauge the functionality of other crystallized proteins with the molecular resolution of vibrational spectroscopy.     

Folding patterns of porin and bacteriorhodopsin.

Porin spans the outer membrane of Escherichia coli with most of the protein embedded within the membrane. It lacks pronounced hydrophobic domains and consists predominantly of beta-pleated sheet. These observations require the accommodation of polar and ionizable residues in an environment that has a low dielectric constant. Owing to a currently limited understanding of the constraints governing membrane protein structure, a minimal approach to structure prediction is proposed that identifies segments causing polypeptides to reverse their direction (turn identification). The application of this procedure avoids hydrophobicity parameters and yields a model of porin which is in good agreement with all experimental data available. The presence of polar and ionizable residues within membrane boundaries implies a dense (saturating) network of hydrogen bond donor and acceptor groups. Application to a paradigm of hydrophobic membrane proteins, bacteriorhodopsin, reveals a pattern consistent with its alpha-helical folding. The postulated structure includes significantly more polar residues in the membrane domain than have been assumed previously, suggesting that there are also hydrogen bonding networks in bacteriorhodopsin. Extensive networks permeating protein interior and surfaces would explain the extraordinary stability and the tight interactions between functional units in the formation of crystalline arrays of both proteins.     

A novel measure characterized by a polar energy surface approximation for recognition and classification of transmembrane protein structures

A new theoretical method has been developed for recognition and classification of membrane proteins. The method is based on computation of a polar energy surface that can reveal characteristic interaction patterns for individual helices even if crystal or NMR structure coordinates are not available. A protein with N transmembrane helices is described as a set of N vectors that are derived from a Fourier analysis of this polar energy surface computed for each helix. We then derive a polarity difference score (PDS) for any two proteins computed as the root mean square deviation between the respective vector coordinate sets. The score was found to correlate with the degree of structural similarity between the following three protein families for which tertiary structures have been determined: bacteriorhodopsin, rhodopsin, and the cytochrome c oxidase III subunit.     

Comparison of NMR and crystal structures of membrane proteins and computational refinement to improve model quality

Membrane proteins are challenging to study and restraints for structure determination are typically sparse or of low resolution because the membrane environment that surrounds them leads to a variety of experimental challenges. When membrane protein structures are determined by different techniques in different environments, a natural question is “which structure is most biologically relevant?” Towards answering this question, we compiled a dataset of membrane proteins with known structures determined by both solution NMR and X‐ray crystallography. By investigating differences between the structures, we found that RMSDs between crystal and NMR structures are below 5 Å in the membrane region, NMR ensembles have a higher convergence in the membrane region, crystal structures typically have a straighter transmembrane region, have higher stereo‐chemical correctness, and are more tightly packed. After quantifying these differences, we used high‐resolution refinement of the NMR structures to mitigate them, which paves the way for identifying and improving the structural quality of membrane proteins.     

Homonuclear three-dimensional NOE-NOE nuclear magnetic resonance/spectra for structure determination of proteins in solution.

The solution structures of two proteins (CMTI-I, a trypsin inhibitor from Cucurbita maxima, and hisactophilin, an actin binding protein of 118 amino acids) have been determined based on the NOE data derived solely from the homonuclear 3D NOE-NOE magnetic resonance spectroscopy. Two different approaches for extraction of the structural information from the 3D NOE-NOE experiment were tested. One approach was based on the transformation of the 3D intensities into distance constraints. In the second, and more robust approach, the 3D NOE intensities were used directly in structure calculations, without the need to transform them into distance constraints. A new 2D potential function representing the 3D NOE-NOE intensity was developed and used in the simulated annealing protocol. For CMTI-I, a comparison between structures determined with the 3D NOE-NOE method and various 2D NOE approaches was carried out. The 3D data set allowed better definition of the structures than was previously possible with the 2D NOE procedures that used the isolated two-spin approximation to derive distance information.     

High-resolution solution structure of reduced French bean plastocyanin and comparison with the crystal structure of poplar plastocyanin

The three-dimensional solution structure of reduced (CuI) plastocyanin from French bean leaves has been determined by distance geometry and restrained molecular dynamics methods using constraints obtained from 1H n.m.r. (nuclear magnetic resonance) spectroscopy. A total of 1244 experimental constraints were used, including 1120 distance constraints, 103 dihedral angle constraints and 21 hydrogen bond constraints. Stereospecific assignments were made for 26 methylene groups and the methyls of 11 valines. Additional constraints on copper co-ordination were included in the restrained dynamics calculations. The structures are well defined with average atomic root-mean-square deviations from the mean of 0.45 A for all backbone heavy atoms and 1.08 A for side-chain heavy atoms. French bean plastocyanin adopts a beta-sandwich structure in solution that is similar to the X-ray structure of reduced poplar plastocyanin; the average atomic root-mean-square difference between 16 n.m.r. structures and the X-ray structure is 0.76 A for all backbone heavy atoms. The conformations of the side-chains that constitute the hydrophobic core of French bean plastocyanin are very well defined. Of 47 conserved residues that populate a single chi 1 angle in solution, 43 have the same rotamer in the X-ray structure. Many surface side-chains adopt highly preferred conformations in solution, although the 3J alpha beta coupling constants often indicate some degree of conformational averaging. Some surface side-chains are disordered in both the solution and crystal structures of plastocyanin. There is a striking correlation between measures of side-chain disorder in solution and side-chain temperature factors in the X-ray structure. Side-chains that form a distinctive acidic surface region, believed to be important in binding other electron transfer proteins, appear to be disordered. Fifty backbone amide protons form hydrogen bonds to carbonyls in more than 60% of the n.m.r. structures; 45 of these amide protons exchange slowly with solvent deuterons. Ten hydrogen bonds are formed between side-chain and backbone atoms, eight of which are correlated with decreased proton exchange. Of the 60 hydrogen bonds formed in French bean plastocyanin, 56 occur in the X-ray structure of the poplar protein; two of the missing hydrogen bonds are absent as a result of mutations. It appears that molecular dynamics refinement of highly constrained n.m.r. structures allows accurate prediction of the pattern of hydrogen bonding.     

Positioning membrane proteins by novel protein engineering and biophysical approaches

Membrane proteins are unique, in that they can function properly only when they are bound to cellular membranes in a distinct manner. Therefore, positioning of membrane proteins with respect to the membrane is required in addition to the three-dimensional structures in order to understand their detailed molecular mechanisms. Atomic-resolution structures of membrane proteins that have been determined to date provide the atom coordinates in arbitrary coordinate systems with no relation to the membrane and therefore provide little or no information on how the protein would interact with the membrane. This is especially true for peripheral membrane proteins, because they, unlike integral proteins, are devoid of well-defined hydrophobic transmembrane domains. Here, we present a novel technique for determination of the configuration of a protein-membrane complex that involves protein ligation, segmental isotope labeling, polarized infrared spectroscopy, membrane depth-dependent fluorescence quenching, and analytical geometry algorithms. We have applied this approach to determine the structure of a membrane-bound phospholipase A2. Our results provide an unprecedented structure of a membrane-bound protein in which the z-coordinate of each atom is the distance from the membrane center and therefore allows precise location of each amino acid relative to the membrane. Given the functional significance of the orientation and location of membrane-bound proteins with respect to the membrane, we propose to specify this structural feature as the "quinary" structure of membrane proteins.     

NMR studies of defensin antimicrobial peptides. 2. Three-dimensional structures of rabbit NP-2 and human HNP-1.

The solution structure of two homologous naturally occurring antimicrobial peptides, rabbit defensin NP-2 and human defensin HNP-1, have been determined by two-dimensional nuclear magnetic resonance spectroscopy, distance geometry, and restrained molecular dynamics calculations. The structure of these defensins consists of an antiparallel beta-sheet in a hairpin conformation, a short region of triple-stranded beta-sheet, several tight turns, and a loop region that has a well-defined local structure but with a global orientation that is not well-defined with respect to the rest of the molecule. The solution structures of these two peptides are compared with the solution and crystal structures of two other homologous defensins. The structures for the defensins are also compared with known structures of other naturally occurring antimicrobial peptides.