Solution structure of a lipid transfer protein extracted from rice seeds. Comparison with homologous proteins.

Nuclear magnetic resonance (NMR) spectroscopy was used to determine the three dimensional structure of rice nonspecific lipid transfer protein (ns-LTP), a 91 amino acid residue protein belonging to the broad family of plant ns-LTP. Sequence specific assignment was obtained for all but three HN backbone 1H resonances and for more than 95% of the 1H side-chain resonances using a combination of 1H 2D NOESY; TOCSY and COSY experiments at 293 K. The structure was calculated on the basis of four disulfide bridge restraints, 1259 distance constraints derived from 1H-1H Overhauser effects, 72 phi angle restraints and 32 hydrogen-bond restraints. The final solution structure involves four helices (H1: Cys3-Arg18, H2: Ala25-Ala37, H3: Thr41-Ala54 and H4: Ala66-Cys73) followed by a long C-terminal tail (T) with no observable regular structure. N-capping residues (Thr2, Ser24, Thr40), whose side-chain oxygen atoms are involved in hydrogen bonds with i + 3 amide proton additionally stabilize the N termini of the first three helices. The fourth helix involving Pro residues display a mixture of alpha and 3(10) conformation. The rms deviation of 14 final structures with respect to the average structure is 1.14 +/- 0.16 A for all heavy atoms (C, N, O and S) and 0.72 +/- 0.01 A for the backbone atoms. The global fold of rice ns-LTP is close to the previously published structures of wheat, barley and maize ns-LTPs exhibiting nearly identical pattern of the numerous sequence specific interactions. As reported previously for different four-helix topology proteins, hydrophobic, hydrogen bonding and electrostatic mechanisms of fold stabilization were found for the rice ns-LTP. The sequential alignment of 36 ns-LTP primary structures strongly suggests that there is a uniform pattern of specific long-range interactions (in terms of sequence), which stabilize the fold of all plant ns-LTPs.     

Three-dimensional structure of acyl carrier protein in solution determined by nuclear magnetic resonance and the combined use of dynamical simulated annealing and distance geometry

The solution conformation of acyl carrier protein from Escherichia coli (77 residues) has been determined on the basis of 423 interproton-distance restraints and 32 hydrogen-bonding restraints derived from NMR measurements. A total of nine structures were computed using a hybrid approach combining metric matrix distance geometry and dynamic simulated annealing. The polypeptide fold is well defined with an average backbone atomic root-mean-square difference of 0.20 +/- 0.03 nm between the final nine converged structures and the mean structure obtained by averaging their coordinates. The principal structural motif is composed of three helices: 1 (residues 3-12), 2 (residues 37-47) and 4 (residues 65-75) which line a hydrophobic cavity. Helices 2 and 4 are approximately parallel to each other and anti-parallel at an angle of approximately equal to 150 degrees to helix 1. The smaller helix 3 (residues 56-63) is at an angle of approximately equal to 100 degrees to helix 4.     

The death-domain fold of the ASC PYRIN domain, presenting a basis for PYRIN/PYRIN recognition

The PYRIN domain is a conserved sequence motif identified in more than 20 human proteins with putative functions in apoptotic and inflammatory signalling pathways. The three-dimensional structure of the PYRIN domain from human ASC was determined by NMR spectroscopy. The structure determination reveals close structural similarity to death domains, death effector domains, and caspase activation and recruitment domains, although the structural alignment with these other members of the death-domain superfamily differs from previously predicted amino acid sequence alignments. Two highly positively and negatively charged surfaces in the PYRIN domain of ASC result in a strong electrostatic dipole moment that is predicted to be present also in related PYRIN domains. These results suggest that electrostatic interactions play an important role for the binding between PYRIN domains. Consequently, the previously reported binding between the PYRIN domains of ASC and ASC2/POP1 or between the zebrafish PYRIN domains of zAsc and Caspy is proposed to involve interactions between helices 2 and 3 of one PYRIN domain with helices 1 and 4 of the other PYRIN domain, in analogy to previously reported homophilic interactions between caspase activation and recruitment domains.     

Human interleukin 4. The solution structure of a four-helix bundle protein.

