Computational and experimental studies of the electronic excitation spectra of EDTA and DTPA substituted tetraphenylporphyrins and their Lu complexes

Ethylendiaminetetraacetic acid (EDTA) substituted and diethylenetriaminopentaacetic acid (DTPA) substituted aminated free-base tetraphenylporphyrins (H₂ATPP) and the corresponding lutetium(III) complexes have been studied computationally at the density functional theory (DFT) and second-order algebraic diagrammatic construction (ADC(2)) levels using triple-ξ basis sets augmented with polarization functions. The molecular structures were optimized using Becke's three-parameter hybrid functional (B3LYP). The electronic excitation spectra in the range of 400–700 nm were calculated using the ADC(2) and the linear-response time-dependent DFT methods. The calculated spectra are compared to those measured in ethanol solution. The calculated excitation energies agree well with those deduced from the experimental spectra. The excitation energies for the Q_x band calculated at the B3LYP and ADC(2) level are 0.20-0.25 eV larger than the experimental values. The excitation energies for the Q_y band calculated at the B3LYP level are 0.10-0.20 eV smaller than the ADC(2) ones and are thus in good agreement with experiment. The calculated excitation energies corresponding to the B_x and B_y bands are 0.10-0.30 eV larger than the experimental values. The excitation energies of the B_x and B_y bands calculated at the B3LYP level are in somewhat better agreement with experiment than the ADC(2) ones. The calculated and measured band strengths largely agree.

The chiroptical properties of proteins. II. Near-ultraviolet circular dichroism of lysozyme

A study of the near-uv CD spectrum of lysozyme was carried out in the presence and absence of the inhibitor tri-N-acetylglucosamine, and theoretical chiroptical calculations based on the tetragonal crystal structure of the enzyme and the enzyme-inhibitor complex were performed. The results of these calculations indicate that the near-uv CD spectrum of lysozyme can be adequately explained in terms of negative rotatory strengths arising from the tryptophan ¹L_a (293–300 nm) and the disulfide n-σ* bands (250 rm), and positive rotatory strength contributions from the tryptophan ¹L_b bands (291 nm) and the tyrosine ¹L_b bands (275 nm). Contributions to the rotatory strength of each band were approximated in terms of specific interactions between chromophores. It was found that the rotatory strength of most of the near-uv transitions arises primarily from coupling interactions involving other side-chain chromophores and amide groups which are in close proximity. Changes which are observed in the lysozyme CD spectrum on binding of tri-N-acetylglucosamine may be explained in terms of changes in the rotatory strength which result from interactions of the ¹L_a transitions of the active-site tryptophans with the acetamide groups of the inhibitor. The reasonable agreement which is found between the experimental and calculated rotatory strengths implies that the crystal conformation of lysozyme must resemble the solution conformation.     

All-valence-shell molecular orbital calculations on the chiroptical properties of substituted and unsubstituted 2,5-diketopiperazine structures

