Biomedical subjects

Jan H Jensen

Publications and source records attributed to Jan H Jensen.

17 recordsLinked to original sources

PubMed2006

Exploring the role of the active site cysteine in human muscle creatine kinase.

All known guanidino kinases contain a conserved cysteine residue that interacts with the non-nucleophilic eta1-nitrogen of the guanidino substrate. Site-directed mutagenesis studies have shown that this cysteine is important, but not essential for activity. In human muscle creatine kinase (HMCK) this residue, Cys283, forms part of a conserved cysteine-proline-serine (CPS) motif and has a pKa about 3 pH units below that of a regular cysteine residue. Here we employ a computational approach to predict the contribution of residues in this motif to the unusually low cysteine pKa. We calculate that hydrogen bonds to the hydroxyl and to the backbone amide of Ser285 would both contribute approximately 1 pH unit, while the presence of Pro284 in the motif lowers the pKa of Cys283 by a further 1.2 pH units. Using UV difference spectroscopy the pKa of the active site cysteine in WT HMCK and in the P284A, S285A, and C283S/S285C mutants was determined experimentally. The pKa values, although consistently about 0.5 pH unit lower, were in broad agreement with those predicted. The effect of each of these mutations on the pH-rate profile was also examined. The results show conclusively that, contrary to a previous report (Wang et al. (2001) Biochemistry 40, 11698-11705), Cys283 is not responsible for the pKa of 5.4 observed in the WT V/K(creatine) pH profile. Finally we use molecular dynamics simulations to demonstrate that, in order to maintain the linear alignment necessary for associative inline transfer of a phosphoryl group, Cys283 needs to be ionized.

Adenosine Diphosphate↗

PubMed2006

Cooperative hydrogen bonding effects are key determinants of backbone amide proton chemical shifts in proteins.

A computational methodology for backbone amide proton chemical shift (delta(H)) predictions based on ab initio quantum mechanical treatment of part of the protein is presented. The method is used to predict and interpret 13 delta(H) values in protein G and ubiquitin. The predicted amide-amide delta(H) values are within 0.6 ppm of experiment, with a root-mean-square deviation (RMSD) of 0.3 ppm. We show that while the hydrogen bond geometry is the most important delta(H)-determinant, longer-range cooperative effects of extended hydrogen networks make significant contributions to delta(H). We present a simple model that accurately relates the protein structure to delta(H).

Amino Acid Sequence↗

PubMed2006

FTIR spectroscopy combined with isotope labeling and quantum chemical calculations to investigate adsorbed bicarbonate formation following reaction of carbon dioxide with surface hydroxyl groups on Fe2O3 and Al2O3.

FTIR spectroscopy combined with isotope labeling experiments and quantum chemical calculations is used to investigate the adsorption of carbon dioxide on hydroxylated metal oxide surfaces. In particular, transmission FTIR spectra following CO2 adsorption on hydroxylated nanoparticulate Fe2O3, alpha-Al2O3, and gamma-Al2O3 particles at 296 K are reported. As expected, reaction of CO2 with these surfaces results in the formation of adsorbed bicarbonate and carbonate. In this study, the vibrational spectrum of the bicarbonate product is analyzed in detail through the use of isotope labeling experiments and quantum chemical calculations. The experimental and calculated vibrational frequencies of adsorbed HC16O3-, DC16O3-, HC18O3-, HC16O18O2-, and HC18O16O2- indicate that bicarbonate bonds to the surface in a bridged structure. There is some evidence from the mixed isotope experiments that following initial nucleophilic attack of OH, the formation of the final bicarbonate structure involves a proton transfer. On the basis of energetic considerations, the proton transfer mechanism most likely occurs through an intermolecular process involving either coadsorbed hydroxyl groups and/or carbonate.

Journal Article↗

PubMed2006

Charge transfer interaction in the effective fragment potential method.

An approximate formula is derived and implemented in the general effective fragment potential (EFP2) method to model the intermolecular charge transfer interaction. This formula is based on second order intermolecular perturbation theory and utilizes canonical molecular orbitals and Fock matrices obtained with preparative self-consistent field calculations. It predicts charge transfer energies that are in reasonable agreement with the reduced variational space energy decomposition analysis. The formulas for the charge transfer gradients with respect to EFP translational and rotational displacements are also derived and implemented.

Chemistry, Physical↗

PubMed2006

