Search PubMed⌕ Search

Biomedical subjects

Youngshang Pak

Publications and source records attributed to Youngshang Pak.

7 recordsLinked to original sources

Direct folding simulation of alpha-helices and beta-hairpins based on a single all-atom force field with an implicit solvation model.

Recently, we have shown that a modified energy model based on the param99 force field with the generalized Born (GB) solvation model produces reliable free energy landscapes of mini-proteins with a betabetaalpha motif (BBA5, 1FSD, and 1PSV), with the native structures of the mini-proteins located in their lowest free energy minimum states. One of the main features in the modified energy model is a significant improvement for more balanced treatments of alpha and beta strands in proteins. In this study, using the replica exchange molecular dynamics (REMD) simulation method with this new force field, we have carried out extensive ab initio folding studies of several well-known peptides with alpha or beta strands (C-peptide, EK-peptide, le0q, and gbl). Starting from fully extended conformations as the initial conditions, all of the native-like structures of the target peptides were successfully identified by REMD, with reasonable representations of free energy surfaces. The present simulation results with the modified energy model are consistent with experiments, demonstrating an extended applicability of the energy model to folding studies of a variety of alpha-helices, beta-strands, and alpha/beta proteins.

Amino Acid Sequence↗

Free energy surfaces of miniproteins with a betabetaalpha motif: replica exchange molecular dynamics simulation with an implicit solvation model.

Designed miniproteins with a betabetaalpha motif, such as BBA5, 1FSD, and 1PSV can serve as a benchmark set to test the validity of all-atom force fields with computer simulation, because they contain all the basic structural elements in protein folding. Unfortunately, it was found that the standard all-atom force fields with the generalized Born (GB) implicit solvation model tend to produce distorted free energy surfaces for the betabetaalpha proteins, not only because energetically those proteins need to be described by more balanced weights of the alpha- and beta-strands, but also because the GB implicit solvation model suffers from overestimated salt bridge effects. In an attempt to resolve these problems, we have modified one of the standard all-atom force fields in conjunction with the GB model, such that each native state of the betabetaalpha proteins is in its free energy minimum state with reasonable energy barriers separating local minima. With this modified energy model, the free energy contour map in each protein was constructed from the replica exchange molecular dynamics REMD simulation. The resulting free energy surfaces are significantly improved in comparison with previous simulation results and consistent with general views on small protein folding behaviors with realistic topology and energetics of all three proteins.

Calorimetry↗

Relationship between the ratio of ligand to metal and the coordinating ability of anions. Synthesis and structural properties of AgX-bearing bis(4-pyridyl)dimethylsilane (X- = NO2-, NO3-, CF3SO3-, and PF6-).

Studies of the anion effects on the molecular construction of a series of AgX complexes with bis(4-pyridyl)dimethylsilane (L) (X- = NO2-, NO3-, CF3SO3-, and PF6-) have been carried out. Formation of the skeletal bonds appears to be primarily associated with a suitable combination of bidentate N-donors of L and a variety of coordination geometries of Ag(I) ions. The L:Ag(I) ratios of the products are dependent on the nature of the polyatomic anions. The 1:1 adduct Ag(I)-L for NO2-, 3:4 adduct for NO3-, 2:3 adduct for CF3SO3-, and 1:2 adduct for PF6- have been obtained. A linear relationship between the ratio of ligand to metal and the coordinating ability of anions was observed. [Ag(NO2)(L)] has a unique sheet structure consisting of double helices, and [Ag3(L)4](NO3)3 is a 2 nm thick interwoven sheet structure consisting of nanotubes. The compound [Ag2(L)3](CF3SO3)2 affords a characteristic ladder-type channel structure, and [Ag(L)2](PF6) is a simple 2D grid structure.

Journal Article↗

Ab initio folding of helix bundle proteins using molecular dynamics simulations.

We have demonstrated that ab initio fast folding simulations at 400 K using a GB implicit solvent model with an all-atom based force field can describe the spontaneous formation of nativelike structures for the 36-residue villin headpiece and the 46-residue fragment B of Staphylococcal protein A. An implicit solvent model combined with high-temperature MD makes it possible to perform direct folding simulations of small- to medium-sized proteins by reducing the computational requirements tremendously. In the early stage of folding of the villin headpiece and protein A, initial hydrophobic collapse and rapid formation of helices were found to play important roles. For protein A, the third helix forms first in the early stage of folding and exhibits higher stability. The free energy profiles calculated from the folding simulations suggested that both of the helix-bundle proteins show a two-state thermodynamic behavior and protein A exhibits rather broad native basins.

Computer Simulation↗

Replica-exchange method using the generalized effective potential.

We propose an effective scheme for fast conformational searches by combining the replica exchange method (REM) with the generalized effect potential concept. The present method introduces the "q" value from the effective potential as a coupling parameter. It is found that the new method not only requires a much smaller number of replicas than the conventional REM, but also makes it possible to perform effective conformational sampling of complex systems with correct distributions maintained. The advantage of the present method has been demonstrated with in vacuo alanine dipeptide using a molecular dynamics simulation.

Computer Simulation↗

Molecular dynamics study of peptides in implicit water: ab initio folding of beta-hairpin, beta-sheet, and beta beta alpha-motif.

In this communication, we have demonstrated that molecular dynamics simulations using a GB implicit solvation model with the all-atom based force field (CHARMM19) can describe the spontaneous folding of small peptides in aqueous solution. The native structures of peptides with various structural motifs (beta-hairpin, beta-sheet, and betabetaalpha-moiety) were successfully predicted within reasonable time scales by MD simulations at moderately elevated temperatures. It is expected that the present simulations provide further insight into mechanism/pathways of the peptide folding.

Amino Acid Sequence↗