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.     

Interactions of Streptomyces griseus aminopeptidase with amino acid reaction products and their implications toward a catalytic mechanism

Streptomyces griseus aminopeptidase (SGAP) is a double-zinc exopeptidase with a high preference toward large hydrophobic amino-terminus residues. It is a monomer of a relatively low molecular weight (30 kDa), it is heat stable, it displays a high and efficient catalytic turnover, and its activity is modulated by calcium ions. The small size, high activity, and heat stability make SGAP a very attractive enzyme for various biotechnological applications, among which is the processing of recombinant DNA proteins and fusion protein products. Several free amino acids, such as phenylalanine, leucine, and methionine, were found to act as weak inhibitors of SGAP and hence were chosen for structural studies. These inhibitors can potentially be regarded as product analogs because one of the products obtained in a normal enzymatic reaction is the cleaved amino terminal amino acid of the substrate. The current study includes the X-ray crystallographic analysis of the SGAP complexes with methionine (1.53 A resolution), leucine (1.70 A resolution), and phenylalanine (1.80 A resolution). These three high-resolution structures have been used to fully characterize the SGAP active site and to identify some of the functional groups of the enzyme that are involved in enzyme-substrate and enzyme-product interactions. A unique binding site for the terminal amine group of the substrate (including the side chains of Glu131 and Asp160, as well as the carbonyl group of Arg202) is indicated to play an important role in the binding and orientation of both the substrate and the product of the catalytic reaction. These studies also suggest that Glu131 and Tyr246 are directly involved in the catalytic mechanism of the enzyme. Both of these residues seem to be important for substrate binding and orientation, as well as the stabilization of the tetrahedral transition state of the enzyme-substrate complex. Glu131 is specifically suggested to function as a general base during catalysis by promoting the nucleophilic attack of the zinc-bound water/hydroxide on the substrate carbonyl carbon. The structures of the three SGAP complexes are compared with recent structures of three related aminopeptidases: Aeromonas proteolytica aminopeptidase (AAP), leucine aminopeptidase (LAP), and methionine aminopeptidase (MAP) and their complexes with corresponding inhibitors and analogs. These structural results have been used for the simulation of several species along the reaction coordinate and for the suggestion of a general scheme for the proteolytic reaction catalyzed by SGAP.     

Influence of rotational energy barriers to the conformational search of protein loops in molecular dynamics and ranking the conformations.

An adjustable-barrier dihedral angle potential was added as an extension to a novel, previously presented soft-core potential to study its contribution to the efficacy of the search of the conformational space in molecular dynamics. As opposed to the conventional soft-core potential functions, the leading principle in the design of the new soft-core potential, as well as of its extension, the soft-core and adjustable-barrier dihedral angle (SCADA) potential (referred as the SCADA potential), was to maintain the main equilibrium properties of the original force field. This qualifies the methods for a variety of a priori modeling problems without need for additional restraints typically required with the conventional soft-core potentials. In the present study, the different potential energy functions are applied to the problem of predicting loop conformations in proteins. Comparison of the performance of the soft-core and SCADA potential showed that the main hurdles for the efficient sampling of the conformational space of (loops in) proteins are related to the high-energy barriers caused by the Lennard-Jones and Coulombic energy terms, and not to the rotational barriers, although the conformational search can be further enhanced by lowering the rotational barriers of the dihedral angles. Finally, different evaluation methods were studied and a few promising criteria found to distinguish the near-native loop conformations from the wrong ones.     

Substrate flow in catalases deduced from the crystal structures of active site variants of HPII from Escherichia coli.

The active site of heme catalases is buried deep inside a structurally highly conserved homotetramer. Channels leading to the active site have been identified as potential routes for substrate flow and product release, although evidence in support of this model is limited. To investigate further the role of protein structure and molecular channels in catalysis, the crystal structures of four active site variants of catalase HPII from Escherichia coli (His128Ala, His128Asn, Asn201Ala, and Asn201His) have been determined at approximately 2.0-A resolution. The solvent organization shows major rearrangements with respect to native HPII, not only in the vicinity of the replaced residues but also in the main molecular channel leading to the heme distal pocket. In the two inactive His128 variants, continuous chains of hydrogen bonded water molecules extend from the molecular surface to the heme distal pocket filling the main channel. The differences in continuity of solvent molecules between the native and variant structures illustrate how sensitive the solvent matrix is to subtle changes in structure. It is hypothesized that the slightly larger H(2)O(2) passing through the channel of the native enzyme will promote the formation of a continuous chain of solvent and peroxide. The structure of the His128Asn variant complexed with hydrogen peroxide has also been determined at 2.3-A resolution, revealing the existence of hydrogen peroxide binding sites both in the heme distal pocket and in the main channel. Unexpectedly, the largest changes in protein structure resulting from peroxide binding are clustered on the heme proximal side and mainly involve residues in only two subunits, leading to a departure from the 222-point group symmetry of the native enzyme. An active role for channels in the selective flow of substrates through the catalase molecule is proposed as an integral feature of the catalytic mechanism. The Asn201His variant of HPII was found to contain unoxidized heme b in combination with the proximal side His-Tyr bond suggesting that the mechanistic pathways of the two reactions can be uncoupled.     

