Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

Generation of deviation parameters for amino acid singlets, doublets and triplets from three-dimensional structures of proteins and its implications for secondary structure prediction from amino acid sequences

We present a new method, secondary structure prediction by deviation parameter (SSPDP) for predicting the secondary structure of proteins from amino acid sequence. Deviation parameters (DP) for amino acid singlets, doublets and triplets were computed with respect to secondary structural elements of proteins based on the dictionary of secondary structure prediction (DSSP)-generated secondary structure for 408 selected nonhomologous proteins. To the amino acid triplets which are not found in the selected dataset, a DP value of zero is assigned with respect to the secondary structural elements of proteins. The total number of parameters generated is 15,432, in the possible parameters of 25,260. Deviation parameter is complete with respect to amino acid singlets, doublets, and partially complete with respect to amino acid triplets. These generated parameters were used to predict secondary structural elements from amino acid sequence. The secondary structure predicted by our method (SSPDP) was compared with that of single sequence (NNPREDICT) and multiple sequence (PHD) methods. The average value of the percentage of prediction accuracy for αhelix by SSPDP, NNPREDICT and PHD methods was found to be 57%, 44% and 69% respectively for the proteins in the selected dataset. For Β-strand the prediction accuracy is found to be 69%, 21% and 53% respectively by SSPDP, NNPREDICT and PHD methods. This clearly indicates that the secondary structure prediction by our method is as good as PHD method but much better than NNPREDICT method.     

Prediction algorithm for amino acid types with their secondary structure in proteins (PLATON) using chemical shifts

The algorithm PLATON is able to assign sets of chemical shifts derived from a single residue to amino acid types with its secondary structure (amino acid species). A subsequent ranking procedure using optionally two different penalty functions yields predictions for possible amino acid species for the given set of chemical shifts. This was demonstrated in the case of the -spectrin SH3 domain and applied to 9 further protein data sets taken from the BioMagRes database. A database consisting of reference chemical shift patterns (reference CSPs) was generated from assigned chemical shifts of proteins with known 3D-structure. This reference CSP database is used in our approach for extracting distributions of amino acid types with their most likely secondary structure elements (namely -helix, -sheet, and coil) for single amino acids by comparison with query CSPs. Results obtained for the 10 investigated proteins indicates that the percentage of correct amino acid species in the first three positions in the ranking list, ranges from 71.4% to 93.2% for the more favorable penalty function. Where only the top result of the ranking list for these 10 proteins is considered, 36.5% to 83.1% of the amino acid species are correctly predicted. The main advantage of our approach, over other methods that rely on average chemical shift values is the ability to increase database content by incorporating newly derived CSPs, and therefore to improve PLATON's performance over time.     

Generation of deviation parameters for amino acid singlets, doublets and triplets from three-dimentional structures of proteins and its implications for secondary structure prediction from amino acid sequences

We present a new method, secondary structure prediction by deviation parameter (SSPDP) for predicting the secondary structure of proteins from amino acid sequence. Deviation parameters (DP) for amino acid singlets, doublets and triplets were computed with respect to secondary structural elements of proteins based on the dictionary of secondary structure prediction (DSSP)-generated secondary structure for 408 selected non-homologous proteins. To the amino acid triplets which are not found in the selected dataset, a DP value of zero is assigned with respect to the secondary structural elements of proteins. The total number of parameters generated is 15,432, in the possible parameters of 25,260. Deviation parameter is complete with respect to amino acid singlets, doublets, and partially complete with respect to amino acid triplets. These generated parameters were used to predict secondary structural elements from amino acid sequence. The secondary structure predicted by our method (SSPDP) was compared with that of single sequence (NNPREDICT) and multiple sequence (PHD) methods. The average value of the percentage of prediction accuracy for a helix by SSPDP, NNPREDICT and PHD methods was found to be 57%, 44% and 69% respectively for the proteins in the selected dataset. For b-strand the prediction accuracy is found to be 69%, 21% and 53% respectively by SSPDP, NNPREDICT and PHD methods. This clearly indicates that the secondary structure prediction by our method is as good as PHD method but much better than NNPREDICT method.     

