When fold is not important: a common structural framework for adenine and AMP binding in 12 unrelated protein families

ATP is a ligand common to many proteins, yet it is unclear whether common recognition patterns do exist among the many different folds that bind ATP. Previously, it was shown that cAMP-dependent protein kinase, D-Ala:D-Ala ligase and the alpha-subunit of the alpha 2 beta 2 ribonucleotide reductase do share extensive common structural elements for ATP recognition although their folds are different. Here, we have made a survey of structures that bind ATP and compared them with the key features seen in these three proteins. Our survey shows that 12 different fold types share a specific recognition pattern for the adenine moiety, and 8 of these folds have a common structural framework for recognition of the AMP moiety of the ligand. The common framework consists of a tripeptide segment plus three additional residues, which provides similar polar and hydrophobic interactions between the protein and mononucleotide. Consensus interactions are represented by four key hydrogen bonds present in each fold type. Two of these four hydrogen bonds, together with three aliphatic residues, form a specific recognition pattern for the adenine moiety in all 12 folds. These similarities point to a structural-functional requirement shared by these different mononucleotide-binding proteins that represent at this time 28% of the adenine mononucleotide complexes found in the Brookhaven Protein Data Bank.     

Diversity in hapten recognition: structural study of an anti-cocaine antibody M82G2

Antibodies against cocaine and other drugs of abuse are the basis for diagnostic tests for the presence of those drugs in human serum. The 1.7A resolution crystal structure of the anti-cocaine monoclonal antibody M82G2 in complex with cocaine is presented. This structure determination was undertaken to establish the stereochemical features in the antibody binding site that confer specificity for cocaine, and as part of an ongoing project to understand the rules that govern molecular recognition. The cocaine-binding site can be characterized topologically as a narrow groove on the protein surface. The antibody utilizes water-mediated hydrogen bonding, and cation-pi and stacking (pi-pi) interactions to provide specificity. Comparison with the previously published structure of the anti-cocaine antibody GNC92H2 shows that binding of a small ligand can be achieved in diverse ways, both in terms of a binding site structure/topology and protein-ligand interactions.     

Solution structure and base perturbation studies reveal a novel mode of alkylated base recognition by 3-methyladenine DNA glycosylase I

The specific recognition mechanisms of DNA repair glycosylases that remove cationic alkylpurine bases in DNA are not well understood partly due to the absence of structures of these enzymes with their cognate bases. Here we report the solution structure of 3-methyladenine DNA glycosylase I (TAG) in complex with its 3-methyladenine (3-MeA) cognate base, and we have used chemical perturbation of the base in combination with mutagenesis of the enzyme to evaluate the role of hydrogen bonding and pi-cation interactions in alkylated base recognition by this DNA repair enzyme. We find that TAG uses hydrogen bonding with heteroatoms on the base, van der Waals interactions with the 3-Me group, and conventional pi-pi stacking with a conserved Trp side chain to selectively bind neutral 3-MeA over the cationic form of the base. Discrimination against binding of the normal base adenine is derived from direct sensing of the 3-methyl group, leading to an induced-fit conformational change that engulfs the base in a box defined by five aromatic side chains. These findings indicate that base specific recognition by TAG does not involve strong pi-cation interactions, and suggest a novel mechanism for alkylated base recognition and removal.     

Binding of the neuroleptic drug haloperidol to a monoclonal antibody: refinement of the binding site molecular model using canonical structures

The ligand binding site of a monoclonal antibody (185), which binds the neuroleptic drug haloperidol, has been modelled using canonical structures and energy minimization techniques. This refined modelling protocol has allowed us to predict the variable region loop conformations. Three key residues, H:50(W), H:100a(D) and L:96(Y) appear to create the basis of the electrostatic, pi-pi stacking interactions and hydrogen bonding required for the high affinity binding site characteristics present in this antibody. The use of computer-aided graphics techniques and appropriate three-dimensional modelling permits inspection of the predicted molecular recognition features of the ligand binding site.     

Probing the “fingers” domain binding pocket of Hepatitis C virus NS5B RdRp and D559G resistance mutation via molecular docking,molecular dynamics simulation and binding free energy calculations

