3D Modeling for TM Receptors: Algorithms and Validations

This paper outlines the basic strategy to build 3D models of transmembrane G-protein coupled receptors (GPCRs) starting from their amino acid sequences in a block-by-block manner: (i) locate possible TM helical fragments in the GPCR sequence; (ii) build 3D structures for these helices; (iii) arrange isolated helices across the membrane; (iv) calculate all pairwise helix-helix interactions; (v) assemble helical bundle(s); (vi) restore interhelical loops and N- and C-termini; and (vii) refine the entire 3D structure(s). Computer algorithms and preliminary results for most of the steps are discussed.     

Tertiary Motifs Revealed in Analyses of Higher-Order RNA Junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Relationship between the Interhelical Packing Angles and the Length of α-Helices in Proteins

Abstract—Mutual arrangement, or packing, of α-helices in proteins depends on several factors, but, tight packing and the chemical nature of the polypeptide chain are the most important. This study shows, for the first time, that the torsion packing angles between axes of α-helices depend on their length. A database of helical pairs formed by two connected and juxtaposed α-helices has been compiled using the Protein Data Bank. These helical pairs have been subdivided into four types: (1) 10474 pairs formed by long helices; (2) 3665 pairs in which the first α-helix is long and the second is short; (3) 3648 pairs in which the first α‑helix is short and the second is long; 4) 1895 pairs in which both helices are short. Analysis of the database showed that most helical pairs in which both the helices are long form α-hairpins having interhelical packing angles of Ω ≈ 20°. Most helical pairs in which one α-helix is long and the other is short or both helices are short form αα-corners having orthogonal (Ω ≈ –70°…–90°) or slanted (Ω ≈ –50°) packing of α-helices. The possible reasons for this relationship between interhelical angles (Ω) and the length of α-helices are discussed. These results are of great importance in protein modeling and prediction since they enable the determination of the mutual arrangement of α-helices in protein molecules.     

A comparative study of peptide models of the alpha-domain of alpha-lactalbumin, lysozyme, and alpha-lactalbumin/lysozyme chimeras allows the elucidation of critical factors that contribute to the ability to form stable partially folded states

alpha-Lactalbumin (alpha LA) forms a well-populated equilibrium molten globule state, while the homologous protein hen lysozyme does not. alpha LA is a two-domain protein and the alpha-domain is more structured in the molten globule state than is the beta-domain. Peptide models derived from the alpha-subdomain that contain the A, B, D, and 3(10) helices of alpha LA are capable of forming a molten globule state in the absence of the remainder of the protein. Here we report comparative studies of a peptide model derived from the same region of hen lysozyme and a set of chimeric alpha-lactalbumin--lysozyme constructs. Circular dichroism, dynamic light scattering, sedimentation equilibrium, and fluorescence experiments indicate that the lysozyme construct does not fold. Chimeric constructs were prepared to probe the origins of the difference in the ability of the two isolated subdomains to fold. The first consists of the A and B helices of alpha LA cross-linked to the D and C-terminal 3(10) helices of lysozyme. This construct is highly helical, while a second construct that contains the A and B helices of lysozyme cross-linked to the D and 3(10) helices of alpha LA does not fold. Furthermore, the disulfide cross-linked homodimer of the alpha LA AB peptide is helical, while the homodimer of the lysozyme AB peptide is unstructured. Thus, the AB helix region of alpha LA appears to have an intrinsic ability to form structure as long as some relatively nonspecific interactions can be made with other regions of the protein. Our studies show that the A and B helices plays a key role in the ability of the respective alpha-subdomains to fold.     

The global structure of the VS ribozyme

The VS ribozyme comprises five helical segments (II-VI) in a formal H shape, organized by two three-way junctions. It interacts with its stem-loop substrate (I) by tertiary interactions. We have determined the global shape of the 3-4-5 junction (relating helices III-V) by electrophoresis and FRET. Estimation of the dihedral angle between helices II and V electrophoretically has allowed us to build a model for the global structure of the complete ribozyme. We propose that the substrate is docked into a cleft between helices II and VI, with its loop making a tertiary interaction with that of helix V. This is consistent with the dependence of activity on the length of helix III. The scissile phosphate is well placed to interact with the probable active site of the ribozyme, the loop containing A730.     

