Structural basis for the unusual specificity of Escherichia coli aminopeptidase N

Aminopeptidase N from Escherichia coli is a M1 class aminopeptidase with the active-site region related to that of thermolysin. The enzyme has unusual specificity, cleaving adjacent to the large, nonpolar amino acids Phe and Tyr but also cleaving next to the polar residues Lys and Arg. To try to understand the structural basis for this pattern of hydrolysis, the structure of the enzyme was determined in complex with the amino acids L-arginine, L-lysine, L-phenylalanine, L-tryptophan, and L-tyrosine. These amino acids all bind with their backbone atoms close to the active-site zinc ion and their side chain occupying the S1 subsite. This subsite is in the form of a cylinder, about 10 A in cross-section and 12 A in length. The bottom of the cylinder includes the zinc ion and a number of polar side chains that make multiple hydrogen-bonding and other interactions with the alpha-amino group and the alpha-carboxylate of the bound amino acid. The walls of the S1 cylinder are hydrophobic and accommodate the nonpolar or largely nonpolar side chains of Phe and Tyr. The top of the cylinder is polar in character and includes bound water molecules. The epsilon-amino group of the bound lysine side chain and the guanidinium group of arginine both make multiple hydrogen bonds to this part of the S1 site. At the same time, the hydrocarbon part of the lysine and arginine side chains is accommodated within the nonpolar walls of the S1 cylinder. This combination of hydrophobic and hydrophilic binding surfaces explains the ability of ePepN to cleave Lys, Arg, Phe, and Tyr. Another favored substrate has Ala at the P1 position. The short, nonpolar side chain of this residue can clearly be bound within the hydrophobic part of the S1 cylinder, but the reason for its facile hydrolysis remains uncertain.     

Amino acid type-selective backbone 1H-15N-correlations for Arg and Lys

Four novel amino acid type-selective triple resonance experiments to identify the backbone amino proton and nitrogen resonances of Arg and Lys and of their sequential neighbors in (¹³C,¹⁵N)-labeled proteins are presented: the R(i+1)-HSQC and R(i,i+1)-HSQC select signals originating from Arg side chains, the K(i+1)-HSQC and K(i,i+1)-HSQC select signals originating from Lys side chains. The selection is based on exploiting the characteristic chemical shifts of a pair of carbon atoms in Arg and Lys side chains using selective 90° pulses. The new experiments are recorded as two-dimensional ¹H-¹⁵N-correlations and their performance is demonstrated with the application to a protein domain of 83 amino acids.     

Investigating the cationic side chains of the antimicrobial peptide tritrpticin: hydrogen bonding properties govern its membrane-disruptive activities

The positively charged side chains of cationic antimicrobial peptides are generally thought to provide the initial long-range electrostatic attractive forces that guide them towards the negatively charged bacterial membranes. Peptide analogs were designed to examine the role of the four Arg side chains in the cathelicidin peptide tritrpticin (VRRFPWWWPFLRR). The analogs include several noncoded Arg and Lys derivatives that offer small variations in side chain length and methylation state. The peptides were tested for bactericidal and hemolytic activities, and their membrane insertion and permeabilization properties were characterized by leakage assays and fluorescence spectroscopy. A net charge of +5 for most of the analogs maintains their high antimicrobial activity and directs them towards preferential insertion into model bacterial membrane systems with a similar extent of burial of the Trp side chains. However the peptides exhibit significant functional differences. Analogs with methylated cationic side chains cause lower levels of membrane leakage and are associated with lower hemolytic activities, making them potentially attractive pharmaceutical candidates. Analogs containing the Arg guanidinium groups cause more membrane disruption than those containing the Lys amino groups. Peptides in the latter group with shorter side chains have increased membrane activity and conversely, elongating the Arg residue causes slightly higher membrane activity. Altogether, the potential for strong hydrogen bonding between the four positive Arg side chains with the phospholipid head groups seems to be a determinant for the membrane disruptive properties of tritrpticin and many related cationic antimicrobial peptides.     

