Microtubule stability studied by three-dimensional molecular theory of solvation

We study microtubular supramolecular architectures of tubulin dimers self-assembling into linear protofilaments, in turn forming a closed tube, which is an important component of the cytoskeleton. We identify the protofilament arrangements with the lowest free energy using molecular dynamics to optimize tubulin conformations. We then use the three-dimensional molecular theory of solvation to obtain the hydration structure of protofilaments built of optimized tubulins and the solvent-mediated effective potential between them. The latter theoretical method, based on first principles of statistical mechanics, is capable of predicting the structure and thermodynamics of solvation of supramolecular architectures. We obtained a set of profiles of the potential of mean force between protofilaments in a periodic two-dimensional sheet in aqueous solution. The profiles were calculated for a number of amino acid sequences, tubulin conformations, and spatial arrangements of protofilaments. The results indicate that the effective interaction between protofilaments in aqueous solution depends little on the isotypes studied; however, it strongly depends on the M loop conformation of beta-tubulin. Based on the analysis of the potential of mean force between adjacent protofilaments, we found the optimal arrangement of protofilaments, which is in good agreement with other studies. We also decomposed the potential of mean force into its energetic and entropic components, and found that both are considerable in the free-energy balance for the stabilized protofilament arrangements.     

The HET-s prion protein of the filamentous fungus Podospora anserina aggregates in vitro into amyloid-like fibrils.

The HET-s protein of Podospora anserina is a fungal prion. This protein behaves as an infectious cytoplasmic element that is transmitted horizontally from one strain to another. Under the prion form, the HET-s protein forms aggregates in vivo. The specificity of this prion model compared with the yeast prions resides in the fact that under the prion form HET-s causes a growth inhibition and cell death reaction when co-expressed with the HET-S protein from which it differs by 13 residues. Herein we describe the purification and initial characterization of recombinant HET-s protein expressed in Escherichia coli. The HET-s protein self-associates over time into high molecular weight aggregates. These aggregates greatly accelerate precipitation of the soluble form. HET-s aggregates appear as amyloid-like fibrils using electron microscopy. They bind Congo Red and show birefringence under polarized light. In the aggregated form, a HET-s fragment of approximately 7 kDa is resistant to proteinase K digestion. CD and FTIR analyses indicate that upon transition to the aggregated state, the HET-s protein undergoes a structural rearrangement characterized by an increase in antiparallel beta-sheet structure content. These results suggest that the [Het-s] prion element propagates in vivo as an infectious amyloid.     

Prion and non-prion amyloids of the HET-s prion forming domain

HET-s is a prion protein of the fungus Podospora anserina. A plausible structural model for the infectious amyloid fold of the HET-s prion-forming domain, HET-s(218-289), makes it an attractive system to study structure-function relationships in amyloid assembly and prion propagation. Here, we report on the diversity of HET-s(218-289) amyloids formed in vitro. We distinguish two types formed at pH 7 from fibrils formed at pH 2, on morphological grounds. Unlike pH 7 fibrils, the pH 2 fibrils show very little if any prion infectivity. They also differ in ThT-binding, resistance to denaturants, assembly kinetics, secondary structure, and intrinsic fluorescence. Both contain 5 nm fibrils, either bundled or disordered (pH 7) or as tightly twisted protofibrils (pH 2). We show that electrostatic interactions are critical for the formation and stability of the infectious prion fold given in the current model. The altered properties of the amyloid assembled at pH 2 may arise from a perturbation in the subunit fold or fibrillar stacking.     

Hydration structure,thermodynamics, and functions of protein studied by the 3D-RISM theory

The three-dimensional reference interaction site model (3D-RISM) theory is a molecular theory of solvation, which has been developed recently. Unlike the conventional statistical-mechanical theories of liquids, the applications of which are limited to rather simple systems, the new theory has been demonstrating great success in predicting the hydration structure and thermodynamic quantities of protein. Here, we review the recent applications of the 3D-RISM theory to some subjects concerning protein hydration: detecting water molecules in protein cavities, the pressure effects on protein conformation in terms of the partial molar volume, and the pressure reversal of anesthesia. These results demonstrate the outstanding capability of the theory for investigating various issues in biochemistry and biophysics.     

Correlations among morphology, beta-sheet stability, and molecular structure in prion peptide aggregates

The misfolding of proteins into beta-sheets and the subsequent aggregation of these sheets into fibrous networks underlies many diseases. In this paper, the role of peptide structure in determining the ordering of beta-sheet aggregates and the morphology of fibrils and protofibrils is dissected. Using a series of peptides based on residues 109-122 of the Syrian hamster prion protein (H1) with a range of substitutions at position 117, the link between side chain interactions and beta-sheet thermal stability has been investigated. The thermal stability of beta-sheets is associated with the peptides' ability to adopt the same alignment as wild-type H1, with residue 117 in register across all beta-strands [Silva, R. A. G. D., Barber-Armstrong, W., and Decatur, S. M. (2003) J. Am. Chem. Soc. 125, 13674-13675]. These aligned strands are capable of forming long, rigid, and twisted fibrils (as visualized by atomic force microscopy) which are thermostable. Peptides which do not adopt this strand alignment aggregate to form thin, flexible, and smooth protofibrils. The ability to form ordered aggregates, and thus to form twisted fibrils, is modulated by the structure of the side chain of residue 117.     

Locating missing water molecules in protein cavities by the three-dimensional reference interaction site model theory of molecular solvation

Water molecules confined in protein cavities are of great importance in understanding the protein structure and functions. However, it is a nontrivial task to locate such water molecules in protein by the ordinary molecular simulation and modeling techniques as well as experimental methods. The present study proves that the three-dimensional reference interaction site model (3D-RISM) theory, a recently developed statistical-mechanical theory of molecular solvation, has an outstanding advantage in locating such water molecules. In this paper, we demonstrate that the 3D-RISM theory is able to reproduce the structure and the number of water molecules in cavities of hen egg-white lysozyme observed commonly in the X-ray structures of different resolutions and conditions. Furthermore, we show that the theory successfully identified a water molecule in a cavity, the existence of which has been ambiguous even from the X-ray results. In contrast, we confirmed that molecular dynamics simulation is helpless at present to find such water molecules because the results substantially depend on the initial coordinates of water molecules. Possible applications of the theory to problems in the fields of biochemistry and biophysics are also discussed.     

Domain organization and structure-function relationship of the HET-s prion protein of Podospora anserina

The [Het-s] infectious element of the fungus Podospora anserina is a prion protein involved in a genetically controlled cell death reaction termed heterokaryon incompatibility. Previous analyses indicate that [Het-s] propagates as a self-perpetuating amyloid aggregate. The HET-s protein is 289 amino acids in length. Herein, we identify the region of the HET-s protein that is responsible for amyloid formation and prion propagation. The region of HET-s spanning residues 218-289 forms amyloid fibers in vitro and allows prion propagation in vivo. Conversely, a C-terminal deletion in HET-s prevents amyloid aggregation in vitro and prion propagation in vivo, and abolishes the incompatibility function. In the soluble form of HET-s, the region from residue 1 to 227 forms a well-folded domain while the C-terminal region is highly flexible. Together, our data establish a domain structure-function relationship for HET-s amyloid formation, prion propagation and incompatibility activity.     

Non-infectious aggregates of the prion protein react with several PrPSc-directed antibodies

The key event in the pathogenesis of prion diseases is the conformational conversion of the normal prion protein (PrP) (PrP^C) into an infectious, aggregated isoform (PrP^Sc) that has a high content of β-sheet. Historically, a great deal of effort has been devoted to developing antibodies that specifically recognize PrP^Sc but not PrP^C, as such antibodies would have enormous diagnostic and experimental value. A mouse monoclonal IgM antibody (designated 15B3) and three PrP motif-grafted monoclonal antibodies (referred to as IgG 19–33, 89–112, and 136–158) have been previously reported to react specifically with infectious PrP^Sc but not PrP^C. In this study, we extend the characterization of these four antibodies by testing their ability to immunoprecipitate and immunostain infectious and non-infectious aggregates of wild-type, mutant, and recombinant PrP. We find that 15B3 as well as the motif-grafted antibodies recognize multiple types of aggregated PrP, both infectious and non-infectious, including forms found in brain, in transfected cells, and induced in vitro from purified recombinant protein. These antibodies are exquisitely selective for aggregated PrP, and do not react with soluble PrP even when present in vast excess. Our results suggest that 15B3 and the motif-grafted antibodies recognize structural features common to both infectious and non-infectious aggregates of PrP. Our study extends the utility of these antibodies for diagnostic and experimental purposes, and it provides new insight into the structural changes that accompany PrP oligomerization and prion propagation.     

Free-energy analysis of protein solvation with all-atom molecular dynamics simulation combined with a theory of solutions

1. Download : Download high-res image (172KB)
2. Download : Download full-size image

     

Recognition complex between the HMG domain of LEF-1 and its cognate DNA studied by molecular dynamics simulations with explicit solvation

Molecular dynamics simulations of the complex formed between the HMG box of the lymphoid enhancer-binding factor (LEF-1) and its cognate DNA duplex were carried out with explicit inclusion of water. The simulation started with an NMR-based model (pdb code 2LEF) and the dynamics was pursued for 10 nanoseconds without constraints. It revealed that water intervenes in many ionic/polar interactions, establishing in particular local equilibria between direct and water-mediated hydrogen bonds, and thus increasing the entropy of the complex. Quite unexpectedly, the simulation indicated that a binding pocket for a specific water molecule may be reversibly formed at the apex of the bend induced in the DNA helix by LEF-1 binding, where a methionine side chain intercalates between two destacked adenines. We observed that the specific water molecule can temporarily replace the intercalated S-CH(3) group, acting as a sort of "extension" of the side chain. The residence time of this water molecule was about 3.5 ns. Simulations of the cognate DNA alone showed that this sequence has no intrinsic tendency to bend; therefore, the bending occurs solely as a consequence of the recognition, following the "induced-fit" mechanism.     

Superoxide dismutase: fluctuations in the structure and solvation of the active site channel studied by molecular dynamics simulation

The molecular dynamics (MD) simulation of superoxide dismutase (SOD) in water is carried out for a total of 23 ps. The simulation system is a 26 A sphere centered at the active site of SOD, including 1602 atoms from SOD and 1761 water molecules. There is no gross deviation from the x-ray structure for the average MD structure. The structure and potential fluctuations around the active site are examined. The results provide new insight to the interactions between SOD and its substrate superoxide.     

Application of a universal solvation model to nucleic acid bases: comparison of semiempirical molecular orbital theory, ab initio Hartree-Fock theory, and density functional theory

The free energies of solvation of six nucleic acid bases (adenine, cytosine, hypoxanthine, guanine, thymine, and uracil) in water and chloroform are calculated using CM2 class IV charges and SM5.42R atomic surface tensions. Using any of three approximations to the electronic wave function (AM1, Hartree-Fock, or DFT), we obtain good agreement with experiment for five cases where the experimental results are known for the partition coefficients between the two solvents. Decomposition of the solvation effects into bulk electrostatic contributions and first-solvation-shell effects shows that the partitioning is dominated by the former, and this illustrates the importance of using accurate partial atomic charges for modeling these molecules in aqueous solution.     

Role of lipid rafts in the processing of the pathogenic prion and Alzheimer's amyloid-beta proteins

The conformational conversion of the cellular form of the prion protein (PrP C) into the infectious form (PrP Sc) and the proteolytic processing of the amyloid-beta (Abeta) peptide are central pathogenetic events in the prion diseases and Alzheimer's disease, respectively. Cholesterol- and sphingolipid-rich lipid rafts have emerged as important sites for the conversion of PrP C into PrP Sc, and for the proteolytic production, degradation and aggregation of Abeta. Here, we discuss these findings and their implications for our understanding of these disease processes. In addition, the potential for rafts as sites for therapeutic intervention in prion diseases and Alzheimer's disease is considered.     

Hydration and packing effects on prion folding and beta-sheet conversion. High pressure spectroscopy and pressure perturbation calorimetry studies

The main hypothesis for prion diseases proposes that the cellular protein (PrP(C)) can be altered into a misfolded, beta-sheet-rich isoform (PrP(Sc)), which undergoes aggregation and triggers the onset of transmissible spongiform encephalopathies. Here, we compare the stability against pressure and the thermomechanical properties of the alpha-helical and beta-sheet conformations of recombinant murine prion protein, designated as alpha-rPrP and beta-rPrP, respectively. High temperature induces aggregates and a large gain in intermolecular antiparallel beta-sheet (beta-rPrP), a conformation that shares structural similarity with PrP(Sc). alpha-rPrP is highly stable, and only pressures above 5 kilobars (1 kilobar = 100 MegaPascals) cause reversible denaturation, a process that leads to a random and turnrich conformation with concomitant loss of alpha-helix, as measured by Fourier transform infrared spectroscopy. In contrast, aggregates of beta-rPrP are very sensitive to pressure, undergoing transition into a dissociated species that differs from the denatured form derived from alpha-rPrP. The higher susceptibility to pressure of beta-rPrP can be explained by its less hydrated structure. Pressure perturbation calorimetry supports the view that the accessible surface area of alpha-rPrP is much higher than that of beta-rPrP, which explains the lower degree of hydration of beta-rPrP. Our findings shed new light on the mechanism of prion conversion and show how water plays a prominent role. Our results allow us to propose a volume and free energy diagram of the different species involved in the conversion and aggregation. The existence of different folded conformations as well as different denatured states of PrP may explain the elusive character of its conversion into a pathogenic form.     

Heat shock proteins, molecular chaperones and the prion encephalopathies

No Abstract Available     

Studies of base pair sequence effects on DNA solvation based on all-atom molecular dynamics simulations

Detailed analyses of the sequence-dependent solvation and ion atmosphere of DNA are presented based on molecular dynamics (MD) simulations on all the 136 unique tetranucleotide steps obtained by the ABC consortium using the AMBER suite of programs. Significant sequence effects on solvation and ion localization were observed in these simulations. The results were compared to essentially all known experimental data on the subject. Proximity analysis was employed to highlight the sequence dependent differences in solvation and ion localization properties in the grooves of DNA. Comparison of the MD-calculated DNA structure with canonical A- and B-forms supports the idea that the G/C-rich sequences are closer to canonical A- than B-form structures, while the reverse is true for the poly A sequences, with the exception of the alternating ATAT sequence. Analysis of hydration density maps reveals that the flexibility of solute molecule has a significant effect on the nature of observed hydration. Energetic analysis of solute-solvent interactions based on proximity analysis of solvent reveals that the GC or CG base pairs interact more strongly with water molecules in the minor groove of DNA that the AT or TA base pairs, while the interactions of the AT or TA pairs in the major groove are stronger than those of the GC or CG pairs. Computation of solvent-accessible surface area of the nucleotide units in the simulated trajectories reveals that the similarity with results derived from analysis of a database of crystallographic structures is excellent. The MD trajectories tend to follow Manning's counterion condensation theory, presenting a region of condensed counterions within a radius of about 17 A from the DNA surface independent of sequence. The GC and CG pairs tend to associate with cations in the major groove of the DNA structure to a greater extent than the AT and TA pairs. Cation association is more frequent in the minor groove of AT than the GC pairs. In general, the observed water and ion atmosphere around the DNA sequences is the MD simulation is in good agreement with experimental observations.     

Designing drugs to stop the formation of prion aggregates and other amyloids

Amyloid protein aggregates are implicated in many neurodegenerative diseases, including Alzheimer's disease and the prion diseases. Therapeutics to block amyloid formation are often tested in vitro, but it is not clear how to extrapolate from these experiments to a clinical setting, where the effective drug dose may be much lower. Here we address this question using a theoretical kinetic model to calculate the growth rate of protein aggregates as a function of the dose of each of three categories of drug. We find that therapeutics which block the growing ends of amyloids are the most promising, as alternative strategies may be ineffective or even accelerate amyloid formation at low drug concentrations. Our mathematical model can be used to identify and optimise an end-blocking drug in vitro. Our model also suggests an alternative explanation for data previously thought to prove the existence of an entity known as protein X.     

Solution three-dimensional structure of surfactin: A cyclic lipopeptide studied by 1H-nmr,distance geometry,and molecular dynamics

The solution three-dimensional structure of the protonated [Leu7]-surfactin, an hepta-peptide extracted from Bacillus subtilis, has been determined from two-dimensional ¹Hnmr performed in ²H₆-dimethylsulfoxide and combined with molecular modeling. Experimental data included 9 coupling constants, 61 nuclear Overhauser effect derived distances, NH temperature coefficients, and ¹³C relaxation times. Two distance geometry (DISMAN) protocols converged toward models of the structure and the best of them were refined by restrained and unrestrained molecular dynamics (GROMOS). Two structures in accord with the set of experimental constraints are presented. Both are characterized by a “horse saddle” topology for ring atoms on which are attached the two polar Glu and Asp side chains showing an orientation clearly opposite to that of the C_11–13 aliphatic chain. Amphipathic and surface properties of surfactin are certainly related to the existence of such minor polar and a major hydrophobic domains. The particular “claw” configuration of acidic residues observed in surfactin gives important clues for the understanding of its cation binding and transporting ability. © 1994 John Wiley & Sons, Inc.