Search PubMed⌕ Search

Biomedical subjects

L G Pedersen

Publications and source records attributed to L G Pedersen.

At least 55 records · Page 3Linked to original sources

The first solvation shell of magnesium ion in a model protein environment with formate, water, and X-NH3, H2S, imidazole, formaldehyde, and chloride as ligands: an Ab initio study.

The first coordination shell of an Mg(II) ion in a model protein environment is studied. Complexes containing a model carboxylate, an Mg(II) ion, various ligands (NH3, H2S, imidazole, and formaldehyde) and water of hydration about the divalent metal ion were geometry optimized. We find that for complexes with the same coordination number, the unidentate carboxylate-Mg(II) ion is greater than 10 kcal mol-1 more stable than the bidentate orientation. Imidazole was found to be the most stable ligand, followed in order by NH3, formaldehyde, H2O, and H2S.

Cations↗

The first solvation shell of magnesium and calcium ions in a model nucleic acid environment: an ab initio study.

The interaction of organophosphate anions with divalent metal ions is central to many biological catalytic events. While experimental structural studies can give insight into the likely geometries that can be adopted, quantum mechanics allows for a more complete exploration of the competing forms. Ab initio quantum mechanical calculations have been performed on a series of complexes comprised of dimethyl phosphate, a divalent metal ion (either Mg(II) or Ca(II)) and water of hydration. An additional series of complexes were studied that included a Cl(I) ion to provide for charge neutrality. The most stable orientation of the hydrated metal ion complexed with the phosphate anion occurs when the metal ion is in a unidentate, rather than bidentate, orientation. The question of whether the divalent metal ion is located in the phosphinyl (-PO2(-)-) plane depends on the identity of the divalent metal ion and on the charge state of the complex.

Binding Sites↗

Is Ca(II) ion binding to prothrombin fragment 1 intrinsically cooperative, or is the cooperative binding accounted for by self-association?

The recent suggestion that the apparent cooperativity seen in the binding of Ca(II) ions to prothrombin fragment 1 is due to protein aggregation is evaluated. Since (1) we find that the Ca(II) ion binding is not dependent upon protein concentration, (2) the analytical expression for the equilibrium constant of the aggregation model is unrealistically large when evaluated at realistic Ca(II) ion concentrations, and (3) a very simple allosteric cooperative binding model (Monod) can be shown to fit the experimental data, we conclude that the aggregation explanation for the apparent cooperativity in the Ca(II) ion binding by prothrombin fragment 1 is not correct.

Animals↗

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↗

Simulation of the solution structure of the H-ras p21-GTP complex.

An unconstrained simulation of the GTP-bound form of the H-ras protein p21 is performed in an aqueous environment with charge-neutralizing counterions. The simulation is compared to the 1.35-A structure of Pai et al. [(1990) EMBO J. 9, 2351] and a proposed alternate structure, in which the loop at residues 60-65 is modeled into a form which may activate a water molecule for the GTP hydrolysis. The simulation suggests that some protein intermolecular H-bond contacts which are present in the crystal structure are lost in the solvation process and this loss may lead to localized refolding of the molecule. For instance, we find that the gamma-phosphate of the GTP has somewhat weaker contact with the protein in the simulation structure. The antiparallel beta-sheet (residues 38-57) partially melts. The 60-65 loop, which is hypervariable in the X-ray study, is initially relatively distant from the gamma-phosphate region. However, this loop moves so as to sample the space around the gamma-phosphate. For a significant fraction of the simulation time, forms similar to the alternate structure are observed, and a water molecule is localized near the hydrolytic site. The molecular dynamics simulations of p21-GTP in solution support a postulated hydrolysis mechanism for the biological inactivation of the nucleotide complex based on crystallographic data.

Amino Acid Sequence↗

A metal ion-binding site in the kringle region of bovine prothrombin fragment 1.

45Ca(II) binding studies (equilibrium dialysis) on the kringle domain of bovine prothrombin fragment 1 were conducted using a mixture of peptides (residues 43-156 and 46-156) resulting from limited alpha-chymotryptic hydrolysis of fragment 1. Analysis of the Scatchard plot of these data indicates a single, low affinity Ca(II)-binding site to be present. Similar results were obtained from studies on the decarboxylated fragment 1 derivative, 10-gamma-MGlu-fragment 1. Acetylation of bovine fragment 1 in the absence of Ca(II) or Mg(II) ions results in the loss of the metal ion-promoted quenching of the intrinsic Trp fluorescence of the protein and the Ca(II)-mediated binding to phosphatidylserine/phosphatidylcholine (PS/PC) vesicles. The acetylation of the NH2 alpha-group of Ala-1 has been shown (Welsch, D. J., and Nelsestuen, G. L. (1988) Biochemistry 27, 4946-4952) to abolish the PS/PC binding property of fragment 1. The present study demonstrates that acetylation of a second site possibly Ser-79 or Thr-81 using the conditions described in the preceding paper results in loss of both the fluorescence transition and the Ca(II)-mediated PS/PC binding of the resulting protein derivative. Removal of the O-acetyl group at the Ser-79/Thr-81 site is accomplished by aminolysis with 0.2 M hydroxylamine, pH 10, 50 degrees C; the fluorescence transition is partially restored. PS/PC binding is partially restored if the NH2 alpha-group of Ala-1 is trinitrophenylated but is not restored if the NH2 alpha-group of Ala-1 is acetylated. We conclude that the Ser-79/Thr-81 site may represent a portion of the metal ion-binding site within the kringle domain of fragment 1. Occupancy of this site by a Ca(II) ion appears to be important in the binding of the protein to PS/PC vesicles.

