Highly polar environments catalyze the unfolding of PrP C helix 1

The first alpha-helix (H1) likely plays an important role in the conversion of the cellular prion protein (PrP(C)) into its pathogenic isoform (PrP(Sc)). In this conversion, H1 may either have to unfold or may represent a site of intermolecular contact. A recent molecular dynamics simulation suggested that H1 can unfold if it is detached from the protein core (Hirschberger et al. in Biophys J 90:3908, 2006). It has been hypothesized that the high dielectric constant epsilon (S) of the bulk water environment facilitates the unfolding of H1. To check this hypothesis, we performed a number of replica exchange molecular dynamics simulations of an H1 peptide in solvents of different epsilon (S). We found that the equilibrium helix fraction in water is less than 40%, in agreement with previous experimental findings, and that the helix unfolds much faster in water than in less polar solvents. The kinetically stabilizing effect of the organic solvents is largely unspecific and correlates well with their dielectric constant epsilon (S).

Substituents at the c(13) position of retinal and their influence on the function of bacteriorhodopsin

Retinal analogues in which the 13-methyl group is replaced by H, C(2)H(5), CF(3), and OCH(3) residues are studied by means of quantumchemical modified neglect of diatomic overlap-correlated version (MNDOC) calculations. The analogues are suitable to test the stereochemical mechanism of proton pumping in bacteriorhodopsin. The results explain the proton-pumping activities of bacterio-opsin reconstituted with these analogues and elucidate the decisive role of retinal's ground-state intramolecular properties in the pump cycle of bacteriorhodopsin.

Structural investigation of bacteriorhodopsin and some of its photoproducts by polarized Fourier transform infrared spectroscopic methods-difference spectroscopy and photoselection

The direction of selected IR-transition moments of the retinal chromophore of bacteriorhodopsin (BR) and functional active amino acid residues are determined for light- and dark-adapted BR and for the intermediates K and L of the photocycle. Torsions around single bonds of the chromophore are found to be present in all the investigated BR states. The number of twisted single bonds and the magnitude of these torsions decreases in the order K, L, light-adapted BR, dark-adapted BR. In the last, only the C(14)-C(15) single bond is twisted. The orientation of molecular planes and chemical bonds of such protein side chains, which are perturbed during the transition of light-adapted BR to the respective intermediates, are deduced and the results compared with the current three dimensional model of BR. Trp 86 and Trp 185 are found to form a rigid part of the protein, whereas Asp 96 and Asp 115 perform molecular rearrangements upon formation of the L-intermediate.

Evidence for a 13,14-cis cycle in bacteriorhodopsin

We discuss to what extent the vibrational spectra of bacteriorhodopsin that have been observed and assigned by Smith et al. (1, 2) by means of resonance Raman and by Gerwert and Siebert (EMBO (Eur. Mol. Biol. Organ.) J. In press) by means of infrared absorption experiments are in agreement with a photo-cycle of bacteriorhodopsin that involves the sequence BR, IO(all-trans) --> K(13,14-cis) --> L(13,14-cis) --> M(13-cis) --> N(13-cis) --> O(all-trans). Our discussion is based on a quantumchemical modified neglect of diatomic overlap [MNDO] calculation of the vibrational spectra of the relevant isomers of the protonated retinal Schiff base. In particular, we investigated in these calculations the effects of different charge environments on the frequencies of the relevant C-C single bond stretching vibrations of these isomers.

A polarizable force field for computing the infrared spectra of the polypeptide backbone

The shapes of the amide bands in the infrared (IR) spectra of proteins and peptides are caused by electrostatically coupled vibrations within the polypeptide backbone and code the structures of these biopolymers. A structural decoding of the amide bands has to resort to simplified models because the huge size of these macromolecules prevents the application of accurate quantum mechanical methods such as density functional theory (DFT). Previous models employed transition-dipole coupling methods that are of limited accuracy. Here we propose a concept for the computation of protein IR spectra, which describes the molecular mechanics (MM) of polypeptide backbones by a polarizable force field of "type II". By extending the concepts of conventional polarizable MM force fields, such a PMM/II approach employs field-dependent parameters not only for the electrostatic signatures of the molecular components but also for the local potentials modeling the stiffness of chemical bonds with respect to elongations, angle deformations, and torsions. Using a PMM/II force field, the IR spectra of the polypeptide backbone can be efficiently calculated from the time dependence of the backbone's dipole moment during a short (e.g., 100 ps) MD simulation by Fourier transformation. PMM/II parameters are derived for harmonic bonding potentials of amide groups in polypeptides from a series of DFT calculations on the model molecule N-methylacetamide (NMA) exposed to homogeneous external electric fields. The amide force constants are shown to vary by as much as 20% for relevant field strengths. As a proof of principle, it is shown that the large solvatochromic effects observed in the IR spectra of NMA upon transfer from the gas phase into aqueous solution are not only excellently reproduced by DFT/MM simulations but are also nicely modeled by the PMM/II approach. The tasks remaining for a proof of practice are specified.

Light-triggered beta-hairpin folding and unfolding

