Search PubMed⌕ Search

Biomedical subjects

Ruhong Zhou

Publications and source records attributed to Ruhong Zhou.

21 records · Page 2Linked to original sources

Can a continuum solvent model reproduce the free energy landscape of a beta -hairpin folding in water?

The folding free energy landscape of the C-terminal beta-hairpin of protein G is explored using the surface-generalized Born (SGB) implicit solvent model, and the results are compared with the landscape from an earlier study with explicit solvent model. The OPLSAA force field is used for the beta-hairpin in both implicit and explicit solvent simulations, and the conformational space sampling is carried out with a highly parallel replica-exchange method. Surprisingly, we find from exhaustive conformation space sampling that the free energy landscape from the implicit solvent model is quite different from that of the explicit solvent model. In the implicit solvent model some nonnative states are heavily overweighted, and more importantly, the lowest free energy state is no longer the native beta-strand structure. An overly strong salt-bridge effect between charged residues (E42, D46, D47, E56, and K50) is found to be responsible for this behavior in the implicit solvent model. Despite this, we find that the OPLSAA/SGB energies of all the nonnative structures are higher than that of the native structure; thus the OPLSAA/SGB energy is still a good scoring function for structure prediction for this beta-hairpin. Furthermore, the beta-hairpin population at 282 K is found to be less than 40% from the implicit solvent model, which is much smaller than the 72% from the explicit solvent model and approximately equal 80% from experiment. On the other hand, both implicit and explicit solvent simulations with the OPLSAA force field exhibit no meaningful helical content during the folding process, which is in contrast to some very recent studies using other force fields.

Animals↗

Development of a polarizable force field for proteins via ab initio quantum chemistry: first generation model and gas phase tests.

We present results of developing a methodology suitable for producing molecular mechanics force fields with explicit treatment of electrostatic polarization for proteins and other molecular system of biological interest. The technique allows simulation of realistic-size systems. Employing high-level ab initio data as a target for fitting allows us to avoid the problem of the lack of detailed experimental data. Using the fast and reliable quantum mechanical methods supplies robust fitting data for the resulting parameter sets. As a result, gas-phase many-body effects for dipeptides are captured within the average RMSD of 0.22 kcal/mol from their ab initio values, and conformational energies for the di- and tetrapeptides are reproduced within the average RMSD of 0.43 kcal/mol from their quantum mechanical counterparts. The latter is achieved in part because of application of a novel torsional fitting technique recently developed in our group, which has already been used to greatly improve accuracy of the peptide conformational equilibrium prediction with the OPLS-AA force field.1 Finally, we have employed the newly developed first-generation model in computing gas-phase conformations of real proteins, as well as in molecular dynamics studies of the systems. The results show that, although the overall accuracy is no better than what can be achieved with a fixed-charges model, the methodology produces robust results, permits reasonably low computational cost, and avoids other computational problems typical for polarizable force fields. It can be considered as a solid basis for building a more accurate and complete second-generation model.

Algorithms↗

Detecting native protein folds among large decoy sets with hydrophobic moment profiling.

A new hydrophobic score will be presented in this paper for detecting native-like folds from a large set of decoy structures. A recent paper (B. D. Silverman, PNAS 98, 4996, 2001) had revealed that for globular proteins there exists a relatively universal constant of 0.74 for a hydrophobic ratio, which is defined as the ratio of radii from the protein centroid at which the second order hydrophobic moment and the zero order moment vanishes. This paper further defines a new hydrophobic score which will be used to examine protein decoys, in particular, the Holm & Sander, Park & Levitt and Baker decoy sets. It will be shown that the hydrophobic score and profile shapes can provide useful information that should be complementary to the information provided by other procedures, such as free energy calculations.

Amino Acids↗