Induction of adiponectin by natural and synthetic phenolamides in mouse and human preadipocytes and its enhancement by docosahexaenoic acid.

Adiponectin, the adipose-derived cytokine, plays an important role in preventing metabolic syndromes. To develop new adiponectin inducers, eight species of ferulic esters and amides, and five related compounds were synthesized and tested on the stimulation of adiponectin production in mouse 3T3-L1 and normal human preadipocytes. The ferulamides with an aromatic ring in the N-substituent are very active in inducing adiponectin as compared with the known active compounds, curcumin, [6]-gingerol, and capsaicin, and furthermore the activities of these ferulamides are remarkably stronger than those of the corresponding esters or the straight chain octylamide. The most active compound, N-(2-phenylethyl)ferulamide (7), was found to activate the PPAR (peroxisome proliferator-activated receptor) gamma-RXR (retinoid X receptor) alpha heterodimeric complex in the PPRE (PPAR-responsive element)-driven luciferase reporter assay. The adiponectin production by 7 is synergistically enhanced by coaddition of a PPARgamma-specific agonist, pioglitazone (PGZ), or another PPARgamma agonist, docosahexaenoic acid (DHA), in cultured preadipocytes. The compound 7 alone did not show a statistically significant effect on the plasma adiponectin level in KK-A(y)/Ta mice, while 1% 7 in the diets significantly lowered the blood glucose and triglyceride levels and 0.3% 7 mixed with DHA oil in the diets significantly increased the adiponectin level as compared with the control. These results suggest that the present ferulamides would be useful lead compounds in developing more potent agents for treatment of metabolic syndromes through promoting the endogenous adiponectin production, and that such an activity is possibly enhanced by the coadministration with DHA.

The fragment molecular orbital method for geometry optimizations of polypeptides and proteins.

The fragment molecular orbital method (FMO) has been used with a large number of wave functions for single-point calculations, and its high accuracy in comparison to ab initio methods has been well established. We have developed the analytic derivative of the electrostatic interaction between far separated fragments and performed a number of restricted Hartree-Fock (RHF) geometry optimizations using FMO and ab initio methods. In particular, the alpha-helix, beta-turn, and extended conformers of a 10-residue polyalanine were studied and the good FMO accuracy was established (the rms deviations for the former two forms were about 0.2 A and for the latter structure about 0.001 A). Met-enkephalin dimer was used as a model for the polypeptide binding and computed at the 3-21G and 6-31G* levels with a similar accuracy achieved; the error in the binding energy predictions (FMO vs ab initio) was 1-3 kcal/mol. Chignolin (PDB: 1uao) and an agonist polypeptide of the erythropoietin receptor protein (emp1) were optimized at the 3-21(+)G level, with the rms deviation from ab initio of about 0.2 A, or 0.5 degrees in terms of bond angles. The effect of solvation on the structure optimization was studied in chignolin and the Trp-cage miniprotein construct (PDB:1l2y), by describing water with TIP3P. The computed structures in gas phase and solution are compared to each other and experiment.

Ab initio fragment molecular orbital (FMO) method applied to analysis of the ligand-protein interaction in a pheromone-binding protein.

Full quantum computation of the electronic state of proteins has recently become possible by the advent of the ab initio fragment molecular orbital (FMO) method. We applied this method to the analysis of the interaction between the Bombyx mori pheromone-binding protein and its ligand, bombykol. The protein-ligand interaction of this molecular complex was minutely analyzed by the FMO method, and the analysis revealed several important interactions between the ligand and amino acid residues.

Ab initio quantum mechanical study of the binding energies of human estrogen receptor alpha with its ligands: an application of fragment molecular orbital method.

We have theoretically examined the relative binding affinities (RBA) of typical ligands, 17beta-estradiol (EST), 17alpha-estradiol (ESTA), genistein (GEN), raloxifene (RAL), 4-hydroxytamoxifen (OHT), tamoxifen (TAM), clomifene (CLO), 4-hydroxyclomifene (OHC), diethylstilbestrol (DES), bisphenol A (BISA), and bisphenol F (BISF), to the alpha-subtype of the human estrogen receptor ligand-binding domain (hERalpha LBD), by calculating their binding energies. The ab initio fragment molecular orbital (FMO) method, which we have recently proposed for the calculations of macromolecules such as proteins, was applied at the HF/STO-3G level. The receptor protein was primarily modeled by 50 amino acid residues surrounding the ligand. The number of atoms in these model complexes is about 850, including hydrogen atoms. For the complexes with EST, RAL, OHT, and DES, the binding energies were calculated again with the entire ERalphaLBD consisting of 241 residues or about 4000 atoms. No significant difference was found in the calculated binding energies between the model and the real protein complexes. This indicates that the binding between the protein and its ligands is well characterized by the model protein with the 50 residues. The calculated binding energies relative to EST were very well correlated with the experimental RBA (the correlation coefficient r=0.837) for the ligands studied in this work. We also found that the charge transfer between ER and ligands is significant on ER-ligand binding. To our knowledge, this is the first achievement of ab initio quantum mechanical calculations of large molecules such as the entire ERalphaLBD protein.

Molecular dynamics simulations revealed Ca(2+)-dependent conformational change of Calmodulin.

Molecular dynamics simulations were performed to simulate Ca(2+)-dependent conformational change of calmodulin (CaM). Simulations of the fully Ca(2+)-bound form of CaM (Holo-CaM) and the Ca(2+)-free form (Apo-CaM) were performed in solution for 4 ns starting from the X-ray crystal structure of Holo-CaM. A striking difference was observed between the trajectories of Holo-CaM and Apo-CaM: the central helix remained straight in the former but became largely bent in the latter. Also, the flexibility of Apo-CaM was higher than that of Holo-CaM. The results indicated that the bound Ca(2+) ions harden the structure of CaM.

Electrostatic Interactions That Determine the Rate of Pseudorotation Processes in Oxyphosphorane Intermediates: Implications with Respect to the Roles of Metal Ions in the Enzymatic Cleavage of RNA.

The enzymatic cleavage of RNA takes place via a cyclic pentacoordinate oxyphosphorane intermediate/transition state. We carried out ab initio investigations on the neutral cyclic oxyphosphorane, which exists as a stable intermediate. As a consequence of the conformational preferences of the pentacoordinate trigonal bipyramidal intermediates, the rotation of the P-OH bonds is strongly coupled with the reaction coordinate for the pseudorotation process. In addition, the neutral PF(4)OH species has a higher barrier to pseudorotation than the corresponding anionic species PF(4)O(-). These findings are related to the positive charge of the hydrogen atoms on the equatorial oxygens in the trigonal bipyramidal structures: the hydrogen atoms preferably adopt eclipsed positions relative to the axial ligands. Fixing the cationic species in these regions causes an increase in the barrier heights for pseudorotation processes and, thus, prevents isomerization by pseudorotation. Consequently, metal coordination in the double-metal ion mechanism for enzymatic cleavage of RNA should serve to exclusively stabilize the trigonal bipyramidal intermediate/transition state for the in-line attack and departure process.