Anne M. Brown

Comparing atomistic molecular mechanics force fields for a difficult target: a case study on the Alzheimer’s amyloid β-peptide

Macromolecular function arises from structure, and many diseases are associated with misfolding of proteins. Molecular simulation methods can augment experimental techniques to understand misfolding and aggregation pathways with atomistic resolution, but the reliability of these predictions is a function of the parameters used for the simulation. There are many biomolecular force fields available, but most are validated using stably folded structures. Here, we present the results of molecular dynamics simulations on the intrinsically disordered amyloid β-peptide (Aβ), whose misfolding and aggregation give rise to the symptoms of Alzheimer’s disease. Because of the link between secondary structure changes and pathology, being able to accurately model the structure of Aβ would greatly improve our understanding of this disease, and it may facilitate application of modeling approaches to other protein misfolding disorders. To this end, we compared five popular atomistic force fields (AMBER03, CHARMM22 + CMAP, GROMOS96 53A6, GROMOS96 54A7, and OPLS-AA) to determine which could best model the structure of Aβ. By comparing secondary structure content, NMR shifts, and radius of gyration to available experimental data, we conclude that AMBER03 and CHARMM22 + CMAP over-stabilize helical structure within Aβ, with CHARMM22 + CMAP also producing elongated Aβ structures, in conflict with experimental findings. OPLS-AA, GROMOS96 53A6, and GROMOS96 54A7 produce very similar results in terms of helical and β-strand content, calculated NMR shifts, and radii of gyration that agree well with experimental data.

Simulations of monomeric amyloid β-peptide (1-40) with varying solution conditions and oxidation state of Met35: implications for aggregation.

The amyloid β-peptide (Aβ) is a 40-42 residue peptide that is the principal toxic species in Alzheimers disease (AD). The oxidation of methionine-35 (Met35) to the sulfoxide form (Met35(ox)) has been identified as potential modulator of Aβ aggregation. The role Met35(ox) plays in Aβ neurotoxicity differs among experimental studies, which may be due to inconsistent solution conditions (pH, buffer, temperature). We applied atomistic molecular dynamics (MD) simulations as a means to probe the dynamics of the monomeric 40-residue alloform of Aβ (Aβ40) containing Met35 or Met35(ox) in an effort to resolve the conflicting experimental results. We found that Met35 oxidation decreases the β-strand content of the C-terminal hydrophobic region (residues 29-40), with a specific effect on the secondary structure of residues 33-35, thus potentially impeding aggregation. Further, there is an important interplay between oxidation state and solution conditions, with pH and salt concentration augmenting the effects of oxidation. The results presented here serve to rationalize the conflicting results seen in experimental studies and provide a fundamental biophysical characterization of monomeric Aβ40 dynamics in both reduced and oxidized forms, providing insight into the biochemical mechanism of Aβ40 and oxidative stress related to AD.

Molecular Dynamics Simulations of Amyloid β-Peptide (1-42): Tetramer Formation and Membrane Interactions

The aggregation cascade and peptide-membrane interactions of the amyloid β-peptide (Aβ) have been implicated as toxic events in the development and progression of Alzheimers disease. Aβ42 forms oligomers and ultimately plaques, and it has been hypothesized that these oligomeric species are the main toxic species contributing to neuronal cell death. To better understand oligomerization events and subsequent oligomer-membrane interactions of Aβ42, we performed atomistic molecular-dynamics (MD) simulations to characterize both interpeptide interactions and perturbation of model membranes by the peptides. MD simulations were utilized to first show the formation of a tetramer unit by four separate Aβ42 peptides. Aβ42 tetramers adopted an oblate ellipsoid shape and showed a significant increase in β-strand formation in the final tetramer unit relative to the monomers, indicative of on-pathway events for fibril formation. The Aβ42 tetramer unit that formed in the initial simulations was used in subsequent MD simulations in the presence of a pure POPC or cholesterol-rich raft model membrane. Tetramer-membrane simulations resulted in elongation of the tetramer in the presence of both model membranes, with tetramer-raft interactions giving rise to the rearrangement of key hydrophobic regions in the tetramer and the formation of a more rod-like structure indicative of a fibril-seeding aggregate. Membrane perturbation by the tetramer was manifested in the form of more ordered, rigid membranes, with the pure POPC being affected to a greater extent than the raft membrane. These results provide critical atomistic insight into the aggregation pathway of Aβ42 and a putative toxic mechanism in the pathogenesis of Alzheimers disease.

Purine biosynthesis in archaea: variations on a theme

