Excitonic couplings and electronic coherence in bridged naphthalene dimers

The electronic excitations of naphthalene and a family of bridged naphthalene dimers are calculated and analyzed by using the Collective Electronic Oscillator method combined with the oblique Lanczos algorithm. All experimentally observed trends in absorption profiles and radiative lifetimes are reproduced. Each electronic excitation is linked to the corresponding real-space transition density matrix, which represents the motions of electrons and holes created in the molecule by photon absorption. Two-dimensional plots of these matrices help visualize the degree of exciton localization and explain the dependence of the electronic interaction between chromophores on their separation.

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.

Structural factors controlling ligand binding to myoglobin: A kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as “kinetic hole burning.” The study of these correlations gives direct insights into the structure–function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Ab initio calculation of electronic coupling in the photosynthetic reaction center

We have carried out an ab initio electronic structure calculations of electron transfer couplings between chromophores in the bacterial photosynthetic reaction center. The couplings agree remarkably well with parameters obtained from recent quantum dynamical modeling of experimental data assuming an explicit intermediate mechanism. We also have computed couplings on the M-side of the reaction center and have found that the interaction of the primary donor to the M-side intermediate bacteriochlorophyll is quite small because of destructive interference of the two localized coupling matrix elements. This may explain the slow rate of electron transfer down the M-side of the reaction center.

Solid and liquid phase 59Co NMR studies of cobalamins and their derivatives

We describe the application of 59Co NMR to the study of naturally occurring cobalamins. Targets of these investigations included vitamin B12, the B12 coenzyme, methylcobalamin, and dicyanocobyrinic acid heptamethylester. These measurements were carried out on solutions and powders of different origins, and repeated at a variety of magnetic field strengths. Particularly informative were the solid-state central transition NMR spectra, which when combined with numerical line shape analyses provided a clear description of the cobalt coupling parameters. These parameters showed a high sensitivity to the type of ligands attached to the metal and to the crystallization history of the sample. 59Co NMR determinations also were carried out on synthetic cobaloximes possessing alkyl, cyanide, aquo, and nitrogenated axial groups, substituents that paralleled the coordination of the natural compounds. These analogs displayed coupling anisotropies comparable to those of the cobalamins, as well as systematic up-field shifts that can be rationalized in terms of their stronger binding affinity to the cobalt atom. Cobaloximes also displayed a higher regularity in the relative orientations of their quadrupole and shielding coupling tensors, reflecting a higher symmetry in their in-plane coordination. For the cobalamines, poor correlations were observed between the values measured for the quadrupole couplings in the solid and the line widths observed in the corresponding solution 59Co NMR resonances.

Femtosecond dynamics of the forbidden carotenoid S1 state in light-harvesting complexes of purple bacteria observed after two-photon excitation

Time-resolved excited-state absorption intensities after direct two-photon excitation of the carotenoid S1 state are reported for light-harvesting complexes of purple bacteria. Direct excitation of the carotenoid S1 state enables the measurement of subsequent dynamics on a fs time scale without interference from higher excited states, such as the optically allowed S2 state or the recently discovered dark state situated between S1 and S2. The lifetimes of the carotenoid S1 states in the B800-B850 complex and B800-B820 complex of Rhodopseudomonas acidophila are 7 ± 0.5 ps and 6 ± 0.5 ps, respectively, and in the light-harvesting complex 2 of Rhodobacter sphaeroides ≈1.9 ± 0.5 ps. These results explain the differences in the carotenoid to bacteriochlorophyll energy transfer efficiency after S2 excitation. In Rps. acidophila the carotenoid S1 to bacteriochlorophyll energy transfer is found to be quite inefficient (φET1 <28%) whereas in Rb. sphaeroides this energy transfer is very efficient (φET1 ≈80%). The results are rationalized by calculations of the ensemble averaged time constants. We find that the Car S1 → B800 electronic energy transfer (EET) pathway (≈85%) dominates over Car S1 → B850 EET (≈15%) in Rb. sphaeroides, whereas in Rps. acidophila the Car S1 → B850 EET (≈60%) is more efficient than the Car S1 → B800 EET (≈40%). The individual electronic couplings for the Car S1 → BChl energy transfer are estimated to be approximately 5–26 cm−1. A major contribution to the difference between the energy transfer efficiencies can be explained by different Car S1 energy gaps in the two species.

