Structural context of exons in protein domains: implications for protein modelling and design

Intron boundaries were extracted from genomic data and mapped onto single-domain human and murine protein structures taken from the Protein Data Bank. A first analysis of this set of proteins shows that intron boundaries prefer to be in non-regular secondary structure elements, while avoiding alpha-helices and beta-strands. This fact alone suggests an evolutionary model in which introns are constrained by protein structure, particularly by tertiary structure contacts. In addition, in silico recombination experiments of a subset of these proteins together with their homologues, including those in different species, show that introns have a tendency to occur away from artificial crossover hot spots. Altogether, these findings support a model in which genes can preferentially harbour introns in less constrained regions of the protein fold they code for. In the light of these findings, we discuss some implications for protein modelling and design.     

Structural and functional analysis of APOA5 mutations identified in patients with severe hypertriglyceridemia

During the diagnosis of three unrelated patients with severe hypertriglyceridemia, three APOA5 mutations [p.(Ser232_Leu235)del, p.Leu253Pro, and p.Asp332ValfsX4] were found without evidence of concomitant LPL, APOC2, or GPIHBP1 mutations. The molecular mechanisms by which APOA5 mutations result in severe hypertriglyceridemia remain poorly understood, and the functional impairment/s induced by these specific mutations was not obvious. Therefore, we performed a thorough structural and functional analysis that included follow-up of patients and their closest relatives, measurement of apoA-V serum concentrations, and sequencing of the APOA5 gene in 200 nonhyperlipidemic controls. Further, we cloned, overexpressed, and purified both wild-type and mutant apoA-V variants and characterized their capacity to activate LPL. The interactions of recombinant wild-type and mutated apoA-V variants with liposomes of different composition, heparin, LRP1, sortilin, and SorLA/LR11 were also analyzed. Finally, to explore the possible structural consequences of these mutations, we developed a three-dimensional model of full-length, lipid-free human apoA-V. A complex, wide array of impairments was found in each of the three mutants, suggesting that the specific residues affected are critical structural determinants for apoA-V function in lipoprotein metabolism and, therefore, that these APOA5 mutations are a direct cause of hypertriglyceridemia.     

Homology modeling of histidine-containing phosphocarrier protein and eosinophil-derived neurotoxin: Construction of models and comparison with experiment

Homology modeling methods have been used to construct models of two proteins—the histidine-containing phosphocarrier protein (HPr) from Mycoplasma capricolum and human eosinophil-derived neurotoxin (EDN). Comparison of the models with the subsequently determined X-ray crystal structures indicates that the core regions of both proteins are reasonably well reproduced, although the template structures are closer to the X-ray structures in these regions—possible enhancements are discussed. The conformations of most of the side chains in the core of HPr are well reproduced in the modeled structure. As expected, the conformations of surface side chains in this protein differ significantly from the X-ray structure. The loop regions of EDN were incorrectly modeled—reasons for this and possible enhancements are discussed. © 1995 Wiley-Liss, Inc.     

In silico structural homology modeling of nif A protein of rhizobial strains in selective legume plants

Symbiosis is a complex genetic regulatory biological evolution which is highly specific pertaining to plant species and microbial strains. Biological nitrogen fixation in legumes is a functional combination of nodulation by nod genes and regulation by nif, fix genes. Three rhizobial strains (Rhizobium leguminosarum, Bradyrhizobium japonicum, and Mesorhizobium ciceri) that we considered for in silico analysis of nif A are proved to be the best isolates with respect to N₂ fixing for ground nut, chick pea and soya bean (in vitro) out of 47 forest soil samples. An attempt has been made to understand the structural characteristics and variations of nif genes that may reveal the factors influencing the nitrogen fixation. The primary, secondary and tertiary structure of nif A protein was analyzed by using multiple bioinformatics tools such as chou-Fasman, GOR, ExPasy ProtParam tools, Prosa -web. Literature shows that the homology modeling of nif A protein have not been explored yet which insisted the immediate development for better understanding of nif A structure and its influence on biological nitrogen fixation. In the present predicted 3D structure, the nif A protein was analyzed by three different software tools (Phyre2, Swiss model, Modeller) and validated accordingly which can be considered as an acceptable model. However further in silico studies are suggested to determine the specific factors responsible for nitrogen fixing in the present three rhizobial strains.     