BackgroundThe ability to perform de novo biosynthesis of purines is present in organisms in all three domains of life, reflecting the essentiality of these molecules to life. Although the pathway is quite similar in eukaryotes and bacteria, the archaeal pathway is more variable. A careful manual curation of genes in this pathway demonstrates the value of manual curation in archaea, even in pathways that have been well-studied in other domains.ResultsWe searched the Integrated Microbial Genome system (IMG) for the 17 distinct genes involved in the 11 steps of de novo purine biosynthesis in 65 sequenced archaea, finding 738 predicted proteins with sequence similarity to known purine biosynthesis enzymes. Each sequence was manually inspected for the presence of active site residues and other residues known or suspected to be required for function.Many apparently purine-biosynthesizing archaea lack evidence for a single enzyme, either glycinamide ribonucleotide formyltransferase or inosine monophosphate cyclohydrolase, suggesting that there are at least two more gene variants in the purine biosynthetic pathway to discover. Variations in domain arrangement of formylglycinamidine ribonucleotide synthetase and substantial problems in aminoimidazole carboxamide ribonucleotide formyltransferase and inosine monophosphate cyclohydrolase assignments were also identified.Manual curation revealed some overly specific annotations in the IMG gene product name, with predicted proteins without essential active site residues assigned product names implying enzymatic activity (21 proteins, 2.8% of proteins inspected) or Enzyme Commission (E. C.) numbers (57 proteins, 7.7%). There were also 57 proteins (7.7%) assigned overly generic names and 78 proteins (10.6%) without E.C. numbers as part of the assigned name when a specific enzyme name and E. C. number were well-justified.ConclusionsThe patchy distribution of purine biosynthetic genes in archaea is consistent with a pathway that has been shaped by horizontal gene transfer, duplication, and gene loss. Our results indicate that manual curation can improve upon automated annotation for a small number of automatically-annotated proteins and can reveal a need to identify further pathway components even in well-studied pathways.ReviewersThis article was reviewed by Dr. Céline Brochier-Armanet, Dr Kira S Makarova (nominated by Dr. Eugene Koonin), and Dr. Michael Galperin.

Transforming Sphingosine Kinase 1 Inhibitors into Dual and Sphingosine Kinase 2 Selective Inhibitors: Design, Synthesis, and in Vivo Activity

Sphingosine 1-phosphate (S1P) is a pleiotropic signaling molecule that interacts with its five G-protein coupled receptors (S1P1-5) to regulate cell growth and survival and has been implicated in a variety of diseases including cancer and sickle cell disease. As the key mediators in the synthesis of S1P, sphingosine kinase (SphK) isoforms 1 and 2 have attracted attention as viable targets for pharmaceutical inhibition. In this article, we describe the design, synthesis, and biological evaluation of aminothiazole-based guanidine inhibitors of SphK. Surprisingly, combining features of reported SphK1 inhibitors generated SphK1/2 dual inhibitor 20l (SLC4011540) (hSphK1 Ki = 120 nM, hSphK2 Ki = 90 nM) and SphK2 inhibitor 20dd (SLC4101431) (Ki = 90 nM, 100-fold SphK2 selectivity). These compounds effectively decrease S1P levels in vitro. In vivo administration of 20dd validated that inhibition of SphK2 increases blood S1P levels.

Structure-Activity Relationship Studies and Molecular Modeling of Naphthalene-Based Sphingosine Kinase 2 Inhibitors.

The two isoforms of sphingosine kinase (SphK1 and SphK2) are the only enzymes that phosphorylate sphingosine to sphingosine-1-phosphate (S1P), which is a pleiotropic lipid mediator involved in a broad range of cellular processes including migration, proliferation, and inflammation. SphKs are targets for various diseases such as cancer, fibrosis, and Alzheimers and sickle cell disease. Herein, we disclose the structure-activity profile of naphthalene-containing SphK inhibitors and molecular modeling studies that reveal a key molecular switch that controls SphK selectivity.

Development of a structured undergraduate research experience: Framework and implications

Participating in undergraduate research can be a pivotal experience for students in life science disciplines. Development of critical thinking skills, in addition to conveying scientific ideas in oral and written formats, is essential to ensuring that students develop a greater understanding of basic scientific knowledge and the research process. Modernizing the current life sciences research environment to accommodate the growing demand by students for experiential learning is needed. By developing and implementing a structured, theory‐based approach to undergraduate research in the life sciences, specifically biochemistry, it has been successfully shown that more students can be provided with a high‐quality, high‐impact research experience. The structure of this approach allowed students to develop novel, independent projects in a computational molecular modeling lab. Students engaged in an experience in which career goals, problem‐solving skills, time management skills, and independence in a research lab were developed. After experiencing this approach to undergraduate research, students reported feeling challenged to think critically and prepared for future career paths. The approach allowed for a progressive learning environment where more undergraduate students could participate in publishable research. Future areas for development include implementation in a bench‐top lab and extension to disciplines beyond biochemistry. In this study, it has been shown that utilizing the structured approach to undergraduate research could allow for more students to experience undergraduate research and develop into more confident, independent life scientists well prepared for graduate schools and professional research environments.

