Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.

Spectroscopic properties from hybrid QM, MM molecular dynamics simulation

The central aim of this thesis work is the application and further development of a hybrid quantum mechanical/molecular mechanics (QM/MM) based approach to compute spectroscopic properties of molecules in complex chemical environments from electronic structure theory. In the framework of this thesis, an existing density functional theory implementation of the QM/MM approach is first used to calculate the nuclear magnetic resonance (NMR) solvent shifts of an adenine molecule in aqueous solution. The findings show that the aqueous solvation with its strongly fluctuating hydrogen bond network leads to specific changes in the NMR resonance lines. Besides the absolute values, also the ordering of the NMR lines changes under the influence of the solvating water molecules. Without the QM/MM scheme, a quantum chemical calculation could have led to an incorrect assignment of these lines. The second part of this thesis describes a methodological improvement of the QM/MM method that is designed for cases in which a covalent chemical bond crosses the QM/MM boundary. The development consists in an automatized protocol to optimize a so-called capping potential that saturates the electronic subsystem in the QM region. The optimization scheme is capable of tuning the parameters in such a way that the deviations of the electronic orbitals between the regular and the truncated (and "capped") molecule are minimized. This in turn results in a considerable improvement of the structural and spectroscopic parameters when computed with the new optimized capping potential within the QM/MM technique. This optimization scheme is applied and benchmarked on the example of truncated carbon-carbon bonds in a set of small test molecules. It turns out that the optimized capping potentials yield an excellent agreement of NMR chemical shifts and protonation energies with respect to the corresponding full molecules. These results are very promising, so that the application to larger biological complexes will significantly improve the reliability of the prediction of the related spectroscopic properties.

Ab-initio study of hydrogen technology materials for hydrogen storage and proton conduction

This dissertation deals with two specific aspects of a potential hydrogen-based energy economy, namely the problems of energy storage and energy conversion. In order to contribute to the solution of these problems, the structural and dynamical properties of two promising materials for hydrogen storage (lithium imide/amide) and proton conduction (poly[vinyl phosphonic acid]) are modeled on an atomistic scale by means of first principles molecular dynamics simulation methods.rnrnrnIn the case of the hydrogen storage system lithium amide/imide (LiNH_2/Li_2NH), the focus was on the interplay of structural features and nuclear quantum effects. For these calculations, Path-Integral Molecular Dynamics (PIMD) simulations were used. The structures of these materials at room temperature were elucidated; in collaboration with an experimental group, a very good agreement between calculated and experimental solid-state 1H-NMR chemical shifts was observed. Specifically, the structure of Li_2NH features a disordered arrangement of the Li lattice, which was not reported in previous studies. In addition, a persistent precession of the NH bonds was observed in our simulations. We provide evidence that this precession is the consequence of a toroid-shaped effective potential, in which the protons in the material are immersed. This potential is essentially flat along the torus azimuthal angle, which might lead to important quantum delocalization effects of the protons over the torus.rnrnOn the energy conversion side, the dynamics of protons in a proton conducting polymer (poly[vinyl phosphonic acid], PVPA) was studied by means of a steered ab-initio Molecular Dynamics approach applied on a simplified polymer model. The focus was put on understanding the microscopic proton transport mechanism in polymer membranes, and on characterizing the relevance of the local environment. This covers particularly the effect of water molecules, which participate in the hydrogen bonding network in the material. The results indicate that these water molecules are essential for the effectiveness of proton conduction. A water-mediated Grotthuss mechanism is identified as the main contributor to proton conduction, which agrees with the experimentally observed decay on conductivity for the same material in the absence of water molecules.rnrnThe gain in understanding the microscopic processes and structures present in this materials can help the development of new materials with improved properties, thus contributing to the solution of problems in the implementation of fuel cells.

Metallkomplexe mit sterisch anspruchsvollen mono- und bifunktionellen Triazenidoliganden