The NS5B RdRp polymerase is a prominent enzyme for the replication of Hepatitis C virus (HCV). During the HCV replication, the template RNA binding takes place in the “fingers” sub-domain of NS5B. The “fingers” domain is a new emerging allosteric site for the HCV drug development. The inhibitors of the “fingers” sub-domain adopt a new antiviral mechanism called RNA intervention. The details of essential amino acid residues, binding mode of the ligand, and the active site intermolecular interactions of RNA intervention reflect that this mechanism is ambiguous in the experimental study. To elucidate these details, we performed molecular docking analysis of the fingers domain inhibitor quercetagetin (QGN) with NS5B polymerase. The detailed analysis of QGN-NS5B intermolecular interactions was carried out and found that QGN interacts with the binding pocket amino acid residues Ala97, Ala140, Ile160, Phe162, Gly283, Gly557, and Asp559; and also forms π?π stacking interaction with Phe162 and hydrogen bonding interaction with Gly283. These are found to be the essential interactions for the RNA intervention mechanism. Among the strong hydrogen bonding interactions, the QGN?Ala140 is a newly identified important hydrogen bonding interaction by the present work and this interaction was not resolved by the previously reported crystal structure. Since D559G mutation at the fingers domain was reported for reducing the inhibition percentage of QGN to sevenfold, we carried out molecular dynamics (MD) simulation for wild and D559G mutated complexes to study the stability of protein conformation and intermolecular interactions. At the end of 50?ns MD simulation, the π?π stacking interaction of Phe162 with QGN found in the wild-type complex is altered into T-shaped π stacking interaction, which reduces the inhibition strength. The origin of the D559G resistance mutation was studied using combined MD simulation, binding free energy calculations and principal component analysis. The results were compared with the wild-type complex. The mutation D559G reduces the binding affinity of the QGN molecule to the fingers domain. The free energy decomposition analysis of each residue of wild-type and mutated complexes revealed that the loss of non-polar energy contribution is the origin of the resistance.

Communicated by Ramaswamy H. Sarma     

Structural requirements for the specific recognition of an m7G mRNA cap

7-Methylguanosine (m(7)G), also known as the mRNA "cap", is used as a molecular tag in eukaryotic cells to mark the 5' end of messenger RNAs. The mRNA cap is required for several key events in gene expression in which the m(7)G moiety is specifically recognized by cellular proteins. The configurations of the m(7)G-binding pockets of a cellular (eIF4E) and a viral (VP39) cap-binding protein have been determined by X-ray crystallography. The binding energy has been hypothesized to result from a pi-pi stacking interaction between aromatic residues sandwiching the m(7)G base in addition to hydrogen bonds between the base and acidic protein side chains. To further understand the structural requirements for the specific recognition of an m(7)G mRNA cap, we determined the effects of amino acid substitutions in eIF4E and VP39 cap-binding sites on their affinity for m(7)GDP. The requirements for residues suggested to pi-pi stack and hydrogen bond with the m(7)G base were examined in each protein by measuring their affinities for m(7)GDP by fluorimetry. The results suggest that both eIF4E and VP39 require a complicated pattern of both orientation and identity of the stacking aromatic residues to permit the selective binding of m(7)GDP.     

Probing the energetic and structural role of amino acid/nucleobase cation-pi interactions in protein-ligand complexes

X-ray structures of proteins bound to ligand molecules containing a nucleic acid base were systematically searched for cation-pi interactions between the base and a positively charged or partially charged side chain group located above it, using geometric criteria. Such interactions were found in 38% of the complexes and are thus even more frequent than pi-pi stacking interactions. They are moreover well conserved in families of related proteins. The overwhelming majority of cation-pi contacts involve Ade bases, as these constitute by far the most frequent ligand building block; Arg-Ade is the most frequent cation-pi pair. Ab initio energy calculations at MP2 level were performed on all recorded pairs. Though cation-pi interactions involving the net positive charge carried by Arg or Lys side chains are the most favorable energetically, those involving the partial positive charge of Asn and Gln side chain amino groups (sometimes referred to as amino-pi interactions) are favorable too, owing to the electron correlation energy contribution. Chains of cation-pi interactions with a nucleobase bound simultaneously to two charged groups or a charged group sandwiched between two aromatic moieties are found in several complexes. The systematic association of these motifs with specific ligand molecules in unrelated protein sequences raises the question of their role in protein-ligand structure, stability, and recognition.     

Complementarity of hydrophobic properties in ATP-protein binding: a new criterion to rank docking solutions

