Fold recognition and ab initio structure predictions using hidden markov models and β-strand pair potentials

Protein structure predictions were submitted for 9 of the target sequences in the competition that ran during 1994. Targets sequences were selected that had no known homology with any sequence of known structure and were members of a reasonably sized family of related but divergent sequences. The objective was either to recognize a compatible fold for the target sequence in the database of known structures or to predict ab initio its rough 3D topology. The main tools used were Hidden Markov models (HMM) for fold recognition, a β- strand pair potential to predict β-sheet topology, and the PHD server for secondary structure prediction. Compatible folds were correctly identified in a number of cases and the β-strand pair potential was shown to be a useful tool for ab initio topology prediction. © 1995 Wiley-Liss, Inc.     

Effects of site-directed mutation on the function of the chloroplast ATP synthase ε subunit

The previous work in our lab showed that the spinach chloroplast ATP synthase ε mutant with 3 amino acid residues deleted from the N-terminus had much lower ability to inhibit ATP hydrolysis and block proton leakage in comparison to a mutant with 1 or 2 residues deleted from the N-terminus. The present study aimed at determining whether there is special importance in the structure and function of the N-terminal third residue of the chloroplast ε subunit. The leucine residue at the N-terminal third site (Leu3) of the spinach chloroplast ε subunit was replaced with Ile, Phe, Thr, Arg, Glu or Pro by site-directed mutagenesis, forming mutants εL3I, εL3F, εL3T, εL3R, εL3E and εL3P, respectively. These ε variants all showed lower abilities to inhibit ATP hydrolysis and to block proton leakage, as compared to the wild type ε subunit (εWT). The abilities of mutants εL3I and εL3F to restore the ATP synthesis activity of reconstituted membranes were higher than those of εWT, but the abilities of the other ε variants were lower than that of εWT. These results indicate that the hydrophobic and neutral characteristics of Leu3 of the chloroplast ε subunit are very important for its ability to inhibit ATP hydrolysis and block proton leakage, and for the ATP synthesis ability of ATP synthase.     

Review: Model Peptides and the Physicochemical Approach to β-Amyloids

β-Amyloid peptides are the main protein components of neuritic plaques and may be important in the pathogenesis of Alzheimer's Disease. The determination of the structure of β-amyloid fibrils poses a challenge because of the limited solubility of β-amyloid peptides and the noncrystalline nature of fibrils formed from these peptides. In this paper, we describe several physicochemical approaches which have been used to examine fibrils and the fibrillogenesis of peptide models of β-amyloid. Recent advances in solid state NMR, such as the DRAWS pulse sequence, have made this approach a particularly attractive one for peptides such as β-amyloid, which are not yet amenable to high-resolution solution phase NMR and crystallography. The application of solid state NMR techniques has yielded information on a model peptide comprising residues 10–35 of human β-amyloid and indicates that in fibrils, this peptide assumes a parallel β-strand conformation, with all residues in exact register. In addition, we discuss the use of block copolymers of Aβ peptides and polyethylene glycol as probes for the pathways of fibrillogenesis. These methods can be combined with other new methods, such as high-resolution synchrotron X-ray diffraction and small angle neutron and X-ray scattering, to yield structural data of relevance not only to disease, but to the broader question of protein folding and self-assembly.     

Characterization of Tiki,a New Family of Wnt-specific Metalloproteases

The Wnt family of secreted glycolipoproteins plays pivotal roles in development and human diseases. Tiki family proteins were identified as novel Wnt inhibitors that act by cleaving the Wnt amino-terminal region to inactivate specific Wnt ligands. Tiki represents a new metalloprotease family that is dependent on Mn²⁺/Co²⁺ but lacks known metalloprotease motifs. The Tiki extracellular domain shares homology with bacterial TraB/PrgY proteins, known for their roles in the inhibition of mating pheromones. The TIKI/TraB fold is predicted to be distantly related to structures of additional bacterial proteins and may use a core β-sheet within an α+β-fold to coordinate conserved residues for catalysis. In this study, using assays for Wnt3a cleavage and signaling inhibition, we performed mutagenesis analyses of human TIKI2 to examine the structural prediction and identify the active site residues. We also established an in vitro assay for TIKI2 protease activity using FRET peptide substrates derived from the cleavage motifs of Wnt3a and Xenopus wnt8 (Xwnt8). We further identified two pairs of potential disulfide bonds that reside outside the β-sheet catalytic core but likely assist the folding of the TIKI domain. Finally, we systematically analyzed TIKI2 cleavage of the 19 human WNT proteins, of which we identified 10 as potential TIKI2 substrates, revealing the hydrophobic nature of Tiki cleavage sites. Our study provides insights into the Tiki family of proteases and its Wnt substrates.     

