Phenylisoserine in the gas-phase and water: Ab initio studies on neutral and zwitterion conformers

The conformational landscape of phenylisoserine (PhIS) was studied. Trial structures were generated by allowing for all combinations of single-bond rotamers. Based on the B3LYP/aug-cc-pVDZ calculations 54 conformers were found to be stable in the gas phase. The six most stable conformers were further optimized at the B3LYP/aug-cc-pVTZ and MP2/aug-cc-pVDZ levels for which characteristic intramolecular hydrogen bond types were classified. To estimate the influence of water on PhIS conformation, the IEF-PCM/B3LYP/aug-cc-pVDZ calculations were carried out and showed 51 neutral and six zwitterionic conformers to be stable in water solution. According to DFT calculations, the conformer equilibrium in the gas phase is dominated by one conformer, whereas the MP2 calculations suggest three PhIS structures to be significantly populated. Comparison of DFT and MP2 energies of all 57 structures stable in water indicates that, in practice, one zwitterionic and one neutral conformer determine the equilibrium in water. Based on the AIM calculations, we found that for the neutral conformers in vacuum and in water, d(H...B) is linearly correlated with Laplacian at the H-bond critical point.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Crystal structure and solid state 13C NMR analysis of N-(methyl 3,4,6-tri-O-acetyl-alpha, and beta-D-glucopyranosid-2-yl)-oxamide derivative of p-chloroaniline, N,N-diethylamine, N-methylaniline and N-ethylaniline

The X-ray diffraction analysis of N-(methyl 3,4,6-tri-O-acetyl-alpha-D-glucopyranosid-2-yl)-N'-p-chlorophenyloxamide (1), N-(methyl 3,4,6-tri-O-acetyl-alpha-D-glucopyranosid-2-yl)-N',N'-diethyloxamide (2), N-acetyl, N-(methyl 3,4,6-tri-O-acetyl-beta-D-glucopyranosid-2-yl), N'-methyl, N'-phenyloxamide (3), N-acetyl, N-(methyl 3,4,6-tri-O-acetyl-beta-D-glucopyranosid-2-yl), N'-ethyl, N'-phenyloxamide (4) was performed. It was found that the oxamide group in compounds 1-4 can be characterized as two structurally independent amides because there is no pi conjugation across the oxalyl OC-CO bond. Only the oxamide group of 1 is planar and adopts trans conformation stabilized as two intramolecular N-H...O hydrogen bonds.     

Resonance Raman investigation of lysine and N-acetylmethionine complexes of ferric and ferrous microperoxidase

In order to evaluate the steric and electronic effects of mixed axial ligations on the heme c structure, lysine (Lys) and N-acetylmethionine (AcMet) complexes of ferric and ferrous microperoxidase-8 (MP8(III) and MP8(II), respectively) are characterized by absorption and resonance Raman (RR) spectroscopies. Spectrophotometric titrations establish that MP8(III) binds one molecule of exogenous ligand while MP8(II) forms mono(ligated) and bis(ligated) compounds. The Soret-excited RR spectra of the six-coordinated low-spin MP8(III) complexes show that the macrocycle can adopt different structures between planar and ruffled conformations. The ferriheme c conformation is primarily determined by the ionization state of the His side chain of MP8(III) and, secondarily, by the bonding and nonbonding heme-ligand interactions. As far as the RR spectra of the MP8(II) complexes are concerned, they permit us to conclude that the mixed His/Lys and His/AcMet coordinations induce a nonplanar heme conformation, the extent of deformation again depending on the ionization state of the endogenous His ligand. In contrast, the RR spectra of the bis(Lys) and bis(AcMet) compounds are associated with a planar heme structure. When the His of MP8 is bound to heme c, the stabilization of distorted heme conformations is thus associated with constraints exerted by the Cys-Ala-Gln-Cys-His-peptide on the porphyrin macrocycle. More generally, the spectroscopic data obtained in this study can be used to predict both the axial coordination and the structure of heme in c-type cytochromes. Received: 19 January 1998 / Revised version: 23 March 1998 / Accepted: 27 March 1998     

Anomeric and rotameric preferences of glucopyranose in vacuo, water and organic solvents

Glucopyranose is the most stable form of glucose in solution. Identification of molecular structure of glucopyranose is very important because of its biological and synthetic significance; it is not an easy task because of the large number of possible configurations. Relative energies of exocyclic hydroxymethyl rotamers and α-β anomers of D-glucopyranose have been determined at the reference MP2/6-31G(d,p) level geometry by ab initio calculations at the infinite basis set limit of MP2 approach and with inclusion of CCSD(T) correction term evaluated with the aug-cc-pVDZ basis set in vacuum, water, dimethylsulfoxide, tetrahydrofurane and ethanol. The infinite basis set limit of MP2 level was determined by two point extrapolation using aug-cc-pVTZ and aug-cc-pVQZ basis sets. Solvent effects, relative energies and binding energies have been considered applying explicit calculations and implicit solvent models.     