Photo-cross-linked hydrogels with polysaccharide-poly(amino acid) structure: new biomaterials for pharmaceutical applications

The aim of this work has been the preparation and characterization of novel hydrogels with polysaccharide-poly(amino acid) structure having suitable physicochemical properties for pharmaceutical applications. In the first step, hyaluronic acid (HA) and alpha,beta-poly(N-2-hydroxyethyl)-DL-aspartamide (PHEA) have been derivatized with methacrylic anhydride (AMA), thus obtaining HA-AMA and PHM derivatives, respectively. In the second step, aqueous solutions of both these derivatives have been irradiated at 313 nm to obtain chemical hydrogels. The hydrogel obtained by irradiating for 15 min an aqueous solution containing 4% w/v of HA-AMA and 4% w/v of PHM resulted in the highest yield. Its swelling ability was dependent on the pH and nature of the external medium. Besides, this hydrogel undergoes a partial hydrolysis, especially in the presence of enzymes, such as esterase or hyaluronidase, but the entity of this degradation is lower than that observed for a hydrogel based on HA-AMA alone. The ability of this hydrogel to entrap drug molecules has been evaluated by using thrombin as a model drug. In vitro release studies and a platelet aggregation test demonstrated that the HA-AMA/PHM hydrogel is able to release thrombin in the active form, thus suggesting its suitability for the treatment of hemorrhages.     

Cloning of genes encoding redox proteins of known amino acid sequence from a library of the Desulfovibrio vulgaris (Hildenborough) genome

A library of 900 recombinant phages has been constructed for the genome of Desulfovibrio vulgaris Hildenborough (1.7 x 10(6) bp) by cloning size-fractionated Sau3A fragments (15-20 kb) into the replacement vector lambda-2001. When a hydrogenase gene probe, a 4.7-kb SalI-EcoRI fragment of known nucleotide sequence, was used to screen the plaque lifted library, 23 positive clones were found, which together span 31 kb of D. vulgaris DNA. To facilitate the cloning of genes with oligodeoxynucleotides as probes, DNA was purified for all clones in the library and spotted on a 16 x 16-cm grid of nitrocellulose. This grid was incubated sequentially to identify lambda clones containing the gene for redox proteins of known amino acid sequence: cytochrome c3 (one 18-mer----four clones), flavodoxin (one 17-mer and one 26-mer----one clone) and rubredoxin (one 44-mer----21 clones). The four cyc-positive clones are also recognized by the rubredoxin oligodeoxynucleotide probe. Restriction mapping defines a 35-kb region of the D. vulgaris chromosome in which the rub and cyc loci are separated by 17.5 kb. The nucleotide sequence of the rubredoxin gene was determined and the deduced amino acid sequence found to agree with that determined in Bruschi [Biochim. Biophys. Acta 434 (1976) 4-17] with the exception of Thr-21 which is found to be encoded by GAC, an Asp codon. A plausible ribosome-binding site precedes the N-terminal initiator methionine residue. Rubredoxin does not have an N-terminal signal sequence which is in agreement with the cytoplasmic location of this redox protein.     

An attempt to evaluate the influence of neighboring amino acids (n-1) and (n+1) on the backbone conformation of amino acid (n) in proteins. Use in predicting the three-dimensional structure of the polypeptide backbone of other proteins

