<Emphasis Type="Italic">π</Emphasis>–<Emphasis Type="Italic">π</Emphasis> interaction between aromatic ring and copper-coordinated His81 imidazole regulates the blue copper active-site structure

Noncovalent weak interactions play important roles in biological systems. In particular, such interactions in the second coordination shell of metal ions in proteins may modulate the structure and reactivity of the metal ion site in functionally significant ways. Recently, π–π interactions between metal ion coordinated histidine imidazoles and aromatic amino acids have been recognized as potentially important contributors to the properties of metal ion sites. In this paper we demonstrate that in pseudoazurin (a blue copper protein) the π–π interaction between a coordinated histidine imidazole ring and the side chains of aromatic amino acids in the second coordination sphere, significantly influences the properties of the blue copper site. Electronic absorption and electron paramagnetic resonance spectra indicate that the blue copper electronic structure is perturbed, as is the redox potential, by the introduction of a second coordination shell π–π interaction. We suggest that the π–π interaction with the metal ion coordinated histidine imidazole ring modulates the electron delocalization in the active site, and that such interactions may be functionally important in refining the reactivity of blue copper sites. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

The contribution of proline residues to protein stability is associated with isomerization equilibrium in both unfolded and folded states

It has long been understood that the proline residue has lower configurational entropy than any other amino acid residue due to pyrrolidine ring hindrance. The peptide bond between proline and its preceding amino acid (Xaa-Pro) typically exists as a mixture of cis- and trans-isomers in the unfolded protein. Cis–trans isomerization of Xaa-Pro peptide bonds are infrequent, but still occur in folded proteins. Therefore, the effects of the cis–trans isomerization equilibrium in both unfolded and folded states should be taken into account when estimating the stability contribution of a specific proline residue. In order to study the stability contribution of the four proline residues to the hyperthermophilic protein Ssh10b, in this work, we expressed and purified a series of Pro→Ala mutants of Ssh10b, and performed correlative unfolding experiments in detail. We proposed a new unfolding model including proline isomerization. The model predicts that the contribution of a proline residue to protein stability is associated with the thermodynamic equilibrium between cis- and trans-isomers both in the unfolded and folded states, agreeing well with the experimental results.     

ATP as effector of inorganic pyrophosphatase of <Emphasis Type="Italic">Escherichia coli</Emphasis>. Identification of the binding site for ATP

The interaction of Escherichia coli inorganic pyrophosphatase (E-PPase) with effector ATP has been studied. The E-PPase has been chemically modified with the dialdehyde derivative of ATP. It has been established that in the experiment only one molecule of effector ATP is bound to each subunit of the hexameric enzyme. Tryptic digestion of the adenylated protein followed by isolation of a modified peptide by HPLC and its mass-spectrometric identification has showed that it is an amino group of Lys146 that undergoes modification. Molecular docking of ATP to E-PPase indicates that the binding site for effector ATP is located in a cluster of positively charged amino acid residues proposed earlier on the basis of site-directed mutagenesis to participate in binding of effector pyrophosphate. Molecular docking also reveals several other amino acid residues probably involved in the interaction with effectors. Published in Russian in Biokhimiya, 2007, Vol. 72, No. 1, pp. 110–117.     

Cell cycle and activation of CK2

Microtubule associated protein tau is considered to play roles in some types of human transmissible spongiform encephalopathies (TSE). In this study, the full-length and several truncated human tau proteins were expressed from E. coli and purified. Using GST pull down, co-immunoprecipitation assay and tau-coated ELISA, the molecular interaction between tau protein and PrP was confirmed in the context of the full-length human tau. The N terminus (amino acids 1–91) and tandem repeats region (amino acids 186–283) of tau protein were responsible for the interaction with PrP. The octapeptide repeats within PrP directly affected the binding activity of PrP with tau. GSS-related mutant PrP102L and fCJD- related mutants with two and seven extra octarepeats showed more active binding capacity with tau than wild-type PrP. The molecular interactions between PrP and tau protein highlight a potential role of tau in the biological function of PrP and the pathogenesis of TSE.     

Contribution of bovine lactoferrin inter-lobe region to iron binding stability and antimicrobial activity against <Emphasis Type="Italic">Staphylococcus aureus</Emphasis>