Key role of hydrazine to the interaction between oxaloacetic against phosphoenolpyruvic carboxykinase (PEPCK): ONIOM calculations

The interactions between oxaloacetic (OAA) and phosphoenolpyruvic carboxykinase (PEPCK) binding pocket in the presence and absence of hydrazine were carried out using quantum chemical calculations, based on the two-layered ONIOM (ONIOM2) approach. The complexes were partially optimized by ONIOM2 (B3LYP/6-31G(d):PM6) method while the interaction energies between OAA and individual residues surrounding the pocket were performed at the MP2/6-31G(d,p) level of theory. The calculated interaction energies (INT) indicated that Arg87, Gly237, Ser286, and Arg405 are key residues for binding to OAA with the INT values of ?1.93, ?2.06, ?2.47, and ?3.16 kcal mol^?1, respectively. The interactions are mainly due to the formation of hydrogen bonding interactions with OAA. Moreover, using ONIOM2 (B3LYP/6-31G(d):PM6) applied on the PEPCKHS complex, two proton transfers were observed; first, the proton was transferred from the carboxylic group of OAA to hydrazine while the second one was from Asp311 to Lys244. Such reactions cause the generation of binding strength of OAA to the pocket via electrostatic interaction. The orientations of Lys243, Lys244, His264, Asp311, Phe333, and Arg405 were greatly deviated after hydrazine incorporation. These indicate that hydrazine plays an important role in terms of not only changing the conformation of the binding pocket, but is also tightly bound to OAA resulting in its conformation change in the pocket. The understanding of such interaction can be useful for the design of hydrazine-based inhibitor for antichachexia agents.

Geometry of interaction of the histidine ring with other planar and basic residues

Among the aromatic residues in protein structures, histidine (His) is unique, as it can exist in the neutral or positively charged form at the physiological pH. As such, it can interact with other aromatic residues as well as form hydrogen bonds with polar and charged (both negative and positive) residues. We have analyzed the geometry of interaction of His residues with nine other planar side chains containing aromatic (residues Phe, Tyr, Trp, and His), carboxylate (Asp and Glu), carboxamide (Asn and Gln) and guanidinium (Arg) groups in 432 polypeptide chains. With the exception of the aspartic (Asp) and glutamic (Glu) acid side-chains, all other residues prefer to interact in a face-to-face or offset-face-stacked orientation with the His ring. Such a geometry is different from the edge-to-face relative orientation normally associated with the aromatic-aromatic interaction. His-His pair prefers to interact in a face-to-face orientation; however, when both the residues bind the same metal ion, the interplanar angle is close to 90 degrees. The occurrence of different interactions (including the nonconventional N-H...pi and C-H...pi hydrogen bonds) have been correlated with the relative orientations between the interacting residues. Several structural motifs, mostly involved in binding metal ions, have been identified by considering the cases where His residues are in contact with four other planar moieties. About 10% of His residues used here are also found in sequence patterns in PROSITE database. There are examples of the amino end of the Lys side chain interacting with His residues in such a way that it is located on an arc around a ring nitrogen atom.     

Identification of a molecular switch that selects between two crystals forms of bovine pancreatic trypsin inhibitor.

Two crystals forms of bovine pancreatic trypsin inhibitor are produced between pH 8.39 and 10.13 when crystals are grown at room temperature from solutions of 1.5 M potassium phosphate. Lower pH values favor the form II crystals, whereas higher pH values favor the form III. The transition from one crystal form to the other occurs at pH 9.35. We examined the crystal lattice contacts in both crystal forms and identified an unusual interaction we believe explains these observations. Spanning the crystallographic 2-fold axis in form III crystals, the Lys 41 side-chain amino nitrogens from 2 symmetry-related molecules are only 2.72 A apart, implying they are hydrogen bonded to one another. In form II crystals, the Lys 41 side-chain amino group is protonated and forms a salt bridge with a solvent-derived phosphate group. For the Lys 41 side-chain amino groups to hydrogen bond in form III crystals, at least 1 member of the pair must be deprotonated. The transition that occurs at pH 9.35 marks the pKa for deprotonation. In solution, the pKa for the Lys 41 side chain is around 10.8. The pKa for one of the interacting Lys 41 side chains in form III crystals is therefore shifted downward by about 1.5 pH units. The energy for lowering the pKa value comes from the many additional intermolecular hydrogen bonds that are present in form III crystals: 19 compared to only 8 in form II crystals.     