A light-switchable peptide is transformed with ultrashort pulses from a beta-hairpin to an unfolded hydrophobic cluster and vice versa. The structural changes are monitored by mid-IR probing. Instantaneous normal mode analysis with a Hamiltonian combining density functional theory with molecular mechanics is used to interpret the absorption transients. Illumination of the beta-hairpin state triggers an unfolding reaction that visits several intermediates and reaches the unfolded state within a few nanoseconds. In this unfolding reaction to the equilibrium hydrophobic cluster conformation, the system does not meet significant barriers on the free-energy surface. The reverse folding process takes much longer because it occurs on the time scale of 30 micros. The folded state has a defined structure, and its formation requires an extended search for the correct hydrogen-bond pattern of the beta-strand.

Electrostatics of proteins in dielectric solvent continua. II. First applications in molecular dynamics simulations

In the preceding paper by Stork and Tavan, [J. Chem. Phys. 126, 165105 (2007)], the authors have reformulated an electrostatic theory which treats proteins surrounded by dielectric solvent continua and approximately solves the associated Poisson equation [B. Egwolf and P. Tavan, J. Chem. Phys. 118, 2039 (2003)]. The resulting solution comprises analytical expressions for the electrostatic reaction field (RF) and potential, which are generated within the protein by the polarization of the surrounding continuum. Here the field and potential are represented in terms of Gaussian RF dipole densities localized at the protein atoms. Quite like in a polarizable force field, also the RF dipole at a given protein atom is induced by the partial charges and RF dipoles at the other atoms. Based on the reformulated theory, the authors have suggested expressions for the RF forces, which obey Newton's third law. Previous continuum approaches, which were also built on solutions of the Poisson equation, used to violate the reactio principle required by this law, and thus were inapplicable to molecular dynamics (MD) simulations. In this paper, the authors suggest a set of techniques by which one can surmount the few remaining hurdles still hampering the application of the theory to MD simulations of soluble proteins and peptides. These techniques comprise the treatment of the RF dipoles within an extended Lagrangian approach and the optimization of the atomic RF polarizabilities. Using the well-studied conformational dynamics of alanine dipeptide as the simplest example, the authors demonstrate the remarkable accuracy and efficiency of the resulting RF-MD approach.

Electrostatics of proteins in dielectric solvent continua. I. Newton’s third law marries qE forces

The authors reformulate and revise an electrostatic theory treating proteins surrounded by dielectric solvent continua [B. Egwolf and P. Tavan, J. Chem. Phys. 118, 2039 (2003)] to make the resulting reaction field (RF) forces compatible with Newton's third law. Such a compatibility is required for their use in molecular dynamics (MD) simulations, in which the proteins are modeled by all-atom molecular mechanics force fields. According to the original theory the RF forces, which are due to the electric field generated by the solvent polarization and act on the partial charges of a protein, i.e., the so-called qE forces, can be quite accurately computed from Gaussian RF dipoles localized at the protein atoms. Using a slightly different approximation scheme also the RF energies of given protein configurations are obtained. However, because the qE forces do not account for the dielectric boundary pressure exerted by the solvent continuum on the protein, they do not obey the principle that actio equals reactio as required by Newton's third law. Therefore, their use in MD simulations is severely hampered. An analysis of the original theory has led the authors now to a reformulation removing the main difficulties. By considering the RF energy, which represents the dominant electrostatic contribution to the free energy of solvation for a given protein configuration, they show that its negative configurational gradient yields mean RF forces obeying the reactio principle. Because the evaluation of these mean forces is computationally much more demanding than that of the qE forces, they derive a suggestion how the qE forces can be modified to obey Newton's third law. Various properties of the thus established theory, particularly issues of accuracy and of computational efficiency, are discussed. A sample application to a MD simulation of a peptide in solution is described in the following paper [M. Stork and P. Tavan, J. Chem. Phys., 126, 165106 (2007).

Identification of anti-prion compounds as efficient inhibitors of polyglutamine protein aggregation in a zebrafish model

Several neurodegenerative diseases, including Huntington disease (HD), are associated with aberrant folding and aggregation of polyglutamine (polyQ) expansion proteins. Here we established the zebrafish, Danio rerio, as a vertebrate HD model permitting the screening for chemical suppressors of polyQ aggregation and toxicity. Upon expression in zebrafish embryos, polyQ-expanded fragments of huntingtin (htt) accumulated in large SDS-insoluble inclusions, reproducing a key feature of HD pathology. Real time monitoring of inclusion formation in the living zebrafish indicated that inclusions grow by rapid incorporation of soluble htt species. Expression of mutant htt increased the frequency of embryos with abnormal morphology and the occurrence of apoptosis. Strikingly, apoptotic cells were largely devoid of visible aggregates, suggesting that soluble oligomeric precursors may instead be responsible for toxicity. As in nonvertebrate polyQ disease models, the molecular chaperones, Hsp40 and Hsp70, suppressed both polyQ aggregation and toxicity. Using the newly established zebrafish model, two compounds of the N'-benzylidene-benzohydrazide class directed against mammalian prion proved to be potent inhibitors of polyQ aggregation, consistent with a common structural mechanism of aggregation for prion and polyQ disease proteins.