The investigation of the recombinant bovine lactoferrin-derived antimicrobial protein (rBLfA) demonstrates that the inter-lobe region of bovine lactoferrin contributes to iron binding stability and antimicrobial activity against Staphylococcus aureus. rBLfA containing N-lobe (amino acid residues 1–333) and inter-lobe region (residues 334–344) was expressed in Pichia pastoris at shaking flask and fermentor level. The recombinant intact bovine lactoferrin (rBLf) and N-lobe (rBLfN) were expressed in the same system as control. The physical–chemical parameters of rBLfA, rBLfN and rBLf including amino acid residues, molecular weight, isoelectric point, net positive charge and instability index were computed and compared. The simulated tertiary structure and the calculated surface net charge showed that rBLfA maintained original structure and exhibited a higher cationic feature than rBLf and rBLfN. The three proteins showed different iron binding stability and antimicrobial activity. rBLfA released iron in the pH range of 7.0–3.5, whereas rBLfN lost its iron over the pH range of 7.0–4.0 and iron release from rBLf occurred in the pH range of 5.5–3.0. However, the minimum inhibition concentration of rBLfA against S. aureus ATCC25923 was 6.5 μmol/L, compared with 12.5 and 25 μmol/L that of rBLfN and rBLf, respectively. These results revealed that S. aureus was more sensitive to rBLfA than rBLfN and rBLf. It appeared that the strong cationic character of inter-lobe region related positively to the higher anti-S. aureus activity.     

Stacking efficiency and flexibility analysis of aromatic amino acids in cap-binding proteins

Recognition of the ribonucleic acid 5' termini (RNA 5' cap) by a wide class of cap-binding proteins is largely accomplished by cation-pi stacking that involves the positively charged 7-methylguanine ring and aromatic amino acid side chains. Quantum calculations of the stacking energy were performed by means of MP2 perturbation method for binary and ternary associates composed of the 7-methylguanine moiety and tryptophan, tyrosine, or phenylalanine, in their spatial orientations known from the crystalline cap-protein complexes. The results clearly pointed to an enhancement of the stacking energy due to a net positive charge in the cap guanine moiety and allowed analysis of a role of various amino acids in stabilization of the complexes. Conformational flexibility of the aromatic amino acids taking part in binding ligands to a wide class of RNA-recognizing proteins, including the cap-binding proteins, was determined by regional order neural network (RONN) algorithm that provides results close to those of the crystallographic B-factors analysis. Interestingly, some of the tyrosines that are classified in general as "rigid" showed high flexibility when engaged in binding the cap to nuclear cap-binding protein complex CBC and to viral methyltransferase VP39. Parallel analyses of the binding energy and flexibility of the protein fragments engaged in the binding leads to understanding differences in molecular mechanisms of the cap recognition by various proteins, CBC compared with the eukaryotic initiation factor eIF4E, and enzymes vs. other protein factors.     

Binding of phosphorylated peptides and inhibition of their interaction with disease‐relevant human proteins by synthetic metal‐chelate receptors

The modulation of biological signal transduction pathways by masking phosphorylated amino acid residues represents a viable route toward pharmacologic protein regulation. Binding of phosphorylated amino acid residues has been achieved with synthetic metal‐chelate receptors. The affinity and selectivity of such receptors can be enhanced if combined with a second binding site. We demonstrate this principle with a series of synthetic ditopic metal‐chelate receptors, which were synthesized and investigated for their binding affinity to phosphorylated short peptides under conditions of physiological pH. The compounds showing highest affinity were subsequently used to inhibit the interaction of the human STAT1 protein to a peptide derived from the interferon‐γ receptor, and between the checkpoint kinase Chk2 and its preferred binding motif. Two of the investigated ditopic synthetic receptors show a significant increase in inhibition activity. The results show that regulation of protein function by binding to phosphorylated amino acids is possible. The introduction of additional binding sites into the synthetic receptors increases their affinity, but the flexibility of the structures investigated so far prohibited stringent amino acid sequence selectivity in peptide binding. Copyright © 2009 John Wiley & Sons, Ltd.     

Stacking interaction and its role in kynurenic acid binding to glutamate ionotropic receptors