Animals↗

Modifications of bovine prothrombin fragment 1 in the presence and absence of Ca(II) ions. Loss of positive cooperativity in Ca(II) ion binding for the modified proteins.

Chemical modification of bovine prothrombin fragment 1 according to the procedure of D. J. Welsch and G. L. Nelsestuen (1988) [Biochemistry 27, 4946-4952 and ealier papers] provided a series of fragment 1 derivatives in which various nitrogen-containing side chains were N-acetylated and/or N-2,4,6-trinitrophenylated. In addition the des-[Ala-1,Asn-2]- and des-[Ala-1,Asn-2,Lys-3]-fragment 1 derivatives were prepared by limited enzymatic hydrolysis of fragment 1 using cathepsin C and plasmin, respectively. Quantitative studies on the Ca(II) binding of these proteins have been accomplished using 45Ca(II) equilibrium dialysis. Binding of these fragment 1 derivatives to phosphatidylserine/phosphatidylcholine (PS/PC) vesicles (25:75) in the presence of Ca(II) ions has been studied using the light-scattering technique. Acylation of the 5 lysine residues of fragment 1 by the action of acetic anhydride (500-fold molar excess) in the presence of 75 mM Ca(II), pH 8.0, results in loss of positive cooperativity in Ca(II) binding (Scatchard plot) and an increase in the number of Ca(II) ions bound. The Ca(II)-dependent PS/PC binding of the acylated protein is reduced. Removal of 2 and 3 residues from the amino terminus likewise leads to loss of positive cooperativity in Ca(II) binding and reduced binding affinity to PS/PC vesicles. The important role of the amino-terminal 1-10 sequence is discussed. We conclude that positive cooperativity in Ca(II) binding is not a prerequisite for the Ca(II)-dependent binding of bovine prothrombin fragment 1 to PS/PC vesicles.

Animals↗

Predicted secondary structure of bovine prothrombin fragment 1 and related proteins in different environments by circular dichroism spectroscopy.

Circular dichroism spectroscopy was used to investigate the structure of bovine prothrombin fragment 1 (BF1) and related proteins in several environments. The conformational change induced in BF1 by the addition of Mg[II] ions was found to be different from that induced by Ca[II] or Sr[II]. The Ca[II] and Sr[II] conformations appear to differ only slightly from the apo-metal conformation. The conformation of the 1-45 fragment of prothrombin, however, is markedly different than the conformation of the same fragment in the presence of either Ca[II] of Mg[II]; both of the latter structures differ substantially from one another. The presence of phospholipids has almost no effect on the structure of either BF1 or the 1-45 fragment; in the presence of both phospholipids and Ca[II] a structural change is seen for the 1-45 fragment but not BF1 (relative to the protein alone). The addition of phospholipids to the Mg[II]/BF1 structure did not induce a CD-detectable conformational change, while the addition of phospholipids to the Ca[II]/BF1 or Sr[II]/BF1 structures induced a change to a conformation similar in secondary structure composition to the relative apometal structures.

Animals↗

Determination of strontium binding to macromolecules.

An equilibrium dialysis technique for determining the binding of strontium to macromolecules is described. The major difficulty to be overcome is that 90Sr has a decay product, 90Y, which is also a beta-emitter. The described protocol is used to determine the Sr binding isotherm to bovine prothrombin fragment 1. The binding is found to be cooperative, somewhat weaker than Ca binding, and to involve approximately nine strontium sites. The stoichiometric equilibrium constants are determined by nonlinear regression. The procedure should be of great utility for many macromolecules that show strontium affinity.

Animals↗

Expression and characterization of human factor IX. Factor IXthr-397 and factor IXval-397.

