Biomedical subjects

Jooyoung Lee

Publications and source records attributed to Jooyoung Lee.

16 recordsLinked to original sources

PubMed2006

Application of the multiensemble sampling to the equilibrium folding of proteins.

MOTIVATION: Conventional Monte Carlo and molecular dynamics simulations of proteins in the canonical ensemble are of little use, because they tend to get trapped in states of energy local minima at low temperatures. One way to surmount this difficulty is to use a non-Boltzmann sampling method in which conformations are sampled upon a general weighting function instead of the conventional Boltzmann weighting function. The multiensemble sampling (MES) method is a non-Boltzmann sampling method that was originally developed to estimate free energy differences between systems with different potential energies and/or at different thermodynamic states. The method has not yet been applied to studies of complex molecular systems such as proteins. RESULTS: MES Monte Carlo simulations of small proteins have been carried out using a united-residue force field. The proteins at several temperatures from the unfolded to the folded states were simulated in a single MC run at a time and their equilibrium thermodynamic properties were calculated correctly. The distributions of sampled conformations clearly indicate that, when going through states of energy local minima, the MES simulation did not get trapped in them but escaped from them so quickly that all the relevant parts of conformation space could be sampled properly. A two-step folding process consisting of a collapse transition followed by a folding transition is observed. This study demonstrates that the use of MES alleviates the multiple-minima problem greatly. AVAILABILITY: Available on request from the authors.

Algorithms↗

PubMed2005

Study of protein-protein interaction using conformational space annealing.

We apply conformational space annealing (CSA), an efficient global optimization method, to the study of protein-protein interaction. The CSA is incorporated into the Tinker molecular modeling package along with a B-spline method for CAPRI Round 5 experiments. We have used an energy function for the protein-protein interaction that consists of electrostatic interaction, van der Waals interaction, and solvation energy terms represented by the occupancy desolvation method. The parameters of the AMBER94 all-atom empirical force field are used. Each energy term is calculated by precalculated grid potentials and B-spline method approximation. The ligand protein is placed inside a sphere of 50 A radius centered at an appropriate location, and the CSA rigid docking studies are carried out to find stable complexes. Up to 10 complexes are selected using the K-mean clustering method and biological information when available. These complexes are energy-minimized for further refinement by considering the flexibility of interacting proteins. The results show that the CSA method has a potential for the study of protein-protein interaction.

Algorithms↗

PubMed2005

Structure optimization by conformational space annealing in an off-lattice protein model.

The optimization results by conformational space annealing are presented for an off-lattice protein model consisting of hydrophobic and hydrophilic residues in Fibonacci sequences. The ground-state energies found are lower than those reported in the literature. In addition, the ground-state conformations in three dimensions exhibit the important aspect of forming a single hydrophobic core in real proteins. The energy landscape for the population of local minima is also investigated.

Algorithms↗

PubMed2005

PPRODO: prediction of protein domain boundaries using neural networks.

Successful prediction of protein domain boundaries provides valuable information not only for the computational structure prediction of multidomain proteins but also for the experimental structure determination. Since protein sequences of multiple domains may contain much information regarding evolutionary processes such as gene-exon shuffling, this information can be detected by analyzing the position-specific scoring matrix (PSSM) generated by PSI-BLAST. We have presented a method, PPRODO (Prediction of PROtein DOmain boundaries) that predicts domain boundaries of proteins from sequence information by a neural network. The network is trained and tested using the values obtained from the PSSM generated by PSI-BLAST. A 10-fold cross-validation technique is performed to obtain the parameters of neural networks using a nonredundant set of 522 proteins containing 2 contiguous domains. PPRODO provides good and consistent results for the prediction of domain boundaries, with accuracy of about 66% using the +/-20 residue criterion. The PPRODO source code, as well as all data sets used in this work, are available from http://gene.kias.re.kr/ approximately jlee/pprodo/.

Amino Acids↗

PubMed2005

Topological determinants of protein unfolding rates.

For proteins that fold by two-state kinetics, the folding and unfolding processes are believed to be closely related to their native structures. In particular, folding and unfolding rates are influenced by the native structures of proteins. Thus, we focus on finding important topological quantities from a protein structure that determine its unfolding rate. After constructing graphs from protein native structures, we investigate the relationships between unfolding rates and various topological quantities of the graphs. First, we find that the correlation between the unfolding rate and the contact order is not as prominent as in the case of the folding rate and the contact order. Next, we investigate the correlation between the unfolding rate and the clustering coefficient of the graph of a protein native structure, and observe no correlation between them. Finally, we find that a newly introduced quantity, the impact of edge removal per residue, has a good overall correlation with protein unfolding rates. The impact of edge removal is defined as the ratio of the change of the average path length to the edge removal probability. From these facts, we conclude that the protein unfolding process is closely related to the protein native structure.

Binding Sites↗

PubMed2005

An efficient molecular docking using conformational space annealing.