Key Residues in the Nicotinic Acetylcholine Receptor β2 Subunit Contribute to α-Conotoxin LvIA Binding

α-Conotoxin LvIA (α-CTx LvIA) is a small peptide from the venom of the carnivorous marine gastropod Conus lividus and is the most selective inhibitor of α3β2 nicotinic acetylcholine receptors (nAChRs) known to date. It can distinguish the α3β2 nAChR subtype from the α6β2* (* indicates the other subunit) and α3β4 nAChR subtypes. In this study, we performed mutational studies to assess the influence of residues of the β2 subunit versus those of the β4 subunit on the binding of α-CTx LvIA. Although two β2 mutations, α3β2[F119Q] and α3β2[T59K], strongly enhanced the affinity of LvIA, the β2 mutation α3β2[V111I] substantially reduced the binding of LvIA. Increased activity of LvIA was also observed when the β2-T59L mutant was combined with the α3 subunit. There were no significant difference in inhibition of α3β2[T59I], α3β2[Q34A], and α3β2[K79A] nAChRs when compared with wild-type α3β2 nAChR. α-CTx LvIA displayed slower off-rate kinetics at α3β2[F119Q] and α3β2[T59K] than at the wild-type receptor, with the latter mutant having the most pronounced effect. Taken together, these data provide evidence that the β2 subunit contributes to α-CTx LvIA binding and selectivity. The results demonstrate that Val¹¹¹ is critical and facilitates LvIA binding; this position has not previously been identified as important to binding of other 4/7 framework α-conotoxins. Thr⁵⁹ and Phe¹¹⁹ of the β2 subunit appear to interfere with LvIA binding, and their replacement by the corresponding residues of the β4 subunit leads to increased affinity.     

The solution structure of the amino-terminal domain of human DNA polymerase epsilon subunit B is homologous to C-domains of AAA+ proteins

DNA polymerases α, δ and ε are large multisubunit complexes that replicate the bulk of the DNA in the eukaryotic cell. In addition to the homologous catalytic subunits, these enzymes possess structurally related B subunits, characterized by a carboxyterminal calcineurin-like and an aminoproximal oligonucleotide/oligosaccharide binding-fold domain. The B subunits also share homology with the exonuclease subunit of archaeal DNA polymerases D. Here, we describe a novel domain specific to the N-terminus of the B subunit of eukaryotic DNA polymerases ε. The N-terminal domain of human DNA polymerases ε (Dpoe2NT) expressed in Escherichia coli was characterized. Circular dichroism studies demonstrated that Dpoe2NT forms a stable, predominantly α-helical structure. The solution structure of Dpoe2NT revealed a domain that consists of a left-handed superhelical bundle. Four helices are arranged in two hairpins and the connecting loops contain short β-strand segments that form a short parallel sheet. DALI searches demonstrated a striking structural similarity of the Dpoe2NT with the α-helical subdomains of ATPase associated with various cellular activity (AAA+) proteins (the C-domain). Like C-domains, Dpoe2NT is rich in charged amino acids. The biased distribution of the charged residues is reflected by a polarization and a considerable dipole moment across the Dpoe2NT. Dpoe2NT represents the first C-domain fold not associated with an AAA+ protein.     

Structure of the Lectin Mannose 6-Phosphate Receptor Homology (MRH) Domain of Glucosidase II,an Enzyme That Regulates Glycoprotein Folding Quality Control in the Endoplasmic Reticulum