Monoanionische, stickstoffhaltige Chelatliganden sind sehr populär und werden oft als Baueinheiten zur Synthese von Metallkomplexen eingesetzt. Im Gegensatz zu den bekannten Amidinat- und β-Diketiminatliganden wurde den Triazeniden bisher weniger Aufmerksamkeit geschenkt. Dies liegt hauptsächlich daran, dass wenige, sterisch anspruchsvolle Vertreter dieser Gruppe bekannt waren. Mit der Synthese neuer biphenyl- und terphenylsubstituierter Triazenide sind diese stärker in den Fokus der Forschung gerückt. Durch die Sperrigkeit der Liganden werden unerwünschte Ligandenumverteilungsreaktionen unterdrückt und der Weg zu neuen Metallkomplexen mit bemerkenswerten Eigenschaften geebnet. Diese Arbeit beschäftigt sich mit Triazenidokomplexen der Alkalimetalle, des Magnesiums, des Zinks und des Aluminiums. Zusätzlich zu den monofunktionellen wurden bifunktionelle Triazene, welche aus einer Biphenylgruppe aufgebaut sind und zwei terphenylsubstituierte Triazengruppen tragen, synthetisiert und deren Funktion als Komplexliganden untersucht. Durch das hohe Volumen der Liganden konnten Metallkomplexe hergestellt werden, die bemerkenswerte strukturelle Eigenschaften aufweisen. So bilden beispielsweise alle dargestellten Verbindungen Kristalle, die aus monomeren Baueinheiten aufgebaut sind. Eine homoleptische Magnesiumverbindung mit zwei Triazenidliganden weist eine ungewöhnliche, planare Koordinationssphäre auf, eine Triazenidodiphenylaluminiumverbindung zeigt eine interessante Ausrichtung der Phenylgruppen, die auf π-Stapeleffekte zurückgeführt wird und die Alkalimetalltriazenide weisen sekundäre Metall-π-Aren-Wechselwirkungen auf, womit eine bessere koordinative Absättigung des Koordinationszentrums erreicht wird. Zusätzlich zur strukturellen Charakterisierung mittels Röntgendiffraktometrie wurden, wo es möglich war 1H und 13C NMR-Messungen durchgeführt, Schmelzpunkte bestimmt sowie IR- und UV-VIS-Spektren aufgenommen. Durch die hohe Dichte an aromatischen Systemen innerhalb der Verbindungen, ist es bei einer dreikernigen Zinktriazenidverbindung gelungen mittels 1H NMR-Messungen die starken Einflüsse des Ringstroms nachzuweisen. Außerdem wurden DFT-Berechnungen durchgeführt, um herauszufinden, wie sehr die räumliche Ausrichtung der Bi- bzw. Terphenylsubstituenten gegenüber der N3-Baugruppe der Triazenide die Verteilung der negativen Ladung innerhalb eines monoanionischen Triazenids und somit dessen Koordinationseigenschaften beeinflusst. Diese Effekte konnten durch experimentelle Ergebnisse unterstützt werden.

Studying Electrophoretic Separations Using Cholate Micelles: Effects Of Ph, Temperature, And Composition

Micelle-forming bile salts have previously been shown to be effective pseudo-stationary phases for separating the chiral isomers of binaphthyl compounds with micellar electrokinetic capillary chromatography (MEKC). Here, cholate micelles are systematically investigated via electrophoretic separations and NMR using R, S-1, 1¿- binaphthyl- 2, 2¿-diylhydrogenphosphate (BNDHP) as a model chiral analyte. The pH, temperature, and concentration of BNDHP were systematically varied while monitoring the chiral resolution obtained with MEKC and the chemical shift of various protons in NMR. NMR data for each proton on BNDHP is monitored as a function of cholate concentration: as cholate monomers begin to aggregate and the analyte molecules begin to sample the micelle aggregate we observe changes in the cholate methyl and S-BNDHP proton chemical shifts. From such NMR data, the apparent CMC of cholate at pH 12 is found to be about 13-14 mM, but this value decreases at higher pH, suggesting that more extreme pHs may give rise to more effective separations. In general, CMCs increase with temperature indicating that one may be able to obtain better separations at lower temperatures. S-BNDHP concentrations ranging from 50 ¿M to 400 ¿M (pH 12.8) gave rise to apparent cholate CMC values from 10 mM to 8 mM, respectively, indicating that S-BNDHP, the chiral analyte molecule, may play an active role in stabilizing cholate aggregates. In all, these data show that NMR can be used to systematically investigate a complex multi-variable landscape of potential optimizations of chiral separations.

Understanding Bile Micelle Chiral Interactions