A nicked duplex decamer DNA with a PEG6 tether

A dumbbell double-stranded DNA decamer tethered with a hexaethylene glycol linker moiety (DDSDPEG), with a nick in the centre of one strand, has been synthesised. The standard NMR methods, E.COSY, TOCSY, NOESY and HMQC, were used to measure 1H, 31P and T1 spectral parameters. Molecular modelling using rMD-simulated annealing was used to compute the structure. Scalar couplings and dipolar contacts show that the molecule adopts a right-handed B-DNA helix in 38 mM phosphate buffer at pH 7. Its high melting temperature confirms the good base stacking and stability of the duplex. This is partly attributed to the presence of the PEG6 linker at both ends of the duplex that restricts the dynamics of the stem pentamers and thus stabilises the oligonucleotide. The inspection of the global parameters shows that the linker does not distort the B-DNA geometry. The computed structure suggests that the presence of the nick is not disturbing the overall tertiary structure, base pair geometry or duplex base pairing to a substantial extent. The nick has, however, a noticeable impact on the local geometry at the nick site, indicated clearly by NMR analysis and reflected in the conformational parameters of the computed structure. The 1H spectra also show much sharper resonances in the presence of K+ indicating that conformational heterogeneity of DDSDPEG is reduced in the presence of potassium as compared to sodium or caesium ions. At the same time the 1H resonances have longer T1 times. This parameter is suggested as a sensitive gauge of stabilisation.

Symmetry in chemistry from the hydrogen atom to proteins

The last 2 decades have seen discoveries in highly excited states of atoms and molecules of phenomena that are qualitatively different from the “planetary” model of the atom, and the near-rigid model of molecules, characteristic of these systems in their low-energy states. A unified view is emerging in terms of approximate dynamical symmetry principles. Highly excited states of two-electron atoms display “molecular” behavior of a nonrigid linear structure undergoing collective rotation and vibration. Highly excited states of molecules described in the “standard molecular model” display normal mode couplings, which induce bifurcations on the route to molecular chaos. New approaches such as rigid–nonrigid correlation, vibrons, and quantum groups suggest a unified view of collective electronic motion in atoms and nuclear motion in molecules.

New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary 13C shifts, 1H shifts, T1/T2 ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.

Functional dynamics in the active site of the ribonuclease binase

Binase, a member of a family of microbial guanyl-specific ribonucleases, catalyzes the endonucleotic cleavage of single-stranded RNA. It shares 82% amino acid identity with the well-studied protein barnase. We used NMR spectroscopy to study the millisecond dynamics of this small enzyme, using several methods including the measurement of residual dipolar couplings in solution. Our data show that the active site of binase is flanked by loops that are flexible at the 300-μs time scale. One of the catalytic residues, His-101, is located on such a flexible loop. In contrast, the other catalytic residue, Glu-72, is located on a β-sheet, and is static. The residues Phe-55, part of the guanine base recognition site, and Tyr-102, stabilizing the base, are the most dynamic. Our findings suggest that binase possesses an active site that has a well-defined bottom, but which has sides that are flexible to facilitate substrate access/egress, and to deliver one of the catalytic residues. The motion in these loops does not change on complexation with the inhibitor d(CGAG) and compares well with the maximum kcat (1,500 s−1) of these ribonucleases. This observation indicates that the NMR-measured loop motions reflect the opening necessary for product release, which is apparently rate limiting for the overall turnover.

Molecular cloning and sequence analysis of the complestatin biosynthetic gene cluster