Structural changes in bacteriorhodopsin during in vitro refolding from a partially denatured state

We report on the formation of the secondary and tertiary structure of bacteriorhodopsin during its in vitro refolding from an SDS-denatured state. We used the mobility of single spin labels in seven samples, attached at various locations to six of the seven helical segments to engineered cysteine residues, to follow coil-to-helix formation. Distance measurements obtained by spin dipolar quenching in six samples labeled at either the cytoplasmic or extracellular ends of pairs of helices revealed the time dependence of the recovery of the transmembrane helical bundle. The secondary structure in the majority of the helical segments refolds with a time constant of <100–140 ms. Recovery of the tertiary structure is achieved by sequential association of the helices and occurs in at least three distinct steps with time constants of 1), well below 1 s; 2), 3–4 s; and 3), 60–130 s (the latter depending on the helical pair). The slowest of these processes occurs in concert with recovery of the retinal chromophore.     

Analysis of Four-Way Junctions in RNA Structures

RNA secondary structures can be divided into helical regions composed of canonical Watson-Crick and related base pairs, as well as single-stranded regions such as hairpin loops, internal loops, and junctions. These elements function as building blocks in the design of diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze existing RNA four-way junctions in terms of base-pair interactions and 3D configurations. Specifically, we identify nine broad junction families according to coaxial stacking patterns and helical configurations. We find that helices within junctions tend to arrange in roughly parallel and perpendicular patterns and stabilize their conformations using common tertiary motifs such as coaxial stacking, loop-helix interaction, and helix packing interaction. Our analysis also reveals a number of highly conserved base-pair interaction patterns and novel tertiary motifs such as A-minor-coaxial stacking combinations and sarcin/ricin motif variants. Such analyses of RNA building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Comprehensive evaluation of solution nuclear magnetic resonance spectroscopy sample preparation for helical integral membrane proteins

The preparation of high quality samples is a critical challenge for the structural characterization of helical integral membrane proteins. Solving the structures of this diverse class of proteins by solution nuclear magnetic resonance spectroscopy (NMR) requires that well-resolved 2D ¹H/¹⁵N chemical shift correlation spectra be obtained. Acquiring these spectra demands the production of samples with high levels of purity and excellent homogeneity throughout the sample. In addition, high yields of isotopically enriched protein and efficient purification protocols are required. We describe two robust sample preparation methods for preparing high quality, homogeneous samples of helical integral membrane proteins. These sample preparation protocols have been combined with screens for detergents and sample conditions leading to the efficient production of samples suitable for solution NMR spectroscopy. We have examined 18 helical integral membrane proteins, ranging in size from approximately 9 kDa to 29 kDa with 1–4 transmembrane helices, originating from a number of bacterial and viral genomes. 2D ¹H/¹⁵N chemical shift correlation spectra acquired for each protein demonstrate well-resolved resonances, and >90% detection of the predicted resonances. These results indicate that with proper sample preparation, high quality solution NMR spectra of helical integral membrane proteins can be obtained greatly enhancing the probability for structural characterization of these important proteins.     

Nonnative helical motif in a chaperone-bound protein fragment

The effect of cotranslationally active chaperones on the conformation of incomplete protein chains is poorly understood. The secondary structure of a 77-residue chaperone-bound N-terminal protein fragment corresponding to the first five helices (A-E) of apomyoglobin (apoMb_1-77) is investigated here at the residue-specific level by multidimensional NMR. The substrate-binding domain of DnaK, DnaK-β, is employed as a chaperone model. By taking advantage of the improved spectral quality resulting from chaperone deuteration, we find that DnaK-β-bound apoMb_1-77 displays a region of nonnative helicity at residues away from the main chaperone binding site. The nonnative structural motif comprises portions of the native D and E helices and has similar characteristics to the reported nonnative DE helical region of acid-unfolded full-length apoMb. Upon incorporation of the missing C-terminal amino acids, a structural kink develops between residues 56 and 57, and two separate native D and E helices are generated. This work highlights, for the first time to our knowledge, the presence of a nonnative helical motif in a large chaperone-bound protein fragment under physiologically relevant solution conditions.     

