PubMed HealthSearch

Biomedical subjects

L G Pedersen

Publications and source records attributed to L G Pedersen.

At least 19 recordsLinked to original sources

The roles of individual gamma-carboxyglutamate residues in the solution structure and cation-dependent properties of conantokin-T.

The solution structure of the Ca2+-loaded conantokin-T (con-T), a gamma-carboxyglutamate (Gla)-containing 21-residue peptide (NH2-G1EgammagammaY5QKMLgamma10NLRgammaA15EVKKN20A-CONH2,gam ma = Gla), has been elucidated by use of distance geometry calculations with experimental distances derived from two-dimensional 1H NMR spectroscopy. An end-to-end alpha-helix was the dominant conformation in solution, similar to that of apo-con-T, except that reorientation of several side chains occurred in the Ca2+-coordinated complex. The most notable examples of this were those of Gla10 and Gla14, which were more optimally positioned for complexation with Ca2+. In addition to the stabilization offered to the alpha-helix by Ca2+ binding, hydrophobic clustering of the side chains of Tyr5, Met8, Leu9, and Leu12, and ionic interactions between Lys7 and Gla3/Gla10 and between Arg13 and Gla14, along with hydrogen bonding between Gln6 and Gla10, were among the side chain interactions likely playing a significant role in maintenance of the alpha-helical conformation. Docking of Ca2+ in the con-T structure was accomplished using genetic algorithm-molecular dynamics simulation approaches. The results showed that one Ca2+ ion is most likely coordinated by four side chain oxygen atoms, two each from Gla10 and Gla14. Another bound Ca2+ ion has as its donor sites three oxygen atoms, two from Gla3 and one from Gln6. To examine the functional roles of the individual Gla residues, a series of variant peptides have been synthesized with Ala substituted for each Gla residue, and several properties of the resulting variants have been examined. The data obtained demonstrated the importance of Gla10 and Gla14 in stabilizing binding of the highest affinity Ca2+ site and in governing the conformational change induced by Ca2+. The critical nature of Gla3 and Gla4 in inhibition of the spermine-induced potentiation of the binding of MK-801 to open ion channels of the N-methyl-D-aspartate receptor was established, as well as the role of Gla4 in stabilizing the apo-con-T alpha-helical conformation.

1-Carboxyglutamic Acid

Role of gamma-carboxyglutamic acid in the calcium-induced structural transition of conantokin G, a conotoxin from the marine snail Conus geographus.

Conantokin G is a gamma-carboxyglutamic acid- (Gla-) containing conotoxin isolated from the venom of the marine cone snail Conus geographus. This 17-residue polypeptide, which contains five gamma-carboxyglutamic acid residues, is a N-methyl-d-aspartate- (NMDA-) type glutamate receptor antagonist. To investigate the role of gamma-carboxyglutamic acid in the calcium-induced structural transition of conantokin G, we determined the three-dimensional structure of the conantokin G/Ca2+ complex by two-dimensional 1H NMR spectroscopy and compared it to the high-resolution structure of conantokin G in the absence of metal ions [Rigby et al. (1997) Biochemistry 36, 6906]. Complete resonance assignments were made by two dimensional 1H NMR spectroscopy at pH 5.6 in the presence of saturating amounts of Ca2+. Distance geometry and simulated annealing methods were used to derive 23 convergent structures from a set of 302 interproton distance restraints and two torsion angle measurements. A high-resolution structure, with the backbone root mean square deviation to the geometric average of the 23 structures of 0.6 +/- 0.1 A, contains a linear alpha-helix from Gla 3 to Lys 15. Gla residues 3, 7, 10, and 14 are aligned in a linear array on one face of the helix. A genetic algorithm was applied to determine the calcium positions in conantokin G, and the conantokin G/Ca2+ complex refined by molecular simulation. Upon binding of Ca2+ to gamma-carboxyglutamic acid, conantokin G undergoes a conformational transition from a distorted curvilinear 310 helix to a linear alpha-helix. Occupancy of the metal binding sites, defined by gamma-carboxyglutamic acids, results in formation of a calcium-carboxylate network that linearizes the helix and exposes the hydrophobic amino acids on the opposite face of the helix.

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

Prediction of solution structures of the Ca2+-bound gamma-carboxyglutamic acid domains of protein S and homolog growth arrest specific protein 6: use of the particle mesh Ewald method.

