PubMed HealthSearch

Biomedical subjects

T A Darden

Publications and source records attributed to T A Darden.

At least 19 recordsLinked to original sources

Vertical-scanning mutagenesis of a critical tryptophan in the minor groove binding track of HIV-1 reverse transcriptase. Molecular nature of polymerase-nucleic acid interactions.

While sequence-specific DNA-binding proteins interact predominantly in the DNA major groove, DNA polymerases bind DNA through interactions in the minor groove that are sequence nonspecific. Through functional analyses of alanine-substituted mutant enzymes that were guided by molecular dynamics modeling of the human immunodeficiency virus type 1-reverse transcriptase and DNA complex, we previously identified a structural element in reverse transcriptase, the minor groove binding track (MGBT). The MGBT is comprised of five residues (Ile94, Gln258, Gly262, Trp266, and Gln269) which interact 2-6 base pairs upstream from the polymerase active site in the DNA minor groove and are important in DNA binding, processivity, and frameshift fidelity. These residues do not contribute equally; functional analysis of alanine mutants suggests that Trp266 contributes the most to binding. To define the molecular interactions between Trp266 and the DNA minor groove, we have analyzed the properties of eight mutants, each with an alternate side chain at this position. A refined molecular dynamics model was used to calculate relative binding free energies based on apolar surface area buried upon complex formation. In general, there was a strong correlation between the relative calculated binding free energies for the alternate residue 266 side chains and the magnitude of the change in the properties which reflect template-primer interactions (template-primer dissociation rate constant, Ki,AZTTP, processivity, and frameshift fidelity). This correlation suggests that hydrophobic interactions make a major contribution to the stability of the polymerase-DNA complex. Additionally, tyrosine and arginine substitutions resulted in mutant enzymes with DNA binding properties better than predicted by buried surface area alone, suggesting that hydrogen bonding could also play a role in DNA binding at this position.

Antiviral Agents

Trans-cis isomerization of proline 22 in bovine prothrombin fragment 1: a surprising result of structural characterization.

The calcium ion-mediated interaction of bovine prothrombin (BF1) with negatively charged phospholipid membranes is assumed to be largely via the Gla domain of BF1 with the fold of the Gla domain essential for binding. It has been reported that Pro22 undergoes classical trans to cis isomerization in the presence of calcium ions with the cis conformation of Pro22 of BF1 responsible for membrane binding [Evans, T. C., Jr., and Nelsestuen, G. L. (1996) Biochemistry 35, 8210-8215]. However, Pro22 was found to be in the trans conformation in the crystal structure of BF1. In the present work, we have used molecular dynamics simulations to investigate the relative importance of the two conformations of Pro22 to the structural and dynamical properties of BF1. The initial trans conformation of Pro22 in BF1 was slowly converted to cis-Pro22 using constrained dynamics. The second-generation AMBER force field in conjunction with the particle mesh Ewald method to accommodate long-range interaction was employed in the trajectory calculations. Comparison of the BF1(trans-Pro22) and BF1(cis-Pro22) equilibrated structures reveals surprisingly that the overall structural changes associated with the trans-cis isomerization is minimal and only minor modifications to the hydrogen bond network and the network of N-terminus Ala1 take place. The calculated electrostatic potential energy surfaces of the two protein structures also appear to be very similar, indicating the near equality of the local interaction site environments in the protein prior to lipid binding.

1-Carboxyglutamic Acid

Reciprocal size-effect relationship of the key residues in determining regio- and stereospecificities of DHEA hydroxylase activity in P450 2a5.

Collectively, the P450 2a4/2a5 system hyrdoxylates DHEA in at least three positions (7alpha, 7beta, and 2alpha). An individual P450, however, exhibits high specificity to one of these products. Using site-directed mutagenesis of mP450 2a5 from the wild mouse Mus minutoides and bacterial expression, we have associated the function of residues 117, 209, and 481 with the respective specificity observed in each P450. Ala at position 117 determines the 7beta-hydroxylase activity, whereas Val at this position defines the 2alpha-hydroxylase activity. Leu at position 209 is essential for high DHEA 7alpha-hydroxylase activity. The substitutions of residue 481 with various hydrophobic amino acids elicited a profound alteration of the specific hydroxylation rates, but did not influence the regio- and stereospecificities at either of the three positions of DHEA. The alterations caused by residue 481 also depended on the residue identity at position 117 or 209. The results indicate that the sizes of several key residues obey a concerted reciprocal relationship whereby the substrate pocket of the P450s adjusts to accommodate DHEA. A limited molecular modeling study successfully correlates DHEA binding to experimental DHEA hydroxylase activities for a series of mutants at key positions.