Synthesis and conformational analysis of a locked analogue of carbovir built on a bicyclo[3.1.0]hex-2-enyl template

The synthesis and biological evaluation of a carbovir analogue (5) built on a bicyclo[3.1.0]hex-2-enyl template is described. A conformational analysis using density functional theory at the B3LYP/6-31G* level has been carried out on the rigid pseudosugar template of 5, the cyclopentene moiety of carbovir and the bicyclo[3.1.0]hex-2-yl pseudosugars of two isomeric carbonucleosides (12 and 13) containing exo- and endo-fused cyclopropane rings. The results show that while the planar configuration of the fused cyclopentane ring of compound 5 helps retain weak anti-HIV activity, the ability of the cyclopentene ring of carbovir to easily adopt a planar or puckered conformation with little energy penalty may prove to be a crucial advantage. The bicyclo[3.1.0]hex-2-yl nucleosides 12 and 13 that were inactive against HIV exhibited stiffer resistance to having a planar, fused cyclopentane moiety.     

Impact of cis-proline analogs on peptide conformation

The beta-turn is a common motif in both proteins and peptides and often a recognition site in protein interactions. A beta-turn of four sequential residues reverses the direction of the peptide chain and is classified by the phi and psi backbone torsional angles of residues i + 1 and i + 2. The type VI turn usually contains a proline with a cis-amide bond at residue i + 2. Cis-proline analogs that constrain the peptide to adopt a type VI turn led to peptidomimetics with enhanced activity or metabolic stability. To compare the impact of different analogs on amide cis-trans isomerism and peptide conformation, the conformational preference for the cis-amide bond and the type VI turn was investigated at the MP2/6-31+G** level of theory in water (polarizable continuum water model). Analogs stabilize the cis-amide conformations through different mechanisms: (1) 5-alkylproline, with bulky hydrocarbon substituent on the C(delta) of proline, increases the cis-amide population through steric hindrance between the alkyl substituent and the N-terminal residues; (2) oxaproline or thioproline, the oxazolidine- or thiazolidine-derived proline analog, favors interactions between the dipole of the heterocyclic ring and the preceding carbonyl oxygen; and (3) azaproline, containing a nitrogen atom in place of the C(alpha) of proline, prefers the cis-amide bond by lone-pair repulsion between the alpha-nitrogen and the preceding carbonyl oxygen. Preference for the cis conformation was augmented by combining different modifications within a single proline. Azaproline and its derivatives are most effective in stabilizing cis-amide bonds without introducing additional steric bulk to compromise receptor interactions.     

Synthesis and structural studies of N-p-toluensulfonyl cyclodipeptides

In view of the chemical and structural interest of cyclopeptides bearing an electron withdrawing substituent directly bonded at the amide nitrogen atom, the two N-p-toluensulfonyl (N-tosyl) derivatives cyclo[-Phe(Tos)-D-Phe-] (I) and cyclo[-Phe(Tos)-D-Pro-] (II) have been synthesized and their stereochemistry defined. The molecular structure of I, as determined by x-ray diffraction analysis, is reported together with 1H-nmr parameters indicating the preferred rotameric conformation in chloroform solution. The N-tosyl group alters the geometry of the cyclodipeptide ring by lengthening both the N-C bonds departing from the tosylated nitrogen and reducing the corresponding ring angle. The 6-membered peptide ring adopts an unusual "sofa" conformation with the Tos-Phe C alpha atom deviating 0.230(3) A out of the mean plane of the other five ring atoms. One of the two S-O bonds forms a planar system that involves the tosylated nitrogen and the corresponding amide carbonyl. In the crystal, both the benzylic side chains are folded over the heterocyclic ring, whereas in chloroform solution, the benzylic side chain of the D-Phe prefers an extended conformation.     

QSAR study of dual cyclooxygenase and 5-lipoxygenase inhibitors 2,6-di-tert-butylphenol derivatives

