Characterizing ATP processing by the AAA+ protein p97 at the atomic level.

The human enzyme p97 regulates various cellular pathways by unfolding hundreds of protein substrates in an ATP-dependent manner, making it an essential component of protein homeostasis and an impactful pharmacological target. The hexameric complex undergoes substantial conformational changes throughout its catalytic cycle. Here we elucidate the molecular motions that occur at the active site in the temporal window immediately before and after ATP hydrolysis by merging cryo-EM, NMR spectroscopy and molecular dynamics simulations. p97 populates a metastable reaction intermediate, the ADP·Pi state, which is poised between hydrolysis and product release. Detailed snapshots reveal that the active site is finely tuned to trap and eventually discharge the cleaved phosphate. Signalling pathways originating at the active site coordinate the action of the hexamer subunits and couple hydrolysis with allosteric conformational changes. Our multidisciplinary approach enables a glimpse into the sophisticated spatial and temporal orchestration of ATP handling by a prototype AAA+ protein.

Structural basis of metabolite transport by the chloroplast outer envelope channel OEP21.

Triose phosphates (TPs) are the primary products of photosynthetic CO2 fixation in chloroplasts, which need to be exported into the cytosol across the chloroplast inner envelope (IE) and outer envelope (OE) membranes to sustain plant growth. While transport across the IE is well understood, the mode of action of the transporters in the OE remains unclear. Here we present the high-resolution nuclear magnetic resonance (NMR) structure of the outer envelope protein 21 (OEP21) from garden pea, the main exit pore for TPs in C3 plants. OEP21 is a cone-shaped ß-barrel pore with a highly positively charged interior that enables binding and translocation of negatively charged metabolites in a competitive manner, up to a size of ~1 kDa. ATP stabilizes the channel and keeps it in an open state. Despite the broad substrate selectivity of OEP21, these results suggest that control of metabolite transport across the OE might be possible.

A Chemical Proteomic Strategy Reveals Inhibitors of Lipoate Salvage in Bacteria and Parasites.

The development of novel anti-infectives requires unprecedented strategies targeting pathways which are solely present in pathogens but absent in humans. Following this principle, we developed inhibitors of lipoic acid (LA) salvage, a crucial pathway for the survival of LA auxotrophic bacteria and parasites but non-essential in human cells. An LA-based probe was selectively transferred onto substrate proteins via lipoate protein ligase (LPL) in intact cells, and their binding sites were determined by mass spectrometry. Probe labeling served as a proxy of LPL activity, enabling in situ screenings for cell-permeable LPL inhibitors. Profiling a focused compound library revealed two substrate analogs (LAMe and C3) as inhibitors, which were further validated by binding studies and co-crystallography. Importantly, LAMe exhibited low toxicity in human cells and achieved killing of Plasmodium falciparum in erythrocytes with an EC50 value of 15 µM, making it the most effective LPL inhibitor reported to date.

The human signal peptidase complex acts as a quality control enzyme for membrane proteins.

Cells need to detect and degrade faulty membrane proteins to maintain homeostasis. In this study, we identify a previously unknown function of the human signal peptidase complex (SPC)-the enzyme that removes endoplasmic reticulum (ER) signal peptides-as a membrane protein quality control factor. We show that the SPC cleaves membrane proteins that fail to correctly fold or assemble into their native complexes at otherwise hidden cleavage sites, which our study reveals to be abundant in the human membrane proteome. This posttranslocational cleavage synergizes with ER-associated degradation to sustain membrane protein homeostasis and contributes to cellular fitness. Cryptic SPC cleavage sites thus serve as predetermined breaking points that, when exposed, help to target misfolded or surplus proteins for degradation, thereby maintaining a healthy membrane proteome.

Amyloid fibril structure from the vascular variant of systemic AA amyloidosis.

Systemic AA amyloidosis is a debilitating protein misfolding disease in humans and animals. In humans, it occurs in two variants that are called 'vascular' and 'glomerular', depending on the main amyloid deposition site in the kidneys. Using cryo electron microscopy, we here show the amyloid fibril structure underlying the vascular disease variant. Fibrils purified from the tissue of such patients are mainly left-hand twisted and contain two non-equal stacks of fibril proteins. They contrast in these properties to the fibrils from the glomerular disease variant which are right-hand twisted and consist of two structurally equal stacks of fibril proteins. Our data demonstrate that the different disease variants in systemic AA amyloidosis are associated with different fibril morphologies.