Here we report for the first time the three-dimensional structure of a mannose 6-phosphate receptor homology (MRH) domain present in a protein with enzymatic activity, glucosidase II (GII). GII is involved in glycoprotein folding in the endoplasmic reticulum. GII removes the two innermost glucose residues from the Glc₃Man₉GlcNAc₂ transferred to nascent proteins and the glucose added by UDP-Glc:glycoprotein glucosyltransferase. GII is composed of a catalytic GIIα subunit and a regulatory GIIβ subunit. GIIβ participates in the endoplasmic reticulum localization of GIIα and mediates in vivo enhancement of N-glycan trimming by GII through its C-terminal MRH domain. We determined the structure of a functional GIIβ MRH domain by NMR spectroscopy. It adopts a β-barrel fold similar to that of other MRH domains, but its binding pocket is the most shallow known to date as it accommodates a single mannose residue. In addition, we identified a conserved residue outside the binding pocket (Trp-409) present in GIIβ but not in other MRHs that influences GII glucose trimming activity.     

Prediction of the structure of the replication initiator protein DnaA

The secondary structure of DnaA protein and its interaction with DNA and ribonucleotides has been predicted using biochemical, biophysical techniques, and prediction methods based on multiple-sequence alignment and neural networks. The core of all proteins from the DnaA family consists of an “open twisted α/β structure,” containing five α-helices alternating with five β-strands. In our proposed structural model the interior of the core is formed by a parallel β-sheet, whereas the α-helices are arranged on the surface of the core. The ATP-binding motif is located within the core, in a loop region following the first β-strand. The N-terminal domain (80 aa) is composed of two α-helices, the first of which contains a potential leucine zipper motif for mediating protein-protein interaction, followed by a β-strand and an additional α-helix. The N-terminal domain and the α/β core region of DnaA are connected by a variable loop (45–70 aa); major parts of the loop region can be deleted without loss of protein activity. The C-terminal DNA-binding domain (94 aa) is mostly α-helical and contains a potential helix-loop-helix motif. DnaA protein does not dimerize in solution; instead, the two longest C-terminal α-helices could interact with each other, forming an internal “coiled coil” and exposing highly basic residues of a small loop region on the surface, probably responsible for DNA backbone contacts. © 1997 Wiley-Liss Inc.     

The Escherichia coli outer membrane protein OmpA acquires secondary structure prior to its integration into the membrane

Almost all proteins that reside in the outer membrane (OM) of Gram-negative bacteria contain a membrane-spanning segment that folds into a unique β barrel structure and inserts into the membrane by an unknown mechanism. To obtain further insight into outer membrane protein (OMP) biogenesis, we revisited the surprising observation reported over 20 years ago that the Escherichia coli OmpA β barrel can be assembled into a native structure in vivo when it is expressed as two noncovalently linked fragments. Here, we show that disulfide bonds between β strand 4 in the N-terminal fragment and β strand 5 in the C-terminal fragment can form in the periplasmic space and greatly increase the efficiency of assembly of “split” OmpA, but only if the cysteine residues are engineered in perfect register (i.e., they are aligned in the fully folded β barrel). In contrast, we observed only weak disulfide bonding between β strand 1 in the N-terminal fragment and β strand 8 in the C-terminal fragment that would form a closed or circularly permutated β barrel. Our results not only demonstrate that β barrels begin to fold into a β-sheet-like structure before they are integrated into the OM but also help to discriminate among the different models of OMP biogenesis that have been proposed.     

A Model of the Reactive Form of Plasminogen Activator Inhibitor-1

A model of the reactive form of plasminogen activator inhibitor-1 (PAI-1) has been constructed using molecular graphics and starting from the known crystal structure of latent PAI-1. The residues P16 to P10′, of which P16-P4 form strand 4 of the β-sheet A (s4A) and P3-P10′ form an extended loop in the latent form, have been removed and remodeled into this structure, based on the structures of ovalbumin and cleaved α1-proteinase inhibitor. Residues P4′-P10′ were remodeled as a β-strand s1C, located on the surface of the molecule and the N-terminal end (P16-P14) of the eliminated loop was rebuilt using appropriate backbone dihedrals. Subsequently, a secondary structure prediction program was applied and further optimization of the model was performed by several molecular dynamics runs. Apparently the β-strand was stabilized by only two hydrogen bonds. Further analysis revealed that, although s4A was removed, s3A and s5A did not approach each other. In this current model it was also found that the large gap between the loop connecting s4C-s3C and the loop connecting s3B-hG remained 11 Å in contrast to the small gap (4Å) at a similar position in other serpins. These observations may explain the ease of a conformational change of the reactive site loop of PAI-1 during transition to the latent and the preinserted form. In addition the current model can be used for the design of stable, functional, PAI-1 mutants. Detailed structural analysis of the latter may facilitate studies on the structure-function relationship in PAI-1 in particular and in other serpins in general.     

