Bayesian Energy Landscape Tilting: Towards Concordant Models of Molecular Ensembles

Predicting biological structure has remained challenging for systems such as disordered proteins that take on myriad conformations. Hybrid simulation/experiment strategies have been undermined by difficulties in evaluating errors from computational model inaccuracies and data uncertainties. Building on recent proposals from maximum entropy theory and nonequilibrium thermodynamics, we address these issues through a Bayesian energy landscape tilting (BELT) scheme for computing Bayesian hyperensembles over conformational ensembles. BELT uses Markov chain Monte Carlo to directly sample maximum-entropy conformational ensembles consistent with a set of input experimental observables. To test this framework, we apply BELT to model trialanine, starting from disagreeing simulations with the force fields ff96, ff99, ff99sbnmr-ildn, CHARMM27, and OPLS-AA. BELT incorporation of limited chemical shift and ³J measurements gives convergent values of the peptide’s α, β, and PP_II conformational populations in all cases. As a test of predictive power, all five BELT hyperensembles recover set-aside measurements not used in the fitting and report accurate errors, even when starting from highly inaccurate simulations. BELT’s principled framework thus enables practical predictions for complex biomolecular systems from discordant simulations and sparse data.     

The Conformational Ensembles of α-Synuclein and Tau: Combining Single-Molecule FRET and Simulations

Intrinsically disordered proteins (IDPs) are increasingly recognized for their important roles in a range of biological contexts, both in normal physiological function and in a variety of devastating human diseases. However, their structural characterization by traditional biophysical methods, for the purposes of understanding their function and dysfunction, has proved challenging. Here, we investigate the model IDPs α-Synuclein (αS) and tau, that are involved in major neurodegenerative conditions including Parkinson’s and Alzheimer’s diseases, using excluded volume Monte Carlo simulations constrained by pairwise distance distributions from single-molecule fluorescence measurements. Using this, to our knowledge, novel approach we find that a relatively small number of intermolecular distance constraints are sufficient to accurately determine the dimensions and polymer conformational statistics of αS and tau in solution. Moreover, this method can detect local changes in αS and tau conformations that correlate with enhanced aggregation. Constrained Monte Carlo simulations produce ensembles that are in excellent agreement both with experimental measurements on αS and tau and with all-atom, explicit solvent molecular dynamics simulations of αS, with much lower configurational sampling requirements and computational expense.Abbreviations used: AAMD, all-atom molecular dynamics; ECMC, experimentally constrained Monte Carlo; ET_eff, energy transfer efficiency; (sm)FRET, (single molecule) Förster resonance energy transfer; IDP, intrinsically disordered protein; LJ, Lennard-Jones; MC, Monte Carlo; MD, molecular dynamics; PRE, paramagnetic relaxation enhancement; SAX(N)S, small-angle x-ray (neutron) scattering; UMC, unconstrained Monte Carlo     

Flexibility of the Bacterial Chaperone Trigger Factor in Microsecond-Timescale Molecular Dynamics Simulations

The bacterial chaperone trigger factor (TF) is the first chaperone to be encountered by a nascent protein chain as it emerges from the ribosome exit tunnel. Experimental results suggest that TF possesses considerable conformational flexibility, and in an attempt to provide an atomic-level view of this flexibility, we have performed independent 1.5-μs molecular dynamics simulations of TF in explicit solvent using two different simulation force fields (OPLS-AA/L and AMBER ff99SB-ILDN). Both simulations indicate that TF possesses tremendous flexibility, with huge excursions from the crystallographic conformation caused by reorientations of the protein’s constituent domains; both simulations also predict the formation of extensive contacts between TF’s PPIase domain and the Arm 1 domain that is involved in nascent-chain binding. In the OPLS simulation, however, TF rapidly settles into a very compact conformation that persists for at least 1 μs, whereas in the AMBER simulation, it remains highly dynamic; additional simulations in which the two force fields were swapped suggest that these differences are at least partly attributable to sampling issues. The simulation results provide potential rationalizations of a number of experimental observations regarding TF’s conformational behavior and have implications for using simulations to model TF’s function on translating ribosomes.     