Streptomyces lavendulae produces complestatin, a cyclic peptide natural product that antagonizes pharmacologically relevant protein–protein interactions including formation of the C4b,2b complex in the complement cascade and gp120-CD4 binding in the HIV life cycle. Complestatin, a member of the vancomycin group of natural products, consists of an α-ketoacyl hexapeptide backbone modified by oxidative phenolic couplings and halogenations. The entire complestatin biosynthetic and regulatory gene cluster spanning ca. 50 kb was cloned and sequenced. It consisted of 16 ORFs, encoding proteins homologous to nonribosomal peptide synthetases, cytochrome P450-related oxidases, ferredoxins, nonheme halogenases, four enzymes involved in 4-hydroxyphenylglycine (Hpg) biosynthesis, transcriptional regulators, and ABC transporters. The nonribosomal peptide synthetase consisted of a priming module, six extending modules, and a terminal thioesterase; their arrangement and domain content was entirely consistent with functions required for the biosynthesis of a heptapeptide or α-ketoacyl hexapeptide backbone. Two oxidase genes were proposed to be responsible for the construction of the unique aryl-ether-aryl-aryl linkage on the linear heptapeptide intermediate. Hpg, 3,5-dichloro-Hpg, and 3,5-dichloro-hydroxybenzoylformate are unusual building blocks that repesent five of the seven requisite monomers in the complestatin peptide. Heterologous expression and biochemical analysis of 4-hydroxyphenylglycine transaminon confirmed its role as an aminotransferase responsible for formation of all three precursors. The close similarity but functional divergence between complestatin and chloroeremomycin biosynthetic genes also presents a unique opportunity for the construction of hybrid vancomycin-type antibiotics.

Studies on the biosynthesis of taxol: the taxane carbon skeleton is not of mevalonoid origin.

A cell culture of Taxus chinensis was established to produce the diterpene 2alpha,5alpha,10beta,14beta-tetra-acetoxy4 ++ +(20),11-taxadiene (taxuyunnanine C) in 2.6% (dry weight) yield. The incorporation of [U-13C6]glucose, [1-13C]glucose, and [1,2-13C2]acetate into this diterpene was analyzed by NMR spectroscopy. Label from [1,2-13C2]acetate was diverted to the four acetyl groups of taxuyunnanine C, but not to the taxane ring system. Label from [1-13C]glucose and [U-13C6]glucose was efficiently incorporated into both the taxane ring system and the acetyl groups. The four isoprenoid moieties of the diterpene showed identical labeling patterns. The analysis of long-range 13C13C couplings in taxuyunnanine C obtained from an experiment with [U-13C6]glucose documents the involvement of an intramolecular rearrangement in the biosynthesis of the isoprenoid precursor. The labeling patterns are inconsistent with the mevalonate pathway. The taxoid data share important features with the alternative pathway of isoprenoid biosynthesis operating in certain eubacteria Rohmer, M., Knani, M., Simonin, P., Sutter, B. & Sahm, H. (1993) Biochem. J. 295, 517-524].

Sample size determination in combinatorial chemistry.

Combinatorial chemistry is gaining wide appeal as a technique for generating molecular diversity. Among the many combinatorial protocols, the split/recombine method is quite popular and particularly efficient at generating large libraries of compounds. In this process, polymer beads are equally divided into a series of pools and each pool is treated with a unique fragment; then the beads are recombined, mixed to uniformity, and redivided equally into a new series of pools for the subsequent couplings. The deviation from the ideal equimolar distribution of the final products is assessed by a special overall relative error, which is shown to be related to the Pearson statistic. Although the split/recombine sampling scheme is quite different from those used in analysis of categorical data, the Pearson statistic is shown to still follow a chi2 distribution. This result allows us to derive the required number of beads such that, with 99% confidence, the overall relative error is controlled to be less than a pregiven tolerable limit L1. In this paper, we also discuss another criterion, which determines the required number of beads so that, with 99% confidence, all individual relative errors are controlled to be less than a pregiven tolerable limit L2 (0 < L2 < 1).

Lattice model for rapidly folding protein-like heteropolymers.

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state. Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude. The evolutionary implications of our findings are discussed.

Nuclear magnetic dipole interactions in field-oriented proteins: information for structure determination in solution.

The measurement of dipolar contributions to the splitting of 15N resonances of 1H-15N amide pairs in multidimensional high-field NMR spectra of field-oriented cyanometmyoglobin is reported. The splittings appear as small field-dependent perturbations of normal scalar couplings. Assignment of more than 90 resonances to specific sequential sites in the protein allows correlation of the dipolar contributions with predictions based on the known susceptibility and known structure of the protein. Implications as an additional source of information for protein structure determination in solution are discussed.