A consensus prediction of the secondary structure for the 6-phospho-β-D-galactosidase superfamily

Two separate unrefined models for the secondary structure of two subfamilies of the 6-phospho-β-D -galactosidase superfamily were independently constructed by examining patterns of variation and conservation within homologous protein sequences, assigning surface, interior, parsing, and active site residues to positions in the alignment, and identifying periodicities in these. A consensus model for the secondary structure of the entire superfamily was then built. The prediction tests the limits of an unrefined prediction made using this approach in a large protein with substantial functional and sequence divergence within the family. The protein belongs to the (α–β class), with the core β strands aligned parallel. The supersecondary structural elements that are readily identified in this model is a parallel β sheet built by strands C, D, and E, with helices 2 and 3 connecting strands (C + D) and (D + E), respectively, and an analogous α–β unit (strand G and helix 7) toward the end of the sequence. The resemblance of the supersecondary model to the tertiary structure formed by 8-fold α–β barrel proteins is almost certainly not coincidental. © 1995 Wiley-Liss, Inc.     

Subunit exchange of lens α-crystallin: a fluorescence energy transfer study with the fluorescent labeled αA-crystallin mutant W9F as a probe

A Trp-free αA-crystallin mutant (W9F) was prepared by site-directed mutation. This mutant appears to be identical to the wild-type in terms of conformation (secondary and tertiary structures). W9F was labeled with a sulfhydryl-specific fluorescent probe, 2-(4′-maleimidylanilino) naphthalene-6-sulfonate (MIANS), and used in a subunit exchange between αA- and αA-crystallins as well as between αA- and αB-crystallins, studied by measurement of fluorescence resonance energy transfer. Energy transfer was observed between Trp (donor, with emission maximum at 336 nm) of wild-type αA- or αB-crystallin and MIANS (acceptor, with absorption maximum at 313 nm) of labeled W9F when subunit exchange occurred. Time-dependent decrease of Trp and increase of MIANS fluorescence were recorded. The exchange was faster at 37°C than at 25°C. The energy transfer efficiency was greater between homogeneous subunits (αA-αA) than between heterogeneous subunits (αA-αB). A previous exchange study with isoelectric focusing indicated a complete but slow exchange between αA and αB subunits. The present study showed that the exchange was a fast process, and the different energy transfer efficiencies between αA-αA and αA-αB indicated that αA- and αB-crystallins were not necessarily structurally equivalent.     

Three-dimensional structure of ribulose-1,5-bisphosphate carboxylase/oxygenase from Rhodospirillum rubrum at 2.9 A resolution

The three-dimensional structure of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) from Rhodospirillum rubrum has been determined at 2.9 Å resolution by X-ray crystallographic methods. The MIR-electron density map was substantially improved by two-fold non-crystallographic symmetry averaging. The polypeptide chains in the dimer were traced using a graphics display system with the help of the BONES option in FRODO. The dimer has approximate dimensions of 50 x 72 x 105 Å. The enzyme subunit is a typical two-domain protein. The smaller, N-terminal domain consists of 137 amino acid residues and forms a central, mixed five-stranded β-sheet with α-helices on both sides of the sheet. The larger C-terminal domain consists of 329 amino acid residues. This domain has an eight-stranded parallel α/β barrel structure as found in triosephosphate isomerase and a number of other functionally non-related proteins. The active site in Rubisco determined by difference Fourier techniques and fitting of active site residues to the electron density map, is located at the carboxy-end of the β-strands in the α/β barrel of the C-terminal domain. There are few domain–domain interactions within the subunit. The interactions at the interface between the two subunits of the dimer are tight and extensive. There are tight contacts between the two C-terminal domains, which build up the core of the molecule. There are also interactions between the N-terminal domain of one subunit and the C-terminal domain of the second subunit, close to the active site.     

Coevolving residues of (beta/alpha)(8)-barrel proteins play roles in stabilizing active site architecture and coordinating protein dynamics

