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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Elucidating protein secondary structures using alpha-carbon recurrence quantifications.

Secondary structures of proteins were studied by recurrence quantification analysis (RQA). High-resolution, 3-dimensional coordinates of alpha-carbon atoms comprising a set of 68 proteins were downloaded from the Protein Data Bank. By fine-tuning four recurrence parameters (radius, line, residue, separation), it was possible to establish excellent agreement between percent contribution of alpha-helix and beta-sheet structures determined independently by RQA and that of the DSSP algorithm (Define Secondary Structure of Proteins). These results indicate that there is an equivalency between these two techniques, which are based upon totally different pattern recognition strategies. RQA enhances qualitative contact maps by quantifying the arrangements of recurrent points of alpha carbons close in 3-dimensional space. For example, the radius was systematically increased, moving the analysis beyond local alpha-carbon neighborhoods in order to capture super-secondary and tertiary structures. However, differences between proteins could only be detected within distances up to about 6-11 A, but not higher. This result underscores the complexity of alpha-carbon spacing when super-secondary structures appear at larger distances. Finally, RQA-defined secondary structures were found to be robust against random displacement of alpha carbons upwards of 1 A. This finding has potential import for the dynamic functions of proteins in motion.     

Combining mutations in HIV-1 protease to understand mechanisms of resistance

HIV-1 develops resistance to protease inhibitors predominantly by selecting mutations in the protease gene. Studies of resistant mutants of HIV-1 protease with single amino acid substitutions have shown a range of independent effects on specificity, inhibition, and stability. Four double mutants, K45I/L90M, K45I/V82S, D30N/V82S, and N88D/L90M were selected for analysis on the basis of observations of increased or decreased stability or enzymatic activity for the respective single mutants. The double mutants were assayed for catalysis, inhibition, and stability. Crystal structures were analyzed for the double mutants at resolutions of 2.2-1.2 A to determine the associated molecular changes. Sequence-dependent changes in protease-inhibitor interactions were observed in the crystal structures. Mutations D30N, K45I, and V82S showed altered interactions with inhibitor residues at P2/P2', P3/P3'/P4/P4', and P1/P1', respectively. One of the conformations of Met90 in K45I/L90M has an unfavorably close contact with the carbonyl oxygen of Asp25, as observed previously in the L90M single mutant. The observed catalytic efficiency and inhibition for the double mutants depended on the specific substrate or inhibitor. In particular, large variation in cleavage of p6(pol)-PR substrate was observed, which is likely to result in defects in the maturation of the protease from the Gag-Pol precursor and hence viral replication. Three of the double mutants showed values for stability that were intermediate between the values observed for the respective single mutants. D30N/V82S mutant showed lower stability than either of the two individual mutations, which is possibly due to concerted changes in the central P2-P2' and S2-S2' sites. The complex effects of combining mutations are discussed.     

Structural and kinetic analysis of drug resistant mutants of HIV-1 protease.

Mutants of HIV-1 protease that are commonly selected on exposure to different drugs, V82S, G48V, N88D and L90M, showed reduced catalytic activity compared to the wild-type protease on cleavage site peptides, CA-p2, p6pol-PR and PR-RT, critical for viral maturation. Mutant V82S is the least active (2-20% of wild-type protease), mutants N88D, R8Q, and L90M exhibit activities ranging from 20 to 40% and G48V from 50 to 80% of the wild-type activity. In contrast, D30N is variable in its activity on different substrates (10-110% of wild-type), with the PR-RT site being the most affected. Mutants K45I and M46L, usually selected in combination with other mutations, showed activities that are similar to (60-110%) or greater than (110-530%) wild-type, respectively. No direct relationship was observed between catalytic activity, inhibition, and structural stability. The mutants D30N and V82S were similar to wild-type protease in their stability toward urea denaturation, while R8Q, G48V, and L90M showed 1.5 to 2.7-fold decreased stability, and N88D and K45I showed 1.6 to 1.7-fold increased stability. The crystal structures of R8Q, K45I and L90M mutants complexed with a CA-p2 analog inhibitor were determined at 2.0, 1.55 and 1.88 A resolution, respectively, and compared to the wild-type structure. The intersubunit hydrophobic contacts observed in the crystal structures are in good agreement with the relative structural stability of the mutant proteases. All these results suggest that viral resistance does not arise by a single mechanism.     

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.     

Resistance mechanism revealed by crystal structures of unliganded nelfinavir-resistant HIV-1 protease non-active site mutants N88D and N88S

Nelfinavir is an inhibitor of HIV-1 protease, and is used for treatment of patients suffering from HIV/AIDS. However, treatment results in drug resistant mutations in HIV-1 protease. N88D and N88S are two such mutations which occur in the non-active site region of the enzyme. We have determined crystal structures of unliganded N88D and N88S mutants of HIV-1 protease to resolution of 1.65 Å and 1.8 Å, respectively. These structures refined against synchrotron data lead to R-factors of 0.1859 and 0.1780, respectively. While structural effects of N88D are very subtle, the mutation N88S has caused a significant conformational change in D30, an active site residue crucial for substrate and inhibitor binding.     

Three-dimensional structure of a monomeric form of a retroviral protease