Proline residues link the active site to transmembrane domain movements in human nucleoside triphosphate diphosphohydrolase 3 (NTPDase3)

The active sites of the membrane-bound nucleoside triphosphate diphosphohydrolases (NTPDases) regulate and are regulated by coordinated and spatially distant movements of their transmembrane helices, modulating enzyme activity, and substrate specificity. Using site-directed mutagenesis, the roles of the conserved proline residues (N-terminal: P52 and P53; C-terminal: P472, P476, P481, P484, and P485) of human NTPDase3, located in the “linker regions” that connect the N- and C-terminal transmembrane helices with the extracellular active site, were examined. Single cysteine substitutions were strategically placed in the transmembrane domain (N-terminal helix: V42C; C-terminal helix: G489C) to serve as cross-linking “sensors” of helical interactions. These “sensor” background mutant proteins (V42C and G489C NTPDase3) are enzymatically active and are cross-linked by copper phenanthroline less efficiently in the presence of adenosine triphosphate (ATP). Proline to alanine substitutions at P53, P481, P484, and P485 in the V42C background, as well as P53, P481, and P484 in the G489C background, exhibited decreased nucleotidase activities. More importantly, alanine substitutions at P53 and P481 in the V42C background and P481 in the G489C background no longer exhibited the ATP-induced decrease in transmembrane cross-linking efficiency. Interestingly, the P485A mutation abolished oxidative cross-linking at G489C both in the presence and absence of ATP. Taken together, these results suggest a role for proline residues 53 and 481 in the linker regions of human NTPDase3 for coupling nucleotide binding at the enzyme active site to movements and/or rearrangements of the transmembrane helices necessary for optimal nucleotide hydrolysis.     

Tuning the beta-turn segment in designed peptide beta-hairpins: construction of a stable type I' beta-turn nucleus and hairpin-helix transition promoting segments

Designed octapeptides Boc-Leu-Val-Val-Aib-(D)Xxx-Leu-Val-Val-OMe ((D)Xxx = (D)Ala, 3a;(D)Val, 3c and (D)Pro, 5a) and Boc-Leu-Phe-Val-Aib-(D)Ala-Leu-Phe-Val-OMe (3b) have been investigated to construct models of a stable type I' beta-turn nucleated hairpin and to generate systems for investigating helix-hairpin conformational transitions. Peptide 5a, which contains a central Aib-(D)Pro segment, is shown to adopt a stable type I' beta-turn nucleated hairpin structure, stabilized by four cross-strand hydrogen bonds. The stability of the structure in diverse solvents is established by the observation of all diagnostic NOEs expected in a beta-hairpin conformation. Replacement of (D)Pro5 by (D)Ala/(D)Val (3a-c) results in sequences that form beta-hairpins in hydrogen bonding solvents like CD(3)OH and DMSO-d(6). However, in CDCl(3) evidence for population of helical conformations is obtained. Peptide 6b (Boc-Leu-Phe-Val-Aib-Aib-Leu-Phe-Val-OMe), which contains a centrally positioned Aib-Aib segment, provides a clear example of a system, which exhibits a helical conformation in CDCl(3) and a significant population of both helices and hairpins in CD(3)OH and DMSO-d(6). The coexistence of multiple conformations is established by the simultaneous observation of diagnostic NOEs. Control over stereochemistry of the central beta-turn permits generation of models for robust beta-hairpins and also for the construction of systems that may be used to probe helix-hairpin conformational transitions.     

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K+-ATPase

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center ofRhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na⁺, K⁺-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na⁺, K⁺-ATPase is discussed in the final part of the work.     

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K+-ATPase

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center ofRhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na⁺, K⁺-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na⁺, K⁺-ATPase is discussed in the final part of the work.     

Identification and location of alpha-helices in mammalian cytochromes P450

