Search

Folylpolyglutamate synthetase (FPGS) catalyzes the sequential addition of several glutamates to folate, forming gamma-linked polyglutamate folates of varying lengths. To understand how this protein is capable of accommodating ligands of different length and net charge, we have performed docking studies for folate substrates and glutamate based on the ternary crystal structure of Lactobacillus casei FPGS. Our results suggest two locations for folate binding, the one seen in the crystal structure and another distinct cavity. According to our model and experimental data, it is highly probable that folate can bind in both sites, and we suggest that the new pocket is especially important for the initial addition of the first glutamate residue. Docking longer substrates, di- and triglutamylated folates, showed how these molecules bind in the same sites. The longer folates also adopted transition-state-like conformations that may help us to understand the ligation reaction in FPGS and influence the design of mechanism-based inhibitors for anticancer or antimicrobial therapy.

The epoxidation of three stereolabeled methyl-substituted chiral allylic alcohols with (1,2)A and/or (1,3)A allylic strain, namely 3-methylbut-3-en-2-ol (1a), pent-3-en-2-ol (1b), and 3-methylpent-3-en-2-ol (1c), have been studied by the density-functional theory method, B3LYP/6-31+G(d,p). For each substrate we calculated the two prereaction complexes with Ti(OH)(4)/MeOOH (the oxidant model for Ti(O-i-Pr)(4)/t-BuOOH), their threo and erythro transition states for oxygen transfer, and the corresponding product complexes. For substrate 1a, the erythro transition state is 0.91 kcal/mol of lower energy than the threo one; for substrates 1b and 1c, the threo compared to the erythro transition states are by 1.05 and 0.21 kcal/mol more favorable, respectively. The threo/erythro product ratios have been estimated from the computed free energies for the competing threo and erythro transition states 3a-c in CH(2)Cl(2) solution to be 12:88 (1a), 92:8 (1b), and 77:23 (1c), which are in good accordance with the experimental values 22:78 (1a), 91:9 (1b), and 83:17 (1c). The diastereoselectivity of this diastereoselective oxyfunctionalization is rationalized in terms of the competition between (1,3)A and (1,2)A strain and the electronic advantage for the spiro transition state. In addition, solvent effects are also play a role for the diastereoselectivity at the same time.

The theoretical investigation of tetramethylammonium (TMA)−imidazole and TMA−furan complexes has been performed to justify the preferred structure of the cation−aromatic complexes predicated by the...

A detailed theoretical investigation of the tetramethylammonium(TMA)-benzene and TMA-pyrrole complexes has been performed to obtain the interaction properties of TMA with aromatics. Diffuse functions have been found to be important in the computational studies of these noncovalent complexes. Adding diffuse functions to the basis set decreases the binding energy by about 10% for the TMA-aromatic systems. Dispersion interactions in the TMA-aromatic systems are very important. They enhance the binding interactions between the TMA and the aromatic ring systems by about 0.5 kcal·mol-1 per interacting atomic pair, which is in agreement with the estimates of Rappe and Bernstein.1 Also, for the TMA-pyrrole complex, the presence of the dispersion interaction leads to a dramatic change in the optimized structure. Because B3LYP cannot handle properly the dispersion in the calculation, use of the Moller−Plesset second-order perturbation or other sophisticated methods should be considered in computational studies of c...

DFT/B3LYP calculations were carried out on complexes formed by NH4+ with aromatics, viz. benzene, phenol, pyrrole, imidazole, pyridine, indole, furane, and thiophene, to characterize the forces involved in such interactions and to gain further insight into the nature and diversity of cation−aromatic interactions. Such calculations may provide valuable information for understanding molecular recognition in biological systems and for force-field development. B3LYP/6-31G** optimization on 35 initial structures resulted in 11 different finally optimized geometries, which could be divided into three types: NH4+−π complexes, protonated heterocyclic−NH3 hydrogen bond complexes, and heterocyclic−NH4+ hydrogen bond complexes. For NH4+−π complexes, NH4+ always tilts toward the carbon−carbon bond rather than toward the heteroatom or the carbon−heteroatom bond. The calculated CHelpG charges suggest that the charge distribution of a free heterocyclic may be used to predict the geometry of its complex. Charge populati...