Stacking interaction is known to play an important role in protein folding, enzyme-substrate and ligand-receptor complex formation. It has been shown to make a contribution into the aromatic antagonists binding with glutamate ionotropic receptors (iGluRs), in particular, the complex of NMDA receptor NR1 subunit with the kynurenic acid (KYNA) derivatives. The specificity of KYNA binding to the glutamate receptors subtypes might partially result from the differences in stacking interaction. We have calculated the optimal geometry and binding energy of KYNA dimers with the four types of aromatic amino acid residues in Rattus and Drosophila ionotropic iGluR subunits. All ab initio quantum chemical calculations were performed taking into account electron correlations at MP2 and MP4 perturbation theory levels. We have also investigated the potential energy surfaces (PES) of stacking and hydrogen bonds (HBs) within the receptor binding site and calculated the free energy of the ligand-receptor complex formation. The energy of stacking interaction depends both on the size of aromatic moieties and the electrostatic effects. The distribution of charges was shown to determine the geometry of polar aromatic ring dimers. Presumably, stacking interaction is important at the first stage of ligand binding when HBs are weak. The freedom of ligand movements and rotation within receptor site provides the precise tuning of the HBs pattern, while the incorrect stacking binding prohibits the ligand-receptor complex formation.     

N-Glycans protect proteins from protease digestion through their binding affinities for aromatic amino acid residues

It was previously revealed [Yamaguchi, H., Nishiyama, T., and Uchida, M. (1999) J. Biochem. 126, 261-265] that N-glycans of both the high-mannose and complex types have binding affinity for aromatic amino acid residues. This study shows that free N-glycans protect proteins from protease digestion through their binding affinities for the aromatic amino acid residues exposed on protein molecules. Protease digestion of bovine pancreatic RNase A and bovine a-lactalbumin was depressed in solutions (1 mM or so) of free N-glycans of both the high-mannose and complex types. The increasing order of the protective effects of the N-glycans paralleled that of their affinities for aromatic amino acid residues; and the presence of aromatic amino acids practically abolished the protective effects of the N-glycans. The N-glycans also depressed the protease digestion of metallothionein, an aromatic amino acid-free protein, in agreement with the observation that the N-glycans also interact with the solvent-exposed aromatic amino acid residues of the proteases. Thus it seems probable that the N-glycans protect proteins from protease digestion by steric hindrance attributable to their binding affinity for the solvent-exposed aromatic amino acid residues of both substrate proteins and proteases.     

Expression of biologically active recombinant antifreeze protein His-MpAFP149 from the desert beetle (<Emphasis Type="Italic">Microdera punctipennis dzungarica</Emphasis>) in <Emphasis Type="Italic">Escherichia coli</Emphasis>

An insect antifreeze protein gene Mpafp149 was cloned by the RT-PCR approach from the desert beetle Microdera punctipennis dzungarica. Sequence analysis revealed that this gene encoding a protein of 120 amid acids and this protein showed 65–76% homology with other insect antifreeze proteins, the deduced amino acid sequence displays very high similarities in those regions that contain tandem the 12-residue repeats (TCTxSxxCxxAx) domain and the TCT motif. Mpafp149 gene was cloned into pET-28a vector and expressed in Escherichia coli. A single-step purification based on specific binding of histidine residues was achieved. The purified His-MpAFP149 was SDS–PAGE analyzed, showing an atypical migration with molecular weight of about 24 kDa. The expression of His-MpAFP149 was confirmed by Western blot with specific binding to anti-GST-MpAFP149 antibody. The thermal hysteresis activity of the purified recombinant protein was 0.915°C at 0.09 mg/ml, and the supercooling point was −9.6°C at 0.03 mg/ml. In vitro antifreeze activity assay by measuring the survival rate of bacteria at −7 and −20°C respectively, with the protection of His-MpAFP149 showed that the His-MpAFP149 fusion protein was able to enhance the freeze resistance of bacteria.     

The interaction of La<Superscript>3+</Superscript> complexes of DOTA/DTPA glycoconjugates with the RCA<Subscript>120</Subscript> lectin: a saturation transfer difference NMR spectroscopic study