Mechanical Resistance in Unstructured Proteins

Single-molecule pulling experiments on unstructured proteins linked to neurodegenerative diseases have measured rupture forces comparable to those for stable folded proteins. To investigate the structural mechanisms of this unexpected force resistance, we perform pulling simulations of the amyloid β-peptide (Aβ) and α-synuclein (αS), starting from simulated conformational ensembles for the free monomers. For both proteins, the simulations yield a set of rupture events that agree well with the experimental data. By analyzing the conformations occurring shortly before rupture in each event, we find that the mechanically resistant structures share a common architecture, with similarities to the folds adopted by Aβ and αS in amyloid fibrils. The disease-linked Arctic mutation of Aβ is found to increase the occurrence of highly force-resistant structures. Our study suggests that the high rupture forces observed in Aβ and αS pulling experiments are caused by structures that might have a key role in amyloid formation.     

Differences in β-strand Populations of Monomeric Aβ40 and Aβ42

Using homonuclear ¹H NOESY spectra, with chemical shifts, ³J_H^N_H^α scalar couplings, residual dipolar couplings, and ¹H-¹⁵N NOEs, we have optimized and validated the conformational ensembles of the amyloid-β 1–40 (Aβ40) and amyloid-β 1–42 (Aβ42) peptides generated by molecular dynamics simulations. We find that both peptides have a diverse set of secondary structure elements including turns, helices, and antiparallel and parallel β-strands. The most significant difference in the structural ensembles of the two peptides is the type of β-hairpins and β-strands they populate. We find that Aβ42 forms a major antiparallel β-hairpin involving the central hydrophobic cluster residues (16–21) with residues 29–36, compatible with known amyloid fibril forming regions, whereas Aβ40 forms an alternative but less populated antiparallel β-hairpin between the central hydrophobic cluster and residues 9–13, that sometimes forms a β-sheet by association with residues 35–37. Furthermore, we show that the two additional C-terminal residues of Aβ42, in particular Ile-41, directly control the differences in the β-strand content found between the Aβ40 and Aβ42 structural ensembles. Integrating the experimental and theoretical evidence accumulated over the last decade, it is now possible to present monomeric structural ensembles of Aβ40 and Aβ42 consistent with available information that produce a plausible molecular basis for why Aβ42 exhibits greater fibrillization rates than Aβ40.     

Nullspace Sampling with Holonomic Constraints Reveals Molecular Mechanisms of Protein Gαs

Proteins perform their function or interact with partners by exchanging between conformational substates on a wide range of spatiotemporal scales. Structurally characterizing these exchanges is challenging, both experimentally and computationally. Large, diffusional motions are often on timescales that are difficult to access with molecular dynamics simulations, especially for large proteins and their complexes. The low frequency modes of normal mode analysis (NMA) report on molecular fluctuations associated with biological activity. However, NMA is limited to a second order expansion about a minimum of the potential energy function, which limits opportunities to observe diffusional motions. By contrast, kino-geometric conformational sampling (KGS) permits large perturbations while maintaining the exact geometry of explicit conformational constraints, such as hydrogen bonds. Here, we extend KGS and show that a conformational ensemble of the α subunit Gαs of heterotrimeric stimulatory protein Gs exhibits structural features implicated in its activation pathway. Activation of protein Gs by G protein-coupled receptors (GPCRs) is associated with GDP release and large conformational changes of its α-helical domain. Our method reveals a coupled α-helical domain opening motion while, simultaneously, Gαs helix α₅ samples an activated conformation. These motions are moderated in the activated state. The motion centers on a dynamic hub near the nucleotide-binding site of Gαs, and radiates to helix α₄. We find that comparative NMA-based ensembles underestimate the amplitudes of the motion. Additionally, the ensembles fall short in predicting the accepted direction of the full activation pathway. Taken together, our findings suggest that nullspace sampling with explicit, holonomic constraints yields ensembles that illuminate molecular mechanisms involved in GDP release and protein Gs activation, and further establish conformational coupling between key structural elements of Gαs.     

Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 μs. Four combinations of protein and water force fields were tested: ff19SB/OPC , ff19SB/TIP4P-D , ff03CMAP/TIP4P-D , and a99SB-disp/TIP4P-disp , with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 μs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and ³J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.     