Structural and thermodynamic studies on a salt-bridge triad in the NADP-binding domain of glutamate dehydrogenase from Thermotoga maritima: cooperativity and electrostatic contribution to stability

Cooperative interactions within ion-pair networks of hyperthermostable proteins are thought to be a major determinant for extreme protein stability. While the favorable thermodynamic contributions of optimized electrostatics in general as well as those of pairwise interactions have been documented, cooperativity between pairwise interactions has not yet been studied thermodynamically in proteins from hyperthermophiles. In this study we use the isolated cofactor binding domain of glutamate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima to analyze pairwise and cooperative interactions within the salt-bridge triad Arg190-Glu231-Lys193. The X-ray structure of the domain was solved at 1.43 A and reveals the salt-bridge network with surrounding solvent molecules in detail. All three participating charges in the network were mutated to alanine in all combinations. The X-ray structure of the variant lacking all three charges reveals that the removal of the side chains has no effect on the overall conformation of the protein. Using solvent denaturation and thermodynamic cycles, the interaction energies between each pair of residues in the network were determined in the presence and in the absence of the third residue. Both the Arg190-Glu231 ion pair and the Lys193-Glu231 salt bridge in the absence of the third residue, contribute favorably to the free energy for unfolding of the domain in urea. Using guanidinium chloride as denaturant reveals a strong cooperativity between the two ion-pair interactions, the presence of the second ion pair converts the first interaction from destabilizing into stabilizing by as much as 1.09 kcal/mol. The different energetics of the salt-bridge triad in urea and GdmCl are discussed with reference to the observed anion binding in the crystal structure at high ionic strength and their possible role in a highly charged, high-temperature environment such as the cytoplasm of hyperthermophiles.     

Molecular mechanics simulation of the interaction of estrogen receptor with estrogen regulatory element.

A model for the interaction of 31 amino acid fragment (protein) from DNA binding domain of human estrogen receptor (hER) with a five base pair DNA sequence 5'GGTCA 3' from estrogen regulatory element (ERE) has been obtained using a step-wise procedure based on structural data on model peptides, DNA binding domain of hER, steric constrains imposed by tetrahedral coordination of the Cys sulphurs with zinc ion and classical secondary structural elements. Structure of the protein as well as its complex with DNA is obtained by energy minimization followed by refinement by molecular mechanics. The complex is stabilized by H-bonds between Lys22, Lys26 and Arg27 with DNA bases G2, T3 and T6. Lys22 also made H-bond with the backbone of G2. The backbone of Cys18 H-bonded with N7 of G1. DNA was in distorted B form and showed evidence of protein-induced conformational changes.     

Saccharomyces cerevisiae phosphoenolpyruvate carboxykinase: relevance of arginine 70 for catalysis

