Solvent effect on lipase enantioselectivity. Evidence for the presence of two thermodynamic states

The enantioselectivity of lipase-catalyzed kinetic resolutions has been measured at various temperatures in binary mixtures of solvents. Varying the solvent composition and temperature had a profound effect on the enantiomeric ratio. The values for delta delta H(R-S)(#) and delta delta S(R-S)(#), calculated from the E values measured at various temperatures, were estimated as a function of the solvent composition. By plotting delta delta H(R-S)(#) versus delta delta S(R-S)(#) as a function of the solvent composition, an extreme was observed. The resulting "hairpin-type" enthalpy-entropy compensation plots can be described by assuming the presence of two thermodynamically distinct physical states, displaying different enantioselectivities, that are in equilibrium with one another. Changing the solvent composition results in a change in the equilibrium constant K(eq) for the two states. The intriguing bell-shaped curves of the enantioselectivity versus solvent composition observed for lipase-catalyzed kinetic resolutions can be described assuming a linear correlation for the logarithm of K(eq) and the solvent composition. Thus, a simulation of the two-state model adequately describes the solvent effects found for lipase-catalyzed kinetic resolutions in binary mixtures of solvents and possibly in series of homologous organic solvents.     

Temperature and pressure effects on conformational equilibria of alanine dipeptide in aqueous solution

We investigated the temperature and pressure effects on conformational equilibria of N-acetyl-L-alanine-N'-methylamide (AAlaMA) in aqueous solution by Raman spectroscopy. Scattering intensities in the skeletal stretching mode of AAlaMA in aqueous solution were decomposed into some component bands by the spectra analysis. Our results indicate that each component band for AAlaMA adopts not only the P(II) and alpha(R) conformations but also the C(7eq) conformation. From temperature and pressure dependencies of the band intensities, we determined the enthalpy differences and the volume differences between the conformers. The C(7eq) conformer is enthalpically most stable due to the intramolecular hydrogen bond. The partial molar volume of the C(7eq) conformer is the smallest through the solvent-exclusion effect rather than the solute-solvent electrostatic interaction effect.     

Conformational investigation of alpha,beta-dehydropeptides. XV: N-acetyl-alpha,beta-dehydroamino acid N 'N '-dimethylamides: conformational properties from infrared and theoretical studies.

The FTIR spectra were analysed in the region of the nu(s)(N-H), AI(C=O) and nu(s)(Calpha=Cbeta) bands for a series of Ac-DeltaXaa-NMe2, where DeltaXaa = DeltaAla, (Z)-DeltaAbu, (Z)-DeltaLeu, (Z)-DeltaPhe and DeltaVal, to determine a predominant solution conformation of these alpha,beta-dehydropeptide-related molecules. Measurements were taken in CCl4, DCM and MeCN solutions. In the same way, spectra of saturated analogues Ac-Xaa-NMe2, where Xaa = Ala, Abu, Leu, Phe and Val, were investigated. To help interpret the spectroscopic results, conformational maps were calculated by the B3LYP/6-31+G** method. Also, the relative energies of all conformers of the dehydro compounds in vacuo as well as in the studied solvents in addition to the theoretical IR frequencies of these conformers were calculated. For comparison, molecules of two saturated analogues, Ac-L-Ala-NMe2 and Ac-L-Phe-NMe2, were calculated in a similar way. Both unsaturated and saturated compounds, which have an aliphatic side chain, occur in CCl4 and DCM mainly as a mixture of extended conformers with the C5 H-bond and open conformers. As solvent polarity increases, participation of the open conformers also increases, and in MeCN, the model amides are almost exclusively in the open form, except Ac-DeltaAla-NMe2, which shows a small amount of the H-bonded conformer. Ac-DeltaAla-NMe2 and Ac-DeltaAbu-NMe2 have stronger C5 hydrogen bonds than those of their saturated counterparts. As the calculations indicate, the open conformation of the unsaturated amides is conformer H/F with phi, psi -44 +/- 5 degrees, 127 +/- 4 degrees. This is the second lowest in energy conformer in vacuo and in CCl4 and the lowest one in more polar solvents. The open conformation of Ac-L-Ala-NMe2 constitutes conformer C with phi, psi -101.5 degrees, 112.7 degrees. For Ac-DeltaAla-NMe2 and Ac-DeltaAbu-NMe2, FTIR also reveals the presence of a third conformer. Calculations indicate that is the semiextended conformer D with the N1-H1...N2 hydrogen bond/contact. In all solvents, Ac-L-Phe-NMe2 and Ac-(Z)-DeltaPhe-NMe2 show only the extended E and the open H/F, respectively. In both there is an amide/pi(Ph) interaction.     

Chiroptical properties of 2,2’‐bioxirane

The two enantiomers of 2,2′‐bioxirane were synthesized, and their chiroptical properties were thoroughly investigated in various solvents by polarimetry, vibrational circular dichroism (VCD), and Raman optical activity (ROA). Density functional theory (DFT) calculations at the B3LYP/aug‐cc‐pVTZ level revealed the presence of three conformers (G₊, G_?, and cis) with Gibbs populations of 51, 44, and 5% for the isolated molecule, respectively. The population ratios of the two main conformers were modified for solvents exhibiting higher dielectric constants (G_? form decreases whereas G₊ form increases). The behavior of the specific optical rotation values with the different solvents was correctly reproduced by time‐dependent DFT calculations using the polarizable continuum model (PCM), except for the benzene for which explicit solvent model should be necessary. Finally, VCD and ROA spectra were perfectly reproduced by the DFT/PCM calculations for the Boltzmann‐averaged G₊ and G_? conformers.     

Theoretical studies on the conformation of saccharides. IV. Solvent effect on the stability of β-maltose conformers

The conformational equilibria of four β-maltose conformers have been studied theoretically in 12 solvents. The stability of the conformers in dilute solution has been compared by using the continuum reaction field method. The solvation energy consists of electrostatic, dispersion, and cavity terms, which have been determined from the calculated properties of β-maltose and physicochemical properties of solvents. The calculated population of four conformers significantly depends on the solvent (e.g., in dioxane, M1:M2:M3:M4 = 53.3: 20.3:17.7:8.7; in dimethyl sulfoxide, 40.1:21.8:22.8:15.3; and in water, 25.7:17.5:26.3:30.5) and was found to correlate with experimentally observed data. The results obtained indicate that the conformation adopted by β-maltose in the crystalline form is not the one preferred in solution. The roles of the individual contributions to the solvation energy have been analyzed. Based on the determined abundance of conformers, averaged residual optical activity and vicinal carbon–proton coupling constants have been calculated and discussed from the point of view of the solution behavior of β-maltose.     

MINDO/3 calculations of the conformation and carcinogenicity of epoxy-metabolites of aromatic hydrocarbons: 7,8-dihydroxy-9,10-oxy-7,8,9,10-tetrahydrobenzo(a)pyrene

Quantum mechanical calculations in the MINDO/3 approximation were performed on the four conformations of the alicyclic moiety of 7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo(a)pyrene.Total charge and frontier orbital densities show that attack by nucleophiles will occur predominantly at Position 10 of all configurations of the dihydroxyepoxybenzo(a)pyrene.The calculations show the cis diastereomer to be more stable than the trans, although no evidence for hydrogen bonding in the (ax, ax′) conformer was found.On the basis of the results obtained for the stability of the various conformers, a model is proposed to explain the higher carcinogenicity of the trans isomer as compared to the cis. Such a model implies the formation of an intercalation complex between the diol epoxide metabolite and nucleic acids.     

Solution conformation of a model hexapeptide containing RGD sequence.

The solution conformation of a model hexapeptide Asp-Arg-Gly-Asp-Ser-Gly (DRGDSG) containing the RGD sequence has been studied in DMSO-d6 as well as in aqueous solution (H2O:D2O/90:10%) by 1H NMR spectroscopy. The unambiguous identification of spin systems of various amino acid residues and sequence specific assignment of all proton resonances was achieved by a combination of two dimensional COSY and NOESY experiments. The temperature coefficient data of the amide proton chemical shifts in conjunction with the vicinal coupling constants, i.e. 3JNH-C alpha H, NOESY and ROESY results indicate that the peptide in both the solvents exists in a blend of conformers with beta-sheet like extended backbone structure and folded conformations. The folded conformers do not appear to be stabilised by intramolecular hydrogen bonding. Our results are consistent with the flexibility of RGD segment observed in the NMR studies on the protein echistatin containing the RGD motif (references 23-25).     

Conformational analysis of genotoxic benzo[a]pyrene-7,8-dione-duplex DNA adducts using a molecular dynamics method (II)

The conformations of the benzo[a]pyrene-7,8-quinone (BPQ) modified oligonucleotide were investigated using molecular dynamic simulation. In the initial structures, the central guanine base was modified with BPQ resulting in the formation of four structurally distinguishable 10-(N2-deoxyguanosyl)-9,10-dihydro-9-hydroxy benzo[a]pyrene-7,8-dione adducts (BPQ-G3,4). Each of the oligonucleotide adduct consisted of two conformers, namely syn and anti conformations, depending on the rotation around the glycosidic bond between BPQ and the guanine base. The results revealed that the BPQ moiety was located in the major groove for all four syn conformers. The relative energies of these conformers were high, and the backbone largely deviated from the B-form. On the other hand, BPQ was located in the minor groove with relatively low energies, and backbone was retained in all of the anti conformer cases. The most conceivable BPQ-modified double stranded oligonucleotide structure was proposed from the energy calculation and the structural analysis.     

Solvent‐dependent inversion of circular dichroism signal in naproxen: An unusual effect!

The electronic circular dichroism (ECD) spectra of naproxen enantiomers were studied as a function of solvents using experimental (circular dichroism) and theoretical (time‐dependent density functional theory) approaches. The (R)‐ and (S)‐naproxen enantiomers presented an unusual inversion in their ECD signals in the presence of ethanol and water when compared with polar aprotic solvents such as acetonitrile. From a practical point of view, these findings deserve great attention because these solvents are widely used for high‐performance liquid chromatography analysis in quality control of chiral pharmaceutical drugs. This is particularly relevant to naproxen because the (S)‐naproxen has anti‐inflammatory properties, whereas (R)‐naproxen is hepatotoxic. A time‐dependent density functional theory computer simulation was conducted to investigate the signal inversion using the solvation model based on density, a reparameterization of polarized continuum model. Electronic circular dichroism signals of conformers were calculated by computer simulation and their contribution to the combined spectra obtained according to Boltzmann weighting. It was found that the experimentally observed ECD signal inversion can be associated with the minor or major contribution of different conformers of naproxen.     

Solvent effect on the stability of mannobiose conformers

The conformational equilibria of seven methyl β-D -mannobioside conformers have been studied theoretically in five solvents. The structure of each individual conformer has been refined by the PCILO quantum-chemical method from the seven distinct low-energy regions determined from (Φ, Ψ) maps calculated by a potential function method. The stability of the conformers in dilute solution has been evaluated by using a method that consists of electrostatic, dispersion, and cavity terms. The calculated abundance of conformers depends on the solvent and results indicate that the preponderant conformer in the solution may not be the one adopted by mannobiose in the crystalline form. Based on the determined abundance of conformers, thermodynamically averaged nmr parameters, dipole moment, and linkage rotation have been calculated. The solvation behavior of methyl β-D -mannobioside is compared to those previously estimated for cellobiose and maltose.     

Crystal and molecular structure of two geometrically restricted chemotactic tripeptides, analogues of formyl-methionine-leucine-phenylalanine

The crystal structures of HCO-Met-Leu-Phe-OC(CH3)3, (CH25H39N3O5S), fMLP-OtBu, and HCO-Met psi [CSNH]-Leu-Phe-OCH3, (C22H33N3O4S2), fMS LP-OMe, have been determined by single crystal X-ray diffraction, and their conformational properties investigated by molecular mechanics energy calculations. Crystals of fMLP-OtBu are monoclinic, space group P2(1), a = 12.027(4), b = 9.492(3), c = 12.660(4) A, beta = 101.99(3) degrees, Z = 2; those of fMS LP-OMe are orthorhombic, space group P2(1)2(1)2(1), a = 7.130(1), b = 12.097(2), c = 31.060(5) A, Z = 4. The first compounds fMLP-OtBu is the t-butyl ester of the tripeptide fMLP that represents one of the most potent compounds in inducing the lysozyme release from human neutrophils that reflects the chemotactic activity. From the crystal structure, it is shown that the orientation of the phenylalanine side chain is largely affected by the presence of the bulky group. fMSLP-OMe was shown to be inactive after thionation of the methionine residue in the original tripeptide. Nevertheless, the crystal structure does not reveal any influence of the presence of the thionated peptidic bond on the backbone conformation. The X-ray results have been used to generate parameters for empirical energy calculations. Subsequently, a strategy based on random generation of conformations followed by energy-minimization was applied to investigate the conformational space of thiopeptides, in comparison with normal peptides. From molecular free energy calculations, it is shown that the main influence of the introduction of a thioamide bond on the molecular structure is to prevent the existence of C7(eq) conformations involving the thiomethionine residue. Consequently, a larger number of conformers are found to form intramolecular hydrogen bonds involving the formyl group, reducing its availability to interact with the receptor. For the first time, the theoretical prediction of the existence of C7eq conformations for fMLP is made. The resulting conformers are compared to previously active structures of these chemotactic agents.     

Photophysical properties of 2-amino-9,10-anthraquinone: evidence for structural changes in the molecule with solvent polarity.

The photophysical properties of 2-amino-9,10-anthraquinone (2AAQ) have been investigated in different solvents and solvent mixtures and correlated with the Lippert-Mataga solvent polarity parameter, Deltaf. In the low solvent polarity region with Deltaf < ca. 0.1, the dye shows unusually high fluorescence quantum yields (Phif) and lifetimes (tauf) in comparison to those in other solvents of medium to high polarities. Similarly, the radiative rate constants (kf) are relatively lower and the non-radiative rate constants (knr) are relatively higher in the low polarity solvents in comparison to those in the medium to high polarity solvents. The current results have been rationalized assuming that the dye adopts different structural forms below and above the Deltaf value of approximately 0.1. It is inferred that in the low solvent polarity region the dye exists in a non-planar structure, with its 2-NH2 plane away from that of the 9,10-anthraquinone moiety in the ground state. In solvents of medium to high polarities, the dye exists in a polar intramolecular charge transfer (ICT) structure, where the amino lone pair of the 2-NH2 group is in strong resonance with the anthraquinone pi-cloud in the ground state. In all the solvents, however the dye is inferred to exist in the ICT structure in its excited (S1) state. Supportive evidence for the above hypothesis has been obtained from the solvent polarity effect on the Stokes' shifts for the dye. Quantum chemical studies on the structures of 2AAQ dye in the gas phase also give qualitative support for the inferences drawn from the photophysical properties of the dye in different solvents.     

CD conformational studies on synthetic peptides encompassing the processing domain of the ocytocin/neurophysin precursor

Synthetic peptides of different size, reproducing the proteolytic processing site of proocytocin, were studied by CD under several experimental conditions in order to ascertain the ability of different solvents to stabilize secondary structural motifs, such as α-helix tracts and β-turns. A combination of deconvolution methods and empirical calculations subtracting the contributions due to unordered structures from the spectra suggests that in solution (a) mainly two distinct families of ordered conformers containing structurally different β-turns are present, (b) the relative stability of the different conformers depends from the nature of the solvent, and (c) in the case of the larger peptides, a population containing an α-helical conformation is also present. From the biological point of view the presence of at least two families of ordered conformers could be in line with current theories assuming that the catalytic effect of the receptor microenvironment may be determinant in shifting the equilibrium toward the active conformation. © 1997 John Wiley & Sons, Inc. Biopoly 41 : 461–479, 1997     

Solvent effect on the physical quenching of singlet molecular oxygen by p-quinones.

Rate constants for the interaction between singlet molecular oxygen [O2(1 delta g)] and the p-quinones 1,4-benzoquinone (BQ), duroquinone (DQ), 9,10-anthraquinone (AQ) and 1,8-dihydroxy-9,10-anthraquinone (OHAQ) are reported for several solvents at room temperature. The solvent effect on the total quenching rate constant (kt) was analysed employing the semiempirical solvatochromic equation proposed by Kamlet and Taft. The higher values of kt (2-7 x 10(7) M(-1) s(-1)) were obtained when the hydrogen-bond donor solvent ability is increased (higher alpha parameter values). The results indicate the importance of specific solvent interactions in governing the rates of the quenching.     

Relaxation kinetics of ribonuclease T1 binding with guanosine and 3'-GMP.

Temperature-jump relaxation kinetic studies were undertaken at 25 degrees C with ribonuclease T1 (RNase T1) alone and in the presence of guanosine (Guo) and 3'-guanylic acid (3'-GMP). No relaxations were observed in the absence of ligands and only one process was observed in their presence which reflected a simple on-off reaction in both cases. Apparent association rate constants, k(on), and dissociation rate constants, k(off), were evaluated at several pH values and their ratios, k(on)/k(off), were contrasted with independently determined values of the equilibrium association constant, Ka(eq). The value of k(on)/k(off) for Guo was significantly greater than Ka(eq), whereas Ka(eq) was significantly greater than k(on)/k(off) for 3'-GMP. The simplest interpretation of the result for Guo is that free RNase T1 undergoes a relatively slow undetected isomerization and Guo can bind only with one isomer. 3'-GMP can be considered to bind with the same preference, but in this case the initial enzyme complex undergoes a relatively slow undetected isomerization. These results are consistent with a recent NMR study which suggested that RNase T1 binding with Guo and 3'-GMP are coupled to slow exchange processes in a ligand dependent manner (Shimada, I. and Inagaki, F. (1990) Biochemistry 29, 757-764). It is tentatively concluded that binding of Guo and 3'-GMP at the active site of RNase T1 is limited to a sub-population of conformers involving the base-recognition site and that the phosphomonoester group of the nucleotide can engage in additional conformationally linked interactions at the adjacent catalytic site.     