Computational study of HIV gp120 as a target for polyanionic entry inhibitors: Exploiting the V3 loop region

Multiple approaches are being utilized to develop therapeutics to treat HIV infection. One approach is designed to inhibit entry of HIV into host cells, with a target being the viral envelope glycoprotein, gp120. Polyanionic compounds have been shown to be effective in inhibiting HIV entry, with a mechanism involving electrostatic interactions with the V3 loop of gp120 being proposed. In this study, we applied computational methods to elucidate molecular interactions between the repeat unit of the precisely alternating polyanion, Poly(4,4′-stilbenedicarboxylate-alt–maleic acid) (DCSti-alt-MA) and the V3 loop of gp120 from strains of HIV against which these polyanions were previously tested (IIIb, BaL, 92UG037, JR-CSF) as well as two strains for which gp120 crystal structures are available (YU2, 2B4C). Homology modeling was used to create models of the gp120 proteins. Using monomers of the gp120 protein, we applied extensive molecular dynamics simulations to obtain dominant morphologies that represent a variety of open-closed states of the V3 loop to examine the interaction of 112 ligands of the repeating units of DCSti-alt-MA docked to the V3 loop and surrounding residues. Using the distance between the V1/V2 and V3 loops of gp120 as a metric, we revealed through MD simulations that gp120 from the lab-adapted strains (BaL and IIIb), which are more susceptible to inhibition by DCSti-alt-MA, clearly transitioned to the closed state in one replicate of each simulation set, whereas none of the replicates from the Tier II strains (92UG037 and JR-CSF) did so. Docking repeat unit microspecies to the gp120 protein before and after MD simulation enabled identification of residues that were key for binding. Notably, only a few residues were found to be important for docking both before and after MD simulation as a result of the conformational heterogeneity provided by the simulations. Consideration of the residues that were consistently involved in interactions with the ligand revealed the importance of both hydrophilic and hydrophobic moieties of the ligand for effective binding. The results also suggest that polymers of DCSti-alt-MA with repeating units of different configurations may have advantages for therapeutic efficacy.

Influence of Sequence and Lipid Type on Membrane Perturbation by Human and Rat Amyloid β-Peptide (1-42)

The hallmark characteristics of plaque formation and neuronal cell death in Alzheimers disease (AD) are caused principally by the amyloid β-peptide (Aβ). Aβ sequence and lipid composition are essential variables to consider when elucidating the impact of biological membranes on Aβ structure and the effect of Aβ on membrane integrity. Atomistic molecular dynamics simulations testing two Aβ sequences, human and rat Aβ (HAβ and RAβ, respectively), and five lipid types were performed to assess the effect of these variables on membrane perturbation and potential link to AD phenotype differences based on differences in sequence. All metrics agree insomuch that monomeric HAβ and RAβ contribute to membrane perturbation by causing a more rigid, gel-like lipid phase. Differences between HAβ and RAβ binding on degree of membrane perturbation were based on lipid headgroup properties. Cholesterol was found to moderate the amount of perturbation caused by HAβ and RAβ in a model raft membrane. The difference in position of an arginine residue between HAβ and RAβ influenced peptide-membrane interactions and was determined to be the mediating factor in observed differences in lipid affinity and degree of membrane disruption. Overall, this work increases our understanding of the influence of sequence and lipid type on Aβ-membrane interactions and their relationship to AD.

Insights into Stabilizing Forces in Amyloid Fibrils of Differing Sizes from Polarizable Molecular Dynamics Simulations

Pathological aggregation of amyloid-forming proteins is a hallmark of a number of human diseases, including Alzheimers, type 2 diabetes, Parkinsons, and more. Despite having very different primary amino acid sequences, these amyloid proteins form similar supramolecular, fibril structures that are highly resilient to physical and chemical denaturation. To better understand the structural stability of disease-related amyloids and to gain a greater understanding of factors that stabilize functional amyloid assemblies, insights into tertiary and quaternary interactions are needed. We performed molecular dynamics simulations on human tau, amyloid-β, and islet amyloid polypeptide fibrils to determine key physicochemical properties that give rise to their unique characteristics and fibril structures. These simulations are the first of their kind in employing a polarizable force field to explore properties of local electric fields on dipole properties and other electrostatic forces that contribute to amyloid stability. Across these different amyloid fibrils, we focused on how the underlying forces stabilize fibrils to elucidate the driving forces behind the protein aggregation. The polarizable model allows for an investigation of how side-chain dipole moments, properties of structured water molecules in the fibril core, and the local environment around salt bridges contribute to the formation of interfaces essential for fibril stability. By systematically studying three amyloidogenic proteins of various fibril sizes for key structural properties and stabilizing forces, we shed light on properties of amyloid structures related to both diseased and functional states at the atomistic level.