The dual or selective ability of 24 derived mono- and 2,6-di-tert-butylphenols (DTBP) to act as inhibitors of cyclooxygenase (COX) and/or 5-lipoxygenase (LOX) enzymes is investigated. Firstly, we explored the conformational variability of the compounds. It is found that dual inhibitors can adopt four minimum energy conformations: cis or trans, depending on the orientation of the carbonyl oxygen atom (localized in the para position) relative to the hydroxyl hydrogen, and alpha or beta, depending on whether the carbonyl oxygen is below or above the phenyl plane. The possible bioactive conformations are selected by molecular superimposition to the optimised structure of tebufelone, a dual inhibitor in the clinical trial phase. From this selected conformation, different molecular parameters were calculated and correlated with both COX/LOX inhibitory activities. The MEP and GRID maps for different probes (hydrogen bond donor/acceptor, hydrophobic and ferric/ferrous interaction areas) are represented and discussed. The results point out the importance of the hydrogen donation and the hydrophobic properties in the COX inhibition whereas, for LOX inhibition, a redox mechanism might be involved. Finally, the predictive ability of the proposed QSAR equations is tested analysing a set of selective COX or LOX inhibitors.     

Accurate micro- and macro- gas phase basicities of hydroxyl-radical-modified pyrimidines estimated by advanced quantum chemistry methods

The gas phase basicities (GPB) of 16 hydroxyl-radical-derived analogues of model pyrimidine nucleosides were analysed at B3LYP/aug-cc-pVDZ and G3MP2B3 levels of theory. All possible tautomeric equilibriums of neutral and protonated forms of pyrimidine analogues were taken into account. The oxidation of pyrimidines usually reduces the GPB values with respect to canonical values. The only exception was observed in the case of 5,6-dihydroxycytidine, which is more basic than model cytidine if both macro- and micro-scopic measures of basicities are used for comparison. In addition, 6-hydroxycytidine is characterised more by the basic character of the O₂ atom compared to the more basic centre of cytidine, despite the fact that the GPB values of both these compounds are almost identical. Although the B3LYP/aug-cc-pVDZ approach seems to be an accurate and robust method for GPB estimation, the microscopic protonation properties are much more sensitive to this method since the difference in energy between some tautomers is often less than 1 kcal/mol with method-dependent succession. In such cases, G3MP2B3 is recommended.     

Structures of nitroreductase in three states: effects of inhibitor binding and reduction.

The crystal structure of the nitroreductase enzyme from Enterobacter cloacae has been determined for the oxidized form in separate complexes with benzoate and acetate inhibitors and for the two-electron reduced form. Nitroreductase is a member of a group of enzymes that reduce a broad range of nitroaromatic compounds and has potential uses in chemotherapy and bioremediation. The monomers of the nitroreductase dimer adopt an alpha+beta fold and together bind two flavin mononucleotide prosthetic groups at the dimer interface. In the oxidized enzyme, the flavin ring system adopts a strongly bent (16 degrees ) conformation, and the bend increases (25 degrees ) in the reduced form of the enzyme, roughly the conformation predicted for reduced flavin free in solution. Because free oxidized flavin is planar, the induced bend in the oxidized enzyme may favor reduction, and it may also account for the characteristic inability of the enzyme to stabilize the one electron-reduced semiquinone flavin, which is also planar. Both inhibitors bind over the pyrimidine and central rings of the flavin in partially overlapping sites. Comparison of the two inhibitor complexes shows that a portion of helix H6 can flex to accommodate the differently sized inhibitors suggesting a mechanism for accommodating varied substrates.     

Can Li2F2 cluster be formed by LiF2/Li2F–Li/F interactions? An ab initio investigation

Superatoms can be considered as potential building blocks not only of non-traditional species but also of many traditional species. This idea is demonstrated by considering Li₂F₂, traditionally a dimer of LiF which can be formed either by LiF₂–Li interaction or by Li₂F–F interaction using MP2/aug-cc-pVDZ method. The existence of superatomic motifs is discussed in the resulting linear as well as square conformers of Li₂F₂ which are compared on the basis of the highest occupied molecular orbitals and vibrational properties. Finally, we have established that formation of Li₂F₂ by these interactions is exothermic at CCSD(T)//MP2/aug-cc-pVDZ level. This study is expected to provide further insights into superatomic interactions and formations of traditional salts.     

1H-NMR studies of receptor-selective substance P analogues reveal distinct predominant conformations in DMSO-d6

