Search PubMed⌕ Search

Biomedical subjects

T A Darden

Publications and source records attributed to T A Darden.

At least 19 recordsLinked to original sources
PubMed2003

An all-atom solution-equilibrated model for human extrinsic blood coagulation complex (sTF-VIIa-Xa): a protein-protein docking and molecular dynamics refinement study.

Tissue factor (TF)-bound factor (F)VIIa plays a critical role in activating FX, an event that rapidly results in blood coagulation. Despite recent advances in the structural information about soluble TF (sTF)-bound VIIa and Xa individually, the atomic details of the ternary complex are not known. As part of our long-term goal to provide a structural understanding of the extrinsic blood coagulation pathway, we built an all atom solution-equilibrated model of the human sTF-VIIa-Xa ternary complex using protein-protein docking and molecular dynamics (MD) simulations. The starting structural coordinates of sTF-VIIa and Xa were derived from dynamically equilibrated solution structures. Due to the flexible nature of the light-chain of the Xa molecule, a three-stage docking approach was employed in which SP (Arg195-Lys448)/EGF2 (Arg86-Arg139), EGF1 (Asp46-Thr85) and GLA (Ala1-Lys45) domains were docked in a sequential manner. The rigid-body docking approach of the FTDOCK method in conjunction with filtering based on biochemical knowledge from experimental site-specific mutagenesis studies provided the strategy. The best complex obtained from the docking experiments was further refined using MD simulations for 3 ns in explicit water. In addition to explaining most of the known experimental site-specific mutagenesis data pertaining to sTF-VIIa, our model also characterizes likely enzyme-binding exosites on FVIIa and Xa that may be involved in the ternary complex formation. According to the equilibrated model, the 140s loop of VIIa serves as the key recognition motif for complex formation. Stable interactions occur between the FVIIa 140s loop and the FXa -strand B2 region near the sodium-binding domain, the 160 s loop and the N-terminal activation loop regions. The helical-hydrophobic stack region that connects the GLA and EGF1 domains of VIIa and Xa appears to play a potential role in the membrane binding region of the ternary complex. The proposed model may serve as a reasonable structural basis for understanding the exosite-mediated substrate recognition of sTF-VIIa and to advance understanding of the TFPI-mediated regulatory pathway of the extrinsic blood coagulation cascade.

Binding Sites↗
PubMed2001

Gene selection for sample classification based on gene expression data: study of sensitivity to choice of parameters of the GA/KNN method.

MOTIVATION: We recently introduced a multivariate approach that selects a subset of predictive genes jointly for sample classification based on expression data. We tested the algorithm on colon and leukemia data sets. As an extension to our earlier work, we systematically examine the sensitivity, reproducibility and stability of gene selection/sample classification to the choice of parameters of the algorithm. METHODS: Our approach combines a Genetic Algorithm (GA) and the k-Nearest Neighbor (KNN) method to identify genes that can jointly discriminate between different classes of samples (e.g. normal versus tumor). The GA/KNN method is a stochastic supervised pattern recognition method. The genes identified are subsequently used to classify independent test set samples. RESULTS: The GA/KNN method is capable of selecting a subset of predictive genes from a large noisy data set for sample classification. It is a multivariate approach that can capture the correlated structure in the data. We find that for a given data set gene selection is highly repeatable in independent runs using the GA/KNN method. In general, however, gene selection may be less robust than classification. AVAILABILITY: The method is available at http://dir.niehs.nih.gov/microarray/datamining CONTACT: LI3@niehs.nih.gov

Algorithms↗
PubMed2001

Biased distribution of inverted and direct Alus in the human genome: implications for insertion, exclusion, and genome stability.