Model for the three-dimensional structure of vitronectin: predictions for the multi-domain protein from threading and docking.

The structure of vitronectin, an adhesive protein that circulates in high concentrations in human plasma, was predicted through a combination of computational methods and experimental approaches. Fold recognition and sequence-structure alignment were performed using the threading program PROSPECT for each of three structural domains, i.e., the N-terminal somatomedin B domain (residues 1-53), the central region that folds into a four-bladed beta-propeller domain (residues 131-342), and the C-terminal heparin-binding domain (residues 347-459). The atomic structure of each domain was generated using MODELLER, based on the alignment obtained from threading. Docking experiments between the central and C-terminal domains were conducted using the program GRAMM, with limits on the degrees of freedom from a known inter-domain disulfide bridge. The docked structure has a large inter-domain contact surface and defines a putative heparin-binding groove at the inter-domain interface. We also docked heparin together with the combined structure of the central and C-terminal domains, using GRAMM. The predictions from the threading and docking experiments are consistent with experimental data on purified plasma vitronectin pertaining to protease sensitivity, ligand-binding sites, and buried cysteines.     

Model of three-dimensional structure of vitamin D receptor and its binding mechanism with 1alpha,25-dihydroxyvitamin D(3).

Comparative modeling of the vitamin D receptor three-dimensional structure and computational docking of 1alpha,25-dihydroxyvitamin D(3) into the putative binding pocket of the two deletion mutant receptors: (207-423) and (120-422, Delta [164-207]) are reported and evaluated in the context of extensive mutagenic analysis and crystal structure of holo hVDR deletion protein published recently. The obtained molecular model agrees well with the experimentally determined structure. Six different conformers of 1alpha,25-dihydroxyvitamin D(3) were used to study flexible docking to the receptor. On the basis of values of conformational energy of various complexes and their consistency with functional activity, it appears that 1alpha,25-dihydroxyvitamin D(3) binds the receptor in its 6-s-trans form. The two lowest energy complexes obtained from docking the hormone into the deletion protein (207-423) differ in conformation of ring A and orientation of the ligand molecule in the VDR pocket. 1alpha,25-Dihydroxyvitamin D(3) possessing the A-ring conformation with axially oriented 1alpha-hydroxy group binds receptor with its 25-hydroxy substituent oriented toward the center of the receptor cavity, whereas ligand possessing equatorial conformation of 1alpha-hydroxy enters the pocket with A ring directed inward. The latter conformation and orientation of the ligand is consistent with the crystal structure of hVDR deletion mutant (118-425, Delta [165-215]). The lattice model of rVDR (120-422, Delta [164-207]) shows excellent agreement with the crystal structure of the hVDR mutant. The complex obtained from docking the hormone into the receptor has lower energy than complexes for which homology modeling was used. Thus, a simple model of vitamin D receptor with the first two helices deleted can be potentially useful for designing a general structure of ligand, whereas the advanced lattice model is suitable for examining binding sites in the pocket.     

Structure of a legume lectin from the bark of Robinia pseudoacacia and its complex with N-acetylgalactosamine

The structure of the bark lectin RPbAI (isoform A4) from Robinia pseudoacacia has been determined by protein crystallography both in the free form and complexed with N-acetylgalactosamine. The free form is refined at 1.80 A resolution to an R-factor of 18.9% whereas the complexed structure has an R-factor of 19.7% at 2.05 A resolution. Both structures are compared to each other and to other available legume lectin structures. The polypeptide chains of the two structures exhibit the characteristic legume lectin tertiary fold. The quaternary structure resembles that of the Phaseolus vulgaris lectin, the soybean agglutinin, and the Dolichos biflorus lectin, but displays some unique features leading to the extreme stability of this lectin.     

Structural implications of drug-resistant mutants of HIV-1 protease: high-resolution crystal structures of the mutant protease/substrate analogue complexes