Protein instability and functional defects caused by mutations of dihydro-orotate dehydrogenase in Miller syndrome patients

Miller syndrome is a recessive inherited disorder characterized by postaxial acrofacial dysostosis. It is caused by dysfunction of the DHODH (dihydroorotate dehydrogenase) gene, which encodes a key enzyme in the pyrimidine de novo biosynthesis pathway and is localized at mitochondria intermembrane space. We investigated the consequence of three missense mutations, G202A, R346W and R135C of DHODH, which were previously identified in patients with Miller syndrome. First, we established HeLa cell lines stably expressing DHODH with Miller syndrome-causative mutations: G202A, R346W and R135C. These three mutant proteins retained the proper mitochondrial localization based on immunohistochemistry and mitochondrial subfractionation studies. The G202A, R346W DHODH proteins showed reduced protein stability. On the other hand, the third one R135C, in which the mutation lies at the ubiquinone-binding site, was stable but possessed no enzymatic activity. In conclusion, the G202A and R346W mutation causes deficient protein stability, and the R135C mutation does not affect stability but impairs the substrate-induced enzymatic activity, suggesting that impairment of DHODH activity is linked to the Miller syndrome phenotype.     

Structural aspects of the functional modules in human protein kinase-Cα deduced from comparative analyses

Three-dimensional models of the five functional modules in human protein kinase Cα (PKCα) have been generated on the basis of known related structures. The catalytic region at the C-terminus of the sequence and the N-terminal auto-inhibitory pseudo-substrate have been modeled using the crystal structure complex of cAMP-dependent protein kinese (cAPK) and PKI peptide. While the N-terminal helix of the catalytic region of PKCα is predicted to be in a different location compared with cAPK, the C-terminal extension is modeled like that in the cAPK. The predicted permissive phosphorylation site of PKCα, Thr 497, is found to be entirely consistent with the mutagenesis studies. Basic Lys and Arg residues in the pseudo-substrate make several specific interactions with acidic residues in the catalytic region and may interact with the permissive phosphorylation site. Models of the two zinc-binding modules of PKCα are based on nuclear magnetic resonance and crystal structures of such modules in other PKC isoforms while the calcium phospholipid binding module (C2) is based on the crystal structure of a repeating unit in synaptotagmin I. Phorbol ester binding regions in zinc-binding modules and the calcium binding region in the C2 domain are similar to those in the basis structures. A hypothetical model of the relative positions of all five modules has the putative lipid binding ends of the C2 and the two zinc-binding domains pointing in the same direction and may serve as a basis for further experiments. © 1996 Wiley-Liss, Inc.     

Structural and functional defects caused by point mutations in the alpha-crystallin domain of a bacterial alpha-heat shock protein

The diverse family of alpha-crystallin-type small heat shock proteins (alpha-Hsps or sHsps) is characterised by a central, moderately conserved alpha-crystallin domain. Oligomerisation followed by dissociation of subparticles is thought to be a prerequisite for chaperone function. We demonstrate that HspH, a bacterial alpha-Hsp from the soybean-symbiont Bradyrhizobium japonicum, assembles into dynamic complexes freely exchanging subunits with homologous and heterologous complexes. The importance of the alpha-crystallin domain for oligomerisation and chaperone activity was tested by site-directed mutagenesis of 12 different residues. In contrast to mammalian alpha-Hsps, the majority of these mutations elicited severe structural and functional defects in HspH. The individual exchange of five amino acid residues throughout the alpha-crystallin domain was found to compromise oligomerisation to various degrees. Assembly defects resulting in complexes of reduced size correlated with greatly decreased or abolished chaperone activity, reinforcing that complete oligomerisation is required for functionality. Mutation of a highly conserved glycine (G114) at the C-terminal end of the alpha-crystallin domain specifically impaired chaperone activity without interfering with oligomerisation properties, indicating that this residue is critical for substrate interaction. The structural and functional importance of this and other residues is discussed in the context of a modeled three-dimensional structure of HspH.     