Infrared spectra of carbonyl hemoglobins: characterization of dynamic heme pocket conformers

The infrared spectra for carbon monoxide complexed to hemoglobins were examined in the C-O stretch region. Deconvolution of the spectra requires four bands and supports the presence of four distinct conformers at the ligand binding site. Most typical hemoglobins exhibit only one predominant conformer for each subunit represented by a band at 1951 cm-1 in contrast to myoglobins, which typically exist in two major conformations. Several hemoglobins with an enlarged heme pocket are shown to shift the C-O frequency into the higher frequency conformer regions. Many factors, including pH, temperature, solvents, and divalent metals, are also shown to be capable of expanding the heme pocket. Only very specific structural changes that can reduce the size of the heme pocket will result in the lower frequency conformers. The weighted averages of the multiple CO vibrational frequencies are linearly related to the single 13CO NMR chemical shift values and to the exponential of fast CO on-rates. Conformer interconversion occurs at a rate greater than 10(4) s-1. The infrared C-O stretch spectra provide qualitative and quantitative information on the structural dynamics, stability, and ligand binding properties of hemoglobins.     

Analysis of the conformational profile of trishomocubane amino acid dipeptide

4-Amino-(D3)-trishomocubane-4-carboxylic acid is a constrained alpha-amino acid residue that exhibits promising conformational characteristics, i.e., helical and beta-turns. As part of the development of conformational guidelines for the design of peptides and protein surrogates, the conformational energy calculations on trishomocubane using molecular mechanics and ab initio methods are presented. The C(alpha) carbon of trishomocubane forms part of the cyclic structure, and consequently a peptidic environment was simulated with an acetyl group on its N-terminus and a methylamide group on its C-terminus. Ramachandran maps computed at the molecular mechanics level using the standard AMBER (parm94) force field libraries compared reasonably well with the corresponding maps computed at the Hartree Fock level, using the 6-31G* basis set. Trishomocubane peptide (Ac-Tris-NHMe) is characterized by four low energy conformers corresponding to the C7ax, C7eq, 3(10), and alpha(L) helical structures.     