Vibrational spectra from atomic fluctuations in dynamics simulations. II. Solvent-induced frequency fluctuations at femtosecond time resolution

The midinfrared (MIR) spectra of molecules in polar solvents exhibit inhomogeneously broadened bands whose spectral positions are shifted as compared to the gas phase. The shifts are caused by interactions with structured solvation shells and the broadenings by fluctuations of these interactions. The MIR spectra can be calculated from hybrid molecular dynamics (MD) simulations, which treat the solute molecule by density functional theory and the solvent by molecular mechanics by the so-called instantaneous normal mode analysis (INMA) or by Fourier transforming the time correlation function (FTTCF) of the molecular dipole moment. In Paper I of this work [M. Schmitz and P. Tavan, J. Chem. Phys. 121, 12233 (2004)] we explored an alternative method based on generalized virial (GV) frequencies noting, however, that GV systematically underestimates frequencies. As shown by us these artifacts are caused by solvent-induced fluctuations of the (i) equilibrium geometry, (ii) force constants, and (iii) normal mode directions as well as by (iv) diagonal and (v) off-diagonal anharmonicities. Here we now show, by analyzing the time scales of fluctuations and sample MD trajectories of formaldehyde in the gas phase and in water, that all these sources of computational artifacts can be made visible by a Fourier analysis of the normal coordinates. Correspondingly, the error sources (i) and (iii)-(v) can be removed by bandpass filtering, as long as the spectral signatures of the respective effects are well separated from the fundamental band. Furthermore, the artifacts arising from effect (ii) can be strongly diminished by a time-resolved version of the GV approach (TF-GV). The TF-GV method then yields for each mode j a trajectory of the vibrational frequency omega(j)(tmid R:tau) at a time resolution tau>tau(j), which is only limited by the corresponding oscillation time tau(j)=2pi/omega(j) and, thus, is in the femtosecond range. A correlation analysis of these trajectories clearly separates the librational motions from the conformational dynamics of the solvation shells and yields the inhomogeneously broadened MIR spectra, if the theory of motional narrowing is properly included. The MIR spectrum of formaldehyde in solution obtained by TF-GV agrees very well with the FTTCF result, if one applies the so-called "harmonic approximation" quantum correction factor and a temperature scaling to the FTTCF intensities. Also for INMA an excellent agreement is achieved if one disregards a slight INMA overestimate of linewidths.

Structural instability of the prion protein upon M205S/R mutations revealed by molecular dynamics simulations

The point mutations M205S and M205R have been demonstrated to severely disturb the folding and maturation process of the cellular prion protein (PrP(C)). These disturbances have been interpreted as consequences of mutation-induced structural changes in PrP, which are suggested to involve helix 1 and its attachment to helix 3, because the mutated residue M205 of helix 3 is located at the interface of these two helices. Furthermore, current models of the prion protein scrapie (PrP(Sc)), which is the pathogenic isoform of PrP(C) in prion diseases, imply that helix 1 disappears during refolding of PrP(C) into PrP(Sc). Based on molecular-dynamics simulations of wild-type and mutant PrP(C) in aqueous solution, we show here that the native PrP(C) structure becomes strongly distorted within a few nanoseconds, once the point mutations M205S and M205R have been applied. In the case of M205R, this distortion is characterized by a motion of helix 1 away from the hydrophobic core into the aqueous environment and a subsequent structural decay. Together with experimental evidence on model peptides, this decay suggests that the hydrophobic attachment of helix 1 to helix 3 at M205 is required for its correct folding into its stable native structure.

Systematic identification of antiprion drugs by high-throughput screening based on scanning for intensely fluorescent targets

Conformational changes and aggregation of specific proteins are hallmarks of a number of diseases, like Alzheimer's disease, Parkinson's disease, and prion diseases. In the case of prion diseases, the prion protein (PrP), a neuronal glycoprotein, undergoes a conformational change from the normal, mainly alpha-helical conformation to a disease-associated, mainly beta-sheeted scrapie isoform (PrP(Sc)), which forms amyloid aggregates. This conversion, which is crucial for disease progression, depends on direct PrP(C)/PrP(Sc) interaction. We developed a high-throughput assay based on scanning for intensely fluorescent targets (SIFT) for the identification of drugs which interfere with this interaction at the molecular level. Screening of a library of 10,000 drug-like compounds yielded 256 primary hits, 80 of which were confirmed by dose response curves with half-maximal inhibitory effects ranging from 0.3 to 60 microM. Among these, six compounds displayed an inhibitory effect on PrP(Sc) propagation in scrapie-infected N2a cells. Four of these candidate drugs share an N'-benzylidene-benzohydrazide core structure. Thus, the combination of high-throughput in vitro assay with the established cell culture system provides a rapid and efficient method to identify new antiprion drugs, which corroborates that interaction of PrP(C) and PrP(Sc) is a crucial molecular step in the propagation of prions. Moreover, SIFT-based screening may facilitate the search for drugs against other diseases linked to protein aggregation.