On the nature of the transition state in catechol O-methyltransferase. A complementary study based on molecular dynamics and potential energy surface explorations

The way in which enzymes influence the rate of chemical processes is still a question of debate. The protein promotes the catalysis of biochemical processes by lowering the free energy barrier in comparison with the reference uncatalyzed reaction in solution. In this article we are reporting static and dynamic aspects of the enzyme catalysis in a bimolecular reaction, namely a methyl transfer from S-adenosylmethionine to the hydroxylate oxygen of a substituted catechol catalyzed by catechol O-methyltransferase. From QM/MM optimizations, we will first analyze the participation of the environment on the transition vector. The study of molecular dynamics trajectories will allow us to estimate the transmission coefficient from a previously localized transition state as the maximum in the potential of mean force profile. The analysis of the reactive and nonreactive trajectories in the enzyme environment and in solution will also allow studying the geometrical and electronic changes, with special attention to the chemical system movements and the coupling with the environment. The main result, coming from both analyses, is the approximation of the magnesium cation to the nucleophilic and the hydroxyl group of the catecholate as a result of a general movement of the protein, stabilizing in this way the transition state. Consequently, the free energy barrier of the enzyme reaction is dramatically decreased with respect to the reaction in solution.     

Computationally-assisted approach to the vibrational spectra of molecular crystals: study of hydrogen-bonding and pseudo-polymorphism.

A new computationally-assisted methodology (PiMM), which accounts for the effects of intermolecular interactions in the crystal, is applied to the complete assignment of the Raman and infrared vibrational spectra of room temperature forms of crystalline caffeine, theobromine, and theophylline. The vibrational shifts due to crystal packing interactions are evaluated from ab initio calculations for a set of suitable molecular pairs, using the B3LYP/6-31G* approach.The proposed methodology provides an answer to the current demand for a reliable assignment of the vibrational spectra of these methyl-xanthines, and clarifies several misleading assignments. The most relevant intermolecular interactions in each system and their effect on the vibrational spectra are considered and discussed. Based on these results, significant insights are obtained for the structure of caffeine in the anhydrous form (stable at room temperature), for which no X-ray structure has been reported. A possible structure based on C((8))--H...N((9)) and C((1,3))--H...O intermolecular interactions is suggested.     

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.     

Characterization of sodiated glycerol phosphatidylcholine phospholipids by mass spectrometry

Fast atom bombardment mass spectrometry in the positive mode was used for the characterization of sodiated glycerol phosphatidylcholines. The relative abundance (RA) of the protonated species is similar to the RA of the sodiated molecular species. The sodiated fragment ion, [M + Na - 59](+), corresponding to the loss of trimethylamine, and other sodiated fragment ions, were also observed. The decomposition of the sodiated molecule is very similar for all the studied glycerol phosphatidylcholines, in which the most abundant ion corresponds to a neutral loss of 59 Da. Upon collision-induced dissociation (CID) of the [M + Na](+) ion informative ions are formed by the losses of the fatty acids in the sn-1 and sn-2 positions. Other major fragment ions of the sodiated molecule result from loss of non-sodiated and sodiated choline phosphate, [M + Na - 183](+), [M + Na - 184](+.) and [M + Na - 205](+), respectively. The main CID fragmentation pathway of the [M + Na - 59](+) ion yields the [M + Na - 183](+) ion, also observed in the CID spectra of the [M + Na](+) molecular ion. Other major fragment ions are [M + Na - 205](+) and the fragment ion at m/z 147. Collisional activation of [M + Na - 205](+) results in charge site remote fragmentation of both fatty acid alkyl chains. The terminal ions of these series of charge remote fragmentations result from loss of part of the R(1) or R(2) alkyl chain. Other major informative ions correspond to acylium ions.     

Constant neutral loss scanning for the characterization of glycerol phosphatidylcholine phospholipids

The product-ion spectra of the sodiated molecules of glycerol phosphatidylcholine phospholipids (GPC) were obtained, using fast-atom bombardment (FAB) as the ionization method, in a tandem mass spectrometer. The product-ion spectra of these sodiated molecules of GPCs were found to differ significantly from those of the protonated GPC molecules. This difference is due to the absence of the ion of m/z 184 (protonated-phosphocholine moiety) and to the presence of an ion resulting from the loss of trimethylamine (m=59 Da) from the polar head group, which is the dominant fragmentation. This characteristic neutral loss provides a means of identification of this class of phospholipids and of differentiation from other phospholipid classes in complex mixtures by performing a constant-neutral-loss scan of 59.     

Chemical composition,antioxidant and cytotoxic activities of extracts from the leaves of Smilax brasiliensis Sprengel (Smilacaceae)

The antioxidant and cytotoxic activities of petroleum ether and methanol extracts, fatty acids and methyl esters from leaves of Smilax brasiliensis were evaluated, and the composition of the extracts was determined. Palmitic, linoleic and linolenic acids were major components of the extracts. For antioxidant activity, all samples exhibited IC₅₀ values lower than BHT (2,6-di-tert-butyl-4-methylphenol). The extracts, fatty acids and methyl esters from S. brasiliensis presented no toxicity to larvae of the brine shrimp, Artemia salina. Among the purified substances, only methyl linolenate showed toxicity (LD₅₀ = 21.47 μg/mL). This study showed, for the first time, the composition of petroleum ether and methanol extracts from S. brasiliensis leaves, as well as the antioxidant and cytotoxic activities of extracts, fatty acids and methyl esters.