Discrete Molecular Dynamics Can Predict Helical Prestructured Motifs in Disordered Proteins

Intrinsically disordered proteins (IDPs) lack a stable tertiary structure, but their short binding regions termed Pre-Structured Motifs (PreSMo) can form transient secondary structure elements in solution. Although disordered proteins are crucial in many biological processes and designing strategies to modulate their function is highly important, both experimental and computational tools to describe their conformational ensembles and the initial steps of folding are sparse. Here we report that discrete molecular dynamics (DMD) simulations combined with replica exchange (RX) method efficiently samples the conformational space and detects regions populating α-helical conformational states in disordered protein regions. While the available computational methods predict secondary structural propensities in IDPs based on the observation of protein-protein interactions, our ab initio method rests on physical principles of protein folding and dynamics. We show that RX-DMD predicts α-PreSMos with high confidence confirmed by comparison to experimental NMR data. Moreover, the method also can dissect α-PreSMos in close vicinity to each other and indicate helix stability. Importantly, simulations with disordered regions forming helices in X-ray structures of complexes indicate that a preformed helix is frequently the binding element itself, while in other cases it may have a role in initiating the binding process. Our results indicate that RX-DMD provides a breakthrough in the structural and dynamical characterization of disordered proteins by generating the structural ensembles of IDPs even when experimental data are not available.     

Evaluating Molecular Mechanical Potentials for Helical Peptides and Proteins

Multiple variants of the AMBER all-atom force field were quantitatively evaluated with respect to their ability to accurately characterize helix-coil equilibria in explicit solvent simulations. Using a global distributed computing network, absolute conformational convergence was achieved for large ensembles of the capped A₂₁ and F_s helical peptides. Further assessment of these AMBER variants was conducted via simulations of a flexible 164-residue five-helix-bundle protein, apolipophorin-III, on the 100 ns timescale. Of the contemporary potentials that had not been assessed previously, the AMBER-99SB force field showed significant helix-destabilizing tendencies, with beta bridge formation occurring in helical peptides, and unfolding of apolipophorin-III occurring on the tens of nanoseconds timescale. The AMBER-03 force field, while showing adequate helical propensities for both peptides and stabilizing apolipophorin-III, (i) predicts an unexpected decrease in helicity with ALA→ARG⁺ substitution, (ii) lacks experimentally observed 3₁₀ helical content, and (iii) deviates strongly from average apolipophorin-III NMR structural properties. As is observed for AMBER-99SB, AMBER-03 significantly overweighs the contribution of extended and polyproline backbone configurations to the conformational equilibrium. In contrast, the AMBER-99φ force field, which was previously shown to best reproduce experimental measurements of the helix-coil transition in model helical peptides, adequately stabilizes apolipophorin-III and yields both an average gyration radius and polar solvent exposed surface area that are in excellent agreement with the NMR ensemble.     

Effects of Hydrophobic Macromolecular Crowders on Amyloid β (16–22) Aggregation