Intramembrane client recognition potentiates the chaperone functions of calnexin.

One-third of the human proteome is comprised of membrane proteins, which are particularly vulnerable to misfolding and often require folding assistance by molecular chaperones. Calnexin (CNX), which engages client proteins via its sugar-binding lectin domain, is one of the most abundant ER chaperones, and plays an important role in membrane protein biogenesis. Based on mass spectrometric analyses, we here show that calnexin interacts with a large number of nonglycosylated membrane proteins, indicative of additional nonlectin binding modes. We find that calnexin preferentially bind misfolded membrane proteins and that it uses its single transmembrane domain (TMD) for client recognition. Combining experimental and computational approaches, we systematically dissect signatures for intramembrane client recognition by calnexin, and identify sequence motifs within the calnexin TMD region that mediate client binding. Building on this, we show that intramembrane client binding potentiates the chaperone functions of calnexin. Together, these data reveal a widespread role of calnexin client recognition in the lipid bilayer, which synergizes with its established lectin-based substrate binding. Molecular chaperones thus can combine different interaction modes to support the biogenesis of the diverse eukaryotic membrane proteome.

Enzyme-substrate interface targeting by imidazole-based γ-secretase modulators activates γ-secretase and stabilizes its interaction with APP.

Alzheimer's disease (AD) pathogenesis has been linked to the accumulation of longer, aggregation-prone amyloid ß (Aß) peptides in the brain. Γ-secretases generate Aß peptides from the amyloid precursor protein (APP). Γ-secretase modulators (GSMs) promote the generation of shorter, less-amyloidogenic Aßs and have therapeutic potential. However, poorly defined drug-target interactions and mechanisms of action have hampered their therapeutic development. Here, we investigate the interactions between the imidazole-based GSM and its target γ-secretase-APP using experimental and in silico approaches. We map the GSM binding site to the enzyme-substrate interface, define a drug-binding mode that is consistent with functional and structural data, and provide molecular insights into the underlying mechanisms of action. In this respect, our analyses show that occupancy of a γ-secretase (sub)pocket, mediating binding of the modulator's imidazole moiety, is sufficient to trigger allosteric rearrangements in γ-secretase as well as stabilize enzyme-substrate interactions. Together, these findings may facilitate the rational design of new modulators of γ-secretase with improved pharmacological properties.

Cryo-EM demonstrates the in vitro proliferation of an ex vivo amyloid fibril morphology by seeding.

Several studies showed that seeding of solutions of monomeric fibril proteins with ex vivo amyloid fibrils accelerated the kinetics of fibril formation in vitro but did not necessarily replicate the seed structure. In this research we use cryo-electron microscopy and other methods to analyze the ability of serum amyloid A (SAA)1.1-derived amyloid fibrils, purified from systemic AA amyloidosis tissue, to seed solutions of recombinant SAA1.1 protein. We show that 98% of the seeded fibrils remodel the full fibril structure of the main ex vivo fibril morphology, which we used for seeding, while they are notably different from unseeded in vitro fibrils. The seeded fibrils show a similar proteinase K resistance as ex vivo fibrils and are substantially more stable to proteolytic digestion than unseeded in vitro fibrils. Our data support the view that the fibril morphology contributes to determining proteolytic stability and that pathogenic amyloid fibrils arise from proteolytic selection.

Molecular mechanism of amyloidogenic mutations in hypervariable regions of antibody light chains.

Systemic light chain (AL) amyloidosis is a fatal protein misfolding disease in which excessive secretion, misfolding, and subsequent aggregation of free antibody light chains eventually lead to deposition of amyloid plaques in various organs. Patient-specific mutations in the antibody VL domain are closely linked to the disease, but the molecular mechanisms by which certain mutations induce misfolding and amyloid aggregation of antibody domains are still poorly understood. Here, we compare a patient VL domain with its nonamyloidogenic germline counterpart and show that, out of the five mutations present, two of them strongly destabilize the protein and induce amyloid fibril formation. Surprisingly, the decisive, disease-causing mutations are located in the highly variable complementarity determining regions (CDRs) but exhibit a strong impact on the dynamics of conserved core regions of the patient VL domain. This effect seems to be based on a deviation from the canonical CDR structures of CDR2 and CDR3 induced by the substitutions. The amyloid-driving mutations are not necessarily involved in propagating fibril formation by providing specific side chain interactions within the fibril structure. Rather, they destabilize the VL domain in a specific way, increasing the dynamics of framework regions, which can then change their conformation to form the fibril core. These findings reveal unexpected influences of CDR-framework interactions on antibody architecture, stability, and amyloid propensity.