Structural and functional studies of conserved nucleotide-binding protein LptB in lipopolysaccharide transport

Lipopolysaccharide (LPS) is the main component of the outer membrane of Gram-negative bacteria, which plays an essential role in protecting the bacteria from harsh conditions and antibiotics. LPS molecules are transported from the inner membrane to the outer membrane by seven LPS transport proteins. LptB is vital in hydrolyzing ATP to provide energy for LPS transport, however this mechanism is not very clear. Here we report wild-type LptB crystal structure in complex with ATP and Mg²⁺, which reveals that its structure is conserved with other nucleotide-binding proteins (NBD). Structural, functional and electron microscopic studies demonstrated that the ATP binding residues, including K42 and T43, are crucial for LptB’s ATPase activity, LPS transport and the vitality of Escherichia coli cells with the exceptions of H195A and Q85A; the H195A mutation does not lower its ATPase activity but impairs LPS transport, and Q85A does not alter ATPase activity but causes cell death. Our data also suggest that two protomers of LptB have to work together for ATP hydrolysis and LPS transport. These results have significant impacts in understanding the LPS transport mechanism and developing new antibiotics.     

Tumor suppressor INK4: refinement of p16INK4A structure and determination of p15INK4B structure by comparative modeling and NMR data

Within the tumor suppressor protein INK4 (inhibitor of cyclin-dependent kinase 4) family, p15INK4B is the smallest and the only one whose structure has not been determined previously, probably due to the protein's conformational flexibility and instability. In this work, multidimensional NMR studies were performed on this protein. The first tertiary structure was built by comparative modeling with p16INK4A as the template, followed by restrained energy minimization with NMR constraints (NOE and H-bonds). For this purpose, the solution structure of pl6INK4A, whose quality was also limited by similar problems, was refined with additional NMR experiments conducted on an 800 MHz spectrometer and by structure-based iterative NOE assignments. The nonhelical regions showed major improvement with root-mean-square deviation (RMSD) improved from 1.23 to 0.68 A for backbone heavy atoms. The completion of p15INK4B coupled with refinement of p16INK4A made it possible to compare the structures of the four INK4 members in depth, and to compare the structures of p16INK4A in the free form and in the p16INK4A-CDK6 complex. This is an important step toward a comprehensive understanding of the precise functional roles of each INK4 member.     

Structural modeling of SoxF protein from Chlorobium tepidum: an approach to understand the molecular basis of thiosulfate oxidation

Microbial redox reactions of inorganic sulfur compounds play a vital role in balancing the turnover of this element in the environment. These vital reactions are carried out by the enzyme system encoded by the sox operon. The central player of the sulfur oxidation biochemistry is the SoxY–Z protein complex. Another protein called SoxF having sulfide dehydrogenase activity has the ability to reactivate the inactivated SoxY–Z protein complex. This SoxF protein is obtained from the sox operon of Chlorobium tepidium. In the present work an attempt has been made to understand the structural details of the activity of SoxF protein. A plausible biochemical mechanism has been predicted regarding the involvement of the SoxF protein in biological sulfur anion oxidation process. Since this is the first report regarding the structural biology of SoxF protein this study may shed light in the hitherto unknown molecular biochemistry of sulfur anion oxidation by sox operon.     

Adrenoleukodystrophy: subcellular localization and degradation of adrenoleukodystrophy protein (ALDP/ABCD1) with naturally occurring missense mutations