Emergence of drug-resistant mutants of HIV-1 protease is an ongoing problem in the fight against AIDS. The mechanisms governing resistance are both complex and varied. We have determined crystal structures of HIV-1 protease mutants, D30N, K45I, N88D, and L90M complexed with peptide inhibitor analogues of CA-p2 and p2-NC cleavage sites in the Gag-pol precursor in order to study the structural mechanisms underlying resistance. The structures were determined at 1.55-1.9-A resolution and compared with the wild-type structure. The conformational disorder seen for most of the hydrophobic side-chains around the inhibitor binding site indicates flexibility of binding. Eight water molecules are conserved in all 9 structures; their location suggests that they are important for catalysis as well as structural stability. Structural differences among the mutants were analyzed in relation to the observed changes in protease activity and stability. Mutant L90M shows steric contacts with the catalytic Asp25 that could destabilize the catalytic loop at the dimer interface, leading to its observed decreased dimer stability and activity. Mutant K45I reduces the mobility of the flap and the inhibitor and contributes to an enhancement in structural stability and activity. The side-chain variations at residue 30 relative to wild-type are the largest in D30N and the changes are consistent with the altered activity observed with peptide substrates. Polar interactions in D30N are maintained, in agreement with the observed urea sensitivity. The side-chains of D30N and N88D are linked through a water molecule suggesting correlated changes at the two sites, as seen with clinical inhibitors. Structural changes seen in N88D are small; however, water molecules that mediate interactions between Asn88 and Thr74/Thr31/Asp30 in other complexes are missing in N88D.     

Solution structure and dynamics of ribonuclease Sa.

We have used NMR methods to characterize the structure and dynamics of ribonuclease Sa in solution. The solution structure of RNase Sa was obtained using the distance constraints provided by 2,276 NOEs and the C6-C96 disulfide bond. The 40 resulting structures are well determined; their mean pairwise RMSD is 0.76 A (backbone) and 1.26 A (heavy atoms). The solution structures are similar to previously determined crystal structures, especially in the secondary structure, but exhibit new features: the loop composed of Pro 45 to Ser 48 adopts distinct conformations and the rings of tyrosines 51, 52, and 55 have reduced flipping rates. Amide protons with greatly reduced exchange rates are found predominantly in interior beta-strands and the alpha-helix, but also in the external 3/10 helix and edge beta-strand linked by the disulfide bond. Analysis of (15)N relaxation experiments (R1, R2, and NOE) at 600 MHz revealed five segments, consisting of residues 1-5, 28-31, 46-50, 60-65, 74-77, retaining flexibility in solution. The change in conformation entropy for RNase SA folding is smaller than previously believed, since the native protein is more flexible in solution than in a crystal.     

Molecular dynamics studies on HIV-1 protease drug resistance and folding pathways

Drug resistance to HIV-1 protease involves the accumulation of multiple mutations in the protein. We investigate the role of these mutations by using molecular dynamics simulations that exploit the influence of the native-state topology in the folding process. Our calculations show that sites contributing to phenotypic resistance of FDA-approved drugs are among the most sensitive positions for the stability of partially folded states and should play a relevant role in the folding process. Furthermore, associations between amino acid sites mutating under drug treatment are shown to be statistically correlated. The striking correlation between clinical data and our calculations suggest a novel approach to the design of drugs tailored to bind regions crucial not only for protein function, but for folding as well.     

Structural basis of oncogenic activation caused by point mutations in the kinase domain of the MET proto-oncogene: modeling studies

Missense mutations in the tyrosine kinase domain of the MET proto-oncogene occur in selected cases of papillary renal carcinoma. In biochemical and biological assays, these mutations produced constitutive activation of the MET kinase and led to tumor formation in nude mice. Some mutations caused transformation of NIH 3T3 cells. To elucidate the mechanism of ligand-independent MET kinase activation by point mutations, we constructed several 3D models of the wild-type and mutated MET catalytic core domains. Analysis of these structures showed that some mutations (e.g., V1110I, Y1248H/D/C, M1268T) directly alter contacts between residues from the activation loop in its inhibitory conformation and those from the main body of the catalytic domain; others (e.g., M1149T, L1213V) increase flexibility at the critical points of the tertiary structure and facilitate subdomain movements. Mutation D1246N plays a role in stabilizing the active form of the enzyme. Mutation M1268T affects the S+1 and S+3 substrate-binding pockets. Models implicate that although these changes do not compromise the affinity toward the C-terminal autophosphorylation site of the MET protein, they allow for binding of the substrate for the c-Abl tyrosine kinase. We provide biochemical data supporting this observation. Mutation L1213V affects the conformation of Tyr1212 in the active form of MET. Several somatic mutations are clustered at the surface of the catalytic domain in close vicinity of the probable location of the MET C-terminal docking site for cytoplasmic effectors.     