Proton nmr parameters are reported for DMSO-d6 solutions of two receptor-selective substance P analogues: Ac[Arg6,Pro9]SP6-11, which is selective for the NK-1 (SP-P) receptor and [pGlu6,N-MePhe8]SP6-11, which selectively activates the NK-3 (SP-N) receptor. Full peak assignments of both analogues were obtained by COSY experiments. The chemical shifts, coupling constants, and temperature coefficients of amide proton chemical shifts as well as NOESY effects and calculated side-chain rotamer populations of Phe side chains are reported for both peptides. Analysis of coupling constants and temperature coefficients together with the nuclear Overhauser enhancement spectroscopy effects suggest that Ac[Arg6,Pro9]SP6-11 has a trans configuration about the Phe8-Pro9 amide bond and the preferred conformation of this analogue has a type I beta-turn. The nmr data for [pGlu6,N-MePhe8]SP6-11 suggest that this peptide exists as a mixture of cis-trans isomers in which the cis isomer can preferably adopt a type VI beta-turn conformation, and the trans isomer can adopt a gamma-turn conformation. There are indications that the two last turns are stabilized by a hydrogen bond between the syn carboxamide proton and the pGlu ring carbonyl.     

High-resolution A-DNA crystal structures of d(AGGGGCCCCT). An A-DNA model of poly(dG) x poly(dC).

A-DNA conformation is favored by guanine-rich sequences, such as (dG)n x (dC)n, or under low-humidity conditions. Earlier A-DNA crystal structures revealed some conformational variations which may be the result of sequence-dependent effects and/or crystal packing forces. Here we report the high-resolution crystal structure of d(AGGGGCCCCT) in two crystal forms (either in the P212121 or the P6122 space group) to gain insights into the conformation and dynamics of the (dG)n x (dC)n sequence. The P212121 form has been analyzed using data to 1.1 A resolution by the anisotropic temperature factor refinement procedure of the SHELX97 program. Such analysis affords us with the detailed geometric, conformational and motional property of an A-DNA structure. The backbone torsional angles fall in a narrow range, except for the alpha/gamma angles which have two distinct combinations (gauche-/gauche+ or trans/trans). An A-DNA model of poly(dG) x poly(dC) has been constructed using the conformational parameters derived from the crystal structure of the P212121 form. In the crystal structure of the P6122 space group, the central eight base pairs of the decamer adopt A-DNA conformation with the two terminal nucleotides flipped out to form base pairs with the neighboring nucleotides. Comparison of the A-DNA structure of the same sequence from two different crystal forms, reinforced the conclusion that molecules crystallized in the same space group have a more similar conformation, whereas the same molecule crystallized in different space groups has different (local) conformations.     

Ab initio conformational analysis of nucleic acid components: intrinsic energetic contributions to nucleic acid structure and dynamics

In recent years, the use of high-level ab initio calculations has allowed for the intrinsic conformational properties of nucleic acid building blocks to be revisited. This has provided new insights into the intrinsic conformational energetics of these compounds and its relationship to nucleic acids structure and dynamics. In this article we review recent developments and present new results. New data include comparison of various levels of theory on conformational properties of nucleic acid building blocks, calculations on the abasic sugar, known to occur in vivo in DNA, on the TA conformation of DNA observed in the complex with the TATA box binding protein, and on inosine. Tests of the Hartree-Fock (HF), second-order M?ller-Plesset (MP2), and Density Functional Theory/Becke3, Lee, Yang and Par (DFT/B3LYP) levels of theory show the overall shape of backbone torsional energy profiles (for gamma, epsilon, and chi) to be similar for the different levels, though some systematic differences are identified between the MP2 and DFT/B3LYP profiles. The east pseudorotation energy barrier in deoxyribonucleosides is also sensitive to the level of theory, with the HF and DFT/B3LYP east barriers being significantly lower (approximately 2.5 kcal/mol) than the MP2 counterpart (approximately 4.0 kcal/mol). Additional calculations at various levels of theory suggest that the east barrier in deoxyribonucleosides is between 3.0 and 4.0 kcal/mol. In the abasic sugar, the west pseudorotation energy barrier is found to be slightly lower than the east barrier and the south pucker is favored more than in standard nucleosides. Results on the TA conformation suggest that, at the nucleoside level, this conformation is significantly destabilized relative to the global energy minimum, or relative to the A- and B-DNA conformations. Deoxyribocytosine would destabilize the TA conformation more than other bases relative to the A-DNA conformation, but not relative to the B-DNA conformation.     

Simulations of dynamical properties of a Michaelis complex: porcine pancreatic elastase and the hexapeptide, Thr-Pro-n Val-Leu-Tyr-Thr