Bile salts are known to aggregate into micelles in biological systems; however, the fundamental structure and dynamics of bile molecule micelle formation are poorly understood. Previous studies have established that the bile salt cholate is capable of performing chirally selective micellar electrokinetic capillary chromatography (MEKC) separations of model racemic binaphthyl compounds 1,1¿-binaphthyl-2,2¿-diyl hydrogen phosphate (R,S-BNDHP) and 1,1¿-bi-2-naphthol (R,S-BN). Nuclear magnetic resonance (NMR) has been established as a complementary technique for understanding chiral selectivity and micelle formation events based on changes in proton chemical shifts of the probe molecules BNDHP and BN as well as of cholate. This work investigated the effects of the probe molecule, the alkali cation identity and temperature on cholate micelle aggregation and MEKC separations of R,S-BN and R,S-BNDHP. The probe molecule was found to mediate micelle formation by MEKC and proton NMR. A low (0.1 mM) concentration of probe was found to have minimal effects on micellization events detected by proton NMR while higher probe concentration (2.5 mM) was found to mediate micellization causing micellization events to occur at lower cholate concentrations. This work also investigated the effects of alkali counterion on chiral separation. Generally, counterions with larger crystal cationic radius were found to cause greater chiral separation power. NMR data suggest that protons near the surface of the cholate micelle are most sensitive to the cation identity, suggesting a model of improved separation based on the cation sterically inhibiting binding of one isomer. Finally, the effect of temperature on MEKC separation was investigated. Separation power of R,S-BN and R,S-BNDHP appeared to increase linearly with temperature for 22.0 mM to 50.0 mM pH 12.0 cholate. In total, these results indicate that cholate aggregation is dependent on multiple conditions. Understanding the roles that these factors play in influencing cholate micellization can inform better separation in MEKC.

Time-dependent interactions of the two porphyrinic compounds chlorin e6 and mono-L-aspartyl-chlorin e6 with phospholipid vesicles probed by NMR spectroscopy

The distribution processes of chlorin e6 (CE) and monoaspartyl-chlorin e6 (MACE) between the outer and inner phospholipid monolayers of 1,2-dioleoyl-phosphatidylcholine (DOPC) vesicles were monitored by 1H NMR spectroscopy through analysis of chemical shifts and line widths of the DOPC vesicle resonances. Chlorin adsorption to the outer vesicle monolayer induced changes in the DOPC 1H NMR spectrum. Most pronounced was a split of the N-methyl choline resonance, allowing for separate analysis of inner and outer vesicle layers. Transbilayer distribution of the chlorin compounds was indicated by time-dependent characteristic spectral changes of the DOPC resonances. Kinetic parameters for the flip-flop processes, that is, half-lives and rate constants, were obtained from the experimental data points. In comparison to CE, MACE transbilayer movement was significantly reduced, with MACE remaining more or less attached to the outer membrane layer. The distribution coefficients for CE and MACE between the vesicular and aqueous phase were determined. Both CE and MACE exhibited a high affinity for the vesicular phase. For CE, a positive correlation was found between transfer rate and increasing molar ratio CE/DOPC. Enhanced membrane rigidity induced by increasing amounts of cholesterol into the model membrane was accompanied by a decrease of CE flip-flop rates across the membrane. The present study shows that the movement of porphyrins across membranes can efficiently be investigated by 1H NMR spectroscopy and that small changes in porphyrin structure can have large effects on membrane kinetics.

Molecular mechanisms of Ca$\sp{2+}$ signal transduction from cardiac troponin C to troponin I

