Biomedical subjects

Jakob P Ulmschneider

Publications and source records attributed to Jakob P Ulmschneider.

3 recordsLinked to original sources

PubMed2007

A generalized born implicit-membrane representation compared to experimental insertion free energies.

An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.

Computer Simulation↗

PubMed2006

Monte Carlo vs molecular dynamics for all-atom polypeptide folding simulations.

An efficient Monte Carlo (MC) algorithm including concerted rotations is directly compared to molecular dynamics (MD) in all-atom statistical mechanics folding simulations of small polypeptides. The previously reported algorithm "concerted rotations with flexible bond angles" (CRA) has been shown to successfully locate the native state of small polypeptides. In this study, the folding of three small polypeptides (trpzip2/H1/Trp-cage) is investigated using MC and MD, for a combined sampling time of approximately 10(11) MC configurations and 8 micros, respectively. Both methods successfully locate the experimentally determined native states of the three systems, but they do so at different speed, with 2-2.5 times faster folding of the MC runs. The comparison reveals that thermodynamic and dynamic properties can reliably be obtained by both and that results from folding simulations do not depend on the algorithm used. Similar to previous comparisons of MC and MD, it is found that one MD integration step of 2 fs corresponds to one MC scan, revealing the good sampling of MC. The simplicity and efficiency of the MC method will enable its future use in folding studies involving larger systems and the combination with replica exchange algorithms.

Computer Simulation↗

PubMed2004

Polypeptide folding using Monte Carlo sampling, concerted rotation, and continuum solvation.

An efficient concerted rotation algorithm for use in Monte Carlo statistical mechanics simulations is applied to fold three polypeptides, U(1-17)T9D, alpha(1), and trpzip2, which exhibit native beta-hairpin and alpha-helix folds. The method includes flexible bond and dihedral angles, and a Gaussian bias is applied with driver bond and dihedral angles to optimize the sampling efficiency. Solvation in water is implemented with the generalized Born (GBSA) model. The computed lowest-energy manifolds for the folded structures of the two beta-hairpins agree closely with the corresponding NMR structures. In the case of the alpha(1) peptide, the folded alpha-helical state, which is observed as oligomers in concentrated solution and crystals, is not stable in isolation. The computed preference for random coil structures is in agreement with NMR experiments at low concentration. The fact that native states can be located on high dimensional energy surfaces starting from extended conformations shows that the present methodology samples all relevant parts of the conformational space. The OPLS-AA force field with the GBSA solvent model was also found to perform well in leading to clear energetic separation of the correctly folded structures from misfolded structures for the two peptides that form beta-turns.

Algorithms↗