Heteronuclear 13C and 15N three-dimensional nuclear magnetic resonance (n.m.r.) techniques have been used to determine the solution structure of human interleukin 4, a four-helix bundle protein. A dynamical simulated annealing protocol was used to calculate an ensemble of structures from an n.m.r. data set of 1735 distance restraints, 101 phi angle restraints and 27 pairs of hydrogen bond restraints. The protein structure has a left-handed up-up-down-down topology for the four helices with the two long overhand loops in the structure being connected by a short section of irregular antiparallel beta-sheet. Analysis of the side-chains in the protein shows a clustering of hydrophobic residues, particularly leucines, in the core of the bundle with the side-chains of charged residues being located on the protein surface. The solution structure has been compared with a recent structure prediction for human interleukin 4 and with crystal structures of other helix bundle proteins.     

Automated method for modeling seven-helix transmembrane receptors from experimental data.

A rule-based automated method is presented for modeling the structures of the seven transmembrane helices of G-protein-coupled receptors. The structures are generated by using a simulated annealing Monte Carlo procedure that positions and orients rigid helices to satisfy structural restraints. The restraints are derived from analysis of experimental information from biophysical studies on native and mutant proteins, from analysis of the sequences of related proteins, and from theoretical considerations of protein structure. Calculations are presented for two systems. The method was validated through calculations using appropriate experimental information for bacteriorhodopsin, which produced a model structure with a root mean square (rms) deviation of 1.87 A from the structure determined by electron microscopy. Calculations are also presented using experimental and theoretical information available for bovine rhodopsin to assign the helices to a projection density map and to produce a model of bovine rhodopsin that can be used as a template for modeling other G-protein-coupled receptors.     

NMR solution structure of KP-TerB, a tellurite-resistance protein from Klebsiella pneumoniae

Klebsiella pneumoniae (KP), a Gram-negative bacterium, is a common cause of hospital-acquired bacterial infections worldwide. Tellurium (Te) compounds, although relatively rare in the environment, have a long history as antimicrobial and therapeutic agents. In bacteria, tellurite (TeO(3) (-2)) resistance is conferred by the ter (Te(r)) operon (terZABCDEF). Here, on the basis of 2593 restraints derived from NMR analysis, we report the NMR structure of TerB protein (151 amino acids) of KP (KP-TerB), which is mainly composed of seven alpha-helices and a 3(10) helix, with helices II to V apparently forming a four-helix bundle. The ensemble of 20 NMR structures was well-defined, with a RMSD of 0.32 +/- 0.06 A for backbone atoms and 1.11 +/- 0.07 A for heavy atoms, respectively. A unique property of the KP-TerB structure is that the positively and negatively charged clusters are formed by the N-terminal positively and C-terminal negatively charged residues, respectively. To the best of our knowledge, the protein sequence and structures of KP-TerB are unique.     

Mapping of POP1-binding site on pyrin domain of ASC

Apoptosis-associated speck-like protein containing a caspase recruitment domain (ASC) is an essential adaptor protein in the formation of a multiprotein complex that activates procaspase-1. ASC is also known as a modulator of NF-kappaB activation pathways. ASC has a bipartite domain structure, consisting of an N-terminal pyrin domain (PYD) and a C-terminal caspase-recruitment domain. The PYD of ASC (ASC_PYD) is known to interact with various PYD-containing intracellular danger signal sensors and PYD-only proteins. Using purified proteins, we characterized the in vitro interaction of ASC_PYD with PYD-only protein 1 (POP1). POP1 specifically interacts with ASC_PYD with a dissociation constant of 4.08 +/- 0.52 microm but does not interact with Cryopyrin. NMR and mutagenesis experiments show that a negative electrostatic potential surface patch (EPSP) on ASC_PYD, consisting of the first (H1) and fourth (H4) helices, is essential in the interaction with POP1. A positive EPSP on POP1, consisting of the second (H2) and third (H3) helices, is a counterpart of this interaction. The interaction between ASC_PYD and POP1 is similar to the interaction between caspase recruitment domains of Apaf-1 and procaspase-9. In addition, we present evidence that conformational changes at the long loop of ASC_PYD between the H2 and H3 helices can affect its interaction with POP1. Based on our observations, we propose that the positive EPSP of ASC_PYD, including the H2 and H3 helices, may be the binding site for Cryopyrin, and the interaction with Cryopyrin may induce the dissociation of POP1 from ASC_PYD.     