Molecular docking falls into the general category of global optimization problems because its main purpose is to find the most stable complex consisting of a receptor and its ligand. Conformational space annealing (CSA), a powerful global optimization method, is incorporated with the Tinker molecular modeling package to perform molecular docking simulations of six receptor-ligand complexes (3PTB, 1ULB, 2CPP, 1STP, 3CPA, and 1PPH) from the Protein Data Bank. In parallel, Monte Carlo with the minimization (MCM) method is also incorporated into the Tinker package for comparison. The energy function, consisting of electrostatic interactions, van der Waals interactions, and torsional energy terms, is calculated using the AMBER94 all-atom empirical force field. Rigid docking simulations for all six complexes and flexible docking simulations for three complexes (1STP, 3CPA, and 1PPH) are carried out using the CSA and the MCM methods. The simulation results show that the docking procedures using the CSA method generally find the most stable complexes as well as the native-like complexes more efficiently and accurately than those using the MCM, demonstrating that CSA is a promising search method for molecular docking problems.

Algorithms↗

PubMed2005

Understanding the structural characteristics of compstatin by conformational space annealing.

The structural characteristics of the 13-residue compstatin molecule are investigated using the conformational space annealing (CSA) method with CHARMM force field and the GBSA continuum solvent model. In order to sample conformations in the energy range of the minimized NMR structures, we have used the stopping criterion to the CSA search when a conformation whose energy is less than -490 kcal/mol is found. With this stopping criterion, a great variety of conformations are generated around experimentally known structures. Twenty independent CSA runs starting from random states find 1000 representative conformations in the energy landscape of the compstatin, which are classified into thirty-one structural families. The majority of the conformations (94.4%) are in the coil state. Other conformers containing a 3(10)-helix, a pi-helix, a beta-hairpin, and an alpha-helix are also found.

Models, Molecular↗

PubMed2005

Folding simulations of small proteins.

Understanding how a protein folds is a long-standing challenge in modern science. We have used an optimized atomistic model (united-residue force field) to simulate folding of small proteins of various structures: HP-36 (alpha protein), protein A (beta), 1fsd (alpha+beta), and betanova (beta). Extensive Monte Carlo folding simulations (ten independent runs with 10(9) Monte Carlo steps at a temperature) starting from non-native conformations are carried out for each protein. In all cases, proteins fold into their native-like conformations at appropriate temperatures, and glassy transitions occur at low temperatures. To investigate early folding trajectories, 200 independent runs with 10(6) Monte Carlo steps are also performed at a fixed temperature for a protein. There are a variety of possible pathways during non-equilibrium early processes (fast process, approximately 10(4) Monte Carlo steps). Finally, these pathways converge to the point unique for each protein. The convergence point of the early folding pathways can be determined only by direct folding simulations. The free energy surface, an equilibrium thermodynamic property, dictates the rest of the folding (slow process, approximately 10(8) Monte Carlo steps).

Computer Simulation↗

PubMed2005

Protein structure prediction based on fragment assembly and parameter optimization.

We propose a novel method for ab-initio prediction of protein tertiary structures based on the fragment assembly and global optimization. Fifteen residue long fragment libraries are constructed using the secondary structure prediction method PREDICT, and fragments in these libraries are assembled to generate full-length chains of a query protein. Tertiary structures of 50 to 100 conformations are obtained by minimizing an energy function for proteins, using the conformational space annealing method that enables one to sample diverse low-lying local minima of the energy. Then in order to enhance the performance of the prediction method, we optimize the linear parameters of the energy function, so that the native-like conformations become energetically more favorable than the non-native ones for proteins with known structures. We test the feasibility of the parameter optimization procedure by applying it to the training set consisting of three proteins: the 10-55 residue fragment of staphylococcal protein A (PDB ID 1bdd), a designed protein betanova, and 1fsd.

Models, Molecular↗

PubMed2004

Prediction of protein tertiary structure using PROFESY, a novel method based on fragment assembly and conformational space annealing.

A novel method for ab initio prediction of protein tertiary structures, PROFESY (PROFile Enumerating SYstem), is proposed. This method utilizes the secondary structure prediction information of a query sequence and the fragment assembly procedure based on global optimization. Fifteen-residue-long fragment libraries are constructed using the secondary structure prediction method PREDICT, and fragments in these libraries are assembled to generate full-length chains of a query protein. Tertiary structures of 50 to 100 conformations are obtained by minimizing an energy function for proteins, using the conformational space annealing method that enables one to sample diverse low-lying local minima of the energy. We apply PROFESY for benchmark tests to proteins with known structures to demonstrate its feasibility. In addition, we participated in CASP5 and applied PROFESY to four new-fold targets for blind prediction. The results are quite promising, despite the fact that PROFESY was in its early stages of development. In particular, PROFESY successfully provided us the best model-one structure for the target T0161.

Caspases↗

PubMed2004

Folding of small proteins using a single continuous potential.