From φ, ψ data on eleven proteins, a 20 × 20 × 20 Table of tripeptides has been computed to evaluate the influence of nearest neighbors (n ? 1) and (n + 1) on the φ, ψ angles of amino acid (n). From this Table, having removed values for horse cytochrome c and using the sequences of 18 cytochromes c and the procedure of Kabat &; Wu (1972), an attempt was made to select a set of φ, ψ angles for positions 2 to 103 of cytochrome c and compare them with the values obtained from the atomic co-ordinates. Agreement was good for 56, intermediate for 29 and poor for 17 residues. Eleven of the 17 with poor agreement were residues contacting the heme or adjacent to a contacting residue. Moreover, 6 of the 17 poor values were in regions of the φ, ψ plot for which no occurrences in the ten known proteins were reported, and for four others known values were minimal so that no basis for selection existed. Frequency distributions on the Ramachandran plot (Ramachandran &; Sasisekharan, 1968) of all φ, ψ values in the eleven known proteins are given as well as a contour plot for such frequencies. The uses and limitations of the procedure are discussed and the need for obtaining accurate estimates of errors in experimentally determined φ, ψ angles is emphasised.     

A cholinesterase genes server (ESTHER): a database of cholinesterase-related sequences for multiple alignments, phylogenetic relationships, mutations and structural data retrieval.

We have built a database of sequences phylogenetically related to cholinesterases (ESTHER) for esterases, alpha/beta hydrolase enzymes and relatives). These sequences define a homogeneous group of enzymes (carboxylesterases, lipases and hormone-sensitive lipases) with some related proteins devoid of enzymatic activity. The purpose of ESTHER is to help comparison and alignment of any new sequence appearing in the field, to favour mutation analysis of structure-function relationships and to allow structural data recovery. ESTHER is a World Wide Web server with the URL http://www.montpellier.inra.fr:70/cholinesterase.     

Selenomethionyl proteins produced for analysis by multiwavelength anomalous diffraction (MAD): a vehicle for direct determination of three-dimensional structure.

An expression system has been established for the incorporation of selenomethionine into recombinant proteins produced from plasmids in Escherichia coli. Replacement of methionine by selenomethionine is demonstrated at the level of 100% for both T4 and E. coli thioredoxins. The natural recombinant proteins and the selenomethionyl variants of both thioredoxins crystallize isomorphously. Anomalous scattering factors were deduced from synchrotron X-ray absorption measurements of crystals of the selenomethionyl proteins. Taken with reference to experience in the structural analysis of selenobiotinyl streptavidin by the method of multiwavelength anomalous diffraction (MAD), these data indicate that recombinant selenomethionyl proteins analyzed by MAD phasing offer a rather general means for the elucidation of atomic structures.     

Novel detection method for human protein C inhibitor (PCI) using immunoblotting with antibodies specific to synthetic peptides corresponding to the N-terminal and C-terminal amino acid sequences of human PCI.

Three types of antibody, which were specific to the purified human protein C inhibitor (PCI) and to two synthetic peptides, which corresponded to the N-terminal 15 amino acid residues (PCI-N) and C-terminal 15 residues (PCI-C) of PCI, were produced. Plasma PCI patterns, which were shown by means of polyacrylamide gel isoelectric focusing followed by immunoblotting with the antibody to purified PCI, were almost identical to those obtained with the antibodies to the peptides, PCI-N and PCI-C. The latter antibodies against the synthetic peptides proved very useful for biochemical or genetic analysis of human plasma PCI, as such peptide antigens could be produced without any of the tedious procedures necessary to obtain purified PCI as an immunogen.     

The effect of zinc on the secondary structure of d(GA.TC)n DNA sequences of different length: a model for the formation *H-DNA.

Alternating d(GA.TC)n DNA sequences are known to undergo transition to *H-DNA in the presence of zinc. Here, the effect of zinc on the secondary DNA structure of d(GA.TC)n sequences of different length (n = 5, 8, 10 and 19) was determined. Short d(GA.TC)n sequences form *H-DNA with a higher difficulty than longer ones. At bacterial negative superhelical density (- sigma = 0.05), zinc still induces transition to the *H-DNA conformation at a d(GA.TC)10 sequence but shorter sequences do not form *H-DNA. Transition to *H-DNA at a d(GA.TC)8 sequence is observed under conditions which destabilize the DNA double helix such as high negative supercoiling or low ionic strength. Our results indicate that a first step in the transition to *H-DNA is the formation of a denaturation bubble at the centre of the repeated DNA sequence, suggesting that the primary role of zinc is to induce a local denaturation of the DNA double helix. Subsequently, zinc might also participate in the stabilization of the altered DNA conformation through its direct interaction with the bases. Based on these results a model for the formation of *H-DNA is proposed.     