Conformational properties of native sperm whale apomyoglobin in solution.

Apomyoglobin from sperm whale is often used for studies of ligand binding, protein folding, and protein stability. In an effort to describe its conformational properties in solution, homonuclear and heteronuclear (13C and 15N) NMR methods were applied to the protein in its native state. Assignments were confirmed for nuclear Overhauser effects (NOEs) involving side chain and backbone protons in the folded regions of the structure. These NOEs were used to derive distance restraints. The shifts induced by the hydrophobic dye 8-anilino-1-naphthalenesulfonic acid (ANS) were inspected in the regions remote from its binding site and served as an indicator of conformational flexibility. 3JalphaH-NH values were obtained to assess dihedral angle averaging and to provide additional restraints. A family of structures was calculated with X-PLOR and an ab initio simulated annealing protocol using holomyoglobin as a template. Where the structure appeared well defined by chemical shift, line width, ANS perturbation, and density of NOEs, the low resolution model of apomyoglobin provides a valid approximation for the structure. The new model offers an improved representation of the folded regions of the protein, which encompass the A, B, E, helices as well as parts of the G and H helices. Regions that are less well defined at this stage of calculations include the CD corner and the end of the H-helix. The EF-F-FG segment remains uncharacterized.     

Three-dimensional structure in solution of acyl-coenzyme A binding protein from bovine liver.

The three-dimensional structure of acyl-coenzyme A binding protein as encoded by the recombinant gene in Escherichia coli has been determined using nuclear magnetic resonance (n.m.r.) spectroscopy. The structure consists of four alpha-helices A1 (residues 3 to 15), A2 (residues 20 to 36), A3 (residues 51 to 60), and A4 (residues 65 to 85). A1 and A4, and A2 and A3, run in parallel pairs. A2 runs anti-parallel to A1 and A4. The three-dimensional structure of the protein is reminiscent of a shallow bowl with a rim. The "rim" is characterized by many polar and charged groups, whereas the inside and outside surface is predominantly hydrophobic with patches of uncharged polar hydroxyl groups of threonyl, serinyl and tyrosyl residues. The inside bottom contains through two epsilon-amino groups of lysine residues (Lys13 and Lys32) suggesting that the binding site for the nucleotide part of the acyl-coenzyme A part of the ligand molecule is at the inside surface of the bowl. The structure determination was done on the basis of measurements of the intensities of nuclear Overhauser effects (NOEs) and coupling constants that were translated into interatom distance restraints for 833 atom pairs, and 87 dihedral angle restraints, of which 23 were in chiral centers. In all, 42 hydrogen bonds were identified by n.m.r. and provided an additional 84 distance restraints. A total of 20 structures were calculated and the structures can be aligned to a root-mean-square deviation of 0.5 A for the backbone atoms of the residues in the four helices. A region of six residues could not be defined by the restraints obtained by n.m.r. The program Pronto was used for the spectrum analysis in general, and especially for the assignment of the individual NOEs, the integration of the cross peaks, and the measurements of the coupling constants. The programs DIANA and X-PLOR have been used in the structure calculations and evaluations.     

Model of the 3-D structure of the GLUT3 glucose transporter and molecular dynamics simulation of glucose transport

A molecular model of the three-dimensional (3-D) structure of the glucose transport protein, GLUT3, has been derived by homology modeling. The model was built on the basis of structural data from the MscL protein, which is a mechanosensitive ion channel, and general insights from aquaporin (a water permeation pore). Structurally conserved regions were defined by amino acid sequence comparisons, optimum interconnecting loops were selected from the protein databank, and amino (N)- and carboxy (C)-terminal ends of the protein were generated as random coil structures. The model was then subjected to energy minimization and molecular dynamics simulations in the presence of bound substrate (D-glucose). In the proposed structure of GLUT3, the 12 transmembrane (TM) helices form a right-hand barrel with a central hydrophilic pore. The pore is shaped like a funnel with dimensions of approximately 5-6 A by 8 A at its narrowest point. A network of polar and aromatic amino acids line the pore region and may facilitate the movement of glucose along the channel. A putative binding site for inhibitory ligands, such as forskolin and cytochalasin B, was identified on an intracellular aspect of the protein. Molecular dynamics studies showed that changes in the tilt and flexibility of key TM helices may modulate the opening of the pore to effect glucose transport. The proposed structure of GLUT3 may prove useful in guiding future experiments aimed at more precisely defining various functional regions of the transporter and may encourage efforts to develop models of other complex membrane proteins.     