Molecular dynamic simulations (30ps) of the Michaelis complex of hexapeptide (Thr-Pro-nVal-Leu-Tyr-Thr) bound to porcine pancreatic elastase (PPE) hydrated by about 2000 water molecules have been performed using the AMBER 3.0 program package. Dynamical properties of the conformation of the active site have been examined. A comparison with previously reported simulations of native PPE shows that after the substrate is bound, the catalytically crucial H-bond between O gamma-H group of (Ser 195) and nitrogen N epsilon (His 57) is more readily formed. These results show, however, that the H-bond does not adopt the most favorable conformation. The O gamma-H group of Ser 195 has a statistical preference for an attractive interaction with the O = C carbonyl (Ser 214) rather than the nitrogen N epsilon (His 57).     

Structure and conformation of peptides containing the sulphonamide junction. I. N-acetyl-tauryl-L-phenylalanine methyl ester

N-acetyl-tauryl-L-phenylalanine methyl ester 1 has been synthesized. The crystal structure and molecular conformation of 1 have been determined. Crystals are monoclinic, space group P2(1) with a = 5.088(2), b = 17.112(17), c = 9.581(6) A, beta = 92.34(4) degrees, Z = 2. The structure has been solved by direct methods and refined to R = 0.043 for 2279 reflections with I greater than 1.5 sigma(I). The sulphonamide junction maintains the peptide backbone folded with Tau and Phe C alpha atoms in a cisoidal arrangement, the torsion angle around the S-N bond being 65.4 degrees. In this conformation the p-orbital of the sulphonamide nitrogen lies in the region of the plane bisecting the O-S-O angle, thus favouring d pi-p pi interactions between nitrogen and sulphur atoms. The S-N bond with a length of 1.618 A has significant pi-bond character. The CO-NH is planar and adopts trans conformation. The Tau residue is extended with the Tau-C1 alpha-Ca beta bond anti-periplanar to the S-N bond. The Phe side chain conformation corresponds to the statistically most favoured g- rotamer and exhibits a chi 1 torsion angle of -67.5 degrees. The packing is characterized by intermolecular H-bonds which the Tau and Phe NH groups form with the acetyl carbonyl and one of the two sulphonamide oxygens, respectively.     

Synthesis and Conformational Analysis of a Locked Analogue of Carbovir Built on a Bicyclo[3.1.0]-hex-2-enyl Template

Abstract

The synthesis and biological evaluation of a carbovir analogue (5) built on a bicyclo[3.1.0]hex-2-enyl template is described. A conformational analysis using density functional theory at the B3LYP/6-31G* level has been carried out on the rigid pseudosugar template of 5, the cyclopentene moiety of carbovir and the bicyclo[3.1.0]hex-2-yl pseudosugars of two isomeric carbonucleosides (12 and 13) containing exo- and endo-fused cyclopropane rings. The results show that while the planar configuration of the fused cyclopentane ring of compound 5 helps retain weak anti-HIV activity, the ability of the cyclopentene ring of carbovir to easily adopt a planar or puckered conformation with little energy penalty may prove to be a crucial advantage. The bicyclo[3.1.0]hex-2-yl nucleosides 12 and 13 that were inactive against HIV exhibited stiffer resistance to having a planar, fused cyclopentane moiety.     

Ab Initio Quantum Mechanics Analysis of Imidazole C-H…O Water Hydrogen Bonding and a Molecular Mechanics Forcefield Correction

Abstract

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but ‘activated’ C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H…O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H…0 and C—O. In general, these calculations indicate that C-H…0 interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ//MP2-FC/6-31++G(2d,2p) level of theory of a water C-H…O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H…O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole C_e?H…O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole C_e?H…O water and imidazole N_d?H…O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H…O water interaction, these forcefields do not adequately account for the imidazole C_e?H…O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H…O hydrogen bonded complex.     

How does modification of adenine by hydroxyl radical influence the stability and the nature of stacking interactions in adenine-cytosine complex?

This study reports on ab initio calculations of adenine - cytosine complexes in two different context alignments appearing in B-DNA. The influence of adenine modification by hydroxyl radical on the stability of the complexes is also discussed. The analysis was performed on over 40 crystallographic structures for each of the sequence contexts. In most cases, modification of adenine by hydroxyl radical leads to less negative intermolecular interaction energies. The issue of the influence of alteration of structural base step parameters on the stability of modified and unmodified adenine - cytosine complexes is also addressed. Analysis of the dependence of intermolecular interaction energy on base step parameters reveals that for twist and shift modification of adenine by hydroxyl radical leads to quite different interaction energy profiles in comparison with unmodified complexes. In order to elucidate the physical origins of this phenomenon, i.e. to analyze how the modification of adenine by hydroxyl radical is reflected in the change of intermolecular interaction energy components, a variational-perturbational decomposition scheme was applied at the MP2/aug-cc-pVDZ level of theory.