Leucyl/phenylalanyl(L/F)-tRNA-protein transferase-mediated aminoacyl transfer of a nonnatural amino acid to the N-terminus of peptides and proteins and subsequent functionalization by bioorthogonal reactions

We report here a new strategy for derivatizing peptides and proteins at the N-terminus. To achieve this, a nonnatural amino acid was charged onto a tRNA and then enzymatically transferred to a lysine (Lys) unit at the N-terminus of a peptide or a protein by using L/F-tRNA-protein transferase. By using the chemoenzymatic technique, beta-(2-quinolyl)-L-alanine, p-azido-L-phenylalanine, and p-acetyl-L-phenylalanine were introduced to the N-terminus. The latter two nonnatural amino acids possess bioorthogonal functional groups to which artificial tags can be introduced. Actually, a biotin tag was coupled to the bioorthogonal ketone group of acetylphenylalanine at the N-terminus of a peptide. N-terminal-specific biotinylation and fluorescence derivatization of the bioorthogonal azido-containing protein or peptide was also carried out based on a [3 + 2] cycloaddition. The enzymatic transfer of a nonnatural amino acid to the N-terminus of target peptides or proteins was also successfully achieved in the presence of other peptides or crude protein mixtures.     

Fibrinopeptides A and B of baboons (Papio anubis, Papio hamadryas, and Theropithecus gelada): their amino acid sequences and evolutionary rates and a molecular phylogeny for the baboons

Amino acid sequences of fibrinopeptides A and B from savanna baboons, Papio anubis and Papio hamadryas, and highland baboon, Theropithecus gelada, were established. The sequences of the fibrinopeptides A from the three baboons were identical: (sequence: see text) The fibrinopeptides B were composed of 9 residues and demonstrated the sequence: (sequence see text) where X3 = Arg in P. anubis, His in P. hamadryas, and Gly in Th. gelada. Position-3 of the B peptides was the only replacement site observed among the 25 amino acid residues in both fibrinopeptides from the baboons. Based on these sequences, a molecular phylogeny for the three species of baboons was deduced. The evolutionary rates of the peptides B of the baboons and macaques were also estimated. It was observed that the fibrinopeptides changed at an uneven rate during the evolution of old-world monkeys, i.e., baboons and macaques.     

Single amino acid of g(alpha16) (Ala(228)) is responsible for the ability of chemoattractant C5a receptor to induce G(alpha16)-mediated inositol phosphate release

Our previous study suggested that the region encompassing residues 220-240 on G(alpha16) is important in coupling with C5a receptor (Lee et al. (1995) Mol. Pharmacol. 47, 218-223). When aligned sequences are compared in the residue 220-240 segment of G(alpha16), there is a block of eight amino acids extending from residue 227 to residue 234 (227-Ile-Ala-Leu-Ile-Tyr-Leu-Ala-Ser-234) in G(alpha16) that is replaced by a heterologous block extending from amino acid residue 224 to residue 231 (224-Thr-Ser-Ile-Met-Phe-Leu-Val-Ala-231) in G(alpha11). In order to identify the specific amino acid residue necessary for coupling to C5a receptor within the extension of eight amino acids in G(alpha16), a series of chimeric G(alpha11)/G(alpha16) cDNA constructs and mutant G(alpha16) cDNAs were expressed. Then the ability of chimeras and mutant proteins to mediate C5a-induced release of inositol phosphate in transfected Cos-7 cells was tested. The results show that single amino acid Ala(228) is responsible for conferring about 40-50% of the activity of G(alpha16) induced by C5a receptor stimulation.     