Extensive Monte Carlo folding simulations for four proteins of various structural classes are carried out, using a single continuous potential (united-residue force field). In all cases, collapse occurs at a very early stage, and proteins fold into their nativelike conformations at appropriate temperatures. We also observe that glassy transitions occur at low temperatures. The simulation results demonstrate that the folding mechanism is controlled not only by thermodynamic factors but also by kinetic factors: The way a protein folds into its native structure is also determined by the convergence point of early folding trajectories, which cannot be obtained by the free energy surface.

Biophysics↗

PubMed2004

Dynamic pathway model for the formation of C(60).

We present a dynamic pathway model for the formation of C(60) using the action-derived molecular dynamics simulations. We propose candidate precursors for dynamic pathway models in which carbons spontaneously aggregate due to favorable energetics and kinetics. Various planar polycyclic models are in a disadvantageous state where they cannot be trapped in the forward reaction due to their high excess internal energies. Our simulation results show that precursors either in the shape of tangled polycyclics or in the shape of open cages are kinetically favored over precursors in the shape of planar hexagonal graphite fragments. Calculated activation energies for the probable precursor models are in good agreement with experiment. Existence of chains in the models of tangled polycyclics and open cages is beneficially for the formation of C(60) molecule. Chains attached to the precursor model are energetically favorable and display lithe movements along the dynamic pathway.

Journal Article↗

PubMed2003

Unbiased global optimization of Lennard-Jones clusters for N < or =201 using the conformational space annealing method.

We apply the conformational space annealing method to the Lennard-Jones clusters and find all known lowest energy configurations up to 201 atoms, without using extra information of the problem such as the structures of the known global energy minima. In addition, the robustness of the algorithm with respect to the randomness of initial conditions of the problem is demonstrated by ten successful independent runs up to 183 atoms. Our results indicate that this method is a general and yet efficient global optimization algorithm applicable to many systems.

Journal Article↗

PubMed2003

Vortex patterns and infinite degeneracy in the uniformly frustrated XY models and lattice Coulomb gas.

We find analytically and numerically that there exist infinite degeneracies in the ground state vortex configurations of a uniformly frustrated XY model on a square lattice for cases of frustrations f=1/6, 1/7, 1/8, and 1/10. More generally, we could also obtain a class of exact solutions for the phase configurations corresponding to the cases of f=1/q=1/(2m) with an integer m>or=3 of which the above cases of m=3, 4, and 5 correspond to the true ground states. These states are analogous to the staircase solutions originally given by Halsey [Phys. Rev. B 31, 5728 (1985)], which, however, are relevant in a different regime of the frustration parameter (or vortex density). In these quasistaircase states, all the gauge-invariant phase differences are found to be integer multiples of pi/2q. The supercurrents are conserved at each node in a trivial manner by separate channels. The infinite ground state degeneracy is preserved in the case of an arbitrary screening length for the corresponding lattice Coulomb gas.

Journal Article↗

PubMed2002

Lignin peroxidase oxidation of veratryl alcohol: effects of the mutants H82A, Q222A, W171A, and F267L.

The site-directed mutations H82A and Q222A (residues near the heme access channel), and W171A and F267L (residues near the surface of the protein) were introduced into the gene encoding lignin peroxidase (LiP) isozyme H8 from Phanerochaete chrysosporium. The variant enzymes were produced by homologous expression in P. chrysosporium, purified to homogeneity, and characterized by kinetic and spectroscopic methods. The molecular masses, the pIs, and the UV-vis absorption spectra of the ferric and oxidized states of these LiP variant enzymes were similar to those of wild-type LiP (wtLiP), suggesting the overall protein and heme environments were not significantly affected by these mutations. The steady-state and transient-state parameters for the oxidation of veratryl alcohol (VA) by the H82A and Q222A variants were very similar to those of wtLiP, demonstrating that these residues are not involved in VA oxidation and that the heme access channel is an unlikely site for VA oxidation. In contrast, the W171A variant was unable to oxidize VA, confirming the apparent essentiality of Trp171 in VA oxidation by LiP. The kinetic rates of spontaneous LiP compound I reduction in the absence of VA were similar for W171A and wild-type LiP, suggesting that there may not be a radical formed on the Trp171 residue of LiP in the absence of VA. For the F267L variant, both the K(m app) value in the steady state and the apparent dissociation constant (K(D)) for compound II reduction were greater than those for wtLiP. These results indicate that the site including W171 and F267, rather than the heme access channel, is the site of VA binding and oxidation in LiP. Whereas Trp171 appears to be essential for VA oxidation, it apparently is not independently responsible for the spontaneous decomposition of oxidized intermediates. The nearby Phe267 apparently is also involved in VA binding.

Benzyl Alcohols↗

PubMed2002

Evolution of physics-based methodology for exploring the conformational energy landscape of proteins.

The evolution of our physics-based computational methods for determining protein conformation without the introduction of secondary-structure predictions, homology modeling, threading, or fragment coupling is described. Initial use of a hard-sphere potential captured much of the structural properties of polypeptide chains, and subsequent more refined force fields, together with efficient methods of global optimization provide indications that progress is being made toward an understanding of the interresidue interactions that underlie protein folding.

Computational Biology↗