Altered Hinge Conformations in APP Transmembrane Helix Mutants May Affect Enzyme-Substrate Interactions of γ-Secretase.

Cleavage of substrates by γ-secretase is an inherently slow process where substrate-enzyme affinities cannot be broken down into specific sequence requirements in contrast to soluble proteases. Nevertheless, despite its apparent sequence tolerance single point mutations in amyloid precursor protein can severely affect cleavage efficiencies and change product line preferences. We have determined by NMR spectroscopy the structures of the transmembrane domain of amyloid precursor protein in TFE/water and compared it to that of four mutants: two FAD mutants, V44M and I45T, and the two diglycine hinge mutants, G38L and G38P. In accordance with previous publications, the transmembrane domain is composed of two helical segments connected by the diglycine hinge. Mutations alter kink angles and structural flexibility. Furthermore, to our surprise, we observe different, but specific mutual orientations of N- and C-terminal helical segments in the four mutants compared to the wildtype. We speculate that the observed orientations for G38L, G38P, V44M, and I45T lead to unfavorable interactions with γ-secretase exosites during substrate movement to the enzyme's active site in presenilin and/or for the accommodation into the substrate-binding cavity of presenilin.

CHARMM-GUI supports the Amber force fields.

As part of our ongoing efforts to support diverse force fields and simulation programs in CHARMM-GUI, this work presents the development of FF-Converter to prepare Amber simulation inputs with various Amber force fields within the current CHARMM-GUI workflow. The currently supported Amber force fields are ff14SB/ff19SB (protein), Bsc1 (DNA), OL3 (RNA), GLYCAM06 (carbohydrate), Lipid17 (lipid), GAFF/GAFF2 (small molecule), TIP3P/TIP4P-EW/OPC (water), and 12-6-4 ions, and more will be added if necessary. The robustness and usefulness of this new CHARMM-GUI extension are demonstrated by two exemplary systems: a protein/N-glycan/ligand/membrane system and a protein/DNA/RNA system. Currently, CHARMM-GUI supports the Amber force fields only for the Amber program, but we will expand the FF-Converter functionality to support other simulation programs that support the Amber force fields.

The dynamics of γ-secretase and its substrates.

γ-Secretase is an intramembrane aspartyl-protease catalyzing the final step in the regulated intramembrane proteolysis of a large number of single-span type-1 transmembrane proteins. The most extensively studied substrates are the amyloid-ß precursor protein (APP) and the NOTCH receptors. An important technique for the characterization of interactions and conformational changes enabling γ-secretase to perform hydrolysis within the confines of the membrane are molecular dynamics simulations on different time and length scales. Here, we review structural and dynamical features of γ-secretase and its substrates including flexibility descriptions from simulations and experiments. We address (1) conformational sampling of apo-enzyme and unbound substrates (exemplified for APP, NOTCH1 and the apparent non-substrate integrin ß1), (2) substrate recruitment pathways, (3) conformational changes associated with the formation of the recognition complex, (4) cleavage-site unfolding upon interaction with the enzyme's active site, (5) substrate processing after endoproteolysis, and (6) inhibition and modulation of γ-secretase. We conclude with still open questions and suggest further investigations in order to advance our understanding on how γ-secretase selects and processes substrates. This knowledge will improve the ability to better target substrates selectively for therapeutic applications.

Uncovering the Binding Mode of γ -Secretase Inhibitors.

