Hydration effects on the HET-s prion and amyloid-beta fibrillous aggregates, studied with three-dimensional molecular theory of solvation

Reaction Entropies in Solution from Analytical Three-Dimensional Reference Interaction Site Model Derivatives with Application to Redox and Spin-Crossover Processes

An analytical approach to compute the excess entropy of solvation at constant pressure in three-dimensional reference interaction site model (3D-RISM) calculations is presented. It includes the changes in the macroscopic dielectric constant of the solvent upon variation of temperature and density. The approach is exact within the framework of force-field descriptions of the solute and gives reasonable results for self-consistently determined electrostatics as used in the 3D-RISM-self-consistent field approach, particularly for entropy differences. The new method is applied to simple examples of reaction entropies of iron complexes in aqueous solution, for which simple gas-phase calculations and many other approaches give unreliable estimates. For both redox half-reactions and spin-crossover processes, (semi)quantitative agreement with experimental reaction entropies can be achieved out of the box.

Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation

Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.

The impact of water on the ambivalent behavior and paradoxical phenomenon of the amyloid-β fibril protein

The crucial role of water in amyloid-β(Aβ) fibril proteins is evaluated in several ways including the water's thermodynamic and kinetic solvation effects. As regards the water's character, its hindered-rotation barriers are also considered. The following protein molecules considered here are: the Aβ40 (PDB ID: 2LMN), Aβ42 (PDB ID: 5KK3 and 2NAO) and the double-layered Aβ17-42 fibril. We discuss: (i) extracellular Aβ40 and Aβ42 fibril monomers exhibit an ambivalent propensity to transform into a helical form toward the N-term region and a β-strand-like form near the C-terminal; (ii) interfacial water molecules play a crucial role in protein-protein interactions, as molecular dynamics simulations have shown a significant impact on the protein-protein binding; (iii) it is shown that the spontaneous dimerization process of the Aβ42 fibril protein in water occurs via a two-step nucleation-accommodation mechanism; (iv) MD simulations of the double-layered Aβ17-42 fibril model show that the C↔C interface appears more energetically favorable than the N↔N interface due to large hydrophobic contacts; (v) the water's role in the HET-s prion and in the Aβ fibrillar aggregates; (vi) it was found that the monomer-oligomer equilibrium spontaneously dissociates into stable monomeric species when they are incubated up to 3 μm for a longer time (>1 week) in a physiological buffer.

How Water's Properties Are Encoded in Its Molecular Structure and Energies

How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties.

Theoretical approaches for dynamical ordering of biomolecular systems

BACKGROUND

Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies.

SCOPE OF REVIEW

We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified.

MAJOR CONCLUSIONS

The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function.

GENERAL SIGNIFICANCE

These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

A closure relation to molecular theory of solvation for macromolecules

We propose a closure to the integral equations of molecular theory of solvation, particularly suitable for polar and charged macromolecules in electrolyte solution. This includes such systems as oligomeric polyelectrolytes at a finite concentration in aqueous and various non-aqueous solutions, as well as drug-like compounds in solution. The new closure by Kobryn, Gusarov, and Kovalenko (KGK closure) imposes the mean spherical approximation (MSA) almost everywhere in the solvation shell but levels out the density distribution function to zero (with the continuity at joint boundaries) inside the repulsive core and in the spatial regions of strong density depletion emerging due to molecular associative interactions. Similarly to MSA, the KGK closure reduces the problem to a linear equation for the direct correlation function which is predefined analytically on most of the solvation shells and has to be determined numerically on a relatively small (three-dimensional) domain of strong depletion, typically within the repulsive core. The KGK closure leads to the solvation free energy in the form of the Gaussian fluctuation (GF) functional. We first test the performance of the KGK closure coupled to the reference interaction site model (RISM) integral equations on the examples of Lennard-Jones liquids, polar and nonpolar molecular solvents, including water, and aqueous solutions of simple ions. The solvation structure, solvation chemical potential, and compressibility obtained from RISM with the KGK closure favorably compare to the results of the hypernetted chain (HNC) and Kovalenko-Hirata (KH) closures, including their combination with the GF solvation free energy. We then use the KGK closure coupled to RISM to obtain the solvation structure and thermodynamics of oligomeric polyelectrolytes and drug-like compounds at a finite concentration in electrolyte solution, for which no convergence is obtained with other closures. For comparison, we calculate their solvation structure from molecular dynamics (MD) simulations. We further couple the 3D-RISM integral equation with the 3D-version of the KGK closure, and solve it for molecular mixtures as well as oligomeric polyelectrolytes and drug-like molecules in electrolyte solutions.