Saccharomyces cerevisiae phosphoenolpyruvate (PEP) carboxykinase is a key enzyme of the gluconeogenic pathway and catalyzes the decarboxylation of oxaloacetate and transfer of the gamma-phosphoryl group of ATP to yield PEP, ADP, and CO2 in the presence of a divalent metal ion. Previous experiments indicate that mutation of amino acid residues at metal site 1 decrease the enzyme catalytic efficiency and the affinity of the protein for PEP, evidencing the relevance of hydrogen-bond interactions between PEP and water molecules of the first coordination sphere of the metal ion for catalysis [Biochemistry 41 (2002) 12763]. To further understand the function of amino acid residues located in the PEP binding site, we have now addressed the catalytic importance of Arg70, whose guanidinium group is close to the PEP carboxyl group. Arg70 mutants of PEP carboxykinase were prepared, and almost unaltered kinetic parameters were found for the Arg70Lys PEP carboxykinase, while a decrease in 4-5 orders of magnitude for the catalytic efficiency was detected for the Arg70Gln and Arg70Met altered enzymes. To evaluate the enzyme interaction with PEP, the phosphopyridoxyl-derivatives of wild type, Arg70Lys, Arg70Gln, and Arg70Met S. cerevisiae PEP carboxykinase were prepared, and the change in the fluorescence emission of the probe upon PEP binding was used to obtain the dissociation equilibrium constant of the corresponding derivatized enzyme-PEP-Mn2+ complex. The titration experiments showed that a loss in 2.1 kcal/mol in PEP binding affinity is produced in the Arg70Met and Arg70Gln mutant enzymes. It is proposed that the electrostatic interaction between the guanidinium group of Arg70 and the carboxyl group of PEP is important for PEP binding and for further steps in catalysis.     

Refinement of the structure of human basic fibroblast growth factor at 1.6 A resolution and analysis of presumed heparin binding sites by selenate substitution.

The three-dimensional structure of human basic fibroblast growth factor has been refined to a crystallographic residual of 16.1% at 1.6 A resolution. The structure has a Kunitz-type fold and is composed of 12 antiparallel beta-strands, 6 of which form a beta-barrel. One bound sulfate ion has been identified in the model, hydrogen bonded to the side chains of Asn 27, Arg 120, and Lys 125. The side chain of Arg 120 has two conformations, both of which permit hydrogen bonds to the sulfate. This sulfate binding site has been suggested as the binding site for heparin (Eriksson, A.E., Cousens, L.S., Weaver, L.H., & Matthews, B.W., 1991, Proc. Natl. Acad. Sci. USA 88, 3441-3445). Two beta-mercaptoethanol (BME) molecules are also included in the model, each forming a disulfide bond to the S gamma atoms of Cys 69 and Cys 92, respectively. The side chain of Cys 92 has two conformations of which only one can bind BME. Therefore the BME molecule is half occupied at this site. The locations of possible sulfate binding sites on the protein were examined by replacing the ammonium sulfate in the crystallization medium with ammonium selenate. Diffraction data were measured to 2.2 A resolution and the structure refined to an R-factor of 13.8%. The binding of the more electron-dense selenate ion was identified at two positions. One position was identical to the sulfate binding site identified previously. The second selenate binding site, which is of lower occupancy, is situated 5.6 A from the first. This ion is hydrogen bonded by the side chain of Lys 135 and Arg 120. Thus the side chain of Arg 120 binds two selenate ions simultaneously. It is suggested that the observed second selenate binding site should also be considered as a possible binding site for heparin, or that both selenate binding sites might simultaneously contribute to the binding of heparin.     

ANS fluorescence: potential to augment the identification of the external binding sites of proteins

8-anilino-1-naphthalenesulfonic acid (ANS) is believed to strongly bind cationic groups of proteins and polyamino acids through ion pair formation. A paucity of data exists on the fluorescent properties of ANS in these interactions. ANS binding to arginine and lysine derivatives was studied by fluorescence and circular dichroism spectroscopies to augment published information attained by isothermal titration calorimetry (ITC). Fluorescence enhancement with a hypsochromic shift results from the interaction of the charged group of lysine and arginine with the sulfonate group of ANS. Ion pairing between Arg (or Lys) and the sulfonate group of ANS reduce the intermolecular charge transfer (CT) rate constant that leads to enhancement of fluorescence. A positive charge near the -NH group of ANS changes the intramolecular CT process producing a blue shift of fluorescence. The Arg side chain compared to that of Lys more effectively interacts with both the -NH and sulfonate groups of ANS. ANS binding also induces a random coil-alpha helix transition in poly-Arg. Our data, in contrast to ITC results, indicate that electrostatic interactions between ANS derivatives and positively charged side chains do not account for binding affinity in the micromolar range. In addition to ion pairing complementary interactions, such as van der Waals, should be considered for high affinity (K(d)<1 mM) external binding sites of proteins.     