ATP is an important substrate of numerous biochemical reactions in living cells. Molecular recognition of this ligand by proteins is very important for understanding enzymatic mechanisms. Considerable insight into the problem may be gained via molecular docking simulations. At the same time, standard docking protocols are often insufficient to predict correct conformations for protein-ATP complexes. Thus, in most cases the native-like solutions can be found among the docking poses, but current scoring functions have only limited ability to discriminate them from false positives. To improve the selection of correct docking solutions obtained with the GOLD software, we developed a new ranking criterion specific for ATP-protein binding. The method is based on detailed analysis of the intermolecular interactions in 40 high-resolution 3D structures of ATP-protein complexes (the training set). We found that the most important factors governing this recognition are hydrogen-bonding, stacking between adenine and aromatic protein residues, and hydrophobic contacts between adenine and protein residues. To address the latter, we applied the formalism of 3D molecular hydrophobicity potential. The results obtained were used to construct an ATP-oriented scoring criterion as a linear combination of the terms describing these intermolecular interactions. The criterion was then validated using the test set of 10 additional ATP-protein complexes. As compared with the standard scoring functions, the new ranking criterion significantly improved the selection of correct docking solutions in both sets and allowed considerable enrichment at the top of the list containing docking poses with correct solutions.     

Insight into the structural role of carotenoids in the photosystem I: a quantum chemical analysis

The structural stabilization role of carotenoids in the formation of photosynthetic pigment-protein complexes is investigated theoretically. The pi-pi stacking and CH-pi interactions between beta-carotenes and their surrounding chlorophylls (and/or aromatic residues) in Photosystem I (PS1) from the cyanobacterium Synechococcus elongatus were studied by means of the supermolecular approach at the level of the second-order M?ller-Plesset perturbation method. PS1 features a core integral antenna system consisting of 22 beta-carotenes intertwined with 90 chlorophyll molecules. The binding environments of all 22 beta-carotenes were systematically analyzed. For 21 out of the 22 cases, one or more chlorophyll molecules exist within van der Waals' contacts of the beta-carotene molecule. The calculated strengths of pi-pi stacking interactions between the conjugated core of beta-carotene and the aromatic tetrapyrrole rings of chlorophyll are substantial, ranging from -3.54 kcal/mol for the perpendicular-positioned BCR4004...CHL1217 pair to -16.01 kcal/mol for the parallel-oriented BCR4007...CHL1122 pair. A strong dependence of the pi-pi stacking interaction energies on the intermolecular configurations of the two interacting pi-planes is observed. The parallel-oriented beta-carotene and chlorophyll pair is energetically much more stable than the perpendicular-positioned pair. The larger the extent of pi-pi overlapping, the stronger the interaction strength. In many cases, the beta-ring ends of beta-carotene molecules are found to interact with the tetrapyrrole rings of chlorophyll via CH-pi interactions. For the latter interactions, the calculated interaction strengths vary from -7.03 to -11.03 kcal/mol, depending on the intermolecular configuration. This work leads to the conclusion that pi-pi stacking and CH-pi interactions between beta-carotene and their surrounding chlorophylls and aromatic residues play an essential role in binding beta-carotenes in PS1 from S. elongatus. Consequently, the molecular basis of the structural stabilization function of carotenoids in formation of the photosynthetic pigment-protein complexes is established.     

Intermolecular interactions in staphylokinase-plasmin(ogen) bimolecular complex: function of His43 and Tyr44

Staphylokinase (SAK) forms a 1:1 stoichiometric complex with human plasmin (Pm) and switches its substrate specificity to generate a plasminogen (Pg) activator complex. Site-directed mutagenesis of SAKHis43 and SAKTyr44 demonstrated the crucial requirement of a positively charged and an aromatic residue, respectively, at these positions for optimal functioning of SAK-Pm activator complex. Molecular modeling studies further revealed the role of these residues in making cation-pi and pi-pi interactions with Trp215 of Pm and thus establishing the crucial intermolecular contacts within the active site cleft of the activator complex for the cofactor activity of SAK.     

Crystal packing driven by metal-ligand interactions, hydrogen bonds, π-π stacking and anion-π stacking

Two polymeric copper(II) compounds have been synthesized with the ligand bis(pyrimidin-2yl)amine (dipm), by means of coordination bonds, Watson-Crick type hydrogen bond interactions and π-π stacking. Both supramolecules hold the same cationic building block, namely [Cu(dipm)₂(H₂O)₂]²⁺. However, their crystal structure significantly differs and this variation apparently arises from the nature of the anion which induces dissimilar crystal packings. The crystal growth is driven by several synergistic intermolecular interactions, i.e., coordination and hydrogen bonds, π-π and anion-π stacks.     