Crystal structure, conformation, and potential energy calculations of the chemotactic peptide N-formyl-L-Met-L-Leu-L-Phe-OMe

The tripeptide N-formyl-L-Met-L-Leu-L-Phe-OMe (FMLP-OMe) crystallizes in the orthorhombic system, space group P 2(1)2(1)2(1), with the following unit-cell parameters: a = 21.727, b = 21.836, c = 5.133 A, Z = 4. The structure has been solved and refined to a final R of 0.068 for 1838 independent reflexions with I greater than 2 omega (I). The peptide backbone is folded at the Leu residue (phi L = -67.7, psi L = -49.1 degrees) without intramolecular hydrogen bonds. Considering each peptide plane, the Leu side-chain is oriented on the same side of that of the Phe residue and on the opposite side of that of the Met residue, respectively. The crystal conformation differs from all the other conformations proposed for FMLP-OMe and the anionic form of N-formyl-L-Met-L-Leu-L-Phe-OH (FMLP) in solution accounts for the amphiphilic character of the peptide, giving rise, through intermolecular hydrogen bonds, to a stacking of molecules which could be maintained in the aggregation states experimentally observed in solvents of low polarity. Intramolecular potential energy calculations have been carried out in order to compare the energies of the various backbone conformers.     

Conformational heterogeneity of quinoxaline peptides in solution.