Indole-3-glycerol phosphate synthase (IGPS) is a representative of (β/α)₈-barrel proteins—the most common enzyme fold in nature. To better understand how the constituent amino-acids work together to define the structure and to facilitate the function, we investigated the evolutionary and dynamical coupling of IGPS residues by combining statistical coupling analysis (SCA) and molecular dynamics (MD) simulations. The coevolving residues identified by the SCA were found to form a network which encloses the active site completely. The MD simulations showed that these coevolving residues are involved in the correlated and anti-correlated motions. The correlated residues are within van der Waals contact and appear to maintain the active site architecture; the anti-correlated residues are mainly distributed on opposite sides of the catalytic cavity and coordinate the motions likely required for the substrate entry and product release. Our findings might have broad implications for proteins with the highly conserved (βα)₈-barrel in assessing the roles of amino-acids that are moderately conserved and not directly involved in the active site of the (β/α)₈-barrel. The results of this study could also provide useful information for further exploring the specific residue motions for the catalysis and protein design based on the (β/α)₈-barrel scaffold.     

Three-dimensional structure of rat acid phosphatase.

The crystal structure of recombinant rat prostatic acid phosphatase was determined to 3 A resolution with protein crystallographic methods. The enzyme subunit is built up of two domains, an alpha/beta domain consisting of a seven-stranded mixed beta-sheet with helices on both sides of the sheet and a smaller alpha domain. Two disulfide bridges between residues 129-340 and 315-319 were found. Electron density at two of the glycosylation sites for parts of the carbohydrate moieties was observed. The dimer of acid phosphatase is formed through two-fold interactions of edge strand 3 from one subunit with strand 3 from the second subunit, thus extending the beta-sheet from seven to 14 strands. Other subunit-subunit interactions involve conserved residues from loops between helices and beta-strands. The fold of the alpha/beta domain is similar to the fold observed in phosphoglycerate mutase. The active site is at the carboxy end of the parallel strands of the alpha/beta domain. There is a strong residual electron density at the phosphate binding site which probably represents a bound chloride ion. Biochemical properties and results from site-directed mutagenesis experiments of acid phosphatase are correlated to the three-dimensional structure.     

Genealogy of the α-crystallin—small heat-shock protein superfamily

Sequences of 40 very diverse representatives of the α-crystallin–small heat-shock protein (α-Hsp) superfamily are compared. Their characteristic C-terminal ‘α-crystallin domain' of 80–100 residues contains short consensus sequences that are highly conserved from prokaryotes to eukaryotes. There are, in addition, some positions that clearly distinguish animal from non-animal α-Hsps. The α-crystallin domain is predicted to consist of two hydrophobic β-sheet motifs, separated by a hydrophilic region which is variable in length. Combination of a conserved α-crystallin domain with a variable N-terminal domain and C-terminal extension probably modulates the properties of the various α-Hsps as stress-protective and structural oligomeric proteins. Phylogeny reconstruction indicates that multiple α-Hsps were already present in the last common ancestor of pro- and eukaryotes. It is suggested that during eukaryote evolution, animal and non-animal α-Hsps originated from different ancestral gene copies. Repeated gene duplications gave rise to the multiple α-Hsps present in most organisms.     

Almost half of the RTX domain is dispensable for complement receptor 3 binding and cell-invasive activity of the Bordetella adenylate cyclase toxin

The whooping cough agent Bordetella pertussis secretes an adenylate cyclase toxin (CyaA) that through its large carboxy-proximal Repeat-in-ToXin (RTX) domain binds the complement receptor 3 (CR3). The RTX domain consists of five blocks (I–V) of characteristic glycine and aspartate-rich nonapeptides that fold into five Ca²⁺-loaded parallel β-rolls. Previous work indicated that the CR3-binding structure comprises the interface of β-rolls II and III. To test if further portions of the RTX domain contribute to CR3 binding, we generated a construct with the RTX block II/III interface (CyaA residues 1132–1294) linked directly to the C-terminal block V fragment bearing the folding scaffold (CyaA residues 1562–1681). Despite deletion of 267 internal residues of the RTX domain, the Ca²⁺-driven folding of the hybrid block III/V β-roll still supported formation of the CR3-binding structure at the interface of β-rolls II and III. Moreover, upon stabilization by N- and C-terminal flanking segments, the block III/V hybrid-comprising constructs competed with CyaA for CR3 binding and induced formation of CyaA toxin-neutralizing antibodies in mice. Finally, a truncated CyaAΔ_1295-1561 toxin bound and penetrated erythrocytes and CR3-expressing cells, showing that the deleted portions of RTX blocks III, IV, and V (residues 1295–1561) were dispensable for CR3 binding and for toxin translocation across the target cell membrane. This suggests that almost a half of the RTX domain of CyaA is not involved in target cell interaction and rather serves the purpose of toxin secretion.     