Recognition templates for predicting adenylate-binding sites in proteins.

Recognition templates encapsulate the structural and energetic features for the specific recognition of a given ligand by a protein active site. These templates identify the major interactions used for specific recognition and may be used to find specific binding sites in proteins of unknown function. We present a grid-based method for deriving recognition templates for adenylate groups from a set of diverse nucleotide-binding proteins. The templates reveal the basis of specific binding of adenylate, including tight shape complementarity, specific hydrogen bonds, and underscoring the importance of a key steric contact for excluding guanylate from adenylate-specific sites. We demonstrate the utility of recognition templates in identifying specific adenylate-binding sites in a diverse set of dinucleotide-binding proteins.     

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.     

Structure-wise discrimination of adenine and guanine by proteins on the basis of their nonbonded interactions

We have analyzed the nonbonded interactions of the structurally similar moieties, adenine and guanine forming complexes with proteins. The results comprise (a) the amino acid–ligand atom preferences, (b) solvent accessibility of ligand atoms before and after complex formation with proteins, and (c) preferred amino acid residue atoms involved in the interactions. We have observed that the amino acid preferences involved in the hydrogen bonding interactions vary for adenine and guanine. The structural variation between the purine atoms is clearly reflected by their burial tendency in the solvent environment. Correlation of the mean amino acid preference values show the variation that exists between adenine and guanine preferences of all the amino acid residues. All our observations provide evidence for the discriminating nature of the proteins in recognizing adenine and guanine.     

Exploring the role of cation-pi interactions in glycoproteins lipid-binding proteins and RNA-binding proteins

We have analyzed and compared the influence of cation-pi interactions in glycoproteins (GPs), lipid-binding proteins (LBPs) and RNA-binding proteins (RBPs) in this study. We observed that all the proteins included in the study had profound cation-pi interactions. There is an average of one energetically significant cation-pi interaction for every 71 residues in GPs, for every 58 residues in LBPs and for every 64 residues in RBPs. Long-range contacts are predominant in all the three types of proteins studied. The pair-wise cation-pi interaction energy between the positively charged and aromatic residues shows that Arg-Trp pair energy was the strongest among all six possible pairs in all the three types of proteins studied. There were considerable differences in the preference of cation-pi interacting residues to different secondary structure elements and ASA and these might contribute to differences in biochemical functions of GPs, LBPs and RBPs. It was interesting to note that all the five residues involved in cation-pi interactions were found to have stabilization centers in GPs, LBPs and RBPs. Majority of the cation-pi interacting residues investigated in the present study had a conservation score of 6, the cutoff value used to identify the stabilizing residues. A small percentage of cation-pi interacting residues were also present as stabilizing residues. The cation-pi interaction-forming residues play an important role in the structural stability of in GPs, LBPs and RBPs. The results obtained in this study will be helpful in further understanding the stability, specificity and differences in the biochemical functions of GPs, LBPs and RBPs.     

Analysis of anisotropic side-chain packing in proteins and application to high-resolution structure prediction

pi-pi, Cation-pi, and hydrophobic packing interactions contribute specificity to protein folding and stability to the native state. As a step towards developing improved models of these interactions in proteins, we compare the side-chain packing arrangements in native proteins to those found in compact decoys produced by the Rosetta de novo structure prediction method. We find enrichments in the native distributions for T-shaped and parallel offset arrangements of aromatic residue pairs, in parallel stacked arrangements of cation-aromatic pairs, in parallel stacked pairs involving proline residues, and in parallel offset arrangements for aliphatic residue pairs. We then investigate the extent to which the distinctive features of native packing can be explained using Lennard-Jones and electrostatics models. Finally, we derive orientation-dependent pi-pi, cation-pi and hydrophobic interaction potentials based on the differences between the native and compact decoy distributions and investigate their efficacy for high-resolution protein structure prediction. Surprisingly, the orientation-dependent potential derived from the packing arrangements of aliphatic side-chain pairs distinguishes the native structure from compact decoys better than the orientation-dependent potentials describing pi-pi and cation-pi interactions.     

Adenine recognition: a motif present in ATP-, CoA-, NAD-, NADP-, and FAD-dependent proteins.