$\rm Ca\sp{2+}$-dependent exposure of an N-terminal hydrophobic region in troponin C (TnC) is thought to be important for the regulation of contraction in striated muscle. To study these conformational changes in cardiac troponin (cTnC), the $\varepsilon$C and $\varepsilon$H chemical shifts for all 10 Met residues in cTnC were sequence-specific assigned on NMR spectra using a combination of two dimensional NMR techniques and site-directed mutagenesis. The assigned methyl-Met chemical shifts were used as structural markers to monitor conformational changes induced by $\rm Ca\sp{2+}.$ The results showed that binding of $\rm Ca\sp{2+}$ to the regulatory site in the N-domain induced large changes in the $\varepsilon$H and $\varepsilon$C chemical shifts of Met 45, Met 80, Met 81 in the predicted N-terminal hydrophobic region, but had no effect on the chemical shifts of Met residues located in the C-domain. These results suggest that the $\rm Ca\sp{2+}$-dependent functions of cTnC are mainly through N-terminal domain of cTnC.^ To further define the molecular mechanism by which TnC regulates muscle contraction, single Cys residues were engineered at positions 45, 81, 84 or 85 in the N-terminal hydrophobic region of cTnC to provide sites for attachment of specific blocking groups. Blocking groups were coupled to these Cys residues in cTnC mutants and the covalent adducts were tested for activity in TnC-extracted myofibrils. Covalent modification of cTnC(C45) had no effect on maximal myofibril ATPase activity. Greatly decreased myofibril ATPase activity resulted when the peptide or biotin was conjugated to residue 81 in cTnC(C81), while less inhibition resulted from covalent modification of cTnC(C84) or cTnC(C85). The results suggest that limited sites of the N-terminal hydrophobic region in cTnC are important for transducing the $\rm Ca\sp{2+}$ signal to troponin I (TnI) and are sensitive to modification, while other regions are less important or can adapt to steric hindrances introduced by bulky blocking groups.^ Although the exposed TnI interaction site in the N-terminal hydrophobic region of TnC is crucial for function of TnC, other regions in the N-domain of TnC may also participate in transducing the $\rm Ca\sp{2+}$ signal and conferring the maximal activation of actomyosin ATPase. The interactions between the B-/C-helices of cTnC and cTnI were characterized using a combination of site-directed mutagenesis, fluorescence and covalent modification. The results suggest that the $\rm Ca\sp{2+}$-dependent interactions of the B-/C-helices of cTnC with TnI may be required for the maximal activation of muscle contraction. ^

NMR structural analysis of cardiac troponin C: Monitoring conformational changes induced by binding calcium, troponin I, and calcium -sensitizing drugs

Contraction of cardiac muscle is regulated through the Ca2+ dependent protein-protein interactions of the troponin complex (Tn). The critical role cardiac troponin C (cTnC) plays as the Ca2+ receptor in this complex makes it an attractive target for positive inotropic compounds. In this study, the ten Met methyl groups in cTnC, [98% 13C ϵ]-Met cTnC, are used as structural markers to monitor conformational changes in cTnC and identify sites of interaction between cTnC and cardiac troponin I (cTnI) responsible for the Ca2+ dependent interactions. In addition the structural consequences that a number of Ca2+-sensitizing compounds have on free cTnC and the cTnC·cTnI complex were characterized. Using heteronuclear NMR experiments and monitoring chemical shift changes in the ten Met methyl 1H-13C correlations in 3Ca2+ cTnC when bound to cTnI revealed an anti-parallel arrangement for the two proteins such that the N-domain of cTnI interacts with the C-domain of cTnC. The large chemical shifts in Mets-81, -120, and -157 identified points of contact between the proteins that include the C-domain hydrophobic surface in cTnC and the A, B, and D helical interface located in the regulatory N-domain of cTnC. TnI association [cTnI(33–80), cTnI(86–211), or cTnI(33–211)] was found also to dramatically reduce flexibility in the D/E central linker of cTnC as monitored by line broadening in the Met 1H- 13C correlations of cTnC induced by a nitroxide spin label, MTSSL, covalently attached to cTnC at Cys 84. TnI association resulted in an extended cTnC that is unlike the compact structure observed for free cTnC. The Met 1H-13C correlations also allowed the binding characteristics of bepridil, TFP, levosimendan, and EMD 57033 to the apo, 2Ca2+, and Ca2+ saturated forms of cTnC to be determined. In addition, the location of drug binding on the 3Ca2+cTnC·cTnI complex was identified for bepridil and TFP. Use of a novel spin-labeled phenothiazine, and detection of isotope filtered NOEs, allowed identification of drug binding sites in the shallow hydrophobic cup in the C-terminal domain, and on two hydrophobic surfaces on N-regulatory domain in free 3Ca2+ cTnC. In contrast, only one N-domain drug binding site exists in 3Ca2+ cTnC·cTnI complex. The methyl groups of Met 45, 60 and 80, which are grouped in a hydrophobic patch near site II in cTnC, showed the greatest change upon titration with bepridil or TFP, suggesting that this is a critical site of drug binding in both free cTnC and when associated with cTnI. The strongest NOEs were seen for Met-60 and -80, which are located on helices C and D, respectively, of Ca2+ binding site II. These results support the conclusion that the small hydrophobic patch which includes Met-45, -60, and -80 constitutes a drug binding site, and that binding drugs to this site will lead to an increase in Ca2+ binding affinity of site II while preserving maximal cTnC activity. Thus, the subregion in cTnC makes a likely target against which to design new and selective Ca2+-sensitizing compounds. ^