Mutation in the X-chromosomal adrenoleukodystrophy gene (ALD; ABCD1) leads to X-linked adrenoleukodystrophy (X-ALD), a severe neurodegenerative disorder. The encoded adrenoleukodystrophy protein (ALDP/ABCD1) is a half-size peroxisomal ATP-binding cassette protein of 745 amino acids in humans. In this study, we chose nine arbitrary mutant human ALDP forms (R104C, G116R, Y174C, S342P, Q544R, S606P, S606L, R617H, and H667D) with naturally occurring missense mutations and examined the intracellular behavior. When expressed in X-ALD fibroblasts lacking ALDP, the expression level of mutant His-ALDPs (S606L, R617H, and H667D) was lower than that of wild type and other mutant ALDPs. Furthermore, mutant ALDP-green fluorescence proteins (S606L and H667D) stably expressed in CHO cells were not detected due to rapid degradation. Interestingly, the wild type ALDP co-expressed in these cells also disappeared. In the case of X-ALD fibroblasts from an ALD patient (R617H), the mutant ALDP was not detected in the cells, but appeared upon incubation with a proteasome inhibitor. When CHO cells expressing mutant ALDP-green fluorescence protein (H667D) were cultured in the presence of a proteasome inhibitor, both the mutant and wild type ALDP reappeared. In addition, mutant His-ALDP (Y174C), which has a mutation between transmembrane domain 2 and 3, did not exhibit peroxisomal localization by immunofluorescense study. These results suggest that mutant ALDPs, which have a mutation in the COOH-terminal half of ALDP, including S606L, R617H, and H667D, were degraded by proteasomes after dimerization. Further, the region between transmembrane domain 2 and 3 is important for the targeting of ALDP to the peroxisome.     

Structural basis for type I and type II deficiencies of antithrombotic plasma protein C: Patterns revealed by three-dimensional molecular modelling of mutations of the protease domain

Familial deficiency of protein C is associated with inherited thrombophilia. To explore how specific missense mutations might cause observed clinical phenotypes, known protein C missense mutations were mapped onto three-dimensional homology models of the protein C protease domain, and the implications for domain folding and structure were evaluated. Most Type I missense mutations either replaced internal hydrophobic residues (I201T, L223F, A259V, A267T, A346T, A346V, G376D) or nearby interacting residues (I403M, T298M, Q184H), thus disrupting the packing of internal hydrophobic side chains, or changed hydrophilic residues, thus disrupting ion pairs (N256D, R178W). Mutations (P168L, R169W) at the activation site destabilized the region containing the activation peptide structure. Most Type II mutations involved solvent-exposed residues and were clustered either in a positively charged region (R147W, R157Q, R229Q, R352W) or were located in or near the active site region (S252N, D359N, G381S, G391S, H211Q). The cluster of arginines 147, 157, 229, and 352 may identify a functionally important exosite. Identification of the spatial relationships of natural mutations in the protein C model is helpful for understanding manifestations of protein C deficiency and for identification of novel, functionally important molecular features and exosites. © 1994 John Wiley & Sons, Inc.     

Molecular modeling of D151Y and M391T mutations in the LDL receptor

The low-density lipoprotein receptor (LDLR) is a key regulator of cholesterol homeostasis, and defects in the function of LDLR result in familial hypercholesterolemia (FH). In the present study, we performed structural analyses of two novel LDLR mutations, D151Y and M391T. Both mutations occurred in conserved residues of LDLR. The D151Y mutation, in the ligand binding domain, caused an elimination of a hydrogen bond in the calcium binding site, higher solvent accessibility and a loss of negative charge in the Y151 residue. On the other hand, the M391T mutation, in the β-propeller of the epidermal growth factor (EGF) precursor homology domain, caused an additional hydrogen bond to form, higher solvent accessibility and a distortion of the β-strand. These data suggest that the irregular structures of the mutated LDLRs are likely to cause the functional defect that contributes to the pathology of FH.     

Glycosylation of the OCTN2 carnitine transporter: Study of natural mutations identified in patients with primary carnitine deficiency