A model of the alpha-helical structure of mammalian cytochromes P450 is proposed. The location and sequence of alpha-helices in mammalian cytochromes P450 were predicted from their homology with those of cytochrome P450cam, and these sequences were generally confirmed as helical in nature by using a secondary structure prediction method. These analyses were applied to 26 sequences in 6 gene families of cytochrome P450. Mammalian cytochromes P450 consist of approximately 100 amino acid residues more than cytochrome P450cam. This difference was accounted for by three major areas of insertion: (1) at the N-terminus, (2) between helices C and D and between helices D and E, and (3) between helices J and K. Insertion 1 has been suggested by others as a membrane anchoring sequence, but the apparent insertions at 2 and 3 are novel observations; it is suggested that they may be involved in the binding of cytochrome P450 reductase. Only the mitochondrial cytochrome P450 family appeared to show a major variation from this pattern, as insertion 2 was absent, replaced by an insertion between helices G and H and between helices H and I. This may reflect the difference in electron donor proteins that bind to members of this cytochrome P450 family. Other than these differences the model of mammalian cytochromes P450 proposed maintains the general structure of cytochrome P450cam as determined by its alpha-helical composition.     

A family of double helices of alternating poly(gamma-benzyl-D-L-glutamate), a stereochemical model for gramicidin A.

Poly(γ-benzyl-d-l-glutamate) with strict alternation of l and d residues is found to exist, in addition to the α_DL and π_DL^4.4 helical structures already described (Heitz et al., 1975a), in four more helical structures. Models based on double helices made of antiparallel strands are proposed for all four structures, based on infrared, X-ray and electron diffraction data. These double helices are, like the single-stranded π_DL helices, specific to polypeptides with a strict stereosequence of alternating l and d residues. The diameter of the helical core of three of these helices appears to depend on the dimensions of the solvent molecules. Conformational angles (located in the β regions) and atomic co-ordinates determined by conformational energy analysis are given for the four structures. Experimental conditions used to obtain these helices, and to induce transconformations between the various helical structures of PBd-lG are described. The present investigations on PBd-lG help to make more precise the structure and geometry of models proposed (Veatch et al., 1974) for the antibiotic gramicidin A.     

Conformational features of a hexapeptide model Ac-TGAAKA-NH2 corresponding to a hydrated alpha helical segment from glyceraldehyde 3-phosphate dehydrogenase: implications for the role of turns in helix folding

Recent analysis of alpha helices in protein crystal structures, available in literature, revealed hydrated alpha helical segments in which, water molecule breaks open helix 5-->1 hydrogen bond by inserting itself, hydrogen bonds to both C=O and NH groups of helix hydrogen bond without disrupting the helix hydrogen bond, and hydrogen bonds to either C=O or NH of helix hydrogen bond. These hydrated segments display a variety of turn conformations and are thought to be 'folding intermediates' trapped during folding-unfolding of alpha helices. A role for reverse turns is implicated in the folding of alpha helices. We considered a hexapeptide model Ac-1TGAAKA6-NH2 from glyceraldehyde 3-phosphate dehydrogenase, corresponding to a hydrated helical segment to assess its role in helix folding. The sequence is a site for two 'folding intermediates'. The conformational features of the model peptide have been investigated by 1H 2D NMR techniques and quantum mechanical perturbative configuration interaction over localized orbitals (PCILO) method. Theoretical modeling largely correlates with experimental observations. Based upon the amide proton temperature coefficients, the observed d alpha n(i, i + 1), d alpha n(i, i + 2), dnn(i, i + 1), d beta n(i, i + 1) NOEs and the results from theoretical modeling, we conclude that the residues of the peptide sample alpha helical and neck regions of the Ramachandran phi, psi map with reduced conformational entropy and there is a potential for turn conformations at N and C terminal ends of the peptide. The role of reduced conformational entropy and turn potential in helix formation have been discussed. We conclude that the peptide sequence can serve as a 'folding intermediate' in the helix folding of glyceraldehyde 3-phosphate dehydrogenase.     

Design,synthesis and conformational analysis of hGM-CSF[13-31]-Gly-Pro-Gly-[103-116]