The conformational heterogeneity of several quinoxaline antibiotics, a class of naturally occurring quinoxaline peptides with antitumor properties, and their synthetic analogues was investigated in polar and nonpolar solvents by high performance liquid chromatography (HPLC) with uv photodiode array detection, uv-absorbance, low-temperature phosphorescence, and nmr techniques. Multiple peak formation and interconversion in the HPLC and 1H-nmr analysis of triostin A, its under-N-methylated synthetic analogues (des-N-tetramethyltriostin A [TANDEM] and [N-MeCys3, N-MeCys7]-TANDEM [MCTAN-DEM]), and echinomycin were examined as a function of temperature, solvent polarity, and residence time in solution prior to analysis. Slow interconversion between HPLC peaks, ascribed to the presence of multiple solution conformers, was exhibited by these peptides although at very different interconversion rates. Among the triostins, the rate of interconversion appeared to vary with the degree of N-methylation of the residues in the cyclic depsipeptide chain. Interconversion of the n and p conformers of triostin A in chloroform occurred on a chromatographic timescale (a few minutes with kn----p calculated to be 0.02 s-1 at 25 degrees C) while the solution conformers of TANDEM in methanol equilibrated very slowly to one preferred conformer over a period of several weeks at ambient temperature. MCTANDEM, a synthetic analogue of triostin A with an intermediate degree of N-methylation of the residues in the peptide ring, consisted of an equilibrium mixture of n and p conformers in methanol that interconverted on a chromatographic time scale. Two additional conformers of MCTANDEM developed within a few weeks' residence time in methanol at ambient temperature. Echinomycin was found to exist in methanol as an interconverting mixture of at least four minor conformers in addition to the major isoform (95% by peak area) of the peptide. The solution conformers of the quinoxaline peptides investigated in this report are most likely a consequence of hindered rotation about the N-methylated peptide bonds in the depsipeptide ring and/or intramolecular hydrogen bonding.