Knowledge of how transition state inhibitors bind to γ-secretase is of major importance for the design of new Alzheimer's disease therapies. On the basis of the known structure of γ-secretase in complex with a fragment of the amyloid precursor protein, we generated a structural model of γ-secretase in complex with the effective L-685,458 transition state inhibitor. The predicted binding mode is in excellent agreement with experimental data, mimicking all enzyme-substrate interactions at the active site and forming the relevant transition state geometry with the active site aspartate residues. The model also indicates the possible location and nature of the amino acid residues forming the proposed binding pockets S1', S2', and S3' near the active site that are occupied by chemical groups of the inhibitor. In addition, we found that the stability of the complex is very likely sensitive to the pH value. Comparative simulations on the binding of L-685,458 and the epimer L682,679 allowed us to explain the strongly reduced affinity of the epimer for γ-secretase. The structural model could form a valuable basis for the design of new or modified γ-secretase inhibitors.

Extracellular interface between APP and Nicastrin regulates Aß length and response to γ-secretase modulators.

γ-Secretase complexes (GSECs) are multimeric membrane proteases involved in a variety of physiological processes and linked to Alzheimer's disease (AD). Presenilin (PSEN, catalytic subunit), Nicastrin (NCT), Presenilin Enhancer 2 (PEN-2), and Anterior Pharynx Defective 1 (APH1) are the essential subunits of GSECs. Mutations in PSEN and the Amyloid Precursor Protein (APP) cause early-onset AD GSECs successively cut APP to generate amyloid-ß (Aß) peptides of various lengths. AD-causing mutations destabilize GSEC-APP/Aßn interactions and thus enhance the production of longer Aßs, which elicit neurotoxic events underlying pathogenesis. Here, we investigated the molecular strategies that anchor GSEC and APP/Aßn during the sequential proteolysis. Our studies reveal that a direct interaction between NCT ectodomain and APPC99 influences the stability of GSEC-Aßn assemblies and thereby modulates Aß length. The data suggest a potential link between single-nucleotide variants in NCSTN and AD risk. Furthermore, our work indicates that an extracellular interface between the protease (NCT, PSEN) and the substrate (APP) represents the target for compounds (GSMs) modulating Aß length. Our findings may guide future rationale-based drug discovery efforts.

γ-Secretase Studied by Atomistic Molecular Dynamics Simulations: Global Dynamics, Enzyme Activation, Water Distribution and Lipid Binding.

γ-secretase, an intramembrane-cleaving aspartyl protease is involved in the cleavage of a large number of intramembrane proteins. The most prominent substrate is the amyloid precursor protein, whose proteolytic processing leads to the production of different amyloid Aß peptides. These peptides are known to form toxic aggregates and may play a key role in Alzheimer's disease (AD). Recently, the three-dimensional structure of γ-secretase has been determined via Cryo-EM, elucidating the spatial geometry of this enzyme complex in different functional states. We have used molecular dynamics (MD) simulations to study the global dynamics and conformational transitions of γ-secretase, as well as the water and lipid distributions in and around the transmembrane domains in atomic detail. Simulations were performed on the full enzyme complex and on the membrane embedded parts alone. The simulations revealed global motions compatible with the experimental enzyme structures and indicated little dependence of the dynamics of the transmembrane domains on the soluble extracellular subunits. During the simulation on the membrane spanning part a transition between an inactive conformation (with catalytic residues far apart) toward a putatively active form (with catalytic residues in close proximity) has been observed. This conformational change is associated with a distinct rearrangement of transmembrane helices, a global compaction of the catalytically active presenilin subunit a change in the water structure near the active site and a rigidification of the protein fold. The observed conformational rearrangement allows the interpretation of the effect of several mutations on the activity of γ-secretase. A number of long-lived lipid binding sites could be identified on the membrane spanning surface of γ-secretase which may coincide with association regions of hydrophobic membrane helices to form putative substrate binding exosites.

The Binding Mode of the Sonic Hedgehog Inhibitor Robotnikinin, a Combined Docking and QM/MM MD Study.

Erroneous activation of the Hedgehog pathway has been linked to a great amount of cancerous diseases and therefore a large number of studies aiming at its inhibition have been carried out. One leverage point for novel therapeutic strategies targeting the proteins involved, is the prevention of complex formation between the extracellular signaling protein Sonic Hedgehog and the transmembrane protein Patched 1. In 2009 robotnikinin, a small molecule capable of binding to and inhibiting the activity of Sonic Hedgehog has been identified, however in the absence of X-ray structures of the Sonic Hedgehog-robotnikinin complex, the binding mode of this inhibitor remains unknown. In order to aid with the identification of novel Sonic Hedgehog inhibitors, the presented investigation elucidates the binding mode of robotnikinin by performing an extensive docking study, including subsequent molecular mechanical as well as quantum mechanical/molecular mechanical molecular dynamics simulations. The attained configurations enabled the identification of a number of key protein-ligand interactions, aiding complex formation and providing stabilizing contributions to the binding of the ligand. The predicted structure of the Sonic Hedgehog-robotnikinin complex is provided via a PDB file as Supplementary Material and can be used for further reference.