E(z), a depth-dependent potential for assessing the energies of insertion of amino acid side-chains into membranes: derivation and applications to determining the orientation of transmembrane and interfacial helices

We have developed an empirical residue-based potential (E(z) potential) for protein insertion in lipid membranes. Propensities for occurrence as a function of depth in the bilayer were calculated for the individual amino acid types from their distribution in known structures of helical membrane proteins. The propensities were then fit to continuous curves and converted to a potential using a reverse-Boltzman relationship. The E(z) potential demonstrated a good correlation with experimental data such as amino acid transfer free energy scales (water to membrane center and water to interface), and it incorporates transmembrane helices of varying composition in the membrane with trends similar to those obtained with translocon-mediated insertion experiments. The potential has a variety of applications in the analysis of natural membrane proteins as well as in the design of new ones. It can help in calculating the propensity of single helices to insert in the bilayer and estimate their tilt angle with respect to the bilayer normal. It can be utilized to discriminate amphiphilic helices that assume a parallel orientation at the membrane interface, such as those of membrane-active peptides. In membrane protein design applications, the potential allows an environment-dependent selection of amino acid identities.     

New general approach for determining the solution structure of a ligand bound weakly to a receptor: structure of a fibrinogen Aalpha-like peptide bound to thrombin (S195A) obtained using NOE distance constraints and an ECEPP/3 flexible docking program.

A new approach incorporating flexible docking simulations and NMR data is presented for calculating the bound conformation of a ligand that interacts weakly with an enzyme. This approach consists of sampling directly the conformation of a flexible ligand inside a receptor active site containing surrounding flexible loops. To make this sampling efficient, a ligand-growing procedure has been adopted. Optimization of the ECEPP/3-plus-NOE constraint function is carried out by using a collective variable Monte Carlo minimization technique. Numerous energy minimizations are made possible for such a large system by using a Bezier splines energy grid technique. This new flexible docking approach was applied to determine the structure of a fibrinogen Aalpha-like peptide (7DFLAEGGGVRGPRV20) bound to an active site mutant of thrombin [thrombin(S195A)]. Structure calculations of the bound ligand, using 2D-transferred NOESY distance constraints in the DIANA program, showed that the N-terminal portion of the peptide (D7-R16) involves a chain reversal, whereas the C-terminal portion (G17-V20) adopts a fold that exists in several different orientations. In addition, the ECEPP/3 flexible docking package was used to assess the conformational variability of the ligand and surrounding 60D-insertion loop of thrombin. Amino acid residues (17-20) of the peptide interact with a region of the enzyme that exhibits broad specificity, with a preferred direction between the 60D-insertion loop and Pro37 of thrombin.     

Accessibility of lysyl residues of Escherichia coli B/r porin (OmpF) to covalent labeling reagents of different sizes. An approach for a three-dimensional structure of a channel-forming protein

The three-dimensional structure of Escherichia coli B/r porin (OmpF) was studied by chemical modification using activated sugars of different size. Galactose and galactosides of different penetration properties through the porin channel were oxidized by galactose oxidase, and the 6-aldehydes formed were linked to amino groups in porin by reduction with NaBH3CN. Tryptic fragments of modified and unmodified porin were separated by reversed-phase high pressure liquid chromatography and identified by amino acid and amino-terminal analysis from the known primary structure of OmpF. Modification of purified native porin trimers in beta-octylglucoside revealed three classes of amino groups: (i) those not modified by any sugars; (ii) those modified only by small sugars that diffuse rapidly through the pore, such as galactose or melibiose; and (iii) those modified by either small or large sugars, the latter including pore-impermeant sugars such as stachyose. The results suggest that the three classes of amino groups correspond, respectively, to groups buried in the trimeric molecule, those in the interior of the pore and those exposed on the surface of porin. In addition modification experiments performed on whole cells suggested that all the reactive groups modified by the pore-impermeant sugars (class iii) are located on the surface of porin exposed on the outside of the outer membrane.     