Solution structure and dynamics of melanoma inhibitory activity protein

Melanoma inhibitory activity (MIA) is a small secreted protein that is implicated in cartilage cell maintenance and melanoma metastasis. It is representative of a recently discovered family of proteins that contain a Src Homologous 3 (SH3) subdomain. While SH3 domains are normally found in intracellular proteins and mediate protein-protein interactions via recognition of polyproline helices, MIA is single-domain extracellular protein, and it probably binds to a different class of ligands.Here we report the assignments, solution structure, and dynamics of human MIA determined by heteronuclear NMR methods. The structures were calculated in a semi-automated manner without manual assignment of NOE crosspeaks, and have a backbone rmsd of 0.38 Å over the ordered regions of the protein. The structure consists of an SH3-like subdomain with N- and C-terminal extensions of approximately 20 amino acids each that together form a novel fold. The rmsd between the solution structure and our recently reported crystal structure is 0.86 Å over the ordered regions of the backbone, and the main differences are localized to the most dynamic regions of the protein. The similarity between the NMR and crystal structures supports the use of automated NOE assignments and ambiguous restraints to accelerate the calculation of NMR structures.     

Use of secondary structural information and C<Emphasis Type="Italic">α</Emphasis>-C<Emphasis Type="Italic">α</Emphasis> distance restraints to model protein structures with MODELLER

Protein secondary structure predictions and amino acid long range contact map predictions from primary sequence of proteins have been explored to aid in modelling protein tertiary structures. In order to evaluate the usefulness of secondary structure and 3D-residue contact prediction methods to model protein structures we have used the known Q3 (alpha-helix, beta-strands and irregular turns/loops) secondary structure information, along with residue-residue contact information as restraints for MODELLER. We present here results of our modelling studies on 30 best resolved single domain protein structures of varied lengths. The results shows that it is very difficult to obtain useful models even with 100% accurate secondary structure predictions and accurate residue contact predictions for up to 30% of residues in a sequence. The best models that we obtained for proteins of lengths 37, 70, 118, 136 and 193 amino acid residues are of RMSDs 4.17, 5.27, 9.12, 7.89 and 9.69, respectively. The results show that one can obtain better models for the proteins which have high percent of alpha-helix content. This analysis further shows that MODELLER restrain optimization program can be useful only if we have truly homologous structure(s) as a template where it derives numerous restraints, almost identical to the templates used. This analysis also clearly indicates that even if we satisfy several true residue-residue contact distances, up to 30% of their sequence length with fully known secondary structural information, we end up predicting model structures much distant from their corresponding native structures.     

Use of secondary structural information and C alpha-C alpha distance restraints to model protein structures with MODELLER

Protein secondary structure predictions and amino acid long range contact map predictions from primary sequence of proteins have been explored to aid in modelling protein tertiary structures. In order to evaluate the usefulness of secondary structure and 3D-residue contact prediction methods to model protein structures we have used the known Q3 (alpha-helix,beta-strands and irregular turns/loops) secondary structure information, along with residue-residue contact information as restraints for MODELLER. We present here results of our modelling studies on 30 best resolved single domain protein structures of varied lengths. The results shows that it is very difficult to obtain useful models even with 100% accurate secondary structure predictions and accurate residue contact predictions for up to 30% of residues in a sequence. The best models that we obtained for proteins of lengths 37, 70, 118, 136 and 193 amino acid residues are of RMSDs 4.17, 5.27, 9.12, 7.89 and 9.69,respectively. The results show that one can obtain better models for the proteins which have high percent of alpha-helix content. This analysis further shows that MODELLER restrain optimization program can be useful only if we have truly homologous structure(s) as a template where it derives numerous restraints, almost identical to the templates used. This analysis also clearly indicates that even if we satisfy several true residue-residue contact distances, up to 30%of their sequence length with fully known secondary structural information, we end up predicting model structures much distant from their corresponding native structures.     

Role of charged and hydrophobic residues in the oligomerization of the PYRIN domain of ASC