In Alzheimer’s disease (AD), the amyloid β (Aβ) peptide aggregates in the brain to form progressively larger oligomers, fibrils, and plaques. The aggregation process is strongly influenced by the presence of other macromolecular species, called crowders, that can exert forces on the proteins. One very common attribute of macromolecular crowders is their hydrophobicity. We examined the effect of hydrophobic crowders on protein aggregation by using discontinuous molecular dynamics (DMD) simulations in combination with an intermediate resolution protein model, PRIME20. The systems considered contained 48 Aβ (16–22) peptides and crowders with diameters of 5 Å, 20 Å, and 40 Å, represented by hard spheres or spheres with square-well/square-shoulder interactions, at a crowder volume fraction of ϕ = 0.10. Results show that low levels of crowder hydrophobicity are capable of increasing the fibrillation lag time and high levels of crowder hydrophobicity can fully prevent the formation of fibrils. The types of structures that remain during the final stages of the simulations are summarized in a global phase diagram that shows fibril, disordered oligomer, or β-sheet phases in the space spanned by crowder size and crowder hydrophobicity. In particular, at high levels of hydrophobicity, simulations with 5 Å crowders result in only disordered oligomers and simulations with 40 Å crowders result in only β-sheets. The presence of hydrophobic crowders reduces the antiparallel β-sheet content of fibrils, whereas hard sphere crowders increase it. Finally, strong hydrophobic crowders alter the secondary structure of the Aβ (16–22) monomers, bending them into a shape that is incapable of forming ordered β-sheets or fibrils. These results qualitatively agree with previous theoretical and experimental work.     

Mechanism of the αβ Conformational Change in F1-ATPase after ATP Hydrolysis: Free-Energy Simulations

One of the motive forces for F₁-ATPase rotation is the conformational change of the catalytically active β subunit due to closing and opening motions caused by ATP binding and hydrolysis, respectively. The closing motion is accomplished in two steps: the hydrogen-bond network around ATP changes and then the entire structure changes via B-helix sliding, as shown in our previous study. Here, we investigated the opening motion induced by ATP hydrolysis using all-atom free-energy simulations, combining the nudged elastic band method and umbrella sampling molecular-dynamics simulations. Because hydrolysis requires residues in the α subunit, the simulations were performed with the αβ dimer. The results indicate that the large-scale opening motion is also achieved by the B-helix sliding (in the reverse direction). However, the sliding mechanism is different from that of ATP binding because sliding is triggered by separation of the hydrolysis products ADP and Pi. We also addressed several important issues: 1), the timing of the product Pi release; 2), the unresolved half-closed β structure; and 3), the ADP release mechanism. These issues are fundamental for motor function; thus, the rotational mechanism of the entire F₁-ATPase is also elucidated through this αβ study. During the conformational change, conserved residues among the ATPase proteins play important roles, suggesting that the obtained mechanism may be shared with other ATPase proteins. When combined with our previous studies, these results provide a comprehensive view of the β-subunit conformational change that drives the ATPase.     

Probing the Conformation of FhaC with Small-Angle Neutron Scattering and Molecular Modeling

Probing the solution structure of membrane proteins represents a formidable challenge, particularly when using small-angle scattering. Detergent molecules often present residual scattering contributions even at their match point in small-angle neutron scattering (SANS) measurements. Here, we studied the conformation of FhaC, the outer-membrane, β-barrel transporter of the Bordetella pertussis filamentous hemagglutinin adhesin. SANS measurements were performed on homogeneous solutions of FhaC solubilized in n-octyl-d17-βD-glucoside and on a variant devoid of the α helix H1, which critically obstructs the FhaC pore, in two solvent conditions corresponding to the match points of the protein and the detergent, respectively. Protein-bound detergent amounted to 142 ± 10 mol/mol as determined by analytical ultracentrifugation. By using molecular modeling and starting from three distinct conformations of FhaC and its variant embedded in lipid bilayers, we generated ensembles of protein-detergent arrangement models with 120–160 detergent molecules. The scattered curves were back-calculated for each model and compared with experimental data. Good fits were obtained for relatively compact, connected detergent belts, which occasionally displayed small detergent-free patches on the outer surface of the β barrel. The combination of SANS and modeling clearly enabled us to infer the solution structure of FhaC, with H1 inside the pore as in the crystal structure. We believe that our strategy of combining explicit atomic detergent modeling with SANS measurements has significant potential for structural studies of other detergent-solubilized membrane proteins.     