The solution structures of the N-terminal domains of protein S, a plasma vitamin K-dependent glycoprotein, and its homolog growth arrest specific protein 6 (Gas6) were predicted by molecular dynamics computer simulations. The initial structures were based on the x-ray crystallographic structure of the corresponding region of bovine prothrombin fragment 1. The subsequent molecular dynamics trajectories were calculated using the second-generation AMBER force field. The long-range electrostatic forces were evaluated by the particle mesh Ewald method. The structures that stabilized over a 400-ps time interval were compared with the corresponding region of the simulated solution structure of bovine prothrombin fragment 1. Structural properties of the gamma-carboxyglutamic acid (Gla) domains obtained from simulations and calcium binding were found to be conserved for all three proteins. Analysis of the predicted solution structure of the Gla domain of Gas6 suggests that this domain should bind with negatively charged phospholipid surfaces analogous to bovine prothrombin fragment 1 and protein S.

1-Carboxyglutamic Acid

Crystal structure of estrogen sulphotransferase.

The structure of estrogen sulphotransferase has been solved in the presence of inactive cofactor PAP and substrate 17 beta-estradiol. This structure reveals structural similarities between cytosolic sulphotransferases and nucleotide kinases.

Amino Acid Sequence

Studies on the peptide corresponding to residues 34-47 of bovine factor X.

Calcium binding studies of a 14-residue peptide corresponding to the 37-46 sequence of bovine factor X were performed using calcium ion selective electrode titrations and equilibrium dialysis. The presence of gamma-carboxyglutamic acid residues at positions 36 and 40 coupled with the assumption that the peptide would bind calcium ions also prompted an investigation of possibly secondary conformational changes in the peptide by use of circular dichroism spectroscopy. Equilibrium dialysis revealed a single relatively weak calcium binding site (log Ka = 2.39); an ion selective electrode experiment confirmed this result (log Ka = 2.17). The peptide maintained a random coil conformation throughout the calcium ion titrations as measured by circular dichroism.

Amino Acid Sequence

Characterization of gamma-carboxyglutamic acid residue 21 of human factor IX.

We investigated the functional role of gamma-carboxyglutamic acid (Gla) residue 21 of human factor IX, using site-directed mutagenesis to change the glutamic acid residue to aspartic acid (FIX21D). FIX21D had reduced activity in an activated partial thromboplastin time (aPTT) assay and was activated by factor XIa more slowly than wild-type factor IX (FIXwt). FIX21D underwent normal, two-stage calcium-dependent intrinsic fluorescence quenching, indicating that a folding event similar to that seen in FIXwt occurred upon the addition of calcium ions. Antibody A-7, which recognizes factor IX-specific residues at positions 33-40, bound FIX21D as well as FIXwt; however, the calcium-specific monoclonal antibody, JK-IX-2, whose epitope includes residues 1 and 22, did not recognize FIX21D. FIX21D bound phosphatidylserine/phosphatidylcholine (PS/PC) vesicles with Kd approximately 10-fold greater than FIXwt, as measured by a fluorescence light scattering assay. Finally, although FIXwt binds endothelial cells with a Kd of 2.8 nM, FIX21D did not bind endothelial cells. Molecular modeling simulations of FIXwt and FIX21D indicate that mutating Gla 21 to Asp causes structural changes in residues 3-5 and 8-10, as well as in two exposed calcium ions, consistent with the reduced function of FIX21D. Immunological and intrinsic fluorescence quenching assays and the molecular dynamics simulations suggest normal folding in the C-terminal region of the Gla domain. Thus we hypothesize the FIX21D has reduced JK-IX-2 and phospholipid and endothelial cell binding due to localized structural changes in residues 3-10 and the exposed calcium ions. Our study suggests that the Gla 21 to Asp mutation disrupts function in the N-terminal region of the Gla domain without affecting structure in the C-terminal Gla domain region.

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

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

Computational studies of human prothrombin fragment 1, the Gla domain of factor IX and several biological interesting mutants.

We have employed molecular dynamics simulations to understand the properties of the gamma-carboxyglutamic acid (Gla) domains of human prothrombin fragment 1 (residues 1-144) and factor IX (residues 1-47). The simulations are based on (1) the crystal structure of bovine prothrombin fragment 1 in the presence of calcium and (2) our subsequent simulation of the solution structure for which we found accommodation of the long range ionic forces critical. In addition, we have estimated the solution structures for key mutant proteins [prothrombin (Gla6 to Asp and Gla16 to Asp) and factor IX (Gly12 to Ala)]. The simulations for the latter two mutants do not stabilize, a result in concert with the known biological data.

1-Carboxyglutamic Acid

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

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

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