1H HR MAS NMR metabolomic and non-destructive 2D NMR relaxometry to assess internal quality in apples.

NMR can be considered a multi-scale multidimensional technology in the sense that it provides both spatial insight at macroscopic (MRI) or microscopic level (relaxometry), together with chemical characterization (HR-MAS). In this study 296 apples (from 4 cultivars) were MRI screened (20 slices per fruit) among which 7 fruits were used for metabolomic study by 1H HR MAS in order to assess various chemical shifts: malic acid, sucrose, glucose, fructose and ethanol. On the first season, tissue samples were taken from the sound and affected apples (near the core, centre and outer part of the mesocarp) belonging to sound and affected locations, while on the second season, tissue samples were focused on the comparison between sound and affected tissue. Beside, MRI and 2D non-destructive relaxometry (on whole fruits, and localized tissue) where performed on 72 and 12 apples respectively in order to compare features at macroscopic (tissue) and microscopic (subcellular) level. HR MAS shows higher content of ?-glucose, ?-glucose, malic acid and aromatic compounds in watercore affected tissues from both seasons, while sound tissue reflects higher sucrose. Microscopic (subcellular) degradation of tissue varies according to disorder development and is in good accordance with macroscopic characterization with MRI.

NMR spectroscopy in studies of light-induced structural changes in mammalian rhodopsin: Applicability of solution 19F NMR

We report high resolution solution 19F NMR spectra of fluorine-labeled rhodopsin mutants in detergent micelles. Single cysteine substitution mutants in the cytoplasmic face of rhodopsin were labeled by attachment of the trifluoroethylthio (TET), CF3-CH2-S, group through a disulfide linkage. TET-labeled cysteine mutants at amino acid positions 67, 140, 245, 248, 311, and 316 in rhodopsin were thus prepared. Purified mutant rhodopsins (6–10 mg), in dodecylmaltoside, were analyzed at 20°C by solution 19F NMR spectroscopy. The spectra recorded in the dark showed the following chemical shifts relative to trifluoroacetate: Cys-67, 9.8 ppm; Cys-140, 10.6 ppm; Cys-245, 9.9 ppm; Cys-248, 9.5 ppm; Cys-311, 9.9 ppm; and Cys-316, 10.0 ppm. Thus, all mutants showed chemical shifts downfield that of free TET (6.5 ppm). On illumination to form metarhodopsin II, upfield changes in chemical shift were observed for 19F labels at positions 67 (−0.2 ppm) and 140 (−0.4 ppm) and downfield changes for positions 248 (+0.1 ppm) and 316 (+0.1 ppm) whereas little or no change was observed at positions 311 and 245. On decay of metarhodopsin II, the chemical shifts reverted largely to those originally observed in the dark. The results demonstrate the applicability of solution 19F NMR spectroscopy to studies of the tertiary structures in the cytoplasmic face of intact rhodopsin in the dark and on light activation.

NMR scalar couplings across Watson–Crick base pair hydrogen bonds in DNA observed by transverse relaxation-optimized spectroscopy

This paper describes the NMR observation of 15N—15N and 1H—15N scalar couplings across the hydrogen bonds in Watson–Crick base pairs in a DNA duplex, hJNN and hJHN. These couplings represent new parameters of interest for both structural studies of DNA and theoretical investigations into the nature of the hydrogen bonds. Two dimensional [15N,1H]-transverse relaxation-optimized spectroscopy (TROSY) with a 15N-labeled 14-mer DNA duplex was used to measure hJNN, which is in the range 6–7 Hz, and the two-dimensional hJNN-correlation-[15N,1H]-TROSY experiment was used to correlate the chemical shifts of pairs of hydrogen bond-related 15N spins and to observe, for the first time, hJHN scalar couplings, with values in the range 2–3.6 Hz. TROSY-based studies of scalar couplings across hydrogen bonds should be applicable for large molecular sizes, including protein-bound nucleic acids.

Fractionation factors and activation energies for exchange of the low barrier hydrogen bonding proton in peptidyl trifluoromethyl ketone complexes of chymotrypsin