Structural and Dynamic Properties of the Human Prion Protein

Prion diseases involve the conformational conversion of the cellular prion protein (PrP^C) to its misfolded pathogenic form (PrP^Sc). To better understand the structural mechanism of this conversion, we performed extensive all-atom, explicit-solvent molecular-dynamics simulations for three structures of the wild-type human PrP (huPrP) at different pH values and temperatures. Residue 129 is polymorphic, being either Met or Val. Two of the three structures have Met in position 129 and the other has Val. Lowering the pH or raising the temperature induced large conformational changes of the C-terminal globular domain and increased exposure of its hydrophobic core. In some simulations, HA and its preceding S1-HA loop underwent large displacements. The C-terminus of HB was unstable and sometimes partially unfolded. Two hydrophobic residues, Phe-198 and Met-134, frequently became exposed to solvent. These conformational changes became more dramatic at lower pH or higher temperature. Furthermore, Tyr-169 and the S2-HB loop, or the X-loop, were different in the starting structures but converged to common conformations in the simulations for the Met-129, but not the Val-129, protein. α-Strands and β-strands formed in the initially unstructured N-terminus. α-Strand propensity in the N-terminus was different between the Met-129 and Val129 proteins, but β-strand propensity was similar. This study reveals detailed structural and dynamic properties of huPrP, providing insight into the mechanism of the conversion of PrP^C to PrP^Sc.     

Na+/K+-ATPase α-subunit (nkaα) Isoforms and Their mRNA Expression Levels,Overall Nkaα Protein Abundance,and Kinetic Properties of Nka in the Skeletal Muscle and Three Electric Organs of the Electric Eel,Electrophorus electricus

This study aimed to obtain the coding cDNA sequences of Na⁺/K⁺-ATPase α (nkaα) isoforms from, and to quantify their mRNA expression in, the skeletal muscle (SM), the main electric organ (EO), the Hunter’s EO and the Sach’s EO of the electric eel, Electrophorus electricus. Four nkaα isoforms (nkaα1c1, nkaα1c2, nkaα2 and nkaα3) were obtained from the SM and the EOs of E. electricus. Based on mRNA expression levels, the major nkaα expressed in the SM and the three EOs of juvenile and adult E. electricus were nkaα1c1 and nkaα2, respectively. Molecular characterization of the deduced Nkaα1c1 and Nkaα2 sequences indicates that they probably have different affinities to Na⁺ and K⁺. Western blotting demonstrated that the protein abundance of Nkaα was barely detectable in the SM, but strongly detected in the main and Hunter’s EOs and weakly in the Sach’s EO of juvenile and adult E. electricus. These results corroborate the fact that the main EO and Hunter’s EO have high densities of Na⁺ channels and produce high voltage discharges while the Sach’s EO produces low voltage discharges. More importantly, there were significant differences in kinetic properties of Nka among the three EOs of juvenile E. electricus. The highest and lowest V _max of Nka were detected in the main EO and the Sach’s EO, respectively, with the Hunter’s EO having a V _max value intermediate between the two, indicating that the metabolic costs of EO discharge could be the highest in the main EO. Furthermore, the Nka from the main EO had the lowest K_m (or highest affinity) for Na⁺ and K⁺ among the three EOs, suggesting that the Nka of the main EO was more effective than those of the other two EOs in maintaining intracellular Na⁺ and K⁺ homeostasis and in clearing extracellular K⁺ after EO discharge.     

Immunoprecipitation of Amyloid Fibrils by the Use of an Antibody that Recognizes a Generic Epitope Common to Amyloid Fibrils