The chiroptical properties of L -3-methyl-2,5-diketopiperazine (L -alanylglycyl anhydride) are examined on a theoretical model in which the electronic wave functions are obtained from semi-empirical all-valence-shell molecular orbital calculations. The INDO molecular orbital model is used to perform SCF-MO calculations on the ground states of six conformation isomers of L -3-methyl-2,5-diketopiperazine and two chiral conformational isomers of unsubstituted 2,5-diketopiperazine. Excited-state wave functions are constructed in the virtual orbital-configuration interaction approximation. The rotatory strengths, dipole strengths, oscillator strengths, and dissymmetry factors of the first eight singlet–singlet transitions for each of the eight structures are calculated and reported. Additionally, ground-state dipole moments, net atomic charges, and the first four ionization potentials (calculated according to Koopman's theorem) are computed for each structure. The signs and the magnitudes of the rotatory strengths are found to be extremely sensitive to the conformation of the piperazine ring as well as to methyl substitution at the α carbon of the ring. Spectra–structure relationships based on the calculations reported here are discussed, and the available experimental CD data on dissymmetric 2,5-diketopiperazine are examined in terms of our theoretical results.     

Density functional theory calculations of optical rotation: employment of ADZP and its comparison with other basis sets

Density function theory calculations of frequency dependent optical rotations ([alpha]omega) for 30 rigid chiral molecules are reported. Calculations have been carried out at the sodium D line frequency, using the augmented double zeta valence quality plus polarization functions (ADZP) basis set and the BP86 nonhybrid and B3LYP hybrid functionals. Gauge-invariant atomic orbitals were used to guarantee origin-independent values of [alpha]D. Comparison between corresponding results obtained with nonhybrid and hybrid functionals as well as with theoretical optical rotations reported in the literature is done. Excited electronic states of three molecules are also discussed in light of circular dichroism spectra and B3LYP and BP86 calculated excitation energies and rotatory strengths. One verifies that the B3LYP/ADZP results are in better agreement with experiment.     

Chiroptical properties of S-proline conformational isomers

The chiroptical properties of S-proline conformational isomers are examined on a theoretical model in which electronic wave functions are obtained from semiempirical molecular orbital calculations. The CNDO/S molecular orbital model is used to perform SCF-MO calculations on ground state electronic structure and excited states are constructed in the virtual orbital-configuration interaction approximation. Electronic rotatory strengths and dipole strengths are calculated directly from the complete (but approximate) molecular electronic wave functions. Zwitterionic, cationic, and anionic S-proline structures are studied twotypes of conformational variables are represented in the calculations: (1) pyrrolidine ring conformation; and (2) rotation about the C^α-COO^? bond. Rotatory strengths are found to be somewhat sensitive to rotational isomerism about the C^α-COO^? bond, but are found to be rather insensitive to conformational changes within the pyrrolidine ring. The CD spectrum of zwitterionic S-proline down to ～160 nm appears to be well accounted for by the theoretically calculated results if conformational preferences with respect to rotation about the C^α-COO^? bond can be assumed to exist in solution media. Furthermore, spectra-structure correlations are offered for the anionic and cationic forms of S-proline in solution.     

The chromophore structures of the Pr States in plant and bacterial phytochromes

The resonance Raman spectra of the Pr state of the N-terminal 65-kDa fragment of plant phytochrome phyA have been measured and analyzed in terms of the configuration and conformation of the tetrapyrroles methine bridges. Spectra were obtained from phyA adducts reconstituted with the natural chromophore phytochromobilin as well as phycocyanobilin and its isotopomers labeled at the terminal methine bridges through (13)C/(12)C and D/H substitution. Upon comparing the resonance Raman spectra of the various phyA adducts, it was possible to identify the bands that originate from normal modes dominated by the stretching coordinates of the terminal methine bridges A-B and C-D. Quantum chemical calculations of the isolated tetrapyrroles reveal that these modes are sensitive indicators for the methine bridge configuration and conformation. For all phyA adducts, the experimental spectra of Pr including this marker band region are well reproduced by the calculated spectra obtained for the ZZZasa configuration. In contrast, there are substantial discrepancies between the experimental spectra and the spectra calculated for the ZZZssa configuration, which has been previously shown to be the chromophore geometry in the Pr state of the bacterial, biliverdin-binding phytochrome from Deinococcus radiodurans (Wagner, J. R., J. S. Brunzelle, K. T. Forest, R. D. Vierstra. 2005. Nature. 438:325-331). The results of this work, therefore, suggest that plant and bacterial (biliverdin-binding) phytochromes exhibit different structures in the parent state although the mechanism of the photoinduced reaction cycle may be quite similar.     

Theoretical CD spectrum calculations of the crown-ether aralkyl-ammonium salt complex

Rotatory strengths of the alpha-(1-naphtyl)-ethylammonium perchlorate (NEA)-phenazino-18-crown-6 ether molecular complex is determined theoretically by the coupled oscillator model and using ab initio random phase approximation (RPA) to describe local excitations on the chromophores. The computational results are compared to the experimental circular dichroism (CD) spectrum published previously. The good qualitative agreement between calculated and measured optical rotatory strengths allows one to assign the CD bands of the complex in a unique manner.     

Critical comparison of the experimental optical activity of helical polypeptides and the predictions of the molecular exciton model

The predictions of the presently accepted molecular exciton model for the optical activity of helical polypeptides are in reasonable agreement with experimental spectra in the accessible wavelength range. However, crucial verification requires the detect of a significant negative rotatory hand just below the accessible range. A computer-oriented method is utilized to obtain information concerning the inaccessible range. Optical rotatory dispersion computed by evaluation of the Kronig-Kramers integral transform from the experimentally determined circular dichroism of several helical homopolypeptides in solution are compared with the experimentally determined optical rotatory dispersion. Computed and experimental curves are congruent within an uncertainty approaching that of the experimental technique, whatever the polypeptide sample lot, side chains, and solvent. It is shown that t his agreement is not a computational or experimental artifact. These results can be interpreted in two ways: (1) that the predicted band does not exist, and (2) that the perturbation of the predicted band is being negated by other inaccessible bands in the vacuum ultraviolet. Arguments are presented to show that the first of these two possibilities is more probable.     

On the eptihermal neutron energy limit for Accelerator-Based Boron Neutron Capture Therapy (AB-BNCT): Study and impact of new energy limits

Background and purpose: Accelerator-Based Boron Neutron Capture Therapy is a radiotherapy based on compact accelerator neutron sources requiring an epithermal neutron field for tumour irradiations. Neutrons of 10 keV are considered as the maximum optimised energy to treat deep-seated tumours. We investigated, by means of Monte Carlo simulations, the epithermal range from 10 eV to 10 keV in order to optimise the maximum epithermal neutron energy as a function of the tumour depth.Methods: A Snyder head phantom was simulated and mono-energetic neutrons with 4 different incident energies were used: 10 eV, 100 eV, 1 keV and 10 keV. ¹⁰B capture rates and absorbed dose composition on every tissue were calculated to describe and compare the effects of lowering the maximum epithermal energy. The Therapeutic Gain (TG) was estimated considering the whole brain volume.Results: For tumours seated at 4 cm depth, 10 eV, 100 eV and 1 keV neutrons provided respectively 54%, 36% and 18% increase on the TG compared to 10 keV neutrons. Neutrons with energies between 10 eV and 1 keV provided higher TG than 10 keV neutrons for tumours seated up to 6.4 cm depth inside the head. The size of the tumour does not change these results.Conclusions: Using lower epithermal energy neutrons for AB-BNCT tumour irradiation could improve treatment efficacy, delivering more therapeutic dose while reducing the dose in healthy tissues. This could lead to new Beam Shape Assembly designs in order to optimise the BNCT irradiation.     

Calculations of the CD spectrum of bovine pancreatic ribonuclease

CD spectra of bovine pancreatic ribonuclease A (RNase A) and its subtilisin-modified form (RNase S) have been calculated, based upon high-resolution structures from x-ray diffraction. All known transitions in the peptide and side-chain groups, especially the aromatic and disulfide groups, have been included. Calculations have been performed with both the matrix method and with first-order perturbation theory. A newly developed method for treating the electrostatic interactions among transition charge densities and between static charge distributions and transition charge densities is used. The effects of local electrostatic fields upon the group transition energies are included for all transitions. Rotational strengths generated by the matrix method were combined with Gaussian band shapes to generate theoretical CD spectra. The calculated spectra reproduce the signs and approximate magnitudes of the near-uv CD bands of both RNase A and S. Agreement is most satisfactory for the negative 275 nm band, dominated by tyrosine contributions. In agreement with two previous studies by other workers, coupling between Tyr 73 and Tyr 115 is the single most important factor in this band. The positive band observed near 240 nm is dominated by disulfide contributions, according to our results. The far-uv CD spectrum is poorly reproduced by the calculations. The observed 208 nm band, characteristic of α-helices, is absent from the calculated spectrum, probably because the helices in RNase are short. A strong positive couplet centered near 190 nm is predicted but not observed. Possible reasons for these incorrect predictions of the current theoretical model in the far-uv are discussed. © 1997 John Wiley & Sons, Inc.     

The electronic structure of cytosine, 5-azacytosine,and 6-azacytosine

The normal and vacuum uv spectra of cytosine, 5-azacytosine, and 6-azacytosine thin layers were determined. The spectra were resolved into a series of Gaussian bands and the experimental band positions and oscillator strengths were compared with the calculated values. The calculations were carried out with the PPP-SCF-CI method. The agreement of the calculations with those of others and with the experiment is discussed.     

Effects of shape,size, and pyrene doping on electronic properties of graphene nanoflakes

Effects of size, shape, and pyrene doping on electronic properties of graphene nanoflakes (GNFs) were theoretically investigated using density functional theory method with PBE, B3PW91, and M06-2X functionals and cc-pVDZ basis set. Two shapes of zigzag GNFs, hexagonal (HGN) and rhomboidal (RGN), were considered. The energy band gap of GNF depends on shape and decreases with size. The HGN has larger band gap energy (1.23–3.96 eV) than the RGN (0.13–2.12 eV). The doping of pyrene and pyrene derivatives on both HGN and RGN was also studied. The adsorption energy of pyrene and pyrene derivatives on GNF does not depend on the shape of GNFs with energies between 21 and 27 kcal mol^?1. The substituent on pyrene enhances the binding to GNF but the strength does not depend on electron withdrawing or donating capability. The doping by pyrene and pyrene derivatives also shifts the HOMO and LUMO energies of GNFs. Both positive (destabilizing) and negative (stabilizing) shifts on HOMO and LUMO of GNFs were seen. The direction and magnitude of the shift do not follow the electron withdrawing and donating capability of pyrene substituents. However, only a slight shift was observed for doped RGN. A shift of 0.19 eV was noticed for HOMO of HGN doped with 1-aminopyrene (pyNH₂) and of 0.04 eV for LUMO of HGN doped with 1-pyrenecarboxylic acid (pyCOOH).

Open image in new window

Graphical Abstract HOMO and LUMO Energies of pyrene/pyrene derivatives doped Graphene Nanoflakes

     

Optical properties of polypeptides in the β-conformation

The rotational strengths and oscillator strengths of the nπ* band and ππ* exciton bands have been calculated for antiparallel and parallel β-structures of varying length and width. The results are compared with experiment and with previous theoretical treatments of β-structures. The generally good agreement of calculations on the antiparallel β-structure with experimental results on poly-L -lysine and poly-L -serine indicates that these systems are indeed in the antiparallel conformation. It is found that the exciton component strongest in absorption shifts to longer wavelengths as the width of an antiparallel structure increases, and it is suggested that the position of the ππ* absorption band may be a useful criterion of sheet width. The results also reconcile the linear dichroism measurements of Rosenheck and Sommer on poly-L -lysine films with an anti-parallel structure. Calculations on parallel β-structures indicate that the CD spectra of this form will be rather similar to that of the antiparallel form. However, the major absorption band in the antiparallel form is associated with a small positive CD band, while in the parallel form it coincides with a large negative CD band. Finally, it is pointed out that the large positive CD bands predicted for single-stranded parallel and antiparallel β-structures at about 200 mμ render unlikely the suggestion that random-coil polypeptides contain a substantial fraction of extended chain.     

Adsorption and dimerisation of thiol molecules on Au(111) using a Z-matrix approach in density functional theory

The adsorption energetics of methanethiolate and benzenethiolate on Au(111) have been calculated using periodic density functional theory (DFT), based on the SIESTA methodology, with an internal coordinates implementation for geometry input and structure optimisation. Both molecules are covalently bound with interaction energies of 1.85 and 1.43 eV for methanethiolate and benzenethiolate, respectively. The preferred binding site is slightly offset from the bridge site in both cases towards the fcc-hollow. The potential energy surfaces (PES) have depths of 0.36 and 0.22 eV, the hollow sites are local maxima in both cases, and there is no barrier to diffusion of the molecule at the bridge site. The corresponding dimers are weakly bound for methanethiolate and benzenethiolate, with binding energies of 0.38 and 0.16 eV, respectively, and the preferred binding geometry is with the two sulphur atoms close to adjacent atop sites. The barrier to dissociation of the dimer dimethyl disulphide is estimated to lie between 0.3 and 0.35 eV.     

Contribution of side-chain chromophores to the optical activity of proteins: Model compound studies. IV. The indole chromophore of yohimbinic acid.

The optical properties of the indole chromophore of the indole alkaloid yohimbinic acid have been investigated as a function of molecular conformation. Theoretical rotatory strengths have been calculated and compared with experimental circular dichroism spectra. Optical data that may be suitable for calculating the chiroptical properties of the near ultraviolet electronic transitions of the indole chromophores, which occur in tryptophan residues of proteins, have been developed. The far ultraviolet transitions of yohimbinic acid have also been investigated.     

Effect of side chains on the conformational energy and rotational strength of the n-pi transition for some alpha-helical poly-alpha-amino acids

Calculations of the dependence of the conformational energy and the rotational strength of the amide n–π* electronic transition (in a series of α-helical polyhel-α- amino acids with different side chains) on conformation have been carried out. The conformational energies were computed by procedures developed in this laboratory; the computation of rotational strengths was carried out by the method of Schellman and Oriel, with a slight modification. Polyamino acids with both nonpolar and polar side chains were considered; in the latter case, it was assumed that the only influence of the polar side chain was on the backbone conformation and on the electrostatic field which perturbs the amide chromophore of the backbone. Only conformations in the range of backbone dihedral angles of the right- and left-handed a-helices were considered, and the assumption of regularity (i.e., uniformity of dihedral angles in every residue) was made. The rotational strength per residue was found to vary markedly with chain length (in oligomers of up to 40 residues long); both the conformational energy per residue and the rotational strength per residue were found to vary significantly with the backbone conformation, which in turn depends on the nature of the side chain. The geometry of the hydrogen bond in the α-helical backbone is the most important factor which influences the dependence of the rotational strength on conformation. The implications of these results, for the interpretation of experimental circular dichroism and optical rotatory dispersion data, are discussed.     

FTIR study of the photoinduced processes of plant phytochrome phyA using isotope-labeled bilins and density functional theory calculations

Fourier transform infrared spectroscopy was used to analyze the chromophore structure in the parent states Pr and Pfr of plant phytochrome phyA and the respective photoproducts lumi-R and lumi-F. The spectra were obtained from phyA adducts assembled with either uniformly or selectively isotope-labeled phytochromobilin and phycocyanobilin. The interpretation of the experimental spectra is based on the spectra of chromophore models calculated by density functional theory. Global ¹³C-labeling of the tetrapyrrole allows for the discrimination between chromophore and protein bands in the Fourier transform infrared difference spectra. All infrared difference spectra display a prominent difference band attributable to a stretching mode with large contributions from the methine bridge between the inner pyrrole rings (B-C stretching). Due to mode coupling, frequencies and isotopic shifts of this mode suggest that the Pr chromophore may adopt a distorted ZZZssa or ZZZasa geometry with a twisted A-B methine bridge. The transition to lumi-R is associated with only minor changes of the amide I bands indicating limited protein structural changes during the isomerization site of the C-D methine bridge. Major protein structural changes occur upon the transition to Pfr in which the chromophore adopts a ZZEssa or ZZEasa-like state. In addition, specific interactions with the protein alter the structure of the B-C methine bridge as concluded from the substantial downshift of the respective stretching mode. These interactions are removed during the photoreaction to lumi-F (ZZE→ZZZ), which involves only small protein structural changes.     

A computer simulation of functional group contributions to free energy in water and a DPPC lipid bilayer

A series of all-atom molecular dynamics simulations has been performed to evaluate the contributions of various functional groups to the free energy of solvation in water and a dipalmitoylphospatidylcholine lipid bilayer membrane and to the free energies of solute transfer (Delta(DeltaG(o))X) from water into the ordered-chain interior of the bilayer. Free energies for mutations of the alpha-H atom in p-toluic acid to six different substituents (-CH3, -Cl, -OCH3, -CN, -OH, -COOH) were calculated by a combined thermodynamic integration and perturbation method and compared to literature results from vapor pressure measurements, partition coefficients, and membrane transport experiments. Convergence of the calculated free energies was indicated by substantial declines in standard deviations for the calculated free energies with increased simulation length, by the independence of the ensemble-averaged Boltzmann factors to simulation length, and the weak dependence of hysteresis effects on simulation length over two different simulation lengths and starting from different initial configurations. Calculated values of Delta(DeltaG(o))X correlate linearly with corresponding values obtained from lipid bilayer transport experiments with a slope of 1.1 and from measurements of partition coefficients between water and hexadecane or decadiene, with slopes of 1.1 and 0.9, respectively. Van der Waals interactions between the functional group of interest and the acyl chains in the ordered chain region account for more than 95% of the overall potential energy of interaction. These results support the view that the ordered chain region within the bilayer interior is the barrier domain for transport and that solvation interactions within this region resemble those occurring in a nonpolar hydrocarbon.     

Optical rotatory dispersion of mononucleosides and 5′-mononucleotides

Optical rotatory dispersions between 200 and 600 mμ are presented for mononucleosides and 5′-mononucleotides in neutral, acid, and in some cases alkaline solutions. All display a single Cotton effect in the ultraviolet region (above 220–240 mμ); the purine ones are negative in sign and the pyrimidine ones positive, with the crossovers (zero rotations) in most cases close to the wavelengths of their respective absorption maxima. The visible rotatory dispersions obey the one-term Drude equation, except for TMP and UMP, which show anomalous dispersions. The rotational strengths of the dichroic bands were estimated from the Cotton effect profiles; cytosine mononucleotides and mononucleotides show the strongest rotational strengths among the compounds studied.