On the basis of the X-ray structure and results from structure--activity relationship studies, the following GM–CSF analogue was designed and synthesized by solid-phase methodology: hGM–CSF[13-31]-Gly-Pro-Gly-[103–116]-NH₂. This analogue was constructed to comprise helices A and D of the native hGM–CSF, covalently linked in an antiparallel orientation by the tripeptide spacer Gly-Pro-Gly, which is known as a turn-inducing sequence. The conformational analysis of the analogue by CD spectroscopy revealed an essentially random structure in water, while α-helix formation was observed upon addition of TFE. In 40% TFE the helix content was ∼45%. By two-dimensional NMR experiments in 1:1 water/trifluoroethanol mixture two helical sequences were identified comprising the segments corresponding to helix A and helix D. In addition to medium-range NOESY connectivities, a long-range cross-peak was found involving the leucine residues at positions 13 and 35. Based on the experimentally derived data (54 NOEs), the structure was refined by restrained molecular dynamics simulations over 120 ps at various temperatures. A representative conformation derived from the computer simulation is mainly characterized by two helical segments connected by a loop region. The overall three-dimensional structure of the analogue is comparable to the X-ray structure of hGM–CSF in that helices A and D are oriented in an antiparallel fashion, forming a two α-helix bundle. Nevertheless, there are small differences in the topology of the helices between the solution structure of the designed analogue and the X-ray structure of hGM–CSF. The possible implications of these conformational features at the effects of biological activity are discussed. © 1997 European Peptide Society and John Wiley & Sons, Ltd. J. Pep. Sci.3: 323–335 No. of Figures: 10. No. of Tables: 5. No. of References: 46     

Different conformational families of pyrimidine.purine.pyrimidine triple helices depending on backbone composition.

Different helical conformations of DNA (D), RNA (R), and DNA.RNA (DR) hybrid double and triple helices have been detected using affinity cleavage analysis. Synthetic methods were developed to attach EDTA.Fe to a single nucleotide on RNA as well as DNA oligonucleotides. Cleavage patterns generated by a localized diffusible oxidant in the major groove on the pyrimidine strand of four purine.pyrimidine double helices consisting of all DNA, all RNA, and the corresponding hybrids reveal that the relative cleavage intensity shifts to the 5' end of the purine strand increasingly in the order: DD < DR < RD < RR. These results are consistent with models derived from structural studies. In six pyrimidine.purine.pyrimidine triple helices, the altered cleavage patterns of the Watson-Crick pyrimidine strands reveal at least two conformational families: (i) D + DD, R + DD, D + DR, and R + DR and (ii) R + RD and R + RR.     

Post-translational cleavage of recombinantly expressed nitrilase from Rhodococcus rhodochrous J1 yields a stable, active helical form

Nitrilases convert nitriles to the corresponding carboxylic acids and ammonia. The nitrilase from Rhodococcus rhodochrous J1 is known to be inactive as a dimer but to become active on oligomerization. The recombinant enzyme undergoes post-translational cleavage at approximately residue 327, resulting in the formation of active, helical homo-oligomers. Determining the 3D structure of these helices using electron microscopy, followed by fitting the stain envelope with a model based on homology with other members of the nitrilase superfamily, enables the interacting surfaces to be identified. This also suggests that the reason for formation of the helices is related to the removal of steric hindrance arising from the 39 C-terminal amino acids from the wild-type protein. The helical form can be generated by expressing only residues 1-327.     

Three-dimensional model for meta-II rhodopsin,an activated G-protein-coupled receptor

A novel approach that iteratively combined the results of energy calculations and experimental data was used to generate a three-dimensional (3D) model of the photoactivated state (R*) of bovine rhodopsin (Rh). The approach started with simplified energy calculations in an effort to find a set of sterically and energetically reasonable options for transmembrane (TM) helix arrangements with all-trans-retinal. Various 3D models of TM helix packing found by computations were then compared to limited site-directed spin-label experimental data regarding the transition of the TM helices of Rh in the inactive state (R) to those in the R* state to identify the most plausible model of the TM helical bundle. At the next step, all non-TM structural elements, such as the non-TM helix 8, the N- and C-terminal fragments, and the loops connecting TM helices, were reconstructed, and after the entire R* structure had been relaxed, all other currently available additional experimental data, both mutational and spectroscopic, on the structure of the meta-II state of rhodopsin were used to validate the resulting 3D model.