Amyloid fibrils are associated with many maladies, including Alzheimer’s disease (AD). The isolation of amyloids from natural materials is very challenging because the extreme structural stability of amyloid fibrils makes it difficult to apply conventional protein science protocols to their purification. A protocol to isolate and detect amyloids is desired for the diagnosis of amyloid diseases and for the identification of new functional amyloids. Our aim was to develop a protocol to purify amyloid from organisms, based on the particular characteristics of the amyloid fold, such as its resistance to proteolysis and its capacity to be recognized by specific conformational antibodies. We used a two-step strategy with proteolytic digestion as the first step followed by immunoprecipitation using the amyloid conformational antibody LOC. We tested the efficacy of this method using as models amyloid fibrils produced in vitro, tissue extracts from C. elegans that overexpress Aβ peptide, and cerebrospinal fluid (CSF) from patients diagnosed with AD. We were able to immunoprecipitate Aβ_1–40 amyloid fibrils, produced in vitro and then added to complex biological extracts, but not α-synuclein and gelsolin fibrils. This method was useful for isolating amyloid fibrils from tissue homogenates from a C. elegans AD model, especially from aged worms. Although we were able to capture picogram quantities of Aβ_1–40 amyloid fibrils produced in vitro when added to complex biological solutions, we could not detect any Aβ amyloid aggregates in CSF from AD patients. Our results show that although immunoprecipitation using the LOC antibody is useful for isolating Aβ_1–40 amyloid fibrils, it fails to capture fibrils of other amyloidogenic proteins, such as α-synuclein and gelsolin. Additional research might be needed to improve the affinity of these amyloid conformational antibodies for an array of amyloid fibrils without compromising their selectivity before application of this protocol to the isolation of amyloids.     

The Cytoprotective Effects of E-α-(4-Methoxyphenyl)-2’,3,4,4'-Tetramethoxychalcone (E-α-p-OMe-C6H4-TMC)—A Novel and Non-Cytotoxic HO-1 Inducer

Cell protection against different noxious stimuli like oxidative stress or chemical toxins plays a central role in the treatment of many diseases. The inducible heme oxygenase isoform, heme oxygenase-1 (HO-1), is known to protect cells against a variety of harmful conditions including apoptosis. Because a number of medium strong electrophiles from a series of α-X-substituted 2’,3,4,4’-tetramethoxychalcones (α-X-TMCs, X = H, F, Cl, Br, I, CN, Me, p-NO₂-C₆H₄, Ph, p-OMe-C₆H₄, NO₂, CF₃, COOEt, COOH) had proven to activate Nrf2 resulting in HO-1 induction and inhibit NF-κB downstream target genes, their protective effect against staurosporine induced apoptosis and reactive oxygen species (ROS) production was investigated. RAW264.7 macrophages treated with 19 different chalcones (15 α-X-TMCs, chalcone, 2’-hydroxychalcone, calythropsin and 2’-hydroxy-3,4,4’-trimethoxychalcone) prior to staurosporine treatment were analyzed for apoptosis and ROS production, as well as HO-1 protein expression and enzyme activity. Additionally, Nrf2 and NF-κB activity was assessed. We found that amongst all tested chalcones only E-α-(4-methoxyphenyl)-2’,3,4,4''-tetramethoxychalcone (E-α-p-OMe-C₆H₄-TMC) demonstrated a distinct, statistically significant antiapoptotic effect in a dose dependent manner, showing no toxic effects, while its double bond isomer Z-α-p-OMe-C₆H₄-TMC displayed no significant activity. Also, E-α-p-OMe-C₆H₄-TMC induced HO-1 protein expression and increased HO-1 activity, whilst inhibition of HO-1 by SnPP-IX abolished its antiapoptotic effect. The only weakly electrophilic chalcone E-α-p-OMe-C₆H₄-TMC reduced the staurosporine triggered formation of ROS, while inducing the translocation of Nrf2 into the nucleus. Furthermore, staurosporine induced NF-κB activity was attenuated following E-α-p-OMe-C₆H₄-TMC treatment. Overall, E-α-p-OMe-C₆H₄-TMC demonstrated its effective cytoprotective potential via a non-toxic induction of HO-1 in RAW264.7 macrophages. The observed cytoprotective effect may partly be related to both, the activation of the Nrf2- and inhibition of the NF-κB pathway.     

