Theoretical modeling of enzyme catalytic power: analysis of "cratic" and electrostatic factors in catechol O-methyltransferase

A comparative theoretical study of a bimolecular reaction in aqueous solution and catalyzed by the enzyme catechol O-methyltransferase (COMT) has been carried out by a combination of two hybrid QM/MM techniques: statistical simulation methods and internal energy minimizations. In contrast to previous studies by other workers, we have located and characterized transition structures for the reaction in the enzyme active site, in water and in a vacuum, and our potential of mean force calculations are based upon reaction coordinates obtained from features of the potential energy surfaces in the condensed media, not from the gas phase. The AM1/CHARMM calculated free energy of activation for the reaction of S-adenosyl methionine (SAM) with catecholate catalyzed by COMT is 15 kcal mol(-1) lower the AM1/TIP3P free-energy barrier for the reaction of the trimethylsulfonium cation with the catecholate anion in water at 300 K, in agreement with previous estimates. The thermodynamically preferred form of the reactants in the uncatalyzed model reaction in water is a solvent-separated ion pair (SSIP). Conversion of the SSIP into a contact ion pair, with a structure resembling that of the Michaelis complex (MC) for the reaction in the COMT active site, is unfavorable by 7 kcal mol(-1), largely due to reorganization of the solvent. We have considered alternative ways to estimate the so-called "cratic" free energy for bringing the reactant species together in the correct orientation for reaction but conclude that direct evaluation of the free energy of association by means of molecular dynamics simulation with a simple standard-state correction is probably the best approach. The latter correction allows for the fact that the size of the unit cell employed with the periodic boundary simulations does not correspond to the standard state concentration of 1 M. Consideration of MC-like species allows a helpful decomposition of the catalytic effect into preorganization and reorganization phases. In the preorganization phase, the substrates are brought together into the MC-like species, either in water or in the enzyme active site. In the reorganization phase, the roles of the enzymic and aqueous environments may be compared directly because reorganization of the substrate is about the same in both cases. Analysis of the electric field along the reaction coordinate demonstrates that in water the TS is destabilized with respect to the MC-like species because the polarity of the solute diminishes and consequently the reaction field is also decreased. In the enzyme, the electric field is mainly a permanent field and consequently there is only a small reorganization of the environment. Therefore, destabilization of the TS is lower than in solution, and the activation barrier is smaller.     

Diheteroarylmethanes. 6.(1) Probing the Electron-Withdrawing Rank of Heteroaryl Groups by Conformational Studies of Push-Pull Ethylenes. Isolation of NH-Enamine Tautomers of alpha,alpha-Diheteroarylacetaldehydes (Heteroaryl = 2-Benzoxazolyl, 2-Benzothiazolyl)

Enamines and enol ethers substituted in the beta position by the 2-benzoxazolyl and 2-benzothiazolyl group have been obtained by condensing dimethylformamide dimethyl acetal and ethyl orthoformate, respectively, with bis(2-benzoxazolyl)methane and bis(2-benzothiazolyl)methane. A dynamic NMR and semiempirical (PM3) investigation on rotational energy barriers has been carried out in order to rank the electron-withdrawing capacity of the heterocyclic rings. The NMR-based evaluation of the energy barriers corresponding to the rotation along the enaminic double bond has shown that the pi-electron-withdrawing power of benzothiazole is larger than that of benzoxazole, in full accord with previously obtained charge demand values based on (13)C and (15)N pi-charge/shift relationships. The NMR and the computational approaches have led to consistent results. The X-ray crystal structure of alpha,alpha-bis(2-benzothiazolyl)-beta-(dimethylamino)ethene shows that only one heteroaryl ring is coplanar with the enaminic double bond, while the second one is twisted by 73 degrees relative to such a plane. Moreover, in this case calculations closely reproduce the experimental results. In the calculated transition states corresponding to the rotational process along the enamine double bond, the two heteroaryl groups become coplanar and conjugation develops between the enaminic nitrogen electron pool and both heterocycles. The lower rotational barrier observed in the case of the 2-benzothiazolyl derivatives, with respect to the 2-benzoxazolyl derivatives, is therefore a direct consequence of the higher electron-withdrawing power of the former group. Furthermore, a stabilizing intramolecular nonbonded S.S interaction in the rotational transition state of the benzothiazolyl derivatives is evidenced by the calculations. An unprecedented isolation of the NH-enamine tautomer involving the benzoxazolyl and benzothiazolyl ring in alpha,alpha-diheteroacetaldehydes has been performed, confirming these heterocycles as strong electron-withdrawing substituents.     

High-frequency and -field electron paramagnetic resonance of high-spin manganese(III) in tetrapyrrole complexes