Solution structure of inhibitor-free human metalloelastase (MMP-12) indicates an internal conformational adjustment

Macrophage metalloelastase or matrix metalloproteinase-12 (MMP-12) appears to exacerbate atherosclerosis, emphysema, aortic aneurysm, rheumatoid arthritis, and inflammatory bowel disease. An inactivating E219A mutation, validated by crystallography and NMR spectra, prevents autolysis of MMP-12 and allows us to determine its NMR structure without an inhibitor. The structural ensemble of the catalytic domain without an inhibitor is based on 2813 nuclear Overhauser effects (NOEs) and has an average RMSD to the mean structure of 0.25 Å for the backbone and 0.61 Å for all heavy atoms for residues Trp109-Gly263. Compared to crystal structures of MMP-12, helix B (hB) at the active site is unexpectedly more deeply recessed under the β-sheet. This opens a pocket between hB and β-strand IV in the active-site cleft. Both hB and an internal cavity are shifted toward β-strand I, β-strand III, and helix A on the back side of the protease. About 25 internal NOE contacts distinguish the inhibitor-free solution structure and indicate hB's greater depth and proximity to the sheet and helix A. Line broadening and multiplicity of amide proton NMR peaks from hB are consistent with hB undergoing a slow conformational exchange among subtly different environments. Inhibitor-binding-induced perturbations of the NMR spectra of MMP-1 and MMP-3 map to similar locations across MMP-12 and encompass the internal conformational adjustments. Evolutionary trace analysis suggests a functionally important network of residues that encompasses most of the locations adjusting in conformation, including 18 residues with NOE contacts unique to inhibitor-free MMP-12. The conformational change, sequence analysis, and inhibitor perturbations of NMR spectra agree on the network they identify between structural scaffold and the active site of MMPs.     

Subcellular dynamics and functional activity of the cleaved intracellular domain of the Na+ channel β1 subunit

The voltage-gated Na⁺ channel β1 subunit, encoded by SCN1B, regulates cell surface expression and gating of α subunits and participates in cell adhesion. β1 is cleaved by α/β and γ-secretases, releasing an extracellular domain and intracellular domain (ICD), respectively. Abnormal SCN1B expression/function is linked to pathologies including epilepsy, cardiac arrhythmia, and cancer. In this study, we sought to determine the effect of secretase cleavage on β1 function in breast cancer cells. Using a series of GFP-tagged β1 constructs, we show that β1-GFP is mainly retained intracellularly, particularly in the endoplasmic reticulum and endolysosomal pathway, and accumulates in the nucleus. Reduction in endosomal β1-GFP levels occurred following γ-secretase inhibition, implicating endosomes and/or the preceding plasma membrane as important sites for secretase processing. Using live-cell imaging, we also report β1ICD-GFP accumulation in the nucleus. Furthermore, β1-GFP and β1ICD-GFP both increased Na⁺ current, whereas β1STOP-GFP, which lacks the ICD, did not, thus highlighting that the β1-ICD is necessary and sufficient to increase Na⁺ current measured at the plasma membrane. Importantly, although the endogenous Na⁺ current expressed in MDA-MB-231 cells is tetrodotoxin (TTX)-resistant (carried by Na_v1.5), the Na⁺ current increased by β1-GFP or β1ICD-GFP was TTX-sensitive. Finally, we found β1-GFP increased mRNA levels of the TTX-sensitive α subunits SCN1A/Na_v1.1 and SCN9A/Na_v1.7. Taken together, this work suggests that the β1-ICD is a critical regulator of α subunit function in cancer cells. Our data further highlight that γ-secretase may play a key role in regulating β1 function in breast cancer.