The study of ligand–receptor interactions using high-resolution NMR techniques, namely the saturation transfer difference (STD), is presented for the recognition process between La(III) complexes of 1,4,7,10-tetrakis(carboxymethyl)-1,4,7,10-tetraazacyclododecane monoamide and diethylenetriaminepentaacetic acid bisamide glycoconjugates and the galactose-specific lectin Ricinus communis agglutinin (RCA₁₂₀). This new class of Gd(III)-based potential targeted MRI contrast agents (CAs), bearing one or two terminal sugar (galactosyl or lactosyl) moieties, has been designed for in vivo binding to the asialoglycoprotein receptor, which is specifically expressed at the surface of liver hepatocytes, with the aim of leading to a new possible diagnosis of liver diseases. The in vitro affinity constants for the affinity of the divalent La(III)–glycoconjugate complexes for RCA₁₂₀, used as a simple, water-soluble receptor model, were higher than those of the monovalent analogues. The combination of the experimental data obtained from the STD NMR experiments with molecular modelling protocols (Autodock 4.1) allowed us to predict the mode of binding of monovalent and divalent forms of these CAs to the galactose 1α binding sites of RCA₁₂₀. The atomic details of the molecular interactions allowed us to corroborate and supported the interaction of both sugar moieties and the linkers with the surface of the protein and, thus, their contribution to the observed interaction stabilities.     

In silico cloning,expression of Rieske-like apoprotein gene and protein subcellular localization in the Pacific oyster, <Emphasis Type="Italic">Crassostrea gigas</Emphasis>

Rieske protein gene in the Pacific oyster Crassostrea gigas was obtained by in silico cloning for the first time, and its expression profiles and subcellular localization were determined, respectively. The full-length cDNA of Cgisp is 985 bp in length and contains a 5′- and 3′-untranslated regions of 35 and 161 bp, respectively, with an open reading frame of 786 bp encoding a protein of 262 amino acids. The predicted molecular weight of 30 kDa of Cgisp protein was verified by prokaryotic expression. Conserved Rieske [2Fe–2S] cluster binding sites and highly matched-pair tertiary structure with 3CWB_E (Gallus gallus) were revealed by homologous analysis and molecular modeling. Eleven putative SNP sites and two conserved hexapeptide sequences, box I (THLGC) and II (PCHGS), were detected by multiple alignments. Real-time PCR analysis showed that Cgisp is expressed in a wide range of tissues, with adductor muscle exhibiting the top expression level, suggesting its biological function of energy transduction. The GFP tagging Cgisp indicated a mitochondrial localization, further confirming its physiological function.     

Synthesis and Structure-Activity Relationship of [Nle10]Neurokinin A (4–10) Analogs with Constraint in the Backbone and at Position Six

Steric requirements of binding [Nle¹⁰]NKA(4–10) to NK-2 receptor were studied by introducing conformationally constrained amino acid analogs into its sequence. Two series of [Nle¹⁰]NKA(4–10) analogs were synthesized to investigate (i) the significance of a putative β-turn in the receptor-ligand interaction by insertion of either (S)- or (R)-Gly⁸{ANC-2}Leu⁹ γ-lactams to mimic a β-turn constraint, and (ii) the effect of hindered rotation in the Φ, χ¹ and χ² dihedral angle space of the crucially important Phe⁶ which was replaced systematically with d-Phe, d- and l-Tyr, as well as with their conformationally constrained analogs, Tic, HOTic and β-MePhe. Competition binding experiments with [³H]NKA were performed using cloned human NK-2 receptors expressed in CHO cells. The analog possessing only an (R)-Gly⁸{ANC-2}Leu⁹ constraint, had the same binding affinity as that of the parent peptide. The rank order of potency of the other analogs showed a cumulative effect of different structural modifications in decreasing the binding affinity, i.e., when changing the configuration of the lactam ring to S, replacing Phe⁶ with constrained analogs, Tic or β-MePhe, changing the configuration of the amino acid at position six to d, and introducing a hydroxyl group on the aromatic ring. Ferenc ?tv?s and Dmitry S. Gembitsky - Made an equal contribution. Abbreviations used for amino acids and peptides follow the recommendations of the IUPAC-IUB Commission of Biochemical Nomenclature, Eur. J. Biochem. (1984) 138, 9–37     