Ab initio Models for Six-Center Multiple Proton Exchange and Ion Pair Formation Assisted by Lewis Acids

High level ab initio and density functional calculations, extrapolated to QCISD(T)/6-311+G(3df,2p)//MP2/6-31+G^**+ZPE, reveal that cyclic ion pairs can form in the hydrogen bonded complexes of haloboric acids BH_nX_3-n–HX, X=F, Cl, with Lewis bases HX, H₂O, CH₃OH, and NH₃, even in isolation (e.g., in the gas phase). The intrinsic acidities (deprotonation energies) required for protonation of these bases with formation of gas phase ion pairs are calculated to be <295 kcal/mol for water, <301 kcal/mol for methanol, and <306 kcal/mol for ammonia; such values are common for acidic sites in zeolites. All gas phase ion pairs prefer symmetric bidentate or tridentate structures. In the other cases where hydrogen bonded complexes prevail, symmetric ion pair-like transition structures for multiple hydrogen exchange are computed.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089400060563     

Effect of lysine side chain length on intra-helical glutamate--lysine ion pairing interactions

Ion-pairing interactions are important for protein stabilization. Despite the apparent electrostatic nature of these interactions, natural positively charged amino acids Lys and Arg have multiple methylenes linking the charged functionality to the backbone. Interestingly, the amino acids Lys and Orn have positively charged side chains that differ by only one methylene. However, only Lys is encoded and incorporated into proteins. To investigate the effect of side chain length of Lys on ion-pairing interactions, a series of 12 monomeric alpha-helical peptides containing potential Glu-Xaa (i, i+3), (i, i+4) and (i, i+5) (Xaa = Lys, Orn, Dab, Dap) interactions were studied by circular dichroism (CD) spectroscopy at pH 7 and 2. At pH 7, no Glu-Xaa (i, i+5) interaction was observed, regardless of the Xaa side chain length. Furthermore, only Lys was capable of supporting Glu-Xaa (i, i+3) interactions, whereas any Xaa side chain length supported Glu-Xaa (i, i+4) interactions. Side chain conformational analysis by molecular mechanics calculations showed that the side chain length of Lys enables the Glu-Xaa (i, i+3) interaction with lower energy conformations compared to residues with side chain lengths shorter than that of Lys. Furthermore, these calculated low energy conformers were consistent with conformations of intra-helical Glu-Lys salt bridges in a non-redundant protein structure database. Importantly, the CD spectra for peptides with Glu-Lys interactions did not alter significantly upon changing the pH because of a greater contribution to these interactions by forces other than electrostatics. Incorporating side chains just one methylene shorter (Orn) resulted in significant pH dependence or lack of interaction, suggesting that nature has chosen Lys to form durable interactions with negatively charged functional groups.     

Arginyl-tRNA synthetase from yeast. Discrimination between 20 amino acids in aminoacylation of tRNA(Arg)-C-C-A and tRNA(Arg)-C-C-A(3'NH2)

For discrimination between arginine and 19 other amino acids in aminoacylation of tRNA(Arg)-C-C-A by arginyl-tRNA synthetase from baker's yeast, discrimination factors (D) have been determined from kcat and Km values. The lowest values were found for Trp, Cys, Lys (D = 800-8500), showing that arginine is 800-8500 times more often incorporated into tRNA(Arg)-C-C-A than noncognate acids at the same amino acid concentrations. The other noncognate amino acids exhibit D values between 10,000 and 60,000. In aminoacylation of tRNA(Arg)-C-C-A(3'NH2) discrimination factors D1 are in the range 10-600. From these values and AMP formation stoichiometry, pretransfer proof-reading factors II1 were determined; from D values and AMP stoichiometry in aminoacylation of tRNA(Arg)-C-C-A, posttransfer proof-reading factors II2 could be calculated, II1 values between 2 and 120 show that pretransfer proof-reading is the main correction step, posttransfer proof-reading (II2 approximately 1-10) plays a marginal role. Initial discrimination factors due to different Gibbs free energies of binding between arginine and the noncognate amino acids were calculated from discrimination and proof-reading factors. According to a two-step binding process, two factors (I1 and I2) were determined. They can be related to hydrophobic interaction forces and hydrogen bonds that are especially formed by the arginine side chain. A hypothetical 'stopper' model of the amino acid recognition site is discussed.     

