Biomedical subjects

John E Straub

Publications and source records attributed to John E Straub.

16 recordsLinked to original sources

PubMed2006

Monomer adds to preformed structured oligomers of Abeta-peptides by a two-stage dock-lock mechanism.

Nonfibrillar soluble oligomers, which are intermediates in the transition from monomers to amyloid fibrils, may be the toxic species in Alzheimer's disease. To monitor the early events that direct assembly of amyloidogenic peptides we probe the dynamics of formation of (Abeta(16-22))(n) by adding a monomer to a preformed (Abeta(16-22))(n-1) (n = 4-6) oligomer in which the peptides are arranged in an antiparallel beta-sheet conformation. All atom molecular dynamics simulations in water and multiple long trajectories, for a cumulative time of 6.9 mus, show that the oligomer grows by a two-stage dock-lock mechanism. The largest conformational change in the added disordered monomer occurs during the rapid ( approximately 50 ns) first dock stage in which the beta-strand content of the monomer increases substantially from a low initial value. In the second slow-lock phase, the monomer rearranges to form in register antiparallel structures. Surprisingly, the mobile structured oligomers undergo large conformational changes in order to accommodate the added monomer. The time needed to incorporate the monomer into the fluid-like oligomer grows even when n = 6, which suggests that the critical nucleus size must exceed six. Stable antiparallel structure formation exceeds hundreds of nanoseconds even though frequent interpeptide collisions occur at elevated monomer concentrations used in the simulations. The dock-lock mechanism should be a generic mechanism for growth of oligomers of amyloidogenic peptides.

Amyloid beta-Peptides↗

PubMed2006

Dynamics of Asp23-Lys28 salt-bridge formation in Abeta10-35 monomers.

In the amyloid fibrils formed from long fragments of the amyloid beta-protein (Abeta-protein), the monomers are arranged in parallel and lie perpendicular to the fibril axis. The structure of the monomers satisfies the amyloid self-organization principle; namely, the low free energy state of the monomer maximizes the number of intra- and interpeptide contacts and salt bridges. The formation of the intramolecular salt bridge between Asp(D)23 and Lys(K)28 ensures that unpaired charges are not buried in the low-dielectric interior. We have investigated, using all-atom molecular dynamics simulations in explicit water, whether the D23-K28 interaction forms spontaneously in the isolated Abeta10-35 monomer. To validate the simulation protocol, we show, using five independent trajectories spanning a total of 100 ns, that the pKa values of the titratable groups are in good agreement with experimental measurements. The computed free energy disconnectvity graph shows that broadly the ensemble of compact random coil conformations can be clustered into four basins that are separated by free energy barriers ranging from 0.3 to 2.7 kcal/mol. There is significant residual structure in the conformation of the peptide in each of the basins. Due to the desolvation penalty, the structural motif with a stable turn involving the residues VGSN and a preformed D23-K28 contact is a minor component of the simulated structures. The extent of solvation of the peptides in the four basins varies greatly, which underscores the dynamical fluctuations in the monomer. Our results suggest that the early event in the oligomerization process must be the expulsion of discrete water molecules that facilitates the formation of interpeptide-interaction-driven stable structures with an intramolecular D23-K28 salt bridge and an intact VGSN turn.

Amyloid beta-Peptides↗

PubMed2006

Statistical-temperature Monte Carlo and molecular dynamics algorithms.

A simulation method is presented that achieves a flat energy distribution by updating the statistical temperature instead of the density of states in Wang-Landau sampling. A novel molecular dynamics algorithm (STMD) applicable to complex systems and a Monte Carlo algorithm are developed from this point of view. Accelerated convergence for large energy bins, essential for large systems, is demonstrated in tests on the Ising model, the Lennard-Jones fluid, and bead models of proteins. STMD shows a superior ability to find local minima in proteins and new global minima are found for the 55 bead AB model in two and three dimensions. Calculations of the occupation probabilities of individual protein inherent structures provide new insights into folding and misfolding.

Journal Article↗

PubMed2006

Efforts toward developing direct probes of protein dynamics.

We report the first IR characterization of a single C-D bond within a protein, methyl-d1 Met80 of horse heart cytochrome c. A comparison was made to methyl-d1/d3 methionine as well as methyl-d3 Met80. We found that for methyl-d1 and the asymmetric stretches of methyl-d3, line widths/line shapes are dominated by inhomogeneous broadening, whereas the symmetric stretch of methyl-d3 has a significant homogeneous component. Vibrational energy relaxation calculations found that a significantly stronger Fermi resonance exists for the symmetric stretch than for the asymmetric stretches, thereby suggesting that a difference in intramolecular vibrational relaxation (IVR) causes the observed line width/line shape difference between the symmetric and asymmetric stretches.

Animals↗

PubMed2006

Time-dependent perturbation theory for vibrational energy relaxation and dephasing in peptides and proteins.

Without invoking the Markov approximation, we derive formulas for vibrational energy relaxation (VER) and dephasing for an anharmonic system oscillator using a time-dependent perturbation theory. The system-bath Hamiltonian contains more than the third order coupling terms since we take a normal mode picture as a zeroth order approximation. When we invoke the Markov approximation, our theory reduces to the Maradudin-Fein formula which is used to describe the VER properties of glass and proteins. When the system anharmonicity and the renormalization effect due to the environment vanishes, our formulas reduce to those derived by and Mikami and Okazaki [J. Chem. Phys. 121, 10052 (2004)] invoking the path-integral influence functional method with the second order cumulant expansion. We apply our formulas to VER of the amide I mode of a small amino-acid like molecule, N-methylacetamide, in heavy water.

Amides↗

PubMed2005

Vibrational energy relaxation in proteins.

An overview of theories related to vibrational energy relaxation (VER) in proteins is presented. VER of a selected mode in cytochrome c is studied by using two theoretical approaches. One approach is the equilibrium simulation approach with quantum correction factors, and the other is the reduced model approach, which describes the protein as an ensemble of normal modes interacting through nonlinear coupling elements. Both methods result in similar estimates of the VER time (subpicoseconds) for a CD stretching mode in the protein at room temperature. The theoretical predictions are in accord with previous experimental data. A perspective on directions for the detailed study of time scales and mechanisms of VER in proteins is presented.

Adenosine Triphosphate↗

PubMed2005

Monte Carlo simulations of polyalanine using a reduced model and statistics-based interaction potentials.

A coarse-grained residue-residue interaction potential derived from a statistical analysis of the Protein Data Bank is used to investigate the coil-to-helix transition for polyalanine. The interaction potentials depend on the radial distance between interaction sites, as well as the relative orientation of the sites. Two types of interaction sites are present in the model: a site representing the amino acid side chain, and a site representing a "virtual backbone," i.e., a site located in the peptide bond which accounts for backbone hydrogen bonding. Two chain lengths are studied and the results for the thermodynamics of the coil-to-helix transition are analyzed in terms of the Zimm-Bragg model. Results agree qualitatively and quantitatively with all-atom Monte Carlo simulations and other reduced-model Monte Carlo simulations.

Computer Simulation↗

PubMed2004

Probing the initial stage of aggregation of the Abeta(10-35)-protein: assessing the propensity for peptide dimerization.

Characterization of the early stages of peptide aggregation is of fundamental importance in elucidating the mechanism of the formation of deposits associated with amyloid disease. The initial step in the pathway of aggregation of the Abeta-protein, whose monomeric NMR structure is known, was studied through the simulation of the structure and stability of the peptide dimer in aqueous solution. A protocol based on shape complementarity was used to generate an assortment of possible dimer structures. The structures generated based on shape complementarity were evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure. The potential of mean force associated with the dimerization of the peptides in aqueous solution was computed for both the hydrophobic and the electrostatic driven forces using umbrella sampling and classical molecular dynamics simulation at constant temperature and pressure with explicit solvent and periodic boundary conditions. The comparison of the two free energy profiles suggests that the structure of the peptide dimer is determined by the favorable desolvation of the hydrophobic residues at the interface. Molecular dynamics trajectories originating from two putative dimer structures indicate that the peptide dimer is stabilized primarily through hydrophobic interactions, while the conformations of the peptide monomers undergo substantial structural reorganization in the dimerization process. The finding that the phi-dimer may constitute the ensemble of stable Abeta(10-35) dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Abeta(16-22) fragments. We conjecture that events prior to the nucleation process themselves might involve crossing free energy barriers which depend on the peptide-peptide and peptide-water interactions. Consistent with existing experimental studies, the peptides within the ensemble of aggregated states show no signs of formation of secondary structure.

Amyloid beta-Peptides↗

PubMed2004

Aqueous urea solution destabilizes Abeta(16-22) oligomers.

We use long multiple trajectories generated by molecular dynamics simulations to probe the stability of oligomers of Abeta(16-22) (KLVFFAE) peptides in aqueous urea solution. High concentration of urea promotes the formation of beta-strand structures in Abeta(16-22) monomers, whereas in water they adopt largely compact random coil structures. The tripeptide system, which forms stable antiparallel beta-sheet structure in water, is destabilized in urea solution. The enhancement of beta-strand content in the monomers and the disruption of oligomeric structure occur largely by direct interaction of urea with the peptide backbone. Our simulations suggest that the oligomer unbinding dynamics is determined by two opposing effects, namely, by the increased propensity of monomers to form beta-strands and the rapid disruption of the oligomers. The qualitative conclusions are affirmed by using two urea models. Because the proposed destabilization mechanism depends largely on hydrogen bond formation between urea and the peptide backbone, we predict that high urea concentration will destabilize oligomers of other amyloidogenic peptides as well.

Amyloid beta-Peptides↗

PubMed2004

Orientational potentials extracted from protein structures improve native fold recognition.

We develop coarse-grained, distance- and orientation-dependent statistical potentials from the growing protein structural databases. For protein structural classes (alpha, beta, and alpha/beta), a substantial number of backbone-backbone and backbone-side-chain contacts stabilize the native folds. By taking into account the importance of backbone interactions with a virtual backbone interaction center as the 21st anisotropic site, we construct a 21 x 21 interaction scheme. The new potentials are studied using spherical harmonics analysis (SHA) and a smooth, continuous version is constructed using spherical harmonic synthesis (SHS). Our approach has the following advantages: (1) The smooth, continuous form of the resulting potentials is more realistic and presents significant advantages for computational simulations, and (2) with SHS, the potential values can be computed efficiently for arbitrary coordinates, requiring only the knowledge of a few spherical harmonic coefficients. The performance of the new orientation-dependent potentials was tested using a standard database of decoy structures. The results show that the ability of the new orientation-dependent potentials to recognize native protein folds from a set of decoy structures is strongly enhanced by the inclusion of anisotropic backbone interaction centers. The anisotropic potentials can be used to develop realistic coarse-grained simulations of proteins, with direct applications to protein design, folding, and aggregation.

Amino Acids↗

PubMed2003

A molecular switch in amyloid assembly: Met35 and amyloid beta-protein oligomerization.

Aberrant protein oligomerization is an important pathogenetic process in vivo. In Alzheimer's disease (AD), the amyloid beta-protein (Abeta) forms neurotoxic oligomers. The predominant in vivo Abeta alloforms, Abeta40 and Abeta42, have distinct oligomerization pathways. Abeta42 monomers oligomerize into pentamer/hexamer units (paranuclei) which self-associate to form larger oligomers. Abeta40 does not form these paranuclei, a fact which may explain the particularly strong linkage of Abeta42 with AD. Here, we sought to determine the structural elements controlling paranucleus formation as a first step toward the development of strategies for treating AD. Because oxidation of Met(35) is associated with altered Abeta assembly, we examined the role of Met(35) in controlling Abeta oligomerization. Oxidation of Met(35) in Abeta42 blocked paranucleus formation and produced oligomers indistinguishable in size and morphology from those produced by Abeta40. Systematic structural alterations of the C(gamma)(35)-substituent group revealed that its electronic nature, rather than its size (van der Waals volume), was the factor controlling oligomerization pathway choice. Preventing assembly of toxic Abeta42 paranuclei through selective oxidation of Met(35) thus represents a potential therapeutic approach for AD.

Amyloid↗

PubMed2003

Structural and dynamical analysis of the hydration of the Alzheimer's beta-amyloid peptide.

An analysis of the water molecules in the first solvation shell obtained from the molecular dynamics simulation of the amyloid beta(10-35)NH2 peptide and the amyloid beta(10-35)NH2E22Q "Dutch" mutant peptide is presented. The structure, energetics, and dynamics of water in the hydration shell have been investigated using a variety of measures, including the hydrogen bond network, the water residence times for all the peptide residues, the diffusion constant, experimentally determined HN amide proton exchange, and the transition probabilities for water to move from one residue to another or into the bulk. The results of the study indicate that: (1) the water molecules at the peptide-solvent interface are organized in an ordered structure similar for the two peptide systems but different from that of the bulk, (2) the peptide structure inhibits diffusion perpendicular to the peptide surface by a factor of 3 to 5 relative to diffusion parallel to the peptide surface, which is comparable to diffusion of bulk water, (3) water in the first solvation shell shows dynamical relaxation on fast (1-2 ps) and slow (10-40 ps) time scales, (4) a novel solvent relaxation master equation is shown to capture the details of the fast relaxation of water in the peptide's first solvation shell, (5) the interaction between the peptide and the solvent is stronger in the wild type than in the E22Q mutant peptide, in agreement with earlier results obtained from computer simulations [Massi, F.; Straub, J. E. Biophys J 2001, 81, 697] correlated with the observed enhanced activity of the E22Q mutant peptide.