Amino Acid Sequence

Refinement of the NMR solution structure of the gamma-carboxyglutamic acid domain of coagulation factor IX using molecular dynamics simulation with initial Ca2+ positions determined by a genetic algorithm.

A genetic algorithm (GA) successfully identified the calcium positions in the crystal structure of bovine prothrombin fragment 1 bound with calcium ions (bf1/Ca). The same protocol was then used to determine the calcium positions in a closely related fragment, the Gla domain of coagulation factor IX, the structure of which had previously been determined by NMR spectroscopy in the presence of calcium ions. The most frequently occurring low-energy structure found by GA was used as the starting structure for a molecular dynamics refinement. The molecular dynamics simulation was performed using explicit water and the Particle-Mesh Ewald method to accommodate the long-range electrostatic forces. While the overall conformation of the NMR structure was preserved, significant refinement is apparent when comparing the simulation average structure with its NMR precursor in terms of the N-terminal (Tyr1-N) network, the total number of hydrogen bonds, the calcium ion coordinations, and the compactness of the structure. It is likely that the placement of calcium ions in the protein is critical for refinement. The calcium ions apparently induce structural changes during the course of the simulation that result in a more compact structure.

1-Carboxyglutamic Acid

Structural flexibility and functional versatility of cytochrome P450 and rapid evolution.

P450 represents a large group of heme-thiolate enzymes that exhibit remarkably diverse activities for the metabolism of numerous endogenous and exogenous chemicals. Recent site-directed mutagenesis studies indicate that a single mutation at any of the key residues can be enough to alter the substrate and/or product specificities in the P450 activities. Molecular modeling predicts that these key residues are located within the substrate heme pocket. Structural elements involved in diversifying P450 activity appear to correspond to the B' helix, the F helix and the F/G interhelical loop in the bacterial P450s. Structures represented by these regions are extremely variable despite the fact that the core of the P450 substrate pocket is well conserved. A mutation within these regions may result in a significant geometrical alteration of the pocket and lead to diversify the P450 activity. Phylogenetical analysis shows a relatively high rate of nonsynonymous substitution within these substrate binding regions. The functional versatility of P450 can thus be largely accounted for in terms of pocket change brought about by rapid mutations.

Animals

The roles of individual amino acids in altering substrate specificity of the P450 2a4/2a5 enzymes.

A single amino acid substitution is sufficient to alter substrate specificity of P450 enzymes. Mouse P450 2a5, for example, has its substrate specificity converted from coumarin 7- to testosterone 15 alpha-hydroxylase activity by the substitution of Phe at position 209 to Leu. Furthermore, placing Asn at this position confers a novel corticosterone 15 alpha-hydroxylase activity to this P450. Recent site-directed mutational studies show the presence of the topologically common residues, each of which can determine the specificities of various mammalian P450s. For instance, residue 209 (in 2a5) corresponds to a residue at position 206 in rat P4502B1 that regulates its steroid hydroxylase activity. High substrate specificity often observed in an individual P450, therefore, can be determined and altered by the identities of a few critical residues. The structural flexibility of the substrate-heme pocket may also provide P450 enzymes with the ability to display a broad range of substrate specificities. Understanding the underlying principles whereby the flexible pocket determines P450 activities may lead us to the prediction of P450 activities based on the identities of key amino acid residues.

Amino Acids

Structural flexibility and functional versatility of mammalian P450 enzymes.

P450 enzymes have evolved into a large superfamily that displays great diversity in substrate and product specificities by fixing the natural amino acid substitutions with high frequency. Site-directed mutagenesis has been used to correlate the substitutions with the diverse specificities in various P450s. As a result, the common residues that determine the specificities of various mammalian P450s have been identified and aligned to the corresponding residues in the substrate-heme pocket of the 3-dimensional structures of bacterial P450s. The substrate-heme pocket appears to be structurally variable so that only a minor substitution (Ala -> -> Val, for example) at the critical positions is enough to define the altered specificity. Thus, the structural variability of the P450s provides the inherent versatility in acquiring a novel activity. Recent mutational studies indicate that the side chain size is the major determining factor of specificity, outweighing other factors such as polarity. Further understanding of the paradoxical characteristics observed may provide us with the underlying principles that determine P450 activities, and may lead to the ability to predict P450 activities based on the types of key amino acid residues.