The influence of metal-ion binding on the structure and surface composition of Sonic Hedgehog: a combined classical and hybrid QM/MM MD study.

In this work, the influence of the metal ions present in vertebrate Sonic Hedgehog was assessed by a series of molecular mechanics molecular dynamics simulations with differing ionic compositions. The obtained data suggest that Ca(ii) binding has a very distinct influence on the composition of the protein surface surrounding the binding site by shaping several ionic interactions with negatively charged sidechains that otherwise would be pointing towards the solvent, repelling potential ligands. Furthermore, the Ca(ii) ions play an important role in the stability of the loop regions where they are coordinated. In contrast, the removal of the Zn(ii) ion results in no noticeable destabilization of its chemical surrounding, however, it is shown that the destabilizing effect of removed Ca(ii) ions is amplified if Zn(ii) is absent as well. Furthermore, a quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation of Sonic Hedgehog with special focus on the Zn(ii) binding site has been conducted. The results indicate that QM/MM in contrast to pure MM accurately reproduces structural features also found by experimental studies and therefore is able to provide credible predictions not only of the dynamical properties of the studied system but also of protein-ligand interactions at the metal ion binding site.

Probing the range of applicability of structure- and energy-adjusted QM/MM link bonds.

The hydrogen-capping method is one of the most popular and widely used coupling-schemes for quantum mechanics/molecular mechanics (QM/MM)-molecular dynamics simulations of macromolecular systems. This is mostly due to the fact that it is fairly convenient to implement and parametrize, thus providing an excellent compromise between accuracy and computational effort. In this work, a viable and straight-forward approach to optimize the placing of the link atom on a suitable distance ratio between the frontier atoms is discussed. To further increase the accuracy, instead of global parameters for all amino acids, different parameter sets for each type of amino acid are derived. The dependency of the link bond parameters on the chemical environment and the used QM-method is probed to assess the range of applicability of the parametrization. Suitable sets of parameters for RI-MP2, B3LYP, (RI)-B3LYP-D3, and RI-BLYP-D3 at triple-zeta level for all relevant proteinogenic amino acids are presented. Furthermore, the scope and range of the perturbation, stemming from the introduction of link bonds is evaluated through application of the presented QM/MM scheme in calculations of the active site of 15S-lipoxygenase.

Tuning the reactivity of a dissociative force field: proton transfer properties of aqueous H3O(+) and their dependence on the three-body interaction.

The proton transfer properties of the dissociative water potential developed by Garofalini et al. were closely examined by carefully analyzing the pairwise screening functions of the three-body interaction. It was shown that a simultaneous adjustment of the exponential screening factor and the three-body cutoff distance enables a selective adjustment of the diffusive properties of an excess proton, while at the same time structural and other dynamical data remain unaffected to a large extend. To investigate proton transfer properties without the influence of nuclear quantum effects, deuterated systems have been investigated in addition to their hydrogen counterparts. It was shown that the suggested parameter set A leads to significantly improved diffusion coefficients and proton hopping rates. Comparison of proton transfer correlation functions to simulation data obtained from Car-Parrinello molecular dynamics simulations confirms the improved performance of the adjusted parametrization.

Solvation properties and behaviour of lutetium(III) in aqueous solution--a quantum mechanical charge field (QMCF) study.

This work presents the first ab initio molecular dynamics study of trivalent lutetium in aqueous solution. The hybrid quantum and molecular mechanics simulation has been carried out on Hartree-Fock level and the results were compared to extended X-ray absorption fine structure and X-ray diffraction data. In addition to the structural characterisation via radial and angular distribution functions, the influence of the ion on the surrounding solvent was further investigated by local-density-corrected three-body distribution functions and frequency calculations. The obtained results for the mean Lu-O bond distance and force constant were in very good agreement with the literature. Furthermore, deeper insight into the dynamics and geometry of the solvation shell and the number of involved solvent molecules was obtained.