NMR investigations have been carried out of complexes between bovine chymotrypsin Aα and a series of four peptidyl trifluoromethyl ketones, listed here in order of increasing affinity for chymotrypsin: N-Acetyl-l-Phe-CF3, N-Acetyl-Gly-l-Phe-CF3, N-Acetyl-l-Val-l-Phe-CF3, and N-Acetyl-l-Leu-l-Phe-CF3. The D/H fractionation factors (φ) for the hydrogen in the H-bond between His 57 and Asp 102 (His 57-Hδ1) in these four complexes at 5°C were in the range φ = 0.32–0.43, expected for a low-barrier hydrogen bond. For this series of complexes, measurements also were made of the chemical shifts of His 57-Hɛ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 8.97–9.18), the exchange rate of the His 57-Hδ1 proton with bulk water protons (284–12.4 s−1), and the activation enthalpies for this hydrogen exchange (14.7–19.4 kcal⋅mol−1). It was found that the previously noted correlations between the inhibition constants (Ki 170–1.2 μM) and the chemical shifts of His 57-Hδ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 18.61–18.95) for this series of peptidyl trifluoromethyl ketones with chymotrypsin [Lin, J., Cassidy, C. S. & Frey, P. A. (1998) Biochemistry 37, 11940–11948] could be extended to include the fractionation factors, hydrogen exchange rates, and hydrogen exchange activation enthalpies. The results support the proposal of low barrier hydrogen bond-facilitated general base catalysis in the addition of Ser 195 to the peptidyl carbonyl group of substrates in the mechanism of chymotrypsin-catalyzed peptide hydrolysis. Trends in the enthalpies for hydrogen exchange and the fractionation factors are consistent with a strong, double-minimum or single-well potential hydrogen bond in the strongest complexes. The lifetimes of His 57-Hδ1, which is solvent shielded in these complexes, track the strength of the hydrogen bond. Because these lifetimes are orders of magnitude shorter than those of the complexes themselves, the enzyme must have a pathway for hydrogen exchange at this site that is independent of dissociation of the complexes.

NMR reveals hydrogen bonds between oxygen and distal histidines in oxyhemoglobin

Compared with free heme, the proteins hemoglobin (Hb) and myoglobin (Mb) exhibit greatly enhanced affinity for oxygen relative to carbon monoxide. This physiologically vital property has been attributed to either steric hindrance of CO or stabilization of O2 binding by a hydrogen bond with the distal histidine. We report here the first direct evidence of such a hydrogen bond in both α- and β-chains of oxyhemoglobin, as revealed by heteronuclear NMR spectra of chain-selectively labeled samples. Using these spectra, we have assigned the imidazole ring 1H and 15N chemical shifts of the proximal and distal histidines in both carbonmonoxy- and oxy-Hb. Because of their proximity to the heme, these chemical shifts are extremely sensitive to the heme pocket conformation. Comparison of the measured chemical shifts with values predicted from x-ray structures suggests differences between the solution and crystal structures of oxy-Hb. The chemical shift discrepancies could be accounted for by very small displacements of the proximal and distal histidines. This suggests that NMR could be used to obtain very high-resolution heme pocket structures of Hb, Mb, and other heme proteins.

Role of the J-domain in the cooperation of Hsp40 with Hsp70

The Escherichia coli Hsp40 DnaJ and Hsp70 DnaK cooperate in the binding of proteins at intermediate stages of folding, assembly, and translocation across membranes. Binding of protein substrates to the DnaK C-terminal domain is controlled by ATP binding and hydrolysis in the N-terminal ATPase domain. The interaction of DnaJ with DnaK is mediated at least in part by the highly conserved N-terminal J-domain of DnaJ that includes residues 2–75. Heteronuclear NMR experiments with uniformly 15N-enriched DnaJ2–75 indicate that the chemical environment of residues located in helix II and the flanking loops is perturbed on interaction with DnaK or a truncated DnaK molecule, DnaK2–388. NMR signals corresponding to these residues broaden and exhibit changes in chemical shifts in the presence of DnaK(MgADP). Addition of MgATP largely reversed the broadening, indicating that NMR signals of DnaJ2–75 respond to ATP-dependent changes in DnaK. The J-domain interaction is localized to the ATPase domain of DnaK and is likely to be dominated by electrostatic interactions. The results suggest that the J-domain tethers DnaK to DnaJ-bound substrates, which DnaK then binds with its C-terminal peptide-binding domain.