Chemically accurate protein structures: validation of protein NMR structures by comparison of measured and predicted pKa values.

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pKa values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program ( http://www.propka.chem.uiowa.edu) allows the user to test the ability of a given structure to reproduce known pKa values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both NMR and X-ray structures, as well as pKa values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.

Aspartic Acid↗

PubMed2006

Hydrogen bonding is the prime determinant of carboxyl pKa values at the N-termini of alpha-helices.

Experimentally determined mean pK(a) values of carboxyl residues located at the N-termini of alpha-helices are lower than their overall mean values. Here, we perform three types of analyses to account for this phenomenon. We estimate the magnitude of the helix macrodipole to determine its potential role in lowering carboxyl pK(a) values at the N-termini. No correlation between the magnitude of the macrodipole and the pK(a) values is observed. Using the pK(a) program propKa we compare the molecular surroundings of 18 N-termini carboxyl residues versus 233 protein carboxyl groups from a previously studied database. Although pK(a) lowering interactions at the N-termini are similar in nature to those encountered in other protein regions, pK(a) lowering backbone and side-chain hydrogen bonds appear in greater number at the N-termini. For both Asp and Glu, there are about 0.5 more hydrogen bonds per residue at the N-termini than in other protein regions, which can be used to explain their lower than average pK(a) values. Using a QM-based pK(a) prediction model, we investigate the chemical environment of the two lowest Asp and the two lowest Glu pK(a) values at the N-termini so as to quantify the effect of various pK(a) determinants. We show that local interactions suffice to account for the acidity of carboxyl residues at the N-termini. The effect of the helix dipole on carboxyl pK(a) values, if any, is marginal. Backbone amide hydrogen bonds constitute the single biggest contributor to the lowest carboxyl pK(a) values at the N-termini. Their estimated pK(a) lowering effects range from about 1.0 to 1.9 pK(a) units.

Carbon Dioxide↗

PubMed2006

Hydride transfer versus hydrogen radical transfer in thymidylate synthase.

The nature of a H-transfer in the thymidylate synthase catalyzed reaction was investigated by comparison of the wild-type enzyme with the W80M mutant. The nature of the H-transfer was not affected, as indicated by intrinsic isotope effects and their temperature dependence. These findings support a single-step hydride transfer instead of a two-step radical transfer.

Hydrogen↗

PubMed2006

The polarizable continuum model (PCM) interfaced with the fragment molecular orbital method (FMO).

The polarizable continuum model (PCM) for the description of solvent effects is combined with the fragment molecular orbital (FMO) method at several levels of theory, using a many-body expansion of the electron density and the corresponding electrostatic potential, thereby determining solute (FMO)-solvent (PCM) interactions. The resulting method, denoted FMO/PCM, is applied to a set of model systems, including alpha-helices and beta-strands of alanine consisting of 10, 20, and 40 residues and their mutants to charged arginine and glutamate residues. The FMO/PCM error in reproducing the PCM solvation energy for a full system is found to be below 1 kcal/mol in all cases if a two-body expansion of the electron density is used in the PCM potential calculation and two residues are assigned to each fragment. The scaling of the FMO/PCM method is demonstrated to be nearly linear at all levels for polyalanine systems. A study of the relative stabilities of alpha-helices and beta-strands is performed, and the magnitude of the contributing factors is determined. The method is applied to three proteins consisting of 20, 129, and 245 residues, and the solvation energy and computational efficiency are discussed.

Journal Article↗

PubMed2005

Very fast empirical prediction and rationalization of protein pKa values.

A very fast empirical method is presented for structure-based protein pKa prediction and rationalization. The desolvation effects and intra-protein interactions, which cause variations in pKa values of protein ionizable groups, are empirically related to the positions and chemical nature of the groups proximate to the pKa sites. A computer program is written to automatically predict pKa values based on these empirical relationships within a couple of seconds. Unusual pKa values at buried active sites, which are among the most interesting protein pKa values, are predicted very well with the empirical method. A test on 233 carboxyl, 12 cysteine, 45 histidine, and 24 lysine pKa values in various proteins shows a root-mean-square deviation (RMSD) of 0.89 from experimental values. Removal of the 29 pKa values that are upper or lower limits results in an RMSD = 0.79 for the remaining 285 pKa values.

Amino Acids↗

PubMed2005

Prediction and rationalization of protein pKa values using QM and QM/MM methods.

We describe the development and application of a computational method for the prediction and rationalization of pKa values of ionizable residues in proteins, based on ab initio quantum mechanics (QM) and the effective fragment potential (EFPs) method (a hybrid QM/MM method). The theoretical developments include (1) a covalent boundary method based on frozen localized orbitals, (2) divide-and-conquer methods for the ab initio computation of protein EFPs consisting of multipoles up to octupoles and dipole polarizability tensors, (3) a method for computing vibrational free energies for a localized molecular region, and (4) solutions of the polarized continuum model of bulk solvation equations for protein-sized systems. The QM-based pKa prediction method is one of the most accurate methods currently available and can be used in cases where other pKa prediction methods fail. Preliminary analysis of the computed results indicate that many pKa values (1) are primarily determined by hydrogen bonds rather than long-range charge-charge interactions and (2) are relatively insensitive to large-scale dynamical fluctuations of the protein structure.

Algorithms↗

PubMed2004

Determinants of cysteine pKa values in creatine kinase and alpha1-antitrypsin.

The structural determinants of the unusually low pK(a) values of Cys282 in human creatine kinase and Cys232 in alpha1-antitrypsin were studied computationally. We have demonstrated that hydrogen bonding to the cysteine residue is the prime determinant for both proteins. In the case of creatine kinase, the hydrogen bond donors are a serine side chain and an amide NH-group, while in alpha1-antitrypsin the donor is an amide NH. Each hydrogen bond lowers the pK(a) by between 0.8 and 1.5 pH units. The 1.1-unit lowering due to the Ser284-Cys282 hydrogen bond is in good agreement with the 1.2-unit difference between the Cys282 pK(a) value of wild-type and the S284A mutant of creatine kinase.

Creatine Kinase↗

PubMed2004

Determinants of the relative reduction potentials of type-1 copper sites in proteins.

The relative Cu(2+)/Cu(+) reduction potentials of six type-1 copper sites (cucumber stellacyanin, P. aeruginosa azurin, poplar plastocyanin, C. cinereus laccase, T. ferrooxidans rusticyanin, and human ceruloplasmin), which lie in a reduction potential range from 260 mV to over 1000 mV, have been studied by quantum mechanical calculations. The range and relative orderings of the reduction potentials are reproduced very well compared to experimental values. The study suggests that the main structural determinants of the relative reduction potentials of the blue copper sites are located within 6 A of the Cu atoms. Further analysis suggests that the reduction potential differences of type-1 copper sites are caused by axial ligand interactions, hydrogen bonding to the S(Cys), and protein constraint on the inner sphere ligand orientations. The low reduction potential of cucumber stellacyanin is due mainly to a glutamine ligand at the axial position, rather than a methionine or a hydrophobic residue as in the other proteins. A stronger interaction with a backbone carbonyl group is a prime contributor to the lower reduction potential of P. aeruginosa azurin as compared to poplar plastocyanin, whereas the reverse is true for C. cinereus laccase and T. ferrooxidans rusticyanin. The lack of an axial methonine ligand also contributes significantly to the increased reduction potentials of C. cinereus laccase and human ceruloplasmin. However, in the case of C. cinereus laccase, this increase is attenuated by the presence of only one amide NH hydrogen bond to the S(Cys) rather than two in the other proteins. In human ceruloplasmin the reduction potential is further increased by the structural distortion of the equatorial ligand orientation.

Azurin↗

PubMed2004

The determinants of carboxyl pKa values in turkey ovomucoid third domain.

A computational methodology for protein pK(a) predictions, based on ab initio quantum mechanical treatment of part of the protein and linear Poisson-Boltzmann equation treatment of the bulk solvent, is presented. The method is used to predict and interpret the pK(a) values of the five carboxyl residues (Asp7, Glu10, Glu19, Asp27, and Glu43) in the serine protease inhibitor turkey ovomucoid third domain. All the predicted pK(a) values are within 0.5 pH units of experiment, with a root-mean-square deviation of 0.31 pH units. We show that the decreased pK(a) values observed for some of the residues are primarily due to hydrogen bonds to the carboxyl oxygens. Hydrogen bonds involving amide protons are shown to be particularly important, and the effect of hydrogen bonding is shown to be nonadditive. Hydrophobic effects are also shown to be important in raising the pK(a). Interactions with charged residues are shown to have relatively little effect on the carboxyl pK(a) values in this protein, in general agreement with experiment.

Animals↗

PubMed2004

Improving the efficiency and convergence of geometry optimization with the polarizable continuum model: new energy gradients and molecular surface tessellation.

New equations are derived and implemented for efficient and accurate computation of solvation energy derivatives for the conductor-like polarizable continuum model (C-PCM) and the isotropic integral equation formalism polarizable continuum model (IEF-PCM). Two new molecular surface tessellation procedures GEPOL-RT and GEPOL-AS that generate near continuous potential energy surfaces are proposed for PCM geometry optimization. The combined use of these new techniques leads to efficient and convergent geometry optimizations with the PCMs.

Journal Article↗

PubMed2003

A planar tetracoordinate carbon and unusual bonding in an organodimetallic propynylidene complex arising from double C-H activation of an allene ligand.

Reduction of the organoditantalum allene complex (eta-C5Me4R)2Ta2(mu-X)X3(mu-eta1,eta3-C3H4) (R = Me (Cp*), Et; X = Cl, Br) with sodium amalgam leads to the propynylidene complex (eta-C5Me4R)2Ta2(mu-H)2X2(mu-HCCCH) by a formal double 1,3-C-H activation of the allene ligand. The solid-state molecular structure contains a planar HCCCH ligand bridging, in parallel coordination mode, the two tantalum atoms, with the HCCCH and Ta atoms coplanar. Key structural features are a Ta-Ta distance of 2.8817(7) A, propynylidene C-C-C angle of 153.7(13) degrees , C-C distance of 1.370(8) A, Ta-C(central) distance of 2.194(9) A, and Ta-C(terminal) distance of 1.970(9) A. Molecular orbital calculations on the complex at the RHF/SBK(d) and B3LYP/LanL2dz levels of theory demonstrate that the propynylidene ligand is best viewed formally as an allenediylidene(4-) ligand bonded to two d0 tantalum atoms via two Ta=C(terminal) double bonds and an unusual three-center, two-electron bridge bond involving both tantalum atoms and a lone pair on the planar, tetracoordinate central carbon. There is no net Ta-Ta bonding based on the orbital analysis.

Journal Article↗

PubMed2003

Intraprotein electrostatics derived from first principles: divide-and-conquer approaches for QM/MM calculations.

Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of N(zeta) of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pK(a) of Lys55, and compared to other computational methodologies in which force field charges are employed.

Journal Article↗

PubMed2002

Boundary conditions for the Swain-Schaad relationship as a criterion for hydrogen tunneling.

Hydrogen quantum mechanical tunneling has been suggested to play a role in a wide variety of hydrogen-transfer reactions in chemistry and enzymology. An important experimental criterion for tunneling is based on the breakdown of the semiclassical prediction for the relationship among the rates of the three isotopes of hydrogen (hydrogen -H, deuterium -D, and tritium -T). This is denoted the Swain-Schaad relationship. This study examines the breakdown of the Swain-Schaad relationship as criterion for tunneling. The semiclassical (no tunneling) limit used hereto (e.g., 3.34, for H/T to D/T kinetic isotope effects), was based on simple theoretical considerations of a diatomic cleavage of a stable covalent bond, for example, a C-H bond. Yet, most experimental evidence for a tunneling contribution has come from breakdown of those relationship for a secondary hydrogen, that is, not the hydrogen whose bond is being cleaved but its geminal neighbor. Furthermore, many of the reported experiments have been mixed-labeling experiments, in which a secondary H/T kinetic isotope effect was measured for C-H cleavage, while the D/T secondary effect accompanied C-D cleavage. In experiments of this type, the breakdown of the Swain-Schaad relationship indicates both tunneling and the degree of coupled motion between the primary and secondary hydrogens. We found a new semiclassical limit (e.g., 4.8 for H/T to D/T kinetic isotope effects), whose breakdown can serve as a more reliable experimental evidence for tunneling in this common mixed-labeling experiment. We study the tunneling contribution to C-H bond activation, for which many relevant experimental and theoretical data are available. However, these studies can be applied to any hydrogen-transfer reaction. First, an extension of the original approach was applied, and then vibrational analysis studies were carried out for a model system (the enzyme alcohol dehydrogenase). Finally, the effect of complex kinetics on the observed Swain-Schaad relationship was examined. All three methods yield a new semiclassical limit (4.8), above which tunneling must be considered. Yet, it was found that for many cases the original, localized limit (3.34), holds fairly well. For experimental results that are between the original and new limits (within statistical errors), several methods are suggested that can support or exclude tunneling. These new and clearer criteria provide a basis for future applications of the Swain-Schaad relationship to demonstrate tunneling in complex systems.

Alcohol Dehydrogenase↗