Dataset Information

Molecular dynamics simulation of tryptophan hydroxylase-1: binding modes and free energy analysis to phenylalanine derivative inhibitors.

ABSTRACT: Serotonin is a neurotransmitter that modulates many central and peripheral functions. Tryptophan hydroxylase-1 (TPH1) is a key enzyme of serotonin synthesis. In the current study, the interaction mechanism of phenylalanine derivative TPH1 inhibitors was investigated using molecular dynamics (MD) simulations, free energy calculations, free energy decomposition analysis and computational alanine scanning. The predicted binding free energies of these complexes are consistent with the experimental data. The analysis of the individual energy terms indicates that although the van der Waals and electrostatics interaction contributions are important in distinguishing the binding affinities of these inhibitors, the electrostatic contribution plays a more crucial role in that. Moreover, it is observed that different configurations of the naphthalene substituent could form different binding patterns with protein, yet lead to similar inhibitory potency. The combination of different molecular modeling techniques is an efficient way to interpret the interaction mechanism of inhibitors and our work could provide valuable information for the TPH1 inhibitor design in the future.

SUBMITTER: Zhong H

PROVIDER: S-EPMC3676822 | biostudies-literature | 2013

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Molecular dynamics simulation of tryptophan hydroxylase-1: binding modes and free energy analysis to phenylalanine derivative inhibitors.

Zhong Hao H Huang Wei W He Gu G Peng Cheng C Wu Fengbo F Ouyang Liang L

International journal of molecular sciences 20130510 5

Serotonin is a neurotransmitter that modulates many central and peripheral functions. Tryptophan hydroxylase-1 (TPH1) is a key enzyme of serotonin synthesis. In the current study, the interaction mechanism of phenylalanine derivative TPH1 inhibitors was investigated using molecular dynamics (MD) simulations, free energy calculations, free energy decomposition analysis and computational alanine scanning. The predicted binding free energies of these complexes are consistent with the experimental d ...[more]

PMID: 23665899

Similar Datasets

Project description:RNA hairpin loops are one of the most abundant secondary structure elements and participate in RNA folding and protein-RNA recognition. To characterize the free energy surface of RNA hairpin folding at an atomic level, we calculated the potential of mean force (PMF) as a function of the end-to-end distance, by using umbrella sampling simulations in explicit solvent. Two RNA hairpins containing tetraloop cUUCGg and cUUUUg are studied with AMBER ff99 and CHARMM27 force fields. Experimentally, the UUCG hairpin is known to be significantly more stable than UUUU. In this study, the calculations using AMBER force field give a qualitatively correct description for the folding of two RNA hairpins, as the calculated PMF confirms the global stability of the folded structures and the resulting relative folding free energy is in quantitative agreement with the experimental result. The hairpin stabilities are also correctly differentiated by the more rapid molecular mechanics-Poisson Boltzmann-surface area approach, but the relative free energy estimated from this method is overestimated. The free energy profile shows that the native state basin and the unfolded state plateau are separated by a wide shoulder region, which samples a variety of native-like structures with frayed terminal basepair. The calculated PMF lacks major barriers that are expected near the transition regions, and this is attributed to the limitation of the 1-D reaction coordinate. The PMF results are compared with other studies of small RNA hairpins using kinetics method and coarse grained models. The two RNA hairpins described by CHARMM27 are significantly more deformable than those represented by AMBER. Compared with the AMBER results, the CHARMM27 calculated DeltaG(fold) for the UUUU tetraloop is in better agreement with the experimental results. However, the CHARMM27 calculation does not confirm the global stability of the experimental UUCG structure; instead, the extended conformations are predicted to be thermodynamically stable in solution. This finding is further supported by separate unrestrained CHARMM27 simulations, in which the UUCG hairpin unfolds spontaneously within 10 ns. The instability of the UUCG hairpin originates from the loop region, and propagates to the stem. The results of this study provide a molecular picture of RNA hairpin unfolding and reveal problems in the force field descriptions for the conformational energy of certain RNA hairpin.

Project description:Phenylalanine hydroxylase (PheH) and tryptophan hydroxylase (TrpH) catalyze the aromatic hydroxylation of phenylalanine and tryptophan, forming tyrosine and 5-hydroxytryptophan, respectively. The reactions of PheH and TrpH have been investigated with [4-(2)H]-, [3,5-(2)H(2)]-, and (2)H(5)-phenylalanine as substrates. All (D)k(cat) values are normal with Delta117PheH, the catalytic core of rat phenylalanine hydroxylase, ranging from 1.12-1.41. In contrast, for Delta117PheH V379D, a mutant protein in which the stoichiometry between tetrahydropterin oxidation and amino acid hydroxylation is altered, the (D)k(cat) value with [4-(2)H]-phenylalanine is 0.92 but is normal with [3,5-(2)H(2)]-phenylalanine. The ratio of tetrahydropterin oxidation to amino acid hydroxylation for Delta117PheH V379D shows a similar inverse isotope effect with [4-(2)H]-phenylalanine. Intramolecular isotope effects, determined from the deuterium contents of the tyrosine formed from [4-(2)H]-and [3,5(2)H(2)]-phenylalanine, are identical for Delta117PheH and Delta117PheH V379D, suggesting that steps subsequent to oxygen addition are unaffected in the mutant protein. The inverse effects are consistent with the reaction of an activated ferryl-oxo species at the para position of the side chain of the amino acid to form a cationic intermediate. The normal effects on the (D)k(cat) value for the wild-type enzyme are attributed to an isotope effect of 5.1 on the tautomerization of a dienone intermediate to tyrosine with a rate constant 6- to7-fold that for hydroxylation. In addition, there is a slight ( approximately 34%) preference for the loss of the hydrogen originally at C4 of phenylalanine. With (2)H(5)-indole-tryptophan as a substrate for Delta117PheH, the (D)k(cat) value is 0.89, consistent with hydroxylation being rate-limiting in this case. When deuterated phenylalanines are used as substrates for TrpH, the (D)k(cat) values are within error of those for Delta117PheH V379D. Overall, these results are consistent with the aromatic amino acid hydroxylases all sharing the same chemical mechanism, but with the isotope effect for hydroxylation by PheH being masked by tautomerization of an enedione intermediate to tyrosine.

Dataset Information

Molecular dynamics simulation of tryptophan hydroxylase-1: binding modes and free energy analysis to phenylalanine derivative inhibitors.

Publications

Molecular dynamics simulation of tryptophan hydroxylase-1: binding modes and free energy analysis to phenylalanine derivative inhibitors.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure