PubMed Health⌕ Search

Biomedical subjects

Hugh Nymeyer

Publications and source records attributed to Hugh Nymeyer.

6 recordsLinked to original sources
PubMed2006

Indole localization in lipid membranes revealed by molecular simulation.

It is commonly known that the amino acid residue tryptophan and its side-chain analogs, e.g., indole, are strongly attracted to the interfacial region of lipid bilayers. Phenylalanine and its side-chain analogs, e.g., benzene, do not localize in the interface but are distributed throughout the lipid bilayer. We use molecular dynamics to investigate the details of indole and benzene localization and orientation within a POPC bilayer and the factors that lead to their different properties. We identify three sites in the bilayer at which indole is localized: 1), a site in the interface near the glycerol moiety; 2), a weakly bound site in the interface near the choline moiety; and 3), a weakly bound site in the center of the bilayer's hydrocarbon core. Benzene is localized in the same three positions, but the most stable position is the hydrocarbon core followed by the site near the glycerol moiety. Transfer of indole from water to the hydrocarbon core shows a classic hydrophobic effect. In contrast, interfacial binding is strongly enthalpy driven. We use several different sets of partial charges to investigate the factors that contribute to indole's and benzene's orientational and spatial distribution. Our simulations show that a number of electrostatic interactions appear to contribute to localization, including hydrogen bonding to the lipid carbonyl groups, cation-pi interactions, interactions between the indole dipole and the lipid bilayer's strong interfacial electric field, and nonspecific electrostatic stabilization due to a mismatch in the variation of the nonpolar forces and local dielectric with position in the bilayer.

Benzene↗
PubMed2005

Folding is not required for bilayer insertion: replica exchange simulations of an alpha-helical peptide with an explicit lipid bilayer.

We implement the replica exchange molecular dynamics algorithm to study the interactions of a model peptide (WALP-16) with an explicitly represented DPPC membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16 into the DPPC bilayer and its folding into an alpha-helix with a transbilayer orientation. The free energy surface suggests that the insertion of the peptide into the DPPC bilayer precedes secondary structure formation. Although the peptide has some propensity to form a partially helical structure in the interfacial region of the DPPC/water system, this state is not a productive intermediate but rather an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that the observed insertion/folding pathway mirrors the potential of mean force (PMF). Calculation of the enthalpic and entropic contributions to this PMF show that the surface bound conformation of WALP-16 is significantly lower in energy than other conformations, and that the insertion of WALP-16 into the bilayer without regular secondary structure is enthalpically unfavorable by 5-10 kcal/mol/residue. The observed insertion/folding pathway disagrees with the dominant conceptual model, which is that a surface-bound helix is an obligatory intermediate for the insertion of alpha-helical peptides into lipid bilayers. In our simulations, the observed insertion/folding pathway is favored because of a large (>100 kcal/mol) increase in system entropy that occurs when the unstructured WALP-16 peptide enters the lipid bilayer interior. The insertion/folding pathway that is lowest in free energy depends sensitively on the near cancellation of large enthalpic and entropic terms. This suggests the possibility that intrinsic membrane peptides may have a diversity of insertion/folding behaviors depending on the exact system of peptide and lipid under consideration.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2003

Simulation of the folding equilibrium of alpha-helical peptides: a comparison of the generalized Born approximation with explicit solvent.

We compare simulations using the generalized Born/surface area (GB/SA) implicit solvent model with simulations using explicit solvent (transferable intermolecular potential 3 point, TIP3P) to test the GB/SA algorithm. We use the replica exchange molecular dynamics method to sample the conformational phase space of two alpha-helical peptides, A21 and the Fs, by using two different classical potentials and both water models. We find that when using GB/SA: (i) A21 is predicted to be more helical than the Fs peptide at all temperatures; (ii) the native structure of the Fs peptide is predicted to be a helical bundle instead of a single helix; and (iii) the persistence length and most probable end-to-end distance are too large in the unfolded state when compared against the explicit solvent simulations. We find that the potential of mean force in the phi(psi) plane is markedly different in the two solvents, making the two simulated peptides respond differently when the backbone torsions are perturbed. A fit of the temperature melting curves obtained in these simulations to a Lifson-Roig model finds that the GB/SA model has an unphysically large nucleation parameter, whereas the explicit solvent model produces values similar to experiment.

Algorithms↗
PubMed2003

Peptide folding simulations.

Developments in the design of small peptides that mimic proteins in complexity, recent advances in nanosecond time-resolved spectroscopy methods to study peptides and the development of modern, highly parallel simulation algorithms have come together to give us a detailed picture of peptide folding dynamics. Two newly implemented simulation techniques, parallel replica dynamics and replica exchange molecular dynamics, can now describe directly from simulations the kinetics and thermodynamics of peptide formation, respectively. Given these developments, the simulation community now has the tools to verify and validate simulation protocols and models (forcefields).

Computer Simulation↗
PubMed2002

Specific and nonspecific collapse in protein folding funnels.

Experiments with fast folding proteins are beginning to address the relationship between collapse and folding. We investigate how different scenarios for folding can arise depending on whether the folding and collapse transitions are concurrent or whether a nonspecific collapse precedes folding. Many earlier studies have focused on the limit in which collapse is fast compared to the folding time; in this work we focus on the opposite limit where, at the folding temperature, collapse and folding occur simultaneously. Real proteins exist in both of these limits. The folding mechanism varies substantially in these two regimes. In the regime of concurrent folding and collapse, nonspecific collapse now occurs at a temperature below the folding temperature (but slightly above the glass transition temperature).

Amino Acid Sequence↗