One Dog,but Which Dog? How Researchers Guide Participants to Select Dogs in Surveys of Human–Dog Relationships

Surveys and questionnaires are regularly used in studies of human–animal relationships. However, little attention has been given to understanding how survey participants are provided with instructions for the selection of a single animal within a multi-pet household, let alone the implications for reporting and interpreting data. We reviewed the instructions for the selection of an individual animal in studies addressing emotional or psychological attachment between people and dogs. By searching multidisciplinary journals from the year 2000 onwards, we identified a total of 128 papers, of which 63 met the inclusion criteria. Where selection criteria/instructions were not clear, authors were contacted. One in five studies (21%, or n = 13) did not report their instructions. When provided, instructions varied considerably. The most commonly provided direction was “favorite/closest relationship” (n = 12, or 19%). The remainder (n = 38, or 60%) were spread across eight different categories. Around half of the studies used a validated questionnaire that already contained an instruction, though a similar proportion of studies implemented author-designed instruments. Overall, the common absence and inconsistency of instructions for individual dog selection is taken to imply that there is no standard expectation or approach for instructions to be reported in studies of human relationships with dogs, or human–animal relationships more generally. We recommend further research on the implication of selection methods to ensure that instructions can be matched with specific research aims.     

Threading with chemostructural restrictions method for predicting fold and functionally significant residues: application to dipeptidylpeptidase IV (DPP-IV)

We present a new method for more accurate modeling of protein structure, called threading with chemostructural restrictions. This method addresses those cases in which a target sequence has only remote homologues of known structure for which sequence comparison methods cannot provide accurate alignments. Although remote homologues cannot provide an accurate model for the whole chain, they can be used in constructing practically useful models for the most conserved-and often the most interesting-part of the structure. For many proteins of interest, one can suggest certain chemostructural patterns for the native structure based on the available information on the structural superfamily of the protein, the type of activity, the sequence location of the functionally significant residues, and other factors. We use such patterns to restrict (1) a number of possible templates, and (2) a number of allowed chain conformations on a template. The latter restrictions are imposed in the form of additional template potentials (including terms acting as sequence anchors) that act on certain residues. This approach is tested on remote homologues of alpha/beta-hydrolases that have significant structural similarity in the positions of their catalytic triads. The study shows that, in spite of significant deviations between the model and the native structures, the surroundings of the catalytic triad (positions of C(alpha) atoms of 20-30 nearby residues) can be reproduced with accuracy of 2-3 A. We then apply the approach to predict the structure of dipeptidylpeptidase IV (DPP-IV). Using experimentally available data identifying the catalytic triad residues of DPP-IV (David et al., J Biol Chem 1993;268:17247-17252); we predict a model structure of the catalytic domain of DPP-IV based on the 3D fold of prolyl oligopeptidase (Fulop et al., Cell 1998;94:161-170) and use this structure for modeling the interaction of DPP-IV with inhibitor.     

Mutation of amino acids in the alpha 1,3-fucosyltransferase motif affects enzyme activity and Km for donor and acceptor substrates

Alpha 1,3-fucosyltransferases (FucT) share a conserved amino acid sequence designated the alpha 1,3 FucT motif that has been proposed to be important for nucleotide sugar binding. To evaluate the importance of the amino acids in this motif, each of the alpha 1,3 FucT motif amino acids was replaced with alanine (alanine scanning mutagenesis) in human FucT VI, and the resulting mutant proteins were analyzed for enzyme activity and kinetically characterized in those cases in which the mutant protein had sufficient activity. Two of the mutant proteins were inactive, six had less than 1% of wild-type activity, and four had approximately 10-50% of wild-type enzyme activity. Three of the mutant proteins with significant enzyme activity had substantially larger Km (5 to 15 times) for GDP-fucose than FucT VI wild-type enzyme. The fourth mutant protein with significant enzyme activity (S249A) had a Km at least 10 times larger than wild-type FucT VI for the acceptor substrate, with only a slightly larger (2-3 times) Km for GDP-fucose. Thus mutation of any of the amino acids within the alpha 1,3 FucT motif to Ala affects alpha 1,3-FucT activity, and substitution of Ala for some of the alpha 1,3 FucT motif amino acids results in proteins with altered kinetic constants for both the acceptor and donor substrates. Secondary structure prediction suggests a helix-loop-helix fold for the alpha 1,3 FucT motif, which can be used to rationalize the effects of mutations in terms of 3D structure.     