Primary carnitine deficiency is caused by impaired activity of the Na⁺-dependent OCTN2 carnitine/organic cation transporter. Carnitine is essential for entry of long-chain fatty acids into mitochondria and its deficiency impairs fatty acid oxidation. Most missense mutations identified in patients with primary carnitine deficiency affect putative transmembrane or intracellular domains of the transporter. Exceptions are the substitutions P46S and R83L located in an extracellular loop close to putative glycosylation sites (N57, N64, and N91) of OCTN2. P46S and R83L impaired glycosylation and maturation of OCTN2 transporters to the plasma membrane. We tested whether glycosylation was essential for the maturation of OCTN2 transporters to the plasma membrane. Substitution of each of the three asparagine (N) glycosylation sites with glutamine (Q) decreased carnitine transport. Substitution of two sites at a time caused a further decline in carnitine transport that was fully abolished when all three glycosylation sites were substituted by glutamine (N57Q/N64Q/N91Q). Kinetic analysis of carnitine and sodium-stimulated carnitine transport indicated that all substitutions decreased the Vmax for carnitine transport, but N64Q/N91Q also significantly increased the Km toward carnitine, indicating that these two substitutions affected regions of the transporter important for substrate recognition. Western blot analysis confirmed increased mobility of OCTN2 transporters with progressive substitutions of asparagines 57, 64 and/or 91 with glutamine. Confocal microscopy indicated that glutamine substitutions caused progressive retention of OCTN2 transporters in the cytoplasm, up to full retention (such as that observed with R83L) when all three glycosylation sites were substituted. Tunicamycin prevented OCTN2 glycosylation, but it did not impair maturation to the plasma membrane. These results indicate that OCTN2 is physiologically glycosylated and that the P46S and R83L substitutions impair this process. Glycosylation does not affect maturation of OCTN2 transporters to the plasma membrane, but the 3 asparagines that are normally glycosylated are located in a region important for substrate recognition and turnover rate.     

Modeling of peptides and proteins in a membrane environment: II. Structural and energetic aspects of glycophorin A in a lipid bilayer

The conformational space of a hydrophobic peptide fragment of glycophorin A in a lipid membrane was studied with the Monte Carlo method using the solvation model described in the first communication of this series. The simulation was performed for various starting orientations of the peptide relative the membrane bilayer: outside, inside, partially immersed, and transbilayer. We showed that the membrane substantially stabilizes the α-helical conformation of the central hydrophobic part of the glycophorin A molecule, which for the most part is immersed in the apolar core of the bilayer. For various conformational states, energy values were calculated and the orientations of the peptide relative to the membrane were characterized. Depending on the thickness of the bilayer, either an entirely α-helical conformation in transbilayer orientation or a conformation with a kink in the central part of the helix with theN- andC-termini exposed on one side of the membrane corresponds to the minimal-energy structure. The transmembrane orientation of glycophorin A is energetically advantageous when the membrane thickness is close to the length of its hydrophobic helical portion, which is consistent with the effect ofhydrophobic match observed experimentally. The prospects for further refinement of the model are discussed. For communication I, see [1].     

On the structure of hisH: protein structure prediction in the context of structural and functional genomics

We predict a structure of the glutamine amidotransferase subunit (hisH) of imidazole glycerol phosphate synthase (IGPS) which catalyzes the fifth step of the histidine biosynthesis in Escherichia coli. The model is constructed using an energy-based threading program augmented by a multiple sequence to structure profile analysis. In developing our model we identified a conserved core region within hisH and a variable domain which is the likely site of interaction with the synthase subunit (hisF) of IGPS. Information available from structural and functional genomics studies was used to improve the structure prediction, to discuss parallels between histidine biosynthesis and other amino acid and nucleotide metabolic pathways, and to better understand the protein-protein interactions between the hisH and hisF domains of IGPS. This work allows us to develop a preliminary model for the structure of the entire IGPS holoenzyme.