Single crystal EPR study of electronic structure and exchange interactions for copper(II)(L-arginine)2(SO4).(H2O)6: a model system to study exchange interactions between unpaired spins in proteins

We report EPR measurements at 9.77 and 34.1 GHz in powder and single crystal samples of the ternary copper amino acid complex Cu(L-arginine)(2)(SO(4)).(H(2)O)(6). The single crystal Electron Paramagnetic Resonance spectra display a single resonance for all magnetic field orientations in the ca and cb crystal planes. In the ab plane they display two resonances for most orientations of the magnetic field, and only one resonance for orientations close to the crystal axes. This behavior is a result of the selective collapse of the resonances corresponding to the four copper sites in the unit cell produced by the exchange interactions between copper ions. From the characteristics of the collapse and the angular dependences of the position and width of the resonances we evaluate the g-tensors of the copper molecules and estimate exchange interactions |J(1)/k(B)|=0.9 K and |J(2)/k(B)|=0.009 K between copper neighbors at 5.908 A and at 15.684 A, respectively. J(1) is assigned to a syn-anti equatorial-apical carboxylate bridge with a total bond length of 7.133 A. J(2) is assigned to a long bridge of 12 atoms with a total bond length of 19.789 A, that includes two hydrogen bonds. The results are discussed in terms of the crystal and electronic structure of Cu(L-arginine)(2)(SO(4)).(H(2)O)(6). We show that J(2) is in excellent agreement with the observed magnetic interaction between the reduced quinone acceptors in the photosynthetic reaction center protein of the bacterium Rb. sphaeroides, which is transmitted along a similar chemical path containing two hydrogen bonds. Our findings indicate that it is valid to estimate values for the exchange interactions between redox centers in proteins transmitted along long chemical paths containing sigma and H-bonds, from data obtained in model systems, and emphasize the importance of measuring exchange interactions in biologically relevant model systems.     

Functional importance of polar and charged amino acid residues in transmembrane helix 14 of multidrug resistance protein 1 (MRP1/ABCC1): identification of an aspartate residue critical for conversion from a high to low affinity substrate binding state

Human multidrug resistance protein 1 (MRP1) confers resistance to many chemotherapeutic agents and transports diverse conjugated organic anions. We previously demonstrated that Glu1089 in transmembrane (TM) 14 is critical for the protein to confer anthracycline resistance. We have now assessed the functional importance of all polar and charged amino acids in this TM helix. Asn1100, Ser1097, and Lys1092, which are all predicted to be on the same face of the helix as to Glu1089, are involved in determining the substrate specificity of the protein. Notably, elimination of the positively charged side chain of Lys1092, increased resistance to the cationic drugs vincristine and doxorubicin, but not the electroneutral drug etoposide (VP-16). In addition, mutations S1097A and N1100A selectively decreased transport of 17beta-estradiol 17-(beta-d-glucuronide) (E217betaG) but not cysteinyl leukotriene 4 (LTC4), demonstrating the importance of multiple residues in this helix in determining substrate specificity. In contrast, mutations of Asp1084 that eliminate the carboxylate side chain markedly decreased resistance to all drugs tested, as well as transport of both E217betaG and LTC4, despite the fact that LTC4 binding was unaffected. We show that these mutations prevent the ATP-dependent transition of the protein from a high to low affinity substrate binding state and drastically diminish ADP trapping at nucleotide binding domain 2. Based on results presented here and crystal structures of prokaryotic ATP binding cassette transporters, Asp1084 may be critical for interaction between the cytoplasmic loop connecting TM13 and TM14 and a region of nucleotide binding domain 2 between the conserved Walker A and ABC signature motifs.