Apoptosis-associated speck-like protein containing a caspase recruitment domain (ASC) is an adaptor protein composed of two homophilic protein-protein interaction domains, a PYRIN domain (PYD) and a caspase recruitment domain. PYD-dependent oligomerization of ASC is thought to play a crucial role in formation of a molecular platform, the inflammasome, which activates caspase-1. When expressed in cells, the PYD of ASC was shown to form cytoplasmic filaments through self-association. Over 70 single point mutants were analyzed for filament formation in cells expressing the mutant proteins. The set of mutations comprised every single amino acid residue with a charged side chain (Arg, Lys, Asp, and Glu) and a large hydrophobic side chain (Ile, Leu, Met, Phe, Pro, and Val). Filament formation of the ASC PYD was prevented by mutation of Lys21, Leu25, Lys26, Pro40, Arg41, Asp48, and Asp51 of helices 2, 3, and 4. These data identify a coherent interaction surface, establishing a molecular model of PYD-PYD complexes with an important role for charge-charge interactions.     

The three-dimensional structure of guanine-specific ribonuclease F1 in solution determined by NMR spectroscopy and distance geometry.

Two-dimensional 1H-NMR studies have been performed on ribonuclease F1 (RNase F1), which contains 106 amino acid residues. Sequence-specific resonance assignments were accomplished for the backbone protons of 99 amino acid residues and for most of their side-chain protons. The three-dimensional structures were constructed on the basis of 820 interproton-distance restraints derived from NOE, 64 distance restraints for 32 hydrogen bonds and 33 phi torsion-angle restraints. A total of 40 structures were obtained by distance geometry and simulated-annealing calculations. The average root-mean-square deviation (residues 1-106) between the 40 converged structures and the mean structure obtained by averaging their coordinates was 0.116 +/- 0.018 nm for the backbone atoms and 0.182 +/- 0.015 nm for all atoms including the hydrogen atoms. RNase F1 was determined to be an alpha/beta-type protein. A well-defined structure constitutes the core region, which consists of a small N-terminal beta-sheet (beta 1, beta 2) and a central five-stranded beta-sheet (beta 3-beta 7) packed on a long helix. The structure of RNase F1 has been compared with that of RNase T1, which was determined by X-ray crystallography. Both belong to the same family of microbial ribonucleases. The polypeptide backbone fold of RNase F1 is basically identical to that of RNase T1. The conformation-dependent chemical shifts of the C alpha protons are well conserved between RNase F1 and RNase T1. The residues implicated in catalysis are all located on the central beta-sheet in a geometry similar to that of RNase T1.     

人EEN蛋白的结构预测和功能分析

采用基于神经网络的算法预测了我们自行克隆的新的白血病相关蛋白ＥＥＮ（ｅｘｔｒａ ｅｌｅｖｅｎｎｉｎｅｔｅｅｎ, ＥＥＮ）全长分子的二级结构,结果表明：ＥＥＮ 蛋白可能有三个结构域,Ｎ 端由三段α螺旋和短β折叠组成,中间为四段α螺旋组成的四螺旋结构,Ｃ端为ＳＨ３结构域,类似于在受体酪氨酸激酶信号传导途径中起重要作用的ＳＥＭ－５／ＧＲＢ２ Ｃ端ＳＨ３结构域;利用同源蛋白结构模拟的方法,模拟了ＥＥＮ ＳＨ３结构域的三维结构,结果表明：ＥＥＮ ＳＨ３结构域与ＳＥＭ－５／ＧＲＢ２ ＳＨ３结构域具有相近的结构,构成脯氨酸结合区的氨基酸非常保守．上述结果提示：ＥＥＮ 蛋白可能为新的信号蛋白,可能涉及新的信号传导途径或新的信号传导旁路,ＳＨ３结构域是其功能区域．     

Folding determinants of disulfide bond forming protein B explored by solution nuclear magnetic resonance spectroscopy