New actinoporins from sea anemone <Emphasis Type="Italic">Heteractis crispa</Emphasis>: Cloning and functional expression

A new actinoporin Hct-S4 (molecular mass 19,414 ± 10 Da) belonging to the sphingomyelin-inhibited α-pore forming toxin (α-PFT) family was isolated from the tropical sea anemone Heteractis crispa (also called Radianthus macrodactylus) and purified by methods of protein chemistry. The N-terminal nucleotide sequence (encoding 20 amino acid residues) of actinoporin Hct-S4 was determined. Genes encoding 18 new isoforms of H. crispa actinoporins were cloned and sequenced. These genes form a multigene Hct-S family characterized by presence of N-terminal serine in the mature proteins. Highly conserved residues comprising the aromatic phosphorylcholine-binding site and significant structure-function changes in the N-terminal segment (10–27 amino acid residues) of actinoporins were established. Two expressed recombinant actinoporins (rHct-S5 and rHct-S6) were one order less hemolytically active than native actinoporins.     

Quantum mechanical origin of the conformational preferences of 4-thiaproline and its <Emphasis Type="Italic">S</Emphasis>-oxides

The saturated ring and secondary amine of proline spawn equilibria between pyrrolidine ring puckers as well as peptide bond isomers. These conformational equilibria can be modulated by alterations to the chemical architecture of proline. For example, C^γ in the pyrrolidine ring can be replaced with sulfur, which can be oxidized either stereoselectively to yield diastereomeric S-oxides or completely to yield a sulfone. Here, the thiazolidine ring and peptide bond conformations of 4-thiaproline and its S-oxides were analyzed in an Ac–Xaa–OMe system using NMR spectroscopy, X-ray crystallography, and hybrid density functional theory. The results indicate that the ring pucker of the S-oxides is governed by the gauche effect, and the prolyl peptide bond conformation is determined by the strength of the n → π* interaction between the amide oxygen and the ester carbonyl group. These findings, which are consistent with those of isologous 4-hydroxyprolines and 4-fluoroprolines, substantiate the importance of electron delocalization in amino acid conformation.     

Investigation of the stereoselectivity of an anti-amino acid antibody using molecular modeling and ligand docking

The structure of the binding site of the stereoselective anti-D-amino acid antibody 67.36 was modeled utilizing web antibody modeling (WAM) and SWISS-MODEL. Although docking experiments performed with an aromatic amino acid as model ligand were unsuccessful with the WAM structure, ligand binding was achieved with the SWISS-MODEL structure. Incorporation of side-chain flexibility within the binding site resulted in a protein structure that stereoselectively binds to the D-enantiomer of the model ligand. In addition to four hydrogen bonds that are formed between amino acid residues in the binding site and the ligand, a number of hydrophobic interactions are involved in the formation of the antibody-ligand complex. The aromatic side chain of the ligand interacts with a tryptophan and a tyrosine residue in the binding site through pi-pi stacking. Fluorescence spectroscopic investigations also suggest the presence of tryptophan residues in the binding site, as ligand binding causes an enhancement of the antibody's intrinsic fluorescence at an emission wavelength of 350 nm. Based on the modeled antibody structure, the L-enantiomer of the model ligand cannot access the binding site due to steric hindrance. Additional docking experiments performed with D-phenylalanine and D-norvaline showed that these ligands are bound to the antibody in a way analogous to the D-enantiomer of the model ligand.     

Phenylalanine-Binding RNAs and Genetic Code Evolution