The assembly of Mason-Pfizer monkey virus Gag polyproteins into immature capsids and their cleavage by the encoded protease are temporally and spatially separated processes, making the virus a particularly useful model for investigation of protease activation. Here we present a high resolution NMR structure of a fully folded monomer of a 12 kDa M-PMV protease (wt 12 PR) and of a Cys7Ala/Asp26Asn/Cys106Ala mutant (12 PR(D26N/C7A/C106A)). The overall structures of both wt 12 PR and 12 PR(D26N/C7A/C106A) follow the conservative structural motif of other retroviral proteases. The most prominent difference from the canonical fold of retroviral proteases is the absence of the interfacial beta-sheet, which leads to the loss of the principal force stabilizing the dimer of M-PMV PR. The monomer-dimer equilibrium can be shifted in favor of the dimer by adding a substrate or an inhibitor, partially compensating for the missing role of the beta-sheet. We also show that cysteines C7 and C106 play a crucial role in stabilizing the dimer and consequently increasing the proteolytic activity of M-PMV PR. This is consistent with the role of reversible oxidative modification of the cysteine residues in the regulation of the maturation of assembled M-PMV capsids in the cytoplasm.     

Comprehensive mutagenesis of HIV-1 protease: a computational geometry approach

A computational geometry technique based on Delaunay tessellation of protein structure, represented by C(alpha) atoms, is used to study effects of single residue mutations on sequence-structure compatibility in HIV-1 protease. Profiles of residue scores derived from the four-body statistical potential are constructed for all 1881 mutants of the HIV-1 protease monomer and compared with the profile of the wild-type protein. The profiles for an isolated monomer of HIV-1 protease and the identical monomer in a dimeric state with an inhibitor are analyzed to elucidate changes to structural stability. Protease residues shown to undergo the greatest impact are those forming the dimer interface and flap region, as well as those known to be involved in inhibitor binding.     

Molecular mechanism for dimerization to regulate the catalytic activity of human cytomegalovirus protease

Biochemical studies indicate that dimerization is required for the catalytic activity of herpesvirus proteases, whereas structural studies show a complete active site in each monomer, away from the dimer interface. Here we report kinetic, biophysical and crystallographic characterizations of structure-based mutants in the dimer interface of human cytomegalovirus (HCMV) protease. Such mutations can produce a 1,700-fold reduction in the kcat while having minimal effects on the K(m). Dimer stability is not affected by these mutations, suggesting that dimerization itself is insufficient for activity. There are large changes in monomer conformation and dimer organization of the apo S225Y mutant enzyme. However, binding of an activated peptidomimetic inhibitor induced a conformation remarkably similar to the wild type protease. Our studies suggest that appropriate dimer formation may be required to indirectly stabilize the protease oxyanion hole, revealing a novel mechanism for dimerization to regulate enzyme activity.     

Drug resistance mutations can effect dimer stability of HIV-1 protease at neutral pH.

The monomer-dimer equilibrium for the human immunodeficiency virus type 1 (HIV-1) protease has been investigated under physiological conditions. Dimer dissociation at pH 7.0 was correlated with a loss in beta-sheet structure and a lower degree of ANS binding. An autolysis-resistant mutant, Q7K/L33I/L63I, was used to facilitate sedimentation equilibrium studies at neutral pH where the wild-type enzyme is typically unstable in the absence of bound inhibitor. The dimer dissociation constant (KD) of the triple mutant was 5.8 microM at pH 7.0 and was below the limit of measurement (approximately 100 nM) at pH 4.5. Similar studies using the catalytically inactive D25N mutant yielded a KD value of 1.0 microM at pH 7.0. These values differ significantly from a previously reported value of 23 nM obtained indirectly from inhibitor binding measurements (Darke et al., 1994). We show that the discrepancy may result from the thermodynamic linkage between the monomer-dimer and inhibitor binding equilibria. Under conditions where a significant degree of monomer is present, both substrates and competitive inhibitors will shift the equilibrium toward the dimer, resulting in apparent increases in dimer stability and decreases in ligand binding affinity. Sedimentation equilibrium studies were also carried out on several drug-resistant HIV-1 protease mutants: V82F, V82F/I84V, V82T/I84V, and L90M. All four mutants exhibited reduced dimer stability relative to the autolysis-resistant mutant at pH 7.0. Our results indicate that reductions in drug affinity may be due to the combined effects of mutations on both dimer stability and inhibitor binding.     

Effect of the active site D25N mutation on the structure, stability, and ligand binding of the mature HIV-1 protease

All aspartic proteases, including retroviral proteases, share the triplet DTG critical for the active site geometry and catalytic function. These residues interact closely in the active, dimeric structure of HIV-1 protease (PR). We have systematically assessed the effect of the D25N mutation on the structure and stability of the mature PR monomer and dimer. The D25N mutation (PR(D25N)) increases the equilibrium dimer dissociation constant by a factor >100-fold (1.3 +/- 0.09 microm) relative to PR. In the absence of inhibitor, NMR studies reveal clear structural differences between PR and PR(D25N) in the relatively mobile P1 loop (residues 79-83) and flap regions, and differential scanning calorimetric analyses show that the mutation lowers the stabilities of both the monomer and dimer folds by 5 and 7.3 degrees C, respectively. Only minimal differences are observed in high resolution crystal structures of PR(D25N) complexed to darunavir (DRV), a potent clinical inhibitor, or a non-hydrolyzable substrate analogue, Ac-Thr-Ile-Nle-r-Nle-Gln-Arg-NH(2) (RPB), as compared with PR.DRV and PR.RPB complexes. Although complexation with RPB stabilizes both dimers, the effect on their T(m) is smaller for PR(D25N) (6.2 degrees C) than for PR (8.7 degrees C). The T(m) of PR(D25N).DRV increases by only 3 degrees C relative to free PR(D25N), as compared with a 22 degrees C increase for PR.DRV, and the mutation increases the ligand dissociation constant of PR(D25N).DRV by a factor of approximately 10(6) relative to PR.DRV. These results suggest that interactions mediated by the catalytic Asp residues make a major contribution to the tight binding of DRV to PR.