Effects of Ceramics Irradiation on Dissociation State of L-Amino Acids

Non-thermal effects of ceramics irradiation on dissociation state of twenty L-amino acids have been investigated. Dissociation constants of the amino acids other than His and Glu varied by a 3-h irradiation under cooling. pK’s of α-carboxyl group of amino acids having longer side chains on the α-carbon were decreased by the irradiation. Although pK’s of α-amino group of Arg, Lys, Asp, and CySH were decreased by the irradiation and pK of Tyr was increased, pK’s of the other fifteen amino acids were not affected. Although the isoelectric points of Lys, Arg, Trp, Asp, and CySH were decreased by the irradiation, those of the other fifteen amino acids were not affected. It was suggested that various changes in pH of amino acids in aqueous solution and dissociation state of the functional groups will be caused from stimulation by the irradiation and stabilization of the hydration layer around amino acids.     

Relative reactivity of lysine and other peptide-bound amino acids to oxidation by hypochlorite

Antibacterial and inflammatory responses of neutrophils and macrophages produce hypochlorite as a major oxidant. Numerous side chains of amino acids found in extracellular proteins can be modified by hypochlorite, including His, Arg, Tyr, Lys, Trp, and Met. We studied the relative reactivity of each of these amino acid residues in short N-blocked peptides, where other residues in the peptide were highly resistant to hypochlorite attack. Hypochlorite treatment led to modified peptides in each case, which were detected by changes in retention on reversed-phase HPLC. A distinct single product, consuming two equivalents of hypochlorite per equivalent of peptide, was obtained from the Lys-containing peptides. UV spectroscopy, nuclear magnetic resonance (NMR), and electrospray/mass spectroscopy identified this product as the dichloramine at the epsilon-amino group of the Lys side chain. The dichloramine at Lys did not decompose to form a detectable amount of carbonyl reactive with dinitrophenylhydrazine. The dichloramine at Lys did however quantitatively revert back to Lys during HCl digestion of the tetrapeptide for amino acid analysis, with simultaneous modification of the adjacent Phe residue. The formation of the dichloramine at Lys was not blocked by peptides or acetylated amino acids that contained Tyr, His, or Arg. In contrast, the presence of equimolar Met-containing peptide, or N-Acetyl-Trp, both inhibited the formation of the dichloramine at Lys. Thus, Met and Trp side chains of proteins might be able to protect Lys from chloramine formation under some circumstances, but this interpretation must consider that Met and Trp are typically found in relatively inaccessible hydrophobic sites, whereas lysine is typically exposed on the protein surface. The hierarchy of amino acid reactivities examined here will aid in the prediction of residues in biological samples most likely to be modified by hypochlorite.     

Stereochemical aspects of interaction of DNA binding domain of human progesterone receptor with d(AGGTCATGCT)2.

Structures of (i) 66 amino-acid fragment (residues 567-633) from DNA binding domain of human progesterone receptor (hPR), (ii) a ten base pair DNA sequence d(AGGTCATGCT)2 from hormone responsive element (HRE) and (iii) a complex of these two are optimised by computer modelling and molecular mechanics technique using extensive steric constraints from secondary structure predictions, comparison with the structures of known metalloproteins, geometric constraints imposed by tetrahedral coordination with the zinc ion and comparison with structures of DNA binding domains of human glucocorticoid and estrogen receptors (hGR and hER). Structure of the complex was obtained using genetic modification data on steroid receptors and general consensus about protein-DNA interaction. DNA is in distorted B conformation. Sequence dependent as well as protein-induced conformation changes are noticed. There is change in propeller twist, buckle and angle between glycosyl bonds. However, H-bonding network is preserved. The complex is stabilized with eighteen hydrogen-bonds, mainly between peptide side-chains and backbone phosphate. There are five specific H-bonds between basic amino acid side chains, Lys 22, Lys 26 and Arg 27, and DNA bases, A1, G3, G16 and A17. Gly 19, Ser 20 and Val 23 are in close proximity of DNA.     