Amino Acid Sequence

Selected new developments in computational chemistry.

Molecular dynamics is a general technique for simulating the time-dependent properties of molecules and their environments. Quantum mechanics, as applied to molecules or clusters of molecules, provides a prescription for predicting properties exactly (in principle). It is reasonable to expect that both will have a profound effect on our understanding of environmental chemistry in the future. In this review, we consider several recent advances and applications in computational chemistry.

Chemical Phenomena

Reduced frameshift fidelity and processivity of HIV-1 reverse transcriptase mutants containing alanine substitutions in helix H of the thumb subdomain.

We have analyzed two human immunodeficiency virus (HIV-1) reverse transcriptase mutants of helix H in the thumb subdomain suggested by x-ray crystallography to interact with the primer strand of the template-primer. These enzymes, G262A and W266A, were previously shown to have greatly elevated dissociation rate constants for template-primer and to be much less sensitive to inhibition by 3'-azidodeoxythymidine 5'-triphosphate. Here we describe their processivity and error specificity. The results reveal that: (i) both enzymes have reduced processivity and lower fidelity for template-primer slippage errors, (ii) they differ from each other in sequence-dependent termination of processive synthesis and in error specificity, and (iii) the magnitude of the mutator effect relative to wild-type enzyme for deletions in homopolymeric sequences decreases as the length of the run increases. Thus amino acid substitutions in a subdomain thought to interact with the duplex template-primer confer a strand slippage mutator phenotype to a replicative DNA polymerase. This suggests that interactions between specific amino acids and the primer stem at positions well removed from the active site are critical determinants of processivity and fidelity. These effects, obtained in aqueous solution during catalytic cycling, are consistent with and support the existing crystallographic structural model.

Alanine

Altering the regiospecificity of androstenedione hydroxylase activity in P450s 2a-4/5 by a mutation of the residue at position 481.

Mouse P450 2a-5 (coumarin 7-hydroxylase) acquires androstenedione (AD) hydroxylase activity by substituting Phe at position 209 with Asn. However, this mutant P450 2a-5 (F209N) and the corresponding mutant P450 2a-4 (L209N) exhibit different regiospecificites of androstenedione (AD) hydroxylase activity. While the former mutant catalyzes both AD 15 alpha- and 7 alpha-hydroxylase activities at similar rates, the latter mutant maintains the original high specificity of AD 15 alpha-hydroxylase activity. The AD hydroxylase activities in chimeric enzymes of the mutants L209N and F209N show that the regiospecificites are determined by the carboxy-terminal halves of the P450 molecules. Mutations at each of the four different residues within the carboxy-terminal halves indicate that the differences in regiospecificity are determined by the Val/Ala mutation at position 481. As the size of the hydrophobic amino acid at position 481 becomes larger (Ala < Val < Ile), the regiospecificities toward the C15 position of the AD molecule are dramatically increased. The regiospecificity is also increased by placing positively-charged Arg at position 481, although the remaining 15 alpha-hydroxylase activity in this mutant is considerably lower than the other P450s. The results indicate that the size of the residue at position 481 is a key factor in regulating the regiospecificity of AD hydroxylase activity in the P450s. Modeling AD in the substrate-heme pocket of bacterial P450 101A provided further support that residue 481 may reside near the steroid molecule so as to possibly affect the AD hydroxylase activity.

Amino Acid Sequence

Multiple steroid-binding orientations: alteration of regiospecificity of dehydroepiandrosterone 2- and 7-hydroxylase activities of cytochrome P-450 2a-5 by mutation of residue 209.