Alu sequences, the most abundant class of large dispersed DNA repeats in human chromosomes, contribute to human genome dynamics. Recently we reported that long inverted repeats, including human Alus, can be strong initiators of genetic change in yeast. We proposed that the potential for interactions between adjacent, closely related Alus would influence their stability and this would be reflected in their distribution. We have undertaken an extensive computational analysis of all Alus (the database is at http://dir.niehs.nih.gov/ALU) to better understand their distribution and circumstances under which Alu sequences might affect genome stability. Alus separated by <650 bp were categorized according to orientation, length of regions sharing high sequence identity, distance between highly identical regions, and extent of sequence identity. Nearly 50% of all Alu pairs have long alignable regions (>275 bp), corresponding to nearly full-length Alus, regardless of orientation. There are dramatic differences in the distributions and character of Alu pairs with closely spaced, nearly identical regions. For Alu pairs that are directly repetitive, approximately 30% have highly identical regions separated by <20 bp, but only when the alignments correspond to near full-size or half-size Alus. The opposite is found for the distribution of inverted repeats: Alu pairs with aligned regions separated by <20 bp are rare. Furthermore, closely spaced direct and inverted Alus differ in their truncation patterns, suggesting differences in the mechanisms of insertion. At larger distances, the direct and inverted Alu pairs have similar distributions. We propose that sequence identity, orientation, and distance are important factors determining insertion of adjacent Alus, the frequency and spectrum of Alu-associated changes in the genome, and the contribution of Alu pairs to genome instability. Based on results in model systems and the present analysis, closely spaced inverted Alu pairs with long regions of alignment are likely at-risk motifs (ARMs) for genome instability.

Alu Elements↗
PubMed2001

Modeling human zymogen factor IX.

Modern theoretical techniques are employed to provide complete three dimensional structure for the zymogen and activated forms of human coagulation factors IX and IXa. These structures are fully calcium bound and equilibrated in an electrically neutral aqueous environment. The relationship of structure to mutational data is examined. We find that a substantial relative orientational change of the catalytic domain occurs on activation. Also, we find that the electrostatistically dipolar nature of the catalytic domain is substantially modified upon activation, with cleavage of the negatively charged activation peptide leaving behind a largely hydrophobic face in factor IXa. While the backbone atoms of the catalytic residues have little relative movement, nearby loops are found that do move. The presence or absence of these changes likely defines specificity.

Amino Acid Substitution↗
PubMed2000

Heparan/chondroitin sulfate biosynthesis. Structure and mechanism of human glucuronyltransferase I.

Human beta1,3-glucuronyltransferase I (GlcAT-I) is a central enzyme in the initial steps of proteoglycan synthesis. GlcAT-I transfers a glucuronic acid moiety from the uridine diphosphate-glucuronic acid (UDP-GlcUA) to the common linkage region trisaccharide Gal beta 1-3Gal beta 1-4Xyl covalently bound to a Ser residue at the glycosaminylglycan attachment site of proteoglycans. We have now determined the crystal structure of GlcAT-1 at 2.3 A in the presence of the donor substrate product UDP, the catalytic Mn(2+) ion, and the acceptor substrate analog Gal beta 1-3Gal beta 1-4Xyl. The enzyme is a alpha/beta protein with two subdomains that constitute the donor and acceptor substrate binding site. The active site residues lie in a cleft extending across both subdomains in which the trisaccharide molecule is oriented perpendicular to the UDP. Residues Glu(227), Asp(252), and Glu(281) dictate the binding orientation of the terminal Gal-2 moiety. Residue Glu(281) is in position to function as a catalytic base by deprotonating the incoming 3-hydroxyl group of the acceptor. The conserved DXD motif (Asp(194), Asp(195), Asp(196)) has direct interaction with the ribose of the UDP molecule as well as with the Mn(2+) ion. The key residues involved in substrate binding and catalysis are conserved in the glucuronyltransferase family as well as other glycosyltransferases.

Amino Acid Sequence↗
PubMed2000

Heparan sulfate biosynthesis: a theoretical study of the initial sulfation step by N-deacetylase/N-sulfotransferase.

Heparan sulfate N-deacetylase/N-sulfotransferase (NDST) catalyzes the deacetylation and sulfation of N-acetyl-D-glucosamine residues of heparan sulfate, a key step in its biosynthesis. Recent crystallographic and mutational studies have identified several potentially catalytic residues of the sulfotransferase domain of this enzyme (, J. Biol. Chem. 274:10673-10676). We have used the x-ray crystal structure of heparan sulfate N-sulfotransferase with 3'-phosphoadenosine 5'-phosphate to build a solution model with cofactor 3'-phosphoadenosine 5'-phosphosulfate (PAPS) and a model heparan sulfate ligand bound, and subsequently performed a 2-ns dynamics solution simulation. The simulation results confirm the importance of residues Glu(642), Lys(614), and Lys(833), with the possible involvement of Thr(617) and Thr(618), in binding PAPS. Additionally, Lys(676) is found in close proximity to the reaction site in our solvated structure. This study illustrates for the first time the possible involvement of water in the catalysis. Three water molecules were found in the binding site, where they are coordinated to PAPS, heparan sulfate, and the catalytic residues.

Amidohydrolases↗
PubMed2000

Modeling zymogen protein C.

A solution structure for the complete zymogen form of human coagulation protein C is modeled. The initial core structure is based on the x-ray crystallographic structure of the gamma-carboxyglutamic acid (Gla)-domainless activated form. The Gla domain (residues 1-48) is modeled from the x-ray crystal coordinates of the factor VII(a)/tissue factor complex and oriented with the epidermal growth factor-1 domain to yield an initial orientation consistent with the x-ray crystal structure of porcine factor IX(a). The missing C-terminal residues in the light chain (residues 147-157) and the activation peptide residues 158-169 were introduced using homology modeling so that the activation peptide residues directly interact with the residues in the calcium binding loop. Molecular dynamics simulations (Amber-particle-mesh-Ewald) are used to obtain the complete calcium-complexed solution structure. The individual domain structures of protein C in solution are largely unaffected by solvation, whereas the Gla-epidermal growth factor-1 orientation evolves to a form different from both factors VII(a) and IX(a). The solution structure of the zymogen protein C is compared with the crystal structures of the existing zymogen serine proteases: chymotrypsinogen, proproteinase, and prethrombin-2. Calculated electrostatic potential surfaces support the involvement of the serine protease calcium ion binding loop in providing a suitable electrostatic environment around the scissile bond for II(a)/thrombomodulin interaction.

1-Carboxyglutamic Acid↗
PubMed2000

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer: influence of the crystal environment.

Molecular dynamics (MD) simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure. Two crystal simulations were carried out; one consisted of one unit cell containing two duplexes, and the other of two unit cells containing four duplexes. Two solution simulations were also carried out, one starting from canonical B-DNA and the other starting from the crystal structure. For many helicoidal parameters, the results from the crystal and solution simulations were essentially identical. However, for other parameters, in particular, alpha, gamma, delta, (epsilon - zeta), phase, and helical twist, differences between crystal and solution simulations were apparent. Notably, during crystal simulations, values of helical twist remained comparable to those in the crystal structure, to include the sequence-dependent differences among base steps, in which values ranged from 20 degrees to 50 degrees per base step. However, in the solution simulations, not only did the average values of helical twist decrease to approximately 30 degrees per base step, but every base step was approximately 30 degrees, suggesting that the sequence-dependent information may be lost. This study reveals that MD simulations of the crystal environment complement solution simulations in validating the applicability of MD to the analysis of DNA structure.

Computer Simulation↗
PubMed1999

Residues in the alphaH and alphaI helices of the HIV-1 reverse transcriptase thumb subdomain required for the specificity of RNase H-catalyzed removal of the polypurine tract primer.

During retrovirus replication, reverse transcriptase (RT) must specifically interact with the polypurine tract (PPT) to generate and subsequently remove the RNA primer for plus-strand DNA synthesis. We have investigated the role that human immunodeficiency virus-1 RT residues in the alphaH and alphaI helices in the thumb subdomain play in specific RNase H cleavage at the 3'-end of the PPT; an in vitro assay modeling the primer removal step was used. Analysis of alanine-scanning mutants revealed that a subgroup exhibits an unusual phenotype in which the PPT is cleaved up to seven bases from its 3'-end. Further analysis of alphaH mutants (G262A, K263A, N265A, and W266A) with changes in residues in or near a structural motif known as the minor groove binding track showed that the RNase H activity of these mutants is more dramatically affected with PPT substrates than with non-PPT substrates. Vertical scan mutants at position 266 were all defective in specific RNase H cleavage, consistent with conservation of tryptophan at this position among lentiviral RTs. Our results indicate that residues in the thumb subdomain and the minor groove binding track in particular, are crucial for unique interactions between RT and the PPT required for correct positioning and precise RNase H cleavage.

Alanine↗
PubMed1999

Probing the structural changes in the light chain of human coagulation factor VIIa due to tissue factor association.

The crystallographic structure of human coagulation factor VIIa/tissue factor complex bound with calcium ions was used to model the solution structure of the light chain of factor VIIa (residues 1-142) in the absence of tissue factor. The Amber force field in conjunction with the particle mesh Ewald summation method to accommodate long-range electrostatic interactions was used in the trajectory calculations. The estimated TF-free solution structure was then compared with the crystal structure of factor VIIa/tissue factor complex to estimate the restructuring of factor VIIa due to tissue factor binding. The solution structure of the light chain of factor VIIa in the absence of tissue factor is predicted to be an extended domain structure similar to that of the tissue factor-bound crystal. Removal of the EGF1-bound calcium ion is shown by simulation to lead to minor structural changes within the EGF1 domain, but also leads to substantial relative reorientation of the Gla and EGF1 domains.

Amino Acid Sequence↗
PubMed1999

An atomic model for the pleated beta-sheet structure of Abeta amyloid protofilaments.

Synchrotron x-ray studies on amyloid fibrils have suggested that the stacked pleated beta-sheets are twisted so that a repeating unit of 24 beta-strands forms a helical turn around the fibril axis (. J. Mol. Biol. 273:729-739). Based on this morphological study, we have constructed an atomic model for the twisted pleated beta-sheet of human Abeta amyloid protofilament. In the model, 48 monomers of Abeta 12-42 stack (four per layer) to form a helical turn of beta-sheet. Each monomer is in an antiparallel beta-sheet conformation with a turn located at residues 25-28. Residues 17-21 and 31-36 form a hydrophobic core along the fibril axis. The hydrophobic core should play a critical role in initializing Abeta aggregation and in stabilizing the aggregates. The model was tested using molecular dynamics simulations in explicit aqueous solution, with the particle mesh Ewald (PME) method employed to accommodate long-range electrostatic forces. Based on the molecular dynamics simulations, we hypothesize that an isolated protofilament, if it exists, may not be twisted, as it appears to be when in the fibril environment. The twisted nature of the protofilaments in amyloid fibrils is likely the result of stabilizing packing interactions of the protofilaments. The model also provides a binding mode for Congo red on Abeta amyloid fibrils. The model may be useful for the design of Abeta aggregation inhibitors.

Amino Acid Sequence↗
PubMed1999

Molecular dynamics simulations of biomolecules: long-range electrostatic effects.

Current computer simulations of biomolecules typically make use of classical molecular dynamics methods, as a very large number (tens to hundreds of thousands) of atoms are involved over timescales of many nanoseconds. The methodology for treating short-range bonded and van der Waals interactions has matured. However, long-range electrostatic interactions still represent a bottleneck in simulations. In this article, we introduce the basic issues for an accurate representation of the relevant electrostatic interactions. In spite of the huge computational time demanded by most biomolecular systems, it is no longer necessary to resort to uncontrolled approximations such as the use of cutoffs. In particular, we discuss the Ewald summation methods, the fast particle mesh methods, and the fast multipole methods. We also review recent efforts to understand the role of boundary conditions in systems with long-range interactions, and conclude with a short perspective on future trends.

Computer Simulation↗
PubMed1998

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↗
PubMed1998

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↗
PubMed1997

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↗
PubMed1997

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↗