PubMed Health⌕ Search

Biomedical subjects

N Yamaotsu

Publications and source records attributed to N Yamaotsu.

5 recordsLinked to original sources

Molecular dynamics simulation of the calmodulin-trifluoperazine complex in aqueous solution.

Trifluoperazine (TFP) has been widely studied in relation to its mode of binding and its inactivation of calmodulin (CaM). Most studies in solution have indicated that CaM has two high-affinity binding sites for TFP. The crystal structure of the 1:4 CaM-TFP complex (CaM-4TFP) shows that three TFP molecules bind to the C-domain of CaM, and that one TFP molecule binds to the N-domain. In contrast, the crystal structure of the 1:1 CaM-TFP complex (CaM-1TFP) shows that one TFP molecule binds to the C-domain. It has been thought that the binding of one TFP molecule to the C-domain is followed by binding to the N-domain. The crystal structure of the 1:2 CaM-TFP complex (CaM-2TFP), moreover, has recently been determined, showing that two TFP molecules bind to the C-domain. In order to determine the structure of the CaM-TFP complex and to clarify the interaction between CaM and TFP in solution, we performed a molecular dynamics simulation of the CaM-TFP complex in aqueous solution starting from the CaM-4TFP crystal structure. The obtained solution structure is very similar to the CaM-2TFP crystal structure. The computer simulation showed that the binding ability of the secondary binding site of the C-domain is higher than that of the primary binding site of the N-domain.

Binding Sites↗

Rational procedure for 3D-QSAR analysis using TRNOE experiments and computational methods: application to thermolysin inhibitors.

The success or failure of 3D QSAR, particularly CoMFA, is most strongly dependent, especially for flexible compounds, on the conformation of the molecule used in the analysis, and on the orientation of the molecule relative to the other molecules in 3D space (i.e., alignment). The present study suggests a rational procedure for the estimation of binding conformation that uses the transferred nuclear Overhauser effect (TRNOE) experiment in combination with conformational analysis using CAMDAS (Conformational Analyzer with Molecular Dynamics And Sampling) program that is developed in our laboratory. In the next step the TRNOE-obtained conformation can be used as a reference template in order to obtain alignment of other ligands, that have a common binding site. In this step we used the SUPERPOSE program created in our laboratory, in order to estimate the binding conformation of other compounds, and to simultaneously obtain the alignment of compounds for CoMFA. The resulting CoMFA models could be expected to closely reproduce the interaction mode with protein represented by the reported X-ray results. In order to confirm the validity of our procedure described above, we show its application in obtaining CoMFA models of thermolysin inhibitors. We obtained twenty CoMFA models, and that with the highest q2 value (q2 = 0.701) was found to provide an interaction mode very similar to that represented by the X-ray results.

Molecular Conformation↗

Analysis of affinities of penicillins for a class C beta-lactamase by molecular dynamics simulations.

We present a calculation for the binding free energy difference between two complexes of the class C beta-lactamase from Enterobacter cloacae with foramidocillin (FOPC) and with piperacillin (PIPC). The calculation was carried out by means of the thermodynamic integration (TI) method implemented with molecular dynamics (MD). By use of the available crystal structure of the class C beta-lactamase from E. cloacae, the structures of the beta-lactamase-FOPC (FOPC complex) and beta-lactamase-PIPC (PIPC complex) complexes were built by molecular modeling and equilibrated with MD simulations. FOPC were gradually converted into PIPC in both the solution and the enzyme system by means of MD/TI methods during the MD simulation. The structure of the PIPC complex as derived by the MD/TI simulation was similar to that of the PIPC complex obtained from molecular modeling. The calculated difference in the free energy of binding (deltadeltaGbind) was -2.2 kcal/mol. This compares well with the experimental value of -1.5 kcal/mol. The results indicate that the binding affinity of FOPC is lower than that of PIPC because of the greater difficulty of desolvation for FOPC upon binding to the enzyme. This calculation suggests that the desolvation of the ligand, as well as its interaction with the beta-lactamase, is important in understanding the relative affinity of the ligands with beta-lactamase.

Crystallography↗

Estimation of stabilities of staphylococcal nuclease mutants (Met32-->Ala and Met32-->Leu) using molecular dynamics/free energy perturbation.

We performed molecular dynamics (MD)/free energy perturbation (FEP) calculations to reproduce the experimental free energy difference of denaturation for staphylococcal nuclease mutant Met32-->Ala (M32A) and to predict the stability of the mutant Met32-->Leu (M32L). The calculated free energy difference of denaturation for the M32A of -1.9 kcal/mol was in good agreement with the experimental value. In the M32A, a small hydrophobic core formed by three aromatic rings (Tyr-27, Phe-34, Phe-76) in a wild-type crumbled as a result of exposure to water. The van der Waals interactions in the native state of the M32A were weaker than those of the wild-type, which strongly suggests that the Met-32 is important for the stability of the enzyme. The M32L has not been available yet, but is expected to retain the small hydrophobic core. The free energy difference of denaturation for the M32L was calculated to be 1.6 kcal/mol. The MD/FEP simulation showed that the native state structure of the M32L was only slightly changed when compared with that of the wild-type. It was suggested that the M32L is more stable than the wild-type because the electrostatic interactions in the denatured state are more disadvantageous than those in the native state.

Enzyme Stability↗

Molecular dynamics study of the stability of staphylococcal nuclease mutants: component analysis of the free energy difference of denaturation.

The stability of two mutants G88V (Gly-88-->Val) and A69T (Ala-69-->Thr) of staphylococcal nuclease was analyzed by molecular dynamics simulations. The calculated free energy differences of denaturation for G88V and A69T were -1.1 and -2.8 kcal/mol, respectively. These values are in good agreement with the experimental values. The free energy differences divided into electrostatic and van der Waals components were analyzed. These two mutants are mainly destabilized due to van der Waals interactions. There is little difference between the electrostatic contribution to the free energy change in the native state and that in the denatured state. In each mutant structure, a small cavity appears in the vicinity of the perturbed residue. It is suggested that intramolecular van der Waals interactions of the mutants are weaker than those of the wild-type. Furthermore, analyses of the contributions of each residue near the perturbed residue and of water to the free energy difference of denaturation suggest that the interaction between water and the perturbed residue plays a very important role in the stability of staphylococcal nuclease, and that a small hydrophobic core consisting of the three aromatic rings (Tyr-27, Phe-34, Phe-76) and the side chain of Met-32 is also important for the stability.

Alanine↗