Factor IXLong Beach has a single amino acid substitution at 397 (Ile to Thr) in the catalytic domain which results in severe hemophilia B. Recent investigations have shown that the substitution of threonine for isoleucine at 397 may affect a part of the macromolecular substrate binding site. Because threonine has a hydroxyl group in its side chain, it is possible that this hydroxyl group makes new hydrogen bonds and disturbs the substrate binding site. We used three techniques: molecular biology, which includes site-directed mutagenesis and recombinant protein expression in tissue culture; computer-aided kinetic data analysis; and molecular modeling to study this mutation site. We have produced two mutant factor IX molecules that have isoleucine 397 replaced by valine or threonine. Factor IXwild type and the two mutants (factor IXVal and factor IXThr) were expressed in human kidney cells and purified using a conformation-specific monoclonal antibody column. After the activation by factor XIa, these three molecules were able to bind p-aminobenzamidine and increase its fluorescence intensity in a similar manner. Factor IXVal and factor IXwild type had indistinguishable activities in an activated partial thromboplastin time (aPTT) assay and similar kinetic parameters with factor X as a substrate. Factor IXThr had only 5% clotting activity compared with normal factor IX, a slightly lower Km and significantly reduced kcat, using factor X as a substrate. We developed energy-refined (AMBER v.3.1) computer models of the three factor IX molecules based on previous work. Three factor IXa models (Ile, Val, or Thr at 397) with a fragment of the factor X activation site were used to predict the effect of the mutation at 397 and evaluate the significance of the new hydrogen bond thought to form between the side chain hydroxyl group of threonine 397 and the carbonyl oxygen of tryptophan 385. This new hydrogen bond would affect the position of an amide proton of adjacent glycine 386 which has been proposed to make a hydrogen bond with a backbone carbonyl oxygen of the P3 residue of factor X. In addition to the new hydrogen bond, there is significant movement in the side chain of tryptophan 385 between the factor IXawild type-factor X model and the factor IXaThr-factor X model that could interfere with substrate binding. This movement could be caused by the change in the molecular volume, the orientation of the side chain at 397, and the new hydrogen bond.

Base Sequence↗

Solution conformations of the gamma-carboxyglutamic acid domain of bovine prothrombin fragment 1, residues 1-65.

Molecular dynamics simulations have been performed (AMBER version 3.1) on solvated residues 1-65 of bovine prothrombin fragment 1 (BF1) by using the 2.8-A resolution crystallographic coordinates as the starting conformation for understanding calcium ion-induced conformational changes that precede experimentally observable phospholipid binding. Simulations were performed on the non-metal-bound crystal structure, the form resulting from addition of eight calcium ions to the 1-65 region of the crystal structure, the form resulting from removal of calcium ions after 107 ps and continuing the simulation, and an isolated hexapeptide loop (residues 18-23). In all cases, the 100-ps time scale seemed adequate to sample an ensemble of solution conformers within a particular region of conformation space. The non-metal-containing BF1 did not unfold appreciably during a 106-ps simulation starting from the crystallographic geometry. The calcium ion-containing structure (Ca-BF1) underwent an interesting conformational reorganization during its evolution from the crystal structure: during the time course of a 107-ps simulation, Ca-BF1 experienced a trans----cis isomerization of the gamma-carboxyglutamic acid-21 (Gla-21)-Pro-22 peptide bond. Removal of the calcium ions from this structure followed by 114 ps of additional molecular dynamics showed significant unfolding relative to the final 20-ps average structure of the 107-ps simulation; however, the Gla-21-Pro-22 peptide bond remained cis. A 265-ps simulation on the termini-protected hexapeptide loop (Cys-18 to Cys-23) containing two calcium ions also did not undergo a trans----cis isomerization. It is believed that the necessary activation energy for the transitional event observed in the Ca-BF1 simulation was largely supplied by global conformational events with a possible assist from relief of intermolecular crystal packing forces. The presence of a Gla preceding Pro-22, the inclusion of Pro-22 in a highly strained loop structure, and the formation of two long-lived salt bridges prior to isomerization may all contribute to this finding.

1-Carboxyglutamic Acid↗

Theoretical and experimental measures of DNA helix stability and their relation to sequence specific repair of O6-ethylguanine lesions.

Recent work (Breslauer et al. (1986) Proc. Natl. Acad. Sci. (U.S.A.), 83, 3746) has provided a method for calculating empirical thermodynamic quantities for helix to coil transitions from the base sequence of any oligomer. It is shown in this work that the DNA helix binding energy, calculated with the AMBER force field, for 9-mers of the type 5'-GGGXGeYGGG-3', where X and Y are any base and the central Ge is O6-ethylguanine, correlates well with the empirical delta G for helix to strand transitions. The mutation spectrum of ethane methylsulfonate (EMS) in the lacI gene of Escherichia coli can be modeled using the calculated local binding energy but the empirical free energies, enthalpies and melting temperatures predict these levels of repair less well. The relation of the binding energy to the mutation spectrum can be somewhat improved by including entropic effects in a theoretical free energy of binding as given by delta G theoretical identical to delta E binding - T delta S.

Base Sequence↗