Adenosine triphosphate (ATP) plays an essential role in energy transfer within the cell. In the form of NAD, adenine participates in multiple redox reactions. Phosphorylation and ATP-hydrolysis reactions have key roles in signal transduction and regulation of many proteins, especially enzymes. In each cell, proteins with many different functions use adenine and its derivatives as ligands; adenine, of course, is present in DNA and RNA. We show that an adenine binding motif, which differs according to the backbone chain direction of a loop that binds adenine (and in one variant by the participation of an aspartate side-chain), is common to many proteins; it was found from an analysis of all adenylate-containing protein structures from the Protein Data Bank. Indeed, 224 protein-ligand complexes (86 different proteins) from a total of 645 protein structure files bind ATP, CoA, NAD, NADP, FAD, or other adenine-containing ligands, and use the same structural elements to recognize adenine, regardless of whether the ligand is a coenzyme, cofactor, substrate, or an allosteric effector. The common adenine-binding motif shown in this study is simple to construct. It uses only (1) backbone polar interactions that are not dependent on the protein sequence or particular properties of amino acid side-chains, and (2) nonspecific hydrophobic interactions. This is probably why so many different proteins with different functions use this motif to bind an adenylate-containing ligand. The adenylate-binding motif reported is present in "ancient proteins" common to all living organisms, suggesting that adenine-containing ligands and the common motif for binding them were exploited very early in evolution. The geometry of adenine binding by this motif mimics almost exactly the geometry of adenine base-pairing seen in DNA and RNA.     

NADP-Dependent enzymes. I: Conserved stereochemistry of cofactor binding

The ubiquitous redox cofactors nicotinamide adenine dinucleotides [NAD and NADP] are very similar molecules, despite their participation in substantially different biochemical processes. NADP differs from NAD in only the presence of an additional phosphate group esterified to the 2′-hydroxyl group of the ribose at the adenine end and yet NADP is confined with few exceptions to the reactions of reductive biosynthesis, whereas NAD is used almost exclusively in oxidative degradations. The discrimination between NAD and NADP is therefore an impressive example of the power of molecular recognition by proteins. The many known tertiary structures of NADP complexes affords the possibility for an analysis of their discrimination. A systematic analysis of several crystal structures of NAD(P)-protein complexes show that: 1) the NADP coenzymes are more flexible in conformation than those of NAD; 2) although the protein-cofactor interactions are largely conserved in the NAD complexes, they are quite variable in those of NADP; and 3) in both cases the pocket around the nicotinamide moiety is substrate dependent. The conserved and variable interactions between protein and cofactors in the respective binding pockets are reported in detail. Discrimination between NAD and NADP is essentially a consequence of the overall pocket and not of a few residues. A clear fingerprint in NAD complexes is a carboxylate side chain that chelates the diol group at the ribose near the adenine, whereas in NADP complexes an arginine side chain faces the adenine plane and interacts with the phosphomonoester. The latter type of interaction might be a general feature of recognition of nucleotides by proteins. Other features such as strand-like hydrogen bonding between the NADP diphosphate moeties and the protein are also significant. The NADP binding pocket properties should prove useful in protein engineering and design. © 1997 Wiley-Liss Inc.     

Structural analysis of cation-pi interactions in DNA binding proteins

Cation-pi interactions play an important role in the stability of protein structures. In this work, we have analyzed the influence of cation-pi interactions in DNA binding proteins. We observed cation-pi interactions in 45 out of 62 DNA binding proteins and there is no significant correlation between the number of amino acid residues and number of cation-pi interactions. These interactions are mainly formed by long-range contacts, and the role of short and medium-range contacts is minimal. The preference of Arg is higher than Lys to form cation-pi interactions. The pair-wise cation-pi interaction energy between aromatic and positively charged residues shows that Arg-Tyr energy is the strongest among the possible six pairs. The structural analysis of cation-pi interaction forming residues shows that Lys, Trp, and Tyr prefer to be in the binding site of protein-DNA complexes. Further, the accessible surface areas of cation-pi interaction forming cationic residues are significantly less than that of other residues. The preference of cation-pi interaction forming residues in different secondary structures shows that Lys prefers to be in strand and Phe prefers to be in turn regions. The results obtained in the present study will be useful in understanding the contribution of cation-pi interactions to the stability and specificity of protein-DNA complexes.