The mutation of Ala-117 to Val conferred dehydroepiandrosterone (DHEA) hydroxylase activity on cytochrome P-450 2a-4, with the production of both 2 alpha- and 7 alpha-hydroxyDHEA at similar rates. P-450 2a-5 which has Val at position 117, acquired high DHEA hydroxylase activity by mutation of Phe-209. Mutant F209L of P-450 2a-5 exhibited strong regiospecificity at the 2-position of the DHEA molecule with the production of 2 alpha-hydroxy DHEA as the major metabolite. On the other hand, mutant F209V of P-450 2a-5 showed the 7-position to be the major hydroxylation site, 7 beta-hydroxyDHEA and 7 alpha-OHDHEA being produced. Therefore the regiospecificity of DHEA hydroxylase activity of P-450 2a-5 is altered between the 2- and 7-position depending on the amino acid at position 209. Modelling of the DHEA molecule in the pocket of bacterial P-450cam showed that the steroid can be accommodated in at least two orientations for which the 2- or 7- position is near the sixth axial position of the haem. Moreover, these two orientations, which are of similar energy, can be interconverted by a 180 degrees rotation of the steroid molecule around its long axis. These results support the hypothesis that the steroid molecule in the pocket is in dynamic equilibrium with multiple binding orientations and that the equilibrium is apparently determined by a few critical residues including those at positions 117 and 209.

Aryl Hydrocarbon Hydroxylases

Atomic-level accuracy in simulations of large protein crystals.

Proper treatment of long-range Coulombic forces presents a major obstacle to providing realistic molecular dynamics simulations of macromolecules. Traditional approximations made to lessen computational cost ultimately lead to unrealistic behavior. The particle mesh Ewald method accommodates long-range Coulombic forces accurately and efficiently by use of fast Fourier transform techniques. We report a 1-ns simulation of bovine pancreatic trypsin inhibitor in a crystal unit cell using the particle mesh Ewald methodology. We find an rms backbone deviation from the x-ray structure (0.33 A) that is lower than that observed between bovine pancreatic trypsin inhibitor in different crystal forms and much lower than those of previous simulations. These results bridge the gap between structures obtained from molecular simulation and those from experiment.

Animals

Inherent versatility of P-450 oxygenase. Conferring dehydroepiandrosterone hydroxylase activity to P-450 2a-4 by a single amino acid mutation at position 117.

Mouse steroid 15 alpha-hydroxylase P-450 2a-4 is restricted in its substrate specificity to the delta 4, 3-ketone steroids such as androstenedione. As a result, the P-450 exhibits little hydroxylase activity toward delta 5, 3-hydroxysteroids including dehydroepiandrosterone (DHEA). A single amino acid mutation of Ala at position 117 to Val, however, is enough to confer a high DHEA hydroxylase activity to P-450 2a-4 with 7 alpha-OH DHEA as one of the two major hydroxylated metabolites. Mouse coumarin 7-hydroxylase P-450 2a-5 contains Val at position 117, but it exhibits very low DHEA hydroxylase activity. P-450 2a-5 acquires high DHEA hydroxylase activity, however, by a mutation of Phe-209 to Asn. Moreover, the mutant P-450 2a-5 loses its activity when Val is replaced by Ala at position 117. The residue at position 117, therefore, plays the principal role in the determination of the DHEA hydroxylase activity of the P-450s. Conversely, mutations at residue 117 have little effect on the androstenedione hydroxylase activities of the P-450s. Further modeling of the DHEA binding orientation in the substrate-heme pocket of bacterial P-450cam (Iwasaki, M., Darden, T., Pedersen, L., Davis, D. G., Juvonen, R. O., Sueyoshi, T., and Negishi, M. (1993) J. Biol. Chem. 268, 759-762) provides support for the hypothesis that the type of residue at position 117 determines the steroid-substrate specificity of the P-450 depending on the substituent at the C3 position of steroid molecule.

Androstenedione

Role of glutamine-61 in the hydrolysis of GTP by p21H-ras: an experimental and theoretical study.