High-field and -frequency electron paramagnetic resonance (HFEPR) spectroscopy has been used to study three complexes of high spin Manganese(III), 3d4, S = 2. The complexes studied were tetraphenylporphyrinatomanganese(III) chloride (MnTPPCI), phthalocyanatomanganese(III) chloride (MnPcCl), and (8,12-diethyl-2,3,7,13,17,18-hexamethylcorrolato)manganese(III) (MnCor). We demonstrate the ability to obtain both field-oriented (single-crystal like) spectra and true powder pattern HFEPR spectra of solid samples. The latter are obtained by immobilizing the powder, either in an n-eicosane mull or KBr pellet. We can also obtain frozen solution HFEPR spectra with good signal-to-noise, and yielding the expected true powder pattern. Frozen solution spectra are described for MnTPPCl in 2:3 (v/v) toluene/CH2Cl2 solution and for MnCor in neat pyridine (py) solution. All of the HFEPR spectra have been fully analyzed using spectral simulation software and a complete set of spin Hamiltonian parameters has been determined for each complex in each medium. Both porphyrinic complexes (MnTPPCl and MnPcCl) are rigorously axial systems, with similar axial zero-field splitting (zfs): D approximately -2.3 cm(-1), and g values quite close to 2.00. In contrast, the corrole complex, MnCor, exhibits slightly larger magnitude, rhombic zfs: D approximtely -2.6 cm(-1), absolute value(E) approximately 0.015 cm(-1), also with g values quite close to 2.00. These results are discussed in terms of the molecular structures of these complexes and their electronic structure. We propose that there is a significant mixing of the triplet (S = 1) excited state with the quintet (S= 2) ground state in Mn(III) complexes with porphyrinic ligands, which is even more pronounced for corroles.     

Variants of Solid-State and Solution Structures of (η3-Allyl)- {2-[2′-(diphenylphosphino)phenyl]-4,5-dihydrooxazole-P,N}palladium(II) hexafluorophosphates and tetraphenylborates

Crystal and solution structures of the enantiomerically pure and the racemic pairs of (η³-allyl) {2-[2′-(diphenylphosphino)phenyl]-4,5-dihydro-4-phenyloxazole}palladium(II) hexafluorophosphates ( 1 , and rac- 1 , resp.) and tetraphenylborates ( 2 , and rac- 2 , resp.) as well as (η³-allyl){2-[2′-(diphenylphosphino)phenyl]-4,5-dihydro-4-isopropyloxazole}palladium(II) tetraphenylborate ( 3 ) were characterized by X-ray crystallography and ¹H-NMR spectroscopy. In the solid state, rac- 1 and rac- 2 proved to be disordered with both diastereoisomeric complexes in the crystal. The complexes 2 and 3 exist only in the ‘exo’ form. The X-ray structures show that the [Pd^II(η³-allyl)] moiety may adopt different configurations between a nearly symmetrical three-electron Pd^II(π³-allyl) system and an asymmetrical allyl group with a η¹- and a η²-bonding to the metal center. The [Pd^II(η³-allyl)] system of rac- 1 and of ‘endo’ rac- 2 is closer to the former, and that of 2 , ‘exo’-rac- 2 , and 3 closer to the later geometry. The ¹H-NMR spectra of the hexafluorophosphates 1 and rac- 1 show two sets of signals of the allylic protons in an ‘exo’/‘endo’ ratio of 2:3. The tetraphenylborates 2, rac- 2 , and 3 give only one set of broad signals of the allylic protons.     

Design, synthesis, and photophysical studies of a porphyrin-fullerene dyad with parachute topology; charge recombination in the marcus inverted region

As part of a continuing investigation of the topological control of intramolecular electron transfer (ET) in donor-acceptor systems, a symmetrical parachute-shaped octaethylporphyrin-fullerene dyad has been synthesized. A symmetrical strap, attached to ortho positions of phenyl groups at opposing meso positions of the porphyrin, was linked to [60]-fullerene in the final step of the synthesis. The dyad structures were confirmed by (1)H, (13)C, and (3)He NMR, and MALDI-TOF mass spectra. The free-base and Zn-containing dyads were subjected to extensive spectroscopic, electrochemical and photophysical studies. UV-vis spectra of the dyads are superimposable on the sum of the spectra of appropriate model systems, indicating that there is no significant ground-state electronic interaction between the component chromophores. Molecular modeling studies reveal that the lowest energy conformation of the dyad is not the C(2)(v)() symmetrical structure, but rather one in which the porphyrin moves over to the side of the fullerene sphere, bringing the two pi-systems into close proximity, which enhances van der Waals attractive forces. To account for the NMR data, it is proposed that the dyad is conformationally mobile at room temperature, with the porphyrin swinging back and forth from one side of the fullerene to the other. The extensive fluorescence quenching in both the free base and Zn dyads is associated with an extremely rapid photoinduced electron-transfer process, k(ET) approximately 10(11) s(-)(1), generating porphyrin radical cations and C(60) radical anions, detected by transient absorption spectroscopy. Back electron transfer (BET) is slower than charge separation by up to 2 orders of magnitude in these systems. The BET rate is slower in nonpolar than in polar solvents, indicating that BET occurs in the Marcus inverted region, where the rate decreases as the thermodynamic driving force for BET increases. Transient absorption and singlet molecular oxygen sensitization data show that fullerene triplets are formed only with the free base dyad in toluene, where triplet formation from the charge-separated state is competitive with decay to the ground state. The photophysical properties of the P-C(60) dyads with parachute topology are very similar to those of structurally related rigid pi-stacked P-C(60) dyads, with the exception that there is no detectable charge-transfer absorption in the parachute systems, attributed to their conformational flexibility. It is concluded that charge separation in these hybrid systems occurs through space in unsymmetrical conformations, where the center-to-center distance between the component pi-systems is minimized. Analysis of the BET data using Marcus theory gives reorganization energies for these systems between 0.6 and 0.8 eV and electronic coupling matrix elements between 4.8 and 5.6 cm(-)(1).     