Nonsequence-specific arginine interactions in the nucleosome core particle

We have analyzed the relative orientation of basic amino acid side chains towards DNA in the nucleosome core particle. The electric field created by DNA phosphates has no apparent preferential orientation: no favored orientation of the arginine guanidinium group is found. Arginine may be either directly hydrogen bonded to a phosphate oxygen or stabilized in the minor groove by van der Waals contacts and the local negative electric field. On the other hand, the phosphate oxygen atoms hydrogen bonded to arginines are always found close to the plane defined by the guanidinium group. Thus it can be concluded that the interactions of arginine are strongly directional, those of phosphate are not. We also find that a highly charged fragment of histone H2B, which is placed between two DNA turns, has a very variable conformation. An increase in protein positive charge density apparently allows multiple nonspecific protein conformations when interacting with DNA.     

Exploring the environmental preference of weak interactions in (alpha/beta)8 barrel proteins

The environmental preference for the occurrence of noncanonical hydrogen bonding and cation-pi interactions, in a data set containing 71 nonredundant (alpha/beta)(8) barrel proteins, with respect to amino acid type, secondary structure, solvent accessibility, and stabilizing residues has been performed. Our analysis reveals some important findings, which include (a) higher contribution of weak interactions mediated by main-chain atoms irrespective of the amino acids involved; (b) domination of the aromatic amino acids among interactions involving side-chain atoms; (c) involvement of strands as the principal secondary structural unit, accommodating cross strand ion pair interaction and clustering of aromatic amino acid residues; (d) significant contribution to weak interactions occur in the solvent exposed areas of the protein; (e) majority of the interactions involve long-range contacts; (f) the preference of Arg is higher than Lys to form cation-pi interaction; and (g) probability of theoretically predicted stabilizing amino acid residues involved in weak interaction is higher for polar amino acids such as Trp, Glu, and Gln. On the whole, the present study reveals that the weak interactions contribute to the global stability of (alpha/beta)(8) TIM-barrel proteins in an environment-specific manner, which can possibly be exploited for protein engineering applications.     

Theoretical studies on protein-nucleic acid interactions. II. Hydrogen bonding of amino acid side chains with bases and base pairs of nucleic acids

With a view to understanding the role of hydrogen bonds in the recognition of nucleic acids by proteins, hydrogen bonding between the bases and base pairs of nucleic acids and the amino acids (Asn, Gln, Asp and Glu, and charged residues Arg⁺, Glu^?, and Asp^?) has been studied by a second-order perturbation theory. Binding energies have been calculated for all possible configurations involving a pair of hydrogen bonds between the base (or base pair) and the amino acid residue. Our results show that the hydrogen bonding in these cases has a large contribution from electrostatic interaction. In general, the charged amino acids, compared to the uncharged ones, form more stable complexes with bases or base pairs. The hydrogen-bond energies are an order of magnitude smaller than the Coulombic interaction energies between basic amino acids (Lys⁺, Arg⁺, and His⁺) and the phosphate groups of nucleic acids. The stabilities of the complexes of amino acids Asn, Gln, Asp, and Glu with bases are in the order: G–X > C–X > A–X U–X or T–X, and G · C–X > A · T(U)–X, where X is one of these amino acid residues. It has been shown that Glu^? and Asp^? can recognize guanine in single-stranded nucleic acids; Arg⁺ can recognize G · C base pairs from A · T base pairs in double-stranded structures.