The active GTP-bound form of p21ras is converted to the biologically inactive GDP-bound form by enzymatic hydrolysis and this function serves to regulate the wild-type ras protein. The side chain of the amino acid at position 61 may play a key role in this hydrolysis of GTP by p21. Experimental studies that define properties of the Q61E mutant of p21H-ras are presented along with supporting molecular dynamics simulations. We find that under saturating concentrations of GTP the Q61E mutant of p21H-ras has a 20-fold greater rate of intrinsic hydrolysis (kcat = 0.57 min-1) than the wild type. The affinity of the Q61E variant for GTP (Kd = 115 microM) is much lower than that of the wild type. GTPase activating protein does not activate the variant. From molecular dynamics simulations, we find that both the wild type and Q61E mutant have the residue 61 side chain in transient contact with a water molecule that is well-positioned for hydrolytic attack on the gamma phosphate. Thr-35 also is found to form a transient hydrogen bond with this critical water. These elements may define the catalytic complex for hydrolysis of the GTP [Pai et al. (1990) EMBO J. 9, 2351]. Similarly, the G12P mutant, which also has an intrinsic hydrolysis rate similar to the wild type, is found to form the same complex in simulation. In contrast, molecular dynamics analysis of the mutants G12R, G12V, and Q61L, which have much lower intrinsic rates than the wild-type p21, do not show this complex.(ABSTRACT TRUNCATED AT 250 WORDS)

Binding Sites

Molecular dynamics simulation of HIV-1 protease in a crystalline environment and in solution.

Simulations of the unbound form of the human immunodeficiency virus type 1 protease have been carried out to 200 ps in a crystalline environment and in solution. Solution simulations were performed with and without charge-balancing counterions. The results are compared with the 2.8-A crystallographic structure of Wlodawer et al. [(1989) Science 245, 616], and a proposed model for the solution structure which involves local refolding of the flap regions is presented. The simulations suggest the crystal packing environment of the protease dimer stabilizes the flaps in an extended conformation. Solvation of the dimer leads to local refolding of the flaps which contract toward the active site, forming increased overlap and stronger intersubunit hydrogn bonding at the tips. The degree to which the flaps overlap in solution is observed to depend on the charge state of the system.

Amino Acid Sequence

Engineering mouse P450coh to a novel corticosterone 15 alpha-hydroxylase and modeling steroid-binding orientation in the substrate pocket.

The F209L mutation alters specificity of P450coh from coumarin 7-hydroxylation to 15 alpha-hydroxylation of 11-deoxysteroids such as testosterone and 11-deoxycorticosterone. Neither the wild-type nor F209L exhibits activity toward 11 beta-hydroxysteroids including corticosterone. Mutation of Phe-209 to Asn, however, confers on mutant F209N a high corticosterone 15 alpha-hydroxylase activity. F209V also exhibits low corticosterone 15 alpha-hydroxylase activity; Km and Vmax are 10-fold higher and lower, respectively, than for F209N. The results are consistent with the hypothesis that direct interaction of Asn-209 with 11OH is responsible for high corticosterone 15 alpha-hydroxylase activity. To support this hypothesis, a possible steroid-binding orientation is modeled in the substrate pocket of P450cam. Our weighted homology and constrained alignments map residue 209 of P450coh to Met-184 and Met-191 of P450cam. Energy minimization of corticosterone in the substrate pocket results in the 11OH of the steroid directed toward Met-184 (7 A) and Met-191 (16 A), and in C15 located near the sixth axial position of the heme. The steroid-binding model suggests that the P450cam's substrate pocket may be conserved in the mammalian P450 and can accommodate a steroid molecule, and that residue 209 appears to be located at the critical site that determines the steroid-substrate specificity of a P450 depending on the type of group at the 11-position of steroid molecule.

Amino Acid Sequence

Molecular modeling studies suggest that zinc ions inhibit HIV-1 protease by binding at catalytic aspartates.

Human immunodeficiency virus type 1 protease is inhibited in vitro by zinc ions at neutral pH. The binding site of these ions is not known; however, experimental data suggest that binding may occur in the active site. To examine the possibility of zinc binding in the active site, molecular dynamics simulations in the presence and absence of zinc have been carried out to 200 psec. The results are compared with the 2.8-A crystallographic structures of a synthetic HIV-1 protease, and a zinc binding site at the catalytic aspartate residues (Asp-25, Asp-25') is proposed. Molecular dynamics simulations show that the zinc ion remains stably bound in this region, coordinating the carboxylate side chains of both aspartate residues. Interaction with zinc does not disrupt the dimeric structure of the protein or significantly alter the structure of the active site. These data are consistent with experimental studies of HIV-1 protease inhibition by zinc and give strong evidence that this is the binding site that leads to inactivation.

Aspartic Acid