(Perhalogenmethylthio)heterocyclen. X. Säurekatalysierte Substitutionen an (Perchlorfluormethylthio)pyrrolen und deren agrobiologische Wirkung

(Perhalomethylthio)heterocycles. X

1 IX. Mitt.: s. [1].

. Acid-catalyzed substitutions on (perchlorofluoromethylthio)pyrroles and their agro-biological activities In the presence of C₄F₉SO₃H the (perhalomethylthio)pyrroles 1a–c react with Cl_3?nF_nCSCl (n = 1–3) to give mixtures of the 2,5- and 2,4-disubstituted pyrroles 2a–f and 3a–h . 2a and 3a react with CF₃SCl (catalyst CF₃SO₃H) yielding 2,3,5-tris (trifloromethylthio)pyrrole ( 4a ), which under similar conditions reacts further to give 2,3,4,5-tetrakis (trifluoromethylthio)pyrrole ( 5 ). As a by-product during the conversion of 3a to 4a 2,3,4-tris (trifluoromethylthio)pyrrole ( 4b ) is formed. The pyrroles 2a , 4a and 5 form the mercury salts 6a–c ; compound 5 yields also a silver salt 7 . The ¹H- and ¹⁹F-NMR. spectra are discussed and the agro-biological properties of the compounds investigated.     

Designing generalized statistical ensembles for numerical simulations of biopolymers

Conformational properties of polymers, such as average dihedral angles or molecular alpha-helicity, display a rather weak dependence on the detailed arrangement of the elementary constituents (atoms). We propose a computer simulation method to explore the polymer phase space using a variant of the standard multicanonical method, in which the density of states associated to suitably chosen configurational variables is considered in place of the standard energy density of states. This configurational density of states is used in the Metropolis acceptance/rejection test when configurations are generated with the help of a hybrid Monte Carlo algorithm. The resulting configurational probability distribution is then modulated by exponential factors derived from the general principle of the maximal constrained entropy by requiring that certain average configurational quantities take preassigned (possibly temperature dependent) values. Thermal averages of other configurational quantities can be computed by using the probability distributions obtained in this way. Moments of the energy distribution require an extra canonical sampling of the system phase space at the desired temperature, in order to locally thermalize the configurational degrees of freedom. As an application of these ideas we present the study of the structural properties of two simple models: a bead-and-spring model of polyethylene with independent hindered torsions and an all-atom model of alanine and glycine oligomers with 12 amino acids in vacuum.     

A new veratramine alkaloid from the bulbs of Fritillaria hupehensis

A new highly conjugated alkaloid of veratramine type, 225,25S,5α-veratramine-7(8),12(14)-diene-3β,13β,23β-triol-6-one (1), was isolated from the bulbs of Fritillaria hupehensis Hsiao et K.C. Hsia. Its structure was determined on the basis of spectroscopic evidences.     

Synthesis of novel 1,4‐dihydropyrido[3′,2′:5,6]thiopyrano[4,3‐c]‐pyrazoles and 5H‐pyrido[3′,2′:5,6]thiopyrano[4,3‐d]pyrimidines as potential antiproliferative agents

The synthesis of new planar derivatives characterized by the presence of a pyridothiopyranopyrazole or pyridothiopyranopyrimidine nucleus, carrying a substituted aryl group, is reported. The novel 1,4‐dihydropyrido[3′,2′:5,6]thiopyrano[4,3‐c]pyrazole derivatives were obtained by condensation of 2,3‐dihydro‐3‐hydroxymethylenethiopyrano[2,3‐b]pyridin‐4(4H)‐ones with appropriate hydrazines. The preparation of 2‐substituted pyrido[3′,2′:5,6]thiopyrano[4,3‐d]pyrimidines was accomplished from the intermediate 2,3‐dihy‐dro‐3‐dimethylaminomethylenethiopyrano[2,3‐b]pyridin‐4(4H)‐ones by reaction with the appropriate binucleophile amidines. The antiproliferative activity of some new products was tested by an in vitro assay on human tumour cell lines (HL‐60 and HeLa), but none of them showed any significant effects in the tests performed. Accordingly, linear flow dichroism measurements indicated their inability to form a molecular complex with DNA.     

BINOL catalyzed enantioselective addition of titanium phenylacetylide to aromatic ketones

An enantioselective addition of titanium phenylacetylide to ketones, promoted by BINOL, is described; this new enantioselective protocol gives high enantiomeric excess (up to 90% ee) with aromatic ketones using a simple procedure without pyrophoric or expensive reagents.