The two-stage model for membrane protein folding postulates that individual helices form first and are subsequently packed against each other. To probe the two-stage model, the structures of peptides representing individual transmembrane helices of the disulfide bond forming protein B have been studied in trifluoroethanol solution as well as in detergent micelles using nuclear magnetic resonance (NMR) and circular dichroism spectroscopy. In TFE solution, peptides showed well-defined α-helical structures. Peptide structures in TFE were compared to the structures of full-length protein obtained by X-ray crystallography and NMR. The extent of α-helical secondary structure coincided well, lending support for the two-stage model for membrane protein folding. However, the conformation of some amino acid side chains differs between the structures of peptide and full-length protein. In micellar solution, the peptides also adopted a helical structure, albeit of reduced definition. Using measurements of paramagnetic relaxation enhancement, peptides were confirmed to be embedded in micelles. These observations may indicate that in the native protein, tertiary interactions additionally stabilize the secondary structure of the individual transmembrane helices.     

Recoverable one-dimensional encoding of three-dimensional protein structures

One-dimensional (1D) structures of proteins such as secondary structure and contact number provide intuitive pictures to understand how the native three-dimensional (3D) structure of a protein is encoded in the amino acid sequence. However, it is still not clear whether a given set of 1D structures contains sufficient information for recovering the underlying 3D structure. Here we show that the 3D structure of a protein can be recovered from a set of three types of 1D structures, namely, secondary structure, contact number and residue-wise contact order which is introduced here for the first time. Using simulated annealing molecular dynamics simulations, the structures satisfying the given native 1D structural restraints were sought for 16 proteins of various structural classes and of sizes ranging from 56 to 146 residues. By selecting the structures best satisfying the restraints, all the proteins showed a coordinate RMS deviation of <4 A from the native structure, and, for most of them, the deviation was even <2 A. The present result opens a new possibility to protein structure prediction and our understanding of the sequence-structure relationship.     

Determination of the solution structures of domains II and III of protein G from Streptococcus by 1H nuclear magnetic resonance.

We have used 1H nuclear magnetic resonance spectroscopy to determine the solution structures of two small (61 and 64 residue) immunoglobulin G (IgG)-binding domains from protein G, a cell-surface protein from Streptococcus strain G148. The two domains differ in sequence by four amino acid substitutions, and differ in their affinity for some subclasses of IgG. The structure of domain II was determined using a total of 478 distance restraints, 31 phi and 9 chi 1 dihedral angle restraints; that of domain III was determined using a total of 445 distance restraints, 31 phi and 9 chi 1 dihedral angle restraints. A protocol which involved distance geometry, simulated annealing and restrained molecular dynamics was used to determine ensembles of 40 structures consistent with these restraints. The structures are found to consist of an alpha-helix packed against a four-stranded antiparallel-parallel-antiparallel beta-sheet. The structures of the two domains are compared to each other and to the reported structure of a similar domain from a protein G from a different strain of Streptococcus. We conclude that the difference in affinity of domains II and III for IgG is due to local changes in amino acid side-chains, rather than a more extensive change in conformation, suggesting that one or more of the residues which differ between them are directly involved in interaction with IgG.     

Three-dimensional solution structure of beta cryptogein, a beta elicitin secreted by a phytopathogenic fungus Phytophthora cryptogea.

Cryptogein belongs to a new family of 10-kDa proteins called elicitins. Elicitins are necrotic and signaling proteins secreted by Phytophthora spp. responsible for the incompatible reaction and systemic hypersensitive-like necroses of diverse plant species leading to resistance against fungal or bacterial plant pathogens. The solution structure of beta cryptogein from Phytophthora cryptogea fungus was determined by using multidimensional heteronuclear nuclear magnetic resonance spectroscopy. A set of 18 structures was calculated using 1360 NOE-derived distance restraints and 40 dihedral angle restraints obtained from 3JHNH alpha couplings. The RMS deviation from the mean structure is 0.87 +/- 0.14 A for backbone atoms and 1.34 +/- 0.14 A for all the non-hydrogen atoms of residues 2 to 98. The structure of beta cryptogein reveals a novel protein fold, with five helices and a double-stranded beta-sheet facing an omega-loop. One edge of the beta-sheet and the adjacent face of the omega-loop form a hydrophobic cavity. This cavity made of highly conserved residues represents a plausible binding site. Residue 13, which has been identified from directed mutagenesis and natural sequence comparison studies as a key amino acid involved in the differential control of necrosis, is surface exposed and could contribute to the binding to a ligand or a receptor. The solution structure is close to the X-ray structure, with slight differences lightly due to the crystal packing.