Alzheimer Disease↗

PubMed2003

Vibrational frequency shifts and relaxation rates for a selected vibrational mode in cytochrome C.

The vibrational energy relaxation of a selected vibrational mode in cytochrome c--a C-D stretch in the terminal methyl group of Met80--has been studied using equilibrium molecular dynamics simulation and normal mode analysis methods. As demonstrated in the pioneering work of Romesberg and co-workers, isotopic labeling of the C-H (to C-D) stretch in alkyl side chains shifts the stretching frequency to the transparent region of the protein's density of states, making it an effective and versatile probe of protein structure and dynamics. Molecular dynamics trajectories of solvated cytochrome c were run at 300 K, and vibrational population relaxation times were estimated using the classical Landau-Teller-Zwanzig model and a number of semiclassical theories of resonant and two-phonon vibrational relaxation processes. The C-D stretch vibrational population relaxation time is estimated to be T(1) = 14-40 ps; the relatively close agreement between various semiclassical estimates of T(1) lends support to the applicability of those expressions. Normal mode calculations were used to identify the dominant coupling between the protein and C-D oscillator. All bath modes strongly coupled to the C-D stretch are in close proximity. Angle bending modes in the terminal methyl group of Met80 appear to be the most likely acceptor modes defining the mechanism of population relaxation of the C-D vibration.

Animals↗

PubMed2002

Gravitational smoothing as a global optimization strategy.

An optimization scheme for atomic cluster structures, based on exaggerating the importance of the gravitational force, is introduced. Results are presented for calculations on Lennard-Jones clusters of 13, 38, and 55 atoms, and the 13-atom Morse cluster.

Journal Article↗

PubMed2002

Long time dynamic simulations: exploring the folding pathways of an Alzheimer's amyloid Abeta-peptide.

We describe the MaxFlux algorithm for the computation of likely pathways for global macromolecular conformational transitions. The algorithm assumes an overdamped diffusive dynamics for the biomolecule, appropriate to large scale conformational changes. As an application of the MaxFlux method, we explore conformational transitions between alpha-helical, collapsed coil, and beta-sheet conformations of an amyloid Abeta-peptide. The resulting transition pathways are analyzed in terms of the mechanism of conformational transition and the progression of the peptide energetics in both an aqueous and a membrane-mimicking nonpolar solvent.

Amyloid beta-Peptides↗

PubMed2002

Charge states rather than propensity for beta-structure determine enhanced fibrillogenesis in wild-type Alzheimer's beta-amyloid peptide compared to E22Q Dutch mutant.

The activity of the Alzheimer's amyloid beta-peptide is a sensitive function of the peptide's sequence. Increased fibril elongation rate of the E22Q Dutch mutant of the Alzheimer's amyloid beta-peptide relative to that of the wild-type peptide has been observed. The increased activity has been attributed to a larger propensity for the formation of beta structure in the monomeric E22Q mutant peptide in solution relative to the WT peptide. That hypothesis is tested using four nanosecond timescale simulations of the WT and Dutch mutant forms of the Abeta(10-35)-peptide in aqueous solution. The simulation results indicate that the propensity for formation of beta-structure is no greater in the E22Q mutant peptide than in the WT peptide. A significant measure of "flickering" of helical structure in the central hydrophobic cluster region of both the WT and mutant peptides is observed. The simulation results argue against the hypothesis that the Dutch mutation leads to a higher probability of formation of beta-structure in the monomeric peptide in aqueous solution. We propose that the greater stability of the solvated WT peptide relative to the E22Q mutant peptide leads to decreased fibril elongation rate in the former. Stability difference is due to the differing charge state of the two peptides. The other proposal leads to the prediction that the fibril elongation rates for the WT and the mutant E22Q should be similar under acid conditions.

Amino Acid Substitution↗