• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.125-137
• /
• 2013

2012년 G 단백질 연결 수용체(G-Protein Coupled Receptors ; GPCR) 연구가 노벨 화학상을 받았다. 상당히 많은 병과 관련되어 있어 잠재력이 크고, 많은 연구가 진행 중이다. 현재 리간드와 단백질간의 정전기적 포텐셜 연구를 통한 예측 연구가 진행되고 있지만, GPCR과 리간드 간의 연구에서 아직 리간드의 전하를 통한 단백질과 리간드간의 상호작용 예측 연구가 되어 있지 않다. 그렇기 때문에 이번 연구에서는 8가지 방법으로 전하(charge)를 띠게 하여서 단백질과 리간드의 상호작용을 계산을 통하여 예측하여 보았다.

• Bulletin of the Korean Chemical Society
• /
• v.31 no.10
• /
• pp.2903-2908
• /
• 2010

Protein S-nitrosation is common in cells under nitrosative stress. In order to model proteins with S-nitrosocysteine (CysSNO) residues, we first developed an Amber force field for S-nitrosoethanethiol (EtSNO) and then transferred it to CysSNO. Partial atomic charges for EtSNO and CysSNO were obtained by a restrained electrostatic potential approach to be compatible with the Amber-99 force field. The force field parameters for bonds and angles in EtSNO were obtained from a generalized Amber force field (GAFF) by running the Antechamber module of the Amber software package. The GAFF parameters for the CC-SN and CS-NO dihedrals were not accurate and thus determined anew. The CC-SN and CS-NO torsional energy profiles of EtSNO were calculated quantum mechanically at the level of B3LYP/cc-pVTZ//HF/6-$31G^*$. Torsional force constants were obtained by fitting the theoretical torsional energies with those obtained from molecular mechanics energy minimization. These parameters for EtSNO reproduced, to a reasonable accuracy, the corresponding torsional energy profiles of the capped tripeptide ACE-CysSNO-NME as well as their structures obtained from quantum mechanical geometry optimization. A molecular dynamics simulation of myoglobin with a CysSNO residue produced a well-behaved trajectory demonstrating that the parameters may be used in modeling other S-nitrosated proteins.

• Journal of the Korean Magnetic Resonance Society
• /
• v.22 no.4
• /
• pp.139-143
• /
• 2018

Short antimicrobial peptide, A4W, have been studied by molecular dynamics (MD) simulation in an explicit dodecylphosphocholine (DPC) micelle. Peptide was aligned with DPC micelle and transferred new peptide-micelle coordinates within the same solvent box using specific micelle topology parameters. After initial energy minimization and equilibration, the conformation and orientation of the peptide were analyzed from trajectories obtained from the RMD (restrained molecular dynamics) or the subsequent free MD. Also, the information of solvation in the backbone and the side chain of the peptide, hydrogen bonding, and the properties of the dynamics were obtained. The results showed that the backbone residues of peptide are either solvated using water or in other case, they relate to hydrogen bonding. These properties could be a critical factor against the insertion mode of interaction. Most of the peptide-micelle interactions come from the hydrophobic interaction between the side chains of peptide and the structural interior of micelle system. The interaction of peptide-micelle, electrostatic potential and hydrogen bonding, between the terminal residues of peptide and the headgroups in micelle were observed. These interactions could be effect on the structure and flexibility of the peptide terminus.

• Bulletin of the Korean Chemical Society
• /
• v.31 no.9
• /
• pp.2581-2587
• /
• 2010

Nitration of tyrosine and tryptophan residues is common in cells under nitrative stress. However, physiological consequences of protein nitration are not well characterized on a molecular level due to limited availability of the 3D structures of nitrated proteins. Molecular dynamics (MD) simulation can be an alternative tool to probe the structural perturbations induced by nitration. In this study we developed molecular mechanics parameters for 3-nitrotyrosine (NIY) and 6-nitrotryptophan (NIW) that are compatible with the AMBER-99 force field. Partial atomic charges were derived by using a multi-conformational restrained electrostatic potential (RESP) methodology that included the geometry optimized structures of both $\alpha$- and $\beta$-conformers of a capped tripeptide ACE-NIY-NME or ACE-NIW-NME. Force constants for bonds and angles were adopted from the generalized AMBER force field. Torsional force constants for the proper dihedral C-C-N-O and improper dihedral C-O-N-O of the nitro group in NIY were determined by fitting the torsional energy profiles obtained from quantum mechanical (QM) geometry optimization with those from molecular mechanical (MM) energy minimization. Force field parameters obtained for NIY were transferable to NIW so that they reproduced the QM torsional energy profiles of ACE-NIW-NME accurately. Moreover, the QM optimized structures of the tripeptides containing NIY and NIW were almost identical to the corresponding structures obtained from MM energy minimization, attesting the validity of the current parameter set. Molecular dynamics simulations of thioredoxin nitrated at the single tyrosine and tryptophan yielded well-behaved trajectories suggesting that the parameters are suitable for molecular dynamics simulations of a nitrated protein.