Biomedical subjects

Bosco K Ho

Publications and source records attributed to Bosco K Ho.

5 recordsLinked to original sources

PubMed2006

Folding very short peptides using molecular dynamics.

Peptides often have conformational preferences. We simulated 133 peptide 8-mer fragments from six different proteins, sampled by replica-exchange molecular dynamics using Amber7 with a GB/SA (generalized-Born/solvent-accessible electrostatic approximation to water) implicit solvent. We found that 85 of the peptides have no preferred structure, while 48 of them converge to a preferred structure. In 85% of the converged cases (41 peptides), the structures found by the simulations bear some resemblance to their native structures, based on a coarse-grained backbone description. In particular, all seven of the beta hairpins in the native structures contain a fragment in the turn that is highly structured. In the eight cases where the bioinformatics-based I-sites library picks out native-like structures, the present simulations are largely in agreement. Such physics-based modeling may be useful for identifying early nuclei in folding kinetics and for assisting in protein-structure prediction methods that utilize the assembly of peptide fragments.

Biophysical Phenomena↗

PubMed2005

The Ramachandran plots of glycine and pre-proline.

BACKGROUND: The Ramachandran plot is a fundamental tool in the analysis of protein structures. Of the 4 basic types of Ramachandran plots, the interactions that determine the generic and proline Ramachandran plots are well understood. The interactions of the glycine and pre-proline Ramachandran plots are not. RESULTS: In glycine, the psi angle is typically clustered at psi = 180 degrees and psi = 0 degrees. We show that these clusters correspond to conformations where either the N(i+1) or O atom is sandwiched between the two Halpha atoms of glycine. We show that the shape of the 5 distinct regions of density (the alpha, alphaL, betaS, betaP and betaPR regions) can be reproduced with electrostatic dipole-dipole interactions. In pre-proline, we analyse the origin of the zeta region of the Ramachandran plot, a region unique to pre-proline. We show that it is stabilized by a CO(i-1)...CdeltaHdelta(i+1) weak hydrogen bond. This is analogous to the CO(i-1)...NH(i+1) hydrogen bond that stabilizes the gamma region in the generic Ramachandran plot. CONCLUSION: We have identified the specific interactions that affect the backbone of glycine and pre-proline. Knowledge of these interactions will improve current force-fields, and help understand structural motifs containing these residues.

Amino Acid Motifs↗

PubMed2005

The flexibility in the proline ring couples to the protein backbone.

In proteins, the proline ring exists predominantly in two discrete states. However, there is also a small but significant amount of flexibility in the proline ring of high-resolution protein structures. We have found that this side-chain flexibility is coupled to the backbone conformation. To study this coupling, we have developed a model that is simply based on geometric and steric factors and not on energetics. We show that the coupling between phi and chi1 torsions in the proline ring can be described by an analytic equation that was developed by Bricard in 1897, and we describe a computer algorithm that implements the equation. The model predicts the observed coupling very well. The strain in the C(gamma)-C(delta)-N angle appears to be the principal barrier between the UP and DOWN pucker. This strain is relaxed to allow the proline ring to flatten in the rare PLANAR conformation.

Algorithms↗

PubMed2003

Revisiting the Ramachandran plot: hard-sphere repulsion, electrostatics, and H-bonding in the alpha-helix.

What determines the shape of the allowed regions in the Ramachandran plot? Although Ramachandran explained these regions in terms of 1-4 hard-sphere repulsions, there are discrepancies with the data where, in particular, the alphaR, alphaL, and beta-strand regions are diagonal. The alphaR-region also varies along the alpha-helix where it is constrained at the center and the amino terminus but diffuse at the carboxyl terminus. By analyzing a high-resolution database of protein structures, we find that certain 1-4 hard-sphere repulsions in the standard steric map of Ramachandran do not affect the statistical distributions. By ignoring these steric clashes (NH(i+1) and O(i-1)C), we identify a revised set of steric clashes (CbetaO, O(i-1)N(i+1), CbetaN(i+1), O(i-1)Cbeta, and O(i-1)O) that produce a better match with the data. We also find that the strictly forbidden region in the Ramachandran plot is excluded by multiple steric clashes, whereas the outlier region is excluded by only one significant steric clash. However, steric clashes alone do not account for the diagonal regions. Using electrostatics to analyze the conformational dependence of specific interatomic interactions, we find that the diagonal shape of the alphaR and alphaL-regions also depends on the optimization of the NH(i+1) and O(i-1)C interactions, and the diagonal beta-strand region is due to the alignment of the CO and NH dipoles. Finally, we reproduce the variation of the Ramachandran plot along the alpha-helix in a simple model that uses only H-bonding constraints. This allows us to rationalize the difference between the amino terminus and the carboxyl terminus of the alpha-helix in terms of backbone entropy.

Amino Acids↗

PubMed2002

Twist and shear in beta-sheets and beta-ribbons.

The structures of the beta-sheets and the beta-ribbons have been analysed using high-resolution protein structure data. Systematic asymmetries measured in both parallel and antiparallel beta-structures include the sheet twist and the strand shear. In order to determine the origin of these asymmetries, numerous interactions and correlations were examined. The strongest correlations are observed for residues in antiparallel beta-sheets and beta-ribbons that form non-H-bonded pairs. For these residues, the sheet twist is correlated to the backbone phi angle but not to the psi angle. Our analysis supports the existence of an inter-strand C(alpha)H(alpha)...O weak H-bond, which, together with the CO...HN H-bond, constitutes a bifurcated H-bond that links neighbouring beta-strands. Residues of beta-sheets and beta-ribbons in high-resolution protein structures form a distinct region of the Ramachandran plot, which is determined by the formation of the bifurcated H-bond, the formation of an intra-strand O...H(alpha) non-bonded polar interaction, and an intra-strand O...C(beta) steric clash. Using beta-strands parameterised by phi-psi values from the allowed beta-sheet region of the Ramachandran plot, the shear and the right-hand twist can be reproduced in a simple model of the antiparallel and parallel beta-ribbon that models the bifurcated H-bonds specifically. The conformations of interior residues of beta-sheets are shown to be subsets of the conformations of residues of beta-ribbons.

Databases, Protein↗