Effects of force fields on the conformational and dynamic properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1‐40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER‐ff99sb‐ILDN, AMBER‐ff99sb*‐ILDN, AMBER‐ff99sb‐NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER‐ff99sb‐ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α‐helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER‐ff99sb‐NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER‐ff99sb‐NMR force field, the others tended to under estimate the expected amount of β‐sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER‐ff99sb‐NMR, reproduce a theoretically expected β‐sheet‐turn‐β‐sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C‐terminal hydrophobic cores from residues 17‐21 and 30‐36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.     

Neuronal Cx3cr1 Deficiency Protects against Amyloid β-Induced Neurotoxicity

Cx3cr1, the receptor for the chemokine Cx3cl1 (fractalkine), has been implicated in the progression and severity of Alzheimer’s disease-like pathology in mice, but the underlying mechanisms remain unclear. A complicating factor is that Cx3cr1 has been demonstrated in both neurons and microglia. Here, we have dissected the differences between neuronal and microglial Cx3cr1, specifically by comparing direct amyloid-β-induced toxicity in cultured, mature, microglia-depleted hippocampal neurons from wild-type and Cx3cr1^-/- mice. Wild-type neurons expressed both Cx3cl1 and Cx3cr1 and released Cx3cl1 in response to amyloid-β. Knockout of neuronal Cx3cr1 abated amyloid-β-induced lactate dehydrogenase release. Furthermore, amyloid-β differentially induced depression of pre- and postsynaptic components of miniature excitatory postsynaptic currents, in a peptide conformation-dependent manner. Knockout of neuronal Cx3cr1 abated effects of both amyloid-β conformational states, which were differentiable by aggregation kinetics and peptide morphology. We obtained similar results after both acute and chronic treatment of cultured neurons with the Cx3cr1 antagonist F1. Thus, neuronal Cx3cr1 may impact Alzheimer’s disease-like pathology by modulating conformational state-dependent amyloid-β-induced synaptotoxicity.     

Combined application of rapamycin and atorvastatin improves lipid metabolism in apolipoprotein E-deficient mice with chronic kidney disease

Atherosclerosis arising from the pro-inflammatory conditions associated with chronic kidney disease (CKD) increases major cardiovascular morbidity and mortality. Rapamycin (RAPA) is known to inhibit atherosclerosis under CKD and non-CKD conditions, but it can cause dyslipidemia; thus, the co-application of lipid-lowering agents is recommended. Atorvastatin (ATV) has been widely used to reduce serum lipids levels, but its synergistic effect with RAPA in CKD remains unclear. Here, we analyzed the effect of their combined treatment on atherosclerosis stimulated by CKD in apolipoprotein E-deficient (ApoE^−/−) mice. Oil Red O staining revealed that treatment with RAPA and RAPA＋ ATV, but not ATV alone, significantly decreased the atherosclerotic lesions in the aorta and aortic sinus, compared to those seen in the control (CKD) group. The co-administration of RAPA and ATV improved the serum lipid profile and raised the expression levels of proteins involved in reverse cholesterol transport (LXRα, CYP7A1, ABCG1, PPARγ, ApoA1) in the liver. The CKD group showed increased levels of various genes encoding atherosclerosis-promoting cytokines in the spleen (Tnf-α, Il-6 and Il-1β) and aorta (Tnf-α and Il-4), and these increases were attenuated by RAPA treatment. ATV and RAPA＋ATV decreased the levels of Tnf-α and Il-1β in the spleen, but not in the aorta. Together, these results indicate that, in CKD-induced ApoE^−/− mice, RAPA significantly reduces the development of atherosclerosis by regulating the expression of inflammatory cytokines and the co-application of ATV improves lipid metabolism.