We isolated RNAs by selection–amplification, selecting for affinity to Phe–Sepharose and elution with free l-phenylalanine. Constant sequences did not contain Phe condons or anticodons, to avoid any possible confounding influence on initially randomized sequences. We examined the eight most frequent Phe-binding RNAs for inclusion of coding triplets. Binding sites were defined by nucleotide conservation, protection, and interference data. Together these RNAs comprise 70% of the 105 sequenced RNAs. The K _D for the strongest sites is ≈50 μM free amino acid, with strong stereoselectivity. One site strongly distinguishes free Phe from Trp and Tyr, a specificity not observed previously. In these eight Phe-binding RNAs, Phe codons are not significantly associated with Phe binding sites. However, among 21 characterized RNAs binding Phe, Tyr, Arg, and Ile, containing 1342 total nucleotides, codons are 2.7-fold more frequent within binding sites than in surrounding sequences in the same molecules. If triplets were not specifically related to binding sites, the probability of this distribution would be 4.8 × 10⁻¹¹. Therefore, triplet concentration within amino acid binding sites taken together is highly likely. In binding sites for Arg, Tyr, and Ile cognate codons are overrepresented. Thus Arg, Tyr, and Ile may be amino acids whose codons were assigned during an era of direct RNA–amino acid affinity. In contrast, Phe codons arguably were assigned by another criterion, perhaps during later code evolution.     

Identification and characterization of a novel legume-like lectin cDNA sequence from the red marine algae <Emphasis Type="Italic">Gracilaria fisheri</Emphasis>

A legume-type lectin (L-lectin) gene of the red algae Gracilaria fisheri (GFL) was cloned by rapid amplification of cDNA ends (RACE). The full-length cDNA of GFL was 1714 bp and contained a 1542 bp open reading frame encoding 513 amino acids with a predicted molecular mass of 56.5 kDa. Analysis of the putative amino acid sequence with NCBI-BLAST revealed a high homology (30–68%) with legume-type lectins (L-lectin) from Griffithsia japonica, Clavispora lusitaniae, Acyrthosiphon pisum, Tetraodon nigroviridis and Xenopus tropicalis. Phylogenetic relationship analysis showed the highest sequence identity to a glycoprotein of the red algae Griffithsia japonica (68%) (GenBank number AAM93989). Conserved Domain Database analysis detected an N-terminal carbohydrate recognition domain (CRD), the characteristic of L-lectins, which contained two sugar binding sites and a metal binding site. The secondary structure prediction of GFL showed a β-sheet structure, connected with turn and coil. The most abundant structural element of GFL was the random coil, while the α-helixes were distributed at the N- and C-termini, and 21 β-sheets were distributed in the CRD. Computer analysis of three-dimensional structure showed a common feature of L-lectins of GFL, which included an overall globular shape that was composed of a β-sandwich of two anti-parallel β-sheets, monosaccharide binding sites, were on the top of the structure and in proximity with a metal binding site. Northern blot analysis using a DIG-labelled probe derived from a partial GFL sequence revealed a hybridization signal of ~1.7 kb consistent with the length of the full-length GFL cDNA identified by RACE. No detectable band was observed from control total RNA extracted from filamentous green algae.     

3D modeling and molecular dynamics simulation of an immune-regulatory cytokine,interleukin-10, from the Indian major carp, <Emphasis Type="Italic">Catla catla</Emphasis>

Interleukin-10 (IL-10) is a pleiotropic immune-regulatory cytokine that is expressed in various species of fish and higher vertebrates, and is activated during infection. In spite of its important role, IL-10 has not been well characterized either functionally or structurally in fish. To analyze its properties and function, we constructed a 3D model of IL-10 in the Indian major carp, the catla (Catla catla), which is a highly preferred fish species and the most commercially important one in the Indian subcontinent. The catla IL-10 model was constructed by comparative modeling using human IL-10 (2ILK) as the template, and a 5 ns molecular dynamics (MD) simulation was carried out to characterize its structural and dynamical features, which was validated by the SAVES, WHAT IF and MolProbity servers. Analysis using the VAST server revealed a comparatively low level of homology between catla and human IL-10 amino acids at the N-terminal (22.7%) compared to the C-terminal (38.29%). Six conserved domains (A–F) were predicted in catla that threaded well with human IL-10, but their putative interaction sites varied significantly. The amino acid residues in helices A and F differed in length between catla and human IL-10, which may lead to the differences in the IL-10/IL-10R complexes of these two species. The existence of two highly conserved amino acid residues (Cys5 and Cys10) in fish IL-10 but not in higher vertebrate (including human) IL-10 was analyzed in this 3D model. CastP, cons-PPISP and InterProSurf server identified several binding pockets with various probe radii, but Cys5 and Cys10 did not form any significant bonds relating to structural stabilization or protein–protein interactions.