• 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.

• BMB Reports
• /
• v.31 no.1
• /
• pp.89-94
• /
• 1998

The structures of the self-complementary decamer duplexes, $d(ACGTATACGT)_2$ (TATA-duplex) and $d(ACGTTAACGT)_2$, (TTAA-duplex) has been obtained in solution by proton NMR spectroscopy and restrained molecular dynamics. The duplexes are essentially B-type, with distortions apparent at the TATA and TTAA steps. Theses distortions and their effects on dynamics have been investigated by the measurement of imino proton exchange time of the base-pairs. The unusual opening kinetics of central A T base-pairs could be correlated to the abnormal structural properties of the corresponding sequences.

• 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.

• Coupled systems mechanics
• /
• v.1 no.3
• /
• pp.269-283
• /
• 2012

In the present paper, a multilevel approach for the local nanobuckling analysis of carbon nanotube (CNT) based composite materials is proposed and described. The approach comprises four levels, all of them at nanoscale. The first level aims to propose the potential that describes the interatomic forces between carbon atoms. In the second level, molecular dynamics simulations are performed to extract the elastic properties of the CNT. The third level aims to determine the stiffness of the material that surrounds the CNT (matrix), using the annular membrane analysis. In the fourth level, finite strip analysis of the CNT elastically restrained by the matrix is performed to calculate the critical strain at which the CNT buckles locally. In order to achieve accurate results and take the CNT-matrix interaction into account, the $3^{rd}$ and $4^{th}$ steps may be repeated iteratively until convergence is achieved. The proposed multilevel approach is applied to several CNTs embedded in a cylindrical representative volume element and illustrated in detail. It shows that (i) the interaction between the CNT and the matrix should be taken into account and (ii) the buckling at nanoscale is sensitive to several types of local buckling modes.

• Proceedings of the Korean Biophysical Society Conference
• /
• 1996.07a
• /
• pp.29-29
• /
• 1996

The structure of an antimicrobial peptide buforin II has been studied by $^1$H NMR and CD spectroscopy. Buforin II is flexible and random structure in H$_2$O but the parts of buforin II become helical structure in trifluoroethanol (TFE)/H$_2$O (1:1, v/v) solution. From the restrained molecular dynamics calculation using NMR data, we obtained the possible conformation of buforin II in TFE/H$_2$O solution. (omitted)

• Archives of Pharmacal Research
• /
• v.25 no.1
• /
• pp.11-24
• /
• 2002

The antitumor antibiotic (+)-CC-1065 can alkylate N3 of guanine in certain sequences. A previous high-field $^1H$ NMR study on the$(+)-CC-1065d[GCGCAATTG*CGC]_2$ adduct ($^*$ indicates the drug alkylation site) showed that drag modification on N3 of guanine results in protonation of the cross-strand cytosine [Park, H-J.; Hurley, L. H. J. Am. Chem. Soc.1997, 119,629]. In this contribution we describe a further analysis of the NMR data sets together with restrained molecular dynamics. This study provides not only a solution structure of the (+)-CC-1065(N3- guanine) DNA duplex adduct but also new insight into the molecular basis for the sequence- specific interaction between (+)-CC-1065 and N3-guanine in the DNA duplex. On the basis of NOESY data, we propose that the narrow minor groove at the 7T8T step and conformational kinks at the junctions of 16C17A and 18A19T are both related to DNA bending in the drugDNA adduct. Analysis of the one-dimensional $^1H$ NMR (in $H_2O$) data and rMD trajectories strongly suggests that hydrogen bonding linkages between the 8-OH group of the (+)-CC-1065 A-sub-unit and the 9G10C phosphate via a water molecule are present. All the phenomena observed here in the (+)-CC-1065(N3-guanine) adduct at 5'$-AATTG^*$are reminiscent of those obtained from the studies on the (+)-CC-1065(N3-adenine) adduct at $5'-AGTTA^*$, suggesting that (+)-CC-1065 takes advantage of the conformational flexibility of the 5'-TPu step to entrap the bent structure required for the covalent bonding reaction. This study reveals a common molecular basis for (+)-CC-1065 alkylation at both $5'-TTG^*$ and $5'-TTA^*$, which involves a trapping out of sequence-dependent DNA conformational flexibility as well as sequence-dependent general acid and general base catalysis by duplex DNA.

• Journal of the Korean Magnetic Resonance Society
• /
• v.26 no.1
• /
• pp.1-9
• /
• 2022

There are two distinct approaches to improving the quality of protein NMR structures during refinement: all-atom force fields and accumulated knowledge-assisted methods that include Rosetta. Mao et al. reported that, for 40 proteins, Rosetta increased the accuracies of their NMR-determined structures with respect to the X-ray crystal structures (Mao et al., J. Am. Chem. Soc. 136, 1893 (2014)). In this study, we calculated 32 structures of those studied by Mao et al. using all-atom force field and implicit solvent model, and we compared the results with those obtained from Rosetta. For a single protein, using only the experimental NOE-derived distances and backbone torsion angle restraints, 20 of the lowest energy structures were extracted as an ensemble from 100 generated structures. Restrained simulated annealing by molecular dynamics simulation searched conformational spaces with a total time step of 1-ns. The use of GPU-accelerated AMBER code allowed the calculations to be completed in hours using a single GPU computer-even for proteins larger than 20 kDa. Remarkably, statistical analyses indicated that the structures determined in this way showed overall higher accuracies to their X-ray structures compared to those refined by Rosetta (p-value < 0.01). Our data demonstrate the capability of sophisticated atomistic force fields in refining NMR structures, particularly when they are coupled with the latest GPU-based calculations. The straightforwardness of the protocol allows its use to be extended to all NMR structures.

• Archives of Pharmacal Research
• /
• v.24 no.5
• /
• pp.455-465
• /
• 2001

The highly potent cytotoxic DNA-DNA cross-linker consists of two cyclopropa[c]pyrrolo[3,4-3]indol-4(5H)-ones insoles [(+)-CPI-I] joined by a bisamido pyrrole (abbreviated to "Pyrrole"). The Pyrrole is a synthetic analog of Bizelesin, which is currently in phase II clinical trials due to its excellent in vivo antitumor activity. The Pyrrole has 10 times more potent cytotoxicity than Bizelesin and mostly form DNA-DNA interstrand cross-links through the N3 of adenines spaced 7 bp apart. The Pyrrole requires a centrally positioned GC base pair for high cross-linking reactivity (i.e., $5^1$-T$AT_2$A*-$3^1$), while Bizelesin prefers purely AT-rich sequences (i.e., $5^1$-T$AT_4$A*-$3^1$, where /(equation omitted) represents the cross-strand adenine alkylation and A* represents an adenine alkylation) (Park et al., 1996). In this study, the high-field $^1$H-NMR and rMD studies are conducted on the 1 1-mer DNA duplex adduct of the Pyrrole where the 5′(equation omitted)TAGTTA*-3′sequence is cross-linked by the drug. A severe structural perturbation is observed in the intervening sequences of cross-linking site, while a normal B-DNA structure is maintained in the region next to the drug-modified adenines. Based upon these observations, we propose that the interplay between the bisamido pyrrole unit of the drug and central C/C base pair (hydrogen-bonding interactions) is involved in the process of cross-linking reaction, and sequence specificity is the outcome of those interactions. This study suggests a mechanism for the sequence specific cross-linking reaction of the Pyrrole, and provides a further insight to develop new DNA sequence selective and distortive cross-linking agents.