Comparative proteomic analysis of PAI-1 and TNF-alpha-derived endothelial microparticles

Endothelium-derived microparticles (EMPs) are small vesicles released from endothelial cells in response to cell injury, apoptosis, or activation. Elevated concentrations of EMPs have been associated with many inflammatory and vascular diseases. EMPs also mediate long range signaling and alter downstream cell function. Unfortunately, the molecular and cellular basis of microparticle production and downstream cell function is poorly understood. We hypothesize that EMPs generated by different agonists will produce distinct populations of EMPs with unique protein compositions. To test this hypothesis, different EMP populations were generated from human umbilical vein endothelial cells by stimulation with plasminogen activator inhibitor type 1 (PAI-1) or tumor necrosis factor-alpha (TNF-alpha) and subjected to proteomic analysis by LC/MS. We identified 432 common proteins in all EMP populations studied. Also identified were 231 proteins unique to control EMPs, 104 proteins unique to PAI-1 EMPs and 70 proteins unique to TNF-alpha EMPs. Interestingly, variations in protein abundance were found among many of the common EMP proteins, suggesting that differences exist between EMPs on a relative scale. Finally, gene ontology (GO) and KEGG pathway analysis revealed many functional similarities and few differences between the EMP populations studied. In summary, our results clearly indicate that EMPs generated by PAI-1 and TNF-alpha produce EMPs with overlapping but distinct protein compositions. These observations provide fundamental insight into the mechanisms regulating the production of these particles and their physiological role in numerous diseases.     

3',5' Cyclic nucleotide phosphodiesterases class III: members, structure, and catalytic mechanism

3',5' Cyclic nucleotide phosphodiesterases (PDEs) comprise a superfamily of enzymes that were previously divided by their primary structure into two major classes: PDE class I and II. The 3',5' cyclic AMP phosphodiesterase from Escherichia coli encoded by the cpdA gene does not show any homology to either PDE class I or class II enzymes and, therefore, represents a new, third class of PDEs. Previously, information about essential structural elements, substrate and cofactor binding sites, and the mechanism of catalysis was unknown for this enzyme. The present study shows by computational analysis that the enzyme encoded by the E. coli cpdA gene belongs to a family of phosphodiesterases that closely resembles the catalytic machinery known from purple acid phosphatases and several other dimetallophosphoesterases. They share both the conserved sequence motif, D-(X)(n) GD-(X)(n)-GNH[E/D]-(X)(n)-H-(X)(n)-GHXH, which contains the invariant residues forming the active site of purple acid phosphatases, a binuclear Fe(3+)-Me(2+)-containing center, as well as a beta(alpha)beta(alpha)beta motif as a typical secondary structure signature. Furthermore, the known biochemical properties of the bacterial phosphodiesterase encoded by the cpdA gene, such as the requirement of iron ions and a reductant for maintaining its catalytic activity, support this hypothesis developed by computational analysis. In addition, the availability of atomic coordinates for several purple acid phosphatases and related proteins allowed the generation of a three-dimensional model for class III cyclic nucleotide phosphodiesterases.     

Neuronal nicotinic receptor deficits in Alzheimer patients with the Swedish amyloid precursor protein 670/671 mutation

The influence of beta-amyloid on cholinergic neurotransmission was studied by measuring alterations in nicotinic acetylcholine receptors (nAChRs) in autopsy brain tissue from subjects carrying the Swedish amyloid precursor protein (APP) 670/671 mutation. Significant reductions in numbers of nAChRs were observed in various cortical regions of the Swedish 670/671 APP mutation family subjects (-73 to -87%) as well as in sporadic Alzheimer's disease (AD) cases (-37 to -57%) using the nicotinic agonists [3H]epibatidine and [3H]nicotine, which bind with high affinity to both alpha3 and alpha4 and to alpha4 nAChR subtypes, respectively. Saturation binding studies with [3H]epibatidine revealed two binding sites in the parietal cortex of AD subjects and controls. A significant decrease in Bmax (-82%) for the high-affinity site was observed in APP 670/671 subjects with no change in K(D) compared with controls (0.018 nM APP 670/671; 0.036 nM control). The highest load of neuronal plaques (NPs) was observed in the parietal cortex of APP 670/671 brains, whereas the number of [3H]nicotine binding sites was less impaired compared with other cortical brain regions. Except for a positive significant correlation between the number of [3H]nicotine binding sites and number of NPs in the parietal cortex, no strict correlation was observed between nAChR deficits and the presence of NPs and neurofibrillary tangles, suggesting that these different processes may be closely related but not strictly dependent on each other.