Small molecule hydration energy and entropy from 3D-RISM

Implicit solvent models offer an attractive way to estimate the effects of a solvent environment on the properties of small or large solutes without the complications of explicit simulations. One common test of accuracy is to compute the free energy of transfer from gas to liquid for a variety of small molecules, since many of these values have been measured. Studies of the temperature dependence of these values (i.e. solvation enthalpies and entropies) can provide additional insights into the performance of implicit solvent models. Here, we show how to compute temperature derivatives of hydration free energies for the 3D-RISM integral equation approach. We have computed hydration free energies of 1123 small drug-like molecules (both neutral and charged). Temperature derivatives were also used to calculate hydration energies and entropies of 74 of these molecules (both neutral and charged) for which experimental data is available. While direct results have rather poor agreement with experiment, we have found that several previously proposed linear hydration free energy correction schemes give good agreement with experiment. These corrections also provide good agreement for hydration energies and entropies though simple extensions are required in some cases.

SAMPL5: 3D-RISM partition coefficient calculations with partial molar volume corrections and solute conformational sampling

Implicit solvent methods for classical molecular modeling are frequently used to provide fast, physics-based hydration free energies of macromolecules. Less commonly considered is the transferability of these methods to other solvents. The Statistical Assessment of Modeling of Proteins and Ligands 5 (SAMPL5) distribution coefficient dataset and the accompanying explicit solvent partition coefficient reference calculations provide a direct test of solvent model transferability. Here we use the 3D reference interaction site model (3D-RISM) statistical-mechanical solvation theory, with a well tested water model and a new united atom cyclohexane model, to calculate partition coefficients for the SAMPL5 dataset. The cyclohexane model performed well in training and testing ([Formula: see text] for amino acid neutral side chain analogues) but only if a parameterized solvation free energy correction was used. In contrast, the same protocol, using single solute conformations, performed poorly on the SAMPL5 dataset, obtaining [Formula: see text] compared to the reference partition coefficients, likely due to the much larger solute sizes. Including solute conformational sampling through molecular dynamics coupled with 3D-RISM (MD/3D-RISM) improved agreement with the reference calculation to [Formula: see text]. Since our initial calculations only considered partition coefficients and not distribution coefficients, solute sampling provided little benefit comparing against experiment, where ionized and tautomer states are more important. Applying a simple [Formula: see text] correction improved agreement with experiment from [Formula: see text] to [Formula: see text], despite a small number of outliers. Better agreement is possible by accounting for tautomers and improving the ionization correction.

Combined effects of solvation and aggregation propensity on the final supramolecular structures adopted by hydrophobic, glycine-rich, elastin-like polypeptides

Previous work on elastin-like polypeptides (ELPs) made of hydrophobic amino acids of the type XxxGlyGlyZzzGly (Xxx, Zzz = Val, Leu) has consistently shown that differing dominant supramolecular structures were formed when the suspending media were varied: helical, amyloid-like fibers when suspended in water and globules evolving into "string of bead" structures, poly(ValGlyGlyValGly), or cigar-like bundles, poly(ValGlyGlyLeuGly), when suspended in methyl alcohol. Comparative experiments with poly(LeuGlyGlyValGly) have further indicated that the interface energy plays a significant role and that solvation effects act in concomitance with the intrinsic aggregation propensity of the repeat sequence. Continuing our investigation on ELPs using surface (X-ray photoelectron spectroscopy, atomic force microscopy) and bulk (circular dichroism, Fourier transform infrared spectroscopy) techniques for their characterization, here we have compared the effect of suspending solvents (H(2)O, dimethylsulfoxide, ethylene glycol, and MeOH) on poly(ValGlyGlyValGly), the polypeptide most inclined to form long and well-refined helical fibers in water, searching for the signature of intermolecular interactions occurring between the polypeptide chains in the given suspension. The influence of sequence specificities has been studied by comparing poly(ValGlyGlyValGly) and poly(LeuGlyGlyValGly) with a similar degree of polymerization. Deposits on substrates of the polypeptides were characterized taking into account the differing evaporation rate of solvents, and tests on their stability in ultra high vacuum were performed. Finally, combining experimental and computational studies, we have revaluated the three-dimensional modeling previously proposed for the supramolecular assembly in water of poly(ValGlyGlyValGly). The results were discussed and rationalized also in the light of published data.

Non-Amyloid-β Component of Human α-Synuclein Oligomers Induces Formation of New Aβ Oligomers: Insight into the Mechanisms That Link Parkinson's and Alzheimer's Diseases

Parkinson's disease (PD) is characterized by the formation of Lewy bodies (LBs), of which their major component is the non-amyloid-β component (NAC) of α-synuclein (AS). Clinical studies have identified a link between PD and Alzheimer's disease (AD), but the question of why PD patients are at risk to develop various types of dementia, such as AD, is still elusive. In vivo studies have shown that Aβ can act as a seed for NAC/AS aggregation, promoting NAC/AS aggregation and thus contributing to the etiology of PD. However, the mechanisms by which NAC/AS oligomers interact with Aβ oligomers are still elusive. This work presents the interactions between NAC oligomers and Aβ oligomers at atomic resolution by applying extensive molecular dynamics simulations for an ensemble of cross-seeded NAC-Aβ(1-42) oligomers. The main conclusions of this study are as follows: first, the cross-seeded NAC-Aβ(1-42) oligomers represent polymorphic states, yet NAC oligomers prefer to interact with Aβ(1-42) oligomers to form double-layer over single-layer conformations due to electrostatic/hydrophobic interactions; second, among the single-layer conformations, the NAC oligomers induce formation of new β-strands in Aβ(1-42) oligomers, thus leading to new Aβ oligomer structures; and third, NAC oligomers stabilize the cross-β structure of Aβ oligomers, i.e., yielding compact Aβ fibril-like structures.

Octanol-Water Partition Coefficient from 3D-RISM-KH Molecular Theory of Solvation with Partial Molar Volume Correction

The octanol-water partition coefficient is an important physical-chemical characteristic widely used to describe hydrophobic/hydrophilic properties of chemical compounds. The partition coefficient is related to the transfer free energy of a compound from water to octanol. Here, we introduce a new protocol for prediction of the partition coefficient based on the statistical-mechanical, 3D-RISM-KH molecular theory of solvation. It was shown recently that with the compound-solvent correlation functions obtained from the 3D-RISM-KH molecular theory of solvation, the free energy functional supplemented with the correction linearly related to the partial molar volume obtained from the Kirkwood-Buff/3D-RISM theory, also called the "universal correction" (UC), provides accurate prediction of the hydration free energy of small compounds, compared to explicit solvent molecular dynamics [ Palmer , D. S. ; J. Phys.: Condens. Matter 2010 , 22 , 492101 ]. Here we report that with the UC reparametrized accordingly this theory also provides an excellent agreement with the experimental data for the solvation free energy in nonpolar solvent (1-octanol) and so accurately predicts the octanol-water partition coefficient. The performance of the Kovalenko-Hirata (KH) and Gaussian fluctuation (GF) functionals of the solvation free energy, with and without UC, is tested on a large library of small compounds with diverse functional groups. The best agreement with the experimental data for octanol-water partition coefficients is obtained with the KH-UC solvation free energy functional.

Negative thermal expansibility change for dissociation of lysozyme variant amyloid protofibril

A disulfide-deficient variant of hen lysozyme, 0SS, is known to form an amyloid protofibril spontaneously, and to dissociate into monomers at high hydrostatic pressure. We carried out native PAGE at various temperatures (20-35°C) and pressures (0.1-200 MPa), to characterize the dissociation equilibrium of disulfide-deficient variant of hen lysozyme amyloid protofibril. Based on the density profiles, the partial molar volume and thermal expansibility changes for dissociation, ΔvD and ΔeD , were obtained to be -74 cm(3) /mol at 25°C and -2.3 cm(3) mol(-1) K(-1) , respectively. The dissociation of amyloid fibril destroys the cross β-structure, and such conformational destruction in native protein fold rarely accompanies negative thermal expansibility change. We discussed the negative thermal expansibility change in terms of hydration and structural packing of the amyloid protofibril.

Role of water in ligand binding to maltose-binding protein: insight from a new docking protocol based on the 3D-RISM-KH molecular theory of solvation

Maltose-binding protein is a periplasmic binding protein responsible for transport of maltooligosaccarides through the periplasmic space of Gram-negative bacteria, as a part of the ABC transport system. The molecular mechanisms of the initial ligand binding and induced large scale motion of the protein's domains still remain elusive. In this study, we use a new docking protocol that combines a recently proposed explicit water placement algorithm based on the 3D-RISM-KH molecular theory of solvation and conventional docking software (AutoDock Vina) to explain the mechanisms of maltotriose binding to the apo-open state of a maltose-binding protein. We confirm the predictions of previous NMR spectroscopic experiments on binding modes of the ligand. We provide the molecular details on the binding mode that was not previously observed in the X-ray experiments. We show that this mode, which is defined by the fine balance between the protein-ligand direct interactions and solvation effects, can trigger the protein's domain motion resulting in the holo-closed structure of the maltose-binding protein with the maltotriose ligand in excellent agreement with the experimental data. We also discuss the role of water in blocking unfavorable binding sites and water-mediated interactions contributing to the stability of observable binding modes of maltotriose.

Multiple time step molecular dynamics in the optimized isokinetic ensemble steered with the molecular theory of solvation: accelerating with advanced extrapolation of effective solvation forces

We develop efficient handling of solvation forces in the multiscale method of multiple time step molecular dynamics (MTS-MD) of a biomolecule steered by the solvation free energy (effective solvation forces) obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model complemented with the Kovalenko-Hirata closure approximation). To reduce the computational expenses, we calculate the effective solvation forces acting on the biomolecule by using advanced solvation force extrapolation (ASFE) at inner time steps while converging the 3D-RISM-KH integral equations only at large outer time steps. The idea of ASFE consists in developing a discrete non-Eckart rotational transformation of atomic coordinates that minimizes the distances between the atomic positions of the biomolecule at different time moments. The effective solvation forces for the biomolecule in a current conformation at an inner time step are then extrapolated in the transformed subspace of those at outer time steps by using a modified least square fit approach applied to a relatively small number of the best force-coordinate pairs. The latter are selected from an extended set collecting the effective solvation forces obtained from 3D-RISM-KH at outer time steps over a broad time interval. The MTS-MD integration with effective solvation forces obtained by converging 3D-RISM-KH at outer time steps and applying ASFE at inner time steps is stabilized by employing the optimized isokinetic Nosé-Hoover chain (OIN) ensemble. Compared to the previous extrapolation schemes used in combination with the Langevin thermostat, the ASFE approach substantially improves the accuracy of evaluation of effective solvation forces and in combination with the OIN thermostat enables a dramatic increase of outer time steps. We demonstrate on a fully flexible model of alanine dipeptide in aqueous solution that the MTS-MD/OIN/ASFE/3D-RISM-KH multiscale method of molecular dynamics steered by effective solvation forces allows huge outer time steps up to tens of picoseconds without affecting the equilibrium and conformational properties, and thus provides a 100- to 500-fold effective speedup in comparison to conventional MD with explicit solvent. With the statistical-mechanical 3D-RISM-KH account for effective solvation forces, the method provides efficient sampling of biomolecular processes with slow and/or rare solvation events such as conformational transitions of hydrated alanine dipeptide with the mean life times ranging from 30 ps up to 10 ns for "flip-flop" conformations, and is particularly beneficial for biomolecular systems with exchange and localization of solvent and ions, ligand binding, and molecular recognition.

Multiscale modeling of solvation in chemical and biological nanosystems and in nanoporous materials

Statistical–mechanical, 3D-RISM-KH molecular theory of solvation (3D reference interaction site model with the Kovalenko–Hirata closure) is promising as an essential part of multiscale methodology for chemical and biomolecular nanosystems in solution. 3D-RISM-KH explains the molecular mechanisms of self-assembly and conformational stability of synthetic organic rosette nanotubes (RNTs), aggregation of prion proteins and β-sheet amyloid oligomers, protein-ligand binding, and function-related solvation properties of complexes as large as the Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) and GroEL/ES chaperone. Molecular mechanics/Poisson–Boltzmann (generalized Born) surface area [MM/PB(GB)SA] post-processing of molecular dynamics (MD) trajectories involving SA empirical nonpolar terms is replaced with MM/3D-RISM-KH statistical–mechanical evaluation of the solvation thermodynamics. 3D-RISM-KH has been coupled with multiple time-step (MTS) MD of the solute biomolecule driven by effective solvation forces, which are obtained analytically by converging the 3D-RISM-KH integral equations at outer time-steps and are calculated in between by using solvation force coordinate extrapolation (SFCE) in the subspace of previous solutions to 3D-RISM-KH. The procedure is stabilized by the optimized isokinetic Nosé–Hoover (OIN) chain thermostatting, which enables gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties. The multiscale OIN/SFCE/3D-RISM-KH algorithm is implemented in the Amber package and illustrated on a fully flexible model of alanine dipeptide in aqueous solution, exhibiting the computational rate of solvent sampling 20 times faster than standard MD with explicit solvent. Further substantial acceleration can be achieved with 3D-RISM-KH efficiently sampling essential events with rare statistics such as exchange and localization of solvent, ions, and ligands at binding sites and pockets of the biomolecule. 3D-RISM-KH was coupled with ab initio complete active space self-consistent field (CASSCF) and orbital-free embedding (OFE) Kohn–Sham (KS) density functional theory (DFT) quantum chemistry methods in an SCF description of electronic structure, optimized geometry, and chemical reactions in solution. The (OFE)KS-DFT/3D-RISM-KH multi-scale method is implemented in the Amsterdam Density Functional (ADF) package and extensively validated against experiment for solvation thermochemistry, photochemistry, conformational equilibria, and activation barriers of various nanosystems in solvents and ionic liquids (ILs). Finally, the replica RISM-KH-VM molecular theory for the solvation structure, thermodynamics, and electrochemistry of electrolyte solutions sorbed in nanoporous materials reveals the molecular mechanisms of sorption and supercapacitance in nanoporous carbon electrodes, which is drastically different from a planar electrical double layer.

Aβ(39-42) modulates Aβ oligomerization but not fibril formation

Recently, certain C-terminal fragments (CTFs) of Aβ42 have been shown to be effective inhibitors of Aβ42 toxicity. Here, we examine the interactions between the shortest CTF in the original series, Aβ(39-42), and full-length Aβ. Mass spectrometry results indicate that Aβ(39-42) binds directly to Aβ monomers and to the n = 2, 4, and 6 oligomers. The Aβ42:Aβ(39-42) complex is further probed using molecular dynamics simulations. Although the CTF was expected to bind to the hydrophobic C-terminus of Aβ42, the simulations show that Aβ(39-42) binds at several locations on Aβ42, including the C-terminus, other hydrophobic regions, and preferentially in the N-terminus. Ion mobility-mass spectrometry (IM-MS) and electron microscopy experiments indicate that Aβ(39-42) disrupts the early assembly of full-length Aβ. Specifically, the ion-mobility results show that Aβ(39-42) prevents the formation of large decamer/dodecamer Aβ42 species and, moreover, can remove these structures from solution. At the same time, thioflavin T fluorescence and electron microscopy results show that the CTF does not inhibit fibril formation, lending strong support to the hypothesis that oligomers and not amyloid fibrils are the Aβ form responsible for toxicity. The results emphasize the role of small, soluble assemblies in Aβ-induced toxicity and suggest that Aβ(39-42) inhibits Aβ-induced toxicity by a unique mechanism, modulating early assembly into nontoxic hetero-oligomers, without preventing fibril formation.

A peptide's perspective of water dynamics

This Perspective is focused on amide groups of peptides interacting with water. The 2D IR spectroscopy has already enabled structural aspects of the peptide backbone to be determined through its ability to measure the coupling between different amide-I modes. Here we describe why nonlinear IR is emerging as the method of choice to examine the fast components of the water dynamics near peptides and how isotopically edited peptide links can be used to probe the local water at a residue level in proteins. This type of research necessarily involves an intimate mix of theory and experiment. The description of the results is underpinned by relatively well established quantum-statistical theories that describe the important manifestations of peptide vibrational frequency fluctuations.

Influence of amino acid specificities on the molecular and supramolecular organization of glycine-rich elastin-like polypeptides in water

Elastin-like polypeptides adopt complex supramolecular structures, showing either a hydrophobic or a hydrophilic surface, depending on their surrounding environment and the supporting substrate. The preferred organization is important in many situations ranging from biocompatibility to bio-function. Here we compare the n-repeat pentamer LeuGlyGlyValGly (n = 7) with the analogue ValGlyGlyValGly (n = 5), as water suspensions and as deposits on silicon substrates. These sequences contain the repeat XxxGlyGlyZzzGly (Xxx, Zzz = Val, Leu) motif belonging to the hydrophobic glycine-rich domain of elastin and represent a simplified model from which to obtain information on molecular interactions functional to elastin itself. The compounds studied differ only by the presence of the -CH(2)- spacer in the Leu moiety and thus the work was aimed at revealing the influence of this spacer element on self assembly. Both polypeptides were studied under identical conditions, using combined techniques, to identify differences in their conformational states both at molecular (CD, FTIR) and supramolecular (XPS, AFM) levels. By these means, together with a Congo Red spectroscopic assay of β-sheet formation in water, a clear correlation between amino acid sequences (sequence specificity) and their kinetics and ordering of aggregation has emerged. The novel outcomes of this work are from the supplementary measurements, made to augment the AFM and XPS studies, showing that the significant step in the self assembly of both polypeptides takes place in the liquid phase and from the finding that the substitution of Val by Leu in the first position of the pentapeptide effectively inhibits the formation of amyloidal fibers.

Comparative studies of nontoxic and toxic amyloids interacting with membrane models at the air-water interface

Many in vitro studies have pointed out the interaction between amyloids and membranes, and their potential involvement in amyloid toxicity. In a previous study, we generated a yeast toxic mutant (M8) of the harmless model amyloid protein HET-s((218-289)). In this study, we compared the self-assembling process of the nontoxic wild-type (WT) and toxic (M8) protein at the air-water interface and in interaction with various phospholipid monolayers (DOPE, DOPC, DOPI, DOPS and DOPG). We first demonstrate using ellipsometry measurements and polarization-modulated infrared reflection absorption spectroscopy (PMIRRAS) that the air-water interface promotes and modifies the assembly of WT since an amyloid-like film was instantaneously formed at the interface with an antiparallel β-sheet structuration instead of the parallel β-sheet commonly observed for amyloid fibers generated in solution. The toxic mutant (M8) behaves in a similar manner at the air-water interface or in bulk, with a fast self-assembling and an antiparallel β-sheet organization. The transmission electron microscopy (TEM) images established the fibrillous morphology of the protein films formed at the air-water interface. Second, we demonstrate for the first time that the main driving force between this particular fungus amyloid and membrane interaction is based on electrostatic interactions with negatively charged phospholipids (DOPG, DOPI, DOPS). Interestingly, the toxic mutant (M8) clearly induces perturbations of the negatively charged phospholipid monolayers, leading to a massive surface aggregation, whereas the nontoxic (WT) exhibits a slight effect on the membrane models. This study allows concluding that the toxicity of the M8 mutant could be due to its high propensity to interact with membranes.

An MM/3D-RISM approach for ligand binding affinities

We have modified the popular MM/PBSA or MM/GBSA approaches (molecular mechanics for a biomolecule, combined with a Poisson-Boltzmann or generalized Born electrostatic and surface area nonelectrostatic solvation energy) by employing instead the statistical-mechanical, three-dimensional molecular theory of solvation (also known as 3D reference interaction site model, or 3D-RISM-KH) coupled with molecular mechanics or molecular dynamics ( Blinov , N. ; et al. Biophys. J. 2010 ; Luchko , T. ; et al. J. Chem. Theory Comput. 2010 ). Unlike the PBSA or GBSA semiempirical approaches, the 3D-RISM-KH theory yields a full molecular picture of the solvation structure and thermodynamics from the first principles, with proper account of chemical specificities of both solvent and biomolecules, such as hydrogen bonding, hydrophobic interactions, salt bridges, etc. We test the method on the binding of seven biotin analogues to avidin in aqueous solution and show it to work well in predicting the ligand-binding affinities. We have compared the results of 3D-RISM-KH with four different generalized Born and two Poisson-Boltzmann methods. They give absolute binding energies that differ by up to 208 kJ/mol and mean absolute deviations in the relative affinities of 10-43 kJ/mol.

Molecular modeling of two distinct triangular oligomers in amyloid beta-protein

Amyloid-beta (Abeta) peptides exhibit many distinct structural morphology at the early aggregate stage, some of which are biological relevant to the pathogenesis of Alzheimer's disease (AD). Atomic-resolution structures of the early Abeta aggregates and their conformational changes in amyloid aggregation remain elusive. Here, we perform all-atom molecular modeling and dynamics simulations to obtain two stable triangular-like Abeta structures with the lowest packing energy, one corresponding to the Tycko's model (Paravastu, A.; Leapman, R.; Yau, W.; Tycko, R. Proc. Nat. Acad. Soc. U.S.A. 2008, 105, 18349-18354) (referred to C-WT model) and the other corresponding to computational model (N-WT model). Both models have the same 3-fold symmetry but distinct beta-sheet organizations in which three Abeta hexamers pack together via either C-terminal beta-strand residues or N-terminal beta-strand residues forming distinct hydrophobic cross section. Structural and energetic comparisons of two 3-fold Abeta oligomers, coupled with structural changes upon the mutations occurring at the interacting interfaces, reveal that although hydrophobic interactions are still dominant forces, electrostatic interactions are more favorable in the N-WT model due to the formation of more and stable intersheet salt bridges, while solvation energy is more favorable in the C-WT model due to more exposed hydrophilic residues to solvent. Both models display many common features similar to other amyloid oligomers and therefore are likely to be biologically relevant.

Dynamics and stability of amyloid-like steric zipper assemblies with hydrophobic dry interfaces

Recent seminal investigations have suggested that the basic structural motif of amyloid fibers may be constituted by a tight association of two facing beta-sheets (steric zipper). Although this model has been derived from crystal structures of small peptide models, several theoretical investigations, essentially focused on steric zipper interface containing large polar and/or aromatic side chains, have confirmed the stability of this motif in a crystal-free context. To analyze the general validity of these findings, we carried out molecular dynamics (MD) simulations on aggregates stabilized by steric zipper interfaces made also of small or hydrophobic residues. In particular, we here characterized assemblies formed by the peptides SSTSAA and VQIVYK, whose structures have been recently solved at high resolution. In contrast to previous results obtained for polar/aromatic aggregates of the same size and with similar interface area, steric zipper assemblies composed of a pair of 10-stranded beta-sheets show high fluctuations and significant distortions in the simulation timescales (40-60 ns). Taking into account the crystal packing, the effect of the addition of an extra sheet to the assemblies was also evaluated. The MD results indicate that this addition does not provide extra-stabilization to the pair of sheet models. Although present data do not preclude the possibility that the steric zipper association identified in the crystal structure is the basic motif of SSTSAA and VQIVYK fibers, our findings highlight the importance of the nature of residues directly involved in the motif. Indeed, polar and aromatic residues that may form intrasheet and intersheet interactions likely provide a strong contribution to the steric zipper motif stability. Along this line, assemblies endowed with hydrophobic residues presumably require larger interfaces. In line with this suggestion, MD analysis of the HET-s(218-289) prion models composed of a similar number of strands shows that the assembly is endowed with a remarkable stability.